Richard Redon

Global variation in copy number in the human genome

Copy number variation (CNV) of DNA sequences is functionally significant but has yet to be fully ascertained. We have constructed a first-generation CNV map of the human genome through the study of 270 individuals from four populations with ancestry in Europe, Africa or Asia (the HapMap collection). DNA from these individuals was screened for CNV using two complementary technologies: single-nucleotide polymorphism (SNP) genotyping arrays, and clone-based comparative genomic hybridization. A total of 1,447 copy number variable regions (CNVRs), which can encompass overlapping or adjacent gains or losses, covering 360 megabases (12% of the genome) were identified in these populations. These CNVRs contained hundreds of genes, disease loci, functional elements and segmental duplications. Notably, the CNVRs encompassed more nucleotide content per genome than SNPs, underscoring the importance of CNV in genetic diversity and evolution. The data obtained delineate linkage disequilibrium patterns for many CNVs, and reveal marked variation in copy number among populations. We also demonstrate the utility of this resource for genetic disease studies.

Origins and functional impact of copy number variation in the human genome

Structural variations of DNA greater than 1 kilobase in size account for most bases that vary among human genomes, but are still relatively under-ascertained. Here we use tiling oligonucleotide microarrays, comprising 42 million probes, to generate a comprehensive map of 11,700 copy number variations (CNVs) greater than 443 base pairs, of which most (8,599) have been validated independently. For 4,978 of these CNVs, we generated reference genotypes from 450 individuals of European, African or East Asian ancestry. The predominant mutational mechanisms differ among CNV size classes. Retrotransposition has duplicated and inserted some coding and non-coding DNA segments randomly around the genome. Furthermore, by correlation with known trait-associated single nucleotide polymorphisms (SNPs), we identified 30 loci with CNVs that are candidates for influencing disease susceptibility. Despite this, having assessed the completeness of our map and the patterns of linkage disequilibrium between CNVs and SNPs, we conclude that, for complex traits, the heritability void left by genome-wide association studies will not be accounted for by common CNVs.

Diet and the evolution of human amylase gene copy number variation

Starch consumption is a prominent characteristic of agricultural societies and hunter-gatherers in arid environments. In contrast, rainforest and circum-arctic hunter-gatherers and some pastoralists consume much less starch. This behavioral variation raises the possibility that different selective pressures have acted on amylase, the enzyme responsible for starch hydrolysis. We found that copy number of the salivary amylase gene (AMY1) is correlated positively with salivary amylase protein level and that individuals from populations with high-starch diets have, on average, more AMY1 copies than those with traditionally low-starch diets. Comparisons with other loci in a subset of these populations suggest that the extent of AMY1 copy number differentiation is highly unusual. This example of positive selection on a copy number–variable gene is, to our knowledge, one of the first discovered in the human genome. Higher AMY1 copy numbers and protein levels probably improve the digestion of starchy foods and may buffer against the fitness-reducing effects of intestinal disease.

Microarray based comparative genomic hybridisation (array-CGH) detects submicroscopic chromosomal deletions and duplications in patients with learning disability/mental retardation and dysmorphic features

The underlying causes of learning disability and dysmorphic features in many patients remain unidentified despite extensive investigation. Routine karyotype analysis is not sensitive enough to detect subtle chromosome rearrangements (less than 5 Mb). The presence of subtle DNA copy number changes was investigated by array-CGH in 50 patients with learning disability and dysmorphism, employing a DNA microarray constructed from large insert clones spaced at approximately 1 Mb intervals across the genome. Twelve copy number abnormalities were identified in 12 patients (24% of the total): seven deletions (six apparently de novo and one inherited from a phenotypically normal parent) and five duplications (one de novo and four inherited from phenotypically normal parents). Altered segments ranged in size from those involving a single clone to regions as large as 14 Mb. No recurrent deletion or duplication was identified within this cohort of patients. On the basis of these results, we anticipate that array-CGH will become a routine method of genome-wide screening for imbalanced rearrangements in children with learning disability.

Array-based comparative genomic hybridisation identifies high frequency of cryptic chromosomal rearrangements in patients with syndromic autism spectrum disorders

