Dataset Information

What human sperm RNA-Seq tells us about the microbiome.

ABSTRACT:

Purpose

The study was designed to assess the capacity of human sperm RNA-seq data to gauge the diversity of the associated microbiome within the ejaculate.

Methods

Semen samples were collected, and semen parameters evaluated at time of collection. Sperm RNA was isolated and subjected to RNA-seq. Microbial composition was determined by aligning sequencing reads not mapped to the human genome to the NCBI RefSeq bacterial, viral and archaeal genomes following RNA-Seq. Analysis of microbial assignments utilized phyloseq and vegan.

Results

Microbial composition within each sample was characterized as a function of microbial associated RNAs. Bacteria known to be associated with the male reproductive tract were present at similar levels in all samples representing 11 genera from four phyla with one exception, an outlier. Shannon diversity index (p?ConclusionThese results provide a first look at the microbiome as a component of human sperm RNA sequencing that has sufficient sensitivity to identify contamination or potential pathogenic bacterial colonization at least among the known contributors.

SUBMITTER: Swanson GM

PROVIDER: S-EPMC7056791 | biostudies-literature | 2020 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

What human sperm RNA-Seq tells us about the microbiome.

Swanson Grace M GM Moskovtsev Sergey S Librach Clifford C Pilsner J Richard JR Goodrich Robert R Krawetz Stephen A SA

Journal of assisted reproduction and genetics 20200104 2

<h4>Purpose</h4>The study was designed to assess the capacity of human sperm RNA-seq data to gauge the diversity of the associated microbiome within the ejaculate.<h4>Methods</h4>Semen samples were collected, and semen parameters evaluated at time of collection. Sperm RNA was isolated and subjected to RNA-seq. Microbial composition was determined by aligning sequencing reads not mapped to the human genome to the NCBI RefSeq bacterial, viral and archaeal genomes following RNA-Seq. Analysis of mic ...[more]

PMID: 31902104

Similar Datasets

Project description:A significant challenge to the effective application of RNA-seq to the complete transcript analysis of low quantity and/or degraded samples is the amplification of minimal input RNA to enable sequencing library construction. Several strategies have been commercialized in order to facilitate this goal. However, each strategy has its own specific protocols and methodology, and each may introduce unique bias and in some cases show specific preference for a collection of sequences. Our wider investigation of human spermatozoal RNAs was able to reveal their complexity despite being generally characterized by low quantity and high fragmentation. In this study, the following four commercially available RNA-seq amplification and library protocols for the preparation of low quantity/highly fragmented samples, SMARTer™ Ultra Low RNA (SU) for Illumina® Sequencing, SeqPlex RNA Amplification (SP), Ovation® RNA-Seq System V2 (OR), and Ovation® RNA-Seq Formalin Fixed Paraffin Embedded System (FFPES) were assessed using human sperm RNAs. Further investigation analyzed the effects on the end results of two different library preparation methods, Encore NGS Multiplex System I (Enc) and Ovation Ultralow Library Systems (UL), that appeared best suited to this type of RNA, along with other potential confounding factors such as FFPE preservation. Our results indicate that for each library preparation protocol, the differences in the initial amount of input RNA and choice of RNA purification step do not generate marked differences in terms of RNA profiling. However, substantial disparity is introduced by individual amplification methods prior to library construction. These significant differences may be caused by the different priming methods or amplification strategies used in each of the four different protocols examined. The observation of intra-sample variation introduced by the choice of protocol highlights the role that external factors play in planning and subsequent reliable interpretation of results of any RNA-seq experiment.

Project description:We report transcriptomes from 430 single glioblastoma cells isolated from 5 individual tumors and 102 single cells from gliomasphere cells lines generated using SMART-seq. In addition, we report population RNA-seq from the five tumors as well as RNA-seq from cell lines derived from 3 tumors (MGH26, MGH28, MGH31) cultured under serum free (CSC) and differentiated (FCS) conditions. This dataset highlights intratumoral heterogeneity with regards to the expression of de novo derived transcriptional modules and established subtype classifiers. Operative specimens from five glioblastoma patients (MGH26, MGH28, MGH29, MGH30, MGH31) were acutely dissociated, depleted for CD45+ inflammatory cells and then sorted as single cells (576 samples). Population controls for each tumor were isolated by sorting 2000-10000 cells and processed in parallel (5 population control samples). Single cells from two established cell lines, GBM6 and GBM8, were also sorted as single cells (192 samples). SMART-seq protocol was implemented to generate single cell full length transcriptomes (modified from Shalek, et al Nature 2013) and sequenced using 25 bp paired end reads. Single cell cDNA libraries for MGH30 were resequenced using 100 bp paired end reads to allow for isoform and splice junction reconstruction (96 samples, annotated MGH30L). Cells were also cultured in serum free conditions to generate gliomasphere cell lines for MGH26, MGH28, and MGH31 (CSC) which were then differentiated using 10% serum (FCS). Population RNA-seq was performed on these samples (3 CSC, 3 FCS, 6 total). The initial dataset included 875 RNA-seq libraries (576 single glioblastoma cells, 96 resequenced MGH30L, 192 single gliomasphere cells, 5 tumor population controls, 6 population libraries from CSC and FCS samples). Data was processed as described below using RSEM for quantification of gene expression. 5,948 genes with the highest composite expression either across all single cells combined (average log2(TPM)>4.5) or within a single tumor (average log2(TPM)>6 in at least one tumor) were included. Cells expressing less than 2,000 of these 5,948 genes were excluded. The final processed dataset then included 430 primary single cell glioblastoma transcriptomes, 102 single cell transcriptomes from cell lines(GBM6,GBM8), 5 population controls (1 for each tumor), and 6 population libraries from cell lines derived from the tumors (CSC and FCS for MGH26, MGH28 and MGH31). The final matrix (GBM_data_matrix.txt) therefore contains 5948 rows (genes) quantified in 543 samples (columns). Please note that the samples which are not included in the data processing are indicated in the sample description field.

Dataset Information

What human sperm RNA-Seq tells us about the microbiome.

Purpose

Methods

Results

Publications

What human sperm RNA-Seq tells us about the microbiome.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets