Systematic analysis of variants escaping nonsense-mediated decay uncovers candidate Mendelian diseases.

Protein-truncating variants (PTVs) near the 3' end of genes may escape nonsense-mediated decay (NMD). PTVs in the NMD-escape region (PTVescs) can cause Mendelian disease but are difficult to interpret given their varying impact on protein function. Previously, PTVesc burden was assessed in an epilepsy cohort, but no large-scale analysis has systematically evaluated these variants in rare disease. We performed a retrospective analysis of 29,031 neurodevelopmental disorder (NDD) parent-offspring trios referred for clinical exome sequencing to identify PTVesc de novo mutations (DNMs). We identified 1,376 PTVesc DNMs and 133 genes that were significantly enriched (binomial p < 0.001). The PTVesc-enriched genes included those with PTVescs previously described to cause dominant Mendelian disease (e.g., SEMA6B, PPM1D, and DAGLA). We annotated ClinVar variants for PTVescs and identified 948 genes with at least one high-confidence pathogenic variant. Twenty-two known Mendelian PTVesc-enriched genes had no prior evidence of PTVesc-associated disease. We found 22 additional PTVesc-enriched genes that are not well established to be associated with Mendelian disease, several of which showed phenotypic similarity between individuals harboring PTVesc variants in the same gene. Four individuals with PTVesc mutations in RAB1A had similar phenotypes including NDD and spasticity. PTVesc mutations in IRF2BP1 were found in two individuals who each had severe immunodeficiency manifesting in NDD. Three individuals with PTVesc mutations in LDB1 all had NDD and multiple congenital anomalies. Using a large-scale, systematic analysis of DNMs, we extend the mutation spectrum for known Mendelian disease-associated genes and identify potentially novel disease-associated genes.

Genetic modifiers of body mass index in individuals with cystic fibrosis.

To identify modifier loci underlying variation in body mass index (BMI) in persons with cystic fibrosis (pwCF), we performed a genome-wide association study (GWAS). Utilizing longitudinal height and weight data, along with demographic information and covariates from 4,393 pwCF, we calculated AvgBMIz representing the average of per-quarter BMI Z scores. The GWAS incorporated 9.8M single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) > 0.005 extracted from whole-genome sequencing (WGS) of each study subject. We observed genome-wide significant association with a variant in FTO (FaT mass and Obesity-associated gene; rs28567725; p value = 1.21e-08; MAF = 0.41, ß = 0.106; n = 4,393 individuals) and a variant within ADAMTS5 (A Disintegrin And Metalloproteinase with ThromboSpondin motifs 5; rs162500; p value = 2.11e-10; MAF = 0.005, ß = -0.768; n = 4,085 pancreatic-insufficient individuals). Notably, BMI-associated variants in ADAMTS5 occur on a haplotype that is much more common in African (AFR, MAF = 0.183) than European (EUR, MAF = 0.006) populations (1000 Genomes project). A polygenic risk score (PRS) calculated using 924 SNPs (excluding 17 in FTO) showed significant association with AvgBMIz (p value = 2.2e-16; r2 = 0.03). Association between variants in FTO and the PRS correlation reveals similarities in the genetic architecture of BMI in CF and the general population. Inclusion of Black individuals in whom the single-gene disorder CF is much less common but genomic diversity is greater facilitated detection of association with variants that are in LD with functional SNPs in ADAMTS5. Our results illustrate the importance of population diversity, particularly when attempting to identify variants that manifest only under certain physiologic conditions.

Contribution and therapeutic implications of retroelement insertions in ataxia telangiectasia.

Certain classes of genetic variation still escape detection in clinical sequencing analysis. One such class is retroelement insertion, which has been reported as a cause of Mendelian diseases and may offer unique therapeutic implications. Here, we conducted retroelement profiling on whole-genome sequencing data from a cohort of 237 individuals with ataxia telangiectasia (A-T). We found 15 individuals carrying retroelement insertions in ATM, all but one of which integrated in noncoding regions. Systematic functional characterization via RNA sequencing, RT-PCR, and/or minigene splicing assays showed that 12 out of 14 intronic insertions led or contributed to ATM loss of function by exon skipping or activating cryptic splice sites. We also present proof-of-concept antisense oligonucleotides that suppress cryptic exonization caused by a deep intronic retroelement insertion. These results provide an initial systematic estimate of the contribution of retroelements to the genetic architecture of recessive Mendelian disorders as ∼2.1%-5.5%. Our study highlights the importance of retroelement insertions as causal variants and therapeutic targets in genetic diseases.

