Small RNA Transcriptome Research Articles

Profiling the cell-specific small non-coding RNA transcriptome of the human placenta

Profiling the cell-specific small non-coding RNA transcriptome of the human placenta

Sponging of five tumour suppressor miRNAs by lncRNA-KCNQ1OT1 activates BMPR1A/BMPR1B-ACVR2A/ACVR2B signalling and promotes chemoresistance in hepatocellular carcinoma

Diverse mechanisms have been established to understand the chemoresistance of hepatocellular carcinoma (HCC), but the contribution of non-coding RNAs is not surveyed well. Here, we aimed to explore the lncRNA-miRNA axis in Hepatitis C and B virus (HCV and HBV) infected HCC to investigate the molecular mechanism of chemoresistance and to identify a potential therapeutic target for HCC. The small RNA transcriptome analysis followed by qRT-PCR validation with the liver tissues of both HCV and HBV infected HCC patients revealed that miR-424-5p, miR-136-3p, miR-139-5p, miR-223-3p, and miR-375-3p were the most downregulated miRNAs in HCC compared to normal (log2 fold change ≤−1.5, Padj ≤ 0.05). In silico pathway analysis with the validated targets of each miRNA revealed that the signalling pathway regulating pluripotency of stem cells is commonly targeted by these five miRNAs. Subsequent validation by 3′UTR-luciferase assay and western blot analysis unveiled that these five miRNAs impeded either same or diverse genes, but all linked to BMP signalling pathway such as BMPR1A/BMPR1B by miR-139-5p, miR-136-3p, and miR-375-3p, and ACVR2A/ACVR2B by miR-424-5p and miR-223-3p. Furthermore, restoration of each miRNA in Huh7/SNU449 cells inhibited phosphorylation of downstream SMAD1/5 and ERK1/2, and attenuated Epithelial-mesenchymal transition, stemness, spheroid formation, chemoresistance, invasion and migration of cells. To investigate the mechanism of suppression of these miRNAs, “DIANA” tool was employed and lncRNA-KCNQ1OT1 was retrieved as interacting partner of all the five miRNAs. In vitro RNA pull-down assay revealed that lncRNA-KCNQ1OT1 physically interacted and sequestered these five miRNAs in the cytoplasm. Hence, KCNQ1OT1 was suppressed in Huh7/SNU449 cells using CRISPR technology and observed regression of oncogenic properties with enhanced chemosensitivity and reduced metastasis in cancer cells. Shrinkage of tumour size and volume in NOD-SCID mice injected with KCNQ1OT1-sgRNA cells further strengthened our observations. Thus, lncRNA-KCNQ1OT1 is the main regulator, which reduces the level of beneficiary miRNAs in the tumour milieu and modulates BMP signalling pathway to promote chemoresistance in HCC, suggesting lncRNA-KCNQ1OT1 might have robust potential to be a therapeutic target in HCC.

Genome of tropical bed bug Cimex hemipterus (Cimicidae, Hemiptera) reveals tetraspanin expanded in bed bug ancestor.

Cimex species are ectoparasites that exclusively feed on warm-blooded animals such as birds and mammals. Three cimicid species are known to be persistent pests for humans, including the tropical bed bug Cimex hemipterus, common bed bug Cimex lectularius, and Eastern bat bug Leptocimex boueti. To date, genomic information is restricted to the common bed bug C. lectularius, which limits understanding their biology and to provide controls of bed bug infestations. Here, a chromosomal-level genome assembly of C. hemipterus (495 Mb [megabase pairs]) contained on 16 pseudochromosomes (scaffold N50 = 34 Mb), together with 9 messenger RNA and small RNA transcriptomes were obtained. In comparison between hemipteran genomes, we found that the tetraspanin superfamily was expanded in the Cimex ancestor. This study provides the first genome assembly for the tropical bed bug C. hemipterus, and offers an unprecedented opportunity to address questions relating to bed bug infestations, as well as genomic evolution to hemipterans more widely.

TsMS combined with EA promotes functional recovery and axonal regeneration via mediating the miR-539-5p/Sema3A/PlexinA1 signalling axis in sciatic nerve-injured rats

Enhancing axonal regeneration is one of the most important processes in treating nerve injuries. Both magnetic and electrical stimulation have the effect of promoting nerve axon regeneration. But few study has investigated the effects of trans-spinal magnetic stimulation (TsMS) combined with electroacupuncture (EA) on nerve regeneration in rats with sciatic nerve injury. In this study, we compared the improvement of neurological function in rats with sciatic nerve crush injuries after 4 weeks of different interventions (EA, TsMS, or TsMS combined with EA). We further explored the morphological and molecular biological alterations following sciatic nerve injury by HE, Masson, RT-PCR, western blotting, immunofluorescence staining and small RNA transcriptome sequencing. The results showed that TsMS combined with EA treatment significantly promoted axonal regeneration, increased the survival rate of neurons, and suppressed denervation atrophy of the gastrocnemius muscle. Subsequent experiments suggested that the combination treatment may play an active role by mediating the miR-539-5p/Sema3A/PlexinA1 signaling axis.

Viral and host small RNA transcriptome analysis of SARS-CoV-1 and SARS-CoV-2-infected human cells reveals novel viral short RNAs

RNA viruses have been shown to express various short RNAs, some of which have regulatory roles during replication, transcription, and translation of viral genomes. However, short viral RNAs generated from SARS-CoV-1 and SARS-CoV-2 genomic RNAs remained largely unexplored, possibly due limitations of the widely used library preparation methods for small RNA deep sequencing and corresponding data processing. By analyzing publicly available small RNA sequencing datasets, we observed that human Calu-3 cells infected by SARS-CoV-1 or SARS-CoV-2 accumulate multiple previously unreported short viral RNAs. In addition, we verified the presence of the five most abundant SARS-CoV-2 short viral RNAs in SARS-CoV-2-infected human lung adenocarcinoma cells by quantitative PCR. Interestingly, the copy number of the observed SARS-CoV-2 short viral RNAs dramatically exceeded the expression of previously reported viral microRNAs in the same cells. We hypothesize that the reported SARS-CoV-2 short viral RNAs could serve as biomarkers for early infection stages due to their high abundance. Furthermore, unlike SARS-CoV-1, the SARS-CoV-2 infection induced significant (Benjamini-Hochberg-corrected p-value <0.05) deregulation of Y-RNA, transfer RNA, vault RNA, as well as more than 300 endogenous short RNAs that aligned predominantly to human protein-coding and long noncoding RNA transcripts. In particular, more than 20-fold upregulation of reads derived from Y-RNA (and several transfer RNAs) have been documented in RNA-seq datasets from SARS-CoV-2 infected cells. Finally, a significant proportion of short RNAs derived from full-length viral genomes also aligned to various human genome (hg38) sequences, suggesting opportunities to investigate regulatory roles of short viral RNAs during infection. Further characterization of the small RNA landscape of both viral and host genomes is clearly warranted to improve our understanding of molecular events related to infection and to design more efficient strategies for therapeutic interventions as well as early diagnosis.

Unraveling the regulatory network of miRNA expression in Potato Y virus-infected of Nicotiana benthamiana using integrated small RNA and transcriptome sequencing.

Potato virus Y (PVY) disease is a global problem that causes significant damage to crop quality and yield. As traditional chemical control methods are ineffective against PVY, it is crucial to explore new control strategies. MicroRNAs (miRNAs) play a crucial role in plant and animal defense responses to biotic and abiotic stresses. These endogenous miRNAs act as a link between antiviral gene pathways and host immunity. Several miRNAs target plant immune genes and are involved in the virus infection process. In this study, we conducted small RNA sequencing and transcriptome sequencing on healthy and PVY-infected N. benthamiana tissues (roots, stems, and leaves). Through bioinformatics analysis, we predicted potential targets of differentially expressed miRNAs using the N. benthamiana reference genome and the PVY genome. We then compared the identified differentially expressed mRNAs with the predicted target genes to uncover the complex relationships between miRNAs and their targets. This study successfully constructed a miRNA-mRNA network through the joint analysis of Small RNA sequencing and transcriptome sequencing, which unveiled potential miRNA targets and identified potential binding sites of miRNAs on the PVY genome. This miRNA-mRNA regulatory network suggests the involvement of miRNAs in the virus infection process.

