Group II Intron RNPs and Reverse Transcriptases: From Retroelements to Research Tools

Cryo-EM: A window into the dynamic world of RNA molecules

RNAs are critical for complex cellular functions, characterized by their structural versatility and ability to undergo conformational transitions in response to cellular cues. The elusive structures of RNAs are being unraveled with unprecedented clarity, thanks to the technological advancements in structural biology, including nuclear magnetic resonance (NMR), X-ray crystallography, cryo-electron microscopy (cryo-EM) etc. This review focuses on examining the revolutionary impact of cryo-EM on our comprehension of RNA structural dynamics, underscoring the technique's contributions to structural biology and envisioning the future trajectory of this rapidly evolving field.

Orientia tsutsugamushi: comprehensive analysis of the mobilome of a highly fragmented and repetitive genome reveals the capacity for ongoing lateral gene transfer in an obligate intracellular bacterium

IMPORTANCE

Obligate intracellular bacteria, or those only capable of growth inside other living cells, have limited opportunities for horizontal gene transfer with other microbes due to their isolated replicative niche. The human pathogen Ot, an obligate intracellular bacterium causing scrub typhus, encodes an unusually high copy number of a ~40 gene mobile genetic element that typically facilitates genetic transfer across microbes. This proliferated element is heavily degraded in Ot and previously assumed to be inactive. Here, we conducted a detailed analysis of this element in eight Ot strains and discovered two strains with at least one intact copy. This implies that the element is still capable of moving across Ot populations and suggests that the genome of this bacterium may be even more dynamic than previously appreciated. Our work raises questions about intracellular microbial evolution and sounds an alarm for gene-based efforts focused on diagnosing and combatting scrub typhus.

Genome-Wide Profiling of tRNA Using an Unexplored Reverse Transcriptase with High Processivity

Monitoring the dynamic changes of cellular tRNA pools is challenging, due to the extensive post-transcriptional modifications of individual species. The most critical component in tRNAseq is a processive reverse transcriptase (RT) that can read through each modification with high efficiency. Here we show that the recently developed group-II intron RT Induro has the processivity and efficiency necessary to profile tRNA dynamics. Using our Induro-tRNAseq, simpler and more comprehensive than the best methods to date, we show that Induro progressively increases readthrough of tRNA over time and that the mechanism of increase is selective removal of RT stops, without altering the misincorporation frequency. We provide a parallel dataset of the misincorporation profile of Induro relative to the related TGIRT RT to facilitate the prediction of non-annotated modifications. We report an unexpected modification profile among human proline isoacceptors, absent from mouse and lower eukaryotes, that indicates new biology of decoding proline codons.

Structural insights into intron catalysis and dynamics during splicing

The group II intron ribonucleoprotein is an archetypal splicing system with numerous mechanistic parallels to the spliceosome, including excision of lariat introns1,2. Despite the importance of branching in RNA metabolism, structural understanding of this process has remained elusive. Here we present a comprehensive analysis of three single-particle cryogenic electron microscopy structures captured along the splicing pathway. They reveal the network of molecular interactions that specifies the branchpoint adenosine and positions key functional groups to catalyse lariat formation and coordinate exon ligation. The structures also reveal conformational rearrangements of the branch helix and the mechanism of splice site exchange that facilitate the transition from branching to ligation. These findings shed light on the evolution of splicing and highlight the conservation of structural components, catalytic mechanism and dynamical strategies retained through time in premessenger RNA splicing machines.

The mitogenomes of Leptographium aureum, Leptographium sp., and Grosmannia fruticeta: expansion by introns

Introduction

Many members of the Ophiostomatales are of economic importance as they are bark-beetle associates and causative agents for blue stain on timber and in some instances contribute towards tree mortality. The taxonomy of these fungi has been challenging due to the convergent evolution of many traits associated with insect dispersal and a limited number of morphological characters that happen to be highly pleomorphic. This study examines the mitochondrial genomes for three members of Leptographium sensu lato [Leptographium aureum (also known as Grosmannia aurea), Grosmannia fruticeta (also known as Leptographium fruticetum), and Leptographium sp. WIN(M)1376)].

