A copy number variation morbidity map of developmental delay

To understand the genetic heterogeneity underlying developmental delay, we compared copy number variants (CNVs) in 15,767 children with intellectual disability and various congenital defects (cases) to CNVs in 8,329 unaffected adult controls. We estimate that ～14.2% of disease in these children is caused by CNVs >400 kb. We observed a greater enrichment of CNVs in individuals with craniofacial anomalies and cardiovascular defects compared to those with epilepsy or autism. We identified 59 pathogenic CNVs, including 14 new or previously weakly supported candidates, refined the critical interval for several genomic disorders, such as the 17q21.31 microdeletion syndrome, and identified 940 candidate dosage-sensitive genes. We also developed methods to opportunistically discover small, disruptive CNVs within the large and growing diagnostic array datasets. This evolving CNV morbidity map, combined with exome and genome sequencing, will be critical for deciphering the genetic basis of developmental delay, intellectual disability and autism spectrum disorders.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

Progress and prospects in rat genetics: a community view

The rat is an important system for modeling human disease. Four years ago, the rich 150-year history of rat research was transformed by the sequencing of the rat genome, ushering in an era of exceptional opportunity for identifying genes and pathways underlying disease phenotypes. Genome-wide association studies in human populations have recently provided a direct approach for finding robust genetic associations in common diseases, but identifying the precise genes and their mechanisms of action remains problematic. In the context of significant progress in rat genomic resources over the past decade, we outline achievements in rat gene discovery to date, show how these findings have been translated to human disease, and document an increasing pace of discovery of new disease genes, pathways and mechanisms. Finally, we present a set of principles that justify continuing and strengthening genetic studies in the rat model, and further development of genomic infrastructure for rat research.     

Systematic assessment of copy number variant detection via genome-wide SNP genotyping

SNP genotyping has emerged as a technology to incorporate copy number variants (CNVs) into genetic analyses of human traits. However, the extent to which SNP platforms accurately capture CNVs remains unclear. Using independent, sequence-based CNV maps, we find that commonly used SNP platforms have limited or no probe coverage for a large fraction of CNVs. Despite this, in 9 samples we inferred 368 CNVs using Illumina SNP genotyping data and experimentally validated over two-thirds of these. We also developed a method (SNP-Conditional Mixture Modeling, SCIMM) to robustly genotype deletions using as few as two SNP probes. We find that HapMap SNPs are strongly correlated with 82% of common deletions, but the newest SNP platforms effectively tag about 50%. We conclude that currently available genome-wide SNP assays can capture CNVs accurately, but improvements in array designs, particularly in duplicated sequences, are necessary to facilitate more comprehensive analyses of genomic variation.     

Integrated detection and population-genetic analysis of SNPs and copy number variation

Dissecting the genetic basis of disease risk requires measuring all forms of genetic variation, including SNPs and copy number variants (CNVs), and is enabled by accurate maps of their locations, frequencies and population-genetic properties. We designed a hybrid genotyping array (Affymetrix SNP 6.0) to simultaneously measure 906,600 SNPs and copy number at 1.8 million genomic locations. By characterizing 270 HapMap samples, we developed a map of human CNV (at 2-kb breakpoint resolution) informed by integer genotypes for 1,320 copy number polymorphisms (CNPs) that segregate at an allele frequency >1%. More than 80% of the sequence in previously reported CNV regions fell outside our estimated CNV boundaries, indicating that large (>100 kb) CNVs affect much less of the genome than initially reported. Approximately 80% of observed copy number differences between pairs of individuals were due to common CNPs with an allele frequency >5%, and more than 99% derived from inheritance rather than new mutation. Most common, diallelic CNPs were in strong linkage disequilibrium with SNPs, and most low-frequency CNVs segregated on specific SNP haplotypes.     

Plasmodium cynomolgi genome sequences provide insight into Plasmodium vivax and the monkey malaria clade

P. cynomolgi, a malaria-causing parasite of Asian Old World monkeys, is the sister taxon of P. vivax, the most prevalent malaria-causing species in humans outside of Africa. Because P. cynomolgi shares many phenotypic, biological and genetic characteristics with P. vivax, we generated draft genome sequences for three P. cynomolgi strains and performed genomic analysis comparing them with the P. vivax genome, as well as with the genome of a third previously sequenced simian parasite, Plasmodium knowlesi. Here, we show that genomes of the monkey malaria clade can be characterized by copy-number variants (CNVs) in multigene families involved in evasion of the human immune system and invasion of host erythrocytes. We identify genome-wide SNPs, microsatellites and CNVs in the P. cynomolgi genome, providing a map of genetic variation that can be used to map parasite traits and study parasite populations. The sequencing of the P. cynomolgi genome is a critical step in developing a model system for P. vivax research and in counteracting the neglect of P. vivax.     

