Evaluation of methodologies for microRNA biomarker detection by next generation sequencing

MicroRNAs as Bile-based biomarkers in pancreaticobiliary cancers (MIRABILE): a cohort study

BACKGROUND

Biliary obstruction can be due to both malignant and benign pancreaticobiliary disease. Currently, there are no biomarkers that can accurately help make this distinction. MicroRNAs (miRNAs) are stable molecules in tissue and biofluids that are commonly deregulated in cancer. The MIRABILE study aimed to identify miRNAs in bile that can differentiate malignant from benign pancreaticobiliary disease.

MATERIALS AND METHODS

There were 111 patients recruited prospectively at endoscopic retrograde cholangiopancreatography (ERCP) or percutaneous transhepatic cholangiography (PTC) for obstructive jaundice, and bile was aspirated for cell-free RNA (cfRNA) extraction and analysis. In a discovery cohort of 78 patients (27 with pancreatic ductal adenocarcinoma (PDAC), 14 cholangiocarcinoma (CCA), 37 benign disease), cfRNA was subjected to small-RNA sequencing. LASSO regression was used to define bile miRNA signatures, and NormFinder to identify endogenous controls. In a second cohort of 87 patients (34 PDAC, 14 CCA, 39 benign disease), RT-qPCR was used for validation.

RESULTS

LASSO regression identified 14 differentially-expressed bile miRNAs of which 6 were selected for validation. When comparing malignant and benign pancreaticobiliary disease, bile miR-340 and miR-182 were validated and significantly differentially expressed (P<0.05 and P<0.001, respectively). This generated an AUC of 0.79 (95%CI 0.70-0.88, sensitivity 65%; specificity 82%) in predicting malignant disease.

CONCLUSION

Bile collected during biliary drainage contains miRNAs able to differentiate benign from malignant pancreaticobiliary diseases in patients with obstructive jaundice. These bile miRNAs have the potential to increase diagnostic accuracy.

Detection and analysis of plasma lncRNA, miRNA and mRNA profile in preterm birth with intraventricular hemorrhage

Intraventricular hemorrhage (IVH) is a cause of morbidity and mortality in preterm infants and is strongly associated with adverse neurological outcomes. The incidence of severe IVH (grade 3 or 4) has persisted despite the overall decline in IVH. IVH has been attributed to changes in cerebral blood flow to the immature germinal matrix microvasculature. The cascade of adverse events following IVH includes inflammation, white matter injury, and delayed oligodendrial maturation. In this study, we aimed to identify long non-coding RNA (lncRNA), microRNA (miRNA), and messenger RNA (mRNA) expression in the peripheral blood of preterm infants with IVH compared to normal controls, resulting in the finding of novel biomarkers for IVH. We conducted transcriptome sequencing and small RNA sequencing for identifying differential expression of RNA in preterm infants with IVH. We identified differentially expressed 47 lncRNAs, 95 miRNAs, and 1,370 mRNAs in preterm infants with IVH compared to normal control. Particularly, lncRNA H19 exhibited significantly high expression in preterm infants with IVH. The functional analysis revealed that differentially expressed RNAs in preterm infants with IVH were associated with ferroptosis, heme metabolism, and immune response such as lymphocyte activation and interferon response. In conclusion, these results demonstrate the potential of lncRNA, miRNA, mRNA as possible diagnostic and prognostic biomarkers for IVH.

Spitting in the wind?-The challenges of RNA sequencing for biomarker discovery from saliva

Forensic trace contextualization, i.e., assessing information beyond who deposited a biological stain, has become an issue of great and steadily growing importance in forensic genetic casework and research. The human transcriptome encodes a wide variety of information and thus has received increasing interest for the identification of biomarkers for different aspects of forensic trace contextualization over the past years. Massively parallel sequencing of reverse-transcribed RNA ("RNA sequencing") has emerged as the gold standard technology to characterize the transcriptome in its entirety and identify RNA markers showing significant expression differences not only between different forensically relevant body fluids but also within a single body fluid between forensically relevant conditions of interest. Here, we analyze the quality and composition of four RNA sequencing datasets (whole transcriptome as well as miRNA sequencing) from two different research projects (the RNAgE project and the TrACES project), aiming at identifying contextualizing forensic biomarker from the forensically relevant body fluid saliva. We describe and characterize challenges of RNA sequencing of saliva samples arising from the presence of oral bacteria, the heterogeneity of sample composition, and the confounding factor of degradation. Based on these observations, we formulate recommendations that might help to improve RNA biomarker discovery from the challenging but forensically relevant body fluid saliva.

Harnessing microRNA-enriched extracellular vesicles for liquid biopsy

Extracellular microRNAs (miRNAs) can be detected in body fluids and hold great potential as cancer biomarkers. Extracellular miRNAs are protected from degradation by binding various proteins and through their packaging into extracellular vesicles (EVs). There is evidence that the diagnostic performance of cancer-associated extracellular miRNAs can be improved by assaying EV-miRNA instead of total cell-free miRNA, but several challenges have hampered the advancement of EV-miRNA in liquid biopsy. Because almost all types of cells release EVs, cancer cell-derived EVs might constitute only a minor fraction of EVs in body fluids of cancer patients with low volume disease. Furthermore, a given cell type can release several subpopulations of EVs that vary in their cargo, and there is evidence that the majority of EVs contain low copy numbers of miRNAs. In this mini-review, we discuss the potential of several candidate EV membrane proteins such as CD147 to define cancer cell-derived EVs, and approaches by which subpopulations of miRNA-rich EVs in body fluids might be identified. By integrating these insights, we discuss strategies by which EVs that are both cancer cell-derived and miRNA-rich could be isolated to enhance the diagnostic performance of extracellular miRNAs.

MicroRNAs in Parkinson's disease: a systematic review and diagnostic accuracy meta-analysis

