RNA Mapping Research Articles

Absolute quantitative and base-resolution sequencing reveals comprehensive landscape of pseudouridine across the human transcriptome.

Pseudouridine (Ψ) is one of the most abundant modifications in cellular RNA. However, its function remains elusive, mainly due to the lack of highly sensitive and accurate detection methods. Here, we introduced 2-bromoacrylamide-assisted cyclization sequencing (BACS), which enables Ψ-to-C transitions, for quantitative profiling of Ψ at single-base resolution. BACS allowed the precise identification of Ψ positions, especially in densely modified Ψ regions and consecutive uridine sequences. BACS detected all known Ψ sites in human rRNA and spliceosomal small nuclear RNAs and generated the quantitative Ψ map of human small nucleolar RNA and tRNA. Furthermore, BACS simultaneously detected adenosine-to-inosine editing sites and N1-methyladenosine. Depletion of pseudouridine synthases TRUB1, PUS7 and PUS1 elucidated their targets and sequence motifs. We further identified a highly abundant Ψ114 site in Epstein-Barr virus-encoded small RNA EBER2. Surprisingly, applying BACS to a panel of RNA viruses demonstrated the absence of Ψ in their viral transcripts or genomes, shedding light on differences in pseudouridylation across virus families.

Analysis of miRNAs in milk of four livestock species

BackgroundMilk is essential for mammalian nutrition because it provides vital nutrients for growth and development. Milk composition, which is influenced by genetic and environmental factors, supports lactation, a complex process crucial for milk production and quality. Recent research has focused on noncoding RNAs, particularly microRNAs (miRNAs), which are present in body fluids and regulate gene expression post-transcriptionally. This study comprehensively characterizes miRNAs in milk of four livestock species, namely Bubalus bubalis, Capra hircus, Equus asinus, and Ovis aries and identifies potential target genes.ResultsHigh-throughput sequencing of milk RNA resulted in distinct read counts across species: B. bubalis (8,790,441 reads), C. hircus (12,976,275 reads), E. asinus (9,385,067 reads), and O. aries (7,295,297 reads). E. asinus had the highest RNA mapping rate (94.6%) and O. aries the lowest (84.8%). A substantially greater proportion of miRNAs over other small RNAs was observed for the donkey milk sample (7.74%) compared to buffalo (0.87%), goat (1.57%), and sheep (1.12%). Shared miRNAs, which included miR-200a, miR-200b, miR-200c, and miR-23a among others, showed varying expression levels across species, confirmed by qPCR analysis. Functional annotation of predicted miRNA target genes highlighted diverse roles, with an enrichment in functions linked to metabolism and immunity. Pathway analysis identified immune response pathways as significant, with several miRNAs targeting specific genes across species, suggesting their regulatory function in milk.ConclusionsBoth conserved and species-specific miRNAs were detected in milk of the investigated species. The identified target genes of these miRNAs have important roles in neonatal development, adaptation, growth, and immune response. Furthermore, they influence milk and meat production traits in livestock.

The RNA interference response to alphanodavirus replication in Phlebotomus papatasi sand fly cells.

In this study, we identified and assembled a strain of American nodavirus (ANV) in the Phlebotomus papatasi-derived PP9ad cell line. This strain most closely resembles Flock House virus and ANV identified in the Drosophila melanogaster S2/S2R cell line. Through small RNA sequencing and analysis, we demonstrate that ANV replication in PP9ad cells is primarily targeted by the exogenous small interfering RNA (exo-siRNA) pathway, with minimal engagement from the PIWI-interacting RNA (piRNA) pathway. In mosquitoes such as Aedes and Culex, the PIWI pathway is expanded and specialised, which actively limits virus replication. This is unlike in Drosophila spp., where the piRNA pathway does not restrict viral replication. In Lutzomyia sandflies (family Psychodidae), close relatives of Phlebotomus species and Drosophila, there appears to be an absence of virus-derived piRNAs. To investigate whether this absence is due to a lack of PIWI pathway proteins, we analysed the piRNA and siRNA diversity and repertoire in PP9ad cells. Previous assemblies of P. papatasi genome (Ppap_1.0) have revealed a patchy repertoire of the siRNA and piRNA pathways. Our analysis of the updated P. papatasi genome (Ppap_2.1) has shown no PIWI protein expansion in sandflies. We found that both siRNA and piRNA pathways are transcriptionally active in PP9ad cells, with genomic mapping of small RNAs generating typical piRNA signatures. Our results suggest that the piRNA pathway may not respond to virus replication in these cells, but an antiviral response is mounted via the exo-siRNA pathway.

