Improving small RNA-seq by using a synthetic spike-in set for size-range quality control together with a set for data normalization

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.

Protocol for a multicentre cross-sectional, longitudinal ambulatory clinical trial in rheumatoid arthritis and Parkinson's disease patients analysing the relation between the gut microbiome, fasting and immune status in Germany (ExpoBiome)

INTRODUCTION

Chronic inflammatory diseases like rheumatoid arthritis (RA) and neurodegenerative disorders like Parkinson's disease (PD) have recently been associated with a decreased diversity in the gut microbiome, emerging as key driver of various diseases. The specific interactions between gut-borne microorganisms and host pathophysiology remain largely unclear. The microbiome can be modulated by interventions comprising nutrition.The aim of our clinical study is to (1) examine effects of prolonged fasting (PF) and time-restricted eating (TRE) on the outcome parameters and the immunophenotypes of RA and PD with (2) special consideration of microbial taxa and molecules associated with changes expected in (1), and (3) identify factors impacting the disease course and treatment by in-depth screening of microorganisms and molecules in personalised HuMiX gut-on-chip models, to identify novel targets for anti-inflammatory therapy.

METHODS AND ANALYSIS

This trial is an open-label, multicentre, controlled clinical trial consisting of a cross-sectional and a longitudinal study. A total of 180 patients is recruited. For the cross-sectional study, 60 patients with PD, 60 patients with RA and 60 healthy controls are recruited at two different, specialised clinical sites. For the longitudinal part, 30 patients with PD and 30 patients with RA undergo 5-7 days of PF followed by TRE (16:8) for a period of 12 months. One baseline visit takes place before the PF intervention and 10 follow-up visits will follow over a period of 12 months (April 2021 to November 2023).

ETHICS AND DISSEMINATION

Ethical approval was obtained to plan and conduct the trial from the institutional review board of the Charité-Universitätsmedizin Berlin (EA1/204/19), the ethics committee of the state medical association (Landesärztekammer) of Hessen (2021-2230-zvBO) and the Ethics Review Panel (ERP) of the University of Luxembourg (ERP 21-001 A ExpoBiome). The results of this study will be disseminated through peer-reviewed publications, scientific presentations and social media.

TRIAL REGISTRATION NUMBER

NCT04847011.

Prostate cancer bone metastases biology and clinical management (Review)

Prostate cancer (PCa) is one of the most prominent causes of cancer-related mortality in the male population. A highly impactful prognostic factor for patients diagnosed with PCa is the presence or absence of bone metastases. The formation of secondary tumours at the bone is the most commonly observed site for the establishment of PCa metastases and is associated with reduced survival of patients in addition to a cohort of life-debilitating symptoms, including mobility issues and chronic pain. Despite the prevalence of this disease presentation and the high medical relevance of bone metastases, the mechanisms underlying the formation of metastases to the bone and the understanding of what drives the osteotropism exhibited by prostate tumours remain to be fully elucidated. This lack of in-depth understanding manifests in limited effective treatment options for patients with advanced metastatic PCa and culminates in the low rate of survival observed for this sub-set of patients. The present review aims to summarise the most recent promising advances in the understanding of how and why prostate tumours metastasise to the bone, with the ultimate aim of highlighting novel treatment and prognostic targets, which may provide the opportunity to improve the diagnosis and treatment of patients with PCa with bone metastases.

Cloning and Sequencing Eukaryotic Small RNAs

Small RNAs are ubiquitous regulators of gene expression that participate in nearly all aspects of physiology in a wide range of organisms. There are many different classes of eukaryotic small RNAs that play regulatory roles at every level of gene expression, including transcription, RNA stability, and translation. While eukaryotic small RNAs display diverse functions across and within classes, they are generally grouped functionally based on the machinery required for their biogenesis, the effector proteins they associate with, and their molecular characteristics. The development of techniques to clone and sequence small RNAs has been critical for their identification, yet the ligation-dependent addition of RNA adapters and the use of reverse transcriptase to generate cDNA in traditional library preparation protocols can be unsuitable to detect certain small RNA subtypes. In particular, 3' or 5' chemical modifications that are characteristic of specific types of small RNAs can impede the ligation-dependent addition of RNA adapters, while internal RNA modifications can interfere with accurate reverse transcription. The inability to clone certain small RNA subtypes with traditional protocols results in an inaccurate assessment of small RNA abundance and diversity, where some RNAs appear over-represented and others are not detected. This overview aims to guide users on how to design small RNA cloning workflows in eukaryotes to more accurately capture specific small RNAs of interest. Hence, we discuss the molecular biology underlying the identification and quantitation of small RNAs, explore the limitations of commonly used protocols, and detail the alternative approaches that can be used to enrich specific small RNA classes. © 2022 Wiley Periodicals LLC.

