unitas: the universal tool for annotation of small RNAs

piRNAs in the human retina and retinal pigment epithelium reveal a potential role in intracellular trafficking and oxidative stress

Long considered active only in the germline, the PIWI/piRNA pathway is now known to play a significant role in somatic cells, especially neurons. In this study, piRNAs were profiled in the human retina and retinal pigment epithelium (RPE). Furthermore, RNA immunoprecipitation with HIWI2 (PIWIL4) in ARPE19 cells yielded 261 piRNAs, and the expression of selective piRNAs in donor eyes was assessed by qRT-PCR. Intriguingly, computational analysis revealed complete and partial seed sequence similarity between piR-hsa-26131 and the sensory organ specific miR-183/96/182 cluster. Furthermore, the expression of retina-enriched piR-hsa-26131 was positively correlated with miR-182 in HIWI2-silenced Y79 cells. In addition, the lnc-ZNF169 sequence matched with two miRNAs of the let-7 family, and piRNAs, piR-hsa-11361 and piR-hsa-11360, which could modulate the regulatory network of retinal differentiation. Interestingly, we annotated four enriched motifs among the piRNAs and found that the piRNAs containing CACAATG and CTCATCAKYG motifs were snoRNA-derived piRNAs, which are significantly associated with developmental functions. However, piRNAs consisting of ACCACTANACCAC and AKCACGYTCSC motifs were mainly tRNA-derived fragments linked to stress response and sensory perception. Additionally, co-expression network analysis revealed cell cycle control, intracellular transport and stress response as the important biological functions regulated by piRNAs in the retina. Moreover, loss of piRNAs in HIWI2 knockdown ARPE19 confirmed altered expression of targets implicated in intracellular transport, circadian clock, and retinal degeneration. Moreover, piRNAs were dysregulated under oxidative stress conditions, indicating their potential role in retinal pathology. Therefore, we postulate that piRNAs, miRNAs, and lncRNAs might have a functional interplay during retinal development and functions to regulate retinal homeostasis.

In silico approach for the identification of tRNA-derived small non-coding RNAs in SARS-CoV infection

tsRNAs (tRNA-derived small non-coding RNAs), including tRNA halves (tiRNAs) and tRNA fragments (tRFs), have been implicated in some viral infections, such as respiratory viral infections. However, their involvement in SARS-CoV infection is completely unknown. A comprehensive analysis was performed to determine tsRNA populations in a mouse model of SARS-CoV-infected samples containing the wild-type and attenuated viruses. Data from the Gene Expression Omnibus (GEO) dataset at NCBI (accession ID GSE90624 ) was used for this study. A count matrix was generated for the tRNAs. Differentially expressed tRNAs, followed by tsRNAs derived from each significant tRNAs at different conditions and time points between the two groups WT(SARS-CoV-MA15-WT) vs Mock and ΔE (SARS-CoV-MA15-ΔE) vs Mock were identified. Notably, significantly differentially expressed tRNAs at 2dpi but not at 4dpi. The tsRNAs originating from differentially expressed tRNAs across all the samples belonging to each condition (WT, ΔE, and Mock) were identified. Intriguingly, tRFs (tRNA-derived RNA fragments) exhibited higher levels compared to tiRNAs (tRNA-derived stress-induced RNAs) across all samples associated with WT SARS-CoV strain compared to ΔE and mock-infected samples. This discrepancy suggests a non-random formation of tsRNAs, hinting at a possible involvement of tsRNAs in SARS-CoV viral infection.

SARS-CoV-2 remodels the landscape of small non-coding RNAs with infection time and symptom severity

The COVID-19 pandemic caused by the coronavirus SARS-CoV-2 has significantly impacted global health, stressing the necessity of basic understanding of the host response to this viral infection. In this study, we investigated how SARS-CoV-2 remodels the landscape of small non-coding RNAs (sncRNA) from a large collection of nasopharyngeal swab samples taken at various time points from patients with distinct symptom severity. High-throughput RNA sequencing analysis revealed a global alteration of the sncRNA landscape, with abundance peaks related to species of 21-23 and 32-33 nucleotides. Host-derived sncRNAs, including microRNAs (miRNAs), transfer RNA-derived small RNAs (tsRNAs), and small nucleolar RNA-derived small RNAs (sdRNAs) exhibited significant differential expression in infected patients compared to controls. Importantly, miRNA expression was predominantly down-regulated in response to SARS-CoV-2 infection, especially in patients with severe symptoms. Furthermore, we identified specific tsRNAs derived from Glu- and Gly-tRNAs as major altered elements upon infection, with 5' tRNA halves being the most abundant species and suggesting their potential as biomarkers for viral presence and disease severity prediction. Additionally, down-regulation of C/D-box sdRNAs and altered expression of tinyRNAs (tyRNAs) were observed in infected patients. These findings provide valuable insights into the host sncRNA response to SARS-CoV-2 infection and may contribute to the development of further diagnostic and therapeutic strategies in the clinic.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug's responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

sRNA-Effector: A tool to expedite discovery of small RNA regulators

microRNAs (miRNAs) are small regulatory RNAs that repress target mRNA transcripts through base pairing. Although the mechanisms of miRNA production and function are clearly established, new insights into miRNA regulation or miRNA-mediated gene silencing are still emerging. In order to facilitate the discovery of miRNA regulators or effectors, we have developed sRNA-Effector, a machine learning algorithm trained on enhanced crosslinking and immunoprecipitation sequencing and RNA sequencing data following knockdown of specific genes. sRNA-Effector can accurately identify known miRNA biogenesis and effector proteins and identifies 9 putative regulators of miRNA function, including serine/threonine kinase STK33, splicing factor SFPQ, and proto-oncogene BMI1. We validated the role of STK33, SFPQ, and BMI1 in miRNA regulation, showing that sRNA-Effector is useful for identifying new players in small RNA biology. sRNA-Effector will be a web tool available for all researchers to identify potential miRNA regulators in any cell line of interest.

Quantitative and qualitative detection of tRNAs, tRNA halves and tRFs in human cancer samples: Molecular grounds for biomarker development and clinical perspectives

Transfer RNAs (tRNAs) are small non-coding RNAs playing a central role during protein synthesis. Besides translation, growing evidence suggests that in many contexts, precursor or mature tRNAs can also be processed into smaller fragments playing many non-canonical regulatory roles in different biological pathways with oncogenic relevance. Depending on the source, these molecules can be classified as tRNA halves (also known as tiRNAs) or tRNA-derived fragments (tRFs), and furtherly divided into 5'-tRNA and 3'-tRNA halves, or tRF-1, tRF-2, tRF-3, tRF-5, and i-tRF, respectively. Unlike DNA and mRNA, high-throughput sequencing of tRNAs is challenging, because of technical limitations of currently developed sequencing methods. In recent years, different sequencing approaches have been proposed allowing the quantification and identification of an increasing number of tRNA fragments with critical functions in distinct physiological and pathophysiological processes. In the present review, we discussed pros and cons of recent advances in different sequencing methods, also introducing the expanding repertoire of bioinformatics tool and resources specifically focused on tRNA research and discussing current issues in the study of these small RNA molecules. Furthermore, we discussed the potential value of tRNA fragments as diagnostic and prognostic biomarkers for different types of cancers.

DNMT3L interacts with Piwi and modulates the expression of piRNAs in transgenic Drosophila

Aim: To explore the role of Piwi protein and piRNAs in DNMT3L-mediated epigenetic inheritance. Materials & methods: Transgenic Drosophila were used to examine the effect of ectopically expressed DNMT3L on the profile of piRNAs by sequencing of small RNAs. Results & conclusion: Our previous work showed accumulation and inheritance of epimutations across multiple generations in transgenic DNMT3L Drosophila. Here, we show interaction of DNMT3L with Piwi and a significant alteration in the piRNA profile across multiple generations in transgenic Drosophila. In the light of its interaction with histone H1, we propose that in addition to its role of modulating core histone modifications, DNMT3L allows for inheritance of epigenetic information through its collaboration with Piwi, piRNAs and histone H1.

miR-Blood - a small RNA atlas of human blood components

