Refining clinically relevant parameters for mis-splicing risk in shortened introns with donor-to-branchpoint space constraint.

Intronic deletions that critically shorten donor-to-branchpoint (D-BP) distance of a precursor mRNA impose biophysical space constraint on assembly of the U1/U2 spliceosomal complex, leading to canonical splicing failure. Here we use a series of ß-globin (HBB) gene constructs with intron 1 deletions to define D-BP lengths that present low/no risk of mis-splicing and lengths which are critically short and likely elicit clinically relevant mis-splicing. We extend our previous observation in EMD intron 5 of 46 nt as the minimal productive D-BP length, demonstrating spliceosome assembly constraint persists at D-BP lengths of 47-56 nt. We exploit the common HBB exon 1 ß-thalassemia variant that strengthens a cryptic donor (NM_000518.5(HBB):c.79G > A) to provide a simple barometer for the earliest signs of space constraint, via cryptic donor activation. For clinical evaluation of intronic deletions, we assert D-BP lengths > 60 nt present low mis-splicing risk while space constraint increases exponentially with D-BP lengths < 55 nt, with critical risk and profound splicing abnormalities with D-BP lengths < 50 nt.

De novo variants in the non-coding spliceosomal snRNA gene RNU4-2 are a frequent cause of syndromic neurodevelopmental disorders.

Around 60% of individuals with neurodevelopmental disorders (NDD) remain undiagnosed after comprehensive genetic testing, primarily of protein-coding genes1. Increasingly, large genome-sequenced cohorts are improving our ability to discover new diagnoses in the non-coding genome. Here, we identify the non-coding RNA RNU4-2 as a novel syndromic NDD gene. RNU4-2 encodes the U4 small nuclear RNA (snRNA), which is a critical component of the U4/U6.U5 tri-snRNP complex of the major spliceosome2. We identify an 18 bp region of RNU4-2 mapping to two structural elements in the U4/U6 snRNA duplex (the T-loop and Stem III) that is severely depleted of variation in the general population, but in which we identify heterozygous variants in 119 individuals with NDD. The vast majority of individuals (77.3%) have the same highly recurrent single base-pair insertion (n.64_65insT). We estimate that variants in this region explain 0.41% of individuals with NDD. We demonstrate that RNU4-2 is highly expressed in the developing human brain, in contrast to its contiguous counterpart RNU4-1 and other U4 homologs, supporting RNU4-2's role as the primary U4 transcript in the brain. Overall, this work underscores the importance of non-coding genes in rare disorders. It will provide a diagnosis to thousands of individuals with NDD worldwide and pave the way for the development of effective treatments for these individuals.

Systematic identification of disease-causing promoter and untranslated region variants in 8,040 undiagnosed individuals with rare disease.

Background: Both promoters and untranslated regions (UTRs) have critical regulatory roles, yet variants in these regions are largely excluded from clinical genetic testing due to difficulty in interpreting pathogenicity. The extent to which these regions may harbour diagnoses for individuals with rare disease is currently unknown. Methods: We present a framework for the identification and annotation of potentially deleterious proximal promoter and UTR variants in known dominant disease genes. We use this framework to annotate de novo variants (DNVs) in 8,040 undiagnosed individuals in the Genomics England 100,000 genomes project, which were subject to strict region-based filtering, clinical review, and validation studies where possible. In addition, we performed region and variant annotation-based burden testing in 7,862 unrelated probands against matched unaffected controls. Results: We prioritised eleven DNVs and identified an additional variant overlapping one of the eleven. Ten of these twelve variants (82%) are in genes that are a strong match to the individual's phenotype and six had not previously been identified. Through burden testing, we did not observe a significant enrichment of potentially deleterious promoter and/or UTR variants in individuals with rare disease collectively across any of our region or variant annotations. Conclusions: Overall, we demonstrate the value of screening promoters and UTRs to uncover additional diagnoses for previously undiagnosed individuals with rare disease and provide a framework for doing so without dramatically increasing interpretation burden.

SpliceVault predicts the precise nature of variant-associated mis-splicing.

