Analysis of mouse germ-cell transcriptome at different stages of spermatogenesis by SAGE: Biological significance

Embryonic stem cell derived germ cells induce spermatogenesis after transplantation into the testes of an adult mouse azoospermia model

The present study aimed to: (i) identify the exogenous factors that allow in vitro differentiation of mouse spermatogonial stem cells (SSCs) from embryonic stem cells (ESCs); (ii) evaluate the effects of Sertoli cells in SSC enrichment; and (iii) assess the success of transplantation using in vitro differentiated SSCs in a mouse busulfan-treated azoospermia model. A 1-day-old embryoid body (EB) received 5 ng/ml of bone morphogenetic protein 4 (BMP4) for 4 days, 3 µM retinoic acid (RA) in a SIM mouse embryo-derived thioguanine and ouabain resistant (STO) co-culture system for 7 days, and was subsequently co-cultured for 2 days with Sertoli cells in the presence or absence of a leukaemia inhibitory factor (LIF), basic fibroblast growth factor (bFGF) and RA composition, and in the presence of these factors in simple culture medium. Higher viability, proliferation and germ cell gene expression were seen in the presence of the LIF, bFGF and RA composition, on top of Sertoli cells. Immunocytochemistry results showed higher CDH1 expression in this group. Sertoli co-culture had no effects on SSC proliferation. Eight weeks after transplantation, injected cells were observed at the base of the seminiferous tubules and in the recipient testes. The number of spermatogonia and the mass of the testes were higher in transplanted testes relative to the control group. It seems that transplantation of these cells can be useful in infertility treatment.

Male Infertility Is Responsible for Nearly Half of the Extinction Observed in the Mouse Collaborative Cross

The goal of the Collaborative Cross (CC) project was to generate and distribute over 1000 independent mouse recombinant inbred strains derived from eight inbred founders. With inbreeding nearly complete, we estimated the extinction rate among CC lines at a remarkable 95%, which is substantially higher than in the derivation of other mouse recombinant inbred populations. Here, we report genome-wide allele frequencies in 347 extinct CC lines. Contrary to expectations, autosomes had equal allelic contributions from the eight founders, but chromosome X had significantly lower allelic contributions from the two inbred founders with underrepresented subspecific origins (PWK/PhJ and CAST/EiJ). By comparing extinct CC lines to living CC strains, we conclude that a complex genetic architecture is driving extinction, and selection pressures are different on the autosomes and chromosome X Male infertility played a large role in extinction as 47% of extinct lines had males that were infertile. Males from extinct lines had high variability in reproductive organ size, low sperm counts, low sperm motility, and a high rate of vacuolization of seminiferous tubules. We performed QTL mapping and identified nine genomic regions associated with male fertility and reproductive phenotypes. Many of the allelic effects in the QTL were driven by the two founders with underrepresented subspecific origins, including a QTL on chromosome X for infertility that was driven by the PWK/PhJ haplotype. We also performed the first example of cross validation using complementary CC resources to verify the effect of sperm curvilinear velocity from the PWK/PhJ haplotype on chromosome 2 in an independent population across multiple generations. While selection typically constrains the examination of reproductive traits toward the more fertile alleles, the CC extinct lines provided a unique opportunity to study the genetic architecture of fertility in a widely genetically variable population. We hypothesize that incompatibilities between alleles with different subspecific origins is a key driver of infertility. These results help clarify the factors that drove strain extinction in the CC, reveal the genetic regions associated with poor fertility in the CC, and serve as a resource to further study mammalian infertility.

Methods for the Study of Gonadal Development

Current knowledge on gonadal development and sex determination is the product of many decades of research involving a variety of scientific methods from different biological disciplines such as histology, genetics, biochemistry, and molecular biology. The earliest embryological investigations, followed by the invention of microscopy and staining methods, were based on histological examinations. The most robust development of histological staining techniques occurred in the second half of the nineteenth century and resulted in structural descriptions of gonadogenesis. These first studies on gonadal development were conducted on domesticated animals; however, currently the mouse is the most extensively studied species. The next key point in the study of gonadogenesis was the advancement of methods allowing for the in vitro culture of fetal gonads. For instance, this led to the description of the origin of cell lines forming the gonads. Protein detection using antibodies and immunolabeling methods and the use of reporter genes were also invaluable for developmental studies, enabling the visualization of the formation of gonadal structure. Recently, genetic and molecular biology techniques, especially gene expression analysis, have revolutionized studies on gonadogenesis and have provided insight into the molecular mechanisms that govern this process. The successive invention of new methods is reflected in the progress of research on gonadal development.

The signaling module cAMP/Epac/Rap1/PLCε/IP3 mobilizes acrosomal calcium during sperm exocytosis

Exocytosis of the sperm's single secretory granule, or acrosome, is a regulated exocytosis triggered by components of the egg's investments. In addition to external calcium, sperm exocytosis (termed the acrosome reaction) requires cAMP synthesized endogenously and calcium mobilized from the acrosome through IP3-sensitive channels. The relevant cAMP target is Epac. In the first part of this paper, we present a novel tool (the TAT-cAMP sponge) to investigate cAMP-related signaling pathways in response to progesterone as acrosome reaction trigger. The TAT-cAMP sponge consists of the cAMP-binding sites of protein kinase A regulatory subunit RIβ fused to the protein transduction domain TAT of the human immunodeficiency virus-1. The sponge permeated into sperm, sequestered endogenous cAMP, and blocked exocytosis. Progesterone increased the population of sperm with Rap1-GTP, Rab3-GTP, and Rab27-GTP in the acrosomal region; pretreatment with the TAT-cAMP sponge prevented the activation of all three GTPases. In the second part of this manuscript, we show that phospholipase Cε (PLCε) is required for the acrosome reaction downstream of Rap1 and upstream of intra-acrosomal calcium mobilization. Last, we present direct evidence that cAMP, Epac, Rap1, and PLCε are necessary for calcium mobilization from sperm's secretory granule. In summary, we describe here a pathway that connects cAMP to calcium mobilization from the acrosome during sperm exocytosis. Never before had direct evidence for each step of the cascade been put together in the same study.

Expression of steroidogenesis-related genes in murine male germ cells

