RNaseIII and T4 polynucleotide Kinase sequence biases and solutions during RNA-seq library construction

Transcriptomic analysis reveals that non-forage or forage fiber source promotes rumen development through different metabolic processes in lambs

Dietary fiber supplementation can stimulate rumen development in lambs during the pre-weaning period. However, it is unclear whether different sources of fiber have varying effects on rumen development. This study aimed to investigate the molecular mechanism of rumen morphological and functional development based on non-forage or forage as a starter dietary fiber source. Twenty-four male Hu lambs with similar body weights (BW, 3.67 ± 0.08 kg) were selected and divided into two groups that received diets supplemented with either alfalfa hay (AH) or soybean hull (SH). At the age of 70 days, six lambs were slaughtered from each treatment group for rumen fermentation and morphological analyses. Three samples of the rumen tissue from the ventral sac were collected for transcriptomic analysis. The results identified 633 differentially expressed genes (DEGs), of which 210 were upregulated and 423 were downregulated in the SH group compared with those in the AH group. The upregulated DEGs were most enriched in the immune function and proteolysis pathways, whereas the downregulated DEGs were mainly involved in cell proliferation, apoptosis, and differentiation pathways. These findings indicated that non-forage as a starter dietary fiber source improved immune function and enhanced nitrogen utilization, whereas forage facilitated rumen morphological development.

Mining Candidate Genes Related to Heavy Metals in Mature Melon (Cucumis melo L.) Peel and Pulp Using WGCNA

The content of metal ions in fruits is inseparable from plant intake of trace elements and health effects in the human body. To understand metal ion content in the fruit and pericarp of melon (Cucumis melo L.) and the candidate genes responsible for controlling this process, we analyzed the metal ion content in distinct parts of melon fruit and pericarp and performed RNA-seq. The results showed that the content of metal ions in melon fruit tissue was significantly higher than that in the pericarp. Based on transcriptome expression profiling, we found that the fruit and pericarp contained elevated levels of DEGs. GO functional annotations included cell surface receptor signaling, signal transduction, organic substance metabolism, carbohydrate derivative binding, and hormone-mediated signaling pathways. KEGG pathways included pectate lyase, pentose and glucuronate interconversions, H⁺-transporting ATPase, oxidative phosphorylation, plant hormone signal transduction, and MAPK signaling pathways. We also analyzed the expression patterns of genes and transcription factors involved in hormone biosynthesis and signal transduction. Using weighted gene co-expression network analysis (WGCNA), a co-expression network was constructed to identify a specific module that was significantly correlated with the content of metal ions in melon, after which the gene expression in the module was measured. Connectivity and qRT–PCR identified five candidate melon genes, LOC103501427, LOC103501539, LOC103503694, LOC103504124, and LOC107990281, associated with metal ion content. This study provides a theoretical basis for further understanding the molecular mechanism of heavy metal ion content in melon fruit and peel and provides new genetic resources for the study of heavy metal ion content in plant tissues.

Early feeding strategies in lambs affect rumen development and growth performance, with advantages persisting for two weeks after the transition to fattening diets

This study aimed to explore the effects of early feeding strategies on the growth and rumen development of lambs from pre-weaning to the transition to fattening diets. Ninety-six newborn, male lambs with similar body weights were randomly assigned to three treatments: fed starter at 42 days old + weaned at 56 days old (Ctrl, n = 36), fed starter at 7 days old + weaned at 56 days old (ES, n = 36), and fed starter at 7 days old + weaned at 28 days old (ES + EW, n = 24). The fattening diets of all lambs were gradually replaced from 60 to 70 days of age. Six randomly selected lambs from each treatment were slaughtered at 14, 28, 42, 56, 70, and 84 days of age. The results showed that the richness and diversity of rumen microbiota of lambs in the Ctrl group were distinct from those of lambs in the other groups at 42 days of age. Moreover, transcriptome analysis revealed 407, 219, and 1,211 unique differentially expressed genes (DEGs) in the rumen tissue of ES vs. Ctrl, ES vs. ES + EW, and ES + EW vs. Ctrl groups, respectively, at 42 days of age. Different early feeding strategies resulted in differences in ruminal anatomy, morphology, and fermentation in lambs from 42 to 84 days of age (P < 0.05). Lambs in the ES + EW group had a higher average starter diet intake than those in the other groups (P < 0.05) from 28 to 56 days of age, which affected their growth performance. After 42 days of age, the body and carcass weights of lambs in the ES and ES + EW groups were higher than those in the Ctrl group (P < 0.05). These findings demonstrate that feeding lambs with a starter diet at 7 days of age and weaning them at 28 days of age can promote rumen development and improve growth performance, and this advantage persists for up to 2 weeks after transition to the fattening diet.

Phospho-RNA sequencing with circAID-p-seq

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

Phospho-RNA-seq: a modified small RNA-seq method that reveals circulating mRNA and lncRNA fragments as potential biomarkers in human plasma

Extracellular RNAs (exRNAs) in biofluids have attracted great interest as potential biomarkers. Although extracellular microRNAs in blood plasma are extensively characterized, extracellular messenger RNA (mRNA) and long non-coding RNA (lncRNA) studies are limited. We report that plasma contains fragmented mRNAs and lncRNAs that are missed by standard small RNA-seq protocols due to lack of 5' phosphate or presence of 3' phosphate. These fragments were revealed using a modified protocol ("phospho-RNA-seq") incorporating RNA treatment with T4-polynucleotide kinase, which we compared with standard small RNA-seq for sequencing synthetic RNAs with varied 5' and 3' ends, as well as human plasma exRNA Analyzing phospho-RNA-seq data using a custom, high-stringency bioinformatic pipeline, we identified mRNA/lncRNA transcriptome fingerprints in plasma, including tissue-specific gene sets. In a longitudinal study of hematopoietic stem cell transplant patients, bone marrow- and liver-enriched exRNA genes were tracked with bone marrow recovery and liver injury, respectively, providing proof-of-concept validation as a biomarker approach. By enabling access to an unexplored realm of mRNA and lncRNA fragments, phospho-RNA-seq opens up new possibilities for plasma transcriptomic biomarker development.

Small RNA-seq: The RNA 5'-end adapter ligation problem and how to circumvent it

The preparation of small RNA cDNA sequencing libraries depends on the unbiased ligation of adapters to the RNA ends. Small RNA with 5' recessed ends are poor substrates for enzymatic adapter ligation, but this 5' adapter ligation problem can go undetected if the library preparation steps are not monitored. Here we illustrate the severity of the 5' RNA end ligation problem using several pre-miRNA-like hairpins that allow us to expand the definition of the problem to include 5' ends close to a hairpin stem, whether recessed or in a short extension. The ribosome profiling method can avoid a difficult 5' adapter ligation, but the enzyme typically used to circularize the cDNA has been reported to be biased, calling into question the benefit of this workaround. Using the TS2126 RNA ligase 1 (a.k.a. CircLigase) as the circularizing enzyme, we devised a bias test for the circularization of first strand cDNA. All possible dinucleotides were circle-ligated with similar efficiency. To re-linearize the first strand cDNA in the ribosome profiling approach, we introduce an improved method wherein a single ribonucleotide is placed between the sequencing primer binding sites in the reverse transcriptase primer, which later serves as the point of re-linearization by RNase A. We incorporate this step into the ribosomal profiling method and describe a complete improved library preparation method, Coligo-seq, for the sequencing of small RNA with secondary structure close to the 5' end. This method accepts a variety of 5' modified RNA, including 5' monophosphorylated RNA, as demonstrated by the construction of a HeLa cell microRNA cDNA library.

Effects of dietary physical or nutritional factors on morphology of rumen papillae and transcriptome changes in lactating dairy cows based on three different forage-based diets

BACKGROUND

Rumen epithelial tissue plays an important role in nutrient absorption and rumen health. However, whether forage quality and particle size impact the rumen epithelial morphology is unclear. The current study was conducted to elucidate the effects of forage quality and forage particle size on rumen epithelial morphology and to identify potential underlying molecular mechanisms by analyzing the transcriptome of the rumen epithelium (RE). To achieve these objectives, 18 mid-lactation dairy cows were allocated to three groups (6 cows per group), and were fed with one of three different forage-based diets, alfalfa hay (AH), corn stover (CS), and rice straw (RS) for 14 weeks, respectively. Ruminal volatile fatty acids (VFAs) and epithelial thickness were determined, and RNA-sequencing was conducted to identify the transcriptomic changes of rumen epithelial under different forage-based diets.

RESULTS

