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.

Alteration of RNA modification signature in human sperm correlates with sperm motility.

RNA modifications, which are introduced post-transcriptionally, have recently been assigned pivotal roles in the regulation of spermatogenesis and embryonic development. However, the RNA modification landscape in human sperm is poorly characterized, hampering our understanding about the potential role played by RNA modification in sperm. Through our recently developed high-throughput RNA modification detection platform based on liquid chromatography with tandem mass spectroscopy, we are the first to have characterized the RNA modification signature in human sperm. The RNA modification signature was generated on the basis of 49 samples from participants, including 13 healthy controls, 21 patients with asthenozoospermia (AZS) and 15 patients with teratozoospermia (TZS). In total, we identified 13 types of RNA modification marks on the total RNA in sperm, and 16 types of RNA modification marks on sperm RNA fragments of different sizes. The levels of these RNA modifications on the RNA of patients with AZS or TZS were altered, compared to controls, especially on sperm RNA fragments > 80 nt. A few types of RNA modifications, such as m1G, m5C, m2G and m1A, showed clear co-expression patterns as well as high linear correlations with clinical sperm motility. In conclusion, we characterized the RNA modification signature of human sperm and identified its correlation with sperm motility, providing promising candidates for use in clinical sperm quality assessment and new research insights for exploring the underlying pathological mechanisms in human male infertility syndromes.

Epigenetic regulations in mammalian spermatogenesis: RNA-m6A modification and beyond.

Emerging evidence shows that m6A, one of the most abundant RNA modifications in mammals, is involved in the entire process of spermatogenesis, including mitosis, meiosis, and spermiogenesis. "Writers" catalyze m6A formation on stage-specific transcripts during male germline development, while "erasers" remove m6A modification to maintain a balance between methylation and demethylation. The different functions of RNA-m6A transcripts depend on their recognition by "readers". m6A modification mediates RNA metabolism, including mRNA splicing, translation, and degradation, as well as the maturity and biosynthesis of non-coding RNAs. Sperm RNA profiles are easily affected by environmental exposure and can even be inherited for several generations, similar to epigenetic inheritance. Here, we review and summarize the critical role of m6A in different developmental stages of male germ cells, to understand of the mechanisms and epigenetic regulation of m6A modifications. In addition, we also outline and discuss the important role of non-coding RNAs in spermatogenesis and RNA modifications in epigenetic inheritance.

Maternal overnutrition programs hedonic and metabolic phenotypes across generations through sperm tsRNAs.

There is a growing body of evidence linking maternal overnutrition to obesity and psychopathology that can be conserved across multiple generations. Recently, we demonstrated in a maternal high-fat diet (HFD; MHFD) mouse model that MHFD induced enhanced hedonic behaviors and obesogenic phenotypes that were conserved across three generations via the paternal lineage, which was independent of sperm methylome changes. Here, we show that sperm tRNA-derived small RNAs (tsRNAs) partly contribute to the transmission of such phenotypes. We observe increased expression of sperm tsRNAs in the F1 male offspring born to HFD-exposed dams. Microinjection of sperm tsRNAs from the F1-HFD male into normal zygotes reproduces obesogenic phenotypes and addictive-like behaviors, such as increased preference of palatable foods and enhanced sensitivity to drugs of abuse in the resultant offspring. The expression of several of the differentially expressed sperm tsRNAs predicted targets such as CHRNA2 and GRIN3A, which have been implicated in addiction pathology, are altered in the mesolimbic reward brain regions of the F1-HFD father and the resultant HFD-tsRNA offspring. Together, our findings demonstrate that sperm tsRNA is a potential vector that contributes to the transmission of MHFD-induced addictive-like behaviors and obesogenic phenotypes across generations, thereby emphasizing its role in diverse pathological outcomes.

Comparison of spermatozoal RNA extraction methods in goats.