Current clinical tests for Parkinson's disease (PD) provide insufficient diagnostic accuracy leading to an urgent need for improved diagnostic biomarkers. As microRNAs (miRNAs) are promising biomarkers of various diseases, including PD, this systematic review and meta-analysis aimed to assess the diagnostic accuracy of biofluid miRNAs in PD. All studies reporting data on miRNAs expression in PD patients compared to controls were included. Gene targets and significant pathways associated with miRNAs expressed in more than 3 biofluid studies with the same direction of change were analyzed using target prediction and enrichment analysis. A bivariate model was used to calculate sensitivity, specificity, likelihood ratios, and diagnostic odds ratio. While miR-24-3p and miR-214-3p were the most reported miRNA (7 each), miR-331-5p was found to be consistently up regulated in 4 different biofluids. Importantly, miR-19b-3p, miR-24-3p, miR-146a-5p, and miR-221-3p were reported in multiple studies without conflicting directions of change in serum and bioinformatic analysis found the targets of these miRNAs to be associated with pathways important in PD pathology. Of the 102 studies from the systematic review, 15 studies reported sensitivity and specificity data on combinations of miRNAs and were pooled for meta-analysis. Studies (17) reporting sensitivity and specificity data on single microRNA were pooled in a separate meta-analysis. Meta-analysis of the combinations of miRNAs (15 studies) showed that biofluid miRNAs can discriminate between PD patients and controls with good diagnostic accuracy (sensitivity = 0.82, 95% CI 0.76-0.87; specificity = 0.80, 95% CI 0.74-0.84; AUC = 0.87, 95% CI 0.83-0.89). However, we found multiple studies included more males with PD than any other group therefore possibly introducing a sex-related selection bias. Overall, our study captures key miRNAs which may represent a point of focus for future studies and the development of diagnostic panels whilst also highlighting the importance of appropriate study design to develop representative biomarker panels for the diagnosis of PD.

Discovery of the major 15-30 nt mammalian small RNAs, their biogenesis and function

Small RNAs (sRNAs) within 15-30 nt such as miRNA, tsRNA, srRNA with 3'-OH have been identified. However, whether these sRNAs are the major 15-30 nt sRNAs is still unknown. Here we show about 90% mammalian sRNAs within 15-30 nt end with 2',3'-cyclic phosphate (3'-cP). TANT-seq was developed to simultaneously profile sRNAs with 3'-cP (sRNA-cPs) and sRNA-OHs, and huge amount of sRNA-cPs were detected. Surprisingly, sRNA-cPs and sRNA-OHs usually have distinct sequences. The data from TANT-seq were validated by a novel method termed TE-qPCR, and Northern blot. Furthermore, we found that Angiogenin and RNase 4 contribute to the biogenesis of sRNA-cPs. Moreover, much more sRNA-cPs than sRNA-OHs bind to Ago2, and can regulate gene expression. Particularly, snR-2-cP regulates Bcl2 by targeting to its 3'UTR dependent on Ago2, and subsequently regulates apoptosis. In addition, sRNA-cPs can guide the cleavage of target RNAs in Ago2 complex as miRNAs without the requirement of 3'-cP. Our discovery greatly expands the repertoire of mammalian sRNAs, and provides strategies and powerful tools towards further investigation of sRNA-cPs.

microRNA-based predictor for diagnosis of frontotemporal dementia

AIMS

This study aimed to explore the non-linear relationships between cell-free microRNAs (miRNAs) and their contribution to prediction of Frontotemporal dementia (FTD), an early onset dementia that is clinically heterogeneous, and too often suffers from delayed diagnosis.

METHODS

We initially studied a training cohort of 219 subjects (135 FTD and 84 non-neurodegenerative controls) and then validated the results in a cohort of 74 subjects (33 FTD and 41 controls).

RESULTS

On the basis of cell-free plasma miRNA profiling by next generation sequencing and machine learning approaches, we develop a non-linear prediction model that accurately distinguishes FTD from non-neurodegenerative controls in ~90% of cases.

CONCLUSIONS

The fascinating potential of diagnostic miRNA biomarkers might enable early-stage detection and a cost-effective screening approach for clinical trials that can facilitate drug development.

Circulating microRNAs as promising testicular translatable safety biomarkers: current state and future perspectives

Drug-induced testicular injury (DITI) is one of the often-observed and challenging safety issues seen during drug development. Semen analysis and circulating hormones currently utilized have significant gaps in their ability to detect testicular damage accurately. In addition, no biomarkers enable a mechanistic understanding of the damage to the different regions of the testis, such as seminiferous tubules, Sertoli, and Leydig cells. MicroRNAs (miRNAs) are a class of non-coding RNAs that modulate gene expression post-transcriptionally and have been indicated to regulate a wide range of biological pathways. Circulating miRNAs can be measured in the body fluids due to tissue-specific cell injury/damage or toxicant exposure. Therefore, these circulating miRNAs have become attractive and promising non-invasive biomarkers for assessing drug-induced testicular injury, with several reports on their use as safety biomarkers for monitoring testicular damage in preclinical species. Leveraging emerging tools such as 'organs-on-chips' that can emulate the human organ's physiological environment and function is starting to enable biomarker discovery, validation, and clinical translation for regulatory qualification and implementation in drug development.

Identification of extremely GC-rich micro RNAs for RT-qPCR data normalization in human plasma

We aimed at extending the repertoire of high-quality miRNA normalizers for reverse transcription-quantitative PCR (RT-qPCR) of human plasma with special emphasis on the extremely guanine-cytosine-rich portion of the miRNome. For high-throughput selection of stable candidates, microarray technology was preferred over small-RNA sequencing (sRNA-seq) since the latter underrepresented miRNAs with a guanine-cytosine (GC) content of at least 75% (p = 0.0002, n = 2). miRNA abundances measured on the microarray were ranked for consistency and uniformity using nine normalization approaches. The eleven most stable sequences included miRNAs of moderate, but also extreme GC content (45%-65%: miR-320d, miR-425-5p, miR-185-5p, miR-486-5p; 80%-95%: miR-1915-3p, miR-3656-5p, miR-3665-5p, miR-3960-5p, miR-4488-5p, miR-4497 and miR-4787-5p). In contrast, the seven extremely GC-rich miRNAs were not found in the two plasma miRNomes screened by sRNA-seq. Stem-loop RT-qPCR was employed for stability verification in 32 plasma samples of healthy male Caucasians (age range: 18-55 years). In general, inter-individual variance of miRNA abundance was low or very low as indicated by coefficient of variation (CV) values of 0.6%-8.2%. miR-3665 and miR-1915-3p outperformed in this analysis (CVs: 0.6 and 2.4%, respectively). The eight most stable sequences included four extremely GC-rich miRNAs (miR-1915-3p, miR-3665, miR-4787-5p and miR-4497). The best-performing duo normalization factor (NF) for the condition of human plasma, miR-320d and miR-4787-5p, also included a GC-extreme miRNA. In summary, the identification of extremely guanine-cytosine-rich plasma normalizers will help to increase accuracy of PCR-based miRNA quantification, thus raise the potential that miRNAs become markers for psychological stress reactions or early and precise diagnosis of clinical phenotypes. The novel miRNAs might also be useful for orthologous contexts considering their conservation in related animal genomes.

Potential of Non-Coding RNA as Biomarkers for Progressive Supranuclear Palsy

