Systematic comparison of small RNA library preparation protocols for next-generation sequencing

De novo mutations promote inflammation in children with STAT3 gain-of-function syndrome by affecting IL-1β expression

STAT3 gain-of-function syndrome, characterized by early-onset autoimmunity and primary immune regulatory disorder, remains poorly understood in terms of its immunological mechanisms. We employed whole-genome sequencing of familial trios to elucidate the pivotal role of de novo mutations in genetic diseases. We identified 37 high-risk pathogenic loci affecting 23 genes, including a novel STAT3 c.508G>A mutation. We also observed significant down-regulation of pathogenic genes in affected individuals, potentially associated with inflammatory responses regulated by PTPN14 via miR378c. These findings enhance our understanding of the pathogenesis of STAT3 gain-of-function syndrome and suggest potential therapeutic strategies. Notably, combined JAK inhibitors and IL-6R antagonists may offer promising treatment avenues for mitigating the severity of STAT3 gain-of-function syndrome.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the diversity in RNA size and chemical modification makes it difficult to capture all RNAs in a cell. We developed a method that combines quasi-random priming with template switching to construct sequencing libraries from RNA molecules of any length and with any type of 3' modifications, allowing for the sequencing of virtually all RNA species. Our ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to identify and quantify all classes of coding and non-coding RNAs. With LIDAR, we comprehensively characterized the transcriptomes of mouse embryonic stem cells, neural progenitor cells, mouse tissues, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared with traditional ligation-dependent sequencing methods and uncovered tDRs with blocked 3' ends that had previously escaped detection. Therefore, LIDAR can capture all RNAs in a sample and uncover RNA species with potential regulatory functions.

Guidelines for mitochondrial RNA analysis

Mitochondria are the energy-producing organelles of mammalian cells with critical involvement in metabolism and signaling. Studying their regulation in pathological conditions may lead to the discovery of novel drugs to treat, for instance, cardiovascular or neurological diseases, which affect high-energy-consuming cells such as cardiomyocytes, hepatocytes, or neurons. Mitochondria possess both protein-coding and noncoding RNAs, such as microRNAs, long noncoding RNAs, circular RNAs, and piwi-interacting RNAs, encoded by the mitochondria or the nuclear genome. Mitochondrial RNAs are involved in anterograde-retrograde communication between the nucleus and mitochondria and play an important role in physiological and pathological conditions. Despite accumulating evidence on the presence and biogenesis of mitochondrial RNAs, their study continues to pose significant challenges. Currently, there are no standardized protocols and guidelines to conduct deep functional characterization and expression profiling of mitochondrial RNAs. To overcome major obstacles in this emerging field, the EU-CardioRNA and AtheroNET COST Action networks summarize currently available techniques and emphasize critical points that may constitute sources of variability and explain discrepancies between published results. Standardized methods and adherence to guidelines to quantify and study mitochondrial RNAs in normal and disease states will improve research outputs, their reproducibility, and translation potential to clinical application.

SingmiR: a single-cell miRNA alignment and analysis tool

Single-cell RNA sequencing (RNA-seq) has revolutionized our understanding of cell biology, developmental and pathophysiological molecular processes, paving the way toward novel diagnostic and therapeutic approaches. However, most of the gene regulatory processes on the single-cell level are still unknown, including post-transcriptional control conferred by microRNAs (miRNAs). Like the established single-cell gene expression analysis, advanced computational expertise is required to comprehensively process newly emerging single-cell miRNA-seq datasets. A web server providing a workflow tailored for single-cell miRNA-seq data with a self-explanatory interface is currently not available. Here, we present SingmiR, enabling the rapid (pre-)processing and quantification of human miRNAs from noncoding single-cell samples. It performs read trimming for different library preparation protocols, generates automated quality control reports and provides feature-normalized count files. Numerous standard and advanced analyses such as dimension reduction, clustered feature heatmaps, sample correlation heatmaps and differential expression statistics are implemented. We aim to speed up the prototyping pipeline for biologists developing single-cell miRNA-seq protocols on small to medium-sized datasets. SingmiR is freely available to all users without the need for a login at https://www.ccb.uni-saarland.de/singmir.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Characterisation and engineering of a thermophilic RNA ligase from Palaeococcus pacificus

