Phenotypic expansion illuminates multilocus pathogenic variation

Biallelic Pathogenic Variants in TNNT3 Associated With Congenital Myopathy

OBJECTIVE

Pathogenic variants in TNNT3, the gene encoding fast skeletal muscle troponin T, were first described in autosomal dominant distal arthrogryposis type 2B2. Recently, a homozygous splice site variant, c.681+1G>A, was identified in a patient with nemaline myopathy and distal arthrogryposis. Here, we describe the second individual with congenital myopathy associated with biallelic TNNT3 variants.

METHODS

Clinical exome sequencing data from a patient with molecularly undiagnosed congenital myopathy underwent research reanalysis. Clinical and histopathologic data were collected and compared with the single reported patient with TNNT3-related congenital myopathy.

RESULTS

A homozygous TNNT3 variant, c.481-1G>A, was identified. This variant alters a consensus splice acceptor and is predicted to affect splicing by multiple in silico prediction tools. Both the patient reported here and the previously published patient exhibited limb, bulbar, and respiratory muscle weakness from birth, which improved over time. Other shared features include history of polyhydramnios, hypotonia, scoliosis, and high-arched palate. Distal arthrogryposis and nemaline rods, findings reported in the first patient with TNNT3-related congenital myopathy, were not observed in the patient reported here.

CONCLUSIONS

This report provides further evidence for the association of biallelic TNNT3 variants with severe recessive congenital myopathy with or without nemaline rods and distal arthrogryposis. TNNT3 sequencing and copy number analysis should be incorporated into the workup of congenital myopathies.

A novel homozygous SLC13A5 whole-gene deletion generated by Alu/Alu-mediated rearrangement in an Iraqi family with epileptic encephalopathy

Biallelic loss-of-function (LoF) of SLC13A5 (solute carrier family 13, member 5) induced deficiency in sodium/citrate transporter (NaCT) causes autosomal recessive developmental epileptic encephalopathy 25 with hypoplastic amelogenesis imperfecta (DEE25; MIM #615905). Many pathogenic SLC13A5 single nucleotide variants (SNVs) and small indels have been described; however, no cases with copy number variants (CNVs) have been sufficiently investigated. We describe a consanguineous Iraqi family harboring an 88.5 kb homozygous deletion including SLC13A5 in Chr17p13.1. The three affected male siblings exhibit neonatal-onset epilepsy with fever-sensitivity, recurrent status epilepticus, global developmental delay/intellectual disability (GDD/ID), and other variable neurological findings as shared phenotypical features of DEE25. Two of the three affected subjects exhibit hypoplastic amelogenesis imperfecta (AI), while the proband shows no evidence of dental abnormalities or AI at 2 years of age with apparently unaffected primary dentition. Characterization of the genomic architecture at this locus revealed evidence for genomic instability generated by an Alu/Alu-mediated rearrangement; confirmed by break-point junction Sanger sequencing. This multiplex family from a distinct population elucidates the phenotypic consequence of complete LoF of SLC13A5 and illustrates the importance of read-depth-based CNV detection in comprehensive exome sequencing analysis to solve cases that otherwise remain molecularly unsolved.

Recessive multiple epiphyseal dysplasia and Stargardt disease in two sisters

BACKGROUND

The rapid spread of genome-wide next-generation sequencing in the molecular diagnosis of rare genetic disorders has produced increasing evidence of multilocus genomic variations in cases with a previously well-characterized molecular diagnosis. Here, we describe two patients with a rare combination of skeletal abnormalities and retinal dystrophy caused by variants in the SLC26A2 and ABCA4 genes, respectively, in a family with parental consanguinity.

METHODS

Next-generation sequencing and Sanger sequencing were performed to obtain a molecular diagnosis for the retinal and skeletal phenotypes, respectively.

RESULTS

Genetic testing revealed that the sisters were homozygous for the p.(Cys653Ser) variant in SLC26A2 and heterozygous for the missense p.(Pro68Leu) and splice donor c.6386+2C>G variants in ABCA4. Segregation analysis confirmed the carrier status of the parents.

CONCLUSION

Despite low frequency of occurrence, the detection of multilocus genomic variations in a single disease gene-oriented approach can provide accurate diagnosis even in cases with high phenotypic complexity. A targeted sequencing approach can detect relationships between observed phenotypes and underlying genotypes, useful for clinical management.

Is Gene-Size an Issue for the Diagnosis of Skeletal Muscle Disorders?

Human genes have a variable length. Those having a coding sequence of extraordinary length and a high number of exons were almost impossible to sequence using the traditional Sanger-based gene-by-gene approach. High-throughput sequencing has partly overcome the size-related technical issues, enabling a straightforward, rapid and relatively inexpensive analysis of large genes.

Several large genes (e.g. TTN, NEB, RYR1, DMD) are recognized as disease-causing in patients with skeletal muscle diseases. However, because of their sheer size, the clinical interpretation of variants in these genes is probably the most challenging aspect of the high-throughput genetic investigation in the field of skeletal muscle diseases.

The main aim of this review is to discuss the technical and interpretative issues related to the diagnostic investigation of large genes and to reflect upon the current state of the art and the future advancements in the field.

Rare and low frequency genomic variants impacting neuronal functions modify the Dup7q11.23 phenotype

BACKGROUND

7q11.23 duplication (Dup7) is one of the most frequent recurrent copy number variants (CNVs) in individuals with autism spectrum disorder (ASD), but based on gold-standard assessments, only 19% of Dup7 carriers have ASD, suggesting that additional genetic factors are necessary to manifest the ASD phenotype. To assess the contribution of additional genetic variants to the Dup7 phenotype, we conducted whole-genome sequencing analysis of 20 Dup7 carriers: nine with ASD (Dup7-ASD) and 11 without ASD (Dup7-non-ASD).

RESULTS

We identified three rare variants of potential clinical relevance for ASD: a 1q21.1 microdeletion (Dup7-non-ASD) and two deletions which disrupted IMMP2L (one Dup7-ASD, one Dup7-non-ASD). There were no significant differences in gene-set or pathway variant burden between the Dup7-ASD and Dup7-non-ASD groups. However, overall intellectual ability negatively correlated with the number of rare loss-of-function variants present in nervous system development and membrane component pathways, and adaptive behaviour standard scores negatively correlated with the number of low-frequency likely-damaging missense variants found in genes expressed in the prenatal human brain. ASD severity positively correlated with the number of low frequency loss-of-function variants impacting genes expressed at low levels in the brain, and genes with a low level of intolerance.

CONCLUSIONS

Our study suggests that in the presence of the same pathogenic Dup7 variant, rare and low frequency genetic variants act additively to contribute to components of the overall Dup7 phenotype.

Targeted Next-Generation Sequencing Indicates a Frequent Oligogenic Involvement in Primary Ovarian Insufficiency Onset

Primary ovarian insufficiency (POI) is one of the major causes of female infertility associated with the premature loss of ovarian function in about 3.7% of women before the age of 40. This disorder is highly heterogeneous and can manifest with a wide range of clinical phenotypes, ranging from ovarian dysgenesis and primary amenorrhea to post-pubertal secondary amenorrhea, with elevated serum gonadotropins and hypoestrogenism. The ovarian defect still remains idiopathic in some cases; however, a strong genetic component has been demonstrated by the next-generation sequencing (NGS) approach of familiar and sporadic POI cases. As recent evidence suggested an oligogenic architecture for POI, we developed a target NGS panel with 295 genes including known candidates and novel genetic determinants potentially involved in POI pathogenesis. Sixty-four patients with early onset POI (range: 10-25 years) of our cohort have been screened with 90% of target coverage at 50×. Here, we report 48 analyzed patients with at least one genetic variant (75%) in the selected candidate genes. In particular, we found the following: 11/64 patients (17%) with two variants, 9/64 (14%) with three variants, 9/64 (14%) with four variants, 3/64 (5%) with five variants, and 2/64 (3%) with six variants. The most severe phenotypes were associated with either the major number of variations or a worse prediction in pathogenicity of variants. Bioinformatic gene ontology analysis identified the following major pathways likely affected by gene variants: 1) cell cycle, meiosis, and DNA repair; 2) extracellular matrix remodeling; 3) reproduction; 4) cell metabolism; 5) cell proliferation; 6) calcium homeostasis; 7) NOTCH signaling; 8) signal transduction; 9) WNT signaling; 10) cell death; and 11) ubiquitin modifications. Consistently, the identified pathways have been described in other studies dissecting the mechanisms of folliculogenesis in animal models of altered fertility. In conclusion, our results contribute to define POI as an oligogenic disease and suggest novel candidates to be investigated in patients with POI.

