Assessing small RNA profiles in potato diploid hybrid and its resynthesized allopolyploid reveals conserved abundance with distinct genomic distribution.

KEY MESSAGE: The shock produced by the allopolyploidization process on a potato interspecific diploid hybrid displays a non-random remobilization of the small RNAs profile on a variety of genomic features. Allopolyploidy, a complex process involving interspecific hybridization and whole genome duplication, significantly impacts plant evolution, leading to the emergence of novel phenotypes. Polyploids often present phenotypic nuances that enhance adaptability, enabling them to compete better and occasionally to colonize new habitats. Whole-genome duplication represents a genomic "shock" that can trigger genetic and epigenetic changes that yield novel expression patterns. In this work, we investigate the polyploidization effect on a diploid interspecific hybrid obtained through the cross between the cultivated potato Solanum tuberosum and the wild potato Solanum kurtzianum, by assessing the small RNAs (sRNAs) profile of the parental diploid hybrid and its derived allopolyploid. Small RNAs are key components of the epigenetic mechanisms involved in silencing by RNA-directed DNA Methylation (RdDM). A sRNA sequencing (sRNA-Seq) analysis was performed to individually profile the 21 to 22 nucleotide (21 to 22-nt) and 24-nt sRNA size classes due to their unique mechanism of biogenesis and mode of function. The composition and distribution of different genomic features and differentially accumulated (DA) sRNAs were evaluated throughout the potato genome. We selected a subset of genes associated with DA sRNAs for messenger RNA (mRNA) expression analysis to assess potential impacts on the transcriptome. Interestingly, we noted that 24-nt DA sRNAs that exclusively mapped to exons were correlated with differentially expressed mRNAs between genotypes, while this behavior was not observed when 24-nt DA sRNAs were mapped to intronic regions. These findings collectively emphasize the nonstochastic nature of sRNA remobilization in response to the genomic shock induced by allopolyploidization.

The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming.

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

Dynamics and Function of sRNA/mRNAs Under the Scrutiny of Computational Simulation Methods.

Molecular dynamics simulations have proved extremely useful in investigating the functioning of proteins with atomic-scale resolution. Many applications to the study of RNA also exist, and their number increases by the day. However, implementing MD simulations for RNA molecules in solution faces challenges that the MD practitioner must be aware of for the appropriate use of this tool. In this chapter, we present the fundamentals of MD simulations, in general, and the peculiarities of RNA simulations, in particular. We discuss the strengths and limitations of the technique and provide examples of its application to elucidate small RNA's performance.

Differential expression of miRNAs associated with pectoral myopathies in young broilers: insights from a comparative transcriptome analysis.

INTRODUCTION: White Striping (WS) and Wooden Breast (WB) pectoral myopathies are relevant disorders for contemporary broiler production worldwide. Several studies aimed to elucidate the genetic components associated with the occurrence of these myopathies. However, epigenetic factors that trigger or differentiate these two conditions are still unclear. The aim of this study was to identify miRNAs differentially expressed (DE) between normal and WS and WB-affected broilers, and to verify the possible role of these miRNAs in metabolic pathways related to the manifestation of these pectoral myopathies in 28-day-old broilers. RESULTS: Five miRNAs were DE in the WS vs control (gga-miR-375, gga-miR-200b-3p, gga-miR-429-3p, gga-miR-1769-5p, gga-miR-200a-3p), 82 between WB vs control and 62 between WB vs WS. Several known miRNAs were associated with WB, such as gga-miR-155, gga-miR-146b, gga-miR-222, gga-miR-146-5p, gga-miR- 29, gga-miR-21-5p, gga-miR-133a-3p and gga-miR-133b. Most of them had not previously been associated with the development of this myopathy in broilers. We also have predicted 17 new miRNAs expressed in the broilers pectoral muscle. DE miRNA target gene ontology analysis enriched 6 common pathways for WS and WB compared to control: autophagy, insulin signaling, FoxO signaling, endocytosis, and metabolic pathways. The WS vs control contrast had two unique pathways, ERBB signaling and the mTOR signaling, while WB vs control had 14 unique pathways, with ubiquitin-mediated proteolysis and endoplasmic reticulum protein processing being the most significant. CONCLUSIONS: We found miRNAs DE between normal broilers and those affected with breast myopathies at 28 days of age. Our results also provide novel evidence of the miRNAs role on the regulation of WS and in the differentiation of both WS and WB myopathies. Overall, our study provides insights into miRNA-mediated and pathways involved in the occurrence of WS and WB helping to better understand these chicken growth disorders in an early age. These findings can help developing new approaches to reduce these complex issues in poultry production possibly by adjustments in nutrition and management conditions. Moreover, the miRNAs and target genes associated with the initial stages of WS and WB development could be potential biomarkers to be used in selection to reduce the occurrence of these myopathies in broiler production.

