Gene expression programs in mammalian spermatogenesis.

Mammalian spermatogenesis, probably the most complex of all cellular developmental processes, is an ideal model both for studying the specific mechanism of gametogenesis and for understanding the basic rules governing all developmental processes, as it entails both cell type-specific and housekeeping molecular processes. Spermatogenesis can be viewed as a mission with many tasks to accomplish, and its success is genetically programmed and ensured by the collaboration of a large number of genes. Here, I present an overview of mammalian spermatogenesis and the mechanisms underlying each step in the process, covering the cellular and molecular activities that occur at each developmental stage and emphasizing their gene regulation in light of recent studies.

Chromodomain Protein CDYL Acts as a Crotonyl-CoA Hydratase to Regulate Histone Crotonylation and Spermatogenesis.

Lysine crotonylation (Kcr) is a newly identified histone modification that is associated with active transcription in mammalian cells. Here we report that the chromodomain Y-like transcription corepressor CDYL negatively regulates histone Kcr by acting as a crotonyl-CoA hydratase to convert crotonyl-CoA to ß-hydroxybutyryl-CoA. We showed that the negative regulation of histone Kcr by CDYL is intrinsically linked to its transcription repression activity and functionally implemented in the reactivation of sex chromosome-linked genes in round spermatids and genome-wide histone replacement in elongating spermatids. Significantly, Cdyl transgenic mice manifest dysregulation of histone Kcr and reduction of male fertility with a decreased epididymal sperm count and sperm cell motility. Our study uncovers a biochemical pathway in the regulation of histone Kcr and implicates CDYL-regulated histone Kcr in spermatogenesis, adding to the understanding of the physiology of male reproduction and the mechanism of the spermatogenic failure in AZFc (Azoospermia Factor c)-deleted infertile men.

The microRNA miR-202 prevents precocious spermatogonial differentiation and meiotic initiation during mouse spermatogenesis.

Spermatogonial differentiation and meiotic initiation during spermatogenesis are tightly regulated by a number of genes, including those encoding enzymes for miRNA biogenesis. However, whether and how single miRNAs regulate these processes remain unclear. Here, we report that miR-202, a member of the let-7 family, prevents precocious spermatogonial differentiation and meiotic initiation in spermatogenesis by regulating the timely expression of many genes, including those for key regulators such as STRA8 and DMRT6. In miR-202 knockout (KO) mice, the undifferentiated spermatogonial pool is reduced, accompanied by age-dependent decline of fertility. In KO mice, SYCP3, STRA8 and DMRT6 are expressed earlier than in wild-type littermates, and Dmrt6 mRNA is a direct target of miR-202-5p. Moreover, the precocious spermatogonial differentiation and meiotic initiation were also observed in KO spermatogonial stem cells when cultured and induced in vitro, and could be partially rescued by the knockdown of Dmrt6. Therefore, we have not only shown that miR-202 is a regulator of meiotic initiation but also identified a previously unknown module in the underlying regulatory network.

MicroRNA-202 safeguards meiotic progression by preventing premature SEPARASE-mediated REC8 cleavage.

MicroRNAs (miRNAs) are believed to play important roles in mammalian spermatogenesis but the in vivo functions of single miRNAs in this highly complex developmental process remain unclear. Here, we report that miR-202, a member of the let-7 family, plays an important role in spermatogenesis by phenotypic evaluation of miR-202 knockout (KO) mice. Loss of miR-202 results in spermatocyte apoptosis and perturbation of the zygonema-to-pachynema transition. Multiple processes during meiosis prophase I including synapsis and crossover formation are disrupted, and inter-sister chromatid synapses are detected. Moreover, we demonstrate that Separase mRNA is a miR-202 direct target and provides evidence that miR-202 upregulates REC8 by repressing Separase expression. Therefore, we have identified miR-202 as a new regulating noncoding gene that acts on the established SEPARASE-REC8 axis in meiosis.

The male germline-specific protein MAPS is indispensable for pachynema progression and fertility.

Meiosis is a specialized cell division that creates haploid germ cells from diploid progenitors. Through differential RNA expression analyses, we previously identified a number of mouse genes that were dramatically elevated in spermatocytes, relative to their very low expression in spermatogonia and somatic organs. Here, we investigated in detail 1700102P08Rik, one of these genes, and independently conclude that it encodes a male germline-specific protein, in agreement with a recent report. We demonstrated that it is essential for pachynema progression in spermatocytes and named it male pachynema-specific (MAPS) protein. Mice lacking Maps (Maps-/- ) suffered from pachytene arrest and spermatocyte death, leading to male infertility, whereas female fertility was not affected. Interestingly, pubertal Maps-/- spermatocytes were arrested at early pachytene stage, accompanied by defects in DNA double-strand break (DSB) repair, crossover formation, and XY body formation. In contrast, adult Maps-/- spermatocytes only exhibited partially defective crossover but nonetheless were delayed or failed in progression from early to mid- and late pachytene stage, resulting in cell death. Furthermore, we report a significant transcriptional dysregulation in autosomes and XY chromosomes in both pubertal and adult Maps-/- pachytene spermatocytes, including failed meiotic sex chromosome inactivation (MSCI). Further experiments revealed that MAPS overexpression in vitro dramatically decreased the ubiquitination levels of cellular proteins. Conversely, in Maps-/- pachytene cells, protein ubiquitination was dramatically increased, likely contributing to the large-scale disruption in gene expression in pachytene cells. Thus, MAPS is a protein essential for pachynema progression in male mice, possibly in mammals in general.

The function of Foxo1 in spermatogonia development is independent of PI3K/PTEN signaling.

Spermatogenesis is a highly coordinated process that initiates shortly after birth and continues throughout the lifespan of male animals. Foxo1 is a transcription factor and is involved in many biological processes. It has been reported that the inactivation of Foxo1 in gonocytes during the embryonic stage causes the defects of spermatogenesis. In the present study, we found that the inactivation of Foxo1 in spermatogonia after birth also caused germ cell loss and male infertility. We found that the initiation of meiosis was not affected; however, the germ cell development was arrested after meiosis and lack of mature spermatozoa in the cauda epididymis. We also found that the proliferation of Foxo1-deficient spermatogonia stem cells was significantly reduced under in vitro conditions. Further study revealed that inactivation of Pten in postnatal spermatogonia using Stra8-Cre did not affect germ cell development and the subcellular location of FOXO1 in Pten-deficient spermatogonia. This study demonstrated that Foxo1 was involved in the development of spermatogonia after birth and the function of Foxo1 was probably not regulated by PI3K/PTEN signaling.

The transcription factor SOX30 is a key regulator of mouse spermiogenesis.

The postmeiotic development of male germ cells, also known as spermiogenesis, features the coordinated expression of a large number of spermatid-specific genes. However, only a limited number of key transcription factors have been identified and the underlying regulatory mechanisms remain largely unknown. Here, we report that SOX30, the most-divergent member of the Sry-related high-motility group box (SOX) family of transcription factors, is essential for mouse spermiogenesis. The SOX30 protein was predominantly expressed in spermatids, while its transcription was regulated by retinoic acid and by MYBL1 before and during meiosis. Sox30 knockout mice arrested spermiogenesis at step 3 round spermatids, which underwent apoptosis and abnormal chromocenter formation. We also determined that SOX30 regulated the expression of hundreds of spermatid-specific protein-coding and long non-coding RNA genes. SOX30 bound to the proximal promoter of its own gene and activated its transcription. These results reveal SOX30 as a novel key regulator of spermiogenesis that regulates its own transcription to enforce and activate this meiotic regulatory pathway.

Wt1 directs the lineage specification of sertoli and granulosa cells by repressing Sf1 expression.

Supporting cells (Sertoli and granulosa) and steroidogenic cells (Leydig and theca-interstitium) are two major somatic cell types in mammalian gonads, but the mechanisms that control their differentiation during gonad development remain elusive. In this study, we found that deletion of Wt1 in the ovary after sex determination caused ectopic development of steroidogenic cells at the embryonic stage. Furthermore, differentiation of both Sertoli and granulosa cells was blocked when Wt1 was deleted before sex determination and most genital ridge somatic cells differentiated into steroidogenic cells in both male and female gonads. Further studies revealed that WT1 repressed Sf1 expression by directly binding to the Sf1 promoter region, and the repressive function was completely abolished when WT1 binding sites were mutated. This study demonstrates that Wt1 is required for the lineage specification of both Sertoli and granulosa cells by repressing Sf1 expression. Without Wt1, the expression of Sf1 was upregulated and the somatic cells differentiated into steroidogenic cells instead of supporting cells. Our study uncovers a novel mechanism of somatic cell differentiation during gonad development.

Speedy A-Cdk2 binding mediates initial telomere-nuclear envelope attachment during meiotic prophase I independent of Cdk2 activation.

