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.

Invariance of SCL-90-R dimensions of symptom distress in patients with peri partum pelvic pain (PPPP) syndrome.

OBJECTIVES: There are no studies available that have examined the factorial invariance of dimensions underlying the Symptom Checklist-90-Revised (SCL-90-R) across at least three distinct samples. In the following study, we wished to determine whether a dimensional model comprising eight primary factors previously identified in psychiatric out-patients, phobics and the general population (Arrindell & Ettema, 2003) could be extended to a homogeneous sample of pain patients comprising females suffering from peri partum pelvic pain (PPPP) syndrome (N = 413). The internal consistency and discriminant validity of the dimensions were also examined. METHOD: The SCL-90-R and measures of disability, pain-related fear, pain intensity and fatigue were administered to the participants. The multiple group method was used to determine factorial invariance. Pearson correlations were determined between the SCL-90-R and aforementioned measures. RESULTS: The factorial invariance of an 8-dimensional model of primary factors underlying the SCL-90-R, namely, agoraphobia, anxiety, depression, somatization, cognitive-performance deficits, interpersonal sensitivity-mistrust, acting-out hostility and sleep difficulties, was extended with success to the present sample of PPPP patients. In spite of substantial correlations between the internally consistent SCL-90-R symptom dimensions, some evidence of discriminant validity was reported in that specific subscales showed different patterns of correlations with measures of disability, pain-related fear, pain intensity and fatigue. CONCLUSIONS: The 8-dimensional system based on the work of Arrindell and Ettema (2003) was invariant across psychiatric patients, phobics, the general population and pain patients. The invariance of the SCL-90-R hostility dimensions may have implications for a re-formulation of Watson and Clark's tripartite model of general distress, specific anxiety and specific depression.

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.

Molecular cytogenetic analysis of clustered sporadic and familial renal cell carcinoma-associated 3q13 approximately q22 breakpoints.

We describe several relatives within one renal cell cancer (RCC) family sharing a constitutional t(2;3) (q35;q21). Based on molecular studies on several independent primary tumors in this family, a causative role for this translocation in tumor development was suggested. Subsequent positional cloning of the 3q21 chromosomal breakpoint revealed that this breakpoint disrupts a novel gene, DIRC2 (disrupted in renal cancer 2). This gene encodes an evolutionary conserved transmembrane protein and represents a novel member of the MFS superfamily of transporters. To evaluate whether DIRC2 is also targeted in sporadic RCC cases with cytogenetically defined 3q21 breakpoints, fluorescence in situ hybridization analysis was performed on metaphase spreads and/or interphase nuclei of 12 primary sporadic RCC using genomic clones from a 3q21 breakpoint-spanning contig as probes. Three breakpoints were mapped proximal to the familial breakpoint and nine breakpoints were mapped distal to this breakpoint. Two of the latter breakpoints were mapped in the contig within 1 Mb distance from the familial breakpoint. Because these clustered 3q21 breakpoints do not coincide with the familial 3q21 breakpoint, they most likely affect genes distinct from DIRC2.

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.

Non-invasive prenatal diagnosis of fetal aneuploidies using massively parallel sequencing-by-ligation and evidence that cell-free fetal DNA in the maternal plasma originates from cytotrophoblastic cells.

Blood plasma of pregnant women contains circulating cell-free fetal DNA (ccffDNA), originating from the placenta. The use of this DNA for non-invasive detection of fetal aneuploidies using massively parallel sequencing (MPS)-by-synthesis has been proven previously. Sequence performance may, however, depend on the MPS platform and therefore we have explored the possibility for multiplex MPS-by-ligation, using the Applied Biosystems SOLiD(™) 4 system. DNA isolated from plasma samples from 52 pregnant women, carrying normal or aneuploid fetuses, was sequenced in multiplex runs of 4, 8 or 16 samples simultaneously. The sequence reads were mapped to the human reference genome and quantified according to their genomic location. In case of a fetal aneuploidy, the number of reads of the aberrant chromosome is expected to be higher or lower than in normal reference samples. To statistically determine this, Z-scores per chromosome were calculated as described previously, with thresholds for aneuploidies set at > +3.0 and < -3.0 for chromosomal over- or underrepresentation, respectively. All samples from fetal aneuploidies yielded Z-scores outside the thresholds for the aberrant chromosomes, with no false negative or positive results. Full-blown fetal aneuploidies can thus be reliably detected in maternal plasma using a multiplex MPS-by-ligation approach. Furthermore, the results obtained with a sample from a pregnancy with 45,X in the cytotrophoblastic cell layer and 46,XX in the mesenchymal core cells show that ccffDNA originates from the cytotrophoblastic cell layer. Discrepancies between the genetic constitution of this cell layer and the fetus itself are well known, and therefore, care should be taken when translating results to the fetus itself.