RNA-sequencing of paired tape-strips and skin biopsies in atopic dermatitis reveals key differences.

Skin tape-strips and biopsies are widely used methods for investigating the skin in atopic dermatitis (AD). Biopsies are more commonly used but can cause scarring and pain, whereas tape-strips are noninvasive but sample less tissue. The study evaluated the performance of skin tape-strips and biopsies for studying AD. Whole-transcriptome RNA-sequencing was performed on paired tape-strips and biopsies collected from lesional and non-lesional skin from AD patients (n = 7) and non-AD controls (n = 5). RNA yield, mapping efficiency, and differentially expressed genes (DEGs) for the two methods (tape-strip/biopsy) and presence of AD (AD/non-AD) were compared. Tape-strips demonstrated a lower RNA yield (22 vs. 4596 ng) and mapping efficiency to known genes (28% vs. 93%) than biopsies. Gene-expression profiles of paired tape-strips and biopsies demonstrated a medium correlation (R2 = 0.431). Tape-strips and biopsies demonstrated systematic differences in measured expression levels of 6483 genes across both AD and non-AD samples. Tape-strips preferentially detected many itch (CCL3/CCL4/OSM) and immune-response (CXCL8/IL4/IL5/IL22) genes as well as markers of epidermal dendritic cells (CD1a/CD207), while certain cytokines (IL18/IL37), skin-barrier genes (KRT2/FLG2), and dermal fibroblasts markers (COL1A/COL3A) were preferentially detected by biopsies. Tape-strips identified more DEGs between AD and non-AD (3157 DEGs) then biopsies (44 DEGs). Tape-strips also detected higher levels of bacterial mRNA than biopsies. This study concludes that tape-strips and biopsies each demonstrate respective advantages for measuring gene-expression changes in AD. Thus, the specific skin layers and genes of interest should be considered before selecting either method.

Abstract PO4-23-11: Resistance Mechanisms to CDK4/6 Inhibitors and/or Tamoxifen Using Comprehensive Thermostable Group II Intron Reverse Transcriptase Sequencing

