Small RNA Sequencing Data Research Articles

Post-Transcriptional Modifications to miRNAs Undergo Widespread Alterations, Creating a Unique Lung Adenocarcinoma IsomiRome.

MicroRNAs (miRNAs) modulate the expression of oncogenes and tumor suppressor genes, functioning as significant epigenetic regulators in cancer. IsomiRs are miRNA molecules that have undergone small modifications during miRNA processing. These modifications can alter an isomiR's binding stability with mRNA targets, and certain isomiRs have been implicated in the development of specific cancers. Still, the isomiRomes of many tissues, including the lung, have not been characterized; Methods: In this study, we analyzed small RNA sequencing data for three cohorts of lung adenocarcinoma (LUAD) and adult non-malignant lung (ANL) samples. We quantified isomiR expression and found 16 A-to-I edited isomiRs expressed in multiple cohorts, as well as 213 5' isomiRs, 128 3' adenylated isomiRs, and 100 3' uridylated isomiRs. Rates of A-to-I editing at editing hotspots correlated with mRNA expression of the editing enzymes ADAR and ADARB1, which were both observed to be deregulated in LUAD. LUAD samples displayed lower overall rates of A-to-I editing and 3' adenylation than ANL samples. Support vector machines and random forest models were trained on one cohort to distinguish ANL and stage I/II LUAD samples using reads per million (RPM) and frequency data for different types of isomiRs. Models trained on A-to-I editing rates at editing hotspots displayed high accuracy when tested on the other two cohorts and compared favorably to classifiers trained on miRNA expression alone; Conclusions: We have identified isomiRs in the human lung and found that their expression differs between non-malignant and tumor tissues, suggesting they hold potential as cancer biomarkers.

The Integrated Analysis of miRNome and Degradome Sequencing Reveals the Regulatory Mechanisms of Seed Development and Oil Biosynthesis in Pecan (Carya illinoinensis).

Pecan seed oil is a valuable source of essential fatty acids and various bioactive compounds; however, the functions of microRNAs and their targets in oil biosynthesis during seed development are still unknown. Here, we found that the oil content increased rapidly in the three early stages in three cultivars, and that oleic acid was the predominant fatty acid component in the mature pecan embryos. We identified, analyzed, and validated the expression levels of miRNAs related to seed development and oil biosynthesis, as well as their potential target genes, using small RNA sequencing data from three stages (120, 135, and 150 days after flowering). During the seed development process, 365 known and 321 novel miRNAs were discovered. In total, 91 known and 181 novel miRNAs were found to be differentially expressed, and 633 target genes were further investigated. The expression trend analysis revealed that the 91 known miRNAs were classified into eight groups, approximately two-thirds of which were up-regulated, whereas most novel miRNAs were down-regulated. The qRT-PCR and degradome sequencing data were used to identify five miRNA- target pairs. Overall, our study provides valuable insights into the molecular regulation of oil biosynthesis in pecan seeds.

MicroRNAs implicated in canine diffuse large B-cell lymphoma prognosis.

Diffuse large B-cell lymphoma (DLBCL) is the most prevalent subtype of non-Hodgkin lymphoma (NHL) in domestic dogs, with many similarities to its human counterpart. The progression of the disease is rapid, and treatment must be initiated early to achieve cancer remission and extend life. This study examined the relationship between progression-free survival (PFS) and microRNA (miRNA) expression in dogs with DLBCL. miRNAs are small non-coding RNA molecules that typically regulate gene expression post-transcriptionally. They are involved in several pathophysiological processes, including the growth and progression of cancer. Based on the analysis of small RNA sequencing (sRNA-seq) data, we validated a group of miRNAs in lymph nodes from 44 DLBCL-affected dogs with known outcomes. We used quantitative PCR to quantify their expression and report a specific subset of miRNAs is associated with decreased PFS in dogs with DLBCL. The miR-192-5p and miR-16-5p expression were significantly downregulated in dogs with increased PFS. These results indicate that miRNA profiling may potentially identify dogs with DLBCL that will experience poor outcomes following treatment. Identifying specific miRNAs that correlate with the progression of canine DLBCL could aid the development of individualized treatment regimens for dogs.