miR-Blood is a high-quality, small RNA expression atlas for the major components of human peripheral blood (plasma, erythrocytes, thrombocytes, monocytes, neutrophils, eosinophils, basophils, natural killer cells, CD4+ T cells, CD8+ T cells, and B cells). Based on the purified blood components from 52 individuals, the dataset provides a comprehensive repository for the expression of 4971 small RNAs from eight non-coding RNA classes.

Depletion of the m1A writer TRMT6/TRMT61A reduces proliferation and resistance against cellular stress in bladder cancer

Background

Bladder cancer (BLCA) is a common and deadly disease that results in a reduced quality of life for the patients and a significant economic burden on society. A better understanding of tumorigenesis is needed to improve clinical outcomes. Recent evidence places the RNA modification m1A and its regulatory proteins TRMT6/TRMT61A and ALKBH3 in BLCA pathogenesis.

Methods

TRMT6/TRMT61A, ALKBH1, and ALKBH3 expression was examined in human BLCA cell lines and a normal urinary tract epithelium cell line through qRT-PCR and western blot analysis. Prestoblue Cell Viability Reagent, wound-healing assay, and live-cell imaging-based cell displacement analysis, were conducted to assess proliferation, migration, and displacement of this BLCA cell line panel. Cell survival was assessed after inducing cellular stress and activating the unfolded protein response (UPR) with tunicamycin. Moreover, siRNA-mediated gene silencing in two BLCA cell lines (5637 and HT1197) was conducted to investigate the biological roles of TRMT6/TRMT61A.

Results

Heterogeneous morphology, proliferation, displacement, tunicamycin sensitivity, and expression levels of m1A regulators were observed among the panel of cell lines examined. In general, TRMT61A expression was increased in BLCA cell lines when compared to SV-HUC-1. Depletion of TRMT6/TRMT61A reduced proliferation capacity in both 5637 and HT1197 cell lines. The average cell displacement of 5637 was also reduced upon TRMT6/TRMT61A depletion. Interestingly, TRMT6/TRMT61A depletion decreased mRNA expression of targets associated with the ATF6-branch of the UPR in 5637 but not in HT1197. Moreover, cell survival after induction of cellular stress was compromised after TRMT6/TRMT61A knockdown in 5637 but not in HT1197 cells.

Conclusion

The findings suggest that TRMT6/TRMT61A plays an oncogenic role in BLCA and is involved in desensitizing BLCA cells against cellular stress. Further investigation into the regulation of TRMT6/TRMT61A expression and its impact on cellular stress tolerance may provide insights for future BLCA treatment.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

Transcriptional and epigenetic changes during tomato yellow leaf curl virus infection in tomato

BACKGROUND

Geminiviruses are DNA plant viruses that cause highly damaging diseases affecting crops worldwide. During the infection, geminiviruses hijack cellular processes, suppress plant defenses, and cause a massive reprogramming of the infected cells leading to major changes in the whole plant homeostasis. The advances in sequencing technologies allow the simultaneous analysis of multiple aspects of viral infection at a large scale, generating new insights into the molecular mechanisms underlying plant-virus interactions. However, an integrative study of the changes in the host transcriptome, small RNA profile and methylome during a geminivirus infection has not been performed yet. Using a time-scale approach, we aim to decipher the gene regulation in tomato in response to the infection with the geminivirus, tomato yellow leaf curl virus (TYLCV).

RESULTS

We showed that tomato undergoes substantial transcriptional and post-transcriptional changes upon TYLCV infection and identified the main altered regulatory pathways. Interestingly, although the principal plant defense-related processes, gene silencing and the immune response were induced, this cannot prevent the establishment of the infection. Moreover, we identified extra- and intracellular immune receptors as targets for the deregulated microRNAs (miRNAs) and established a network for those that also produced phased secondary small interfering RNAs (phasiRNAs). On the other hand, there were no significant genome-wide changes in tomato methylome at 14 days post infection, the time point at which the symptoms were general, and the amount of viral DNA had reached its maximum level, but we were able to identify differentially methylated regions that could be involved in the transcriptional regulation of some of the differentially expressed genes.

CONCLUSION

We have conducted a comprehensive and reliable study on the changes at transcriptional, post-transcriptional and epigenetic levels in tomato throughout TYLCV infection. The generated genomic information is substantial for understanding the genetic, molecular and physiological changes caused by TYLCV infection in tomato.

Human sperm heads harbor modified YsRNA as transgenerationally inherited non-coding RNAs

Most epigenetic information is reprogrammed during gametogenesis and early development. However, some epigenetic information persists and can be inherited, a phenomenon that is common in plants. On the other hand, there are increasing examples of epigenetic inheritance in metazoans, especially for small non-coding RNAs. The presence of regulatory important RNAs in oocytes is undisputed, whereas the corresponding RNA payload in spermatozoa and its regulatory influence in the zygote and early embryogenesis is largely enigmatic. For humans, we herein describe small YRNA fragments (YsRNA) as a paternal contribution to the zygote. First, we trace the biogenesis of these YsRNAs from the source YRNAs with respect to the 5' and 3' modifications. Both the length and modifications make these YsRNAs reminiscent of canonical piRNAs that are not derived from piRNA clusters. Second, from the early stages of spermatogenesis to maturation in the epididymis, we observe distinct YsRNA profile dynamics in the male germline. We detected YsRNAs exclusively in mature sperm heads, the precursor of the male pronucleus in the zygote, suggesting an important role of the epididymis as a site for transmitting and modification of epigenetic information in the form of YsRNA between soma and germline in humans. Since this YsRNA-based epigenetic mechanism is effective across generations, we wondered whether this phenomenon of epigenetic inheritance has an adaptive value. Full-length YRNAs bind to Ro60, an RNA chaperone that additionally binds to non-coding RNAs. We described the profiles of non-coding RNAs bound to Ro60 in the human sperm head and detected specific binding profiles of RNA to Ro60 but no YRNA bound to Ro60. We hypothesize that the sperm head Ro60 system is functional. An adaptive phenotype mediated by the presence of a large amount of YsRNA in the sperm head, and thus as a paternal contribution in the zygote, might be related to an association of YsRNA with YRNA that prevents the adoption of a YRNA secondary structure capable of binding to Ro60. We hypothesize that preventing YRNAs from acting as Ro60-associated gatekeepers for misfolded RNAs in the zygote and early development may enhance RNA chaperoning and, thus, represent the adaptive molecular phenotype.

Long noncoding RNAs emerge from transposon-derived antisense sequences and may contribute to infection stage-specific transposon regulation in a fungal phytopathogen

BACKGROUND

The genome of the obligate biotrophic phytopathogenic barley powdery mildew fungus Blumeria hordei is inflated due to highly abundant and possibly active transposable elements (TEs). In the absence of the otherwise common repeat-induced point mutation transposon defense mechanism, noncoding RNAs could be key for regulating the activity of TEs and coding genes during the pathogenic life cycle.

RESULTS

We performed time-course whole-transcriptome shotgun sequencing (RNA-seq) of total RNA derived from infected barley leaf epidermis at various stages of fungal pathogenesis and observed significant transcript accumulation and time point-dependent regulation of TEs in B. hordei. Using a manually curated consensus database of 344 TEs, we discovered phased small RNAs mapping to 104 consensus transposons, suggesting that RNA interference contributes significantly to their regulation. Further, we identified 5,127 long noncoding RNAs (lncRNAs) genome-wide in B. hordei, of which 823 originated from the antisense strand of a TE. Co-expression network analysis of lncRNAs, TEs, and coding genes throughout the asexual life cycle of B. hordei points at extensive positive and negative co-regulation of lncRNAs, subsets of TEs and coding genes.

CONCLUSIONS

Our work suggests that similar to mammals and plants, fungal lncRNAs support the dynamic modulation of transcript levels, including TEs, during pivotal stages of host infection. The lncRNAs may support transcriptional diversity and plasticity amid loss of coding genes in powdery mildew fungi and may give rise to novel regulatory elements and virulence peptides, thus representing key drivers of rapid evolutionary adaptation to promote pathogenicity and overcome host defense.

PhasiHunter: a robust phased siRNA regulatory cascade mining tool based on multiple reference sequences

SUMMARY

