FANCM branchpoint translocase: Master of traverse, reverse and adverse DNA repair.

FANCM is a multifunctional DNA repair enzyme that acts as a sensor and coordinator of replication stress responses, especially interstrand crosslink (ICL) repair mediated by the Fanconi anaemia (FA) pathway. Its specialised ability to bind and remodel branched DNA structures enables diverse genome maintenance activities. Through ATP-powered "branchpoint translocation", FANCM can promote fork reversal, facilitate replication traverse of ICLs, resolve deleterious R-loop structures, and restrain recombination. These remodelling functions also support a role as sensor of perturbed replication, eliciting checkpoint signalling and recruitment of downstream repair factors like the Fanconi anaemia FANCI:FANCD2 complex. Accordingly, FANCM deficiency causes chromosome fragility and cancer susceptibility. Other recent advances link FANCM to roles in gene editing efficiency and meiotic recombination, along with emerging synthetic lethal relationships, and targeting opportunities in ALT-positive cancers. Here we review key properties of FANCM's biochemical activities, with a particular focus on branchpoint translocation as a distinguishing characteristic.

The splicing regulators RBM5 and RBM10 are subunits of the U2 snRNP engaged with intron branch sites on chromatin.

Understanding the mechanisms of pre-mRNA splicing is limited by the technical challenges to examining spliceosomes in vivo. Here, we report the isolation of RNP complexes derived from precatalytic A or B-like spliceosomes solubilized from the chromatin pellet of mammalian cell nuclei. We found that these complexes contain U2 snRNP proteins and a portion of the U2 snRNA bound with protected RNA fragments that precisely map to intronic branch sites across the transcriptome. These U2 complexes also contained the splicing regulators RBM5 and RBM10. We found RBM5 and RBM10 bound to nearly all branch site complexes and not simply those at regulated exons. The deletion of a conserved RBM5/RBM10 peptide sequence, including a zinc finger motif, disrupted U2 interaction and rendered the proteins inactive for the repression of many alternative exons. We propose a model where RBM5 and RBM10 regulate splicing as components of the U2 snRNP complex following branch site base pairing.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites.

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

ABCA4 c.6480-35A>G, a novel branchpoint variant associated with Stargardt disease.

Introduction: Inherited retinal dystrophies (IRDs) can be caused by variants in more than 280 genes. The ATP-binding cassette transporter type A4 (ABCA4) gene is one of these genes and has been linked to Stargardt disease type 1 (STGD1), fundus flavimaculatus, cone-rod dystrophy (CRD), and pan-retinal CRD. Approximately 25% of the reported ABCA4 variants affect RNA splicing. In most cases, it is necessary to perform a functional assay to determine the effect of these variants. Methods: Whole genome sequencing (WGS) was performed in one Spanish proband with Stargardt disease. The putative pathogenicity of c.6480-35A>G on splicing was investigated both in silico and in vitro. The in silico approach was based on the deep-learning tool SpliceAI. For the in vitro approach we used a midigene splice assay in HEK293T cells, based on a previously established wild-type midigene (BA29) containing ABCA4 exons 46 to 48. Results: Through the analysis of WGS data, we identified two candidate variants in ABCA4 in one proband: a previously described deletion, c.699_768+342del (p.(Gln234Phefs*5)), and a novel branchpoint variant, c.6480-35A>G. Segregation analysis confirmed that the variants were in trans. For the branchpoint variant, SpliceAI predicted an acceptor gain with a high score (0.47) at position c.6480-47. A midigene splice assay in HEK293T cells revealed the inclusion of the last 47 nucleotides of intron 47 creating a premature stop codon and allowed to categorize the variant as moderately severe. Subsequent analysis revealed the presence of this variant as a second allele besides c.1958G>A p.(Arg653His) in an additional Spanish proband in a large cohort of IRD cases. Conclusion: A splice-altering effect of the branchpoint variant, confirmed by the midigene splice assay, along with the identification of this variant in a second unrelated individual affected with STGD, provides sufficient evidence to classify the variant as likely pathogenic. In addition, this research highlights the importance of studying non-coding regions and performing functional assays to provide a conclusive molecular diagnosis.

Branchpoints as potential targets of exon-skipping therapies for genetic disorders.