RNA ligases are important enzymes in molecular biology and are highly useful for the manipulation and analysis of nucleic acids, including adapter ligation in next-generation sequencing of microRNAs. Thermophilic RNA ligases belonging to the RNA ligase 3 family are gaining attention for their use in molecular biology, for example a thermophilic RNA ligase from Methanobacterium thermoautotrophicum is commercially available for the adenylation of nucleic acids. Here we extensively characterise a newly identified RNA ligase from the thermophilic archaeon Palaeococcus pacificus (PpaRnl). PpaRnl exhibited significant substrate adenylation activity but low ligation activity across a range of oligonucleotide substrates. Mutation of Lys92 in motif I to alanine, resulted in an enzyme that lacked adenylation activity, but demonstrated improved ligation activity with pre-adenylated substrates (ATP-independent ligation). Subsequent structural characterisation revealed that in this mutant enzyme Lys238 was found in two alternate positions for coordination of the phosphate tail of ATP. In contrast mutation of Lys238 in motif V to glycine via structure-guided engineering enhanced ATP-dependent ligation activity via an arginine residue compensating for the absence of Lys238. Ligation activity for both mutations was higher than the wild-type, with activity observed across a range of oligonucleotide substrates with varying sequence and secondary structure.

Single-Cell Transcriptome Sequencing Reveals Molecular Expression Differences and Marker Genes in Testes during the Sexual Maturation of Mongolian Horses

This study aimed to investigate differences in testicular tissue morphology, gene expression, and marker genes between sexually immature (1-year-old) and sexually mature (10-year-old) Mongolian horses. The purposes of our research were to provide insights into the reproductive physiology of male Mongolian horses and to identify potential markers for sexual maturity. The methods we applied included the transcriptomic profiling of testicular cells using single-cell sequencing techniques. Our results revealed significant differences in tissue morphology and gene expression patterns between the two age groups. Specifically, 25 cell clusters and 10 cell types were identified, including spermatogonial and somatic cells. Differential gene expression analysis highlighted distinct patterns related to cellular infrastructure in sexually immature horses and spermatogenesis in sexually mature horses. Marker genes specific to each stage were also identified, including APOA1, AMH, TAC3, INHA, SPARC, and SOX9 for the sexually immature stage, and PRM1, PRM2, LOC100051500, PRSS37, HMGB4, and H1-9 for the sexually mature stage. These findings contribute to a deeper understanding of testicular development and spermatogenesis in Mongolian horses and have potential applications in equine reproductive biology and breeding programs. In conclusion, this study provides valuable insights into the molecular mechanisms underlying sexual maturity in Mongolian horses.

Small RNA structural biochemistry in a post-sequencing era

High-throughput sequencing has had an enormous impact on small RNA research during the past decade. However, sequencing only offers a one-dimensional view of the transcriptome and is often highly biased. Additionally, the 'sequence, map and annotate' approach, used widely in small RNA research, can lead to flawed interpretations of the data, lacking biological plausibility, due in part to database issues. Even in the absence of technical biases, the loss of three-dimensional information is a major limitation to understanding RNA stability, turnover and function. For example, noncoding RNA-derived fragments seem to exist mainly as dimers, tetramers or as nicked forms of their parental RNAs, contrary to widespread assumptions. In this perspective, we will discuss main sources of bias during small RNA-sequencing, present several useful bias-reducing strategies and provide guidance on the interpretation of small RNA-sequencing results, with emphasis on RNA fragmentomics. As sequencing offers a one-dimensional projection of a four-dimensional reality, prior structure-level knowledge is often needed to make sense of the data. Consequently, while less-biased sequencing methods are welcomed, integration of orthologous experimental techniques is also strongly recommended.

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.

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.

Comparison of Extracellular Vesicle Isolation Methods for miRNA Sequencing

MicroRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) are potential diagnostic and prognostic biomarkers. However, discrepancies in miRNA patterns and their validation are still frequent due to differences in sample origin, EV isolation, and miRNA sequencing methods. The aim of the present study is to find a reliable EV isolation method for miRNA sequencing, adequate for clinical application. To this aim, two comparative studies were performed in parallel with the same human plasma sample: (i) isolation and characterization of EVs obtained using three procedures: size exclusion chromatography (SEC), iodixanol gradient (GRAD), and its combination (SEC+GRAD) and (ii) evaluation of the yield of miRNA sequences obtained using NextSeq 500 (Illumina) and three miRNA library preparation protocols: NEBNext, NEXTFlex, and SMARTer smRNA-seq. The conclusion of comparison (i) is that recovery of the largest amount of EVs and reproducibility were attained with SEC, but GRAD and SEC+GRAD yielded purer EV preparations. The conclusion of (ii) is that the NEBNext library showed the highest reproducibility in the number of miRNAs recovered and the highest diversity of miRNAs. These results render the combination of GRAD EV isolation and NEBNext library preparation for miRNA retrieval as adequate for clinical applications using plasma samples.

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.

Non-coding RNAs from seminal plasma extracellular vesicles and success of live birth among couples undergoing fertility treatment