Untying the regulatory roles of miRNAs in CuO-NPs stress response mechanism in maize: A genome-wide sRNA transcriptome analysis

MicroRNAs (miRNAs), a group of tiny non-coding RNAs play pivotal role in plant responses to environmental stress. The present small RNA transcriptome study aims to untie the role of miRNAs in CuO-NPs stress adaptation in maize seedlings. Restricted seedling growth, enhanced ROS generation and higher membrane damage were recorded under CuO-NPs [<50 nm, 8 mM] treatment. Deep sequencing reveals 7 up- and 36 down-regulated known miRNAs from CuO-NPs challenged leaves. Gene ontology study demonstrates involvement of CuO-NPs responsive miRNAs in a variety of biological processes including plant growth (miR159a, miR159b), redox homeostasis (miR156e, miR395a), detoxification of heavy metals (miR156e, miR827), signal transduction (miR156e, miR156d), and cell signalling (miR167b-3p, miR393a). Enhanced transcriptional abundance of ABC transporter G family member 41 isoform X2 and HM-associated isoprenylated plant protein 45 isoform X1 might be involved in sequestration and detoxification of excess Cu, essential for metal homeostasis in maize. The miR528-5p mediated up-regulation of superoxide dismutase does not give much protection against CuO-NPs induced oxidative stress damages as evident after histochemical staining with NBT. Moreover, CuO-NPs stress mediated down regulation of miR396 could be an underlying cause of the restricted seedling growth. Taken together, our findings provide insights into the miRNA-guided stress regulatory networks involved in plant's adaptive responses to CuO-NPs stress.

Integrated analysis of miRNA profiles and gut bacterial changes in Altica viridicyanea following antibiotic treatment.

The gut bacteria involves in insect homeostasis by playing essential roles in host physiology, metabolism, innate immunity, and so forth. microRNAs (miRNAs) are endogenous small noncoding RNAs that posttranscriptionally regulate gene expression to affect immune or metabolic processes in insects. For several non-model insects, the available knowledge on the relationship between changes in the gut bacteria and miRNA profiles is limited. In this study, we investigated the gut bacterial diversity, composition, and function from Altica viridicyanea feeding on normal- and antibiotic-treated host plants using 16S rRNA amplicon sequencing; antibiotics have been shown to affect the body weight and development time in A. viridicyanea, suggesting that the gut bacteria of the normal sample were more diverse and abundant than those of the antibiotic-fed group, and most of them were involved in various physical functions by enrichment analysis. Furthermore, we executed small RNA transcriptome sequencing using the two experimental groups to obtain numerous sRNAs, such as piRNAs, siRNAs, and known and novel miRNAs, by data mapping and quality control, and furthermore, a total of 224 miRNAs were identified as significantly differentially expressed miRNAs, of which some DEMs and their target genes participated in immune- and metabolism-related pathways based on GO and KEGG annotation. Besides, regarding the regulatory roles of miRNA and target genes, a interaction network of DEM-target gene pairs from eight immune- or metabolism-related signaling pathways were constructed. Finally, we discovered that DEMs from above pathways were significantly positively or negatively correlated with gut bacterial alterations following antibiotic treatment. Collectively, the observations of this study expand our understanding of how the disturbance of gut bacteria affects miRNA profiles in A. viridicyanea and provide new valuable resources from extreme ranges for future studies on the adaptive evolution in insects.

Long-term iTBS improves neural functional recovery by reducing the inflammatory response and inhibiting neuronal apoptosis via miR-34c-5p/p53/Bax signaling pathway in cerebral ischemic rats.

