Mutation screening in 86 known X-linked mental retardation genes by droplet-based multiplex PCR and massive parallel sequencing

Locus heterogeneity in two siblings presenting with developmental delay, intellectual disability and autism spectrum disorder

Correct diagnosis of children presenting with developmental delay and intellectual disability remains challenging due to the complex and heterogeneous etiology. High throughput sequencing technologies like exome sequencing have become more commonly available and are significantly improving genetic testing. We present two siblings - a 14-year old male and an 8-year old female patient - with a similar clinical phenotype that was characterized by combined developmental delay primarily affecting speech, mild to moderate intellectual disability, behavioral abnormalities, and autism spectrum disorder, but with no congenital anomalies. The sister showed additional muscular hypotonia and more pronounced dysmorphic features compared to her brother. Both parents had psychiatric disorders and mild to moderate intellectual disability. A common genetic etiology in the siblings was suspected. Metabolic, psychological and neuroradiological examinations were complemented by basic genetic testing including chromosome analysis and array comparative genomics hybridization analysis (CGH), followed by exome sequencing and combined data analysis of the family. Exome sequencing identified two different underlying genetic conditions: in the sister, a maternally inherited pathogenic variant c.1661C > T, p.Pro554Leu in SLC6A8 (NM_005629.4) was identified causing cerebral creatine deficiency syndrome 1 (MIM #300352) which was confirmed by MR spectroscopy and treated accordingly. In the brother, a paternally inherited 16p13.11 duplication was identified by exome sequencing and considered to be likely associated with his and possibly his father's phenotype. The 16p13.11 duplication had been previously identified in an array CGH but had not been prioritized due to the lack of segregation in the siblings. In conclusion, we report a case of intra-familial locus heterogeneity of developmental delay in two siblings. We advocate for the need of unbiased and comprehensive genetic testing to provide accurate diagnosis despite locus heterogeneity.

Assorted dysfunctions of endosomal alkali cation/proton exchanger SLC9A6 variants linked to Christianson syndrome

Genetic screening has identified numerous variants of the endosomal solute carrier family 9 member A6 (SLC9A6)/(Na+,K+)/H+ exchanger 6 (NHE6) gene that cause Christianson syndrome, a debilitating X-linked developmental disorder associated with a range of neurological, somatic, and behavioral symptoms. Many of these variants cause complete loss of NHE6 expression, but how subtler missense substitutions or nonsense mutations that partially truncate its C-terminal cytoplasmic regulatory domain impair NHE6 activity and endosomal function are poorly understood. Here, we describe the molecular and cellular consequences of six unique mutations located in the N-terminal cytoplasmic segment (A9S), the membrane ion translocation domain (L188P and G383D), and the C-terminal regulatory domain (E547*, R568Q, and W570*) of human NHE6 that purportedly cause disease. Using a heterologous NHE6-deficient cell expression system, we show that the biochemical, catalytic, and cellular properties of the A9S and R568Q variants were largely indistinguishable from those of the WT transporter, which obscured their disease significance. By contrast, the L188P, G383D, E547*, and W570* mutants exhibited variable deficiencies in biosynthetic post-translational maturation, membrane sorting, pH homeostasis in recycling endosomes, and cargo trafficking, and they also triggered apoptosis. These findings broaden our understanding of the molecular dysfunctions of distinct NHE6 variants associated with Christianson syndrome.

Thyroid Hormone Transporters

Thyroid hormone transporters at the plasma membrane govern intracellular bioavailability of thyroid hormone. Monocarboxylate transporter (MCT) 8 and MCT10, organic anion transporting polypeptide (OATP) 1C1, and SLC17A4 are currently known as transporters displaying the highest specificity toward thyroid hormones. Structure-function studies using homology modeling and mutational screens have led to better understanding of the molecular basis of thyroid hormone transport. Mutations in MCT8 and in OATP1C1 have been associated with clinical disorders. Different animal models have provided insight into the functional role of thyroid hormone transporters, in particular MCT8. Different treatment strategies for MCT8 deficiency have been explored, of which thyroid hormone analogue therapy is currently applied in patients. Future studies may reveal the identity of as-yet-undiscovered thyroid hormone transporters. Complementary studies employing animal and human models will provide further insight into the role of transporters in health and disease. (Endocrine Reviews 41: 1 - 55, 2020).

The complete loss of function of the SMS gene results in a severe form of Snyder-Robinson syndrome

Snyder-Robinson syndrome (SRS) is an X-linked syndromic intellectual disability condition caused by variants in the spermine synthase gene (SMS). The syndrome is characterized by facial dysmorphism, thin body build, kyphoscoliosis, osteoporosis, hypotonia, developmental delay and associated neurological features (seizures, unsteady gait, abnormal speech). Until now, only missense variants with a functionally characterized partial loss of function (LoF) have been described. Here we describe the first complete LoF variant, Met303Lysfs*, in a male patient with a severe form of Snyder-Robinson syndrome. He presented with multiple malformations and severly delayed development, and died at 4 months of age. Functional in vitro assays showed a complete absence of functional SMS protein. Taken together, our findings and those of previously reported patients confirm that pathogenic variants of SMS are indeed LoF and that there might exist a genotype-phenotype correlation between the type of variant and the severity of the syndrome.

Amyloid clearance defect in ApoE4 astrocytes is reversed by epigenetic correction of endosomal pH

Endosomes have emerged as a central hub and pathogenic driver of Alzheimer's disease (AD). The earliest brain cytopathology in neurodegeneration, occurring decades before amyloid plaques and cognitive decline, is an expansion in the size and number of endosomal compartments. The strongest genetic risk factor for sporadic AD is the ε4 allele of Apolipoprotein E (ApoE4). Previous studies have shown that ApoE4 potentiates presymptomatic endosomal dysfunction and defective endocytic clearance of amyloid beta (Aβ), although how these two pathways are linked at a cellular and mechanistic level has been unclear. Here, we show that aberrant endosomal acidification in ApoE4 astrocytes traps the low-density lipoprotein receptor-related protein (LRP1) within intracellular compartments, leading to loss of surface expression and Aβ clearance. Pathological endosome acidification is caused by ε4 risk allele-selective down-regulation of the Na⁺/H⁺ exchanger isoform NHE6, which functions as a critical leak pathway for endosomal protons. In vivo, the NHE6 knockout (NHE6^KO) mouse model showed elevated Aβ in the brain, consistent with a causal effect. Increased nuclear translocation of histone deacetylase 4 (HDAC4) in ApoE4 astrocytes, compared with the nonpathogenic ApoE3 allele, suggested a mechanistic basis for transcriptional down-regulation of NHE6. HDAC inhibitors that restored NHE6 expression normalized ApoE4-specific defects in endosomal pH, LRP1 trafficking, and amyloid clearance. Thus, NHE6 is a downstream effector of ApoE4 and emerges as a promising therapeutic target in AD. These observations have prognostic implications for patients who have Christianson syndrome with loss of function mutations in NHE6 and exhibit prominent glial pathology and progressive hallmarks of neurodegeneration.

De novo and inherited mutations in the X-linked gene CLCN4 are associated with syndromic intellectual disability and behavior and seizure disorders in males and females

Variants in CLCN4, which encodes the chloride/hydrogen ion exchanger CIC-4 prominently expressed in brain, were recently described to cause X-linked intellectual disability and epilepsy. We present detailed phenotypic information on 52 individuals from 16 families with CLCN4-related disorder: 5 affected females and 2 affected males with a de novo variant in CLCN4 (6 individuals previously unreported) and 27 affected males, 3 affected females and 15 asymptomatic female carriers from 9 families with inherited CLCN4 variants (4 families previously unreported). Intellectual disability ranged from borderline to profound. Behavioral and psychiatric disorders were common in both child- and adulthood, and included autistic features, mood disorders, obsessive-compulsive behaviors and hetero- and autoaggression. Epilepsy was common, with severity ranging from epileptic encephalopathy to well-controlled seizures. Several affected individuals showed white matter changes on cerebral neuroimaging and progressive neurological symptoms, including movement disorders and spasticity. Heterozygous females can be as severely affected as males. The variability of symptoms in females is not correlated with the X inactivation pattern studied in their blood. The mutation spectrum includes frameshift, missense and splice site variants and one single-exon deletion. All missense variants were predicted to affect CLCN4's function based on in silico tools and either segregated with the phenotype in the family or were de novo. Pathogenicity of all previously unreported missense variants was further supported by electrophysiological studies in Xenopus laevis oocytes. We compare CLCN4-related disorder with conditions related to dysfunction of other members of the CLC family.

Disorder of thyroid hormone transport into the tissues

Transport of thyroid hormone (TH) across the plasma membrane is essential for intracellular TH metabolism and action, and this is mediated by specific transporter proteins. During the last two decades several transporters capable of transporting TH have been identified, including monocarboxylate transporter 8 (MCT8), MCT10 and organic anion transporting polypeptide 1C1 (OATP1C1). In particular MCT8 and OATP1C1 are important for the regulation of local TH activity in the brain and thus for brain development. MCT8 is a protein containing 12 transmembrane domains, and is encoded by the SLC16A2 gene located on the X chromosome. It facilitates both TH uptake and efflux across the cell membrane. Male subjects with hemizygous mutations in MCT8 are afflicted with severe intellectual and motor disability, also known as the Allan-Herndon-Dudley syndrome (AHDS), which goes together with low serum T4 and high T3 levels. This review concerns molecular and clinical aspects of MCT8 function.

Mutations in two large pedigrees highlight the role of ZNF711 in X-linked intellectual disability

Intellectual disability (ID) affects approximately 1-2% of the general population and is characterized by impaired cognitive abilities. ID is both clinically as well as genetically heterogeneous, up to 2000 genes are estimated to be involved in the emergence of the disease with various clinical presentations. For many genes, only a few patients have been reported and causality of some genes has been questioned upon the discovery of apparent loss-of-function mutations in healthy controls. Description of additional patients strengthens the evidence for the involvement of a gene in the disease and can clarify the clinical phenotype associated with mutations in a particular gene. Here, we present two large four-generation families with a total of 11 males affected with ID caused by mutations in ZNF711, thereby expanding the total number of families with ID and a ZNF711 mutation to four. Patients with mutations in ZNF711 all present with mild to moderate ID and poor speech accompanied by additional features in some patients, including autistic features and mild facial dysmorphisms, suggesting that ZNF711 mutations cause non-syndromic ID.

A Christianson syndrome-linked deletion mutation (∆(287)ES(288)) in SLC9A6 disrupts recycling endosomal function and elicits neurodegeneration and cell death

Background

Christianson Syndrome, a recently identified X-linked neurodevelopmental disorder, is caused by mutations in the human gene SLC9A6 encoding the recycling endosomal alkali cation/proton exchanger NHE6. The patients have pronounced limitations in cognitive ability, motor skills and adaptive behaviour. However, the mechanistic basis for this disorder is poorly understood as few of the more than 20 mutations identified thus far have been studied in detail.

Methods

Here, we examined the molecular and cellular consequences of a 6 base-pair deletion of amino acids Glu²⁸⁷ and Ser²⁸⁸ (∆ES) in the predicted seventh transmembrane helix of human NHE6 expressed in established cell lines (CHO/AP-1, HeLa and neuroblastoma SH-SY5Y) and primary cultures of mouse hippocampal neurons by measuring levels of protein expression, stability, membrane trafficking, endosomal function and cell viability.

Results

In the cell lines, immunoblot analyses showed that the nascent mutant protein was properly synthesized and assembled as a homodimer, but its oligosaccharide maturation and half-life were markedly reduced compared to wild-type (WT) and correlated with enhanced ubiquitination leading to both proteasomal and lysosomal degradation. Despite this instability, a measurable fraction of the transporter was correctly sorted to the plasma membrane. However, the rates of clathrin-mediated endocytosis of the ∆ES mutant as well as uptake of companion vesicular cargo, such as the ligand-bound transferrin receptor, were significantly reduced and correlated with excessive endosomal acidification. Notably, ectopic expression of ∆ES but not WT induced apoptosis when examined in AP-1 cells. Similarly, in transfected primary cultures of mouse hippocampal neurons, membrane trafficking of the ∆ES mutant was impaired and elicited marked reductions in total dendritic length, area and arborization, and triggered apoptotic cell death.

Conclusions

These results suggest that loss-of-function mutations in NHE6 disrupt recycling endosomal function and trafficking of cargo which ultimately leads to neuronal degeneration and cell death in Christianson Syndrome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13024-016-0129-9) contains supplementary material, which is available to authorized users.