FlaHMM: unistrand flamenco-like piRNA cluster prediction in Drosophila species using hidden Markov models.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs that are essential for transposon control in animal gonads. In Drosophila ovarian somatic cells, piRNAs are transcribed from large genomic regions called piRNA clusters, which are enriched for transposon fragments and act as a memory of past invasions. Despite being widely present across Drosophila species, somatic piRNA clusters are difficult to identify and study due to their lack of sequence conservation and limited synteny. Current identification methods rely on either extensive manual curation or availability of high-throughput small RNA sequencing data, limiting large-scale comparative studies. We now present FlaHMM, a hidden Markov model developed to automate genomic annotation of flamenco-like unistrand piRNA clusters in Drosophila species, requiring only a genome assembly and transposon annotations. FlaHMM uses transposable element content across 5- or 10-kb bins, which can be calculated from genome sequence alone, and is thus able to detect candidate piRNA clusters without the need to obtain flies and experimentally perform small RNA sequencing. We show that FlaHMM performs on par with piRNA-guided or manual methods, and thus provides a scalable and efficient approach to piRNA cluster annotation in new genome assemblies. FlaHMM is freely available at https://github.com/Hannon-lab/FlaHMM under an MIT licence.

Analysis of Amblyomma americanum microRNAs in response to Ehrlichia chaffeensis infection and their potential role in vectorial capacity.

MicroRNAs (miRNAs) represent a subset of small noncoding RNAs and carry tremendous potential for regulating gene expression at the post-transcriptional level. They play pivotal roles in distinct cellular mechanisms including inhibition of bacterial, parasitic, and viral infections via immune response pathways. Intriguingly, pathogens have developed strategies to manipulate the host's miRNA profile, fostering environments conducive to successful infection. Therefore, changes in an arthropod host's miRNA profile in response to pathogen invasion could be critical in understanding host-pathogen dynamics. Additionally, this area of study could provide insights into discovering new targets for disease control and prevention. The main objective of the present study is to investigate the functional role of differentially expressed miRNAs upon Ehrlichia chaffeensis, a tick-borne pathogen, infection in tick vector, Amblyomma americanum. Small RNA libraries from uninfected and E. chaffeensis-infected Am. americanum midgut and salivary gland tissues were prepared using the Illumina Truseq kit. Small RNA sequencing data was analyzed using miRDeep2 and sRNAtoolbox to identify novel and known miRNAs. The differentially expressed miRNAs were validated using a quantitative PCR assay. Furthermore, a miRNA inhibitor approach was used to determine the functional role of selected miRNA candidates. The sequencing of small RNA libraries generated >147 million raw reads in all four libraries and identified a total of >250 miRNAs across the four libraries. We identified 23 and 14 differentially expressed miRNAs in salivary glands, and midgut tissues infected with E. chaffeensis, respectively. Three differentially expressed miRNAs (miR-87, miR-750, and miR-275) were further characterized to determine their roles in pathogen infection. Inhibition of target miRNAs significantly decreased the E. chaffeensis load in tick tissues, which warrants more in-depth mechanistic studies. The current study identified known and novel miRNAs and suggests that interfering with these miRNAs may impact the vectorial capacity of ticks to harbor Ehrlichia. This study identified several new miRNAs for future analysis of their functions in tick biology and tick-pathogen interaction studies.

Integrative Bioinformatics Analysis Reveals a Transcription Factor EB-Driven MicroRNA Regulatory Network in Endothelial Cells.

