Selective depletion of rRNA enables whole transcriptome profiling of archival fixed tissue

The impact of microplastics on tissue-specific gene expression in the tropical house cricket, Gryllodes sigillatus

Microplastics are ubiquitous in our environment, resulting in animal exposure and consumption via food, water, and air. Animals that consume microplastics may suffer from physiological effects like immunotoxicity or mitochondrial dysfunction, but how specific tissues may differentially respond to plastic consumption is poorly understood, particularly in terrestrial insects. Here, we measured transcriptomic responses of tissues (midgut, hindgut, fat body and ovaries) to microplastic consumption in a generalist ground-dwelling insect, the tropical house cricket, Gryllodes sigillatus. Using this approach, we provide insights on how microplastics may impact specific organ systems. We generated a de novo transcriptome, a useful resource for further studies on this emerging model insect, that we then used to infer differential gene expression due to microplastic consumption in individual organs. Ingestion of microplastics elicited unique changes in gene expression depending on the tissue of focus, with notable differentially-expressed genes related to survival and stress pathways as well as those related to metabolism, immunity, and cancer.

Direct mRNA-to-sgRNA conversion generates design-free ultra-dense CRISPRi libraries for systematic phenotypic screening

CRISPR interference (CRISPRi) is a versatile tool for high-throughput phenotypic screening. However, rational design and synthesis of the single-guide RNA (sgRNA) library required for each genome-wide CRISPRi application is time-consuming, expensive, and unfeasible if the target organisms lack comprehensive sequencing and characterization. We developed an ultra-dense random sgRNA library generation method applicable to any organism, including those that are not well-characterized. Our method converts transcriptome-wide mRNA into 20 nt of sgRNA spacer sequences through enzymatic reactions. The generated sgRNA library selectively binds to the non-template strand of the coding sequence, leading to more efficient repression compared to binding the template strand. We then generated a genome-scale library for Escherichia coli by applying this method and identified essential and auxotrophic genes through phenotypic screening. Furthermore, we tuned the production levels of lycopene and violacein and identified new repression targets for violacein production. Our results demonstrated that a genome-scale sgRNA library can be generated without rational design and can be utilized simultaneously in a range of phenotypic screenings.

Mosquito-Based Detection of Endogenous Jaagsiekte Sheep Retrovirus in Senegal: Expanding the Scope of Xenosurveillance

Background

Mosquitoes are well-known vectors for arthropod-borne viruses, yet their role as passive carriers of non-arthropod-borne viruses remains underexplored. Xenosurveillance, a method that utilizes blood-feeding arthropods to sample host and pathogen genetic material, has emerged as a valuable tool in viral ecology. In this study, we report the first identification of Jaagsiekte Sheep Retrovirus (JSRV)-related sequences in blood-fed mosquitoes collected in Senegal. JSRV, a betaretrovirus responsible for ovine pulmonary adenocarcinoma, is typically found in sheep, but its genetic trace in mosquitoes offers a novel perspective on host-vector contact and surveillance. Our study aimed to investigate whether mosquitoes can serve as sentinels for detecting both pathogens and host-derived markers in complex ecosystems.

Methods

Mosquitoes were collected between 2016 and 2019 from three ecologically significant regions in Senegal (Louga, Barkedji, and Kedougou). Blood-fed mosquitoes were pooled and subjected to RNA extraction and metagenomic sequencing using Illumina NextSeq550. Sequencing data were analyzed with CZ-ID and BLAST for viral identification. RT-qPCR assays were designed to validate the presence of JSRV-related sequences, targeting conserved regions of the envelope gene and 3' untranslated region. Phylogenetic analysis was conducted using MAFFT and IQ-TREE to compare the detected sequence with global exogenous and endogenous JSRV references.

Results

A diverse array of viruses across mosquito species, including both arboviruses and non-arthropod-borne viruses. A JSRV-related sequence was detected in a single blood-fed mosquito pool collected in Barkedji (2019). The RT-qPCR assay confirmed JSRV presence, validating the sequencing results. Phylogenetic analysis revealed strong similarity to known endogenous JSRV (enJSRV) sequences integrated in the sheep genome, indicating that the detected material likely originated from host DNA ingested during blood feeding.

Discussion

This study presents the first report of endogenous retroviral sequences detected in mosquitoes, alongside the identification of actively circulating viruses, highlighting the broader potential of mosquitoes as environmental sentinels. While mosquitoes are not biological vectors for JSRV, their ability to capture both host-derived retroviral material and pathogenic viral genomes through bloodmeals reinforces the value of xenosurveillance for monitoring livestock-vector-environment interactions. These findings contribute to broader efforts in integrated disease surveillance and underscore the utility of combining metagenomics with molecular diagnostics to detect diverse viral signals in high-risk ecological settings.

Chloroplast precursor protein preClpD overaccumulation triggers multilevel reprogramming of gene expression and a heat shock-like response

Thousands of nucleus-encoded chloroplast proteins are synthesized as precursors on cytosolic ribosomes and posttranslationally imported into chloroplasts. Cytosolic accumulation of unfolded chloroplast precursor proteins (e.g., under stress conditions) is hazardous to the cell. The global cellular responses and regulatory pathways involved in triggering appropriate responses are largely unknown. Here, by inducible and constitutive overexpression of ClpD-GFP to result in precursor protein overaccumulation, we present a comprehensive picture of multilevel reprogramming of gene expression in response to chloroplast precursor overaccumulation stress (cPOS), reveal a critical role of translational activation in the expression of cytosolic chaperones (heat-shock proteins, HSPs), and demonstrate that chloroplast-derived reactive oxygen species act as retrograde signal for the transcriptional activation of small HSPs. Furthermore, we reveal an important role of the chaperone ClpB1/HOT1 in maintaining cellular proteostasis upon cPOS. Together, our observations uncover a cytosolic heat shock-like response to cPOS and provide insights into the underlying molecular mechanisms.

A Csde1-Strap complex regulates plasma cell differentiation by coupling mRNA translation and decay

Upon encountering antigens, B cells may undergo multiple differentiation paths, including becoming plasma cells and memory B cells. Although it is well-known that transcription factors govern gene expression programs underpinning these fate decisions in transcriptional level, the role of post-transcriptional regulators, with a focus on RNA-binding proteins, in the fate determination are lesser known. Here we find by RNA interactome capture-coupled CRISPR/Cas9 functional screening that the Csde1-Strap complex plays an important role in plasma cell differentiation. Mechanistically, the Csde1-Strap complex establishes the expression kinetics of Bach2, a key regulator of plasma cell differentiation. Bach2 expression is rapidly induced to promote B cell expansion and then decreased to initiate plasma cell differentiation. The Csde1-Strap interaction is critical for their binding to Bach2 mRNA to couple its decay with translation to restrain the magnitude and duration of Bach2 protein expression. In the absence of Csde1 or Strap, Bach2 translation is de-coupled from mRNA decay, leading to elevated and prolonged expression of Bach2 protein and impaired plasma cell differentiation. This study thus establishes the functional RBP landscape in B cells and illustrates the fundamental importance of controlling protein expression kinetics in cell fate determination.

