Comprehensive assessment of multiple biases in small RNA sequencing reveals significant differences in the performance of widely used methods

Characteristics of MicroRNA Expression Depending on the Presence or Absence of Meningioma in Patients with Neurofibromatosis Type 2: A Secondary Analysis

Meningioma is the second most frequent tumor in patients with neurofibromatosis type 2 (NF2). The presence of meningioma is believed to be a negative prognostic marker in these patients. However, the molecular mechanisms involved in the tumorigenesis of NF2-associated meningioma are not well characterized. Epigenetic regulation, including microRNAs (miRNAs), may be involved in the development of different tumor types in patients with NF2. The objective of this study is to explore the different characteristics of serum miRNA expression depending on the presence or absence of meningioma in patients with NF2. Nine patients with NF2 who were treated at the Department of Neurosurgery, Hiroshima University Hospital, were included. Total RNA (including small RNAs) was extracted from serum samples for the preparation of a small RNA library for next-generation sequencing analysis. Differentially expressed miRNAs (DEMs) were analyzed using the DESeq2 package to compare the characteristic miRNA expression profiles of patients with and without meningioma. In small RNA sequencing analysis, out of a total of 1,879 miRNAs registered in the database, the expressions of 657 miRNAs were observed. In DEM analysis, the expressions of four miRNAs, namely, hsa-miR-664b, hsa-miR-7706, hsa-miR-590, and hsa-miR-6513, were downregulated in patients with NF2 with meningioma compared with patients with NF2 without meningioma. Hsa-miR-193a was identified as the only upregulated miRNA in patients with NF2 with meningioma. In conclusion, we identified different circulating miRNA expression characteristics depending on the presence or absence of meningioma in patients with NF2.

Recommendations for detection, validation, and evaluation of RNA editing events in cardiovascular and neurological/neurodegenerative diseases

RNA editing, a common and potentially highly functional form of RNA modification, encompasses two different RNA modifications, namely adenosine to inosine (A-to-I) and cytidine to uridine (C-to-U) editing. As inosines are interpreted as guanosines by the cellular machinery, both A-to-I and C-to-U editing change the nucleotide sequence of the RNA. Editing events in coding sequences have the potential to change the amino acid sequence of proteins, whereas editing events in noncoding RNAs can, for example, affect microRNA target binding. With advancing RNA sequencing technology, more RNA editing events are being discovered, studied, and reported. However, RNA editing events are still often overlooked or discarded as sequence read quality defects. With this position paper, we aim to provide guidelines and recommendations for the detection, validation, and follow-up experiments to study RNA editing, taking examples from the fields of cardiovascular and brain disease. We discuss all steps, from sample collection, storage, and preparation, to different strategies for RNA sequencing and editing-sensitive data analysis strategies, to validation and follow-up experiments, as well as potential pitfalls and gaps in the available technologies. This paper may be used as an experimental guideline for RNA editing studies in any disease context.

Highly Multiplexed, Efficient, and Automated Single-Cell MicroRNA Sequencing with Digital Microfluidics

Single-cell microRNA (miRNA) sequencing has allowed for comprehensively studying the abundance and complex networks of miRNAs, which provides insights beyond single-cell heterogeneity into the dynamic regulation of cellular events. Current benchtop-based technologies for single-cell miRNA sequencing are low throughput, limited reaction efficiency, tedious manual operations, and high reagent costs. Here, a highly multiplexed, efficient, integrated, and automated sample preparation platform is introduced for single-cell miRNA sequencing based on digital microfluidics (DMF), named Hiper-seq. The platform integrates major steps and automates the iterative operations of miRNA sequencing library construction by digital control of addressable droplets on the DMF chip. Based on the design of hydrophilic micro-structures and the capability of handling droplets of DMF, multiple single cells can be selectively isolated and subject to sample processing in a highly parallel way, thus increasing the throughput and efficiency for single-cell miRNA measurement. The nanoliter reaction volume of this platform enables a much higher miRNA detection ability and lower reagent cost compared to benchtop methods. It is further applied Hiper-seq to explore miRNAs involved in the ossification of mouse skeletal stem cells after bone fracture and discovered unreported miRNAs that regulate bone repairing.

Mesenchymal stromal cell derived extracellular vesicles as a therapeutic tool: immune regulation, MSC priming, and applications to SLE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by a dysfunction of the immune system. Mesenchymal stromal cell (MSCs) derived extracellular vesicles (EVs) are nanometer-sized particles carrying a diverse range of bioactive molecules, such as proteins, miRNAs, and lipids. Despite the methodological disparities, recent works on MSC-EVs have highlighted their broad immunosuppressive effect, thus driving forwards the potential of MSC-EVs in the treatment of chronic diseases. Nonetheless, their mechanism of action is still unclear, and better understanding is needed for clinical application. Therefore, we describe in this review the diverse range of bioactive molecules mediating their immunomodulatory effect, the techniques and possibilities for enhancing their immune activity, and finally the potential application to SLE.

Isothermal amplification methods in cancer-related miRNA detection; a new paradigm in study of cancer pathology

MicroRNAs (miRNAs) are short, non-coding RNA molecules that regulate gene expression. They are involved in a wide range of biological processes, including development, differentiation, cell cycle regulation, and response to stress. Numerous studies have demonstrated that miRNAs are present in different bodily fluids, which could serve as an important biomarker. The advancement of techniques and strategies for the identification of cancer-associated miRNAs in human specimens offers a novel opportunity to diagnose cancer in early stages, predict patient prognosis and evaluate response to treatment. Isothermal techniques including loop-mediated isothermal amplification (LAMP), rolling circle amplification (RCA), or recombinase polymerase amplification (RPA) offer simplicity, efficiency, and rapidity in miRNA detection processes. In contrast to traditional PCR (polymerase chain reaction), these techniques analysis and quantify miRNA molecules in specimens using a single constant temperature. In this comprehensive review, we summarized the recent advances in cancer-related miRNA detection via highly sensitive isothermal amplification methods by more focusing on the involved mechanism.

Single Cell Micro RNA Sequencing Library Preparation

Micro RNAs represent important post-transcriptional regulators in health and are involved in the onset of many diseases. Therefore, the further characterization of physiological miRNA functions is an important basic research question, and miRNAs even have high potential as biomarkers both for prognosis and diagnosis. In order to exploit this potential, it is mandatory to accurately quantify the miRNA expression not only in bulk but also on the single-cell level. Here, we describe a protocol, which facilitates miRNA sequencing library preparation of very low input samples, single cells, and even clinical samples such as circulating tumor cells. The protocol can be combined with different single-cell isolation methods (e.g., micromanipulation and FACS sorting). After cell lysis, sequencing adapters are ligated to the miRNAs, other ncRNA species, and adapter dimers are reduced by exonuclease digest, the miRNA library is reverse transcribed, amplified, and purified. Furthermore, quality controls are described to select only high-quality samples for sequencing.

The Typical tRNA Co-Expresses Multiple 5' tRNA Halves Whose Sequences and Abundances Depend on Isodecoder and Isoacceptor and Change with Tissue Type, Cell Type, and Disease