Using miRNAs as diagnostic biomarkers for male infertility: opportunities and challenges

The non-coding genome has been extensively studied for its role in human development and diseases. MicroRNAs (miRNAs) are small non-coding RNAs, which can regulate the expression of hundreds of genes at the post-transcriptional level. Therefore, any defects in miRNA biogenesis or processing can affect the genes and have been linked to several diseases. Male infertility is a clinical disorder with a significant number of cases being idiopathic. Problems in spermatogenesis and epididymal maturation, testicular development, sperm maturation or migration contribute to male infertility, and many of these idiopathic cases are related to issues with the miRNAs which tightly regulate these processes. This review summarizes the recent research on various such miRNAs and puts together the candidate miRNAs that may be used as biomarkers for diagnosis. The development of strategies for male infertility treatment using anti-miRs or miRNA mimics is also discussed. Although promising, the development of miRNA diagnostics and therapeutics is challenging, and ways to overcome some of these challenges are also reviewed.

Profile of Small RNAs, vDNA Forms and Viral Integrations in Late Chikungunya Virus Infection of Aedes albopictus Mosquitoes

The Asian tiger mosquito Aedes albopictus is contributing to the (re)-emergence of Chikungunya virus (CHIKV). To gain insights into the molecular underpinning of viral persistence, which renders a mosquito a life-long vector, we coupled small RNA and whole genome sequencing approaches on carcasses and ovaries of mosquitoes sampled 14 days post CHIKV infection and investigated the profile of small RNAs and the presence of vDNA fragments. Since Aedes genomes harbor nonretroviral Endogenous Viral Elements (nrEVEs) which confers tolerance to cognate viral infections in ovaries, we also tested whether nrEVEs are formed after CHIKV infection. We show that while small interfering (si)RNAs are evenly distributed along the full viral genome, PIWI-interacting (pi)RNAs mostly arise from a ~1000 bp window, from which a unique vDNA fragment is identified. CHIKV infection does not result in the formation of new nrEVEs, but piRNAs derived from existing nrEVEs correlate with differential expression of an endogenous transcript. These results demonstrate that all three RNAi pathways contribute to the homeostasis during the late stage of CHIKV infection, but in different ways, ranging from directly targeting the viral sequence to regulating the expression of mosquito transcripts and expand the role of nrEVEs beyond immunity against cognate viruses.

WIND (Workflow for pIRNAs aNd beyonD): a strategy for in-depth analysis of small RNA-seq data

Current bioinformatics workflows for PIWI-interacting RNA (piRNA) analysis focus primarily on germline-derived piRNAs and piRNA-clusters. Frequently, they suffer from outdated piRNA databases, questionable quantification methods, and lack of reproducibility. Often, pipelines specific to miRNA analysis are used for the piRNA research in silico. Furthermore, the absence of a well-established database for piRNA annotation, as for miRNA, leads to uniformity issues between studies and generates confusion for data analysts and biologists.

For these reasons, we have developed WIND ( Workflow for p IRNAs a Nd beyon D), a bioinformatics workflow that addresses the crucial issue of piRNA annotation, thereby allowing a reliable analysis of small RNA sequencing data for the identification of piRNAs and other small non-coding RNAs (sncRNAs) that in the past have been incorrectly classified as piRNAs. WIND allows the creation of a comprehensive annotation track of sncRNAs combining information available in RNAcentral, with piRNA sequences from piRNABank, the first database dedicated to piRNA annotation. WIND was built with Docker containers for reproducibility and integrates widely used bioinformatics tools for sequence alignment and quantification. In addition, it includes Bioconductor packages for exploratory data and differential expression analysis. Moreover, WIND implements a "dual" approach for the evaluation of sncRNAs expression level quantifying the aligned reads to the annotated genome and carrying out an alignment-free transcript quantification using reads mapped to the transcriptome. Therefore, a broader range of piRNAs can be annotated, improving their quantification and easing the subsequent downstream analysis. WIND performance has been tested with several small RNA-seq datasets, demonstrating how our approach can be a useful and comprehensive resource to analyse piRNAs and other classes of sncRNAs.

