Base-specific mutational intolerance near splice sites clarifies the role of nonessential splice nucleotides

Modulation of prion protein expression through cryptic splice site manipulation

Lowering expression of prion protein (PrP) is a well-validated therapeutic strategy in prion disease, but additional modalities are urgently needed. In other diseases, small molecules have proven capable of modulating pre-mRNA splicing, sometimes by forcing inclusion of cryptic exons that reduce gene expression. Here, we characterize a cryptic exon located in human PRNP's sole intron and evaluate its potential to reduce PrP expression through incorporation into the 5' untranslated region. This exon is homologous to exon 2 in nonprimate species but contains a start codon that would yield an upstream open reading frame with a stop codon prior to a splice site if included in PRNP mRNA, potentially downregulating PrP expression through translational repression or nonsense-mediated decay. We establish a minigene transfection system and test a panel of splice site alterations, identifying mutants that reduce PrP expression by as much as 78%. Our findings nominate a new therapeutic target for lowering PrP.

Achieving an optimal pregnancy outcome through the combined utilization of micro-TESE and ICSI in cryptorchidism associated with a non-canonical splicing variant in RXFP2

PURPOSE

To identify the genetic cause of a cryptorchidism patient carrying a non-canonical splicing variant highlighted by SPCards platform in RXFP2 and to provide a comprehensive overview of RXFP2 variants with cryptorchidism correlation.

METHODS

We identified a homozygous non-canonical splicing variant by whole-exome sequencing and Sanger sequencing in a case with cryptorchidism and non-obstructive azoospermia (NOA). As the pathogenicity of this non-canonical splicing variant remained unclear, we initially utilized the SPCards platform to predict its pathogenicity. Subsequently, we employed a minigene splicing assay to further evaluate the influence of the identified splicing variant. Microdissection testicular sperm extraction (micro-TESE) combined with intracytoplasmic sperm injection (ICSI) was performed. PubMed and Human Genome Variant Database (HGMD) were queried to search for RXFP2 variants.

RESULTS

We identified a homozygous non-canonical splicing variant (NM_130806: c.1376-12A > G) in RXFP2, and confirmed this variant caused aberrant splicing of exons 15 and 16 of the RXFP2 gene: 11 bases were added in front of exon 16, leading to an abnormal transcript initiation and a frameshift. Fortunately, the patient successfully obtained his biological offspring through micro-TESE combined with ICSI. Four cryptorchidism-associated variants in RXFP2 from 90 patients with cryptorchidism were identified through a literature search in PubMed and HGMD, with different inheritance patterns.

CONCLUSION

This is the first cryptorchidism case carrying a novel causative non-canonical splicing RXFP2 variant. The combined approach of micro-TESE and ICSI contributed to an optimal pregnancy outcome. Our literature review demonstrated that RXFP2 variants caused cryptorchidism in a recessive inheritance pattern, rather than a dominant pattern.

A benchmarked, high-efficiency prime editing platform for multiplexed dropout screening

Prime editing installs precise edits into the genome with minimal unwanted byproducts, but low and variable editing efficiencies have complicated application of the approach to high-throughput functional genomics. Leveraging several recent advances, we assembled a prime editing platform capable of high-efficiency substitution editing across a set of engineered prime editing guide RNAs (epegRNAs) and corresponding target sequences (80% median intended editing). Then, using a custom library of 240,000 epegRNAs targeting >17,000 codons with 175 different substitution types, we benchmarked our platform for functional interrogation of small substitution variants (1-3 nucleotides) targeted to essential genes. Resulting data identified negative growth phenotypes for nonsense mutations targeted to ~8,000 codons, and comparing those phenotypes to results from controls demonstrated high specificity. We also observed phenotypes for synonymous mutations that disrupted splice site motifs at 3' exon boundaries. Altogether, we establish and benchmark a high-throughput prime editing approach for functional characterization of genetic variants with simple readouts from multiplexed experiments.

Recurring homozygous ACTN2 variant (p.Arg506Gly) causes a recessive myopathy

OBJECTIVE

ACTN2, encoding alpha-actinin-2, is essential for cardiac and skeletal muscle sarcomeric function. ACTN2 variants are a known cause of cardiomyopathy without skeletal muscle involvement. Recently, specific dominant monoallelic variants were reported as a rare cause of core myopathy of variable clinical onset, although the pathomechanism remains to be elucidated. The possibility of a recessively inherited ACTN2-myopathy has also been proposed in a single series.

METHODS

We provide clinical, imaging, and histological characterization of a series of patients with a novel biallelic ACTN2 variant.

RESULTS

We report seven patients from five families with a recurring biallelic variant in ACTN2: c.1516A>G (p.Arg506Gly), all manifesting with a consistent phenotype of asymmetric, progressive, proximal, and distal lower extremity predominant muscle weakness. None of the patients have cardiomyopathy or respiratory insufficiency. Notably, all patients report Palestinian ethnicity, suggesting a possible founder ACTN2 variant, which was confirmed through haplotype analysis in two families. Muscle biopsies reveal an underlying myopathic process with disruption of the intermyofibrillar architecture, Type I fiber predominance and atrophy. MRI of the lower extremities demonstrate a distinct pattern of asymmetric muscle involvement with selective involvement of the hamstrings and adductors in the thigh, and anterior tibial group and soleus in the lower leg. Using an in vitro splicing assay, we show that c.1516A>G ACTN2 does not impair normal splicing.