Abstract Background Understanding the mechanisms of resistance to CDK4/6 inhibitors (CDK4/6i) and endocrine therapy (ET) is pivotal in exploring new therapeutic strategies for hormone receptor-positive (HR+), HER2-negative metastatic breast cancer. To decipher these resistance mechanisms, we analyzed the alterations in comprehensive RNA-seq (coding and non-coding RNAs) of HR+ HER2-negative BC cell lines via thermostable group II intron reverse transcriptase sequencing (TGIRT-seq). Methods We established Tamoxifen-resistant (TMR), Abemaciclib-resistant (ACR), Palbociclib-resistant (PCR), Tamoxifen/Abemaciclib double-resistant (TMR-ACR), and Tamoxifen/Palbociclib double-resistant (TMR-PCR) BC cell lines from MCF7 and T47D HR+/HER2- BC cell lines through stepwise dose-escalation continuous drug exposure. We performed a TGIRT-seq transcriptomic analysis using a protocol that allows the sequencing of both long and short non-coding and protein-coding RNAs in a single library. All libraries were sequenced using paired-end 150 bp on the Novaseq platform, resulting in an average of 50 million reads per library. For analysis, raw reads underwent adapter trimming, small RNA mapping, whole genome mapping, and the generation of annotated genes' read count. We used raw counts to detect differentially expressed genes (DEGs) with DESeq2 in R using the cut-off (log2 [fold change] > 1, FDR < 1e-3). We then used the selected DEGs for Gene Set Enrichment Analysis (GSEA) and Kaplan-Meier survival analysis from the TCGA database. Results The TGIRT-seq analysis identified 1171 to 3472 DEGs in different drug resistance and cell line combinations. Most DEGs in resistant cells were consistently downregulated compared to parent cells. Seventy-five percent of DEGs were protein-coding genes, with the rest being non-coding RNAs, such as small nucleolar RNAs, microRNAs, long non-coding RNAs, and transposable element RNAs, which conventional RNA-seq poorly detects. In principal component analysis (PCA) plots, replicates per drug resistance clustered clearly in both cell line backgrounds. PCA also indicated distinctive pathways for acquiring drug resistance in MCF7 and T47D cells. The estrogen receptor (ER) expression decreased, while HER2 increased in resistant cell lines. DEGs identified were enriched in ESR1-related pathways in all MCF7 and T47D cell lines resistant models. Several ER-regulating genes like IL-1R1 and RET were upregulated, while ADCY1 was consistently downregulated across different resistance types. A significant overlap with single-resistance DEGs was observed among the DEGs in double-resistance cell lines, but 20-50% of up/down-regulated DEGs in double-resistance cell lines were uniquely altered. Several previously identified DEGs (e.g., SALL4, TOP2A) and 21 novel candidate genes correlated with poor survival outcomes. Conclusion The analysis identified unique DEGs in double resistance cell lines, suggesting that double resistance mechanisms may not merely be a cumulative effect of single resistance mechanisms, necessitating further investigation and validation. CDK4/6i and/or ET-resistant BC cell lines displayed significant transcriptome reprogramming during the development of drug resistance. In previously published research, ESR1-related pathway alterations were proposed in tamoxifen resistance cell lines. Here we find that CDK4/6i-resistant cells also modify these ESR1-related pathways. Additionally, we identified several targetable genes, such as IL-1R1 and RET, involved in ESR1-related pathways that could pave the way for developing new treatment strategies. Citation Format: Toshiaki Iwase, Hengyi Xu, Nakyung Oh, Jangsoon Lee, Alan M Lambowitz, Naoto Ueno. Resistance Mechanisms to CDK4/6 Inhibitors and/or Tamoxifen Using Comprehensive Thermostable Group II Intron Reverse Transcriptase Sequencing [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO4-23-11.

Co-localization and analysis of miR477b with fiber length quantitative trait loci in cotton.

Cotton is an important cash crop for the textile industry. However, the understanding of natural genetic variation of fiber elongation in relation to miRNA is lacking. A miRNA gene (miR477b) was found to co-localize with a previously mapped fiber length (FL) quantitative trait locus (QTL). The miR477b was differentially expressed during fiber elongation between two backcross inbred lines (BILs) differing in FL and its precursor sequences. Bioinformatics and qRT-PCR analysis were further used to analyse the miRNA genes, which could produce mature miR477b. Cotton plants with virus-induced gene silencing (VIGS) constructs to over-express the allele of miR477b from the BIL with longer fibers had significantly longer fibers as compared with negative control plants, while the VIGS plants with suppressed miRNA expression had significantly shorter fibers. The expression level of the target gene (DELLA) and related genes (RDL1 and EXPA1 for DELLA through HOX3 protein) in the two BILs and/or the VIGS plants were generally congruent, as expected. This report represents one of the first comprehensive studies to integrate QTL linkage mapping and physical mapping of small RNAs with both small and mRNA transcriptome analysis, followed by VIGS, to identify candidate small RNA genes affecting the natural variation of fiber elongation in cotton.

Bioorthogonal masked acylating agents for proximity-dependent RNA labelling.

