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.

Detection of clinically relevant copy-number variants by exome sequencing in a large cohort of genetic disorders.

PURPOSE: Copy-number variation is a common source of genomic variation and an important genetic cause of disease. Microarray-based analysis of copy-number variants (CNVs) has become a first-tier diagnostic test for patients with neurodevelopmental disorders, with a diagnostic yield of 10-20%. However, for most other genetic disorders, the role of CNVs is less clear and most diagnostic genetic studies are generally limited to the study of single-nucleotide variants (SNVs) and other small variants. With the introduction of exome and genome sequencing, it is now possible to detect both SNVs and CNVs using an exome- or genome-wide approach with a single test. METHODS: We performed exome-based read-depth CNV screening on data from 2,603 patients affected by a range of genetic disorders for which exome sequencing was performed in a diagnostic setting. RESULTS: In total, 123 clinically relevant CNVs ranging in size from 727 bp to 15.3 Mb were detected, which resulted in 51 conclusive diagnoses and an overall increase in diagnostic yield of ~2% (ranging from 0 to -5.8% per disorder). CONCLUSIONS: This study shows that CNVs play an important role in a broad range of genetic disorders and that detection via exome-based CNV profiling results in an increase in the diagnostic yield without additional testing, bringing us closer to single-test genomics.Genet Med advance online publication 27 October 2016.

Mutations in DDHD2, encoding an intracellular phospholipase A(1), cause a recessive form of complex hereditary spastic paraplegia.

We report on four families affected by a clinical presentation of complex hereditary spastic paraplegia (HSP) due to recessive mutations in DDHD2, encoding one of the three mammalian intracellular phospholipases A(1) (iPLA(1)). The core phenotype of this HSP syndrome consists of very early-onset (<2 years) spastic paraplegia, intellectual disability, and a specific pattern of brain abnormalities on cerebral imaging. An essential role for DDHD2 in the human CNS, and perhaps more specifically in synaptic functioning, is supported by a reduced number of active zones at synaptic terminals in Ddhd-knockdown Drosophila models. All identified mutations affect the protein's DDHD domain, which is vital for its phospholipase activity. In line with the function of DDHD2 in lipid metabolism and its role in the CNS, an abnormal lipid peak indicating accumulation of lipids was detected with cerebral magnetic resonance spectroscopy, which provides an applicable diagnostic biomarker that can distinguish the DDHD2 phenotype from other complex HSP phenotypes. We show that mutations in DDHD2 cause a specific complex HSP subtype (SPG54), thereby linking a member of the PLA(1) family to human neurologic disease.

Identification of pathogenic gene variants in small families with intellectually disabled siblings by exome sequencing.

BACKGROUND: Intellectual disability (ID) is a common neurodevelopmental disorder affecting 1-3% of the general population. Mutations in more than 10% of all human genes are considered to be involved in this disorder, although the majority of these genes are still unknown. OBJECTIVES: We investigated 19 small non-consanguineous families with two to five affected siblings in order to identify pathogenic gene variants in known, novel and potential ID candidate genes. Non-consanguineous families have been largely ignored in gene identification studies as small family size precludes prior mapping of the genetic defect. METHODS AND RESULTS: Using exome sequencing, we identified pathogenic mutations in three genes, DDHD2, SLC6A8, and SLC9A6, of which the latter two have previously been implicated in X-linked ID phenotypes. In addition, we identified potentially pathogenic mutations in BCORL1 on the X-chromosome and in MCM3AP, PTPRT, SYNE1, and ZNF528 on autosomes. CONCLUSIONS: We show that potentially pathogenic gene variants can be identified in small, non-consanguineous families with as few as two affected siblings, thus emphasising their value in the identification of syndromic and non-syndromic ID genes.

Heterozygous mutations of FREM1 are associated with an increased risk of isolated metopic craniosynostosis in humans and mice.