In recent years, phased small interfering RNA has been found to play crucial roles in many biological processes in plants. However, efficiently predicting phasiRNA regulatory cascades with computational methods is still challenging. Here, we introduce PhasiHunter, a phasiRNA regulatory network prediction tool that has several distinctive features compared to existing tools: (i) PhasiHunter employs two major phasiRNA prediction algorithms, namely phase score and hypergeometric distribution-based methods, to ensure the integrity and accuracy of prediction; (ii) PhasiHunter can identify phasiRNAs and their regulatory networks based on multiple reference sequences and the predicted results can be automatically integrated; (iii) PhasiHunter can efficiently identify the phasiRNAs generated through alternative splicing events; and (iv) the excellent data structure and parallel computing architecture allow PhasiHunter to predict phasiRNAs and their regulatory pathways with high efficiency.

AVAILABILITY AND IMPLEMENTATION

PhasiHunter is an open-source tool that is available at https://github.com/HuangLab-CBI/PhasiHunter.

Early Detection of Lung Cancer Using Small RNAs

INTRODUCTION

Lung cancer remains the deadliest cancer in the world, and lung cancer survival is heavily dependent on tumor stage at the time of detection. Low-dose computed tomography screening can reduce mortality; however, annual screening is limited by low adherence in the United States of America and still not broadly implemented in Europe. As a result, less than 10% of lung cancers are detected through existing programs. Thus, there is a great need for additional screening tests, such as a blood test, that could be deployed in the primary care setting.

METHODS

We prospectively recruited 1384 individuals meeting the National Lung Screening Trial demographic eligibility criteria for lung cancer and collected stabilized whole blood to enable the pipetting-free collection of material, thus minimizing preanalytical noise. Ultra-deep small RNA sequencing (20 million reads per sample) was performed with the addition of a method to remove highly abundant erythroid RNAs, and thus open bandwidth for the detection of less abundant species originating from the plasma or the immune cellular compartment. We used 100 random data splits to train and evaluate an ensemble of logistic regression classifiers using small RNA expression of 943 individuals, discovered an 18-small RNA feature consensus signature (miLung), and validated this signature in an independent cohort (441 individuals). Blood cell sorting and tumor tissue sequencing were performed to deconvolve small RNAs into their source of origin.

RESULTS

We generated diagnostic models and report a median receiver-operating characteristic area under the curve of 0.86 (95% confidence interval [CI]: 0.84-0.86) in the discovery cohort and generalized performance of 0.83 in the validation cohort. Diagnostic performance increased in a stage-dependent manner ranging from 0.73 (95% CI: 0.71-0.76) for stage I to 0.90 (95% CI: 0.89-0.90) for stage IV in the discovery cohort and from 0.76 to 0.86 in the validation cohort. We identified a tumor-shed, plasma-bound ribosomal RNA fragment of the L1 stalk as a dominant predictor of lung cancer. The fragment is decreased after surgery with curative intent. In additional experiments, results of dried blood spot collection and sequencing revealed that small RNA analysis could potentially be conducted through home sampling.

CONCLUSIONS

These data suggest the potential of a small RNA-based blood test as a viable alternative to low-dose computed tomography screening for early detection of smoking-associated lung cancer.

Integrated analysis of inflammatory mRNAs, miRNAs, and lncRNAs elucidates the molecular interactome behind bovine mastitis

Mastitis is known as intramammary inflammation, which has a multifactorial complex phenotype. However, the underlying molecular pathogenesis of mastitis remains poorly understood. In this study, we utilized a combination of RNA-seq and miRNA-seq techniques, along with computational systems biology approaches, to gain a deeper understanding of the molecular interactome involved in mastitis. We retrieved and processed one hundred transcriptomic libraries, consisting of 50 RNA-seq and 50 matched miRNA-seq data, obtained from milk-isolated monocytes of Holstein-Friesian cows, both infected with Streptococcus uberis and non-infected controls. Using the weighted gene co-expression network analysis (WGCNA) approach, we constructed co-expressed RNA-seq-based and miRNA-seq-based modules separately. Module-trait relationship analysis was then performed on the RNA-seq-based modules to identify highly-correlated modules associated with clinical traits of mastitis. Functional enrichment analysis was conducted to understand the functional behavior of these modules. Additionally, we assigned the RNA-seq-based modules to the miRNA-seq-based modules and constructed an integrated regulatory network based on the modules of interest. To enhance the reliability of our findings, we conducted further analyses, including hub RNA detection, protein-protein interaction (PPI) network construction, screening of hub-hub RNAs, and target prediction analysis on the detected modules. We identified a total of 17 RNA-seq-based modules and 3 miRNA-seq-based modules. Among the significant highly-correlated RNA-seq-based modules, six modules showed strong associations with clinical characteristics of mastitis. Functional enrichment analysis revealed that the turquoise module was directly related to inflammation persistence and mastitis development. Furthermore, module assignment analysis demonstrated that the blue miRNA-seq-based module post-transcriptionally regulates the turquoise RNA-seq-based module. We also identified a set of different RNAs, including hub-hub genes, hub-hub TFs (transcription factors), hub-hub lncRNAs (long non-coding RNAs), and hub miRNAs within the modules of interest, indicating their central role in the molecular interactome underlying the pathogenic mechanisms of S. uberis infection. This study provides a comprehensive insight into the molecular crosstalk between immunoregulatory mRNAs, miRNAs, and lncRNAs during S. uberis infection. These findings offer valuable directions for the development of molecular diagnosis and biological therapies for mastitis.

Circulating microRNA Profiles for Premature Cardiovascular Death in Patients with Kidney Failure with Replacement Therapy

INTRODUCTION

Patients with kidney failure with replacement therapy (KFRT) suffer from a disproportionately high cardiovascular disease burden. Circulating small non-coding RNAs (c-sncRNAs) have emerged as novel epigenetic regulators and are suggested as novel biomarkers and therapeutic targets for cardiovascular disease; however, little is known about the associations of c-sncRNAs with premature cardiovascular death in KFRT.

METHODS

In a pilot case-control study of 50 hemodialysis patients who died of cardiovascular events as cases, and 50 matched hemodialysis controls who remained alive during a median follow-up of 2.0 years, we performed c-sncRNAs profiles using next-generation sequencing to identify differentially expressed circulating microRNAs (c-miRNAs) between the plasma of cases and that of controls. mRNA target prediction and pathway enrichment analysis were performed to examine the functional relevance of differentially expressed c-miRNAs to cardiovascular pathophysiology. The association of differentially expressed c-miRNAs with cardiovascular mortality was examined using multivariable conditional logistic regression.

RESULTS

The patient characteristics were similar between cases and controls, with a mean age of 63 years, 48% male, and 54% African American in both groups. We detected a total of 613 miRNAs in the plasma, among which five miRNAs (i.e., miR-129-1-5p, miR-500b-3p, miR-125b-1-3p, miR-3648-2-5p, and miR-3150b-3p) were identified to be differentially expressed between cases and controls with cut-offs of p < 0.05 and log2 fold-change (log2FC) > 1. When using more stringent cut-offs of p-adjusted < 0.05 and log2FC > 1, only miR-129-1-5p remained significantly differentially expressed, with higher levels of miR-129-1-5p in the cases than in the controls. The pathway enrichment analysis using predicted miR-129-1-5p mRNA targets demonstrated enrichment in adrenergic signaling in cardiomyocytes, arrhythmogenic right ventricular cardiomyopathy, and oxytocin signaling pathways. In parallel, the circulating miR-129-1-5p levels were significantly associated with the risk of cardiovascular death (adjusted OR [95% CI], 1.68 [1.01-2.81] for one increase in log-transformed miR-129-1-5p counts), independent of potential confounders.

CONCLUSIONS

Circulating miR-129-1-5p may serve as a novel biomarker for premature cardiovascular death in KFRT.

Obesity and PCOS radically alters the snRNA composition of follicular fluid extracellular vesicles

Introduction

The ovarian follicle consists of the oocyte, somatic cells, and follicular fluid (FF). Proper signalling between these compartments is required for optimal folliculogenesis. The association between polycystic ovarian syndrome (PCOS) and extracellular vesicular small non-coding RNAs (snRNAs) signatures in follicular fluid (FF) and how this relates to adiposity is unknown. The purpose of this study was to determine whether FF extracellular vesicle (FFEV)-derived snRNAs are differentially expressed (DE) between PCOS and non-PCOS subjects; and if these differences are vesicle-specific and/or adiposity-dependent.

