Small non-coding RNA profiling in human biofluids and surrogate tissues from healthy individuals: description of the diverse and most represented species

Circulating serum profile of small non-coding RNAs in patients with anaphylaxis beyond microRNAs

Introduction

Anaphylaxis is the most severe manifestation of allergic disorders. Currently, an increasing number of cells, pathways and molecules involved in the etiopathogenesis of anaphylaxis are being discovered. However, there are no conclusive biomarkers to confirm its diagnosis. Small non-coding RNAs (sncRNAs) are 18-200 nucleotide molecules that can be divided into: microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), small nucleolar RNAs (snoRNAs), small nuclear RNAs (snRNAs), transference RNA derived fragments (tRFs) and YRNA derived fragments (YRFs). These molecules participate in cell-cell communication modulating various physiological processes and have been postulated as non-invasive biomarkers of several pathologies. Therefore, in this study we characterized the serum circulating profile of other sncRNA beyond miRNAs in two populations of 5 adults and 5 children with drug- and food-mediated anaphylaxis, respectively.

Methods

Samples were obtained from each patient under two different conditions: during anaphylaxis and 14 days after the reaction (control). The sncRNA analysis was carried out by Next Generation Sequencing (NGS).

Results

A total of 671 sncRNAs (3 piRNAs, 74 snoRNAs, 54 snRNAs, 348 tRFs and 192 YRFs) were identified in adults with drug-induced anaphylaxis, while 612 sncRNAs (2 piRNAs, 73 snoRNAs, 52 snRNAs, 321 tRFs and 164 YRFs) were characterized in children with food-mediated anaphylaxis. However, only 33 (1 piRNA, 4 snoRNAs, 1 snRNAs, 7 tRFs and 20 YRFs) and 80 (4 snoRNAs, 6 snRNAs, 54 tRFs and 16 YRFs) of them were statistically different between both conditions, respectively. Among them, only three (Y_RNA.394, Y_RNA.781 and SCARNA2) were common to both adults and children analysis.

Discussion

This study provides a differential profile of circulating serum sncRNAs beyond miRNAs in patients with anaphylaxis, postulating them as candidate biomarkers for this pathological event and as novel mediators of the reaction.

Characterization of the internal working-life exposome using minimally and non-invasive sampling methods - a narrative review

During recent years, we are moving away from the 'one exposure, one disease'-approach in occupational settings and towards a more comprehensive approach, taking into account the totality of exposures during a life course by using an exposome approach. Taking an exposome approach however is accompanied by many challenges, one of which, for example, relates to the collection of biological samples. Methods used for sample collection in occupational exposome studies should ideally be minimally invasive, while at the same time sensitive, and enable meaningful repeated sampling in a large population and over a longer time period. This might be hampered in specific situations e.g., people working in remote areas, during pandemics or with flexible work hours. In these situations, using self-sampling techniques might offer a solution. Therefore, our aim was to identify existing self-sampling techniques and to evaluate the applicability of these techniques in an occupational exposome context by conducting a literature review. We here present an overview of current self-sampling methodologies used to characterize the internal exposome. In addition, the use of different biological matrices was evaluated and subdivided based on their level of invasiveness and applicability in an occupational exposome context. In conclusion, this review and the overview of self-sampling techniques presented herein can serve as a guide in the design of future (occupational) exposome studies while circumventing sample collection challenges associated with exposome studies.

Profiling small RNAs in fecal immunochemical tests: is it possible?

Fecal microRNAs represent promising molecules with potential clinical interest as non-invasive diagnostic and prognostic biomarkers. Colorectal cancer (CRC) screening based on the fecal immunochemical test (FIT) is an effective tool for prevention of cancer development. However, due to the poor sensitivity of FIT especially for premalignant lesions, there is a need for implementation of complementary tests. Improving the identification of individuals who would benefit from further investigation with colonoscopy using molecular analysis, such as miRNA profiling of FIT samples, would be ideal due to their widespread use. In the present study, we assessed the feasibility of applying small RNA sequencing to measure human miRNAs in FIT leftover buffer in samples from two European screening populations. We showed robust detection of miRNAs with profiles similar to those obtained from specimens sampled using the established protocol of RNA stabilizing buffers, or in long-term archived samples. Detected miRNAs exhibited differential abundances for CRC, advanced adenoma, and control samples that were consistent for FIT and RNA-stabilizing buffers. Interestingly, the sequencing data also allowed for concomitant evaluation of small RNA-based microbial profiles. We demonstrated that it is possible to explore the human miRNome in FIT leftover samples across populations and envision that the analysis of small RNA biomarkers can complement the FIT in large scale screening settings.

A Fecal MicroRNA Signature by Small RNA Sequencing Accurately Distinguishes Colorectal Cancers: Results From a Multicenter Study

BACKGROUND & AIMS

Fecal tests currently used for colorectal cancer (CRC) screening show limited accuracy in detecting early tumors or precancerous lesions. In this respect, we comprehensively evaluated stool microRNA (miRNA) profiles as biomarkers for noninvasive CRC diagnosis.

METHODS

A total of 1273 small RNA sequencing experiments were performed in multiple biospecimens. In a cross-sectional study, miRNA profiles were investigated in fecal samples from an Italian and a Czech cohort (155 CRCs, 87 adenomas, 96 other intestinal diseases, 141 colonoscopy-negative controls). A predictive miRNA signature for cancer detection was defined by a machine learning strategy and tested in additional fecal samples from 141 CRC patients and 80 healthy volunteers. miRNA profiles were compared with those of 132 tumors/adenomas paired with adjacent mucosa, 210 plasma extracellular vesicle samples, and 185 fecal immunochemical test leftover samples.

RESULTS

Twenty-five miRNAs showed altered levels in the stool of CRC patients in both cohorts (adjusted P < .05). A 5-miRNA signature, including miR-149-3p, miR-607-5p, miR-1246, miR-4488, and miR-6777-5p, distinguished patients from control individuals (area under the curve [AUC], 0.86; 95% confidence interval [CI], 0.79-0.94) and was validated in an independent cohort (AUC, 0.96; 95% CI, 0.92-1.00). The signature classified control individuals from patients with low-/high-stage tumors and advanced adenomas (AUC, 0.82; 95% CI, 0.71-0.97). Tissue miRNA profiles mirrored those of stool samples, and fecal profiles of different gastrointestinal diseases highlighted miRNAs specifically dysregulated in CRC. miRNA profiles in fecal immunochemical test leftover samples showed good correlation with those of stool collected in preservative buffer, and their alterations could be detected in adenoma or CRC patients.

CONCLUSIONS

Our comprehensive fecal miRNome analysis identified a signature accurately discriminating cancer aimed at improving noninvasive diagnosis and screening strategies.

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.

Gut Microbiome and Small RNA Integrative-Omic Perspective of Meconium and Milk-FED Infant Stool Samples

The human gut microbiome plays an important role in health, and its initial development is conditioned by many factors, such as feeding. It has also been claimed that this colonization is guided by bacterial populations, the dynamic virome, and transkingdom interactions between host and microbial cells, partially mediated by epigenetic signaling. In this article, we characterized the bacteriome, virome, and smallRNome and their interaction in the meconium and stool samples from infants. Bacterial and viral DNA and RNA were extracted from the meconium and stool samples of 2- to 4-month-old milk-fed infants. The bacteriome, DNA and RNA virome, and smallRNome were assessed using 16S rRNA V4 sequencing, viral enrichment sequencing, and small RNA sequencing protocols, respectively. Data pathway analysis and integration were performed using the R package mixOmics. Our findings showed that the bacteriome differed among the three groups, while the virome and smallRNome presented significant differences, mainly between the meconium and stool of milk-fed infants. The gut environment is rapidly acquired after birth, and it is highly adaptable due to the interaction of environmental factors. Additionally, transkingdom interactions between viruses and bacteria can influence host and smallRNome profiles. However, virome characterization has several protocol limitations that must be considered.