Empirical estimation of sequencing error rates using smoothing splines

Background

Next-generation sequencing has been used by investigators to address a diverse range of biological problems through, for example, polymorphism and mutation discovery and microRNA profiling. However, compared to conventional sequencing, the error rates for next-generation sequencing are often higher, which impacts the downstream genomic analysis. Recently, Wang et al. (BMC Bioinformatics 13:185, 2012) proposed a shadow regression approach to estimate the error rates for next-generation sequencing data based on the assumption of a linear relationship between the number of reads sequenced and the number of reads containing errors (denoted as shadows). However, this linear read-shadow relationship may not be appropriate for all types of sequence data. Therefore, it is necessary to estimate the error rates in a more reliable way without assuming linearity. We proposed an empirical error rate estimation approach that employs cubic and robust smoothing splines to model the relationship between the number of reads sequenced and the number of shadows.

Results

We performed simulation studies using a frequency-based approach to generate the read and shadow counts directly, which can mimic the real sequence counts data structure. Using simulation, we investigated the performance of the proposed approach and compared it to that of shadow linear regression. The proposed approach provided more accurate error rate estimations than the shadow linear regression approach for all the scenarios tested. We also applied the proposed approach to assess the error rates for the sequence data from the MicroArray Quality Control project, a mutation screening study, the Encyclopedia of DNA Elements project, and bacteriophage PhiX DNA samples.

