Exome copy number variant detection, analysis, and classification in a large cohort of families with undiagnosed rare genetic disease.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs.

Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Erythromelalgia caused by the missense mutation p.Arg220Pro in an alternatively spliced exon of SCN9A (NaV1.7).

Erythromelalgia (EM), is a familial pain syndrome characterized by episodic 'burning' pain, warmth, and erythema. EM is caused by monoallelic variants in SCN9A, which encodes the voltage-gated sodium channel (NaV) NaV1.7. Over 25 different SCN9A mutations attributed to EM have been described to date, all identified in the SCN9A transcript utilizing exon 6N. Here we report a novel SCN9A missense variant identified in seven related individuals with stereotypic episodes of bilateral lower limb pain presenting in childhood. The variant, XM_011511617.3:c.659G>C;p.(Arg220Pro), resides in the exon 6A of SCN9A, an exon previously shown to be selectively incorporated by developmentally regulated alternative splicing. The mutation is located in the voltage-sensing S4 segment of domain I, which is important for regulating channel activation. Functional analysis showed the p.Arg220Pro mutation altered voltage-dependent activation and delayed channel inactivation, consistent with a NaV1.7 gain-of-function molecular phenotype. These results demonstrate that alternatively spliced isoforms of SCN9A should be included in all genomic testing of EM.

Bi-allelic variants in HMGCR cause an autosomal-recessive progressive limb-girdle muscular dystrophy.

Statins are a mainstay intervention for cardiovascular disease prevention, yet their use can cause rare severe myopathy. HMG-CoA reductase, an essential enzyme in the mevalonate pathway, is the target of statins. We identified nine individuals from five unrelated families with unexplained limb-girdle like muscular dystrophy and bi-allelic variants in HMGCR via clinical and research exome sequencing. The clinical features resembled other genetic causes of muscular dystrophy with incidental high CPK levels (>1,000 U/L), proximal muscle weakness, variable age of onset, and progression leading to impaired ambulation. Muscle biopsies in most affected individuals showed non-specific dystrophic changes with non-diagnostic immunohistochemistry. Molecular modeling analyses revealed variants to be destabilizing and affecting protein oligomerization. Protein activity studies using three variants (p.Asp623Asn, p.Tyr792Cys, and p.Arg443Gln) identified in affected individuals confirmed decreased enzymatic activity and reduced protein stability. In summary, we showed that individuals with bi-allelic amorphic (i.e., null and/or hypomorphic) variants in HMGCR display phenotypes that resemble non-genetic causes of myopathy involving this reductase. This study expands our knowledge regarding the mechanisms leading to muscular dystrophy through dysregulation of the mevalonate pathway, autoimmune myopathy, and statin-induced myopathy.

Bi-allelic variants in INTS11 are associated with a complex neurological disorder.

The Integrator complex is a multi-subunit protein complex that regulates the processing of nascent RNAs transcribed by RNA polymerase II (RNAPII), including small nuclear RNAs, enhancer RNAs, telomeric RNAs, viral RNAs, and protein-coding mRNAs. Integrator subunit 11 (INTS11) is the catalytic subunit that cleaves nascent RNAs, but, to date, mutations in this subunit have not been linked to human disease. Here, we describe 15 individuals from 10 unrelated families with bi-allelic variants in INTS11 who present with global developmental and language delay, intellectual disability, impaired motor development, and brain atrophy. Consistent with human observations, we find that the fly ortholog of INTS11, dIntS11, is essential and expressed in the central nervous systems in a subset of neurons and most glia in larval and adult stages. Using Drosophila as a model, we investigated the effect of seven variants. We found that two (p.Arg17Leu and p.His414Tyr) fail to rescue the lethality of null mutants, indicating that they are strong loss-of-function variants. Furthermore, we found that five variants (p.Gly55Ser, p.Leu138Phe, p.Lys396Glu, p.Val517Met, and p.Ile553Glu) rescue lethality but cause a shortened lifespan and bang sensitivity and affect locomotor activity, indicating that they are partial loss-of-function variants. Altogether, our results provide compelling evidence that integrity of the Integrator RNA endonuclease is critical for brain development.