Genetic Testing in Various Neurodevelopmental Disorders Which Manifest as Cerebral Palsy: A Case Study From Iran

Purpose: Cerebral palsy (CP) is a heterogeneous permanent disorder impacting movement and posture. Investigations aimed at diagnosing this disorder are expensive and time-consuming and can eventually inconclusive. This study aimed to determine the diagnostic yield of next generation sequencing in patients with atypical CP (ACP). Methods: Patient eligibility criteria included impaired motor function with onset at birth or within the first year of life, and one or more of the following conditions: severe intellectual disability, positive family history, brain imaging findings not typical for cerebral palsy, abnormal neurometabolic profile, intractable seizure, normal neuroimaging despite severe psychomotor disability, after pediatric neurologist assessment including neuroimaging and biochemical-metabolic study offered for genetic study. Results: Exome sequencing was done for 66 patients which revealed pathogenic, likely pathogenic, and variants of unknown significance in 36.2, 9, and 43.9%, respectively. We also found 10 new mutations and were able to suggest specific and personalized treatments for nine patients. We also found three different mutations with different phenotypical spectrum in one gene that have not been reported for cerebral palsy. Conclusion: An accurate history and physical examination and determination of patients with atypical cerebral palsy for doing exome sequencing result in improved genetic counseling and personalized management.

Comprehensive clinical, biochemical, radiological and genetic analysis of 28 Turkish cases with suspected metachromatic leukodystrophy and their relatives

Metachromatic leukodystrophy (MLD) is a glycosphingolipid storage disease caused by deficiency of the lysosomal enzyme arylsulfatase A (ASA) or its activator protein saposin B. MLD can affect all age groups in severity varying from a severe fatal form to milder adult onset forms. Diagnosis is usually made by measuring leukocyte ASA activity. However, this test can give false negative or false positive laboratory results due to pseudodeficiency of ASA and saposin B deficiency, respectively. Therefore, we aimed to evaluate patients with suspected MLD in a Turkish population by comprehensive clinical, biochemical, radiological, and genetic analyses for molecular and phenotypic characterization. We analyzed 28 suspected MLD patients and 41 relatives from 24 families. ASA activity was found to be decreased in 21 of 28 patients. Sixteen patients were diagnosed as MLD (11 late infantile, 2 juvenile and 3 adult types), 2 MSD, 2 pseudodeficiency (PD) and the remaining 8 patients were diagnosed as having other leukodystrophies. Enzyme analysis showed that the age of onset of MLD did not correlate with residual ASA activity. Sequence analysis showed 11 mutations in ARSA, of which 4 were novel (p.Trp195GlyfsTer5, p.Gly298Asp, p.Arg301Leu, and p.Gly311Asp), and 2 mutations in SUMF1 causing multiple sulfatase deficiency, and confirmed the diagnosis of MLD in 2 presymptomatic relatives. All individuals with confirmed mutations had low ASA activity and urinary sulfatide excretion. Intra- and inter-familial variability was high for the same ARSA missense genotypes, indicating the contribution of other factors to disease expression. Imaging findings were evaluated through a modified brain MRI scoring system which indicated patients with protein-truncating mutations had more severe MRI findings and late-infantile disease onset. MRI findings were not specific for the diagnosis. Anti-sulfatide IgM was similar to control subjects, and IgG, elevated in multiple sulfatase deficiency. In conclusion, the knowledge on the biochemical, clinical and genetic basis of MLD was expanded, a modified diagnostic laboratory algorithm for MLD based on integrated evaluation of ASA activity, urinary sulfatide excretion and genetic tests was devised.

The role of acid-labile subunit (ALS) in the modulation of GH-IGF-I action

Acid-labile subunit (ALS) deficiency (ACLSD) constitutes the first monogenic defect involving a member of the Insulin-like Growth Factor (IGF) binding protein system. The lack of ALS completely disrupts the circulating IGF system. Autocrine/paracrine action of local produced IGF-I could explain the mild effect on growth. In the present work we have revised the more relevant clinical and biochemical consequences of complete ACLSD in 61 reported subjects from 31 families. Low birth weight and/or length, reduced head circumference, height between -2 and -3 SD, pubertal delay and insulin resistance are commonly observed. Partial ACLSD could be present in children initially labeled as idiopathic short stature, presenting low IGF-I levels, suggesting that one functional IGFALS allele is insufficient to stabilize ternary complexes. Dysfunction of the GH-IGF axis observed in ACLSD may eventually result in increased risk for type-2 diabetes and tumor progression. Consequently, long term surveillance is recommended in these patients.

Parental exome analysis identifies shared carrier status for a second recessive disorder in couples with an affected child

Consanguinity, commonplace in many regions around the globe, is associated with an increased risk of autosomal recessive (AR) genetic disorders. Consequently, consanguineous couples undergoing preimplantation genetic diagnosis (PGD) for one Mendelian disorder may be at increased risk for a child with a second, unrelated AR genetic disorder. We examined the yield of exome analysis for carrier screening of additional AR disorders, beyond the primary diagnosis, amongst consanguineous vs. non-consanguineous populations. Parental samples from trio exomes of 102 consanguineous families and 105 non-consanguineous controls were evaluated for shared carrier status, after disregarding the primary molecular diagnosis. Results were sub-classified according to disease severity. Secondary shared carrier status for pathogenic and likely pathogenic variants leading to AR disorders of moderate to profound severity was identified in 10/102 (9.8%) consanguineous couples, as compared to 1/105 (0.95%) non-consanguineous couples (χ² = 8.0565, p value < 0.005). Higher inbreeding coefficient values, calculated from individual exomes, correlated with secondary shared carrier status for diseases of moderate to profound severity (r = 0.17, p value < 0.0125). Our results indicate that consanguineous couples undergoing PGD are at increased risk for a second genetic disease of moderate to profound severity. This study represents an underestimate of the rate of secondary shared carrier status due to inability to detect deep intronic variants, no assessment of copy number variants, and false negative results stemming from stringent variant interpretation. False positive results may result from inaccuracies in public databases. Additional studies in consanguineous populations will determine whether exome-based carrier screening should be recommended to all couples undergoing PGD.

De novo SMARCA2 variants clustered outside the helicase domain cause a new recognizable syndrome with intellectual disability and blepharophimosis distinct from Nicolaides-Baraitser syndrome

PURPOSE

Nontruncating variants in SMARCA2, encoding a catalytic subunit of SWI/SNF chromatin remodeling complex, cause Nicolaides-Baraitser syndrome (NCBRS), a condition with intellectual disability and multiple congenital anomalies. Other disorders due to SMARCA2 are unknown.

METHODS

By next-generation sequencing, we identified candidate variants in SMARCA2 in 20 individuals from 18 families with a syndromic neurodevelopmental disorder not consistent with NCBRS. To stratify variant interpretation, we functionally analyzed SMARCA2 variants in yeasts and performed transcriptomic and genome methylation analyses on blood leukocytes.