Objective markers for the neurodegenerative disorder progressive supranuclear palsy (PSP) are needed to provide a timely diagnosis with greater certainty. Non-coding RNA (ncRNA), including microRNA, piwi-interacting RNA, and transfer RNA, are good candidate markers in other neurodegenerative diseases, but have not been investigated in PSP. Therefore, as proof of principle, we sought to identify whether they were dysregulated in matched serum and cerebrospinal fluid (CSF) samples of patients with PSP. Small RNA-seq was undertaken on serum and CSF samples from healthy controls (n = 20) and patients with PSP (n = 31) in two cohorts, with reverse transcription-quantitative PCR (RT-qPCR) to confirm their dysregulation. Using RT-qPCR, we found in serum significant down-regulation in hsa-miR-92a-3p, hsa-miR-626, hsa-piR-31068, and tRNA-ValCAC. In CSF, both hsa-let-7a-5p and hsa-piR-31068 showed significant up-regulation, consistent with their changes observed in the RNA-seq results. Interestingly, we saw no correlation in the expression of hsa-piR-31068 within our matched serum and CSF samples, suggesting there is no common dysregulatory mechanism between the two biofluids. While these changes were in a small cohort of samples, we have provided novel evidence that ncRNA in biofluids could be possible diagnostic biomarkers for PSP and further work will help to expand this potential.

A Bioinformatics Approach to MicroRNA-Sequencing Analysis Based on Human Saliva Samples of Patients with Endometriosis

Endometriosis, defined by the presence of endometrium-like tissue outside the uterus, affects 2-10% of the female population, i.e., around 190 million women, worldwide. The aim of the prospective ENDO-miRNA study was to develop a bioinformatics approach for microRNA-sequencing analysis of 200 saliva samples for miRNAome expression and to test its diagnostic accuracy for endometriosis. Among the 200 patients, 76.5% (n = 153) had confirmed endometriosis and 23.5% (n = 47) had no endometriosis (controls). Small RNA-seq of 200 saliva samples yielded ~4642 M raw sequencing reads (from ~13.7 M to ~39.3 M reads/sample). The number of expressed miRNAs ranged from 1250 (outlier) to 2561 per sample. Some 2561 miRNAs were found to be differentially expressed in the saliva samples of patients with endometriosis compared with the control patients. Among these, 1.17% (n = 30) were up- or downregulated. Among these, the F1-score, sensitivity, specificity, and AUC ranged from 11-86.8%, 5.8-97.4%, 10.6-100%, and 39.3-69.2%, respectively. Here, we report a bioinformatic approach to saliva miRNA sequencing and analysis. We underline the advantages of using saliva over blood in terms of ease of collection, reproducibility, stability, safety, non-invasiveness. This report describes the whole saliva transcriptome to make miRNA quantification a validated, standardized, and reliable technique for routine use. The methodology could be applied to build a saliva signature of endometriosis.

Computational methods and translational applications for targeted next-generation sequencing platforms

During the past decade, next-generation sequencing (NGS) technologies have become widely adopted in cancer research and clinical care. Common applications within the clinical setting include patient stratification into relevant molecular subtypes, identification of biomarkers of response and resistance to targeted and systemic therapies, assessment of heritable cancer risk based on known pathogenic variants, and longitudinal monitoring of treatment response. The need for efficient downstream processing and reliable interpretation of sequencing data has led to the development of novel algorithms and computational pipelines, as well as structured knowledge bases that link genomic alterations to currently available drugs and ongoing clinical trials. Cancer centers around the world use different types of targeted solid-tissue and blood based NGS assays to analyze the genomic and transcriptomic profile of patients as part of their routine clinical care. Recently, cross-institutional collaborations have led to the creation of large pooled datasets that can offer valuable insights into the genomics of rare cancers.

Global MicroRNA Profiling of Vascular Endothelial Cells

MicroRNA sequencing (miRNA-seq) enables the detection and characterization of the cell miRNome, including miRNA isoforms (isomiRs) and novel miRNA species. In roughly half of the cases, the most abundant isomiR in the cells is not the reference miRNA given in miRBase, which highlights the importance of isomiR-specific analysis. Here, we describe a gel-free protocol for global miRNA profiling in vascular endothelial cells and the main steps of the subsequent data analysis with two alternative analysis methods. In addition to endothelial cells, the protocol is suitable for other cell and tissue types and has been successfully used to obtain miRNA-seq data from human cardiac tissue, plasma, pericardial fluid, and biofluid exosomes.

Circulating miR-181 is a prognostic biomarker for amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a relentless neurodegenerative disease of the human motor neuron system, where variability in progression rate limits clinical trial efficacy. Therefore, better prognostication will facilitate therapeutic progress. In this study, we investigated the potential of plasma cell-free microRNAs (miRNAs) as ALS prognostication biomarkers in 252 patients with detailed clinical phenotyping. First, we identified, in a longitudinal cohort, miRNAs whose plasma levels remain stable over the course of disease. Next, we showed that high levels of miR-181, a miRNA enriched in neurons, predicts a greater than two-fold risk of death in independent discovery and replication cohorts (126 and 122 patients, respectively). miR-181 performance is similar to neurofilament light chain (NfL), and when combined together, miR-181 + NfL establish a novel RNA-protein biomarker pair with superior prognostication capacity. Therefore, plasma miR-181 alone and a novel miRNA-protein biomarker approach, based on miR-181 + NfL, boost precision of patient stratification. miR-181-based ALS biomarkers encourage additional validation and might enhance the power of clinical trials.

Plasma Exosome-Enriched Extracellular Vesicles From Lactating Mothers With Type 1 Diabetes Contain Aberrant Levels of miRNAs During the Postpartum Period

Type 1 diabetes is an immune-driven disease, where the insulin-producing beta cells from the pancreatic islets of Langerhans becomes target of immune-mediated destruction. Several studies have highlighted the implication of circulating and exosomal microRNAs (miRNAs) in type 1 diabetes, underlining its biomarker value and novel therapeutic potential. Recently, we discovered that exosome-enriched extracellular vesicles carry altered levels of both known and novel miRNAs in breast milk from lactating mothers with type 1 diabetes. In this study, we aimed to characterize exosomal miRNAs in the circulation of lactating mothers with and without type 1 diabetes, hypothesizing that differences in type 1 diabetes risk in offspring from these groups are reflected in the circulating miRNA profile. We performed small RNA sequencing on exosome-enriched extracellular vesicles extracted from plasma of 52 lactating mothers around 5 weeks postpartum (26 with type 1 diabetes and 26 age-matched controls), and found a total of 2,289 miRNAs in vesicles from type 1 diabetes and control libraries. Of these, 176 were differentially expressed in plasma from mothers with type 1 diabetes (167 upregulated; 9 downregulated, using a cut-off of abs(log2FC) >1 and FDR adjusted p-value <0.05). Extracellular vesicles were verified by nanoparticle tracking analysis, transmission electron microscopy and immunoblotting. Five candidate miRNAs were selected based on their involvement in diabetes and immune modulation/beta-cell functions: hsa-miR-127-3p, hsa-miR-146a-5p, hsa-miR-26a-5p, hsa-miR-24-3p and hsa-miR-30d-5p. Real-time qPCR validation confirmed that hsa-miR-146a-5p, hsa-miR-26a-5p, hsa-miR-24-3p, and hsa-miR-30d-5p were significantly upregulated in lactating mothers with type 1 diabetes as compared to lactating healthy mothers. To determine possible target genes and affected pathways of the 5 miRNA candidates, computational network-based analyses were carried out with TargetScan, mirTarBase, QIAGEN Ingenuity Pathway Analysis and PantherDB database. The candidates showed significant association with inflammatory response and cytokine and chemokine mediated signaling pathways. With this study, we detect aberrant levels of miRNAs within plasma extracellular vesicles from lactating mothers with type 1 diabetes during the postpartum period, including miRNAs with associations to disease pathogenesis and inflammatory responses.