GABBR1 monoallelic de novo variants linked to neurodevelopmental delay and epilepsy.

GABAB receptors are obligatory heterodimers responsible for prolonged neuronal inhibition in the central nervous system. The two receptor subunits are encoded by GABBR1 and GABBR2. Variants in GABBR2 have been associated with a Rett-like phenotype (MIM: 617903), epileptic encephalopathy (MIM: 617904), and milder forms of developmental delay with absence epilepsy. To date, however, no phenotypes associated with pathogenic variants of GABBR1 have been established. Through GeneMatcher, we have ascertained four individuals who each have a monoallelic GABBR1 de novo non-synonymous variant; these individuals exhibit motor and/or language delay, ranging from mild to severe, and in one case, epilepsy. Further phenotypic features include varying degrees of intellectual disability, learning difficulties, autism, ADHD, ODD, sleep disorders, and muscular hypotonia. We functionally characterized the four de novo GABBR1 variants, p.Glu368Asp, p.Ala397Val, p.Ala535Thr, and p.Gly673Asp, in transfected HEK293 cells. GABA fails to efficiently activate the variant receptors, most likely leading to an increase in the excitation/inhibition balance in the central nervous system. Variant p.Gly673Asp in transmembrane domain 3 (TMD3) renders the receptor completely inactive, consistent with failure of the receptor to reach the cell surface. p.Glu368Asp is located near the orthosteric binding site and reduces GABA potency and efficacy at the receptor. GABA exhibits normal potency but decreased efficacy at the p.Ala397Val and p.Ala535Thr variants. Functional characterization of GABBR1-related variants provides a rationale for understanding the severity of disease phenotypes and points to possible therapeutic strategies.

A spectrum of recessiveness among Mendelian disease variants in UK Biobank.

Recent work has found increasing evidence of mitigated, incompletely penetrant phenotypes in heterozygous carriers of recessive Mendelian disease variants. We leveraged whole-exome imputation within the full UK Biobank cohort (n ∼ 500K) to extend such analyses to 3,475 rare variants curated from ClinVar and OMIM. Testing these variants for association with 58 quantitative traits yielded 102 significant associations involving variants previously implicated in 34 different diseases. Notable examples included a POR missense variant implicated in Antley-Bixler syndrome that associated with a 1.76 (SE 0.27) cm increase in height and an ABCA3 missense variant implicated in interstitial lung disease that associated with reduced FEV1/FVC ratio. Association analyses with 1,134 disease traits yielded five additional variant-disease associations. We also observed contrasting levels of recessiveness between two more-common, classical Mendelian diseases. Carriers of cystic fibrosis variants exhibited increased risk of several mitigated disease phenotypes, whereas carriers of spinal muscular atrophy alleles showed no evidence of altered phenotypes. Incomplete penetrance of cystic fibrosis carrier phenotypes did not appear to be mediated by common allelic variation on the functional haplotype. Our results show that many disease-associated recessive variants can produce mitigated phenotypes in heterozygous carriers and motivate further work exploring penetrance mechanisms.

MRSD: A quantitative approach for assessing suitability of RNA-seq in the investigation of mis-splicing in Mendelian disease.