Methods

Illumina sequencing combined with gene and intron annotations and phylogenetic analysis were performed.

Results

Sequence analysis showed that gene content and gene synteny are conserved but mitochondrial genome sizes were variable: G. fruticeta at 63,821 bp, Leptographium sp. WIN(M)1376 at 81,823 bp and L. aureum at 104,547 bp. The variation in size is due to the number of introns and intron-associated open reading frames. Phylogenetic analysis of currently available mitochondrial genomes for members of the Ophiostomatales supports currently accepted generic arrangements within this order and specifically supports the separation of members with Leptographium-like conidiophores into two genera, with L. aureum grouping with Leptographium and G. fruticeta aligning with Grosmannia.

Discussion

Mitochondrial genomes are promising sequences for resolving evolutionary relationships within the Ophiostomatales.

Inflammatory breast cancer biomarker identification by simultaneous TGIRT-seq profiling of coding and non-coding RNAs in tumors and blood

Inflammatory breast cancer (IBC) is the most aggressive and lethal breast cancer subtype, but lags in biomarker identification. Here, we used an improved Thermostable Group II Intron Reverse Transcriptase RNA sequencing (TGIRT-seq) method to simultaneously profile coding and non-coding RNAs from tumors, PBMCs, and plasma of IBC and non-IBC patients and healthy donors. Besides RNAs from known IBC-relevant genes, we identified hundreds of other overexpressed coding and non-coding RNAs (p≤0.001) in IBC tumors and PBMCs, including higher proportions with elevated intron-exon depth ratios (IDRs), likely reflecting enhanced transcription resulting in accumulation of intronic RNAs. As a consequence, differentially represented protein-coding gene RNAs in IBC plasma were largely intron RNA fragments, whereas those in healthy donor and non-IBC plasma were largely fragmented mRNAs. Potential IBC biomarkers in plasma included T-cell receptor pre-mRNA fragments traced to IBC tumors and PBMCs; intron RNA fragments correlated with high IDR genes; and LINE-1 and other retroelement RNAs that we found globally up-regulated in IBC and preferentially enriched in plasma. Our findings provide new insights into IBC and demonstrate advantages of broadly analyzing transcriptomes for biomarker identification. The RNA-seq and data analysis methods developed for this study may be broadly applicable to other diseases.

Orientia tsutsugamushi: analysis of the mobilome of a highly fragmented and repetitive genome reveals ongoing lateral gene transfer in an obligate intracellular bacterium

The rickettsial human pathogen Orientia tsutsugamushi (Ot) is an obligate intracellular Gram-negative bacterium with one of the most highly fragmented and repetitive genomes of any organism. Around 50% of its ~2.3 Mb genome is comprised of repetitive DNA that is derived from the highly proliferated Rickettsiales amplified genetic element (RAGE). RAGE is an integrative and conjugative element (ICE) that is present in a single Ot genome in up to 92 copies, most of which are partially or heavily degraded. In this report, we analysed RAGEs in eight fully sequenced Ot genomes and manually curated and reannotated all RAGE-associated genes, including those encoding DNA mobilisation proteins, P-type (vir) and F-type (tra) type IV secretion system (T4SS) components, Ankyrin repeat- and tetratricopeptide repeat-containing effectors, and other piggybacking cargo. Originally, the heavily degraded Ot RAGEs led to speculation that they are remnants of historical ICEs that are no longer active. Our analysis, however, identified two Ot genomes harbouring one or more intact RAGEs with complete F-T4SS genes essential for mediating ICE DNA transfer. As similar ICEs have been identified in unrelated rickettsial species, we assert that RAGEs play an ongoing role in lateral gene transfer within the Rickettsiales. Remarkably, we also identified in several Ot genomes remnants of prophages with no similarity to other rickettsial prophages. Together these findings indicate that, despite their obligate intracellular lifestyle and host range restricted to mites, rodents and humans, Ot genomes are highly dynamic and shaped through ongoing invasions by mobile genetic elements and viruses.