Conclusions

The proposed empirical error rate estimation approach does not assume a linear relationship between the error-free read and shadow counts and provides more accurate estimations of error rates for next-generation, short-read sequencing data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1052-3) contains supplementary material, which is available to authorized users.

X-exome sequencing of 405 unresolved families identifies seven novel intellectual disability genes

X-linked intellectual disability (XLID) is a clinically and genetically heterogeneous disorder. During the past two decades in excess of 100 X-chromosome ID genes have been identified. Yet, a large number of families mapping to the X-chromosome remained unresolved suggesting that more XLID genes or loci are yet to be identified. Here, we have investigated 405 unresolved families with XLID. We employed massively parallel sequencing of all X-chromosome exons in the index males. The majority of these males were previously tested negative for copy number variations and for mutations in a subset of known XLID genes by Sanger sequencing. In total, 745 X-chromosomal genes were screened. After stringent filtering, a total of 1297 non-recurrent exonic variants remained for prioritization. Co-segregation analysis of potential clinically relevant changes revealed that 80 families (20%) carried pathogenic variants in established XLID genes. In 19 families, we detected likely causative protein truncating and missense variants in 7 novel and validated XLID genes (CLCN4, CNKSR2, FRMPD4, KLHL15, LAS1L, RLIM and USP27X) and potentially deleterious variants in 2 novel candidate XLID genes (CDK16 and TAF1). We show that the CLCN4 and CNKSR2 variants impair protein functions as indicated by electrophysiological studies and altered differentiation of cultured primary neurons from Clcn4(-/-) mice or after mRNA knock-down. The newly identified and candidate XLID proteins belong to pathways and networks with established roles in cognitive function and intellectual disability in particular. We suggest that systematic sequencing of all X-chromosomal genes in a cohort of patients with genetic evidence for X-chromosome locus involvement may resolve up to 58% of Fragile X-negative cases.

Next-generation sequencing using a pre-designed gene panel for the molecular diagnosis of congenital disorders in pediatric patients

BACKGROUND

Next-generation sequencing (NGS) has revolutionized genetic research and offers enormous potential for clinical application. Sequencing the exome has the advantage of casting the net wide for all known coding regions while targeted gene panel sequencing provides enhanced sequencing depths and can be designed to avoid incidental findings in adult-onset conditions. A HaloPlex panel consisting of 180 genes within commonly altered chromosomal regions is available for use on both the Ion Personal Genome Machine (PGM) and MiSeq platforms to screen for causative mutations in these genes.

METHODS

We used this Haloplex ICCG panel for targeted sequencing of 15 patients with clinical presentations indicative of an abnormality in one of the 180 genes. Sequencing runs were done using the Ion 318 Chips on the Ion Torrent PGM. Variants were filtered for known polymorphisms and analysis was done to identify possible disease-causing variants before validation by Sanger sequencing. When possible, segregation of variants with phenotype in family members was performed to ascertain the pathogenicity of the variant.

RESULTS

More than 97% of the target bases were covered at >20×. There was an average of 9.6 novel variants per patient. Pathogenic mutations were identified in five genes for six patients, with two novel variants. There were another five likely pathogenic variants, some of which were unreported novel variants.

CONCLUSIONS

In a cohort of 15 patients, we were able to identify a likely genetic etiology in six patients (40%). Another five patients had candidate variants for which further evaluation and segregation analysis are ongoing. Our results indicate that the HaloPlex ICCG panel is useful as a rapid, high-throughput and cost-effective screening tool for 170 of the 180 genes. There is low coverage for some regions in several genes which might have to be supplemented by Sanger sequencing. However, comparing the cost, ease of analysis, and shorter turnaround time, it is a good alternative to exome sequencing for patients whose features are suggestive of a genetic etiology involving one of the genes in the panel.