Impact of commercial RNA extraction methods on the recovery of human RNA sequence data from archival fixed tissues

Archival fixed tissues hold key insights into the evolutionary history of RNA viruses and the associated host immune response, yet access to the RNA sequence data is limited by a lack of robust methods for RNA extraction and sequence retrieval from these tissue types. Here we compared three commercial RNA extraction techniques (bead, column, and phase-based) on five fixed human brain tissues done in triplicate, that have been stored for up to 43 years. We found that for this sample set, bead-based extractions captured longer molecules and yielded a greater proportion of unique reads when aligned to the human genome, than did column and phase-based extraction methods. Via the incorporation of multiple extraction replicates, we quantified the variability in sequencing metrics resulting from tissue sample and extraction technique heterogeneity. Additionally, we compared pre- and post-sequencing metrics and found that the former poorly predicted post-sequencing on-target success. Our findings help inform future research on the recovery of RNA from archival fixed tissues.

A probe-based capture enrichment method for detection of A-to-I editing in low abundance transcripts

Exactly two decades ago, the ability to use high-throughput RNA sequencing technology to identify sites of editing by ADARs was employed for the first time. Since that time, RNA sequencing has become a standard tool for researchers studying RNA biology and led to the discovery of RNA editing sites present in a multitude of organisms, across tissue types, and in disease. However, transcriptome-wide sequencing is not without limitations. Most notably, RNA sequencing depth of a given transcript is correlated with expression, and sequencing depth impacts the ability to robustly detect RNA editing events. This chapter focuses on a method for enrichment of low-abundance transcripts that can facilitate more efficient sequencing and detection of RNA editing events. An important note is that while we describe aspects of the protocol important for capturing intron-containing transcripts, this probe-based enrichment method could be easily modified to assess editing within any low-abundance transcript. We also provide some perspectives on the current limitations as well as important future directions for expanding this technology to gain more insights into how RNA editing can impact transcript diversity.

Ribosomal RNA Depletion for Poly(A)-Tail-Independent Quantification of Genome Activation

High-throughput RNA sequencing (RNA-seq) is commonly used to quantify gene expression transcriptome-wide. While usually paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA) to the exclusion of highly abundant ribosomal RNA (rRNA) in the cell, this strategy will under-quantify mRNA with short or absent poly(A) tails and can conflate changes in poly(A) tail length with changes in RNA level. This is notably an issue during early development, when cytoplasmic polyadenylation of maternal mRNA over time can be mistaken for genome activation in poly(A) + RNA-seq time courses. Here, we present a method to perform total RNA-seq using a streamlined rRNA depletion strategy customizable to any taxon. Antisense DNA oligos are designed with the aid of our Oligo-ASST web tool to sparsely tile the length of the rRNA, which are used with thermostable RNaseH to digest rRNA from a total RNA sample. After column cleanup, the mRNA-enriched sample is ready for sequencing library construction.

Adaptable and comprehensive approaches for long-read nanopore sequencing of polyadenylated and non-polyadenylated RNAs

The advent of long-read (LR) sequencing technologies has provided a direct opportunity to determine the structure of transcripts with potential for end-to-end sequencing of full-length RNAs. LR methods that have been described to date include commercial offerings from Oxford Nanopore Technologies (ONT) and Pacific Biosciences. These kits are based on selection of polyadenylated (polyA+) RNAs and/or oligo-dT priming of reverse transcription. Thus, these approaches do not allow comprehensive interrogation of the transcriptome due to their exclusion of non-polyadenylated (polyA-) RNAs. In addition, polyA + specificity also results in 3'-biased measurements of PolyA+ RNAs especially when the RNA input is partially degraded. To address these limitations of current LR protocols, we modified rRNA depletion protocols that have been used in short-read sequencing: one approach representing a ligation-based method and the other a template-switch cDNA synthesis-based method to append ONT-specific adaptor sequences and by removing any deliberate fragmentation/shearing of RNA/cDNA. Here, we present comparisons with poly+ RNA-specific versions of the two approaches including the ONT PCR-cDNA Barcoding kit. The rRNA depletion protocols displayed higher proportions (30%-50%) of intronic content compared to that of the polyA-specific protocols (5%-8%). In addition, the rRNA depletion protocols enabled ∼20-50% higher detection of expressed genes. Other metrics that were favourable to the rRNA depletion protocols include better coverage of long transcripts, and higher accuracy and reproducibility of expression measurements. Overall, these results indicate that the rRNA depletion-based protocols described here allow the comprehensive characterization of polyadenylated and non-polyadenylated RNAs. While the resulting reads are long enough to help decipher transcript structures, future endeavors are warranted to improve the proportion of individual reads representing end-to-end spanning of transcripts.

A novel N4,N4-dimethylcytidine in the archaeal ribosome enhances hyperthermophily

Ribosome structure and activity are challenged at high temperatures, often demanding modifications to ribosomal RNAs (rRNAs) to retain translation fidelity. LC-MS/MS, bisulfite-sequencing, and high-resolution cryo-EM structures of the archaeal ribosome identified an RNA modification, N4,N4-dimethylcytidine (m42C), at the universally conserved C918 in the 16S rRNA helix 31 loop. Here, we characterize and structurally resolve a class of RNA methyltransferase that generates m42C whose function is critical for hyperthermophilic growth. m42C is synthesized by the activity of a unique family of RNA methyltransferase containing a Rossman-fold that targets only intact ribosomes. The phylogenetic distribution of the newly identified m42C synthase family implies that m42C is biologically relevant in each domain. Resistance of m42C to bisulfite-driven deamination suggests that efforts to capture m5C profiles via bisulfite sequencing are also capturing m42C.

A germline PAF1 paralog complex ensures cell type-specific gene expression

Animal germline development and fertility rely on paralogs of general transcription factors that recruit RNA polymerase II to ensure cell type-specific gene expression. It remains unclear whether gene expression processes downstream from such paralog-based transcription is distinct from that of canonical RNA polymerase II genes. In Drosophila, the testis-specific TBP-associated factors (tTAFs) activate over a thousand spermatocyte-specific gene promoters to enable meiosis and germ cell differentiation. Here, we show that efficient termination of tTAF-activated transcription relies on testis-specific paralogs of canonical polymerase-associated factor 1 complex (PAF1C) proteins, which form a testis-specific PAF1C (tPAF). Consequently, tPAF mutants show aberrant expression of hundreds of downstream genes due to read-in transcription. Furthermore, tPAF facilitates expression of Y-linked male fertility factor genes and thus serves to maintain spermatocyte-specific gene expression. Consistently, tPAF is required for the segregation of meiotic chromosomes and male fertility. Supported by comparative in vivo protein interaction assays, we provide a mechanistic model for the functional divergence of tPAF and the PAF1C and identify transcription termination as a developmentally regulated process required for germline-specific gene expression.