INTERPRETATION

This series further establishes ACTN2 as a muscle disease gene, now also including variants with a recessive inheritance mode, and expands the clinical spectrum of actinopathies to adult-onset progressive muscle disease.

Prioritizing de novo potential non-canonical splicing variants in neurodevelopmental disorders

BACKGROUND

Genomic variants outside of the canonical splicing site (±2) may generate abnormal mRNA splicing, which are defined as non-canonical splicing variants (NCSVs). However, the clinical interpretation of NCSVs in neurodevelopmental disorders (NDDs) is largely unknown.

METHODS

We investigated the contribution of NCSVs to NDDs from 345,787 de novo variants (DNVs) in 47,574 patients with NDDs. We performed functional enrichment and protein-protein interaction analysis to assess the association between genes carrying prioritised NCSVs and NDDs. Minigene was used to validate the impact of NCSVs on mRNA splicing.

FINDINGS

We observed significantly more NCSVs (p = 0.02, odds ratio [OR] = 2.05) among patients with NDD than in controls. Both canonical splicing variants (CSVs) and NCSVs contributed to an equal proportion of patients with NDD (0.76% vs. 0.82%). The candidate genes carrying NCSVs were associated with glutamatergic synapse and chromatin remodelling. Minigene successfully validated 59 of 79 (74.68%) NCSVs that led to abnormal splicing in 40 candidate genes, and 9 of the genes (ARID1B, KAT6B, TCF4, SMARCA2, SHANK3, PDHA1, WDR45, SCN2A, SYNGAP1) harboured recurrent NCSVs with the same variant present in more than two unrelated patients with NDD. Moreover, 36 of 59 (61.02%) NCSVs are novel clinically relevant variants, including 34 unreported and 2 clinically conflicting interpretations or of uncertain significance NCSVs in the ClinVar database.

INTERPRETATION

This study highlights the common pathology and clinical importance of NCSVs in unsolved patients with NDD.

FUNDING

The present study was funded by grants from the National Natural Science Foundation of China, China Postdoctoral Science Foundation, the Hunan Youth Science and Technology Innovation Talent Project, the Provincial Natural Science Foundation of Hunan, The Scientific Research Program of FuRong laboratory, and the Natural Science Project of the University of Anhui Province.

Modulation of prion protein expression through cryptic splice site manipulation