Methods

FF and granulosa cells (GC) were collected from 35 patients matched by demographic and stimulation parameters. FFEVs were isolated and snRNA libraries were constructed, sequenced, and analyzed.

Results

miRNAs were the most abundant biotype present, with specific enrichment in exosomes (EX), whereas in GCs long non-coding RNAs were the most abundant biotype. In obese PCOS vs. lean PCOS, pathway analysis revealed target genes involved in cell survival and apoptosis, leukocyte differentiation and migration, JAK/STAT, and MAPK signalling. In obese PCOS FFEVs were selectively enriched (FFEVs vs. GCs) for miRNAs targeting p53 signalling, cell survival and apoptosis, FOXO, Hippo, TNF, and MAPK signalling.

Discussion

We provide comprehensive profiling of snRNAs in FFEVs and GCs of PCOS and non-PCOS patients, highlighting the effect of adiposity on these findings. We hypothesize that the selective packaging and release of miRNAs specifically targeting anti-apoptotic genes into the FF may be an attempt by the follicle to reduce the apoptotic pressure of the GCs and stave off premature apoptosis of the follicle observed in PCOS.

Site-specific analysis reveals candidate cross-kingdom small RNAs, tRNA and rRNA fragments, and signs of fungal RNA phasing in the barley-powdery mildew interaction

The establishment of host-microbe interactions requires molecular communication between both partners, which may involve the mutual transfer of noncoding small RNAs. Previous evidence suggests that this is also true for powdery mildew disease in barley, which is caused by the fungal pathogen Blumeria hordei. However, previous studies lacked spatial resolution regarding the accumulation of small RNAs upon host infection by B. hordei. Here, we analysed site-specific small RNA repertoires in the context of the barley-B. hordei interaction. To this end, we dissected infected leaves into separate fractions representing different sites that are key to the pathogenic process: epiphytic fungal mycelium, infected plant epidermis, isolated haustoria, a vesicle-enriched fraction from infected epidermis, and extracellular vesicles. Unexpectedly, we discovered enrichment of specific 31-33-base 5'-terminal fragments of barley 5.8S ribosomal RNA in extracellular vesicles and infected epidermis, as well as particular B. hordei transfer RNA fragments in haustoria. We describe canonical small RNAs from both the plant host and the fungal pathogen that may confer cross-kingdom RNA interference activity. Interestingly, we found first evidence of phased small interfering RNAs in B. hordei, a feature usually attributed to plants, which may be associated with the posttranscriptional control of fungal coding genes, pseudogenes, and transposable elements. Our data suggest a key and possibly site-specific role for cross-kingdom RNA interference and noncoding RNA fragments in the host-pathogen communication between B. hordei and its host barley.

MiSiPi-Rna: an integrated tool for characterizing small regulatory RNA processing

RNA interference (RNAi) is mediated by small (20-30 nucleotide) RNAs that are produced by complex processing pathways. In animals, three main classes are recognized: microRNAs (miRNAs), small-interfering RNAs (siRNAs) and piwi-interacting RNAs (piRNAs). Understanding of small RNA pathways has benefited from genetic models where key enzymatic events were identified that lead to stereotypical positioning of small RNAs relative to precursor transcripts. Increasingly there is interest in using RNAi in non-model systems due to ease of generating synthetic small RNA precursors for research and biotechnology. Unfortunately, small RNAs are often rapidly evolving, requiring investigation of a species' endogenous small RNAs prior to deploying an RNAi approach. This can be accomplished through small non-coding RNA sequencing followed by applying various computational tools; however, the complexity and separately maintained packages lead to significant challenges for annotating global small RNA populations. To address this need, we developed a simple and efficient R package (MiSiPi-Rna) which can be used to characterize pre-selected loci with plots and statistics, aiding researchers understanding RNAi biology specific to their target species. Additionally, MiSiPi-Rna pioneers several computational approaches to identifying Dicer processing to assist annotation of miRNA and siRNA.

Exosome-derived Small RNAs in mouse Sertoli cells inhibit spermatogonial apoptosis

Spermatogenesis is a highly complicated biological process that occurs in the epithelium of the seminiferous tubules. It is regulated by a complex network of endocrine and paracrine factors. Sertoli cells (SCs) play a key role in spermatogenesis due to their production of trophic, differentiation, and immune-modulating factors. However, many of the molecular pathways of SC action remain controversial and unclear. Recently, many studies have focused on exosomes as an important mechanism of intercellular communication. We found that the exosomes derived from mouse SCs inhibited the apoptosis of primary spermatogonia. A total of 1016 miRNAs in SCs and 556 miRNAs in exosomes were detected using miRNA high-throughput sequencing. A total of 294 miRNAs were differentially expressed between SCs and exosomes. Furthermore, 19 tsRNA families appeared in SCs, while 6 tsRNA families appeared in exosomes. A total of 57 and 1 miRNAs (RPM >4) and 14 and 1 tsRNAs were exclusively expressed in SCs and exosomes, respectively. MiR-10b is one of the top ten exosomes with a relatively large enrichment of miRNA. Overexpression of miR-10b downregulates the expression of the target KLF4 to reduce spermatogonial apoptosis in primary spermatogonia or the C18-4 cell line.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

stepRNA: Identification of Dicer cleavage signatures and passenger strand lengths in small RNA sequences

The enzyme Dicer is a component of many small RNA (sRNA) pathways involved in RNA processing for post-transcriptional regulation, anti-viral response and control of transposable elements. Cleavage of double-stranded RNA by Dicer produces a signature overhanging sequence at the 3’ end of the sRNA sequence relative to a complementary passenger strand in a RNA duplex. There is a need for reliable tools to computationally search for Dicer cleavage signatures to help characterise families of sRNAs. This is increasingly important due to the rising popularity of sRNA sequencing, especially in non-model organisms. Here, we present stepRNA, a fast, local tool that identifies (i) overhang signatures strongly indicative of Dicer cleavage in RNA sequences, and (ii) the length of the passenger strand in sRNAs duplexes. We demonstrate the use of stepRNA with simulated and biological datasets to detect Dicer cleavage signatures in experimentally validated examples. Compared to currently available tools, stepRNA is more accurate, requires only sRNA sequence data rather than a reference genome, and provides information about other important features such as passenger strand length. stepRNA is freely available at https://github.com/Vicky-Hunt-Lab/stepRNA and is easily installable.

Investigation of the transability of dietary small non-coding RNAs to animals

Our daily diet not only provides essential nutrients needed for survival and growth but also supplies bioactive ingredients to promote health and prevent disease. Recent studies have shown that exogenous microRNAs (miRNAs), xenomiRs, may enter the consumer’s body through dietary intake and regulate gene expression. This fascinating phenomenon suggests that xenomiRs can act as a new class of bioactive substances associated with mammalian systems. In contrast, several studies have failed to detect xenomiRs in consumers and reported that the observed diet-derived miRNAs in the previous studies can be related to the false positive effects of experiments. This discrepancy can be attributed to the potential artifacts related to the process of experiments, small sample size, and inefficient bioinformatics pipeline. Since this hypothesis is not generally accepted yet, more studies are required. Here, a stringent and reliable bioinformatics pipeline was used to analyze 133 miRNA sequencing data from seven different studies to investigate this phenomenon. Generally, our results do not support the transfer of diet-derived miRNAs into the animal/human tissues in every situation. Briefly, xenomiRs were absent from most samples, and also, their expressions were very low in the samples where they were present, which is unlikely to be sufficient to regulate cell transcripts. Furthermore, this study showed that the possibility of miRNAs being absorbed through animals’ diets and thus influencing gene expression during specific periods of biological development is not inconceivable. In this context, our results were in agreement with the theory of the transfer of small RNAs under certain conditions and periods as xenomiRs were found in colostrum which may modulate infants’ immune systems via post-transcriptional regulation. These findings provide evidence for the selective absorption of diet-derived small RNAs, which need to be investigated in future studies to shed light on the mechanisms underlying the transference of diet-derived miRNAs.

A multiomics disease progression signature of low-risk ccRCC