The 8q24 region hosts miRNAs altered in biospecimens of colorectal and bladder cancer patients

BACKGROUND

The 8q24 locus is enriched in cancer-associated polymorphisms and, despite containing relatively few protein-coding genes, it hosts the MYC oncogene and other genetic elements connected to tumorigenesis, including microRNAs (miRNAs). Research on miRNAs may provide insights into the transcriptomic regulation of this multiple cancer-associated region.

MATERIAL AND METHODS

We profiled all miRNAs located in the 8q24 region in 120 colorectal cancer (CRC) patients and 80 controls. miRNA profiling was performed on cancer/non-malignant adjacent mucosa, stool, and plasma extracellular vesicles (EVs), and the results validated with The Cancer Genome Atlas (TCGA) data. To verify if the 8q24-annotated miRNAs altered in CRC were dysregulated in other cancers and biofluids, we evaluated their levels in bladder cancer (BC) cases from the TCGA dataset and in urine and plasma EVs from a set of BC cases and healthy controls.

RESULTS

Among the detected mature miRNAs in the region, 12 were altered between CRC and adjacent mucosa (adj. p < 0.05). Five and four miRNAs were confirmed as dysregulated in the CRC and BC TCGA dataset, respectively. A co-expression analysis of tumor/adjacent tissue data from the CRC group revealed a correlation between the dysregulated miRNAs and CRC-related genes (PVT1 and MYC) annotated in 8q24 region. miR-30d-5p and miR-151a-3p, altered in CRC tissue, were also dysregulated in stool of CRC patients and urine of BC cases, respectively. Functional enrichment of dysregulated miRNA target genes highlighted terms related to TP53-mediated cell cycle control.

CONCLUSIONS

Altered expression of 8q24-annotated miRNAs may be relevant for the initiation and/or progression of cancer.

Insight on Non-Coding RNAs from Biofluids in Ovarian Tumors

Ovarian tumors are the most frequent adnexal mass, raising diagnostic and therapeutic issues linked to a large spectrum of tumors, with a continuum from benign to malignant. Thus far, none of the available diagnostic tools have proven efficient in deciding strategy, and no consensus exists on the best strategy between "single test", "dual testing", "sequential testing", "multiple testing options" and "no testing". In addition, there is a need for prognostic tools such as biological markers of recurrence and theragnostic tools to detect women not responding to chemotherapy in order to adapt therapies. Non-coding RNAs are classified as small or long based on their nucleotide count. Non-coding RNAs have multiple biological functions such as a role in tumorigenesis, gene regulation and genome protection. These ncRNAs emerge as new potential tools to differentiate benign from malignant tumors and to evaluate prognostic and theragnostic factors. In the specific setting of ovarian tumors, the goal of the present work is to offer an insight into the contribution of biofluid non-coding RNAs (ncRNA) expression.

Serum Extracellular Vesicle-Derived microRNAs as Potential Biomarkers for Pleural Mesothelioma in a European Prospective Study

Malignant pleural mesothelioma (MPM) is an aggressive cancer with a dismal prognosis. Early therapeutic interventions could improve patient outcomes. We aimed to identify a pattern of microRNAs (miRNAs) as potential early non-invasive markers of MPM. In a case-control study nested in the European Prospective Investigation into Cancer and Nutrition cohort, we screened the whole miRNome in serum extracellular vesicles (EVs) of preclinical MPM cases. In a subgroup of 20 preclinical samples collected five years prior MPM diagnosis, we observed an upregulation of miR-11400 (fold change (FC) = 2.6, adjusted p-value = 0.01), miR-148a-3p (FC = 1.5, p-value = 0.001), and miR-409-3p (FC = 1.5, p-value = 0.04) relative to matched controls. The 3-miRNA panel showed a good classification capacity with an area under the receiver operating characteristic curve (AUC) of 0.81 (specificity = 0.75, sensitivity = 0.70). The diagnostic ability of the model was also evaluated in an independent retrospective cohort, yielding a higher predictive power (AUC = 0.86). A signature of EV miRNA can be detected up to five years before MPM; moreover, the identified miRNAs could provide functional insights into the molecular changes related to the late carcinogenic process, preceding MPM development.

A systematic review of EV-piRNA in human body fluid and its role in disease progression

The state of host cells is reflected in the cargo carried by their extracellular vesicles (EVs). This makes EV a potential source of biomarkers for human diseases. Piwi-interacting RNA (piRNA) regulates gene expression through epigenetic regulation and post-transcriptional gene silencing. Thus, piRNA profiles in EVs derived from human clinical samples could identify markers that characterize disease stages, and unveil their roles in disease pathology. This review aimed to report the expression profiles of EV-derived piRNA (EV-piRNA) in various human samples, as well as their role in each pathology. A systematic review was conducted to collate the findings of human EV-piRNA from original research articles published in indexed scientific journals up to 16th Feb 2022. Article searches were performed in PubMed, Web of Science, and Scopus databases, using a combination of keywords including 'EV' and 'piRNA'. A total of 775 non-redundant original articles were identified. After subjecting articles to inclusion and exclusion criteria, 34 articles were accepted for this review. The piRNA expression levels among the small RNA profiles of human-derived EVs range from 0.09% to 43.84%, with the lowest expression level reported in urine-derived EVs and the highest percentage in plasma-derived EVs. Differentially expressed EV-piRNAs have been identified in patients with specific disease conditions compared to their counterparts (healthy control), suggesting an association between piRNA and progression in various diseases. Seven articles identified piRNA putative target genes and/or the pathway enrichment of piRNA target genes, and one study demonstrated a direct tole of piRNA candidates in disease pathology. In conclusion, EV-piRNA has been isolated successfully from various human body fluids. EV-piRNA is a new research niche in human disease pathology. The expression profiles of EV-piRNA in various tissue types and disease conditions remain largely unexplored. Furthermore, there is currently a lack of guidelines on piRNA bioinformatics analysis, which could lead to inconsistent results and thus hinder the progression of piRNA discoveries. Lastly, the lack of published scientific evidence on the role of EV-piRNA supports the need for future research to focus on the functional analysis of EV-piRNA as part of the route in piRNA discoveries.

Stool microRNA profiles reflect different dietary and gut microbiome patterns in healthy individuals

OBJECTIVES

MicroRNA (miRNA) profiles have been evaluated in several biospecimens in relation to common diseases for which diet may have a considerable impact. We aimed at characterising how specific diets are associated with the miRNome in stool of vegans, vegetarians and omnivores and how this is reflected in the gut microbial composition, as this is still poorly explored.

DESIGN

We performed small RNA and shotgun metagenomic sequencing in faecal samples and dietary recording from 120 healthy volunteers, equally distributed for the different diets and matched for sex and age.

RESULTS

We found 49 miRNAs differentially expressed among vegans, vegetarians and omnivores (adj. p <0.05) and confirmed trends of expression levels of such miRNAs in vegans and vegetarians compared with an independent cohort of 45 omnivores. Two miRNAs related to lipid metabolism, miR-636 and miR-4739, were inversely correlated to the non-omnivorous diet duration, independently of subject age. Seventeen miRNAs correlated (|rho|>0.22, adj. p <0.05) with the estimated intake of nutrients, particularly animal proteins, phosphorus and, interestingly, lipids. In omnivores, higher Prevotella and Roseburia and lower Bacteroides abundances than in vegans and vegetarians were observed. Lipid metabolism-related miR-425-3p and miR-638 expression levels were associated with increased abundances of microbial species, such as Roseburia sp. CAG 182 and Akkermansia muciniphila, specific of different diets. An integrated analysis identified 25 miRNAs, 25 taxa and 7 dietary nutrients that clearly discriminated (area under the receiver operating characteristic curve=0.89) the three diets.