Fukutin (FKTN) c.647+2084G>T creates a pseudo-exon with a premature stop codon, which causes Fukuyama congenital muscular dystrophy (FCMD). We aimed to ameliorate aberrant splicing of FKTN caused by this variant. We screened compounds focusing on splicing regulation using the c.647+2084G>T splicing reporter and discovered that the branchpoint, which is essential for splicing reactions, could be a potential therapeutic target. To confirm the effectiveness of branchpoints as targets for exon skipping, we designed branchpoint-targeted antisense oligonucleotides (BP-AONs). This restored normal FKTN mRNA and protein production in FCMD patient myotubes. We identified a functional BP by detecting splicing intermediates and creating BP mutations in the FKTN reporter gene; this BP was non-redundant and sufficiently blocked by BP-AONs. Next, a BP-AON was designed for a different FCMD-causing variant, which induces pathogenic exon trapping by a common SINE-VNTR-Alu-type retrotransposon. Notably, this BP-AON also restored normal FKTN mRNA and protein production in FCMD patient myotubes. Our findings suggest that BPs could be potential targets in exon-skipping therapeutic strategies for genetic disorders.

Evolutionary responses of energy metabolism, development, and reproduction to artificial selection for increasing heat tolerance in Drosophila subobscura.

Adaptation to warming conditions involves increased heat tolerance and metabolic changes to reduce maintenance costs and maximize biological functions close to fitness. Evidence shows that energy metabolism evolves in response to warming conditions, but we know little about how heat stress intensity determines the evolutionary responses of metabolism and life history traits. Here, we evaluated the evolutionary responses of energy metabolism and life-history traits to artificial selection for increasing heat tolerance in Drosophila subobscura, using 2 protocols to measure and select heat tolerance: slow and fast ramping protocols. We found that the increase in heat tolerance was associated with reduced activity of the enzymes involved in the glucose-6-phosphate branchpoint but no changes of the metabolic rate in selected lines. We also found that the evolution of increased heat tolerance increased the early fecundity in selected lines and increased the egg-to-adult viability only in the slow-ramping selected lines. This work shows heat tolerance can evolve under different thermal scenarios but with different evolutionary outcomes on associated traits depending on the heat stress intensity. Therefore, spatial and temporal variability of thermal stress intensity should be taken into account to understand and predict the adaptive response to ongoing and future climatic conditions.

Genome-wide detection of human variants that disrupt intronic branchpoints.

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Investigation of tracer gas transport in a new numerical model of lung acini.

Obstructive pulmonary diseases are associated with considerable morbidity. For an early diagnosis of these diseases, inert gas washouts can potentially be used. However, the complex interaction between lung anatomy and gas transport mechanisms complicates data analysis. In order to investigate this interaction, a numerical model, based on the finite difference method, consisting of two lung units connected in parallel, was developed to simulate the tracer gas transport within the human acinus. Firstly, the geometries of the units were varied and the diffusion coefficients (D) were kept constant. Secondly, D was changed and the geometry was kept constant. Furthermore, simple monoexponential growth functions were applied to evaluate the simulated data. In 109 of the 112 analyzed curves, monoexponential function matched simulated data with an accuracy of over 90%, potentially representing a suitable numerical tool to predict transport processes in further model extensions. For total flows greater than 5 × 10-4 ml/s, the exponential growth constants increased linearly with linear increasing flow to an accuracy of over 95%. The slopes of these linear trend lines of 1.23 µl-1 (D = 0.6 cm2/s), 1.69 µl-1 (D = 0.3 cm2/s), and 2.25 µl-1 (D = 0.1 cm2/s) indicated that gases with low D are more sensitive to changes in flows than gases with high D.

The splicing effect of variants at branchpoint elements in cancer genes.