Clear cell renal cell carcinoma (ccRCC) is the most common renal cancer. Identification of ccRCC likely to progress, despite an apparent low risk at the time of surgery, represents a key clinical issue. From a cohort of adult ccRCC patients (n = 443), we selected low-risk tumors progressing within a 5-years average follow-up (progressors: P, n = 8) and non-progressing (NP) tumors (n = 16). Transcriptome sequencing, miRNA sequencing and proteomics were performed on tissues obtained at surgery. We identified 151 proteins, 1167 mRNAs and 63 miRNAs differentially expressed in P compared to NP low-risk tumors. Pathway analysis demonstrated overrepresentation of proteins related to “LXR/RXR and FXR/RXR Activation”, “Acute Phase Response Signaling” in NP compared to P samples. Integrating mRNA, miRNA and proteomic data, we developed a 10-component classifier including two proteins, three genes and five miRNAs, effectively differentiating P and NP ccRCC and capturing underlying biological differences, potentially useful to identify “low-risk” patients requiring closer surveillance and treatment adjustments. Key results were validated by immunohistochemistry, qPCR and data from publicly available databases. Our work suggests that LXR, FXR and macrophage activation pathways could be critically involved in the inhibition of the progression of low-risk ccRCC. Furthermore, a 10-component classifier could support an early identification of apparently low-risk ccRCC patients.

Characterization of novel small non-coding RNAs and their modifications in bladder cancer using an updated small RNA-seq workflow

Background: Bladder cancer (BLCA) is one of the most common cancer types worldwide. The disease is responsible for about 200,000 deaths annually, thus improved diagnostics and therapy is needed. A large body of evidence reveal that small RNAs of less than 40 nucleotides may act as tumor suppressors, oncogenes, and disease biomarkers, with a major focus on microRNAs. However, the role of other families of small RNAs is not yet deciphered. Recent results suggest that small RNAs and their modification status, play a role in BLCA development and are promising biomarkers due to their high abundance in the exomes and body fluids (including urine). Moreover, free modified nucleosides have been detected at elevated levels from the urine of BLCA patients. A genome-wide view of small RNAs, and their modifications, will help pinpoint the molecules that could be used as biomarker or has important biology in BLCA development. Methods: BLCA tumor tissue specimens were obtained from 12 patients undergoing transurethral resection of non-muscle invasive papillary urothelial carcinomas. Genome-wide profiling of small RNAs less than 40 bases long was performed by a modified protocol with TGIRT (thermostable group II reverse transcriptase) to identify novel small RNAs and their modification status. Results: Comprehensive analysis identified not only microRNAs. Intriguingly, 57 ± 15% (mean ± S.D.) of sequencing reads mapped to non-microRNA-small RNAs including tRNA-derived fragments (tRFs), ribosomal RNA-derived fragments (rRFs) and YRNA-derived fragments (YRFs). Misincorporation (mismatch) sites identified potential base modification positions on the small RNAs, especially on tRFs, corresponding to m1A (N1-methyladenosine), m1G (N1-methylguanosine) and m2 2G (N2, N2-dimethylguanosine). We also detected mismatch sites on rRFs corresponding to known modifications on 28 and 18S rRNA. Conclusion: We found abundant non-microRNA-small RNAs in BLCA tumor samples. Small RNAs, especially tRFs and rRFs, contain modifications that can be captured as mismatch by TGIRT sequencing. Both the modifications and the non-microRNA-small RNAs should be explored as a biomarker for BLCA detection or follow-up.

Dynamics of Small Non-coding RNA Profiles and the Intestinal Microbiome of High and Low Weight Chickens

The host and its symbiotic bacteria form a biological entity, holobiont, in which they share a dynamic connection characterized by symbiosis, co-metabolism, and coevolution. However, how these collaborative relationships were maintained over evolutionary time remains unclear. In this research, the small non-coding RNA (sncRNA) profiles of cecum and their bacteria contents were measured from lines of chickens that have undergone long-term selection for high (HWS) or low (LWS) 56-day body weight. The results from these lines that originated from a common founder population and maintained under the same husbandry showed an association between host intestinal sncRNA expression profile (miRNA, lncRNA fragment, mRNA fragment, snoRNA, and snRNA) and intestinal microbiota. Correlation analyses suggested that some central miRNAs and mRNA fragments had interactions with the abundance of intestinal microbial species and microbiota functions. miR-6622-3p, a significantly differentially expressed (DE) miRNA was correlated with a body weight gain related bacterium, Alistipes putredinis. Our results showed that host sncRNAs may be mediators of interaction between the host and its intestinal microbiome. This provides additional clue for holobiont concepts.

piRNA-like small RNAs target transposable elements in a Clade IV parasitic nematode

The small RNA (sRNA) pathways identified in the model organism Caenorhabditis elegans are not widely conserved across nematodes. For example, the PIWI pathway and PIWI-interacting RNAs (piRNAs) are involved in regulating and silencing transposable elements (TE) in most animals but have been lost in nematodes outside of the C. elegans group (Clade V), and little is known about how nematodes regulate TEs in the absence of the PIWI pathway. Here, we investigated the role of sRNAs in the Clade IV parasitic nematode Strongyloides ratti by comparing two genetically identical adult stages (the parasitic female and free-living female). We identified putative small-interfering RNAs, microRNAs and tRNA-derived sRNA fragments that are differentially expressed between the two adult stages. Two classes of sRNAs were predicted to regulate TE activity including (i) a parasite-associated class of 21-22 nt long sRNAs with a 5' uridine (21-22Us) and a 5' monophosphate, and (ii) 27 nt long sRNAs with a 5' guanine/adenine (27GAs) and a 5' modification. The 21-22Us show striking resemblance to the 21U PIWI-interacting RNAs found in C. elegans, including an AT rich upstream sequence, overlapping loci and physical clustering in the genome. Overall, we have shown that an alternative class of sRNAs compensate for the loss of piRNAs and regulate TE activity in nematodes outside of Clade V.

A transcriptomic atlas underlying developmental plasticity of seasonal forms of Bicyclus anynana butterflies

Organisms residing in regions with alternating seasons often develop different phenotypes, or forms, in each season. These forms are often adaptations to each season and result from an altered developmental response to specific environmental cues such as temperature. While multiple studies have examined form-specific gene expression profiles in a diversity of species, little is known about how environments and developmental transitions, cued by hormone pulses, alter post-transcriptional patterns. In this study, we examine how gene expression, alternative splicing, and miRNA-mediated gene silencing in Bicyclus anynana butterfly hindwing tissue, varies across two rearing temperatures at four developmental timepoints. These timepoints flank two temperature-sensitive periods that coincide with two pulses of the insect hormone 20E. Our results suggest that developmental transitions, coincident with 20E pulses, elicit a greater impact on all these transcriptomic patterns than rearing temperatures per se. More similar transcriptomic patterns are observed pre-20E pulses than those observed post-20E pulses. We also found functionally distinct sets of differentially expressed and differentially spliced genes in the seasonal forms. Furthermore, around 10% of differentially expressed genes are predicted to be direct targets of, and regulated by, differentially expressed miRNAs between the seasonal forms. Many differentially expressed genes, miRNAs, or differentially spliced genes potentially regulate eyespot size plasticity, and we validated the differential splicing pattern of one such gene, daughterless. We present a comprehensive and interactive transcriptomic atlas of the hindwing tissue of both seasonal forms of B. anynana throughout development, a model organism of seasonal plasticity.

Characteristics of miRNAs Present in Bovine Sperm and Associations With Differences in Fertility

Small non-coding RNAs have been linked to different phenotypes in bovine sperm, however attempts to identify sperm-borne molecular biomarkers of male fertility have thus far failed to identify a robust profile of expressed miRNAs related to fertility. We hypothesized that some differences in bull fertility may be reflected in the levels of different miRNAs in sperm. To explore such differences in fertility that are not due to differences in visible metrics of sperm quality, we employed Next Generation Sequencing to compare the miRNA populations in Bos taurus sperm from bulls with comparable motility and morphology but varying Sire Conception Rates. We identified the most abundant miRNAs in both populations (miRs -34b-3p; -100-5p; -191-5p; -30d-4p; -21-5p) and evaluated differences in the overall levels and specific patterns of isomiR expression. We also explored correlations between specific pairs of miRNAs in each population and identified 10 distinct pairs of miRNAs that were positively correlated in bulls with higher fertility and negatively correlated in comparatively less fertile individuals. Furthermore, 8 additional miRNA pairs demonstrated the opposite trend; negatively correlated in high fertility animals and positively correlated in less fertile bulls. Finally, we performed pathway analysis to identify potential roles of miRNAs present in bull sperm in the regulation of specific genes that impact spermatogenesis and embryo development. Together, these results present a comprehensive picture of the bovine sperm miRNAome that suggests multiple potential roles in fertility.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