Diagnostic Approach to Hereditary Colorectal Cancer Syndromes

Approximately 5 to 10% of colorectal cancers develop within a known hereditary syndrome. Specific underlying genetic mutations drive the clinical phenotype and it is imperative to determine the genetic etiology to provide meaningful surveillance and intervention. Recognizing potential patients and families with a hereditary predisposition is the first step in management. Syndromes can be categorized according to polyp burden as polyposis or nonpolyposis. Clinical assessment should start with a personal and family medical history, physical examination, and evaluation for the presence and type of colorectal polyps or cancers. Key information is gained from these simple steps and should guide the specific genetic analysis for diagnosis. Genetic counseling is a critical component to any hereditary colorectal cancer program and should be conducted before genetic testing to provide education about the implications of test results. This review focuses on the thought process that drives initial clinical evaluation and guides genetic testing for patients with suspected hereditary colorectal cancer syndromes.

Integrated sequence analysis pipeline provides one-stop solution for identifying disease-causing mutations

Next-generation sequencing has greatly accelerated the search for disease-causing defects, but even for experts the data analysis can be a major challenge. To facilitate the data processing in a clinical setting, we have developed a novel medical resequencing analysis pipeline (MERAP). MERAP assesses the quality of sequencing, and has optimized capacity for calling variants, including single-nucleotide variants, insertions and deletions, copy-number variation, and other structural variants. MERAP identifies polymorphic and known causal variants by filtering against public domain databases, and flags nonsynonymous and splice-site changes. MERAP uses a logistic model to estimate the causal likelihood of a given missense variant. MERAP considers the relevant information such as phenotype and interaction with known disease-causing genes. MERAP compares favorably with GATK, one of the widely used tools, because of its higher sensitivity for detecting indels, its easy installation, and its economical use of computational resources. Upon testing more than 1,200 individuals with mutations in known and novel disease genes, MERAP proved highly reliable, as illustrated here for five families with disease-causing variants. We believe that the clinical implementation of MERAP will expedite the diagnostic process of many disease-causing defects.

Mung bean nuclease treatment increases capture specificity of microdroplet-PCR based targeted DNA enrichment

Targeted DNA enrichment coupled with next generation sequencing has been increasingly used for interrogation of select sub-genomic regions at high depth of coverage in a cost effective manner. Specificity measured by on-target efficiency is a key performance metric for target enrichment. Non-specific capture leads to off-target reads, resulting in waste of sequencing throughput on irrelevant regions. Microdroplet-PCR allows simultaneous amplification of up to thousands of regions in the genome and is among the most commonly used strategies for target enrichment. Here we show that carryover of single-stranded template genomic DNA from microdroplet-PCR constitutes a major contributing factor for off-target reads in the resultant libraries. Moreover, treatment of microdroplet-PCR enrichment products with a nuclease specific to single-stranded DNA alleviates off-target load and improves enrichment specificity. We propose that nuclease treatment of enrichment products should be incorporated in the workflow of targeted sequencing using microdroplet-PCR for target capture. These findings may have a broad impact on other PCR based applications for which removal of template DNA is beneficial.

An inside job: how endosomal Na(+)/H(+) exchangers link to autism and neurological disease

Autism imposes a major impediment to childhood development and a huge emotional and financial burden on society. In recent years, there has been rapidly accumulating genetic evidence that links the eNHE, a subset of Na(+)/H(+) exchangers that localize to intracellular vesicles, to a variety of neurological conditions including autism, attention deficit hyperactivity disorder (ADHD), intellectual disability, and epilepsy. By providing a leak pathway for protons pumped by the V-ATPase, eNHE determine luminal pH and regulate cation (Na(+), K(+)) content in early and recycling endosomal compartments. Loss-of-function mutations in eNHE cause hyperacidification of endosomal lumen, as a result of imbalance in pump and leak pathways. Two isoforms, NHE6 and NHE9 are highly expressed in brain, including hippocampus and cortex. Here, we summarize evidence for the importance of luminal cation content and pH on processing, delivery and fate of cargo. Drawing upon insights from model organisms and mammalian cells we show how eNHE affect surface expression and function of membrane receptors and neurotransmitter transporters. These studies lead to cellular models of eNHE activity in pre- and post-synaptic neurons and astrocytes, where they could impact synapse development and plasticity. The study of eNHE has provided new insight on the mechanism of autism and other debilitating neurological disorders and opened up new possibilities for therapeutic intervention.

X-exome sequencing in Finnish families with intellectual disability--four novel mutations and two novel syndromic phenotypes

Background

X-linked intellectual disability (XLID) is a group of genetically heterogeneous disorders characterized by substantial impairment in cognitive abilities, social and behavioral adaptive skills. Next generation sequencing technologies have become a powerful approach for identifying molecular gene mutations relevant for diagnosis.

Methods & objectives

Enrichment of X-chromosome specific exons and massively parallel sequencing was performed for identifying the causative mutations in 14 Finnish families, each of them having several males affected with intellectual disability of unknown cause.

Results

We found four novel mutations in known XLID genes. Two mutations; one previously reported missense mutation (c.1111C > T), and one novel frameshift mutation (c. 990_991insGCTGC) were identified in SLC16A2, a gene that has been linked to Allan-Herndon-Dudley syndrome (AHDS). One novel missense mutation (c.1888G > C) was found in GRIA3 and two novel splice donor site mutations (c.357 + 1G > C and c.985 + 1G > C) were identified in the DLG3 gene. One missense mutation (c.1321C > T) was identified in the candidate gene ZMYM3 in three affected males with a previously unrecognized syndrome characterized by unique facial features, aortic stenosis and hypospadia was detected. All of the identified mutations segregated in the corresponding families and were absent in > 100 Finnish controls and in the publicly available databases. In addition, a previously reported benign variant (c.877G > A) in SYP was identified in a large family with nine affected males in three generations, who have a syndromic phenotype.

Conclusions

