scholarly journals Ribocutter: Cas9-mediated rRNA depletion from multiplexed riboseq libraries

2021 ◽  
Author(s):  
Oscar G Wilkins ◽  
Jernej Ule

RNA sequencing libraries produced from fragmented RNA, especially Ribo-seq libraries, contain high proportions of reads from abundant non-coding RNAs. Here, we describe a streamlined Cas9-based protocol for removing abundant rRNA/ncRNA contaminants from Ribo-seq (or other small RNA-seq) libraries and an easy-to-use software tool, ribocutter, for designing sgRNA templates. Following sgRNA template design, the pool of templates is in vitro transcribed using a 1-step commercial kit, which produces enough sgRNAs for multiple treatments. A single multiplexed sequencing library is then treated with Cas9/sgRNAs, followed by a short PCR program, which can increase the fraction of useful reads by more than 3-fold. Comparison of samples before and after depletion demonstrates that Cas9 produces minimal off-target effects and preserves key features (eg. footprint length, periodicity) of Ribo-seq libraries. The method is thus highly effective, costs <€1.00 per sample, and minimises non-specific depletion and technical variation between samples.

2019 ◽  
Author(s):  
Amelie J. Kraus ◽  
Benedikt G. Brink ◽  
T. Nicolai Siegel

SummaryIn most organisms, ribosomal RNA (rRNA) contributes to >85% of total RNA. Thus, to obtain useful information from RNA-sequencing (RNA-seq) analyses at reasonable sequencing depth, typically, mature polyadenylated transcripts are enriched or rRNA molecules are depleted. Targeted depletion of rRNA or other highly abundant transcripts is particularly useful when studying transcripts lacking a poly(A) tail, such as some non-coding RNAs (ncRNAs), most bacterial RNAs and partially degraded or immature transcripts. While several commercially available kits allow effective rRNA depletion, their efficiency relies on a high degree of sequence homology between oligonucleotide probes and the target RNA. This restricts the use of such kits to a limited number of organisms with conserved rRNA sequences.In this study we describe the use of biotinylated oligos and streptavidin-coated paramagnetic beads for the efficient and specific depletion of trypanosomal rRNA. Our approach reduces the levels of the most abundant rRNA transcripts to less than 5% with minimal off-target effects.By adjusting the sequence of the oligonucleotide probes, our approach can be used to deplete rRNAs or other abundant transcripts independent of species. Thus, our protocol provides a useful alternative for rRNA removal where enrichment of polyadenylated transcripts is not an option and commercial kits for rRNA are not available.


BMC Biology ◽  
2021 ◽  
Vol 19 (1) ◽  
Author(s):  
Milda Mickutė ◽  
Kotryna Kvederavičiūtė ◽  
Aleksandr Osipenko ◽  
Raminta Mineikaitė ◽  
Saulius Klimašauskas ◽  
...  

Abstract Background Targeted installation of designer chemical moieties on biopolymers provides an orthogonal means for their visualisation, manipulation and sequence analysis. Although high-throughput RNA sequencing is a widely used method for transcriptome analysis, certain steps, such as 3′ adapter ligation in strand-specific RNA sequencing, remain challenging due to structure- and sequence-related biases introduced by RNA ligases, leading to misrepresentation of particular RNA species. Here, we remedy this limitation by adapting two RNA 2′-O-methyltransferases from the Hen1 family for orthogonal chemo-enzymatic click tethering of a 3′ sequencing adapter that supports cDNA production by reverse transcription of the tagged RNA. Results We showed that the ssRNA-specific DmHen1 and dsRNA-specific AtHEN1 can be used to efficiently append an oligonucleotide adapter to the 3′ end of target RNA for sequencing library preparation. Using this new chemo-enzymatic approach, we identified miRNAs and prokaryotic small non-coding sRNAs in probiotic Lactobacillus casei BL23. We found that compared to a reference conventional RNA library preparation, methyltransferase-Directed Orthogonal Tagging and RNA sequencing, mDOT-seq, avoids misdetection of unspecific highly-structured RNA species, thus providing better accuracy in identifying the groups of transcripts analysed. Our results suggest that mDOT-seq has the potential to advance analysis of eukaryotic and prokaryotic ssRNAs. Conclusions Our findings provide a valuable resource for studies of the RNA-centred regulatory networks in Lactobacilli and pave the way to developing novel transcriptome and epitranscriptome profiling approaches in vitro and inside living cells. As RNA methyltransferases share the structure of the AdoMet-binding domain and several specific cofactor binding features, the basic principles of our approach could be easily translated to other AdoMet-dependent enzymes for the development of modification-specific RNA-seq techniques.