RESULTS

Of 20 individuals, 14 showed a recognizable phenotype with recurrent features including epicanthal folds, blepharophimosis, and downturned nasal tip along with variable degree of intellectual disability (or blepharophimosis intellectual disability syndrome [BIS]). In contrast to most NCBRS variants, all SMARCA2 variants associated with BIS are localized outside the helicase domains. Yeast phenotype assays differentiated NCBRS from non-NCBRS SMARCA2 variants. Transcriptomic and DNA methylation signatures differentiated NCBRS from BIS and those with nonspecific phenotype. In the remaining six individuals with nonspecific dysmorphic features, clinical and molecular data did not permit variant reclassification.

CONCLUSION

We identified a novel recognizable syndrome named BIS associated with clustered de novo SMARCA2 variants outside the helicase domains, phenotypically and molecularly distinct from NCBRS.

Update of the genotype and phenotype of KMT2D and KDM6A by genetic screening of 100 patients with clinically suspected Kabuki syndrome

Kabuki syndrome is characterized by a variable degree of intellectual disability, characteristic facial features, and complications in various organs. Many variants have been identified in two causative genes, that is, lysine methyltransferase 2D (KMT2D) and lysine demethylase 6A (KDM6A). In this study, we present the results of genetic screening of 100 patients with a suspected diagnosis of Kabuki syndrome in our center from July 2010 to June 2018. We identified 76 variants (43 novel) in KMT2D and 4 variants (3 novel) in KDM6A as pathogenic or likely pathogenic. Rare variants included a deep splicing variant (c.14000-8C>G) confirmed by RNA sequencing and an 18% mosaicism level for a KMT2D mutation. We also characterized a case with a blended phenotype consisting of Kabuki syndrome, osteogenesis imperfecta, and 16p13.11 microdeletion. We summarized the clinical phenotypes of 44 patients including a patient who developed cervical cancer of unknown origin at 16 years of age. This study presents important details of patients with Kabuki syndrome including rare clinical cases and expands our genetic understanding of this syndrome, which will help clinicians and researchers better manage and understand patients with Kabuki syndrome they may encounter.

Novel Compound Heterozygous Mutations in CRTAP Cause Rare Autosomal Recessive Osteogenesis Imperfecta

Whole-exome sequencing (WES) has advantages over the traditional molecular test by screening 20,000 genes simultaneously and has become an invaluable tool for genetic diagnosis in clinical practice. Here, we reported a family with a child and a fetus presenting undiagnosed skeletal dysplasia phenotypes, while the parents were asymptomatic. WES was applied to the parents and affected fetus to identify the genetic cause of the phenotypes. We identified novel compound heterozygous mutations consisting of a single-nucleotide variant (SNV) and a large deletion in the CRTAP gene (NM_006371.4:c.1153-3C > G/hg19 chr3:g.32398837_34210906del). Genetic alterations of CRTAP are known to cause osteogenesis imperfecta (OI) in an autosomal recessive manner. Further examination of the proband’s elder sibling who was diagnosed as OI after birth found that she shares the inherited compound heterozygous mutations of CRTAP; thus, the findings support the disease-causing role of CRTAP mutations. Through the in vitro molecular test and in silico analysis, the deleterious effects of the splicing-altering SNV in CRTAP (c.1153-3C > G) on gene product were confirmed. Collectively, our WES-based pathogenic variant discovery pipeline identifies the SNVs and copy number variation to delineate the genetic cause on the proband affected with OI. The data not only extend the knowledge of mutation spectrum in patients with skeletal dysplasia but also demonstrate that WES holds great promise for genetic screening of rare diseases in clinical settings.

Neurodevelopmental Disorders: From Genetics to Functional Pathways

Neurodevelopmental disorders (NDDs) are a class of disorders affecting brain development and function and are characterized by wide genetic and clinical variability. In this review, we discuss the multiple factors that influence the clinical presentation of NDDs, with particular attention to gene vulnerability, mutational load, and the two-hit model. Despite the complex architecture of mutational events associated with NDDs, the various proteins involved appear to converge on common pathways, such as synaptic plasticity/function, chromatin remodelers and the mammalian target of rapamycin (mTOR) pathway. A thorough understanding of the mechanisms behind these pathways will hopefully lead to the identification of candidates that could be targeted for treatment approaches.

Functional biology of the Steel syndrome founder allele and evidence for clan genomics derivation of COL27A1 pathogenic alleles worldwide

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

Biallelic GRM7 variants cause epilepsy, microcephaly, and cerebral atrophy

Objective

Defects in ion channels and neurotransmitter receptors are implicated in developmental and epileptic encephalopathy (DEE). Metabotropic glutamate receptor 7 (mGluR7), encoded by GRM7, is a presynaptic G‐protein‐coupled glutamate receptor critical for synaptic transmission. We previously proposed GRM7 as a candidate disease gene in two families with neurodevelopmental disorders (NDDs). One additional family has been published since. Here, we describe three additional families with GRM7 biallelic variants and deeply characterize the associated clinical neurological and electrophysiological phenotype and molecular data in 11 affected individuals from six unrelated families.

Methods

Exome sequencing and family‐based rare variant analyses on a cohort of 220 consanguineous families with NDDs revealed three families with GRM7 biallelic variants; three additional families were identified through literature search and collaboration with a clinical molecular laboratory.

Results

We compared the observed clinical features and variants of 11 affected individuals from the six unrelated families. Identified novel deleterious variants included two homozygous missense variants (c.2671G>A:p.Glu891Lys and c.1973G>A:p.Arg685Gln) and one homozygous stop‐gain variant (c.1975C>T:p.Arg659Ter). Developmental delay, neonatal‐ or infantile‐onset epilepsy, and microcephaly were universal. Three individuals had hypothalamic–pituitary–axis dysfunction without pituitary structural abnormality. Neuroimaging showed cerebral atrophy and hypomyelination in a majority of cases. Two siblings demonstrated progressive loss of myelination by 2 years in both and an acquired microcephaly pattern in one. Five individuals died in early or late childhood.

Conclusion

Detailed clinical characterization of 11 individuals from six unrelated families demonstrates that rare biallelic GRM7 pathogenic variants can cause DEEs, microcephaly, hypomyelination, and cerebral atrophy.

Expanding the phenotype of thrombocytopenia absent radius syndrome with hypospadias

Rare genetic diseases and syndromes may appear with unique features in some patients. In genetically-solved cases, this situation indicates a phenotypic expansion of the syndrome with additional features (i.e. the disease-associated gene gives rise to unusual clinical presentation). However, this situation can also hide a multilocus pathogenic variation that cannot be solved genetically except by a massive sequencing approach, such as exome sequencing. Here we describe the case of a child with bilateral radial aplasia, transient thrombocytopenia and anemia, cow's milk intolerance, hypospadias, facial dysmorphism, mild hypothyroidism and umbilical and inguinal hernia. Bilaterally absent radius, presence of thumbs and low platelet count are pathognomonic of thrombocytopenia absent radius (TAR) syndrome, but the child also showed other features beyond those reported in the literature. Since various diseases resembling the proband's phenotype required differential diagnosis, clinical exome sequencing was performed. The results showed compound heterozygous mutations in the RBM8A gene, confirming the suspicion of TAR syndrome. A truncating heterozygous variant in the DUOX2 gene, known to be associated with transient thyroid dyshormonogenesis type 6 (TDH6), was also detected and may explain the proband's mild hypothyroidism.

Double homozygosity in CEP57 and DYNC2H1 genes detected by WES: Composite or expanded phenotype?

BACKGROUND