For decades, only few tissues and cell types were defined as steroidogenic, capable of de novo steroid synthesis from cholesterol. However, with the refinement of detection methods, several tissues have now been added to the list of steroidogenic tissues. Besides their critical role as long-range acting hormones, steroids are also playing more discreet roles as local mediators and signaling molecules within the tissues they are produced. In testis, steroidogenesis is carried out by the Leydig cells through a broad network of proteins, mediating cholesterol delivery to CYP11A1, the first cytochrome of the steroidogenic cascade, and the sequential action of enzymes insuring the production of active steroids, the main one being testosterone. The knowledge that male germ cells can be directly regulated by steroids and that they express several steroidogenesis-related proteins led us to hypothesize that germ cells could produce steroids, acting as autocrine, intracrine and juxtacrine modulators, as a way to insure synchronized progression within spermatogenic cycles, and preventing inappropriate cell behaviors between neighboring cells. Gene expression and protein analyses of mouse and rat germ cells from neonatal gonocytes to spermatozoa showed that most steroidogenesis-associated genes are expressed in germ cells, showing cell type-, spermatogenic cycle-, and age-specific expression profiles. Highly expressed genes included genes involved in steroidogenesis and other cell functions, such as Acbd1 and 3, Tspo and Vdac1-3, and genes involved in fatty acids metabolism or synthesis, including Hsb17b4 10 and 12, implying broader roles than steroid synthesis in germ cells. These results support the possibility of an additional level of regulation of spermatogenesis exerted between adjacent germ cells.

The Founder Strains of the Collaborative Cross Express a Complex Combination of Advantageous and Deleterious Traits for Male Reproduction

Surveys of inbred strains of mice are standard approaches to determine the heritability and range of phenotypic variation for biomedical traits. In addition, they may lead to the identification of novel phenotypes and models of human disease. Surprisingly, male reproductive phenotypes are among the least-represented traits in the Mouse Phenome Database. Here we report the results of a broad survey of the eight founder inbred strains of both the Collaborative Cross (CC) and the Diversity Outbred populations, two new mouse resources that are being used as platforms for systems genetics and sources of mouse models of human diseases. Our survey includes representatives of the three main subspecies of the house mice and a mix of classical and wild-derived inbred strains. In addition to standard staples of male reproductive phenotyping such as reproductive organ weights, sperm counts, and sperm morphology, our survey includes sperm motility and the first detailed survey of testis histology. As expected for such a broad survey, heritability varies widely among traits. We conclude that although all eight inbred strains are fertile, most display a mix of advantageous and deleterious male reproductive traits. The CAST/EiJ strain is an outlier, with an unusual combination of deleterious male reproductive traits including low sperm counts, high levels of morphologically abnormal sperm, and poor motility. In contrast, sperm from the PWK/PhJ and WSB/EiJ strains had the greatest percentages of normal morphology and vigorous motility. Finally, we report an abnormal testis phenotype that is highly heritable and restricted to the WSB/EiJ strain. This phenotype is characterized by the presence of a large, but variable, number of vacuoles in at least 10% of the seminiferous tubules. The onset of the phenotype between 2 and 3 wk of age is temporally correlated with the formation of the blood-testis barrier. We speculate that this phenotype may play a role in high rates of extinction in the CC project and in the phenotypes associated with speciation in genetic crosses that use the WSB/EiJ strain as representative of the Mus muculus domesticus subspecies.

A gene expression resource generated by genome-wide lacZ profiling in the mouse

Knowledge of the expression profile of a gene is a critical piece of information required to build an understanding of the normal and essential functions of that gene and any role it may play in the development or progression of disease. High-throughput, large-scale efforts are on-going internationally to characterise reporter-tagged knockout mouse lines. As part of that effort, we report an open access adult mouse expression resource, in which the expression profile of 424 genes has been assessed in up to 47 different organs, tissues and sub-structures using a lacZ reporter gene. Many specific and informative expression patterns were noted. Expression was most commonly observed in the testis and brain and was most restricted in white adipose tissue and mammary gland. Over half of the assessed genes presented with an absent or localised expression pattern (categorised as 0-10 positive structures). A link between complexity of expression profile and viability of homozygous null animals was observed; inactivation of genes expressed in ≥ 21 structures was more likely to result in reduced viability by postnatal day 14 compared with more restricted expression profiles. For validation purposes, this mouse expression resource was compared with Bgee, a federated composite of RNA-based expression data sets. Strong agreement was observed, indicating a high degree of specificity in our data. Furthermore, there were 1207 observations of expression of a particular gene in an anatomical structure where Bgee had no data, indicating a large amount of novelty in our data set. Examples of expression data corroborating and extending genotype-phenotype associations and supporting disease gene candidacy are presented to demonstrate the potential of this powerful resource.

GermlncRNA: a unique catalogue of long non-coding RNAs and associated regulations in male germ cell development

Spermatogenic failure is a major cause of male infertility, which affects millions of couples worldwide. Recent discovery of long non-coding RNAs (lncRNAs) as critical regulators in normal and disease development provides new clues for delineating the molecular regulation in male germ cell development. However, few functional lncRNAs have been characterized to date. A major limitation in studying lncRNA in male germ cell development is the absence of germ cell-specific lncRNA annotation. Current lncRNA annotations are assembled by transcriptome data from heterogeneous tissue sources; specific germ cell transcript information of various developmental stages is therefore under-represented, which may lead to biased prediction or fail to identity important germ cell-specific lncRNAs. GermlncRNA provides the first comprehensive web-based and open-access lncRNA catalogue for three key male germ cell stages, including type A spermatogonia, pachytene spermatocytes and round spermatids. This information has been developed by integrating male germ transcriptome resources derived from RNA-Seq, tiling microarray and GermSAGE. Characterizations on lncRNA-associated regulatory features, potential coding gene and microRNA targets are also provided. Search results from GermlncRNA can be exported to Galaxy for downstream analysis or downloaded locally. Taken together, GermlncRNA offers a new avenue to better understand the role of lncRNAs and associated targets during spermatogenesis.

Database URL: http://germlncrna.cbiit.cuhk.edu.hk/

The proteins of exocytosis: lessons from the sperm model