Exosome-derived small non-coding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca)

Exosomes, a subset of small extracellular vesicles, carry various nucleic acids, proteins, lipids, amino acids and metabolites. They function as a mode of intercellular communication and molecular transfer. Exosome cargo molecules, including small non-coding RNAs (sncRNAs), are involved in the immune response in various organisms. However, the role of exosome-derived sncRNAs in immune responses in molluscs remains unclear. Here, we aimed to reveal the sncRNAs involved in the immune response during grafting transplantation by the pearl oyster Pinctada fucata. Exosomes were successfully extracted from the P. fucata haemolymph during graft transplantation. Abundant microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) were simultaneously discovered in P. fucata exosomes by small RNA sequencing. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to immune response upon grafting transplantation, whereas piRNAs may also be associated with transposon silencing by targeting with genome transposon elements. This work provides the basis for a functional understanding of exosome-derived sncRNAs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in molluscs.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

The Multiverse of Plant Small RNAs: How Can We Explore It?

Plant small RNAs (sRNAs) are a heterogeneous group of noncoding RNAs with a length of 20–24 nucleotides that are widely studied due to their importance as major regulators in various biological processes. sRNAs are divided into two main classes—microRNAs (miRNAs) and small interfering RNAs (siRNAs)—which differ in their biogenesis and functional pathways. Their identification and enrichment with new structural variants would not be possible without the use of various high-throughput sequencing (NGS) techniques, allowing for the detection of the total population of sRNAs in plants. Classifying sRNAs and predicting their functional role based on such high-performance datasets is a nontrivial bioinformatics task, as plants can generate millions of sRNAs from a variety of biosynthetic pathways. Over the years, many computing tools have been developed to meet this challenge. Here, we review more than 35 tools developed specifically for plant sRNAs over the past few years and explore some of their basic algorithms for performing tasks related to predicting, identifying, categorizing, and quantifying individual sRNAs in plant samples, as well as visualizing the results of these analyzes. We believe that this review will be practical for biologists who want to analyze their plant sRNA datasets but are overwhelmed by the number of tools available, thus answering the basic question of how to choose the right one for a particular study.

Activation of DNA Transposons and Evolution of piRNA Genes Through Interspecific Hybridization in Xenopus Frogs

Interspecific hybridization between two closely related species sometimes resulted in a new species with allotetraploid genomes. Many clawed frog species belonging to the Xenopus genus have diverged from the allotetraploid ancestor created by the hybridization of two closely related species with the predicted L and S genomes. There are species-specific repeated sequences including transposable elements in each genome of organisms that reproduce sexually. To understand what happened on and after the hybridization of the two distinct systems consisting of repeated sequences and their corresponding piRNAs, we isolated small RNAs from ovaries and testes of three Xenopus species consisting of allotetraploid X. laevis and X. borealis and diploid X. tropicalis as controls. After a comprehensive sequencing and selection of piRNAs, comparison of their sequences showed that most piRNA sequences were different between the ovaries and testes in all three species. We compared piRNA and genome sequences and specified gene clusters for piRNA expression in each genome. The synteny and homology analyses showed many distinct piRNA clusters among the three species and even between the two L and/or S subgenomes, indicating that most clusters of the two allotetraploid species changed after hybridization. Moreover, evolutionary analysis showed that DNA transposons including Kolobok superfamily might get activated just after hybridization and then gradually inactivated. These findings suggest that some DNA transposons and their piRNAs might greatly influence allotetraploid genome evolution after hybridization.

piRNAs and PIWI Proteins as Diagnostic and Prognostic Markers of Genitourinary Cancers

piRNAs (PIWI-interacting RNAs) are small non-coding RNAs capable of regulation of transposon and gene expression. piRNAs utilise multiple mechanisms to affect gene expression, which makes them potentially more powerful regulators than microRNAs. The mechanisms by which piRNAs regulate transposon and gene expression include DNA methylation, histone modifications, and mRNA degradation. Genitourinary cancers (GC) are a large group of neoplasms that differ by their incidence, clinical course, biology, and prognosis for patients. Regardless of the GC type, metastatic disease remains a key therapeutic challenge, largely affecting patients’ survival rates. Recent studies indicate that piRNAs could serve as potentially useful biomarkers allowing for early cancer detection and therapeutic interventions at the stage of non-advanced tumour, improving patient’s outcomes. Furthermore, studies in prostate cancer show that piRNAs contribute to cancer progression by affecting key oncogenic pathways such as PI3K/AKT. Here, we discuss recent findings on biogenesis, mechanisms of action and the role of piRNAs and the associated PIWI proteins in GC. We also present tools that may be useful for studies on the functioning of piRNAs in cancers.

Methodologies for Discovery and Quantitative Profiling of sRNAs in Potato

Small RNAs (sRNAs) are short noncoding RNAs involved in the regulation of a wide range of biological processes in plants. Advances in high-throughput sequencing and development of new computational tools had facilitated the discovery of different classes of sRNAs, their quantification, and elucidation of their functional role in gene expression regulation by target transcript predictions. The workflow presented here allows identification of different sRNA species: known and novel potato miRNAs, and their sequence variants (isomiRs), as well as identification of phased small interfering RNAs (phasiRNAs). Moreover, it includes steps for differential expression analysis to search for regulated sRNAs across different tested biological conditions. In addition, it describes two different methods for predicting sRNA targets, in silico prediction, and degradome sequencing data analysis. All steps of the workflow are written in a clear and user-friendly way; thus they can be followed also by the users with minimal bioinformatics knowledge. We also included several in-house scripts together with valuable notes to facilitate data (pre)processing steps and to reduce the analysis time.

Distinct small non-coding RNA landscape in the axons and released extracellular vesicles of developing primary cortical neurons and the axoplasm of adult nerves

Neurons have highlighted the needs for decentralized gene expression and specific RNA function in somato-dendritic and axonal compartments, as well as in intercellular communication via extracellular vesicles (EVs). Despite advances in miRNA biology, the identity and regulatory capacity of other small non-coding RNAs (sncRNAs) in neuronal models and local subdomains has been largely unexplored.We identified a highly complex and differentially localized content of sncRNAs in axons and EVs during early neuronal development of cortical primary neurons and in adult axons in vivo. This content goes far beyond miRNAs and includes most known sncRNAs and precisely processed fragments from tRNAs, sno/snRNAs, Y RNAs and vtRNAs. Although miRNAs are the major sncRNA biotype in whole-cell samples, their relative abundance is significantly decreased in axons and neuronal EVs, where specific tRNA fragments (tRFs and tRHs/tiRNAs) mainly derived from tRNAs Gly-GCC, Val-CAC and Val-AAC predominate. Notably, although 5'-tRHs compose the great majority of tRNA-derived fragments observed in vitro, a shift to 3'-tRNAs is observed in mature axons in vivo.The existence of these complex sncRNA populations that are specific to distinct neuronal subdomains and selectively incorporated into EVs, equip neurons with key molecular tools for spatiotemporal functional control and cell-to-cell communication.

Large Drosophila germline piRNA clusters are evolutionarily labile and dispensable for transposon regulation

PIWI proteins and their guiding Piwi-interacting small RNAs (piRNAs) are crucial for fertility and transposon defense in the animal germline. In most species, the majority of piRNAs are produced from distinct large genomic loci, called piRNA clusters. It is assumed that germline-expressed piRNA clusters, particularly in Drosophila, act as principal regulators to control transposons dispersed across the genome. Here, using synteny analysis, we show that large clusters are evolutionarily labile, arise at loci characterized by recurrent chromosomal rearrangements, and are mostly species-specific across the Drosophila genus. By engineering chromosomal deletions in D. melanogaster, we demonstrate that the three largest germline clusters, which account for the accumulation of >40% of all transposon-targeting piRNAs in ovaries, are neither required for fertility nor for transposon regulation in trans. We provide further evidence that dispersed elements, rather than the regulatory action of large Drosophila germline clusters in trans, may be central for transposon defense.