Genes ◽  
2021 ◽  
Vol 12 (5) ◽  
pp. 630
Author(s):  
Yongqing Lan ◽  
Meng Li ◽  
Shuangli Mi

Hematopoietic differentiation is a well-orchestrated process by many regulators such as transcription factor and long non-coding RNAs (lncRNAs). However, due to the large number of lncRNAs and the difficulty in determining their roles, the study of lncRNAs is a considerable challenge in hematopoietic differentiation. Here, through gene co-expression network analysis over RNA-seq data generated from representative types of mouse myeloid cells, we obtained a catalog of potential key lncRNAs in the context of mouse myeloid differentiation. Then, employing a widely used in vitro cell model, we screened a novel lncRNA, named Gdal1 (Granulocytic differentiation associated lncRNA 1), from this list and demonstrated that Gdal1 was required for granulocytic differentiation. Furthermore, knockdown of Cebpe, a principal transcription factor of granulocytic differentiation regulation, led to down-regulation of Gdal1, but not vice versa. In addition, expression of genes involved in myeloid differentiation and its regulation, such as Cebpa, were influenced in Gdal1 knockdown cells with differentiation blockage. We thus systematically identified myeloid differentiation associated lncRNAs and substantiated the identification by investigation of one of these lncRNAs on cellular phenotype and gene regulation levels. This study promotes our understanding of the regulation of myeloid differentiation and the characterization of roles of lncRNAs in hematopoietic system.


2019 ◽  
Vol 9 (1) ◽  
Author(s):  
Xia Tang ◽  
Delong Feng ◽  
Min Li ◽  
Jinxue Zhou ◽  
Xiaoyuan Li ◽  
...  

Abstract Fully elucidating the molecular mechanisms of non-coding RNAs (ncRNAs), including micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs), underlying hepatocarcinogenesis is challenging. We characterized the expression profiles of ncRNAs and constructed a regulatory mRNA-lncRNA-miRNA (MLMI) network based on transcriptome sequencing (RNA-seq) of hepatocellular carcinoma (HCC, n = 9) patients. Of the identified miRNAs (n = 203) and lncRNAs (n = 1,090), we found 16 significantly differentially expressed (DE) miRNAs and three DE lncRNAs. The DE RNAs were highly enriched in 21 functional pathways implicated in HCC (p < 0.05), including p53, MAPK, and NAFLD signaling. Potential pairwise interactions between DE ncRNAs and mRNAs were fully characterized using in silico prediction and experimentally-validated evidence. We for the first time constructed a MLMI network of reciprocal interactions for 16 miRNAs, three lncRNAs, and 253 mRNAs in HCC. The predominant role of MEG3 in the MLMI network was validated by its overexpression in vitro that the expression levels of a proportion of MEG3-targeted miRNAs and mRNAs was changed significantly. Our results suggested that the comprehensive MLMI network synergistically modulated carcinogenesis, and the crosstalk of the network provides a new avenue to accurately describe the molecular mechanisms of hepatocarcinogenesis.


2020 ◽  
Vol 6 (3) ◽  
pp. 32 ◽  
Author(s):  
Anna R. Dahlgren ◽  
Erica Y. Scott ◽  
Tamer Mansour ◽  
Erin N. Hales ◽  
Pablo J. Ross ◽  
...  