RNA localization is highly regulated, with subcellular organization driving context-dependent cell physiology. Although proximity-based labelling technologies that use highly reactive radicals or carbenes provide a powerful method for unbiased mapping of protein organization within a cell, methods for unbiased RNA mapping are scarce and comparably less robust. Here we develop α-alkoxy thioenol and chloroenol esters that function as potent acylating agents upon controlled ester unmasking. We pair these probes with subcellular-localized expression of a bioorthogonal esterase to establish a platform for spatial analysis of RNA: bioorthogonal acylating agents for proximity labelling and sequencing (BAP-seq). We demonstrate that, by selectively unmasking the enol probe in a locale of interest, we can map RNA distribution in membrane-bound and membrane-less organelles. The controlled-release acylating agent chemistry and corresponding BAP-seq method expand the scope of proximity labelling technologies and provide a powerful approach to interrogate the cellular organization of RNAs.

Challenges to mapping and defining m6A function in viral RNA.

Viral RNA molecules contain multiple layers of regulatory information. This includes features beyond the primary sequence, such as RNA structures and RNA modifications, including N6-methyladenosine (m6A). Many recent studies have identified the presence and location of m6A in viral RNA and have found diverse regulatory roles for this modification during viral infection. However, to date, viral m6A mapping strategies have limitations that prevent a complete understanding of the function of m6A on individual viral RNA molecules. While m6A sites have been profiled on bulk RNA from many viruses, the resulting m6A maps of viral RNAs described to date present a composite picture of m6A across viral RNA molecules in the infected cell. Thus, for most viruses, it is unknown if unique viral m6A profiles exist throughout infection, nor if they regulate specific viral life cycle stages. Here, we describe several challenges to defining the function of m6A in viral RNA molecules and provide a framework for future studies to help in the understanding of how m6A regulates viral infection.

Abstract 4351: Differential dependency mapping of chimeric RNAs across cancer reveals a new landscape of functional fusion transcripts