PURPOSE: Branchpoint elements are required for intron removal, and variants at these elements can result in aberrant splicing. We aimed to assess the value of branchpoint annotations generated from recent large-scale studies to select branchpoint-abrogating variants, using hereditary cancer genes as model. METHODS: We identified branchpoint elements in 119 genes associated with hereditary cancer from 3 genome-wide experimentally-inferred and 2 predicted branchpoint data sets. We then identified variants that occur within branchpoint elements from public databases. We compared conservation, unique variant observations, and population frequencies at different nucleotides within branchpoint motifs. Finally, selected minigene assays were performed to assess the splicing effect of variants at branchpoint elements within mismatch repair genes. RESULTS: There was poor overlap between predicted and experimentally-inferred branchpoints. Our analysis of cancer genes suggested that variants at -2 nucleotide, -1 nucleotide, and branchpoint positions in experimentally-inferred canonical motifs are more likely to be clinically relevant. Minigene assay data showed the -2 nucleotide to be more important to branchpoint motif integrity but also showed fluidity in branchpoint usage. CONCLUSION: Data from cancer gene analysis suggest that there are few high-risk alleles that severely impact function via branchpoint abrogation. Results of this study inform a general scheme to prioritize branchpoint motif variants for further study.

Splicing efficiency of minor introns in a mouse model of SMA predominantly depends on their branchpoint sequence and can involve the contribution of major spliceosome components.

Spinal muscular atrophy (SMA) is a devastating neurodegenerative disease caused by reduced amounts of the ubiquitously expressed Survival of Motor Neuron (SMN) protein. In agreement with its crucial role in the biogenesis of spliceosomal snRNPs, SMN-deficiency is correlated to numerous splicing alterations in patient cells and various tissues of SMA mouse models. Among the snRNPs whose assembly is impacted by SMN-deficiency, those involved in the minor spliceosome are particularly affected. Importantly, splicing of several, but not all U12-dependent introns has been shown to be affected in different SMA models. Here, we have investigated the molecular determinants of this differential splicing in spinal cords from SMA mice. We show that the branchpoint sequence (BPS) is a key element controlling splicing efficiency of minor introns. Unexpectedly, splicing of several minor introns with suboptimal BPS is not affected in SMA mice. Using in vitro splicing experiments and oligonucleotides targeting minor or major snRNAs, we show for the first time that splicing of these introns involves both the minor and major machineries. Our results strongly suggest that splicing of a subset of minor introns is not affected in SMA mice because components of the major spliceosome compensate for the loss of minor splicing activity.

Characterization of the aberrant splicing of DVL2 induced by cancer-associated SF3B1 mutation.

SF3B1, an essential component of the U2 snRNP, is frequently mutated in cancers. Cancer-associated SF3B1 mutation causes aberrant RNA splicing, mostly at 3' splice sites (3'ss). RNA splicing of DVL2, a regulator of Notch signaling, is affected by SF3B1 mutation. Here, we report that the mutated SF3B1 use an alternative branchpoint sequence (BPS) for the aberrant splicing of DVL2, which has a higher affinity to U2 snRNA than the BPS for the canonical splicing of DVL2. Swapping the position of the alternative BPS with the position of the canonical BPS decreased the aberrant splicing of DVL2, suggesting that the mutated SF3B1 prefers to use BPS with high affinity to U2 snRNA for splicing. Additionally, swapping the positions of two BPSs associated with the canonical splicing of DVL2 demonstrated that both the affinity to the U2 snRNA and the distance to the 3'ss are important to the selection of BPS. Importantly, the aberrant splicing of DVL2 does not require the canonical 3'ss and the canonical polypyrimidine tract, which reveals a novel type of aberrant splicing induced by SF3B1 mutation. These findings provide a more comprehensive understanding of the mechanisms underlying aberrant splicing induced by SF3B1 mutation in cancer.

Rewards of divergence in sequences, 3-D structures and dynamics of yeast and human spliceosome SF3b complexes.