RNA sequencing (RNAseq) has divulged newer role of spermatozoal RNA in male fertility. This study aimed to evaluate different sperm purification and RNA extraction methods for long-read RNA sequencing of poly(A) transcriptome in goat spermatozoa. Sperm samples were purified by swim-up separation using different purification medium. Spermatozoal RNA was extracted by seven different methods with additional supplementation of reducing agents in lysis buffer. poly(A) selected RNA was used for cDNA library preparation and long-read RNAseq in Nanopore sequencer. Sperm purification by 1 h swim-up resulted in higher recovery (89.20 ± 1.15%), motility (82.33 ± 1.53%), viability (88.10 ± 5.03%) and plasma membrane integrity (71.33 ± 4.51%) in sperm TALP (sp-TL) medium. A monophasic solution of GITC with phenol and DTT resulted in the highest yield of large sized RNAs (3.89 ± 0.46 ng/million cells) suitable for long-read RNAseq of poly(A) transcripts. RNAseq resulted in reads of length, ranging from 500bp to 2 Kb. A total of 123 transcripts were identified in goat spermatozoa by de novo assembly and included sperm-specific transcripts such as CATSPERG, PRM2, CYLC2, SPATA6, PLCZ1 etc. This study is the first report of long-read RNAseq of poly(A) transcriptome in goat spermatozoa.

Post-ejaculation thermal stress causes changes to the RNA profile of sperm in an external fertilizer.

Sperm cells experience considerable post-ejaculation environmental variation. However, little is known about whether this affects their molecular composition, probably owing to the assumption that sperm are transcriptionally quiescent. Nevertheless, recent evidence shows sperm have distinct RNA profiles that affect fertilization and embryo viability. Moreover, RNAs are expected to be highly sensitive to extracellular changes. One such group of RNAs are heat shock protein (hsp) transcripts, which function in stress responses and are enriched in sperm. Here, we exploit the experimental tractability of the mussel Mytilus galloprovincialis by exposing paired samples of ejaculated sperm to ambient (19°C) and increased (25°C) temperatures, then measure (i) sperm motility phenotypes, and (ii) messenger RNA (mRNA) levels of two target genes (hsp70 and hsp90) and several putative reference genes. We find no phenotypic changes in motility, but reduced mRNA levels for hsp90 and the putative reference gene gapdh at 25°C. This could reflect either decay of specific RNAs, or changes in translation and degradation rates of transcripts to maintain sperm function under stress. These findings represent, to our knowledge, the first evidence for changes in sperm RNA profiles owing to post-ejaculation environments, and suggest that sperm may be more vulnerable to stress from rising temperatures than currently thought.

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

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 < 0.001) and beta diversity (unweighted UniFrac distances, p = 9.99e-4; beta dispersion, p = 0.006) indicated the outlier was significantly different from all other samples. The outlier sample exhibited a dramatic increase in Streptococcus. Multiple testing indicated two operational taxonomic units, S. agalactiae and S. dysgalactiae (p = 0.009), were present. CONCLUSION: These 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.

A novel combination of silane-coated silica colloid with hybrid RNA extraction protocol and RNA enrichment for downstream applications of spermatozoal RNA.

Spermatozoa are specialised cells with low RNA content as compared to somatic cells. The suitable sperm RNA extraction and enrichment protocols for downstream applications are available for human, cattle, stallion and mouse but not for buffalo spermatozoa. Therefore, the present work was conducted to find out suitable colloidal solution for sperm purification and appropriate protocol for sperm RNA extraction and enrichment/amplification of RNA. For purification, we used PVP-coated silica colloidal solution (PVP-Si), silane-coated silica colloidal solution (Silane-Si) and iodixanol. Sperm recovery rate, total sperm motility and progressive sperm motility were significantly improved after separation by Silane-Si and iodixanol compared to PVA-Si method. The combined guanidinium thiocyanate-phenol-chloroform (GTPC) with silica matrix (SM)-based RNA extraction yielded more quantity of RNA in compared to individual method. The hybrid of SM and GTPC into a single protocol yielded 360-450 ng RNA from 30 million buffalo spermatozoa. For the first time, we adopted new way to enrich sperm RNA that increased the RNA concentration 4-5 times that was sufficient for downstream applications. The linear amplification of sperm RNA increased RNA concentration around 27-45 times. In summary, Silane-Si colloid for sperm separation, hybrid SM and GTPC protocol for sperm RNA extraction followed by enrichment or amplification of RNA was found suitable for high-throughput analyses of buffalo sperm RNA.