Background: Infertility remains a global health problem with male-factor infertility accounting for around 50% of cases. Understanding the molecular markers for the male contribution of live birth success has been limited. Here, we evaluated the expression levels of seminal plasma extracellular vesicle (spEV) non-coding RNAs (ncRNAs) in men of couples in relation with those with and without a successful live birth after infertility treatment. Method: Sperm-free spEV small RNA profiles were generated from 91 semen samples collected from male participants of couples undergoing assisted reproductive technology (ART) treatment. Couples were classified into two groups based on successful live birth (yes, n = 28) and (no, n = 63). Mapping of reads to human transcriptomes followed the order: miRNA > tRNA > piRNA > rRNA> "other" RNA > circRNA > lncRNA. Differential expression analysis of biotype-specific normalized read counts between groups were assessed using EdgeR (FDR<0.05). Result: We found a total of 12 differentially expressed spEV ncRNAs which included 10 circRNAs and two piRNAs between the live birth groups. Most (n = 8) of the identified circRNAs were downregulated in the no live birth group and targeted genes related to ontology terms such as negative reproductive system and head development, tissue morphogenesis, embryo development ending in birth or egg hatching, and vesicle-mediated transport. The differentially upregulated piRNAs overlapped with genomic regions including coding PID1 genes previously known to play a role in mitochondrion morphogenesis, signal transduction and cellular proliferation. Conclusion: This study identified novel ncRNAs profiles of spEVs differentiating men of couples with and without live birth and emphasizes the role of the male partner for ART success.

A ligation-independent sequencing method reveals tRNA-derived RNAs with blocked 3' termini

Despite the numerous sequencing methods available, the vast diversity in size and chemical modifications of RNA molecules makes the capture of the full spectrum of cellular RNAs a difficult task. By combining quasi-random hexamer priming with a custom template switching strategy, we developed a method to construct sequencing libraries from RNA molecules of any length and with any type of 3' terminal modification, allowing the sequencing and analysis of virtually all RNA species. Ligation-independent detection of all types of RNA (LIDAR) is a simple, effective tool to comprehensively characterize changes in small non-coding RNAs and mRNAs simultaneously, with performance comparable to separate dedicated methods. With LIDAR, we comprehensively characterized the coding and non-coding transcriptome of mouse embryonic stem cells, neural progenitor cells, and sperm. LIDAR detected a much larger variety of tRNA-derived RNAs (tDRs) compared to traditional ligation-dependent sequencing methods, and uncovered the presence of tDRs with blocked 3' ends that had previously escaped detection. Our findings highlight the potential of LIDAR to systematically detect all RNAs in a sample and uncover new RNA species with potential regulatory functions.

Adapter dimer contamination in sRNA-sequencing datasets predicts sequencing failure and batch effects and hampers extracellular vesicle-sRNA analysis

Small RNA (sRNA) profiling of Extracellular Vesicles (EVs) by Next-Generation Sequencing (NGS) often delivers poor outcomes, independently of reagents, platforms or pipelines used, which contributes to poor reproducibility of studies. Here we analysed pre/post-sequencing quality controls (QC) to predict issues potentially biasing biological sRNA-sequencing results from purified human milk EVs, human and mouse EV-enriched plasma and human paraffin-embedded tissues. Although different RNA isolation protocols and NGS platforms were used in these experiments, all datasets had samples characterized by a marked removal of reads after pre-processing. The extent of read loss between individual samples within a dataset did not correlate with isolated RNA quantity or sequenced base quality. Rather, cDNA electropherograms revealed the presence of a constant peak whose intensity correlated with the degree of read loss and, remarkably, with the percentage of adapter dimers, which were found to be overrepresented sequences in high read-loss samples. The analysis through a QC pipeline, which allowed us to monitor quality parameters in a step-by-step manner, provided compelling evidence that adapter dimer contamination was the main factor causing batch effects. We concluded this study by summarising peer-reviewed published workflows that perform consistently well in avoiding adapter dimer contamination towards a greater likelihood of sequencing success.

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.

Diversity and Differential Expression of MicroRNAs in the Human Skeletal Muscle with Distinct Fiber Type Composition