Circulating non-coding RNAs as a diagnostic and management biomarker for breast cancer: current insights

Cancer biomarkers can be used to determine the molecular status of a tumor or its metastases, which either release them directly into body fluids or indirectly through disruption of tumor/metastatic tissue. New minimally invasive and repeatable sample collection methods, such as liquid biopsy, have been developed in the last decade to apply cancer knowledge and track its progression. Circulating non-coding RNAs, which include microRNAs, long non-coding RNAs, and PIWI-interacting RNAs, are increasingly being recognized as potential cancer biomarkers. The growing understanding of cancer's molecular pathogenesis, combined with the rapid development of new molecular techniques, encourages the study of early molecular alterations associated with cancer development in body fluids. Specific genetic and epigenetic changes in circulating free RNA (cf-RNA) in plasma, serum, and urine could be used as diagnostic biomarkers for a variety of cancers. Only a subset of these cf-RNAs have been studied in breast cancer, with the most extensive research focusing on cf-miRNA in plasma. These findings pave the way for immediate use of selected cf-RNAs as biomarkers in breast cancer liquid biopsy, as well as additional research into other cf-RNAs to advance.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase

Retrotransposons are noninfectious, mobile genetic elements that proliferate in host genomes via an RNA intermediate that is copied into DNA by a reverse transcriptase (RT) enzyme. RTs are important for biotechnological applications involving information capture from RNA since RNA is first converted into complementary DNA for detection or sequencing. Here, we biochemically characterized RTs from two retroelements and uncovered several activities that allowed us to design a streamlined, efficient workflow for determining the inventory of RNA sequences in processed RNA pools. The unique properties of nonretroviral RT activities obviate many technical issues associated with current methods of RNA sequence analysis, with wide applications in research, biotechnology, and diagnostics.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale. Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3′ ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

High-Throughput Sequencing of Circulating MicroRNAs in Plasma and Serum during Pregnancy Progression

Although circulating microRNAs (miRNAs) in maternal blood may play an important role in regulation of pregnancy progression and serve as non-invasive biomarkers for different gestation complications, little is known about their profile in blood during normally developing pregnancy. In this study we evaluated the miRNA profiles in paired plasma and serum samples from pregnant women without health or gestational abnormalities at three time points using high-throughput sequencing technology. Sequencing revealed that the percentage of miRNA reads in plasma and serum decreased by a third compared to first and second trimesters. We found two miRNAs in plasma (hsa-miR-7853-5p and hsa-miR-200c-3p) and 10 miRNAs in serum (hsa-miR-203a-5p, hsa-miR-495-3p, hsa-miR-4435, hsa-miR-340-5p, hsa-miR-4417, hsa-miR-1266-5p, hsa-miR-4494, hsa-miR-134-3p, hsa-miR-5008-5p, and hsa-miR-6756-5p), that exhibit level changes during pregnancy (p-value adjusted < 0.05). In addition, we observed differences for 36 miRNAs between plasma and serum (p-value adjusted < 0.05), which should be taken into consideration when comparing the results between studies performed using different biosample types. The results were verified by analysis of three miRNAs using qRT-PCR (p < 0.05). The present study confirms that the circulating miRNA profile in blood changes during gestation. Our results set the basis for further investigation of molecular mechanisms, involved in regulation of pregnancy, and the search for biomarkers of gestation abnormalities.

GH Overexpression Alters Spermatic Cells MicroRNAome Profile in Transgenic Zebrafish

Overexpression of growth hormone (GH) in gh-transgenic zebrafish of a highly studied lineage F0104 has earlier been reported to cause increased muscle growth. In addition to this, GH affects a broad range of cellular processes in transgenic fish, such as morphology, physiology, and behavior. Reports show changes such as decreased sperm quality and reduced reproductive performance in transgenic males. It is hypothesized that microRNAs are directly involved in the regulation of fertility potential during spermatogenesis. The primary aim of our study was to verify whether gh overexpression disturbs the sperm miRNA profile and influences the sperm quality in transgenic zebrafish. We report a significant increase in body weight of gh-transgenic males along with associated reduced sperm motility and other kinetic parameters in comparison to the non-transgenic group. MicroRNA transcriptome sequencing of gh-transgenic zebrafish sperms revealed expressions of 186 miRNAs, among which six miRNA were up-regulated (miR-146b, miR-200a-5p, miR-146a, miR-726, miR-184, and miR-738) and sixteen were down-regulated (miR-19d-3p, miR-126a-5p, miR-126b-5p, miR-22a-5p, miR-16c-5p, miR-20a-5p, miR-126b-3p, miR-107a-3p, miR-93, miR-2189, miR-202–5p, miR-221–3p, miR-125a, miR-125b-5p, miR-126a-3p, and miR-30c-5p) in comparison to non-transgenic zebrafish. Some of the dysregulated miRNAs were previously reported to be related to abnormalities in sperm quality and reduced reproduction ability in other species. In this study, an average of 134 differentially expressed miRNAs-targeted genes were predicted using the in silico approach. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis demonstrated that the genes of affected pathways were primarily related to spermatogenesis, sperm motility, and cell apoptosis. Our results suggested that excess GH caused a detrimental effect on sperm microRNAome, consequently reducing the sperm quality and reproductive potential of zebrafish males.

Quantitative mapping of the cellular small RNA landscape with AQRNA-seq

