Digital gene expression tag profiling analysis of the gene expression patterns regulating the early stage of mouse spermatogenesis

The Role of Placental DNA Methylation at Reproduction-Related Genes in Associations between Prenatal Bisphenol Analogues Exposure and the Digit Ratio in Children at Age 4: A Birth Cohort Study

Placental DNA methylation (DNAm) may be a potential mechanism underlying the effects of prenatal bisphenol analogues (BPs) exposure on reproductive health. Based on the Shanghai-Minhang Birth Cohort Study (S-MBCS), this study investigated associations of placental DNAm at reproduction-related genes with prenatal BPs exposure and children's digit ratios at age 4 using multiple linear regression models, and mediation analysis was further used to examine the mediating role of placental DNAm in the associations between prenatal BPs exposure and digit ratios among 345 mother-child pairs. Prenatal exposure to bisphenol A (BPA) was associated with hypermethylation at Protocadherin 8 (PCDH8), RBMX Like 2 (RBMXL2), and Sperm Acrosome Associated 1 (SPACA1), while bisphenol F (BPF) exposure was associated with higher methylation levels of Fibroblast Growth Factor 13 (FGF13). Consistent patterns were found in associations between higher DNAm at the 4 genes and increased digit ratios. Further mediation analysis showed that about 15% of the effect of BPF exposure on increased digit ratios was mediated by placental FGF13 methylation. In conclusion, the altered placental DNAm status might be a mediator underlying the feminizing effect of prenatal BPs exposure.

Identification of Key Genes Related to Postnatal Testicular Development Based on Transcriptomic Data of Testis in Hu Sheep

The selection of testis size can improve the reproductive capacity of livestock used for artificial insemination and has been considered as an important strategy for accelerating the breeding process. Although much work has been done to investigate the mechanisms of testis development in various species, there is little information available in regard to the differences in transcriptomic profiling of sheep testes at different developmental stages. In this work, we aimed to identify differentially expressed genes (DEGs) by RNA-Seq in sheep during different growth stages, including 0 month old (infant, M0), 3 months old (puberty, M3), 6 months old (sexual maturity, M6) and 12 months old (body maturity, M12). A total of 4,606 (2,381 up and 2,225 down), 7,500 (4,368 up and 3,132 down), 15 (8 up and seven down) DEGs were identified in M3_vs_M0, M6_vs_M3, and M12_vs_M6 comparison, respectively. Of which, a number of genes were continuously up-regulated and down-regulated with testicular development, including ODF3, ZPBP1, PKDREJ, MYBL1, PDGFA, IGF1, LH, INSL3, VIM, AMH, INHBA, COL1A1, COL1A2, and INHA. GO analysis illustrated that DEGs were mainly involved in testis development and spermatogenesis. KEGG analysis identified several important pathways and verified several reproduction-associated DEGs such as COL1A1, COL1A2, PDGFA, and IGF1. In addition, two gene modules highly associated with testis development and core genes with testis size were identified using weighted gene co-expression network analysis (WGCNA), including hub genes positively associated with testis size such as RANBP9, DNAH17, SPATA4, CIB4 and SPEM1, and those negatively associated with testis size such as CD81, CSK, PDGFA, VIM, and INHBA. This study comprehensively identified key genes related to sheep testicular development, which may provide potential insights for understanding male fertility and better guide in animal breeding.

Transcriptomics of Meiosis in the Male Mouse

Molecular studies of meiosis in mammals have been long relegated due to some intrinsic obstacles, namely the impossibility to reproduce the process in vitro, and the difficulty to obtain highly pure isolated cells of the different meiotic stages. In the recent years, some technical advances, from the improvement of flow cytometry sorting protocols to single-cell RNAseq, are enabling to profile the transcriptome and its fluctuations along the meiotic process. In this mini-review we will outline the diverse methodological approaches that have been employed, and some of the main findings that have started to arise from these studies. As for practical reasons most studies have been carried out in males, and mostly using mouse as a model, our focus will be on murine male meiosis, although also including specific comments about humans. Particularly, we will center on the controversy about gene expression during early meiotic prophase; the widespread existing gap between transcription and translation in meiotic cells; the expression patterns and potential roles of meiotic long non-coding RNAs; and the visualization of meiotic sex chromosome inactivation from the RNAseq perspective.