The ratio of fast- and slow-twitch fibers in human skeletal muscle is variable and largely determined by genetic factors. In this study, we investigated the contribution of microRNA (miRNA) in skeletal muscle fiber type composition. The study involved biopsy samples of the vastus lateralis muscle from 24 male participants with distinct fiber type ratios. The miRNA study included samples from five endurance athletes and five power athletes with the predominance of slow-twitch (61.6-72.8%) and fast-twitch (69.3-80.7%) fibers, respectively. Total and small RNA were extracted from tissue samples. Total RNA sequencing (N = 24) revealed 352 differentially expressed genes between the groups with the predominance of fast- and slow-twitch muscle fibers. Small RNA sequencing showed upregulation of miR-206, miR-501-3p and miR-185-5p, and downregulation of miR-499a-5p and miR-208-5p in the group of power athletes with fast-twitch fiber predominance. Two miRtronic miRNAs, miR-208b-3p and miR-499a-5p, had strong correlations in expression with their host genes (MYH7 and MYH7B, respectively). Correlations between the expression of miRNAs and their experimentally validated messenger RNA (mRNA) targets were calculated, and 11 miRNA-mRNA interactions with strong negative correlations were identified. Two of them belonged to miR-208b-3p and miR-499a-5p, indicating their regulatory links with the expression of CDKN1A and FOXO4, respectively.

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.

The Impact of the Anticoagulant Type in Blood Collection Tubes on Circulating Extracellular Plasma MicroRNA Profiles Revealed by Small RNA Sequencing

Pre-analytical factors have a significant influence on circulating microRNA (miRNA) profiling. The aim of this study was a comprehensive assessment of the impact of the anticoagulant type in blood collection tubes on circulating plasma miRNA profiles using small RNA sequencing. Blood from ten healthy participants (five males and five females from 25 to 40 years old) was taken in collection tubes with four different anticoagulants: acid citrate dextrose (ACD-B), sodium citrate, citrate-theophylline-adenosine-dipyridamole (CTAD) and dipotassium-ethylenediaminetetraacetic acid (K2 EDTA). Platelet-free plasma samples were obtained by double centrifugation. EDTA plasma samples had elevated levels of hemolysis compared to samples obtained using other anticoagulants. Small RNA was extracted from plasma samples and small RNA sequencing was performed on the Illumina NextSeq 500 system. A total of 30 samples had been successfully sequenced starting from ~1 M reads mapped to miRNAs, allowing us to analyze their diversity and isoform content. The principal component analysis showed that the EDTA samples have distinct circulating plasma miRNA profiles compared to samples obtained using other anticoagulants. We selected 50 miRNA species that were differentially expressed between the sample groups based on the type of anticoagulant. We found that the EDTA samples had elevated levels of miRNAs which are abundant in red blood cells (RBC) and associated with hemolysis, while the levels of some platelet-specific miRNAs in these samples were lowered. The ratio between RBC-derived and platelet-derived miRNAs differed between the EDTA samples and other sample groups, which was validated by quantitative PCR. This study provides full plasma miRNA profiles of 10 healthy adults, compares them with previous studies and shows that the profile of circulating miRNAs in the EDTA plasma samples is altered primarily due to an increased level of hemolysis.

Small but strong: Pivotal roles and potential applications of snoRNAs in hematopoietic malignancies

Small nucleolar RNAs (snoRNAs) belong to a family of noncoding RNAs that are 60-300 nucleotides in length, and they are classified into two classes according to their structure and function: C/D box snoRNAs, playing an essential role in 2’-O-methylation modification on ribosomal RNA; H/ACA box snoRNAs, involved in the pseudouridylation of rRNA. SnoRNAs with unclear functions, no predictable targets, and unusual subcellular locations are called orphan snoRNAs. Recent studies have revealed abnormal expression and demonstrated the pivotal roles of snoRNAs and their host genes in various types of hematological malignancies. This review discusses recent discoveries concerning snoRNAs in a variety of hematological malignancies, including multiple myeloma, lymphoma and leukemia, and sheds light on the application of snoRNAs as diagnostic and prognostic markers as well as therapeutic targets of hematological malignancies in the future.

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

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

The PPR domain of mitochondrial RNA polymerase is an exoribonuclease required for mtDNA replication in Drosophila melanogaster

Mitochondrial DNA (mtDNA) replication and transcription are of paramount importance to cellular energy metabolism. Mitochondrial RNA polymerase is thought to be the primase for mtDNA replication. However, it is unclear how this enzyme, which normally transcribes long polycistronic RNAs, can produce short RNA oligonucleotides to initiate mtDNA replication. We show that the PPR domain of Drosophila mitochondrial RNA polymerase (PolrMT) has 3'-to-5' exoribonuclease activity, which is indispensable for PolrMT to synthesize short RNA oligonucleotides and prime DNA replication in vitro. An exoribonuclease-deficient mutant, PolrMT^E423P, partially restores mitochondrial transcription but fails to support mtDNA replication when expressed in PolrMT-mutant flies, indicating that the exoribonuclease activity is necessary for mtDNA replication. In addition, overexpression of PolrMT^E423P in adult flies leads to severe neuromuscular defects and a marked increase in mtDNA transcript errors, suggesting that exoribonuclease activity may contribute to the proofreading of mtDNA transcription.