Identification of novel small RNAs in extracellular vesicles produced by Actinobacillus pleuropneumoniae.

Extracellular vesicle (EV) production by bacteria is an important mechanism for microbial communication and host-pathogen interaction. EVs of some bacterial species have been reported to contain nucleic acids. However, the role of small RNAs (sRNAs) packaged in EVs is poorly understood. Here, we report on the RNA cargo of EVs produced by the pig pathogen Actinobacillus pleuropneumoniae, the causal agent of porcine pleuropneumonia, a disease which causes substantial economic losses to the swine industry worldwide. The EVs produced by aerobically and anaerobically grown bacteria were only slightly different in size and distribution. Total cell and outer membrane protein profiles and lipid composition of A. pleuropneumoniae whole cell extracts and EVs were similar, although EVs contained rough lipopolysaccharide compared to the smooth form in whole cells. Approximately 50% of Galleria mellonella larvae died after the injection of EVs. RNAseq, RT-PCR, protection from nuclease degradation, and database searching identified previously described and 13 novel A. pleuropneumoniae sRNAs in EVs, some of which were enriched compared to whole cell content. We conclude that A. pleuropneumoniae EVs contain sRNAs, including those known to be involved in virulence, and some with homologs in other Pasteurellaceae and/or non-Pasteurellaceae. Further work will establish whether the novel sRNAs in A. pleuropneumoniae EVs play any role in pathogenesis.

Computational Prediction of RNA-RNA Interactions between Small RNA Tracks from Betacoronavirus Nonstructural Protein 3 and Neurotrophin Genes during Infection of an Epithelial Lung Cancer Cell Line: Potential Role of Novel Small Regulatory RNA.

Whether RNA-RNA interactions of cytoplasmic RNA viruses, such as Betacoronavirus, might end in the biogenesis of putative virus-derived small RNAs as miRNA-like molecules has been controversial. Even more, whether RNA-RNA interactions of wild animal viruses may act as virus-derived small RNAs is unknown. Here, we address these issues in four ways. First, we use conserved RNA structures undergoing negative selection in the genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 circulating in different bat species, intermediate animals, and human hosts. Second, a systematic literature review was conducted to identify Betacoronavirus-targeting hsa-miRNAs involved in lung cell infection. Third, we employed sophisticated long-range RNA-RNA interactions to refine the seed sequence homology of hsa-miRNAs with conserved RNA structures. Fourth, we used high-throughput RNA sequencing of a Betacoronavirus-infected epithelial lung cancer cell line (Calu-3) to validate the results. We proposed nine potential virus-derived small RNAs: two vsRNAs in SARS-CoV (Bats: SB-vsRNA-ORF1a-3p; SB-vsRNA-S-5p), one vsRNA in MERS-CoV (Bats: MB-vsRNA-ORF1b-3p), and six vsRNAs in SARS-CoV-2 (Bats: S2B-vsRNA-ORF1a-5p; intermediate animals: S2I-vsRNA-ORF1a-5p; and humans: S2H-vsRNA-ORF1a-5p, S2H-vsRNA-ORF1a-3p, S2H-vsRNA-ORF1b-3p, S2H-vsRNA-ORF3a-3p), mainly encoded by nonstructural protein 3. Notably, Betacoronavirus-derived small RNAs targeted 74 differentially expressed genes in infected human cells, of which 55 upregulate the molecular mechanisms underlying acute respiratory distress syndrome (ARDS), and the 19 downregulated genes might be implicated in neurotrophin signaling impairment. These results reveal a novel small RNA-based regulatory mechanism involved in neuropathogenesis that must be further studied to validate its therapeutic use.

The Exosome-like Vesicles of Giardia Assemblages A, B, and E Are Involved in the Delivering of Distinct Small RNA from Parasite to Parasite.

The genetically related assemblages of the intestinal protozoa parasite Giardia lamblia are morphologically indistinguishable and are often derived from specific hosts. The Giardia assemblages are separated by large genetic distances, which might account for their relevant biological and pathogenic differences. In this work, we analyzed the RNAs cargo released into exosomal-like vesicles (ElVs) by the assemblages A and B, which differentially infect humans, and the assemblage E, which infects hoofed animals. The RNA sequencing analysis revealed that the ElVs of each assemblage contained distinct small RNA (sRNA) biotypes, suggesting a preference for specific packaging in each assemblage. These sRNAs were classified into three categories, ribosomal-small RNAs (rsRNAs), messenger-small RNAs (msRNAs), and transfer-small RNAs (tsRNAs), which may play a regulatory role in parasite communication and contribute to host-specificity and pathogenesis. Uptake experiments showed, for the first time, that ElVs were successfully internalized by the parasite trophozoites. Furthermore, we observed that the sRNAs contained inside these ElVs were first located below the plasma membrane but then distributed along the cytoplasm. Overall, the study provides new insights into the molecular mechanisms underlying the host-specificity and pathogenesis of G. lamblia and highlights the potential role of sRNAs in parasite communication and regulation.

Modulation of Pseudomonas aeruginosa quorum sensing by ajoene through direct competition with small RNAs for binding at the proximal site of Hfq - A structure-based perspective.

Bacteria can communicate to each other via quorum sensing, a cell density-dependent gene regulation system that stimulates the expression of virulence factors in the neighboring cells. Although the interaction of the natural product ajoene with the Hfq protein has been associated with the disruption of the quorum sensing system in Pseudomonas aeruginosa, there is no information concerning the corresponding ligand-target interaction process. Herein we observed a strong correlation (p < 0.00001) between the estimated affinities for the binding of 23 ajoene analogues at the proximal site of the Hfq protein of P. aeruginosa and their corresponding IC50 values, which reflect the reduction in the transcription of a virulence factor after quorum sensing inhibition. In this concern, our analyses reinforces previous propositions suggesting that ajoene could target the Hfq protein and affects its interaction with RNAs. Based on docking simulations, we tried to elucidate the binding mode of ajoene into the proximal Hfq site and we also established the minimum set of groups that would be necessary for a good interaction at this site, which includes a single hydrogen bond acceptor feature surrounded by groups that interact via π-sulfur (i.e., disulfide sulfurs) and/or π-alkyl/π-π stacking interactions (e.g., vinyl or small aryl/heteroaryl/heterocyclic groups). Because of the widespread role of Hfq as a matchmaker between messenger and small regulatory RNAs in Gram-negatives, we believe the discussion here provided for P. aeruginosa could be extrapolated for Gram-negatives in general, while the interaction of ajoene over the Hfq protein of Gram-positives would still remain more controversial.

Differences in Bacterial Small RNAs in Stool Samples from Hypercholesterolemic and Normocholesterolemic Subjects.

Cholesterol metabolism is important at the physiological level as well as in several diseases, with small RNA being an element to consider in terms of its epigenetic control. Thus, the aim of this study was to identify differences between bacterial small RNAs present at the gut level in hypercholesterolemic and normocholesterolemic individuals. Twenty stool samples were collected from hypercholesterolemic and normocholesterolemic subjects. RNA extraction and small RNA sequencing were performed, followed by bioinformatics analyses with BrumiR, Bowtie 2, BLASTn, DESeq2, and IntaRNA, after the filtering of the reads with fastp. In addition, the prediction of secondary structures was obtained with RNAfold WebServer. Most of the small RNAs were of bacterial origin and presented a greater number of readings in normocholesterolemic participants. The upregulation of small RNA ID 2909606 associated with Coprococcus eutactus (family Lachnospiraceae) was presented in hypercholesterolemic subjects. In addition, a positive correlation was established between small RNA ID 2149569 from the species Blautia wexlerae and hypercholesterolemic subjects. Other bacterial and archaeal small RNAs that interacted with the LDL receptor (LDLR) were identified. For these sequences, the prediction of secondary structures was also obtained. There were significant differences in bacterial small RNAs associated with cholesterol metabolism in hypercholesterolemic and normocholesterolemic participants.

Keeping up with the miRNAs: current paradigms of the biogenesis pathway.