Variable levels of gene expression between tissues complicates the use of RNA sequencing of patient biosamples to delineate the impact of genomic variants. Here, we describe a gene- and tissue-specific metric to inform the feasibility of RNA sequencing. This overcomes limitations of using expression values alone as a metric to predict RNA-sequencing utility. We have derived a metric, minimum required sequencing depth (MRSD), that estimates the depth of sequencing required from RNA sequencing to achieve user-specified sequencing coverage of a gene, transcript, or group of genes. We applied MRSD across four human biosamples: whole blood, lymphoblastoid cell lines (LCLs), skeletal muscle, and cultured fibroblasts. MRSD has high precision (90.1%-98.2%) and overcomes transcript region-specific sequencing biases. Applying MRSD scoring to established disease gene panels shows that fibroblasts, of these four biosamples, are the optimum source of RNA for 63.1% of gene panels. Using this approach, up to 67.8% of the variants of uncertain significance in ClinVar that are predicted to impact splicing could be assayed by RNA sequencing in at least one of the biosamples. We demonstrate the utility and benefits of MRSD as a metric to inform functional assessment of splicing aberrations, in particular in the context of Mendelian genetic disorders to improve diagnostic yield.

An autosomal recessive variant in PYGM causes myophosphorylase deficiency in Red Angus composite cattle.

BACKGROUND: Between 2020 and 2022, eight calves in a Nebraska herd (composite Simmental, Red Angus, Gelbvieh) displayed exercise intolerance during forced activity. In some cases, the calves collapsed and did not recover. Available sire pedigrees contained a paternal ancestor within 2-4 generations in all affected calves. Pedigrees of the calves' dams were unavailable, however, the cows were ranch-raised and retained from prior breeding seasons, where bulls used for breeding occasionally had a common ancestor. Therefore, it was hypothesized that a de novo autosomal recessive variant was causative of exercise intolerance in these calves. RESULTS: A genome-wide association analysis utilizing SNP data from 6 affected calves and 715 herd mates, followed by whole-genome sequencing of 2 affected calves led to the identification of a variant in the gene PYGM (BTA29:g.42989581G > A). The variant, confirmed to be present in the skeletal muscle transcriptome, was predicted to produce a premature stop codon (p.Arg650*). The protein product of PYGM, myophosphorylase, breaks down glycogen in skeletal muscle. Glycogen concentrations were fluorometrically assayed as glucose residues demonstrating significantly elevated glycogen concentrations in affected calves compared to cattle carrying the variant and to wild-type controls. The absence of the PYGM protein product in skeletal muscle was confirmed by immunohistochemistry and label-free quantitative proteomics analysis; muscle degeneration was confirmed in biopsy and necropsy samples. Elevated skeletal muscle glycogen persisted after harvest, resulting in a high pH and dark-cutting beef, which is negatively perceived by consumers and results in an economic loss to the industry. Carriers of the variant did not exhibit differences in meat quality or any measures of animal well-being. CONCLUSIONS: Myophosphorylase deficiency poses welfare concerns for affected animals and negatively impacts the final product. The association of the recessive genotype with dark-cutting beef further demonstrates the importance of genetics to not only animal health but to the quality of their product. Although cattle heterozygous for the variant may not immediately affect the beef industry, identifying carriers will enable selection and breeding strategies to prevent the production of affected calves.

Mendelian pathway analysis of laboratory traits reveals distinct roles for ciliary subcompartments in common disease pathogenesis.

Rare monogenic disorders of the primary cilium, termed ciliopathies, are characterized by extreme presentations of otherwise common diseases, such as diabetes, hepatic fibrosis, and kidney failure. However, despite a recent revolution in our understanding of the cilium's role in rare disease pathogenesis, the organelle's contribution to common disease remains largely unknown. Hypothesizing that common genetic variants within Mendelian ciliopathy genes might contribute to common complex diseases pathogenesis, we performed association studies of 16,874 common genetic variants across 122 ciliary genes with 12 quantitative laboratory traits characteristic of ciliopathy syndromes in 452,593 individuals in the UK Biobank. We incorporated tissue-specific gene expression analysis, expression quantitative trait loci, and Mendelian disease phenotype information into our analysis and replicated our findings in meta-analysis. 101 statistically significant associations were identified across 42 of the 122 examined ciliary genes (including eight novel replicating associations). These ciliary genes were widely expressed in tissues relevant to the phenotypes being studied, and eQTL analysis revealed strong evidence for correlation between ciliary gene expression levels and laboratory traits. Perhaps most interestingly, our analysis identified different ciliary subcompartments as being specifically associated with distinct sets of phenotypes. Taken together, our data demonstrate the utility of a Mendelian pathway-based approach to genomic association studies, challenge the widely held belief that the cilium is an organelle important mainly in development and in rare syndromic disease pathogenesis, and provide a framework for the continued integration of common and rare disease genetics to provide insight into the pathophysiology of human diseases of immense public health burden.

