Optimized design of antisense oligomers for targeted rRNA depletion

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

RIP-PEN-seq identifies a class of kink-turn RNAs as splicing regulators

A kink-turn (K-turn) is a three-dimensional RNA structure that exists in all three primary phylogenetic domains. In this study, we developed the RIP-PEN-seq method to identify the full-length sequences of RNAs bound by the K-turn binding protein 15.5K and discovered a previously uncharacterized class of RNAs with backward K-turn motifs (bktRNAs) in humans and mice. All bktRNAs share two consensus sequence motifs at their fixed terminal position and have complex folding properties, expression and evolution patterns. We found that a highly conserved bktRNA1 guides the methyltransferase fibrillarin to install RNA methylation of U12 small nuclear RNA in humans. Depletion of bktRNA1 causes global splicing dysregulation of U12-type introns by impairing the recruitment of ZCRB1 to the minor spliceosome. Most bktRNAs regulate the splicing of local introns by interacting with the 15.5K protein. Taken together, our findings characterize a class of small RNAs and uncover another layer of gene expression regulation that involves crosstalk among bktRNAs, RNA splicing and RNA methylation.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The mammalian mitochondrial genome encodes thirteen oxidative phosphorylation system proteins, crucial in aerobic energy transduction. These proteins are translated from 9 monocistronic and 2 bicistronic transcripts, whose native structures remain unexplored, leaving fundamental molecular determinants of mitochondrial gene expression unknown. To address this gap, we developed a mitoDMS-MaPseq approach and used DREEM clustering to resolve the native human mitochondrial mt-mRNA structurome. We gained insights into mt-mRNA biology and translation regulatory mechanisms, including a unique programmed ribosomal frameshifting for the ATP8/ATP6 transcript. Furthermore, absence of the mt-mRNA maintenance factor LRPPRC led to a mitochondrial transcriptome structured differently, with specific mRNA regions exhibiting increased or decreased structuredness. This highlights the role of LRPPRC in maintaining mRNA folding to promote mt-mRNA stabilization and efficient translation. In conclusion, our mt-mRNA folding maps reveal novel mitochondrial gene expression mechanisms, serving as a detailed reference and tool for studying them in different physiological and pathological contexts.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus. Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities similar to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far particular to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

CRISPR/Cas9-based depletion of 16S ribosomal RNA improves library complexity of single-cell RNA-sequencing

Background

Single-cell RNA-sequencing (scRNA-seq) relies on PCR amplification to retrieve information from vanishingly small amounts of starting material. To selectively enrich mRNA from abundant non-polyadenylated transcripts, poly(A) selection is a key step during library preparation. However, some transcripts, such as mitochondrial genes, can escape this elimination and overwhelm libraries. Often, these transcripts are removed in silico, but whether physical depletion improves detection of rare transcripts in single cells is unclear.

Results

We find that a single 16S ribosomal RNA is widely enriched in planarian scRNA-seq datasets, independent of the library preparation method. To deplete this transcript from scRNA-seq libraries, we design 30 single-guide RNAs spanning its length. To evaluate the effects of depletion, we perform a side-by-side comparison of the effects of eliminating the 16S transcript and find a substantial increase in the number of genes detected per cell, coupled with virtually complete loss of the 16S RNA. Moreover, we systematically determine that library complexity increases with a limited number of PCR cycles following CRISPR treatment. When compared to in silico depletion of 16S, physically removing it reduces dropout rates, retrieves more clusters, and reveals more differentially-expressed genes.

Conclusions

Our results show that abundant transcripts reduce the retrieval of informative transcripts in scRNA-seq and distort the analysis. Physical removal of these contaminants enables the detection of rare transcripts at lower sequencing depth, and also outperforms in silico depletion. Importantly, this method can be easily customized to deplete any abundant transcript from scRNA-seq libraries.

Parallels and contrasts between the cnidarian and bilaterian maternal-to-zygotic transition are revealed in Hydractinia embryos

Embryogenesis requires coordinated gene regulatory activities early on that establish the trajectory of subsequent development, during a period called the maternal-to-zygotic transition (MZT). The MZT comprises transcriptional activation of the embryonic genome and post-transcriptional regulation of egg-inherited maternal mRNA. Investigation into the MZT in animals has focused almost exclusively on bilaterians, which include all classical models such as flies, worms, sea urchin, and vertebrates, thus limiting our capacity to understand the gene regulatory paradigms uniting the MZT across all animals. Here, we elucidate the MZT of a non-bilaterian, the cnidarian Hydractinia symbiolongicarpus . Using parallel poly(A)-selected and non poly(A)-dependent RNA-seq approaches, we find that the Hydractinia MZT is composed of regulatory activities analogous to many bilaterians, including cytoplasmic readenylation of maternally contributed mRNA, delayed genome activation, and separate phases of maternal mRNA deadenylation and degradation that likely depend on both maternally and zygotically encoded clearance factors, including microRNAs. But we also observe massive upregulation of histone genes and an expanded repertoire of predicted H4K20 methyltransferases, aspects thus far unique to the Hydractinia MZT and potentially underlying a novel mode of early embryonic chromatin regulation. Thus, similar regulatory strategies with taxon-specific elaboration underlie the MZT in both bilaterian and non-bilaterian embryos, providing insight into how an essential developmental transition may have arisen in ancestral animals.

Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito rRNA datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNA-seq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Efficient and specific DNA oligonucleotide rRNA probe-based rRNA removal in Talaromyces marneffei

Emerging evidence showed that lncRNAs play important roles in a wide range of biological processes of fungi such as Saccharomyces cerevisiae. However, systemic identification of lncRNAs in non-model fungi is a challenging task as the efficiency of rRNA removal has been proved to be affected by mismatches of universal rRNA-targeting probes of commercial kits, which forces deeper sequencing depth and increases costs. Here, we developed a low-cost and simple rRNA depletion method (rProbe) that could efficiently remove more than 99% rRNA in both yeast and mycelium samples of Talaromyces marneffei. The efficiency and robustness of rProbe were demonstrated to outperform the Illumina Ribo-Zero kit. Using rProbe RNA-seq, we identified 115 differentially expressed lncRNAs and constructed lncRNA-mRNA co-expression network related to dimorphic switch of T. marneffei. Our rRNA removal method has the potential to be a useful tool to explore non-coding transcriptomes of non-model fungi by adjusting rRNA probe sequences species specifically.