An efficient and sensitive method for preparing cDNA libraries from scarce biological samples

Integrated transcriptome and microbiome analyses of residual feed intake in ducks during high production period

Residual feed intake (RFI) is a crucial parameter for assessing the feeding efficiency of poultry. Minimizing RFI can enhance feed utilization and reduce costs. In this study, 315 healthy female ducks were individually housed in cages. Growth performance was monitored during the high laying period, from 290 to 325 d of age. The cecal transcriptome and microbiome of 12 ducks with high RFI and 12 with low residual feed intake (LRFI) were analyzed. Regarding growth performance, the LRFI group exhibited significantly lower RFI, feed conversion ratio (FCR), and feed intake (Fi) compared to the HRFI group (p < 0.01). However, there were no significant differences observed in body weight (BW), body weight gain (BWG), and egg mass (EML) between the groups (p > 0.05). Microbiome analysis demonstrated that RFI impacted gut microbial abundance, particularly affecting metabolism and disease-related microorganisms such as Romboutsia, Enterococcus, and Megamonas funiformis. Transcriptome analysis revealed that varying RFI changed the expression of genes related to glucose metabolism and lipid metabolism, including APOA1, G6PC1, PCK1, and PLIN1. The integrated analysis indicated that host genes were closely linked to the microbiota and primarily function in lipid metabolism, which may enhance feeding efficiency by influencing metabolism and maintaining gut homeostasis.

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.

Integrated transcriptome and metabolome profiling of Camellia reticulata reveal mechanisms of flower color differentiation

Camellia reticulata (Lindl.) is an important ornamental plant in China. Long-term natural or artificial selections have resulted in diverse phenotypes, especially for flower colors. Modulating flower colors can enhance the visual appeal and economic value in ornamental plants. In this study, we investigated the molecular mechanisms underlying flower color differentiation in C. reticulata. We performed a combined transcriptome and metabolome analysis of the petals of a popular variety C. reticulata (HHYC) (red), and its two cultivars "Xuejiao" (XJ) (pink) and "Tongzimian" (TZM) (white). Targeted metabolome profiling identified 310 flavonoid compounds of which 18 anthocyanins were differentially accumulated among the three samples with an accumulation pattern of HHYC > XJ > TZM. Likewise, transcriptome analysis showed that carotenoid and anthocyanin biosynthetic structural genes were mostly expressed in order of HHYC > XJ > TZM. Two genes (gene-LOC114287745765 and gene-LOC114289234) encoding for anthocyanidin 3-O-glucosyltransferase are predicted to be responsible for red coloration in HHYC and XJ. We also detected 42 MYB and 29 bHLH transcription factors as key regulators of anthocyanin-structural genes. Overall, this work showed that flavonoids, particularly anthocyanins contents are the major determinants of flower color differentiation among the 3 C. reticulata samples. In addition, the main regulatory and structural genes modulating anthocyanin contents in C. reticulata have been unveiled. Our results will help in the development of Camellia varieties with specific flower color and quality.

Interactions of Muscovy duck reovirus, gut microbiota, and host innate immunity: Transcriptome and gut microbiota analysis

It has been shown that Muscovy duck reovirus (MDRV) infection causes severe intestinal barrier damage and intestinal mucosal immune suppression. The health and balance of gut microbes is essential for the progression of intestinal infectious diseases. To investigate the interaction of MDRV, intestinal bacteria with host intestinal innate immunity, an MDRV contact-infection model was established in this study. High-throughput sequencing technology was used to sequence 16S rDNA and transcripts in ileal samples from experimental Muscovy ducklings. Our results suggest that intestinal opportunistic pathogens such as Streptococcus and Corynebacterium proliferated massively in MDRV-infected Muscovy ducklings. The body initiates antiviral and antibacterial immunity and actively fights the infection of gut microbes. The synthesis of peptidoglycan, lipopolysaccharide, and flagellin by intestinal bacteria activates the Toll-like receptor signaling pathway resulting in increased secretion of IFN-β, IL-1β, and IL-8. The RIG-I-like receptor signaling pathway is an important signaling pathway for the interaction between MDRV and the host. At the same time, we also observed that multiple genes in the JAK-STAT signaling pathway were significantly different. These genes are important targets for studying the immunosuppression caused by MDRV. In conclusion, we analyzed the interaction of MDRV, intestinal flora and host immune system during MDRV infection, which provides a basis for the further study on the mechanism of intestinal immunosuppression caused by MDRV.

