Structural and Functional Characterization of a Testicular Long Non-coding RNA (4930463O16Rik) Identified in the Meiotic Arrest of the Mouse Topaz1–/– Testes

Spermatogenesis involves coordinated processes, including meiosis, to produce functional gametes. We previously reported Topaz1 as a germ cell-specific gene highly conserved in vertebrates. Topaz1 knockout males are sterile with testes that lack haploid germ cells because of meiotic arrest after prophase I. To better characterize Topaz1^–/– testes, we used RNA-sequencing analyses at two different developmental stages (P16 and P18). The absence of TOPAZ1 disturbed the expression of genes involved in microtubule and/or cilium mobility, biological processes required for spermatogenesis. Moreover, a quarter of P18 dysregulated genes are long non-coding RNAs (lncRNAs), and three of them are testis-specific and located in spermatocytes, their expression starting between P11 and P15. The suppression of one of them, 4939463O16Rik, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-4939463O16Rik^–/– testes was altered and the molecular pathways affected included microtubule-based processes, the regulation of cilium movement and spermatogenesis. The absence of TOPAZ1 protein or 4930463O16Rik produced the same enrichment clusters in mutant testes despite a contrasted phenotype on male fertility. In conclusion, although Topaz1 is essential for the meiosis in male germ cells and regulate the expression of numerous lncRNAs, these studies have identified a Topaz1 regulated lncRNA (4930463O16Rik) that is key for both sperm production and motility.

Introduction

In mammals, an organism derives from two parental haploid gametes, a maternal oocyte and paternal sperm. Meiosis is a highly specialized event that leads to the production of these haploid germ cells (Kleckner, 1996). In females, meiosis is initiated during fetal life while male germ cells are involved in the meiosis process around puberty. In males, meiosis is essential during spermatogenesis that involves mitotic division and the multiplication of spermatogonia, the segregation of homologous chromosomes and the spermiogenesis of haploid germ cells. This complex process of spermatogenesis, which progresses through precisely timed and highly organized cycles, is primordial for male fertility. All these different events are highly regulated and associated with the controlled expression of several testis-enriched genes. A previous study had demonstrated the essential role of Topaz1 during meiosis in male mice (Luangpraseuth-Prosper et al., 2015). Topaz1 is a highly conserved gene in vertebrates (Baillet et al., 2011). Its expression is germ cell-specific in mice (Baillet et al., 2011). The suppression of Topaz1 in mice (Topaz1^–/–) results in azoospermia (Luangpraseuth-Prosper et al., 2015). Male meiotic blockage occurs without a deregulation of chromosome alignment and TOPAZ1 is not involved in formation of the XY body or the maintenance of MSCI (Meiotic Sex Chromosome Inactivation). Topaz1 depletion increases the apoptosis of mouse adult male pachytene cells and triggers chromosome misalignment at the metaphase I plate in mouse testes (Luangpraseuth-Prosper et al., 2015). This misalignment leads to an arrest at the prophase to metaphase transition during the first meiosis division (Luangpraseuth-Prosper et al., 2015). Microarray-based gene expression profiling of Topaz1^–/– mouse testes revealed that TOPAZ1 influences the expression of one hundred transcripts, including several long non-coding RNAs (lncRNAs) and unknown genes, at postnatal day 20 (P20) (Luangpraseuth-Prosper et al., 2015).

Since discovery of the maternal H19 lncRNA (Brannan et al., 1990) and the Xist (Brockdorff et al., 1992) genes that regulate the structure of chromosomes and mediate gene repression during X chromosome inactivation, interest in studying the role of non-coding RNAs (ncRNAs) has grown considerably. Non-coding RNAs are present in many organisms, from bacteria to humans, where only 1.2% of the human genome codes for functional proteins (Carninci et al., 2005; ENCODE Project Consortium, Birney et al., 2007; Gil and Latorre, 2012). While much remains to be discovered about the functions of ncRNAs and their molecular interactions, accumulated evidence suggests that ncRNAs participate in various biological processes such as cell differentiation, development, proliferation, apoptosis, and cancers.

They are divided into two groups: small and long non-coding RNAs (sncRNAs and lncRNAs, respectively). sncRNAs contain transcripts smaller than 200 nucleotides (nt) (Guttman et al., 2009). They include microRNAs (miRNAs, 20–25 nt), small interfering RNAs (siRNAs), PIWI-interacting RNAs (piRNAs, 26–31 nt) and circular RNAs (cricRNA). These sncRNAs are essential for several functions such as the regulation of gene expression and genome protection (Ref in Morris and Mattick, 2014) as well as during mammalian spermatogenesis (Yadav and Kotaja, 2014; Bie et al., 2018; Quan and Li, 2018). The second group, lncRNAs, contains transcripts longer than 200 nt without a significant open reading frame. Advances in high-throughput sequencing have enabled the identification of new transcripts, including lncRNAs, most of which are transcribed by RNA polymerase II and possess a 5′ cap and polyadenylated tail (ref in Jarroux et al., 2017). They are classified according to their length, location in the genome (e.g., surrounding regulatory elements) or functions.

Several studies have pointed out that testes contain a very high proportion of lncRNA compared to other organs (Necsulea and Kaessmann, 2014; Sarropoulos et al., 2019). However, this high testicular expression is only observed in the adult organ, as the level of lncRNAs in the developing testis is comparable to that seen in somatic organs (Sarropoulos et al., 2019). In mice, some testis-expressed lncRNAs were functionally characterized during spermatogenesis. Thus, the lncRNA Mrhl repressed Wnt signaling in the GC1-Spg spermatogonial cell line, suggesting a role in spermatocyte differentiation (Arun et al., 2012). Expression of the Testis-specific X-linked gene was specific to, and highly abundant in, mouse pachytene-stage spermatocytes and could regulate germ cell progression in meiosis (Anguera et al., 2011). Moreover, in male germ cells, it has been shown that the Dmrt1-related gene negatively regulates Dmrt1 (doublesex and mab-3 related transcription factor 1) and that this regulation might be involved in the switch between mitosis and meiosis in spermatogenesis (Zhang et al., 2010). Lastly, in a non-mammalian model as Drosophila, Wen et al. (2016) produced mutant fly lines by deleting 105 testis-specific lncRNAs and demonstrated the essential role of 33 of them in spermatogenesis and/or male fertility.

In a previous study, we performed comparative microarray analyses of wild-type and Topaz1^–/– testis RNAs at P15 and P20 (Luangpraseuth-Prosper et al., 2015). Here, we have performed bulk RNA-sequencing (RNA-seq) of testes collected at P16 and P18 in order to refine the developmental stages that display transcriptional differences between the two mouse lines and to obtain an exhaustive list of dysregulated transcripts, especially non-coding ones. Since the proportion of deregulated lncRNAs represented about a quarter of the differentially expressed genes (DEGs), we studied the testicular localization of three of them. In order to approach the role of testicular lncRNAs, we created a mouse line in which one of them was deleted (4930463O16Rik). These knockout mice displayed normal fertility in both sexes, but the male mutants produced half as much sperm as wild-type controls.

Materials and Methods

Ethics Statement

All animal experiments were performed in strict accordance with the guidelines of the Code for Methods and Welfare Considerations in Behavioral Research with Animals (Directive 2016/63/UE). All experiments were approved by the INRAE Ethical Committee for Animal Experimentation covering Jouy-en-Josas (COMETHEA, no. 18-12) and authorized by the French Ministry for Higher Education, Research and Innovation (No. 815-2015073014516635).

Mice

The generation and preliminary analysis of Topaz1-null transgenic mouse line has been described previously (Luangpraseuth-Prosper et al., 2015).

Generation of the 4630493O16Rik-null transgenic mouse line was achieved using CrispR-Cas9 genome editing technology. The RNA mix comprised an mRNA encoding for SpCas9-HF1 nuclease and the four sgRNA (Supplementary Table 1) targeting the 4930463o16Rik gene (NC_000076: 84324157-84333540). These sgRNAs were chosen according to CRISPOR software¹ in order to remove the four exons and introns of the 4930463o16Rik gene. Cas9-encoding mRNA and the four sgRNAs were injected at a rate of 100 ng/μL each into one cell fertilized C57Bl/6N mouse eggs (Henao-Mejia et al., 2016). The surviving injected eggs were transferred into pseudo-pregnant recipient mice. Tail-DNA analysis of the resulting live pups was performed using PCR with genotyping oligonucleotides (Supplementary Table 1) and the Takara Ex Taq^® DNA Polymerase kit. The PCR conditions were 94°C 30 s, 60°C 30 s, and 72°C 30 s, with 35 amplification cycles. Two transgenic founder mice were then crossed with wild-type C57Bl/6N mice to establish transgenic lines. F1 heterozygote mice were crossed together in each line to obtain F2 homozygote mice, thus establishing the 4630493O16Rik^–/– mouse lines. Both mouse lines were fertile and the number of pups was equivalent, so we worked with one mouse line.

All mice were fed ad libitum and were housed at a temperature of 25°C under a 12 h/12 h light/dark cycle at the UE0907 unit (INRAE, Jouy-en-Josas, France). The animals were placed in an enriched environment in order to improve their receptiveness while respecting the 3R. All mice were then sacrificed by cervical dislocation. Tissues at different developmental stages were dissected and fixed as indicated below, or flash frozen immediately in liquid nitrogen before storage at -80°C. The frozen tissues were used for the molecular biology experiments described below.

Histological and Immunohistochemical Analyses

For histological studies, fresh tissues from 8-week-old mice were fixed in 4% paraformaldehyde (Electron Microscopy Sciences reference 50-980-495) in phosphate buffer saline (PBS) at 4°C. After rinsing the tissues in PBS, they were stored in 70% ethanol at 4°C. Paraffin inclusions were then performed using a Citadel automat (Thermo Scientific Shandon Citadel 1000) according to a standard protocol. Tissues included in paraffin blocks were sectioned at 4 μm and organized on Superfrost Plus Slides (reference J1800AMNZ). Once dry, the slides were stored at 4°C. On the day of the experiment, these slides of sectioned tissues were deparaffinized and rehydrated in successive baths of xylene and ethanol at room temperature. For histology, testes sections were stained with hematoxylin and eosin (HE) by the @Bridge platform (INRAE, Jouy-en-Josas) using an automatic Varistain Slide Stainer (Thermo Fisher Scientific). Periodic acid-Schiff staining (PAS) was used to determine seminiferous epithelium stages.

In situ hybridization experiments were performed using the RNAscope^® system (ACB, Bio-Techne SAS, Rennes, France). Briefly, probes (around 1000 nt long) for Topaz1 (NM_001199736.1), 4930463o16Rik (NR_108059.1), Gm21269 (NR_102375.1), and 4921513H07Rik (NR_153846.1) were designed by ACB and referenced with the catalog numbers 402321, 431411, 549421 and 549441, respectively. Negative (dapB, Bacillus subtilis dihydrodipicolinate reductase) and positive (PPIB, Mus musculus peptidylprolyl isomerase B, and UBC, Homo sapiens ubiquitin C) controls were ordered from ACD. Hybridization was performed according to the manufacturer’s instructions using a labeling kit (RNAscope^® 2.5HD assay-brown). Slides were counterstained according to a PAS staining protocol and then observed for visible signals. Hybridization was considered to be positive when at least one dot was observed in a cell. Stained sections were scanned using a 3DHISTECH panoramic scanner at the @Bridge platform (INRAE, Jouy-en-Josas) and analyzed with Case Viewer software (3DHISTECH). In situ hybridization with the control probes was performed using tissue sections of 2-month-old WT testes (Supplementary Figure 1). We also used the RNAscope^® 2.5HD assay-red kit in combination with immunofluorescence in order to achieve the simultaneous visualization of RNA and protein on the same slide. The ISH protocol was thus stopped by immersion in water before hematoxylin counterstaining. Instead, the slides were washed in PBS at room temperature. The Mouse on mouse (M.O.M.) kit (BMK-2202, Vector laboratories) was used and slides were incubated for 1 h in Blocking Reagent, 5 min in Working solution and 2 h with a primary antibody: DDX4 (ab13840, Abcam) or γH2AX(Ser139) (Merck), diluted at 1:200 in Blocking Reagent. Detection was ensured using secondary antibody conjugated to DyLight 488 (green, KPL). Diluted DAPI (1:1000 in PBS) was then applied to the slides for 8 min. The slides were then mounted with Vectashield Hard Set Mounting Medium for fluorescence H-1400 and images were captured at the MIMA2 platform²,³ using an inverted ZEISS AxioObserver Z1 microscope equipped with an ApoTome slider, a Colibri light source and Axiocam MRm camera. Images were analyzed using Axiovision software 4.8.2 (Carl Zeiss, Germany).