Current next-generation RNA-sequencing (RNA-seq) methods do not provide accurate quantification of small RNAs within a sample, due to sequence-dependent biases in capture, ligation and amplification during library preparation. We present a method, absolute quantification RNA-sequencing (AQRNA-seq), that minimizes biases and provides a direct, linear correlation between sequencing read count and copy number for all small RNAs in a sample. Library preparation and data processing were optimized and validated using a 963-member microRNA reference library, oligonucleotide standards of varying length, and RNA blots. Application of AQRNA-seq to a panel of human cancer cells revealed >800 detectable miRNAs that varied during cancer progression, while application to bacterial transfer RNA pools, with the challenges of secondary structure and abundant modifications, revealed 80-fold variation in tRNA isoacceptor levels, stress-induced site-specific tRNA fragmentation, quantitative modification maps, and evidence for stress-induced, tRNA-driven, codon-biased translation. AQRNA-seq thus provides a versatile means to quantitatively map the small RNA landscape in cells.

Small RNA-Sequencing: Approaches and Considerations for miRNA Analysis

MicroRNAs (miRNAs) are a class of small RNA molecules that have an important regulatory role in multiple physiological and pathological processes. Their disease-specific profiles and presence in biofluids are properties that enable miRNAs to be employed as non-invasive biomarkers. In the past decades, several methods have been developed for miRNA analysis, including small RNA sequencing (RNA-seq). Small RNA-seq enables genome-wide profiling and analysis of known, as well as novel, miRNA variants. Moreover, its high sensitivity allows for profiling of low input samples such as liquid biopsies, which have now found applications in diagnostics and prognostics. Still, due to technical bias and the limited ability to capture the true miRNA representation, its potential remains unfulfilled. The introduction of many new small RNA-seq approaches that tried to minimize this bias, has led to the existence of the many small RNA-seq protocols seen today. Here, we review all current approaches to cDNA library construction used during the small RNA-seq workflow, with particular focus on their implementation in commercially available protocols. We provide an overview of each protocol and discuss their applicability. We also review recent benchmarking studies comparing each protocol's performance and summarize the major conclusions that can be gathered from their usage. The result documents variable performance of the protocols and highlights their different applications in miRNA research. Taken together, our review provides a comprehensive overview of all the current small RNA-seq approaches, summarizes their strengths and weaknesses, and provides guidelines for their applications in miRNA research.

Size-Exclusion Chromatography Separation Reveals That Vesicular and Non-Vesicular Small RNA Profiles Differ in Cell Free Urine

Urinary extracellular vesicles (EVs) and their RNA cargo are a novel source of biomarkers for various diseases. We aimed to identify the optimal method for isolating small (<200 nm) EVs from human urine prior to small RNA analysis. EVs from filtered healthy volunteer urine were concentrated using three methods: ultracentrifugation (UC); a precipitation-based kit (PR); and ultrafiltration (UF). EVs were further purified by size-exclusion chromatography (SEC). EV preparations were analysed with transmission electron microscopy (TEM), Western blotting, nanoparticle tracking analysis (NTA) and an Agilent Bioanalyzer Small RNA kit. UF yielded the highest number of particles both before and after SEC. Small RNA analysis from UF-concentrated urine identified two major peaks at 10–40 nucleotides (nt) and 40–80 nt. In contrast, EV preparations obtained after UC, PR or SEC combined with any concentrating method, contained predominantly 40–80 nt sized small RNA. Protein fractions from UF+SEC contained small RNA of 10–40 nt in size (consistent with miRNAs). These data indicate that most of the microRNA-sized RNAs in filtered urine are not associated with small-sized EVs, and highlights the importance of removing non-vesicular proteins and RNA from urine EV preparations prior to small RNA analysis.

Methods to Investigate miRNA Function: Focus on Platelet Reactivity

MicroRNAs (miRNAs) are small noncoding RNAs modulating protein production. They are key players in regulation of cell function and are considered as biomarkers in several diseases. The identification of the proteins they regulate, and their impact on cell physiology, may delineate their role as diagnostic or prognostic markers and identify new therapeutic strategies. During the last 3 decades, development of a large panel of techniques has given rise to multiple models dedicated to the study of miRNAs. Since plasma samples are easily accessible, circulating miRNAs can be studied in clinical trials. To quantify miRNAs in numerous plasma samples, the choice of extraction and purification techniques, as well as normalization procedures, are important for comparisons of miRNA levels in populations and over time. Recent advances in bioinformatics provide tools to identify putative miRNAs targets that can then be validated with dedicated assays. In vitro and in vivo approaches aim to functionally validate candidate miRNAs from correlations and to understand their impact on cellular processes. This review describes the advantages and pitfalls of the available techniques for translational research to study miRNAs with a focus on their role in regulating platelet reactivity.

Human milk extracellular vesicle miRNA expression and associations with maternal characteristics in a population-based cohort from the Faroe Islands

Human milk plays a critical role in infant development and health, particularly in cognitive, immune, and cardiometabolic functions. Milk contains extracellular vesicles (EVs) that can transport biologically relevant cargo from mother to infant, including microRNAs (miRNAs). We aimed to characterize milk EV-miRNA profiles in a human population cohort, assess potential pathways and ontology, and investigate associations with maternal characteristics. We conducted the first study to describe the EV miRNA profile of human milk in 364 mothers from a population-based mother-infant cohort in the Faroe Islands using small RNA sequencing. We detected 1523 miRNAs with ≥ one read in 70% of samples. Using hierarchical clustering, we determined five EV-miRNA clusters, the top three consisting of 15, 27 and 67 miRNAs. Correlation coefficients indicated that the expression of many miRNAs within the top three clusters was highly correlated. Top-cluster human milk EV-miRNAs were involved in pathways enriched for the endocrine system, cellular community, neurodevelopment, and cancers. miRNA expression was associated with time to milk collection post-delivery, maternal body mass index, and maternal smoking, but not maternal parity. Future studies investigating determinants of human EV-miRNAs and associated health outcomes are needed to elucidate the role of human milk EV-miRNAs in health and disease.

MicroRNA-1291 Is Associated With Locoregional Metastases in Patients With Early-Stage Breast Cancer

Evidence that microRNAs (miRNAs) regulate the various steps of metastasis is increasing. Several studies have looked at the miRNA expression profile in primary breast tumors but few have compared primary tumor and sentinel lymph node (SLN) metastasis. We correlated the expression of miRNAs with the SLN status and the outcome of axillary lymph node dissection (ALND) in 60 patients with early breast cancer. We profiled the expression of miRNAs in paired breast tumor samples and SLNs using the NextSeq500 Illumina platform and key findings were validated by qPCR. MultiMiR Bioconductor and Reactome pathways analysis were performed to identify target genes and signaling pathways affected by altered expressed miRNAs. Our results show that nine miRNAs were differentially expressed in tumor tissues (q ≤ 0.05). In tumor samples, a 13.5-fold up-regulation of miR-7641-2 (q < 0.001) and a 2.9-fold down-regulation of miR-1291 (q < 0.001) were associated with tumors with positive SLNs. However, only down-regulation of miR-1291 (q = 0.048) remained significant in paired SLNs samples. Interestingly, a 10.5 up-regulation of miR-1291 in SLNs samples was associated with additional axillary lymph node involvement (q < 0.001). The enrichment analyses showed that canonical and non-canonical WNT pathways and negative regulation of various receptor tyrosine kinases signaling pathways were targets of miR-1291 and supports the role of miR-1291 as a tumor suppressor gene (TSG). Further studies are warranted to investigate the use of miR-1291 as a surrogate biomarker of SLN node metastasis in patients with early-stage breast cancer.