Making Invisible RNA Visible: Discriminative Sequencing Methods for RNA Molecules with Specific Terminal Formations

Next generation sequencing of RNA molecules (RNA-seq) has become a common tool to characterize the expression profiles of RNAs and their regulations in normal physiological processes and diseases. Although increasingly accumulating RNA-seq data are widely available through publicly accessible sites, most of the data for short non-coding RNAs (sncRNAs) have been obtained for microRNA (miRNA) analyses by standard RNA-seq, which only capture the sncRNAs with 5'-phosphate (5'-P) and 3'-hydroxyl (3'-OH) ends. The sncRNAs with other terminal formations such as those with a 5'-hydroxyl end (5'-OH), a 3'-phosphate (3'-P) end, or a 2',3'-cyclic phosphate end (2',3'-cP) cannot be efficiently amplified and sequenced by standard RNA-seq. Due to the invisibility in standard RNA-seq data, these non-miRNA-sncRNAs have been a hidden component in the transcriptome. However, as the functional significances of these sncRNAs have become increasingly apparent, specific RNA-seq methods compatible with various terminal formations of sncRNAs have been developed and started shedding light on the previously unrecognized sncRNAs that lack 5'-P/3'-OH ends. In this review, we summarize the expanding world of sncRNAs with various terminal formations and the strategic approaches of specific RNA-seq methods to distinctively characterize their expression profiles.

Deciphering the RNA universe in sperm in its role as a vertical information carrier

The inheritance of neurophysiologic and neuropsychologic complex diseases can only partly be explained by the Mendelian concept of genetic inheritance. Previous research showed that both psychological disorders like post-traumatic stress disorder and metabolic diseases are more prevalent in the progeny of affected parents. This could suggest an epigenetic mode of transmission. Human studies give first insight into the scope of intergenerational influence of stressors but are limited in exploring the underlying mechanisms. Animal models have elucidated the mechanistic underpinnings of epigenetic transmission. In this review, we summarize progress on the mechanisms of paternal intergenerational transmission by means of sperm RNA in mouse models. We discuss relevant details for the modelling of RNA-mediated transmission, point towards currently unanswered questions and propose experimental considerations for tackling these questions.

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.

Identification of high-confidence human poly(A) RNA isoform scaffolds using nanopore sequencing

Nanopore sequencing devices read individual RNA strands directly. This facilitates identification of exon linkages and nucleotide modifications; however, using conventional direct RNA nanopore sequencing, the 5' and 3' ends of poly(A) RNA cannot be identified unambiguously. This is due in part to RNA degradation in vivo and in vitro that can obscure transcription start and end sites. In this study, we aimed to identify individual full-length human RNA isoforms among ∼4 million nanopore poly(A)-selected RNA reads. First, to identify RNA strands bearing 5' m⁷G caps, we exchanged the biological cap for a modified cap attached to a 45-nt oligomer. This oligomer adaptation method improved 5' end sequencing and ensured correct identification of the 5' m⁷G capped ends. Second, among these 5'-capped nanopore reads, we screened for features consistent with a 3' polyadenylation site. Combining these two steps, we identified 294,107 individual high-confidence full-length RNA scaffolds from human GM12878 cells, most of which (257,721) aligned to protein-coding genes. Of these, 4876 scaffolds indicated unannotated isoforms that were often internal to longer, previously identified RNA isoforms. Orthogonal data for m⁷G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

Recent Major Transcriptomics and Epitranscriptomics Contributions toward Personalized and Precision Medicine

With the advent of genome-wide screening methods-beginning with microarray technologies and moving onto next generation sequencing methods-the era of precision and personalized medicine was born. Genomics led the way, and its contributions are well recognized. However, "other-omics" fields have rapidly emerged and are becoming as important toward defining disease causes and exploring therapeutic benefits. In this review, we focus on the impacts of transcriptomics, and its extension-epitranscriptomics-on personalized and precision medicine efforts. There has been an explosion of transcriptomic studies particularly in the last decade, along with a growing number of recent epitranscriptomic studies in several disease areas. Here, we summarize and overview major efforts for cancer, cardiovascular disease, and neurodevelopmental disorders (including autism spectrum disorder and intellectual disability) for transcriptomics/epitranscriptomics in precision and personalized medicine. We show that leading advances are being made in both diagnostics, and in investigative and landscaping disease pathophysiological studies. As transcriptomics/epitranscriptomics screens become more widespread, it is certain that they will yield vital and transformative precision and personalized medicine contributions in ways that will significantly further genomics gains.

MGcount: a total RNA-seq quantification tool to address multi-mapping and multi-overlapping alignments ambiguity in non-coding transcripts

