Inborn Errors of RNA Lariat Metabolism in Humans with Brainstem Viral Infection.

Viruses that are typically benign sometimes invade the brainstem in otherwise healthy children. We report bi-allelic DBR1 mutations in unrelated patients from different ethnicities, each of whom had brainstem infection due to herpes simplex virus 1 (HSV1), influenza virus, or norovirus. DBR1 encodes the only known RNA lariat debranching enzyme. We show that DBR1 expression is ubiquitous, but strongest in the spinal cord and brainstem. We also show that all DBR1 mutant alleles are severely hypomorphic, in terms of expression and function. The fibroblasts of DBR1-mutated patients contain higher RNA lariat levels than control cells, this difference becoming even more marked during HSV1 infection. Finally, we show that the patients' fibroblasts are highly susceptible to HSV1. RNA lariat accumulation and viral susceptibility are rescued by wild-type DBR1. Autosomal recessive, partial DBR1 deficiency underlies viral infection of the brainstem in humans through the disruption of tissue-specific and cell-intrinsic immunity to viruses.

A functional link between lariat debranching enzyme and the intron-binding complex is defective in non-photosensitive trichothiodystrophy.

The pre-mRNA life cycle requires intron processing; yet, how intron-processing defects influence splicing and gene expression is unclear. Here, we find that TTDN1/MPLKIP, which is encoded by a gene implicated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to spliceosomal function. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, whereas its N-terminal intrinsically disordered region (IDR) binds the intron-binding complex (IBC). TTDN1 loss, or a mutated IDR, causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. A Ttdn1-deficient mouse model recapitulates intron-processing defects and certain neurodevelopmental phenotypes seen in NP-TTD. Fusing DBR1 to the TTDN1 IDR is sufficient to recruit DBR1 to the IBC and circumvents the functional requirement for TTDN1. Collectively, our findings link RNA lariat processing with splicing outcomes by revealing the molecular function of TTDN1.

Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans.

Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.

Massively parallel reporter assays discover de novo exonic splicing mutants in paralogs of Autism genes.

To determine the contribution of defective splicing in Autism Spectrum Disorders (ASD), the most common neurodevelopmental disorder, a high throughput Massively Parallel Splicing Assay (MaPSY) was employed and identified 42 exonic splicing mutants out of 725 coding de novo variants discovered in the sequencing of ASD families. A redesign of the minigene constructs in MaPSY revealed that upstream exons with strong 5' splice sites increase the magnitude of skipping phenotypes observed in downstream exons. Select hits were validated by RT-PCR and amplicon sequencing in patient cell lines. Exonic splicing mutants were enriched in probands relative to unaffected siblings -especially synonymous variants (7.5% vs 3.5%, respectively). Of the 26 genes disrupted by exonic splicing mutations, 6 were in known ASD genes and 3 were in paralogs of known ASD genes. Of particular interest was a synonymous variant in TNRC6C - an ASD gene paralog with interactions with other ASD genes. Clinical records of 3 ASD patients with TNRC6C variant revealed respiratory issues consistent with phenotypes observed in TNRC6 depleted mice. Overall, this study highlights the need for splicing analysis in determining variant pathogenicity, especially as it relates to ASD.

Activation of human RNA lariat debranching enzyme Dbr1 by binding protein TTDN1 occurs though an intrinsically disordered C-terminal domain.

In eukaryotic cells, the introns are excised from pre-mRNA by the spliceosome. These introns typically have a lariat configuration due to the 2'-5' phosphodiester bond between an internal branched residue and the 5' terminus of the RNA. The only enzyme known to selectively hydrolyze the 2'-5' linkage of these lariats is the RNA lariat debranching enzyme Dbr1. In humans, Dbr1 is involved in processes such as class-switch recombination of immunoglobulin genes, and its dysfunction is implicated in viral encephalitis, HIV, ALS, and cancer. However, mechanistic details of precisely how Dbr1 affects these processes are missing. Here we show that human Dbr1 contains a disordered C-terminal domain through sequence analysis and nuclear magnetic resonance. This domain stabilizes Dbr1 in vitro by reducing aggregation but is dispensable for debranching activity. We establish that Dbr1 requires Fe2+ for efficient catalysis and demonstrate that the noncatalytic protein Drn1 and the uncharacterized protein trichothiodystrophy nonphotosensitive 1 directly bind to Dbr1. We demonstrate addition of trichothiodystrophy nonphotosensitive 1 to in vitro debranching reactions increases the catalytic efficiency of human Dbr1 19-fold but has no effect on the activity of Dbr1 from the amoeba Entamoeba histolytica, which lacks a disordered C-terminal domain. Finally, we systematically examine how the identity of the branchpoint nucleotide affects debranching rates. These findings describe new aspects of Dbr1 function in humans and further clarify how Dbr1 contributes to human health and disease.

