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

Manon Chadourne, Elodie Poumerol, Luc Jouneau, Bruno Passet, Johan

Castille, Eli Sellem, Eric Pailhoux, Béatrice Mandon-Pépin

► To cite this version:

Manon Chadourne, Elodie Poumerol, Luc Jouneau, Bruno Passet, Johan Castille, et al.. Structural and Functional Characterization of a Testicular Long Non-coding RNA (4930463O16Rik) Identified in the Meiotic Arrest of the Mouse Topaz1-/- Testes. Frontiers in Cell and Developmental Biology, 2021, 9, 10.3389/fcell.2021.700290. [https://doi.org/10.32988/iv2](https://doi.org/10.32988/iv2)

HAL Id: hal-03298848

https://hal.science/hal-03298848v2

Submitted on 31 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Structural and Functional Characterization of a Testicular Long Non-coding RNA (4930463016Rik) Identified in the Meiotic Arrest of the Mouse Topaz1<sup>-/-</sup> Testes

Manon Chadourne¹, Elodie Poumerol¹, Luc Jouneau¹, Bruno Passet², Johan Castille², Eli Sellem³, Eric Pailhoux¹ and Béatrice Mandon-Pépin¹*

¹ UVSQ, INRAE, BREED, Université Paris-Saclay, Jouy-en-Josas, France, ² INRAE, AgroParisTech, GABI, Université Paris-Saclay, Jouy-en-Josas, France, ³ R&D Department, ALLICE, Paris, France

OPEN ACCESS

**Edited by:**
Serge Nef,
Université de Genève, Switzerland

**Reviewed by:**
Atsushi P. Kimura,
Hokkaido University, Japan
Mary Ann Handel,
Jackson Laboratory, United States

*Correspondence:
Béatrice Mandon-Pépin
beatrice.mandon-pepin@inrae.fr

**Specialty section:**
This article was submitted to Molecular and Cellular Reproduction, a section of the journal Frontiers in Cell and Developmental Biology

**Received:** 25 April 2021
**Accepted:** 14 June 2021
**Published:** 01 July 2021

**Citation:**
Chadourne M, Poumerol E, Jouneau L, Passet B, Castille J, Sellem E, Pailhoux E and Mandon-Pépin B (2021) Structural and Functional Characterization of a Testicular Long Non-coding RNA (4930463016Rik) Identified in the Meiotic Arrest of the Mouse Topaz1<sup>-/-</sup> Testes. Front. Cell Dev. Biol. 9:700290. doi: 10.3389/fcell.2021.700290

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*<sup>-/-</sup> 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, *4939463016Rik*, did not alter fertility although sperm parameters were disturbed and sperm concentration fell. The transcriptome of P18-*4939463016Rik*<sup>-/-</sup> 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 *4930463016Rik* 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 (*4930463016Rik*) that is key for both sperm production and motility.

**Keywords:** TOPAZ1, spermatogenesis, meiosis, germ cells, mouse, lncRNAs

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*<sup>-/-</sup>) 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*<sup>-/-</sup> 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*<sup>-/-</sup> 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 Jouyen-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<fnref n="1" /> 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<fnref n="2" />,<fnref n="3" /> 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™ 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)].

<footnote n="1">¹http://crispor.tefor.net/</footnote>

<footnote n="2">²https://www6.jouy.inrae.fr/mima2/</footnote>

<footnote n="3">³https://doi.org/10.15454/1.5572348210007727E12</footnote>
The primer sequences used for RT-qPCR are shown in **Supplementary Table 1**.

## RNA-Sequencing

Total RNA was extracted from WT, *Topaz1*<sup>-/-</sup> or *1630463O16Rik*<sup>-/-</sup> 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<fnref n="4" />, BioProject ID PRJNA698440.

## Biotype Determination of DEGs

Data available on the NCBI, MGI<fnref n="5" /> and Ensembl<fnref n="6" /> websites were used simultaneously to determine the DEG biotypes. For this purpose, information on the mouse genome was obtained by ftp from NCBI<fnref n="7" />; the annotation BioMart file from Ensembl<fnref n="8" /> (Ensembl genes 95, Mouse genes GRCm28.p6) and feature types from MGI<fnref n="9" /> (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<fnref n="11" />.

## Sperm Analysis

Evaluations of the concentrations and motility of sperm in WT and *4930463O16Rik*<sup>-/-</sup> 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*<sup>-/-</sup> 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

<footnote n="4">⁴https://www.ncbi.nlm.nih.gov/bioproject/PRJNA698440</footnote>

<footnote n="5">⁵http://www.informatics.jax.org</footnote>

<footnote n="6">⁶https://www.ensembl.org/</footnote>

<footnote n="7">⁷ftp://ftp.ncbi.nih.gov/gene/DATA/GENE_INFO/Mammalia/Mus_musculus.gene_info.gz</footnote>

<footnote n="8">⁸http://www.ensembl.org/biomart/martview</footnote>

<footnote n="9">⁹http://www.informatics.jax.org/marker/</footnote>

<footnote n="10">¹⁰https://david.ncifcrf.gov/</footnote>

<footnote n="11">¹¹https://maayanlab.cloud/Enrichr/</footnote>
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**).

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
<figure><img src="image_3.png" /></figure>