Chronic stress promotes EMT-mediated metastasis through activation of STAT3 signaling pathway by miR-337-3p in breast cancer

Chronic stress could induce cancer metastasis by constant activation of the sympathetic nervous system, while cellular mechanism remains obscure. The aim of this research is to explore the metastasis associated negative effect of chronic stress. The analysis of transcriptome sequencing implied that activation of STAT3 signaling pathway by downregulated miR-337-3p might be a potential mechanism to induce epithelial to mesenchymal transition (EMT) of cancer cell and promote metastasis under chronic stress. We also verified this biological process in further experiments. Downregulation of miR-337-3p could downregulate E-cadherin expression and upregulate vimentin expression in vitro and in vivo. STAT3, related signal pathways of which are involved in metastasis regulation, was directly targeted by miR-337-3p. In conclusion, the above results denoted that activation of miR-337-3p/STAT3 axis might be a potential pathway for the increasing metastasis of breast cancer under chronic stress.

Can miRNAs Be Considered as Diagnostic and Therapeutic Molecules in Ischemic Stroke Pathogenesis?-Current Status

Ischemic stroke is one of the leading causes of death worldwide. Clinical manifestations of stroke are long-lasting and causing economic burden on the patients and society. Current therapeutic modalities to treat ischemic stroke (IS) are unsatisfactory due to the intricate pathophysiology and poor functional recovery of brain cellular compartment. MicroRNAs (miRNA) are endogenously expressed small non-coding RNA molecules, which can act as translation inhibitors and play a pivotal role in the pathophysiology associated with IS. Moreover, miRNAs may be used as potential diagnostic and therapeutic tools in clinical practice; yet, the complete role of miRNAs is enigmatic during IS. In this review, we explored the role of miRNAs in the regulation of stroke risk factors viz., arterial hypertension, metabolic disorders, and atherosclerosis. Furthermore, the role of miRNAs were reviewed during IS pathogenesis accompanied by excitotoxicity, oxidative stress, inflammation, apoptosis, angiogenesis, neurogenesis, and Alzheimer's disease. The functional role of miRNAs is a double-edged sword effect in cerebral ischemia as they could modulate pathological mechanisms associated with risk factors of IS. miRNAs pertaining to IS pathogenesis could be potential biomarkers for stroke; they could help researchers to identify a particular stroke type and enable medical professionals to evaluate the severity of brain injury. Thus, ascertaining the role of miRNAs may be useful in deciphering their diagnostic role consequently it is plausible to envisage a suitable therapeutic modality against IS.

Flåm S, Leithaug M, Dalland M, Rayner S, Lie BA, Gilfillan GD. An extension to: Systematic assessment of commercially available low-input miRNA library preparation kits

High-throughput sequencing has emerged as the favoured method to study microRNA (miRNA) expression, but biases introduced during library preparation have been reported. We recently compared the performance (sensitivity, reliability, titration response and differential expression) of six commercially-available kits on synthetic miRNAs and human RNA, where library preparation was performed by the vendors. We hereby supplement this study with data from two further commonly used kits (NEBNext, NEXTflex) whose manufacturers initially declined to participate. NEXTflex demonstrated the highest sensitivity, which may reflect its use of partially-randomized adapter sequences, but overall performance was lower than the QIAseq and TailorMix kits. NEBNext showed intermediate performance. We reaffirm that biases are kit specific, complicating the comparison of miRNA datasets generated using different kits.

Deep Sequencing of Small RNAs from Neurosurgical Extracellular Vesicles Substantiates miR-486-3p as a Circulating Biomarker that Distinguishes Glioblastoma from Lower-Grade Astrocytoma Patients

Extracellular vesicles (EVs) play key roles in glioblastoma (GBM; astrocytoma grade IV) biology and are novel sources of biomarkers. EVs released from GBM tumors can cross the blood-brain-barrier into the periphery carrying GBM molecules, including small non-coding RNA (sncRNA). Biomarkers cargoed in circulating EVs have shown great promise for assessing the molecular state of brain tumors in situ. Neurosurgical aspirate fluids captured during tumor resections are a rich source of GBM-EVs isolated directly from tumor microenvironments. Using density gradient ultracentrifugation, EVs were purified from cavitron ultrasonic surgical aspirate (CUSA) washings from GBM (n = 12) and astrocytoma II-III (GII-III, n = 5) surgeries. The sncRNA contents of surgically captured EVs were profiled using the Illumina^® NextSeq^TM 500 NGS System. Differential expression analysis identified 27 miRNA and 10 piRNA species in GBM relative to GII-III CUSA-EVs. Resolved CUSA-EV sncRNAs could discriminate serum-EV sncRNA profiles from GBM and GII-III patients and healthy controls and 14 miRNAs (including miR-486-3p and miR-106b-3p) and cancer-associated piRNAs (piR_016658, _016659, _020829 and _204090) were also significantly expressed in serum-EVs. Circulating EV markers that correlate with histological, neuroradiographic and clinical parameters will provide objective measures of tumor activity and improve the accuracy of GBM tumor surveillance.

Identification of a potential non-coding RNA biomarker signature for amyotrophic lateral sclerosis