Expression and role of HEPIS in breast cancer

Human embryo lung cellular protein interacting with severe acute respiratory syndrome-coronavirus nonstructural protein-10 (HEPIS) is expressed at varying levels in multiple organs and breast cancer cell lines. However, its expression and function in breast cancer cells has yet to be studied. Therefore, RNA in situ hybridization was used to detect the expression of HEPIS in breast cancer and cancer-adjacent normal breast tissue. HEPIS was expressed at lower levels in breast cancer compared with that in adjacent normal tissue. Subcellular localization of HEPIS was mainly found in the cytoplasm of HeLa cells. Cell Counting Kit-8 and 5-ethynyl-2′-deoxyuridine cell proliferation assays were used to investigate the role of HEPIS in cancer cell proliferation. Ectopic expression of HEPIS in MCF-7 cells was found to significantly inhibit cell proliferation. In contrast, knockdown of HEPIS by RNA interference exhibited the opposite effect. Furthermore, a dual-luciferase reporter assay was performed and HEPIS overexpression specifically inhibited the activity of the NF-κB reporter gene. Results of the gene chip assay revealed that 2,231 genes were differentially expressed in HEPIS-overexpressing cells. Results of the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that these genes were enriched in the ‘mitogen-activated protein kinase signaling pathway’, ‘JAK-STAT signaling pathway’ and ‘focal adhesion’. Reverse transcription-quantitative PCR was used to confirm the expression levels of the differentially expressed genes interleukin 2 receptor subunit α (IL2RA), interferon α and β receptor subunit 2 (IFNAR2) and IFα8 (IFNA8). In conclusion, the results of the present study indicated that HEPIS may function as a potential repressor of breast cancer.

SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data

Quality control (QC) and preprocessing are essential steps for sequencing data analysis to ensure the accuracy of results. However, existing tools cannot provide a satisfying solution with integrated comprehensive functions, proper architectures, and highly scalable acceleration. In this article, we demonstrate SOAPnuke as a tool with abundant functions for a “QC-Preprocess-QC” workflow and MapReduce acceleration framework. Four modules with different preprocessing functions are designed for processing datasets from genomic, small RNA, Digital Gene Expression, and metagenomic experiments, respectively. As a workflow-like tool, SOAPnuke centralizes processing functions into 1 executable and predefines their order to avoid the necessity of reformatting different files when switching tools. Furthermore, the MapReduce framework enables large scalability to distribute all the processing works to an entire compute cluster.

We conducted a benchmarking where SOAPnuke and other tools are used to preprocess a ∼30× NA12878 dataset published by GIAB. The standalone operation of SOAPnuke struck a balance between resource occupancy and performance. When accelerated on 16 working nodes with MapReduce, SOAPnuke achieved ∼5.7 times the fastest speed of other tools.

Identification and comparative analysis of the ovarian microRNAs of prolific and non-prolific goats during the follicular phase using high-throughput sequencing

The kidding rate is one of the most important economic traits for goat production, but the genetic mechanism that is associated with ovulation rate is poorly understood. Recently, increasing evidence has suggested that microRNAs (miRNAs) influence ovarian biological processes. The present study provides the first comparison of the ovarian miRNAs of prolific Jintang black goats (JTGs) and non-prolific Tibetan goats (TBGs) during the follicular phase using RNA-Seq technology. We generated 11.19 million (M) and 11.34 M clean reads from the TBG and JTG libraries, respectively, from which a total of 389 known miRNAs were identified and 142 novel miRNAs were predicted. A total of 191 miRNAs were differentially expressed between the two breeds. Among the 10 most abundant miRNAs, miR-21-5p was defined as differentially expressed miRNA with a higher level in the JTG library than in the TBG library, but the other miRNAs were not different between the breeds. The predicted miRNA-targeted genes were further analyzed by Gene Ontology and KEGG pathway analyses. The results revealed that miR-21, miR-99a, miRNA-143, let-7f, miR-493 and miR-200b may affect follicular development. These findings will increase the current understanding of the role of ovarian miRNAs in the regulation of ovulation rate in goats.

Molecular and Cellular Mechanisms of Apoptosis during Dissociated Spermatogenesis