BACKGROUND

Total-RNA sequencing (total-RNA-seq) allows the simultaneous study of both the coding and the non-coding transcriptome. Yet, computational pipelines have traditionally focused on particular biotypes, making assumptions that are not fullfilled by total-RNA-seq datasets. Transcripts from distinct RNA biotypes vary in length, biogenesis, and function, can overlap in a genomic region, and may be present in the genome with a high copy number. Consequently, reads from total-RNA-seq libraries may cause ambiguous genomic alignments, demanding for flexible quantification approaches.

RESULTS

Here we present Multi-Graph count (MGcount), a total-RNA-seq quantification tool combining two strategies for handling ambiguous alignments. First, MGcount assigns reads hierarchically to small-RNA and long-RNA features to account for length disparity when transcripts overlap in the same genomic position. Next, MGcount aggregates RNA products with similar sequences where reads systematically multi-map using a graph-based approach. MGcount outputs a transcriptomic count matrix compatible with RNA-sequencing downstream analysis pipelines, with both bulk and single-cell resolution, and the graphs that model repeated transcript structures for different biotypes. The software can be used as a python module or as a single-file executable program.

CONCLUSIONS

MGcount is a flexible total-RNA-seq quantification tool that successfully integrates reads that align to multiple genomic locations or that overlap with multiple gene features. Its approach is suitable for the simultaneous estimation of protein-coding, long non-coding and small non-coding transcript concentration, in both precursor and processed forms. Both source code and compiled software are available at https://github.com/hitaandrea/MGcount .

Rolling circle reverse transcription enables high fidelity nanopore sequencing of small RNA

Small RNAs (sRNAs) are an important group of non-coding RNAs that have great potential as diagnostic and prognostic biomarkers for treatment of a wide variety of diseases. The portability and affordability of nanopore sequencing technology makes it ideal for point of care and low resource settings. Currently sRNAs can't be reliably sequenced on the nanopore platform due to the short size of sRNAs and high error rate of the nanopore sequencer. Here, we developed a highly efficient nanopore-based sequencing strategy for sRNAs (SR-Cat-Seq) in which sRNAs are ligated to an adapter, circularized, and undergo rolling circle reverse transcription to generate concatemeric cDNA. After sequencing, the resulting tandem repeat sequences within the individual cDNA can be aligned to generate highly accurate consensus sequences. We compared our sequencing strategy with other sRNA sequencing methods on a short-read sequencing platform and demonstrated that SR-Cat-Seq can obtain low bias and highly accurate sRNA transcriptomes. Therefore, our method could enable nanopore sequencing for sRNA-based diagnostics and other applications.

PIWI-interacting RNAs and human testicular function

Small noncoding RNAs (sncRNAs) are pieces of RNA with a length below 200 bp and represent a diverse group of RNAs having many different biological functions. The best described subtype is the microRNAs which primarily function in posttranscriptional gene regulation and appear essential for most physiological processes. Of particular interest for the germline is the PIWI-interacting RNAs (piRNAs) which are a class of sncRNA of 21-35 bp in length that are almost exclusively found in germ cells. Recently, it has become clear that piRNAs are essential for testicular function, and in this perspective, we outline the current knowledge of piRNAs in humans. Although piRNAs appear unique to germ cells, they have also been described in various somatic cancers and biofluids. Here, we discuss the potential function of piRNAs in somatic tissues and whether detection in biofluids may be used as a biomarker for testicular function. This article is categorized under: Reproductive System Diseases > Genetics/Genomics/Epigenetics Reproductive System Diseases > Molecular and Cellular Physiology.

Small Noncoding RNAs in Reproduction and Infertility

Infertility has been reported as one of the most common reproductive impairments, affecting nearly one in six couples worldwide. A large proportion of infertility cases are diagnosed as idiopathic, signifying a deficit in information surrounding the pathology of infertility and necessity of medical intervention such as assisted reproductive therapy. Small noncoding RNAs (sncRNAs) are well-established regulators of mammalian reproduction. Advanced technologies have revealed the dynamic expression and diverse functions of sncRNAs during mammalian germ cell development. Mounting evidence indicates sncRNAs in sperm, especially microRNAs (miRNAs) and transfer RNA (tRNA)-derived small RNAs (tsRNAs), are sensitive to environmental changes and mediate the inheritance of paternally acquired metabolic and mental traits. Here, we review the critical roles of sncRNAs in mammalian germ cell development. Furthermore, we highlight the functions of sperm-borne sncRNAs in epigenetic inheritance. We also discuss evidence supporting sncRNAs as promising biomarkers for fertility and embryo quality in addition to the present limitations of using sncRNAs for infertility diagnosis and treatment.