CONCLUSION

Stool miRNA profiles are associated with specific diets and support the role of lipids as a driver of epigenetic changes and host-microbial molecular interactions in the gut.

Small and intermediate size structural RNAs in the unicellular parasite Cryptosporidium parvum as revealed by sRNA-seq and comparative genomics

Small and intermediate-size noncoding RNAs (sRNAs and is-ncRNAs) have been shown to play important regulatory roles in the development of several eukaryotic organisms. However, they have not been thoroughly explored in Cryptosporidium parvum, an obligate zoonotic protist parasite responsible for the diarrhoeal disease cryptosporidiosis. Using Illumina sequencing of a small RNA library, a systematic identification of novel small and is-ncRNAs was performed in C. parvum excysted sporozoites. A total of 79 novel is-ncRNA candidates, including antisense, intergenic and intronic is-ncRNAs, were identified, including 7 new small nucleolar RNAs (snoRNAs). Expression of select novel is-ncRNAs was confirmed by RT-PCR. Phylogenetic conservation was analysed using covariance models (CMs) in related Cryptosporidium and apicomplexan parasite genome sequences. A potential new type of small ncRNA derived from tRNA fragments was observed. Overall, a deep profiling analysis of novel is-ncRNAs in C. parvum and related species revealed structural features and conservation of these novel is-ncRNAs. Covariance models can be used to detect is-ncRNA genes in other closely related parasites. These findings provide important new sequences for additional functional characterization of novel is-ncRNAs in the protist pathogen C. parvum.

Combined miRNA and SERS urine liquid biopsy for the point-of-care diagnosis and molecular stratification of bladder cancer

BACKGROUND

Bladder cancer (BC) has the highest per-patient cost of all cancer types. Hence, we aim to develop a non-invasive, point-of-care tool for the diagnostic and molecular stratification of patients with BC based on combined microRNAs (miRNAs) and surface-enhanced Raman spectroscopy (SERS) profiling of urine.

METHODS

Next-generation sequencing of the whole miRNome and SERS profiling were performed on urine samples collected from 15 patients with BC and 16 control subjects (CTRLs). A retrospective cohort (BC = 66 and CTRL = 50) and RT-qPCR were used to confirm the selected differently expressed miRNAs. Diagnostic accuracy was assessed using machine learning algorithms (logistic regression, naïve Bayes, and random forest), which were trained to discriminate between BC and CTRL, using as input either miRNAs, SERS, or both. The molecular stratification of BC based on miRNA and SERS profiling was performed to discriminate between high-grade and low-grade tumors and between luminal and basal types.

RESULTS

Combining SERS data with three differentially expressed miRNAs (miR-34a-5p, miR-205-3p, miR-210-3p) yielded an Area Under the Curve (AUC) of 0.92 ± 0.06 in discriminating between BC and CTRL, an accuracy which was superior either to miRNAs (AUC = 0.84 ± 0.03) or SERS data (AUC = 0.84 ± 0.05) individually. When evaluating the classification accuracy for luminal and basal BC, the combination of miRNAs and SERS profiling averaged an AUC of 0.95 ± 0.03 across the three machine learning algorithms, again better than miRNA (AUC = 0.89 ± 0.04) or SERS (AUC = 0.92 ± 0.05) individually, although SERS alone performed better in terms of classification accuracy.

CONCLUSION

miRNA profiling synergizes with SERS profiling for point-of-care diagnostic and molecular stratification of BC. By combining the two liquid biopsy methods, a clinically relevant tool that can aid BC patients is envisaged.

The Clinical Significance of Transfer RNAs Present in Extracellular Vesicles

Extracellular vesicles (EVs) are important for intercellular signalling in multi-cellular organisms. However, the role of mature transfer RNAs (tRNAs) and tRNA fragments in EVs has yet to be characterised. This systematic review aimed to identify up-to-date literature on tRNAs present within human EVs and explores their potential clinical significance in health and disease. A comprehensive and systematic literature search was performed, and the study was conducted in accordance with PRISMA guidelines. Electronic databases MEDLINE and EMBASE were searched up until 1 January 2022. From 685 papers, 60 studies were identified for analysis. The majority of papers reviewed focussed on the role of EV tRNAs in cancers (31.7%), with numerous other conditions represented. Blood and cell lines were the most common EV sources, representing 85.9% of protocols used. EV isolation methods included most known methods, precipitation being the most common (49.3%). The proportion of EV tRNAs was highly variable, ranging between 0.04% to >95% depending on tissue source. EV tRNAs are present in a multitude of sources and show promise as disease markers in breast cancer, gastrointestinal cancers, and other diseases. EV tRNA research is an emerging field, with increasing numbers of papers highlighting novel methodologies for tRNA and tRNA fragment discovery.

A Novel Tissue Atlas and Online Tool for the Interrogation of Small RNA Expression in Human Tissues and Biofluids

One promising goal for utilizing the molecular information circulating in biofluids is the discovery of clinically useful biomarkers. Extracellular RNAs (exRNAs) are one of the most diverse classes of molecular cargo, easily assayed by sequencing and with expressions that rapidly change in response to subject status. Despite diverse exRNA cargo, most evaluations from biofluids have focused on small RNA sequencing and analysis, specifically on microRNAs (miRNAs). Another goal of characterizing circulating molecular information, is to correlate expression to injuries associated with specific tissues of origin. Biomarker candidates are often described as being specific, enriched in a particular tissue or associated with a disease process. Likewise, miRNA data is often reported to be specific, enriched for a tissue, without rigorous testing to support the claim. Here we provide a tissue atlas of small RNAs from 30 different tissues and three different blood cell types. We analyzed the tissues for enrichment of small RNA sequences and assessed their expression in biofluids: plasma, cerebrospinal fluid, urine, and saliva. We employed published data sets representing physiological (resting vs. acute exercise) and pathologic states (early- vs. late-stage liver fibrosis, and differential subtypes of stroke) to determine differential tissue-enriched small RNAs. We also developed an online tool that provides information about exRNA sequences found in different biofluids and tissues. The data can be used to better understand the various types of small RNA sequences in different tissues as well as their potential release into biofluids, which should help in the validation or design of biomarker studies.

Faecal miRNA profiles associated with age, sex, BMI, and lifestyle habits in healthy individuals

For their stability and detectability faecal microRNAs represent promising molecules with potential clinical interest as non-invasive diagnostic and prognostic biomarkers. However, there is no evidence on how stool miRNA profiles change according to an individual’s age, sex, and body mass index (BMI) or how lifestyle habits influence the expression levels of these molecules. We explored the relationship between the stool miRNA levels and common traits (sex, age, BMI, and menopausal status) or lifestyle habits (physical activity, smoking status, coffee, and alcohol consumption) as derived by a self-reported questionnaire, using small RNA-sequencing data of samples from 335 healthy subjects. We detected 151 differentially expressed miRNAs associated with one variable and 52 associated with at least two. Differences in miR-638 levels were associated with age, sex, BMI, and smoking status. The highest number of differentially expressed miRNAs was associated with BMI (n = 92) and smoking status (n = 84), with several miRNAs shared between them. Functional enrichment analyses revealed the involvement of the miRNA target genes in pathways coherent with the analysed variables. Our findings suggest that miRNA profiles in stool may reflect common traits and lifestyle habits and should be considered in relation to disease and association studies based on faecal miRNA expression.

Extracellular microRNA 3' end modification across diverse body fluids