Total RNA Extraction and Quantitative RT-PCR (RT-qPCR)

Total RNAs from post-natal mouse testes or other organs were isolated using Trizol reagent. The RNAs were purified using the RNeasy Mini kit (Qiagen) following the manufacturer’s instructions and then DNAse-treated (Qiagen). The quantification of total RNAs was achieved with a Qbit^® Fluorometric Quantitation. Maxima First-Strand cDNA Synthesis Kit (Thermo Scientific) was used to reverse transcribe RNA into cDNA. The Step One system with Fast SYBR^TM Green Master Mix (Applied Biosystems, ThermoFisher France) was used for qPCR, which was performed in duplicate for all tested genes and the results were normalized with qBase⁺ software (Biogazelle) (Hellemans et al., 2007). Gapdh, Ywhaz, and Mapk1 were used as the reference genes. For each experiment, median values were plotted using GraphPad Prism, and statistical analyses were performed with KrusKall–Wallis tests under R software [Rcmdr package (p-value < 0.05)]. The primer sequences used for RT-qPCR are shown in Supplementary Table 1.

RNA-Sequencing

Total RNA was extracted from WT, Topaz1^–/– or 1630463O16Rik^–/– mouse testes at P16 and P18 (n = 3 for each mouse phenotype and each developmental stage). Total RNA quality was verified on an Agilent 2100 Bioanalyser (Matriks, Norway) and samples with a RIN > 9 were made available for RNA-sequencing. This work benefited from the facilities and expertise of the I2BC High-throughput Sequencing Platform (Gif-sur-Yvette, Université Paris-Saclay, France) for oriented library preparation (Illumina Truseq RNA Sample Preparation Kit) and sequencing (Paired-end 75 bp; NextSeq). More than 38 million 75 bp paired-end reads per sample were generated. Demultiplexing was done (bcl2fastq2-2.18.12) and adapters were removed (Cutadapt1.15) at the I2BC High-throughput Sequencing Platform. Only reads longer than 10 pb were used for analysis. Quality control of raw RNA-Seq data were processed by FastQC v0.11.5.

Transcriptomic Analysis

Sequence libraries were aligned with the Ensembl 95 genome using TopHat (Trapnell et al., 2009), and gene table counts were obtained by applying featureCounts to these alignments (Liao et al., 2014). Data normalization and single-gene level analyses of differential expression were performed using DESeq2 (Love et al., 2014). Since some samples were sequenced several months apart, a batch effect after computation of hierarchical clustering was observed. In order to take this effect into account, we introduced the batch number into the DESeq2 model, as well as the study conditions. Differences were considered to be significant for Benjamini–Hochberg adjusted p-values < 0.05, and absolute fold changes > 2 (absolute Log2FC > 1) (Benjamini and Hochberg, 1995). Raw RNA-seq data were deposited via the SRA Submission portal⁴, BioProject ID PRJNA698440.

Biotype Determination of DEGs

Data available on the NCBI, MGI⁵ and Ensembl⁶ websites were used simultaneously to determine the DEG biotypes. For this purpose, information on the mouse genome was obtained by ftp from NCBI⁷; the annotation BioMart file from Ensembl⁸ (Ensembl genes 95, Mouse genes GRCm28.p6) and feature types from MGI⁹ (with the protein coding gene, non-coding RNA gene, unclassified gene and pseudogenic region). Only data corresponding to the DEGs were conserved. The files from these three databases were therefore cross-referenced to determine DEG biotypes. When the biotype of a gene differed between databases, the annotation was then listed as genes with a “biotype conflict.”

Gene Ontology Enrichment

The mouse DEGS thus identified were analyzed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway membership with Database performed using the DAVID Bioinformatic Database 6.8¹⁰. These analyses and pathways were considered to be significant for a Benjamini-corrected enrichment p-value of less than 0.05. The Mouse Atlas Genome of differentially expressed genes extracted from this study was performed via the Enrichr website¹¹.

Sperm Analysis

Evaluations of the concentrations and motility of sperm in WT and 4930463O16Rik^–/– 8-week-old mice were performed using the IVOS II Computer Assisted Sperm Analysis (CASA) system (Hamilton Thorne, Beverly, MA, United States). The two fresh cauda epididymes from each individual were removed and plunged into 200 μL TCF buffer (Tris, citrate and fructose buffer) where they were chopped up with small scissors. For sperm release, the samples were incubated for 10 min at 37°C. A 4 μl aliquot was placed in a standardized four-chamber Leja counting slide (Leja Products B.V., Nieuw-Vennep, Netherlands). Ten microscope fields were analyzed using the predetermined starting position within each chamber with an automated stage. Statistical analyses were performed using the mean of the 10 analyzed fields containing at least 300 cells. The IVOS settings chosen were those defined for mouse sperm cell analysis (by Hamilton Thorne). The principal parameters were fixed as follows: 45 frames were captured at 60 Hz; for cell detection, the camera considered a signal as a spermatozoon when the elongation percentage was between 70 (maximum) and 2 (minimum); the minimal brightness of the head at 186, and the minimum and maximum size of the head at 7 and 100 μm², respectively. The kinematic thresholds applied were: cell travel max at 10 μm, progressive STR at 45%, progressive VAP at 45 μm/s, slow VAP at 20 μm/s, slow VSL at 30 μm/s, static VAP at 4 μm/S and static VSL at 1 μm/s. The full settings used are listed in Supplementary Table 2. The CASA parameters thus recorded included the average path velocity (VAP in μm/s), straight line velocity (VSL in μm/s), curvilinear velocity (VCL in μm/s), amplitude of lateral head displacement (ALH in μm), motility (percentage), and sperm concentration (0.10⁶/mL). Slow cells were recorded as static. Median and interquartile ranges were plotted with GraphPad. To compare the sperm parameters between WT and 4930463O16Rik^–/– mice, statistical analyses were performed using the Kruskal–Wallis non-parametric test.

Results

Topaz1 Mutant Testes Experience a Deregulated Transcriptome as Early as P16

To expand on the previous comparative microarray analyses of wild-type and mutant testes performed at P15 and P20 during the first wave of spermatogenesis (Luangpraseuth-Prosper et al., 2015), transcriptomic analyses by RNA-seq were performed on WT and Topaz1^–/– mouse testes at two developmental stages: P16 and P18. The P15 stage previously chosen for microarray analyses was too early, with only one gene differentially expressed (Topaz1). At P15, seminiferous tubules contain spermatocytes that have advanced to mid and late-pachytene. Therefore, we chose a one day later, P16, development stage to perform RNA-seq analysis. At P16, seminiferous tubules contain spermatocyte cells that have progressed from the end-pachytene to early diplotene of meiosis I. The other developmental stage for the RNA-seq analysis was P18. This is just after the first meiosis I division. At P18, late-pachytene spermatocytes are abundant and the very first spermatocytes II appear (Drumond et al., 2011). The P20 stage chosen for the previous microarray analyses was, therefore, too late since the first round spermatids had already appeared (Supplementary Figure 2). Therefore, the P16 and P18 stages chosen for this study include the just before and just after the first meiosis I division of spermatogenesis. Only 10% of the differential genes found by microarray at P20 correspond to lncRNAs, which may be underrepresented on microarray. Therefore, we re-analyzed the transcriptome of Topaz1 KO testes by RNA-sequencing technology in order to obtain an exhaustive list of dysregulated transcripts, especially non-coding transcripts.

Differential analyses of RNA-seq results revealed that 205 and 2748 genes were significantly deregulated in Topaz1^–/– testes compared to WT at P16 and P18, respectively [adjusted p-value (Benjamini–Hochberg) < 0.05 and absolute Log2 Fold Change > 1 (Log2FC > 1)] (Figure 1A and Supplementary Table 3). At P16, out of the 205 DEGs, 97 genes were significantly down-regulated (Log2FC < −1 or FC < 0.5) and 108 were up-regulated (Log2FC > 1 or FC > 2). However, at P18, down-regulated DEGs accounted for 91% (2491 genes) and up-regulated genes for only 9% (257 genes). Among all these DEGs, 120 were common to both P16 and P18 (Figure 1A). According to the mouse gene atlas, the 2748 DEGs at developmental stage P18 were largely testis-enriched DEGs in mouse testis-specific genes (Figure 1B).

FIGURE 1

Figure 1. WT vs. Topaz1^–/– deregulated gene analysis from mouse testes. (A) Venn diagram showing the overlap of differentially expressed genes between P16 and P18 Topaz1^–/– mouse testes [adjusted p < 0.05 and down-regulated FC < 0.5 (log2FC < −1) or up-regulated FC > 2 (log2FC > 1)]. (B) Clustergrammer was generated by the Enrichr website. Top 10 enriched tissues are the columns, input genes (2748 DEGs of P18 Topaz1^–/– compared to normal testes) are the rows, and cells in the matrix indicate whether a DEG is associated with a tissue from the Mouse Gene Atlas. (C,D) Biotype of DEGs in Topaz1^–/– testes at (C) P16 and (D) P18. Around half of them are protein-coding genes whereas around one quarter is ncRNA at both developmental stages. (E,F) Biotype of DEGs in Topaz1^–/– testes at (E) P16 and (F) P18, depending on whether they were up- or down-regulated.

The validation of several DEGs was achieved using RT-qPCR. Two randomly selected genes up-regulated at P16 (B3galt2 and Hp) and three at P18 (B3galt2, Afm, and Cx3cr1), four genes down-regulated at P16 and P18 (Gstt2, 4930463O16Rik, 4921513H07Rik, and Gm21269) and two non-differential genes (Cdc25c and Nop10) were analyzed (Supplementary Figure 3). The results confirmed those obtained using RNA-seq.

The biotypes of the differential transcripts (protein-coding, non-coding RNAs, etc.) were determined from the annotations of the NCBI, MGI and Ensembl databases. Two major deregulated groups were highlighted at both stages. The protein-coding gene biotype accounted for half of the deregulated genes (51.2 and 57.2% at P16 and P18, respectively) (Figures 1C,D). A quarter of Topaz1^–/– DEGs; 24.9 and 27.9% at P16 and P18, respectively, was found to belong to the second ncRNA group. Among the latter, the major biotype was lncRNAs at both stages, being 23.4 and 27.1% at P16 and P18, respectively. This significant proportion of deregulated lncRNA thus raised the question of their potential involvement in spermatogenesis.

Pathway and Functional Analysis of DEGs

To further understand the biological functions and pathways involved, these DEGs were functionally annotated based on GO terms and KEGG pathway or on InterPro databases through the ontological Database for Annotation, Visualization and Integrated Discovery (DAVID v6.8, see footnote 10) using the default criteria (Huang et al., 2009a,b).