Metal content and kinetic properties of yeast RNA lariat debranching enzyme Dbr1.

In eukaryotic cells, intron lariats produced by the spliceosome contain a 2'5' phosphodiester linkage. The RNA lariat debranching enzyme, Dbr1, is the only enzyme known to hydrolyze this bond. Dbr1 is a member of the metallophosphoesterase (MPE) family of enzymes, and recent X-ray crystal structures and biochemistry data demonstrate that Dbr1 from Entamoeba histolytica uses combinations of Mn2+, Zn2+, and Fe2+ as enzymatic cofactors. Here, we examine the kinetic properties and metal dependence of the Dbr1 homolog from Saccharomyces cerevisiae (yDbr1). Elemental analysis measured stoichiometric quantities of Fe and Zn in yDbr1 purified following heterologous expression E. coli We analyzed the ability of Fe2+, Zn2+, and Mn2+ to reconstitute activity in metal-free apoenzyme. Purified yDbr1 was highly active, turning over substrate at 5.6 sec-1, and apo-yDbr1 reconstituted with Fe2+ was the most active species, turning over at 9.2 sec-1 We treated human lymphoblastoid cells with the iron-chelator deferoxamine and measured a twofold increase in cellular lariats. These data suggest that Fe is an important biological cofactor for Dbr1 enzymes.

Crystal Structure of the RNA Lariat Debranching Enzyme Dbr1 with Hydrolyzed Phosphorothioate RNA Product.

The RNA lariat debranching enzyme is the sole enzyme responsible for hydrolyzing the 2'-5' phosphodiester bond in RNA lariats produced by the spliceosome. Here, we test the ability of Dbr1 to hydrolyze branched RNAs (bRNAs) that contain a 2'-5'-phosphorothioate linkage, a modification commonly used to resist degradation. We attempted to cocrystallize a phosphorothioate-branched RNA (PS-bRNA) with wild-type Entamoeba histolytica Dbr1 (EhDbr1) but observed in-crystal hydrolysis of the phosphorothioate bond. The crystal structure revealed EhDbr1 in a product-bound state, with the hydrolyzed 2'-5' fragment of the PS-bRNA mimicking the binding mode of the native bRNA substrate. These findings suggest that product inhibition may contribute to the kinetic mechanism of Dbr1. We show that Dbr1 enzymes cleave phosphorothioate linkages at rates ∼10,000-fold more slowly than native phosphate linkages. This new product-bound crystal structure offers atomic details, which can aid inhibitor design. Dbr1 inhibitors could be therapeutic or investigative compounds for human diseases such as human immunodeficiency virus (HIV), amyotrophic lateral sclerosis (ALS), cancer, and viral encephalitis.

Hereditary cancer genes are highly susceptible to splicing mutations.

Substitutions that disrupt pre-mRNA splicing are a common cause of genetic disease. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations mapping to the canonical 5' and 3' splice sites. However, splicing mutations present in exons and deeper intronic positions are vastly underreported. A recent re-analysis of coding mutations in exon 10 of the Lynch Syndrome gene, MLH1, revealed an extremely high rate (77%) of mutations that lead to defective splicing. This finding is confirmed by extending the sampling to five other exons in the MLH1 gene. Further analysis suggests a more general phenomenon of defective splicing driving Lynch Syndrome. Of the 36 mutations tested, 11 disrupted splicing. Furthermore, analyzing past reports suggest that MLH1 mutations in canonical splice sites also occupy a much higher fraction (36%) of total mutations than expected. When performing a comprehensive analysis of splicing mutations in human disease genes, we found that three main causal genes of Lynch Syndrome, MLH1, MSH2, and PMS2, belonged to a class of 86 disease genes which are enriched for splicing mutations. Other cancer genes were also enriched in the 86 susceptible genes. The enrichment of splicing mutations in hereditary cancers strongly argues for additional priority in interpreting clinical sequencing data in relation to cancer and splicing.

Large-scale analysis of branchpoint usage across species and cell lines.