To investigate intermittent theta-burst stimulation (iTBS) effect on ischemic stroke and the underlying mechanism of neurorehabilitation, we developed an ischemia/reperfusion (I/R) injury model in Sprague-Dawley (SD) rats using the middle cerebral artery occlusion/reperfusion (MCAO/r) method. Next, using different behavioral studies, we compared the improvement of the whole organism with and without iTBS administration for 28 days. We further explored the morphological and molecular biological alterations associated with neuronal apoptosis and neuroinflammation by TTC staining, HE staining, Nissl staining, immunofluorescence staining, ELISA, small RNA sequencing, RT-PCR, and western blot assays. The results showed that iTBS significantly protected against neurological deficits and neurological damage induced by cerebral I/R injury. iTBS also significantly decreased brain infarct volume and increased the number of surviving neurons after 28 days. Additionally, it was observed that iTBS decreased synaptic loss, suppressed activation of astrocytes and M1-polarized microglia, and simultaneously promoted M2-polarized microglial activation. Furthermore, iTBS intervention inhibited neuronal apoptosis and exerted a positive impact on the neuronal microenvironment by reducing neuroinflammation in cerebral I/R injured rats. To further investigate the iTBS mechanism, this study was conducted using small RNA transcriptome sequencing of various groups of peri-infarcted tissues. Bioinformatics analysis and RT-PCR discovered the possible involvement of miR-34c-5p in the mechanism of action. The target genes prediction and detection of dual-luciferase reporter genes confirmed that miR-34c-5p could inhibit neuronal apoptosis in cerebral I/R injured rats by regulating the p53/Bax signaling pathway. We also confirmed by RT-PCR and western blotting that miR-34c-5p inhibited Bax expression. In conclusion, our study supports that iTBS is vital in inhibiting neuronal apoptosis in cerebral I/R injured rats by mediating the miR-34c-5p involvement in regulating the p53/Bax signaling pathway.

Multiomics studies with co-transformation reveal microRNAs via miRNA-TF-mRNA network participating in wood formation in Hevea brasiliensis.

MicroRNAs (miRNAs) are small endogenous non-coding RNAs that play an important role in wood formation in plants. However, the significance of the link between miRNAs and their target transcripts in wood formation remains unclear in rubber tree (Hevea brasiliensis). In this study, we induced the formation of reaction wood by artificially bending rubber trees for 300 days and performed small RNA sequencing and transcriptome deep sequencing (RNA-seq) to describe the complement of miRNAs and their targets contributing to this process. We identified 5, 11, and 2 differentially abundant miRNAs in normal wood (NW) compared to tension wood (TW), in NW relative to opposite wood (OW), and between TW and OW, respectively. We also identified 12 novel miRNAs and 39 potential miRNA-mRNA pairs with different accumulation patterns in NW, TW, and OW. We noticed that many miRNAs targeted transcription factor genes, which were enriched in KEGG pathways associated with phenylpropanoid biosynthesis, phenylalanine metabolism, and pyruvate metabolism. Thus, miRNA-TF-mRNA network involved in wood formation via tension wood model were constructed. We validated the differential accumulation of miRNAs and their targets by RT-qPCR analysis and overexpressed miRNA in Nicotiana benthamiana with its potential target gene. These results will provide a reference for a deep exploration of growth and development in rubber tree.

Time-resolved small RNA transcriptomics of the ichthyosporean Sphaeroforma arctica.

Ichthyosporea, a clade of holozoans, represent a clade closely related to animals, and thus hold a key phylogenetic position for understanding the origin of animals. We have previously discovered that an ichthyosporean, Sphaeroforma arctica, contains microRNAs (miRNAs) as well as the miRNA processing machinery. This was the first discovery of miRNAs among the closest single-celled relatives of animals and raised intriguing questions about the roles of regulatory small RNAs in cell development and differentiation in unicellular eukaryotes. Like many ichthyosporeans, S. arctica also undergoes a transient multicellular developmental life cycle. As miRNAs are, among other roles, key regulators of gene expression during development in animals, we wanted to investigate the dynamics of miRNAs during the developmental cycle in S. arctica. Here we have therefore collected a comprehensive time-resolved small RNA transcriptome linked to specific life stages with a substantially higher sequencing depth than before, which can enable further discovery of functionally relevant small RNAs. The data consists of Illumina-sequenced small RNA libraries from two independent biological replicates of the entire life cycle of S. arctica with high temporal resolution. The dataset is directly linked and comes from the same samples as a previously published mRNA-seq dataset, thus enabling direct cross-functional analyses.