For many years we have studied the processes involved in producing miRNAs in plants and the numerous differences from their metazoan counterpart. A well-defined catalytic process, mostly carried out by the RNase III enzyme DICER-LIKE1 (DCL1), it was identified early after the discovery of RNAi and was followed by the isolation of a plethora of miRNA biogenesis cofactors. The production of miRNAs, which later are loaded in ARGONAUTE (AGO) proteins to perform their RNA silencing functions both within the cell and non-cell autonomously, appears to be a highly regulated and dynamic process. Many regulatory events during miRNA biogenesis require the action of specific proteins. However, in recent years, many post-transcriptional modifications, structural features, and coupling with other cellular processing emerged as critical elements controlling the production of miRNA and, thus, a plant's physiology. This review discusses new evidence that has changed the way we understand how miRNAs are produced in plants. We also provide an updated view of the miRNA biogenesis pathways, focusing on the gaps in our knowledge and the most compelling questions that remain open.

Trichoderma atroviride hyphal regeneration and conidiation depend on cell-signaling processes regulated by a microRNA-like RNA.

The ability to respond to injury is essential for the survival of an organism and involves analogous mechanisms in animals and plants. Such mechanisms integrate coordinated genetic and metabolic reprogramming events requiring regulation by small RNAs for adequate healing of the wounded area. We have previously reported that the response to injury of the filamentous fungus Trichoderma atroviride involves molecular mechanisms closely resembling those of plants and animals that lead to the formation of new hyphae (regeneration) and the development of asexual reproduction structures (conidiophores). However, the involvement of microRNAs in this process has not been investigated in fungi. In this work, we explore the participation of microRNA-like RNAs (milRNAs) molecules by sequencing messenger and small RNAs during the injury response of the WT strain and RNAi mutants. We found that Dcr2 appears to play an important role in hyphal regeneration and is required to produce the majority of sRNAs in T. atroviride. We also determined that the three main milRNAs produced via Dcr2 are induced during the damage-triggered developmental process. Importantly, elimination of a single milRNA phenocopied the main defects observed in the dcr2 mutant. Our results demonstrate the essential role of milRNAs in hyphal regeneration and asexual development by post-transcriptionally regulating cellular signalling processes involving phosphorylation events. These observations allow us to conclude that fungi, like plants and animals, in response to damage activate fine-tuning regulatory mechanisms.

Altered RNome expression in Murine Gastrocnemius Muscle following Exposure to Jararhagin, a Metalloproteinase from Bothrops jararaca Venom.