Telomere attachment to the nuclear envelope (NE) is a prerequisite for chromosome movement during meiotic prophase I that is required for pairing of homologous chromosomes, synapsis, and homologous recombination. Here we show that Speedy A, a noncanonical activator of cyclin-dependent kinases (Cdks), is specifically localized to telomeres in prophase I male and female germ cells in mice, and plays an essential role in the telomere-NE attachment. Deletion of Spdya in mice disrupts telomere-NE attachment, and this impairs homologous pairing and synapsis and leads to zygotene arrest in male and female germ cells. In addition, we have identified a telomere localization domain on Speedy A covering the distal N terminus and the Cdk2-binding Ringo domain, and this domain is essential for the localization of Speedy A to telomeres. Furthermore, we found that the binding of Cdk2 to Speedy A is indispensable for Cdk2's localization on telomeres, suggesting that Speedy A and Cdk2 might be the initial components that are recruited to the NE for forming the meiotic telomere complex. However, Speedy A-Cdk2-mediated telomere-NE attachment is independent of Cdk2 activation. Our results thus indicate that Speedy A and Cdk2 might mediate the initial telomere-NE attachment for the efficient assembly of the telomere complex that is essential for meiotic prophase I progression.

Non-canonical RNA polyadenylation polymerase FAM46C is essential for fastening sperm head and flagellum in mice†.

Family with sequence similarity 46, member C (FAM46C) is a highly conserved non-canonical RNA polyadenylation polymerase that is abundantly expressed in human and mouse testes and is frequently mutated in patients with multiple myeloma. However, its physiological role remains largely unknown. In this study, we found that FAM46C is specifically localized to the manchette of spermatids in mouse testes, a transient microtubule-based structure mainly involved in nuclear shaping and intra-flagellar protein traffic. Gene knockout of FAM46C in mice resulted in male sterility, characterized by the production of headless spermatozoa in testes. Sperm heads were intermittently found in the epididymides of FAM46C knockout mice, but their fertilization ability was severely compromised based on the results of intracytoplasmic sperm injection assays. Interestingly, our RNA-sequencing analyses of FAM46C knockout testes revealed that mRNA levels of only nine genes were significantly altered compared to wild-type ones (q < 0.05). When considering alternate activities for FAM46C, in vitro assays demonstrated that FAM46C does not exhibit protein kinase or AMPylation activity against general substrates. Together, our data show that FAM46C in spermatids is a novel component in fastening the sperm head and flagellum.

MicroRNA-202 maintains spermatogonial stem cells by inhibiting cell cycle regulators and RNA binding proteins.

miRNAs play important roles during mammalian spermatogenesis. However, the function of most miRNAs in spermatogenesis and the underlying mechanisms remain unknown. Here, we report that miR-202 is highly expressed in mouse spermatogonial stem cells (SSCs), and is oppositely regulated by Glial cell-Derived Neurotrophic Factor (GDNF) and retinoic acid (RA), two key factors for SSC self-renewal and differentiation. We used inducible CRISPR-Cas9 to knockout miR-202 in cultured SSCs, and found that the knockout SSCs initiated premature differentiation accompanied by reduced stem cell activity and increased mitosis and apoptosis. Target genes were identified with iTRAQ-based proteomic analysis and RNA sequencing, and are enriched with cell cycle regulators and RNA-binding proteins. Rbfox2 and Cpeb1 were found to be direct targets of miR-202 and Rbfox2 but not Cpeb1, is essential for the differentiation of SSCs into meiotic cells. Accordingly, an SSC fate-regulatory network composed of signaling molecules of GDNF and RA, miR-202 and diverse downstream effectors has been identified.

Diversity analysis of intestinal microflora between healthy and diarrheal neonatal piglets from the same litter in different regions.