All of the mutations found in this study are being reported for the first time in Finnish families with several affected male patients whose etiological diagnoses have remained unknown to us, in some families, for more than 30 years. This study illustrates the impact of X-exome sequencing to identify rare gene mutations and the challenges of interpreting the results. Further functional studies are required to confirm the cause of the syndromic phenotypes associated with ZMYM3 and SYP in this study.

Target enrichment using parallel nanoliter quantitative PCR amplification

BACKGROUND

Next generation targeted resequencing is replacing Sanger sequencing at high pace in routine genetic diagnosis. The need for well validated, high quality enrichment platforms to complement the bench-top next generation sequencing devices is high.

RESULTS

We used the WaferGen Smartchip platform to perform highly parallelized PCR based target enrichment for a set of known cancer genes in a well characterized set of cancer cell lines from the NCI60 panel. Optimization of PCR assay design and cycling conditions resulted in a high enrichment efficiency. We provide proof of a high mutation rediscovery rate and have included technical replicates to enable SNP calling validation demonstrating the high reproducibility of our enrichment platform.

CONCLUSIONS

Here we present our custom developed quantitative PCR based target enrichment platform. Using highly parallel nanoliter singleplex PCR reactions makes this a flexible and efficient platform. The high mutation validation rate shows this platform's promise as a targeted resequencing method for multi-gene routine sequencing diagnostics.

Next generation sequencing in nonsyndromic intellectual disability: from a negative molecular karyotype to a possible causative mutation detection

The identification of causes underlying intellectual disability (ID) is one of the most demanding challenges for clinical Geneticists and Researchers. Despite molecular diagnostics improvements, the vast majority of patients still remain without genetic diagnosis. Here, we report the results obtained using Whole Exome and Target Sequencing on nine patients affected by isolated ID without pathological copy number variations, which were accurately selected from an initial cohort of 236 patients. Three patterns of inheritance were used to search for: (1) de novo, (2) X-linked, and (3) autosomal recessive variants. In three of the nine proband-parent trios analyzed, we identified and validated two de novo and one X-linked potentially causative mutations located in three ID-related genes. We proposed three genes as ID candidate, carrying one de novo and three X-linked mutations. Overall, this systematic proband-parent trio approach using next generation sequencing could explain a consistent percentage of patients with isolated ID, thus increasing our knowledge on the molecular bases of this disease and opening new perspectives for a better diagnosis, counseling, and treatment.

Detection of a novel intragenic rearrangement in the creatine transporter gene by next generation sequencing

Deficiency caused by mutations in the creatine transporter gene (SLC6A8/CT1) is an X-linked form of intellectual disability. The presence of highly homologous pseudogenes and high GC content of SLC6A8 genomic sequence complicates the molecular diagnosis of this disorder. To minimize the pseudogene interference, exons 2 to 13 of SLC6A8 were amplified as a single PCR product using gene-specific long-range PCR (LR-PCR) primers. The GC-rich exon 1 and its flanking intronic sequences were amplified separately in a short fragment under GC-rich conditions and a touchdown PCR program. Traditional Sanger sequence analysis of all coding exons of SLC6A8 from a 3-year-old boy with creatine transporter deficiency did not detect deleterious mutations. The long-range PCR product was used as template followed by massively parallel sequencing (MPS) on HiSeq2000. We were able to detect a tandem duplication involving part of exons 11 and 12 in the SLC6A8 gene. The deduced c.1592_1639dup133 mutation was confirmed to be a hemizygous insertion by targeted genomic DNA and cDNA Sanger sequencing. Combination of deep sequencing technology with long-range PCR revealed a novel intragenic duplication in the SLC6A8 gene, providing a definitive molecular diagnosis of creatine transporter deficiency in a male patient.

Massively parallel sequencing for diagnosing clinically and genetically heterogeneous disorders

Massively parallel sequencing technology has great potential to identify the precise molecular causes of human disorders, which were once considered difficult to diagnose for either inherited Mendelian or oncological diseases. Besides being a valuable tool for gene discovery, this technology has also gained tremendous momentum in clinical molecular diagnostic laboratories. There are a wide variety of clinical applications recently developed to meet the clinical demands for personalized medicine. This article discusses these up-to-date massively parallel sequencing clinical applications and emphasizes various disorders that can be targeted, as well as some challenges faced in the process of implementing these assays in the clinical setting.

A 3-base pair deletion, c.9711_9713del, in DMD results in intellectual disability without muscular dystrophy

We have identified a deletion of 3 base pairs in the dystrophin gene (DMD), c.9711_9713del, in a family with nonspecific X-linked intellectual disability (ID) by sequencing of the exons of 86 known X-linked ID genes. This in-frame deletion results in the deletion of a single-amino-acid residue, Leu3238, in the brain-specific isoform Dp71 of dystrophin. Linkage analysis supported causality as the mutation was present in the 7.6 cM linkage interval on Xp22.11-Xp21.1 with a maximum positive LOD score of 2.41 (MRX85 locus). Molecular modeling predicts that the p.(Leu3238del) deletion results in the destabilization of the C-terminal domain of dystrophin and hence reduces the ability to interact with β-dystroglycan. Correspondingly, Dp71 protein levels in lymphoblastoid cells from the index patient are 6.7-fold lower than those in control cell lines (P=0.08). Subsequent determination of the creatine kinase levels in blood of the index patient showed a mild but significant elevation in serum creatine kinase, which is in line with impaired dystrophin function. In conclusion, we have identified the first DMD mutation in Dp71 that results in ID without muscular dystrophy.

A systematic evaluation of hybridization-based mouse exome capture system

BACKGROUND

Exome sequencing is increasingly used to search for phenotypically-relevant sequence variants in the mouse genome. All of the current hybridization-based mouse exome capture systems are designed based on the genome reference sequences of the C57BL/6 J strain. Given that the substantial sequence divergence exists between C57BL/6 J and other distantly-related strains, the impact of sequence divergence on the efficiency of such capture systems needs to be systematically evaluated before they can be widely applied to the study of those strains.