Transfer RNA-derived fragments (tRFs) are noncoding RNAs that arise from either mature transfer RNAs (tRNAs) or their precursors. One important category of tRFs comprises the tRNA halves, which are generated through cleavage at the anticodon. A given tRNA typically gives rise to several co-expressed 5'-tRNA halves (5'-tRHs) that differ in the location of their 3' ends. These 5'-tRHs, even though distinct, have traditionally been treated as indistinguishable from one another due to their near-identical sequences and lengths. We focused on co-expressed 5'-tRHs that arise from the same tRNA and systematically examined their exact sequences and abundances across 10 different human tissues. To this end, we manually curated and analyzed several hundred human RNA-seq datasets from NCBI's Sequence Run Archive (SRA). We grouped datasets from the same tissue into their own collection and examined each group separately. We found that a given tRNA produces different groups of co-expressed 5'-tRHs in different tissues, different cell lines, and different diseases. Importantly, the co-expressed 5'-tRHs differ in their sequences, absolute abundances, and relative abundances, even among tRNAs with near-identical sequences from the same isodecoder or isoacceptor group. The findings suggest that co-expressed 5'-tRHs that are produced from the same tRNA or closely related tRNAs have distinct, context-dependent roles. Moreover, our analyses show that cell lines modeling the same tissue type and disease may not be interchangeable when it comes to experimenting with tRFs.

Preparation and processing of dried blood spots for microRNA sequencing

Dried blood spots (DBS) are biological samples commonly collected from newborns and in geographic areas distanced from laboratory settings for the purposes of disease testing and identification. MicroRNAs (miRNAs)-small non-coding RNAs that regulate gene activity at the post-transcriptional level-are emerging as critical markers and mediators of disease, including cancer, infectious diseases, and mental disorders. This protocol describes optimized procedural steps for utilizing DBS as a reliable source of biological material for obtaining peripheral miRNA expression profiles. We outline key practices, such as the method of DBS rehydration that maximizes RNA extraction yield, and the use of degenerate oligonucleotide adapters to mitigate ligase-dependent biases that are associated with small RNA sequencing. The standardization of miRNA readout from DBS offers numerous benefits: cost-effectiveness in sample collection and processing, enhanced reliability and consistency of miRNA profiling, and minimal invasiveness that facilitates repeated testing and retention of participants. The use of DBS-based miRNA sequencing is a promising method to investigate disease mechanisms and to advance personalized medicine.

Urinary phthalate metabolites and small non-coding RNAs from seminal plasma extracellular vesicles among men undergoing infertility treatment

Non-coding RNA (ncRNA) cargo of extracellular vesicles (EVs) in the male reproductive tract play critical roles in semen quality and emerging evidence suggests their susceptibility to environmental factors. Male phthalate exposures have been linked to poor semen quality, sperm DNA methylation profiles and embryo development; however, there is limited evidence on their potential impact on EV ncRNAs profiles. We evaluated the association between urinary phthalate metabolites and small ncRNAs (sncRNAs) of seminal plasma EVs (spEV) among men receiving clinical infertility care. We conducted sncRNA sequencing of EVs in 96 seminal plasma samples collected from the Sperm Environmental Epigenetics and Development Study (SEEDS). Sequencing reads were mapped to human transcriptome databases using STAR. Urinary metabolite concentrations of thirteen phthalates and two DiNCH, a phthalate alternative, were measured via tandem mass spectrometry. Associations with normalized counts were assessed using EdgeR (FDR<0.05) adjusting for urinary dilution via specific gravity, age, BMI, batch, and biotype-specific total counts. Select metabolites, MEOHP, MECPP, ∑DEHP, MCPP, MCNP, MCOP, were negatively (p < 0.05) correlated with miRNA relative abundance. Similarly, nine metabolites including MEOHP, MECPP, MEHP, MCPP, MHBP, MHiNCH, MiBP, MEHHP, MCOP and ∑DEHP were associated (q < 0.05) with normalized counts from 23 unique ncRNA transcripts (7 miRNAs (pre & mature); 6 tRFs; and 10 piRNAs), most (78%) of which displayed increased expression patterns. miRNA and tRFs gene targets were enriched in vesicle-mediated transport and developmental-related ontology terms, such as tyrosine kinase, head development, and cell morphogenesis. Six genes (MAPK1, BMPR1A/2, PTEN, TGFBR2, TP53 and APP) were present in all the ontology terms and predicted to form protein association networks. piRNAs were annotated to pseudogenes of genes important in EV cargo transfer and embryonic development. This is the first study to associate phthalate exposures to altered spEV sncRNA profiles. Future studies are needed to determine their impact on reproductive outcomes.

A systematic review of non-coding RNA genes with differential expression profiles associated with autism spectrum disorders

AIMS

To identify differential expression of shorter non-coding RNA (ncRNA) genes associated with autism spectrum disorders (ASD).

BACKGROUND

ncRNA are functional molecules that derive from non-translated DNA sequence. The HUGO Gene Nomenclature Committee (HGNC) have approved ncRNA gene classes with alignment to the reference human genome. One subset is microRNA (miRNA), which are highly conserved, short RNA molecules that regulate gene expression by direct post-transcriptional repression of messenger RNA. Several miRNA genes are implicated in the development and regulation of the nervous system. Expression of miRNA genes in ASD cohorts have been examined by multiple research groups. Other shorter classes of ncRNA have been examined less. A comprehensive systematic review examining expression of shorter ncRNA gene classes in ASD is timely to inform the direction of research.

METHODS

We extracted data from studies examining ncRNA gene expression in ASD compared with non-ASD controls. We included studies on miRNA, piwi-interacting RNA (piRNA), small NF90 (ILF3) associated RNA (snaR), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), transfer RNA (tRNA), vault RNA (vtRNA) and Y RNA. The following electronic databases were searched: Cochrane Library, EMBASE, PubMed, Web of Science, PsycINFO, ERIC, AMED and CINAHL for papers published from January 2000 to May 2022. Studies were screened by two independent investigators with a third resolving discrepancies. Data was extracted from eligible papers.

RESULTS

Forty-eight eligible studies were included in our systematic review with the majority examining miRNA gene expression alone. Sixty-four miRNA genes had differential expression in ASD compared to controls as reported in two or more studies, but often in opposing directions. Four miRNA genes had differential expression in the same direction in the same tissue type in at least 3 separate studies. Increased expression was reported in miR-106b-5p, miR-155-5p and miR-146a-5p in blood, post-mortem brain, and across several tissue types, respectively. Decreased expression was reported in miR-328-3p in bloods samples. Seven studies examined differential expression from other classes of ncRNA, including piRNA, snRNA, snoRNA and Y RNA. No individual ncRNA genes were reported in more than one study. Six studies reported differentially expressed snoRNA genes in ASD. A meta-analysis was not possible because of inconsistent methodologies, disparate tissue types examined, and varying forms of data presented.

CONCLUSION

There is limited but promising evidence associating the expression of certain miRNA genes and ASD, although the studies are of variable methodological quality and the results are largely inconsistent. There is emerging evidence associating differential expression of snoRNA genes in ASD. It is not currently possible to say whether the reports of differential expression in ncRNA may relate to ASD aetiology, a response to shared environmental factors linked to ASD such as sleep and nutrition, other molecular functions, human diversity, or chance findings. To improve our understanding of any potential association, we recommend improved and standardised methodologies and reporting of raw data. Further high-quality research is required to shine a light on possible associations, which may yet yield important information.