Exocytosis is a highly regulated process that consists of multiple functionally, kinetically and/or morphologically definable stages such as recruitment, targeting, tethering and docking of secretory vesicles with the plasma membrane, priming of the fusion machinery and calcium-triggered membrane fusion. After fusion, the membrane around the secretory vesicle is incorporated into the plasma membrane and the granule releases its contents. The proteins involved in these processes belong to several highly conserved families: Rab GTPases, SNAREs (soluble NSF-attachment protein receptors), α-SNAP (α-NSF attachment protein), NSF (N-ethylmaleimide-sensitive factor), Munc13 and -18, complexins and synaptotagmins. In the present article, the molecules of exocytosis are reviewed, using human sperm as a model system. Sperm exocytosis is driven by isoforms of the same proteinaceous fusion machinery mentioned above, with their functions orchestrated in a hierarchically organized and unidirectional signalling cascade. In addition to the universal exocytosis regulator calcium, this cascade includes other second messengers such as diacylglycerol, inositol 1,4,5-trisphosphate and cAMP, as well as the enzymes that synthesize them and their target proteins. Of special interest is the cAMP-binding protein Epac (exchange protein directly activated by cAMP) due in part to its enzymatic activity towards Rap. The activation of Epac and Rap leads to a highly localized calcium signal which, together with assembly of the SNARE complex, governs the final stages of exocytosis. The source of this releasable calcium is the secretory granule itself.

A lacZ reporter gene expression atlas for 313 adult KOMP mutant mouse lines

Expression of the bacterial beta-galactosidase reporter gene (lacZ) in the vector used for the Knockout Mouse Project (KOMP) is driven by the endogenous promoter of the target gene. In tissues from KOMP mice, histochemical staining for LacZ enzyme activity can be used to determine gene expression patterns. With this technique, we have produced a comprehensive resource of gene expression using both whole mount (WM) and frozen section (FS) LacZ staining in 313 unique KOMP mutant mouse lines. Of these, ∼80% of mutants showed specific staining in one or more tissues, while ∼20% showed no specific staining, ∼13% had staining in only one tissue, and ∼25% had staining in >6 tissues. The highest frequency of specific staining occurred in the brain (∼50%), male gonads (42%), and kidney (39%). The WM method was useful for rapidly identifying whole organ and some substructure staining, while the FS method often revealed substructure and cellular staining specificity. Both staining methods had >90% repeatability in biological replicates. Nonspecific LacZ staining occurs in some tissues due to the presence of bacteria or endogenous enzyme activity. However, this can be effectively distinguished from reporter gene activity by the combination of the WM and FS methods. After careful annotation, LacZ staining patterns in a high percentage of mutants revealed a unique structure-function not previously reported for many of these genes. The validation of methods for LacZ staining, annotation, and expression analysis reported here provides unique insights into the function of genes for which little is currently known.

Long noncoding RNAs in spermatogenesis: insights from recent high-throughput transcriptome studies

Spermatogenesis is a complex developmental process in which undifferentiated spermatogonia are differentiated into spermatocytes and spermatids through two rounds of meiotic division and finally giving rise to mature spermatozoa (sperm). These processes involve many testis- or male germ cell-specific gene products that undergo strict developmental regulations. As a result, identifying critical, regulatory genes controlling spermatogenesis provide the clues not only to the regulatory mechanism of spermatogenesis at the molecular level, but also to the identification of candidate genes for infertility or contraceptives development. Despite the biological importance in male germ cell development, the underlying mechanisms of stage-specific gene regulation and cellular transition during spermatogenesis remain largely elusive. Previous genomic studies on transcriptome profiling were largely limited to protein-coding genes. Importantly, protein-coding genes only account for a small percentage of transcriptome; the majority are noncoding transcripts that do not translate into proteins. Although small noncoding RNAs (ncRNAs) such as microRNAs, siRNAs, and Piwi-interacting RNAs are extensively investigated in male germ cell development, the role of long ncRNAs (lncRNAs), commonly defined as ncRNAs longer than 200 bp, is relatively unexplored. Herein, we summarize recent transcriptome studies on spermatogenesis and show examples that a subset of noncoding transcript population, known as lncRNAs, constitutes a novel regulatory target in spermatogenesis.

Transcriptome profiling of the murine testis during the first wave of spermatogenesis

Correct gene expression patterns form the basis for male germ cell differentiation and male fertility. Although previous studies have elucidated the importance of testis specific gene expression, the exact transcripts and comprehensive gene expression patterns remain unknown. Large scale sequencing techniques have enabled cost effective analysis of gene expression and isoform studies. Using the SOLiD 4 next-generation sequencing platform we have investigated the gene expression patterns at five different time points during the first wave on murine spermatogenesis. Our results highlight the upregulation of spermatogenesis related biological processes and associated cellular components. Elucidation of differential gene expression at important time points during the sperm development emphasizes the importance of correct timing of gene expression within biological processes. Differential gene level expression was analyzed with R/Bioconductor’s Limma package and isoform analysis was conducted with the Cufflinks pipeline. At gene level total of 2494 differentially expressed genes were identified and Cufflinks characterized over 160 000 gene isoforms, of which 29% were novel transcripts assigned to known genes. Isoforms were detected for 57% of expressed genes and in a total over 26 000 genes were expressed in the testis. Differential promoter and transcription start site usage appears also to play a role in regulation of gene expression during spermatogenesis. Furthermore, we identified 947 upregulated long non-coding RNAs during the first wave of spermatogenesis. These RNAs appeared to be highly specific to different time points. Transcriptomic analysis of testis tissue samples is highly informative due to the large number of expressed genes and identified isoforms. Our study provides a very valuable basis for investigation of gene isoforms and regulation and factors contributing to male fertility.

Transcriptome profiling of the developing postnatal mouse testis using next-generation sequencing

Mammalian testis development is a complex and highly sophisticated process. To study the dynamic change of normal testis development at the transcriptional level, we investigated mouse testes at three postnatal ages: 6 days postnatal, 4 weeks old, and 10 weeks old, representing infant (PN1), juvenile (PN2), and adult (PN3) stages, respectively. Using ultra high-throughput RNA sequencing (RNA-seq) technology, we obtained 211 million reads with a length of 35 bp. We identified 18837 genes that were expressed in mouse testes, and found that genes expressed at the highest level were involved in spermatogenesis. The gene expression pattern in PN1 was distinct from that in PN2 and PN3, which indicates that spermatogenesis has commenced in PN2. We analyzed a large number of genes related to spermatogenesis and somatic development of the testis, which play important roles at different developmental stages. We also found that the MAPK, Hedgehog, and Wnt signaling pathways were significantly involved at different developmental stages. These findings further our understanding of the molecular mechanisms that regulate testis development. Our study also demonstrates significant advantages of RNA-seq technology for studying transcriptome during development.