RESULTS

Using the Agilent SureSelect mouse exome capture system, we performed exome sequencing on F1 generation hybrid mice that were derived by crossing two divergent strains, C57BL/6 J and SPRET/EiJ. Our results showed that the C57BL/6 J-based probes captured the sequences derived from C57BL/6 J alleles more efficiently and that the bias was higher for the target regions with greater sequence divergence. At low sequencing depths, the bias also affected the efficiency of variant detection. However, the effects became negligible when sufficient sequencing depth was achieved.

CONCLUSION

Sufficient sequence depth needs to be planned to match the sequence divergence between C57BL/6 J and the strain to be studied, when the C57BL/6 J-based Agilent SureSelect exome capture system is to be used.

A Y328C missense mutation in spermine synthase causes a mild form of Snyder-Robinson syndrome

Snyder-Robinson syndrome (SRS, OMIM: 309583) is an X-linked intellectual disability (XLID) syndrome, characterized by a collection of clinical features including facial asymmetry, marfanoid habitus, hypertonia, osteoporosis and unsteady gait. It is caused by a significant decrease or loss of spermine synthase (SMS) activity. Here, we report a new missense mutation, p.Y328C (c.1084A>G), in SMS in a family with XLID. The affected males available for evaluation had mild ID, speech and global delay, an asthenic build, short stature with long fingers and mild kyphosis. The spermine/spermidine ratio in lymphoblasts was 0.53, significantly reduced compared with normal (1.87 average). Activity analysis of SMS in the index patient failed to detect any activity above background. In silico modeling demonstrated that the Y328C mutation has a significant effect on SMS stability, resulting in decreased folding free energy and larger structural fluctuations compared with those of wild-type SMS. The loss of activity was attributed to the increase in conformational dynamics in the mutant which affects the active site geometry, rather than preventing dimer formation. Taken together, the biochemical and in silico studies confirm the p.Y328C mutation in SMS is responsible for the patients having a mild form of SRS and reveal yet another molecular mechanism resulting in a non-functional SMS causing SRS.

ZC4H2 mutations are associated with arthrogryposis multiplex congenita and intellectual disability through impairment of central and peripheral synaptic plasticity

Arthrogryposis multiplex congenita (AMC) is caused by heterogeneous pathologies leading to multiple antenatal joint contractures through fetal akinesia. Understanding the pathophysiology of this disorder is important for clinical care of the affected individuals and genetic counseling of the families. We thus aimed to establish the genetic basis of an AMC subtype that is associated with multiple dysmorphic features and intellectual disability (ID). We used haplotype analysis, next-generation sequencing, array comparative genomic hybridization, and chromosome breakpoint mapping to identify the pathogenic mutations in families and simplex cases. Suspected disease variants were verified by cosegregation analysis. We identified disease-causing mutations in the zinc-finger gene ZC4H2 in four families affected by X-linked AMC plus ID and one family affected by cerebral palsy. Several heterozygous females were also affected, but to a lesser degree. Furthermore, we found two ZC4H2 deletions and one rearrangement in two female and one male unrelated simplex cases, respectively. In mouse primary hippocampal neurons, transiently produced ZC4H2 localized to the postsynaptic compartment of excitatory synapses, and the altered protein influenced dendritic spine density. In zebrafish, antisense-morpholino-mediated zc4h2 knockdown caused abnormal swimming and impaired α-motoneuron development. All missense mutations identified herein failed to rescue the swimming defect of zebrafish morphants. We conclude that ZC4H2 point mutations, rearrangements, and small deletions cause a clinically variable broad-spectrum neurodevelopmental disorder of the central and peripheral nervous systems in both familial and simplex cases of both sexes. Our results highlight the importance of ZC4H2 for genetic testing of individuals presenting with ID plus muscle weakness and minor or major forms of AMC.

Loss of function of KIAA2022 causes mild to severe intellectual disability with an autism spectrum disorder and impairs neurite outgrowth

Existence of a discrete new X-linked intellectual disability (XLID) syndrome due to KIAA2022 deficiency was questioned by disruption of KIAA2022 by an X-chromosome pericentric inversion in a XLID family we reported in 2004. Three additional families with likely pathogenic KIAA2022 mutations were discovered within the frame of systematic parallel sequencing of familial cases of XLID or in the context of routine array-CGH evaluation of sporadic intellectual deficiency (ID) cases. The c.186delC and c.3597dupA KIAA2022 truncating mutations were identified by X-chromosome exome sequencing, while array CGH discovered a 70 kb microduplication encompassing KIAA2022 exon 1 in the third family. This duplication decreased KIAA2022 mRNA level in patients' lymphocytes by 60%. Detailed clinical examination of all patients, including the two initially reported, indicated moderate-to-severe ID with autistic features, strabismus in all patients, with no specific dysmorphic features other than a round face in infancy and no structural brain abnormalities on magnetic resonance imaging (MRI). Interestingly, the patient with decreased KIAA2022 expression had only mild ID with severe language delay and repetitive behaviors falling in the range of an autism spectrum disorder (ASD). Since little is known about KIAA2022 function, we conducted morphometric studies in cultured rat hippocampal neurons. We found that siRNA-mediated KIAA2022 knockdown resulted in marked impairment in neurite outgrowth including both the dendrites and the axons, suggesting a major role for KIAA2022 in neuron development and brain function.

Comprehensive mutation analysis for congenital muscular dystrophy: a clinical PCR-based enrichment and next-generation sequencing panel