The coding sequence of each human pre-mRNA is interrupted, on average, by 11 introns that must be spliced out for proper gene expression. Each intron contains three obligate signals: a 5' splice site, a branch site, and a 3' splice site. Splice site usage has been mapped exhaustively across different species, cell types, and cellular states. In contrast, only a small fraction of branch sites have been identified even once. The few reported annotations of branch site are imprecise as reverse transcriptase skips several nucleotides while traversing a 2-5 linkage. Here, we report large-scale mapping of the branchpoints from deep sequencing data in three different species and in the SF3B1 K700E oncogenic mutant background. We have developed a novel method whereby raw lariat reads are refined by U2snRNP/pre-mRNA base-pairing models to return the largest current data set of branchpoint sequences with quality metrics. This analysis discovers novel modes of U2snRNA:pre-mRNA base-pairing conserved in yeast and provides insight into the biogenesis of intron circles. Finally, matching branch site usage with isoform selection across the extensive panel of ENCODE RNA-seq data sets offers insight into the mechanisms by which branchpoint usage drives alternative splicing.

Future directions for high-throughput splicing assays in precision medicine.

Classification of variants of unknown significance is a challenging technical problem in clinical genetics. As up to one-third of disease-causing mutations are thought to affect pre-mRNA splicing, it is important to accurately classify splicing mutations in patient sequencing data. Several consortia and healthcare systems have conducted large-scale patient sequencing studies, which discover novel variants faster than they can be classified. Here, we compare the advantages and limitations of several high-throughput splicing assays aimed at mitigating this bottleneck, and describe a data set of ~5,000 variants that we analyzed using our Massively Parallel Splicing Assay (MaPSy). The Critical Assessment of Genome Interpretation group (CAGI) organized a challenge, in which participants submitted machine learning models to predict the splicing effects of variants in this data set. We discuss the winning submission of the challenge (MMSplice) which outperformed existing software. Finally, we highlight methods to overcome the limitations of MaPSy and similar assays, such as tissue-specific splicing, the effect of surrounding sequence context, classifying intronic variants, synthesizing large exons, and amplifying complex libraries of minigene species. Further development of these assays will greatly benefit the field of clinical genetics, which lack high-throughput methods for variant interpretation.

Assessing predictions of the impact of variants on splicing in CAGI5.

Precision medicine and sequence-based clinical diagnostics seek to predict disease risk or to identify causative variants from sequencing data. The Critical Assessment of Genome Interpretation (CAGI) is a community experiment consisting of genotype-phenotype prediction challenges; participants build models, undergo assessment, and share key findings. In the past, few CAGI challenges have addressed the impact of sequence variants on splicing. In CAGI5, two challenges (Vex-seq and MaPSY) involved prediction of the effect of variants, primarily single-nucleotide changes, on splicing. Although there are significant differences between these two challenges, both involved prediction of results from high-throughput exon inclusion assays. Here, we discuss the methods used to predict the impact of these variants on splicing, their performance, strengths, and weaknesses, and prospects for predicting the impact of sequence variation on splicing and disease phenotypes.

RNA structure replaces the need for U2AF2 in splicing.

RNA secondary structure plays an integral role in catalytic, ribosomal, small nuclear, micro, and transfer RNAs. Discovering a prevalent role for secondary structure in pre-mRNAs has proven more elusive. By utilizing a variety of computational and biochemical approaches, we present evidence for a class of nuclear introns that relies upon secondary structure for correct splicing. These introns are defined by simple repeat expansions of complementary AC and GT dimers that co-occur at opposite boundaries of an intron to form a bridging structure that enforces correct splice site pairing. Remarkably, this class of introns does not require U2AF2, a core component of the spliceosome, for its processing. Phylogenetic analysis suggests that this mechanism was present in the ancestral vertebrate lineage prior to the divergence of tetrapods from teleosts. While largely lost from land dwelling vertebrates, this class of introns is found in 10% of all zebrafish genes.

Widespread intra-dependencies in the removal of introns from human transcripts.

Research into the problem of splice site selection has followed a reductionist approach focused on how individual splice sites are recognized. Early applications of information theory uncovered an inconsistency. Human splice signals do not contain enough information to explain the observed fidelity of splicing. Here, we conclude that introns do not necessarily contain 'missing' information but rather may require definition from neighboring processing events. For example, there are known cases where an intronic mutation disrupts the splicing of not only the local intron but also adjacent introns. We present a genome-wide measurement of the order of splicing within human transcripts. The observed order of splicing cannot be explained by a simple kinetic model. Simulations reveal a bias toward a particular, transcript-specific order of intron removal in human genes. We validate an extreme class of intron that can only splice in a multi-intron context. Special categories of splicing such as exon circularization, first and last intron processing, alternative 5 and 3'ss usage and exon skipping are marked by distinct patterns of ordered intron removal. Excessive intronic length and silencer density tend to delay splicing. Shorter introns that contain enhancers splice early.