The extensive m5C epitranscriptome of Thermococcus kodakarensis is generated by a suite of RNA methyltransferases that support thermophily

RNAs are often modified to invoke new activities. While many modifications are limited in frequency, restricted to non-coding RNAs, or present only in select organisms, 5-methylcytidine (m5C) is abundant across diverse RNAs and fitness-relevant across Domains of life, but the synthesis and impacts of m5C have yet to be fully investigated. Here, we map m5C in the model hyperthermophile, Thermococcus kodakarensis. We demonstrate that m5C is ~25x more abundant in T. kodakarensis than human cells, and the m5C epitranscriptome includes ~10% of unique transcripts. T. kodakarensis rRNAs harbor tenfold more m5C compared to Eukarya or Bacteria. We identify at least five RNA m5C methyltransferases (R5CMTs), and strains deleted for individual R5CMTs lack site-specific m5C modifications that limit hyperthermophilic growth. We show that m5C is likely generated through partial redundancy in target sites among R5CMTs. The complexity of the m5C epitranscriptome in T. kodakarensis argues that m5C supports life in the extremes.

CleanUpRNAseq: An R/Bioconductor Package for Detecting and Correcting DNA Contamination in RNA-Seq Data

RNA sequencing (RNA-seq) has become a standard method for profiling gene expression, yet genomic DNA (gDNA) contamination carried over to the sequencing library poses a significant challenge to data integrity. Detecting and correcting this contamination is vital for accurate downstream analyses. Particularly, when RNA samples are scarce and invaluable, it becomes essential not only to identify but also to correct gDNA contamination to maximize the data's utility. However, existing tools capable of correcting gDNA contamination are limited and lack thorough evaluation. To fill the gap, we developed CleanUpRNAseq, which offers a comprehensive set of functionalities for identifying and correcting gDNA-contaminated RNA-seq data. Our package offers three correction methods for unstranded RNA-seq data and a dedicated approach for stranded data. Through rigorous validation on published RNA-seq datasets with known levels of gDNA contamination and real-world RNA-seq data, we demonstrate CleanUpRNAseq's efficacy in detecting and correcting detrimental levels of gDNA contamination across diverse library protocols. CleanUpRNAseq thus serves as a valuable tool for post-alignment quality assessment of RNA-seq data and should be integrated into routine workflows for RNA-seq data analysis. Its incorporation into OneStopRNAseq should significantly bolster the accuracy of gene expression quantification and differential expression analysis of RNA-seq data.

Dissecting infectious bronchitis virus-induced host shutoff at the translation level

Viruses have evolved a range of strategies to utilize or manipulate the host's cellular translational machinery for efficient infection, although the mechanisms by which infectious bronchitis virus (IBV) manipulates the host translation machinery remain unclear. In this study, we firstly demonstrate that IBV infection causes host shutoff, although viral protein synthesis is not affected. We then screened 23 viral proteins, and identified that more than one viral protein is responsible for IBV-induced host shutoff, the inhibitory effects of proteins Nsp15 were particularly pronounced. Ribosome profiling was used to draw the landscape of viral mRNA and cellular genes expression model, and the results showed that IBV mRNAs gradually dominated the cellular mRNA pool, the translation efficiency of the viral mRNAs was lower than the median efficiency (about 1) of cellular mRNAs. In the analysis of viral transcription and translation, higher densities of RNA sequencing (RNA-seq) and ribosome profiling (Ribo-seq) reads were observed for structural proteins and 5' untranslated regions, which conformed to the typical transcriptional characteristics of nested viruses. Translational halt events and the number of host genes increased significantly after viral infection. The translationally paused genes were enriched in translation, unfolded-protein-related response, and activation of immune response pathways. Immune- and inflammation-related mRNAs were inefficiently translated in infected cells, and IBV infection delayed the production of IFN-β and IFN-λ. Our results describe the translational landscape of IBV-infected cells and demonstrate new strategies by which IBV induces host gene shutoff to promote its replication.

IMPORTANCE

Infectious bronchitis virus (IBV) is a γ-coronavirus that causes huge economic losses to the poultry industry. Understanding how the virus manipulates cellular biological processes to facilitate its replication is critical for controlling viral infections. Here, we used Ribo-seq to determine how IBV infection remodels the host's biological processes and identified multiple viral proteins involved in host gene shutoff. Immune- and inflammation-related mRNAs were inefficiently translated, the translation halt of unfolded proteins and immune activation-related genes increased significantly, benefitting IBV replication. These data provide new insights into how IBV modulates its host's antiviral responses.

A hidden proteome encoded by circRNAs in human placentas: Implications for uncovering preeclampsia pathogenesis

BACKGROUND

CircRNA-encoded proteins (CEPs) are emerging as new players in health and disease, and function as baits for the common partners of their cognate linear-spliced RNA encoded proteins (LEPs). However, their prevalence across human tissues and biological roles remain largely unexplored. The placenta is an ideal model for identifying CEPs due to its considerable protein diversity that is required to sustain fetal development during pregnancy. The aim of this study was to evaluate circRNA translation in the human placenta, and the potential roles of the CEPs in placental development and dysfunction.

METHODS

Multiomics approaches, including RNA sequencing, ribosome profiling, and LC-MS/MS analysis, were utilised to identify novel translational events of circRNAs in human placentas. Bioinformatics methods and the protein bait hypothesis were employed to evaluate the roles of these newly discovered CEPs in placentation and associated disorders. The pathogenic role of a recently identified CEP circPRKCB119aa in preeclampsia was investigated through qRT-PCR, Western blotting, immunofluorescence imaging and phenotypic analyses.

RESULTS

We found that 528 placental circRNAs bound to ribosomes with active translational elongation, and 139 were translated to proteins. The CEPs showed considerable structural homology with their cognate LEPs, but are more stable, hydrophobic and have a lower molecular-weight than the latter, all of which are conducive to their function as baits. On this basis, CEPs are deduced to be closely involved in placental function. Furthermore, we focused on a novel CEP circPRKCB119aa, and illuminated its pathogenic role in preeclampsia; it enhanced trophoblast autophagy by acting as a bait to inhibit phosphorylation of the cognate linear isoform PKCβ.

CONCLUSIONS

We discovered a hidden circRNA-encoded proteome in the human placenta, which offers new insights into the mechanisms underlying placental development, as well as placental disorders such as preeclampsia. Key points A hidden circRNA-encoded proteome in the human placenta was extensively identified and systematically characterised. The circRNA-encoded proteins (CEPs) are potentially related to placental development and associated disorders. A novel conserved CEP circPRKCB119aa enhanced trophoblast autophagy by inhibiting phosphorylation of its cognate linear-spliced isoform protein kinase C (PKC) β in preeclampsia.

Comparison of RNA-Sequencing Methods for Degraded RNA