Functional characterization of bitter-taste receptors expressed in mammalian testis

Mammalian spermatogenesis and sperm maturation are susceptible to the effects of internal and external factors. However, how male germ cells interact with and respond to these elements including those potentially toxic substances is poorly understood. Here, we show that many bitter-taste receptors (T2rs), which are believed to function as gatekeepers in the oral cavity to detect and innately prevent the ingestion of poisonous bitter-tasting compounds, are expressed in mouse seminiferous tubules. Our in situ hybridization results indicate that Tas2r transcripts are expressed postmeiotically. Functional analysis showed that mouse spermatids and spermatozoa responded to both naturally occurring and synthetic bitter-tasting compounds by increasing intracellular free calcium concentrations, and individual male germ cells exhibited different ligand-activation profiles, indicating that each cell may express a unique subset of T2r receptors. These calcium responses could be suppressed by a specific bitter-tastant blocker or abolished by the knockout of the gene for the G protein subunit α-gustducin. Taken together, our data strongly suggest that male germ cells, like taste bud cells in the oral cavity and solitary chemosensory cells in the airway, utilize T2r receptors to sense chemicals in the milieu that may affect sperm behavior and fertilization.

Identification of novel long noncoding RNA transcripts in male germ cells

Emerging evidence from these studies suggested that the male germ cell transcriptome is more complex than previously envisioned. In addition to protein-coding genes, the transcriptome also encodes a significant number of nonprotein-coding transcripts. These noncoding (nc) RNAs appear to be involved in a variety of cellular activities, ranging from simple housekeeping to complex regulatory functions. A class of ncRNAs known as long ncRNAs (lncRNAs) were recently shown to be expressed in a developmentally regulated manner during brain and embryonic stem cell development. This protocol aims to predict and identify potential lncRNA candidates using Serial Analysis of Gene Expression (SAGE) data. We also illustrate how to validate the potential lncRNAs by expression analyses using real-time PCR and Northern Blot. Potential lncRNA candidates in male germ cells are identified using our previously established male germ cell SAGE database (GermSAGE).

Mouse models in male fertility research

Limited knowledge of the genetic causes of male infertility has resulted in few treatment and targeted therapeutic options. Although the ideal approach to identify infertility causing mutations is to conduct studies in the human population, this approach has progressed slowly due to the limitations described herein. Given the complexity of male fertility, the entire process cannot be modeled in vitro. As such, animal models, in particular mouse models, provide a valuable alternative for gene identification and experimentation. Since the introduction of molecular biology and recent advances in animal model production, there has been a substantial acceleration in the identification and characterization of genes associated with many diseases, including infertility. Three major types of mouse models are commonly used in biomedical research, including knockout/knockin/gene-trapped, transgenic and chemical-induced point mutant mice. Using these mouse models, over 400 genes essential for male fertility have been revealed. It has, however, been estimated that thousands of genes are involved in the regulation of the complex process of male fertility, as many such genes remain to be characterized. The current review is by no means a comprehensive list of these mouse models, rather it contains examples of how mouse models have advanced our knowledge of post-natal germ cell development and male fertility regulation.

Testicular postgenomics: targeting the regulation of spermatogenesis

Sperm are, arguably, the most differentiated cells produced within the body of any given species. This is owing to the fact that spermatogenesis is an intricate and highly specialized process evolved to suit the individual particularities of each sexual species. Despite a vast diversity in method, the aim of spermatogenesis is always the same, the idealized transmission of genetic patrimony. Towards this goal certain requirements must always be met, such as a relative twofold reduction in ploidy, repackaging of the chromatin for transport and specialized enhancements for cell motility, recognition and fusion. In the past 20 years, the study of molecular networks coordinating male germ cell development, particularly in mammals, has become more and more facilitated thanks to large-scale analyses of genome expression. Such postgenomic endeavors have generated landscapes of data for both fundamental and clinical reproductive biology. Continuous, large-scale integration analyses of these datasets are undertaken which provide access to very precise information on a myriad of biomolecules. This review presents commonly used transcriptomic and proteomic workflows applied to various testicular germ cell studies. We will also provide a general overview of the technical possibilities available to reproductive genomic biologists, noting the advantages and drawbacks of each technique.

Stage-specific gene expression during spermatogenesis in the dogfish (Scyliorhinus canicula)

In the dogfish testis, the cystic arrangement and polarization of germ cell stages make it possible to observe all stages of spermatogenesis in a single transverse section. By taking advantage of the zonation of this organ, we have used suppressive subtractive libraries construction, real-time PCR, and in situ hybridization to identify 32 dogfish genes showing differential expressions during spermatogenesis. These include homologs of genes already known to be expressed in the vertebrate testis, but found here to be specifically expressed either in pre-meiotic and/or meiotic zones (ribosomal protein S8, high-mobility group box 3, ubiquitin carboxyl-terminal esterase L3, 20beta-hydroxysteroid dehydrogenase, or cyclophilin B) or in post-meiotic zone (speriolin, Soggy, zinc finger protein 474, calreticulin, or phospholipase c-zeta). We also report, for the first time, testis-specific expression patterns for dogfish genes coding for A-kinase anchor protein 5, ring finger protein 152, or F-box only protein 7. Finally, the study highlights the differential expression of new sequences whose identity remains to be assessed. This study provides the first molecular characterization of spermatogenesis in a chondrichthyan, a key species to gain insight into the evolution of this process in gnathostomes.

Age affects gene expression in mouse spermatogonial stem/progenitor cells