Apoptosis is a tightly controlled process by which tissues eliminate unwanted cells. Spontaneous germ cell apoptosis in testis has been broadly investigated in mammals that have an associated spermatogenesis pattern. However, the mechanism of germ cell apoptosis in seasonally breeding reptiles following a dissociated spermatogenesis has remained enigmatic. In the present study, morphological evidence has clearly confirmed the dissociated spermatogenesis pattern in Pelodiscus sinensis. TUNEL and TEM analyses presented dynamic changes and ultrastructural characteristics of apoptotic germ cells during seasonal spermatogenesis, implying that apoptosis might be one of the key mechanisms to clear degraded germ cells. Furthermore, using RNA-Seq and digital gene expression (DGE) profiling, a large number of apoptosis-related differentially expressed genes (DEGs) at different phases of spermatogenesis were identified and characterized in the testis. DGE and RT-qPCR analysis revealed that the critical anti-apoptosis genes, such as Bcl-2, BAG1, and BAG5, showed up-regulated patterns during intermediate and late spermatogenesis. Moreover, the increases in mitochondrial transmembrane potential in July and October were detected by JC-1 staining. Notably, the low protein levels of pro-apoptotic cleaved caspase-3 and CytC in cytoplasm were detected by immunohistochemistry and western blot analyses, indicating that the CytC-Caspase model might be responsible for the effects of germ cell apoptosis on seasonal spermatogenesis. These results facilitate understanding the regulatory mechanisms of apoptosis during spermatogenesis and uncovering the biological process of the dissociated spermatogenesis system in reptiles.

Complete depletion of primordial germ cells in an All-female fish leads to Sex-biased gene expression alteration and sterile All-male occurrence

Background

Gynogenesis is one of unisexual reproduction modes in vertebrates, and produces all-female individuals with identical genetic background. In sexual reproduction vertebrates, the roles of primordial germ cells on sexual dimorphism and gonadal differentiation have been largely studied, and two distinct functional models have been proposed. However, the role of primordial germ cells remains unknown in unisexual animals, and it is also unclear whether the functional models in sexual reproduction animals are common in unisexual animals.

Results

To solve these puzzles, we attempt to utilize the gynogenetic superiority of polyploid Carassius gibelio to create a complete germ cell-depleted gonad model by a similar morpholino-mediated knockdown approach used in other examined sexual reproduction fishes. Through the germ cell-depleted gonad model, we have performed comprehensive and comparative transcriptome analysis, and revealed a complete alteration of sex-biased gene expression. Moreover, the expression alteration leads to up-regulation of testis-biased genes and down-regulation of ovary-biased genes, and results in the occurrence of sterile all-males with testis-like gonads and secondary sex characteristics in the germ cell-depleted gynogenetic Carassius gibelio.

Conclusions

Our current results have demonstrated that unisexual gynogenetic embryos remain keeping male sex determination information in the genome, and the complete depletion of primordial germ cells in the all-female fish leads to sex-biased gene expression alteration and sterile all-male occurrence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2130-z) contains supplementary material, which is available to authorized users.

The Analysis of Gene Expression on Fertility Decline in Caenorhabditis elegans after the Treatment with 5-Fluorouracil

BACKGROUND

5-Fluorouracil could lead to a decline in fertility in Caenorhabditis elegans. The aim of this study was to describe the mechanisms underlying such an altered fertility phenotype and to illustrate the specific genes and pathways that are involved in the related phenotypic changes in C. elegans.

METHODS

We isolated total RNA from the samples and used a new method called Digital Gene Expression (DGE), which can rapidly identify genes with altered transcript levels. The random genes were confirmed by real-time RT-PCR.

RESULTS

We analyzed the results of two methods to draw conclusions based on a comparison between C. elegans and other harmful parasites. Compared with controls, 1147 genes were up-regulated, and 1067 were down-regulated. Overall, 101 up-regulated genes had a log2 ratio higher than 8, whereas the log2 ratio of 141 down-regulated genes was higher than 8. After mapping to the reference database, 4 pathways were confirmed to be involved in this phenomenon, with statistically significant participation from 19 genes.

CONCLUSION

For the first time, the transcript sequence of 5-Fu-treated worms and controls was detected. We found that 4 possible pathways, i.e., ECM-receptor interaction pathway, TGF-beta signaling pathway, Focal adhesion and Hypertrophic cardiomyopathy, may be involved in the number decline in the embryos of C. elegans. Specifically, the ECM-receptor interaction pathway and Focal adhesion may be very important pathways that alter the reproduction of C. elegans.