In the last few years trio-whole exome sequencing (WES) analysis has demonstrated its potential in obtaining genetic diagnoses even in nonspecific clinical pictures and in atypical presentations of known diseases. Moreover WES allows the detection of variants in multiple genes causing different genetic conditions in a single patient, in about 5% of cases. The resulting phenotype may be clinically discerned as variability in the expression of a known phenotype, or as a new unreported syndromic condition.

METHODS

Trio-WES was performed on a 4-month-old baby with a complex clinical presentation characterized by skeletal anomalies, congenital heart malformation, congenital hypothyroidism, generalized venous and arterial hypoplasia, and recurrent infections.

RESULTS

WES detected two different homozygous variants, one in CEP57, the gene responsible for mosaic variegated aneuploidy syndrome 2, the other in DYNC2H1, the main gene associated with short-rib thoracic dysplasia.

CONCLUSION

The contribution of these two different genetic causes in determining the phenotype of our patient is discussed, including some clinical signs not explained by the detected variants. The report then highlights the role of WES in providing complete and fast diagnosis in patients with complex presentations of rare genetic syndromes, with important implications in the assessment of recurrence risk.

Deficiencies in vesicular transport mediated by TRAPPC4 are associated with severe syndromic intellectual disability

The conserved transport protein particle (TRAPP) complexes regulate key trafficking events and are required for autophagy. TRAPPC4, like its yeast Trs23 orthologue, is a core component of the TRAPP complexes and one of the essential subunits for guanine nucleotide exchange factor activity for Rab1 GTPase. Pathogenic variants in specific TRAPP subunits are associated with neurological disorders. We undertook exome sequencing in three unrelated families of Caucasian, Turkish and French-Canadian ethnicities with seven affected children that showed features of early-onset seizures, developmental delay, microcephaly, sensorineural deafness, spastic quadriparesis and progressive cortical and cerebellar atrophy in an effort to determine the genetic aetiology underlying neurodevelopmental disorders. All seven affected subjects shared the same identical rare, homozygous, potentially pathogenic variant in a non-canonical, well-conserved splice site within TRAPPC4 (hg19:chr11:g.118890966A>G; TRAPPC4: NM_016146.5; c.454+3A>G). Single nucleotide polymorphism array analysis revealed there was no haplotype shared between the tested Turkish and Caucasian families suggestive of a variant hotspot region rather than a founder effect. In silico analysis predicted the variant to cause aberrant splicing. Consistent with this, experimental evidence showed both a reduction in full-length transcript levels and an increase in levels of a shorter transcript missing exon 3, suggestive of an incompletely penetrant splice defect. TRAPPC4 protein levels were significantly reduced whilst levels of other TRAPP complex subunits remained unaffected. Native polyacrylamide gel electrophoresis and size exclusion chromatography demonstrated a defect in TRAPP complex assembly and/or stability. Intracellular trafficking through the Golgi using the marker protein VSVG-GFP-ts045 demonstrated significantly delayed entry into and exit from the Golgi in fibroblasts derived from one of the affected subjects. Lentiviral expression of wild-type TRAPPC4 in these fibroblasts restored trafficking, suggesting that the trafficking defect was due to reduced TRAPPC4 levels. Consistent with the recent association of the TRAPP complex with autophagy, we found that the fibroblasts had a basal autophagy defect and a delay in autophagic flux, possibly due to unsealed autophagosomes. These results were validated using a yeast trs23 temperature sensitive variant that exhibits constitutive and stress-induced autophagic defects at permissive temperature and a secretory defect at restrictive temperature. In summary we provide strong evidence for pathogenicity of this variant in a member of the core TRAPP subunit, TRAPPC4 that associates with vesicular trafficking and autophagy defects. This is the first report of a TRAPPC4 variant, and our findings add to the growing number of TRAPP-associated neurological disorders.

Mutations in ANKLE2, a ZIKA Virus Target, Disrupt an Asymmetric Cell Division Pathway in Drosophila Neuroblasts to Cause Microcephaly

The apical Par complex, which contains atypical protein kinase C (aPKC), Bazooka (Par-3), and Par-6, is required for establishing polarity during asymmetric division of neuroblasts in Drosophila, and its activity depends on L(2)gl. We show that loss of Ankle2, a protein associated with microcephaly in humans and known to interact with Zika protein NS4A, reduces brain volume in flies and impacts the function of the Par complex. Reducing Ankle2 levels disrupts endoplasmic reticulum (ER) and nuclear envelope morphology, releasing the kinase Ballchen-VRK1 into the cytosol. These defects are associated with reduced phosphorylation of aPKC, disruption of Par-complex localization, and spindle alignment defects. Importantly, removal of one copy of ballchen or l(2)gl suppresses Ankle2 mutant phenotypes and restores viability and brain size. Human mutational studies implicate the above-mentioned genes in microcephaly and motor neuron disease. We suggest that NS4A, ANKLE2, VRK1, and LLGL1 define a pathway impinging on asymmetric determinants of neural stem cell division.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

A Genocentric Approach to Discovery of Mendelian Disorders

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Bi-allelic Pathogenic Variants in TUBGCP2 Cause Microcephaly and Lissencephaly Spectrum Disorders

Lissencephaly comprises a spectrum of malformations of cortical development. This spectrum includes agyria, pachygyria, and subcortical band heterotopia; each represents anatomical malformations of brain cortical development caused by neuronal migration defects. The molecular etiologies of neuronal migration anomalies are highly enriched for genes encoding microtubules and microtubule-associated proteins, and this enrichment highlights the critical role for these genes in cortical growth and gyrification. Using exome sequencing and family based rare variant analyses, we identified a homozygous variant (c.997C>T [p.Arg333Cys]) in TUBGCP2, encoding gamma-tubulin complex protein 2 (GCP2), in two individuals from a consanguineous family; both individuals presented with microcephaly and developmental delay. GCP2 forms the multiprotein γ-tubulin ring complex (γ-TuRC) together with γ-tubulin and other GCPs to regulate the assembly of microtubules. By querying clinical exome sequencing cases and through GeneMatcher-facilitated collaborations, we found three additional families with bi-allelic variation and similarly affected phenotypes including a homozygous variant (c.1843G>C [p.Ala615Pro]) in two families and compound heterozygous variants consisting of one missense variant (c.889C>T [p.Arg297Cys]) and one splice variant (c.2025-2A>G) in another family. Brain imaging from all five affected individuals revealed varying degrees of cortical malformations including pachygyria and subcortical band heterotopia, presumably caused by disruption of neuronal migration. Our data demonstrate that pathogenic variants in TUBGCP2 cause an autosomal recessive neurodevelopmental trait consisting of a neuronal migration disorder, and our data implicate GCP2 as a core component of γ-TuRC in neuronal migrating cells.

The genomic and clinical landscape of fetal akinesia

PURPOSE

Fetal akinesia has multiple clinical subtypes with over 160 gene associations, but the genetic etiology is not yet completely understood.

METHODS

In this study, 51 patients from 47 unrelated families were analyzed using next-generation sequencing (NGS) techniques aiming to decipher the genomic landscape of fetal akinesia (FA).

RESULTS

We have identified likely pathogenic gene variants in 37 cases and report 41 novel variants. Additionally, we report putative pathogenic variants in eight cases including nine novel variants. Our work identified 14 novel disease-gene associations for fetal akinesia: ADSSL1, ASAH1, ASPM, ATP2B3, EARS2, FBLN1, PRG4, PRICKLE1, ROR2, SETBP1, SCN5A, SCN8A, and ZEB2. Furthermore, a sibling pair harbored a homozygous copy-number variant in TNNT1, an ultrarare congenital myopathy gene that has been linked to arthrogryposis via Gene Ontology analysis.

CONCLUSION