Sperm RNA preparation for transcriptomic analysis: Review of the techniques and personal experience.

In the last 10 years, several approaches, including microarrays, have been applied to investigate sperm transcript levels. However, success using microarray profiling is highly dependent of the quality of the RNA obtained. Therefore, the development of methods that deliver highly purified and intact RNA is of utmost importance. The three steps used to achieve this goal, purification of spermatozoa, RNA extraction and evaluation of RNA quality, are reviewed. Following that review and preliminary experiments, we processed sperm samples from seven normozoospermic men with a combination of gradient centrifugation and somatic cell lysis. RNA was extracted using the NucleoSpin RNA XS kit (Macherey-Nagel) and its purity checked using the BioAnalyzer. Hybridisation was done on Agilent SurePrint G3 Human GE 8 × 60K V2 microarrays. We identified 900 transcripts among the 1000 high abundance sperm transcripts reported in the literature. These genes are known to be involved in several biological processes, notably spermatogenesis, transcription regulation, cell growth and differentiation, sperm motility and capacitation, fertilisation, and embryogenesis. Therefore, our methodology is highly suitable for sperm transcriptomic analyses and can be used, notably, to compare mRNA profiles between fertile and infertile males.

Chicken sperm transcriptome profiling by microarray analysis.

It has been confirmed that mammalian sperm contain thousands of functional RNAs, and some of them have vital roles in fertilization and early embryonic development. Therefore, we attempted to characterize transcriptome of the sperm of fertile chickens using microarray analysis. Spermatozoal RNA was pooled from 10 fertile males and used for RNA preparation. Prior to performing the microarray, RNA quality was assessed using a bioanalyzer, and gDNA and somatic cell RNA contamination was assessed by CD4 and PTPRC gene amplification. The chicken sperm transcriptome was cross-examined by analysing sperm and testes RNA on a 4 × 44K chicken array, and results were verified by RT-PCR. Microarray analysis identified 21,639 predominantly nuclear-encoded transcripts in chicken sperm. The majority (66.55%) of the sperm transcripts were shared with the testes, while surprisingly, 33.45% transcripts were detected (raw signal intensity greater than 50) only in the sperm and not in the testes. The greatest proportion of up-regulated transcripts were responsible for signal transduction (63.20%) followed by embryonic development (56.76%) and cell structure (56.25%). Of the 20 most abundant transcripts, 18 remain uncharacterized, whereas the least abundant genes were mostly associated with the ribosome. These findings lay a foundation for more detailed investigations on sperm RNAs in chickens to identify sperm-based biomarkers for fertility.

piRNA expression patterns in high vs. low fertility bovine sperm.

PIWI-interacting RNAs (piRNAs) are 24-32 nucleotide RNA sequences primarily expressed in germ cells and developing embryos that suppress transposable element expression to protect genomic integrity during epigenetic reprogramming events. We characterized the expression of piRNA sequences and their encoding clusters in sperm samples from an idiopathic fertility model of Holstein bulls with high and low Sire Conception Rates. The piRNA populations were determined to be mostly similar between fertility conditions when investigated by principal component and differential expression analysis, suggesting that a high degree of conservation in the piRNA system is likely necessary for the production of viable sperm. Both fertility conditions demonstrated evidence of 'ping-pong' activity - a secondary biogenesis pathway associated with active transposable element targeting and suppression. Most sperm-borne piRNAs were between 29-30 nucleotides in length and originated from 226 clusters across the genome, with the exception of chromosome 20. Mapping analysis revealed abundant targeting of several transposable element families, suggesting a suppressive function of sperm piRNAs consistent with their established roles. Expression of genes targeted by sperm-borne piRNAs is significantly reduced throughout early embryogenesis compared to the mRNA population. Limited transposable element expression is known to be essential for spermatogenesis, thus epigenetic regulation of this pathway is likely to influence sperm quality and fertilizing capacity.

Phthalates impact on the epigenetic factors contributed specifically by the father at fertilization.