RIPiT-Seq: A tandem immunoprecipitation approach to reveal global binding landscape of multisubunit ribonucleoproteins

RNA-binding proteins (RBPs) regulate all aspects of RNA metabolism. The ability to identify RNA targets bound by RBPs is critical for understanding RBP function. While powerful techniques are available to identify binding sites of individual RBPs at high resolution, it remains challenging to unravel binding sites of multicomponent ribonucleoproteins (RNPs) where multiple RBPs or proteins function cooperatively to bind to target RNAs. To fill this gap, we have previously developed RNA Immunoprecipitation in Tandem followed by high-throughput sequencing (RIPiT-seq) to characterize RNA targets of compositionally distinct RNP complexes by sequentially immunoprecipitating two proteins from the same RNP and sequencing the co-purifying RNA footprints. Here, we provide an updated and improved protocol for RIPiT-seq. In this protocol, we have used CRISPR-Cas9 to introduce affinity tag to endogenous protein of interest to capture a more representative state of an RNP complex. We present a modified protocol for library preparation for high-throughput sequencing so that it exclusively uses equipment and reagents available in a standard molecular biology lab. This updated custom library preparation protocol is compatible with commercial PCR multiplexing systems for Illumina sequencing platform for simultaneous and cost-effective analysis of large number of samples.

Extracellular microRNAs in human circulation are associated with miRISC complexes that are accessible to anti-AGO2 antibody and can bind target mimic oligonucleotides

MicroRNAs (miRNAs) function cell-intrinsically to regulate gene expression by base-pairing to complementary mRNA targets while in association with Argonaute, the effector protein of the miRNA-mediated silencing complex (miRISC). A relatively dilute population of miRNAs can be found extracellularly in body fluids such as human blood plasma and cerebrospinal fluid (CSF). The remarkable stability of circulating miRNAs in such harsh extracellular environments can be attributed to their association with protective macromolecular complexes, including extracellular vesicles (EVs), proteins such as Argonaut 2 (AGO2), or high-density lipoproteins. The precise origins and the potential biological significance of various forms of miRNA-containing extracellular complexes are poorly understood. It is also not known whether extracellular miRNAs in their native state may retain the capacity for miRISC-mediated target RNA binding. To explore the potential functionality of circulating extracellular miRNAs, we comprehensively investigated the association between circulating miRNAs and the miRISC Argonaute AGO2. Using AGO2 immunoprecipitation (IP) followed by small-RNA sequencing, we find that miRNAs in circulation are primarily associated with antibody-accessible miRISC/AGO2 complexes. Moreover, we show that circulating miRNAs can base-pair with a target mimic in a seed-based manner, and that the target-bound AGO2 can be recovered from blood plasma in an ∼1:1 ratio with the respective miRNA. Our findings suggest that miRNAs in circulation are largely contained in functional miRISC/AGO2 complexes under normal physiological conditions. However, we find that, in human CSF, the assortment of certain extracellular miRNAs into free miRISC/AGO2 complexes can be affected by pathological conditions such as amyotrophic lateral sclerosis.

Transcriptome analysis of the effects of Hericium erinaceus polysaccharide on the lymphocyte homing in Muscovy duck reovirus-infected ducklings