Though the problem of neonatal piglet diarrhea is well known, the differences in the bacterial diversity and community composition between healthy and diarrheal piglets are still unknown. We investigated these differences in neonatal piglets from Jiangxi Province, China. Healthy (H, n = 20) and diarrheal (D, n = 20) piglets were selected from six regions. The fecal microbial communities were analyzed by sequencing the V3V4 region of 16S rRNA gene. We found the ratio of major phyla (Fusobacteria, Bacteroidetes, Firmicutes and Proteobacteria) was >99% of 7 phyla. The overall alpha diversity indices, such as chao, sobs, coverage and Shannon, were not significantly different. Moreover, the relative abundance of the predicted functions was highly similar in the two groups. Our results indicated that Clostridium was divided into two major groups: Clostridium sensu stricto_1 and stricto_2. Sensu stricto_2 was highly abundant in the D group and low abundance in the H group, whereas the results of sensu stricto_1 were opposite. Comparative analyses within the H or D groups showed that Escherichia-Shigella and Streptococcus at the genus level and unclassified Lactobacillus at the species level were significant difference. Comparative analyses of the two groups showed that unclassified Prevotellaceae at the genus level and Fusobacterium mortifierum were significantly different and had high linear discriminant analysis (LDA) scores. The significantly different microbes composition results also existed in the same litter, based on excluding regions influence. These results suggested that piglet diarrhea was closely associated with these microbes. This study provides insights into gut microbial interactions and prevention of piglet diarrhea.

Expression dynamics, relationships, and transcriptional regulations of diverse transcripts in mouse spermatogenic cells.

Among all tissues of the metazoa, the transcritpome of testis displays the highest diversity and specificity. However, its composition and dynamics during spermatogenesis have not been fully understood. Here, we have identified 20,639 message RNAs (mRNAs), 7,168 long non-coding RNAs (lncRNAs) and 15,101 circular RNAs (circRNAs) in mouse spermatogenic cells, and found many of them were specifically expressed in testes. lncRNAs are significantly more testis-specific than mRNAs. At all stages, mRNAs are generally more abundant than lncRNAs, and linear transcripts are more abundant than circRNAs. We showed that the productions of circRNAs and piRNAs were highly regulated instead of random processes. Based on the results of a small-scale functional screening experiment using cultured mouse spermatogonial stem cells, many evolutionarily conserved lncRNAs are likely to play roles in spermatogenesis. Typical classes of transcription factor binding sites are enriched in the promoters of testis-specific m/lncRNA genes. Target genes of CREM and RFX2, 2 key TFs for spermatogenesis, were further validated by using ChIP-chip assays and RNA-seq on RFX2-knockout spermatogenic cells. Our results contribute to the current understanding of the transcriptomic complexity of spermatogenic cells and provide a valuable resource from which many candidate genes may be selected for further functional studies.

Genome-wide loss of 5-hmC is a novel epigenetic feature of Huntington's disease.

5-Hydroxymethylcytosine (5-hmC) may represent a new epigenetic modification of cytosine. While the dynamics of 5-hmC during neurodevelopment have recently been reported, little is known about its genomic distribution and function(s) in neurodegenerative diseases such as Huntington's disease (HD). We here observed a marked reduction of the 5-hmC signal in YAC128 (yeast artificial chromosome transgene with 128 CAG repeats) HD mouse brain tissues when compared with age-matched wild-type (WT) mice, suggesting a deficiency of 5-hmC reconstruction in HD brains during postnatal development. Genome-wide distribution analysis of 5-hmC further confirmed the diminishment of the 5-hmC signal in striatum and cortex in YAC128 HD mice. General genomic features of 5-hmC are highly conserved, not being affected by either disease or brain regions. Intriguingly, we have identified disease-specific (YAC128 versus WT) differentially hydroxymethylated regions (DhMRs), and found that acquisition of DhmRs in gene body is a positive epigenetic regulator for gene expression. Ingenuity pathway analysis (IPA) of genotype-specific DhMR-annotated genes revealed that alternation of a number of canonical pathways involving neuronal development/differentiation (Wnt/ß-catenin/Sox pathway, axonal guidance signaling pathway) and neuronal function/survival (glutamate receptor/calcium/CREB, GABA receptor signaling, dopamine-DARPP32 feedback pathway, etc.) could be important for the onset of HD. Our results indicate that loss of the 5-hmC marker is a novel epigenetic feature in HD, and that this aberrant epigenetic regulation may impair the neurogenesis, neuronal function and survival in HD brain. Our study also opens a new avenue for HD treatment; re-establishing the native 5-hmC landscape may have the potential to slow/halt the progression of HD.

Integrative proteomic and transcriptomic analyses reveal multiple post-transcriptional regulatory mechanisms of mouse spermatogenesis.