microRNAs (miRNAs) are small non-coding RNAs that play critical roles in gene regulation. The presence of miRNAs in extracellular biofluids is increasingly recognized. However, most previous characterization of extracellular miRNAs focused on their overall expression levels. Alternative sequence isoforms and modifications of miRNAs were rarely considered in the extracellular space. Here, we developed a highly accurate bioinformatic method, called miNTA, to identify 3' non-templated additions (NTAs) of miRNAs using small RNA-sequencing data. Using miNTA, we conducted an in-depth analysis of miRNA 3' NTA profiles in 1047 extracellular RNA-sequencing data sets of 4 types of biofluids. This analysis identified hundreds of miRNAs with 3' uridylation or adenylation, with the former being more prevalent. Among these miRNAs, up to 53% (22%) had an average 3' uridylation (adenylation) level of at least 10% in a specific biofluid. Strikingly, we found that 3' uridylation levels enabled segregation of different types of biofluids, more effectively than overall miRNA expression levels. This observation suggests that 3' NTA levels possess fluid-specific information relatively robust to batch effects. In addition, we observed that extracellular miRNAs with 3' uridylations are enriched in processes related to angiogenesis, apoptosis, and inflammatory response, and this type of modification may stabilize base-pairing between miRNAs and their target genes. Together, our study provides a comprehensive landscape of miRNA NTAs in human biofluids, which paves way for further biomarker discoveries. The insights generated in our work built a foundation for future functional, mechanistic, and translational discoveries.

The miRNA content of circulating exosomes in DLBCL patients and in vitro influence of DLBCL-derived exosomes on miRNA expression of healthy B cells from peripheral blood

BACKGROUND

Due to the heterogeneous nature of Diffuse Large B-cell Lymphoma (DLBCL), the mechanisms underlying tumor development and progression have not yet been fully elucidated.

OBJECTIVE

This study aimed to compare the characteristics of plasma exosomes of DLBCL patients and healthy individuals and to evaluate the exosomal interactions between DLBCL cell lines and normal B-cells.

METHODS

Exosome isolation was performed using an ultracentrifugation-based protocol from plasma of 20 patients with DLBCL and 20 controls. The expression of miRNAs from exosome samples was analyzed using a miRNA expression microarray. The presence of exosome-mediated communication between the lymphoma cells and normal B-cells was determined by the co-culture model.

RESULTS

A significant increase in plasma exosome concentrations of DLBCL patients was observed. There was also a significant decrease in the expression of 33 miRNAs in plasma exosomes of DLBCL patients. It was determined that normal B-cells internalize DLBCL-derived exosomes and then miRNA expression differences observed in normal B-cells are specific to lymphoma-subtypes.

CONCLUSIONS

MiR-3960, miR-6089 and miR-939-5p can be used as the miRNA signature in DLBCL diagnosis. We suppose that the exosomes changed the molecular signature of the target cells depending on the genomic characterization of the lymphoma cells they have originated.

Identification of novel circulating microRNAs in advanced heart failure by next-generation sequencing

Aims

Risk stratification in patients with advanced chronic heart failure (HF) is an unmet need. Circulating microRNA (miRNA) levels have been proposed as diagnostic and prognostic biomarkers in several diseases including HF. The aims of the present study were to characterize HF‐specific miRNA expression profiles and to identify miRNAs with prognostic value in HF patients.

Methods and results

We performed a global miRNome analysis using next‐generation sequencing in the plasma of 30 advanced chronic HF patients and of matched healthy controls. A small subset of miRNAs was validated by real‐time PCR (P < 0.0008). Pearson's correlation analysis was computed between miRNA expression levels and common HF markers. Multivariate prediction models were exploited to evaluate miRNA profiles' prognostic role. Thirty‐two miRNAs were found to be dysregulated between the two groups. Six miRNAs (miR‐210‐3p, miR‐22‐5p, miR‐22‐3p, miR‐21‐3p, miR‐339‐3p, and miR‐125a‐5p) significantly correlated with HF biomarkers, among which N‐terminal prohormone of brain natriuretic peptide. Inside the cohort of advanced HF population, we identified three miRNAs (miR‐125a‐5p, miR‐10b‐5p, and miR‐9‐5p) altered in HF patients experiencing the primary endpoint of cardiac death, heart transplantation, or mechanical circulatory support implantation when compared with those without clinical events. The three miRNAs added substantial prognostic power to Barcelona Bio‐HF score, a multiparametric and validated risk stratification tool for HF (from area under the curve = 0.72 to area under the curve = 0.82).

Conclusions

This discovery study has characterized, for the first time, the advanced chronic HF‐specific miRNA expression pattern. We identified a few miRNAs able to improve the prognostic stratification of HF patients based on common clinical and laboratory values. Further studies are needed to validate our results in larger populations.

Adipocyte-derived extracellular vesicles regulate survival and function of pancreatic β cells

Extracellular vesicles (EVs) are implicated in the crosstalk between adipocytes and other metabolic organs, and an altered biological cargo has been observed in EVs from human obese adipose tissue (AT). Yet, the role of adipocyte-derived EVs in pancreatic β cells remains to be determined. Here, we explored the effects of EVs released from adipocytes isolated from both rodents and humans and human AT explants on survival and function of pancreatic β cells and human pancreatic islets. EVs from healthy 3T3-L1 adipocytes increased survival and proliferation and promoted insulin secretion in INS-1E β cells and human pancreatic islets, both those untreated or exposed to cytokines or glucolipotoxicity, whereas EVs from inflamed adipocytes caused β cell death and dysfunction. Human lean adipocyte-derived EVs produced similar beneficial effects, whereas EVs from obese AT explants were harmful for human EndoC-βH3 β cells. We observed differential expression of miRNAs in EVs from healthy and inflamed adipocytes, as well as alteration in signaling pathways and expression of β cell genes, adipokines, and cytokines in recipient β cells. These in vitro results suggest that, depending on the physiopathological state of AT, adipocyte-derived EVs may influence β cell fate and function.

Exosomal microRNAs and other non-coding RNAs as colorectal cancer biomarkers: a review

The circulating human transcriptome, which includes both coding and non-coding RNA (ncRNA) molecules, represents a rich source of potential biomarkers for colorectal cancer (CRC) that has only recently been explored. In particular, the release of RNA-containing extracellular vesicles (EVs), in a multitude of different in vitro cell systems and in a variety of body fluids, has attracted wide interest. The role of RNA species in EVs is still not fully understood, but their capacity to act as a form of distant communication between cells and their higher abundance in association with cancer demonstrated their relevance. In this review, we report the evidence from both in vitro and human studies on microRNAs (miRNAs) and other ncRNA profiles analysed in EVs in relation to CRC as diagnostic, prognostic and predictive markers. The studies so far highlighted that, in exosomes, the most studied category of EVs, several miRNAs are able to accurately discriminate CRC cases from controls as well as to describe the progression of the disease and its prognosis. Most of the time, the in vitro findings support the miRNA profiles detected in human exosomes. The expression profiles measured in exosomes and other EVs differ and, interestingly, there is a variability of expression also among different subsets of exosomes according to their proteic profile. On the other hand, evidence is still limited for what concerns exosome miRNAs as early diagnostic and predictive markers of treatment. Several other ncRNAs that are carried by exosomes, mostly long ncRNAs and circular RNAs, seem also to be dysregulated in CRC. Besides various technical challenges, such as the standardisation of EVs isolation methods and the optimisation of methodologies to characterise the whole spectrum of RNA molecules in exosomes, further studies are needed in order to elucidate their relevance as CRC markers.

Intake of Natural Compounds and Circulating microRNA Expression Levels: Their Relationship Investigated in Healthy Subjects With Different Dietary Habits