Small RNA transcriptome analysis using parallel single-cell small RNA sequencing

miRNA and other forms of small RNAs are known to regulate many biological processes. Single-cell small RNA sequencing can be used to profile small RNAs of individual cells; however, limitations of efficiency and scale prevent its widespread application. Here, we developed parallel single-cell small RNA sequencing (PSCSR-seq), which can overcome the limitations of existing methods and enable high-throughput small RNA expression profiling of individual cells. Analysis of PSCSR-seq data indicated that diverse cell types could be identified based on patterns of miRNA expression, and showed that miRNA content in nuclei is informative (for example, cell type marker miRNAs can be detected in isolated nuclei). PSCSR-seq is very sensitive: analysis of only 732 peripheral blood mononuclear cells (PBMCs) detected 774 miRNAs, whereas bulk small RNA analysis would require input RNA from approximately 106 cells to detect as many miRNAs. We identified 42 miRNAs as markers for PBMC subpopulations. Moreover, we analyzed the miRNA profiles of 9,533 cells from lung cancer biopsies, and by dissecting cell subpopulations, we identified potentially diagnostic and therapeutic miRNAs for lung cancers. Our study demonstrates that PSCSR-seq is highly sensitive and reproducible, thus making it an advanced tool for miRNA analysis in cancer and life science research.

NGS-based profiling identifies miRNAs and pathways dysregulated in cisplatin-resistant esophageal cancer cells.

Esophageal cancer (EC) incidence remains to be on a global rise supported by an unchanged recurrence and 5-year survival rate owing to the development of chemoresistance. Resistance to cisplatin, one of the majorly used chemotherapeutic drugs in EC, is a major nuisance. This study sheds light on miRNA dysregulation and its inverse relation with dysregulated mRNAs to guide pathways into the manifestation of cisplatin resistance in EC. A cisplatin-resistant version of an EC cell line was established and comparative profiling by NGS with the parental cell line was employed to identify dysregulation in miRNA and mRNA levels. Protein-protein interaction network analysis was done using Cytoscape, followed by Funrich pathway analysis. Furthermore, selective significant miRNAs were validated using qRT-PCR. miRNA-mRNA integrated analysis was carried out using the Ingenuity Pathway Analysis (IPA) tool. Expression of various established resistance markers supported the successful establishment of cisplatin-resistant cell line. Whole-cell small RNA sequencing and transcriptome sequencing identified 261 miRNAs and 1892 genes to be significantly differentially expressed (DE), respectively. Pathway analysis indicated enrichment of EMT signaling, supported by NOTCH, mTOR, TNF receptor, and PI3K-mediated AKT signaling pathways, in chemoresistant cells. Validation by qRT-PCR confirmed upregulation of miR-10a-5p, miR-618, miR-99a-5p, and miR-935 and downregulation of miR-335-3p, miR-205-5p, miR-944, miR-130a-3p, and miR-429 in resistant cells. Pathway analysis that followed IPA analysis indicated that the dysregulation of these miRNAs and their target genes may be instrumental in the development and regulation of chemoresistance via p53 signaling, xenobiotic metabolism, and NRF2-mediated oxidative stress. This study concludes the interplay between miRNA and mRNA as an important aspect and occurrence in guiding the regulation, acquisition, and maintenance of chemoresistance in esophageal cancer in vitro.

Comparison of Methods for MicroRNA Isolation from Extracellular Vesicles Obtained from Ascitic Fluids.