Leveraging Clinical Intuition to Improve Accuracy of Phenotype-Driven Prioritization.

PURPOSE: Clinical intuition is commonly incorporated into the differential diagnosis as an assessment of the likelihood of candidate diagnoses based either on the patient population being seen in a specific clinic or on the signs and symptoms of the initial presentation. Algorithms to support diagnostic sequencing in individuals with a suspected rare genetic disease do not yet incorporate intuition and instead assume that each Mendelian disease has an equal pretest probability. METHODS: The LIRICAL algorithm calculates the likelihood ratio of clinical manifestations represented by Human Phenotype Ontology (HPO) terms to rank candidate diagnoses. The initial version of LIRICAL assumed an equal pretest probability for each disease in its calculation of the posttest probability (where the test is diagnostic exome or genome sequencing). We introduce Clinical Intuition for Likelihood Ratios (ClintLR), an extension of the LIRICAL algorithm that boosts the pretest probability of groups of related diseases deemed to be more likely. RESULTS: The average rank of the correct diagnosis in simulations using ClintLR showed a statistically significant improvement over a range of adjustment factors. CONCLUSION: ClintLR successfully encodes clinical intuition to improve ranking of rare diseases in diagnostic sequencing. ClintLR is freely available at https://github.com/TheJacksonLaboratory/ClintLR.

A heterozygous SPRY4 variant identified in female infertility characterized by reduced oocyte potential and early embryonic arrest.

STUDY QUESTION: Can novel genetic factors contributing to early embryonic arrest in infertile patients be identified, along with the underlying mechanisms of the pathogenic variant? SUMMARY ANSWER: We identified a heterozygous variant in the SPRY4 (sprouty RTK signaling antagonist 4) in infertile patients and conducted in vitro and in vivo studies to investigate the effects of the variant/deletion, highlighting its critical role in female reproductive health. WHAT IS KNOWN ALREADY: SPRY4 acts as a negative regulator of receptor tyrosine kinases (RTKs) and functions as a tumor suppressor. Its abnormal expression can lead to recurrent miscarriage by affecting trophoblast function. In mice, Spry4 knockout (KO) leads to craniofacial anomalies and growth defects. A human study links the SPRY4 variant to a male patient with isolated hypogonadotropic hypogonadism (IHH), hypothetically impacting gonadotropin-releasing hormone (GnRH) neurons, and causing reproductive dysfunctions. SPRY4 is thus potentially integral in regulating endocrine homeostasis and reproductive function. To date, no study has reported SPRY4 variants associated with female fertility, and a causal relationship has not been established with functional evidence. STUDY DESIGN, SIZE, DURATION: Whole-exome sequencing (WES) was performed in 392 infertile women who suffered from primary infertility of unknown reason, and the heterozygous SPRY4 variant were identified in one independent family. The infertile patients presenting were recruited from July 2017 to November 2023. PARTICIPANTS/MATERIALS, SETTING, METHODS: Women diagnosed with primary infertility were recruited from the Reproduction Center of Zhongshan Hospital, Fudan University. Genomic DNA was extracted from peripheral blood for WES analysis. The SPRY4 variant were identified through WES, in silico analysis, and variant screening. All variants were confirmed by Sanger sequencing. The effects of the variants were investigated in human embryonic kidney (HEK) 293T (HEK293T) cells via western blotting, and in mouse oocytes and embryos through complementary RNA (cRNA) injection, RNA sequencing, fluorescence, absorbance, and RT-qPCR assays. Gene function was further examined in Spry4 KO mice via histology, western blotting, ELISA, and RT-qPCR assays. MAIN RESULTS AND THE ROLE OF CHANCE: We identified a missense heterozygous pathogenic variant in SPRY4 (GRCh38, GenBank: NM_030964.5, c.157C>T p.(Arg53Trp), rs200531302) that reduces SPRY4 protein levels in HEK293T cells and disrupts the redox system and mitochondrial function in mouse oocyte, and perturbs developmental potential in mouse embryos. These phenotypes could be partially reversed by the exogenous addition of Nrf1 cRNA. Additionally, Spry4-/- mice exhibit ovarian oxidative stress and decreased ovarian function. LIMITATIONS, REASONS FOR CAUTION: Due to the limited WES data and population, we identified only one family with a SPRY4 mutation. The deeper mechanism and therapeutic strategy should be further investigated through mutant mice and recovery experiment. WIDER IMPLICATIONS OF THE FINDINGS: Our study has identified a pathogenic variant in SPRY4 associated with early embryonic arrest in humans. These findings enhance our understanding of the role of SPRY4 in early embryonic development and present a new genetic marker for female infertility. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the National Natural Science Foundation of China (82071643 and 82171655) and Natural Science Foundation of Shanghai (22ZR1456200). None of the authors have any competing interests. TRIAL REGISTRATION NUMBER: N/A.