The premature fusion of the paired frontal bones results in metopic craniosynostosis (MC) and gives rise to the clinical phenotype of trigonocephaly. Deletions of chromosome 9p22.3 are well described as a cause of MC with variably penetrant midface hypoplasia. In order to identify the gene responsible for the trigonocephaly component of the 9p22.3 syndrome, a cohort of 109 patients were assessed by high-resolution arrays and MLPA for copy number variations (CNVs) involving 9p22. Five CNVs involving FREM1, all of which were de novo variants, were identified by array-based analyses. The remaining 104 patients with MC were then subjected to targeted FREM1 gene re-sequencing, which identified 3 further mutant alleles, one of which was de novo. Consistent with a pathogenic role, mouse Frem1 mRNA and protein expression was demonstrated in the metopic suture as well as in the pericranium and dura mater. Micro-computed tomography based analyses of the mouse posterior frontal (PF) suture, the human metopic suture equivalent, revealed advanced fusion in all mice homozygous for either of two different Frem1 mutant alleles, while heterozygotes exhibited variably penetrant PF suture anomalies. Gene dosage-related penetrance of midfacial hypoplasia was also evident in the Frem1 mutants. These data suggest that CNVs and mutations involving FREM1 can be identified in a significant percentage of people with MC with or without midface hypoplasia. Furthermore, we present Frem1 mutant mice as the first bona fide mouse model of human metopic craniosynostosis and a new model for midfacial hypoplasia.

Mutations in a new member of the chromodomain gene family cause CHARGE syndrome.

CHARGE syndrome is a common cause of congenital anomalies affecting several tissues in a nonrandom fashion. We report a 2.3-Mb de novo overlapping microdeletion on chromosome 8q12 identified by array comparative genomic hybridization in two individuals with CHARGE syndrome. Sequence analysis of genes located in this region detected mutations in the gene CHD7 in 10 of 17 individuals with CHARGE syndrome without microdeletions, accounting for the disease in most affected individuals.

Rare pathogenic microdeletions and tandem duplications are microhomology-mediated and stimulated by local genomic architecture.

Genomic copy number variation (CNV) plays a major role in various human diseases as well as in normal phenotypic variability. For some recurrent disease-causing CNVs that convey genomic disorders, the causative mechanism is meiotic, non-allelic, homologous recombination between breakpoint regions exhibiting extensive sequence homology (e.g. low-copy repeats). For the majority of recently identified rare pathogenic CNVs, however, the mechanism is unknown. Recently, a model for CNV formation implicated mitotic replication-based mechanisms, such as (alternative) non-homologous end joining and fork stalling and template switching, in the etiology of human pathogenic CNVs. The extent to which such mitotic mechanisms contribute to rare pathogenic CNVs remains to be determined. In addition, it is unexplored whether genomic architectural features such as repetitive elements or sequence motifs associated with DNA breakage stimulate the formation of rare pathogenic CNVs. To this end, we have sequenced breakpoint junctions of 30 rare pathogenic microdeletions and eight tandem duplications, representing the largest series of such CNVs examined to date in this much detail. Our results demonstrate the presence of (micro)homology ranging from 2 to over 75 bp, in 79% of the breakpoint junctions. This indicates that microhomology-mediated repair mechanisms, including the recently reported fork stalling and template switching and/or microhomology-mediated break-induced replication, prevail in rare pathogenic CNVs. In addition, we found that the vast majority of all breakpoints (81%) were associated with at least one of the genomic architectural features evaluated. Moreover, 75% of tandem duplication breakpoints were associated with the presence of one of two novel sequence motifs. These data suggest that rare pathogenic microdeletions and tandem duplications do not occur at random genome sequences, but are stimulated and potentially catalyzed by various genomic architectural features.

Effects of developmental plasticity on heat tolerance may be mediated by changes in cell size in Drosophila melanogaster.

There is a growing interest in the physiology underpinning heat tolerance of ectotherms and their responses to the ongoing rise in temperature. However, there is no consensus about the underlying physiological mechanisms. According to "the maintain aerobic scope and regulate oxygen supply" hypothesis, responses to warming at different organizational levels contribute to the ability to safeguard energy metabolism via aerobic pathways. At the cellular level, a decrease in cell size increases the capacity for the uptake of resources (e.g., food and oxygen), but the maintenance of electrochemical gradients across cellular membranes implies greater energetic costs in small cells. In this study, we investigated how different rearing temperatures affected cell size and heat tolerance in the fruit fly Drosophila melanogaster. We tested the hypothesis that smaller-celled flies are more tolerant to acute, intense heat stress whereas larger-celled flies are more tolerant to chronic, mild heat stress. We used the thermal tolerance landscape framework, which incorporates the intensity and duration of thermal challenge. Rearing temperatures strongly affected both cell size and survival times. We found different effects of developmental plasticity on tolerance to either chronic or acute heat stress. Warm-reared flies had both smaller cells and exhibited higher survival times under acute, intense heat stress when compared to cold-reared flies. However, under chronic, mild heat stress, the situation was reversed and cold-reared flies, consisting of larger cells, showed better survival. These differences in heat tolerance could have resulted from direct effects of rearing temperature or they may be mediated by the correlated changes in cell size. Notably, our results are consistent with the idea that a smaller cell size may confer tolerance to acute temperatures via enhanced oxygen supply, while a larger cell may confer greater tolerance to chronic and less intense heat stress via more efficient use of resources.

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.

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.