Transposable elements and Piwi‑interacting RNAs in hemato‑oncology with a focus on myelodysplastic syndrome (Review)

Our current understanding of hematopoietic stem cell differentiation and the abnormalities that lead to leukemogenesis originates from the accumulation of knowledge regarding protein‑coding genes. However, the possible impact of transposable element (TE) mobilization and the expression of P‑element‑induced WImpy testis‑interacting RNAs (piRNAs) on leukemogenesis has been beyond the scope of scientific interest to date. The expression profiles of these molecules and their importance for human health have only been characterized recently due to the rapid progress of high‑throughput sequencing technology development. In the present review, current knowledge on the expression profile and function of TEs and piRNAs was summarized, with specific focus on their reported involvement in leukemogenesis and pathogenesis of myelodysplastic syndrome.

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.

A New Specific and Sensitive RT-qPCR Method Based on Splinted 5' Ligation for the Quantitative Detection of RNA Species Shorter than microRNAs

Recently, we discovered a new family of unusually short RNAs mapping to 5.8S ribosomal RNA (rRNA) and which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. To confirm these small RNA-sequencing (RNA-Seq) data, validate the existence of the two overly abundant doRNAs-the minimal core 12-nt doRNA sequence and its + 1-nt variant bearing a 5' Cytosine, C-doRNA-and streamline their analysis, we developed a new specific and sensitive splinted 5' ligation reverse transcription (RT)-quantitative polymerase chain reaction (qPCR) method. This method is based on a splint-assisted ligation of an adapter to the 5' end of doRNAs, followed by RT-qPCR amplification and quantitation. Our optimized protocol, which may discriminate between doRNA, C-doRNA, mutated and precursor sequences, can accurately detect as low as 240 copies and is quantitatively linear over a range of 7 logs. This method provides a unique tool to expand and facilitate studies exploring the molecular and cellular biology of RNA species shorter than microRNAs.

Identification of Abundant and Functional dodecaRNAs (doRNAs) Derived from Ribosomal RNA