The congenital muscular dystrophies (CMDs) comprise a heterogeneous group of heritable muscle disorders with often difficult to interpret muscle pathology, making them challenging to diagnose. Serial Sanger sequencing of suspected CMD genes, while the current molecular diagnostic method of choice, can be slow and expensive. A comprehensive panel test for simultaneous screening of mutations in all known CMD-associated genes would be a more effective diagnostic strategy. Thus, the CMDs are a model disorder group for development and validation of next-generation sequencing (NGS) strategies for diagnostic and clinical care applications. Using a highly multiplexed PCR-based target enrichment method (RainDance) in conjunction with NGS, we performed mutation detection in all CMD genes of 26 samples and compared the results with Sanger sequencing. The RainDance NGS panel showed great consistency in coverage depth, on-target efficiency, versatility of mutation detection, and genotype concordance with Sanger sequencing, demonstrating the test's appropriateness for clinical use. Compared to single tests, a higher diagnostic yield was observed by panel implementation. The panel's limitation is the amplification failure of select gene-specific exons which require Sanger sequencing for test completion. Successful validation and application of the CMD NGS panel to improve the diagnostic yield in a clinical laboratory was shown.

Clinical application of massively parallel sequencing in the molecular diagnosis of glycogen storage diseases of genetically heterogeneous origin

PURPOSE

Glycogen storage diseases are a group of inborn errors of glycogen synthesis or catabolism. The outcome for untreated patients can be devastating. Given the genetic heterogeneity and the limited availability of enzyme study data, the definitive diagnosis of glycogen storage diseases is made on the basis of sequence analysis of selected potentially causative genes.

METHODS

A massively parallel sequencing test was developed for simultaneous sequencing of 16 genes known to cause muscle and liver forms of glycogen storage diseases: GYS2, GYS1, G6PC, SLC37A4, GAA, AGL, GBE1, PYGM, PYGL, PFKM, PHKA2, PHKB, PHKG2, PHKA1, PGAM2, and PGM1. All the nucleotides in the coding regions of these 16 genes have been enriched with sufficient coverage in an unbiased manner.

RESULTS

Massively parallel sequencing demonstrated 100% sensitivity and specificity as compared with Sanger sequencing. Massively parallel sequencing correctly identified all types of mutations, including single-nucleotide substitutions, small deletions and duplications, and large deletions involving one or more exons. In addition, we have confirmed the molecular diagnosis in 11 of 17 patients in whom glycogen storage diseases were suspected.

CONCLUSION

This report demonstrates the clinical utility of massively parallel sequencing technology in the diagnostic testing of a group of clinically and genetically heterogeneous disorders such as glycogen storage diseases, in a cost- and time-efficient manner.

Estimation of sequencing error rates in short reads

BACKGROUND

Short-read data from next-generation sequencing technologies are now being generated across a range of research projects. The fidelity of this data can be affected by several factors and it is important to have simple and reliable approaches for monitoring it at the level of individual experiments.

RESULTS

We developed a fast, scalable and accurate approach to estimating error rates in short reads, which has the added advantage of not requiring a reference genome. We build on the fundamental observation that there is a linear relationship between the copy number for a given read and the number of erroneous reads that differ from the read of interest by one or two bases. The slope of this relationship can be transformed to give an estimate of the error rate, both by read and by position. We present simulation studies as well as analyses of real data sets illustrating the precision and accuracy of this method, and we show that it is more accurate than alternatives that count the difference between the sample of interest and a reference genome. We show how this methodology led to the detection of mutations in the genome of the PhiX strain used for calibration of Illumina data. The proposed method is implemented in an R package, which can be downloaded from http://bcb.dfci.harvard.edu/∼vwang/shadowRegression.html.

CONCLUSIONS

The proposed method can be used to monitor the quality of sequencing pipelines at the level of individual experiments without the use of reference genomes. Furthermore, having an estimate of the error rates gives one the opportunity to improve analyses and inferences in many applications of next-generation sequencing data.

ColoSeq provides comprehensive lynch and polyposis syndrome mutational analysis using massively parallel sequencing

Lynch syndrome (hereditary nonpolyposis colon cancer) and adenomatous polyposis syndromes frequently have overlapping clinical features. Current approaches for molecular genetic testing are often stepwise, taking a best-candidate gene approach with testing of additional genes if initial results are negative. We report a comprehensive assay called ColoSeq that detects all classes of mutations in Lynch and polyposis syndrome genes using targeted capture and massively parallel next-generation sequencing on the Illumina HiSeq2000 instrument. In blinded specimens and colon cancer cell lines with defined mutations, ColoSeq correctly identified 28/28 (100%) pathogenic mutations in MLH1, MSH2, MSH6, PMS2, EPCAM, APC, and MUTYH, including single nucleotide variants (SNVs), small insertions and deletions, and large copy number variants. There was 100% reproducibility of detection mutation between independent runs. The assay correctly identified 222 of 224 heterozygous SNVs (99.4%) in HapMap samples, demonstrating high sensitivity of calling all variants across each captured gene. Average coverage was greater than 320 reads per base pair when the maximum of 96 index samples with barcodes were pooled. In a specificity study of 19 control patients without cancer from different ethnic backgrounds, we did not find any pathogenic mutations but detected two variants of uncertain significance. ColoSeq offers a powerful, cost-effective means of genetic testing for Lynch and polyposis syndromes that eliminates the need for stepwise testing and multiple follow-up clinical visits.

On the future of genetic risk assessment

Next-generation sequencing (NGS) techniques have greatly accelerated the molecular elucidation of Mendelian disorders, and affordable NGS-based diagnostic tests are around the corner that promise to detect or rule out mutations in specific subsets of the known disease genes. Whole exome sequencing and shortly afterwards whole genome sequencing (WGS) will become an even more comprehensive alternative to such targeted tests. In view of the current enthusiasm to implement these methods, but also given their rapidly dropping costs, it is quite possible that WGS will soon be adopted as universal intake test in Clinical Genetics. Central databases and large-scale genotype-phenotype comparisons will be required to progressively identify the clinically relevant sequence variants and to distinguish them from neutral polymorphisms in the human genome, and these databases will become indispensable for the interpretation of individual genome sequences. In this scenario, there will be massively growing demand for genetic counselling, but the need for experienced syndromologists will not increase proportionally, as the success of the diagnostic process will become far less dependent on the ability of clinical geneticists to reliably recognize genetic syndromes.

