Xavier Estivill

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.

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project—the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

Mutations in the cystic fibrosis gene in patients with congenital absence of the vas deferens

BACKGROUND Congenital bilateral absence of the vas deferens (CBAVD) is a form of male infertility in which mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene have been identified. The molecular basis of CBAVD is not completely understood. Although patients with cystic fibrosis have mutations in both copies of the CFTR gene, most patients with CBAVD have mutations in only one copy of the gene. METHODS To investigate CBAVD at the molecular level, we have characterized the mutations in the CFTR gene in 102 patients with this condition. None had clinical manifestations of cystic fibrosis. We also analyzed a DNA variant (the 5T allele) in a noncoding region of CFTR that causes reduced levels of the normal CFTR protein. Parents of patients with cystic fibrosis, patients with types of infertility other than CBAVD, and normal subjects were studied as controls. RESULTS Nineteen of the 102 patients with CBAVD had mutations in both copies of the CFTR gene, and none of them had the 5T allele. Fifty-four patients had a mutation in one copy of CFTR, and 34 of them (63 percent) had the 5T allele in the other CFTR gene. In 29 patients no CFTR mutations were found, but 7 of them (24 percent) had the 5T allele. In contrast, the frequency of this allele in the general population was about 5 percent. CONCLUSIONS Most patients with CBAVD have mutations in the CFTR gene. The combination of the 5T allele in one copy of the CFTR gene with a cystic fibrosis mutation in the other copy is the most common cause of CBAVD: The 5T allele mutation has a wide range of clinical presentations, occurring in patients with CBAVD or moderate forms of cystic fibrosis and in fertile men.

Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.

A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertook a genome-wide association study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified associations at eight previously unreported genomic loci. Seven loci harbored genes with recognized immune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These associations were replicated in 9,079 European samples (six loci with a combined P < 5 × 10−8 and two loci with a combined P < 5 × 10−7). We also report compelling evidence for an interaction between the HLA-C and ERAP1 loci (combined P = 6.95 × 10−6). ERAP1 plays an important role in MHC class I peptide processing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk allele. Our findings implicate pathways that integrate epidermal barrier dysfunction with innate and adaptive immune dysregulation in psoriasis pathogenesis.

Connexin-26 mutations in sporadic and inherited sensorineural deafness

BACKGROUND Hearing impairment affects one infant in 1000 and 4% of people aged younger than 45 years. Congenital deafness is inherited or apparently sporadic. We have shown previously that DFNB1 on chromosome 13 is a major locus for recessive deafness in about 80% of Mediterranean families and that the connexin-26 gene gap junction protein beta2 (GJB2) is mutated in DFNB1 families. We investigated mutations in the GJB2 gene in familial and sporadic cases of deafness. METHODS We obtained DNA samples from 82 families from Italy and Spain with recessive non-syndromic deafness and from 54 unrelated participants with apparently sporadic congenital deafness. We analysed the coding region of the GJB2 gene for mutations. We also tested 280 unrelated people from the general populations of Italy and Spain for the frameshift mutation 35delG. FINDINGS 49% of participants with recessive deafness and 37% of sporadic cases had mutations in the GJB2 gene. The 35delG mutation accounted for 85% of GJB2 mutations, six other mutations accounted for 6% of alleles, and no changes in the coding region of GJB2 were detected in 9% of DFNB1 alleles. The carrier frequency of mutation 35delG among people from the general population was one in 31 (95% CI one in 19 to one in 87). INTERPRETATION Mutations in the GJB2 gene are a major cause of inherited and apparently sporadic congenital deafness. Mutation 35delG is the most common mutation for sensorineural deafness. Identification of 35delG and other mutations in the GJB2 gene should facilitate diagnosis and counselling for the most common genetic form of deafness.

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.

Dating the Origin of the CCR5-Δ32 AIDS-Resistance Allele by the Coalescence of Haplotypes

The CCR5-Delta32 deletion obliterates the CCR5 chemokine and the human immunodeficiency virus (HIV)-1 coreceptor on lymphoid cells, leading to strong resistance against HIV-1 infection and AIDS. A genotype survey of 4,166 individuals revealed a cline of CCR5-Delta32 allele frequencies of 0%-14% across Eurasia, whereas the variant is absent among native African, American Indian, and East Asian ethnic groups. Haplotype analysis of 192 Caucasian chromosomes revealed strong linkage disequilibrium between CCR5 and two microsatellite loci. By use of coalescence theory to interpret modern haplotype genealogy, we estimate the origin of the CCR5-Delta32-containing ancestral haplotype to be approximately 700 years ago, with an estimated range of 275-1,875 years. The geographic cline of CCR5-Delta32 frequencies and its recent emergence are consistent with a historic strong selective event (e.g. , an epidemic of a pathogen that, like HIV-1, utilizes CCR5), driving its frequency upward in ancestral Caucasian populations.