De novo nonsense mutations in ASXL1 cause Bohring-Opitz syndrome.

Bohring-Opitz syndrome is characterized by severe intellectual disability, distinctive facial features and multiple congenital malformations. We sequenced the exomes of three individuals with Bohring-Opitz syndrome and in each identified heterozygous de novo nonsense mutations in ASXL1, which is required for maintenance of both activation and silencing of Hox genes. In total, 7 out of 13 subjects with a Bohring-Opitz phenotype had de novo ASXL1 mutations, suggesting that the syndrome is genetically heterogeneous.

A de novo paradigm for mental retardation.

The per-generation mutation rate in humans is high. De novo mutations may compensate for allele loss due to severely reduced fecundity in common neurodevelopmental and psychiatric diseases, explaining a major paradox in evolutionary genetic theory. Here we used a family based exome sequencing approach to test this de novo mutation hypothesis in ten individuals with unexplained mental retardation. We identified and validated unique non-synonymous de novo mutations in nine genes. Six of these, identified in six different individuals, are likely to be pathogenic based on gene function, evolutionary conservation and mutation impact. Our findings provide strong experimental support for a de novo paradigm for mental retardation. Together with de novo copy number variation, de novo point mutations of large effect could explain the majority of all mental retardation cases in the population.

Identification of tumor-specific molecular signatures in intracranial ependymoma and association with clinical characteristics.

PURPOSE: To delineate clinically relevant molecular signatures of intracranial ependymoma. MATERIALS AND METHODS: We analyzed 24 primary intracranial ependymomas. For genomic profiling, microarray-based comparative genomic hybridization (CGH) was used and results were validated by fluorescent in situ hybridization and loss of heterozygosity mapping. We performed gene expression profiling using microarrays, real-time quantitative reverse transcriptase polymerase chain reaction, and methylation analysis of selected genes. We applied class comparison analyses to compare both genomic and expression profiling data with clinical characteristics. RESULTS: A variable number of genomic imbalances were detected by array CGH, revealing multiple regions of recurrent gain (including 2q23, 7p21, 12p, 13q21.1, and 20p12) and loss (including 5q31, 6q26, 7q36, 15q21.1, 16q24, 17p13.3, 19p13.2, and 22q13.3). An ependymoma-specific gene expression signature was characterized by the concurrent abnormal expression of developmental and differentiation pathways, including NOTCH and sonic hedgehog signaling. We identified specific differentially imbalanced genomic clones and gene expression signatures significantly associated with tumor location, patient age at disease onset, and retrospective risk for relapse. Integrated genomic and expression profiling allowed us to identify genes of which the expression is deregulated in intracranial ependymoma, such as overexpression of the putative proto-oncogene YAP1 (located at 11q22) and downregulation of the SULT4A1 gene (at 22q13.3). CONCLUSION: The present exploratory molecular profiling study allowed us to refine previously reported intervals of genomic imbalance, to identify novel restricted regions of gain and loss, and to identify molecular signatures correlating with various clinical variables. Validation of these results on independent data sets represents the next step before translation into the clinical setting.

Identification of recurrent chromosomal aberrations in germ cell tumors of neonates and infants using genomewide array-based comparative genomic hybridization.

Human germ cell tumors (GCTs) of neonates and infants comprise a heterogeneous group of neoplasms, including teratomas and yolk sac tumors with distinct clinical and epidemiologic features. As yet, little is known about the cytogenetic constitution of these tumors. We applied the recently developed genomewide array-based comparative genomic hybridization (array CGH) technology to 24 GCTs derived from patients under the age of 5 years. In addition, we included seven tumors derived from children and adolescents older than 5 years. In the series from those under the age of 5 years, most teratomas displayed normal profiles, except for some minor recurrent aberrations. In contrast, the yolk sac tumors displayed recurrent losses of 1p35-pter and gains of 3p21-pter and of 20q13. In the GCTs of patients older than 5 years, the main recurrent anomalies included gains of 12p and of whole chromosomes 7 and 8. In addition, gains of the 1q32-qter region and losses of the 6q24-qter and 18q21-qter regions were frequent in GCTs of varied histology, independent of age. We concluded that array CGH is a highly suitable method for identifying recurrent chromosomal anomalies in GCTs of neonates and infants. The recurrent anomalies observed point to chromosomal regions that may harbor novel diagnostic/prognostic identifiers and genes relevant to the development of these neoplasms.