A New EBS2b-IBS2b Base Paring (A-8/T-8) Improved the Gene-Targeting Efficiency of Thermotargetron in Escherichia coli

Thermophilic group II intron is one type of retrotransposon composed of intron RNA and intron-encoded protein (IEP), which can be utilized in gene targeting by harnessing their novel ribozyme-based DNA integration mechanism termed "retrohoming." It is mediated by a ribonucleoprotein (RNP) complex that contains the excised intron lariat RNA and an IEP with reverse transcriptase (RT) activity. The RNP recognizes targeting sites by exon-binding sequences 2 (EBS2)/intron-binding sequences 2 (IBS2), EBS1/IBS1, and EBS3/IBS3 bases pairing. Previously, we developed the TeI3c/4c intron as a thermophilic gene targeting system-Thermotargetron (TMT). However, we found that the targeting efficiency of TMT varies significantly at different targeting sites, which leads to a relatively low success rate. To further improve the success rate and gene-targeting efficiency of TMT, we constructed a Random Gene-targeting Plasmids Pool (RGPP) to analyze the sequence recognition preference of TMT. A new base pairing, located at the -8 site between EBS2/IBS2 and EBS1/IBS1 (named EBS2b-IBS2b), increased the success rate (2.45- to 5.07-fold) and significantly improved gene-targeting efficiency of TMT. A computer algorithm (TMT 1.0), based on the newly discovered sequence recognition roles, was also developed to facilitate the design of TMT gene-targeting primers. The present work could essentially expand the practicalities of TMT in the genome engineering of heat-tolerance mesophilic and thermophilic bacteria. IMPORTANCE The randomized base pairing in the interval of IBS2 and IBS1 of Tel3c/4c intron (-8 and -7 sites) in Thermotargetron (TMT) results in a low success rate and gene-targeting efficiency in bacteria. In the present work, we constructed a randomized gene-targeting plasmids pool (RGPP) to study whether there is a base preference in target sequences. Among all the successful "retrohoming" targets, we found that a new EBS2b-IBS2b base paring (A-8/T-8) significantly increased TMT's gene-targeting efficiency, and the concept is also applicable to other gene targets in redesigned gene-targeting plasmids pool in E. coli. The improved TMT is a promising tool for the genetic engineering of bacteria and could promote metabolic engineering and synthetic biology research in valuable microbes that recalcitrance for genetic manipulation.

Pitfalls of using sequence databases for heterologous expression studies - a technical review

Synthesis of DNA fragments based on gene sequences that are available in public resources has become an efficient and affordable method that has gradually replaced traditional cloning efforts such as PCR cloning from cDNA. However, database entries based on genome sequencing results are prone to errors which can lead to false sequence information and, ultimately, errors in functional characterisation of proteins such as ion channels and transporters in heterologous expression systems. We have identified five common problems that repeatedly appear in public resources: (1) Not every gene has yet been annotated; (2) not all gene annotations are necessarily correct; (3) transcripts may contain automated corrections; (4) there are mismatches between gene, mRNA and protein sequences; and (5) splicing patterns often lack experimental validation. This technical review highlights and provides a strategy to bypass these issues in order to avoid critical mistakes that could impact future studies of any gene/protein of interest in heterologous expression systems.

The first two mitochondrial genomes for the genus Ramaria reveal mitochondrial genome evolution of Ramaria and phylogeny of Basidiomycota

