Jayne Y. Hehir-Kwa

Genome sequencing identifies major causes of severe intellectual disability

Severe intellectual disability (ID) occurs in 0.5% of newborns and is thought to be largely genetic in origin. The extensive genetic heterogeneity of this disorder requires a genome-wide detection of all types of genetic variation. Microarray studies and, more recently, exome sequencing have demonstrated the importance of de novo copy number variations (CNVs) and single-nucleotide variations (SNVs) in ID, but the majority of cases remain undiagnosed. Here we applied whole-genome sequencing to 50 patients with severe ID and their unaffected parents. All patients included had not received a molecular diagnosis after extensive genetic prescreening, including microarray-based CNV studies and exome sequencing. Notwithstanding this prescreening, 84 de novo SNVs affecting the coding region were identified, which showed a statistically significant enrichment of loss-of-function mutations as well as an enrichment for genes previously implicated in ID-related disorders. In addition, we identified eight de novo CNVs, including single-exon and intra-exonic deletions, as well as interchromosomal duplications. These CNVs affected known ID genes more frequently than expected. On the basis of diagnostic interpretation of all de novo variants, a conclusive genetic diagnosis was reached in 20 patients. Together with one compound heterozygous CNV causing disease in a recessive mode, this results in a diagnostic yield of 42% in this extensively studied cohort, and 62% as a cumulative estimate in an unselected cohort. These results suggest that de novo SNVs and CNVs affecting the coding region are a major cause of severe ID. Genome sequencing can be applied as a single genetic test to reliably identify and characterize the comprehensive spectrum of genetic variation, providing a genetic diagnosis in the majority of patients with severe ID.

Whole-genome sequence variation, population structure and demographic history of the Dutch population

Whole-genome sequencing enables complete characterization of genetic variation, but geographic clustering of rare alleles demands many diverse populations be studied. Here we describe the Genome of the Netherlands (GoNL) Project, in which we sequenced the whole genomes of 250 Dutch parent-offspring families and constructed a haplotype map of 20.4 million single-nucleotide variants and 1.2 million insertions and deletions. The intermediate coverage (∼13×) and trio design enabled extensive characterization of structural variation, including midsize events (30–500 bp) previously poorly catalogued and de novo mutations. We demonstrate that the quality of the haplotypes boosts imputation accuracy in independent samples, especially for lower frequency alleles. Population genetic analyses demonstrate fine-scale structure across the country and support multiple ancient migrations, consistent with historical changes in sea level and flooding. The GoNL Project illustrates how single-population whole-genome sequencing can provide detailed characterization of genetic variation and may guide the design of future population studies.

High-resolution genomic profiling of childhood ALL reveals novel recurrent genetic lesions affecting pathways involved in lymphocyte differentiation and cell cycle progression

Gross cytogenetic anomalies are traditionally being used as diagnostic, prognostic and therapeutic markers in the clinical management of cancer, including childhood acute lymphoblastic leukemia (ALL). Recently, it has become increasingly clear that genetic lesions driving tumorigenesis frequently occur at the submicroscopic level and, consequently, escape standard cytogenetic observations. Therefore, we profiled the genomes of 40 childhood ALLs at high resolution. We detected multiple de novo genetic lesions, including gross aneuploidies and segmental gains and losses, some of which were subtle and affected single genes. Many of these lesions involved recurrent (partially) overlapping deletions and duplications, containing various established leukemia-associated genes, such as ETV6, RUNX1 and MLL. Importantly, the most frequently affected genes were those controlling G1/S cell cycle progression (e.g. CDKN2A, CDKN1B and RB1), followed by genes associated with B-cell development. The latter group includes microdeletions of the B-lineage transcription factors PAX5, EBF, E2-2 and IKZF1 (Ikaros), as well as genes with other established roles in B-cell development, that is RAG1 and RAG2, FYN, PBEF1 or CBP/PAG. The fact that we frequently encountered multiple lesions affecting genes involved in cell cycle regulation and B-cell differentiation strongly suggests that both these processes need to be targeted independently and simultaneously to trigger ALL development.

Refining analyses of copy number variation identifies specific genes associated with developmental delay