Spermatogenesis in man starts with spermatogonial stem cells (SSCs), and leads to the production of sperm in approximately 64 days, common to old and young men. Sperm from elderly men are functional and able to fertilize eggs and produce offspring, even though daily sperm production is more than 50% lower and damage to sperm DNA is significantly higher in older men than in those who are younger. Our hypothesis is that the SSC/spermatogonial progenitors themselves age. To test this hypothesis, we studied the gene expression profile of mouse SSC/progenitor cells at several ages using microarrays. After sequential enzyme dispersion, we purified the SSC/progenitors with immunomagnetic cell sorting using an antibody to GFRA1, a known SSC/progenitor cell marker. RNA was isolated and used for the in vitro synthesis of amplified and labeled cRNAs that were hybridized to the Affymetrix mouse genome microarrays. The experiments were repeated twice with different cell preparations, and statistically significant results are presented. Quantitative RT-PCR analysis was used to confirm the microarray results. Comparison of four age groups (6 days, 21 days, 60 days, and 8 months old) showed a number of genes that were expressed specifically in the older mice. Two of them (i.e. Icam1 and Selp) have also been shown to mark aging hematopoietic stem cells. On the other hand, the expression levels of the genes encoding the SSC markers Gfra1 and Plzf did not seem to be significantly altered by age, indicating that age affects only certain SSC/progenitor properties.

Screening of genes involved in chromosome segregation during meiosis I: toward the identification of genes responsible for infertility in humans

Prophase I of male meiosis during early spermatogenesis involves dynamic chromosome segregation processes, including synapsis, meiotic recombination and cohesion. Genetic defects in the genes that participate in these processes consistently cause reproduction failure in mice. To identify candidate genes responsible for infertility in humans, we performed gene expression profiling of mouse spermatogenic cells undergoing meiotic prophase I. Cell fractions enriched in spermatogonia, leptotene/zygotene spermatocytes or pachytene spermatocytes from developing mouse testis were separately isolated by density gradient sedimentation and subjected to microarray analysis. A total of 726 genes were identified that were upregulated in leptotene/zygotene spermatocytes. To evaluate the screening efficiency for meiosis-specific genes, we randomly selected 12 genes from this gene set and characterized each gene product using reverse transcription (RT)-PCR of RNA from gonadal tissues, in situ hybridization on testicular tissue sections and subcellular localization analysis of the encoded protein. Four of the 12 genes were confirmed as genes expressed in meiotic stage and 2 of these 4 genes were novel, previously uncharacterized genes. Among the three confirmation methods that were used, RT-PCR appeared to be the most efficient method for further screening. These 726 candidates for human infertility genes might serve as a useful resource for next-generation sequencing combined with exon capture by microarray.

GonadSAGE: a comprehensive SAGE database for transcript discovery on male embryonic gonad development

SUMMARY

Serial analysis of gene expression (SAGE) provides an alternative, with additional advantages, to microarray gene expression studies. GonadSAGE is the first publicly available web-based SAGE database on male gonad development that covers six male mouse embryonic gonad stages, including E10.5, E11.5, E12.5, E13.5, E15.5 and E17.5. The sequence coverage of each SAGE library is beyond 150K, 'which is the most extensive sequence-based male gonadal transcriptome to date'. An interactive web interface with customizable parameters is provided for analyzing male gonad transcriptome information. Furthermore, the data can be visualized and analyzed with the other genomic features in the UCSC genome browser. It represents an integrated platform that leads to a better understanding of male gonad development, and allows discovery of related novel targets and regulatory pathways.

Histone gene expression profile during spermatogenesis

Gene bank search and reverse transcription-polymerase chain reaction were used to analyze the expression profile for histone genes during spermatogenesis. The objective of this study was to provide a systematic view of histone gene expression.

Modified PCR methods for 3' end amplification from serial analysis of gene expression (SAGE) tags

Serial analysis of gene expression (SAGE) is a powerful technique to study gene expression at the genome level. However, a disadvantage of the shortness of SAGE tags is that it prevents further study of SAGE library data, thus limiting extensive application of the SAGE method in gene expression studies. However, this problem can be solved by extension of the SAGE tags to 3' cDNAs. Therefore, several methods based on PCR have been developed to generate a 3' longer fragment cDNA corresponding to a SAGE tag. The list of modified methods is extensive, and includes rapid RT-PCR analysis of unknown SAGE tags (RAST-PCR), generation of longer cDNA fragments from SAGE tags for gene identification (GLGI), a high-throughput GLGI procedure, reverse SAGE (rSAGE), two-step analysis of unknown SAGE tags (TSAT-PCR), etc. These procedures are constantly being updated because they have characteristics and advantages that can be shared. Development of these methods has promoted the widespread use of the SAGE technique, and has accelerated the speed of studies of large-scale gene expression.

Developmental staging of male murine embryonic gonad by SAGE analysis

Despite the identification of key genes such as Sry integral to embryonic gonadal development, the genomic classification and identification of chromosomal activation of this process is still poorly understood. To better understand the genetic regulation of gonadal development, we performed Serial Analysis of Gene Expression (SAGE) to profile the genes and novel transcripts, and an average of 152,000 tags from male embryonic gonads at E10.5 (embryonic day 10.5), E11.5, E12.5, E13.5, E15.5 and E17.5 were analyzed. A total of 275,583 non-singleton tags that do not map to any annotated sequence were identified in the six gonad libraries, and 47,255 tags were mapped to 24,975 annotated sequences, among which 987 sequences were uncharacterized. Utilizing an unsupervised pattern identification technique, we established molecular staging of male gonadal development. Rather than providing a static descriptive analysis, we developed algorithms to cluster the SAGE data and assign SAGE tags to a corresponding chromosomal position; these data are displayed in chromosome graphic format. A prominent increase in global genomic activity from E10.5 to E17.5 was observed. Important chromosomal regions related to the developmental processes were identified and validated based on established mouse models with developmental disorders. These regions may represent markers for early diagnosis for disorders of male gonad development as well as potential treatment targets.

Genomic landscape of developing male germ cells

Spermatogenesis is a highly orchestrated developmental process by which spermatogonia develop into mature spermatozoa. This process involves many testis- or male germ cell-specific gene products whose expressions are strictly regulated. In the past decade the advent of high-throughput gene expression analytical techniques has made functional genomic studies of this process, particularly in model animals such as mice and rats, feasible and practical. These studies have just begun to reveal the complexity of the genomic landscape of the developing male germ cells. Over 50% of the mouse and rat genome are expressed during testicular development. Among transcripts present in germ cells, 40% - 60% are uncharacterized. A number of genes, and consequently their associated biological pathways, are differentially expressed at different stages of spermatogenesis. Developing male germ cells present a rich repertoire of genetic processes. Tissue-specific as well as spermatogenesis stage-specific alternative splicing of genes exemplifies the complexity of genome expression. In addition to this layer of control, discoveries of abundant presence of antisense transcripts, expressed psuedogenes, non-coding RNAs (ncRNA) including long ncRNAs, microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), and retrogenes all point to the presence of multiple layers of expression and functional regulation in male germ cells. It is anticipated that application of systems biology approaches will further our understanding of the regulatory mechanism of spermatogenesis.