Abstract Chimeric RNAs are RNA transcripts containing sequences originating from multiple distinct genetic loci and can form through a variety of mechanisms such as DNA rearrangement, cis-splicing of adjacent genes, or RNA trans-splicing. Chimeras have long been known to be a hallmark of cancer and due to their unique properties are promising targets for precision medicine. Appropriately, fusion proteins transcribed from chimeras such as BCR-ABL1 and TPM3-NTRK1 have been shown to be effective targets. Due to these past successes, there exists an opportunity to identify new chimeras as therapeutic targets. While modern chimera prediction software allows for the fast and accurate identification of chimeric RNAs from RNA sequencing data, investigations separating therapeutically relevant transcripts from “transcriptional noise” remain lacking. In this study we performed an in-silico functional screen of chimeras across cell lines representing a wide variety of cancers. We used state-of-the-art chimeric RNA prediction software to create a database of chimeric RNAs across 1017 cell lines from The Cancer Cell Line Encyclopedia. For each chimera we assessed factors such as frequency, recurrence, breakpoint coordinates, frame, and coding potential. To identify specific functional transcripts, we integrated publicly available shRNA knockdown data with our predictions and developed an in-silico functional analysis pipeline comparing differential knockdown effects between chimera and corresponding parental transcripts. For each transcript we assessed the average fold change of chimera-mapping probes, a function score defined as the difference between the fold change of chimera-mapping and parent-specific probes, and a p-value assessing the confidence of our functional score. From our initial screen of 127,819 transcripts, we identified 1088 high-confidence functional chimeras. We successfully identified nearly all known functional chimeras screened including PAX3-FOXO1, EWSR1-FLI1, and TCF3-PBX1. We also identified previously unknown chimeras that that we predict have a function in cell growth and proliferation. Overall, the results of our study reveal a new landscape of functionally relevant chimeric transcripts in cancer. Follow-up studies can further investigate these transcripts to determine their potential as therapeutic targets. Citation Format: Samir Lalani, Sehajroop Gadh, Justin Elfman, Sandeep Singh, Hui Li. Differential dependency mapping of chimeric RNAs across cancer reveals a new landscape of functional fusion transcripts [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4351.

Hydroxypropyl-Beta Cyclodextrin Barrier Prevents Respiratory Viral Infections: A Preclinical Study.

The emergence and mutation of pathogenic viruses have been occurring at an unprecedented rate in recent decades. The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has developed into a global public health crisis due to extensive viral transmission. In situ RNA mapping has revealed angiotensin-converting enzyme 2 (ACE2) expression to be highest in the nose and lower in the lung, pointing to nasal susceptibility as a predominant route for infection and the cause of subsequent pulmonary effects. By blocking viral attachment and entry at the nasal airway using a cyclodextrin-based formulation, a preventative therapy can be developed to reduce viral infection at the site of entry. Here, we assess the safety and antiviral efficacy of cyclodextrin-based formulations. From these studies, hydroxypropyl beta-cyclodextrin (HPBCD) and hydroxypropyl gamma-cyclodextrin (HPGCD) were then further evaluated for antiviral effects using SARS-CoV-2 pseudotypes. Efficacy findings were confirmed with SARS-CoV-2 Delta variant infection of Calu-3 cells and using a K18-hACE2 murine model. Intranasal pre-treatment with HPBCD-based formulations reduced viral load and inflammatory signaling in the lung. In vitro efficacy studies were further conducted using lentiviruses, murine hepatitis virus (MHV), and influenza A virus subtype H1N1. These findings suggest HPBCD may be used as an agnostic barrier against transmissible pathogens, including but not limited to SARS-CoV-2.

Chemical diversity of reagents that modify RNA 2'-OH in water: a review.

Electrophilic water-soluble compounds have proven versatile in reacting selectively with 2'-OH groups in RNA, enabling structure mapping, probing, caging, labeling, crosslinking, and conjugation of RNAs in vitro and in living cells. While early work focused on one or two types of reagents with limited properties, recent studies have greatly diversified the structure, properties, and applications of these reagents. Here we review the scope of documented RNA hydroxyl-reactive species reported to date, with an eye to the effects of chemical structure on reactivity with RNA and other useful properties. Multiple forms of carbonyl electrophiles are now known to react at the 2'-OH, and recently, sulfonyl and aryl electrophiles have also been documented to form bonds there in high yields as well. In addition to electrophilicity, data also point to significant effects of reagent stability, steric bulk, and chirality on reaction yields and selectivity. Finally, we outline reagent properties and principles that define utility in applications with RNA, with an eye to the design of future reagents.

Circular RNA mapping reveals CircCWC22 as a MiR-3059-x sponge in yak fat deposition by regulating HMGCL

The regulatory mechanisms and functions of circular RNAs (circRNAs) in yak intramuscular fat (IMF) deposition remain unclear. This study aimed to investigate yak circRNAs with high and low IMF content using high-throughput sequencing. A total of 270 differentially expressed circRNAs were identified, of which 129 were upregulated and 141 were downregulated. Among these circRNAs, circCWC22, derived from the yak CWC22 gene, was further studied to understand its functions and regulatory mechanisms. Sequencing and RNase R processing confirmed the circular nature of circCWC22. By constructing a circRNA-miRNA-mRNA co-expression network, the potential regulatory pathway of circCWC22/miR-3059-x/HMGCL was identified. To investigate the roles of circCWC22, miR-3059-x, and HMGCL in the deposition of yak intramuscular preadipocytes (YIMAs), CCK-8, EdU, BODIPY, triglyceride content, and qRT-PCR analyses were performed. The results demonstrated that circCWC22, miR-3059-x, and HMGCL promoted the differentiation and inhibited the proliferation of YIMAs. Using the dual-luciferase reporter system and qRT-PCR, we confirmed that circCWC22 adsorbed miR-3059-x, and HMGCL was identified as a target gene of miR-3059-x. In conclusion, this study uncovered a large number of potential circRNAs involved in IMF deposition and highlighted the significant role of circCWC22 in yak IMF deposition via the circCWC22/miR-3059-x/HMGCL axis.

Genetic basis of the oil biosynthesis in ultra-high-oil maize grains with an oil content exceeding 20.

Vegetable oil is an important part of the human diet and has multiple industrial uses. The rapid increase in vegetable oil consumption has necessitated the development of viable methods for optimizing the oil content of plants. The key genes regulating the biosynthesis of maize grain oil remain mostly uncharacterized. In this study, by analyzing oil contents and performing bulked segregant RNA sequencing and mapping analyses, we determined that su1 and sh2-R mediate the shrinkage of ultra-high-oil maize grains and contribute to the increase in the grain oil content. Functional kompetitive allele-specific PCR (KASP) markers developed for su1 and sh2-R detected su1su1Sh2Sh2, Su1Su1sh2sh2, and su1su1sh2sh2 mutants among 183 sweet maize inbred lines. An RNA sequencing (RNA-seq) analysis indicated that genes differentially expressed between two conventional sweet maize lines and two ultra-high-oil maize lines were significantly associated with linoleic acid metabolism, cyanoamino acid metabolism, glutathione metabolism, alanine, aspartate, and glutamate metabolism, and nitrogen metabolism. A bulk segregant analysis and sequencing (BSA-seq) analysis identified another 88 genomic intervals related to grain oil content, 16 of which overlapped previously reported maize grain oil-related QTLs. The combined analysis of BSA-seq and RNA-seq data enabled the identification of candidate genes. The KASP markers for GRMZM2G176998 (putative WD40-like beta propeller repeat family protein), GRMZM2G021339 (homeobox-transcription factor 115), and GRMZM2G167438 (3-ketoacyl-CoA synthase) were significantly related to maize grain oil content. Another candidate gene, GRMZM2G099802 (GDSL-like lipase/acylhydrolase), catalyzes the final step of the triacylglycerol synthesis pathway and was expressed at significantly higher levels in the two ultra-high-oil maize lines than in the two conventional sweet maize lines. These novel findings will help clarify the genetic basis of the increased oil production in ultra-high-oil maize lines with grain oil contents exceeding 20%. The KASP markers developed in this study may be useful for breeding new high-oil sweet maize varieties.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy.

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

In vivo secondary structural analysis of Influenza A virus genomic RNA

Influenza A virus (IAV) is a respiratory virus that causes epidemics and pandemics. Knowledge of IAV RNA secondary structure in vivo is crucial for a better understanding of virus biology. Moreover, it is a fundament for the development of new RNA-targeting antivirals. Chemical RNA mapping using selective 2’-hydroxyl acylation analyzed by primer extension (SHAPE) coupled with Mutational Profiling (MaP) allows for the thorough examination of secondary structures in low-abundance RNAs in their biological context. So far, the method has been used for analyzing the RNA secondary structures of several viruses including SARS-CoV-2 in virio and in cellulo. Here, we used SHAPE-MaP and dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) for genome-wide secondary structure analysis of viral RNA (vRNA) of the pandemic influenza A/California/04/2009 (H1N1) strain in both in virio and in cellulo environments. Experimental data allowed the prediction of the secondary structures of all eight vRNA segments in virio and, for the first time, the structures of vRNA5, 7, and 8 in cellulo. We conducted a comprehensive structural analysis of the proposed vRNA structures to reveal the motifs predicted with the highest accuracy. We also performed a base-pairs conservation analysis of the predicted vRNA structures and revealed many highly conserved vRNA motifs among the IAVs. The structural motifs presented herein are potential candidates for new IAV antiviral strategies.