Charting Extracellular Transcriptomes in The Human Biofluid RNA Atlas

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

Targets, pitfalls and reference materials for liquid biopsy tests in cancer diagnostics

Assessment of cell free DNA (cfDNA) and RNA (cfRNA), circulating tumor cells (CTC) and extracellular vesicles (EV) in blood or other bodily fluids can enable early cancer detection, tumor dynamics assessment, minimal residual disease detection and therapy monitoring. However, few liquid biopsy tests progress towards clinical application because results are often discordant and challenging to reproduce. Reproducibility can be enhanced by the development and implementation of standard operating procedures and reference materials to identify and correct for pre-analytical variables. In this review we elaborate on the technological considerations, pre-analytical variables and the use and availability of reference materials for the assessment of liquid biopsy targets in blood and highlight initiatives towards the standardization of liquid biopsy testing.

Molecular and Functional Characterization of the Somatic PIWIL1/piRNA Pathway in Colorectal Cancer Cells

PIWI-like (PIWIL) proteins and small non-coding piRNAs, involved in genome regulation in germline cells, are found aberrantly expressed in human tumors. Gene expression data from The Cancer Genome Atlas (TCGA), the Genotype-Tissue Expression (GTEx) project, and the European Genome-Phenome Archive (EGA) indicate that the PIWIL1 gene is ectopically activated in a significant fraction of colorectal cancers (CRCs), where this is accompanied by promoter demethylation, together with germline factors required for piRNA production. Starting from this observation, the PIWIL/piRNA pathway was studied in detail in COLO 205 CRC cells, which express significant levels of this protein, to investigate role and significance of ectopic PIWIL1 expression in human tumors. RNA sequencing and cell and computational biology led to the demonstration that PIWIL1 localizes in a nuage-like structure located in the perinuclear region of the cell and that a significant fraction of the piRNAs expressed in these cells are methylated, and, therefore, present in an active form. This was further supported by RNA immunoprecipitation, which revealed how several piRNAs can be found loaded into PIWIL1 to form complexes also comprising their target mRNAs. The mature transcripts associated with the PIWIL-piRNA complex encode key regulatory proteins involved in the molecular mechanisms sustaining colorectal carcinogenesis, suggesting that the PIWI/piRNA pathway may actively contribute to the establishment and/or maintenance of clinico-pathological features of CRCs.

Bias-minimized quantification of microRNA reveals widespread alternative processing and 3' end modification

MicroRNAs (miRNAs) modulate diverse biological and pathological processes via post-transcriptional gene silencing. High-throughput small RNA sequencing (sRNA-seq) has been widely adopted to investigate the functions and regulatory mechanisms of miRNAs. However, accurate quantification of miRNAs has been limited owing to the severe ligation bias in conventional sRNA-seq methods. Here, we quantify miRNAs and their variants (known as isomiRs) by an improved sRNA-seq protocol, termed AQ-seq (accurate quantification by sequencing), that utilizes adapters with terminal degenerate sequences and a high concentration of polyethylene glycol (PEG), which minimize the ligation bias during library preparation. Measurement using AQ-seq allows us to correct the previously misannotated 5' end usage and strand preference in public databases. Importantly, the analysis of 5' terminal heterogeneity reveals widespread alternative processing events which have been underestimated. We also identify highly uridylated miRNAs originating from the 3p strands, indicating regulations mediated by terminal uridylyl transferases at the pre-miRNA stage. Taken together, our study reveals the complexity of the miRNA isoform landscape, allowing us to refine miRNA annotation and to advance our understanding of miRNA regulation. Furthermore, AQ-seq can be adopted to improve other ligation-based sequencing methods including crosslinking-immunoprecipitation-sequencing (CLIP-seq) and ribosome profiling (Ribo-seq).

Characterization of Chinese hamster ovary cells with disparate chromosome numbers: Reduction of the amount of mRNA relative to total protein