BACKGROUND: Preconception exposure to phthalates such as the anti-androgenic dibutyl-phthalate (DBP) impacts both male and female reproduction, yet how this occurs largely remains unknown. Previously we defined a series of RNAs expressly provided by sperm at fertilization and separately, and in parallel, those that responded to high DBP exposure. Utilizing both populations of RNAs, we now begin to unravel the impact of high-DBP exposure on those RNAs specifically delivered by the father. RESULTS: Enrichment of RNAs altered by DBP exposure within the Molecular Signature Database highlighted cellular stress, cell cycle, apoptosis, DNA damage response, and gene regulation pathways. Overlap within each of these five pathways identified those RNAs that were specifically (≥ fivefold enriched) or primarily (≥ twofold enriched) provided as part of the paternal contribution compared to the oocyte at fertilization. Key RNAs consistently altered by DBP, including CAMTA2 and PSME4, were delivered by sperm reflective of these pathways. The majority (64/103) of overlapping enriched gene sets were related to gene regulation. Many of these RNAs (45 RNAs) corresponded to key interconnected CRREWs (Chromatin remodeler cofactors, RNA interactors, Readers, Erasers, and Writers). Modeling suggests that CUL2, PHF10, and SMARCC1 may coordinate and mechanistically modulate the phthalate response. CONCLUSIONS: Mediated through a CRREW regulatory network, the cell responded to exposure presenting stressed-induced changes in the cell cycle-DNA damage-apoptosis. Interestingly, the majority of these DBP-responsive epigenetic mediators' direct acetylation or deacetylation, impacting the sperm's cargo delivered at fertilization and that of the embryo.

Small RNA shuffling between murine sperm and their cytoplasmic droplets during epididymal maturation.

Reports that mouse sperm gain small RNAs from the epididymosomes secreted by epididymal epithelial cells and that these "foreign" small RNAs act as an epigenetic information carrier mediating the transmission of acquired paternal traits have drawn great attention because the findings suggest that heritable information can flow from soma to germ line, thus invalidating the long-standing Weismann's barrier theory on heritable information flow. Using small RNA sequencing (sRNA-seq), northern blots, sRNA in situ hybridization, and immunofluorescence, we detected substantial changes in the small RNA profile in murine caput epididymal sperm (sperm in the head of the epididymis), and we further determined that the changes resulted from sperm exchanging small RNAs, mainly tsRNAs and rsRNAs, with cytoplasmic droplets rather than the epididymosomes. Moreover, the murine sperm-borne small RNAs were mainly derived from the nuclear small RNAs in late spermatids. Thus, caution is needed regarding sperm gaining foreign small RNAs as an underlying mechanism of epigenetic inheritance.

Heat shock protein family D member 1 in boar spermatozoa is strongly related to the litter size of inseminated sows.

BACKGROUND: Sperm quality evaluation is the logical first step in increasing field fertility. Spermatozoa contain cytoplasmic organelles and biomolecules known as sperm-intrinsic factors, which play key roles in sperm maturation, sperm-oocyte fusion, and embryo development. In particular, sperm membrane proteins [e.g., arginine vasopressin receptor 2, beta-actin, prohibitin, and heat shock protein family D member 1 (HSPD1)] and RNA could be used as functional indicators of male fertility. We sought to clarify the effects of differential mRNA expression of selected genes on several fertilisation parameters, including sperm motility, motion kinematics, capacitation, and litter size, in a porcine model. RESULTS: Our results demonstrated that HSPD1 expression was significantly correlated with male fertility, as measured by the litter size of inseminated sows. The expression of HSPD1 mRNA was linked to sperm motility and other motion kinematic characteristics. Furthermore, HSPD1 had a 66.7% overall accuracy in detecting male fertility, and the high-litter size group which was selected with the HSPD1 marker had a 1.34 greater litter size than the low-litter size group. CONCLUSIONS: Our findings indicate that HSPD1 might be a helpful biomarker for superior boar selection for artificial insemination, which could boost field fertility.

Establishment of a male fertility prediction model with sperm RNA markers in pigs as a translational animal model.