Various human diseases are triggered by molecular alterations influencing the fine-tuned expression and activity of transcription factors, usually due to imbalances in targets including protein-coding genes and non-coding RNAs, such as microRNAs (miRNAs). The transcription factor EB (TFEB) modulates human cellular networks, overseeing lysosomal biogenesis and function, plasma-membrane trafficking, autophagic flux, and cell cycle progression. In endothelial cells (ECs), TFEB is essential for the maintenance of endothelial integrity and function, ensuring vascular health. However, the comprehensive regulatory network orchestrated by TFEB remains poorly understood. Here, we provide novel mechanistic insights into how TFEB regulates the transcriptional landscape in primary human umbilical vein ECs (HUVECs), using an integrated approach combining high-throughput experimental data with dedicated bioinformatics analysis. By analyzing HUVECs ectopically expressing TFEB using ChIP-seq and examining both polyadenylated mRNA and small RNA sequencing data from TFEB-silenced HUVECs, we have developed a bioinformatics pipeline mapping the different gene regulatory interactions driven by TFEB. We show that TFEB directly regulates multiple miRNAs, which in turn post-transcriptionally modulate a broad network of target genes, significantly expanding the repertoire of gene programs influenced by this transcription factor. These insights may have significant implications for vascular biology and the development of novel therapeutics for vascular disease.

Restricting Colorectal Cancer Cell Metabolism with Metformin: An Integrated Transcriptomics Study.

Metformin is a first-line therapy for type 2 diabetes as it disrupts cellular metabolism. Despite the association between metformin and lower cancer incidence, the anti-tumour activity of the drug in colorectal cancer (CRC) is incompletely understood. This study identifies underlying molecular mechanisms by which metformin slows colorectal cancer cell proliferation by investigating metformin-associated microRNA (miRNA) and target gene pairs implicated in signalling pathways. The present study analysed changes in miRNAs and the coding transcriptome in CRC cells treated with a sublethal dose of metformin, followed by the contextual validation of potential miRNA-target gene pairs. Analyses of small RNA and transcriptome sequencing data revealed 104 miRNAs and 1221 mRNAs to be differentially expressed in CRC cells treated with metformin for 72 h. Interaction networks between differentially expressed miRNAs and putative target mRNAs were identified. Differentially expressed genes were mainly implicated in metabolism and signalling processes, such as the PI3K-Akt and MAPK/ERK pathways. Further validation of potential miRNA-target mRNA pairs revealed that metformin induced miR-2110 and miR-132-3p to target PIK3R3 and, consequently, regulate CRC cell proliferation, cell cycle progression and the PI3K-Akt signalling pathway. Metformin also induced miR-222-3p and miR-589-3p, which directly target STMN1 to inhibit CRC cell proliferation and cell cycle progression. This study identified novel changes in the coding transcriptome and small non-coding RNAs associated with metformin treatment of CRC cells. Integration of these datasets highlighted underlying mechanisms by which metformin impedes cell proliferation in CRC. Importantly, it identified the post-transcriptional regulation of specific genes that impact both metabolism and cell proliferation.

Inhibition of miR-199b-5p reduces pathological alterations in osteoarthritis by potentially targeting Fzd6 and Gcnt2

Osteoarthritis (OA) is a degenerative disease with a high prevalence in the elderly population, but our understanding of its mechanisms remains incomplete. Analysis of serum exosomal small RNA sequencing data from clinical patients and gene expression data from OA patient serum and cartilage obtained from the GEO database revealed a common dysregulated miRNA, miR-199b-5p. In vitro cell experiments demonstrated that miR-199b-5p inhibits chondrocyte vitality and promotes extracellular matrix degradation. Conversely, inhibition of miR-199b-5p under inflammatory conditions exhibited protective effects against damage. Local viral injection of miR-199b-5p into mice induced a decrease in pain threshold and OA-like changes. In an OA model, inhibition of miR-199b-5p alleviated the pathological progression of OA. Furthermore, bioinformatics analysis and experimental validation identified Gcnt2 and Fzd6 as potential target genes of MiR-199b-5p. Thus, these results indicated that MiR-199b-5p/Gcnt2 and Fzd6 axis might be a novel therapeutic target for the treatment of OA.

Inhibition of miR-199b-5p reduces pathological alterations in osteoarthritis by potentially targeting Fzd6 and Gcnt2.