Small RNAs (sRNAs) and microRNAs (miRNAs) are small endogenous noncoding single-stranded RNAs that regulate gene expression in eukaryotes. Experiments in mice and humans have revealed that a typical small RNA can affect the expression of a wide range of genes, implying that small RNAs function as global regulators. Here, we used small RNA deep sequencing to investigate how jararhagin, a metalloproteinase toxin produced from the venom of Bothrops jararaca, affected mmu-miRNAs expression in mice 2 hours (Jar 2hrs) and 24 hours (Jar 24hrs) after injection compared to PBS control. The findings revealed that seven mmu-miRNAs were substantially differentially expressed (p value (p (Corr) cut-off 0.05, fold change ≥ 2) at 2 hrs after jararhagin exposure and that the majority of them were upregulated when compared to PBS. In contrast to these findings, a comparison of Jar 24hrs vs. PBS 24hrs demonstrated that the majority of identified mmu-miRNAs were downregulated. Furthermore, the studies demonstrated that mmu-miRNAs can target the expression of several genes involved in the MAPK signaling pathway. The steady antithetical regulation of mmu-miRNAs may correlate with the expression of genes that trigger apoptosis via MAPK in the early stages, and this effect intensifies with time. The findings expand our understanding of the effects of jararhagin on local tissue lesions at the molecular level.

Comparative Profiling of Circulating Exosomal Small RNAs Derived From Peruvian Patients With Tuberculosis and Pulmonary Adenocarcinoma.

Tuberculosis (TB) is one of the most fatal infectious diseases, caused by the aerobic bacteria Mycobacterium tuberculosis. It is estimated that one-third of the world's population is infected with the latent (LTB) version of this disease, with only 5-10% of infected individuals developing its active (ATB) form. Pulmonary adenocarcinoma (PA) is the most common and diverse form of primary lung carcinoma. The simultaneous or sequential occurrence of TB and lung cancer in patients has been widely reported and is known to be an issue for diagnosis and surgical treatment. Raising evidence shows that patients cured of TB represent a group at risk for developing PA. In this work, using sRNA-sequencing, we evaluated the expression patterns of circulating small RNAs available in exosomes extracted from blood samples of Peruvian patients affected by latent tuberculosis, active tuberculosis, or pulmonary adenocarcinoma. Differential expression analysis revealed a set of 24 microRNAs perturbed in these diseases, revealing potential biomarker candidates for the Peruvian population. Most of these miRNAs are normally expressed in healthy lung tissue and are potential regulators of different shared and unique KEGG pathways related to cancers, infectious diseases, and immunology.

Altered RNome expression in murine gastrocnemius muscle following exposure to Jararhagin, a metalloproteinase from Bothrops jararaca venom

Small RNAs (sRNAs) and microRNAs (miRNAs) are small endogenous noncoding single-stranded RNAs that regulate gene expression in eukaryotes. Experiments in mice and humans have revealed that a typical small RNA can affect the expression of a wide range of genes, implying that small RNAs function as global regulators. Here, we used small RNA deep sequencing to investigate how jararhagin, a metalloproteinase toxin produced from the venom of Bothrops jararaca, affected mmu-miRNAs expression in mice 2 hours (Jar 2hrs) and 24 hours (Jar 24hrs) after injection compared to PBS control. The findings revealed that seven mmu-miRNAs were substantially differentially expressed (p value (p (Corr) cut-off 0.05, fold change ≥ 2) at 2 hrs after jararhagin exposure and that the majority of them were upregulated when compared to PBS. In contrast to these findings, a comparison of Jar 24hrs vs. PBS 24hrs demonstrated that the majority of identified mmu-miRNAs were downregulated. Furthermore, the studies demonstrated that mmu-miRNAs can target the expression of several genes involved in the MAPK signaling pathway. The steady antithetical regulation of mmu-miRNAs may correlate with the expression of genes that trigger apoptosis via MAPK in the early stages, and this effect intensifies with time. The findings expand our understanding of the effects of jararhagin on local tissue lesions at the molecular level.

Interspecies RNA Interactome of Pathogen and Host in a Heritable Defensive Strategy.

Communication with bacteria deeply impacts the life history traits of their hosts. Through specific molecules and metabolites, bacteria can promote short- and long-term phenotypic and behavioral changes in the nematode Caenorhabditis elegans. The chronic exposure of C. elegans to pathogens promotes the adaptive behavior in the host's progeny called pathogen-induced diapause formation (PIDF). PIDF is a pathogen avoidance strategy induced in the second generation of animals infected and can be recalled transgenerationally. This behavior requires the RNA interference machinery and specific nematode and bacteria small RNAs (sRNAs). In this work, we assume that RNAs from both species co-exist and can interact with each other. Under this principle, we explore the potential interspecies RNA interactions during PIDF-triggering conditions, using transcriptomic data from the holobiont. We study two transcriptomics datasets: first, the dual sRNA expression of Pseudomonas aeruginosa PAO1 and C. elegans in a transgenerational paradigm for six generations and second, the simultaneous expression of sRNAs and mRNA in intergenerational PIDF. We focus on those bacterial sRNAs that are systematically overexpressed in the intestines of animals compared with sRNAs expressed in host-naïve bacteria. We selected diverse in silico methods that represent putative mechanisms of RNA-mediated interspecies interaction. These interactions are as follows: heterologous perfect and incomplete pairing between bacterial RNA and host mRNA; sRNAs of similar sequence expressed in both species that could mimic each other; and known or predicted eukaryotic motifs present in bacterial transcripts. We conclude that a broad spectrum of tools can be applied for the identification of potential sRNA and mRNA targets of the interspecies RNA interaction that can be subsequently tested experimentally.

Genome-wide identification of miRNAs and target regulatory network in the invasive ectoparasitic mite Varroa destructor.

Varroa destructor is an ectoparasite mite that attacks bees leading to colony disorders worldwide. microRNAs (miRNAs) are key molecules used by eukaryotes to post-transcriptional control of gene expression. Nevertheless, still lack information aboutV. destructor miRNAs and its regulatory networks. Here, we used an integrative strategy to characterize the miRNAs in the V. destructor mite. We identified 310 precursors that give rise to 500 mature miRNAs, which 257 are likely mite-specific elements. miRNAs showed canonical length ranging between 18 and 25 nucleotides and 5' uracil preference. Top 10 elements concentrated over 80% of total miRNA expression, with bantam alone representing ~50%. We also detected non-templated bases in precursor-derived small RNAs, indicative of miRNA post-transcriptional regulatory mechanisms. Finally, we note that conserved miRNAs control similar processes in different organisms, suggesting a conservative role. Altogether, our findings contribute to the better understanding of the mite biology that can assist future studies on varroosis control.

The Heritability of Behaviors Associated With the Host Gut Microbiota.

What defines whether the interaction between environment and organism creates a genetic memory able to be transferred to subsequent generations? Bacteria and the products of their metabolism are the most ubiquitous biotic environments to which every living organism is exposed. Both microbiota and host establish a framework where environmental and genetic factors are integrated to produce adaptive life traits, some of which can be inherited. Thus, the interplay between host and microbe is a powerful model to study how phenotypic plasticity is inherited. Communication between host and microbe can occur through diverse molecules such as small RNAs (sRNAs) and the RNA interference machinery, which have emerged as mediators and carriers of heritable environmentally induced responses. Notwithstanding, it is still unclear how the organism integrates sRNA signaling between different tissues to orchestrate a systemic bacterially induced response that can be inherited. Here we discuss current evidence of heritability produced by the intestinal microbiota from several species. Neurons and gut are the sensing systems involved in transmitting changes through transcriptional and post-transcriptional modifications to the gonads. Germ cells express inflammatory receptors, and their development and function are regulated by host and bacterial metabolites and sRNAs thus suggesting that the dynamic interplay between host and microbe underlies the host's capacity to transmit heritable behaviors. We discuss how the host detects changes in the microbiota that can modulate germ cells genomic functions. We also explore the nature of the interactions that leave permanent or long-term memory in the host and propose mechanisms by which the microbiota can regulate the development and epigenetic reprogramming of germ cells, thus influencing the inheritance of the host. We highlight the vast contribution of the bacterivore nematode C. elegans and its commensal and pathogenic bacteria to the understanding on how behavioral adaptations can be inter and transgenerational inherited.

Network Analyses Predict Small RNAs That Might Modulate Gene Expression in the Testis and Epididymis of Bos indicus Bulls.

Spermatogenesis relies on complex molecular mechanisms, essential for the genesis and differentiation of the male gamete. Germ cell differentiation starts at the testicular parenchyma and finishes in the epididymis, which has three main regions: head, body, and tail. RNA-sequencing data of the testicular parenchyma (TP), head epididymis (HE), and tail epididymis (TE) from four bulls (three biopsies per bull: 12 samples) were subjected to differential expression analyses, functional enrichment analyses, and co-expression analyses. The aim was to investigate the co-expression and infer possible regulatory roles for transcripts involved in the spermatogenesis of Bos indicus bulls. Across the three pairwise comparisons, 3,826 differentially expressed (DE) transcripts were identified, of which 384 are small RNAs. Functional enrichment analysis pointed to gene ontology (GO) terms related to ion channel activity, detoxification of copper, neuroactive receptors, and spermatogenesis. Using the regulatory impact factor (RIF) algorithm, we detected 70 DE small RNAs likely to regulate the DE transcripts considering all pairwise comparisons among tissues. The pattern of small RNA co-expression suggested that these elements are involved in spermatogenesis regulation. The 3,826 DE transcripts (mRNAs and small RNAs) were further subjected to co-expression analyses using the partial correlation and information theory (PCIT) algorithm for network prediction. Significant correlations underpinned the co-expression network, which had 2,216 transcripts connected by 158,807 predicted interactions. The larger network cluster was enriched for male gamete generation and had 15 miRNAs with significant RIF. The miRNA bta-mir-2886 showed the highest number of connections (601) and was predicted to down-regulate ELOVL3, FEZF2, and HOXA13 (negative co-expression correlations and confirmed with TargetScan). In short, we suggest that bta-mir-2886 and other small RNAs might modulate gene expression in the testis and epididymis, in Bos indicus cattle.

When junk DNA turns functional: transposon-derived non-coding RNAs in plants.

Transposable elements (TEs) are major contributors to genome complexity in eukaryotes. TE mobilization may cause genome instability, although it can also drive genome diversity throughout evolution. TE transposition may influence the transcriptional activity of neighboring genes by modulating the epigenomic profile of the region or by altering the relative position of regulatory elements. Notably, TEs have emerged in the last few years as an important source of functional long and small non-coding RNAs. A plethora of small RNAs derived from TEs have been linked to the trans regulation of gene activity at the transcriptional and post-transcriptional levels. Furthermore, TE-derived long non-coding RNAs have been shown to modulate gene expression by interacting with protein partners, sequestering active small RNAs, and forming duplexes with DNA or other RNA molecules. In this review, we summarize our current knowledge of the functional and mechanistic paradigms of TE-derived long and small non-coding RNAs and discuss their role in plant development and evolution.