Even for essential splice-site variants that are almost guaranteed to alter mRNA splicing, no current method can reliably predict whether exon-skipping, cryptic activation or multiple events will result, greatly complicating clinical interpretation of pathogenicity. Strikingly, ranking the four most common unannotated splicing events across 335,663 reference RNA-sequencing (RNA-seq) samples (300K-RNA Top-4) predicts the nature of variant-associated mis-splicing with 92% sensitivity. The 300K-RNA Top-4 events correctly identify 96% of exon-skipping events and 86% of cryptic splice sites for 140 clinical cases subject to RNA testing, showing higher sensitivity and positive predictive value than SpliceAI. Notably, RNA re-analyses showed we had missed 300K-RNA Top-4 events for several clinical cases tested before the development of this empirical predictive method. Simply, mis-splicing events that happen around a splice site in RNA-seq data are those most likely to be activated by a splice-site variant. The SpliceVault web portal allows users easy access to 300K-RNA for informed splice-site variant interpretation and classification.

Prevalence, parameters, and pathogenic mechanisms for splice-altering acceptor variants that disrupt the AG exclusion zone.

Predicting the pathogenicity of acceptor splice-site variants outside the essential AG is challenging, due to high sequence diversity of the extended splice-site region. Critical analysis of 24,445 intronic extended acceptor splice-site variants reported in ClinVar and the Leiden Open Variation Database (LOVD) demonstrates 41.9% of pathogenic variants create an AG dinucleotide between the predicted branchpoint and acceptor (AG-creating variants in the AG exclusion zone), 28.4% result in loss of a pyrimidine at the -3 position, and 15.1% result in loss of one or more pyrimidines in the polypyrimidine tract. Pathogenicity of AG-creating variants was highly influenced by their position. We define a high-risk zone for pathogenicity: > 6 nucleotides downstream of the predicted branchpoint and >5 nucleotides upstream from the acceptor, where 93.1% of pathogenic AG-creating variants arise and where naturally occurring AG dinucleotides are concordantly depleted (5.8% of natural AGs). SpliceAI effectively predicts pathogenicity of AG-creating variants, achieving 95% sensitivity and 69% specificity. We highlight clinical examples showing contrasting mechanisms for mis-splicing arising from AG variants: (1) cryptic acceptor created; (2) splicing silencer created: an introduced AG silences the acceptor, resulting in exon skipping, intron retention, and/or use of an alternative existing cryptic acceptor; and (3) splicing silencer disrupted: loss of a deep intronic AG activates inclusion of a pseudo-exon. In conclusion, we establish AG-creating variants as a common class of pathogenic extended acceptor variant and outline factors conferring critical risk for mis-splicing for AG-creating variants in the AG exclusion zone, between the branchpoint and acceptor.

Empirical prediction of variant-activated cryptic splice donors using population-based RNA-Seq data.

Predicting which cryptic-donors may be activated by a splicing variant in patient DNA is notoriously difficult. Through analysis of 5145 cryptic-donors (versus 86,963 decoy-donors not used; any GT or GC), we define an empirical method predicting cryptic-donor activation with 87% sensitivity and 95% specificity. Strength (according to four algorithms) and proximity to the annotated-donor appear important determinants of cryptic-donor activation. However, other factors such as splicing regulatory elements, which are difficult to identify, play an important role and are likely responsible for current prediction inaccuracies. We find that the most frequently recurring natural mis-splicing events at each exon-intron junction, summarised over 40,233 RNA-sequencing samples (40K-RNA), predict with accuracy which cryptic-donor will be activated in rare disease. 40K-RNA provides an accurate, evidence-based method to predict variant-activated cryptic-donors in genetic disorders, assisting pathology consideration of possible consequences of a variant for the encoded protein and RNA diagnostic testing strategies.

Standardized practices for RNA diagnostics using clinically accessible specimens reclassifies 75% of putative splicing variants.

PURPOSE: Genetic variants causing aberrant premessenger RNA splicing are increasingly being recognized as causal variants in genetic disorders. In this study, we devise standardized practices for polymerase chain reaction (PCR)-based RNA diagnostics using clinically accessible specimens (blood, fibroblasts, urothelia, biopsy). METHODS: A total of 74 families with diverse monogenic conditions (31% prenatal-congenital onset, 47% early childhood, and 22% teenage-adult onset) were triaged into PCR-based RNA testing, with comparative RNA sequencing for 19 cases. RESULTS: Informative RNA assay data were obtained for 96% of cases, enabling variant reclassification for 75% variants that can be used for genetic counseling (71%), to inform clinical care (32%) and prenatal counseling (41%). Variant-associated mis-splicing was highly reproducible for 28 cases with samples from ≥2 affected individuals or heterozygotes and 10 cases with ≥2 biospecimens. PCR amplicons encompassing another segregated heterozygous variant was vital for clinical interpretation of 22 of 79 variants to phase RNA splicing events and discern complete from partial mis-splicing. CONCLUSION: RNA diagnostics enabled provision of a genetic diagnosis for 64% of recruited cases. PCR-based RNA diagnostics has capacity to analyze 81.3% of clinically significant genes, with long amplicons providing an advantage over RNA sequencing to phase RNA splicing events. The Australasian Consortium for RNA Diagnostics (SpliceACORD) provide clinically-endorsed, standardized protocols and recommendations for interpreting RNA assay data.