Bioinformatics and Machine Learning Approaches to Understand the Regulation of Mobile Genetic Elements

Simple Summary

Transposable elements (TEs) are DNA sequences that are, or were, able to move (transpose) within the genome of a single cell. They were first discovered by Barbara McClintock while working on maize, and they make up a large fraction of the genome. Transpositions can result in mutations and they can alter the genome size. Cells regulate the activity of TEs using a variety of mechanisms, such as chemical modifications of DNA and small RNAs. Machine learning (ML) is an interdisciplinary subject that studies computer algorithms that can improve through experience and by the use of data. ML has been successfully applied to a variety of problems in bioinformatics and has exhibited favorable precision and speed. Here, we provide a systematic and guided review on the ML and bioinformatic methods and tools that are used for the analysis of the regulation of TEs.

Abstract

Transposable elements (TEs, or mobile genetic elements, MGEs) are ubiquitous genetic elements that make up a substantial proportion of the genome of many species. The recent growing interest in understanding the evolution and function of TEs has revealed that TEs play a dual role in genome evolution, development, disease, and drug resistance. Cells regulate TE expression against uncontrolled activity that can lead to developmental defects and disease, using multiple strategies, such as DNA chemical modification, small RNA (sRNA) silencing, chromatin modification, as well as sequence-specific repressors. Advancements in bioinformatics and machine learning approaches are increasingly contributing to the analysis of the regulation mechanisms. A plethora of tools and machine learning approaches have been developed for prediction, annotation, and expression profiling of sRNAs, for methylation analysis of TEs, as well as for genome-wide methylation analysis through bisulfite sequencing data. In this review, we provide a guided overview of the bioinformatic and machine learning state of the art of fields closely associated with TE regulation and function.

An integrated analysis of mRNAs, lncRNAs, and miRNAs based on weighted gene co-expression network analysis involved in bovine endometritis

In dairy cattle, endometritis is a severe infectious disease that occurs following parturition. It is clear that genetic factors are involved in the etiology of endometritis, however, the molecular pathogenesis of endometritis is not entirely understood. In this study, a system biology approach was used to better understand the molecular mechanisms underlying the development of endometritis. Forty transcriptomic datasets comprising of 20 RNA-Seq (GSE66825) and 20 miRNA-Seq (GSE66826) were obtained from the GEO database. Next, the co-expressed modules were constructed based on RNA-Seq (Rb-modules) and miRNA-Seq (mb-modules) data, separately, using a weighted gene co-expression network analysis (WGCNA) approach. Preservation analysis was used to find the non-preserved Rb-modules in endometritis samples. Afterward, the non-preserved Rb-modules were assigned to the mb-modules to construct the integrated regulatory networks. Just highly connected genes (hubs) in the networks were considered and functional enrichment analysis was used to identify the biological pathways associated with the development of the disease. Furthermore, additional bioinformatic analysis including protein-protein interactions network and miRNA target prediction were applied to enhance the reliability of the results. Thirty-five Rb-modules and 10 mb-modules were identified and 19 and 10 modules were non-preserved, respectively, which were enriched in biological pathways related to endometritis like inflammation and ciliogenesis. Two non-preserved Rb-modules were significantly assigned to three mb-modules and three and two important sub-networks in the Rb-modules were identified, respectively, including important mRNAs, lncRNAs and miRNAs genes like IRAK1, CASP3, CCDC40, CCDC39, ZMYND10, FOXJ1, TLR4, IL10, STAT3, FN1, AKT1, CD68, ENSBTAG00000049936, ENSBTAG00000050527, ENSBTAG00000051242, ENSBTAG00000049287, bta-miR-449, bta-miR-484, bta-miR-149, bta-miR-30b and bta-miR-423. The potential roles of these genes have been previously demonstrated in endometritis or related pathways, which reinforced putative functions of the suggested integrated regulatory networks in the endometritis pathogenesis. These findings may help further elucidate the underlying mechanisms of bovine endometritis.

Genome-wide transcriptomic analysis of highly virulent African swine fever virus infection reveals complex and unique virus host interaction

African swine fever virus (ASFV), one of the most devastating emerging swine pathogens in China, causes nearly 100 % mortality in naive herds. Here, whole-transcriptome RNA-seq analysis was conducted in porcine alveolar macrophages (PAMs) infected with Pig/Heilongjiang/2018 (Pig/HLJ/18) ASFV at different time points. Our data suggested that ASFV genes expression demonstrated a time-depended pattern and ASFV early genes were involved in antagonizing host innate immunity. Moreover, viral small RNA (vsRNA) was generated as well. Meanwhile, transcriptome analysis of host genes suggested a strong inhibition host immunity-related genes by ASFV infection in PAMs, while enhanced chemokine-mediated signaling pathways and neutrophil chemotaxis were observed in ASFV infected PAMs. Furthermore, ASFV infection also down-regulated host microRNAs (miRNAs) that putatively targeted viral genes, while also triggering dysregulation of host metabolism that promoted virus replication at transcription level. Most importantly, infection of PAMs with ASFV induced a different transcriptome pattern from that of highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV), which is known to trigger a host cytokine storm. In conclusion, our transcriptome data implied that ASFV infection in PAMs appeared to be associated with strong inhibition of host immune responses, dysregulation of host chemokine axis and metabolic pathways.

Specific PIWI-Interacting RNAs and Related Small Noncoding RNAs Are Associated With Ovarian Aging in Ames Dwarf (df/df) Mice

The Ames dwarf (df/df) mouse is a well-established model for delayed aging. MicroRNAs (miRNAs), the most studied small noncoding RNAs (sncRNAs), may regulate ovarian aging to maintain a younger ovarian phenotype in df/df mice. In this study, we profile other types of ovarian sncRNAs, PIWI-interacting RNAs (piRNAs) and piRNA-like RNAs (piLRNAs), in young and aged df/df and normal mice. Half of the piRNAs derive from transfer RNA fragments (tRF-piRNAs). Aging and dwarfism alter the ovarian expression of these novel sncRNAs. Specific tRF-piRNAs that increased with age might target and decrease the expression of the breast cancer antiestrogen resistance protein 3 (BCAR3) gene in the ovaries of old df/df mice. A set of piLRNAs that decreased with age and map to D10Wsu102e mRNA may have trans-regulatory functions. Other piLRNAs that decreased with age potentially target and may de-repress transposable elements, leading to a beneficial impact on ovarian aging in df/df mice. These results identify unique responses in ovarian tissues with regard to aging and dwarfism. Overall, our findings highlight the complexity of the aging effects on gene expression and suggest that, in addition to miRNAs, piRNAs, piLRNAs, tRF-piRNAs, and their potential targets can be central players in the maintenance of a younger ovarian phenotype in df/df mice.

piRNAclusterDB 2.0: update and expansion of the piRNA cluster database