Ultrasensitive and multiplexed miRNA detection system with DNA-PAINT

MiRNAs hold great potential as biomarkers for the early detection and monitoring of diseases based on their differential expression profiles. Therefore, the sensitive, specific and accurate detection of miRNAs represents an emerging new tool to improve diagnosis and treatment of several diseases, cancer in particular. DNA origami-based miRNA detection is particularly advantageous as it allows to incorporate multiple attachment sites to capture different target miRNAs at the nanoscale. In this work, we present a DNA origami nanoarray system providing distance-dependent recognition of miRNAs by applying super-resolution microscopy technique; DNA-PAINT (point accumulation for imaging in nanoscale topography). The sensor can detect up to 4 miRNAs either separately or in combination based on the relative distance to the boundary markers on the structure using a single imager strand. The detection is highly sensitive, with a limit of detection down to the low femtomolar range (11 fM - 388 fM) and has a large dynamic range up to 10 nM without need for amplification. Moreover, our detection system can discriminate single base mismatches with low false positive rates. Using our strategy, we demonstrate the detection of endogenous miRNAs from cell extracts of cancer cell lines and plasma from breast cancer patients. Overall, we developed an ultrasensitive and amplification-free, DNA-PAINT imaging-based miRNA detection method using DNA origami nanoarray system for the detection of breast-cancer associated miRNAs which potentially provides a sensitive and accurate alternative to the current multiplexed diagnostic technologies.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.

Epigenetic Mechanisms and Nephrotic Syndrome: A Systematic Review

A small subset of people with nephrotic syndrome (NS) have genetically driven disease. However, the disease mechanisms for the remaining majority are unknown. Epigenetic marks are reversible but stable regulators of gene expression with utility as biomarkers and therapeutic targets. We aimed to identify and assess all published human studies of epigenetic mechanisms in NS. PubMed (MEDLINE) and Embase were searched for original research articles examining any epigenetic mechanism in samples collected from people with steroid resistant NS, steroid sensitive NS, focal segmental glomerulosclerosis or minimal change disease. Study quality was assessed by using the Joanna Briggs Institute critical appraisal tools. Forty-nine studies met our inclusion criteria. The majority of these examined micro-RNAs (n = 35, 71%). Study quality was low, with only 23 deemed higher quality, and most of these included fewer than 100 patients and failed to validate findings in a second cohort. However, there were some promising concordant results between the studies; higher levels of serum miR-191 and miR-30c, and urinary miR-23b-3p and miR-30a-5p were observed in NS compared to controls. We have identified that the epigenome, particularly DNA methylation and histone modifications, has been understudied in NS. Large clinical studies, which utilise the latest high-throughput technologies and analytical pipelines, should focus on addressing this critical gap in the literature.

Transcriptional and Epigenetic Regulation of Context-Dependent Plasticity in T-Helper Lineages

Th cell lineage determination and functional specialization are tightly linked to the activation of lineage-determining transcription factors (TFs) that bind cis-regulatory elements. These lineage-determining TFs act in concert with multiple layers of transcriptional regulators to alter the epigenetic landscape, including DNA methylation, histone modification and three-dimensional chromosome architecture, in order to facilitate the specific Th gene expression programs that allow for phenotypic diversification. Accumulating evidence indicates that Th cell differentiation is not as rigid as classically held; rather, extensive phenotypic plasticity is an inherent feature of T cell lineages. Recent studies have begun to uncover the epigenetic programs that mechanistically govern T cell subset specification and immunological memory. Advances in next generation sequencing technologies have allowed global transcriptomic and epigenomic interrogation of CD4+ Th cells that extends previous findings focusing on individual loci. In this review, we provide an overview of recent genome-wide insights into the transcriptional and epigenetic regulation of CD4+ T cell-mediated adaptive immunity and discuss the implications for disease as well as immunotherapies.

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.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.

Embryonic microRNAs are essential for bovine preimplantation embryo development

MicroRNAs (miRNAs) are small, noncoding RNAs that regulate gene expression after transcription. miRNAs are present in transcriptionally quiescent full-grown oocytes and preimplantation embryos that display a low level of transcription prior to embryonic genome activation. The role of miRNAs, if any, in preimplantation development is not known. The temporal pattern of expression of miRNAs during bovine preimplantation development was determined by small RNA-sequencing using eggs and preimplantation embryos (1-cell, 2-cell, 4-cell, 8-cell, 16-cell, morula, and blastocyst). Embryos cultured in the presence of α-amanitin, which permitted the distinguishing of maternal miRNAs from embryonic miRNAs, indicated that embryonic miRNA expression was first detected at the two-cell stage but dramatically increased during the morula and blastocyst stages. Targeting DGCR8 by a small-interfering RNA/morpholino approach revealed a role for miRNAs in the morula-to-blastocyst transition. Knockdown of DGCR8 not only inhibited expression of embryonically expressed miRNAs but also inhibited the morula-to-blastocyst transition. In addition, RNA-sequencing identified an increased relative abundance of messenger RNAs potentially targeted by embryonic miRNAs in DGCR8-knockdown embryos when compared with controls. Results from these experiments implicate an essential role for miRNAs in bovine preimplantation embryo development.

Circulating microRNAs in cellular and antibody-mediated heart transplant rejection

BACKGROUND

Noninvasive monitoring of heart allograft health is important to improve clinical outcomes. MicroRNAs (miRs) are promising biomarkers of cardiovascular disease and limited studies suggest they can be used to noninvasively diagnose acute heart transplant rejection.

METHODS

The Genomic Research Alliance for Transplantation (GRAfT) is a multicenter prospective cohort study that phenotyped heart transplant patients from 5 mid-Atlantic centers. Patients who had no history of rejection after transplant were compared to patients with acute cellular rejection (ACR) or antibody-mediated rejection (AMR). Small RNA sequencing was performed on plasma samples collected at the time of an endomyocardial biopsy. Differential miR expression was performed with adjustment for clinical covariates. Regression was used to develop miR panels with high diagnostic accuracy for ACR and AMR. These panels were then validated in independent samples from GRAfT and Stanford University. Receiver operating characteristic curves were generated and area under the curve (AUC) statistics calculated. Distinct ACR and AMR clinical scores were developed to translate miR expression data for clinical use.

RESULTS

The GRAfT cohort had a median age of 52 years, with 35% females and 45% Black patients. Between GRAfT and Stanford, we included 157 heart transplant patients: 108 controls and 49 with rejection (50 ACR and 38 AMR episodes). After differential miR expression and regression analysis, we identified 12 miRs that accurately discriminate ACR and 17 miRs in AMR. Independent validation of the miR panels within GRAfT led to an ACR AUC 0.92 (95% confidence interval [CI]: 0.86-0.98) and AMR AUC 0.82 (95% CI: 0.74-0.90). The externally validated ACR AUC was 0.72 (95% CI: 0.59-0.82). We developed distinct ACR and AMR miR clinical scores (range 0-100), a score ≥ 65, identified ACR with 86% sensitivity, 76% specificity, and 98% negative predictive value, for AMR score performance was 82%, 84% and 97%, respectively.