Secreted extracellular vesicles (EVs) contain active biomolecules, including miRNAs, composition of which reflects epigenetic changes occurring in cells during pathological processes, in particular, malignant transformation. The accumulated pool of data on the role of EVs in carcinogenesis has stimulated investigations of the EV-derived cancer markers. The most important factor limiting development of this scientific direction is lack of "gold standards" both for methods of EV isolation from biological fluids and for analyzing their molecular content, including composition of miRNAs. Here we first examined efficacy of various methods for small RNA isolation from EVs contained in ascitic fluid for subsequent miRNA analysis. Comparison of different commercial kits showed advantages of the methods based on phenol-chloroform extraction: Total Exosome RNA & Protein Isolation Kit and miRNeasy Serum/Plasma Kit. Analysis of the small RNA transcriptome showed presence of various classes of molecules in the EVs, among which proportion of miRNAs averaged 6% and reaching 10% with the Total Exosome RNA & Protein Isolation Kit. The PureLink miRNA Isolation Kit demonstrated the lowest efficiency. The miRNeasy Advanced Serum/Plasma Kit showed the highest concentration of the small RNA fraction, miRNA proportion of which, however, did not exceed that obtained with the miRNeasy Serum/Plasma Kit and Total Exosome RNA & Protein Isolation Kit. Moreover, RT-PCR analysis of the individual molecules showed lower levels of each of investigated miRNAs (miR-1246, miR-200b-5p, miR-200c-3p, and miR-23a-3p) when using the miRNeasy Advanced Serum/Plasma Kit. In conclusion, Total Exosome RNA & Protein Isolation Kit and miRNeasy Serum/Plasma Kit can be considered as optimal kits in terms of performance based on combination of the studied characteristics, including small RNA concentration, percentage of microRNA according to bioanalyzer and sequencing results, and levels of individual miRNAs detected by RT-PCR.

Comparative transcriptomics reveals small RNA composition and differential microRNA responses underlying interferon-mediated antiviral regulation in porcine alveolar macrophages.

Previous studies have shown that interferon-mediated antiviral activity is subtype-dependent. Using a whole transcriptome procedure, we aimed to characterize the small RNA transcriptome (sRNA-Seq) and specifically the differential microRNA (miRNA) responses in porcine alveolar macrophages (PAMs) upon antiviral activation during viral infection and interferon (IFN) stimulation. Data showed that near 90% of the qualified reads of sRNA were miRNAs, and about 10% of the other sRNAs included rRNA, snoRNA, snRNA, and tRNA in order of enrichment. As the majority of sRNA (>98%) were commonly detected in all PAM samples under different treatments, about 2% sRNA were differentially expressed between the different antiviral treatments. Focusing on miRNA, 386 miRNA were profiled, including 331 known and 55 novel miRNA sequences, of which most were ascribed to miRNA families conserved among vertebrates, particularly mammalian species. Of the miRNA profiles comparably generated across the different treatments, in general, significantly differentially expressed miRNA (SEM) demonstrated that: (1) the wild-type and vaccine strains of a porcine arterivirus (a.k.a., PRRSV) induced nearly reversed patterns of up- or down-regulated SEMs; (2) similar SEM patterns were found among the treatments by the vaccine strain and antiviral IFN-α1/-ω5 subtypes; and (3) the weak antiviral IFN-ω1, however, remarked a suppressive SEM pattern as to SEMs upregulated in the antiviral treatments by the vaccine and IFN-α1/-ω5 subtypes. Further articulation identified SEMs commonly or uniquely expressed in different treatments, and experimentally validated that some SEMs including miR-10b and particularly miR-9-1 acted significantly in regulation of differential antiviral reactions stimulated by different IFN subtypes. Therefore, this study provides a general picture of porcine sRNA composition and pinpoints key SEMs underlying antiviral regulation in PAMs correlated to a typical respiratory RNA virus in pigs.

Screening and identification of miRNAs regulating Tbx4/5 genes of Pampus argenteus.