Long non-coding RNAs (lncRNAs) are untranslated regulatory transcripts longer than 200 nucleotides that can play a role in transcriptional, post-translational, and epigenetic regulation. Traditionally, RNA-sequencing (RNA-seq) libraries have been created by isolating transcriptomic RNA via poly-A+ selection. In the past 10 years, methods to perform ribosomal RNA (rRNA) depletion of total RNA have been developed as an alternative, aiming for better coverage of whole transcriptomic RNA, both polyadenylated and non-polyadenylated transcripts. The purpose of this study was to determine which library preparation method is optimal for lncRNA investigations in the horse. Using liver and cerebral parietal lobe tissues from two healthy Thoroughbred mares, RNA-seq libraries were prepared using standard poly-A+ selection and rRNA-depletion methods. Averaging the two biologic replicates, poly-A+ selection yielded 327 and 773 more unique lncRNA transcripts for liver and parietal lobe, respectively. More lncRNA were found to be unique to poly-A+ selected libraries, and rRNA-depletion identified small nucleolar RNA (snoRNA) to have a higher relative expression than in the poly-A+ selected libraries. Overall, poly-A+ selection provides a more thorough identification of total lncRNA in equine tissues while rRNA-depletion may allow for easier detection of snoRNAs.


2019 ◽  
Author(s):  
Hua Su ◽  
Lin Zhou ◽  
Na Li ◽  
Guowen Wang ◽  
Lingfang Wu ◽  
...  

AbstractHypoxic pulmonary hypertension (HPH) is a lethal disease. CircRNAs and m6A circRNAs have been reported to be associated with cancer progression, but the expression profiling of m6A circRNAs has not been identified in HPH. This study was to investigate the transcriptome-wide map of m6A circRNAs in HPH. In this study, hypoxia-induced PH rat model was established. Total RNA was extracted and purified from lungs of rats, then circRNAs were detected and annotated by RNA-seq analysis. m6A RNA Immunoprecipitation (MeRIP) was performed following rRNA depletion, and RNA-seq library was constructed. CircRNA–miRNA–mRNA co-expression network was also constructed. In vitro, total m6A was measured. m6A circXpo6 and m6A circTmtc3 were detected in pulmonary artery smooth muscle cells (PASMCs) and pulmonary artery endothelial cells (PAECs) exposed to 21% O2 and 1% O2 for 48 h, respectively. m6A abundance in 166 circRNAs was significantly upregulated and m6A abundance in 191 circRNAs was significantly downregulated in lungs of HPH rats. m6A abundance in circRNAs was significantly reduced in hypoxia in vitro. m6A circRNAs were mainly derived from single exons of protein-coding genes. m6A influenced the circRNA–miRNA–mRNA co-expression network in hypoxia. m6A circXpo6 and m6A circTmtc3 were downregulated in hypoxia. In general, our study firstly identified the transcriptome-wide map of m6A circRNAs in HPH. m6A level in circRNAs was decreased in lungs of HPH rats and in PASMCs and PAECs exposed to hypoxia. Downregulated or upregulated m6A level influenced circRNA–miRNA–mRNA co-expression network in HPH. Moreover, we firstly identified two downregulated m6A circRNAs in HPH: circXpo6 and circTmtc3. We suggested that m6A circRNAs may be used as a potential diagnostic marker or therapy target in the future.Author summaryHPH is a disease with great morbidity and mortality. It is often caused by chronic hypoxic lung diseases, such as chronic obstructive pulmonary disease and interstitial lung diseases. It lacks effective therapy methods so far. CircRNAs are a type of non-coding RNAs and can be used as biomarkers because they are differentially enriched in specific cell types or tissues and not easily degraded. m6A is identified as the most universal modification on non-coding RNAs in eukaryotes. CircRNAs can be modified by m6A. m6A circRNAs in HPH is not well understood yet. Here we identify the transcriptome-wide map of m6A circRNAs in HPH. We elucidate that m6A level in circRNAs is decreased in lungs of HPH rats and in PASMCs and PAECs exposed to hypoxia. We find that downregulated or upregulated m6A level influences circRNA– miRNA–mRNA co-expression network in HPH. Moreover, we are the first to identify two downregulated m6A circRNAs in HPH: circXpo6 and circTmtc3. We suggest that m6A circRNAs may be used as a potential diagnostic marker or therapy target in the future.


2021 ◽  
Author(s):  
Yichi Chen ◽  
Haitao Shang ◽  
Chunyue Wang ◽  
Jiaqi Zeng ◽  
Shentao Zhang ◽  
...  