RNA sequencing (RNA-Seq) is a powerful technique and is increasingly being used in clinical research and drug development. Currently, several RNA-Seq methods have been developed. However, the relative advantage of each method for degraded RNA and low-input RNA, such as RNA samples collected in the field of clinical setting, has remained unknown. The Standard method of RNA-Seq captures mRNA by poly(A) capturing using Oligo dT beads, which is not suitable for degraded RNA. Here, we used three commercially available RNA-Seq library preparation kits (SMART-Seq, xGen Broad-range, and RamDA-Seq) using random primer instead of Oligo dT beads. To evaluate the performance of these methods, we compared the correlation, the number of detected expressing genes, and the expression levels with the Standard RNA-Seq method. Although the performance of RamDA-Seq was similar to that of Standard RNA-Seq, the performance for low-input RNA and degraded RNA has decreased. The performance of SMART-Seq was better than xGen and RamDA-Seq in low-input RNA and degraded RNA. Furthermore, the depletion of ribosomal RNA (rRNA) improved the performance of SMART-Seq and xGen due to increased expression levels. SMART-Seq with rRNA depletion has relative advantages for RNA-Seq using low-input and degraded RNA.

A scalable and cost-efficient rRNA depletion approach to enrich RNAs for molecular biology investigations

Transcriptomics analyses play pivotal roles in understanding the complex regulatory networks that govern cellular processes. The abundance of rRNAs, which account for 80%-90% of total RNA in eukaryotes, limits the detection and investigation of other transcripts. While mRNAs and long noncoding RNAs have poly(A) tails that are often used for positive selection, investigations of poly(A)- RNAs, such as circular RNAs, histone mRNAs, and small RNAs, typically require the removal of the abundant rRNAs for enrichment. Current approaches to deplete rRNAs for downstream molecular biology investigations are hampered by restrictive RNA input masses and high costs. To address these challenges, we developed rRNA Removal by RNaseH (rRRR), a method to efficiently deplete rRNAs from a wide range of human, mouse, and rat RNA inputs and of varying qualities at a cost 10- to 20-fold cheaper than other approaches. We used probe-based hybridization and enzymatic digestion to selectively target and remove rRNA molecules while preserving the integrity of non-rRNA transcripts. Comparison of rRRR to two commercially available approaches showed similar rRNA depletion efficiencies and comparable off-target effects. Our developed method provides researchers with a valuable tool for investigating gene expression and regulatory mechanisms across a wide range of biological systems at an affordable price that increases the accessibility for researchers to enter the field, ultimately advancing our understanding of cellular processes.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Thylakoid protein FPB1 synergistically cooperates with PAM68 to promote CP47 biogenesis and Photosystem II assembly

In chloroplasts, insertion of proteins with multiple transmembrane domains (TMDs) into thylakoid membranes usually occurs in a co-translational manner. Here, we have characterized a thylakoid protein designated FPB1 (Facilitator of PsbB biogenesis1) which together with a previously reported factor PAM68 (Photosynthesis Affected Mutant68) is involved in assisting the biogenesis of CP47, a subunit of the Photosystem II (PSII) core. Analysis by ribosome profiling reveals increased ribosome stalling when the last TMD segment of CP47 emerges from the ribosomal tunnel in fpb1 and pam68. FPB1 interacts with PAM68 and both proteins coimmunoprecipitate with SecY/E and Alb3 as well as with some ribosomal components. Thus, our data indicate that, in coordination with the SecY/E translocon and the Alb3 integrase, FPB1 synergistically cooperates with PAM68 to facilitate the co-translational integration of the last two CP47 TMDs and the large loop between them into thylakoids and the PSII core complex.

AGO2 Protects Against Diabetic Cardiomyopathy by Activating Mitochondrial Gene Translation

BACKGROUND

Diabetes is associated with cardiovascular complications. microRNAs translocate into subcellular organelles to modify genes involved in diabetic cardiomyopathy. However, functional properties of subcellular AGO2 (Argonaute2), a core member of miRNA machinery, remain elusive.

METHODS

We elucidated the function and mechanism of subcellular localized AGO2 on mouse models for diabetes and diabetic cardiomyopathy. Recombinant adeno-associated virus type 9 was used to deliver AGO2 to mice through the tail vein. Cardiac structure and functions were assessed by echocardiography and catheter manometer system.

RESULTS

AGO2 was decreased in mitochondria of diabetic cardiomyocytes. Overexpression of mitochondrial AGO2 attenuated diabetes-induced cardiac dysfunction. AGO2 recruited TUFM, a mitochondria translation elongation factor, to activate translation of electron transport chain subunits and decrease reactive oxygen species. Malonylation, a posttranslational modification of AGO2, reduced the importing of AGO2 into mitochondria in diabetic cardiomyopathy. AGO2 malonylation was regulated by a cytoplasmic-localized short isoform of SIRT3 through a previously unknown demalonylase function.

CONCLUSIONS

Our findings reveal that the SIRT3-AGO2-CYTB axis links glucotoxicity to cardiac electron transport chain imbalance, providing new mechanistic insights and the basis to develop mitochondria targeting therapies for diabetic cardiomyopathy.

Two distinct waves of transcriptome and translatome changes drive Drosophila germline stem cell differentiation

The tight control of fate transitions during stem cell differentiation is essential for proper tissue development and maintenance. However, the challenges in studying sparsely distributed adult stem cells in a systematic manner have hindered efforts to identify how the multilayered regulation of gene expression programs orchestrates stem cell differentiation in vivo. Here, we synchronised Drosophila female germline stem cell (GSC) differentiation in vivo to perform in-depth transcriptome and translatome analyses at high temporal resolution. This characterisation revealed widespread and dynamic changes in mRNA level, promoter usage, exon inclusion, and translation efficiency. Transient expression of the master regulator, Bam, drives a first wave of expression changes, primarily modifying the cell cycle program. Surprisingly, as Bam levels recede, differentiating cells return to a remarkably stem cell-like transcription and translation program, with a few crucial changes feeding into a second phase driving terminal differentiation to form the oocyte. Altogether, these findings reveal that rather than a unidirectional accumulation of changes, the in vivo differentiation of stem cells relies on distinctly regulated and developmentally sequential waves.

NAP-seq reveals multiple classes of structured noncoding RNAs with regulatory functions

Up to 80% of the human genome produces "dark matter" RNAs, most of which are noncapped RNAs (napRNAs) that frequently act as noncoding RNAs (ncRNAs) to modulate gene expression. Here, by developing a method, NAP-seq, to globally profile the full-length sequences of napRNAs with various terminal modifications at single-nucleotide resolution, we reveal diverse classes of structured ncRNAs. We discover stably expressed linear intron RNAs (sliRNAs), a class of snoRNA-intron RNAs (snotrons), a class of RNAs embedded in miRNA spacers (misRNAs) and thousands of previously uncharacterized structured napRNAs in humans and mice. These napRNAs undergo dynamic changes in response to various stimuli and differentiation stages. Importantly, we show that a structured napRNA regulates myoblast differentiation and a napRNA DINAP interacts with dyskerin pseudouridine synthase 1 (DKC1) to promote cell proliferation by maintaining DKC1 protein stability. Our approach establishes a paradigm for discovering various classes of ncRNAs with regulatory functions.