Our analysis indicates that genetic defects leading to primary skeletal muscle diseases might have been underdiagnosed, especially pathogenic variants in RYR1. We discuss three novel putative fetal akinesia genes: GCN1, IQSEC3 and RYR3. Of those, IQSEC3, and RYR3 had been proposed as neuromuscular disease-associated genes recently, and our findings endorse them as FA candidate genes. By combining NGS with deep clinical phenotyping, we achieved a 73% success rate of solved cases.

Missense variants in TAF1 and developmental phenotypes: challenges of determining pathogenicity

We recently described a new neurodevelopmental syndrome (TAF1/MRXS33 intellectual disability syndrome) (MIM# 300966) caused by pathogenic variants involving the X-linked gene TAF1, which participates in RNA polymerase II transcription. The initial study reported eleven families, and the syndrome was defined as presenting early in life with hypotonia, facial dysmorphia, and developmental delay that evolved into intellectual disability (ID) and/or autism spectrum disorder (ASD). We have now identified an additional 27 families through a genotype-first approach. Familial segregation analysis, clinical phenotyping, and bioinformatics were capitalized on to assess potential variant pathogenicity, and molecular modelling was performed for those variants falling within structurally characterized domains of TAF1. A novel phenotypic clustering approach was also applied, in which the phenotypes of affected individuals were classified using 51 standardized Human Phenotype Ontology (HPO) terms. Phenotypes associated with TAF1 variants show considerable pleiotropy and clinical variability, but prominent among previously unreported effects were brain morphological abnormalities, seizures, hearing loss, and heart malformations. Our allelic series broadens the phenotypic spectrum of TAF1/MRXS33 intellectual disability syndrome and the range of TAF1 molecular defects in humans. It also illustrates the challenges for determining the pathogenicity of inherited missense variants, particularly for genes mapping to chromosome X. This article is protected by copyright. All rights reserved.

Novel IQCE variations confirm its role in postaxial polydactyly and cause ciliary defect phenotype in zebrafish

Polydactyly is one of the most frequent inherited defects of the limbs characterized by supernumerary digits and high-genetic heterogeneity. Among the many genes involved, either in isolated or syndromic forms, eight have been implicated in postaxial polydactyly (PAP). Among those, IQCE has been recently identified in a single consanguineous family. Using whole-exome sequencing in patients with uncharacterized ciliopathies, including PAP, we identified three families with biallelic pathogenic variations in IQCE. Interestingly, the c.895_904del (p.Val301Serfs*8) was found in all families without sharing a common haplotype, suggesting a recurrent mechanism. Moreover, in two families, the systemic phenotype could be explained by additional pathogenic variants in known genes (TULP1, ATP6V1B1). RNA expression analysis on patients' fibroblasts confirms that the dysfunction of IQCE leads to the dysregulation of genes associated with the hedgehog-signaling pathway, and zebrafish experiments demonstrate a full spectrum of phenotypes linked to defective cilia: Body curvature, kidney cysts, left-right asymmetry, misdirected cilia in the pronephric duct, and retinal defects. In conclusion, we identified three additional families confirming IQCE as a nonsyndromic PAP gene. Our data emphasize the importance of taking into account the complete set of variations of each individual, as each clinical presentation could finally be explained by multiple genes.

Defective tubulin detyrosination causes structural brain abnormalities with cognitive deficiency in humans and mice

Reversible detyrosination of tubulin, the building block of microtubules, is crucial for neuronal physiology. Enzymes responsible for detyrosination were recently identified as complexes of vasohibins (VASHs) one or two with small VASH-binding protein (SVBP). Here we report three consanguineous families, each containing multiple individuals with biallelic inactivation of SVBP caused by truncating variants (p.Q28* and p.K13Nfs*18). Affected individuals show brain abnormalities with microcephaly, intellectual disability and delayed gross motor and speech development. Immunoblot testing in cells with pathogenic SVBP variants demonstrated that the encoded proteins were unstable and non-functional, resulting in a complete loss of VASH detyrosination activity. Svbp knockout mice exhibit drastic accumulation of tyrosinated tubulin and a reduction of detyrosinated tubulin in brain tissue. Similar alterations in tubulin tyrosination levels were observed in cultured neurons and associated with defects in axonal differentiation and architecture. Morphological analysis of the Svbp knockout mouse brains by anatomical magnetic resonance imaging showed a broad impact of SVBP loss, with a 7% brain volume decrease, numerous structural defects and a 30% reduction of some white matter tracts. Svbp knockout mice display behavioural defects, including mild hyperactivity, lower anxiety and impaired social behaviour. They do not, however, show prominent memory defects. Thus, SVBP-deficient mice recapitulate several features observed in human patients. Altogether, our data demonstrate that deleterious variants in SVBP cause this neurodevelopmental pathology, by leading to a major change in brain tubulin tyrosination and alteration of microtubule dynamics and neuron physiology.

Genetic modifiers and non-Mendelian aspects of CMT

Charcot-Marie-Tooth (CMT) neuropathies are amongst the most common inherited diseases in neurology. While great strides have been made to identify the genesis of these diseases, a diagnostic gap of 30-60% remains. Classic models of genetic causation may be limited to fully close this gap and, thus, we review the current state and future role of alternative, non-Mendelian forms of genetics in CMT. Promising synergies exist to further define the full genetic architecture of inherited neuropathies, including affordable whole-genome sequencing, increased data aggregation and clinical collaboration, improved bioinformatics and statistical methodology, and vastly improved computational resources. Given the recent advances in genetic therapies for rare diseases, it becomes a matter of urgency to diagnose CMT patients with great fidelity. Otherwise, they will not be able to benefit from such therapeutic options, or worse, suffer harm when pathogenicity of genetic variation is falsely evaluated. In addition, the newly identified modifier and risk genes may offer alternative targets for pharmacotherapy of inherited and, potentially, even acquired forms of neuropathies.

Homozygous variants in MAPRE2 and CDON in individual with skin folds, growth delay, retinal coloboma, and pyloric stenosis

Cases with multiple molecular diagnoses are challenging to diagnose clinically, yet may be resolved by unbiased exome sequencing analysis. We report an infant with developmental delay, severe growth delay, dysmorphic features, and multiple congenital anomalies including retinal coloboma, congenital pyloric stenosis, and circumferential skin creases. Exome sequencing identified a homozygous missense variant in MAPRE2 and a homozygous stopgain (nonsense) variant in CDON. Variants in MAPRE2, encoding a regulator of microtubule dynamics, lead to congenital symmetric circumferential skin creases type 2, with associated dysmorphism, small growth parameters, and congenital cardiac and genital anomalies. Monoallelic variants in CDON, encoding a coreceptor for sonic hedgehog, have been associated with autosomal dominant pituitary stalk interruption syndrome and holoprosencephaly. Cdon-/- mice have multiple eye defects including coloboma, consistent with the observed human phenotype. Thus, the complex phenotypic presentation of the infant may potentially be attributed to a dual molecular diagnosis. Furthermore, we present CDON as a candidate gene for coloboma formation in addition to the known holoprosencephaly phenotype, and propose to expand the allelic spectrum of CDON to variants associated with autosomal recessive inheritance in addition to dominant inheritance.

A Patient with Berardinelli-Seip Syndrome, Novel AGPAT2 Splicesite Mutation and Concomitant Development of Non-diabetic Polyneuropathy

Primary polyneuropathy in the context of Seip-Berardinelli type 1 seipinopathy, or congenital generalized lipodystrophy type 1 (CGL1) has not been previously reported. We report the case history of a 27 year old female CGL1 patient presenting with an unusual additional development of non-diabetic peripheral neuropathy and learning disabilities in early adolescence. Whole exome sequencing (WES) of the patient genome identified a novel variant, homozygous for a 52 bp intronic deletion in the AGPAT2 locus, coding for 1-acylglycerol-3-phosphate O-acyltransferase 2, which is uniquely associated with CGL1 seipinopathies, with no molecular evidence for dual diagnosis. Functional studies using RNA isolated from patient peripheral blood leucocytes showed abnormal RNA splicing resulting in the loss of 25 amino acids from the patient AGPAT2 protein coding sequence. Stability and transcription levels for the misspliced AGPAT2 mRNA in our patient nonetheless remained normal. Any AGPAT2 protein produced in our patient is therefore likely to be dysfunctional. However, formal linkage of this deletion to the neuropathy observed remains to be shown. The classical clinical presentation of a patient with AGPAT2-associated lipodystrophy shows normal cognition and no development of polyneuropathy. Cognitive disabilities and polyneuropathy are features associated exclusively with clinical CGL type 2 arising from seipin (BSCL2) gene mutations. This case study suggests that in some genetic contexts, AGPAT2 mutations can also produce phenotypes with primary polyneuropathy.

Exome Sequencing of a Primary Ovarian Insufficiency Cohort Reveals Common Molecular Etiologies for a Spectrum of Disease

CONTEXT

Primary ovarian insufficiency (POI) encompasses a spectrum of premature menopause, including both primary and secondary amenorrhea. For 75% to 90% of individuals with hypergonadotropic hypogonadism presenting as POI, the molecular etiology is unknown. Common etiologies include chromosomal abnormalities, environmental factors, and congenital disorders affecting ovarian development and function, as well as syndromic and nonsyndromic single gene disorders suggesting POI represents a complex trait.

OBJECTIVE

To characterize the contribution of known disease genes to POI and identify molecular etiologies and biological underpinnings of POI.

DESIGN, SETTING, AND PARTICIPANTS

We applied exome sequencing (ES) and family-based genomics to 42 affected female individuals from 36 unrelated Turkish families, including 31 with reported parental consanguinity.

RESULTS

This analysis identified likely damaging, potentially contributing variants and molecular diagnoses in 16 families (44%), including 11 families with likely damaging variants in known genes and five families with predicted deleterious variants in disease genes (IGSF10, MND1, MRPS22, and SOHLH1) not previously associated with POI. Of the 16 families, 2 (13%) had evidence for potentially pathogenic variants at more than one locus. Absence of heterozygosity consistent with identity-by-descent mediated recessive disease burden contributes to molecular diagnosis in 15 of 16 (94%) families. GeneMatcher allowed identification of additional families from diverse genetic backgrounds.

CONCLUSIONS

ES analysis of a POI cohort further characterized locus heterogeneity, reaffirmed the association of genes integral to meiotic recombination, demonstrated the likely contribution of genes involved in hypothalamic development, and documented multilocus pathogenic variation suggesting the potential for oligogenic inheritance contributing to the development of POI.

Paralog Studies Augment Gene Discovery: DDX and DHX Genes

Members of a paralogous gene family in which variation in one gene is known to cause disease are eight times more likely to also be associated with human disease. Recent studies have elucidated DHX30 and DDX3X as genes for which pathogenic variant alleles are involved in neurodevelopmental disorders. We hypothesized that variants in paralogous genes encoding members of the DExD/H-box RNA helicase superfamily might also underlie developmental delay and/or intellectual disability (DD and/or ID) disease phenotypes. Here we describe 15 unrelated individuals who have DD and/or ID, central nervous system (CNS) dysfunction, vertebral anomalies, and dysmorphic features and were found to have probably damaging variants in DExD/H-box RNA helicase genes. In addition, these individuals exhibit a variety of other tissue and organ system involvement including ocular, outer ear, hearing, cardiac, and kidney tissues. Five individuals with homozygous (one), compound-heterozygous (two), or de novo (two) missense variants in DHX37 were identified by exome sequencing. We identified ten total individuals with missense variants in three other DDX/DHX paralogs: DHX16 (four individuals), DDX54 (three individuals), and DHX34 (three individuals). Most identified variants are rare, predicted to be damaging, and occur at conserved amino acid residues. Taken together, these 15 individuals implicate the DExD/H-box helicases in both dominantly and recessively inherited neurodevelopmental phenotypes and highlight the potential for more than one disease mechanism underlying these disorders.

Biallelic CACNA2D2 variants in epileptic encephalopathy and cerebellar atrophy

OBJECTIVE

To characterize the molecular and clinical phenotypic basis of developmental and epileptic encephalopathies caused by rare biallelic variants in CACNA2D2.

METHODS

Two affected individuals from a family with clinical features of early onset epileptic encephalopathy were recruited for exome sequencing at the Centers for Mendelian Genomics to identify their molecular diagnosis. GeneMatcher facilitated identification of a second family with a shared candidate disease gene identified through clinical gene panel-based testing.

RESULTS

Rare biallelic CACNA2D2 variants have been previously reported in three families with developmental and epileptic encephalopathy, and one family with congenital ataxia. We identified three individuals in two unrelated families with novel homozygous rare variants in CACNA2D2 with clinical features of developmental and epileptic encephalopathy and cerebellar atrophy. Family 1 includes two affected siblings with a likely damaging homozygous rare missense variant c.1778G>C; p.(Arg593Pro) in CACNA2D2. Family 2 includes a proband with a homozygous rare nonsense variant c.485_486del; p.(Tyr162Ter) in CACNA2D2. We compared clinical and molecular findings from all nine individuals reported to date and note that cerebellar atrophy is shared among all.

INTERPRETATION

Our study supports the candidacy of CACNA2D2 as a disease gene associated with a phenotypic spectrum of neurological disease that include features of developmental and epileptic encephalopathy, ataxia, and cerebellar atrophy. Age at presentation may affect apparent penetrance of neurogenetic trait manifestations and of a particular clinical neurological endophenotype, for example, seizures or ataxia.

The Genomics of Arthrogryposis, a Complex Trait: Candidate Genes and Further Evidence for Oligogenic Inheritance

Arthrogryposis is a clinical finding that is present either as a feature of a neuromuscular condition or as part of a systemic disease in over 400 Mendelian conditions. The underlying molecular etiology remains largely unknown because of genetic and phenotypic heterogeneity. We applied exome sequencing (ES) in a cohort of 89 families with the clinical sign of arthrogryposis. Additional molecular techniques including array comparative genomic hybridization (aCGH) and Droplet Digital PCR (ddPCR) were performed on individuals who were found to have pathogenic copy number variants (CNVs) and mosaicism, respectively. A molecular diagnosis was established in 65.2% (58/89) of families. Eleven out of 58 families (19.0%) showed evidence for potential involvement of pathogenic variation at more than one locus, probably driven by absence of heterozygosity (AOH) burden due to identity-by-descent (IBD). RYR3, MYOM2, ERGIC1, SPTBN4, and ABCA7 represent genes, identified in two or more families, for which mutations are probably causative for arthrogryposis. We also provide evidence for the involvement of CNVs in the etiology of arthrogryposis and for the idea that both mono-allelic and bi-allelic variants in the same gene cause either similar or distinct syndromes. We were able to identify the molecular etiology in nine out of 20 families who underwent reanalysis. In summary, our data from family-based ES further delineate the molecular etiology of arthrogryposis, yielded several candidate disease-associated genes, and provide evidence for mutational burden in a biological pathway or network. Our study also highlights the importance of reanalysis of individuals with unsolved diagnoses in conjunction with sequencing extended family members.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.

Copy number variant and runs of homozygosity detection by microarrays enabled more precise molecular diagnoses in 11,020 clinical exome cases

Background

Exome sequencing (ES) has been successfully applied in clinical detection of single nucleotide variants (SNVs) and small indels. However, identification of copy number variants (CNVs) using ES data remains challenging. The purpose of this study is to understand the contribution of CNVs and copy neutral runs of homozygosity (ROH) in molecular diagnosis of patients referred for ES.

Methods

In a cohort of 11,020 consecutive ES patients, an Illumina SNP array analysis interrogating mostly coding SNPs was performed as a quality control (QC) measurement and for CNV/ROH detection. Among these patients, clinical chromosomal microarray analysis (CMA) was performed at Baylor Genetics (BG) on 3229 patients, either before, concurrently, or after ES. We retrospectively analyzed the findings from CMA and the QC array.

Results

The QC array can detect ~ 70% of pathogenic/likely pathogenic CNVs (PCNVs) detectable by CMA. Out of the 11,020 ES cases, the QC array identified PCNVs in 327 patients and uniparental disomy (UPD) disorder-related ROH in 10 patients. The overall PCNV/UPD detection rate was 5.9% in the 3229 ES patients who also had CMA at BG; PCNV/UPD detection rate was higher in concurrent ES and CMA than in ES with prior CMA (7.2% vs 4.6%). The PCNVs/UPD contributed to the molecular diagnoses in 17.4% (189/1089) of molecularly diagnosed ES cases with CMA and were estimated to contribute in 10.6% of all molecularly diagnosed ES cases. Dual diagnoses with both PCNVs and SNVs were detected in 38 patients. PCNVs affecting single recessive disorder genes in a compound heterozygous state with SNVs were detected in 4 patients, and homozygous deletions (mostly exonic deletions) were detected in 17 patients. A higher PCNV detection rate was observed for patients with syndromic phenotypes and/or cardiovascular abnormalities.

Conclusions

Our clinical genomics study demonstrates that detection of PCNV/UPD through the QC array or CMA increases ES diagnostic rate, provides more precise molecular diagnosis for dominant as well as recessive traits, and enables more complete genetic diagnoses in patients with dual or multiple molecular diagnoses. Concurrent ES and CMA using an array with exonic coverage for disease genes enables most effective detection of both CNVs and SNVs and therefore is recommended especially in time-sensitive clinical situations.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0639-5) contains supplementary material, which is available to authorized users.