Chromosomes in Chinese hamster ovary (CHO) cells are labile. We have shown that high-chromosome-number CHO cells have greater potential to become robust producers of recombinant proteins. One explanation being the increase in transgene integration sites. However, high-chromosome-number cell clones produce more IgG3 following culture of single-cell clones, even under conditions that yield the same number of integrations as cells with normal chromosome numbers. Here, we characterized high-chromosome-number cells by transcriptome analysis. RNA standards were used to normalize transcriptomes of cells that had different chromosome numbers. Our results demonstrate that the mRNA ratio of β-actin and many other genes in high-chromosome-number cells to that in normal-chromosome-number cells per cell (normalized to RNA standards) was smaller than the equivalent genomic size and cell volume ratios. Many genes encoding membrane proteins are more highly expressed in high-chromosome-number cells, probably due to differences in cell size caused by the increase in chromosomes. In addition, genes related to histone modification and lipid metabolism are differentially expressed. The reduced transcript level required per protein produced in total and the different intracellular signal transductions might be key factors for antibody production.

Identifying and characterizing functional 3' nucleotide addition in the miRNA pathway

Over the past decade, modifications to microRNAs (miRNAs) via 3' end nucleotide addition have gone from a deep-sequencing curiosity to experimentally confirmed drivers of a range of regulatory activities. Here we overview the methods that have been deployed by researchers seeking to untangle these diverse functional roles and include characterizing not only the nucleotidyl transferases catalyzing the additions but also the nucleotides being added, and the timing of their addition during the miRNA pathway. These methods and their further development are key to clarifying the diverse and sometimes contradictory functional findings presently attributed to these nucleotide additions.

miRge 2.0 for comprehensive analysis of microRNA sequencing data

Background

miRNAs play important roles in the regulation of gene expression. The rapidly developing field of microRNA sequencing (miRNA-seq; small RNA-seq) needs comprehensive, robust, user-friendly and standardized bioinformatics tools to analyze these large datasets. We present miRge 2.0, in which multiple enhancements were made towards these goals.

Results

miRge 2.0 has become more comprehensive with increased functionality including a novel miRNA detection method, A-to-I editing analysis, integrated standardized GFF3 isomiR reporting, and improved alignment to miRNAs. The novel miRNA detection method uniquely uses both miRNA hairpin sequence structure and composition of isomiRs resulting in higher specificity for potential miRNA identification. Using known miRNA data, our support vector machine (SVM) model predicted miRNAs with an average Matthews correlation coefficient (MCC) of 0.939 over 32 human cell datasets and outperformed miRDeep2 and miRAnalyzer regarding phylogenetic conservation. The A-to-I editing detection strongly correlated with a reference dataset with adjusted R² = 0.96. miRge 2.0 is the most up-to-date aligner with custom libraries to both miRBase v22 and MirGeneDB v2.0 for 6 species: human, mouse, rat, fruit fly, nematode and zebrafish; and has a tool to create custom libraries. For user-friendliness, miRge 2.0 is incorporated into bcbio-nextgen and implementable through Bioconda.

Conclusions

miRge 2.0 is a redesigned, leading miRNA RNA-seq aligner with several improvements and novel utilities. miRge 2.0 is freely available at: https://github.com/mhalushka/miRge.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2287-y) contains supplementary material, which is available to authorized users.

Importance of miRNA stability and alternative primary miRNA isoforms in gene regulation during Drosophila development

Mature microRNAs (miRNAs) are processed from primary transcripts (pri-miRNAs), and their expression is controlled at transcriptional and post-transcriptional levels. However, how regulation at multiple levels achieves precise control remains elusive. Using published and new datasets, we profile a time course of mature and pri-miRNAs in Drosophila embryos and reveal the dynamics of miRNA production and degradation as well as dynamic changes in pri-miRNA isoform selection. We found that 5' nucleotides influence stability of mature miRNAs. Furthermore, distinct half-lives of miRNAs from the mir-309 cluster shape their temporal expression patterns, and the importance of rapid degradation of the miRNAs in gene regulation is detected as distinct evolutionary signatures at the target sites in the transcriptome. Finally, we show that rapid degradation of miR-3/-309 may be important for regulation of the planar cell polarity pathway component Vang. Altogether, the results suggest that complex mechanisms regulate miRNA expression to support normal development.