High expression circRALGPS2 in atretic follicle induces chicken granulosa cell apoptosis and autophagy via encoding a new protein

BACKGROUND

The reproductive performance of chickens mainly depends on the development of follicles. Abnormal follicle development can lead to decreased reproductive performance and even ovarian disease among chickens. Chicken is the only non-human animal with a high incidence of spontaneous ovarian cancer. In recent years, the involvement of circRNAs in follicle development and atresia regulation has been confirmed.

RESULTS

In the present study, we used healthy and atretic chicken follicles for circRNA RNC-seq. The results showed differential expression of circRALGPS2. It was then confirmed that circRALGPS2 can translate into a protein, named circRALGPS2-212aa, which has IRES activity. Next, we found that circRALGPS2-212aa promotes apoptosis and autophagy in chicken granulosa cells by forming a complex with PARP1 and HMGB1.

CONCLUSIONS

Our results revealed that circRALGPS2 can regulate chicken granulosa cell apoptosis and autophagy through the circRALGPS2-212aa/PARP1/HMGB1 axis.

TREX reveals proteins that bind to specific RNA regions in living cells

Different regions of RNA molecules can often engage in specific interactions with distinct RNA-binding proteins (RBPs), giving rise to diverse modalities of RNA regulation and function. However, there are currently no methods for unbiased identification of RBPs that interact with specific RNA regions in living cells and under endogenous settings. Here we introduce TREX (targeted RNase H-mediated extraction of crosslinked RBPs)-a highly sensitive approach for identifying proteins that directly bind to specific RNA regions in living cells. We demonstrate that TREX outperforms existing methods in identifying known interactors of U1 snRNA, and reveals endogenous region-specific interactors of NORAD long noncoding RNA. Using TREX, we generated a comprehensive region-by-region interactome for 45S rRNA, uncovering both established and previously unknown interactions that regulate ribosome biogenesis. With its applicability to different cell types, TREX is an RNA-centric tool for unbiased positional mapping of endogenous RNA-protein interactions in living cells.

Spitting in the wind?-The challenges of RNA sequencing for biomarker discovery from saliva

Forensic trace contextualization, i.e., assessing information beyond who deposited a biological stain, has become an issue of great and steadily growing importance in forensic genetic casework and research. The human transcriptome encodes a wide variety of information and thus has received increasing interest for the identification of biomarkers for different aspects of forensic trace contextualization over the past years. Massively parallel sequencing of reverse-transcribed RNA ("RNA sequencing") has emerged as the gold standard technology to characterize the transcriptome in its entirety and identify RNA markers showing significant expression differences not only between different forensically relevant body fluids but also within a single body fluid between forensically relevant conditions of interest. Here, we analyze the quality and composition of four RNA sequencing datasets (whole transcriptome as well as miRNA sequencing) from two different research projects (the RNAgE project and the TrACES project), aiming at identifying contextualizing forensic biomarker from the forensically relevant body fluid saliva. We describe and characterize challenges of RNA sequencing of saliva samples arising from the presence of oral bacteria, the heterogeneity of sample composition, and the confounding factor of degradation. Based on these observations, we formulate recommendations that might help to improve RNA biomarker discovery from the challenging but forensically relevant body fluid saliva.

The 28S rRNA RT-qPCR assay for host depletion evaluation to enhance avian virus detection in Illumina and Nanopore sequencing

Abundant host and bacterial sequences can obscure the detection of less prevalent viruses in untargeted next-generation sequencing (NGS). Efficient removal of these non-targeted sequences is vital for accurate viral detection. This study presents a novel 28S ribosomal RNA (rRNA) RT-qPCR assay designed to assess the efficiency of avian rRNA depletion before conducting costly NGS for the detection of avian RNA viruses. The comprehensive evaluation of this 28S-test focuses on substituting DNase I with alternative DNases in our established depletion protocols and finetuning essential parameters for reliable host rRNA depletion. To validate the effectiveness of the 28S-test, we compared its performance with NGS results obtained from both Illumina and Nanopore sequencing platforms. This evaluation utilized swab samples from chickens infected with highly pathogenic avian influenza virus, subjected to established and modified depletion protocols. Both methods significantly reduced host rRNA levels, but using the alternative DNase had superior performance. Additionally, utilizing the 28S-test, we explored cost- and time-effective strategies, such as reduced probe concentrations and other alternative DNase usage, assessed the impact of filtration pre-treatment, and evaluated various experimental parameters to further optimize the depletion protocol. Our findings underscore the value of the 28S-test in optimizing depletion methods for advancing improvements in avian disease research through NGS.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

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 IGF2BP3-COPS7B Axis Facilitates mRNA Translation to Drive Colorectal Cancer Progression

Many studies have provided valuable information about genomic and transcriptomic changes that occur in colorectal cancer. However, protein abundance cannot be reliably predicted by DNA alteration or mRNA expression, which can be partially attributed to posttranscriptional and/or translational regulation of gene expression. In this study, we identified increased translational efficiency (TE) as a hallmark of colorectal cancer by evaluating the transcriptomic and proteomic features of patients with colorectal cancer, along with comparative transcriptomic and ribosome-protected mRNA analysis in colon epithelial cells and colon cancer cells. COP9 signalosome subunit 7B (COPS7B) was among the key genes that consistently showed both significant TE increase and protein elevation without transcriptional alteration in colorectal cancer. Insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) enhanced the TE of COPS7B mRNA to promote colorectal cancer growth and metastasis. COPS7B was found to be a component of the ribo-interactome that interacted with ribosomes to facilitate ribosome biogenesis and mRNA translation initiation. Collectively, this study revealed the proteomic features of colorectal cancer and highlighted elevated mRNA translation as a hallmark of colorectal cancer. The identification of the IGF2BP3-COPS7B axis underlying the increased protein synthesis rate in colorectal cancer provided a promising therapeutic target to treat this aggressive disease.

SIGNIFICANCE

Increased expression of COPS7B mediated by IGF2BP3 elevates the translational efficiency of genes enriched in mRNA translation and ribosome biogenesis pathways, promoting protein synthesis and driving progression in colorectal cancer.

Alternative probe hybridization buffers for target RNA depletion and viral sequence recovery in NGS for poultry samples

Non-targeted next generation sequencing (NGS) is widely applied to identify the diversity of pathogens in field samples. However, abundance of host RNA (especially rRNA) and other environmental nucleic acids can reduce the abundance of pathogen specific reads of interest, reduce depth of coverage and increase surveillance costs. We presently deplete chicken- and selected bacterial-specific rRNAs in poultry field RNA samples with complementary DNA probes in a commercially available probe hybridization buffer followed by digestion of the RNA:DNA hybrids with RNase H. Because the current buffer is an expensive special order reagent of proprietary composition, we tested in-house and other commercially available buffers and identified a viable alternative that yields equivalent host rRNA depletion and viral-specific reads in poultry samples as the current special order reagent but at a reduced cost.