Using a modified RNA-sequencing (RNA-seq) approach, we discovered a new family of unusually short RNAs mapping to ribosomal RNA 5.8S, which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. Using a new quantitative detection method that we developed, we confirmed our RNA-seq data and determined that the minimal core doRNA sequence and its 13-nt variant C-doRNA (doRNA with a 5' Cytosine) are the two most abundant doRNAs, which, together, may outnumber microRNAs. The C-doRNA/doRNA ratio is stable within species but differed between species. doRNA and C-doRNA are mainly cytoplasmic and interact with heterogeneous nuclear ribonucleoproteins (hnRNP) A0, A1 and A2B1, but not Argonaute 2. Reporter gene activity assays suggest that C-doRNA may function as a regulator of Annexin II receptor (AXIIR) expression. doRNAs are differentially expressed in prostate cancer cells/tissues and may control cell migration. These findings suggest that unusually short RNAs may be more abundant and important than previously thought.

pre-piRNA trimming and 2'-O-methylation protect piRNAs from 3' tailing and degradation in C. elegans

The Piwi-interacting RNA (piRNA) pathway suppresses transposable elements and promotes fertility in diverse organisms. Maturation of piRNAs involves pre-piRNA trimming followed by 2'-O-methylation at their 3' termini. Here, we report that the 3' termini of Caenorhabditis elegans piRNAs are subject to nontemplated nucleotide addition, and piRNAs with 3' addition exhibit extensive base-pairing interaction with their target RNAs. Animals deficient for PARN-1 (pre-piRNA trimmer) and HENN-1 (2'-O-methyltransferase) accumulate piRNAs with 3' nontemplated nucleotides. In henn-1 mutants, piRNAs are shortened prior to 3' addition, whereas long isoforms of untrimmed piRNAs are preferentially modified in parn-1 mutant animals. Loss of either PARN-1 or HENN-1 results in modest reduction in steady-state levels of piRNAs. Deletion of both enzymes leads to depletion of piRNAs, desilenced piRNA targets, and impaired fecundity. Together, our findings suggest that pre-piRNA trimming and 2'-O-methylation act collaboratively to protect piRNAs from tailing and degradation.

Single paternal dexamethasone challenge programs offspring metabolism and reveals multiple candidates in RNA-mediated inheritance

Single traumatic events that elicit an exaggerated stress response can lead to the development of neuropsychiatric conditions. Rodent studies suggested germline RNA as a mediator of effects of chronic environmental exposures to the progeny. The effects of an acute paternal stress exposure on the germline and their potential consequences on offspring remain to be seen. We find that acute administration of an agonist for the stress-sensitive Glucocorticoid receptor, using the common corticosteroid dexamethasone, affects the RNA payload of mature sperm as soon as 3 hr after exposure. It further impacts early embryonic transcriptional trajectories, as determined by single-embryo sequencing, and metabolism in the offspring. We show persistent regulation of tRNA fragments in sperm and descendant 2-cell embryos, suggesting transmission from sperm to embryo. Lastly, we unravel environmentally induced alterations in sperm circRNAs and their targets in the early embryo, highlighting this class as an additional candidate in RNA-mediated inheritance of disease risk.

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.

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.

Distinct Extracellular RNA Profiles in Different Plasma Components

Circulating extracellular RNAs (exRNAs) have great potential to serve as biomarkers for a wide range of diagnostic, therapeutic, and prognostic applications. So far, knowledge of the difference among different sources of exRNAs is limited. To address this issue, we performed a sequential physical and biochemical precipitation to collect four fractions (platelets and cell debris, the thrombin-induced precipitates, extracellular vesicles, and supernatant) from each of 10 plasma samples. From total RNAs of the 40 fractions, we prepared ligation-free libraries to profile full spectrum of all RNA species, without size selection and rRNA reduction. Due to complicated RNA composition in these libraries, we utilized a successive stepwise alignment strategy to map the RNA sequences to different RNA categories, including miRNAs, piwi-interacting RNAs, tRNAs, rRNAs, lincRNAs, snoRNAs, snRNAs, other ncRNAs, protein coding RNAs, and circRNAs. Our data showed that each plasma fraction had its own unique distribution of RNA species. Hierarchical cluster analyses using transcript abundance demonstrated similarities in the same plasma fraction and significant differences between different fractions. In addition, we observed various unique transcripts, and novel predicted miRNAs among these plasma fractions. These results demonstrate that the distribution of RNA species and functional RNA transcripts is plasma fraction-dependent. Appropriate plasma preparation and thorough inspection of different plasma fractions are necessary for an exRNA-based biomarker study.

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.

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).

Systematic evaluation of multiple qPCR platforms, NanoString and miRNA-Seq for microRNA biomarker discovery in human biofluids

Aberrant miRNA expression has been associated with many diseases, and extracellular miRNAs that circulate in the bloodstream are remarkably stable. Recently, there has been growing interest in identifying cell-free circulating miRNAs that can serve as non-invasive biomarkers for early detection of disease or selection of treatment options. However, quantifying miRNA levels in biofluids is technically challenging due to their low abundance. Using reference samples, we performed a cross-platform evaluation in which miRNA profiling was performed on four different qPCR platforms (MiRXES, Qiagen, Applied Biosystems, Exiqon), nCounter technology (NanoString), and miRNA-Seq. Overall, our results suggest that using miRNA-Seq for discovery and targeted qPCR for validation is a rational strategy for miRNA biomarker development in clinical samples that involve limited amounts of biofluids.

Measurements Methods for the Development of MicroRNA-Based Tests for Cancer Diagnosis

Studies investigating microRNAs as potential biomarkers for cancer, immune-related diseases, or cardiac pathogenic diseases, among others, have exponentially increased in the last years. In particular, altered expression of specific miRNAs correlates with the occurrence of several diseases, making these molecules potential molecular tools for non-invasive diagnosis, prognosis, and response to therapy. Nonetheless, microRNAs are not in clinical use yet, due to inconsistencies in the literature regarding the specific miRNAs identified as biomarkers for a specific disease, which in turn can be attributed to several reasons, including lack of assay standardization and reproducibility. Technological limitations in circulating microRNAs measurement have been, to date, the biggest challenge for using these molecules in clinical settings. In this review we will discuss pre-analytical, analytical, and post-analytical challenges to address the potential technical biases and patient-related parameters that can have an influence and should be improved to translate miRNA biomarkers to the clinical stage. Moreover, we will describe the currently available methods for circulating miRNA expression profiling and measurement, underlining their advantages and potential pitfalls.

A Small RNA-Seq Protocol with Less Bias and Improved Capture of 2'-O-Methyl RNAs

The study of small RNAs (sRNAs) by next-generation sequencing (NGS) is challenged by bias issues during library preparation. Several types of sRNAs such as plant microRNAs (miRNAs) carry a 2'-O-methyl (2'-OMe) modification at their 3' terminal nucleotide. This modification adds another level of difficulty as it inhibits 3' adapter ligation. We previously demonstrated that modified versions of the "TruSeq (TS)" protocol have less bias and an improved detection of 2'-OMe RNAs. Here we describe in detail protocol "TS5," which showed the best overall performance. We also provide guidelines for bioinformatics analysis of the sequencing data.

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.