Diet has a strong influence on many physiological processes, which in turn have important implications on a variety of pathological conditions. In this respect, microRNAs (miRNAs), a class of small non-coding RNAs playing a relevant epigenetic role in controlling gene expression, may represent mediators between the dietary intake and the healthy status. Despite great advances in the field of nutri-epigenomics, it remains unclear how miRNA expression is modulated by the diet and, specifically, the intake of specific nutrients. We investigated the whole circulating miRNome by small RNA-sequencing performed on plasma samples of 120 healthy volunteers with different dietary habits (vegans, vegetarians, and omnivores). Dietary intakes of specific nutrients were estimated for each subject from the information reported in the food-frequency questionnaire previously validated in the EPIC study. We focused hereby on the intake of 23 natural compounds (NCs) of the classes of lipids, micro-elements, and vitamins. We identified 78 significant correlations (rho > 0.300, p-value < 0.05) among the estimated daily intake of 13 NCs and the expression levels of 58 plasma miRNAs. Overall, vitamin D, sodium, and vitamin E correlated with the largest number of miRNAs. All the identified correlations were consistent among the three dietary groups and 22 of them were confirmed as significant (p-value < 0.05) by age-, gender-, and body-mass index-adjusted Generalized Linear regression Model analysis. miR-23a-3p expression levels were related with different NCs including a significant positive correlation with sodium (rho = 0.377) and significant negative correlations with lipid-related NCs and vitamin E. Conversely, the estimated intake of vitamin D was negatively correlated with the expression of the highest number of circulating miRNAs, particularly miR-1277-5p (rho = −0.393) and miR-144-3p (rho = −0.393). Functional analysis of the targets of sodium intake-correlated miRNAs highlighted terms related to cardiac development. A similar approach on targets of those miRNAs correlated with vitamin D intake showed an enrichment in genes involved in hormone metabolisms, while the response to chronic inflammation was among the top enriched processes involving targets of miRNAs negatively related with vitamin E intake. Our findings show that nutrients through the habitual diet influence circulating miRNA profiles and highlight that this aspect must be considered in the nutri-epigenomic research.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas

Extracellular RNAs present in biofluids have emerged as potential biomarkers for disease. Where most studies focus on blood-derived fluids, other biofluids may be more informative. We present an atlas of messenger, circular, and small RNA transcriptomes of a comprehensive collection of 20 human biofluids. By means of synthetic spike-in controls, we compare RNA content across biofluids, revealing a 10,000-fold difference in concentration. The circular RNA fraction is increased in most biofluids compared to tissues. Each biofluid transcriptome is enriched for RNA molecules derived from specific tissues and cell types. Our atlas enables an informed selection of the most relevant biofluid to monitor particular diseases. To verify the biomarker potential in these biofluids, four validation cohorts representing a broad spectrum of diseases were profiled, revealing numerous differential RNAs between case and control subjects. Spike-normalized data are publicly available in the R2 web portal for further exploration.

The Inhibitory Role of miR-486-5p on CSC Phenotype Has Diagnostic and Prognostic Potential in Colorectal Cancer

Colorectal cancer (CRC) is the third most frequent cancer worldwide and the second cause of cancer deaths. Increasing evidences supports the idea that the poor prognosis of patients is related to the presence of cancer stem cells (CSCs), a cell population able to drive cancer recurrence and metastasis. The deregulation of microRNAs (miRNAs) plays a role in the formation of CSC. We investigated the role of hsa-miR-486-5p (miR-486-5p) in CRC, CSCs, and metastasis, in order to reach a better understanding of the biomolecular and epigenetic mechanisms mir-486-5p-related. The expression of miR-486-5p was investigated in three different matrices from CRC patients and controls and in CSCs obtained from the CRC cell lines HCT-116, HT-29, and T-84. In the human study, miR-486-5p was up-regulated in serum and stool of CRC patients in comparison with healthy controls but down-regulated in tumor tissue when compared with normal mucosa. miR-486-5p was also down-regulated in the sera of metastatic patients. In vitro, miR-486-5p was down-regulated in CSC models and it induced an inhibitory effect on stem factors and oncogenes in the main pathways of CSCs. Our results provide a step forward in understanding the role of mir-486-5p in CRC and CSC, and suggest that further studies are needed to investigate its diagnostic and prognostic power, possibly in combination with other biomarkers.

Non-coding RNA biomarkers in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, which is usually diagnosed at an advanced stage. The late disease diagnosis, the limited availability of effective therapeutic interventions and lack of robust diagnostic biomarkers, are some of the primary reasons for the dismal 5-year survival rates (∼8%) in patients with PDAC. The pancreatic cancer develops through accumulation of a series of genomic and epigenomic alterations which lead to the transformation of normal pancreatic epithelium into an invasive carcinoma - a process that can take up to 15-20 years to develop, from the occurrence of first initiating mutational event. These facts highlight a unique window of opportunity for the earlier detection of PDAC, which could allow timely disease interception and improvement in the overall survival outcomes in patients suffering from this fatal malignancy. Non-coding RNAs (ncRNAs) have been recognized to play a central role in PDAC pathogenesis and are emerging as attractive candidates for biomarker development in various cancers, including PDAC. More specifically, the ncRNAs play a pivotal role in PDAC biology as they affect tumor growth, migration, and invasion by regulating cellular processes including cell cycle, apoptosis, and epithelial-mesenchymal transition. In this review, we focus on three types of well-established ncRNAs - microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) - and discuss their potential as diagnostic, prognostic and predictive biomarkers in PDAC.

Mast Cells, microRNAs and Others: The Role of Translational Research on Colorectal Cancer in the Forthcoming Era of Precision Medicine

Colorectal cancer (CRC) is a heterogeneous disease, molecularly and anatomically, that develops in a multi-step process requiring the accumulation of several genetic or epigenetic mutations that lead to the gradual transformation of normal mucosa into cancer. In fact, tumorigenesis is extremely complex, with many immunologic and non-immunologic factors present in the tumor microenvironment that can influence tumorigenesis. In the last few years, a role for mast cells (MCs), microRNAs (miRNAs), Kirsten rat sarcoma (KRAS) and v-raf murine sarcoma viral oncogene homologue B (BRAF) in cancer development and progression has been suggested, and numerous efforts have been made to thoroughly assess their correlation with CRC to improve patient survival and quality of life. The identification of easily measurable, non-invasive and cost-effective biomarkers, the so-called "ideal biomarkers", for CRC screening and treatment remains a high priority. The aim of this review is to discuss the emerging role of mast cells (MCs), microRNAs (miRNAs), KRAS and BRAF as diagnostic and prognostic biomarkers for CRC, evaluating their influence as potential therapy targets in the forthcoming era of precision medicine.

Circulating MicroRNAs from Serum Exosomes May Serve as a Putative Biomarker in the Diagnosis and Treatment of Patients with Focal Cortical Dysplasia