SRSF1 haploinsufficiency is responsible for a syndromic developmental disorder associated with intellectual disability.

SRSF1 (also known as ASF/SF2) is a non-small nuclear ribonucleoprotein (non-snRNP) that belongs to the arginine/serine (R/S) domain family. It recognizes and binds to mRNA, regulating both constitutive and alternative splicing. The complete loss of this proto-oncogene in mice is embryonically lethal. Through international data sharing, we identified 17 individuals (10 females and 7 males) with a neurodevelopmental disorder (NDD) with heterozygous germline SRSF1 variants, mostly de novo, including three frameshift variants, three nonsense variants, seven missense variants, and two microdeletions within region 17q22 encompassing SRSF1. Only in one family, the de novo origin could not be established. All individuals featured a recurrent phenotype including developmental delay and intellectual disability (DD/ID), hypotonia, neurobehavioral problems, with variable skeletal (66.7%) and cardiac (46%) anomalies. To investigate the functional consequences of SRSF1 variants, we performed in silico structural modeling, developed an in vivo splicing assay in Drosophila, and carried out episignature analysis in blood-derived DNA from affected individuals. We found that all loss-of-function and 5 out of 7 missense variants were pathogenic, leading to a loss of SRSF1 splicing activity in Drosophila, correlating with a detectable and specific DNA methylation episignature. In addition, our orthogonal in silico, in vivo, and epigenetics analyses enabled the separation of clearly pathogenic missense variants from those with uncertain significance. Overall, these results indicated that haploinsufficiency of SRSF1 is responsible for a syndromic NDD with ID due to a partial loss of SRSF1-mediated splicing activity.

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder.

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Pathogenic variants in nucleoporin TPR (translocated promoter region, nuclear basket protein) cause severe intellectual disability in humans.

The nuclear pore complex (NPC) is a multi-protein complex that regulates the trafficking of macromolecules between the nucleus and cytoplasm. Genetic variants in components of the NPC have been shown to cause a range of neurological disorders, including intellectual disability and microcephaly. Translocated promoter region, nuclear basket protein (TPR) is a critical scaffolding element of the nuclear facing interior of the NPC. Here, we present two siblings with biallelic variants in TPR who present with a phenotype of microcephaly, ataxia and severe intellectual disability. The variants result in a premature truncation variant, and a splice variant leading to a 12-amino acid deletion respectively. Functional analyses in patient fibroblasts demonstrate significantly reduced TPR levels, and decreased TPR-containing NPC density. A compensatory increase in total NPC levels was observed, and decreased global RNA intensity in the nucleus. The discovery of variants that partly disable TPR function provide valuable insight into this essential protein in human disease, and our findings suggest that TPR variants are the cause of the siblings' neurological disorder.

ABHD16A deficiency causes a complicated form of hereditary spastic paraplegia associated with intellectual disability and cerebral anomalies.

ABHD16A (abhydrolase domain-containing protein 16A, phospholipase) encodes the major phosphatidylserine (PS) lipase in the brain. PS lipase synthesizes lysophosphatidylserine, an important signaling lipid that functions in the mammalian central nervous system. ABHD16A has not yet been associated with a human disease. In this report, we present a cohort of 11 affected individuals from six unrelated families with a complicated form of hereditary spastic paraplegia (HSP) who carry bi-allelic deleterious variants in ABHD16A. Affected individuals present with a similar phenotype consisting of global developmental delay/intellectual disability, progressive spasticity affecting the upper and lower limbs, and corpus callosum and white matter anomalies. Immunoblot analysis on extracts from fibroblasts from four affected individuals demonstrated little to no ABHD16A protein levels compared to controls. Our findings add ABHD16A to the growing list of lipid genes in which dysregulation can cause complicated forms of HSP and begin to describe the molecular etiology of this condition.