Abstract Background The combination of sonodynamic therapy and oxygenation strategy is widely used in cancer treatment. However, due to the complexity, heterogeneity and irreversible hypoxic environment produced by Hepatocellular carcinoma (HCC) tissues, oxygen-enhancing sonodynamic therapy (SDT) have failed to achieve the desired results. With the emergence of ferroptosis with reactive oxygen species (ROS) cytotoxicity, this novel cell death method has attracted widespread attention. Methods In this study, nanobubbles (NBs) were connected with the sonosensitizer Indocyanine green (ICG) to construct a 2-in-1 nanoplatform loaded with RAS-selective lethal (RSL3, ferroptosis promoter) (RSL3@O2-ICG NBs), combined with oxygen-enhanced SDT and potent ferroptosis. In addition, nanobubbles (NBs) combined with low-frequency ultrasound (LFUS) are called ultrasound-targeted nanobubble destruction (UTND) to ensure specific drug release and improve safety. Results MDA/GSH and other related experimental results show that RSL3@O2-ICG NBs can enhance SDT and ferroptosis. Through RNA sequencing (RNA-seq), the differential expression of LncRNA and mRNA before and after synergistic treatment was identified, and then GO and KEGG pathways were used to enrich and analyze target genes and pathways related ferroptosis sensitivity. Then we searched for the expression of differentially expressed genes in the TCGA Hepatocellular carcinoma cohort. At the same time, we evaluated the proportion of immune cell infiltration and the identification of co-expression network modules and related prognostic analysis. We found that it was significantly related to the tumor microenvironment of hepatocellular carcinoma. The prognostic risk genes "SLC37A2" and "ITGB7" may represent new hepatocellular carcinoma ferroptosis-inducing markers and have guiding significance for treating hepatocellular carcinoma. Conclusions The therapeutic effect of the in vitro synergistic treatment has been proven to be significant, revealing the prospect of 2-in-1 nanobubbles combined with SDT and ferroptosis in treating HCC.


2014 ◽  
Author(s):  
Nicholas F Lahens ◽  
Ibrahim Halil Kavakli ◽  
Ray Zhang ◽  
Katharina Hayer ◽  
Michael B Black ◽  
...  

Background RNA sequencing (RNA-seq) is a powerful technique for identifying and quantifying transcription and splicing events, both known and novel. However, given its recent development and the proliferation of library construction methods, understanding the bias it introduces is incomplete but critical to realizing its value. Results Here we present a method, in vitro transcription sequencing (IVT-seq), for identifying and assessing the technical biases in RNA-seq library generation and sequencing at scale. We created a pool of > 1000 in vitro transcribed (IVT) RNAs from a full-length human cDNA library and sequenced them with poly-A and total RNA-seq, the most common protocols. Because each cDNA is full length and we show IVT is incredibly processive, each base in each transcript should be equivalently represented. However, with common RNA-seq applications and platforms, we find ~50% of transcripts have > 2-fold and ~10% have > 10-fold differences in within-transcript sequence coverage. Strikingly, we also find > 6% of transcripts have regions of high, unpredictable sequencing coverage, where the same transcript varies dramatically in coverage between samples, confounding accurate determination of their expression. To get at causal factors, we used a combination of experimental and computational approaches to show that rRNA depletion is responsible for the most significant variability in coverage and that several sequence determinants also strongly influence representation. Conclusions In sum, these results show the utility of IVT-seq in promoting better understanding of bias introduced by RNA-seq and suggest caution in its interpretation. Furthermore, we find that rRNA-depletion is responsible for substantial, unappreciated biases in coverage. Perhaps most importantly, these coverage biases introduced during library preparation suggest exon level expression analysis may be inadvisable.


2020 ◽  
Author(s):  
Edward Chen

The SARS-CoV-2 virus invades human cells via the angiotensin-converting enzyme 2 (ACE2)receptor. ACE2 expression in the lung cells creates a target for the SARS-CoV-2 virus, andconsequently infected individuals often see respiratory symptoms. Here, RNA sequencing(RNA-seq) datasets from in vitro human bronchial organoids (hBO) are analyzed to identifyunderlying gene expression changes as a result of an infection by SARS-CoV-2. Throughdifferential expression analysis, it was found that ACE2 expression levels did not changesignificantly after SARS-CoV-2 infection, but other genes that encode for major immunesignaling pathways were upregulated. This study can lend insight into potential targets for drugdevelopment and treatment for the COVID-19 disease.