Focal cortical dysplasia (FCD) is a congenital malformation of cortical development where the cortical neurons located in the brain area fail to migrate in the proper formation. Epilepsy, particularly medically refractory epilepsy, is the most common clinical presentation for all types of FCD. This study aimed to explore the expression change of circulating miRNAs in patients with FCD from serum exosomes. A total of nine patients with FCD and four healthy volunteers were enrolled in this study. The serum exosomes were isolated from the peripheral blood of the subjects. Transmission electron microscopy (TEM) was used to identify the exosomes. Both exosomal markers and neuronal markers were detected by Western blotting analysis to prove that we could obtain central nervous system-derived exosomes from the circulation. The expression profiles of circulating exosomal miRNAs were assessed using next-generation sequencing analysis (NGS). We obtained a total of 107 miRNAs with dominant fold change (>2-fold) from both the annotated 5p-arm and 3p-arm of 2780 mature miRNAs. Based on the integrated platform of HMDD v3.2, miRway DB and DIANA-miRPath v3.0 online tools, and confirmed by MiRBase analysis, four potentially predicted miRNAs from serum exosomes in patients with FCD were identified, including miR194-2-5p, miR15a-5p, miR-132-3p, and miR-145-5p. All four miRNAs presented upregulated expression in patients with FCD compared with controls. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and pathway category of four target miRNAs, we found eight possible signaling pathways that may be related to FCD. Among them, we suggest that the mTOR signaling pathway, PI3K-Akt signaling pathway, p53 signaling pathway, and cell cycle regulation and TGF-beta signaling pathway are high-risk pathways that play a crucial role in the pathogenesis of FCD and refractory epilepsy. Our results suggest that the circulating miRNAs from exosomes may provide a potential biomarker for diagnostic, prognostic, and therapeutic adjuncts in patients with FCD and refractory epilepsy.

Small Non-Coding RNA Profiling in Plasma Extracellular Vesicles of Bladder Cancer Patients by Next-Generation Sequencing: Expression Levels of miR-126-3p and piR-5936 Increase with Higher Histologic Grades

Bladder cancer (BC) is the tenth most frequent cancer worldwide. Due to the need for recurrent cystoscopies and the lack of non-invasive biomarkers, BC is associated with a high management burden. In this respect, small non-coding RNAs (sncRNAs) have been investigated in urine as possible biomarkers for BC, but in plasma their potential has not yet been defined. The expression levels of sncRNAs contained in plasma extracellular vesicles (EVs) from 47 men with BC and 46 healthy controls were assessed by next-generation sequencing. The sncRNA profiles were compared with urinary profiles from the same subjects. miR-4508 resulted downregulated in plasma EVs of muscle-invasive BC patients, compared to controls (adj-p = 0.04). In World Health Organization (WHO) grade 3 (G3) BC, miR-126-3p was upregulated both in plasma EVs and urine, when compared to controls (for both, adj-p < 0.05). Interestingly, two sncRNAs were associated with the risk class: miR-4508 with a downward trend going from controls to high risk BC, and piR-hsa-5936 with an upward trend (adj-p = 0.04 and adj-p = 0.05, respectively). Additionally, BC cases with low expression of miR-185-5p and miR-106a-5p or high expression of miR-10b-5p showed shorter survival (adj-p = 0.0013, adj-p = 0.039 and adj-p = 0.047, respectively). SncRNAs from plasma EVs could be diagnostic biomarkers for BC, especially in advanced grade.

MIR1246 in body fluids as a biomarker for pancreatic cancer

Pancreatic cancer is an aggressive tumor associated with poor survival, and early detection is important to improve patient outcomes. In the present study, we examined MIR1246 expression as a biomarker of pancreatic cancer. Total RNA was extracted from serum, urine and saliva samples from healthy subjects (n = 30) and patients with pancreatic cancer (n = 41, stage 0–IV). The MIR1246 level in each fluid was analyzed by quantitative reverse transcription-polymerase chain reaction. Significantly higher MIR1246 expression in serum and urine was observed in patients with cancer than in healthy controls. A significant positive correlation was found between serum and urine MIR1246 expression (r = 0.34). Receiver operating characteristic curves were constructed for MIR1246 in all three body fluids. The area under the curve for serum MIR1246 was 0.87 (sensitivity, 92.3%; specificity, 73.3%), and that for urine MIR1246 was 0.90 (sensitivity, 90.2%; specificity, 83.3%). With a cut-off of the control group’s mean plus twice the standard deviation, the sensitivities of MIR1246 in serum and urine for pancreatic cancer were 60.9 and 58.5%, respectively. Combining both serum and urine MIR1246 expression yielded a sensitivity of 85%. These results indicate that MIR246 may be a useful diagnostic biomarker for pancreatic cancer.

Aberrant expression of miR-21 in patients with inflammatory bowel disease: A protocol for systematic review and meta analysis

BACKGROUND

microRNAs have drawn more attention due to their function on the inflammatory process. The association between microRNA-21 (miR-21) expression and risk of inflammatory bowel diseases (IBD) remain inconclusive. This study was aimed to acquire a more exact estimation of this relationship.

METHODS

Relevant studies were identified through searching PubMed, Embase, Wanfang, and China National Knowledge Infrastructure database. Pooled standardized mean difference and 95% confidence intervals were calculated using a random-effect model. Publication bias test, sensitivity analysis and subgroup analysis were carried out.

RESULTS

A total of 20 relevant articles comprising 540 patients with ulcerative colitis (UC), 459 patients with Crohn disease (CD) and 511 non-IBD controls were included in this analysis. The expression of miR-21 was significantly increased in colon tissue of both UC and CD patients compared with non-IBD controls. However, there were no significant differences between patients with UC and CD. Moreover, increased miR-21 expression was associated with disease activity status in UC patients, but not in CD patients.

CONCLUSIONS

This meta-analysis demonstrates that the higher miR-21 expression in colon tissue is positively associated with the development of UC and CD, and miR-21 might serve as a disease marker of IBD.

The complex interplay between DNA methylation and miRNAs in gene expression regulation

The short, non-coding RNAs, also called microRNAs (miRNAs) can bind complementary sequences on cellular mRNAs. The consequence of this binding is generally the degradation of mRNA and the inhibition of its translation. For this reason, miRNAs are included among the epigenetic factors acting as a modulator of gene expression. How miRNAs expression is, in turn, regulated is still the object of active investigation, but DNA methylation, another epigenetic modification, seems to play a central role in this sense. The "one-carbon" metabolism is responsible for the metabolic regulation of trans-methylation reactions and, therefore, DNA methylation. For this reason, to investigate the possible correlations between alterations of the one-carbon metabolism and differential DNA methylation sounds interesting. Moreover, recent evidence indicates that, vice-versa, miRNAs are associated with DNA methylation modulation, in a mutual cross-talk. The present review will discuss the interplay between miRNAs and DNA methylation and its fall-out on gene expression regulation.

tRNA-derived fragments and microRNAs in the maternal-fetal interface of a mouse maternal-immune-activation autism model

tRNA-derived small fragments (tRFs) and tRNA halves have emerging functions in different biological pathways, such as regulating gene expression, protein translation, retrotransposon activity, transgenerational epigenetic changes and response to environmental stress. However, small RNAs like tRFs and microRNAs in the maternal-fetal interface during gestation have not been studied extensively. Here we investigated the small RNA composition of mouse placenta/decidua, which represents the interface where the mother communicates with the foetus, to determine whether there are specific differences in tRFs and microRNAs during fetal development and in response to maternal immune activation (MIA). Global tRF expression pattern, just like microRNAs, can distinguish tissue types among placenta/decidua, fetal brain and fetal liver. In particular, 5' tRNA halves from tRNA^Gly, tRNA^Glu, tRNA^Val and tRNA^Lys are abundantly expressed in the normal mouse placenta/decidua. Moreover, tRF and microRNA levels in the maternal-fetal interface change dynamically over the course of embryonic development. To see if stress alters non-coding RNA expression at the maternal-fetal interface, we treated pregnant mice with a viral infection mimetic, which has been shown to promote autism-related phenotypes in the offspring. Acute changes in the levels of specific tRFs and microRNAs were observed 3-6 h after MIA and are suppressed thereafter. A group of 5' tRNA halves is down-regulated by MIA, whereas a group of 18-nucleotide tRF-3a is up-regulated. In conclusion, tRFs show tissue-specificity, developmental changes and acute response to environmental stress, opening the possibility of them having a role in the fetal response to MIA.