In the present study, we assembled and analyzed the mitogenomes of two Ramaria species. The assembled mitogenomes of Ramaria cfr. rubripermanens and R. rubella were circularized, with sizes of 126,497 bp and 143,271 bp, respectively. Comparative mitogenome analysis showed that intron region contributed the most (contribution rate, 43.74%) to the size variations of Ramaria mitogenomes. The genetic contents, gene length, tRNAs, and codon usages of the two Ramaria mitogenomes varied greatly. In addition, the evolutionary rates of different core protein coding genes (PCGs) in Phallomycetidae mitogenomes varied. We detected large-scale gene rearrangements between Phallomycetidae mitogenomes, including gene displacement and tRNA doubling. A total of 4499 bp and 7746 bp aligned fragments were detected between the mitochondrial and nuclear genomes of R. cfr. rubripermanens and R. rubella, respectively, indicating possible gene transferring events. We further found frequent intron loss/gain and potential intron transfer events in Phallomycetidae mitogenomes during the evolution, and the mitogenomes of R. rubella contained a novel intron P44. Phylogenetic analyses using both Bayesian inference (BI) and Maximum Likelihood (ML) methods based on a combined mitochondrial gene dataset obtained an identical and well-supported phylogenetic tree for Basidiomycota, wherein R. cfr. rubripermanens and Turbinellus floccosus are sister species. This study served as the first report on mitogenomes from the genus Ramaria, which provides a basis for understanding the evolution, genetics, and taxonomy of this important fungal group.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

Identification of Group II Intron RmInt1 Binding Sites in a Bacterial Genome

RmInt1 is a group II intron encoding a reverse transcriptase protein (IEP) lacking the C-terminal endonuclease domain. RmInt1 is an efficient mobile retroelement that predominantly reverse splices into the transient single-stranded DNA at the template for lagging strand DNA synthesis during host replication, a process facilitated by the interaction of the RmInt1 IEP with DnaN at the replication fork. It has been suggested that group II intron ribonucleoprotein particles bind DNA nonspecifically, and then scan for their correct target site. In this study, we investigated RmInt1 binding sites throughout the Sinorhizobium meliloti genome, by chromatin-immunoprecipitation coupled with next-generation sequencing. We found that RmInt1 binding sites cluster around the bidirectional replication origin of each of the three replicons comprising the S. meliloti genome. Our results provide new evidence linking group II intron mobility to host DNA replication.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

A new RNA-DNA interaction required for integration of group II intron retrotransposons into DNA targets

Mobile group II introns are site-specific retrotransposable elements abundant in bacterial and organellar genomes. They are composed of a large and highly structured ribozyme and an intron-encoded reverse transcriptase that binds tightly to its intron to yield a ribonucleoprotein (RNP) particle. During the first stage of the mobility pathway, the intron RNA catalyses its own insertion directly into the DNA target site. Recognition of the proper target rests primarily on multiple base-pairing interactions between the intron RNA and the target DNA, while the protein makes contacts with only a few target positions by yet-unidentified mechanisms. Using a combination of comparative sequence analyses and in vivo mobility assays we demonstrate the existence of a new base-pairing interaction named EBS2a–IBS2a between the intron RNA and its DNA target site. This pairing adopts a Watson–Crick geometry and is essential for intron mobility, most probably by driving unwinding of the DNA duplex. Importantly, formation of EBS2a–IBS2a also requires the reverse transcriptase enzyme which stabilizes the pairing in a non-sequence-specific manner. In addition to bringing to light a new structural device that allows subgroup IIB1 and IIB2 introns to invade their targets with high efficiency and specificity our work has important implications for the biotechnological applications of group II introns in bacterial gene targeting.

Prokaryotic reverse transcriptases: from retroelements to specialized defense systems

Reverse transcriptases (RTs) catalyze the polymerization of DNA from an RNA template. These enzymes were first discovered in RNA tumor viruses in 1970, but it was not until 1989 that they were found in prokaryotes as a key component of retrons. Apart from RTs encoded by the 'selfish' mobile retroelements known as group II introns, prokaryotic RTs are extraordinarily diverse, but their function has remained elusive. However, recent studies have revealed that different lineages of prokaryotic RTs, including retrons, those associated with CRISPR-Cas systems, Abi-like RTs and other yet uncharacterized RTs, are key components of different lines of defense against phages and other mobile genetic elements. Prokaryotic RTs participate in various antiviral strategies, including abortive infection (Abi), in which the infected cell is induced to commit suicide to protect the host population, adaptive immunity, in which a memory of previous infection is used to build an efficient defense, and other as yet unidentified mechanisms. These prokaryotic enzymes are attracting considerable attention, both for use in cutting-edge technologies, such as genome editing, and as an emerging research topic. In this review, we discuss what is known about prokaryotic RTs, and the exciting evidence for their domestication from retroelements to create specialized defense systems.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Bacterial Group II Intron Genomic Neighborhoods Reflect Survival Strategies: Hiding and Hijacking