Changes in the process of alternative RNA splicing results in soluble B and T lymphocyte attenuator with biological and clinical implications in critical illness.

BACKGROUND: Critically ill patients with sepsis and acute respiratory distress syndrome have severely altered physiology and immune system modifications. RNA splicing is a basic molecular mechanism influenced by physiologic alterations. Immune checkpoint inhibitors, such as B and T Lymphocyte Attenuator (BTLA) have previously been shown to influence outcomes in critical illness. We hypothesize altered physiology in critical illness results in alternative RNA splicing of the immune checkpoint protein, BTLA, resulting in a soluble form with biologic and clinical significance. METHODS: Samples were collected from critically ill humans and mice. Levels soluble BTLA (sBTLA) were measured. Ex vivo experiments assessing for cellular proliferation and cytokine production were done using splenocytes from critically ill mice cultured with sBTLA. Deep RNA sequencing was done to look for alternative splicing of BTLA. sBTLA levels were fitted to models to predict sepsis diagnosis. RESULTS: sBTLA is increased in the blood of critically ill humans and mice and can predict a sepsis diagnosis on hospital day 0 in humans. Alternative RNA splicing results in a premature stop codon that results in the soluble form. sBTLA has a clinically relevant impact as splenocytes from mice with critical illness cultured with soluble BTLA have increased cellular proliferation. CONCLUSION: sBTLA is produced as a result of alternative RNA splicing. This isoform of BTLA has biological significance through changes in cellular proliferation and can predict the diagnosis of sepsis.

The effects of structure on pre-mRNA processing and stability.

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

Defective splicing of the RB1 transcript is the dominant cause of retinoblastomas.

Defective splicing is a common cause of genetic diseases. On average, 13.4% of all hereditary disease alleles are classified as splicing mutations with most mapping to the critical GT or AG nucleotides within the 5' and 3' splice sites. However, splicing mutations are underreported and the fraction of splicing mutations that compose all disease alleles varies greatly across disease gene. For example, there is a great excess (46%; ~threefold) of hereditary disease alleles that map to splice sites in RB1 that cause retinoblastoma. Furthermore, mutations in the exons and deeper intronic position may also affect splicing. We recently developed a high-throughput method that assays reported disease mutations for their ability to disrupt pre-mRNA splicing. Surprisingly, 27% of RB1-coding mutations tested also disrupt splicing. High-throughput in vitro spliceosomal assembly assay reveals heterogeneity in which stage of spliceosomal assembly is affected by splicing mutations. 58% of exonic splicing mutations were primarily blocked at the A complex in transition to the B complex and 33% were blocked at the B complex. Several mutants appear to reduce more than one step in the assembly. As RB1 splicing mutants are enriched in retinoblastoma disease alleles, additional priority should be allocated to this class of allele while interpreting clinical sequencing experiments. Analysis of the spectrum of RB1 variants observed in 60,706 exomes identifies 197 variants that have enough potential to disrupt splicing to warrant further consideration.

A general mechanism for transcription regulation by Oct1 and Oct4 in response to genotoxic and oxidative stress.

Oct1 and Oct4 are homologous transcription factors with similar DNA-binding specificities. Here we show that Oct1 is dynamically phosphorylated in vivo following exposure of cells to oxidative and genotoxic stress. We further show that stress regulates the selectivity of both proteins for specific DNA sequences. Mutation of conserved phosphorylation target DNA-binding domain residues in Oct1, and Oct4 confirms their role in regulating binding selectivity. Using chromatin immunoprecipitation, we show that association of Oct4 and Oct1 with a distinct group of in vivo targets is inducible by stress, and that Oct1 is essential for a normal post-stress transcriptional response. Finally, using an unbiased Oct1 target screen we identify a large number of genes targeted by Oct1 specifically under conditions of stress, and show that several of these inducible Oct1 targets are also inducibly bound by Oct4 in embryonic stem cells following stress exposure.

Soluble programmed cell death receptor-1 (sPD-1): a potential biomarker with anti-inflammatory properties in human and experimental acute respiratory distress syndrome (ARDS).