Angiogenin generates specific stress-induced tRNA halves and is not involved in tRF-3-mediated gene silencing

tRNA fragments (tRFs) and tRNA halves have been implicated in various cellular processes, including gene silencing, translation, stress granule assembly, cell differentiation, retrotransposon activity, symbiosis, apoptosis, and more. Overexpressed angiogenin (ANG) cleaves tRNA anticodons and produces tRNA halves similar to those produced in response to stress. However, it is not clear whether endogenous ANG is essential for producing the stress-induced tRNA halves. It is also not clear whether smaller tRFs are generated from the tRNA halves. Here, using global short RNA-Seq approach, we found that ANG overexpression selectively cleaves a subset of tRNAs, including tRNA^Glu, tRNA^Gly, tRNA^Lys, tRNA^Val, tRNA^His, tRNA^Asp, and tRNA^SeC to produce tRNA halves and tRF-5s that are 26-30 bases long. Surprisingly, ANG knockout revealed that the majority of stress-induced tRNA halves, except for the 5' half from tRNA^HisGTG and the 3' half from tRNA^AspGTC, are ANG independent, suggesting there are other RNases that produce tRNA halves. We also found that the 17-25 bases-long tRF-3s and tRF-5s that could enter into Argonaute complexes are not induced by ANG overexpression, suggesting that they are generated independently from tRNA halves. Consistent with this, ANG knockout did not decrease tRF-3 levels or gene-silencing activity. We conclude that ANG cleaves specific tRNAs and is not the only RNase that creates tRNA halves and that the shorter tRFs are not generated from the tRNA halves or from independent tRNA cleavage by ANG.

Altered Fecal Small RNA Profiles in Colorectal Cancer Reflect Gut Microbiome Composition in Stool Samples

The characteristics of microbial small RNA transcription are largely unknown, while it is of primary importance for a better identification of molecules with functional activities in the gut niche under both healthy and disease conditions. By performing combined analyses of metagenomic and small RNA sequencing (sRNA-Seq) data, we characterized both the human and microbial small RNA contents of stool samples from healthy individuals and from patients with colorectal carcinoma or adenoma. With the integrative analyses of metagenomic and sRNA-Seq data, we identified a human and microbial small RNA signature which can be used to improve diagnosis of the disease. Our analysis of human and gut microbiome small RNA expression is relevant to generation of the first hypotheses about the potential molecular interactions occurring in the gut of CRC patients, and it can be the basis for further mechanistic studies and clinical tests.

Dysbiotic configurations of the human gut microbiota have been linked to colorectal cancer (CRC). Human small noncoding RNAs are also implicated in CRC, and recent findings suggest that their release in the gut lumen contributes to shape the gut microbiota. Bacterial small RNAs (bsRNAs) may also play a role in carcinogenesis, but their role has been less extensively explored. Here, we performed small RNA and shotgun sequencing on 80 stool specimens from patients with CRC or with adenomas and from healthy subjects collected in a cross-sectional study to evaluate their combined use as a predictive tool for disease detection. We observed considerable overlap and a correlation between metagenomic and bsRNA quantitative taxonomic profiles obtained from the two approaches. We identified a combined predictive signature composed of 32 features from human and microbial small RNAs and DNA-based microbiome able to accurately classify CRC samples separately from healthy and adenoma samples (area under the curve [AUC] = 0.87). In the present study, we report evidence that host-microbiome dysbiosis in CRC can also be observed by examination of altered small RNA stool profiles. Integrated analyses of the microbiome and small RNAs in the human stool may provide insights for designing more-accurate tools for diagnostic purposes.

IMPORTANCE The characteristics of microbial small RNA transcription are largely unknown, while it is of primary importance for a better identification of molecules with functional activities in the gut niche under both healthy and disease conditions. By performing combined analyses of metagenomic and small RNA sequencing (sRNA-Seq) data, we characterized both the human and microbial small RNA contents of stool samples from healthy individuals and from patients with colorectal carcinoma or adenoma. With the integrative analyses of metagenomic and sRNA-Seq data, we identified a human and microbial small RNA signature which can be used to improve diagnosis of the disease. Our analysis of human and gut microbiome small RNA expression is relevant to generation of the first hypotheses about the potential molecular interactions occurring in the gut of CRC patients, and it can be the basis for further mechanistic studies and clinical tests.

Noncoding RNAs in Extracellular Fluids as Cancer Biomarkers: The New Frontier of Liquid Biopsies

The last two decades of cancer research have been devoted in two directions: (1) understanding the mechanism of carcinogenesis for an effective treatment, and (2) improving cancer prevention and screening for early detection of the disease. This last aspect has been developed, especially for certain types of cancers, thanks also to the introduction of new concepts such as liquid biopsies and precision medicine. In this context, there is a growing interest in the application of alternative and noninvasive methodologies to search for cancer biomarkers. The new frontiers of the research lead to a search for RNA molecules circulating in body fluids. Searching for biomarkers in extracellular body fluids represents a better option for patients because they are easier to access, less painful, and potentially more economical. Moreover, the possibility for these types of samples to be taken repeatedly, allows a better monitoring of the disease progression or treatment efficacy for a better intervention and dynamic treatment of the patient, which is the fundamental basis of personalized medicine. RNA molecules, freely circulating in body fluids or packed in microvesicles, have all the characteristics of the ideal biomarkers owing to their high stability under storage and handling conditions and being able to be sampled several times for monitoring. Moreover, as demonstrated for many cancers, their plasma/serum levels mirror those in the primary tumor. There are a large variety of RNA species noncoding for proteins that could be used as cancer biomarkers in liquid biopsies. Among them, the most studied are microRNAs, but recently the attention of the researcher has been also directed towards Piwi-interacting RNAs, circular RNAs, and other small noncoding RNAs. Another class of RNA species, the long noncoding RNAs, is larger than microRNAs and represents a very versatile and promising group of molecules which, apart from their use as biomarkers, have also a possible therapeutic role. In this review, we will give an overview of the most common noncoding RNA species detectable in extracellular fluids and will provide an update concerning the situation of the research on these molecules as cancer biomarkers.

Comparison of RNA isolation and library preparation methods for small RNA sequencing of canine biofluids

BACKGROUND

Small RNA sequencing (RNA-seq) of biofluids is challenging due to the relative scarcity of microRNAs (miRNAs), limited sample volumes, and the lack of a gold standard isolation method. Additionally, few comparisons exist for the RNA isolation and sequencing methods of biofluids.

OBJECTIVES

We aimed to compare the performance of six commercial RNA isolation kits and two library preparation methods for small RNA-seq using canine serum and urine.

METHODS

Serum and urine were collected from seven dogs with protein-losing nephropathy, and the samples were pooled. Total RNA from serum (2 mL) and urine (10 mL) was isolated in triplicate using three methods each for serum (Zymo Direct-zol, mirVana PARIS, miRCURY Biofluids) and urine (Qiagen exoRNeasy, Norgen Urine Exosome, miRCURY Exosome). For each sample type, the two kits yielding the highest RNA concentration were selected, and small RNA-seq was performed using TruSeq and NEXTflex library preparations. Data were analyzed by CPSS 2.0 and DESeq2.

RESULTS

For serum, Zymo Direct-zol combined with NEXTflex was the only combination that enabled successful library preparation, while for urine, Qiagen exoRNeasy combined with NEXTflex outperformed other combinations for detecting miRNAs. The total number of miRNAs detected in serum and urine was 198 and up to 115, respectively. miRNA expression in serum was distinct from urine. Furthermore, the library preparation method introduced a higher variation of urine results than the RNA isolation method.

CONCLUSIONS