CONCLUSIONS

We identified novel miRs that had excellent performance to noninvasively diagnose acute rejection after heart transplantation. Once rigorously validated, the unique clinical ACR and AMR scores usher in an era whereby genomic biomarkers can be used to screen and diagnose the subtype of rejection. These novel biomarkers may potentially alleviate the need for an endomyocardial biopsy while facilitating the initiation of targeted therapy based on the noninvasive diagnosis of ACR or AMR.

A curated human cellular microRNAome based on 196 primary cell types

Background

An incomplete picture of the expression distribution of microRNAs (miRNAs) across human cell types has long hindered our understanding of this important regulatory class of RNA. With the continued increase in available public small RNA sequencing datasets, there is an opportunity to more fully understand the general distribution of miRNAs at the cell level.

Results

From the NCBI Sequence Read Archive, we obtained 6,054 human primary cell datasets and processed 4,184 of them through the miRge3.0 small RNA sequencing alignment software. This dataset was curated down, through shared miRNA expression patterns, to 2,077 samples from 196 unique cell types derived from 175 separate studies. Of 2,731 putative miRNAs listed in miRBase (v22.1), 2,452 (89.8%) were detected. Among reasonably expressed miRNAs, 108 were designated as cell specific/near specific, 59 as infrequent, 52 as frequent, 54 as near ubiquitous, and 50 as ubiquitous. The complexity of cellular microRNA expression estimates recapitulates tissue expression patterns and informs on the miRNA composition of plasma.

Conclusions