A loss-of-function AGTR1 variant in a critically-ill infant with renal tubular dysgenesis: case presentation and literature review.

BACKGROUND: Renal tubular dysgenesis (RTD) is a severe disorder with poor prognosis significantly impacting the proximal tubules of the kidney while maintaining an anatomically normal gross structure. The genetic origin of RTD, involving variants in the ACE, REN, AGT, and AGTR1 genes, affects various enzymes or receptors within the Renin angiotensin system (RAS). This condition manifests prenatally with oligohydramninos and postnatally with persistent anuria, severe refractory hypotension, and defects in skull ossification. CASE PRESENTATION: In this report, we describe a case of a female patient who, despite receiving multi vasopressor treatment, experienced persistent hypotension, ultimately resulting in early death at five days of age. While there was a history of parental consanguinity, no reported family history of renal disease existed. Blood samples from the parents and the remaining DNA sample of the patient underwent Whole Genome Sequencing (WGS). The genetic analysis revealed a rare homozygous loss of function variant (NM_000685.5; c.415C > T; p.Arg139*) in the Angiotensin II Receptor Type 1 (AGTR1) gene. CONCLUSION: This case highlights the consequence of loss-of-function variants in AGTR1 gene leading to RTD, which is characterized by high mortality rate at birth or during the neonatal period. Furthermore, we provide a comprehensive review of previously reported variants in the AGTR1 gene, which is the least encountered genetic cause of RTD, along with their associated clinical features.

Attention challenges in Kabuki syndrome.

BACKGROUND: Understanding the specific neurobehavioural profile of rare genetic diseases enables clinicians to provide the best possible care for patients and families, including prognostic and treatment advisement. Previous studies suggested that a subset of individuals with Kabuki syndrome (KS), a genetic disorder causing intellectual disability and other neurodevelopmental phenotypes, have attentional deficits. However, these studies looked at relatively small numbers of molecularly confirmed cases and often used retrospective clinical data instead of standardised assessments. METHODS: Fifty-five individuals or caregivers of individuals with molecularly confirmed KS completed assessments to investigate behaviour and adaptive function. Additionally, information was collected on 23 unaffected biological siblings as controls. RESULTS: Attention Problems in children was the only behavioural category that, when averaged, was clinically significant, with the individual scores of nearly 50% of the children with KS falling in the problematic range. Children with KS scored significantly higher than their unaffected sibling on nearly all behavioural categories. A significant correlation was found between Attention Problems scores and adaptive function scores (P = 0.032), which was not explained by lower general cognitive ability. CONCLUSIONS: We found that the rates of children with attentional deficits are much more elevated than would be expected in the general population, and that attention challenges are negatively correlated with adaptive function. When averaged across KS participants, none of the behavioural categories were in the clinically significant range except Attention Problems for children, which underscores the importance of clinicians screening for attention deficit hyperactivity disorder (ADHD) in children with KS.

Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease.