BACKGROUND: Male infertility is an important issue that causes low production in the animal industry. To solve the male fertility crisis in the animal industry, the prediction of sperm quality is the most important step. Sperm RNA is the potential marker for male fertility prediction. We hypothesized that the expression of functional genes related to fertilization will be the best target for male fertility prediction markers. To investigate optimum male fertility prediction marker, we compared target genes expression level and a wide range of field data acquired from artificial insemination of boar semen. RESULTS: Among the genes related to acrosomal vesicle exocytosis and sperm-oocyte fusion, equatorin (EQTN), zona pellucida sperm-binding protein 4 (ZP4), and sperm acrosome membrane-associated protein 3 exhibited high accuracy (70%, 90%, and 70%, respectively) as markers to evaluate male fertility. Combinations of EQTN-ZP4, ZP4-protein unc-13 homolog B, and ZP4-regulating synaptic membrane exocytosis protein 1 (RIMS1) showed the highest prediction value, and all these markers are involved in the acrosome reaction. CONCLUSION: The EQTN-ZP4 model was efficient in clustering the high-fertility group and may be useful for selection of animal that has superior fertility in the livestock industry. Compared to the EQTN-ZP4 model, the ZP4-RIMS1 model was more efficient in clustering the low-fertility group and may be useful in the diagnosis of male infertility in humans and other animals. The appointed translational animal model and established biomarker combination can be widely used in various scientific fields such as biomedical science.

Evolutionary history of sexual selection affects microRNA profiles in Drosophila sperm.

The presence of small RNAs in sperm is a relatively recent discovery and little is currently known about their importance and functions. Environmental changes including social conditions and dietary manipulations are known to affect the composition and expression of some small RNAs in sperm and may elicit a physiological stress response resulting in an associated change in gamete miRNA profiles. Here, we tested how microRNA profiles in sperm are affected by variation in both sexual selection and dietary regimes in Drosophila melanogaster selection lines. The selection lines were exposed to standard versus low yeast diet treatments and three different population sex ratios (male-biased, female-biased, or equal sex) in a full-factorial design. After 38 generations of selection, all males were maintained on their selected diet and in a common garden male-only environment prior to sperm sampling. We performed transcriptome analyses on miRNAs in purified sperm samples. We found 11 differentially expressed miRNAs with the majority showing differences between male- and female-biased lines. Dietary treatment only had a significant effect on miRNA expression levels in interaction with sex ratio. Our findings suggest that long-term adaptation may affect miRNA profiles in sperm and that these may show varied interactions with short-term environmental changes.

Characterization of the Impact of Density Gradient Centrifugation on the Profile of the Pig Sperm Transcriptome by RNA-Seq.

RNA-Seq data from human semen suggests that the study of the sperm transcriptome requires the previous elimination from the ejaculates of somatic cells carrying a larger load of RNA. Semen purification is also carried to study the sperm transcriptome in other species including swine and it is often done by density gradient centrifugation to obtain viable spermatozoa from fresh ejaculates or artificial insemination doses, thereby limiting the throughput and remoteness of the samples that can be processed in one study. The aim of this work was to evaluate the impact of purification with density gradient centrifugation by BoviPureTM on porcine sperm. Four boar ejaculates were purified with BoviPureTM and their transcriptome sequenced by RNA-Seq was compared with the RNA-Seq profiles of their paired non-purified sample. Seven thousand five hundred and nineteen protein coding genes were identified. Correlation, cluster, and principal component analysis indicated high-although not complete-similarity between the purified and the paired non-purified ejaculates. 372 genes displayed differentially abundant RNA levels between treatments. Most of these genes had lower abundances after purification and were mostly related to translation, transcription and metabolic processes. We detected a significant change in the proportion of genes of epididymal origin within the differentially abundant genes (1.3%) when compared with the catalog of unaltered genes (0.2%). In contrast, the proportion of testis-specific genes was higher in the group of unaltered genes (4%) when compared to the list of differentially abundant genes (0%). No proportion differences were identified for prostate, white blood, lymph node, tonsil, duodenum, skeletal muscle, liver, and mammary gland. Altogether, these results suggest that the purification impacts on the RNA levels of a small number of genes which are most likely caused by the removal of epididymal epithelial cells but also premature germinal cells, immature or abnormal spermatozoa or seminal exosomes with a distinct load of RNAs.