Characterizing expression changes in noncoding RNAs during aging and heterochronic parabiosis across mouse tissues

Molecular mechanisms of organismal and cell aging remain incompletely understood. We, therefore, generated a body-wide map of noncoding RNA (ncRNA) expression in aging (16 organs at ten timepoints from 1 to 27 months) and rejuvenated mice. We found molecular aging trajectories are largely tissue-specific except for eight broadly deregulated microRNAs (miRNAs). Their individual abundance mirrors their presence in circulating plasma and extracellular vesicles (EVs) whereas tissue-specific ncRNAs were less present. For miR-29c-3p, we observe the largest correlation with aging in solid organs, plasma and EVs. In mice rejuvenated by heterochronic parabiosis, miR-29c-3p was the most prominent miRNA restored to similar levels found in young liver. miR-29c-3p targets the extracellular matrix and secretion pathways, known to be implicated in aging. We provide a map of organism-wide expression of ncRNAs with aging and rejuvenation and identify a set of broadly deregulated miRNAs, which may function as systemic regulators of aging via plasma and EVs.

Observation of structural switch in nascent SAM-VI riboswitch during transcription at single-nucleotide and single-molecule resolution

Growing RNAs fold differently as they are transcribed, which modulates their finally adopted structures. Riboswitches regulate gene expression by structural change, which are sensitive to co-transcriptionally structural biology. Here we develop a strategy to track the structural change of RNAs during transcription at single-nucleotide and single-molecule resolution and use it to monitor individual transcripts of the SAM-VI riboswitch (riboSAM) as transcription proceeds, observing co-existence of five states in riboSAM. We report a bifurcated helix in one newly identified state from NMR and single-molecule FRET (smFRET) results, and its presence directs the translation inhibition in our cellular translation experiments. A model is proposed to illustrate the distinct switch patterns and gene-regulatory outcome of riboSAM when SAM is present or absent. Our strategy enables the precise mapping of RNAs’ conformational landscape during transcription, and may combine with detection methods other than smFRET for structural studies of RNAs in general.

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