This study represents the most complete reference, to date, of miRNA expression patterns by primary cell type. The data are available through the human cellular microRNAome track at the UCSC Genome Browser (https://genome.ucsc.edu/cgi-bin/hgHubConnect) and an R/Bioconductor package (https://bioconductor.org/packages/microRNAome/).

Characterization of presence and activity of microRNAs in the rumen of cattle hints at possible host-microbiota cross-talk mechanism

MicroRNAs (miRNAs), as important post-transcriptional regulators, are ubiquitous in various tissues. The aim of this exploratory study was to determine the presence of miRNAs in rumen fluid, and to investigate the possibility of miRNA-mediated cross-talk within the ruminal ecosystem. Rumen fluid samples from four cannulated Holstein cows were collected during two feeding regimes (forage and high-grain diet) and DNA and RNA were extracted for amplicon and small RNA sequencing. Epithelial biopsies were simultaneously collected to investigate the co-expression of miRNAs in papillae and rumen fluid. We identified 377 miRNAs in rumen fluid and 638 in rumen papillae, of which 373 were shared. Analysis of microbiota revealed 20 genera to be differentially abundant between the two feeding regimes, whereas no difference in miRNAs expression was detected. Correlations with at least one genus were found for 170 miRNAs, of which, 39 were highly significant (r > |0.7| and P < 0.01). Both hierarchical clustering of the correlation matrix and WGCNA analysis identified two main miRNA groups. Putative target and functional prediction analysis for the two groups revealed shared pathways with the predicted metabolic activities of the microbiota. Hence, our study supports the hypothesis of a cross-talk within the rumen at least partly mediated by miRNAs.

Implementation of GA-VirReport, a Web-Based Bioinformatics Toolkit for Post-Entry Quarantine Screening of Virus and Viroids in Plants

High-throughput sequencing (HTS) of host plant small RNA (sRNA) is a popular approach for plant virus and viroid detection. The major bottlenecks for implementing this approach in routine virus screening of plants in quarantine include lack of computational resources and/or expertise in command-line environments and limited availability of curated plant virus and viroid databases. We developed: (1) virus and viroid report web-based bioinformatics workflows on Galaxy Australia called GA-VirReport and GA-VirReport-Stats for detecting viruses and viroids from host plant sRNA extracts and (2) a curated higher plant virus and viroid database (PVirDB). We implemented sRNA sequencing with unique dual indexing on a set of plants with known viruses. Sequencing data were analyzed using GA-VirReport and PVirDB to validate these resources. We detected all known viruses in this pilot study with no cross-sample contamination. We then conducted a large-scale diagnosis of 105 imported plants processed at the post-entry quarantine facility (PEQ), Australia. We detected various pathogens in 14 imported plants and discovered that de novo assembly using 21–22 nt sRNA fraction and the megablast algorithm yielded better sensitivity and specificity. This study reports the successful, large-scale implementation of HTS and a user-friendly bioinformatics workflow for virus and viroid screening of imported plants at the PEQ.

An Expanded Landscape of Unusually Short RNAs in 11 Samples from Six Eukaryotic Organisms

Small RNA sequencing (sRNA-Seq) approaches unveiled sequences derived from longer non-coding RNAs, such as transfer RNA (tRNA) and ribosomal RNA (rRNA) fragments, known as tRFs and rRFs, respectively. However, rRNAs and RNAs shorter than 16 nt are often depleted from library preparations/sequencing analyses, although they may be functional. Here, we sought to obtain a complete repertoire of small RNAs by sequencing the total RNA from 11 samples of 6 different eukaryotic organisms, from yeasts to human, in an extended 8- to 30-nt window of RNA length. The 8- to 15-nt window essentially contained fragments of longer non-coding RNAs, such as microRNAs, PIWI-associated RNAs (piRNAs), small nucleolar RNAs (snoRNAs), tRNAs and rRNAs. Notably, unusually short RNAs < 16 nt were more abundant than those >16 nt in bilaterian organisms. A new RT-qPCR method confirmed that two unusually short rRFs of 12 and 13 nt were more overly abundant (~3-log difference) than two microRNAs. We propose to not deplete rRNA and to reduce the lower threshold of RNA length to include unusually short RNAs in sRNA-Seq analyses and datasets, as their abundance and diversity support their potential role and importance as biomarkers of disease and/or mediators of cellular function.

Small RNA-Sequencing for Analysis of Circulating miRNAs: Benchmark Study

Small RNA-sequencing (RNA-Seq) is being increasingly used for profiling of circulating microRNAs (miRNAs), a new group of promising biomarkers. Unfortunately, small RNA-Seq protocols are prone to biases limiting quantification accuracy, which motivated development of several novel methods. Here, we present comparison of all small RNA-Seq library preparation approaches that are commercially available for quantification of miRNAs in biofluids. Using synthetic and human plasma samples, we compared performance of traditional two-adaptor ligation protocols (Lexogen, Norgen), as well as methods using randomized adaptors (NEXTflex), polyadenylation (SMARTer), circularization (RealSeq), capture probes (EdgeSeq), or unique molecular identifiers (QIAseq). There was no single protocol outperforming others across all metrics. Limited overlap of measured miRNA profiles was documented between methods largely owing to protocol-specific biases. Methods designed to minimize bias largely differ in their performance, and contributing factors were identified. Usage of unique molecular identifiers has rather negligible effect and, if designed incorrectly, can even introduce spurious results. Together, these results identify strengths and weaknesses of all current methods and provide guidelines for applications of small RNA-Seq in biomarker research.

Exosomal-microRNA transcriptome profiling of Parental and CSC-like MDA-MB-231 cells in response to cisplatin treatment

Triple negative breast cancer (TNBC) is the most aggressive breast cancer subtype with higher risk of metastasis and cancer reoccurrence. Cisplatin is one of the potential anticancer drugs for treating TNBC, where its effectiveness remains challenged by frequent occurrence of cisplatin resistance. Since acquirement of drug resistance often being associated with presence of cancer stem cells (CSCs), investigation has been conducted, suggesting CSC-like subpopulation to be more resistant to cisplatin than their parental counterpart. On the other hand, plethora evidences showed the transmission of exosomal-miRNAs are capable of promoting drug resistance in breast cancers. In this study, we aim to elucidate the differential expression of exosomal-microRNAs profile and reveal the potential target genes in correlation to cisplatin resistance associated with CSC-like subpopulation by using TNBC cell line (MDA-MB-231). Utilizing next generation sequencing and Nanostring techniques, cisplatin-induced dysregulation of exosomal-miRNAs were evaluated in maximal for CSC-like subpopulation as compared to parental cells. Intriguingly, more oncogenic exosomal-miRNAs profile was detected from treated CSC-like subpopulation, which may correlate to enhancement of drug resistance and maintenance of CSCs. In treated CSC-like subpopulation, unique clusters of exosomal-miRNAs namely miR-221-3p, miR-196a-5p, miR-17-5p and miR-126-3p were predicted to target on six genes (ATXN1, LATS1, GSK3β, ITGA6, JAG1 and MYC), aligned with previous finding which demonstrated dysregulation of these genes in treated CSC-like subpopulation. Our results highlight the potential correlation of exosomal-miRNAs and their target genes as well as novel perspectives of the corresponding pathways that may be essential to contribute to the attenuated cytotoxicity of cisplatin in CSC-like subpopulation.

A comparative analysis of single cell small RNA sequencing data reveals heterogeneous isomiR expression and regulation

MicroRNAs (miRNAs) are non-coding small RNAs which play a critical role in the regulation of gene expression in cells. It is known that miRNAs are often expressed as multiple isoforms, called isomiRs, which may have alternative regulatory functions. Despite the recent development of several single cell small RNA sequencing protocols, these methods have not been leveraged to investigate isomiR expression and regulation to better understand their role on a single cell level. Here we integrate sequencing data from three independent studies and find substantial differences in isomiR composition that suggest that cell autonomous mechanisms may drive isomiR processing. We also find evidence of altered regulatory functions of different classes of isomiRs, when compared to their respective wild-type miRNA, which supports a biological role for many of the isomiRs that are expressed.

A MicroRNA Next-Generation-Sequencing Discovery Assay (miND) for Genome-Scale Analysis and Absolute Quantitation of Circulating MicroRNA Biomarkers

The plasma levels of tissue-specific microRNAs can be used as diagnostic, disease severity and prognostic biomarkers for chronic and acute diseases and drug-induced injury. Thereby, the combination of diverse microRNAs into biomarker signatures using multivariate statistics seems especially powerful from the perspective of tissue and condition specific microRNA shedding into the plasma. Although next-generation sequencing (NGS) technology enables one to analyse circulating microRNAs on a genome-scale level, it suffers from potential biases (e.g., adapter ligation bias) and lacks absolute transcript quantitation as well as tailor-made quality controls. In order to develop a robust NGS discovery assay for genome-scale quantitation of circulating microRNAs, we first evaluated the sensitivity, repeatability and ligation bias of four commercially available small RNA library preparation protocols. The protocol from RealSeq Biosciences was selected based on its performance and usability and coupled with a novel panel of exogenous small RNA spike-in controls to enable quality control and absolute quantitation, thus ensuring comparability of data across independent NGS experiments. The established microRNA Next-Generation-Sequencing Discovery Assay (miND) was validated for its relative accuracy, precision, analytical measurement range and sequencing bias and was considered fit-for-purpose for microRNA biomarker discovery. Summarized, all these criteria were met, and thus, our analytical platform is considered fit-for-purpose for microRNA biomarker discovery from biofluids in the setting of any diagnostic, prognostic or patient stratification need. The established miND assay was tested on serum, cerebrospinal fluid (CSF), synovial fluid (SF) and extracellular vesicles (EV) extracted from cell culture medium of primary cells and proved its potential to be used across different sample types.

Quality control and evaluation of plant epigenomics data

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

DNA Methylation and Non-Coding RNAs during Tissue-Injury Associated Pain

While about half of the population experience persistent pain associated with tissue damages during their lifetime, current symptom-based approaches often fail to reduce such pain to a satisfactory level. To provide better patient care, mechanism-based analgesic approaches must be developed, which necessitates a comprehensive understanding of the nociceptive mechanism leading to tissue injury-associated persistent pain. Epigenetic events leading the altered transcription in the nervous system are pivotal in the maintenance of pain in tissue injury. However, the mechanisms through which those events contribute to the persistence of pain are not fully understood. This review provides a summary and critical evaluation of two epigenetic mechanisms, DNA methylation and non-coding RNA expression, on transcriptional modulation in nociceptive pathways during the development of tissue injury-associated pain. We assess the pre-clinical data and their translational implication and evaluate the potential of controlling DNA methylation and non-coding RNA expression as novel analgesic approaches and/or biomarkers of persistent pain.

Phospho-RNA sequencing with circAID-p-seq

Most RNA footprinting approaches that require ribonuclease cleavage generate RNA fragments bearing a phosphate or cyclic phosphate group at their 3′ end. Unfortunately, current library preparation protocols rely only on a 3′ hydroxyl group for adaptor ligation or poly-A tailing. Here, we developed circAID-p-seq, a PCR-free library preparation for selective 3′ phospho-RNA sequencing. As a proof of concept, we applied circAID-p-seq to ribosome profiling, which is based on sequencing of RNA fragments protected by ribosomes after endonuclease digestion. CircAID-p-seq, combined with the dedicated computational pipeline circAidMe, facilitates accurate, fast and highly efficient sequencing of phospho-RNA fragments from eukaryotic cells and tissues. We used circAID-p-seq to portray ribosome occupancy in transcripts, providing a versatile and PCR-free strategy to possibly unravel any endogenous 3′-phospho RNA molecules.

Role of microRNAs in Disorders of Gut-Brain Interactions: Clinical Insights and Therapeutic Alternatives

Disorders of gut–brain interactions (DGBIs) are heterogeneous in nature and intertwine with diverse pathophysiological mechanisms. Regular functioning of the gut requires complex coordinated interplay between a variety of gastrointestinal (GI) cell types and their functions are regulated by multiple mechanisms at the transcriptional, post-transcriptional, translational, and post-translational levels. MicroRNAs (miRNAs) are small non-coding RNA molecules that post-transcriptionally regulate gene expression by binding to specific mRNA targets to repress their translation and/or promote the target mRNA degradation. Dysregulation of miRNAs might impair gut physiological functions leading to DGBIs and gut motility disorders. Studies have shown miRNAs regulate gut functions such as visceral sensation, gut immune response, GI barrier function, enteric neuronal development, and GI motility. These biological processes are highly relevant to the gut where neuroimmune interactions are key contributors in controlling gut homeostasis and functional defects lead to DGBIs. Although extensive research has explored the pathophysiology of DGBIs, further research is warranted to bolster the molecular mechanisms behind these disorders. The therapeutic targeting of miRNAs represents an attractive approach for the treatment of DGBIs because they offer new insights into disease mechanisms and have great potential to be used in the clinic as diagnostic markers and therapeutic targets. Here, we review recent advances regarding the regulation of miRNAs in GI pacemaking cells, immune cells, and enteric neurons modulating pathophysiological mechanisms of DGBIs. This review aims to assess the impacts of miRNAs on the pathophysiological mechanisms of DGBIs, including GI dysmotility, impaired intestinal barrier function, gut immune dysfunction, and visceral hypersensitivity. We also summarize the therapeutic alternatives for gut microbial dysbiosis in DGBIs, highlighting the clinical insights and areas for further exploration. We further discuss the challenges in miRNA therapeutics and promising emerging approaches.

Chemical capping improves template switching and enhances sequencing of small RNAs

Template-switching reverse transcription is widely used in RNA sequencing for low-input and low-quality samples, including RNA from single cells or formalin-fixed paraffin-embedded (FFPE) tissues. Previously, we identified the native eukaryotic mRNA 5′ cap as a key structural element for enhancing template switching efficiency. Here, we introduce CapTS-seq, a new strategy for sequencing small RNAs that combines chemical capping and template switching. We probed a variety of non-native synthetic cap structures and found that an unmethylated guanosine triphosphate cap led to the lowest bias and highest efficiency for template switching. Through cross-examination of different nucleotides at the cap position, our data provided unequivocal evidence that the 5′ cap acts as a template for the first nucleotide in reverse transcriptase-mediated post-templated addition to the emerging cDNA—a key feature to propel template switching. We deployed CapTS-seq for sequencing synthetic miRNAs, human total brain and liver FFPE RNA, and demonstrated that it consistently improves library quality for miRNAs in comparison with a gold standard template switching-based small RNA-seq kit.

miRge3.0: a comprehensive microRNA and tRF sequencing analysis pipeline

MicroRNAs and tRFs are classes of small non-coding RNAs, known for their roles in translational regulation of genes. Advances in next-generation sequencing (NGS) have enabled high-throughput small RNA-seq studies, which require robust alignment pipelines. Our laboratory previously developed miRge and miRge2.0, as flexible tools to process sequencing data for annotation of miRNAs and other small-RNA species and further predict novel miRNAs using a support vector machine approach. Although miRge2.0 is a leading analysis tool in terms of speed with unique quantifying and annotation features, it has a few limitations. We present miRge3.0 that provides additional features along with compatibility to newer versions of Cutadapt and Python. The revisions of the tool include the ability to process Unique Molecular Identifiers (UMIs) to account for PCR duplicates while quantifying miRNAs in the datasets, correct erroneous single base substitutions in miRNAs with miREC and an accurate mirGFF3 formatted isomiR tool. miRge3.0 also has speed improvements benchmarked to miRge2.0, Chimira and sRNAbench. Finally, miRge3.0 output integrates into other packages for a streamlined analysis process and provides a cross-platform Graphical User Interface (GUI). In conclusion miRge3.0 is our third generation small RNA-seq aligner with improvements in speed, versatility and functionality over earlier iterations.

An essential role for the piRNA pathway in regulating the ribosomal RNA pool in C. elegans

Piwi-interacting RNAs (piRNAs) are RNA effectors with key roles in maintaining genome integrity and promoting fertility in metazoans. In Caenorhabditis elegans loss of piRNAs leads to a transgenerational sterility phenotype. The plethora of piRNAs and their ability to silence transcripts with imperfect complementarity have raised several (non-exclusive) models for the underlying drivers of sterility. Here, we report the extranuclear and transferable nature of the sterility driver, its suppression via mutations disrupting the endogenous RNAi and poly-uridylation machinery, and copy-number amplification at the ribosomal DNA locus. In piRNA-deficient animals, several small interfering RNA (siRNA) populations become increasingly overabundant in the generations preceding loss of germline function, including ribosomal siRNAs (risiRNAs). A concomitant increase in uridylated sense rRNA fragments suggests that poly-uridylation may potentiate RNAi-mediated gene silencing of rRNAs. We conclude that loss of the piRNA machinery allows for unchecked amplification of siRNA populations, originating from abundant highly structured RNAs, to deleterious levels.

Single-cell microRNA sequencing method comparison and application to cell lines and circulating lung tumor cells

Molecular single cell analyses provide insights into physiological and pathological processes. Here, in a stepwise approach, we first evaluate 19 protocols for single cell small RNA sequencing on MCF7 cells spiked with 1 pg of 1,006 miRNAs. Second, we analyze MCF7 single cell equivalents of the eight best protocols. Third, we sequence single cells from eight different cell lines and 67 circulating tumor cells (CTCs) from seven SCLC patients. Altogether, we analyze 244 different samples. We observe high reproducibility within protocols and reads covered a broad spectrum of RNAs. For the 67 CTCs, we detect a median of 68 miRNAs, with 10 miRNAs being expressed in 90% of tested cells. Enrichment analysis suggested the lung as the most likely organ of origin and enrichment of cancer-related categories. Even the identification of non-annotated candidate miRNAs was feasible, underlining the potential of single cell small RNA sequencing.

Technologies for small non-coding RNA sequencing at the single-cell level are less mature than for sequencing mRNAs. Here the authors evaluate available protocols for analysis of circulating lung cancer tumour cells.

Effects of non-inherited ancestral genotypes on offspring phenotypes

It is well established that environmental exposures can modify the profile of heritable factors in an individual's germ cells, ultimately affecting the inheritance of phenotypes in descendants. Similar to exposures, an ancestor's genotype can also affect the inheritance of phenotypes across generations, sometimes in offspring who do not inherit the genetic aberration. This can occur via a variety of prenatal, in utero, or postnatal mechanisms. In this review, we discuss the evidence for this process in mammals, with a focus on examples that are potentially mediated through the germline, while also considering alternate routes of inheritance. Non-inherited ancestral genotypes may influence descendant's disease risk to a much greater extent than currently appreciated, and focused evaluation of this phenomenon may reveal novel mechanisms of inheritance.

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.

Paired-end small RNA sequencing reveals a possible overestimation in the isomiR sequence repertoire previously reported from conventional single read data analysis

BACKGROUND

Next generation sequencing has allowed the discovery of miRNA isoforms, termed isomiRs. Some isomiRs are derived from imprecise processing of pre-miRNA precursors, leading to length variants. Additional variability is introduced by non-templated addition of bases at the ends or editing of internal bases, resulting in base differences relative to the template DNA sequence. We hypothesized that some component of the isomiR variation reported so far could be due to systematic technical noise and not real.

RESULTS

We have developed the XICRA pipeline to analyze small RNA sequencing data at the isomiR level. We exploited its ability to use single or merged reads to compare isomiR results derived from paired-end (PE) reads with those from single reads (SR) to address whether detectable sequence differences relative to canonical miRNAs found in isomiRs are true biological variations or the result of errors in sequencing. We have detected non-negligible systematic differences between SR and PE data which primarily affect putative internally edited isomiRs, and at a much smaller frequency terminal length changing isomiRs. This is relevant for the identification of true isomiRs in small RNA sequencing datasets.

CONCLUSIONS

We conclude that potential artifacts derived from sequencing errors and/or data processing could result in an overestimation of abundance and diversity of miRNA isoforms. Efforts in annotating the isomiRnome should take this into account.

Commercially Available Blocking Oligonucleotides Effectively Suppress Unwanted Hemolysis-Related miRNAs in a Large Whole-Blood RNA Cohort

When sequencing small RNA libraries derived from whole blood, the most abundant microRNAs (miRs) detected are often miR-486-5p, miR-451a, and miR-92a-3p. These highly expressed erythropoietic miRs are released into the sample from red blood cell hemolysis. Next-generation sequencing of these unwanted miRs leads to a waste in sequencing cost and diminished detection of lowly expressed miRNAs, including many potential miRNA biomarkers. Previous work has developed a method to reduce targeted miRNAs using oligonucleotides that bind their target miRNA and prevent its ligation during library construction, although the extent to which oligonucleotides can be multiplexed and their effect on larger cohorts has not been thoroughly explored. We present a method for suppressing detection of three highly abundant heme miRs in a single multiplexed blocking oligonucleotide reaction. In a small paired-sample pilot (n = 8) and a large cohort of samples (n = 901), multiplexed oligos reduced detection of their target miRNAs by approximately 70%, allowing for an approximately 10-fold increase in reads mapping to nonheme miRs and increased detection of very lowly expressed miRs, with minimal off-target effects. By removing all three highly expressed erythropoietic miRNAs from next-generational sequencing libraries, this commercially available multiplexed blocking oligonucleotide method allows for greater detection of lowly expressed biomarkers, improving the efficacy, cost-efficiency, and sensitivity of biomarker studies and diagnostic tests.

EZcount: An all-in-one software for microRNA expression quantification from NGS sequencing data

MicroRNAs (miRNAs) are short endogenous molecules of RNA that influence cell regulation by suppressing genes. Their ubiquity throughout all branches of the tree of life has suggested their central role in many cellular functions. Nowadays, several personalized medicine applications rely on miRNAs as biomarkers for diagnoses, prognoses, and prediction of drug response. The increasing ease of sequencing miRNAs contrasts with the difficulty of accurately quantifying their concentration. The use of general purpose aligners is only a partial solution as they have limited possibilities to accurately solve ambiguous mapping due to the short length of these sequences. We developed EZcount, an all-in-one software that, with a single command, performs the entire quantification process: from raw fastq files to read counts. Experiments show that EZcount is more sensitive and accurate than methods based on sequence alignment, independently of the library preparation protocol and sequencing machine. The parallel architecture of EZcount makes it fast enough to process a sample in minutes using a standard workstation. EZcount runs on all of the most common operating systems (Linux, Windows and MacOS) and is freely available for download at https://gitlab.com/BioAlgo/miR-pipe. A detailed description of the datasets, the raw experimental results, and all the scripts used for testing are available as supplementary material.

An integrated workflow for biomarker development using microRNAs in extracellular vesicles for cancer precision medicine

EV-miRNAs are microRNA (miRNA) molecules encapsulated in extracellular vesicles (EVs), which play crucial roles in tumor pathogenesis, progression, and metastasis. Recent studies about EV-miRNAs have gained novel insights into cancer biology and have demonstrated a great potential to develop novel liquid biopsy assays for various applications. Notably, compared to conventional liquid biomarkers, EV-miRNAs are more advantageous in representing host-cell molecular architecture and exhibiting higher stability and specificity. Despite various available techniques for EV-miRNA separation, concentration, profiling, and data analysis, a standardized approach for EV-miRNA biomarker development is yet lacking. In this review, we performed a substantial literature review and distilled an integrated workflow encompassing important steps for EV-miRNA biomarker development, including sample collection and EV isolation, EV-miRNA extraction and quantification, high-throughput data preprocessing, biomarker prioritization and model construction, functional analysis, as well as validation. With the rapid growth of "big data", we highlight the importance of efficient mining of high-throughput data for the discovery of EV-miRNA biomarkers and integrating multiple independent datasets for in silico and experimental validations to increase the robustness and reproducibility. Furthermore, as an efficient strategy in systems biology, network inference provides insights into the regulatory mechanisms and can be used to select functionally important EV-miRNAs to refine the biomarker candidates. Despite the encouraging development in the field, a number of challenges still hinder the clinical translation. We finally summarize several common challenges in various biomarker studies and discuss potential opportunities emerging in the related fields.

On the impact of batch effect correction in TCGA isomiR expression data

MicroRNAs (miRNAs) are small non-coding RNAs with diverse functions in post-transcriptional regulation of gene expression. Sequence and length variants of miRNAs are called isomiRs and can exert different functions compared to their canonical counterparts. The Cancer Genome Atlas (TCGA) provides isomiR-level expression data for patients of various cancer entities collected in a multi-center approach over several years. However, the impact of batch effects within individual cohorts has not been systematically investigated and corrected for before. Therefore, the aim of this study was to identify relevant cohort-specific batch variables and generate batch-corrected isomiR expression data for 16 TCGA cohorts. The main batch variables included sequencing platform, plate, sample purity and sequencing depth. Platform bias was related to certain length and sequence features of individual recurrently affected isomiRs. Furthermore, significant downregulation of reported tumor suppressive isomiRs in lung tumor tissue compared to normal samples was only observed after batch correction, highlighting the importance of working with corrected data. Batch-corrected datasets for all cohorts including quality control are provided as supplement. In summary, this study reveals that batch effects present in the TCGA dataset might mask biologically relevant effects and provides a valuable resource for research on isomiRs in cancer (accessible through GEO: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164767).

An improved method for specific-target preamplification PCR analysis of single blastocysts useful for embryo sexing and high-throughput gene expression analysis

Gene expression analysis in preimplantation embryos has been used for answering fundamental questions related to development, prediction of pregnancy outcome, and other topics. Limited amounts of mRNA in preimplantation embryos hinders progress in studying the preimplantation embryo. Here, a method was developed involving direct synthesis and specific-target preamplification (STA) of cDNA for gene expression analysis in single blastocysts. Effective cell lysis and genomic DNA removal steps were incorporated into the method. In addition, conditions for real-time PCR of cDNA generated from these processes were improved. By using this system, reliable embryo sexing results based on expression of sex-chromosome linked genes was demonstrated. Calibration curve analysis of PCR results using the Fluidigm Biomark microfluidic platform (Fluidigm, South San Francisco, CA) was performed to evaluate 96 STA cDNA from single blastocysts. In total, 93.75% of the genes were validated. Robust amplification was detected even when STA cDNA from a single blastocyst was diluted 1,024-fold. Further analysis showed that within-assay variation increased when cycle threshold values exceeded 18. Overall, STA quantitative real-time PCR analysis was shown to be useful for analysis of gene expression of multiple specific targets in single blastocysts.

Quantification of purified endogenous miRNAs with high sensitivity and specificity

MicroRNAs (miRNAs) are short (19-24 nt) non-coding RNAs that suppress the expression of protein coding genes at the post-transcriptional level. Differential expression profiles of miRNAs across a range of diseases have emerged as powerful biomarkers, making a reliable yet rapid profiling technique for miRNAs potentially essential in clinics. Here, we report an amplification-free multi-color single-molecule imaging technique that can profile purified endogenous miRNAs with high sensitivity, specificity, and reliability. Compared to previously reported techniques, our technique can discriminate single base mismatches and single-nucleotide 3'-tailing with low false positive rates regardless of their positions on miRNA. By preloading probes in Thermus thermophilus Argonaute (TtAgo), miRNAs detection speed is accelerated by more than 20 times. Finally, by utilizing the well-conserved linearity between single-molecule spot numbers and the target miRNA concentrations, the absolute average copy numbers of endogenous miRNA species in a single cell can be estimated. Thus our technique, Ago-FISH (Argonaute-based Fluorescence In Situ Hybridization), provides a reliable way to accurately profile various endogenous miRNAs on a single miRNA sensing chip.

Small RNA Sequencing to Discover Circulating MicroRNA Biomarkers of Testicular Toxicity in Dogs

Predictive indicators of testicular toxicity could improve drug development by allowing early in-life screening for this adverse effect before it becomes severe. We hypothesized that circulating microRNAs (miRNAs) could serve as testicular toxicity biomarkers in dogs. Herein, we describe the results of an exploratory study conducted to discover biomarkers of drug-induced testicular injury. Following a dose-selection study using the testicular toxicant ethylene glycol monomethyl ether (EGME), we chose a dose of 50 mg/kg/d EGME to avoid systemic toxicity and treated 2 groups of dogs (castrated, non-castrated) for 14 to 28 days. Castrated animals were used as negative controls to identify biomarkers specific for testicular toxicity because EGME can cause toxicity to organ systems in addition to the testis. Blood was collected daily during the dosing period, followed by recovery for 29 to 43 days with less frequent sampling. Dosing was well tolerated, resulting in mild-to-moderate degeneration in testes and epididymides. Global profiling of serum miRNAs at selected dosing and recovery time points was completed by small RNA sequencing. Bioinformatics data analysis using linear modeling demonstrated several circulating miRNAs that were differentially abundant during the dosing period compared with baseline and/or castrated control samples. Confirmatory reverse transcription quantitative polymerase chain reaction data in these animals was unable to detect sustained alterations of miRNAs in serum, except for 1 potential candidate cfa-miR-146b. Taken together, we report the results of a comprehensive exploratory study and suggest future directions for follow-up research to address the challenge of developing diagnostic biomarkers of testicular toxicity.

OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data

Over the past decade, a large amount of RNA sequencing (RNA-seq) data were deposited in public repositories, and more are being produced at an unprecedented rate. However, there are few open source tools with point-and-click interfaces that are versatile and offer streamlined comprehensive analysis of RNA-seq datasets. To maximize the capitalization of these vast public resources and facilitate the analysis of RNA-seq data by biologists, we developed a web application called OneStopRNAseq for the one-stop analysis of RNA-seq data. OneStopRNAseq has user-friendly interfaces and offers workflows for common types of RNA-seq data analyses, such as comprehensive data-quality control, differential analysis of gene expression, exon usage, alternative splicing, transposable element expression, allele-specific gene expression quantification, and gene set enrichment analysis. Users only need to select the desired analyses and genome build, and provide a Gene Expression Omnibus (GEO) accession number or Dropbox links to sequence files, alignment files, gene-expression-count tables, or rank files with the corresponding metadata. Our pipeline facilitates the comprehensive and efficient analysis of private and public RNA-seq data.

A low-bias and sensitive small RNA library preparation method using randomized splint ligation

Small RNAs are important regulators of gene expression and are involved in human development and disease. Next generation sequencing (NGS) allows for scalable, genome-wide studies of small RNA; however, current methods are challenged by low sensitivity and high bias, limiting their ability to capture an accurate representation of the cellular small RNA population. Several studies have shown that this bias primarily arises during the ligation of single-strand adapters during library preparation, and that this ligation bias is magnified by 2′-O-methyl modifications (2′OMe) on the 3′ terminal nucleotide. In this study, we developed a novel library preparation process using randomized splint ligation with a cleavable adapter, a design which resolves previous challenges associated with this ligation strategy. We show that a randomized splint ligation based workflow can reduce bias and increase the sensitivity of small RNA sequencing for a wide variety of small RNAs, including microRNA (miRNA) and tRNA fragments as well as 2′OMe modified RNA, including Piwi-interacting RNA and plant miRNA. Finally, we demonstrate that this workflow detects more differentially expressed miRNA between tumorous and matched normal tissues. Overall, this library preparation process allows for highly accurate small RNA sequencing and will enable studies of 2′OMe modified RNA with new levels of detail.

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API

MOTIVATION

MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods.

RESULTS

To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification.

AVAILABILITY AND IMPLEMENTATION

https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dynamic mRNP Remodeling in Response to Internal and External Stimuli

Signal transduction and the regulation of gene expression are fundamental processes in every cell. RNA-binding proteins (RBPs) play a key role in the post-transcriptional modulation of gene expression in response to both internal and external stimuli. However, how signaling pathways regulate the assembly of RBPs with mRNAs remains largely unknown. Here, we summarize observations showing that the formation and composition of messenger ribonucleoprotein particles (mRNPs) is dynamically remodeled in space and time by specific signaling cascades and the resulting post-translational modifications. The integration of signaling events with gene expression is key to the rapid adaptation of cells to environmental changes and stress. Only a combined approach analyzing the signal transduction pathways and the changes in post-transcriptional gene expression they cause will unravel the mechanisms coordinating these important cellular processes.

Insufficiently complex unique-molecular identifiers (UMIs) distort small RNA sequencing

The attachment of unique molecular identifiers (UMIs) to RNA molecules prior to PCR amplification and sequencing, makes it possible to amplify libraries to a level that is sufficient to identify rare molecules, whilst simultaneously eliminating PCR bias through the identification of duplicated reads. Accurate de-duplication is dependent upon a sufficiently complex pool of UMIs to allow unique labelling. In applications dealing with complex libraries, such as total RNA-seq, only a limited variety of UMIs are required as the variation in molecules to be sequenced is enormous. However, when sequencing a less complex library, such as small RNAs for which there is a more limited range of possible sequences, we find increased variation in UMIs are required, even beyond that provided in a commercial kit specifically designed for the preparation of small RNA libraries for sequencing. We show that a pool of UMIs randomly varying across eight nucleotides is not of sufficient depth to uniquely tag the microRNAs to be sequenced. This results in over de-duplication of reads and the marked under-estimation of expression of the more abundant microRNAs. Whilst still arguing for the utility of UMIs, this work demonstrates the importance of their considered design to avoid errors in the estimation of gene expression in libraries derived from select regions of the transcriptome or small genomes.

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.