Mammalian spermatogenesis consists of many cell types and biological processes and serves as an excellent model for studying gene regulation at transcriptional and post-transcriptional levels. Many key proteins, miRNAs, and perhaps piRNAs have been shown to be involved in post-transcriptional regulation of spermatogenesis. However, a systematic method for assessing the relationship between protein and mRNA expression has not been available for studying mechanisms of post-transcriptional regulation. In the present study, we used the iTRAQ-based quantitative proteomic approach to identify 2008 proteins in mouse type A spermatogonia, pachytene spermatocytes, round spermatids, and elongative spermatids with high confidence. Of these proteins, 1194 made up four dynamically changing clusters, which reflect the mitotic amplification, meiosis, and post-meiotic development of germ cells. We identified five major regulatory mechanisms termed "transcript only," "transcript degradation," "translation repression," "translation de-repression," and "protein degradation" based on changes in protein level relative to changes in mRNA level at the mitosis/meiosis transition and the meiosis/post-meiotic development transition. We found that post-transcriptional regulatory mechanisms are related to the generation of piRNAs and antisense transcripts. Our results provide a valuable inventory of proteins produced during mouse spermatogenesis and contribute to elucidating the mechanisms of the post-transcriptional regulation of gene expression in mammalian spermatogenesis.

Confocal Microscopic Analysis of the Spindle and Chromosome Configurations of in vitro-Matured Oocytes from Different Types of Polycystic Ovary Syndrome Patients.

OBJECTIVE: To evaluate the meiotic spindle and chromosome distribution of in vitro-matured oocytes from infertile nonobese or obese women with polycystic ovary syndrome (PCOS) and to compare in vitro maturation (IVM) rates between groups. DESIGN: Prospective study. SETTING: Hospital-based IVF center. PATIENTS: A total of 99 patients (26 obese women with PCOS, 25 nonobese women with PCOS, and 48 controls) undergoing stimulated cycles for intracytoplasmic sperm injection had immature oocytes for IVM. INTERVENTIONS: Immature oocytes (germinal vesicle and metaphase I stages) were collected from obese and nonobese PCOS patients and controlled infertile patients. The meiotic spindle and chromosome configurations in oocytes matured in vitro were studied by confocal microscopy, with fluorescent labeling techniques for visualization of both microtubules and chromosomes. MAIN OUTCOME MEASURES: Meiotic spindle and associated chromosome configurations. RESULTS: There were no significant differences between the different types of PCOS and the control group with respect to IVM rates (61.8, 63.8, and 63.2%, respectively), the percentage of spindle abnormalities in metaphase II oocytes (40.6, 42.9, and 37.5%, respectively) or chromosome abnormalities in metaphase II oocytes (31.2, 34.3, and 33.3%, respectively). CONCLUSIONS: In vitro-matured oocytes obtained from stimulated cycles had a high ratio of meiotic abnormalities. The different types of PCOS had the same ratio of meiotic abnormalities.

Transcription factor PBX4 regulates limb development and haematopoiesis in mice.

The mammalian Pre-B cell leukaemia transcription factors 1-4 (PBX1-4) constitutes the PBC class of the homeodomain (HD)-containing proteins, which play important roles in diverse developmental processes. The functions and the underlying molecular mechanisms of PBX1-3 but not PBX4 have been extensively studied, and they have been reported to direct essential morphogenetic processes and organogenesis. In the present study, we generated knockin mice of FLAG-tagged PBX4 and the Pbx4 knockout (KO) mice and carried out in-depth characterisation of PBX4 expression and function. PBX4 was initially detected only in the testis among several organs of the adult mice and was expressed in spermatocytes and spermatids. However, no abnormality in spermatogenesis, but growth retardation and premature death after birth were observed in most adult Pbx4 KO mice. These animals were inactive and had shorter hindlimbs and lower numbers of reticulocytes and lymphocytes, probably caused by abnormalities at earlier developmental stages. Pbx4 mRNAs were indeed detected in several embryonic cell types related to limb development by in situ hybridisation and single-cell RNA-sequencing analysis. Pbx4 protein was also detected in the bone marrow of adult mice with a lower level compared with that in the testis. PBX4 preferentially binds to the promoters of a large number of genes including those for other HD-containing proteins and ribosomal proteins whose mutations are related to anaemia. PBX4-binding sites are enriched in motifs similar to those of other HD-containing proteins such as PKNOX1 indicating that PBX4 may also act as a co-transcription factor like other PBC proteins. Together, these results show that PBX4 participates in limb development and haematopoiesis while its function in spermatogenesis has not been revealed by gene KO probably due to the complementary effects of other genes.

Identification of most representative hub-genes for diagnosis, prognosis, and therapies of hepatocellular carcinoma.