Regulation of RNase E during the UV stress response in the cyanobacterium Synechocystis sp. PCC 6803.

Endoribonucleases govern the maturation and degradation of RNA and are indispensable in the posttranscriptional regulation of gene expression. A key endoribonuclease in Gram-negative bacteria is RNase E. To ensure an appropriate supply of RNase E, some bacteria, such as Escherichia coli, feedback-regulate RNase E expression via the rne 5'-untranslated region (5' UTR) in cis. However, the mechanisms involved in the control of RNase E in other bacteria largely remain unknown. Cyanobacteria rely on solar light as an energy source for photosynthesis, despite the inherent ultraviolet (UV) irradiation. Inthisstudy, we first investigated globally the changes in gene expression in the cyanobacterium Synechocystis sp. PCC 6803 after a brief exposure to UV. Among the 407 responding genes 2 h after UV exposure was a prominent upregulation of rnemRNA level. Moreover, the enzymatic activity ofRNaseErapidly increased as well, although the protein stability decreased. This unique response was underpinned by the increased accumulation of full-length rne mRNA caused by the stabilization of its 5'UTR and suppression of premature transcriptional termination, but not by an increased transcription rate. Mapping of RNA 3' ends and in vitro cleavage assays revealed that RNase E cleaves within a stretch of six consecutive uridine residues within the rne 5'UTR, indicating autoregulation. These observations suggest that RNase E in cyanobacteria contributes to reshaping the transcriptome during the UV stress response and that its required activity level is secured atthe RNA level despite the enhanced turnover of the protein.

Selective Chemical Labeling Strategy for Oligonucleotides Determination: A First Application to Full-Range Profiling of Transfer RNA Modifications.

To date, the extremely high polarity and poor signal intensity of macromolecular nucleic acids are greatly impeding the progress of mass spectrometry technology in the quality control of nucleic acid drugs and the characterization of DNA oxidation and RNA modifications. We recently described a general N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBSTFA) labeling method for oligonucleotide determination and applied it to the full-range profiling of tRNA in vitro and in vivo studies for the first time. The primary advantages of this method include strong retention, no observable byproducts, predictable and easily interpreted MS2 data, and the circumvention of instrument harmful reagents that were necessary in previous methods. Selective labeling of N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide to the terminal phosphate groups of oligonucleotides endows it broadly applicable for DNA/RNA profiling. Moreover, the improvement of sequence coverage was achieved in yeast tRNAphe(GAA) analysis owing to this method's good detection capability of 1-12 nucleotides in length. We also extended this strategy to determine the abundance of modified bases and discover new modifications via digesting RNA into single-nucleotide products, promoting the comprehensive mapping of RNA. The easy availability of derivatization reagent and the simple, rapid one-step reaction render it easy to operate for researchers. When applied in characterizing tRNAs in HepG2 cells and rats with nonalcoholic fatty liver disease, a fragment of U[m1G][m2G], specific for tRNAAsn(QUU) in cells, was significantly upregulated, indicating a possible clue to nonalcoholic fatty liver disease pathogenesis.