Multi-omic approach identifies hypoxic tumor-associated myeloid cells that drive immunobiology of high-risk pediatric ependymoma

Ependymoma (EPN) is a devastating childhood brain tumor. Single-cell analyses have illustrated the cellular heterogeneity of EPN tumors, identifying multiple neoplastic cell states including a mesenchymal-differentiated subpopulation which characterizes the PFA1 subtype. Here, we characterize the EPN immune environment, in the context of both tumor subtypes and tumor cell subpopulations using single-cell sequencing (scRNAseq, n = 27), deconvolution of bulk tumor gene expression (n = 299), spatial proteomics (n = 54), and single-cell cytokine release assays (n = 12). We identify eight distinct myeloid-derived subpopulations from which a group of cells, termed hypoxia myeloid cells, demonstrate features of myeloid-derived suppressor cells, including IL6/STAT3 pathway activation and wound healing ontologies. In PFA tumors, hypoxia myeloid cells colocalize with mesenchymal-differentiated cells in necrotic and perivascular niches and secrete IL-8, which we hypothesize amplifies the EPN immunosuppressive microenvironment. This myeloid cell-driven immunosuppression will need to be targeted for immunotherapy to be effective in this difficult-to-cure childhood brain tumor.

Acute and Long-Term Consequences of Co-opted doublesex on the Development of Mimetic Butterfly Color Patterns

Novel phenotypes are increasingly recognized to have evolved by co-option of conserved genes into new developmental contexts, yet the process by which co-opted genes modify existing developmental programs remains obscure. Here, we provide insight into this process by characterizing the role of co-opted doublesex in butterfly wing color pattern development. dsx is the master regulator of insect sex differentiation but has been co-opted to control the switch between discrete nonmimetic and mimetic patterns in Papilio alphenor and its relatives through the evolution of novel mimetic alleles. We found dynamic spatial and temporal expression pattern differences between mimetic and nonmimetic butterflies throughout wing development. A mimetic color pattern program is switched on by a pulse of dsx expression in early pupal development that causes acute and long-term differential gene expression, particularly in Wnt and Hedgehog signaling pathways. RNAi suggested opposing, novel roles for these pathways in mimetic pattern development. Importantly, Dsx co-option caused Engrailed, a primary target of Hedgehog signaling, to gain a novel expression domain early in pupal wing development that is propagated through mid-pupal development to specify novel mimetic patterns despite becoming decoupled from Dsx expression itself. Altogether, our findings provide multiple views into how co-opted genes can both cause and elicit changes to conserved networks and pathways to result in development of novel, adaptive phenotypes.

FACT regulates pluripotency through proximal and distal regulation of gene expression in murine embryonic stem cells

BACKGROUND

The FACT complex is a conserved histone chaperone with critical roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear.

RESULTS

To explore the mechanism for FACT function, we generated AID-tagged murine embryonic cell lines for FACT subunit SPT16 and paired depletion with nascent transcription and chromatin accessibility analyses. We also analyzed SPT16 occupancy using CUT&RUN and found that SPT16 localizes to both promoter and enhancer elements, with a strong overlap in binding with OCT4, SOX2, and NANOG. Over a timecourse of SPT16 depletion, nucleosomes invade new loci, including promoters, regions bound by SPT16, OCT4, SOX2, and NANOG, and TSS-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (encoding OCT4), Sox2, Nanog, and enhancer RNAs produced from these genes' associated enhancers are downregulated.

CONCLUSIONS

We propose that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency.

Genome Organization of Four Brazilian Xanthomonas albilineans Strains Does Not Correlate with Aggressiveness

An integrative approach combining genomics, transcriptomics, and cell biology is presented to address leaf scald disease, a major problem for the sugarcane industry. To gain insight into the biology of the causal agent, the complete genome sequences of four Brazilian Xanthomonas albilineans strains with differing virulence capabilities are presented and compared to the GPEPC73 reference strain and FJ1. Based on the aggressiveness index, different strains were compared: Xa04 and Xa11 are highly aggressive, Xa26 is intermediate, and Xa21 is the least, while, based on genome structure, Xa04 shares most of its genomic features with Xa26, and Xa11 share most of its genomic features with Xa21. In addition to presenting more clustered regularly interspaced short palindromic repeats (CRISPR) clusters, four more novel prophage insertions are present than the previously sequenced GPEPC73 and FJ1 strains. Incorporating the aggressiveness index and in vitro cell biology into these genome features indicates that disease establishment is not a result of a single determinant factor, as in most other Xanthomonas species. The Brazilian strains lack the previously described plasmids but present more prophage regions. In pairs, the most virulent and the least virulent share unique prophages. In vitro transcriptomics shed light on the 54 most highly expressed genes among the 4 strains compared to ribosomal proteins (RPs), of these, 3 outer membrane proteins. Finally, comparative albicidin inhibition rings and in vitro growth curves of the four strains also do not correlate with pathogenicity. In conclusion, the results disclose that leaf scald disease is not associated with a single shared characteristic between the most or the least pathogenic strains. IMPORTANCE An integrative approach is presented which combines genomics, transcriptomics, and cell biology to address leaf scald disease. The results presented here disclose that the disease is not associated with a single shared characteristic between the most pathogenic strains or a unique genomic pattern. Sequence data from four Brazilian strains are presented that differ in pathogenicity index: Xa04 and Xa11 are highly virulent, Xa26 is intermediate, and Xa21 is the least pathogenic strain, while, based on genome structure, Xa04 shares with Xa26, and Xa11 shares with X21 most of the genome features. Other than presenting more CRISPR clusters and prophages than the previously sequenced strains, the integration of aggressiveness and cell biology points out that disease establishment is not a result of a single determinant factor as in other xanthomonads.

Ribosome footprint profiling enables elucidating the systemic regulation of fatty acid accumulation in Acer truncatum

BACKGROUND

The accumulation of fatty acids in plants covers a wide range of functions in plant physiology and thereby affects adaptations and characteristics of species. As the famous woody oilseed crop, Acer truncatum accumulates unsaturated fatty acids and could serve as the model to understand the regulation and trait formation in oil-accumulation crops. Here, we performed Ribosome footprint profiling combing with a multi-omics strategy towards vital time points during seed development, and finally constructed systematic profiling from transcription to proteomes. Additionally, we characterized the small open reading frames (ORFs) and revealed that the translational efficiencies of focused genes were highly influenced by their sequence features.

RESULTS

The comprehensive multi-omics analysis of lipid metabolism was conducted in A. truncatum. We applied the Ribo-seq and RNA-seq techniques, and the analyses of transcriptional and translational profiles of seeds collected at 85 and 115 DAF were compared. Key members of biosynthesis-related structural genes (LACS, FAD2, FAD3, and KCS) were characterized fully. More meaningfully, the regulators (MYB, ABI, bZIP, and Dof) were identified and revealed to affect lipid biosynthesis via post-translational regulations. The translational features results showed that translation efficiency tended to be lower for the genes with a translated uORF than for the genes with a non-translated uORF. They provide new insights into the global mechanisms underlying the developmental regulation of lipid metabolism.