Group II (gII) introns are mobile retroelements that can spread to new DNA sites through retrotransposition, which can be influenced by a variety of host factors. To determine if these host factors bear any relationship to the genomic location of gII introns, we developed a bioinformatic pipeline wherein we focused on the genomic neighborhoods of bacterial gII introns within their native contexts and sought to determine global relationships between introns and their surrounding genes. We found that, although gII introns inhabit diverse regions, these neighborhoods are often functionally enriched for genes that could promote gII intron retention or proliferation. On one hand, we observe that gII introns are frequently found hiding in mobile elements or after transcription terminators. On the other hand, gII introns are enriched in locations in which they could hijack host functions for their movement, potentially timing expression of the intron with genes that produce favorable conditions for retrotransposition. Thus, we propose that gII intron distributions have been shaped by relationships with their surrounding genomic neighbors.

Methylation of rRNA as a host defense against rampant group II intron retrotransposition

BACKGROUND

Group II introns are mobile retroelements, capable of invading new sites in DNA. They are self-splicing ribozymes that complex with an intron-encoded protein to form a ribonucleoprotein that targets DNA after splicing. These molecules can invade DNA site-specifically, through a process known as retrohoming, or can invade ectopic sites through retrotransposition. Retrotransposition, in particular, can be strongly influenced by both environmental and cellular factors.

RESULTS

To investigate host factors that influence retrotransposition, we performed random insertional mutagenesis using the ISS1 transposon to generate a library of over 1000 mutants in Lactococcus lactis, the native host of the Ll.LtrB group II intron. By screening this library, we identified 92 mutants with increased retrotransposition frequencies (RTP-ups). We found that mutations in amino acid transport and metabolism tended to have increased retrotransposition frequencies. We further explored a subset of these RTP-up mutants, the most striking of which is a mutant in the ribosomal RNA methyltransferase rlmH, which exhibited a reproducible 20-fold increase in retrotransposition frequency. In vitro and in vivo experiments revealed that ribosomes in the rlmH mutant were defective in the m3Ψ modification and exhibited reduced binding to the intron RNA.

CONCLUSIONS

Taken together, our results reinforce the importance of the native host organism in regulating group II intron retrotransposition. In particular, the evidence from the rlmH mutant suggests a role for ribosome modification in limiting rampant retrotransposition.

The mitochondrial genome of Ophiostoma himal-ulmi and comparison with other fungi causing Dutch elm disease

The mitochondrial genome of Ophiostoma himal-ulmi, a species endemic to the Western Himalayas and one of the fungi that cause Dutch elm disease, has been sequenced and characterized. The mitochondrial genome was compared with other available genomes for members of the Ophiostomatales, including other agents of Dutch elm disease (Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi, and Ophiostoma novo-ulmi subspecies americana), and it was observed that gene synteny is highly conserved, and variability among members of the fungi that cause Dutch-elm disease is primarily due to the number of intron insertions. Among the fungi that cause Dutch elm disease that we examined, O. himal-ulmi has the largest mitochondrial genomes (ranging from 94 934 to 111 712 bp), owing to the expansion of the number of introns.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Exon and protein positioning in a pre-catalytic group II intron RNP primed for splicing