Assessment of target enrichment platforms using massively parallel sequencing for the mutation detection for congenital muscular dystrophy

Sequencing individual genes by Sanger sequencing is a time-consuming and costly approach to resolve clinically heterogeneous genetic disorders. Panel testing offers the ability to efficiently and cost-effectively screen all of the genes for a particular genetic disorder. We assessed the analytical sensitivity and specificity of two different enrichment technologies, solution-based hybridization and microdroplet-based PCR target enrichment, in conjunction with next-generation sequencing (NGS), to identify mutations in 321 exons representing 12 different genes involved with congenital muscular dystrophies. Congenital muscular dystrophies present diagnostic challenges due to phenotypic variability, lack of standard access to and inherent difficulties with muscle immunohistochemical stains, and a general lack of clinician awareness. NGS results were analyzed across several parameters, including sequencing metrics and genotype concordance with Sanger sequencing. Genotyping data showed that both enrichment technologies produced suitable calls for use in clinical laboratories. However, microdroplet-based PCR target enrichment is more appropriate for a clinical laboratory, due to excellent sequence specificity and uniformity, reproducibility, high coverage of the target exons, and the ability to distinguish the active gene versus known pseudogenes. Regardless of the method, exons with highly repetitive and high GC regions are not well enriched and require Sanger sequencing for completeness. Our study demonstrates the successful application of targeted sequencing in conjunction with NGS to screen for mutations in hundreds of exons in a genetically heterogeneous human disorder.

Next generation sequencing: chemistry, technology and applications

High-throughput next generation sequencing (NGS) has been quickly adapted into many aspects of biomedical research and begun to engage with the clinical practice. The latter aspect will enable the application of genomic knowledge into clinical practice in this and next decades and will profoundly change the diagnosis, prognosis and treatment of many human diseases. It will further demand both philosophical and medical curriculum reforms in the training of our future physicians. However, significant huddles need to be overcome before an ultimate application of NGS in genomic medicine can be practical and fruitful.

Novel GDI1 mutation in a large family with nonsyndromic X-linked intellectual disability

X-linked intellectual disability (XLID) is a heterogeneous disorder, and mutations in more than 90 genes have been associated with XLID to date. We report on a large multi-generational German family in which the affected male family members had nonsyndromic intellectual disability, that is, they had neither abnormal body measurements nor any other significant clinical problems. Molecular genetic analysis revealed a frameshift mutation in GDI1 (c.1185_1186delAG; Ser396ProfsX15) that co-segregated with the disease. GDI1 encodes for the GDP-dissociation inhibitor alpha (αGDI), a protein involved in the regulation of the activity of Rab GTPases. Only three families with GDI1 mutations have been reported so far. The present family supports the lack of additional phenotypic features in patients with GDI1 mutations, rendering a clinical diagnosis of GDI1-associated XLID impossible. Thus, this family not only broadens the spectrum of GDI1 mutations but also emphasizes the need for parallel testing of all known genes associated with ID in patients with an unspecific phenotype.

Carrier testing for severe childhood recessive diseases by next-generation sequencing

Of 7028 disorders with suspected Mendelian inheritance, 1139 are recessive and have an established molecular basis. Although individually uncommon, Mendelian diseases collectively account for ~20% of infant mortality and ~10% of pediatric hospitalizations. Preconception screening, together with genetic counseling of carriers, has resulted in remarkable declines in the incidence of several severe recessive diseases including Tay-Sachs disease and cystic fibrosis. However, extension of preconception screening to most severe disease genes has hitherto been impractical. Here, we report a preconception carrier screen for 448 severe recessive childhood diseases. Rather than costly, complete sequencing of the human genome, 7717 regions from 437 target genes were enriched by hybrid capture or microdroplet polymerase chain reaction, sequenced by next-generation sequencing (NGS) to a depth of up to 2.7 gigabases, and assessed with stringent bioinformatic filters. At a resultant 160x average target coverage, 93% of nucleotides had at least 20x coverage, and mutation detection/genotyping had ~95% sensitivity and ~100% specificity for substitution, insertion/deletion, splicing, and gross deletion mutations and single-nucleotide polymorphisms. In 104 unrelated DNA samples, the average genomic carrier burden for severe pediatric recessive mutations was 2.8 and ranged from 0 to 7. The distribution of mutations among sequenced samples appeared random. Twenty-seven percent of mutations cited in the literature were found to be common polymorphisms or misannotated, underscoring the need for better mutation databases as part of a comprehensive carrier testing strategy. Given the magnitude of carrier burden and the lower cost of testing compared to treating these conditions, carrier screening by NGS made available to the general population may be an economical way to reduce the incidence of and ameliorate suffering associated with severe recessive childhood disorders.

Next generation sequencing in a family with autosomal recessive Kahrizi syndrome (OMIM 612713) reveals a homozygous frameshift mutation in SRD5A3

As part of a large-scale, systematic effort to unravel the molecular causes of autosomal recessive mental retardation, we have previously described a novel syndrome consisting of mental retardation, coloboma, cataract and kyphosis (Kahrizi syndrome, OMIM 612713) and mapped the underlying gene to a 10.4-Mb interval near the centromere on chromosome 4. By combining array-based exon enrichment and next generation sequencing, we have now identified a homozygous frameshift mutation (c.203dupC; p.Phe69LeufsX2) in the gene for steroid 5α-reductase type 3 (SRD5A3) as the disease-causing change in this interval. Recent evidence indicates that this enzyme is required for the conversion of polyprenol to dolichol, a step that is essential for N-linked protein glycosylation. Independently, another group has recently observed SRD5A3 mutations in several families with a type 1 congenital disorder of glycosylation (CDG type Ix, OMIM 212067), mental retardation, cerebellar ataxia and eye disorders. Our results show that Kahrizi syndrome and this CDG Ix subtype are allelic disorders, and they illustrate the potential of next-generation sequencing strategies for the elucidation of single gene defects.