Objective biomarkers for the clinically heterogeneous adult-onset neurodegenerative disorder amyotrophic lateral sclerosis are crucial to facilitate assessing emerging therapeutics and improve the diagnostic pathway in what is a clinically heterogeneous syndrome. With non-coding RNA transcripts including microRNA, piwi-RNA and transfer RNA present in human biofluids, we sought to identify whether non-coding RNA in serum could be biomarkers for amyotrophic lateral sclerosis. Serum samples from our Oxford Study for Biomarkers in motor neurone disease/amyotrophic lateral sclerosis discovery cohort of amyotrophic lateral sclerosis patients (n = 48), disease mimics (n = 16) and age- and sex-matched healthy controls (n = 24) were profiled for non-coding RNA expression using RNA-sequencing, which showed a wide range of non-coding RNA to be dysregulated. We confirmed significant alterations with reverse transcription-quantitative PCR in the expression of hsa-miR-16-5p, hsa-miR-21-5p, hsa-miR-92a-3p, hsa-piR-33151, TRV-AAC4-1.1 and TRA-AGC6-1.1. Furthermore, hsa-miR-206, a previously identified amyotrophic lateral sclerosis biomarker, showed a binary-like pattern of expression in our samples. Using the expression of these non-coding RNA, we were able to discriminate amyotrophic lateral sclerosis samples from healthy controls in our discovery cohort using a random forest analysis with 93.7% accuracy with promise in predicting progression rate of patients. Importantly, cross-validation of this novel signature using a new geographically distinct cohort of samples from the United Kingdom and Germany with both amyotrophic lateral sclerosis and control samples (n = 156) yielded an accuracy of 73.9%. The high prediction accuracy of this non-coding RNA-based biomarker signature, even across heterogeneous cohorts, demonstrates the strength of our approach as a novel platform to identify and stratify amyotrophic lateral sclerosis patients.

Elevated Plasma Levels of Circulating Extracellular miR-320a-3p in Patients with Paroxysmal Atrial Fibrillation

The potential of extracellular circulating microRNAs (miRNAs) as non-invasive biomarkers of atrial fibrillation (AF) has been confirmed by a number of recent studies. However, the current data for some miRNAs are controversial and inconsistent, probably due to pre-analytical and methodological differences. In this work, we attempted to fulfill the basic pre-analytical requirements provided for circulating miRNA studies for application to paroxysmal atrial fibrillation (PAF) research. We used quantitative PCR (qPCR) to determine the relative plasma levels of circulating miRNAs expressed in the heart or associated with atrial remodeling or fibrillation with reported altered plasma/serum levels in AF: miR-146a-5p, miR-150-5p, miR-19a-3p, miR-21-5p, miR-29b-3p, miR-320a-3p, miR-328-3p, miR-375-3p, and miR-409-3p. First, in a cohort of 90 adult outpatient clinic patients, we found that the plasma level of miR-320a-3p was elevated in PAF patients compared to healthy controls and hypertensive patients without AF. We further analyzed the impact of medication therapies on miRNA relative levels and found elevated miR-320a-3p levels in patients receiving angiotensin-converting-enzyme inhibitors (ACEI) therapy. Additionally, we found that miR-320a-3p, miR-21-5p, and miR-146a-5p plasma levels positively correlated with the CHA₂DS₂-Vasc score and were elevated in subjects with CHA₂DS₂-Vasc ≥ 2. Our results indicate that, amongst the analyzed miRNAs, miR-320a-3p may be considered as a potential PAF circulating plasma biomarker, leading to speculation as to whether this miRNA is a marker of platelet state change due to ACEI therapy.

Effects of Cardiac Sympathetic Neurodegeneration and PPARγ Activation on Rhesus Macaque Whole Blood miRNA and mRNA Expression Profiles

Degeneration of sympathetic innervation of the heart occurs in numerous diseases, including diabetes, idiopathic REM sleep disorder, and Parkinson's disease (PD). In PD, cardiac sympathetic denervation occurs in 80-90% of patients and can begin before the onset of motor symptoms. Today, there are no disease-modifying therapies for cardiac sympathetic neurodegeneration, and biomarkers are limited to radioimaging techniques. Analysis of expression levels of coding mRNA and noncoding RNAs, such as microRNAs (miRNAs), can uncover pathways involved in disease, leading to the discovery of biomarkers, pathological mechanisms, and potential drug targets. Whole blood in particular is a clinically relevant source of biomarkers, as blood sampling is inexpensive and simple to perform. Our research group has previously developed a nonhuman primate model of cardiac sympathetic denervation by intravenous administration of the catecholaminergic neurotoxin 6-hydroxydopamine (6-OHDA). In this rhesus macaque (Macaca mulatta) model, imaging with positron emission tomography showed that oral administration of the peroxisome proliferator-activated receptor gamma (PPARγ) agonist pioglitazone (n = 5; 5 mg/kg daily) significantly decreased cardiac inflammation and oxidative stress compared to placebo (n = 5). Here, we report our analysis of miRNA and mRNA expression levels over time in the whole blood of these monkeys. Differential expression of three miRNAs was induced by 6-OHDA (mml-miR-16-2-3p, mml-miR-133d-3p, and mml-miR-1262-5p) and two miRNAs by pioglitazone (mml-miR-204-5p and mml-miR-146b-5p) at 12 weeks posttoxin, while expression of mRNAs involved in inflammatory cytokines and receptors was not significantly affected. Overall, this study contributes to the characterization of rhesus coding and noncoding RNA profiles in normal and disease-like conditions, which may facilitate the identification and clinical translation of biomarkers of cardiac neurodegeneration and neuroprotection.

Serum microRNA in patients undergoing atrial fibrillation ablation

MicroRNAs mediate posttranscriptional gene regulation. The aim of the study was to find a microRNA predictor of successful atrial fibrillation (AF) ablation. A total of 109 patients undergoing first-time AF ablation were included. Nineteen patients were selected to undergo serum microRNA sequencing (study group). The sequencing data were used to select several microRNAs that correlated with 12-month recurrences after AF ablation. Those microRNAs were validated by digital droplet PCR in samples from remaining 90 patients. All patients underwent pulmonary vein isolation (RF ablation, contact force catheter, electroanatomical system). The endpoint of the study was the 12-month AF recurrence rate; the overall recurrence rate was 42.5%. In total, levels of 34 miRNAs were significantly different in sera from patients with AF recurrence compared to patients without AF recurrence. Six microRNAs (miR-183-5p, miR-182-5p, miR-32-5p, miR-107, miR-574-3p, and miR-144-3p) were validated in the whole group. Data from the validation group did not confirm the observations from the study group, as no significant differences were found between miRNAs serum levels in patients with and without recurrences 12 months after AF ablation.

Technological Challenges and Future Issues for the Detection of Circulating MicroRNAs in Patients With Cancer

In the era of precision medicine, the success of clinical trials, notably for patients diagnosed with cancer, strongly relies on biomarkers with pristine clinical value but also on robust and versatile analytical technologies to ensure proper patients' stratification and treatment. In this review, we will first address whether plasmatic and salivary microRNAs can be considered as a reliable source of biomarkers for cancer diagnosis and prognosis. We will then discuss the pre-analytical steps preceding miRNA quantification (from isolation to purification), and how such process could be biased and time-consuming. Next, we will review the most recent tools derived from micro- and nano-technologies for microRNA detection available to date and how they may compete with current standards. This review will prioritize publications using relevant biological samples. The significance of various physical transduction schemes (mechanical, optical, electrical, etc.) for biological detection will be compared, and pros and cons of each method will be widely discussed. Finally, we will debate on how micro and nanotechnologies could widespread the use of biomarkers in modern medicine, to help manage patients with serious diseases such as cancer.