Cloning, characterization, and expression analysis of the novel acetyltransferase retrogene Ard1b in the mouse

N-alpha-terminal acetylation is a modification process that occurs cotranslationally on most eukaryotic proteins. The major enzyme responsible for this process, N-alpha-terminal acetyltransferase, is composed of the catalytic subunit ARD1A and the auxiliary subunit NAT1. We cloned, characterized, and studied the expression pattern of Ard1b (also known as Ard2), a novel homolog of the mouse Ard1a. Comparison of the genomic structures suggests that the autosomal Ard1b is a retroposed copy of the X-linked Ard1a. Expression analyses demonstrated a testis predominance of Ard1b. A reciprocal expression pattern between Ard1a and Ard1b is also observed during spermatogenesis, suggesting that Ard1b is expressed to compensate for the loss of Ard1a starting from meiosis. Both ARD1A and ARD1B can interact with NAT1 to constitute a functional N-alpha-terminal acetyltransferase in vitro. The expression of ARD1B protein can be detected in mouse testes but is delayed until the first appearance of round spermatids. In a cell culture model, the inclusion of the long 3' untranslated region of Ard1b leads to reduction of luciferase reporter activity, which implicates its role in translational repression of Ard1b during spermatogenesis. Our results suggest that ARD1B may have an important role in the later course of the spermatogenic process.

The molecular signature of spermatogonial stem/progenitor cells in the 6-day-old mouse testis

To characterize the molecular phenotype of spermatogonial stem cells (SSCs), we examined genes that are differentially expressed in the stem/progenitor spermatogonia compared to nonstem spermatogonia. We isolated type A spermatogonia (stem and nonstem type A) from 6-day-old mice using sedimentation velocity at unit gravity and further selected the stem/progenitor cell subpopulation by magnetic activated cell sorting with an antibody to GDNF-receptor-alpha-1 (GFRA1). It has been previously shown that GFRA1 is expressed in SSCs and is required for their stemness. The purity of the isolated cells was approximately 95% to 99% as indicated by immunocytochemistry using anti-GFRA1. Comparison of GFRA1-positive and GFRA1-negative spermatogonia by microarray analysis revealed 99 known genes and 12 uncharacterized transcripts that are overexpressed in the former cell population with a >2-fold change. Interestingly, the highest level of overexpression was observed for Csf1r, encoding the receptor for macrophage colony-stimulating factor (M-CSF, official symbol CSF1), which has a well-established role in the regulation of myeloid progenitor cells. Analysis of our microarray data with a bioinformatics software program (Ingenuity Systems) revealed the potential role of various signaling pathways in stem/progenitor spermatogonia and suggested a common pathway for GFRA1 and CSF1R that may lead to their proliferation. Further investigation to test this hypothesis has shown that CSF1 promotes cell proliferation in primary cultures of the isolated type A spermatogonia and in the spermatogonial-derived stem cell line C18-4. Semiquantitative RT-PCR and immunohistochemistry confirmed the previously mentioned microarray data. Collectively, this study provides novel molecular signatures for stem/progenitor spermatogonia and demonstrates a role for CSF1/CSF1R signaling in regulating their proliferation.

3640 unique EST clusters from the medaka testis and their potential use for identifying conserved testicular gene expression in fish and mammals

Background

The fish medaka is the first vertebrate capable of full spermatogenesis in vitro from self-renewing spermatogonial stem cells to motile test-tube sperm. Precise staging and molecular dissection of this process has been hampered by the lack of suitable molecular markers.

Methodology and Principal Findings

We have generated a normalized medaka testis cDNA library and obtained 7040 high quality sequences representing 3641 unique gene clusters. Among these, 1197 unique clusters are homologous to known genes, and 2444 appear to be novel genes. Ontology analysis shows that the 1197 gene products are implicated in diverse molecular and cellular processes. These genes include markers for all major types of testicular somatic and germ cells. Furthermore, markers were identified for major spermatogenic stages ranging from spermatogonial stem cell self-renewal to meiosis entry, progression and completion. Intriguingly, the medaka testis expresses at least 13 homologs of the 33 mouse X-chromosomal genes that are enriched in the testis. More importantly, we show that key components of several signaling pathways known to be important for testicular function in mammals are well represented in the medaka testicular EST collection.

Conclusions/Significance

Medaka exhibits a considerable similarity in testicular gene expression to mammals. The medaka testicular EST collection we obtained has wide range coverage and will not only consolidate our knowledge on the comparative analysis of known genes' functions in the testis but also provide a rich resource to dissect molecular events and mechanism of spermatogenesis in vivo and in vitro in medaka as an excellent vertebrate model.

GermSAGE: a comprehensive SAGE database for transcript discovery on male germ cell development

GermSAGE is a comprehensive web-based database generated by Serial Analysis of Gene Expression (SAGE) representing major stages in mouse male germ cell development, with 150 000 sequence tags in each SAGE library. A total of 452 095 tags derived from type A spermatogonia (Spga), pachytene spermatocytes (Spcy) and round spermatids (Sptd) were included. GermSAGE provides web-based tools for browsing, comparing and searching male germ cell transcriptome data at different stages with customizable searching parameters. The data can be visualized in a tabulated format or further analyzed by aligning with various annotations available in the UCSC genome browser. This flexible platform will be useful for gaining better understanding of the genetic networks that regulate spermatogonial cell renewal and differentiation, and will allow novel gene discovery. GermSAGE is freely available at http://germsage.nichd.nih.gov/

Serial analysis of gene expression