A form of muscular dystrophy associated with pathogenic variants in JAG2.

JAG2 encodes the Notch ligand Jagged2. The conserved Notch signaling pathway contributes to the development and homeostasis of multiple tissues, including skeletal muscle. We studied an international cohort of 23 individuals with genetically unsolved muscular dystrophy from 13 unrelated families. Whole-exome sequencing identified rare homozygous or compound heterozygous JAG2 variants in all 13 families. The identified bi-allelic variants include 10 missense variants that disrupt highly conserved amino acids, a nonsense variant, two frameshift variants, an in-frame deletion, and a microdeletion encompassing JAG2. Onset of muscle weakness occurred from infancy to young adulthood. Serum creatine kinase (CK) levels were normal or mildly elevated. Muscle histology was primarily dystrophic. MRI of the lower extremities revealed a distinct, slightly asymmetric pattern of muscle involvement with cores of preserved and affected muscles in quadriceps and tibialis anterior, in some cases resembling patterns seen in POGLUT1-associated muscular dystrophy. Transcriptome analysis of muscle tissue from two participants suggested misregulation of genes involved in myogenesis, including PAX7. In complementary studies, Jag2 downregulation in murine myoblasts led to downregulation of multiple components of the Notch pathway, including Megf10. Investigations in Drosophila suggested an interaction between Serrate and Drpr, the fly orthologs of JAG1/JAG2 and MEGF10, respectively. In silico analysis predicted that many Jagged2 missense variants are associated with structural changes and protein misfolding. In summary, we describe a muscular dystrophy associated with pathogenic variants in JAG2 and evidence suggests a disease mechanism related to Notch pathway dysfunction.

Bi-allelic variants in IPO8 cause a connective tissue disorder associated with cardiovascular defects, skeletal abnormalities, and immune dysregulation.

Dysregulated transforming growth factor TGF-ß signaling underlies the pathogenesis of genetic disorders affecting the connective tissue such as Loeys-Dietz syndrome. Here, we report 12 individuals with bi-allelic loss-of-function variants in IPO8 who presented with a syndromic association characterized by cardio-vascular anomalies, joint hyperlaxity, and various degree of dysmorphic features and developmental delay as well as immune dysregulation; the individuals were from nine unrelated families. Importin 8 belongs to the karyopherin family of nuclear transport receptors and was previously shown to mediate TGF-ß-dependent SMADs trafficking to the nucleus in vitro. The important in vivo role of IPO8 in pSMAD nuclear translocation was demonstrated by CRISPR/Cas9-mediated inactivation in zebrafish. Consistent with IPO8's role in BMP/TGF-ß signaling, ipo8-/- zebrafish presented mild to severe dorso-ventral patterning defects during early embryonic development. Moreover, ipo8-/- zebrafish displayed severe cardiovascular and skeletal defects that mirrored the human phenotype. Our work thus provides evidence that IPO8 plays a critical and non-redundant role in TGF-ß signaling during development and reinforces the existing link between TGF-ß signaling and connective tissue defects.

Bi-allelic ADARB1 Variants Associated with Microcephaly, Intellectual Disability, and Seizures.

The RNA editing enzyme ADAR2 is essential for the recoding of brain transcripts. Impaired ADAR2 editing leads to early-onset epilepsy and premature death in a mouse model. Here, we report bi-allelic variants in ADARB1, the gene encoding ADAR2, in four unrelated individuals with microcephaly, intellectual disability, and epilepsy. In one individual, a homozygous variant in one of the double-stranded RNA-binding domains (dsRBDs) was identified. In the others, variants were situated in or around the deaminase domain. To evaluate the effects of these variants on ADAR2 enzymatic activity, we performed in vitro assays with recombinant proteins in HEK293T cells and ex vivo assays with fibroblasts derived from one of the individuals. We demonstrate that these ADAR2 variants lead to reduced editing activity on a known ADAR2 substrate. We also demonstrate that one variant leads to changes in splicing of ADARB1 transcript isoforms. These findings reinforce the importance of RNA editing in brain development and introduce ADARB1 as a genetic etiology in individuals with intellectual disability, microcephaly, and epilepsy.