Systematic assessment of commercially available low-input miRNA library preparation kits

High-throughput sequencing is increasingly favoured to assay the presence and abundance of microRNAs (miRNAs) in biological samples, even from low RNA amounts, and a number of commercial vendors now offer kits that allow miRNA sequencing from sub-nanogram (ng) inputs. Although biases introduced during library preparation have been documented, the relative performance of current reagent kits has not been investigated in detail. Here, six commercial kits capable of handling <100ng total RNA input were used for library preparation, performed by kit manufactures, on synthetic miRNAs of known quantities and human total RNA samples. We compared the performance of miRNA detection sensitivity, reliability, titration response and the ability to detect differentially expressed miRNAs. In addition, we assessed the use of unique molecular identifiers (UMI) sequence tags in one kit. We observed differences in detection sensitivity and ability to identify differentially expressed miRNAs between the kits, but none were able to detect the full repertoire of synthetic miRNAs. The reliability within the replicates of all kits was good, while larger differences were observed between the kits, although none could accurately quantify the relative levels of the majority of miRNAs. UMI tags, at least within the input ranges tested, offered little advantage to improve data utility. In conclusion, biases in miRNA abundance are heavily influenced by the kit used for library preparation, suggesting that comparisons of datasets prepared by different procedures should be made with caution. This article is intended to assist researchers select the most appropriate kit for their experimental conditions.

Transcriptomic profiling of cell-free and vesicular microRNAs from matched arterial and venous sera

Extracellular vesicles (EVs) play central physiological and pathophysiological roles in intercellular communication. Biomarker studies addressing disorders such as cardiovascular diseases often focus on circulating microRNAs (miRNAs) and may, depending on the type of disease and clinic routine, utilise patient specimens sampled from arterial or venous blood vessels. Thus, it is essential to test whether circulating miRNA profiles depend on the respective sampling site. We assessed potential differences in arterial and venous cell-free miRNA profiles in a cohort of 20 patients scheduled for cardiac surgery. Prior to surgery, blood was simultaneously sampled from the radial artery and the internal jugular vein. After precipitating crude EVs, we performed small RNA Sequencing, which failed to detect significantly regulated miRNAs using stringent filtering criteria for differential expression analysis. Filtering with less strict criteria, we detected four miRNAs slightly upregulated in arterial samples, one of which could be validated by reverse transcription real-time PCR. The applicability of these findings to purified arterial and venous EVs was subsequently tested in a subset of the initial study population. While an additional clean-up step using size-exclusion chromatography seemed to reduce overall miRNA yield compared to crude EV samples, no miRNAs with differential arteriovenous expression were detected. Unsupervised clustering approaches were unable to correctly classify samples drawn from arteries or veins based on miRNAs in either crude or purified preparations. Particle characterisation of crude preparations as well as characterisation of EV markers in purified EVs resulted in highly similar characteristics for arterial and venous samples. With the exception of specific pathologies (e.g. severe pulmonary disorders), arterial versus venous blood sampling should therefore not represent a likely confounder when studying differentially expressed circulating miRNAs. The use of either arterial or venous serum EV samples should result in highly similar data on miRNA expression profiles for the majority of biomarker studies.

Abbreviations ACE inhibitors: Angiotensin-converting-enzyme inhibitors; ApoA1: Apolipoprotein A1; CNX: Calnexin; Cv: Coefficient of variation; cDNA: Complementary DNA; CABG: Coronary artery bypass graft; DGE: Differential gene expression; DPBS: Dulbecco’s Phosphate Buffered Saline; EVs: Extracellular vesicles; log2FC: Log2 fold change; baseMean: Mean miRNA expression; miRNA: MicroRNA; NTA: Nanoparticle Tracking Analysis; NGS: Next-Generation Sequencing; RT-qPCR: Reverse transcription quantitative real-time PCR; rRNA: Ribosomal RNA; RT: Room temperature; SEC: Size-exclusion chromatography; snoRNA: Small nucleolar RNA; snRNA: Small nuclear RNA; small RNA-Seq: Small RNA Sequencing; SD: Standard deviation; tRNA: Transfer RNA; TEM: Transmission electron microscopy; UA: Uranyl acetate.

Cloudy with a Chance of Insights: Context Dependent Gene Regulation and Implications for Evolutionary Studies

Research in various fields of evolutionary biology has shown that divergence in gene expression is a key driver for phenotypic evolution. An exceptional contribution of cis-regulatory divergence has been found to contribute to morphological diversification. In the light of these findings, the analysis of genome-wide expression data has become one of the central tools to link genotype and phenotype information on a more mechanistic level. However, in many studies, especially if general conclusions are drawn from such data, a key feature of gene regulation is often neglected. With our article, we want to raise awareness that gene regulation and thus gene expression is highly context dependent. Genes show tissue- and stage-specific expression. We argue that the regulatory context must be considered in comparative expression studies.

A comparison of RNA extraction and sequencing protocols for detection of small RNAs in plasma

BACKGROUND

Circulating microRNAs (miRNAs) are attractive non-invasive biomarkers for a variety of conditions due to their stability and altered pathophysiological expression levels. Reliable detection of global expression profiles is required to maximise miRNA biomarker discovery. Although developments in small RNA-Seq technology have improved detection of plasma-based miRNAs, the low RNA content and sequencing bias introduced during library preparation remain challenging. In this study we compare commercially available RNA extraction methods using MagnaZol (Bioo Scientific) or miRNeasy (QIAGEN) and three library preparation methods - CleanTag (TriLink), NEXTflex (Bioo Scientific) and QIAseq (QIAGEN) - which aim to address one or both of these issues.

RESULTS

Different RNA extractions and library preparation protocols result in differential detection of miRNAs. A greater proportion of reads mapped to miRNAs in libraries prepared with MagnaZol RNA than with miRNeasy RNA. Libraries prepared using QIAseq demonstrated the greatest miRNA diversity with many more very low abundance miRNAs detected (~ 2-3 fold more with < 10 reads), whilst CleanTag detected the fewest individual miRNAs and considerably over-represented miR-486-5p. Libraries prepared with QIAseq had the strongest correlation with RT-qPCR quantification. Analysis of unique molecular indices (UMIs) incorporated in the QIAseq protocol indicate that little PCR bias is introduced during small RNA library preparation.

CONCLUSIONS

Small RNAs were consistently detected using all RNA extraction and library preparation protocols tested, but with some miRNAs at significantly different levels. Choice of the most suitable protocol should be informed by the relative importance of minimising the total sequencing required, detection of rare miRNAs or absolute quantification.