Genome assembly comparison identifies structural variants in the human genome

Numerous types of DNA variation exist, ranging from SNPs to larger structural alterations such as copy number variants (CNVs) and inversions. Alignment of DNA sequence from different sources has been used to identify SNPs and intermediate-sized variants (ISVs). However, only a small proportion of total heterogeneity is characterized, and little is known of the characteristics of most smaller-sized (<50 kb) variants. Here we show that genome assembly comparison is a robust approach for identification of all classes of genetic variation. Through comparison of two human assemblies (Celera's R27c compilation and the Build 35 reference sequence), we identified megabases of sequence (in the form of 13,534 putative non-SNP events) that were absent, inverted or polymorphic in one assembly. Database comparison and laboratory experimentation further demonstrated overlap or validation for 240 variable regions and confirmed >1.5 million SNPs. Some differences were simple insertions and deletions, but in regions containing CNVs, segmental duplication and repetitive DNA, they were more complex. Our results uncover substantial undescribed variation in humans, highlighting the need for comprehensive annotation strategies to fully interpret genome scanning and personalized sequencing projects.     

Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs

Accurate and complete measurement of single nucleotide (SNP) and copy number (CNV) variants, both common and rare, will be required to understand the role of genetic variation in disease. We present Birdsuite, a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus (for example, including genotypes such as A-null, AAB and BBB in addition to AA, AB and BB calls). Such genotypes more accurately depict the underlying sequence of each individual, reducing the rate of apparent mendelian inconsistencies. The Birdsuite software is applied here to data from the Affymetrix SNP 6.0 array. Additionally, we describe a method, implemented in PLINK, to utilize these combined SNP and CNV genotypes for association testing with a phenotype.     

A genetic interaction network of five genes for human polycystic kidney and liver diseases defines polycystin-1 as the central determinant of cyst formation

Autosomal dominant polycystic liver disease results from mutations in PRKCSH or SEC63. The respective gene products, glucosidase IIβ and SEC63p, function in protein translocation and quality control pathways in the endoplasmic reticulum. Here we show that glucosidase IIβ and Sec63p are required in mice for adequate expression of a functional complex of the polycystic kidney disease gene products, polycystin-1 and polycystin-2. We find that polycystin-1 is the rate-limiting component of this complex and that there is a dose-response relationship between cystic dilation and levels of functional polycystin-1 following mutation of Prkcsh or Sec63. Reduced expression of polycystin-1 also serves to sensitize the kidney to cyst formation resulting from mutations in Pkhd1, the recessive polycystic kidney disease gene. Finally, we show that proteasome inhibition increases steady-state levels of polycystin-1 in cells lacking glucosidase IIβ and that treatment with a proteasome inhibitor reduces cystic disease in orthologous gene models of human autosomal dominant polycystic liver disease.     

Genome-wide genetic changes during modern breeding of maize

The success of modern maize breeding has been demonstrated by remarkable increases in productivity over the last four decades. However, the underlying genetic changes correlated with these gains remain largely unknown. We report here the sequencing of 278 temperate maize inbred lines from different stages of breeding history, including deep resequencing of 4 lines with known pedigree information. The results show that modern breeding has introduced highly dynamic genetic changes into the maize genome. Artificial selection has affected thousands of targets, including genes and non-genic regions, leading to a reduction in nucleotide diversity and an increase in the proportion of rare alleles. Genetic changes during breeding happen rapidly, with extensive variation (SNPs, indels and copy-number variants (CNVs)) occurring, even within identity-by-descent regions. Our genome-wide assessment of genetic changes during modern maize breeding provides new strategies as well as practical targets for future crop breeding and biotechnology.     

Oral-facial-digital type I protein is required for primary cilia formation and left-right axis specification

The oral-facial-digital type I (OFD1) syndrome (OMIM 311200) is a human developmental disorder; affected individuals have craniofacial and digital abnormalities and, in 15% of cases, polycystic kidney. The disease is inherited as an X-linked dominant male-lethal trait. Using a Cre-loxP system, we generated knockout animals lacking Ofd1 and reproduced the main features of the disease, albeit with increased severity, possibly owing to differences of X inactivation patterns between human and mouse. We found failure of left-right axis specification in mutant male embryos, and ultrastructural analysis showed a lack of cilia in the embryonic node. Formation of cilia was defective in cystic kidneys from heterozygous females, implicating ciliogenesis as a mechanism underlying cyst development. In addition, we found impaired patterning of the neural tube and altered expression of the 5' Hoxa and Hoxd genes in the limb buds of mice lacking Ofd1, suggesting that Ofd1 could have a role beyond primary cilium organization and assembly.     

Using FlyAtlas to identify better Drosophila melanogaster models of human disease

FlyAtlas, a new online resource, provides the most comprehensive view yet of expression in multiple tissues of Drosophila melanogaster. Meta-analysis of the data shows that a significant fraction of the genome is expressed with great tissue specificity in the adult, demonstrating the need for the functional genomic community to embrace a wide range of functional phenotypes. Well-known developmental genes are often reused in surprising tissues in the adult, suggesting new functions. The homologs of many human genetic disease loci show selective expression in the Drosophila tissues analogous to the affected human tissues, providing a useful filter for potential candidate genes. Additionally, the contributions of each tissue to the whole-fly array signal can be calculated, demonstrating the limitations of whole-organism approaches to functional genomics and allowing modeling of a simple tissue fractionation procedure that should improve detection of weak or tissue-specific signals.     

The K-Cl cotransporter KCC3 is mutant in a severe peripheral neuropathy associated with agenesis of the corpus callosum

Peripheral neuropathy associated with agenesis of the corpus callosum (ACCPN) is a severe sensorimotor neuropathy associated with mental retardation, dysmorphic features and complete or partial agenesis of the corpus callosum. ACCPN is transmitted in an autosomal recessive fashion and is found at a high frequency in the province of Quebec, Canada. ACCPN has been previously mapped to chromosome 15q. The gene SLC12A6 (solute carrier family 12, member 6), which encodes the K+-Cl- transporter KCC3 and maps within the ACCPN candidate region, was screened for mutations in individuals with ACCPN. Four distinct protein-truncating mutations were found: two in the French Canadian population and two in non-French Canadian families. The functional consequence of the predominant French Canadian mutation (2436delG, Thr813fsX813) was examined by heterologous expression of wildtype and mutant KCC3 in Xenopus laevis oocytes; the truncated mutant is appropriately glycosylated and expressed at the cellular membrane, where it is non-functional. Mice generated with a targeted deletion of Slc12a6 have a locomotor deficit, peripheral neuropathy and a sensorimotor gating deficit, similar to the human disease. Our findings identify mutations in SLC12A6 as the genetic lesion underlying ACCPN and suggest a critical role for SLC12A6 in the development and maintenance of the nervous system.     

The European dimension for the mouse genome mutagenesis program

The European Mouse Mutagenesis Consortium is the European initiative contributing to the international effort on functional annotation of the mouse genome. Its objectives are to establish and integrate mutagenesis platforms, gene expression resources, phenotyping units, storage and distribution centers and bioinformatics resources. The combined efforts will accelerate our understanding of gene function and of human health and disease.     

Action of BTN1, the yeast orthologue of the gene mutated in Batten disease.

Neuronal ceroid-lipofuscinoses (NCL) are autosomal recessive disorders that form the most common group of progressive neurodegenerative diseases in children, with an incidence as high as 1 in 12,500 live births, and with approximately 440,000 carriers in the United States. Disease progression is characterized by a decline in mental abilities, increased severity of untreatable seizures, blindness, loss of motor skills and premature death. The CLN3 gene, which is responsible for Batten disease, has been positionally cloned. The yeast gene, denoted BTN1, encodes a non-essential protein that is 39% identical and 59% similar to human CLN3. Strains lacking Btn1p, btn1-delta, are resistant to D-(-)-threo-2-amino-1-[p-nitrophenyl]-1,3-propanediol (ANP) in a pH-dependent manner. This phenotype was complemented by expression of human CLN3, demonstrating that yeast Btn1p and human CLN3 share the same function. Here, we report that btn1-delta yeast strains have an abnormally acidic vacuolar pH in the early phases of growth. Furthermore, DNA microarray analysis of BTN1 and btn1-delta strains revealed differential expression of two genes, with at least one, HSP30, involved in pH control. Because Btn1p is located in the vacuole, we suggest that Batten disease is caused by a defect in vacuolar (lysosomal) pH control. Our findings draw parallels between fundamental biological processes in yeast and previously observed characteristics of neurodegeneration in humans.     

Leukaemia disease genes: large-scale cloning and pathway predictions

Retroviral insertional mutagenesis in BXH2 and AKXD recombinant inbred mice induces a high incidence of myeloid or B- and T-cell leukaemia and the proviral integration sites in the leukaemias provide powerful genetic tags for disease gene identification. Some of the disease genes identified by proviral tagging are also associated with human disease, validating this approach for human disease gene identification. Although many leukaemia disease genes have been identified over the years, many more remain to be cloned. Here we describe an inverse PCR (IPCR) method for proviral tagging that makes use of automated DNA sequencing and the genetic tools provided by the Mouse Genome Project, which increases the throughput for disease gene identification. We also use this IPCR method to clone and analyse more than 400 proviral integration sites from AKXD and BXH2 leukaemias and, in the process, identify more than 90 candidate disease genes. Some of these genes function in pathways already implicated in leukaemia, whereas others are likely to define new disease pathways. Our studies underscore the power of the mouse as a tool for gene discovery and functional genomics.