Identification of novel mutations in patients with Leber congenital amaurosis and juvenile RP by genome-wide homozygosity mapping with SNP microarrays.

PURPOSE: Leber congenital amaurosis (LCA) and juvenile retinitis pigmentosa (RP) cause severe visual impairment early in life. Thus far, mutations in 13 genes have been associated with autosomal recessive LCA and juvenile RP. The purpose of this study was to use homozygosity mapping to identify mutations in known LCA and juvenile RP genes. METHODS: The genomes of 93 consanguineous and nonconsanguineous patients with LCA and juvenile RP were analyzed for homozygous chromosomal regions by using SNP microarrays. This patient cohort was highly selected, as mutations in the known genes had been excluded with the LCA mutation chip, or a significant number of LCA genes had been excluded by comprehensive mutation analysis. Known LCA and juvenile RP genes residing in the identified homozygous regions were analyzed by sequencing. Detailed ophthalmic examinations were performed on the genotyped patients. RESULTS: Ten homozygous mutations, including seven novel mutations, were identified in the CRB1, LRAT, RPE65, and TULP1 genes in 12 patients. Ten patients were from consanguineous marriages, but in two patients no consanguinity was reported. In 10 of the 12 patients, the causative mutation was present in the largest or second largest homozygous segment of the patient's genome. CONCLUSIONS: Homozygosity mapping using SNP microarrays identified mutations in a significant proportion (30%) of consanguineous patients with LCA and juvenile RP and in a small number (3%) of nonconsanguineous patients. Significant homozygous regions which did not map to known LCA or juvenile RP genes and may be instrumental in identifying novel disease genes were detected in 33 patients.

Chromosomal breakpoint mapping by arrayCGH using flow-sorted chromosomes.

Despite the recent completion of the human genome project, the mapping of disease-related chromosomal translocation breakpoints and genes has remained laborious. Here, we describe a novel and rapid procedure to map such translocation breakpoints using flow-sorted chromosomes in combination with array-based comparative genomic hybridization (arrayCGH). To test the feasibility of this approach, we used a t(12;15)(q13;q25)-positive cell line with known breakpoint positions as a model. The derivative 12 chromosomes were flow-sorted, labeled, and hybridized to a genome-wide array containing 3648 well-characterized human genomic clones. The exact locations of the breakpoints on both chromosome 12 and 15 could be determined in a single hybridization experiment. In addition, we have tested the minimal amount of material necessary to perform these experiments and show that it is possible to obtain highly reliable profiles using as little as 10,000 flow-sorted chromosomes.

Reliable noninvasive prenatal testing by massively parallel sequencing of circulating cell-free DNA from maternal plasma processed up to 24h after venipuncture.

OBJECTIVES: Circulating cell-free fetal DNA (ccffDNA) in maternal plasma is an attractive source for noninvasive prenatal testing (NIPT). The amount of total cell-free DNA significantly increases 24h after venipuncture, leading to a relative decrease of the ccffDNA fraction in the blood sample. In this study, we evaluated the downstream effects of extended processing times on the reliability of aneuploidy detection by massively parallel sequencing (MPS). DESIGN AND METHODS: Whole blood from pregnant women carrying normal and trisomy 21 (T21) fetuses was collected in regular EDTA anti-coagulated tubes and processed within 6h, 24 and 48h after venipuncture. Samples of all three different time points were further analyzed by MPS using Z-score calculation and the percentage of ccffDNA based on X-chromosome reads. RESULTS: Both T21 samples were correctly identified as such at all time-points. However, after 48h, a higher deviation in Z-scores was noticed. Even though the percentage of ccffDNA in a plasma sample has been shown previously to significantly decrease 24h after venipuncture, the percentages based on MPS results did not show a significant decrease after 6, 24 or 48h. CONCLUSIONS: The quality and quantity of ccffDNA extracted from plasma samples processed up to 24h after venipuncture are sufficiently high for reliable downstream NIPT analysis by MPS. Furthermore, we show that it is important to determine the percentage of ccffDNA in the fraction of the sample that is actually used for NIPT, as downstream procedures might influence the fetal or maternal fraction.