BackgroundSilver pomfret (Pampus argenteus) is one of the most widely distributed and economically important pelagic fish species. However, an unique morphological feature of P. argenteus is the loss of pelvic fins, which can increase the energy requirement during food capture to some extent and is therefore not conducive to artificial culture. Tbx4/5 genes are highly conserved regulatory factors that regulate limb development in vertebrates and are in turn regulated by microRNAs (miRNAs). However, the miRNAs that directly regulate the Tbx4/5 genes in P. argenteus remain to be elucidated.MethodsThe Tbx4/5 genes of P. argenteus were first cloned, and the small RNA transcriptomes were sequenced by high-throughput sequencing during the critical period of the fin development at days 1, 7, and 13 of hatching. The miRNAs regulating the Tbx4/5 genes of P. argenteus were subsequently predicted by bioinformatics analysis, and the related miRNAs were verified in vitro using a dual fluorescence reporter system.ResultsA total of 662 miRNAs were identified, of which 257 were known miRNAs and 405 were novel miRNAs were identified. Compared to day 1, 182 miRNAs were differentially expressed (DE) on day 7, of which 77 and 105 miRNAs were downregulated and upregulated, respectively, while 278 miRNAs were DE on day 13, of which 136 and 142 miRNAs were downregulated and upregulated, respectively. Compared to day 13, four miRNAs were DE on day 7, of which three miRNAs were downregulated and one miRNA was upregulated. The results of hierarchical clustering of the miRNAs revealed that the DE genes were inversely expressed between days 1 and 7, and between days 1 and 13 of larval development, indicating that the larvae were in the peak stage of differentiation. However, the number of DE genes between days 7 and 13 of larval development was relatively small, suggesting the initiation of development. The potential target genes of the DE miRNAs were subsequently predicted, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of target genes were performed. The results suggested that the DE miRNAs were involved in growth, development, and signal transduction pathways, of which the Wnt and Fgfs signaling pathways are known to play important roles in the growth and development of fins. The results of dual fluorescence reporter assays demonstrated that miR-102, miR-301c, and miR-589 had a significant negative regulatory effect on the 3′-UTR of the Tbx4 gene, while miR-187, miR-201, miR-219, and miR-460 had a significant negative regulatory effect on the 3′-UTR of the Tbx5 gene. Altogether, the findings indicated that miRNAs play an important role in regulating the growth and development of pelvic fins in P. argenteus. This study provides a reference for elucidating the interactions between the miRNAs and target genes of P. argenteus in future studies.

The role of transposon inverted repeats in balancing drought tolerance and yield-related traits in maize.

The genomic basis underlying the selection for environmental adaptation and yield-related traits in maize remains poorly understood. Here we carried out genome-wide profiling of the small RNA (sRNA) transcriptome (sRNAome) and transcriptome landscapes of a global maize diversity panel under dry and wet conditions and uncover dozens of environment-specific regulatory hotspots. Transgenic and molecular studies of Drought-Related Environment-specific Super eQTL Hotspot on chromosome 8 (DRESH8) and ZmMYBR38, a target of DRESH8-derived small interfering RNAs, revealed a transposable element-mediated inverted repeats (TE-IR)-derived sRNA- and gene-regulatory network that balances plant drought tolerance with yield-related traits. A genome-wide scan revealed that TE-IRs associate with drought response and yield-related traits that were positively selected and expanded during maize domestication. These results indicate that TE-IR-mediated posttranscriptional regulation is a key molecular mechanism underlying the tradeoff between crop environmental adaptation and yield-related traits, providing potential genomic targets for the breeding of crops with greater stress tolerance but uncompromised yield.

Genome-wide profiling of dysregulated piRNAs and their target genes implicated in oncogenicity of tongue squamous cell carcinoma

PIWI-interacting RNAs (piRNAs) are single-stranded, 23–36 nucleotide long RNAs that regulate gene expression in the germline but are also detected in some cancers. However, there are no reports yet on piRNA expression in tongue squamous cell carcinoma (TSCC), the most common oral cancer (80–90% percent of all oral cancers). We performed small RNA and whole transcriptome sequencing in H357 tongue cancer and HOK cells (GEO database accession numbers: GSE196674 and GSE196688). We also examined nine published sets of gene expression array data of TSCC tissues from the GEO database to decode piRNAs and their putative targets that may be involved in tumorigenesis. We identified a pool of 16,058 and 25,677 piRNAs in H357 and HOK, respectively, among which 406 are differentially expressed. We also found that 2094 protein-coding genes are differentially expressed in either TSCC tissues or cell lines. We performed target predictions for these piRNAs, pathway and disease function (DF) analyses, as well as qRT-PCR validation of piRNA-target pairs. These experiments revealed one up-regulated (FDFT1) and four down-regulated (OGA, BDH1, TAT, HYAL4) target genes that are enriched in 11 canonical pathways (CPs), with postulated roles in the initiation and progression of TSCC. Downregulation of piR-33422 is predicted to upregulate the FDFT1 gene, which encodes a mevalonate/cholesterol-pathway related farnesyl-diphosphate farnesyltransferase. The FDFT1 appears to be involved in the largest number of oncogenesis-related processes and is interacting with statins, which is a classical cancer drug. This study provides the first evidence of the piRNome of TSCC, which could be investigated further to decode piRNA-mediated gene regulations in malignancy and potential drug targets, such as FDFT1.