CONCLUSIONS

We performed Ribosome footprint profiling combing with a multi-omics strategy in A. truncatum seed development, which provides an example of the use of Ribosome footprint profiling in deciphering the complex regulation network and will be useful for elucidating the metabolism of A. truncatum seed oil and the regulatory mechanisms.

The m6A reader PRRC2A is essential for meiosis I completion during spermatogenesis

N6-methyladenosine (m6A) and its reader proteins YTHDC1, YTHDC2, and YTHDF2 have been shown to exert essential functions during spermatogenesis. However, much remains unknown about m6A regulation mechanisms and the functions of specific readers during the meiotic cell cycle. Here, we show that the m6A reader Proline rich coiled-coil 2A (PRRC2A) is essential for male fertility. Germ cell-specific knockout of Prrc2a causes XY asynapsis and impaired meiotic sex chromosome inactivation in late-prophase spermatocytes. Moreover, PRRC2A-null spermatocytes exhibit delayed metaphase entry, chromosome misalignment, and spindle disorganization at metaphase I and are finally arrested at this stage. Sequencing data reveal that PRRC2A decreases the RNA abundance or improves the translation efficiency of targeting transcripts. Specifically, PRRC2A recognizes spermatogonia-specific transcripts and downregulates their RNA abundance to maintain the spermatocyte expression pattern during the meiosis prophase. For genes involved in meiotic cell division, PRRC2A improves the translation efficiency of their transcripts. Further, co-immunoprecipitation data show that PRRC2A interacts with several proteins regulating mRNA metabolism or translation (YBX1, YBX2, PABPC1, FXR1, and EIF4G3). Our study reveals post-transcriptional functions of PRRC2A and demonstrates its critical role in the completion of meiosis I in spermatogenesis.

Blocking Abundant RNA Transcripts by High-Affinity Oligonucleotides during Transcriptome Library Preparation

BACKGROUND

RNA sequencing has become the gold standard for transcriptome analysis but has an inherent limitation of challenging quantification of low-abundant transcripts. In contrast to microarray technology, RNA sequencing reads are proportionally divided in function of transcript abundance. Therefore, low-abundant RNAs compete against highly abundant - and sometimes non-informative - RNA species.

RESULTS

We developed an easy-to-use strategy based on high-affinity RNA-binding oligonucleotides to block reverse transcription and PCR amplification of specific RNA transcripts, thereby substantially reducing their abundance in the final sequencing library. To demonstrate the broad application potential of our method, we applied it to different transcripts and library preparation strategies, including YRNAs in small RNA sequencing of human blood plasma, mitochondrial rRNAs in both 3' end sequencing and long-read sequencing, and MALAT1 in single-cell 3' end sequencing. We demonstrate that the blocking strategy is highly efficient, reproducible, specific, and generally results in better transcriptome coverage and complexity.

CONCLUSION

Our method does not require modifications of the library preparation procedure apart from simply adding blocking oligonucleotides to the RT reaction and can thus be easily integrated into virtually any RNA sequencing library preparation protocol.

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.

Tumor-intrinsic YTHDF1 drives immune evasion and resistance to immune checkpoint inhibitors via promoting MHC-I degradation

The recently described role of RNA methylation in regulating immune cell infiltration into tumors has attracted interest, given its potential impact on immunotherapy response. YTHDF1 is a versatile and powerful m6A reader, but the understanding of its impact on immune evasion is limited. Here, we reveal that tumor-intrinsic YTHDF1 drives immune evasion and immune checkpoint inhibitor (ICI) resistance. Additionally, YTHDF1 deficiency converts cold tumors into responsive hot tumors, which improves ICI efficacy. Mechanistically, YTHDF1 deficiency inhibits the translation of lysosomal genes and limits lysosomal proteolysis of the major histocompatibility complex class I (MHC-I) and antigens, ultimately restoring tumor immune surveillance. In addition, we design a system for exosome-mediated CRISPR/Cas9 delivery to target YTHDF1 in vivo, resulting in YTHDF1 depletion and antitumor activity. Our findings elucidate the role of tumor-intrinsic YTHDF1 in driving immune evasion and its underlying mechanism.

A novel protein RASON encoded by a lncRNA controls oncogenic RAS signaling in KRAS mutant cancers

Mutations of the RAS oncogene are found in around 30% of all human cancers yet direct targeting of RAS is still considered clinically impractical except for the KRAS^G12C mutant. Here we report that RAS-ON (RASON), a novel protein encoded by the long intergenic non-protein coding RNA 00673 (LINC00673), is a positive regulator of oncogenic RAS signaling. RASON is aberrantly overexpressed in pancreatic ductal adenocarcinoma (PDAC) patients, and it promotes proliferation of human PDAC cell lines in vitro and tumor growth in vivo. CRISPR/Cas9-mediated knockout of Rason in mouse embryonic fibroblasts inhibits KRAS-mediated tumor transformation. Genetic deletion of Rason abolishes oncogenic KRAS-driven pancreatic and lung cancer tumorigenesis in LSL-Kras^G12D; Trp53^R172H/+ mice. Mechanistically, RASON directly binds to KRAS^G12D/V and inhibits both intrinsic and GTPase activating protein (GAP)-mediated GTP hydrolysis, thus sustaining KRAS^G12D/V in the GTP-bound hyperactive state. Therapeutically, deprivation of RASON sensitizes KRAS mutant pancreatic cancer cells and patient-derived organoids to EGFR inhibitors. Our findings identify RASON as a critical regulator of oncogenic KRAS signaling and a promising therapeutic target for KRAS mutant cancers.

A-MYB/TCFL5 regulatory architecture ensures the production of pachytene piRNAs in placental mammals

In male mice, the transcription factor A MYB initiates the transcription of pachytene piRNA genes during meiosis. Here, we report that A MYB activates the transcription factor Tcfl5 produced in pachytene spermatocytes. Subsequently, A MYB and TCFL5 reciprocally reinforce their own transcription to establish a positive feedback circuit that triggers pachytene piRNA production. TCFL5 regulates the expression of genes required for piRNA maturation and promotes transcription of evolutionarily young pachytene piRNA genes, whereas A-MYB activates the transcription of older pachytene piRNA genes. Intriguingly, pachytene piRNAs from TCFL5-dependent young loci initiates the production of piRNAs from A-MYB-dependent older loci ensuring the self-propagation of pachytene piRNAs. A MYB and TCFL5 act via a set of incoherent feedforward loops that drive regulation of gene expression by pachytene piRNAs during spermatogenesis. This regulatory architecture is conserved in rhesus macaque, suggesting that it was present in the last common ancestor of placental mammals.

Impact of RNA degradation on next-generation sequencing transcriptome data