Blood ◽  
2013 ◽  
Vol 122 (21) ◽  
pp. 48-48
Author(s):  
Min Luo ◽  
Mira Jeong ◽  
Deqiang Sun ◽  
Hyun-Jung Park ◽  
Liubin Yang ◽  
...  

The mammalian genome encodes a significant number of long non-coding RNAs (lncRNAs). The functions of some lncRNAs have been determined in biological processes, such as cancer progression, cell-cycle regulation and embryonic stem cell (ESC) pluripotency. However, our understanding of the basic function of lncRNAs in hematopoietic stem cell (HSC) is still limited. Here, we aim to identify the full complement of lncRNAs expressed in mouse HSCs and to determine whether they control HSC function. To uncover lncRNAs expressed in HSC across different ages, we performed RNA-seq on highly purified HSCs (SP-KSL-CD150+) from 4 month (m04), 12 month (m12) and 24 month (m24) old mice. With two biological replicates for each age, deep sequencing generated 368, 311 and 293 million mapped reads for m04, m12 and m24 HSC, respectively. After combining these datasets, assembly of over 1 billion mapped reads for the HSC transcriptome reconstructed 3,104 novel transcripts, which do not correspond to any UCSC, RefSeq or Ensemble known genes. Among them, 2,853 transcripts have multiple assembled exons and a total length >200 bp, representing potential lncRNAs. It has been shown that lncRNAs usually exhibit stage- or cell type-specific expression. To identify lncRNAs specifically expressed in HSC, we further performed RNA-seq on differentiated lymphoid lineage B cells (B220+) and myeloid lineage Granulocytes (Gr1+). Comparison of the expression of these 2,853 transcripts in the three cell types revealed that 173 transcripts are specifically expressed in HSC. As epigenetic mechanisms play critical roles to regulate gene transcription, we checked the chromatin map associated with those novel transcripts by ChIP-seq for H3K4me3, H3K27me3 and H3K36me3 in purified HSCs. Like protein-coding genes, these HSC specific novel transcripts typically contain the H3K4me3 mark at their transcriptional start site (TSS) and H3K36me3 along the gene body. Remarkably, one fifth of those 173 transcripts showed altered expression with HSC aging. Given that HSC function declines with aging, we hypothesize that those transcripts contribute to control HSC homeostasis. We selected three of these transcripts for further validation: LncHSC-1, LncHSC-2 and LncHSC-3. RT-PCR confirmed that they were highly expressed in stem and progenitor populations (KSL), but not in differentiated lineages (B220, CD4, CD8, Mac1, Gr1 and Ter119). Next, we generated retrovirally expressed-miRNA constructs to knockdown these transcripts. In vitro methocult colony forming assay showed that knockdown of LncHSC-1 in HSC significantly increased the colony number after second plating. Lineage analysis revealed that the majority of those cells are c-Kit+, and exhibit similar morphology, possibly representing expanded myeloid progenitors. To confirm our in vitro findings, we further examined their functions in vivo by HSC transplantation. Progenitors in which LncHSC-3 was knocked down failed to contribute to long-term hematopoietic reconstitution, as revealed by loss of retrovirally transduced population in the peripheral blood and bone marrow. In contrast, progenitors in which LncHSC-1 was knocked down resulted in augmented myeloid differentiation, consistent with in vitro CFU results that knockdown increased myeloid colony number. To understand the molecular mechanism through which lncRNAs influence hematopoiesis, we checked gene expression changes upon knockdown of specific transcripts in KSL cells. Overall, 80-100 genes were significantly changed after knockdown of specific transcripts, including cell cycle regulators and chromatin modification enzymes. For example, after LncHSC-3 knockdown, cell cycle regulator Cdkn1a (p21) expression increased, possibly contributing to the inhibition of hematopoietic reconstitution. In summary, here we carried out a comprehensive lncRNAs analysis in HSC and determined HSC specific novel transcripts. Loss-of-function experiments demonstrated that these transcripts may play important roles for HSC self-renewal and differentiation. These findings provide a useful resource to study lncRNA functions in normal hematopoiesis and disease progression. Disclosures: No relevant conflicts of interest to declare.


Sign in / Sign up

Export Citation Format

Share Document