Lluís Armengol

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.

Recurrent rearrangements of chromosome 1q21.1 and variable pediatric phenotypes

BACKGROUND Duplications and deletions in the human genome can cause disease or predispose persons to disease. Advances in technologies to detect these changes allow for the routine identification of submicroscopic imbalances in large numbers of patients. METHODS We tested for the presence of microdeletions and microduplications at a specific region of chromosome 1q21.1 in two groups of patients with unexplained mental retardation, autism, or congenital anomalies and in unaffected persons. RESULTS We identified 25 persons with a recurrent 1.35-Mb deletion within 1q21.1 from screening 5218 patients. The microdeletions had arisen de novo in eight patients, were inherited from a mildly affected parent in three patients, were inherited from an apparently unaffected parent in six patients, and were of unknown inheritance in eight patients. The deletion was absent in a series of 4737 control persons (P=1.1x10(-7)). We found considerable variability in the level of phenotypic expression of the microdeletion; phenotypes included mild-to-moderate mental retardation, microcephaly, cardiac abnormalities, and cataracts. The reciprocal duplication was enriched in nine children with mental retardation or autism spectrum disorder and other variable features (P=0.02). We identified three deletions and three duplications of the 1q21.1 region in an independent sample of 788 patients with mental retardation and congenital anomalies. CONCLUSIONS We have identified recurrent molecular lesions that elude syndromic classification and whose disease manifestations must be considered in a broader context of development as opposed to being assigned to a specific disease. Clinical diagnosis in patients with these lesions may be most readily achieved on the basis of genotype rather than phenotype.

SNPassoc: an R package to perform whole genome association studies

UNLABELLED The popularization of large-scale genotyping projects has led to the widespread adoption of genetic association studies as the tool of choice in the search for single nucleotide polymorphisms (SNPs) underlying susceptibility to complex diseases. Although the analysis of individual SNPs is a relatively trivial task, when the number is large and multiple genetic models need to be explored it becomes necessary a tool to automate the analyses. In order to address this issue, we developed SNPassoc, an R package to carry out most common analyses in whole genome association studies. These analyses include descriptive statistics and exploratory analysis of missing values, calculation of Hardy-Weinberg equilibrium, analysis of association based on generalized linear models (either for quantitative or binary traits), and analysis of multiple SNPs (haplotype and epistasis analysis). AVAILABILITY Package SNPassoc is available at CRAN from http://cran.r-project.org. SUPPLEMENTARY INFORMATION A tutorial is available on Bioinformatics online and in http://davinci.crg.es/estivill_lab/snpassoc.

Deletion of the late cornified envelope LCE3B and LCE3C genes as a susceptibility factor for psoriasis.

Psoriasis is a common inflammatory skin disease with a prevalence of 2–3% in individuals of European ancestry. In a genome-wide search for copy number variants (CNV) using a sample pooling approach, we have identified a deletion comprising LCE3B and LCE3C, members of the late cornified envelope (LCE) gene cluster. The absence of LCE3B and LCE3C (LCE3C_LCE3B-del) is significantly associated (P = 1.38E–08) with risk of psoriasis in 2,831 samples from Spain, The Netherlands, Italy and the United States, and in a family-based study (P = 5.4E–04). LCE3C_LCE3B-del is tagged by rs4112788 (r 2 = 0.93), which is also strongly associated with psoriasis (P < 6.6E–09). LCE3C_LCE3B-del shows epistatic effects with the HLA-Cw6 allele on the development of psoriasis in Dutch samples and multiplicative effects in the other samples. LCE expression can be induced in normal epidermis by skin barrier disruption and is strongly expressed in psoriatic lesions, suggesting that compromised skin barrier function has a role in psoriasis susceptibility.

Copy Number Variants and Common Disorders: Filling the Gaps and Exploring Complexity in Genome-Wide Association Studies

Genome-wide association scans (GWASs) using single nucleotide polymorphisms (SNPs) have been completed successfully for several common disorders and have detected over 30 new associations. Considering the large sample sizes and genome-wide SNP coverage of the scans, one might have expected many of the common variants underpinning the genetic component of various disorders to have been identified by now. However, these studies have not evaluated the contribution of other forms of genetic variation, such as structural variation, mainly in the form of copy number variants (CNVs). Known CNVs account for over 15% of the assembled human genome sequence. Since CNVs are not easily tagged by SNPs, might have a wide range of copy number variability, and often fall in genomic regions not well covered by whole-genome arrays or not genotyped by the HapMap project, current GWASs have largely missed the contribution of CNVs to complex disorders. In fact, some CNVs have already been reported to show association with several complex disorders using candidate gene/region approaches, underpinning the importance of regions not investigated in current GWASs. This reveals the need for new generation arrays (some already in the market) and the use of tailored approaches to explore the full dimension of genome variability beyond the single nucleotide scale.

A Polymorphic Genomic Duplication on Human Chromosome 15 Is a Susceptibility Factor for Panic and Phobic Disorders