Epigenetic landscape of testis specific histone H2B variant and its influence on sperm function.

BACKGROUND: Biological relevance of the major testis specific histone H2B variant (TH2B) in sperm is not fully understood. Studies in TH2A/TH2B double knockout male mice indicate its role in chromatin compaction and male fertility. Additionally, the presence of TH2B and TH2A reportedly generates more dynamic nucleosomes, leading to an open chromatin structure characteristic of transcriptionally active genome. Given that mature human sperm are transcriptionally and translationally inactive, the presence of TH2B in mature sperm is intriguing. To address its role in sperm, we investigated the genome-wide localization of TH2B in sperm of fertile men. RESULTS: We have identified the genomic loci associated with TH2B in fertile human sperm by ChIP-seq analysis. Bioinformatic analysis revealed ~ 5% sperm genome and 5527 genes to be associated with TH2B. Out of these 105 (1.9%) and 144 (2.6%) genes showed direct involvement in sperm function and early embryogenesis, respectively. Chromosome wide analysis for TH2B distribution indicated its least distribution on X and Y chromosomes and varied distribution on autosomes. TH2B showed relatively higher percentage of gene association on chromosome 4, 18, 3 and 2. TH2B enrichment was more in promoter and gene body region. Gene Ontology (GO) analysis revealed signal transduction and associated kinase activity as the most enriched biological and molecular function, respectively. We also observed the enrichment of TH2B at developmentally important loci, such as HOXA and HOXD and on genes required for normal sperm function, few of which were validated by ChIP-qPCR. The relative expression of these genes was altered in particular subgroup of infertile men showing abnormal chromatin packaging. Chromatin compaction positively correlated with sperm- motility, concentration, viability and with transcript levels of PRKAG2 and CATSPER B. CONCLUSION: ChIP-seq analysis of TH2B revealed a putative role of TH2B in sperm function and embryo development. Altered expression of TH2B associated genes in infertile individuals with sperm chromatin compaction defects indicates involvement of TH2B in transcriptional regulation of these genes in post meiotic male germ cells. This altered transcriptome may be a consequence or cause of abnormal nuclear remodeling during spermiogenesis.

Transgenerational epigenetic reprogramming of early embryos: a mechanistic model.

The notion that epigenetic information can be transmitted across generations is supported by mounting waves of data, but the underlying mechanisms remain elusive. Here, a model is proposed which combines different lines of experimental evidence. First, it has been shown that somatic tissues exposed to stressing stimuli release circulating RNA-containing extracellular vesicles; second, epididymal spermatozoa can take up, internalize and deliver the RNA-containing extracellular vesicles to oocytes at fertilization; third, early embryos can process RNA-based information. These elements constitute the building blocks upon which the model is built. The model proposes that a continuous stream of epigenetic information flows from parental somatic tissues to the developing embryos. The flow can cross the Weismann barrier, is mediated by circulating vesicles and epididymal spermatozoa, and has the potential to generate epigenetic traits that are then stably acquired in the offspring. In a broader perspective, it emerges that a natural 'assembly line' operates continuously, aiming at passing the parental epigenetic blueprint in growing embryos.

Human chromatin remodeler cofactor, RNA interactor, eraser and writer sperm RNAs responding to obesity.

In the United States almost 33% of adults are considered obese (BMI > 30 kg/m2). Both animal models and to a lesser extent human studies, have associated BMI, a measure of obesity, with alterations in sperm DNA methylation and RNAs. Sperm RNAs from the Assessment of Multiple Gestations from Ovarian Stimulation trial, were isolated and sequenced. A Generalized Linear Model identified 487 BMI associated human sperm RNA elements (short exon-sized sequences). They partitioned into four patterns; a continual increase with BMI, increase once obese (BMI>30 kg/m2); a steady decrease with BMI; and decrease once overweight (BMI 25 - 30 kg/m2). Gene Ontology revealed a unique relationship between BMI and transcripts associated with chromosome organization, adipogenesis, cellular stress and obesity-related inflammation. Coregulatory networks linked by Chromatin remodeler cofactors, RNA interactors, Erasers and Writers (CREWs) were uncovered to reveal a hierarchical epigenetic response pathway.