Osteoarthritis (OA) is a degenerative disease with a high prevalence in the elderly population, but our understanding of its mechanisms remains incomplete. Analysis of serum exosomal small RNA sequencing data from clinical patients and gene expression data from OA patient serum and cartilage obtained from the GEO database revealed a common dysregulated miRNA, miR-199b-5p. In vitro cell experiments demonstrated that miR-199b-5p inhibits chondrocyte vitality and promotes extracellular matrix degradation. Conversely, inhibition of miR-199b-5p under inflammatory conditions exhibited protective effects against damage. Local viral injection of miR-199b-5p into mice induced a decrease in pain threshold and OA-like changes. In an OA model, inhibition of miR-199b-5p alleviated the pathological progression of OA. Furthermore, bioinformatics analysis and experimental validation identified Gcnt2 and Fzd6 as potential target genes of MiR-199b-5p. Thus, these results indicated that MiR-199b-5p/Gcnt2 and Fzd6 axis might be a novel therapeutic target for the treatment of OA.

How microRNAs go awry in cancer: Exploring the underlying factors of miRNA dysregulation

MicroRNAs (miRNA) are short, non-coding molecules (~23nt) that fine-tune mRNA expression at the post-transcriptional level by repressing or cleaving their targets. A single miRNA targets several mRNAs; dysfunctions disrupt biological systems. Unravelling the interactome of regulators and targets is vital for biological networks. Cancer research is one of the hottest fields, due to its high prevalence and financial burden. Despite improvements in prevention and therapy, challenges from disease complexity and precision medicine limit the effciency of targeted approaches. Accurate characterization of all interactions and mechanisms is essential to fight it. In this study, we matched whole genome, RNA and small RNA sequencing data acquired from the International Cancer Genome Consortium (ICGC) for Liver cancer patients, resulting in a dataset of 40 samples with the aim to explain the dysregulation of miRNA expression. DNA variants were readily extracted using the “Broad variant call pipeline” of ICGC. Utilization of DNase-seq data indicated an average promoter region activity of 1500bp upstream and 500bp downstream of the transcription start site. Variations were then intersected to promoter regions. Analysis of transcription factors (TF), which regulate the expression of DNA bound genes, as well as miRNAs, was done using FIMO prediction algorithm (25bp proximal window centered on identified variants, p-value < 0.000001) combined with JASPAR transcription factor motifs to identify mutated transcription factor binding sites. Matching RNA-seq expression was then combined with small RNA-seq expression and correlation analysis was performed so as to identify dysregulated TF::miRNA interacting pairs with mutations on their regulatory elements, on the basis that their correlated expression is potentially disrupted due to overlapping variation events. A total of 172691 variants events were identified in liver promoter regions. Analysis of mutational load using chi-square test in various gene biotypes indicated differences with protein-coding and small RNA (excluding miRNA) categories exhibiting the higher rates. Variants in TF genes generally showed low frequency (3.64% of all TF genes analyzed), necessitating a deeper analysis to discover if the underling mechanisms of TFBS are affected by variation events. Of all TFs utilized in the FIMO analysis (n=39), ZNF263, which has been associated to a variety of tumours, was found to have the highest percentage of mutated binding motif (31%) followed by ZNF384 (28%). Pearson correlation analysis between the expression of miRNAs and TFs with and without mutations overlapping TF binding sites in miRNA promoters identified HNF4A (Hepatocyte Nuclear Factor 4) and hsa-mir-122 having significant difference in their correlated expression (Transcripts per million normalized values) between wildtype samples (n=20, r=0.8, p=0.00002) and mutated individuals (n=20, r=0.25, p=026), with both molecules being associated with liver localized functions in the literature. These results demonstrate the importance of further enriching the interactome towards a more complete network in order to achieve a more effective and accurate precision medicine. Hellenic Foundation for Research and Innovation (H.F.R.I.) under the ‘1st Call for H.F.R.I. Research Projects to support Faculty Members & Researchers and the procurement of high-cost research equipment grant’ [2563]. Funding for open access charge. This is the full abstract presented at the American Physiology Summit 2024 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.