Novel chromosomal imbalances in mantle cell lymphoma detected by genome-wide array-based comparative genomic hybridization.

Mantle cell lymphoma (MCL) is an aggressive, highly proliferative B-cell non-Hodgkin lymphoma, characterized by the specific t(11;14)(q13;q32) translocation. It is well established that this translocation alone is not sufficient to promote MCL development, but that additional genetic changes are essential for malignant transformation. We have identified such additional tumorigenic triggers in MCL tumors, by applying genome-wide array-based comparative genomic hybridization with an 800-kilobase (kb) resolution. This strategy, combined with a newly developed statistical approach, enabled us to confirm previously reported genomic alterations such as loss of 1p, 6q, 11q, 13q and gain of 3q and 8q, but it also facilitated the detection of novel recurrent genomic imbalances, such as gain of 4p12-13 and loss of 20p12.1-12.3, 20q12-13.2, 22q12.1-12.3, and 22q13.31-13.32. Genomic hotspot detection allowed for the identification of small genomic intervals that are frequently affected (57%-93%), resulting in interesting positional candidate genes such as KITLG, GPC5, and ING1. Finally, by assessing multiple biopsies from the same patient, we show that seemingly stable genomes do show subtle genomic changes over time. The follow-up of multiple biopsies of patients with MCL by high-resolution genomic profiling is expected to provide us with new clues regarding the relation between clinical outcome and in vivo cytogenetic evolution.

Diagnostic genome profiling in mental retardation.

Mental retardation (MR) occurs in 2%-3% of the general population. Conventional karyotyping has a resolution of 5-10 million bases and detects chromosomal alterations in approximately 5% of individuals with unexplained MR. The frequency of smaller submicroscopic chromosomal alterations in these patients is unknown. Novel molecular karyotyping methods, such as array-based comparative genomic hybridization (array CGH), can detect submicroscopic chromosome alterations at a resolution of 100 kb. In this study, 100 patients with unexplained MR were analyzed using array CGH for DNA copy-number changes by use of a novel tiling-resolution genomewide microarray containing 32,447 bacterial artificial clones. Alterations were validated by fluorescence in situ hybridization and/or multiplex ligation-dependent probe amplification, and parents were tested to determine de novo occurrence. Reproducible DNA copy-number changes were present in 97% of patients. The majority of these alterations were inherited from phenotypically normal parents, which reflects normal large-scale copy-number variation. In 10% of the patients, de novo alterations considered to be clinically relevant were found: seven deletions and three duplications. These alterations varied in size from 540 kb to 12 Mb and were scattered throughout the genome. Our results indicate that the diagnostic yield of this approach in the general population of patients with MR is at least twice as high as that of standard GTG-banded karyotyping.

Array-based comparative genomic hybridization for the genomewide detection of submicroscopic chromosomal abnormalities.

Microdeletions and microduplications, not visible by routine chromosome analysis, are a major cause of human malformation and mental retardation. Novel high-resolution, whole-genome technologies can improve the diagnostic detection rate of these small chromosomal abnormalities. Array-based comparative genomic hybridization allows such a high-resolution screening by hybridizing differentially labeled test and reference DNAs to arrays consisting of thousands of genomic clones. In this study, we tested the diagnostic capacity of this technology using approximately 3,500 flourescent in situ hybridization-verified clones selected to cover the genome with an average of 1 clone per megabase (Mb). The sensitivity and specificity of the technology were tested in normal-versus-normal control experiments and through the screening of patients with known microdeletion syndromes. Subsequently, a series of 20 cytogenetically normal patients with mental retardation and dysmorphisms suggestive of a chromosomal abnormality were analyzed. In this series, three microdeletions and two microduplications were identified and validated. Two of these genomic changes were identified also in one of the parents, indicating that these are large-scale genomic polymorphisms. Deletions and duplications as small as 1 Mb could be reliably detected by our approach. The percentage of false-positive results was reduced to a minimum by use of a dye-swap-replicate analysis, all but eliminating the need for laborious validation experiments and facilitating implementation in a routine diagnostic setting. This high-resolution assay will facilitate the identification of novel genes involved in human mental retardation and/or malformation syndromes and will provide insight into the flexibility and plasticity of the human genome.