The RS diet exhibited greater particle size but low quality, the AH diet was high nutritional value but small particle size, and CS diet was low quality and small particle size. The ruminal total VFA concentration was greater in AH compared with those in CS or RS. The width of the rumen papillae was greater in RS-fed cows than in cows fed AH or CS. In total, 31, 40, and 28 differentially expressed (DE, fold change > 2, FDR < 0.05) genes were identified via pair-wise comparisons including AH vs. CS, AH vs. RS, and RS vs. CS, respectively. Functional classification analysis of DE genes revealed dynamic changes in ion binding (such as DSG1) between AH and CS, proliferation and apoptotic processes (such as BAG3, HLA-DQA1, and UGT2B17) and complement activation (such as C7) between AH or RS and CS. The expression of HLA-DQA1 was down-regulated in RS compared with AH and CS, and the expression of UGT2B17 was down-regulated in RS compared with CS, with positive (R = 0.94) and negative (R = -0.96) correlation with the width of rumen epithelial papillae (P < 0.05), respectively.

CONCLUSION

Our results suggest that both nutrients (VFAs) and particle sizes can alter expression of genes involved in cell proliferation/apoptosis process and complement complex. Our results suggest that particle size may be more important in regulating rumen epithelial morphology when animals are fed with low-quality forage diets and the identified DE genes may affect the RE nutrient absorption or morphology of RE. Our findings provide insights into the effects of the dietary particle size in the future management of dairy cow feeding, that when cows were fed with low-quality forage (such as rice straw), smaller particle size may be beneficial for nutrients absorption and milk production.

Digital gene expression analysis with sample multiplexing and PCR duplicate detection: A straightforward protocol

Tag-Seq is a high-throughput approach used for discovering SNPs and characterizing gene expression. In comparison to RNA-Seq, Tag-Seq eases data processing and allows detection of rare mRNA species using only one tag per transcript molecule. However, reduced library complexity raises the issue of PCR duplicates, which distort gene expression levels. Here we present a novel Tag-Seq protocol that uses the least biased methods for RNA library preparation combined with a novel approach for joint PCR template and sample labeling. In our protocol, input RNA is fragmented by hydrolysis, and poly(A)-bearing RNAs are selected and directly ligated to mixed DNA-RNA P5 adapters. The P5 adapters contain i5 barcodes composed of sample-specific (moderately) degenerate base regions (mDBRs), which later allow detection of PCR duplicates. The P7 adapter is attached via reverse transcription with individual i7 barcodes added during the amplification step. The resulting libraries can be sequenced on an Illumina sequencer. After sample demultiplexing and PCR duplicate removal with a free software tool we designed, the data are ready for downstream analysis. Our protocol was tested on RNA samples from predator-induced and control Daphnia microcrustaceans.

An optimized protocol for generation and analysis of Ion Proton sequencing reads for RNA-Seq

BACKGROUND

Previous studies compared running cost, time and other performance measures of popular sequencing platforms. However, comprehensive assessment of library construction and analysis protocols for Proton sequencing platform remains unexplored. Unlike Illumina sequencing platforms, Proton reads are heterogeneous in length and quality. When sequencing data from different platforms are combined, this can result in reads with various read length. Whether the performance of the commonly used software for handling such kind of data is satisfactory is unknown.

RESULTS

By using universal human reference RNA as the initial material, RNaseIII and chemical fragmentation methods in library construction showed similar result in gene and junction discovery number and expression level estimated accuracy. In contrast, sequencing quality, read length and the choice of software affected mapping rate to a much larger extent. Unspliced aligner TMAP attained the highest mapping rate (97.27 % to genome, 86.46 % to transcriptome), though 47.83 % of mapped reads were clipped. Long reads could paradoxically reduce mapping in junctions. With reference annotation guide, the mapping rate of TopHat2 significantly increased from 75.79 to 92.09 %, especially for long (>150 bp) reads. Sailfish, a k-mer based gene expression quantifier attained highly consistent results with that of TaqMan array and highest sensitivity.

CONCLUSION

We provided for the first time, the reference statistics of library preparation methods, gene detection and quantification and junction discovery for RNA-Seq by the Ion Proton platform. Chemical fragmentation performed equally well with the enzyme-based one. The optimal Ion Proton sequencing options and analysis software have been evaluated.

RNA-seq of human reference RNA samples using a thermostable group II intron reverse transcriptase