Bi-allelic Variations of SMO in Humans Cause a Broad Spectrum of Developmental Anomalies Due to Abnormal Hedgehog Signaling.

The evolutionarily conserved hedgehog (Hh) pathway is essential for organogenesis and plays critical roles in postnatal tissue maintenance and renewal. A unique feature of the vertebrate Hh pathway is that signal transduction requires the primary cilium (PC) where major pathway components are dynamically enriched. These factors include smoothened (SMO) and patched, which constitute the core reception system for sonic hedgehog (SHH) as well as GLI transcription factors, the key mediators of the pathway. Here, we report bi-allelic loss-of-function variations in SMO in seven individuals from five independent families; these variations cause a wide phenotypic spectrum of developmental anomalies affecting the brain (hypothalamic hamartoma and microcephaly), heart (atrioventricular septal defect), skeleton (postaxial polydactyly, narrow chest, and shortening of long bones), and enteric nervous system (aganglionosis). Cells derived from affected individuals showed normal ciliogenesis but severely altered Hh-signal transduction as a result of either altered PC trafficking or abnormal activation of the pathway downstream of SMO. In addition, Hh-independent GLI2 accumulation at the PC tip in cells from the affected individuals suggests a potential function of SMO in regulating basal ciliary trafficking of GLI2 when the pathway is off. Thus, loss of SMO function results in abnormal PC dynamics of key components of the Hh signaling pathway and leads to a large continuum of malformations in humans.

A mutation affecting laminin alpha 5 polymerisation gives rise to a syndromic developmental disorder.

Laminin alpha 5 (LAMA5) is a member of a large family of proteins that trimerise and then polymerise to form a central component of all basement membranes. Consequently, the protein plays an instrumental role in shaping the normal development of the kidney, skin, neural tube, lung and limb, and many other organs and tissues. Pathogenic mutations in some laminins have been shown to cause a range of largely syndromic conditions affecting the competency of the basement membranes to which they contribute. We report the identification of a mutation in the polymerisation domain of LAMA5 in a patient with a complex syndromic disease characterised by defects in kidney, craniofacial and limb development, and by a range of other congenital defects. Using CRISPR-generated mouse models and biochemical assays, we demonstrate the pathogenicity of this variant, showing that the change results in a failure of the polymerisation of α/ß/γ laminin trimers. Comparing these in vivo phenotypes with those apparent upon gene deletion in mice provides insights into the specific functional importance of laminin polymerisation during development and tissue homeostasis.

De novo PHF5A variants are associated with craniofacial abnormalities, developmental delay, and hypospadias.

PURPOSE: The SF3B splicing complex is composed of SF3B1-6 and PHF5A. We report a developmental disorder caused by de novo variants in PHF5A. METHODS: Clinical, genomic, and functional studies using subject-derived fibroblasts and a heterologous cellular system were performed. RESULTS: We studied 9 subjects with congenital malformations, including preauricular tags and hypospadias, growth abnormalities, and developmental delay who had de novo heterozygous PHF5A variants, including 4 loss-of-function (LOF), 3 missense, 1 splice, and 1 start-loss variant. In subject-derived fibroblasts with PHF5A LOF variants, wild-type and variant PHF5A mRNAs had a 1:1 ratio, and PHF5A mRNA levels were normal. Transcriptome sequencing revealed alternative promoter use and downregulated genes involved in cell-cycle regulation. Subject and control fibroblasts had similar amounts of PHF5A with the predicted wild-type molecular weight and of SF3B1-3 and SF3B6. SF3B complex formation was unaffected in 2 subject cell lines. CONCLUSION: Our data suggest the existence of feedback mechanisms in fibroblasts with PHF5A LOF variants to maintain normal levels of SF3B components. These compensatory mechanisms in subject fibroblasts with PHF5A or SF3B4 LOF variants suggest disturbed autoregulation of mutated splicing factor genes in specific cell types, that is, neural crest cells, during embryonic development rather than haploinsufficiency as pathomechanism.