Approaches to Detect microRNA Expression in T Cell Subsets and T Cell Differentiation

MicroRNAs (miRNAs) are crucial components of the molecular networks regulating differentiation and responses of T lymphocytes in health and disease. It is therefore essential to rely on robust methods of qualitative and quantitative investigation of miRNA expression in T cell subsets, and during T cell activation and differentiation. Here, we focus on different methods for miRNA analysis, including Northern blots, quantitative RT-PCR, and next-generation sequencing, and we discuss advantages and disadvantages of each method. While we mainly focus on the study of miRNA expression in human T lymphocytes, these methods can also be applied to other species and/or cell types.

Novel small RNA spike-in oligonucleotides enable absolute normalization of small RNA-Seq data

Normalization of high-throughput small RNA sequencing (sRNA-Seq) data is required to compare sRNA levels across different samples. Commonly used relative normalization approaches can cause erroneous conclusions due to fluctuating small RNA populations between tissues. We developed a set of sRNA spike-in oligonucleotides (sRNA spike-ins) that enable absolute normalization of sRNA-Seq data across independent experiments, as well as the genome-wide estimation of sRNA:mRNA stoichiometries when used together with mRNA spike-in oligonucleotides.

Formidable challenges to the notion of biologically important roles for dietary small RNAs in ingesting mammals

The notion of uptake of active diet-derived small RNAs (sRNAs) in recipient organisms could have significant implications for our understanding of oral therapeutics and nutrition, for the safe use of RNA interference (RNAi) in agricultural biotechnology, and for ecological relationships. Yet, the transfer and subsequent regulation of gene activity by diet-derived sRNAs in ingesting mammals are still heavily debated. Here, we synthesize current information based on multiple independent studies of mammals, invertebrates, and plants. Rigorous assessment of these data emphasize that uptake of active dietary sRNAs is neither a robust nor a prevalent mechanism to maintain steady-state levels in higher organisms. While disagreement still continues regarding whether such transfer may occur in specialized contexts, concerns about technical difficulties and a lack of consensus on appropriate methods have led to questions regarding the reproducibility and biologic significance of some seemingly positive results. For any continuing investigations, concerted efforts should be made to establish a strong mechanistic basis for potential effects of dietary sRNAs and to agree on methodological guidelines for realizing such proof. Such processes would ensure proper interpretation of studies aiming to prove dietary sRNA activity in mammals and inform potential for application in therapeutics and agriculture.

Tissue dual RNA-seq allows fast discovery of infection-specific functions and riboregulators shaping host-pathogen transcriptomes

Pathogenic bacteria need to rapidly adjust their virulence and fitness program to prevent eradication by the host. So far, underlying adaptation processes that drive pathogenesis have mostly been studied in vitro, neglecting the true complexity of host-induced stimuli acting on the invading pathogen. In this study, we developed an unbiased experimental approach that allows simultaneous monitoring of genome-wide infection-linked transcriptional alterations of the host and colonizing extracellular pathogens. Using this tool for Yersinia pseudotuberculosis-infected lymphatic tissues, we revealed numerous alterations of host transcripts associated with inflammatory and acute-phase responses, coagulative activities, and transition metal ion sequestration, highlighting that the immune response is dominated by infiltrating neutrophils and elicits a mixed T_H17/T_H1 response. In consequence, the pathogen's response is mainly directed to prevent phagocytic attacks. Yersinia up-regulates the gene and expression dose of the antiphagocytic type III secretion system (T3SS) and induces functions counteracting neutrophil-induced ion deprivation, radical stress, and nutritional restraints. Several conserved bacterial riboregulators were identified that impacted this response. The strongest influence on virulence was found for the loss of the carbon storage regulator (Csr) system, which is shown to be essential for the up-regulation of the T3SS on host cell contact. In summary, our established approach provides a powerful tool for the discovery of infection-specific stimuli, induced host and pathogen responses, and underlying regulatory processes.

Cajal body function in genome organization and transcriptome diversity