Multi-omics integration identifies regulatory factors underlying bovine subclinical mastitis

BackgroundMastitis caused by multiple factors remains one of the most common and costly disease of the dairy industry. Multi-omics approaches enable the comprehensive investigation of the complex interactions between multiple layers of information to provide a more holistic view of disease pathogenesis. Therefore, this study investigated the genomic and epigenomic signatures and the possible regulatory mechanisms underlying subclinical mastitis by integrating RNA sequencing data (mRNA and lncRNA), small RNA sequencing data (miRNA) and DNA methylation sequencing data of milk somatic cells from 10 healthy cows and 20 cows with naturally occurring subclinical mastitis caused by Staphylococcus aureus or Staphylococcus chromogenes.ResultsFunctional investigation of the data sets through gene set analysis uncovered 3458 biological process GO terms and 170 KEGG pathways with altered activities during subclinical mastitis, provided further insights into subclinical mastitis and revealed the involvement of multi-omics signatures in the altered immune responses and impaired mammary gland productivity during subclinical mastitis. The abundant genomic and epigenomic signatures with significant alterations related to subclinical mastitis were observed, including 30,846, 2552, 1276 and 57 differential methylation haplotype blocks (dMHBs), differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs), respectively. Next, 5 factors presenting the principal variation of differential multi-omics signatures were identified. The important roles of Factor 1 (DEG, DEM and DEL) and Factor 2 (dMHB and DEM), in the regulation of immune defense and impaired mammary gland functions during subclinical mastitis were revealed. Each of the omics within Factors 1 and 2 explained about 20% of the source of variation in subclinical mastitis. Also, networks of important functional gene sets with the involvement of multi-omics signatures were demonstrated, which contributed to a comprehensive view of the possible regulatory mechanisms underlying subclinical mastitis. Furthermore, multi-omics integration enabled the association of the epigenomic regulatory factors (dMHBs, DELs and DEMs) of altered genes in important pathways, such as ‘Staphylococcus aureus infection pathway’ and ‘natural killer cell mediated cytotoxicity pathway’, etc., which provides further insights into mastitis regulatory mechanisms. Moreover, few multi-omics signatures (14 dMHBs, 25 DEGs, 18 DELs and 5 DEMs) were identified as candidate discriminant signatures with capacity of distinguishing subclinical mastitis cows from healthy cows.ConclusionThe integration of genomic and epigenomic data by multi-omics approaches in this study provided a better understanding of the molecular mechanisms underlying subclinical mastitis and identified multi-omics candidate discriminant signatures for subclinical mastitis, which may ultimately lead to the development of more effective mastitis control and management strategies.

5'-Isoforms of miR-1246 Have Distinct Targets and Stronger Functional Impact Compared with Canonical miR-1246 in Colorectal Cancer Cells In Vitro.

Colorectal cancer (CRC) is a multifactorial disease involving genetic and epigenetic factors, such as miRNAs. Sequencing-based studies have revealed that miRNAs have many isoforms (isomiRs) with modifications at the 3'- and 5'-ends or in the middle, resulting in distinct targetomes and, consequently, functions. In the present study, we aimed to evaluate the putative targets and functional role of miR-1246 and its two 5'-isoforms (ISO-miR-1246_a and ISO-miR-1246_G) in vitro. Commercial Caco-2 cells of CRC origin were analyzed for the expression of WT-miR-1246 and its 5'-isoforms using small RNA sequencing data, and the overabundance of the two miR-1246 isoforms was determined in cells. The transcriptome analysis of Caco-2 cells transfected with WT-miR-1246, ISO-miR-1246_G, and ISO-miR-1246_a indicated the minor overlap of the targetomes between the studied miRNA isoforms. Consequently, an enrichment analysis showed the involvement of the potential targets of the miR-1246 isoforms in distinct signaling pathways. Cancer-related pathways were predominantly more enriched in dysregulated genes in ISO-miR-1246_G and ISO-miR-1246_a, whereas cell cycle pathways were more enriched in WT-miR-1246. The functional analysis of WT-miR-1246 and its two 5'-isoforms revealed that the inhibition of any of these molecules had a tumor-suppressive role (reduced cell viability and migration and promotion of early cell apoptosis) in CRC cells. However, the 5'-isoforms had a stronger effect on viability compared with WT-miR-1246. To conclude, this research shows that WT-miR-1246 and its two 5'-isoforms have different targetomes and are involved in distinct signaling pathways but collectively play an important role in CRC pathogenesis.