Next-generation RNA sequencing (RNA-seq) has revolutionized our ability to analyze transcriptomes. Current RNA-seq methods are highly reproducible, but each has biases resulting from different modes of RNA sample preparation, reverse transcription, and adapter addition, leading to variability between methods. Moreover, the transcriptome cannot be profiled comprehensively because highly structured RNAs, such as tRNAs and snoRNAs, are refractory to conventional RNA-seq methods. Recently, we developed a new method for strand-specific RNA-seq using thermostable group II intron reverse transcriptases (TGIRTs). TGIRT enzymes have higher processivity and fidelity than conventional retroviral reverse transcriptases plus a novel template-switching activity that enables RNA-seq adapter addition during cDNA synthesis without using RNA ligase. Here, we obtained TGIRT-seq data sets for well-characterized human RNA reference samples and compared them to previous data sets obtained for these RNAs by the Illumina TruSeq v2 and v3 methods. We find that TGIRT-seq recapitulates the relative abundance of human transcripts and RNA spike-ins in ribo-depleted, fragmented RNA samples comparably to non-strand-specific TruSeq v2 and better than strand-specific TruSeq v3. Moreover, TGIRT-seq is more strand specific than TruSeq v3 and eliminates sampling biases from random hexamer priming, which are inherent to TruSeq. The TGIRT-seq data sets also show more uniform 5' to 3' gene coverage and identify more splice junctions, particularly near the 5' ends of mRNAs, than do the TruSeq data sets. Finally, TGIRT-seq enables the simultaneous profiling of mRNAs and lncRNAs in the same RNA-seq experiment as structured small ncRNAs, including tRNAs, which are essentially absent with TruSeq.

Meeting report: discussions and preliminary findings on extracellular RNA measurement methods from laboratories in the NIH Extracellular RNA Communication Consortium

Extracellular RNAs (exRNAs) have been identified in all tested biofluids and have been associated with a variety of extracellular vesicles, ribonucleoprotein complexes and lipoprotein complexes. Much of the interest in exRNAs lies in the fact that they may serve as signalling molecules between cells, their potential to serve as biomarkers for prediction and diagnosis of disease and the possibility that exRNAs or the extracellular particles that carry them might be used for therapeutic purposes. Among the most significant bottlenecks to progress in this field is the lack of robust and standardized methods for collection and processing of biofluids, separation of different types of exRNA-containing particles and isolation and analysis of exRNAs. The Sample and Assay Standards Working Group of the Extracellular RNA Communication Consortium is a group of laboratories funded by the U.S. National Institutes of Health to develop such methods. In our first joint endeavour, we held a series of conference calls and in-person meetings to survey the methods used among our members, placed them in the context of the current literature and used our findings to identify areas in which the identification of robust methodologies would promote rapid advancements in the exRNA field.

Altered gamma-aminobutyric acid type B receptor subunit 1 splicing in alcoholics

BACKGROUND

Chronic alcohol exposure can change splice variant expression. The gamma-aminobutyric acid type B (GABAB) receptor undergoes splicing and is an alcoholism treatment target, but there is little information about splicing changes in this receptor in alcoholics. We studied GABAB receptor subunit 1 (GABAB1) splicing in alcoholic postmortem brains.

METHODS

To maximize GABAB1 splice junction identification, we combined gene specific libraries with RNA-seq. Splice junctions and mapped reads were also found from intronic and intergenic regions. We compared GABAB1 splice junctions in prefrontal cortices from 14 alcoholic and 15 control subjects and introduced new strategies, reads per kilobase of splice junction model per million mapped reads and reads per kilobase of gene model per million mapped reads, for quantitating splice junction and gene expression.

RESULTS

Novel splice junction detection indicated that the GABAB1 gene is at least two times longer than the previously reported gene length. GABAB1 exon and intron expression data showed low expression at the 5' end exons and exon grouping. This indicated that there are short splicing variants in addition to GABAB receptor subunit GABAB1a, the longest known major transcript. We found that chronic alcohol altered exon/intron expression and splice junction levels. Decreased expression of the gamma-aminobutyric acid binding site, a transmembrane domain and a microRNA binding site may decrease normal GABAB1 transcript population and thereby decrease normal signal transduction in alcoholics.

CONCLUSIONS

We discovered novel, complex splicing of GABAB1 in human brain and showed that chronic alcohol produces additional splicing complexity.