Serial analysis of gene expression (SAGE) is a method used to obtain comprehensive, unbiased and quantitative gene-expression profiles. Its major advantage over arrays is that it does not require a priori knowledge of the genes to be analyzed and reflects absolute mRNA levels. Since the original SAGE protocol was developed in a short-tag (10-bp) format, several modifications have been made to produce longer SAGE tags for more precise gene identification and to decrease the amount of starting material necessary. Several SAGE-like methods have also been developed for the genome-wide analysis of DNA copy-number changes and methylation patterns, chromatin structure and transcription factor targets. In this protocol, we describe the 17-bp longSAGE method for transcriptome profiling optimized for a small amount of starting material. The generation of such libraries can be completed in 7-10 d, whereas sequencing and data analysis require an additional 2-3 wk.

Altered expression of genes involved in regulation of vitamin A metabolism, solute transportation, and cytoskeletal function in the androgen-insensitive tfm mouse testis

Androgens are essential for the development and maintenance of spermatogenesis, but the underlying mechanisms of androgen action in the testis remain unclear. To help clarify these mechanisms, gene expression was measured in testes of pubertal (20 d old), androgen-insensitive, testicular feminized (Tfm) mice and in normal controls. Using microarrays (Affymetrix chips 430A and 430B), initial data identified a large number of genes down-regulated in the Tfm testis (>4700). These genes were largely of germ cell origin, reflecting the arrest of spermatogenesis that is apparent in the 20-d-old Tfm testis. Subsequent screening in vitro and in silico of this gene set identified 20 genes of a somatic tubular origin that were significantly down-regulated in the Tfm testis and six genes that were significantly up-regulated. Altered expression of these genes was confirmed by real-time PCR, and genes down-regulated in the Tfm testis were shown to be up-regulated in testes of hypogonadal (hpg) mice treated with androgen. In a developmental study using real-time PCR most of the regulated genes showed normal expression during fetal and neonatal development and deviated from control only between 10 and 20 d. In all cases, expression was also reduced in the adult, although interpretation is more complex because of the inherent cryptorchidism in the adult Tfm mouse. Of the total number of somatic genes showing differential expression in the Tfm testis, 50% were associated with three separate groups of genes involved in regulation of vitamin A metabolism, solute transportation, and cytoskeletal function. Thus, effects of androgens on tubular function and spermatogenesis may be mediated in part through regulation of the tubular environment and control of retinoic acid concentrations.

Distribution of menin-occupied regions in chromatin specifies a broad role of menin in transcriptional regulation

Menin is the protein product of the MEN1 tumor-suppressor gene; one allele of MEN1 is inactivated in the germ line of patients with "multiple endocrine neoplasia type 1" (MEN1) cancer syndrome. Menin interacts with several proteins involved in transcriptional regulation. RNA expression analyses have identified several menin-regulated genes that could represent proximal or distal interaction sites for menin. This report presents a substantial and unbiased sampling of menin-occupied chromatin regions using Serial Analysis of Chromatin Occupancy; this method combines chromatin immuno-precipitation with Serial Analysis of Gene Expression. Hundreds of menin-occupied genomic sites were identified in promoter regions (32% of menin-occupied loci), near the 3' end of genes (14%), or inside genes (21%), extending other data about menin recruitments to many sites of transcriptional activity. A large number of menin-occupied sites (33%) were located outside known gene regions. Additional annotation of the human genome could help in identifying genes at these loci, or these might be gene-free regions of the genome where menin occupancy could play some structural or regulatory role. Menin occupancy at many intragenic positions distant from the core promoter reveals an unexpected type of menin target region at many loci in the genome. These unbiased data also suggest that menin could play a broad role in transcriptional regulation.

The complexity of antisense transcription revealed by the study of developing male germ cells

Computational analyses have identified the widespread occurrence of antisense transcripts in the human and the mouse genome. However, the structure and the origin of the majority of the antisense transcripts are unknown. The presence of antisense transcripts for 19 of 64 differentially expressed genes during mouse spermatogenesis was demonstrated with orientation-specific RT-PCR. These antisense transcripts were derived from a wide variety of origins, including processed sense transcripts, intronic and exonic sequences of a single gene or multiple genes, intergenic sequences, and pseudogenes. They underwent normal and alternative splicing, 5' capping, and 3' polyadenylation, similar to the sense transcripts. There were also antisense transcripts that were not capped and/or polyadenylated. The testicular levels of the sense transcripts were higher than those of the antisense transcripts in all cases, while the relative expression in nontesticular tissues was variable. Thus antisense transcripts have complex origins and structures and the sense and antisense transcripts can be regulated independently.

Transcriptome analyses of male germ cells with serial analysis of gene expression (SAGE)

Serial analysis of gene expression (SAGE) provides an alternative with additional advantages to microarrays for studying gene expression during spermatogenesis. The digitized transcriptome provided by SAGE of purified mouse germ cells identified 27,504 species of transcripts expressed in type A spermatogonia, pachytene spermatocytes, and round spermatids. Over 2700 of these transcripts were novel. Computational analyses allowed the identification of clusters of co-regulated genes, cell-specific promoter modules, cell-specific biological processes, as well as "preferential" biological networks in different cell types. These analyses provided potential drug targets for interference of specific pathways at different stages of spermatogenesis. Analyses of the transcriptomes revealed the prominent role of cytochrome c oxidase in germ cells and suggest a novel role for this enzyme in cytochrome c-mediated apoptosis in spermatogonia. A number of genes were shown to undergo differential splicing during spermatogenesis giving rise to cell-specific splice variants.

Expression profiling reveals meiotic male germ cell mRNAs that are translationally up- and down-regulated

Gametes rely heavily on posttranscriptional control mechanisms to regulate their differentiation. In eggs, maternal mRNAs are stored and selectively activated during development. In the male, transcription ceases during spermiogenesis, necessitating the posttranscriptional regulation of many paternal mRNAs required for spermatozoan assembly and function. To date, most of the testicular mRNAs known to be translationally regulated are initially transcribed in postmeiotic cells. Because protein synthesis occurs on polysomes and translationally inactive mRNAs are sequestered as ribonucleoproteins (RNPs), movement of mRNAs between these fractions is indicative of translational up- and down-regulation. Here, we use microarrays to analyze mRNAs in RNPs and polysomes from testis extracts of prepuberal and adult mice to characterize the translation state of individual mRNAs as spermatogenesis proceeds. Consistent with published reports, many of the translationally delayed postmeiotic mRNAs shift from the RNPs into the polysomes, establishing the validity of this approach. In addition, we detect another 742 mouse testicular transcripts that show dramatic shifts between RNPs and polysomes. One subgroup of 35 genes containing the known, translationally delayed phosphoglycerate kinase 2 (Pgk2) is initially transcribed during meiosis and is translated in later-stage cells. Another subgroup of 82 meiotically expressed genes is translationally down-regulated late in spermatogenesis. This high-throughput approach defines the changing translation patterns of populations of genes as male germ cells differentiate and identifies groups of meiotic transcripts that are translationally up- and down-regulated.

Application of transcriptional and biological network analyses in mouse germ-cell transcriptomes

Serial analysis of gene expression (SAGE) provides a global analysis platform for profiling mRNA populations present in cells of interest without the constraint of gene selection and the ambiguous nature of data obtained. However, most of the reports on SAGE and germ cell development are limited to descriptive analyses. Here, we report a series of bioinformatic analyses using recently published SAGE data on the transcriptome of mouse type A spermatogonia (Spga), pachytene spermatocytes (Spcy), and round spermatids (Sptd). Tags with a total count of > or =20 in three SAGE libraries were examined. Our aim was to identify and discover potential transcriptional regulators and pathways involved at different stages of spermatogenesis. Unsupervised hierarchical clustering based on tag expression and Gene Ontology analysis were applied to identify genes and biological processes overrepresented at a particular stage of development. The 5' cis-regulatory elements were examined for common regulators in different functional clusters. Potential biological networks were also constructed to reveal the link between the gene candidates. Biological pathways related to the three germ cell stages were constructed. A number of known transcription regulators in spermatogenesis, including NF-kappaB, SP1, AP-1, and EGR, were identified. Novel promoter elements such as the E box in Spga-specific genes, GATA in Spcy-specific genes, and GKLF in Sptd-specific genes were also observed. Taken together, our approach is reliable and provides a foundation for the generation of novel biological hypotheses for studying spermatogenesis.

Expression of Col1a1, Col1a2 and procollagen I in germ cells of immature and adult mouse testis

The objective of this study was to compare the expression of Col1a1, Col1a2, and procollagen I in the seminiferous tubules of immature and adult mice and to characterize the cellular expression pattern of procollagen I in germ cells during spermatogenesis in order to provide necessary groundwork for further functional studies in the process of spermatogenesis. Microarray analysis demonstrated that Col1a1 and Col1a2 were abundantly expressed in the seminiferous tubules of 6-day-old mice compared with 60-day-old mice, and the expression levels of Col1a1 and Col1a2 mRNA were validated using a semi-quantitative RT-PCR assay. Western blot analysis further confirmed that procollagen I was expressed at a higher level in the seminiferous tubules of 6-day-old mice compared with 60-day-old mice. Immunohistochemical analysis revealed that type A spermatogonia were positive for procollagen I in the testis of 6-day-old mice, whereas Sertoli cells were negative for this protein. The in vivo procollagen I staining in type A spermatogonia was corroborated in spermatogonia exhibiting a high potential for proliferation and the ability to form germ cell colonies in in vitro culture. Moreover, procollagen I was also detected in type A spermatogonia, intermediate spermatogonia, type B spermatogonia, and preleptotene spermatocytes in the adult mouse testes, but positive staining disappeared in more differentiated germ cell lineages detaching from the basement membrane, including leptotene spermatocytes, pachytene spermatocytes, round spermatids and elongated spermatids. These data suggest that Col1a1, Col1a2 and procollagen I are associated with type A spermatogonia and play a potential role in mediating the detachment and migration of germ cells during spermatogenesis.

Expression profiling of purified male germ cells: stage-specific expression patterns related to meiosis and postmeiotic development

Gene expression profiling was performed using the National Institute on Aging 15,000-cDNA microarray to reveal the differential expression pattern of 160 genes between meiotic pachytene spermatocytes and postmeiotic round spermatids of the mouse. Our results indicate that more genes are expressed in spermatids than in spermatocytes. Genes participating in cell cycle regulation and chromatin structure and dynamics are preferentially expressed in spermatocytes, while genes for protein turnover, signal transduction, energy metabolism, and intracellular transport are prevalent in spermatids. This suggests that a switch of functional requirement occurs when meiotic germ cells differentiate into haploid spermatids. Concordant expression patterns were obtained when quantitative real-time polymerase chain reaction was performed to verify the microarray data. Interestingly, the majority of the differentially expressed genes were underrepresented in mitotic type A spermatogonia, and they were preferentially expressed in the testis. Our results suggest that an even higher proportion of the mouse genome is devoted to male gamete development from meiosis than was previously estimated. We also provide evidence that underscores the advantage of using purified germ cells over whole testes in profiling spermatogenic gene expression to identify transcripts that demonstrate stage-specific expression patterns.

Haploinsufficiency of cytochrome P450 17alpha-hydroxylase/17,20 lyase (CYP17) causes infertility in male mice

Cytochrome P450 17alpha-hydroxylase/17,20-lyase (CYP17) is critical in determining cortisol and sex steroid biosynthesis. To investigate how CYP17 functions in vivo, we generated mice with a targeted deletion of CYP17. Although in chimeric mice Leydig cell CYP17 mRNA and intratesticular and circulating testosterone levels were dramatically reduced (80%), the remaining testosterone was sufficient to support spermatogenesis as evidenced by the generation of phenotypical black C57BL/6 mice. However, male chimeras consistently failed to generate heterozygous CYP17 mice and after five matings chimeric mice stopped mating indicating a change in sexual behavior. These results suggested that CYP17 deletion caused a primary phenotype (infertility), probably not due to the anticipated androgen imbalance and a secondary phenotype (change in sexual behavior) due to the androgen imbalance. Surprisingly, CYP17 mRNA was found in mature sperm, and serial analysis of gene expression identified CYP17 mRNA in other testicular germ cells. CYP17 mRNA levels were directly related to percent chimerism. Moreover, more than 50% of the sperm from high-percentage chimeric mice were morphologically abnormal, and half of them failed the swim test. Furthermore, 60% of swimming abnormal sperm was devoid of CYP17. These results suggest that CYP17, in addition to its role in steroidogenesis and androgen formation, is present in germ cells where it is essential for sperm function, and deletion of one allele prevents genetic transmission of mutant and wild-type alleles causing infertility followed by change in sexual behavior due to androgen imbalance.