Background: Autism spectrum disorders (ASD) refer to a broader group of neurobiological conditions, pervasive developmental disorders. They are characterised by a symptomatic triad associated with qualitative changes in social interactions, defect in communication abilities, and repetitive and stereotyped interests and activities. ASD is prevalent in 1 to 3 per 1000 people. Despite several arguments for a strong genetic contribution, the molecular basis of a most cases remains unexplained. About 5% of patients with autism have a chromosome abnormality visible with cytogenetic methods. The most frequent are 15q11–q13 duplication, 2q37 and 22q13.3 deletions. Many other chromosomal imbalances have been described. However, most of them remain undetectable using routine karyotype analysis, thus impeding diagnosis and genetic counselling. Methods and results: 29 patients presenting with syndromic ASD were investigated using a DNA microarray constructed from large insert clones spaced at approximately 1 Mb intervals across the genome. Eight clinically relevant rearrangements were identified in 8 (27.5%) patients: six deletions and two duplications. Altered segments ranged in size from 1.4 to 16 Mb (2–19 clones). No recurrent abnormality was identified. Conclusion: These results clearly show that array comparative genomic hybridisation should be considered to be an essential aspect of the genetic analysis of patients with syndromic ASD. Moreover, besides their importance for diagnosis and genetic counselling, they may allow the delineation of new contiguous gene syndromes associated with ASD. Finally, the detailed molecular analysis of the rearranged regions may pave the way for the identification of new ASD genes.

Common variants at SCN5A-SCN10A and HEY2 are associated with Brugada syndrome, a rare disease with high risk of sudden cardiac death

Brugada syndrome is a rare cardiac arrhythmia disorder, causally related to SCN5A mutations in around 20% of cases. Through a genome-wide association study of 312 individuals with Brugada syndrome and 1,115 controls, we detected 2 significant association signals at the SCN10A locus (rs10428132) and near the HEY2 gene (rs9388451). Independent replication confirmed both signals (meta-analyses: rs10428132, P = 1.0 × 10−68; rs9388451, P = 5.1 × 10−17) and identified one additional signal in SCN5A (at 3p21; rs11708996, P = 1.0 × 10−14). The cumulative effect of the three loci on disease susceptibility was unexpectedly large (Ptrend = 6.1 × 10−81). The association signals at SCN5A-SCN10A demonstrate that genetic polymorphisms modulating cardiac conduction can also influence susceptibility to cardiac arrhythmia. The implication of association with HEY2, supported by new evidence that Hey2 regulates cardiac electrical activity, shows that Brugada syndrome may originate from altered transcriptional programming during cardiac development. Altogether, our findings indicate that common genetic variation can have a strong impact on the predisposition to rare diseases.

Copy number variation and evolution in humans and chimpanzees

Copy number variants (CNVs) underlie many aspects of human phenotypic diversity and provide the raw material for gene duplication and gene family expansion. However, our understanding of their evolutionary significance remains limited. We performed comparative genomic hybridization on a single human microarray platform to identify CNVs among the genomes of 30 humans and 30 chimpanzees as well as fixed copy number differences between species. We found that human and chimpanzee CNVs occur in orthologous genomic regions far more often than expected by chance and are strongly associated with the presence of highly homologous intrachromosomal segmental duplications. By adapting population genetic analyses for use with copy number data, we identified functional categories of genes that have likely evolved under purifying or positive selection for copy number changes. In particular, duplications and deletions of genes with inflammatory response and cell proliferation functions may have been fixed by positive selection and involved in the adaptive phenotypic differentiation of humans and chimpanzees.

A robust statistical method for case-control association testing with copy number variation

Copy number variation (CNV) is pervasive in the human genome and can play a causal role in genetic diseases. The functional impact of CNV cannot be fully captured through linkage disequilibrium with SNPs. These observations motivate the development of statistical methods for performing direct CNV association studies. We show through simulation that current tests for CNV association are prone to false-positive associations in the presence of differential errors between cases and controls, especially if quantitative CNV measurements are noisy. We present a statistical framework for performing case-control CNV association studies that applies likelihood ratio testing of quantitative CNV measurements in cases and controls. We show that our methods are robust to differential errors and noisy data and can achieve maximal theoretical power. We illustrate the power of these methods for testing for association with binary and quantitative traits, and have made this software available as the R package CNVtools.

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 (Celeras 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.

Truncating mutations in the last exon of NOTCH2 cause a rare skeletal disorder with osteoporosis

Hajdu-Cheney syndrome is a rare autosomal dominant skeletal disorder with facial anomalies, osteoporosis and acro-osteolysis. We sequenced the exomes of six unrelated individuals with this syndrome and identified heterozygous nonsense and frameshift mutations in NOTCH2 in five of them. All mutations cluster to the last coding exon of the gene, suggesting that the mutant mRNA products escape nonsense-mediated decay and that the resulting truncated NOTCH2 proteins act in a gain-of-function manner.