BACKGROUND: Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related deaths globally. To reduce HCC-related mortality, early diagnosis and therapeutic improvement are essential. Hub differentially expressed genes (HubGs) may serve as potential diagnostic and prognostic biomarkers, also offering therapeutic targets for precise therapies. Therefore, we aimed to identify top-ranked hub genes for the diagnosis, prognosis, and therapy of HCC. METHODS: Through a systematic literature review, 202 HCC-related HubGs were derived from 59 studies, yet consistent detection across these was lacking. Then, we identified top-ranked HubGs (tHubGs) by integrated bioinformatics analysis, highlighting their functions, pathways, and regulators that might be more representative of the diagnosis, prognosis, and therapies of HCC. RESULTS: In this study, eight HubGs (CDK1, AURKA, CDC20, CCNB2, TOP2A, PLK1, BUB1B, and BIRC5) were identified as the tHubGs through the protein-protein interaction (PPI) network and survival analysis. Their differential expression in different stages of HCC, validated using The Cancer Genome Atlas (TCGA) Program database, suggests their potential as early HCC markers. The enrichment analyses revealed some important roles in HCC-related biological processes (BPs), molecular functions (MFs), cellular components (CCs), and signaling pathways. Moreover, the gene regulatory network analysis highlighted key transcription factors (TFs) and microRNAs (miRNAs) that regulate these tHubGs at transcriptional and post-transcriptional. Finally, we selected three drugs (CD437, avrainvillamide, and LRRK2-IN-1) as candidate drugs for HCC treatment as they showed strong binding with all of our proposed and published protein receptors. CONCLUSIONS: The findings of this study may provide valuable resources for early diagnosis, prognosis, and therapies for HCC.

VEGF-FGF Signaling Activates Quiescent CD63+ Liver Stem Cells to Proliferate and Differentiate.

Understanding the liver stem cells (LSCs) holds great promise for new insights into liver diseases and liver regeneration. However, the heterogenicity and plasticity of liver cells have made it controversial. Here, by employing single-cell RNA-sequencing technology, transcriptome features of Krt19+ bile duct lineage cells isolated from Krt19CreERT; Rosa26R-GFP reporter mouse livers are examined. Distinct biliary epithelial cells which include adult LSCs, as well as their downstream hepatocytes and cholangiocytes are identified. Importantly, a novel cell surface LSCs marker, CD63, as well as CD56, which distinguished active and quiescent LSCs are discovered. Cell expansion and bi-potential differentiation in culture demonstrate the stemness ability of CD63+ cells in vitro. Transplantation and lineage tracing of CD63+ cells confirm their contribution to liver cell mass in vivo upon injury. Moreover, CD63+CD56+ cells are proved to be activated LSCs with vigorous proliferation ability. Further studies confirm that CD63+CD56- quiescent LSCs express VEGFR2 and FGFR1, and they can be activated to proliferation and differentiation through combination of growth factors: VEGF-A and bFGF. These findings define an authentic adult liver stem cells compartment, make a further understanding of fate regulation on LSCs, and highlight its contribution to liver during pathophysiologic processes.

piRNA profiling during specific stages of mouse spermatogenesis.

PIWI-interacting RNAs (piRNAs) are a class of small RNAs abundantly expressed in animal gonads. piRNAs that map to retrotransposons are generated by a "ping-pong" amplification loop to suppress the activity of retrotransposons. However, the biogenesis and function of other categories of piRNAs have yet to be investigated. In this study, we first profiled the expression of small RNAs in type A spermatogonia, pachytene spermatocytes, and round spermatids by deep sequencing. We then focused on the computational analysis of the potential piRNAs generated in the present study as well as other published sets. piRNAs mapping to retrotransposons, mRNAs, and intergenic regions had different length distributions and were differentially regulated in spermatogenesis. piRNA-generating mRNAs (PRMRs), whose expression positively correlated with their piRNA products, constituted one-third of the protein-coding genes and were evolutionarily conserved and enriched with splicing isoforms and antisense transcripts. PRMRs with piRNAs preferentially mapped to CDSs and 3' UTRs partitioned into three clusters differentially expressed during spermatogenesis and enriched with unique sets of functional annotation terms related to housekeeping activities as well as spermatogenesis-specific processes. Intergenic piRNAs were divided into 2992 clusters probably representing novel transcriptional units that have not been reported. The transcripts of a large number of genes involved in spermatogenesis are the precursors of piRNAs, and these genes are intricately regulated by alternative splicing and antisense transcripts. piRNAs, whose regulatory role in gene expression awaits to be identified, are clearly products of a novel regulatory process that needs to be defined.