Lowering expression of prion protein (PrP) is a well-validated therapeutic strategy in prion disease, but additional modalities are urgently needed. In other diseases, small molecules have proven capable of modulating pre-mRNA splicing, sometimes by forcing inclusion of cryptic exons that reduce gene expression. Here, we characterize a cryptic exon located in human PRNP's sole intron and evaluate its potential to reduce PrP expression through incorporation into the 5' untranslated region (5'UTR). This exon is homologous to exon 2 in non-primate species, but contains a start codon that would yield an upstream open reading frame (uORF) with a stop codon prior to a splice site if included in PRNP mRNA, potentially downregulating PrP expression through translational repression or nonsense-mediated decay. We establish a minigene transfection system and test a panel of splice site alterations, identifying mutants that reduce PrP expression by as much as 78%. Our findings nominate a new therapeutic target for lowering PrP.

Analysis of 200 unrelated individuals with a constitutional NF1 deep intronic pathogenic variant reveals that variants flanking the alternatively spliced NF1 exon 31 [23a] cause a classical neurofibromatosis type 1 phenotype while altering predominantly NF1 isoform type II

Neurofibromatosis type 1 results from loss-of-function NF1 pathogenic variants (PVs). Up to 30% of all NF1 PVs disrupt mRNA splicing, including deep intronic variants. Here, we retrospectively investigated the spectrum of NF1 deep intronic PVs in a cohort of 8,090 unrelated individuals from the University of Alabama at Birmingham (UAB) dataset with a molecularly confirmed neurofibromatosis type 1. All variants were identified through their effect on the NF1 transcript, followed by variant characterization at the DNA-level. A total of 68 distinct variants, which were ≥ 20 nucleotides away from the closest exon-intron junction, were identified in 2.5% unrelated individuals with NF1 (200/8,090). Nine different pathogenic splice variants, identified in 20 probands, led to exonization of different parts of intron 30 [23.2] or 31 [23a]. The two major NF1 transcript isoforms, distinguished by the absence (type I) or presence (type II) of the alternatively spliced cassette exon 31 [23a], are equally expressed in blood in control individuals without NF1 or NF1-affected individuals carrying their PV not in the introns flanking exon 31 [23a]. By fragment and cloning analysis we demonstrated that the exonization of intron 31 [23a] sequences due to deep intronic PV predominantly affects the NF1 isoform II. Seven additional (likely) pathogenic NF1 deep intronic variants not observed in the UAB dataset were found by classification of 36 variants identified by a literature search. Hence, the unique list of these 75 deep intronic (likely) PVs should be included in any comprehensive NF1 testing strategy.

A spliceosome-associated gene signature aids in predicting prognosis and tumor microenvironment of hepatocellular carcinoma

Splicing alterations have been shown to be key tumorigenesis drivers. In this study, we identified a novel spliceosome-related genes (SRGs) signature to predict the overall survival (OS) of patients with hepatocellular carcinoma (HCC). A total of 25 SRGs were identified from the GSE14520 dataset (training set). Univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were utilized to construct the signature using genes with predictive significance. We then constructed a risk model using six SRGs (BUB3, IGF2BP3, RBM3, ILF3, ZC3H13, and CCT3). The reliability and predictive power of the gene signature were validated in two validation sets (TCGA and GSE76427 dataset). Patients in training and validation sets were divided into high and low-risk groups based on the gene signature. Patients in high-risk groups exhibited a poorer OS than in low-risk groups both in the training set and two validation sets. Next, risk score, BCLC staging, TNM staging, and multinodular were combined in a nomogram for OS prediction, and the decision curve analysis (DCA) curve exhibited the excellent prediction performance of the nomogram. The functional enrichment analyses demonstrated high-risk score patients were closely related to multiple oncology characteristics and invasive-related pathways, such as Cell cycle, DNA replication, and Spliceosome. Different compositions of the tumor microenvironment and immunocyte infiltration ratio might contribute to the prognostic difference between high and low-risk score groups. In conclusion, a spliceosome-related six-gene signature exhibited good performance for predicting the OS of patients with HCC, which may aid in clinical decision-making for individual treatment.

Mendelian inheritance revisited: dominance and recessiveness in medical genetics

Understanding the consequences of genotype for phenotype (which ranges from molecule-level effects to whole-organism traits) is at the core of genetic diagnostics in medicine. Many measures of the deleteriousness of individual alleles exist, but these have limitations for predicting the clinical consequences. Various mechanisms can protect the organism from the adverse effects of functional variants, especially when the variant is paired with a wild type allele. Understanding why some alleles are harmful in the heterozygous state - representing dominant inheritance - but others only with the biallelic presence of pathogenic variants - representing recessive inheritance - is particularly important when faced with the deluge of rare genetic alterations identified by high throughput DNA sequencing. Both awareness of the specific quantitative and/or qualitative effects of individual variants and the elucidation of allelic and non-allelic interactions are essential to optimize genetic diagnosis and counselling.

Machine learning-based identification of SOX10 as an immune regulator of macrophage in gliomas

Gliomas, originating from the glial cells, are the most lethal type of primary tumors in the central nervous system. Standard treatments like surgery have not significantly improved the prognosis of glioblastoma patients. Recently, immune therapy has become a novel and effective option. As a conserved group of transcriptional regulators, the Sry-type HMG box (SOX) family has been proved to have a correlation with numerous diseases. Based on the large-scale machine learning, we found that the SOX family, with significant immune characteristics and genomic profiles, can be divided into two distinct clusters in gliomas, among which SOX10 was identified as an excellent immune regulator of macrophage in gliomas. The high expression of SOX10 is related to a shorter OS in LGG, HGG, and pan-cancer groups but benefited from the immunotherapy. It turned out in single-cell sequencing that SOX10 is high in neurons, M1 macrophages, and neural stem cells. Also, macrophages are found to be elevated in the SOX10 high-expression group. SOX10 has a positive correlation with macrophage cytokine production and negative regulation of macrophages' chemotaxis and migration. In conclusion, our study demonstrates the outstanding cluster ability of the SOX family, indicating that SOX10 is an immune regulator of macrophage in gliomas, which can be an effective target for glioma immunotherapy.

Selection of Diagnostically Significant Regions of the SLC26A4 Gene Involved in Hearing Loss

Screening pathogenic variants in the SLC26A4 gene is an important part of molecular genetic testing for hearing loss (HL) since they are one of the common causes of hereditary HL in many populations. However, a large size of the SLC26A4 gene (20 coding exons) predetermines the difficulties of its complete mutational analysis, especially in large samples of patients. In addition, the regional or ethno-specific prevalence of SLC26A4 pathogenic variants has not yet been fully elucidated, except variants c.919-2A>G and c.2168A>G (p.His723Arg), which have been proven to be most common in Asian populations. We explored the distribution of currently known pathogenic and likely pathogenic (PLP) variants across the SLC26A4 gene sequence presented in the Deafness Variation Database for the selection of potential diagnostically important parts of this gene. As a result of this bioinformatic analysis, we found that molecular testing ten SLC26A4 exons (4, 6, 10, 11, 13−17 and 19) with flanking intronic regions can provide a diagnostic rate of 61.9% for all PLP variants in the SLC26A4 gene. The primary sequencing of these SLC26A4 regions may be applied as an initial effective diagnostic testing in samples of patients of unknown ethnicity or as a subsequent step after the targeted testing of already-known ethno- or region-specific pathogenic SLC26A4 variants.

A systematic analysis of splicing variants identifies new diagnoses in the 100,000 Genomes Project

BACKGROUND

Genomic variants which disrupt splicing are a major cause of rare genetic diseases. However, variants which lie outside of the canonical splice sites are difficult to interpret clinically. Improving the clinical interpretation of non-canonical splicing variants offers a major opportunity to uplift diagnostic yields from whole genome sequencing data.

METHODS

Here, we examine the landscape of splicing variants in whole-genome sequencing data from 38,688 individuals in the 100,000 Genomes Project and assess the contribution of non-canonical splicing variants to rare genetic diseases. We use a variant-level constraint metric (the mutability-adjusted proportion of singletons) to identify constrained functional variant classes near exon-intron junctions and at putative splicing branchpoints. To identify new diagnoses for individuals with unsolved rare diseases in the 100,000 Genomes Project, we identified individuals with de novo single-nucleotide variants near exon-intron boundaries and at putative splicing branchpoints in known disease genes. We identified candidate diagnostic variants through manual phenotype matching and confirmed new molecular diagnoses through clinical variant interpretation and functional RNA studies.

RESULTS

We show that near-splice positions and splicing branchpoints are highly constrained by purifying selection and harbour potentially damaging non-coding variants which are amenable to systematic analysis in sequencing data. From 258 de novo splicing variants in known rare disease genes, we identify 35 new likely diagnoses in probands with an unsolved rare disease. To date, we have confirmed a new diagnosis for six individuals, including four in whom RNA studies were performed.

CONCLUSIONS

Overall, we demonstrate the clinical value of examining non-canonical splicing variants in individuals with unsolved rare diseases.

Elevated Small Nuclear Ribonucleoprotein Polypeptide an Expression Correlated With Poor Prognosis and Immune Infiltrates in Patients With Hepatocellular Carcinoma

Background

Elevated Small Nuclear Ribonucleoprotein Polypeptide A (SNRPA) can enhance tumor cell growth and proliferation in various cancers. However, rarely studies focus on the comprehensive analysis of SNRPA in hepatocellular carcinoma (HCC).

Methods

TCGA and GEO databases were used to analyze the mRNA expression of SNRPA in HCC. Protein expression of SNAPA was validated using immunohistochemistry. Stably transfected HCC cells were used to investigate the role of SNRPA in the progression of HCC. The functional enrichment analysis was utilized for the biological function prediction. The CIBERSORT and ssGSEA algorithms were used to evaluate the composition of the tumor microenvironment and immunocyte infiltration ratio.

Results

The SNRPA expression was upregulated in HCC and positively correlated with tumor stage and grade. SNRPA overexpression were independent risk factors for poor overall survival (OS) and recurrence-free survival (RFS). In patients with early-stage disease, low alpha-fetoprotein expression, and better differentiation, SNRPA still exhibited the excellent prognostic value. Knockdown of SNRPA inhibited the proliferation and migration while promoting the apoptosis of HCC cells. Higher methylation of the CpG site cg16596691 correlated with longer OS in HCC patients. Genes co-expressed with SNRPA were overexpressed in HCC and correlated with shorter OS. The GO and KEGG enrichment analysis showed that SNRPA expression was related to mRNA splicing, spliceosome signaling. GSEA demonstrated that the main enrichment pathway of SNRPA-related differential genes was spliceosome signaling, cell cycle signaling, P53 signaling pathway, T cell receptor signaling pathway, natural killer cell-mediated signaling. CIBERSORT and ssGSEA algorithm revealed that SNRPA was mainly associated with the higher proportion of CD8+T cells, T cells follicular helper, T cells regulatory, Macrophages M0, and the lower proportion of T cells CD4 memory resting, NK cells resting, Monocytes, and Mast cells resting.

Conclusion

Elevated SNRPA enhances tumor cell proliferation and correlated with poor prognosis and immune infiltrates in patients with HCC.

A Pan-Cancer Analysis Revealing the Dual Roles of Lysine (K)-Specific Demethylase 6B in Tumorigenesis and Immunity

Introduction: Epigenetic-targeted therapy has been increasingly applied in the treatment of cancers. Lysine (K)-specific demethylase 6B (KDM6B) is an epigenetic enzyme involved in the coordinated control between cellular intrinsic regulators and the tissue microenvironment whereas the pan-cancer analysis of KDM6B remains unavailable.

Methods: The dual role of KDM6B in 33 cancers was investigated based on the GEO (Gene Expression Omnibus) and TCGA (The Cancer Genome Atlas) databases. TIMER2 and GEPIA2 were applied to investigate the KDM6B levels in different subtypes or stages of tumors. Besides, the Human Protein Atlas database allowed us to conduct a pan-cancer study of the KDM6B protein levels. GEPIA2 and Kaplan–Meier plotter were used for the prognosis analysis in different cancers. Characterization of genetic modifications of the KDM6B gene was analyzed by the cBioPortal. DNA methylation levels of different KDM6B probes in different TCGA tumors were analyzed by MEXPRESS. TIMER2 was applied to determine the association of the KDM6B expression and immune infiltration and DNA methyltransferases. Spearman correlation analysis was used to assess the association of the KDM6B expression with TMB (tumor mutation burden) and MSI (microsatellite instability). The KEGG (Kyoto encyclopedia of genes and genomes) pathway analysis and GO (Gene ontology) enrichment analysis were used to further investigate the potential mechanism of KDM6B in tumor pathophysiology.

Results: KDM6B was downregulated in 11 cancer types and upregulated across five types. In KIRC (kidney renal clear cell carcinoma) and OV (ovarian serous cystadenocarcinoma), the KDM6B level was significantly associated with the pathological stage. A high level of KDM6B was related to poor OS (overall survival) outcomes for THCA (thyroid carcinoma), while a low level was correlated with poor OS and DFS (disease-free survival) prognosis of KIRC. The KDM6B expression level was associated with TMB, MSI, and immune cell infiltration, particularly cancer-associated fibroblasts, across various cancer types with different correlations. Furthermore, the enrichment analysis revealed the relationship between H3K4 and H3K27 methylation and KDM6B function.

Conclusion: Dysregulation of the DNA methyltransferase activity and methylation levels of H3K4 and H3K27 may involve in the dual role of KDM6B in tumorigenesis and development. Our study offered a relatively comprehensive understanding of KDM6B’s dual role in cancer development and response to immunotherapy.

Mutations at a split codon in the GTPase-encoding domain of OPA1 cause dominant optic atrophy through different molecular mechanisms

Exonic (i.e. coding) variants in genes associated with disease can exert pathogenic effects both at the protein and mRNA level, either by altering the amino acid sequence or by affecting pre-mRNA splicing. The latter is often neglected due to the lack of RNA analyses in genetic diagnostic testing. In this study we considered both pathomechanisms and performed a comprehensive analysis of nine exonic nucleotide changes in OPA1, which is the major gene underlying autosomal dominant optic atrophy (DOA) and is characterized by pronounced allelic heterogeneity. We focused on the GTPase-encoding domain of OPA1, which harbors most of the missense variants associated with DOA. Given that the consensus splice sites extend into the exons, we chose a split codon, namely codon 438, for our analyses. Variants at this codon are the second most common cause of disease in our large cohort of DOA patients harboring disease-causing variants in OPA1. In silico splice predictions, heterologous splice assays, analysis of patient’s RNA when available, and protein modeling revealed different molecular outcomes for variants at codon 438. The wildtype aspartate residue at amino acid position 438 is directly involved in the dimerization of OPA1 monomers. We found that six amino acid substitutions at codon 438 (i.e. all substitutions of the first and second nucleotide of the codon) destabilized dimerization while only substitutions of the first nucleotide of the codon caused exon skipping. Our study highlights the value of combining RNA analysis and protein modeling approaches to accurately assign patients to future precision therapies.

Splicing in the Diagnosis of Rare Disease: Advances and Challenges

Mutations which affect splicing are significant contributors to rare disease, but are frequently overlooked by diagnostic sequencing pipelines. Greater ascertainment of pathogenic splicing variants will increase diagnostic yields, ending the diagnostic odyssey for patients and families affected by rare disorders, and improving treatment and care strategies. Advances in sequencing technologies, predictive modeling, and understanding of the mechanisms of splicing in recent years pave the way for improved detection and interpretation of splice affecting variants, yet several limitations still prohibit their routine ascertainment in diagnostic testing. This review explores some of these advances in the context of clinical application and discusses challenges to be overcome before these variants are comprehensively and routinely recognized in diagnostics.

Disrupting upstream translation in mRNAs is associated with human disease

Ribosome-profiling has uncovered pervasive translation in non-canonical open reading frames, however the biological significance of this phenomenon remains unclear. Using genetic variation from 71,702 human genomes, we assess patterns of selection in translated upstream open reading frames (uORFs) in 5’UTRs. We show that uORF variants introducing new stop codons, or strengthening existing stop codons, are under strong negative selection comparable to protein-coding missense variants. Using these variants, we map and validate gene-disease associations in two independent biobanks containing exome sequencing from 10,900 and 32,268 individuals, respectively, and elucidate their impact on protein expression in human cells. Our results suggest translation disrupting mechanisms relating uORF variation to reduced protein expression, and demonstrate that translation at uORFs is genetically constrained in 50% of human genes.

The significance of translated upstream open reading frames is not well known. Here, the authors investigate genetic variants in these regions, finding that they are under high evolutionary constraint and may contribute to disease.

Massively parallel assessment of human variants with base editor screens

Understanding the functional consequences of single-nucleotide variants is critical to uncovering the genetic underpinnings of diseases, but technologies to characterize variants are limiting. Here, we leverage CRISPR-Cas9 cytosine base editors in pooled screens to scalably assay variants at endogenous loci in mammalian cells. We benchmark the performance of base editors in positive and negative selection screens, identifying known loss-of-function mutations in BRCA1 and BRCA2 with high precision. To demonstrate the utility of base editor screens to probe small molecule-protein interactions, we screen against BH3 mimetics and PARP inhibitors, identifying point mutations that confer drug sensitivity or resistance. We also create a library of single guide RNAs (sgRNAs) predicted to generate 52,034 ClinVar variants in 3,584 genes and conduct screens in the presence of cellular stressors, identifying loss-of-function variants in numerous DNA damage repair genes. We anticipate that this screening approach will be broadly useful to readily and scalably functionalize genetic variants.

RET compound inheritance in Chinese patients with Hirschsprung disease: lack of penetrance from insufficient gene dysfunction

Hirschsprung disease (HSCR) is a neurocristopathy characterized by the absence of enteric ganglia along variable lengths of the intestine. Genetic defects play a major role in HSCR pathogenesis with nearly 50% of patients having a structural or regulatory deficiency in the major susceptibility gene RET. However, complete molecular defects remain poorly characterized in most patients. Here, we performed detailed genetic, molecular, and populational investigations of rare null mutations and modifiers at the RET locus. We first verified the pathogenicity of three RET splice site mutants (c.1879 + 1G > A, c.2607 + 5G > A and c.2608-3C > G) at the RNA level. We also identified significantly higher risk allele (genotype) frequencies, and their over-transmission, from unaffected parents to affected offspring of three functionally independent enhancer variants (rs2506030, rs7069590 and rs2435357, with odd ratios (OR) of 2.09, 2.71 and 7.59, respectively, P < 0.001). These three common variants are in significant (P < 4.64 × 10^-186) linkage disequilibrium in the Han Chinese population with ~ 60% of them carrying at least one copy and > 10% with two copies. We show that RET compound inheritance of rare and common variants prevails in 64% (seven out of 11) of Chinese HSCR families. This study supports the idea that common RET variants can modify the penetrance of rare null RET mutations in HSCR, and the combined high susceptibility allele dosage may constitute the unique raised "risk baseline" among the Chinese population.

Programmed Cell Death 2-Like (Pdcd2l) Is Required for Mouse Embryonic Development

Globozoospermia is a rare form of male infertility where men produce round-headed sperm that are incapable of fertilizing an oocyte naturally. In a previous study where we undertook a whole exome screen to define novel genetic causes of globozoospermia, we identified homozygous mutations in the gene PDCD2L. Two brothers carried a p.(Leu225Val) variant predicted to introduce a novel splice donor site, thus presenting PDCD2L as a potential regulator of male fertility. In this study, we generated a Pdcd2l knockout mouse to test its role in male fertility. Contrary to the phenotype predicted from its testis-enriched expression pattern, Pdcd2l null mice died during embryogenesis. Specifically, we identified that Pdcd2l is essential for post-implantation embryonic development. Pdcd2l^−/− embryos were resorbed at embryonic days 12.5-17.5 and no knockout pups were born, while adult heterozygous Pdcd2l males had comparable fertility to wildtype males. To specifically investigate the role of PDCD2L in germ cells, we employed Drosophila melanogaster as a model system. Consistent with the mouse data, global knockdown of trus, the fly ortholog of PDCD2L, resulted in lethality in flies at the third instar larval stage. However, germ cell-specific knockdown with two germ cell drivers did not affect male fertility. Collectively, these data suggest that PDCD2L is not essential for male fertility. By contrast, our results demonstrate an evolutionarily conserved role of PDCD2L in development.

Pathogenic Variants in the Myosin Chaperone UNC-45B Cause Progressive Myopathy with Eccentric Cores

The myosin-directed chaperone UNC-45B is essential for sarcomeric organization and muscle function from Caenorhabditis elegans to humans. The pathological impact of UNC-45B in muscle disease remained elusive. We report ten individuals with bi-allelic variants in UNC45B who exhibit childhood-onset progressive muscle weakness. We identified a common UNC45B variant that acts as a complex hypomorph splice variant. Purified UNC-45B mutants showed changes in folding and solubility. In situ localization studies further demonstrated reduced expression of mutant UNC-45B in muscle combined with abnormal localization away from the A-band towards the Z-disk of the sarcomere. The physiological relevance of these observations was investigated in C. elegans by transgenic expression of conserved UNC-45 missense variants, which showed impaired myosin binding for one and defective muscle function for three. Together, our results demonstrate that UNC-45B impairment manifests as a chaperonopathy with progressive muscle pathology, which discovers the previously unknown conserved role of UNC-45B in myofibrillar organization.

Mutations primarily alter the inclusion of alternatively spliced exons

Genetic analyses and systematic mutagenesis have revealed that synonymous, non-synonymous and intronic mutations frequently alter the inclusion levels of alternatively spliced exons, consistent with the concept that altered splicing might be a common mechanism by which mutations cause disease. However, most exons expressed in any cell are highly-included in mature mRNAs. Here, by performing deep mutagenesis of highly-included exons and by analysing the association between genome sequence variation and exon inclusion across the transcriptome, we report that mutations only very rarely alter the inclusion of highly-included exons. This is true for both exonic and intronic mutations as well as for perturbations in trans. Therefore, mutations that affect splicing are not evenly distributed across primary transcripts but are focussed in and around alternatively spliced exons with intermediate inclusion levels. These results provide a resource for prioritising synonymous and other variants as disease-causing mutations.

Evaluation of both exonic and intronic variants for effects on RNA splicing allows for accurate assessment of the effectiveness of precision therapies

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.

Genetic variants that impact pre-mRNA splicing are a common cause of genetic disease and have varying downstream molecular consequences. As a result, precision therapies that function at the protein level are not always effective for these variants and thus careful assessment is necessary. Here we evaluate RNA-level effects of 52 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene and show that study of splicing and its consequences allows for more accurate assignment of precision therapies.

Transcriptomic signatures across human tissues identify functional rare genetic variation

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

Recessive MYH3 variants cause "Contractures, pterygia, and variable skeletal fusions syndrome 1B" mimicking Escobar variant multiple pterygium syndrome

The multiple pterygium syndromes (MPS) are rare disorders with disease severity ranging from lethal to milder forms. The nonlethal Escobar variant MPS (EVMPS) is characterized by multiple pterygia and arthrogryposis, as well as various additional features including congenital anomalies. The genetic etiology of EVMPS is heterogeneous and the diagnosis has been based either on the detection of pathogenic CHRNG variants (~23% of patients), or suggestive clinical features. We describe four patients with a clinical suspicion of EVMPS who manifested with multiple pterygia, mild flexion contractures of several joints, and vertebral anomalies. We revealed recessively inherited MYH3 variants as the underlying cause in all patients: two novel variants, c.1053C>G, p.(Tyr351Ter) and c.3102+5G>C, as compound heterozygous with the hypomorphic MYH3 variant c.-9+1G>A. Recessive MYH3 variants have been previously associated with spondylocarpotarsal synostosis syndrome. Our findings now highlight multiple pterygia as an important feature in patients with recessive MYH3 variants. Based on all patients with recessive MYH3 variants reported up to date, we consider that this disease entity should be designated as "Contractures, pterygia, and variable skeletal fusions syndrome 1B," as recently suggested by OMIM. Our findings underline the importance of analyzing MYH3 in the differential diagnosis of EVMPS, particularly as the hypomorphic MYH3 variant might remain undetected by routine exome sequencing.

Pathogenic deep intronic MTM1 variant activates a pseudo-exon encoding a nonsense codon resulting in severe X-linked myotubular myopathy

X-linked myotubular myopathy (XLMTM) is a severe congenital myopathy characterised by generalised weakness and respiratory insufficiency. XLMTM is associated with pathogenic variants in MTM1; a gene encoding the lipid phosphatase myotubularin. Whole genome sequencing (WGS) of an exome-negative male proband with severe hypotonia, respiratory insufficiency and centralised nuclei on muscle biopsy identified a deep intronic MTM1 variant NG_008199.1(NM_000252.2):c.1468-577A>G, which strengthened a cryptic 5' splice site (A>G substitution at the +5 position). Muscle RNA sequencing was non-diagnostic due to low read depth. Reverse transcription PCR (RT-PCR) of muscle RNA confirmed the c.1468-577A>G variant activates inclusion of a pseudo-exon encoding a premature stop codon into all detected MTM1 transcripts. Western blot analysis establishes deficiency of myotubularin protein, consistent with the severe XLMTM phenotype. We expand the genotypic spectrum of XLMTM and highlight benefits of screening non-coding regions of MTM1 in male probands with phenotypically concordant XLMTM who remain undiagnosed following exome sequencing.

Genome-wide association study of cerebral small vessel disease reveals established and novel loci

Intracerebral haemorrhage and small vessel ischaemic stroke (SVS) are the most acute manifestations of cerebral small vessel disease, with no established preventive approaches beyond hypertension management. Combined genome-wide association study (GWAS) of these two correlated diseases may improve statistical power to detect novel genetic factors for cerebral small vessel disease, elucidating underlying disease mechanisms that may form the basis for future treatments. Because intracerebral haemorrhage location is an adequate surrogate for distinct histopathological variants of cerebral small vessel disease (lobar for cerebral amyloid angiopathy and non-lobar for arteriolosclerosis), we performed GWAS of intracerebral haemorrhage by location in 1813 subjects (755 lobar and 1005 non-lobar) and 1711 stroke-free control subjects. Intracerebral haemorrhage GWAS results by location were meta-analysed with GWAS results for SVS from MEGASTROKE, using 'Multi-Trait Analysis of GWAS' (MTAG) to integrate summary data across traits and generate combined effect estimates. After combining intracerebral haemorrhage and SVS datasets, our sample size included 241 024 participants (6255 intracerebral haemorrhage or SVS cases and 233 058 control subjects). Genome-wide significant associations were observed for non-lobar intracerebral haemorrhage enhanced by SVS with rs2758605 [MTAG P-value (P) = 2.6 × 10-8] at 1q22; rs72932727 (P = 1.7 × 10-8) at 2q33; and rs9515201 (P = 5.3 × 10-10) at 13q34. In the GTEx gene expression library, rs2758605 (1q22), rs72932727 (2q33) and rs9515201 (13q34) are significant cis-eQTLs for PMF1 (P = 1 × 10-4 in tibial nerve), NBEAL1, FAM117B and CARF (P < 2.1 × 10-7 in arteries) and COL4A2 and COL4A1 (P < 0.01 in brain putamen), respectively. Leveraging S-PrediXcan for gene-based association testing with the predicted expression models in tissues related with nerve, artery, and non-lobar brain, we found that experiment-wide significant (P < 8.5 × 10-7) associations at three genes at 2q33 including NBEAL1, FAM117B and WDR12 and genome-wide significant associations at two genes including ICA1L at 2q33 and ZCCHC14 at 16q24. Brain cell-type specific expression profiling libraries reveal that SEMA4A, SLC25A44 and PMF1 at 1q22 and COL4A1 and COL4A2 at 13q34 were mainly expressed in endothelial cells, while the genes at 2q33 (FAM117B, CARF and NBEAL1) were expressed in various cell types including astrocytes, oligodendrocytes and neurons. Our cross-phenotype genetic study of intracerebral haemorrhage and SVS demonstrates novel genome-wide associations for non-lobar intracerebral haemorrhage at 2q33 and 13q34. Our replication of the 1q22 locus previous seen in traditional GWAS of intracerebral haemorrhage, as well as the rediscovery of 13q34, which had previously been reported in candidate gene studies with other cerebral small vessel disease-related traits strengthens the credibility of applying this novel genome-wide approach across intracerebral haemorrhage and SVS.

Characterising the loss-of-function impact of 5' untranslated region variants in 15,708 individuals

Upstream open reading frames (uORFs) are tissue-specific cis-regulators of protein translation. Isolated reports have shown that variants that create or disrupt uORFs can cause disease. Here, in a systematic genome-wide study using 15,708 whole genome sequences, we show that variants that create new upstream start codons, and variants disrupting stop sites of existing uORFs, are under strong negative selection. This selection signal is significantly stronger for variants arising upstream of genes intolerant to loss-of-function variants. Furthermore, variants creating uORFs that overlap the coding sequence show signals of selection equivalent to coding missense variants. Finally, we identify specific genes where modification of uORFs likely represents an important disease mechanism, and report a novel uORF frameshift variant upstream of NF2 in neurofibromatosis. Our results highlight uORF-perturbing variants as an under-recognised functional class that contribute to penetrant human disease, and demonstrate the power of large-scale population sequencing data in studying non-coding variant classes.

Clinical impact of splicing in neurodevelopmental disorders

Clinical exome sequencing is frequently used to identify gene-disrupting variants in individuals with neurodevelopmental disorders. While splice-disrupting variants are known to contribute to these disorders, clinical interpretation of cryptic splice variants outside of the canonical splice site has been challenging. Here, we discuss papers that improve such detection.

RNA Splicing Defects in Hypertrophic Cardiomyopathy: Implications for Diagnosis and Therapy

Hypertrophic cardiomyopathy (HCM), the most common inherited heart disease, is predominantly caused by mutations in genes that encode sarcomere-associated proteins. Effective gene-based diagnosis is critical for the accurate clinical management of patients and their family members. However, the introduction of high-throughput DNA sequencing approaches for clinical diagnostics has vastly expanded the number of variants of uncertain significance, leading to many inconclusive results that limit the clinical utility of genetic testing. More recently, developments in RNA analysis have been improving diagnostic outcomes by identifying new variants that interfere with splicing. This review summarizes recent discoveries of RNA mis-splicing in HCM and provides an overview of research that aims to apply the concept of RNA therapeutics to HCM.

One byte at a time: evidencing the quality of clinical service next-generation sequencing for germline and somatic variants

Next-generation sequencing (NGS) is replacing other molecular techniques to become the de facto gene diagnostics approach, transforming the speed of diagnosis for patients and expanding opportunities for precision medicine. Consequently, for accredited laboratories as well as those seeking accreditation, both objective measures of quality and external review of laboratory processes are required. External quality assessment (EQA), or Proficiency Testing (PT), can assess a laboratory’s service through an independent external agency, the EQA provider. The analysis of a growing number of genes and whole exome and genomes is now routine; therefore, an EQA must be delivered to enable all testing laboratories to participate. In this paper, we describe the development of a unique platform and gene target independent EQA scheme for NGS, designed to scale from current to future requirements of clinical diagnostic laboratories testing for germline and somatic variants. The EQA results from three annual rounds indicate that clinical diagnostic laboratories are providing an increasingly high-quality NGS service and variant calling abilities are improving. From an EQA provider perspective, challenges remain regarding delivery and performance criteria, as well as in analysing similar NGS approaches between cohorts with meaningful metrics, sample sourcing and data formats.

Common pathogenic mechanism in patients with dropped head syndrome caused by different mutations in the MYH7 gene

Mutations in the MYH7 gene are the source of an allelic series of diseases, including various cardiomyopathies and skeletal myopathies that usually manifest in adulthood. We observed a 1.5 y.o. male patient with congenital weaknesses of the axial muscles, "dropped head" syndrome, and dilated cardiomyopathy. The clinical evaluation included medical history, an echocardiogram, electromyography, and a histopathological study. The genetic evaluation included whole exome sequencing. Muscle biopsy samples from the proband were used for mRNA extraction. We revealed a novel genetic variant c.5655 + 5G > C in the MYH7 gene. The analysis of the cDNA showed an in-frame skipping of exon 38 (p.1854_1885del). This variant and two previously published mutations (c.5655G > A and c.5655 + 1G > A), also presumably leading to exon 38 skipping, were studied by expression analysis in the HEK293T cell line transfected with 4 plasmids containing the MYH7 minigene (wt, c.5655G > C, c.5655 + 1G > A and c.5655 + 5G > A). A quantitative difference in expression was shown for cell lines with each of the three mutant plasmids. All mutation carriers had a similar phenotype and included congenital axial myopathy and variable cardiac involvement. Prominent dropped head syndrome was mentioned in all patients. Early-onset axial myopathy with a dropped head syndrome is a distinct clinical entity within MYH7-related disorders. We suggest that mutations in the MYH7 gene affecting the C-terminal domain of beta-myosin heavy chain should also be considered as a possible cause in cases of early-onset myopathy with "dropped head" syndrome.