Recent studies report the genetic loss of the lariat debranching enzyme ( DBR1 ) activity increases susceptibility to viral infection. Here, we show that more than 25% of human introns contain large hairpin structures created by the folding of two Alu elements inserted in opposite orientation. In wildtype cells, this large reservoir of endogenous dsRNA is efficiently degraded. In DBR1 -null cells, lariats accumulate in the cytosol and dsRNA becomes enriched. We demonstrate how the chronic exposure to these lariats attenuates the dsRNA sensors, reducing the response of the MDA5, RIG-I, RNase L and PKR sensing pathways. We observe evidence for both attenuation and endogenous dsRNA in anti-viral response and viral evasion. Lariats are transiently elevated during infection (e.g. HSV-1, influenza, KSHV). The HSV-1 genome expresses multiple, stable lariats that may attenuate dsRNA sensors during latency. HIGHLIGHTS:Intronic inverted repeat Alu elements constitute largest source of endogenous dsRNA. In the absence of DBR1 , lariats accumulate in the cytoplasm and form dsRNA. Chronic exposure to endogenous dsRNA in a DBR1 -depleted environment desensitizes the dsRNA sensing pathway. ICP0 intron 1 of HSV-1 has a highly-structured stable lariat. GRAPHICAL ABSTRACT:

The molecular mechanism by which inborn errors of the human RNA lariat-debranching enzyme 1 (DBR1) underlie brainstem viral encephalitis is unknown. We show here that the accumulation of RNA lariats in human DBR1-deficient cells interferes with stress granule (SG) assembly, promoting the proteasome degradation of at least G3BP1 and G3BP2, two key components of SGs. In turn, impaired assembly of SGs, which normally recruit PKR, impairs PKR activation and activity against viruses, including HSV-1. Remarkably, the genetic ablation of PKR abolishes the corresponding antiviral effect of DBR1 in vitro. We also show that Dbr1Y17H/Y17H mice are susceptible to similar viral infections in vivo. Moreover, cells and brain samples from Dbr1Y17H/Y17H mice exhibit decreased G3BP1/2 expression and PKR phosphorylation. Thus, the debranching of RNA lariats by DBR1 permits G3BP1/2- and SG assembly-mediated PKR activation and cell-intrinsic antiviral immunity in mice and humans. DBR1-deficient patients are prone to viral disease because of intracellular lariat accumulation, which impairs G3BP1/2- and SG assembly-dependent PKR activation.

Mutations that impact splicing play a significant role in disease etiology but are not fully understood. To characterize the impact of exonic variants on splicing in 71 clinically-actionable disease genes in asymptomatic people, we analyzed 32,112 exonic mutations from ClinVar and Geisinger MyCode using a minigene reporter assay. We identify 1,733 splice-disrupting mutations, of which the most extreme 1-2% of variants are likely to be deleterious. We report that these variants are not distributed evenly across exons but are mostly concentrated in the ∼8% of exons that are most susceptible to splicing mutations (i.e. hotspot exons). We demonstrate that splicing defects in these exons can be reverted by ASOs targeting the splice sites of either their upstream or downstream flanking exons. This finding supports the feasibility of developing single therapeutic ASOs that could revert all splice-altering variants localized to a particular exon.

BACKGROUND:The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used. RESULTS:Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC. CONCLUSIONS:We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of the first viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectroscopy. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

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.

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.

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.

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.

Summary The pre-mRNA life cycle requires intron processing; yet, how intron processing defects influence splicing and gene expression is unclear. Here, we find TTDN1, which is frequently mutated in non-photosensitive trichothiodystrophy (NP-TTD), functionally links intron lariat processing to the spliceosome. The conserved TTDN1 C-terminal region directly binds lariat debranching enzyme DBR1, while its N-terminal intrinsically disordered region (IDR) binds the intron binding complex (IBC). The IDR forms condensates in vitro and is needed for IBC interaction. TTDN1 loss causes significant intron lariat accumulation, as well as splicing and gene expression defects, mirroring phenotypes observed in NP-TTD patient cells. Ttdn1 Δ/Δ mice recapitulate intron processing defects and neurodevelopmental phenotypes seen in NP-TTD. A DBR1-IDR fusion recruits DBR1 to the IBC and circumvents the requirement for TTDN1, indicating this tethering role as its major molecular function. Collectively, our findings unveil key functional connections between lariat processing, splicing outcomes, and NP-TTD molecular pathology.

Severe acute respiratory syndromeecoronavirus 2 (SARS-CoV-2) is transmitted through airborne particles in exhaled breath, causing severe respiratory disease, coronavirus diseasee2019 (COVID-19), in some patients. Samples for SARS-CoV-2 testing are typically collected by nasopharyngeal swab, with the virus detected by PCR; however, patients can test positive for 3 months after infection. Without the capacity to assay SARS-CoV-2 in breath, it is not possible to understand the risk for transmission from infected individuals. To detect virus in breath, the Bubblerda breathalyzer that reverse-transcribes RNA from SARS-CoV-2 particles into a sample-specific barcoded cDNAdwas developed. In a study of 70 hospitalized patients, the Bubbler was both more predictive of lower respiratory tract involvement (abnormal chest X-ray) and less invasive than alternatives. Samples tested using the Bubbler were threefold more enriched for SARS-CoV-2 RNA than were samples from tongue swabs, implying that virus particles were being directly sampled. The barcode-enabled Bubbler was used for simultaneous diagnosis in large batches of pooled samples at a lower limit of detection of 334 genomic copies per sample. Diagnosis by sequencing can provide additional information, such as viral load and strain identity. The Bubbler was configured to sample nucleic acids in water droplets circulating in air, demonstrating its potential in environmental monitoring and the protective effect of adequate ventilation.

High-throughput splicing assays have demonstrated that many exonic variants can disrupt splicing; however, splice-disrupting variants distribute non-uniformly across genes. We propose the existence of exons that are particularly susceptible to splice-disrupting variants, which we refer to as hotspot exons. Hotspot exons are also more susceptible to splicing perturbation through drug treatment and knock-down of RNA-binding proteins. We develop a classifier for exonic splice-disrupting variants and use it to infer hotspot exons. We estimate that 1400 exons in the human genome are hotspots. Using panels of splicing reporters, we demonstrate how the ability of an exon to tolerate a mutation is inversely proportional to the strength of its neighboring splice sites.

Exonic splicing enhancers (ESEs) are enriched in exons relative to introns and bind splicing activators. This study considers a fundamental question of co-evolution: How did ESE motifs become enriched in exons prior to the evolution of ESE recognition? We hypothesize that the high exon to intron motif ratios necessary for ESE function were created by mutational bias coupled with purifying selection on the protein code. These two forces retain certain coding motifs in exons while passively depleting them from introns. Through the use of simulations, genomic analyses, and high throughput splicing assays, we confirm the key predictions of this hypothesis, including an overlap between protein and splicing information in ESEs. We discuss the implications of mutational bias as an evolutionary driver in other cis-regulatory systems.