Hericium erinaceus polysaccharide (HEP) is a bioactive substance present in the fruiting bodies of H. erinaceus. Previously we have shown that HEP can repair the intestinal injury caused by Muscovy duck reovirus (MDRV) infection in Muscovy ducklings. To examine the effect of HEP on intestine mucosal MDRV immunity and explore its possible mechanisms, an MDRV contact-infection model in the Muscovy ducklings was established. Transcriptome sequencing analysis was then performed to investigate the mechanism of action of HEP on intestine mucosal MDRV immunity. During the infection, the expression levels of genes involved in cellular activities (protein translation and binding, cytokine interaction, and adhesion molecules activities) in the infected ducklings were increased. The expression levels of adhesion molecules (α4β7, LFA-1) and chemotaxis cytokine receptors (CCR7, CCR9, and CCR10) were also significantly upregulated. Following HEP treatment, cellular activities and cytokines upregulated to various degrees play crucial roles in the immune defenses and antiviral activities of Muscovy ducklings. ELISA analysis results were consistent with the results of the transcriptome analysis. Overall, our results provide a basis for further studying the underlying mechanisms of HEP in regulating mucosal immunity and for the clinical application of HEP in controlling MDRV infection in the Muscovy duck industry.

Decreasing miRNA sequencing bias using a single adapter and circularization approach

The ability to accurately quantify all the microRNAs (miRNAs) in a sample is important for understanding miRNA biology and for development of new biomarkers and therapeutic targets. We develop a new method for preparing miRNA sequencing libraries, RealSeq®-AC, that involves ligating the miRNAs with a single adapter and circularizing the ligation products. When compared to other methods, RealSeq®-AC provides greatly reduced miRNA sequencing bias and allows the identification of the largest variety of miRNAs in biological samples. This reduced bias also allows robust quantification of miRNAs present in samples across a wide range of RNA input levels.

CleanTag Adapters Improve Small RNA Next-Generation Sequencing Library Preparation by Reducing Adapter Dimers

Next-generation small RNA sequencing is a valuable tool which is increasing our knowledge regarding small noncoding RNAs and their function in regulating genetic information. Library preparation protocols for small RNA have thus far been restricted due to higher RNA input requirements (>10 ng), long workflows, and tedious manual gel purifications. Small RNA library preparation methods focus largely on the prevention or depletion of a side product known as adapter dimer that tends to dominate the reaction. Adapter dimer is the ligation of two adapters to one another without an intervening library RNA insert or any useful sequencing information. The amplification of this side reaction is favored over the amplification of tagged library since it is shorter. The small size discrepancy between these two species makes separation and purification of the tagged library very difficult. Adapter dimer hinders the use of low input samples and the ability to automate the workflow so we introduce an improved library preparation protocol which uses chemically modified adapters (CleanTag) to significantly reduce the adapter dimer. CleanTag small RNA library preparation workflow decreases adapter dimer to allow for ultra-low input samples (down to approx. 10 pg total RNA), elimination of the gel purification step, and automation. We demonstrate how to carry out this streamlined protocol to improve NGS data quality and allow for the use of sample types with limited RNA material.

Cloning and Identification of Recombinant Argonaute-Bound Small RNAs Using Next-Generation Sequencing

Argonaute proteins (AGOs) are loaded with small RNAs as guides to recognize target mRNAs. Since the target specificity heavily depends on the base complementarity between two strands, it is important to identify small guide and long target RNAs bound to AGOs. For this purpose, next-generation sequencing (NGS) technologies have extended our appreciation truly to the nucleotide level. However, the identification of RNAs via NGS from scarce RNA samples remains a challenge. Further, most commercial and published methods are compatible with either small RNAs or long RNAs, but are not equally applicable to both. Therefore, a single method that yields quantitative, bias-free NGS libraries to identify small and long RNAs from low levels of input will be of wide interest. Here, we introduce such a procedure that is based on several modifications of two published protocols and allows robust, sensitive, and reproducible cloning and sequencing of small amounts of RNAs of variable lengths. The method was applied to the identification of small RNAs bound to a purified eukaryotic AGO. Following ligation of a DNA adapter to RNA 3'-end, the key feature of this method is to use the adapter for priming reverse transcription (RT) wherein biotinylated deoxyribonucleotides specifically incorporated into the extended complementary DNA. Such RT products are enriched on streptavidin beads, circularized while immobilized on beads and directly used for PCR amplification. We provide a stepwise guide to generate RNA-Seq libraries, their purification, quantification, validation, and preparation for next-generation sequencing. We also provide basic steps in post-NGS data analyses using Galaxy, an open-source, web-based platform.