Copy number variants (CNVs) are associated with many neurocognitive disorders; however, these events are typically large, and the underlying causative genes are unclear. We created an expanded CNV morbidity map from 29,085 children with developmental delay in comparison to 19,584 healthy controls, identifying 70 significant CNVs. We resequenced 26 candidate genes in 4,716 additional cases with developmental delay or autism and 2,193 controls. An integrated analysis of CNV and single-nucleotide variant (SNV) data pinpointed 10 genes enriched for putative loss of function. Follow-up of a subset of affected individuals identified new clinical subtypes of pediatric disease and the genes responsible for disease-associated CNVs. These genetic changes include haploinsufficiency of SETBP1 associated with intellectual disability and loss of expressive language and truncations of ZMYND11 in individuals with autism, aggression and complex neuropsychiatric features. This combined CNV and SNV approach facilitates the rapid discovery of new syndromes and genes involved in neuropsychiatric disease despite extensive genetic heterogeneity.

Genomic microarrays in mental retardation: a practical workflow for diagnostic applications.

Microarray‐based copy number analysis has found its way into routine clinical practice, predominantly for the diagnosis of patients with unexplained mental retardation. However, the clinical interpretation of submicroscopic copy number variants (CNVs) is complicated by the fact that many CNVs are also present in the general population. Here we introduce and discuss a workflow that can be used in routine diagnostics to assess the clinical significance of the CNVs identified. We applied this scheme to our cohort of 386 individuals with unexplained mental retardation tested using a genome‐wide tiling‐resolution DNA microarray and to 978 additional patients with mental retardation reported in 15 genome‐wide microarray studies extracted from the literature. In our cohort of 386 patients we identified 25 clinically significant copy number losses (median size 2.6 Mb), nine copy number gains (median size 2.0 Mb), and one mosaic numerical chromosome aberration. Accordingly, the overall diagnostic yield of clinically significant CNVs was 9.1%. Taken together, our cohort and the patients described in the literature include a total of 1,364 analyses of DNA copy number in which a total of 11.2% (71.9% losses, 19.6% gains, 8.5% complex) could be identified, reflecting the overall diagnostic yield of clinically significant CNVs in individuals with unexplained mental retardation. Hum Mutat 0, 1–10, 2008.

A germline homozygous mutation in the base-excision repair gene NTHL1 causes adenomatous polyposis and colorectal cancer

The genetic cause underlying the development of multiple colonic adenomas, the premalignant precursors of colorectal cancer (CRC), frequently remains unresolved in patients with adenomatous polyposis. Here we applied whole-exome sequencing to 51 individuals with multiple colonic adenomas from 48 families. In seven affected individuals from three unrelated families, we identified a homozygous germline nonsense mutation in the base-excision repair (BER) gene NTHL1. This mutation was exclusively found in a heterozygous state in controls (minor allele frequency of 0.0036; n = 2,329). All three families showed recessive inheritance of the adenomatous polyposis phenotype and progression to CRC in at least one member. All three affected women developed an endometrial malignancy or premalignancy. Genetic analysis of three carcinomas and five adenomas from different affected individuals showed a non-hypermutated profile enriched for cytosine-to-thymine transitions. We conclude that a homozygous loss-of-function germline mutation in the NTHL1 gene predisposes to a new subtype of BER-associated adenomatous polyposis and CRC.

The Genome of the Netherlands: design, and project goals

Within the Netherlands a national network of biobanks has been established (Biobanking and Biomolecular Research Infrastructure-Netherlands (BBMRI-NL)) as a national node of the European BBMRI. One of the aims of BBMRI-NL is to enrich biobanks with different types of molecular and phenotype data. Here, we describe the Genome of the Netherlands (GoNL), one of the projects within BBMRI-NL. GoNL is a whole-genome-sequencing project in a representative sample consisting of 250 trio-families from all provinces in the Netherlands, which aims to characterize DNA sequence variation in the Dutch population. The parent–offspring trios include adult individuals ranging in age from 19 to 87 years (mean=53 years; SD=16 years) from birth cohorts 1910–1994. Sequencing was done on blood-derived DNA from uncultured cells and accomplished coverage was 14–15x. The family-based design represents a unique resource to assess the frequency of regional variants, accurately reconstruct haplotypes by family-based phasing, characterize short indels and complex structural variants, and establish the rate of de novo mutational events. GoNL will also serve as a reference panel for imputation in the available genome-wide association studies in Dutch and other cohorts to refine association signals and uncover population-specific variants. GoNL will create a catalog of human genetic variation in this sample that is uniquely characterized with respect to micro-geographic location and a wide range of phenotypes. The resource will be made available to the research and medical community to guide the interpretation of sequencing projects. The present paper summarizes the global characteristics of the project.