Interchromosomal template-switching as a novel molecular mechanism for imprinting perturbations associated with Temple syndrome

Background

Intrachromosomal triplications (TRP) can contribute to disease etiology via gene dosage effects, gene disruption, position effects, or fusion gene formation. Recently, post-zygotic de novo triplications adjacent to copy-number neutral genomic intervals with runs of homozygosity (ROH) have been shown to result in uniparental isodisomy (UPD). The genomic structure of these complex genomic rearrangements (CGRs) shows a consistent pattern of an inverted triplication flanked by duplications (DUP-TRP/INV-DUP) formed by an iterative DNA replisome template-switching mechanism during replicative repair of a single-ended, double-stranded DNA (seDNA), the ROH results from an interhomolog or nonsister chromatid template switch. It has been postulated that these CGRs may lead to genetic abnormalities in carriers due to dosage-sensitive genes mapping within the copy-number variant regions, homozygosity for alleles at a locus causing an autosomal recessive (AR) disease trait within the ROH region, or imprinting-associated diseases.

Methods

Here, we report a family wherein the affected subject carries a de novo 2.2-Mb TRP followed by 42.2 Mb of ROH and manifests clinical features overlapping with those observed in association with chromosome 14 maternal UPD (UPD(14)mat). UPD(14)mat can cause clinical phenotypic features enabling a diagnosis of Temple syndrome. This CGR was then molecularly characterized by high-density custom aCGH, genome-wide single-nucleotide polymorphism (SNP) and methylation arrays, exome sequencing (ES), and the Oxford Nanopore long-read sequencing technology.