Structure-seq2: sensitive and accurate genome-wide profiling of RNA structure in vivo

RNA serves many functions in biology such as splicing, temperature sensing, and innate immunity. These functions are often determined by the structure of RNA. There is thus a pressing need to understand RNA structure and how it changes during diverse biological processes both in vivo and genome-wide. Here, we present Structure-seq2, which provides nucleotide-resolution RNA structural information in vivo and genome-wide. This optimized version of our original Structure-seq method increases sensitivity by at least 4-fold and improves data quality by minimizing formation of a deleterious by-product, reducing ligation bias, and improving read coverage. We also present a variation of Structure-seq2 in which a biotinylated nucleotide is incorporated during reverse transcription, which greatly facilitates the protocol by eliminating two PAGE purification steps. We benchmark Structure-seq2 on both mRNA and rRNA structure in rice (Oryza sativa). We demonstrate that Structure-seq2 can lead to new biological insights. Our Structure-seq2 datasets uncover hidden breaks in chloroplast rRNA and identify a previously unreported N¹-methyladenosine (m¹A) in a nuclear-encoded Oryza sativa rRNA. Overall, Structure-seq2 is a rapid, sensitive, and unbiased method to probe RNA in vivo and genome-wide that facilitates new insights into RNA biology.

The in vivo genetic toolkit for studying expression and functions of Drosophila melanogaster microRNAs

Since the initial reports that a group of small RNAs, now known as microRNAs (miRNAs), regulates gene expression without being translated into proteins, there has been an explosion of studies on these important expression modulators. Drosophila melanogaster has proven to be one of the most amenable animal models for investigations of miRNA biogenesis and gene regulatory activities. Here, we highlight the publicly available genetic tools and strategies for in vivo functional studies of miRNA activity in D. melanogaster. By coupling genetic approaches using available strain libraries with technologies for miRNA expression analysis and target and pathway prediction, researchers' ability to test functional activities of miRNAs in vivo is now greatly enhanced. We also comment on the tools that need to be developed to aid in comprehensive evaluation of Drosophila miRNA activities that impact traits of interest.

TargetLink, a new method for identifying the endogenous target set of a specific microRNA in intact living cells

MicroRNAs are small non-coding RNAs acting as posttranscriptional repressors of gene expression. Identifying mRNA targets of a given miRNA remains an outstanding challenge in the field. We have developed a new experimental approach, TargetLink, that applied locked nucleic acid (LNA) as the affinity probe to enrich target genes of a specific microRNA in intact cells. TargetLink also consists a rigorous and systematic data analysis pipeline to identify target genes by comparing LNA-enriched sequences between experimental and control samples. Using miR-21 as a test microRNA, we identified 12 target genes of miR-21 in a human colorectal cancer cell by this approach. The majority of the identified targets interacted with miR-21 via imperfect seed pairing. Target validation confirmed that miR-21 repressed the expression of the identified targets. The cellular abundance of the identified miR-21 target transcripts varied over a wide range, with some targets expressed at a rather low level, confirming that both abundant and rare transcripts are susceptible to regulation by microRNAs, and that TargetLink is an efficient approach for identifying the target set of a specific microRNA in intact cells. C20orf111, one of the novel targets identified by TargetLink, was found to reside in the nuclear speckle and to be reliably repressed by miR-21 through the interaction at its coding sequence.

Plasma vesicle miRNAs for therapy response monitoring in Hodgkin lymphoma patients

BACKGROUND. Cell-free circulating nucleic acids, including 22-nt microRNAs (miRNAs), represent noninvasive biomarkers for treatment response monitoring of cancer patients. While the majority of plasma miRNA is bound to proteins, a smaller, less well-characterized pool is associated with extracellular vesicles (EVs). Here, we addressed whether EV-associated miRNAs reflect metabolic disease in classical Hodgkin lymphoma (cHL) patients.