Identification of clinically significant, submicroscopic chromosome alterations and UPD in fetuses with ultrasound anomalies using genome-wide 250k SNP array analysis

Background The implementation of microarray analysis in prenatal diagnostics is a topic of discussion, as rare copy number variants with unknown/uncertain clinical consequences are likely to be found. The application of targeted microarrays limits such findings, but the potential disadvantage is that relevant, so far unknown, aberrations might be overlooked. Therefore, we explore the possibilities for the prenatal application of the genome-wide 250k single nucleotide polymorphism array platform. Methods Affymetrix 250k NspI single nucleotide polymorphism array analysis (Affymetrix, Inc., Santa Clara, California, USA) was performed on DNA from 38 prenatally karyotyped fetuses with ultrasound anomalies. Analyses were performed after termination of pregnancy, intrauterine fetal death or birth on DNA isolated from fetal or neonatal material. Results Aberrations were detected in 17 of 38 fetuses, 6 of whom with a previously identified chromosomal abnormality and 11 with previously normal or balanced karyotypes. Of the latter, the detected aberration occurred de novo and was considered of clinical relevance in five cases (16%), inherited from a healthy parent in four cases (12%), and de novo yet with unclear clinical relevance in two cases (6%). The clinically relevant abnormalities either were novel copy number variants (n=3) or concerned a uniparental disomy (n=2). Conclusion In at least 16% of fetuses with ultrasound anomalies and a normal or balanced karyotype, causal (submicroscopic) aberrations were detected, illustrating the importance of the (careful) implementation of microarray analysis in prenatal diagnosis. The fact that the identified, clinically relevant, aberrations would have gone undetected with most targeted approaches underscores the added value of a genome-wide approach.

Genome-wide copy number profiling on high-density bacterial artificial chromosomes, single-nucleotide polymorphisms, and oligonucleotide microarrays: a platform comparison based on statistical power analysis.

Abstract Recently, comparative genomic hybridization onto bacterial artificial chromosome (BAC) arrays (array-based comparative genomic hybridization) has proved to be successful for the detection of submicroscopic DNA copy-number variations in health and disease. Technological improvements to achieve a higher resolution have resulted in the generation of additional microarray platforms encompassing larger numbers of shorter DNA targets (oligonucleotides). Here, we present a novel method to estimate the ability of a microarray to detect genomic copy-number variations of different sizes and types (i.e. deletions or duplications). We applied our method, which is based on statistical power analysis, to four widely used high-density genomic microarray platforms. By doing so, we found that the high-density oligonucleotide platforms are superior to the BAC platform for the genome-wide detection of copy-number variations smaller than 1 Mb. The capacity to reliably detect single copy-number variations below 100 kb, however, appeared to be limited for all platforms tested. In addition, our analysis revealed an unexpected platform-dependent difference in sensitivity to detect a single copy-number loss and a single copy-number gain. These analyses provide a first objective insight into the true capacities and limitations of different genomic microarrays to detect and define DNA copy-number variations.

Clinical significance of de novo and inherited copy-number variation.

Copy‐number variations (CNVs) are a common cause of intellectual disability and/or multiple congenital anomalies (ID/MCA). However, the clinical interpretation of CNVs remains challenging, especially for inherited CNVs. Well‐phenotyped patients (5,531) with ID/MCA were screened for rare CNVs using a 250K single‐nucleotide polymorphism array platform in order to improve the understanding of the contribution of CNVs to a patients phenotype. We detected 1,663 rare CNVs in 1,388 patients (25.1%; range 0–5 per patient) of which 437 occurred de novo and 638 were inherited. The detected CNVs were analyzed for various characteristics, gene content, and genotype–phenotype correlations. Patients with severe phenotypes, including organ malformations, had more de novo CNVs (P < 0.001), whereas patient groups with milder phenotypes, such as facial dysmorphisms, were enriched for both de novo and inherited CNVs (P < 0.001), indicating that not only de novo but also inherited CNVs can be associated with a clinically relevant phenotype. Moreover, patients with multiple CNVs presented with a more severe phenotype than patients with a single CNV (P < 0.001), pointing to a combinatorial effect of the additional CNVs. In addition, we identified 20 de novo single‐gene CNVs that directly indicate novel genes for ID/MCA, including ZFHX4, ANKH, DLG2, MPP7, CEP89, TRIO, ASTN2, and PIK3C3.