MiR-708–5p deficiency involves the degeneration of mandibular condylar chondrocytes via the TLR4/NF-κB pathway

ObjectiveAgeing and aberrant biomechanical stimulation are two major risk factors for osteoarthritis (OA). One of the main characteristics of aged cartilage is cellular senescence. One of the main characteristics of osteoarthritic joints is cartilage degeneration. The cells in the temporomandibular joint (TMJ) cartilage are zonally arranged. The deep zone cells are differentiated from the superficial zone cells (SZCs). The purpose of the present study was to investigate whether degenerative shear stress (SS) stimulates the senescence programme in TMJ SZCs, and to determine which miRNA is involved in this process. MethodSZCs were isolated from the TMJ condyles of 3-week-old rats and treated with continuous passaging or SS. RNA sequencing was conducted to identify miRNA(s) that overlap with those involved in the replication senescence process and the SS-induced degeneration programme. Unilateral anterior crossbite (UAC), which is TMJ-OA inducible, was applied to 2-month-old and 12-month-old mice for 3 weeks. The effect of TMJ local injection of agomiR-708-5p was evaluated histologically. ResultsBoth replication and SS treatment induced SZC senescence. miR-708-5p was identified. Knocking down miR-708-5p in SS treated SZCs led to more severe senescence by alleviating the inhibitory impact of miR-708-5p on the TLR4/NF-κB pathway. miR-708-5p expression in mouse TMJ cartilage decreased with age. UAC induced more severe osteoarthritic cartilage lesions in 12-month-old mice than in 2-month-old mice. Injection of agomiR-708-5p suppressed UAC-induced osteoarthritic cartilage lesions. ConclusionsAge-related miR-708-5p deficiency is involved in the mechanically stimulated OA process. Intra-articular administration of agomiR-708-5p is a promising new strategy for OA treatment. Data availabilityThe small RNA and mRNA sequencing data included in this study are available in the NCBI Gene Expression Omnibus (GEO) database (accession numbers GSE236980, GSE236977, and GSE 236978).

A novel saliva-based miRNA profile to diagnose and predict oral cancer

Oral cancer (OC) is the most common form of head and neck cancer. Despite the high incidence and unfavourable patient outcomes, currently, there are no biomarkers for the early detection of OC. This study aims to discover, develop, and validate a novel saliva-based microRNA signature for early diagnosis and prediction of OC risk in oral potentially malignant disorders (OPMD). The Cancer Genome Atlas (TCGA) miRNA sequencing data and small RNA sequencing data of saliva samples were used to discover differentially expressed miRNAs. Identified miRNAs were validated in saliva samples of OC (n = 50), OPMD (n = 52), and controls (n = 60) using quantitative real-time PCR. Eight differentially expressed miRNAs (miR-7-5p, miR-10b-5p, miR-182-5p, miR-215-5p, miR-431-5p, miR-486-3p, miR-3614-5p, and miR-4707-3p) were identified in the discovery phase and were validated. The efficiency of our eight-miRNA signature to discriminate OC and controls was: area under curve (AUC): 0.954, sensitivity: 86%, specificity: 90%, positive predictive value (PPV): 87.8% and negative predictive value (NPV): 88.5% whereas between OC and OPMD was: AUC: 0.911, sensitivity: 90%, specificity: 82.7%, PPV: 74.2% and NPV: 89.6%. We have developed a risk probability score to predict the presence or risk of OC in OPMD patients. We established a salivary miRNA signature that can aid in diagnosing and predicting OC, revolutionising the management of patients with OPMD. Together, our results shed new light on the management of OC by salivary miRNAs to the clinical utility of using miRNAs derived from saliva samples.