The evolution of homologous and functionally equivalent multiprotein assemblies is intriguing considering sequence divergence of constituent proteins. Here, we studied the implications of protein sequence divergence on the structure, dynamics and function of homologous yeast and human SF3b spliceosomal subcomplexes. Human and yeast SF3b comprise of 7 and 6 proteins respectively, with all yeast proteins homologous to their human counterparts at moderate sequence identity. SF3b6, an additional component in the human SF3b, interacts with the N-terminal extension of SF3b1 while the yeast homologue Hsh155 lacks the equivalent region. Through detailed homology studies, we show that SF3b6 is absent not only in yeast but in multiple lineages of eukaryotes implying that it is critical in specific organisms. We probed for the potential role of SF3b6 in the spliceosome assembled form through structural and flexibility analyses. By analysing normal modes derived from anisotropic network models of SF3b1, we demonstrate that when SF3b1 is bound to SF3b6, similarities in the magnitude of residue motions (0.86) and inter-residue correlated motions (0.94) with Hsh155 are significantly higher than when SF3b1 is considered in isolation (0.21 and 0.89 respectively). We observed that SF3b6 promotes functionally relevant 'open-to-close' transition in SF3b1 by enhancing concerted residue motions. Such motions are found to occur in the Hsh155 without SF3b6. The presence of SF3b6 influences motions of 16 residues that interact with U2 snRNA/branchpoint duplex and supports the participation of its interface residues in long-range communication in the SF3b1. These results advocate that SF3b6 potentially acts as an allosteric regulator of SF3b1 for BPS selection and might play a role in alternative splicing. Furthermore, we observe variability in the relative orientation of SF3b4 and in the local structure of three ß-propeller domains of SF3b3 with reference to their yeast counterparts. Such differences influence the inter-protein interactions of SF3b between these two organisms. Together, our findings highlight features of SF3b evolution and suggests that the human SF3b may have evolved sophisticated mechanisms to fine tune its molecular function.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars.

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

New insights into diosgenin biosynthesis pathway and its regulation in Trigonella foenum-graecum L.

INTRODUCTION: Throughout history, thousands of medicinal and aromatic plants have been widely utilised by people worldwide. Owing to them possessing of valuable compounds with little side effects in comparison with chemical drugs, herbs have been of interest to humans for a number of purposes. Diosgenin, driven from fenugreek, Trigonella foenum-graecum L., has extensively drawn scientist's attention owing to having curable properties and being a precursor of steroid hormones synthesis. Nonetheless, complete knowledge about the biosynthesis pathway of this metabolite is still elusive. OBJECTIVE: In the present research, we isolated the full-length CDS of 14 genes involving in diosgenin formation and measured their expression rate in various genotypes, which had illustrated different amount of diosgenin. METHODOLOGY: The genes were successfully isolated, and functional motifs were also assessed using in silico approaches. RESULTS: Moreover, combining transcript and metabolite analysis revealed that there are many genes playing the role in diosgenin formation, some of which are highly influential. Among them, ∆24 -reductase, which converts cycloartenol to cycloartanol, is the first-committed and rate-limiting enzyme in this pathway. Additionally, no transcripts indicating to the presence or expression of lanosterol synthase were detected, contradicting the previous hypothesis about the biosynthetic pathway of diosgenin in fenugreek. CONCLUSION: Considering all these, therefore, we propose the most possible pathway of diosgenin. This knowledge will then pave the way toward cloning the genes as well as engineering the diosgenin biosynthesis pathway.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity.

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R1C2T3R4A5Y6, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G4/A4 ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G4 increase is consistent with A4 to G4 transitions in evolution. Of the G4/A4-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G3/A3 co-vary with the G4 to some extent. However, corresponding increase of the G4-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G4-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G4/A4 ratio correlates well with the abundance of alternative splicing in the two phyla, and G4 enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

Pathogenic Abnormal Splicing Due to Intronic Deletions that Induce Biophysical Space Constraint for Spliceosome Assembly.