Digital gene expression approach over multiple RNA-Seq data sets to detect neoblast transcriptional changes in Schmidtea mediterranea

BACKGROUND

The freshwater planarian Schmidtea mediterranea is recognised as a valuable model for research into adult stem cells and regeneration. With the advent of the high-throughput sequencing technologies, it has become feasible to undertake detailed transcriptional analysis of its unique stem cell population, the neoblasts. Nonetheless, a reliable reference for this type of studies is still lacking.

RESULTS

Taking advantage of digital gene expression (DGE) sequencing technology we compare all the available transcriptomes for S. mediterranea and improve their annotation. These results are accessible via web for the community of researchers. Using the quantitative nature of DGE, we describe the transcriptional profile of neoblasts and present 42 new neoblast genes, including several cancer-related genes and transcription factors. Furthermore, we describe in detail the Smed-meis-like gene and the three Nuclear Factor Y subunits Smed-nf-YA, Smed-nf-YB-2 and Smed-nf-YC.

CONCLUSIONS

DGE is a valuable tool for gene discovery, quantification and annotation. The application of DGE in S. mediterranea confirms the planarian stem cells or neoblasts as a complex population of pluripotent and multipotent cells regulated by a mixture of transcription factors and cancer-related genes.

Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development

The present study profiled and analyzed gene expression of the maize ear at four key developmental stages. Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Some of the differentially expressed genes (DEGs) were predicted to be potential candidates of maize ear development. Several well-known genes were found with reported mutant analyses, such as, compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). MicroRNAs such as microRNA156 were predicted to target genes involved in maize ear development. Antisense transcripts were widespread throughout all the four stages, and are suspected to play important roles in maize ear development. Thus, identification and characterization of important genes and regulators at all the four developmental stages will contribute to an improved understanding of the molecular mechanisms responsible for maize ear development.

The Wilms tumor gene, Wt1, is critical for mouse spermatogenesis via regulation of sertoli cell polarity and is associated with non-obstructive azoospermia in humans

Azoospermia is one of the major reproductive disorders which cause male infertility in humans; however, the etiology of this disease is largely unknown. In the present study, six missense mutations of WT1 gene were detected in 529 human patients with non-obstructive azoospermia (NOA), indicating a strong association between WT1 mutation and NOA. The Wilms tumor gene, Wt1, is specifically expressed in Sertoli cells (SCs) which support spermatogenesis. To examine the functions of this gene in spermatogenesis, Wt1 was deleted in adult testis using Wt1^flox and Cre-ER^TM mice strains. We found that inactivation of Wt1 resulted in massive germ cell death and only SCs were present in most of the seminiferous tubules which was very similar to NOA in humans. In investigating the potential mechanism for this, histological studies revealed that the blood–testis barrier (BTB) was disrupted in Wt1 deficient testes. In vitro studies demonstrated that Wt1 was essential for cell polarity maintenance in SCs. Further studies found that the expression of cell polarity associated genes (Par6b and E-cadherin) and Wnt signaling genes (Wnt4, Wnt11) were downregulated in Wt1 deficient SCs, and that the expression of Par6b and E-cadherin was regulated by Wnt4. Our findings suggest that Wt1 is important in spermatogenesis by regulating the polarity of SCs via Wnt signaling pathway and that WT1 mutation is one of the genetic causes of NOA in humans.

Infertility is one of the most common health problems, affecting about 15% of the couples in the world. In about half of these couples, infertility is related to male reproductive defect. Azoospermia is one of the major causes of male infertility in humans. Previous studies have found that the mutation or deletion of some genes is associated with azoospermia; however, the genetic cause of this remains largely unknown. In the present study, we detected Wt1 missense mutations in men with non-obstructive azoospermia (NOA). An essential function for WT1 in male spermatogenesis was confirmed by the use of a Wt1 conditional knockout mouse strain. Inactivation of Wt1 resulted in germ cell loss in mice, which was similar to NOA in human patients. Our data indicate that WT1 mutation is one genetic cause of male infertility and suggest that WT1 mutational analysis will be useful for diagnosis in a clinical setting.