Nuclear bodies contribute to non-random organization of the human genome and nuclear function. Using a major prototypical nuclear body, the Cajal body, as an example, we suggest that these structures assemble at specific gene loci located across the genome as a result of high transcriptional activity. Subsequently, target genes are physically clustered in close proximity in Cajal body-containing cells. However, Cajal bodies are observed in only a limited number of human cell types, including neuronal and cancer cells. Ultimately, Cajal body depletion perturbs splicing kinetics by reducing target small nuclear RNA (snRNA) transcription and limiting the levels of spliceosomal snRNPs, including their modification and turnover following each round of RNA splicing. As such, Cajal bodies are capable of shaping the chromatin interaction landscape and the transcriptome by influencing spliceosome kinetics. Future studies should concentrate on characterizing the direct influence of Cajal bodies upon snRNA gene transcriptional dynamics. Also see the video abstract here.

Toward the promise of microRNAs - Enhancing reproducibility and rigor in microRNA research

The fields of applied and translational microRNA research have exploded in recent years as microRNAs have been implicated across a spectrum of diseases. MicroRNA biomarkers, microRNA therapeutics, microRNA regulation of cellular physiology and even xenomiRs have stimulated great interest, which have brought many researchers into the field. Despite many successes in determining general mechanisms of microRNA generation and function, the application of microRNAs in translational areas has not had as much success. It has been a challenge to localize microRNAs to a given cell type within tissues and assay them reliably. At supraphysiologic levels, microRNAs may regulate hosts of genes that are not the physiologic biochemical targets. Thus the applied and translational microRNA literature is filled with pitfalls and claims that are neither scientifically rigorous nor reproducible. This review is focused on increasing awareness of the challenges of working with microRNAs in translational research and recommends better practices in this area of discovery.

Toward reliable biomarker signatures in the age of liquid biopsies - how to standardize the small RNA-Seq workflow

Small RNA-Seq has emerged as a powerful tool in transcriptomics, gene expression profiling and biomarker discovery. Sequencing cell-free nucleic acids, particularly microRNA (miRNA), from liquid biopsies additionally provides exciting possibilities for molecular diagnostics, and might help establish disease-specific biomarker signatures. The complexity of the small RNA-Seq workflow, however, bears challenges and biases that researchers need to be aware of in order to generate high-quality data. Rigorous standardization and extensive validation are required to guarantee reliability, reproducibility and comparability of research findings. Hypotheses based on flawed experimental conditions can be inconsistent and even misleading. Comparable to the well-established MIQE guidelines for qPCR experiments, this work aims at establishing guidelines for experimental design and pre-analytical sample processing, standardization of library preparation and sequencing reactions, as well as facilitating data analysis. We highlight bottlenecks in small RNA-Seq experiments, point out the importance of stringent quality control and validation, and provide a primer for differential expression analysis and biomarker discovery. Following our recommendations will encourage better sequencing practice, increase experimental transparency and lead to more reproducible small RNA-Seq results. This will ultimately enhance the validity of biomarker signatures, and allow reliable and robust clinical predictions.

Small RNAs from Bemisia tabaci Are Transferred to Solanum lycopersicum Phloem during Feeding

The phloem-feeding whitefly Bemisia tabaci is a serious pest to a broad range of host plants, including many economically important crops such as tomato. These insects serve as a vector for various devastating plant viruses. It is known that whiteflies are capable of manipulating host-defense responses, potentially mediated by effector molecules in the whitefly saliva. We hypothesized that, beside putative effector proteins, small RNAs (sRNA) are delivered by B. tabaci into the phloem, where they may play a role in manipulating host plant defenses. There is already evidence to suggest that sRNAs can mediate the host-pathogen dialogue. It has been shown that Botrytis cinerea, the causal agent of gray mold disease, takes advantage of the plant sRNA machinery to selectively silence host genes involved in defense signaling. Here we identified sRNAs originating from B. tabaci in the phloem of tomato plants on which they are feeding. sRNAs were isolated and sequenced from tomato phloem of whitefly-infested and control plants as well as from the nymphs themselves, control leaflets, and from the infested leaflets. Using stem-loop RT-PCR, three whitefly sRNAs have been verified to be present in whitefly-infested leaflets that were also present in the whitefly-infested phloem sample. Our results show that whitefly sRNAs are indeed present in tomato tissues upon feeding, and they appear to be mobile in the phloem. Their role in the host-insect interaction can now be investigated.