A precise genetic diagnosis is the single most important step for families with genetic disorders to enable personalized and preventative medicine. In addition to genetic variants in coding regions (exons) that can change a protein sequence, abnormal pre-mRNA splicing can be devastating for the encoded protein, inducing a frameshift or in-frame deletion/insertion of multiple residues. Non-coding variants that disrupt splicing are extremely challenging to identify. Stemming from an initial clinical discovery in two index Australian families, we define 25 families with genetic disorders caused by a class of pathogenic non-coding splice variant due to intronic deletions. These pathogenic intronic deletions spare all consensus splice motifs, though they critically shorten the minimal distance between the 5' splice-site (5'SS) and branchpoint. The mechanistic basis for abnormal splicing is due to biophysical constraint precluding U1/U2 spliceosome assembly, which stalls in A-complexes (that bridge the 5'SS and branchpoint). Substitution of deleted nucleotides with non-specific sequences restores spliceosome assembly and normal splicing, arguing against loss of an intronic element as the primary causal basis. Incremental lengthening of 5'SS-branchpoint length in our index EMD case subject defines 45-47 nt as the critical elongation enabling (inefficient) spliceosome assembly for EMD intron 5. The 5'SS-branchpoint space constraint mechanism, not currently factored by genomic informatics pipelines, is relevant to diagnosis and precision medicine across the breadth of Mendelian disorders and cancer genomics.

Cus2 enforces the first ATP-dependent step of splicing by binding to yeast SF3b1 through a UHM-ULM interaction.

Stable recognition of the intron branchpoint (BP) by the U2 snRNP to form the pre-spliceosome is the first ATP-dependent step of splicing. Genetic and biochemical data from yeast indicate that Cus2 aids U2 snRNA folding into the stem IIa conformation prior to pre-spliceosome formation. Cus2 must then be removed by an ATP-dependent function of Prp5 before assembly can progress. However, the location from which Cus2 is displaced and the nature of its binding to the U2 snRNP are unknown. Here, we show that Cus2 contains a conserved UHM (U2AF homology motif) that binds Hsh155, the yeast homolog of human SF3b1, through a conserved ULM (U2AF ligand motif). Mutations in either motif block binding and allow pre-spliceosome formation without ATP. A 2.0 Å resolution structure of the Hsh155 ULM in complex with the UHM of Tat-SF1, the human homolog of Cus2, and complementary binding assays show that the interaction is highly similar between yeast and humans. Furthermore, we show that Tat-SF1 can replace Cus2 function by enforcing ATP dependence of pre-spliceosome formation in yeast extracts. Cus2 is removed before pre-spliceosome formation, and both Cus2 and its Hsh155 ULM binding site are absent from available cryo-EM structure models. However, our data are consistent with the apparent location of the disordered Hsh155 ULM between the U2 stem-loop IIa and the HEAT repeats of Hsh155 that interact with Prp5. We propose a model in which Prp5 uses ATP to remove Cus2 from Hsh155 such that extended base-pairing between U2 snRNA and the intron BP can occur.

Dual coding potential of a 2',5'-branched ribonucleotide in DNA.

Branchpoints in RNA templates are highly mutagenic, but it is not known yet whether this also applies to branchpoints in DNA templates. Here, we report how nucleic acid polymerases replicate a 2',5'-branched DNA (bDNA) molecule. We constructed long-chained bDNA templates containing a branch guanosine and T7 promoters at both arms by splinted ligation. Quantitative real-time PCR analysis was used to investigate whether a branchpoint blocks DNA synthesis from the two arms in the same manner. We find that the blocking effect of a branchpoint is arm-specific. DNA synthesis from the 2'-arm is more than 20,000-fold decreased, whereas from the 3'-arm only 15-fold. Our sequence analysis of full-length nucleic acid generated by Taq DNA polymerase, Moloney murine leukemia virus reverse transcriptase, and T7 RNA polymerase from the 2'-arm of bDNA shows that the branched guanine has a dual coding potential and can base-pair with cytosine and guanine. We find that branchpoint templating is influenced by the type of the surrounding nucleic acid and is probably modulated by polymerase and RNase H active sites. We show that the branchpoint bypass by the polymerases from the 3'-arm of bDNA is predominantly error-free, indicating that bDNA is not as highly mutagenic as 2',5'-branched RNA.

Transcriptome-wide Interrogation of the Functional Intronome by Spliceosome Profiling.

Full understanding of eukaryotic transcriptomes and how they respond to different conditions requires deep knowledge of all sites of intron excision. Although RNA sequencing (RNA-seq) provides much of this information, the low abundance of many spliced transcripts (often due to their rapid cytoplasmic decay) limits the ability of RNA-seq alone to reveal the full repertoire of spliced species. Here, we present "spliceosome profiling," a strategy based on deep sequencing of RNAs co-purifying with late-stage spliceosomes. Spliceosome profiling allows for unambiguous mapping of intron ends to single-nucleotide resolution and branchpoint identification at unprecedented depths. Our data reveal hundreds of new introns in S. pombe and numerous others that were previously misannotated. By providing a means to directly interrogate sites of spliceosome assembly and catalysis genome-wide, spliceosome profiling promises to transform our understanding of RNA processing in the nucleus, much as ribosome profiling has transformed our understanding mRNA translation in the cytoplasm.