Different isolation and library preparation methods show significant differences in miRNA results that could affect biomarker discovery. Small RNA-seq provides an unbiased, global assessment to compare these methods in canine biofluids.

MMP23B expression and protein levels in blood and urine are associated with bladder cancer

Urothelial bladder cancer (UBC) represents a public health problem because of its high incidence/relapse rates. At present, there are no suitable biomarkers for early diagnosis or relapse/progression prognosis. To improve diagnostic accuracy and overcome the disadvantages of current diagnostic strategies, the detection of UBC biomarkers in easily accessible biofluids, such as urine, represents a promising approach compared with painful biopsies. We investigated the levels of MMP23 genes (microarray and qPCR) and protein (western blot and enzyme-linked immunosorbent assay) in a set of samples (blood, plasma and urine) from patients with UBC and controls as biomarkers for this cancer. MMP23B and its pseudogene MMP23A resulted downregulated in blood cells from UBC compared with controls (66 cases, 70 controls; adjusted P-value = 0.02 and 0.03, respectively). In contrast, MMP23B protein levels in plasma (53 UBC, 49 controls) and urine (59 UBC, 57 controls) increased in cases, being statistically significant in urine. MMP23B dosage observed in urine samples was related to both tumor risk classification and grading. As the lack of correlation between mRNA and protein levels could be due to a posttranscriptional regulation mediated by microRNAs (miRNAs), we investigated the expression of urinary miRNAs targeting MMP23B. Five miRNAs resulted differentially expressed between cases and controls. We reported the first evidence of MMP23B secretion in plasma and urine, suggesting a role of this poorly characterized metalloproteinase in UBC as a potential non-invasive biomarker for this cancer. Further analyses are needed to elucidate the mechanism of regulation of MMP23B expression by miRNAs.

Reproducible bioinformatics project: a community for reproducible bioinformatics analysis pipelines

Background

Reproducibility of a research is a key element in the modern science and it is mandatory for any industrial application. It represents the ability of replicating an experiment independently by the location and the operator. Therefore, a study can be considered reproducible only if all used data are available and the exploited computational analysis workflow is clearly described. However, today for reproducing a complex bioinformatics analysis, the raw data and the list of tools used in the workflow could be not enough to guarantee the reproducibility of the results obtained. Indeed, different releases of the same tools and/or of the system libraries (exploited by such tools) might lead to sneaky reproducibility issues.

Results

To address this challenge, we established the Reproducible Bioinformatics Project (RBP), which is a non-profit and open-source project, whose aim is to provide a schema and an infrastructure, based on docker images and R package, to provide reproducible results in Bioinformatics. One or more Docker images are then defined for a workflow (typically one for each task), while the workflow implementation is handled via R-functions embedded in a package available at github repository. Thus, a bioinformatician participating to the project has firstly to integrate her/his workflow modules into Docker image(s) exploiting an Ubuntu docker image developed ad hoc by RPB to make easier this task. Secondly, the workflow implementation must be realized in R according to an R-skeleton function made available by RPB to guarantee homogeneity and reusability among different RPB functions. Moreover she/he has to provide the R vignette explaining the package functionality together with an example dataset which can be used to improve the user confidence in the workflow utilization.

Conclusions

Reproducible Bioinformatics Project provides a general schema and an infrastructure to distribute robust and reproducible workflows. Thus, it guarantees to final users the ability to repeat consistently any analysis independently by the used UNIX-like architecture.

Transcriptome Signatures of Canine Mammary Gland Tumors and Its Comparison to Human Breast Cancers

Breast cancer (BC)/mammary gland carcinoma (MGC) is the most frequently diagnosed and leading cause of cancer-related mortality in both women and canines. To better understand both canine MGC and human BC-specific genes, we sequenced RNAs obtained from eight pairs of carcinomas and adjacent normal tissues in dogs. By comprehensive transcriptome analysis, 351 differentially expressed genes (DEGs) were identified in overall canine MGCs. Based on the DEGs, comparative analysis revealed correlation existing among the three histological subtypes of canine MGC (ductal, simple, and complex) and four molecular subtypes of human BC (HER2+, ER+, ER&HER2+, and TNBC). Eight DEGs shared by all three subtypes of canine MGCs had been previously reported as cancer-associated genes in human studies. Gene ontology and pathway analyses using the identified DEGs revealed that the biological processes of cell proliferation, adhesion, and inflammatory responses are enriched in up-regulated MGC DEGs. In contrast, fatty acid homeostasis and transcription regulation involved in cell fate commitment were down-regulated in MGC DEGs. Moreover, correlations are demonstrated between upstream promoter transcripts and DEGs. Canine MGC- and subtype-enriched gene expression allows us to better understand both human BC and canine MGC, yielding new insight into the development of biomarkers and targets for both diseases.

Evaluation of commercially available small RNASeq library preparation kits using low input RNA

Background

Evolving interest in comprehensively profiling the full range of small RNAs present in small tissue biopsies and in circulating biofluids, and how the profile differs with disease, has launched small RNA sequencing (RNASeq) into more frequent use. However, known biases associated with small RNASeq, compounded by low RNA inputs, have been both a significant concern and a hurdle to widespread adoption. As RNASeq is becoming a viable choice for the discovery of small RNAs in low input samples and more labs are employing it, there should be benchmark datasets to test and evaluate the performance of new sequencing protocols and operators. In a recent publication from the National Institute of Standards and Technology, Pine et al., 2018, the investigators used a commercially available set of three tissues and tested performance across labs and platforms.

Results

In this paper, we further tested the performance of low RNA input in three commonly used and commercially available RNASeq library preparation kits; NEB Next, NEXTFlex, and TruSeq small RNA library preparation. We evaluated the performance of the kits at two different sites, using three different tissues (brain, liver, and placenta) with high (1 μg) and low RNA (10 ng) input from tissue samples, or 5.0, 3.0, 2.0, 1.0, 0.5, and 0.2 ml starting volumes of plasma. As there has been a lack of robust validation platforms for differentially expressed miRNAs, we also compared low input RNASeq data with their expression profiles on three different platforms (Abcam Fireplex, HTG EdgeSeq, and Qiagen miRNome).

Conclusions

The concordance of RNASeq results on these three platforms was dependent on the RNA expression level; the higher the expression, the better the reproducibility. The results provide an extensive analysis of small RNASeq kit performance using low RNA input, and replication of these data on three downstream technologies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4726-6) contains supplementary material, which is available to authorized users.

microRNA profiles in urine by next-generation sequencing can stratify bladder cancer subtypes

Bladder cancer (BC) is the most frequent malignancy of the urinary tract with a high incidence in men and smokers. Currently, there are no non-invasive markers useful for BC diagnosis and subtypes classification that could overcome invasive procedures such as cystoscopy. Dysregulated miRNA profiles have been associated with numerous cancers, including BC. Cell-free miRNAs are abundantly present in a variety of biofluids including urine and make them promising candidates in cancer biomarker discovery.

In the present study, the identification of miRNA fingerprints associated with different BC status was performed by next-generation sequencing on urine samples from 66 BC and 48 controls. Three signatures based on dysregulated miRNAs have been identified by regression models, assessing the power to discriminate different BC subtypes. Altered miRNAs according to invasiveness and grade were validated by qPCR on 112 cases and 65 controls (among which 46 cases and 16 controls were an independent group of subjects while the rest were replica samples).

The area under the curve (AUC) computed including three miRNAs (miR-30a-5p, let-7c-5p and miR-486-5p) altered in all BC subtypes showed a significantly increased accuracy in the discrimination of cases and controls (AUC model = 0.70; p-value = 0.01).

In conclusions, the non-invasive detection in urine of a selected number of miRNAs altered in different BC subtypes could lead to an accurate early diagnosis of cancer and stratification of patients.