PIWI-interacting RNAs (piRNAs) and their partnering PIWI proteins defend the animal germline against transposable elements and play a crucial role in fertility. Numerous studies in the past have uncovered many additional functions of the piRNA pathway, including gene regulation, anti-viral defense, and somatic transposon repression. Further, comparative analyses across phylogenetic groups showed that the PIWI/piRNA system evolves rapidly and exhibits great evolutionary plasticity. However, the presence of so-called piRNA clusters as the major source of piRNAs is common to nearly all metazoan species. These genomic piRNA-producing loci are highly divergent across taxa and critically influence piRNA populations in different evolutionary lineages. We launched the initial version of the piRNA cluster database to facilitate research on regulation and evolution of piRNA-producing loci across tissues und species. In recent years the amount of small RNA sequencing data that was generated and the abundance of species that were studied has grown rapidly. To keep up with this recent progress, we have released a major update for the piRNA cluster database (https://www.smallrnagroup.uni-mainz.de/piRNAclusterDB), expanding it from 12 to a total of 51 species with hundreds of new datasets, and revised its overall structure to enable easy navigation through this large amount of data.

tsRNAsearch: A pipeline for the identification of tRNA and ncRNA fragments from small RNA-sequencing data

MOTIVATION

tRNAs were originally considered uni-functional RNA molecules involved in the delivery of amino acids to growing peptide chains on the ribosome. More recently, the liberation of tRNA fragments from tRNAs via specific enzyme cleavage has been characterized. Detection of tRNA fragments in sequencing data is difficult due to tRNA sequence redundancy and the short length of both tRNAs and their fragments.

RESULTS

Here we introduce tsRNAsearch, a Nextflow pipeline for the identification of differentially abundant tRNA fragments and other non-coding RNAs from small RNA-sequencing data. tsRNAsearch is intended for use when comparing two groups of datasets, such as control and treatment groups. tsRNAsearch comparatively searches for tRNAs and ncRNAs with irregular read distribution profiles (a proxy for RNA cleavage) using a combined score made up of four novel methods and a differential expression analysis, and reports the top ranked results in simple PDF and TEXT files. In this study, we used publicly available small RNA-seq data to replicate the identification of tsRNAs from chronic hepatitis-infected liver tissue data. In addition, we applied tsRNAsearch to pancreatic ductal adenocarcinoma (PDAC) and matched healthy pancreatic tissue small RNA-sequencing data. Our results support the identification of miR135b from the original study as a potential biomarker of PDAC and identify other potentially stronger miRNA biomarkers of PDAC.

AVAILABILITY

https://github.com/GiantSpaceRobot/tsRNAsearch.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrative Analyses of Circulating Small RNAs and Kidney Graft Transcriptome in Transplant Glomerulopathy

Transplant glomerulopathy develops through multiple mechanisms, including donor-specific antibodies, T cells and innate immunity. This study investigates circulating small RNA profiles in serum samples of kidney transplant recipients with biopsy-proven transplant glomerulopathy. Among total small RNA population, miRNAs were the most abundant species in the serum of kidney transplant patients. In addition, fragments arising from mature tRNA and rRNA were detected. Most of the tRNA fragments were generated from 5′ ends of mature tRNA and mainly from two parental tRNAs: tRNA-Gly and tRNA-Glu. Moreover, transplant patients with transplant glomerulopathy displayed a novel tRNA fragments signature. Gene expression analysis from allograft tissues demonstrated changes in canonical pathways related to immune activation such as iCos-iCosL signaling pathway in T helper cells, Th1 and Th2 activation pathway, and dendritic cell maturation. mRNA targets of down-regulated miRNAs such as miR-1224-5p, miR-4508, miR-320, miR-378a from serum were globally upregulated in tissue. Integration of serum miRNA profiles with tissue gene expression showed that changes in serum miRNAs support the role of T-cell mediated mechanisms in ongoing allograft injury.

Differentially Expressed Extracellular Vesicle-Contained microRNAs before and after Transurethral Resection of Bladder Tumors

Bladder cancer (BC) is currently diagnosed and monitored by cystoscopy, a costly and invasive procedure. Potential biomarkers in urine, blood, and, more recently, extracellular vesicles (EVs), have been explored as non-invasive alternatives for diagnosis and surveillance of BC. EVs are nanovesicles secreted by most cell types containing diverse molecular cargo, including different types of small RNAs, such as microRNA (miRNA). In this study, we performed next-generation sequencing of EV-contained miRNA isolated from urine and serum of 41 patients with non-muscle invasive BC (27 stage Ta, 14 stage T1) and 15 non-cancer patients (NCP) with benign cystoscopy findings. MiRNA sequencing was also performed on serum supernatant samples for T1 patients. To identify potential BC-specific biomarkers, expression levels of miRNA in presurgery samples were compared to those at postsurgery check-ups, and to NCPs. Results showed that two miRNAs, urinary EV-contained miR-451a and miR-486-5p, were significantly upregulated in presurgery samples from T1 patients compared to postsurgery check-up samples. This was confirmed in a replica EV/RNA isolation and sequencing run of 10 T1 patients from the primary run; however, analyses revealed no differential expression of miRNAs in serum EVs, serum supernatant, or when comparing BC patients to NCPs. This is the first study to investigate EV-containing miRNA sequencing in pre- and postsurgery BC patient samples and our findings suggest that urinary EV-contained miR-451a and miR-486-5p may be potential biomarkers for recurrence-free survival of BC patients with stage T1 disease.

Dynamics of small non-coding RNAs in bovine scNT embryos through the maternal-to-embryonic transition

The efficiency of somatic cell nuclear transfer (scNT) for production of viable offspring is relatively low as compared to in vitro fertilization (IVF), presumably due to deficiencies in epigenetic reprogramming of the donor cell genome. Such defects may also involve the population of small non-coding RNAs (sncRNAs), which are important during early embryonic development. The objective of this study was to examine dynamic changes in relative abundance of sncRNAs during the maternal-to embryonic transition (MET) in bovine embryos produced by scNT as compared to IVF by using RNA sequencing. When comparing populations of miRNA in scNT versus IVF embryos, only miR-2340, miR-345, and miR34a were differentially expressed in morulae, though many more miRNAs were differentially expressed when comparing across developmental stages. Also of interest, distinct populations of piwi-interacting like RNAs (pilRNAs) were identified in bovine embryos prior to and during embryonic genome activation (EGA) as compared bovine embryos post EGA and differentiated cells. Overall, sncRNA sequencing analysis of preimplantation embryos revealed largely similar profiles of sncRNAs for IVF and scNT embryos at the 2-cell, 8-cell, morula and blastocyst stages of development. However, these sncRNA profiles, including miRNA, piRNA and tRNA fragments, were notably distinct prior to and after completion of the MET.

Differential Response of Grapevine to Infection with 'Candidatus Phytoplasma solani' in Early and Late Growing Season through Complex Regulation of mRNA and Small RNA Transcriptomes

Bois noir is the most widespread phytoplasma grapevine disease in Europe. It is associated with ‘Candidatus Phytoplasma solani’, but molecular interactions between the causal pathogen and its host plant are not well understood. In this work, we combined the analysis of high-throughput RNA-Seq and sRNA-Seq data with interaction network analysis for finding new cross-talks among pathways involved in infection of grapevine cv. Zweigelt with ‘Ca. P. solani’ in early and late growing seasons. While the early growing season was very dynamic at the transcriptional level in asymptomatic grapevines, the regulation at the level of small RNAs was more pronounced later in the season when symptoms developed in infected grapevines. Most differentially expressed small RNAs were associated with biotic stress. Our study also exposes the less-studied role of hormones in disease development and shows that hormonal balance was already perturbed before symptoms development in infected grapevines. Analysis at the level of communities of genes and mRNA-microRNA interaction networks revealed several new genes (e.g., expansins and cryptdin) that have not been associated with phytoplasma pathogenicity previously. These novel actors may present a new reference framework for research and diagnostics of phytoplasma diseases of grapevine.

Altered non-coding RNA expression profile in F1 progeny 1 year after parental irradiation is linked to adverse effects in zebrafish

Gamma radiation produces DNA instability and impaired phenotype. Previously, we observed negative effects on phenotype, DNA methylation, and gene expression profiles, in offspring of zebrafish exposed to gamma radiation during gametogenesis. We hypothesize that previously observed effects are accompanied with changes in the expression profile of non-coding RNAs, inherited by next generations. Non-coding RNA expression profile was analysed in F₁ offspring (5.5 h post-fertilization) by high-throughput sequencing 1 year after parental irradiation (8.7 mGy/h, 5.2 Gy total dose). Using our previous F₁-γ genome-wide gene expression data (GSE98539), hundreds of mRNAs were predicted as targets of differentially expressed (DE) miRNAs, involved in pathways such as insulin receptor, NFkB and PTEN signalling, linking to apoptosis and cancer. snRNAs belonging to the five major spliceosomal snRNAs were down-regulated in the F₁-γ group, Indicating transcriptional and post-transcriptional alterations. In addition, DEpiRNA clusters were associated to 9 transposable elements (TEs) (LTR, LINE, and TIR) (p = 0.0024), probable as a response to the activation of these TEs. Moreover, the expression of the lincRNAs malat-1, and several others was altered in the offspring F₁, in concordance with previously observed phenotypical alterations. In conclusion, our results demonstrate diverse gamma radiation-induced alterations in the ncRNA profiles of F₁ offspring observable 1 year after parental irradiation.