(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-Proasper 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<fnref n="12" /> (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;<fnref n="13" />), revealed little or no expression of these three genes in *Topaz1*⁻/⁻ testes (Supplementary Figures 6A-C).

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-Proasper 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.

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).

## 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.

<footnote n="12">¹²https://rgv.genouest.org/</footnote>

<footnote n="13">¹³http://software.broadinstitute.org/software/igv/</footnote>
<figure><img src="image_4.png" /></figure>

**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.

## 4930463O16Rik<sup>-/-</sup> 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<sup>-/-</sup> 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<sup>-/-</sup> 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).

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<sup>-/-</sup> Mouse Testes

Transcriptomic RNA-seq analyses were performed in WT and 4930463O16Rik<sup>-/-</sup> mouse testes at two developmental stages; i.e., at P16 and 18, as in the Topaz1<sup>-/-</sup> mouse line.
<figure><img src="image_5.png" /></figure>

**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.

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.

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
<figure><img src="image_6.png" /></figure>

**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 *4930463O16Rik*⁻/⁻ 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 *4930463O16Rik*⁻/⁻ 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).

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 the*Topaz1* 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
<figure><img src="image_7.png" /></figure>

**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.

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
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**FIGURE 6** | Evaluation of sperm parameters. Comparison of sperm-specific parameters from WT (circle, *n* = 11) and 4930463*O16Rik*⁻/⁻ (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.
**A**
WT vs 4930463O16Rik⁻/⁻

258 DEGs

<table><thead><tr><th>258 DEGs</th></tr></thead><tbody><tr><td>199 down-<br>regulated</td><td>59 up-<br>regulated</td></tr></tbody></table>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

**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.

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 I 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*<sup>-/-</sup> 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*<sup>-/-</sup> 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*<sup>-/-</sup> 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.

## 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* dihydrodi Nicolate 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*<sup>-/-</sup> 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*<sup>-/-</sup>). A Kruskal-Wallis statistical test was performed (*p* < 0.05).

**Supplementary Figure 4** | Reprogenomic data on the dynamic expression of *4930463016Rik*. The dynamic expression of *4930463016Rik* in five different tissues from male and female adult mice (A), and in embryonic primordial germ cells and adult male germ cells (B). *4930463016Rik* 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 *4930463016Rik* (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*<sup>-/-</sup> 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 *4930463016Rik* (A), *Gm21269* (B), and *492151H07Rik* (C) transcripts (red) from Ensembl data (GRCm38) is shown at the top. *4930463016Rik* and *4921513H07Rik* gene transcriptions overlap in 3' or 5', respectively.

**Supplementary Figure 7** | Expression of *Gm21269*, *4930463016Rik*, and *492151H07Rik* mRNAs in testes from 5 days to adulthood. Quantitative RT-PCR analysis of *Gm21269*, *4930463016Rik* and *492151H07Rik* gene expressions at different developmental stages in WT (blue) and *Topaz1*<sup>-/-</sup> (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 *4930463016Rik* (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) *4930463016Rik*, (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 IncRNAs (red) was evident. Scale bar = 20 µm.

**Supplementary Figure 10** | Validation of several DEGs by RT-qPCR (RNA-seq *4930463016Rik*−/− 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: *4930463016Rik*−/−). 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 *4930463016Rik*−/− testes compared to WT. List of deregulated genes in *4930463016Rik* 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 *4930463016Rik*−/− 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 *4930463016Rik*−/− testes.

## 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

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

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

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

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

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

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

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.

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

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/J.B978-0-12-394309-5.00005-5

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

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

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

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

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

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

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/JCell.2011.09.028

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/efLife.04530

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

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

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

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

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

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

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

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

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

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
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

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

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

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

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

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.rna29624

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

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

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

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

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

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

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

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

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/biolreio/ab007

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

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

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

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

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

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

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

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

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.

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

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

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

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.26433.120

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

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

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

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

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

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

Luangpraseuth-Proser, A., Lesueur, E., Jouveau, 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

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

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

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.16084581113

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.05045

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

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

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

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

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

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

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

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

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

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

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

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

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/RD-19246

**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.

Copyright © 2021 Chadourne, Poumerol, Jouveau, Passet, Castille, Sellem, Pailhoux and Mandon-Pépin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.