BACKGROUND: Acute respiratory distress syndrome (ARDS) remains a common organ dysfunction in the critically ill patient. Mechanisms for its development have focused on immune mediated causes, aspects of our understanding are not complete, and we lack biomarkers. DESIGN, SETTING, AND SUBJECTS: Blood and bronchial alveolar lavage fluid (BAL) from humans (n = 10-13) with ARDS and controls (n = 5-10) as well as a murine model of ARDS (n = 5-6) with controls (n = 6-7) were studied. METHODS: ARDS was induced in mice by hemorrhagic shock (day 1) followed by poly-microbial sepsis (day 2). Samples were then collected on the third day after the animals were euthanized. Ex vivo experiments used splenocytes from animals with ARDS cultured with and without soluble programmed death receptor-1 (sPD-1). RESULTS: Levels of sPD-1 are increased in both the serum (11,429.3 pg/mL(SD 2133.3) vs. 8061.4(SD 4187.8), p = 0.036) and bronchial alveolar lavage (BAL) fluid (6,311.1 pg/mL(SD 3758.0) vs. 90.7 pg/mL(SD 202.8), p = 0.002) of humans with ARDS. Similar results are seen in the serum (9396.1 pg/mL(SD 1546.0) vs. 3464.5 pg/mL(SD 2511.8), p = 0.001) and BAL fluid (2891.7 pg/mL(SD 868.1) vs. 1385.9 pg/mL(SD 927.8), p = 0.012) of mice. sPD-1 levels in murine blood (AUC = 1(1-1), p = 0.006), murine BAL fluid (AUC = 0.905(0.717-1.093), p = 0.015), and human BAL (AUC = 1(1-1), p = 0.001) fluid predicted ARDS. To assess the importance of sPD-1 in ARDS, ex vivo experiments were undertaken. BAL fluid from mice with ARDS dampens the TNF-α production compared to cells cultured with BAL lacking sPD-1 (2.7 pg/mL(SD 3.8) vs. 52.38 pg/mL(SD 25.1), p = 0.002). CONCLUSIONS: This suggests sPD-1 is elevated in critical illness and may represent a potential biomarker for ARDS. In addition, sPD-1 has an anti-inflammatory mechanism in conditions of marked stress and aids in the resolution of severe inflammation. sPD-1 could be used to not only diagnose ARDS, but may be a potential therapy.

RNA structure in splicing: An evolutionary perspective.

Pre-mRNA splicing is a key post-transcriptional regulation process in which introns are excised and exons are ligated together. A novel class of structured intron was recently discovered in fish. Simple expansions of complementary AC and GT dimers at opposite boundaries of an intron were found to form a bridging structure, thereby enforcing correct splice site pairing across the intron. In some fish introns, the RNA structures are strong enough to bypass the need of regulatory protein factors for splicing. Here, we discuss the prevalence and potential functions of highly structured introns. In humans, structured introns usually arise through the co-occurrence of C and G-rich repeats at intron boundaries. We explore the potentially instructive example of the HLA receptor genes. In HLA pre-mRNA, structured introns flank the exons that encode the highly polymorphic ß sheet cleft, making the processing of the transcript robust to variants that disrupt splicing factor binding. While selective forces that have shaped HLA receptor are fairly atypical, numerous other highly polymorphic genes that encode receptors contain structured introns. Finally, we discuss how the elevated mutation rate associated with the simple repeats that often compose structured intron can make structured introns themselves rapidly evolving elements.

ISWI contributes to ArsI insulator function in development of the sea urchin.

Insulators are genomic elements that regulate transcriptional activity by forming chromatin boundaries. Various DNA insulators have been identified or are postulated in many organisms, and the paradigmatic CTCF-dependent insulators are perhaps the best understood and most widespread in function. The diversity of DNA insulators is, however, understudied, especially in the context of embryonic development, when many new gene territories undergo transitions in functionality. Here we report the functional analysis of the arylsulfatase insulator (ArsI) derived from the sea urchin, which has conserved insulator activity throughout the many metazoans tested, but for which the molecular mechanism of function is unknown. Using a rapid in vivo assay system and a high-throughput mega-shift assay, we identified a minimal region in ArsI that is responsible for its insulator function. We discovered a small set of proteins specifically bound to the minimal ArsI region, including ISWI, a known chromatin-remodeling protein. During embryogenesis, ISWI was found to interact with select ArsI sites throughout the genome, and when inactivated led to misregulation of select gene expression, loss of insulator activity and aberrant morphogenesis. These studies reveal a mechanistic basis for ArsI function in the gene regulatory network of early development.