Variants in Mitochondrial ATP Synthase Cause Variable Neurologic Phenotypes.

OBJECTIVE: ATP synthase (ATPase) is responsible for the majority of ATP production. Nevertheless, disease phenotypes associated with mutations in ATPase subunits are extremely rare. We aimed at expanding the spectrum of ATPase-related diseases. METHODS: Whole-exome sequencing in cohorts with 2,962 patients diagnosed with mitochondrial disease and/or dystonia and international collaboration were used to identify deleterious variants in ATPase-encoding genes. Findings were complemented by transcriptional and proteomic profiling of patient fibroblasts. ATPase integrity and activity were assayed using cells and tissues from 5 patients. RESULTS: We present 10 total individuals with biallelic or de novo monoallelic variants in nuclear ATPase subunit genes. Three unrelated patients showed the same homozygous missense ATP5F1E mutation (including one published case). An intronic splice-disrupting alteration in compound heterozygosity with a nonsense variant in ATP5PO was found in one patient. Three patients had de novo heterozygous missense variants in ATP5F1A, whereas another 3 were heterozygous for ATP5MC3 de novo missense changes. Bioinformatics methods and populational data supported the variants' pathogenicity. Immunohistochemistry, proteomics, and/or immunoblotting revealed significantly reduced ATPase amounts in association to ATP5F1E and ATP5PO mutations. Diminished activity and/or defective assembly of ATPase was demonstrated by enzymatic assays and/or immunoblotting in patient samples bearing ATP5F1A-p.Arg207His, ATP5MC3-p.Gly79Val, and ATP5MC3-p.Asn106Lys. The associated clinical profiles were heterogeneous, ranging from hypotonia with spontaneous resolution (1/10) to epilepsy with early death (1/10) or variable persistent abnormalities, including movement disorders, developmental delay, intellectual disability, hyperlactatemia, and other neurologic and systemic features. Although potentially reflecting an ascertainment bias, dystonia was common (7/10). INTERPRETATION: Our results establish evidence for a previously unrecognized role of ATPase nuclear-gene defects in phenotypes characterized by neurodevelopmental and neurodegenerative features. ANN NEUROL 2022;91:225-237.

Variants in DTNA cause a mild, dominantly inherited muscular dystrophy.

DTNA encodes α-dystrobrevin, a component of the macromolecular dystrophin-glycoprotein complex (DGC) that binds to dystrophin/utrophin and α-syntrophin. Mice lacking α-dystrobrevin have a muscular dystrophy phenotype, but variants in DTNA have not previously been associated with human skeletal muscle disease. We present 12 individuals from four unrelated families with two different monoallelic DTNA variants affecting the coiled-coil domain of α-dystrobrevin. The five affected individuals from family A harbor a c.1585G > A; p.Glu529Lys variant, while the recurrent c.1567_1587del; p.Gln523_Glu529del DTNA variant was identified in the other three families (family B: four affected individuals, family C: one affected individual, and family D: two affected individuals). Myalgia and exercise intolerance, with variable ages of onset, were reported in 10 of 12 affected individuals. Proximal lower limb weakness with onset in the first decade of life was noted in three individuals. Persistent elevations of serum creatine kinase (CK) levels were detected in 11 of 12 affected individuals, 1 of whom had an episode of rhabdomyolysis at 20 years of age. Autism spectrum disorder or learning disabilities were reported in four individuals with the c.1567_1587 deletion. Muscle biopsies in eight affected individuals showed mixed myopathic and dystrophic findings, characterized by fiber size variability, internalized nuclei, and slightly increased extracellular connective tissue and inflammation. Immunofluorescence analysis of biopsies from five affected individuals showed reduced α-dystrobrevin immunoreactivity and variably reduced immunoreactivity of other DGC proteins: dystrophin, α, ß, δ and γ-sarcoglycans, and α and ß-dystroglycans. The DTNA deletion disrupted an interaction between α-dystrobrevin and syntrophin. Specific variants in the coiled-coil domain of DTNA cause skeletal muscle disease with variable penetrance. Affected individuals show a spectrum of clinical manifestations, with severity ranging from hyperCKemia, myalgias, and exercise intolerance to childhood-onset proximal muscle weakness. Our findings expand the molecular etiologies of both muscular dystrophy and paucisymptomatic hyperCKemia, to now include monoallelic DTNA variants as a novel cause of skeletal muscle disease in humans.