Results

We confirmed the postulated DUP-TRP/INV-DUP structure by multiple orthogonal genomic technologies in the proband. The methylation status of known differentially methylated regions (DMRs) on chromosome 14 revealed that the subject shows the typical methylation pattern of UPD(14)mat. Consistent with these molecular findings, the clinical features overlap with those observed in Temple syndrome, including speech delay.

Conclusions

These data provide experimental evidence that, in humans, triplication can lead to segmental UPD and imprinting disease. Importantly, genotype/phenotype analyses further reveal how a post-zygotically generated complex structural variant, resulting from a replication-based mutational mechanism, contributes to expanding the clinical phenotype of known genetic syndromes. Mechanistically, such events can distort transmission genetics resulting in homozygosity at a locus for which only one parent is a carrier as well as cause imprinting diseases.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0633-y) contains supplementary material, which is available to authorized users.

Insights into genetics, human biology and disease gleaned from family based genomic studies

Identifying genes and variants contributing to rare disease phenotypes and Mendelian conditions informs biology and medicine, yet potential phenotypic consequences for variation of >75% of the ~20,000 annotated genes in the human genome are lacking. Technical advances to assess rare variation genome-wide, particularly exome sequencing (ES), enabled establishment in the United States of the National Institutes of Health (NIH)-supported Centers for Mendelian Genomics (CMGs) and have facilitated collaborative studies resulting in novel "disease gene" discoveries. Pedigree-based genomic studies and rare variant analyses in families with suspected Mendelian conditions have led to the elucidation of hundreds of novel disease genes and highlighted the impact of de novo mutational events, somatic variation underlying nononcologic traits, incompletely penetrant alleles, phenotypes with high locus heterogeneity, and multilocus pathogenic variation. Herein, we highlight CMG collaborative discoveries that have contributed to understanding both rare and common diseases and discuss opportunities for future discovery in single-locus Mendelian disorder genomics. Phenotypic annotation of all human genes; development of bioinformatic tools and analytic methods; exploration of non-Mendelian modes of inheritance including reduced penetrance, multilocus variation, and oligogenic inheritance; construction of allelic series at a locus; enhanced data sharing worldwide; and integration with clinical genomics are explored. Realizing the full contribution of rare disease research to functional annotation of the human genome, and further illuminating human biology and health, will lay the foundation for the Precision Medicine Initiative.

A retrospective review of multiple findings in diagnostic exome sequencing: half are distinct and half are overlapping diagnoses

PURPOSE

We evaluated clinical and genetic features enriched in patients with multiple Mendelian conditions to determine which patients are more likely to have multiple potentially relevant genetic findings (MPRF).

METHODS

Results of the first 7698 patients who underwent exome sequencing at Ambry Genetics were reviewed. Clinical and genetic features were examined and degree of phenotypic overlap between the genetic diagnoses was evaluated.

RESULTS

Among patients referred for exome sequencing, 2% had MPRF. MPRF were more common in patients from consanguineous families and patients with greater clinical complexity. The difference in average number of organ systems affected is small: 4.3 (multiple findings) vs. 3.9 (single finding) and may not be distinguished in clinic.

CONCLUSION

Patients with multiple genetic diagnoses had a slightly higher number of organ systems affected than patients with single genetic diagnoses, largely because the comorbid conditions affected overlapping organ systems. Exome testing may be beneficial for all cases with multiple organ systems affected. The identification of multiple relevant genetic findings in 2% of exome patients highlights the utility of a comprehensive molecular workup and updated interpretation of existing genomic data; a single definitive molecular diagnosis from analysis of a limited number of genes may not be the end of a diagnostic odyssey.

TBX6-associated congenital scoliosis (TACS) as a clinically distinguishable subtype of congenital scoliosis: further evidence supporting the compound inheritance and TBX6 gene dosage model

PURPOSE

To characterize clinically measurable endophenotypes, implicating the TBX6 compound inheritance model.

METHODS