Group II introns are the putative progenitors of nuclear spliceosomal introns and use the same two-step splicing pathway. In the cell, the intron RNA forms a ribonucleoprotein (RNP) complex with the intron-encoded protein (IEP), which is essential for splicing. Although structures of spliced group II intron RNAs and RNP complexes have been characterized, structural insights into the splicing process remain enigmatic due to lack of pre-catalytic structural models. Here, we report two cryo-EM structures of endogenously produced group II intron RNPs trapped in their pre-catalytic state. Comparison of the catalytically activated precursor RNP to its previously reported spliced counterpart allowed identification of key structural rearrangements accompanying splicing, including a remodeled active site and engagement of the exons. Importantly, altered RNA-protein interactions were observed upon splicing among the RNP complexes. Furthermore, analysis of the catalytically inert precursor RNP demonstrated the structural impact of the formation of the active site on RNP architecture. Taken together, our results not only fill a gap in understanding the structural basis of IEP-assisted group II intron splicing, but also provide parallels to evolutionarily related spliceosomal splicing.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Zinc-finger protein CNBP alters the 3-D structure of lncRNA Braveheart in solution

Long non-coding RNAs (lncRNAs) constitute a significant fraction of the transcriptome, playing important roles in development and disease. However, our understanding of structure-function relationships for this emerging class of RNAs has been limited to secondary structures. Here, we report the 3-D atomistic structural study of epigenetic lncRNA, Braveheart (Bvht), and its complex with CNBP (Cellular Nucleic acid Binding Protein). Using small angle X-ray scattering (SAXS), we elucidate the ensemble of Bvht RNA conformations in solution, revealing that Bvht lncRNA has a well-defined, albeit flexible 3-D structure that is remodeled upon CNBP binding. Our study suggests that CNBP binding requires multiple domains of Bvht and the RHT/AGIL RNA motif. We show that RHT/AGIL, previously shown to interact with CNBP, contains a highly flexible loop surrounded by more ordered helices. As one of the largest RNA-only 3-D studies, the work lays the foundation for future structural studies of lncRNA-protein complexes.

Branch site bulge conformations in domain 6 determine functional sugar puckers in group II intron splicing

Although group II intron ribozymes are intensively studied the question how structural dynamics affects splicing catalysis has remained elusive. We report for the first time that the group II intron domain 6 exists in a secondary structure equilibrium between a single- and a two-nucleotide bulge conformation, which is directly linked to a switch between sugar puckers of the branch site adenosine. Our study determined a functional sugar pucker equilibrium between the transesterification active C2'-endo conformation of the branch site adenosine in the 1nt bulge and an inactive C3'-endo state in the 2nt bulge fold, allowing the group II intron to switch its activity from the branching to the exon ligation step. Our detailed NMR spectroscopic investigation identified magnesium (II) ions and the branching reaction as regulators of the equilibrium populations. The tuneable secondary structure/sugar pucker equilibrium supports a conformational selection mechanism to up- and downregulate catalytically active and inactive states of the branch site adenosine to orchestrate the multi-step splicing process. The conformational dynamics of group II intron domain 6 is also proposed to be a key aspect for the directionality selection in reversible splicing.

Improved TGIRT-seq methods for comprehensive transcriptome profiling with decreased adapter dimer formation and bias correction

Thermostable group II intron reverse transcriptases (TGIRTs) with high fidelity and processivity have been used for a variety of RNA sequencing (RNA-seq) applications, including comprehensive profiling of whole-cell, exosomal, and human plasma RNAs; quantitative tRNA-seq based on the ability of TGIRT enzymes to give full-length reads of tRNAs and other structured small ncRNAs; high-throughput mapping of post-transcriptional modifications; and RNA structure mapping. Here, we improved TGIRT-seq methods for comprehensive transcriptome profiling by rationally designing RNA-seq adapters that minimize adapter dimer formation. Additionally, we developed biochemical and computational methods for remediating 5′- and 3′-end biases, the latter based on a random forest regression model that provides insight into the contribution of different factors to these biases. These improvements, some of which may be applicable to other RNA-seq methods, increase the efficiency of TGIRT-seq library construction and improve coverage of very small RNAs, such as miRNAs. Our findings provide insight into the biochemical basis of 5′- and 3′-end biases in RNA-seq and suggest general approaches for remediating biases and decreasing adapter dimer formation.