Anxiety disorders are complex and common psychiatric illnesses associated with considerable morbidity and social cost. We have studied the molecular basis of the cooccurrence of panic and phobic disorders with joint laxity. We have identified an interstitial duplication of human chromosome 15q24-26 (named DUP25), which is significantly associated with panic/agoraphobia/social phobia/joint laxity in families, and with panic disorder in nonfamilial cases. Mosaicism, different forms of DUP25 within the same family, and absence of segregation of 15q24-26 markers with DUP25 and the psychiatric phenotypes suggest a non-Mendelian mechanism of disease-causing mutation. We propose that DUP25, which is present in 7% control subjects, is a susceptibility factor for a clinical phenotype that includes panic and phobic disorders and joint laxity.

Met66 in the brain-derived neurotrophic factor (BDNF) precursor is associated with anorexia nervosa restrictive type

Several lines of evidence support a role for brain-derived neurotrophic factor (BDNF) alterations in the etiology of eating disorders (EDs). BDNF heterozygous knockout mice show alterations in eating behavior, increased body weight and adipocyte hypertrophy. BDNF also regulates the synaptic efficiency through the modulation of key neurotransmitter systems previously known to be involved in ED. These findings, together with the fact that this neurotrophin is expressed in the hypothalamus nuclei associated with weight regulation and feeding control, led us to propose BDNF as a candidate gene for ED. To investigate the possible involvement of this neurotrophin in eating behavior, we screened the BDNF gene in 95 ED patients and identified four sequence variants. Two of them, −374A/T and −256G/A, were found in two patients with anorexia nervosa (AN) and consisted of single-nucleotide mutations within the 5′ untranslated region (5′UTR). The other two polymorphisms resulted in a C to T transition located at the 5′UTR of the BDNF gene and an amino-acid substitution within the BDNF precursor protein (Val66Met). We performed a case–control study for these two Single-nucleotide polymorphisms in a sample of 143 ED patients and 112 unrelated controls and found a strong association of restricting AN (ANR) with the Met allele of the Val66Met BDNF polymorphism (2p=0.002). There was also evidence for a significant effect of this sequence variant on the minimum body mass index (MBMI) (2p=0.006). These results suggest that the BDNF Met66 variant may be a susceptibility factor to ED, mainly to ANR and low MBMI.

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.

Uniparental disomies, homozygous deletions, amplifications, and target genes in mantle cell lymphoma revealed by integrative high-resolution whole-genome profiling

Mantle cell lymphoma (MCL) is genetically characterized by the t(11;14)(q13;q32) translocation and a high number of secondary chromosomal alterations. However, only a limited number of target genes have been identified. We have studied 10 MCL cell lines and 28 primary tumors with a combination of a high-density single-nucleotide polymorphism array and gene expression profiling. We detected highly altered genomes in the majority of the samples with a high number of partial uniparental disomies (UPDs). The UPD at 17p was one of the most common, and it was associated with TP53 gene inactivation. Homozygous deletions targeted 4 known tumor suppressor genes (CDKN2C, BCL2L11, CDKN2A, and RB1) and 6 new genes (FAF1, MAP2, SP100, MOBKL2B, ZNF280A, and PRAME). Gene amplification coupled with overexpression was identified in 35 different regions. The most recurrent amplified regions were 11q13.3-q13.5, 13q31.3, and 18q21.33, which targeted CCND1, C13orf25, and BCL2, respectively. Interestingly, the breakpoints flanking all the genomic alterations, including UPDs, were significantly associated with genomic regions enriched in copy number variants and segmental duplications, suggesting that the recombination at these regions may play a role in the genomic instability of MCL. This integrative genomic analysis has revealed target genes that may be potentially relevant in MCL pathogenesis.

Mosaic Uniparental Disomies and Aneuploidies as Large Structural Variants of the Human Genome

Mosaicism is defined as the coexistence of cells with different genetic composition within an individual, caused by postzygotic somatic mutation. Although somatic mosaicism for chromosomal abnormalities is a well-established cause of developmental and somatic disorders and has also been detected in different tissues, its frequency and extent in the adult normal population are still unknown. We provide here a genome-wide survey of mosaic genomic variation obtained by analyzing Illumina 1M SNP array data from blood or buccal DNA samples of 1991 adult individuals from the Spanish Bladder Cancer/EPICURO genome-wide association study. We found mosaic abnormalities in autosomes in 1.7% of samples, including 23 segmental uniparental disomies, 8 complete trisomies, and 11 large (1.5-37 Mb) copy-number variants. Alterations were observed across the different autosomes with recurrent events in chromosomes 9 and 20. No case-control differences were found in the frequency of events or the percentage of cells affected, thus indicating that most rearrangements found are not central to the development of bladder cancer. However, five out of six events tested were detected in both blood and bladder tissue from the same individual, indicating an early developmental origin. The high cellular frequency of the anomalies detected and their presence in normal adult individuals suggest that this type of mosaicism is a widespread phenomenon in the human genome. Somatic mosaicism should be considered in the expanding repertoire of inter- and intraindividual genetic variation, some of which may cause somatic human diseases but also contribute to modifying inherited disorders and/or late-onset multifactorial traits.