Non-coding transcriptional regulatory elements are critical for controlling the spatiotemporal expression of genes. Here, we demonstrate that the sizes and number of enhancers linked to a gene reflect its disease pathogenicity. Moreover, genes with redundant enhancer domains are depleted of cis-acting genetic variants that disrupt gene expression, and they are buffered against the effects of disruptive non-coding mutations. Our results demonstrate that dosage-sensitive genes have evolved a robustness to the disruptive effects of genetic variation by expanding their regulatory domains. This solves a puzzle about why genes associated with human disease are depleted of cis-eQTLs (cis-expression quantitative trait loci), suggesting that this relationship might complicate gene identification in causal genome-wide association studies (GWASs) using eQTL information, and establishes a framework for identifying non-coding regulatory variation with phenotypic consequences.

Mutations in MYLPF Cause a Novel Segmental Amyoplasia that Manifests as Distal Arthrogryposis.

We identified ten persons in six consanguineous families with distal arthrogryposis (DA) who had congenital contractures, scoliosis, and short stature. Exome sequencing revealed that each affected person was homozygous for one of two different rare variants (c.470G>T [p.Cys157Phe] or c.469T>C [p.Cys157Arg]) affecting the same residue of myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF). In a seventh family, a c.487G>A (p.Gly163Ser) variant in MYLPF arose de novo in a father, who transmitted it to his son. In an eighth family comprised of seven individuals with dominantly inherited DA, a c.98C>T (p.Ala33Val) variant segregated in all four persons tested. Variants in MYLPF underlie both dominant and recessively inherited DA. Mylpf protein models suggest that the residues associated with dominant DA interact with myosin whereas the residues altered in families with recessive DA only indirectly impair this interaction. Pathological and histological exam of a foot amputated from an affected child revealed complete absence of skeletal muscle (i.e., segmental amyoplasia). To investigate the mechanism for this finding, we generated an animal model for partial MYLPF impairment by knocking out zebrafish mylpfa. The mylpfa mutant had reduced trunk contractile force and complete pectoral fin paralysis, demonstrating that mylpf impairment most severely affects limb movement. mylpfa mutant muscle weakness was most pronounced in an appendicular muscle and was explained by reduced myosin activity and fiber degeneration. Collectively, our findings demonstrate that partial loss of MYLPF function can lead to congenital contractures, likely as a result of degeneration of skeletal muscle in the distal limb.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome (XGS; MIM# 615829) is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Monogenic focal segmental glomerulosclerosis: A conceptual framework for identification and management of a heterogeneous disease.

Focal segmental glomerulosclerosis (FSGS) is not a disease, rather a pattern of histological injury occurring from a variety of causes. The exact pathogenesis has yet to be fully elucidated but is likely varied based on the type of injury and the primary target of that injury. However, the approach to treatment is often based on the degree of podocyte foot process effacement and clinical presentation without sufficient attention paid to etiology. In this regard, there are many monogenic causes of FSGS with variable presentation from nephrotic syndrome with histological features of primary podocytopathy to more modest degrees of proteinuria with limited evidence of podocyte foot process injury. It is likely that genetic causes are largely underdiagnosed, as the role and the timing of genetic testing in FSGS is not established and genetic counseling, testing options, and interpretation of genotype in the context of phenotype may be outside the scope of practice for both nephrologists and geneticists. Yet most clinicians believe that a genetic diagnosis can lead to targeted therapy, limit the use of high-dose corticosteroids as a therapeutic trial, and allow the prediction of the natural history and risk for recurrence in the transplanted kidney. In this manuscript, we emphasize that genetic FSGS is not monolithic in its presentation, opine on the importance of genetic testing and provide an algorithmic approach to deployment of genetic testing in a timely fashion when faced with a patient with FSGS.

Reorganization of inter-chromosomal interactions in the 2q37-deletion syndrome.

Chromosomes occupy distinct interphase territories in the three-dimensional nucleus. However, how these chromosome territories are arranged relative to one another is poorly understood. Here, we investigated the inter-chromosomal interactions between chromosomes 2q, 12, and 17 in human mesenchymal stem cells (MSCs) and MSC-derived cell types by DNA-FISH We compared our findings in normal karyotypes with a three-generation family harboring a 2q37-deletion syndrome, featuring a heterozygous partial deletion of histone deacetylase 4 (HDAC4) on chr2q37. In normal karyotypes, we detected stable, recurring arrangements and interactions between the three chromosomal territories with a tissue-specific interaction bias at certain loci. These inter-chromosomal interactions were confirmed by Hi-C. Interestingly, the disease-related HDAC4 deletion resulted in displaced inter-chromosomal arrangements and altered interactions between the deletion-affected chromosome 2 and chromosome 12 and/or 17 in 2q37-deletion syndrome patients. Our findings provide evidence for a direct link between a structural chromosomal aberration and altered interphase architecture that results in a nuclear configuration, supporting a possible molecular pathogenesis.

Expanding the Boundaries of RNA Sequencing as a Diagnostic Tool for Rare Mendelian Disease.

Gene-panel and whole-exome analyses are now standard methodologies for mutation detection in Mendelian disease. However, the diagnostic yield achieved is at best 50%, leaving the genetic basis for disease unsolved in many individuals. New approaches are thus needed to narrow the diagnostic gap. Whole-genome sequencing is one potential strategy, but it currently has variant-interpretation challenges, particularly for non-coding changes. In this study we focus on transcriptome analysis, specifically total RNA sequencing (RNA-seq), by using monogenetic neuromuscular disorders as proof of principle. We examined a cohort of 25 exome and/or panel "negative" cases and provided genetic resolution in 36% (9/25). Causative mutations were identified in coding and non-coding exons, as well as in intronic regions, and the mutational pathomechanisms included transcriptional repression, exon skipping, and intron inclusion. We address a key barrier of transcriptome-based diagnostics: the need for source material with disease-representative expression patterns. We establish that blood-based RNA-seq is not adequate for neuromuscular diagnostics, whereas myotubes generated by transdifferentiation from an individual's fibroblasts accurately reflect the muscle transcriptome and faithfully reveal disease-causing mutations. Our work confirms that RNA-seq can greatly improve diagnostic yield in genetically unresolved cases of Mendelian disease, defines strengths and challenges of the technology, and demonstrates the suitability of cell models for RNA-based diagnostics. Our data set the stage for development of RNA-seq as a powerful clinical diagnostic tool that can be applied to the large population of individuals with undiagnosed, rare diseases and provide a framework for establishing minimally invasive strategies for doing so.

Deleterious Variation in BRSK2 Associates with a Neurodevelopmental Disorder.

Developmental delay and intellectual disability (DD and ID) are heterogeneous phenotypes that arise in many rare monogenic disorders. Because of this rarity, developing cohorts with enough individuals to robustly identify disease-associated genes is challenging. Social-media platforms that facilitate data sharing among sequencing labs can help to address this challenge. Through one such tool, GeneMatcher, we identified nine DD- and/or ID-affected probands with a rare, heterozygous variant in the gene encoding the serine/threonine-protein kinase BRSK2. All probands have a speech delay, and most present with intellectual disability, motor delay, behavioral issues, and autism. Six of the nine variants are predicted to result in loss of function, and computational modeling predicts that the remaining three missense variants are damaging to BRSK2 structure and function. All nine variants are absent from large variant databases, and BRSK2 is, in general, relatively intolerant to protein-altering variation among humans. In all six probands for whom parents were available, the mutations were found to have arisen de novo. Five of these de novo variants were from cohorts with at least 400 sequenced probands; collectively, the cohorts span 3,429 probands, and the observed rate of de novo variation in these cohorts is significantly higher than the estimated background-mutation rate (p = 2.46 × 10-6). We also find that exome sequencing provides lower coverage and appears less sensitive to rare variation in BRSK2 than does genome sequencing; this fact most likely reduces BRSK2's visibility in many clinical and research sequencing efforts. Altogether, our results implicate damaging variation in BRSK2 as a source of neurodevelopmental disease.

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.