Pro-inflammatory dopamine-2 receptor-specific T cells in paediatric movement and psychiatric disorders.

OBJECTIVES: A dysregulated inflammatory response against the dopamine-2 receptor (D2R) has been implicated in movement and psychiatric disorders. D2R antibodies were previously reported in a subset of these patients; however, the role of T cells in these disorders remains unknown. Our objective was to identify and characterise pro-inflammatory D2R-specific T cells in movement and psychiatric disorders. METHODS: Blood from paediatric patients with movement and psychiatric disorders of suspected autoimmune and neurodevelopmental aetiology (n = 24) and controls (n = 16) was cultured in vitro with a human D2R peptide library, and D2R-specific T cells were identified by flow cytometric quantification of CD4+CD25+CD134+ T cells. Cytokine secretion was analysed using a cytometric bead array and ELISA. HLA genotypes were examined in D2R-specific T-cell-positive patients. D2R antibody seropositivity was determined using a flow cytometry live cell-based assay. RESULTS: Three immunodominant regions of D2R, amino acid (aa)121-131, aa171-181 and aa396-416, specifically activated CD4+ T cells in 8/24 patients. Peptides corresponding to these regions were predicted to bind with high affinity to the HLA of the eight positive patients and had also elicited the secretion of pro-inflammatory cytokines IL-2, IFN- γ, TNF, IL-6, IL-17A and IL-17F. All eight patients were seronegative for D2R antibodies. CONCLUSION: Autoreactive D2R-specific T cells and a pro-inflammatory Th1 and Th17 cytokine profile characterise a subset of paediatric patients with movement and psychiatric disorders, further underpinning the theory of immune dysregulation in these disorders. These findings offer new perspectives into the neuroinflammatory mechanisms of movement and psychiatric disorders and can influence patient diagnosis and treatment.

Gene discovery informatics toolkit defines candidate genes for unexplained infertility and prenatal or infantile mortality.

Despite a recent surge in novel gene discovery, genetic causes of prenatal-lethal phenotypes remain poorly defined. To advance gene discovery in prenatal-lethal disorders, we created an easy-to-mine database integrating known human phenotypes with inheritance pattern, scores of genetic constraint, and murine and cellular knockout phenotypes-then critically assessed defining features of known prenatal-lethal genes, among 3187 OMIM genes, and relative to 16,009 non-disease genes. While around one-third (39%) of protein-coding genes are essential for murine development, we curate only 3% (624) of human protein-coding genes linked currently to prenatal/infantile lethal disorders. 75% prenatal-lethal genes are linked to developmental lethality in knockout mice, compared to 54% for all OMIM genes and 34% among non-disease genes. Genetic constraint correlates with inheritance pattern (autosomal recessive <90% of recessive genes show neither missense nor loss-of-function constraint, even for prenatal-lethal genes. Detailed ontology mapping for 624 prenatal-lethal genes shows marked enrichment among dominant genes for nuclear proteins with roles in RNA/DNA biology, with recessive genes enriched in cytoplasmic (mitochondrial) metabolic proteins. We conclude that genes without genetic constraint should not be excluded as potential novel disease genes, and especially for recessive conditions (<10% constrained). Prenatal lethal genes are 5.9-fold more likely to be associated with a lethal murine phenotype than non-disease genes. Cell essential genes are largely a subset of mouse-lethal genes, notably under-represented among known OMIM genes, and strong candidates for gamete/embryo non-viability. We therefore curate 3435 'candidate developmental lethal' human genes: essential for murine development or cellular viability, not yet linked to human disorders, presenting strong candidates for unexplained infertility and prenatal/infantile mortality.