Patients with congenital scoliosis (CS) from China(N = 345, cohort 1), Japan (N = 142, cohort 2), and the United States (N = 10, cohort 3) were studied. Clinically measurable endophenotypes were compared according to the TBX6 genotypes. A mouse model for Tbx6 compound inheritance (N = 52) was investigated by micro computed tomography (micro-CT). A clinical diagnostic algorithm (TACScore) was developed to assist in clinical recognition of TBX6-associated CS (TACS).

RESULTS

In cohort 1, TACS patients (N = 33) were significantly younger at onset than the remaining CS patients (P = 0.02), presented with one or more hemivertebrae/butterfly vertebrae (P = 4.9 × 10^‒8), and exhibited vertebral malformations involving the lower part of the spine (T8-S5, P = 4.4 × 10^‒3); observations were confirmed in two replication cohorts. Simple rib anomalies were prevalent in TACS patients (P = 3.1 × 10^‒7), while intraspinal anomalies were uncommon (P = 7.0 × 10^‒7). A clinically usable TACScore was developed with an area under the curve (AUC) of 0.9 (P = 1.6 × 10^‒15). A Tbx6^{-/mh (mild-hypomorphic)} mouse model supported that a gene dosage effect underlies the TACS phenotype.

CONCLUSION

TACS is a clinically distinguishable entity with consistent clinically measurable endophenotypes. The type and distribution of vertebral column abnormalities in TBX6/Tbx6 compound inheritance implicate subtle perturbations in gene dosage as a cause of spine developmental birth defects responsible for about 10% of CS.

Limb girdle muscular dystrophy D3 HNRNPDL related in a Chinese family with distal muscle weakness caused by a mutation in the prion-like domain

Limb-girdle muscular dystrophies (LGMD) are a group of clinically and genetically heterogeneous diseases characterized by weakness and wasting of the pelvic and shoulder girdle muscles. Twenty-four recessive LGMD (types R1-R24) and five dominant LGMD (types D1-D5) have been identified with characterization of mutations in various genes. To date, LGMD D3 (previously known as LGMD1G) has been characterized in only two families with Brazilian or Uruguayan origin. Each was caused by a distinct mutation at codon 378 in the prion-like domain of HNRNPDL encoding heterogeneous nuclear ribonucleoprotein D like (HNRNPDL), an RNA processing protein. Our study characterized eight patients suffering from LGMD D3 in a Chinese family spanning three generations. Muscle biopsy specimens from two patients showed a myopathy with rimmed vacuoles. Sequencing analysis revealed a heterozygous c.1132G > A (p.D378N) mutation in HNRNPDL that co-segregated with disease phenotype in the family. The same mutation has been identified previously in the Brazilian family with LGMD D3. However, most patients in the current family showed distal as well as proximal limb weakness rather than weakness of toe and finger flexor muscles that were typical features in the other two LGMD D3 families reported previously. The present study indicates that the same mutation in HNRNPDL results in various phenotypes of LGMD D3. That all mutations in three unrelated families with different ethnic background occur at the same position in codon 378 of HNRNPDL gene suggests a mutation hotspot. Acceleration of intrinsic self-aggregation of HNRNPDL caused by mutation of the prior-like domain may contribute to the pathogenesis of the disease.

Clinical exome sequencing for fetuses with ultrasound abnormalities and a suspected Mendelian disorder

Background

Exome sequencing is now being incorporated into clinical care for pediatric and adult populations, but its integration into prenatal diagnosis has been more limited. One reason for this is the paucity of information about the clinical utility of exome sequencing in the prenatal setting.

Methods

We retrospectively reviewed indications, results, time to results (turnaround time, TAT), and impact of exome results for 146 consecutive “fetal exomes” performed in a clinical diagnostic laboratory between March 2012 and November 2017. We define a fetal exome as one performed on a sample obtained from a fetus or a product of conception with at least one structural anomaly detected by prenatal imaging or autopsy. Statistical comparisons were performed using Fisher’s exact test.

Results

Prenatal exome yielded an overall molecular diagnostic rate of 32% (n = 46/146). Of the 46 molecular diagnoses, 50% were autosomal dominant disorders (n = 23/46), 41% were autosomal recessive disorders (n = 19/46), and 9% were X-linked disorders (n = 4/46). The molecular diagnostic rate was highest for fetuses with anomalies affecting multiple organ systems and for fetuses with craniofacial anomalies. Out of 146 cases, a prenatal trio exome option designed for ongoing pregnancies was performed on 62 fetal specimens, resulting in a diagnostic yield of 35% with an average TAT of 14 days for initial reporting (excluding tissue culture time). The molecular diagnoses led to refined recurrence risk estimates, altered medical management, and informed reproductive planning for families.

Conclusion

Exome sequencing is a useful diagnostic tool when fetal structural anomalies suggest a genetic etiology, but other standard prenatal genetic tests did not provide a diagnosis.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0582-x) contains supplementary material, which is available to authorized users.

Molecular genetics of the POMT1-related muscular dystrophy-dystroglycanopathies

Protein O-mannosyltransferase 1 (POMT1) is a critical enzyme participating in the first step of protein O-mannosylation. Mutations in the coding gene, POMT1, have been described to be related to a series of autosomal recessive disorders associated with defective alpha-dystroglycan glycosylation, later termed muscular dystrophy-dystroglycanopathies (MDDGs). MDDGs are characterized by a broad phenotypic spectrum of congenital muscular dystrophy or later-onset limb-girdle muscular dystrophy, accompanied by variable degrees of intellectual disability, brain defects, and ocular abnormalities. To date, at least 76 disease-associated mutations in the POMT1 gene, including missense, nonsense, splicing, deletion, insertion/duplication, and insertion-deletion mutations, have been reported in the literature. In this review, we highlight the present knowledge of the identified disease-associated POMT1 gene mutations and genetic animal models related to the POMT1 gene. This review may help further normative classification of phenotypes, assist in definite clinical and genetic diagnoses, and genetic counseling, and may comprehensively improve our understanding of the basis of complex phenotypes and possible pathogenic mechanisms involved.

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants

Purpose

To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants.

Methods

We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants.

Results

The number of rare likely deleterious variants in functionally intolerant genes (“other hits”) correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes.

Conclusion

Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

Identification of likely pathogenic and known variants in TSPEAR, LAMB3, BCOR, and WNT10A in four Turkish families with tooth agenesis

Tooth agenesis (TA), the failure of development of one or more permanent teeth, is a common craniofacial abnormality observed in different world populations. The genetic etiology of TA is heterogeneous; more than a dozen genes have been associated with isolated or nonsyndromic TA, and more than 80 genes with syndromic forms. In this study, we applied whole exome sequencing (WES) to identify candidate genes contributing to TA in four Turkish families. Likely pathogenic variants with a low allele frequency in the general population were identified in four disease-associated genes, including two distinct variants in TSPEAR, associated with syndromic and isolated TA in one family each; a variant in LAMB3 associated with syndromic TA in one family; and a variant in BCOR plus a disease-associated WNT10A variant in one family with syndromic TA. With the notable exception of WNT10A (Tooth agenesis, selective, 4, MIM #150400), the genotype-phenotype relationships described in the present cohort represent an expansion of the clinical spectrum associated with these genes: TSPEAR (Deafness, autosomal recessive 98, MIM #614861), LAMB3 (Amelogenesis imperfecta, type IA, MIM #104530; Epidermolysis bullosa, junctional, MIMs #226700 and #226650), and BCOR (Microphthalmia, syndromic 2, MIM #300166). We provide evidence supporting the candidacy of these genes with TA, and propose TSPEAR as a novel nonsyndromic TA gene. Our data also suggest potential multilocus genomic variation, or mutational burden, in a single family, involving the BCOR and WNT10A loci, underscoring the complexity of the genotype-phenotype relationship in the common complex trait of TA.