A cryptic pathogenic NDUFV1 variant identified by RNA-seq in a patient with normal complex I activity in muscle and transient magnetic resonance imaging changes.

Mitochondrial respiratory chain disorders (MRC) are amongst the most common group of inborn errors of metabolism. MRC, of which complex I deficiency accounts for approximately a quarter, are very diverse, causing a wide range of clinical problems and can be difficult to diagnose. We report an illustrative MRC case whose diagnosis was elusive. Clinical signs included failure to thrive caused by recurrent vomiting, hypotonia and progressive loss of motor milestones. Initial brain imaging suggested Leigh syndrome but without expected diffusion restriction. Muscle respiratory chain enzymology was unremarkable. Whole-genome sequencing identified a maternally inherited NDUFV1 missense variant [NM_007103.4 (NDUFV1):c.1157G > A; p.(Arg386His)] and a paternally inherited synonymous variant [NM_007103.4 (NDUFV1):c.1080G > A; (p.Ser360=)]. RNA sequencing demonstrated aberrant splicing. This case emphasizes the diagnostic odyssey of a patient in whom a confirmed diagnosis was elusive because of atypical features and normal muscle respiratory chain enzyme (RCE) activities, along with a synonymous variant, which are often filtered out from genomic analyses. It also illustrates the following points: (1) complete resolution of magnetic resonance imaging changes may be part of the picture in mitochondrial disease; (2) analysis for synonymous variants is important for undiagnosed patients; and (3) RNA-seq is a powerful tool to demonstrate pathogenicity of putative splicing variants.

Lessons learnt from multifaceted diagnostic approaches to the first 150 families in Victoria's Undiagnosed Diseases Program.

BACKGROUND: Clinical exome sequencing typically achieves diagnostic yields of 30%-57.5% in individuals with monogenic rare diseases. Undiagnosed diseases programmes implement strategies to improve diagnostic outcomes for these individuals. AIM: We share the lessons learnt from the first 3 years of the Undiagnosed Diseases Program-Victoria, an Australian programme embedded within a clinical genetics service in the state of Victoria with a focus on paediatric rare diseases. METHODS: We enrolled families who remained without a diagnosis after clinical genomic (panel, exome or genome) sequencing between 2016 and 2018. We used family-based exome sequencing (family ES), family-based genome sequencing (family GS), RNA sequencing (RNA-seq) and high-resolution chromosomal microarray (CMA) with research-based analysis. RESULTS: In 150 families, we achieved a diagnosis or strong candidate in 64 (42.7%) (37 in known genes with a consistent phenotype, 3 in known genes with a novel phenotype and 24 in novel disease genes). Fifty-four diagnoses or strong candidates were made by family ES, six by family GS with RNA-seq, two by high-resolution CMA and two by data reanalysis. CONCLUSION: We share our lessons learnt from the programme. Flexible implementation of multiple strategies allowed for scalability and response to the availability of new technologies. Broad implementation of family ES with research-based analysis showed promising yields post a negative clinical singleton ES. RNA-seq offered multiple benefits in family ES-negative populations. International data sharing strategies were critical in facilitating collaborations to establish novel disease-gene associations. Finally, the integrated approach of a multiskilled, multidisciplinary team was fundamental to having diverse perspectives and strategic decision-making.