METHODS. With standardized size-exclusion chromatography (SEC), we isolated EV-associated extracellular RNA (exRNA) fractions and protein-bound miRNA from plasma of cHL patients and healthy subjects. We performed a comprehensive small RNA sequencing analysis and validation by TaqMan qRT-PCR for candidate discovery. Fluorodeoxyglucose-PET (FDG-PET) status before treatment, directly after treatment, and during long-term follow-up was compared directly with EV miRNA levels.

RESULTS. The plasma EV miRNA repertoire was more extensive compared with protein-bound miRNA that was heavily dominated by a few abundant miRNA species and was less informative of disease status. Purified EV fractions of untreated cHL patients and tumor EVs had enriched levels of miR24-3p, miR127-3p, miR21-5p, miR155-5p, and let7a-5p compared with EV fractions from healthy subjects and disease controls. Serial monitoring of EV miRNA levels in patients before treatment, directly after treatment, and during long-term follow-up revealed robust, stable decreases in miRNA levels matching a complete metabolic response, as observed with FDG-PET. Importantly, EV miRNA levels rose again in relapse patients.

CONCLUSION. We conclude that cHL-related miRNA levels in circulating EVs reflect the presence of vital tumor tissue and are suitable for therapy response and relapse monitoring in individual cHL patients.

FUNDING. Cancer Center Amsterdam Foundation (CCA-2013), Dutch Cancer Society (KWF-5510), Technology Foundation STW (STW Perspectief CANCER-ID).

The extracellular RNA repertoire in circulating extracellular vesicles is useful indicator of therapy response and relapse in classical Hodgkin lymphoma patients.

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.

Staufen Negatively Modulates MicroRNA Activity in Caenorhabditis elegans

The double-stranded RNA-binding protein Staufen has been implicated in various posttranscriptional gene regulatory processes. Here, we demonstrate that the Caenorhabditis elegans homolog of Staufen, STAU-1, functionally interacts with microRNAs. Loss-of-function mutations of stau-1 significantly suppress phenotypes of let-7 family microRNA mutants, a hypomorphic allele of dicer, and a lsy-6 microRNA partial loss-of-function mutant. Furthermore, STAU-1 modulates the activity of lin-14, a target of lin-4 and let-7 family microRNAs, and this modulation is abolished when the 3' untranslated region of lin-14 is removed. Deep sequencing of small RNA cDNA libraries reveals no dramatic change in the levels of microRNAs or other small RNA populations between wild-type and stau-1 mutants, with the exception of certain endogenous siRNAs in the WAGO pathway. The modulation of microRNA activity by STAU-1 does not seem to be associated with the previously reported enhanced exogenous RNAi (Eri) phenotype of stau-1 mutants, since eri-1 exhibits the opposite effect on microRNA activity. Altogether, our results suggest that STAU-1 negatively modulates microRNA activity downstream of microRNA biogenesis, possibly by competing with microRNAs for binding on the 3' untranslated region of target mRNAs.

An optimized kit-free method for making strand-specific deep sequencing libraries from RNA fragments

Deep sequencing of strand-specific cDNA libraries is now a ubiquitous tool for identifying and quantifying RNAs in diverse sample types. The accuracy of conclusions drawn from these analyses depends on precise and quantitative conversion of the RNA sample into a DNA library suitable for sequencing. Here, we describe an optimized method of preparing strand-specific RNA deep sequencing libraries from small RNAs and variably sized RNA fragments obtained from ribonucleoprotein particle footprinting experiments or fragmentation of long RNAs. Our approach works across a wide range of input amounts (400 pg to 200 ng), is easy to follow and produces a library in 2–3 days at relatively low reagent cost, all while giving the user complete control over every step. Because all enzymatic reactions were optimized and driven to apparent completion, sequence diversity and species abundance in the input sample are well preserved.