At P16, so therefore before the first meiosis division, out of 205 differentially expressed genes, 32% of down-regulated and 20% of up-regulated genes corresponded to non-coding RNAs with no GO annotation or no pathway affiliation for the vast majority (Figure 1E), leading to less powerful functional annotation clustering (Supplementary Table 4). Five clusters with an enrichment score > 1.3 were obtained (an enrichment score > 1.3 was used for a cluster to be statistically significant, as recommended by Huang et al. (2009a) but the number of genes in each cluster was small except for annotation cluster number 4. In this, an absence of TOPAZ1 appeared to affect the extracellular compartment. The others referred to the antioxidant molecular function and the biological detoxification process, suggesting stressful conditions.

At P18, corresponding to the first transitions from prophase to metaphase, and considering either all DEGs (2748 DEGs; 2404 DAVID IDs) or only down-regulated genes (2491 DEGs; 2164 DAVID IDs) in the P18 Topaz1^–/– versus WT testes, 33 and 23 clusters with an enrichment score > 1.3 were obtained respectively. Some clusters had a strong enrichment score, it was possible to identify five identical clusters with an enrichment score higher than 12 (Figure 1F and Supplementary Table 5). However, the enrichment scores were higher when only down-regulated genes were considered. These clusters include the following GO terms: (i) for cellular components: motile cilium, ciliary part, sperm flagellum, axoneme, acrosomal vesicule; (ii) for biological processes: microtubule-based process, spermatogenesis, germ cell development, spermatid differentiation (Supplementary Table 5).

Finally, using the InterPro database, four clusters with enrichment scores > 1.3 were obtained based on down-regulated genes (Supplementary Table 5) and with up-regulated genes, an absence of TOPAZ1 from mouse testes highlighted the biological pathway of the response to external stimulus or the defense response in the testes. Once again, and as for P16, this suggested stressful conditions in these Topaz1^–/– testes.

These results indicate that an absence of TOPAZ1 induced alterations to the murine transcriptome of the mutant testis transcriptome as early as 16 days after birth. Two days later (P18), these effects were amplified and predominantly involved a down-regulation of genes (91% of DEGs). The loss of TOPAZ1 appeared to disrupt the regulation of genes involved in microtubule and/or cilium mobility, spermatogenesis and first meiotic division during the prophase to metaphase transition. This was in agreement with the Topaz1^–/– phenotype in testes.

Selection of Three Deregulated lncRNA, All Spermatocyte-Specific

The vast majority of deregulated lncRNAs in Topaz1^–/– testes has an unknown function. We decided to study three of the 35 down-regulated lncRNAs that are shared at the P16 and P18 stages. The two first, 4930463O16Rik (ENSMUSG00000020033), that is the most down-regulated gene at P16 with a Log2FC of 11.85, and 4921513H07Rik (ENSMUSG00000107042), were both already highlighted by the previous microarray comparative analyses (Luangpraseuth-Prosper et al., 2015). The third one, Gm21269 (ENSMUSG00000108448), has the lowest adjusted p-value at P18.

We quantified these transcripts by qPCR in several somatic tissues (brain, heart, liver, lung, small intestine, muscle, spleen, kidney, epididymis, and placenta) and in the gonads (testes and ovary). These three lncRNAs were almost exclusively expressed in testes (Figures 2A,C,E). These results were in agreement with RNA-seq data available for 4930463O16Rik and Gm21269 on the ReproGenomics viewer¹² (Supplementary Figures 4A,5A, respectively) (Darde et al., 2015, 2019). Our RNA-seq results, summarized using our read density data (bigwig) and the Integrative Genomics Viewer (IGV;¹³), revealed little or no expression of these three genes in Topaz1^–/– testes (Supplementary Figures 6A–C).

FIGURE 2

Figure 2. Expression analysis of three lncRNAs. (A,C,E) RT-qPCR analysis of three different lncRNAs. (A) 4930463O16Rik; (C) Gm21269; (E) 4921513H07Rik in different 2-month-old tissues from WT mice. The red lines represent the median for each tissue; n = 5 for testes and n = 3 for other organs. Statistical analyses were performed using the non-parametric Kruskal–Wallis test. *p < 0.05 (B,D,F) Schematic representation of the results of (B) 4930463O16Rik, (D) Gm21269, and (F) 49215113H07Rik ISH expression in meiotic and post-meiotic cells of the WT mouse seminiferous epithelial cycle.

Quantification of these transcripts using qPCR from postnatal to adulthood in WT and Topaz1^–/– testes had previously been reported, as for 4930463O16Rik and 4921513H07Rik (Figure 9 in Luangpraseuth-Prosper et al., 2015) or performed for Gm21269 (Supplementary Figure 7, also including the postnatal expression of 4930463O16Rik and 4921513H07Rik). The difference in expression between normal and Topaz1^–/– testes was detected as being significant as early as P15 (detected as insignificant in the previous microarray analysis and Gm21269 was absent from the microarray employed). All showed an absence of expression, or at least an important down-regulation, in mutant testes.

To determine the testicular localization of these lncRNA, in situ hybridization (ISH) on adult WT testes sections was performed (Supplementary Figure 8) and the results summarized (Figures 2B,D,F). These lncRNAs were expressed in spermatocytes and the most intense probe labeling was observed at the pachytene stage. These results were confirmed by data on the ReproGenomics viewer for 4930463O16Rik and Gm21269 (see footnote 12) (Supplementary Figures 4B,5B; Darde et al., 2015, 2019).

To refine the subcellular localization of these transcripts in adult mouse testes, we paired ISH experiments and the IF staining of DDX4 protein (or MVH, Mouse Vasa homolog). DDX4 is a germ cell cytoplasmic marker of germ cells, especially in the testes (Toyooka et al., 2000). Our results showed that the three lncRNAs observed displayed different intensities of expression depending on seminiferous epithelium stages. 4930463O16Rik was expressed in the nucleus of spermatocytes with diffuse fluorescence, surrounded by cytoplasmic DDX4 labeling from the zygotene to the diplotene stages (Figures 3A–C). At the same spermatocyte stages (zygotene to diplotene), a diffuse labeling of Gm21269, similar to that of 4930463O16Rik, was observed but with the addition of dot-shaped labeling that co-localized with DDX4 fluorescence (Figures 3D–F). Gm21269 was therefore localized in the cytoplasm and nuclei of spermatocytes during meiosis. 4921513H07Rik appeared to be cytoplasmic, with fluorescent red dots (ISH) surrounding the nuclei, and located in close proximity to DDX4 (IF) labeling (Figures 3G–I). At other stages, identified by DDX4 staining, ISH labeling of these three lncRNA revealed single dots in a few spermatogonia and in round spermatids. The same experiment was then repeated: ISH was followed by IF staining of γH2Ax to highlight the sex body in spermatocytes (Supplementary Figure 9). No co-localization between the sex body and the three lncRNA was revealed.

FIGURE 3

Figure 3. lncRNA cellular localizations on WT 2-month-old mouse testes. In situ hybridization using (A) 4930463O16Rik, (D) Gm21269, and (G) 4921513H07Rik probes (red). (B,E,H) Immunofluorescence staining with DDX4 antibody was achieved at the same stage of seminiferous epithelium to identify male germ cells (green). (C,F,I) DAPI (blue), visualizing nuclear chromosomes, was merged with ISH (green) and IF (red) signals. Zooms in white squares show spermatocytes during the first meiotic division (zygotene to diplotene stages). Zooms in circles show spermatid cells with one spot of DDX4 staining per cell. Scale bar = 20 μm.

Taken together, these results indicate that these spermatocyte-specific lncRNAs had different subcellular localizations in spermatocytes, suggesting functions in these male germ cells.

Generation of 4930463O16Rik-Deleted Mice

In order to evaluate a potential role in spermatogenesis for the lncRNA with the highest Log2FC at P16, namely 4930463O16Rik (the nuclear expressed gene revealed by HIS), a mouse knockout model was established.

The 4930463O16Rik gene (Chr10: 84,488,293-84,497,435 – GRCm38:CM001003.2) is described in public databases as consisting of four exons spanning approximately 10 kb in an intergenic locus on mouse chromosome 10. Using PCR and sequencing, we confirmed this arrangement (data not shown). In order to understand the role of 4930463O16Rik, a new mouse line depleted of this lncRNA was created using CRISPR/Cas9 technology (Figure 4A). Briefly, multiple single guide RNAs (sgRNAs) were chosen, two sgRNAs in 5′ of exon 1 and two sgRNAs in 3′ of exon 4, so as to target the entire length of this gene (Figures 4A,C) and enhance the efficiency of gene deletion in the mouse (Han et al., 2014). Mice experiencing disruption of the target site were identified after the Sanger sequencing of PCR amplification of the genomic region surrounding the deleted locus (Figure 4D). 4930463O16Rik ^+/− mice were fertile and grew normally. Male and female 4930463O16Rik ^+/− animals were mated to obtain 4930463O16Rik^–/– mice. Once the mouse line had been established, all pups were genotyped with a combination of primers (listed in Supplementary Table 1; Figure 4B).

FIGURE 4

Figure 4. Deletion of the 4930463O16Rik gene in the mouse. (A) Schematic design of CRISPR/Cas9 deletion of the 4930463O16Rik gene with the suppression of four exons and three introns. The white boxes and lines represent exons and introns, respectively. (B) PCR genotyping on DNA of WT and 4930463O161Rik^–/– mice. The primer pairs used 1F/4R (located in 5′ of exon 1 and in 3′ of exon 4 of the 4930463O16Rik gene, respectively) or 2F/2R (located in the exon2 of 4930463O16Rik) to determine the genotypes of the mice. Results showed the following amplicon sizes: (*) 352 bp with the 2F/2R primers in WT (no amplification in mutant mice); (#) 935 bp with the 1F/4R primers in 4930463O161Rik^–/– mice (no amplification in WT mice under the PCR conditions used). (L) DNA ladder. (C) Transcription of the forward strand of chromosome 10 around the 4930463O16Rik gene with RNA-seq coverage (BigWig format representation) in WT (top blue tracks) and 4930463O16Rik^–/– (bottom tracks) mouse P18 testes. A continuous (WT) or very low transcription (4930463O16Rik^–/–) was observed from 4930463O16Rik to E230014E18Rik genes. (D) Electrophoregram of 4930463O16Rik^–/– mouse genomic DNA showing 9997 bp deletion and the insertion of three nucleotides (GTT, highlighted in pink).

The Absence of 4930463O16Rik Does Not Affect Mouse Fertility

Fertility was then investigated in 4930463O16Rik-deficient mice. Eight-week-old male and female 4930463O16Rik^–/– mice were mated, and both sexes were fertile. Their litter sizes (7.5 ± 2.10 pups per litter, n = 28) were similar to those of their WT counterparts (6.9 ± 2.12 pups per litter, n = 20). There were no significant differences in terms of testicular size, testis morphology and histology and cauda and caput epididymis between WT and 4930463O16Rik^–/– adult mice (Figures 5A,B). In addition, the different stages of seminiferous tubules divided into seven groups were quantified between 4930463O16Rik^–/– and WT adult mice. No significant differences were observed between the two genotypes (Figure 5C). These results therefore demonstrated that 4930463O16Rik is not required for mouse fertility.

FIGURE 5

Figure 5. Study of 4930463O16Rik^–/– testicular phenotype. (A) Testis/body weight ratio of 8-week-old mice. No significant difference was observed in the two mouse lines. Median lines are in black. (B) Hematoxylin and eosin (HE) staining of testis and epididymis sections from WT and 4930463O16Rik^–/– 8-week-old mice. Scale bar = 50 μm. Spermatozoa were visible in the lumen of the testes and epididymis of WT and 4930463O16Rik^–/– mice. (C) Quantification of the different seminiferous epithelium stages in WT and 4930463O16Rik^–/– 8-week-old mice. No significant difference was found between WT and 4930463O16Rik^–/– mice.

4930463O16Rik^–/– Mice Display Modified Sperm Parameters

The sperm parameters of 8-week-old 4930463O16Rik-deficient testes were compared to WT testes of the same age. Sperm concentrations obtained from the epididymis of 4930463O16Rik^–/– mice were significantly reduced by 57.2% compared to WT (Figure 6A) despite an unmodified testis/body weight ratio (Figure 5A). Motility parameters such as the percentage of motile spermatozoa, the motile mean expressed as beat cross frequency (bcf) and progressive spermatozoa were significantly higher in 4930463O16Rik^–/– mice compared to WT (Figures 6B–D). From a morphological point of view, however, two parameters were significantly modified in the testes of mutant mice: the distal mid-piece reflex (DMR), a defect developing in the epididymis and indicative of a sperm tail abnormality (Johnson, 1997) and the percentage of spermatozoa with coiled tail (Figures 6E,F). In addition, two kinetic parameters were also significantly reduced in mutant sperm: the motile mean vsl (related to the progressive velocity in a straight line) and the average path velocity vap; both parameters measuring sperm velocity (Figures 6G,H).

FIGURE 6

Figure 6. Evaluation of sperm parameters. Comparison of sperm-specific parameters from WT (circle, n = 11) and 4930463O16Rik^–/– (square, n = 20) mice. Significantly affected sperm parameters were (A) spermatozoa concentration (10⁶/mL), (B) spermatozoa motility (%), (C) motile mean bcf (beat cross frequency), (D) progressive spermatozoa (%), (E) DMR (distal midpiece reflex, abnormality of the sperm tail (%), (F) spermatozoa with coiled tail (%), (G) motile mean VSL (μm/s), and (H) VAP (μm/s). Statistical analyses were performed using the non-parametric Kruskal–Wallis test. ^∗p-val < 0.05, ^∗∗p-val < 0.01.

These results obtained using computer-aided sperm analysis (CASA) thus showed that several sperm parameters; namely concentration, motility, morphology and kinetics were impacted in 4930463O16Rik lncRNA-deficient mice. Some of them might negatively impact fertility, such as the sperm concentration, the DMR, the percentage of coiled tail and sperm velocity, while others would tend to suggest increased fertility, such as the motile mean percentage and bcf, and the progressive spermatozoa. These observations might have explained the normal fertility of 4930463O16Rik lncRNA-deficient mice.

Normal Male Fertility Despite a Modified Transcriptome in 4930463O16Rik^–/– Mouse Testes

Transcriptomic RNA-seq analyses were performed in WT and 4930463O16Rik^–/– mouse testes at two developmental stages; i.e., at P16 and 18, as in the Topaz1^–/– mouse line.

At P16, seven genes were differentially expressed (adjusted p-value < 0.05; absolute Log2FC > 1), including 4930463O16Rik, 1700092E16Rik, and E230014E18Rik (Supplementary Table 6). These latter two Riken cDNAs are in fact situated within the 3′ transcribed RNA of 4930463O16Rik (positioned in Figure 4C) and correspond to a unique locus. The transcriptional activity of this new locus stops toward the 3′ end of the cKap4 gene (cytoskeleton-associated protein 4 or Climp-63). This gene was down-regulated 1.7-fold in both knock-out lines (Topaz1 and 4930463O16Rik), which could suggest a newly discovered positive regulatory role for this lncRNA on the cKap4 gene.

At P18, 258 genes were differentially expressed (199 down-regulated and 59 up-regulated using the same statistical parameters; Supplementary Table 6). Among them, 206 were protein-coding genes accounting for 79.8% of DEGs (Figure 7A). Thus, P18 DEGs highlighted a direct or indirect relationship between the loss of the 4930463O16Rik lncRNA and protein-coding genes. In addition, loss of this lncRNA also resulted in the deregulation of 37 (14.3%) other lncRNAs.

FIGURE 7

Figure 7. Deregulated genes from 4930463O16Rik^–/– mouse testes. (A) Biotype of differentially expressed genes at P18 in 4930463O16Rik^–/– testes. Most of the deregulated genes are coding protein genes [adjusted p < 0.05 and down-regulated FC < 2 (log2FC < –1) or up-regulated FC > 2 (log2FC > 1]. (B) Venn diagram showing the overlap of differentially expressed genes between 4930463O16Rik^–/– and Topaz1^–/– mouse testes at 18 dpp.

The validation of several DEGs was performed by RT-qPCR using WT and 4930463O16Rik^–/– testicular RNAs at both developmental stages (P16 and P18) (Supplementary Figure 10). The qPCR results for the genes thus tested were consistent with those of RNA-seq.

The 4930463O16Rik^–/– DEGs were also analyzed using the DAVID database (Supplementary Table 7). At P18, six functional clusters had an enrichment score > 1.3 (Huang et al., 2009a). As for Topaz1^–/– mouse testes, they included the following GO terms: (i) for cellular components: cilium movement, ciliary part, axoneme, and (ii) for biological processes: microtubule-based process, regulation of cilium movement, spermatogenesis, male gamete generation, spermatid development and differentiation (Supplementary Table 7). An analysis that discriminated up-regulated from down-regulated genes only increased the value of the enrichment scores for the latter. The absence of TOPAZ1 protein or of 4930463O16Rik lncRNA caused the same enrichment clusters in mutant testes despite different outcomes regarding the fertility of male mice. The other clusters from the DAVID analysis referred to the GO terms: cell surface, external side of membrane, defense or immune response and response to external stimulus. These clusters were only found in the DAVID analysis with up-regulated genes.

Therefore, the 4930463O16Rik gene would appear to regulate genes related to spermatogenesis, microtubule or ciliary organization and the cytoskeleton in P18 testes. In the absence of this lncRNA, some genes involved in defense mechanisms or the immune response are also deregulated, suggesting stressful conditions. It should be noted that the majority (228/258 or 88%) of the DEGs from P18-4930463O16Rik^–/– testes was in common with those deregulated in Topaz1^–/– mice (Figure 7B). This led to similar results following ontological analyses of the DEGs of the two mutant lines. In the Topaz1^–/– testes, these 228 genes could be a consequence of the down-regulation of 4930463O16Rik lncRNA. On the other hand, these 228 DEGs alone do not explain meiotic arrest in the Topaz1^–/– testes.

Discussion

Topaz1 was initially reported as a germ-cell specific factor (Baillet et al., 2011) essential for meiotic progression and male fertility in mice (Luangpraseuth-Prosper et al., 2015). The suppression of Topaz1 led to an arrest of meiosis progression at the diplotene-metaphase I transition associated with germ cell apoptosis. Moreover, an initial transcriptomic approach, based on DNA microarrays, enabled the observation of a large but not exhaustive repertoire of deregulated transcripts. Using this technology, 10% of differentially expressed probes were lncRNAs and presented deregulated expression in P20 Topaz1^–/– testes compared to WT.

During this study, we were able to show that the effects of theTopaz1 gene being absent were visible on the mouse testicular transcriptome as early as 16 days post-partum, i.e., before the first meiotic division and the production of haploid germ cells. These effects were amplified at 18 days post-partum, just before or at the very start of the first haploid germ cells appearing. The molecular pathways involved in the absence of TOPAZ1 form part of spermatogenesis and the establishment of the cell cytoskeleton. At these two stages, P16 and P18, about a quarter of the deregulated genes in testes were ncRNAs (mainly lncRNAs) some of which displayed almost no expression in Topaz1^–/– testes. Three of these lncRNAs were specific of spermatocytes and had different subcellular localizations in spermatocytes, suggesting functions in these male germ cells. Suppressing one of them did not prevent the production of haploid spermatids and spermatozoa, but halved the murine sperm concentration. Furthermore, by deleting ∼10 kb corresponding to this 4930463O16Rik lincRNA, we showed that the transcriptional extinction was even longer, encompassing ∼35 kb in total and two other genes [1700092E16Rik (unknown gene type according to Ensembl) and the lincRNA E230014E18Rik]. Indeed, our transcriptional data suggested that these three annotated loci belong to a unique gene. Transcription of this lincRNA ended near the 3′ region of the cKap4 gene, known to be associated with the cytoskeleton (Vedrenne and Hauri, 2006). Remarkably, cKap4 expression was down-regulated 1.7 fold in both types of knockout mouse (Topaz1^–/– and 4930463O16Rik^–/–), suggesting a previously unknown and positive regulatory role of 4930463O16Rik on cKap4.

Topaz1 Ablation Down-Regulates a Significant Number of Cytoskeleton-Related Genes

Meiosis is well-orchestrated sequences of events controlled by different genes. In males, meiosis is a continuous process during the post-natal period, just before puberty, and results in the formation of four male gametes from one spermatocyte. During the first reductional division of meiosis, several checkpoint mechanisms regulate and coordinate prophase I during mammalian meiosis. This has been demonstrated during the past 25 years by the use of a large number of mutant mouse models, mainly gene knockout mice (Morelli and Cohen, 2005; Handel and Schimenti, 2010; Su et al., 2011).

A major checkpoint in males is the synaptic checkpoint that controls zygotene-pachytene transition, highlighted in male mice lacking Sycp3, Dmc1, Spo11, mei1, Msh4-5 or OvoL1 genes (Pittman et al., 1998; Edelmann et al., 1999; Baudat et al., 2000; Kneitz et al., 2000; Romanienko and Camerini-Otero, 2000; Yuan et al., 2000; Li et al., 2005; Reinholdt and Schimenti, 2005). Other gene suppressions have highlighted a second meiotic checkpoint of metaphase I in males (such as those of Mlh1-3, or cyclin A1) due to a misalignment of chromosomes on the spindle (Liu et al., 1998; Eaker et al., 2002; Lipkin et al., 2002). Male mice devoid of the Topaz1 gene belong to the latter models since Topaz1 is essential for the diplotene-metaphase I transition.

Absence of the Topaz1 gene disturbs the transcriptome of murine testes as early as 16 days postnatal. Of the 205 DEGs at P16, 85 were specific to this stage of development compared to P18 (Figure 1A), such as Ptgs1 (Cox1) a marker of peritubular cells (Rey-Ares et al., 2018), and Krt18 a marker of Sertoli cell maturation (Tarulli et al., 2006). Moreover, different genes involved in the TGFβ pathway were also P16-DEGs, such as Bmpr1b, Amh, and Fstl3. The latter, for example, was demonstrated to reduce the number of Sertoli cells in mouse testes and to limit their size (Oldknow et al., 2013). Ptgds (L-Pgds), which plays a role in the PGD2 molecular pathway during mammalian testicular organogenesis, is also deregulated (Moniot et al., 2009). All these genes, specifically deregulated at P16 due to the absence of Topaz1, thus appeared to participate in regulating cell communication. Absence of germ cell-specific genes lead to deregulation of Sertoli cell genes. In mammals, Sertoli cells are the only somatic cells that directly interact with spermatogenic cells (Griswold, 1998). Sertoli cells offer nutritional, structural, and regulatory support for germ cells during their whole development (Hess and Renato de Franca, 2008; Zhang et al., 2011). Therefore, it is not surprising to find deregulated Sertoli genes in case where spermatogenesis failure occurs due to germ cell meiotic arrest. At P18, these genes were no longer differentially expressed but they could be replaced by genes belonging to the same gene families, such as the cadherin or keratin families. It is easily conceivable that, in the near future, studies will define the transcriptomic evolution of Sertoli cells during the different stages of spermatogenesis.

Many of the 205 DEGs at P16 were involved in defense response pathways. For example, Ifit3 and Gbp3 are immune response genes in spermatocyte-derived GC-2spd(ts) cells (Kurihara et al., 2017).

Two days later, at P18, just before the prophase I-metaphase 1 transition, ten times more genes were deregulated. Among the 120 DEGs common to P16 and P18, there was at least one gene that might be involved in meiosis such as Aym1, an activator of yeast meiotic promoters 1. The absence of Topaz1 led to a lack of the testicular expression of Aym1. This gene is germ cell-specific (Malcov et al., 2004). In male mice, Aym1 is expressed from 10 dpp in early meiotic spermatocytes. The small murine AYM1 protein (44 amino acids) is immunolocalized in the nucleus of primary spermatocytes, mainly late pachytene and diplotene, suggesting a nuclear role for AYM1 in germ cells during the first meiotic division (Malcov et al., 2004).

At P18, the testicular transcriptome of Topaz1^–/– mice was largely disturbed when compared to WT animals, and most DEGs were down-regulated (Figure 1F), suggesting that TOPAZ1 promotes gene expression in normal mice. As TOPAZ1 is predicted to be an RNA-binding protein, it is tempting to speculate that its absence disorganized ribonucleic-protein complexes, including their instabilities and degradation. This could partly explain why 90% of DEGs were down-regulated at P18 and included a large proportion of lincRNAs. These down-regulated genes at P18 concerned microtubule-based movement and microtubule-based processes, and cellular components relative to motile cilium, ciliary part, sperm flagellum and axoneme. In addition, DAVID analysis revealed GO terms such as centriole, microtubule and spermatogenesis. All these terms relate to elements of the cytoskeleton that are indispensable for mitotic and/or meiotic divisions, motility and differentiation and are also widely involved in spermiogenesis, as might be expected because most DEGS are testis-specific. The absence of Topaz1 leads to a failure to pass the second meiotic checkpoint, namely the diplotene-metaphase I transition. One working hypothesis is that elements of the centriole, microtubules and/or centrosome are impaired in Topaz1^–/– spermatocytes. Failure in diplotene-metaphase I transition occurs in different mouse knockout models. For example, aberrant prometaphase-like cells were observed in Mlh1- or Meioc-deficient testes (Meioc is down-regulated 1.51-fold in P18 in Topaz1^–/– testes) (Eaker et al., 2002; Abby et al., 2016). These mutant mice have been described as displaying an arrest of male meiosis, and testes devoid of haploid germ cells leading to male sterility like mice lacking the Topaz1 gene. In Topaz1^–/– testes, Mlh1 is not a DEG.

During spermatogenesis, the dysregulation of centrosome proteins may affect meiotic division and genome stability. The expression of centriolar protein-coding genes Cep126, Cep128, Cep63 were down-regulated (FC from 2.1 to 2.7 compared to WT) at P18 in Topaz1^–/– testes. CEP126 is localized with γ-tubulin on the centriole during the mitosis of hTERT-RPE-1 (human telomerase-immortalized retinal pigmented epithelial cells) (Bonavita et al., 2014) but has never been studied in germ cells during meiosis. CEP128 was localized to the mother centriole and required to regulate ciliary signaling in zebrafish (Mönnich et al., 2018). Cep128 deletion decreased the stability of centriolar microtubules in F9 cells (epithelial cells from testicular teratoma of mouse embryo) (Kashihara et al., 2019). Centriole separation normally occurs at the end of prophase I or in early metaphase I, and CEP63 is associated with the mother centrioles. The mouse model devoid of Cep63 leads to male infertility (Marjanović et al., 2015), and in spermatocytes from these mice, the centriole duplication was impaired.

Finally, our ontology analysis of Topaz1^–/– P18-DEGs revealed significant enrichment scores for the several clusters relative to the final structure of spermatozoa such as tetratricopeptide repeat (TPR) and dynein heavy chain (DNAH1) clusters. Dynein chains are macromolecular complexes connecting central or doublet pairs of microtubules together to form the flagellar axoneme, the motility apparatus of spermatozoa (ref in Miyata et al., 2020a). Dynein proteins have also been identified as being involved in the microtubule-based intracellular transport of vesicles, and in both mitosis and meiosis (Mountain and Compton, 2000).

The TPR or PPR (pentatricopeptide repeat) domains consist of several 34 or 36 amino acid repeats that make up αα-hairpin repeat units, respectively (D’Andrea and Regan, 2003). The functions of TPR or PPR proteins were firstly documented in plants and are involved in RNA editing (D’Andrea and Regan, 2003; Schmitz-Linneweber and Small, 2008). In the mouse, Cfap70, a tetratricopeptide repeat-containing gene, was shown to be expressed in the testes (Shamoto et al., 2018), or as Spag1 in late-pachytene spermatocytes or round spermatids (Takaishi and Huh, 1999). Moreover, Ttc21a knockout mice have displayed sperm structural defects of the flagella and the connecting piece. In humans, Ttc21a has been associated with asthenoteratospermia in the Chinese population (Liu et al., 2019). Numerous components of the intraflagellar transport (IFT) complex contain TPR. The regulatory mechanisms of these TPR domain containing proteins remain unclear. Several genes coding for such tetratricopeptide repeat-containing proteins are down-regulated in P18 testes devoid of Topaz1, such as Cfap70, Spag1, Tct21a, and Ift140. Based on TPRpred (Karpenahalli et al., 2007) that predicted TPR- or PPR-containing proteins, the TOPAZ1 protein was predicted to contain such domains; seven in mice (p-val = 7.5E-08, probability of being PPR = 46.80%) and ten in humans (p-val = 3.4E-09, probability of being PPR = 88.76%).

A recent study of single cell-RNA-seq from all types of homogeneous spermatogenetic cells identified clusters of cells at similar developmental stages (Chen et al., 2018). This study shown that most of the genes involved in spermiogenesis start being expressed from the early pachytene stage. This is consistent with our RNA-seq results. Taken together, these data indicate that the absence of Topaz1 down-regulated a significant number of cytoskeleton-related genes as early as 18 days post-natal, leading to the meiotic arrest.

Topaz1 Ablation Deregulates a High Proportion of lncRNAs

Differentially expressed genes between Topaz1^–/– and WT mouse testes also revealed a high proportion of deregulated lncRNAs. We showed that three lincRNAs, whose expression was almost abolished as early as P16 in Topaz1-deficient mouse testes, were testis- and germ cell-specific. We showed that these genes are expressed in spermatocytes and round spermatids, suggesting a role in spermatogenesis. Their functions are still unknown.

Several investigations have revealed that the testis allow the expression of many lncRNAs (Necsulea et al., 2014). In mammals, the testis is the organ with the highest transcription rate (Soumillon et al., 2013). However, during the long stage of prophase I, these levels of transcription are not consistent. Indeed, transcription is markedly reduced or even abolished in the entire nucleus of spermatocytes during the early stages of prophase I. This is accompanied in particular by the nuclear processes of DNA division, the pairing of homologous chromosomes and telomeric rearrangements (Bolcun-Filas and Schimenti, 2012; Baudat et al., 2013; Shibuya and Watanabe, 2014), and also by the appearance of MSCI (meiotic sex chromosome inactivation) markers (Page et al., 2012). These processes are supported by epigenetic changes such as histone modifications and the recruitment of specific histone variants (references in Page et al., 2012). Transcription then takes up an important role in late-pachytene to diplotene spermatocytes (Monesi, 1964). The aforementioned scRNA-seq study of individual spermatogenic cells showed that almost 80% of annotated autosomal lncRNAs were expressed in spermatogenetic cells, mainly in mid-pachytene- to metaphase I-spermatocytes but also in round spermatids (Chen et al., 2018). The three lncRNAs investigated during our study (4930463O16Rik, Gm21269, and 4921513H07Rik) were also expressed at these developmental stages in mouse testes (Chen et al., 2018; Li et al., 2021). In the latter study (Li et al., 2021), the authors identified certain male germline-associated lncRNAs as being potentially important to spermatogenesis in vivo, based on several computational and experimental data sets; these lncRNAs included Gm21269 and 4921513H07Rik. The localization of lncRNAs in cells may be indicative of their potential function (Chen, 2016). 4930463O16Rik is expressed in the nucleus of spermatocytes. As mentioned above, 4930463O16Rik may play a positive role in the expression of cKap4 at the neighboring locus. Some nuclear lncRNA are involved in regulating transcription with a cis-regulatory role, such as Malat1 or Air (Sleutels et al., 2002; Zhang et al., 2012) on a nearby gene. Other nuclear lncRNAs act in trans and regulate gene transcription at another locus, such as HOTAIR (Rinn et al., 2007; Chu et al., 2011). In addition, some cytoplasmic lncRNA have been shown to play a role in miRNA competition, acting as miRNA sponges or decoys (such as linc-MD1 in human myoblasts; Cesana et al., 2011). Gm21269 is localized in the cytoplasm and nuclei of spermatocytes during meiosis. Both cytoplasmic and nuclear lncRNAs may act as a molecular scaffold for the assembly of functional protein complexes, such as HOTAIR or Dali (Tsai et al., 2010; Chalei et al., 2014), regulating protein localization and/or direct protein degradation, or acting as an miRNA precursor (Cai and Cullen, 2007). Finally, multiple other roles can be observed for lncRNAs. For example, the Dali lincRNA locally regulates its neighboring Pou3f3 gene, acts as a molecular scaffold for POU3F3 protein and interacts with DNMT1 in regulating the DNA methylation status of CpG island-associated promoters in trans during neural differentiation (Chalei et al., 2014).

Finally, TOPAZ1 possess a CCCH domain (Baillet et al., 2011). CCCH-type zinc-finger proteins are RNA binding proteins with regulatory functions at all stages of mRNA metabolism, such as mRNA transport, sub-cellular localization and stability/degradation (Brown, 2005; Hall, 2005). In addition, TOPAZ1 contains seven TPR domains. Interactions between TPR domains and RNAs have already been described. For example, some IFIT proteins (interferon-induced proteins with tetratricopeptide repeats) inhibit the infection of many viruses by recognizing viral RNAs (Johnson et al., 2018). The high proportion of deregulated testicular lncRNAs following Topaz1 ablation, and the presence of CCCH and TPR domains in the TOPAZ1 protein, thus strongly suggest a role for TOPAZ1 in RNA recognition and binding.

The Deletion of One lncRNA Alters Sperm Parameters Without Affecting Fertility

To decipher the biological function of an lncRNA affected by Topaz1 invalidation, a mouse model devoid of 4930463O16Rik was produced, with the same genetic background as Topaz1^–/– mice. This knockout mouse model did not exhibit meiosis disruption and the fertility of these mutant mice remained intact under standard laboratory conditions. Using a similar approach, Sox30 is a testis-specific factor that is essential to obtain haploid germ cells during spermatogenesis (Bai et al., 2018). SOX30 regulates Dnajb8 expression, but the deletion of Dnajb8 is not essential for spermatogenesis and male fertility (Wang et al., 2020).

The fact that a gene is either testis-specific or highly testis-enriched is no indication that its deletion will impair male fertility. Some laboratories have recently generated several dozen testis-enriched knockout mouse lines and shown that all these genes are individually dispensable in terms of male fertility in mice (Iwamori et al., 2011; Miyata et al., 2016; Khan et al., 2018; Lu et al., 2019; Chotiner et al., 2020; He et al., 2020; Khan et al., 2020; Jamin et al., 2021). The same is true for lncRNAs since their abundant expression during spermatogenesis has prompted other laboratories to produce knockout mouse models of testis-specific lncRNAs. This was the case for 1700121C10Rik or lncRNA5512 lncRNAs where mutant mice were also fertile without variations in their sperm parameters (Li et al., 2020; Zhu et al., 2020). One working hypothesis might be that some lncRNAs may regulate subsets of functional spermatogenetic-gene expression, in line with their nuclear localization, by binding to their regulatory genomic region.

Nevertheless, although our 4930463O16Rik-knockout mouse model are fertile, several parameters were altered including reduction in epididymal sperm concentrations (by more than half) and sperm motility. In Tslrn1 knockout mice (testis-specific long non-coding RNA 1) the males were fertile and displayed significantly lower sperm levels (–20%) but no reduction in litter size, or major defects in testis histology or variations in sperm motility (Wichman et al., 2017). Moreover, in Kif9-mutant male mice (Kif9, a coding gene), no testes abnormalities were found (Miyata et al., 2020b) but sperm motility was impaired: the VSL and VAP velocity parameters were reduced, as in 4930463O16Rik knockout mice. The authors concluded that Kif9 mutant mice were still fertile and this was probably due to variations in the motility of individual spermatozoa; those with good motility could still fertilize oocytes. The same conclusion applies to 4930463O16Rik^–/– mice.

In our study, the suppression of a gene – in this case 4930463O16Rik lincRNA – whose expression is markedly down-regulated in the testes of sterile Topaz1^–/– mice (FC = 40), has no effect on spermatogenesis. Our data suggest that the expression of 4930463O16Rik is not essential for fertility and meiotic division but adds to the terminal differentiation of male germ cells. Outside the laboratory, in wild reproductive life, one might imagine that biological functions may differ under more natural conditions due to stress and reproductive competition. This has been shown in particular for Pkdrej-deficient male mice which are fertile, whereas the Pkdrej gene (polycystin family receptor for egg jelly), is important to postcopulatory reproductive selection (Sutton et al., 2008; Miyata et al., 2016).

In summary, Topaz1 is a gene that is essential for fertility in male mice. Its absence leads to meiotic arrest before the first division. Our RNA-Seq analyses highlighted that Topaz1 stabilizes the expression of numerous lncRNAs. The suppression of one of them is not essential to mouse fertility but it is necessary during the terminal differentiation of male germ cells to achieve optimal function.

The absence of a specific anti-TOPAZ1 antibody did not enable us to further advance in our understanding of its function during murine spermatogenesis. The creation of a Flag-tagged Topaz1 knockin mouse model will allow us to gain further insights, and Rip-seq experiments will enable the determination of RNA-TOPAZ1 complexes during spermatogenesis. Given the large number of lncRNAs expressed in meiotic testes, one hypothesis may be that the function of several lncRNAs – including 4930463O16Rik- is partly redundant with that of other testicular lncRNAs.

Data Availability Statement

The datasets generated for this study can be found in the sequence read archive at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA698440.

Ethics Statement

The animal study was reviewed and approved by Ethical Committee for Animal Experimentation covering Jouy-en-Josas (COMETHEA).

Author Contributions

EPa and BM-P: conceptualization and design, project designing, and funding acquisition. MC, EPo, LJ, and BM-P: formal analysis; MC, EPo, LJ, BP, JC, and ES: investigation and methodology. MC, EPo, and BM-P: validation and visualization. LJ: software and data curation. MC, EPa, and BM-P: write the original data. BM-P: supervision, project administration, and writing – review and editing. All authors read and approved the final version of the manuscript.

Funding

This work was supported by the INRAE “Biology of Reproduction, Environment, Epigenetics and Development” unit (https://www6.jouy.inrae.fr/breed_eng/) and by the ACI-PHASE, INRAE (FLAG-TOPAZ) (https://www.inrae.fr/en/divisions/phase). MC was supported by a doctoral fellowship from the Ministère délégué à l’Enseignement Supérieur et à la Recherche at Paris-Saclay University (https://www.universite-paris-saclay.fr/en) and received an additional fellowship from the PHASE department of INRAE for the last 2 months of her Ph.D.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Our warmest thanks go to Jean-Luc Vilotte for allowing lncRNA deletion to take place in his laboratory, and also for proofreading this manuscript. We would like to thank the TACGENE facility at U1154-UMR7196, MNHN, Paris, and especially Anne De Cian and Jean-Paul Concordet, for the synthesis of gRNAs and Cas9 mRNA. We are grateful to all members of our mouse experimental unit (IERP, INRA, Jouy en Josas, France). We would like to thank the @Bridge platform for use of the Agilent Bioanalyzer and for histology facilities (UMR 1313 GABI, Jouy-en-Josas, France), and particularly Marthe Vilotte for HE staining. We would also like to thank the MIMA2 platform for providing access to the virtual slide scanner (Panoramic SCAN, 3DHISTECH). Finally, we are grateful to all members of the “Gonad Differentiation and its Disturbances” team for our scientific and technical discussions on Monday mornings, either in person or by videoconference. Vicky Hawken was responsible for the English language editing of this manuscript. We thanks must also be extended to Drs Paul Fowler (University of Aberdeen) and Richard Lea (University of Nottingham) for English corrections of the revised version.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.700290/full#supplementary-material

Supplementary Figure 1 | ISH control probe expression in WT 2-month-old testes. (A) Probe targeting the bacterial gene dapB (dapB, Bacillus subtilis dihydrodipicolinate reductase) was used as a negative control; (B,C) are examples using positive control probes (PPIB, Mus musculus peptidylprolyl isomerase B and UBC, Homo sapiens ubiquitin C). Probes for housekeeping genes (PPIB and UBC, respectively weakly and highly expressed genes) were used as positive controls (brown dots).

Supplementary Figure 2 | Hematoxylin and eosin staining of wild-type testes (A) at P15; (B) P17; and (C) P20 in mice (C57/Bl6). At P15, seminiferous tubules contain spermatocytes that have advanced to early and mid-pachynema (blue arrowhead); at P17, early-mid diplotene cells (red arrowhead) appear. At P20, according to the seminiferous epithelium stages, metaphase I plate in spermatocytes are visible (yellow arrowhead) or first round spermatids (green arrowhead) resulting from the first meiotic division. Scale bar = 20 μm.

Supplementary Figure 3 | Validation of several DEGs by RT-qPCR (RNA-seq Topaz1^–/– vs. WT testes). Validation of several differentially expressed up- or down-regulated genes and of non-DEGs from RNA-seq analysis by qRT-PCR from P16 (A) or P18 (B) mouse testis RNAs. The lines represent the median of each genotype (blue: WT; red: Topaz1^–/–). A Kruskal–Wallis statistical test was performed (^∗p < 0.05).

Supplementary Figure 4 | Reprogenomic data on the dynamic expression of 4930463O16Rik. The dynamic expression of 4930463O16Rik in five different tissues from male and female adult mice (A), and in embryonic primordial germ cells and adult male germ cells (B). 4930463O16Rik is expressed in testes in germ cells during post-natal life. The strongest dynamic expression is found in pachytene spermatocytes. The X-axes represented normalized counts.

Supplementary Figure 5 | Reprogenomic data on the dynamic expression of Gm21269. The dynamic expression of Gm21269 in five different tissues from male and female adult mice (A), in embryonic primordial germ cells and adult male germ cells (B). Gm21269 is expressed in testes in germ cells during post-natal life. The strongest dynamic expression is found in pachytene spermatocytes. The X-axes represented normalized counts.

Supplementary Figure 6 | IGV representation of P18-testis RNA-seq. Expression of 4930463O16Rik (A), Gm21269 (B), and 4921513H07Rik (C) from BigWig files of strand-specific RNA-seq data. The first four tracks represent transcripts of WT testes at P18; the next three tracks represent transcripts of Topaz1^–/– testes at the same developmental stage. Representations of the genes (from mm10 or GRCm38) are shown at the bottom of each graph. A representation of the size of 4930463O16Rik (A), Gm21269 (B), and 492151H07Rik (C) transcripts (red) from Ensembl data (GRCm38) is shown at the top. 4930463O16Rik and 4921513H07Rik gene transcriptions overlap in 3′ or 5′, respectively.

Supplementary Figure 7 | Expression of Gm21269, 4930463O16Rik, and 492151H07Rik mRNAs in testes from 5 days to adulthood. Quantitative RT-PCR analysis of Gm21269, 4930463O16Rik and 492151H07Rik gene expressions at different developmental stages in WT (blue) and Topaz1^–/– (red) testes. The lines represent the median of each genotype. A Kruskal–Wallis statistical test was performed (^∗p < 0.05; ^∗∗p < 0.01).

Supplementary Figure 8 | ISH with PAS counterstained in WT mouse testes. Visualization of 4930463O16Rik (A), Gm21269 (B), and 4921513H07Rik (C) mRNAs, respectively, by ISH (brown dots) at different seminiferous epithelium stages highlighted by PAS staining. Scale bar = 20 μm.

Supplementary Figure 9 | LncRNA cellular localizations in testes from 2-month-old WT mice. ISH using (A) 4930463O16Rik, (D) Gm21269, and (G) 4921513H07Rik probes (red). (B,E,H) Immunofluorescence staining with γH2Ax antibody was performed at the same stage of seminiferous epithelium to identify male germ cells (green). (C,F,I) DAPI (blue), visualizing nuclear chromosomes, was merged with ISH (green) and IF (red) signals. Zooms in white squares show spermatocytes during prophase I. No colocation between the sex body (γH2Ax) and lncRNAs (red) was evident. Scale bar = 20 μm.

Supplementary Figure 10 | Validation of several DEGs by RT-qPCR (RNA-seq 4930463O16Rik^–/– vs. WT testes). Validation of several differentially expressed up- or down-regulated genes and of non-DEGs of RNA-seq analysis by RT-qPCR from P16 (A) or P18 (B) mouse testis RNAs. The lines represent the median of each genotype (blue: WT; red: 4930463O16Rik^–/–). A Kruskal–Wallis statistical test was performed (^∗p < 0.05).

Supplementary Table 1 | List of primers. List of different primers used during this study for genotyping, RT-qPCR and gRNAs.

Supplementary Table 2 | Casa system settings.

Supplementary Table 3 | List of DEGs in Topaz1^–/– testes compared to WT. List of deregulated genes in Topaz1 KO testes at P16 (sheet 1) and P18 (sheet 2) (adjusted p-value < 0.05 and absolute Log2FC > 1).

Supplementary Table 4 | Functional annotation of P16 DEGs (RNA-seq Topaz1^–/– vs. WT testes). DAVID functional Annotation Clustering (DAVID 6.8) analysis (based on GO terms and KEGG pathway) of all P16-differentially expressed genes (sheet 1) or only up-regulated DEGs (sheet 2) or down-regulated DEGs (sheet 3) in Topaz1^–/– testes.

Supplementary Table 5 | Functional annotation of P18 DEGs (RNA-seq Topaz1^–/– vs. WT testes). DAVID functional Annotation Clustering (DAVID 6.8) analysis (based on GO terms and KEGG pathway) of P18-differentially expressed genes (sheet 1) or only up-regulated DEGs (sheet 2) or down-regulated DEGs (sheet 3) in Topaz1^–/– testes. Annotation clusters based on the InterPro database of P18-down-regulated DEGs are mentioned in sheet 4.

Supplementary Table 6 | List of DEGs in 4930463O16Rik^–/– testes compared to WT. List of deregulated genes in 4930463O16Rik KO testes at P16 (sheet 1) and P18 (sheet 2) (adjusted p-value < 0.05 and absolute Log2FC > 1).

Supplementary Table 7 | Functional annotation of P18 DEGs (RNA-seq 4930463O16Rik^–/– vs. WT testes). DAVID functional Annotation Clustering (DAVID 6.8) analysis of P18-differentially expressed genes (sheet 1) or only up-regulated (sheet 2) or down-regulated DEGs (sheet 3) in 4930463O16Rik^–/– testes.

Footnotes

1. ^ http://crispor.tefor.net/
2. ^ https://www6.jouy.inrae.fr/mima2/
3. ^ https://doi.org/10.15454/1.5572348210007727E12
4. ^ https://www.ncbi.nlm.nih.gov/bioproject/PRJNA698440
5. ^ http://www.informatics.jax.org
6. ^ https://www.ensembl.org/
7. ^ ftp://ftp.ncbi.nih.gov/gene/DATA/GENE_INFO/Mammalia/Mus_musculus.gene_info.gz
8. ^ http://www.ensembl.org/biomart/martview
9. ^ http://www.informatics.jax.org/marker/
10. ^ https://david.ncifcrf.gov/
11. ^ https://maayanlab.cloud/Enrichr/
12. ^ https://rgv.genouest.org/
13. ^ http://software.broadinstitute.org/software/igv/

References

Abby, E., Tourpin, S., Ribeiro, J., Daniel, K., Messiaen, S., Moison, D., et al. (2016). Implementation of meiosis prophase I programme requires a conserved retinoid-independent stabilizer of meiotic transcripts. Nat. Commun. 7:10324. doi: 10.1038/ncomms10324

PubMed Abstract | CrossRef Full Text | Google Scholar

Anguera, M. C., Ma, W., Clift, D., Namekawa, S., Kelleher, R. J., and Lee, J. T. (2011). Tsx produces a long noncoding RNA and has general functions in the germline, stem cells, and brain. PLoS Genet. 7:e1002248. doi: 10.1371/journal.pgen.1002248

PubMed Abstract | CrossRef Full Text | Google Scholar

Arun, G., Akhade, V. S., Donakonda, S., and Rao, M. R. S. (2012). mrhl RNA, a long noncoding RNA, negatively regulates Wnt signaling through its protein partner Ddx5/p68 in mouse spermatogonial cells. Mol. Cell. Biol. 32, 3140–3152. doi: 10.1128/MCB.00006-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, S., Fu, K., Yin, H., Cui, Y., Yue, Q., Li, W., et al. (2018). Sox30 initiates transcription of haploid genes during late meiosis and spermiogenesis in mouse testes. Development 145:dev164855. doi: 10.1242/dev.164855

PubMed Abstract | CrossRef Full Text | Google Scholar

Baillet, A., Le Bouffant, R., Volff, J. N., Luangpraseuth, A., Poumerol, E., Thépot, D., et al. (2011). TOPAZ1, a novel germ cell-specific expressed gene conserved during evolution across vertebrates. PLoS One 6:e26950. doi: 10.1371/journal.pone.0026950

PubMed Abstract | CrossRef Full Text | Google Scholar

Baudat, F., Imai, Y., and de Massy, B. (2013). Meiotic recombination in mammals: localization and regulation. Nat. Rev. Genet. 14, 794–806. doi: 10.1038/nrg3573

PubMed Abstract | CrossRef Full Text | Google Scholar

Baudat, F., Manova, K., Yuen, J. P., Jasin, M., and Keeney, S. (2000). Chromosome synapsis defects and sexually dimorphic meiotic progression in mice lacking Spo11. Mol. Cell 6, 989–998. doi: 10.1016/s1097-2765(00)00098-8

CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B (Methodol.) 57, 289–300.

Google Scholar

Bie, B., Wang, Y., Li, L., Fang, H., Liu, L., and Sun, J. (2018). Noncoding RNAs: potential players in the self-renewal of mammalian spermatogonial stem cells. Mol. Reprod. Dev. 85, 720–728. doi: 10.1002/mrd.23041

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolcun-Filas, E., and Schimenti, J. C. (2012). Genetics of meiosis and recombination in mice. Int. Rev. Cell Mol. Biol. 298, 179–227. doi: 10.1016/B978-0-12-394309-5.00005-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonavita, R., Walas, D., Brown, A. K., Luini, A., Stephens, D. J., and Colanzi, A. (2014). Cep126 is required for pericentriolar satellite localisation to the centrosome and for primary cilium formation. Biol. Cell 106, 254–267. doi: 10.1111/boc.201300087

PubMed Abstract | CrossRef Full Text | Google Scholar

Brannan, C. I., Dees, E. C., Ingram, R. S., and Tilghman, S. M. (1990). The product of the H19 gene may function as an RNA. Mol. Cell. Biol. 10, 28–36. doi: 10.1128/mcb.10.1.28

CrossRef Full Text | Google Scholar

Brockdorff, N., Ashworth, A., Kay, G. F., McCabe, V. M., Norris, D. P., Cooper, P. J., et al. (1992). The product of the mouse Xist gene is a 15 kb inactive X-specific transcript containing no conserved ORF and located in the nucleus. Cell 71, 515–526. doi: 10.1016/0092-8674(92)90519-i

CrossRef Full Text | Google Scholar

Brown, R. S. (2005). Zinc finger proteins: getting a grip on RNA. Curr. Opin. Struct. Biol. 15, 94–98. doi: 10.1016/j.sbi.2005.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, X., and Cullen, B. R. (2007). The imprinted H19 noncoding RNA is a primary microRNA precursor. RNA 13, 313–316. doi: 10.1261/rna.351707

PubMed Abstract | CrossRef Full Text | Google Scholar

Carninci, P., Kasukawa, T., Katayama, S., Gough, J., Frith, M. C., Maeda, N., et al. (2005). The transcriptional landscape of the mammalian genome. Science 309, 1559–1563. doi: 10.1126/science.1112014

PubMed Abstract | CrossRef Full Text | Google Scholar

Cesana, M., Cacchiarelli, D., Legnini, I., Santini, T., Sthandier, O., Chinappi, M., et al. (2011). A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell 147, 358–369. doi: 10.1016/j.cell.2011.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Chalei, V., Sansom, S. N., Kong, L., Lee, S., Montiel, J. F., Vance, K. W., et al. (2014). The long non-coding RNA Dali is an epigenetic regulator of neural differentiation. Elife 3:e04530. doi: 10.7554/eLife.04530

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L.-L. (2016). Linking long noncoding RNA localization and function. Trends Biochem. Sci. 41, 761–772. doi: 10.1016/j.tibs.2016.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Zheng, Y., Gao, Y., Lin, Z., Yang, S., Wang, T., et al. (2018). Single-cell RNA-seq uncovers dynamic processes and critical regulators in mouse spermatogenesis. Cell Res. 28, 879–896. doi: 10.1038/s41422-018-0074-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Chotiner, J. Y., Leu, N. A., and Wang, P. J. (2020). FLACC1 is testis-specific but dispensable for fertility in mice. Mol. Reprod. Dev. 87, 1199–1201. doi: 10.1002/mrd.23435

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, C., Qu, K., Zhong, F. L., Artandi, S. E., and Chang, H. Y. (2011). Genomic maps of long noncoding RNA occupancy reveal principles of RNA-chromatin interactions. Mol. Cell 44, 667–678. doi: 10.1016/j.molcel.2011.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

D’Andrea, L. D., and Regan, L. (2003). TPR proteins: the versatile helix. Trends Biochem. Sci. 28, 655–662. doi: 10.1016/j.tibs.2003.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Darde, T. A., Lecluze, E., Lardenois, A., Stévant, I., Alary, N., Tüttelmann, F., et al. (2019). The ReproGenomics Viewer: a multi-omics and cross-species resource compatible with single-cell studies for the reproductive science community. Bioinformatics 35, 3133–3139. doi: 10.1093/bioinformatics/btz047

PubMed Abstract | CrossRef Full Text | Google Scholar

Darde, T. A., Sallou, O., Becker, E., Evrard, B., Monjeaud, C., Le Bras, Y., et al. (2015). The ReproGenomics Viewer: an integrative cross-species toolbox for the reproductive science community. Nucleic Acids Res. 43, W109–W116. doi: 10.1093/nar/gkv345

PubMed Abstract | CrossRef Full Text | Google Scholar

Drumond, A. L., Meistrich, M. L., and Chiarini-Garcia, H. (2011). Spermatogonial morphology and kinetics during testis development in mice: a high-resolution light microscopy approach. Reproduction 142, 145–155. doi: 10.1530/REP-10-0431

PubMed Abstract | CrossRef Full Text | Google Scholar

Eaker, S., Cobb, J., Pyle, A., and Handel, M. A. (2002). Meiotic prophase abnormalities and metaphase cell death in MLH1-deficient mouse spermatocytes: insights into regulation of spermatogenic progress. Dev. Biol. 249, 85–95. doi: 10.1006/dbio.2002.0708

PubMed Abstract | CrossRef Full Text | Google Scholar

Edelmann, W., Cohen, P. E., Kneitz, B., Winand, N., Lia, M., Heyer, J., et al. (1999). Mammalian MutS homologue 5 is required for chromosome pairing in meiosis. Nat. Genet. 21, 123–127. doi: 10.1038/5075

PubMed Abstract | CrossRef Full Text | Google Scholar

ENCODE Project Consortium, Birney, E., Stamatoyannopoulos, J. A., Dutta, A., Guigó, R., Gingeras, T. R., et al. (2007). Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature 447, 799–816. doi: 10.1038/nature05874

PubMed Abstract | CrossRef Full Text | Google Scholar

Gil, R., and Latorre, A. (2012). Factors behind junk DNA in bacteria. Genes (Basel) 3, 634–650. doi: 10.3390/genes3040634

PubMed Abstract | CrossRef Full Text | Google Scholar

Griswold, M. D. (1998). The central role of Sertoli cells in spermatogenesis. Semin. Cell Dev. Biol. 9, 411–416. doi: 10.1006/scdb.1998.0203

PubMed Abstract | CrossRef Full Text | Google Scholar

Guttman, M., Amit, I., Garber, M., French, C., Lin, M. F., Feldser, D., et al. (2009). Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature 458, 223–227. doi: 10.1038/nature07672

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, T. M. T. (2005). Multiple modes of RNA recognition by zinc finger proteins. Curr. Opin. Struct. Biol. 15, 367–373. doi: 10.1016/j.sbi.2005.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, J., Zhang, J., Chen, L., Shen, B., Zhou, J., Hu, B., et al. (2014). Efficient in vivo deletion of a large imprinted lncRNA by CRISPR/Cas9. RNA Biol. 11, 829–835. doi: 10.4161/rna.29624

PubMed Abstract | CrossRef Full Text | Google Scholar

Handel, M. A., and Schimenti, J. C. (2010). Genetics of mammalian meiosis: regulation, dynamics and impact on fertility. Nat. Rev. Genet. 11, 124–136. doi: 10.1038/nrg2723

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., Xie, W., Li, H., Cui, Y., Wang, Y., Guo, X., et al. (2020). The testis-specifically expressed gene Trim69 is not essential for fertility in mice. J. Biomed. Res. 35, 47–60. doi: 10.7555/JBR.34.20200069

PubMed Abstract | CrossRef Full Text | Google Scholar

Hellemans, J., Mortier, G., De Paepe, A., Speleman, F., and Vandesompele, J. (2007). qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol. 8:R19. doi: 10.1186/gb-2007-8-2-r19

PubMed Abstract | CrossRef Full Text | Google Scholar

Henao-Mejia, J., Williams, A., Rongvaux, A., Stein, J., Hughes, C., and Flavell, R. A. (2016). Generation of genetically modified mice using the CRISPR-Cas9 genome-editing system. Cold Spring Harb. Protoc. 2016:pdb.prot090704. doi: 10.1101/pdb.prot090704

PubMed Abstract | CrossRef Full Text | Google Scholar

Hess, R. A., and Renato de Franca, L. (2008). Spermatogenesis and cycle of the seminiferous epithelium. Adv. Exp. Med. Biol. 636, 1–15. doi: 10.1007/978-0-387-09597-4_1

CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009a). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D. W., Sherman, B. T., Zheng, X., Yang, J., Imamichi, T., Stephens, R., et al. (2009b). Extracting biological meaning from large gene lists with DAVID. Curr. Protoc. Bioinformatics Chapter 13, Unit13.11. doi: 10.1002/0471250953.bi1311s27

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwamori, N., Zhao, M., Meistrich, M. L., and Matzuk, M. M. (2011). The testis-enriched histone demethylase, KDM4D, regulates methylation of histone H3 lysine 9 during spermatogenesis in the mouse but is dispensable for fertility. Biol. Reprod. 84, 1225–1234. doi: 10.1095/biolreprod.110.088955

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamin, S. P., Petit, F. G., Demini, L., and Primig, M. (2021). Tex55 encodes a conserved putative A-kinase anchoring protein dispensable for male fertility in the mouse. Biol. Reprod. 104, 731–733. doi: 10.1093/biolre/ioab007

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarroux, J., Morillon, A., and Pinskaya, M. (2017). History, discovery, and classification of lncRNAs. Adv. Exp. Med. Biol. 1008, 1–46. doi: 10.1007/978-981-10-5203-3_1

CrossRef Full Text | Google Scholar

Johnson, B., VanBlargan, L. A., Xu, W., White, J. P., Shan, C., Shi, P.-Y., et al. (2018). Human IFIT3 modulates IFIT1 RNA binding specificity and protein stability. Immunity 48, 487–499.e5. doi: 10.1016/j.immuni.2018.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, W. H. (1997). The significance to bull fertility of morphologically abnormal sperm. Vet. Clin. North Am. Food Anim. Pract. 13, 255–270. doi: 10.1016/s0749-0720(15)30339-x

CrossRef Full Text | Google Scholar

Karpenahalli, M. R., Lupas, A. N., and Söding, J. (2007). TPRpred: a tool for prediction of TPR-, PPR- and SEL1-like repeats from protein sequences. BMC Bioinformatics 8:2. doi: 10.1186/1471-2105-8-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashihara, H., Chiba, S., Kanno, S.-I., Suzuki, K., Yano, T., and Tsukita, S. (2019). Cep128 associates with Odf2 to form the subdistal appendage of the centriole. Genes Cells 24, 231–243. doi: 10.1111/gtc.12668

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, M., Jabeen, N., Khan, T., Hussain, H. M. J., Ali, A., Khan, R., et al. (2018). The evolutionarily conserved genes: Tex37, Ccdc73, Prss55 and Nxt2 are dispensable for fertility in mice. Sci. Rep. 8:4975. doi: 10.1038/s41598-018-23176-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, R., Ye, J., Yousaf, A., Shah, W., Aftab, A., Shah, B., et al. (2020). Evolutionarily conserved and testis-specific gene, 4930524B15Rik, is not essential for mouse spermatogenesis and fertility. Mol. Biol. Rep. 47, 5207–5213. doi: 10.1007/s11033-020-05595-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleckner, N. (1996). Meiosis: how could it work? Proc. Natl. Acad. Sci. U.S.A. 93, 8167–8174. doi: 10.1073/pnas.93.16.8167

PubMed Abstract | CrossRef Full Text | Google Scholar

Kneitz, B., Cohen, P. E., Avdievich, E., Zhu, L., Kane, M. F., Hou, H., et al. (2000). MutS homolog 4 localization to meiotic chromosomes is required for chromosome pairing during meiosis in male and female mice. Genes Dev. 14, 1085–1097.

Google Scholar

Kurihara, M., Otsuka, K., Matsubara, S., Shiraishi, A., Satake, H., and Kimura, A. P. (2017). A testis-specific long non-coding RNA, lncRNA-Tcam1, regulates immune-related genes in mouse male germ cells. Front. Endocrinol. (Lausanne) 8:299. doi: 10.3389/fendo.2017.00299

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Nair, M., Mackay, D. R., Bilanchone, V., Hu, M., Fallahi, M., et al. (2005). Ovol1 regulates meiotic pachytene progression during spermatogenesis by repressing Id2 expression. Development 132, 1463–1473. doi: 10.1242/dev.01658

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Shen, C., Shang, X., Tang, L., Xiong, W., Ge, H., et al. (2020). Two novel testis-specific long noncoding RNAs produced by 1700121C10Rik are dispensable for male fertility in mice. J. Reprod. Dev. 66, 57–65. doi: 10.1262/jrd.2019-104

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Xu, J., Luo, Y., Zou, D., Han, R., Zhong, S., et al. (2021). Panoramic transcriptome analysis and functional screening of long noncoding RNAs in mouse spermatogenesis. Genome Res. 31, 13–26. doi: 10.1101/gr.264333.120

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656

PubMed Abstract | CrossRef Full Text | Google Scholar

Lipkin, S. M., Moens, P. B., Wang, V., Lenzi, M., Shanmugarajah, D., Gilgeous, A., et al. (2002). Meiotic arrest and aneuploidy in MLH3-deficient mice. Nat. Genet. 31, 385–390. doi: 10.1038/ng931

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D., Matzuk, M. M., Sung, W. K., Guo, Q., Wang, P., and Wolgemuth, D. J. (1998). Cyclin A1 is required for meiosis in the male mouse. Nat. Genet. 20, 377–380. doi: 10.1038/3855

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., He, X., Yang, S., Zouari, R., Wang, J., Wu, H., et al. (2019). Bi-allelic mutations in TTC21A induce asthenoteratospermia in humans and mice. Am. J. Hum. Genet. 104, 738–748. doi: 10.1016/j.ajhg.2019.02.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Y., Oura, S., Matsumura, T., Oji, A., Sakurai, N., Fujihara, Y., et al. (2019). CRISPR/Cas9-mediated genome editing reveals 30 testis-enriched genes dispensable for male fertility in mice†. Biol. Reprod. 101, 501–511. doi: 10.1093/biolre/ioz103

PubMed Abstract | CrossRef Full Text | Google Scholar

Luangpraseuth-Prosper, A., Lesueur, E., Jouneau, L., Pailhoux, E., Cotinot, C., and Mandon-Pépin, B. (2015). TOPAZ1, a germ cell specific factor, is essential for male meiotic progression. Dev. Biol. 406, 158–171. doi: 10.1016/j.ydbio.2015.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Malcov, M., Cesarkas, K., Stelzer, G., Shalom, S., Dicken, Y., Naor, Y., et al. (2004). Aym1, a mouse meiotic gene identified by virtue of its ability to activate early meiotic genes in the yeast Saccharomyces cerevisiae. Dev. Biol. 276, 111–123. doi: 10.1016/j.ydbio.2004.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Marjanović, M., Sánchez-Huertas, C., Terré, B., Gómez, R., Scheel, J. F., Pacheco, S., et al. (2015). CEP63 deficiency promotes p53-dependent microcephaly and reveals a role for the centrosome in meiotic recombination. Nat. Commun. 6:7676. doi: 10.1038/ncomms8676

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyata, H., Castaneda, J. M., Fujihara, Y., Yu, Z., Archambeault, D. R., Isotani, A., et al. (2016). Genome engineering uncovers 54 evolutionarily conserved and testis-enriched genes that are not required for male fertility in mice. Proc. Natl. Acad. Sci. U.S.A. 113, 7704–7710. doi: 10.1073/pnas.1608458113

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyata, H., Morohoshi, A., and Ikawa, M. (2020a). Analysis of the sperm flagellar axoneme using gene-modified mice. Exp. Anim. 69, 374–381. doi: 10.1538/expanim.20-0064

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyata, H., Shimada, K., Morohoshi, A., Oura, S., Matsumura, T., Xu, Z., et al. (2020b). Testis-enriched kinesin KIF9 is important for progressive motility in mouse spermatozoa. FASEB J. 34, 5389–5400. doi: 10.1096/fj.201902755R

PubMed Abstract | CrossRef Full Text | Google Scholar

Monesi, V. (1964). Ribonucleic acid synthesis during mitosis and meiosis in the mouse testis. J. Cell Biol. 22, 521–532. doi: 10.1083/jcb.22.3.521

PubMed Abstract | CrossRef Full Text | Google Scholar

Moniot, B., Declosmenil, F., Barrionuevo, F., Scherer, G., Aritake, K., Malki, S., et al. (2009). The PGD2 pathway, independently of FGF9, amplifies SOX9 activity in Sertoli cells during male sexual differentiation. Development 136, 1813–1821. doi: 10.1242/dev.032631

PubMed Abstract | CrossRef Full Text | Google Scholar

Mönnich, M., Borgeskov, L., Breslin, L., Jakobsen, L., Rogowski, M., Doganli, C., et al. (2018). CEP128 localizes to the subdistal appendages of the mother centriole and regulates TGF-β/BMP signaling at the primary cilium. Cell Rep. 22, 2584–2592. doi: 10.1016/j.celrep.2018.02.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Morelli, M. A., and Cohen, P. E. (2005). Not all germ cells are created equal: aspects of sexual dimorphism in mammalian meiosis. Reproduction 130, 761–781. doi: 10.1530/rep.1.00865

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, K. V., and Mattick, J. S. (2014). The rise of regulatory RNA. Nat. Rev. Genet. 15, 423–437. doi: 10.1038/nrg3722

PubMed Abstract | CrossRef Full Text | Google Scholar

Mountain, V., and Compton, D. A. (2000). Dissecting the role of molecular motors in the mitotic spindle. Anat. Rec. 261, 14–24.

Google Scholar

Necsulea, A., and Kaessmann, H. (2014). Evolutionary dynamics of coding and non-coding transcriptomes. Nat. Rev. Genet. 15, 734–748. doi: 10.1038/nrg3802

PubMed Abstract | CrossRef Full Text | Google Scholar

Necsulea, A., Soumillon, M., Warnefors, M., Liechti, A., Daish, T., Zeller, U., et al. (2014). The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature 505, 635–640. doi: 10.1038/nature12943

PubMed Abstract | CrossRef Full Text | Google Scholar

Oldknow, K. J., Seebacher, J., Goswami, T., Villen, J., Pitsillides, A. A., O’Shaughnessy, P. J., et al. (2013). Follistatin-like 3 (FSTL3) mediated silencing of transforming growth factor β (TGFβ) signaling is essential for testicular aging and regulating testis size. Endocrinology 154, 1310–1320. doi: 10.1210/en.2012-1886

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, J., de la Fuente, R., Manterola, M., Parra, M. T., Viera, A., Berríos, S., et al. (2012). Inactivation or non-reactivation: what accounts better for the silence of sex chromosomes during mammalian male meiosis? Chromosoma 121, 307–326. doi: 10.1007/s00412-012-0364-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pittman, D. L., Cobb, J., Schimenti, K. J., Wilson, L. A., Cooper, D. M., Brignull, E., et al. (1998). Meiotic prophase arrest with failure of chromosome synapsis in mice deficient for Dmc1, a germline-specific RecA homolog. Mol. Cell 1, 697–705. doi: 10.1016/s1097-2765(00)80069-6

CrossRef Full Text | Google Scholar

Quan, G., and Li, J. (2018). Circular RNAs: biogenesis, expression and their potential roles in reproduction. J. Ovarian Res. 11:9. doi: 10.1186/s13048-018-0381-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Reinholdt, L. G., and Schimenti, J. C. (2005). Mei1 is epistatic to Dmc1 during mouse meiosis. Chromosoma 114, 127–134. doi: 10.1007/s00412-005-0346-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rey-Ares, V., Rossi, S. P., Dietrich, K.-G., Köhn, F.-M., Schwarzer, J. U., Welter, H., et al. (2018). Prostaglandin E2 (PGE2) is a testicular peritubular cell-derived factor involved in human testicular homeostasis. Mol. Cell. Endocrinol. 473, 217–224. doi: 10.1016/j.mce.2018.01.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Rinn, J. L., Kertesz, M., Wang, J. K., Squazzo, S. L., Xu, X., Brugmann, S. A., et al. (2007). Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell 129, 1311–1323. doi: 10.1016/j.cell.2007.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Romanienko, P. J., and Camerini-Otero, R. D. (2000). The mouse Spo11 gene is required for meiotic chromosome synapsis. Mol. Cell 6, 975–987. doi: 10.1016/s1097-2765(00)00097-6

CrossRef Full Text | Google Scholar

Sarropoulos, I., Marin, R., Cardoso-Moreira, M., and Kaessmann, H. (2019). Developmental dynamics of lncRNAs across mammalian organs and species. Nature 571, 510–514. doi: 10.1038/s41586-019-1341-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitz-Linneweber, C., and Small, I. (2008). Pentatricopeptide repeat proteins: a socket set for organelle gene expression. Trends Plant Sci 13, 663–670. doi: 10.1016/j.tplants.2008.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shamoto, N., Narita, K., Kubo, T., Oda, T., and Takeda, S. (2018). CFAP70 is a novel axoneme-binding protein that localizes at the base of the outer dynein arm and regulates ciliary motility. Cells 7:124. doi: 10.3390/cells7090124

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibuya, H., and Watanabe, Y. (2014). The meiosis-specific modification of mammalian telomeres. Cell Cycle 13, 2024–2028. doi: 10.4161/cc.29350

PubMed Abstract | CrossRef Full Text | Google Scholar

Sleutels, F., Zwart, R., and Barlow, D. P. (2002). The non-coding Air RNA is required for silencing autosomal imprinted genes. Nature 415, 810–813. doi: 10.1038/415810a

PubMed Abstract | CrossRef Full Text | Google Scholar

Soumillon, M., Necsulea, A., Weier, M., Brawand, D., Zhang, X., Gu, H., et al. (2013). Cellular source and mechanisms of high transcriptome complexity in the mammalian testis. Cell Rep. 3, 2179–2190. doi: 10.1016/j.celrep.2013.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, Y., Li, Y., and Ye, P. (2011). Mammalian meiosis is more conserved by sex than by species: conserved co-expression networks of meiotic prophase. Reproduction 142, 675–687. doi: 10.1530/REP-11-0260

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutton, K. A., Jungnickel, M. K., and Florman, H. M. (2008). A polycystin-1 controls postcopulatory reproductive selection in mice. Proc. Natl. Acad. Sci. U.S.A. 105, 8661–8666. doi: 10.1073/pnas.0800603105

PubMed Abstract | CrossRef Full Text | Google Scholar

Takaishi, M., and Huh, N. H. (1999). A tetratricopeptide repeat-containing protein gene, tpis, whose expression is induced with differentiation of spermatogenic cells. Biochem. Biophys. Res. Commun. 264, 81–85. doi: 10.1006/bbrc.1999.1477

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarulli, G. A., Stanton, P. G., Lerchl, A., and Meachem, S. J. (2006). Adult sertoli cells are not terminally differentiated in the Djungarian hamster: effect of FSH on proliferation and junction protein organization. Biol. Reprod. 74, 798–806. doi: 10.1095/biolreprod.105.050450

PubMed Abstract | CrossRef Full Text | Google Scholar

Toyooka, Y., Tsunekawa, N., Takahashi, Y., Matsui, Y., Satoh, M., and Noce, T. (2000). Expression and intracellular localization of mouse Vasa-homologue protein during germ cell development. Mech. Dev. 93, 139–149. doi: 10.1016/s0925-4773(00)00283-5

CrossRef Full Text | Google Scholar

Trapnell, C., Pachter, L., and Salzberg, S. L. (2009). TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–1111. doi: 10.1093/bioinformatics/btp120

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, M.-C., Manor, O., Wan, Y., Mosammaparast, N., Wang, J. K., Lan, F., et al. (2010). Long noncoding RNA as modular scaffold of histone modification complexes. Science 329, 689–693. doi: 10.1126/science.1192002

PubMed Abstract | CrossRef Full Text | Google Scholar

Vedrenne, C., and Hauri, H.-P. (2006). Morphogenesis of the endoplasmic reticulum: beyond active membrane expansion. Traffic 7, 639–646. doi: 10.1111/j.1600-0854.2006.00419.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F., Kong, S., Hu, X., Li, X., Xu, B., Yue, Q., et al. (2020). Dnajb8, a target gene of SOX30, is dispensable for male fertility in mice. PeerJ 8:e10582. doi: 10.7717/peerj.10582

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, K., Yang, L., Xiong, T., Di, C., Ma, D., Wu, M., et al. (2016). Critical roles of long noncoding RNAs in Drosophila spermatogenesis. Genome Res. 26, 1233–1244. doi: 10.1101/gr.199547.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wichman, L., Somasundaram, S., Breindel, C., Valerio, D. M., McCarrey, J. R., Hodges, C. A., et al. (2017). Dynamic expression of long noncoding RNAs reveals their potential roles in spermatogenesis and fertility. Biol. Reprod. 97, 313–323. doi: 10.1093/biolre/iox084

PubMed Abstract | CrossRef Full Text | Google Scholar

Yadav, R. P., and Kotaja, N. (2014). Small RNAs in spermatogenesis. Mol. Cell. Endocrinol. 382, 498–508. doi: 10.1016/j.mce.2013.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, L., Liu, J. G., Zhao, J., Brundell, E., Daneholt, B., and Höög, C. (2000). The murine SCP3 gene is required for synaptonemal complex assembly, chromosome synapsis, and male fertility. Mol. Cell 5, 73–83. doi: 10.1016/s1097-2765(00)80404-9

CrossRef Full Text | Google Scholar

Zhang, B., Arun, G., Mao, Y. S., Lazar, Z., Hung, G., Bhattacharjee, G., et al. (2012). The lncRNA Malat1 is dispensable for mouse development but its transcription plays a cis-regulatory role in the adult. Cell Rep. 2, 111–123. doi: 10.1016/j.celrep.2012.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Eto, K., Honmyou, A., Nakao, K., Kiyonari, H., and Abé, S. (2011). Neuregulins are essential for spermatogonial proliferation and meiotic initiation in neonatal mouse testis. Development 138, 3159–3168. doi: 10.1242/dev.062380

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Lu, H., Xin, D., Cheng, H., and Zhou, R. (2010). A novel ncRNA gene from mouse chromosome 5 trans-splices with Dmrt1 on chromosome 19. Biochem. Biophys. Res. Commun. 400, 696–700. doi: 10.1016/j.bbrc.2010.08.130

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Lin, Y., He, Y., Wang, H., Chen, S., Li, Z., et al. (2020). Deletion of lncRNA5512 has no effect on spermatogenesis and reproduction in mice. Reprod. Fertil. Dev. 32, 706–713. doi: 10.1071/RD19246

PubMed Abstract | CrossRef Full Text | Google Scholar