FindAdapt: A python package for fast and accurate adapter detection in small RNA sequencing.

Adapter trimming is an essential step for analyzing small RNA sequencing data, where reads are generally longer than target RNAs ranging from 18 to 30 bp. Most adapter trimming tools require adapter information as input. However, adapter information is hard to access, specified incorrectly, or not provided with publicly available datasets, hampering their reproducibility and reusability. Manual identification of adapter patterns from raw reads is labor-intensive and error-prone. Moreover, the use of randomized adapters to reduce ligation biases during library preparation makes adapter detection even more challenging. Here, we present FindAdapt, a Python package for fast and accurate detection of adapter patterns without relying on prior information. We demonstrated that FindAdapt was far superior to existing approaches. It identified adapters successfully in 180 simulation datasets with diverse read structures and 3,184 real datasets covering a variety of commercial and customized small RNA library preparation kits. FindAdapt is stand-alone software that can be easily integrated into small RNA sequencing analysis pipelines.

SRNAflow: A Tool for the Analysis of Small RNA-Seq Data.

The analysis of small RNA sequencing data across a range of biofluids is a significant research area, given the diversity of RNA types that hold potential diagnostic, prognostic, and predictive value. The intricate task of segregating the complex mixture of small RNAs from both human and other species, including bacteria, fungi, and viruses, poses one of the most formidable challenges in the analysis of small RNA sequencing data, currently lacking satisfactory solutions. This study introduces sRNAflow, a user-friendly bioinformatic tool with a web interface designed for the analysis of small RNAs obtained from biological fluids. Tailored to the unique requirements of such samples, the proposed pipeline addresses various challenges, including filtering potential RNAs from reagents and environment, classifying small RNA types, managing small RNA annotation overlap, conducting differential expression assays, analysing isomiRs, and presenting an approach to identify the sources of small RNAs within samples. sRNAflow also encompasses an alternative alignment-free analysis of RNA-seq data, featuring clustering and initial RNA source identification using BLAST. This comprehensive approach facilitates meaningful comparisons of results between different analytical methods.

Milk exosomal microRNA profiling identified miR-375 and miR-199-5p for regulation of immune response during subclinical mastitis of crossbred cattle.

The dairy industry has experienced significant economic losses as a result of mastitis, an inflammatory disease of cows, including both subclinical and clinical cases. Milk exosome microRNAs have gained attention due to their stable and selective wrapping nature, offering potential for the prognosis and diagnosis of bovine mastitis, the most common pathological condition of the mammary gland. In the present investigation, the microRNA profile of milk exosomes was explored using high-throughput small RNA sequencing data in sub-clinical mastitic and healthy crossbred Vrindavani cattle. In both groups, 349 microRNAs were identified, with 238 (68.19%) microRNAs co-expressed; however, 35 and 76 distinct microRNAs were found in subclinical mastitic and healthy cattle, respectively. Differential expression analysis revealed 11 microRNAs upregulated, and 18 microRNAs were downregulated in sub-clinical mastitic cattle. The functional annotation of the target genes of differentially expressed known and novel microRNAs including bta-miR-375, bta-miR-199-5p and bta-miR-12030 reveals their involvement in the regulation of immune response and inflammatory mechanisms and could be involved in development of mastitis. The analysis of milk exosomal miRNAs cargos hold great promise as an approach to study the underlying molecular mechanisms associated with mastitis in high milk producing dairy cattle. Concurrently, the significantly downregulated miR-375 may upregulate key target genes, including CTLA4, IHH, IRF1, and IL7R. These genes are negative regulators of immune response pathways, which could be associated with impaired inflammatory mechanisms in mammary cells. According to the findings, bta-miR-375 could be a promising biomarker for the development of mastitis in dairy cattle.