RNA sequencing is an innovative technology to study transcriptomes in both biological and clinical research. However, clinical specimens from patients undergoing surgical operations have a major challenge due to sample degradation. This study replicated the process of RNA degradation by maintaining cells at room temperature to achieve none, slight, middle, and high levels of RNA degradation with decreasing RNA integrity numbers (RIN) of approximately 9.8, 6.7, 4.4, and 2.5, respectively. Next, the differential expression of mRNA and long non-coding RNA (lncRNA) was analyzed in the four degradation groups along with pathway enrichment analysis. The results showed that the similarity of lncRNAs exhibited significant differences even for a slight level of RNA degradation compared with the non-degraded RNA sample. Also, the RNA degradation process was found to be universal, global, and random; the differentially expressed genes increased with an increase in degradation but the pathway enrichment phenomenon was not significantly observed.

Nsun4 and Mettl3 mediated translational reprogramming of Sox9 promotes BMSC chondrogenic differentiation

The chondrogenic differentiation of bone marrow-derived mesenchymal stem cells (BMSCs) has been used in the treatment and repair of cartilage defects; however, the in-depth regulatory mechanisms by which RNA modifications are involved in this process are still poorly understood. Here, we found that Sox9, a critical transcription factor that mediates chondrogenic differentiation, exhibited enhanced translation by ribosome sequencing in chondrogenic pellets, which was accompanied by increased 5-methylcytosine (m⁵C) and N6-methyladenosine (m⁶A) levels. Nsun4-mediated m⁵C and Mettl3-mediated m⁶A modifications were required for Sox9-regulated chondrogenic differentiation. Interestingly, we showed that in the 3’UTR of Sox9 mRNA, Nsun4 catalyzed the m⁵C modification and Mettl3 catalyzed the m⁶A modification. Furthermore, we found that Nsun4 and Mettl3 co-regulated the translational reprogramming of Sox9 via the formation of a complex. Surface plasmon resonance (SPR) assays showed that this complex was assembled along with the recruitment of Ythdf2 and eEF1α-1. Moreover, BMSCs overexpressing Mettl3 and Nsun4 can promote the repair of cartilage defects in vivo. Taken together, our study demonstrates that m⁵C and m⁶A co-regulate the translation of Sox9 during the chondrogenic differentiation of BMSCs, which provides a therapeutic target for clinical implications.

Nsun4-mediated m5C and Mettl3-mediated m6A are found to be required for Sox9-regulated chondrogenic differentiation, whereby Nsun4 and Mettl3 interact with each other and recruit Ythdf2 and eEF1a-1 to form a complex at the 3’UTR of Sox9.

circFL-seq, A Full-length circRNA Sequencing Method

Due to overlapping sequences with linear cognates, identifying internal sequences of circular RNA (circRNA) remains a challenge. Recently, we have developed a full-length circRNA sequencing method (circFL-seq) and computational pipeline, to profile ordinary and fusion circRNA at the isoform level. Compared to short-read RNA-seq, rolling circular reverse transcription and nanopore long-read sequencing of circFL-seq make circRNA reads more than tenfold enriched, and show advantages for identification of both short (<100 nt) and long (>2,000 nt) circRNA transcripts. circFL-seq allows identification of differential alternative splicing suggested wide application prospects for functional studies of internal sequences in circRNAs. In addition, the experimental protocol and computational pipeline of circFL-seq shows better sensitivity and precision for identification of back-splicing junctions than current long-read sequencing methods. Together, the accurate identification and quantification of full-length circRNAs makes circFL-seq a potential tool for large-scale screening of functional circRNAs.

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.

TOX4 facilitates promoter-proximal pausing and C-terminal domain dephosphorylation of RNA polymerase II in human cells

TOX4 is one of the regulatory factors of PP1 phosphatases with poorly understood functions. Here we show that chromatin occupancy pattern of TOX4 resembles that of RNA polymerase II (Pol II), and its loss increases cellular level of C-terminal domain (CTD) phosphorylated Pol II but mainly decreases Pol II occupancy on promoters. In addition, elongation rate analyses by 4sUDRB-seq suggest that TOX4 restricts pause release and early elongation but promotes late elongation. Moreover, TT-seq analyses indicate that TOX4 loss mainly decreases transcriptional output. Mechanistically, TOX4 may restrict pause release through facilitating CTD serine 2 and DSIF dephosphorylation, and promote Pol II recycling and reinitiation through facilitating CTD serines 2 and 5 dephosphorylation. Furthermore, among the PP1 phosphatases, TOX4 preferentially binds PP1α and is capable of facilitating Pol II CTD dephosphorylation in vitro. These results lay the foundation for a better understanding of the role of TOX4 in transcriptional regulation.

As a role of TOX4, one of the regulatory proteins of PP1 phosphatases, in transcriptional regulation, authors here show that TOX4 restricts pause release and early productive elongation, and promotes Pol II recycling and transcriptional reinitiation.

Recent Major Transcriptomics and Epitranscriptomics Contributions toward Personalized and Precision Medicine

With the advent of genome-wide screening methods-beginning with microarray technologies and moving onto next generation sequencing methods-the era of precision and personalized medicine was born. Genomics led the way, and its contributions are well recognized. However, "other-omics" fields have rapidly emerged and are becoming as important toward defining disease causes and exploring therapeutic benefits. In this review, we focus on the impacts of transcriptomics, and its extension-epitranscriptomics-on personalized and precision medicine efforts. There has been an explosion of transcriptomic studies particularly in the last decade, along with a growing number of recent epitranscriptomic studies in several disease areas. Here, we summarize and overview major efforts for cancer, cardiovascular disease, and neurodevelopmental disorders (including autism spectrum disorder and intellectual disability) for transcriptomics/epitranscriptomics in precision and personalized medicine. We show that leading advances are being made in both diagnostics, and in investigative and landscaping disease pathophysiological studies. As transcriptomics/epitranscriptomics screens become more widespread, it is certain that they will yield vital and transformative precision and personalized medicine contributions in ways that will significantly further genomics gains.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.

Rubredoxin 1 Is Required for Formation of the Functional Photosystem II Core Complex in Arabidopsis thaliana

Chloroplast thylakoid protein rubredoxin 1 (RBD1) in Chlamydomonas and its cyanobacterial homolog RubA contain a rubredoxin domain. These proteins have been proposed to participate in the assembly of photosystem II (PSII) at early stages. However, the effects of inactivation of RBD1 on PSII assembly in higher plants are largely unclear. Here, we characterized an Arabidopsis rbd1 mutant in detail. A drastic reduction of intact PSII complex but relatively higher levels of assembly intermediates including PSII RC, pre-CP47, and pre-CP43 were found in rbd1. Polysome association and ribosome profiling revealed that ribosome recruitment of psbA mRNA is specifically reduced. Consistently, in vivo protein pulse-chase labeling showed that the rate of D1/pD1 synthesis is significantly reduced in rbd1 compared with WT. Moreover, newly synthesized mature D1 and pD1/D2 can assemble into the PSII reaction center (RC) complex but further formation of larger PSII complexes is nearly totally blocked in rbd1. Our data imply that RBD1 is not only required for the formation of a functional PSII core complex during the early stages of PSII assembly but may also be involved in the translation of D1 in higher plants.