1
|
miRNA-seq identification and clinical validation of CD138+ and circulating miR-25 in treatment response of multiple myeloma. J Transl Med 2023; 21:245. [PMID: 37024879 PMCID: PMC10080848 DOI: 10.1186/s12967-023-04034-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2023] [Accepted: 03/03/2023] [Indexed: 04/08/2023] Open
Abstract
BACKGROUND Despite significant advancements in multiple myeloma (MM) therapy, the highly heterogenous treatment response hinders reliable prognosis and tailored therapeutics. Herein, we have studied the clinical utility of miRNAs in ameliorating patients' management. METHODS miRNA-seq was performed in bone marrow CD138+ plasma cells (PCs) of 24 MM and smoldering MM (sMM) patients to analyze miRNAs profile. CD138+ and circulating miR-25 levels were quantified using in house RT-qPCR assays in our screening MM/sMM cohort (CD138+ plasma cells n = 167; subcohort of MM peripheral plasma samples n = 69). Two external datasets (Kryukov et al. cohort n = 149; MMRF CoMMpass study n = 760) served as institutional-independent validation cohorts. Patients' mortality and disease progression were assessed as clinical endpoints. Internal validation was performed by bootstrap analysis. Clinical benefit was estimated by decision curve analysis. RESULTS miRNA-seq highlighted miR-25 of CD138+ plasma cells to be upregulated in MM vs. sMM, R-ISS II/III vs. R-ISS I, and in progressed compared to progression-free patients. The analysis of our screening cohort highlighted that CD138+ miR-25 levels were correlated with short-term progression (HR = 2.729; p = 0.009) and poor survival (HR = 4.581; p = 0.004) of the patients; which was confirmed by Kryukov et al. cohort (HR = 1.878; p = 0.005) and MMRF CoMMpass study (HR = 1.414; p = 0.039) validation cohorts. Moreover, multivariate miR-25-fitted models contributed to superior risk-stratification and clinical benefit in MM prognostication. Finally, elevated miR-25 circulating levels were correlated with poor survival of MM patients (HR = 5.435; p = 0.021), serving as a potent non-invasive molecular prognostic tool. CONCLUSIONS Our study identified miR-25 overexpression as a powerful independent predictor of poor treatment outcome and post-treatment progression, aiding towards modern non-invasive disease prognosis and personalized treatment decisions.
Collapse
|
2
|
Characteristics of Human and Microbiome RNA Profiles in Saliva. RNA Biol 2023; 20:398-408. [PMID: 37401028 PMCID: PMC10321209 DOI: 10.1080/15476286.2023.2229596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Revised: 06/08/2023] [Accepted: 06/17/2023] [Indexed: 07/05/2023] Open
Abstract
Saliva is a convenient non-invasive source of liquid biopsy to monitor human health and diagnose diseases. In particular, extracellular vesicles (EVs) in saliva can potentially reveal clinically relevant information for systemic health. Recent studies have shown that RNA in saliva EVs could be exploited as biomarkers for disease diagnosis. However, there is no standardized protocol for profiling RNA in saliva EV nor clear guideline on selecting saliva fractions for biomarker analysis. To address these issues, we established a robust protocol for small RNA profiling from fractionated saliva. With this method, we performed comprehensive small RNA sequencing of four saliva fractions, including cell-free saliva (CFS), EV-depleted saliva (EV-D), exosome (EXO), and microvesicle (MV) from ten healthy volunteers. By comparing the expression profiles of total RNA from these fractions, we found that MV was most enriched in microbiome RNA (76.2% of total reads on average), whereas EV-D was notably enriched in human RNA (70.3% of total reads on average). As for human RNA composition, CFS and EV-D were both enriched in snoRNA and tRNA compared with the two EV fractions (EXO and MV, P < 0.05). Interestingly, EXO and MV had highly correlated expression profiles for various noncoding RNAs such as miRNA, tRNA, and yRNA. Our study revealed unique characteristics of circulating RNAs in various saliva fractions, which provides a guideline on preparing saliva samples to study specific RNA biomarkers of interest.
Collapse
|
3
|
Determining extracellular vesicles properties and miRNA cargo variability in bovine milk from healthy cows and cows undergoing subclinical mastitis. BMC Genomics 2022; 23:189. [PMID: 35255807 PMCID: PMC8903571 DOI: 10.1186/s12864-022-08377-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2021] [Accepted: 12/03/2021] [Indexed: 12/13/2022] Open
Abstract
Background Subclinical mastitis, the inflammation of the mammary gland lacking clinical symptoms, is one of the most prevalent and costly diseases in dairy farming worldwide. Milk microRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) have been proposed as potential biomarkers of different mammary gland conditions, including subclinical mastitis. However, little is known about the robustness of EVs analysis regarding sampling time-point and natural infections. To estimate the reliability of EVs measurements in raw bovine milk, we first evaluated changes in EVs size and concentration using Tunable Resistive Pulse Sensing (TRPS) during three consecutive days of sampling. Then, we analysed daily differences in miRNA cargo using small RNA-seq. Finally, we compared milk EVs differences from naturally infected udder quarters with their healthy adjacent quarters and quarters from uninfected udders, respectively. Results We found that the milk EV miRNA cargo was very stable over the course of three days regardless of the health status of the quarter, and that infected quarters did not induce relevant changes in milk EVs of adjacent healthy quarters. Chronic subclinical mastitis induced changes in milk EV miRNA cargo, but neither in EVs size nor concentration. We observed that the changes in immunoregulatory miRNAs in quarters with chronic subclinical mastitis were cow-individual, however, the most upregulated miRNA was bta-miR-223-3p across all individuals. Conclusions Our results showed that the miRNA profile and particle size characteristics remained constant throughout consecutive days, suggesting that miRNAs packed in EVs are physiological state-specific. In addition, infected quarters were solely affected while adjacent healthy quarters remained unaffected. Finally, the cow-individual miRNA changes pointed towards infection-specific alterations. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08377-z.
Collapse
|
4
|
Identification and expression pattern analysis of miRNAs in pectoral muscle during pigeon ( Columba livia) development. PeerJ 2021; 9:e11438. [PMID: 34221709 PMCID: PMC8234919 DOI: 10.7717/peerj.11438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Accepted: 04/21/2021] [Indexed: 11/20/2022] Open
Abstract
MicroRNAs (miRNAs) are a group of crucial regulators in the process of animal growth and development. However, little is known about the expression and function of miRNAs in pigeon muscles. To identify the miRNAs participating in the rapid development of pigeon pectoral muscles and quantitate their expression levels of pectoral muscles in different age stages, we performed miRNA transcriptome analysis in pigeon pectoral muscles by sequencing small RNAs over three different age stages (1-day old, 28 days old, and 2 years old). Dual-luciferase reporter assay was applied to validate the interaction between miRNA and its target gene. We identified 304 known miRNAs, 201 conserved miRNAs, and 86 novel miRNAs in pigeon pectoral muscles. 189 differentially expressed (DE) miRNAs were screened out during pigeon development. A short time-series expression miner (STEM) analysis indicated 89 DE miRNAs were significantly clustered in a progressively decreasing expression profile, and mainly enriched in biosynthesis-related GO categories and signaling pathways for MAPK and TGF-β. Dual-luciferase reporter assay indicated that a progressively down-regulated miRNA (miR-20b-5p) could directly target Krüppel-like factor 3 (KLF3) gene. To sum-up, our data expand the repertoire of pigeon miRNAs and enhance understanding of the mechanisms underlying rapid development in squabs.
Collapse
|
5
|
Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea. Funct Integr Genomics 2021; 21:283-298. [PMID: 33630193 DOI: 10.1007/s10142-021-00777-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Revised: 01/12/2021] [Accepted: 02/09/2021] [Indexed: 11/30/2022]
Abstract
Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.
Collapse
|
6
|
Computational Identification of ceRNA and Reconstruction of ceRNA Regulatory Network Based on RNA-seq and Small RNA-seq Data in Plants. Methods Mol Biol 2021; 2328:261-275. [PMID: 34251632 DOI: 10.1007/978-1-0716-1534-8_17] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Competing endogenous RNAs (ceRNAs) are transcripts with the ability to competitively titrate microRNAs (miRNAs) against miRNA repressing target genes to post-transcriptionally regulate the expression of corresponding miRNAs. It is a newly discovered gene regulation pattern between longer RNA and miRNA molecules. Recent research has gradually revealed the functional significance of ceRNAs in regulating normal development and stress response processes in plants and animals, as well as in cancer genesis and metastasis. Therefore, ceRNA identification is an important and necessary step to deepen our understanding of the regulation mechanisms of various biological processes. Here, we provide a pipeline used to computationally identify plant ceRNAs and reconstruct ceRNA regulatory networks based on RNA-seq and small RNA-seq data.
Collapse
|
7
|
Abstract
The study of small RNAs (sRNAs) by next-generation sequencing (NGS) is challenged by bias issues during library preparation. Several types of sRNAs such as plant microRNAs (miRNAs) carry a 2'-O-methyl (2'-OMe) modification at their 3' terminal nucleotide. This modification adds another level of difficulty as it inhibits 3' adapter ligation. We previously demonstrated that modified versions of the "TruSeq (TS)" protocol have less bias and an improved detection of 2'-OMe RNAs. Here we describe in detail protocol "TS5," which showed the best overall performance. We also provide guidelines for bioinformatics analysis of the sequencing data.
Collapse
|
8
|
Small double-stranded RNA with anti-HIV activity abundantly produced by Bacillus subtilis MTCC5480 isolated from fermented soybean. Int J Biol Macromol 2020; 161:828-835. [PMID: 32553954 DOI: 10.1016/j.ijbiomac.2020.06.112] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2020] [Revised: 06/10/2020] [Accepted: 06/11/2020] [Indexed: 12/27/2022]
Abstract
Anti-viral RNA therapy is on high demand nowadays due to the emergence of several new viral infections. The small non-coding regulatory RNAs (dsRNA) from the microbial sources are not yet explored for anti-viral activity. In this study, we assessed the anti-HIV activity of the small dsRNA produced by 12 different microbial species isolated from naturally fermented foods of North-East India. For this, we selectively extracted the dsRNA from the microbial culture, confirmed its double-stranded nature by immunoblotting, and deep sequenced the cDNA library using Illumina platform. Further, we used conventional algorithms to predict the potential targets of the dsRNA sequences within the 3'-UTR region of HIV-1. A small dsRNA fragment with 34 bases in size with a sequence of 3'-UUGGUACACGAGAUGGUUCGACUCGAUGAAGGGC-5' produced abundantly (9.17% of the total dsRNA fraction) by Bacillus subtilis MTCC5480 showed a much higher base complementarity values than previously reported miRNAs analysed against HIV-1. We separated the dsRNA fraction and validated the anti-HIV activity against human peripheral blood mononuclear cells (PBMC) infected with JRCSF strain of HIV-1 virus and the EC50 value ranges from 0.2-0.3 μM. This small dsRNA abundantly produced by B. subtilis could be studied further for its application as an anti-viral therapeutic agent.
Collapse
|
9
|
Detection and characterization of microRNA expression profiling and its target genes in response to canine parvovirus in Crandell Reese Feline Kidney cells. PeerJ 2020; 8:e8522. [PMID: 32095352 PMCID: PMC7023829 DOI: 10.7717/peerj.8522] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Accepted: 01/07/2020] [Indexed: 12/15/2022] Open
Abstract
Background MicroRNAs (miRNAs) play an essential role in gene regulators in many biological and molecular phenomena. Unraveling the involvement of miRNA as a key cellular factor during in vitro canine parvovirus (CPV) infection may facilitate the discovery of potential intervention candidates. However, the examination of miRNA expression profiles in CPV in tissue culture systems has not been fully elucidated. Method In the present study, we utilized high-throughput small RNA-seq (sRNA-seq) technology to investigate the altered miRNA profiling in miRNA libraries from uninfected (Control) and CPV-2c infected Crandell Reese Feline Kidney cells. Results We identified five of known miRNAs (miR-222-5p, miR-365-2-5p, miR-1247-3p, miR-322-5p and miR-361-3p) and three novel miRNAs (Novel 137, Novel 141 and Novel 102) by sRNA-seq with differentially expressed genes in the miRNA repertoire of CPV-infected cells over control. We further predicted the potential target genes of the aforementioned miRNAs using sequence homology algorithms. Notably, the targets of miR-1247-3p exhibited a potential function associated with cellular defense and humoral response to CPV. To extend the probing scheme for gene targets of miR-1247-3p, we explored and performed Gene Ontology (GO) enrichment analysis of its target genes. We discovered 229 putative targets from a total of 38 enriched GO terms. The top over-represented GO enrichment in biological process were lymphocyte activation and differentiation, marginal zone B cell differentiation, negative regulation of cytokine production, negative regulation of programed cell death, and negative regulation of signaling. We next constructed a GO biological process network composed of 28 target genes of miR-1247-3p, of which, some genes, namely BCL6, DLL1, GATA3, IL6, LEF1, LFNG and WNT1 were among the genes with obviously intersected in multiple GO terms. Conclusion The miRNA-1247-3p and its cognate target genes suggested their great potential as novel therapeutic targets or diagnostic biomarkers of CPV or other related viruses.
Collapse
|
10
|
Water-deficit responsive microRNAs in the primary root growth zone of maize. BMC PLANT BIOLOGY 2019; 19:447. [PMID: 31651253 PMCID: PMC6814125 DOI: 10.1186/s12870-019-2037-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 09/12/2019] [Indexed: 05/19/2023]
Abstract
BACKGROUND MicroRNA-mediated gene regulatory networks play a significant role in plant growth and development and environmental stress responses. RESULTS We identified 79 microRNAs (miRNAs) and multiple miRNA variants (isomiRs) belonging to 26 miRNA families in the primary root growth zone of maize seedlings grown at one of three water potentials: well-watered (- 0.02 MPa), mild water deficit stress (- 0.3 MPa), and severe water deficit stress (- 1.6 MPa). The abundances of 3 miRNAs (mild stress) and 34 miRNAs representing 17 families (severe stress) were significantly different in water-deficit stressed relative to well-watered controls (FDR < 0.05 and validated by stem loop RT-qPCR). Degradome sequencing revealed 213 miRNA-regulated transcripts and trancriptome profiling revealed that the abundance of 77 (miRNA-regulated) were regulated by water-defecit stress. miR399e,i,j-3p was strongly regulated by water-defcit stress implicating the possibility of nutrient deficiency during stress. CONCLUSIONS We have identified a number of maize miRNAs that respond to specific water deficits applied to the primary root growth zone. We have also identified transcripts that are targets for miRNA regulation in the root growth zone under water-deficit stress. The miR399e,i,j-3p that is known to regulate phosphate uptake in response to nutrient deficiencies responds to water-deficit stress, however, at the seedling stage the seed provides adequate nutrients for root growth thus miR399e,i,j-3p may play a separate role in water-deficit responses. A water-deficit regulated maize transcript, similar to known miR399 target mimics, was identified and we hypothesized that it is another regulatory player, moderating the role of miR399e,i,j-3p, in primary root growth zone water deficit responses.
Collapse
|
11
|
miRkwood: a tool for the reliable identification of microRNAs in plant genomes. BMC Genomics 2019; 20:532. [PMID: 31253093 PMCID: PMC6599362 DOI: 10.1186/s12864-019-5913-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 06/18/2019] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND MicroRNAs (miRNAs) play crucial roles in post-transcriptional regulation of eukaryotic gene expression and are involved in many aspects of plant development. Although several prediction tools are available for metazoan genomes, the number of tools dedicated to plants is relatively limited. RESULTS Here, we present miRkwood, a user-friendly tool for the identification of miRNAs in plant genomes using small RNA sequencing data. Deep-sequencing data of Argonaute associated small RNAs showed that miRkwood is able to identify a large diversity of plant miRNAs and limits false positive predictions. Moreover, it outperforms current tools such as ShortStack and contrary to ShortStack, miRkwood provides a quality score allowing users to rank miRNA predictions. CONCLUSION miRkwood is a very efficient tool for the annotation of miRNAs in plant genomes. It is available as a web server, as a standalone version, as a docker image and as a Galaxy tool: http://bioinfo.cristal.univ-lille.fr/mirkwood.
Collapse
|
12
|
Abstract
BACKGROUND MicroRNAs (miRNAs) play crucial roles in post-transcriptional regulation of eukaryotic gene expression and are involved in many aspects of plant development. Although several prediction tools are available for metazoan genomes, the number of tools dedicated to plants is relatively limited. RESULTS Here, we present miRkwood, a user-friendly tool for the identification of miRNAs in plant genomes using small RNA sequencing data. Deep-sequencing data of Argonaute associated small RNAs showed that miRkwood is able to identify a large diversity of plant miRNAs and limits false positive predictions. Moreover, it outperforms current tools such as ShortStack and contrary to ShortStack, miRkwood provides a quality score allowing users to rank miRNA predictions. CONCLUSION miRkwood is a very efficient tool for the annotation of miRNAs in plant genomes. It is available as a web server, as a standalone version, as a docker image and as a Galaxy tool: http://bioinfo.cristal.univ-lille.fr/mirkwood.
Collapse
|
13
|
Exploration of exosomal microRNA expression profiles in pigeon 'Milk' during the lactation period. BMC Genomics 2018; 19:828. [PMID: 30458711 PMCID: PMC6245878 DOI: 10.1186/s12864-018-5201-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2017] [Accepted: 10/25/2018] [Indexed: 11/10/2022] Open
Abstract
Background Pigeon crop has the unique ability to produce a nutrient rich substance termed pigeon ‘milk’ (PM), which has functional resemblance with the mammalian milk. Previous researches have demonstrated that a large number of exosomes and exosomal miRNAs exist in mammalian milk, and many of them are associated with immunity, growth and development. However, to date, little is known about the exosomes and exosomal miRNAs in PM. Results In this study, we isolated the exosomes from PM and used small RNA sequencing to investigate the distribution and expression profiles of exosomal miRNAs. A total of 301 mature miRNAs including 248 conserved and 53 novel miRNAs were identified in five lactation stages i.e. 1d, 5d, 10d, 15d, and 20d. From these, four top 10 conserved miRNAs (cli-miR-21-5p, cli-miR-148a-3p, cli-miR-10a-5p and cli-miR-26a-5p) were co-expressed in all five stages. We speculate that these miRNAs may have important role in the biosynthesis and metabolism of PM. Moreover, similar to the mammalian milk, a significant proportion of immune and growth-related miRNAs were also present and enriched in PM exosomes. Furthermore, we also identified 41 orthologous miRNAs group (giving rise to 81 mature miRNA) commonly shared with PM, human, bovine and porcine breast milk. Additionally, functional enrichment analysis revealed the role of exosomal miRNAs in organ development and in growth-related pathways including the MAPK, Wnt and insulin pathways. Conclusions To sum-up, this comprehensive analysis will contribute to a better understanding of the underlying functions and regulatory mechanisms of PM in squabs. Electronic supplementary material The online version of this article (10.1186/s12864-018-5201-0) contains supplementary material, which is available to authorized users.
Collapse
|
14
|
Small RNA-seq analysis of single porcine blastocysts revealed that maternal estradiol-17beta exposure does not affect miRNA isoform (isomiR) expression. BMC Genomics 2018; 19:590. [PMID: 30081835 PMCID: PMC6090871 DOI: 10.1186/s12864-018-4954-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Accepted: 07/23/2018] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND The expression of microRNAs (miRNAs) is essential for the proper development of the mammalian embryo. A maternal exposure to endocrine disrupting chemicals during preimplantation bears the potential for transgenerational inheritance of disease through the epigenetic perturbation of the developing embryo. A comprehensive assembly of embryo-specific miRNAs and respective isoforms (isomiR) is lacking to date. We aimed at revealing the sex-specific miRNA expression profile of single porcine blastocysts developing in gilts orally exposed to exogenous estradiol-17 (E2). Therefore we analyzed the miRNA profile specifically focusing on isomiRs and potentially embryo-specific miRNAs. RESULTS Deep sequencing of small RNA (small RNA-seq) result in the detection of miRNA sequences mapping to known and predicted porcine miRNAs as well as novel miRNAs highly conserved in human and cattle. A set of highly abundant miRNAs and a large number of rarely expressed miRNAs were identified by using a small RNA analysis pipeline, which was integrated into a novel Galaxy workflow specifically benefits incompletely annotated species. In particular, orthologue species information increased the total number of annotated miRNAs, while mapping to other non-coding RNAs avoided falsely annotated miRNAs. Neither the low nor the high dose of E2 treatment (10 and 1000 µ E2/kg body weight daily, respectively) affected the miRNA profile in blastocysts despite the distinct differential mRNA expression and DNA methylation found in previous studies. The high number of generated sequence reads enabled a comprehensive analysis of the isomiR repertoire showing various templated and non-templated modifications. Furthermore, potentially blastocyst-specific miRNAs were identified. CONCLUSIONS In pre-implantation embryos, numerous distinct isomiRs were discovered indicating a high complexity of miRNA expression. Neither the sex of the embryo nor a maternal E2 exposure affected the miRNA expression profile of developing porcine blastocysts. The adaptation to the continuous duration of the E2 treatment might explain the lack of an effect.
Collapse
|
15
|
miRge 2.0 for comprehensive analysis of microRNA sequencing data. BMC Bioinformatics 2018; 19:275. [PMID: 30153801 PMCID: PMC6112139 DOI: 10.1186/s12859-018-2287-y] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2018] [Accepted: 07/18/2018] [Indexed: 12/21/2022] Open
Abstract
Background miRNAs play important roles in the regulation of gene expression. The rapidly developing field of microRNA sequencing (miRNA-seq; small RNA-seq) needs comprehensive, robust, user-friendly and standardized bioinformatics tools to analyze these large datasets. We present miRge 2.0, in which multiple enhancements were made towards these goals. Results miRge 2.0 has become more comprehensive with increased functionality including a novel miRNA detection method, A-to-I editing analysis, integrated standardized GFF3 isomiR reporting, and improved alignment to miRNAs. The novel miRNA detection method uniquely uses both miRNA hairpin sequence structure and composition of isomiRs resulting in higher specificity for potential miRNA identification. Using known miRNA data, our support vector machine (SVM) model predicted miRNAs with an average Matthews correlation coefficient (MCC) of 0.939 over 32 human cell datasets and outperformed miRDeep2 and miRAnalyzer regarding phylogenetic conservation. The A-to-I editing detection strongly correlated with a reference dataset with adjusted R2 = 0.96. miRge 2.0 is the most up-to-date aligner with custom libraries to both miRBase v22 and MirGeneDB v2.0 for 6 species: human, mouse, rat, fruit fly, nematode and zebrafish; and has a tool to create custom libraries. For user-friendliness, miRge 2.0 is incorporated into bcbio-nextgen and implementable through Bioconda. Conclusions miRge 2.0 is a redesigned, leading miRNA RNA-seq aligner with several improvements and novel utilities. miRge 2.0 is freely available at: https://github.com/mhalushka/miRge. Electronic supplementary material The online version of this article (10.1186/s12859-018-2287-y) contains supplementary material, which is available to authorized users.
Collapse
|
16
|
Elimination of PCR duplicates in RNA-seq and small RNA-seq using unique molecular identifiers. BMC Genomics 2018; 19:531. [PMID: 30001700 PMCID: PMC6044086 DOI: 10.1186/s12864-018-4933-1] [Citation(s) in RCA: 87] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2018] [Accepted: 07/08/2018] [Indexed: 12/12/2022] Open
Abstract
Background RNA-seq and small RNA-seq are powerful, quantitative tools to study gene regulation and function. Common high-throughput sequencing methods rely on polymerase chain reaction (PCR) to expand the starting material, but not every molecule amplifies equally, causing some to be overrepresented. Unique molecular identifiers (UMIs) can be used to distinguish undesirable PCR duplicates derived from a single molecule and identical but biologically meaningful reads from different molecules. Results We have incorporated UMIs into RNA-seq and small RNA-seq protocols and developed tools to analyze the resulting data. Our UMIs contain stretches of random nucleotides whose lengths sufficiently capture diverse molecule species in both RNA-seq and small RNA-seq libraries generated from mouse testis. Our approach yields high-quality data while allowing unique tagging of all molecules in high-depth libraries. Conclusions Using simulated and real datasets, we demonstrate that our methods increase the reproducibility of RNA-seq and small RNA-seq data. Notably, we find that the amount of starting material and sequencing depth, but not the number of PCR cycles, determine PCR duplicate frequency. Finally, we show that computational removal of PCR duplicates based only on their mapping coordinates introduces substantial bias into data analysis. Electronic supplementary material The online version of this article (10.1186/s12864-018-4933-1) contains supplementary material, which is available to authorized users.
Collapse
|
17
|
miR-MaGiC improves quantification accuracy for small RNA-seq. BMC Res Notes 2018; 11:296. [PMID: 29764489 PMCID: PMC5952827 DOI: 10.1186/s13104-018-3418-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2018] [Accepted: 05/09/2018] [Indexed: 12/17/2022] Open
Abstract
Objective Many tools have been developed to profile microRNA (miRNA) expression from small RNA-seq data. These tools must contend with several issues: the small size of miRNAs, the small number of unique miRNAs, the fact that similar miRNAs can be transcribed from multiple loci, and the presence of miRNA isoforms known as isomiRs. Methods failing to address these issues can return misleading information. We propose a novel quantification method designed to address these concerns. Results We present miR-MaGiC, a novel miRNA quantification method, implemented as a cross-platform tool in Java. miR-MaGiC performs stringent mapping to a core region of each miRNA and defines a meaningful set of target miRNA sequences by collapsing the miRNA space to “functional groups”. We hypothesize that these two features, mapping stringency and collapsing, provide more optimal quantification to a more meaningful unit (i.e., miRNA family). We test miR-MaGiC and several published methods on 210 small RNA-seq libraries, evaluating each method’s ability to accurately reflect global miRNA expression profiles. We define accuracy as total counts close to the total number of input reads originating from miRNAs. We find that miR-MaGiC, which incorporates both stringency and collapsing, provides the most accurate counts. Electronic supplementary material The online version of this article (10.1186/s13104-018-3418-2) contains supplementary material, which is available to authorized users.
Collapse
|
18
|
Deep sequencing and high-throughput analysis of PIWI-associated small RNAs. Methods 2017; 126:66-75. [PMID: 28552266 DOI: 10.1016/j.ymeth.2017.05.020] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2017] [Revised: 05/08/2017] [Accepted: 05/21/2017] [Indexed: 02/07/2023] Open
Abstract
Small RNAs are now known to be major regulatory factors of gene expression. Emerging methods based on deep-sequencing have enabled the analysis of small RNA expression in a high-throughput manner, leading to the identification of large numbers of small RNAs in various species. Moreover, profiling small RNA data together with transcriptome data enables transcriptional and post-transcriptional regulation mediated by small RNAs to be hypothesized. Here, we isolated PIWIL1 (MIWI)-associated small RNAs from mouse testes, and performed small RNA-seq analysis. In addition, directional RNA-seq was performed using Piwil1 mutant mouse testes. Using these data, we describe protocols for analyzing small RNA-seq reads to obtain profiles of small RNAs associated with PIWI proteins. We also present bioinformatic protocols for analyzing RNA-seq reads that aim to annotate expression of piRNA clusters and identify genes regulated by piRNAs.
Collapse
|
19
|
Functional Impact of RNA editing and ADARs on regulation of gene expression: perspectives from deep sequencing studies. Cell Biosci 2014; 4:44. [PMID: 25949793 PMCID: PMC4422215 DOI: 10.1186/2045-3701-4-44] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2014] [Accepted: 07/14/2014] [Indexed: 11/13/2022] Open
Abstract
Cells regulate gene expression at multiple levels leading to a balance between robustness and complexity within their proteome. One core molecular step contributing to this important balance during metazoan gene expression is RNA editing, such as the co-transcriptional recoding of RNA transcripts catalyzed by the adenosine deaminse acting on RNA (ADAR) family of enzymes. Understanding of the adenosine-to-inosine RNA editing process has been broadened considerably by the next generation sequencing (NGS) technology, which allows for in-depth demarcation of an RNA editome at nucleotide resolution. However, critical issues remain unresolved with regard to how RNA editing cooperates with other transcript-associated events to underpin regulated gene expression. Here we review the growing body of evidence, provided by recent NGS-based studies, that links RNA editing to other mechanisms of post-transcriptional RNA processing and gene expression regulation including alternative splicing, transcript stability and localization, and the biogenesis and function of microRNAs (miRNAs). We also discuss the possibility that systematic integration of NGS data may be employed to establish the rules of an “RNA editing code”, which may give us new insights into the functional consequences of RNA editing.
Collapse
|