Identification of Virus-Derived Small Interfering RNAs (vsiRNAs) from Infected sRNA-Seq Samples.

Plants have developed sophisticated defense mechanisms to combat viral infections, prominently utilizing Dicer-like enzymes (DCL) for generating virus-derived small interfering RNAs (vsiRNAs) through RNA interference (RNAi). This intrinsic mechanism effectively impedes virus replication. Exploiting their potential, vsiRNAs have become a major focus area for comprehensive viral investigations in plants, integrating both bioinformatics and experimental strategies. This chapter introduces an up-to-date computational workflow optimized for identifying and comprehensively annotating vsiRNAs with the utilization of small RNA sequencing (sRNA-seq) data collected from virus-infected plants. The workflow detailed in this chapter centers on known plant-targeting viruses, providing step-by-step guidance to enhance vsiRNA analysis, ultimately advancing the comprehension of plant-virus interactions.

New plasma diagnostic markers for colorectal cancer: transporter fragments of glutamate tRNA origin.

Colorectal cancer (CRC) is the second leading cause of cancer-related deaths worldwide. Early diagnosis of the disease can greatly improve the clinical prognosis for patients with CRC. Unfortunately, there are no current simple and effective early diagnostic markers available. The transfer RNA (tRNA)-derived RNA fragments (tRFs) are a class of small non-coding RNAs (sncRNAs), which have been shown to play an important role in the development and prognosis of CRC. However, only a few studies on tRFs as early diagnostic markers in CRC have been conducted. In this study, previously ignored tRFs expression data were extracted from six paired small RNA sequencing data in the Sequence Read Archive (SRA) database using MINTmap. Three i-tRFs, derived from the tRNA that transports glutamate (i-tRF-Glu), were identified and used to construct a random forest diagnostic model. The model performance was evaluated using the receiver operating characteristic (ROC) curve and precision-recall (PR) curve. The area under the curves (AUC) for the ROC and PR was 0.941 and 0.944, respectively. We further verified the differences in expression of the these i-tRF-Glu in the tissue and plasma of both CRC patients and healthy subjects using quantitative real-time PCR (qRT-PCR). We found that the ROC-AUC of the three was greater than traditional plasma tumor markers such as CEA and CA199. Our bioinformatics analysis suggested that the these i-tRF-Glu are associated with cancer development and glutamate (Glu)-glutamine (Gln) metabolism. Overall, our study uncovered these i-tRF-Glu that have early diagnostic significance and therapeutic potential for CRC, this warrants further investigation into the diagnostic and therapeutic potential of these i-tRF-Glu in CRC.

Differential microRNA editing may drive target pathway switching in human temporal lobe epilepsy.

MicroRNAs have emerged as important regulators of the gene expression landscape in temporal lobe epilepsy. The mechanisms that control microRNA levels and influence target choice remain, however, poorly understood. RNA editing is a post-transcriptional mechanism mediated by the adenosine acting on RNA (ADAR) family of proteins that introduces base modification that diversifies the gene expression landscape. RNA editing has been studied for the mRNA landscape but the extent to which microRNA editing occurs in human temporal lobe epilepsy is unknown. Here, we used small RNA-sequencing data to characterize the identity and extent of microRNA editing in human temporal lobe epilepsy brain samples. This detected low-to-high editing in over 40 of the identified microRNAs. Among microRNA exhibiting the highest editing was miR-376a-3p, which was edited in the seed region and this was predicted to significantly change the target pool. The edited form was expressed at lower levels in human temporal lobe epilepsy samples. We modelled the shift in editing levels of miR-376a-3p in human-induced pluripotent stem cell-derived neurons. Reducing levels of the edited form of miR-376a-3p using antisense oligonucleotides resulted in extensive gene expression changes, including upregulation of mitochondrial and metabolism-associated pathways. Together, these results show that differential editing of microRNAs may re-direct targeting and result in altered functions relevant to the pathophysiology of temporal lobe epilepsy and perhaps other disorders of neuronal hyperexcitability.