Small and long RNA transcriptome of whole human cerebrospinal fluid and serum as compared to their extracellular vesicle fractions reveal profound differences in expression patterns and impacts on biological processes

BackgroundNext generation sequencing (NGS) of human specimen is expected to improve prognosis and diagnosis of human diseases, but its sensitivity urges for well-defined sampling and standardized protocols in order to avoid error-prone conclusions.MethodsIn this study, large volumes of pooled human cerebrospinal fluid (CSF) were used to prepare RNA from human CSF-derived extracellular vesicles (EV) and from whole CSF, as well as from whole human serum and serum-derived EV. In all four fractions small and long coding and non-coding RNA expression was analyzed with NGS and transcriptome analyses.ResultsWe show, that the source of sampling has a large impact on the acquired NGS pattern, and differences between small RNA fractions are more distinct than differences between long RNA fractions. The highest percentual discrepancy between small RNA fractions and the second highest difference between long RNA fractions is seen in the comparison of CSF-derived EV and whole CSF. Differences between miR (microRNA) and mRNA fractions of EV and the respective whole body fluid have the potential to affect different cellular and biological processes. I.e. a comparison of miR in both CSF fractions reveals that miR from EV target four transcripts sets involved in neurobiological processes, whereas eight others, also involved in neurobiological processes are targeted by miR found in whole CSF only. Likewise, three mRNAs sets derived from CSF-derived EV are associated with neurobiological and six sets with mitochondrial metabolism, whereas no such mRNA transcript sets are found in the whole CSF fraction. We show that trace amounts of blood-derived contaminations of CSF can bias RNA-based CSF diagnostics.ConclusionsThis study shows that the composition of small and long RNA differ significantly between whole body fluid and its respective EV fraction and thus can affect different cellular and molecular functions. Trace amounts of blood-derived contaminations of CSF can bias CSF analysis. This has to be considered for a meaningful RNA-based diagnostics. Our data imply a transport of EV from serum to CSF across the blood–brain barrier.

Organization of minicircle cassettes and guide RNA genes in Trypanosoma brucei.

Mitochondrial DNA of protists of order Kinetoplastida comprises thousands of interlinked circular molecules arranged in a network. There are two types of molecules called minicircles and maxicircles. Minicircles encode guide RNA (gRNA) genes whose transcripts mediate post-transcriptional editing of maxicircle encoded genes. Minicircles are diverse. The human sleeping sickness parasite Trypanosoma brucei has one of the most diverse sets of minicircle classes of all studied trypanosomatids with hundreds of different classes, each encoding one to four genes mainly within cassettes framed by 18 bp inverted repeats. A third of cassettes have no identifiable gRNA genes even though their sequence structures are similar to cassettes with identifiable genes. Only recently have almost all minicircle classes for some subspecies and isolates of T. brucei been sequenced and annotated with corresponding verification of gRNA expression by small-RNA transcriptome data. These data sets provide a rich resource for understanding the structure of minicircle classes, cassettes and gRNA genes and their transcription. Here, we provide a statistical description of the functionality, expression status, structure and sequence of gRNA genes in a differentiation-competent, laboratory-adapted strain of T. brucei. We obtain a clearer definition of what is a gRNA gene. Our analysis supports the idea that many, if not all, cassettes without an identifiable gRNA gene contain decaying remnants of once functional gRNA genes. Finally, we report several new, unexplained discoveries such as the association between cassette position on the minicircle and gene expression and functionality, and the association between gene initiation sequence and anchor position.