# A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka

Amaury Herpin, Cornelia Schmidt, Susanne Kneitz, Clara Gobe, Martina Regensburger, Aurélie Le Cam, Jérôme Montfort, Mateus C. Adolfi, Christina Lillesaar, Dagmar Wilhelm, et al.

► To cite this version:

Amaury Herpin, Cornelia Schmidt, Susanne Kneitz, Clara Gobe, Martina Regensburger, et al. A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka. PLoS Biology, 2019, 17 (4), pp.1-29. 10.1371/journal.pbio.3000185. [https://doi.org/10.1371/journal.pbio.3000185](https://doi.org/10.1371/journal.pbio.3000185)

HAL Id: hal-02093379

https://hal.science/hal-02093379

Submitted on 8 Apr 2019

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

# A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka

Amaury Herpin<sup>1,2*</sup>, Cornelia Schmidt<sup>2</sup>, Susanne Kneitz<sup>2</sup>, Clara Gobé<sup>1,3</sup>, Martina Regensburger<sup>2</sup>, Aurélie Le Cam<sup>1</sup>, Jérôme Montfort<sup>1</sup>, Mateus C. Adolfi<sup>2</sup>, Christina Lillesaar<sup>2</sup>, Dagmar Wilhelm<sup>4</sup>, Michael Kraeusling<sup>2†</sup>, Brigitte Mourot<sup>1,4</sup>, Béatrice Porcon<sup>1</sup>, Maëlle Pannetier<sup>1,3</sup>, Eric Pailhoux<sup>3</sup>, Laurence Ettwiller<sup>1,5</sup>, Dirk Dolle<sup>5</sup>, Yann Guiguen<sup>1,6</sup>, Manfred Schartl<sup>2,6,7</sup>

1 INRA, UR1037 Fish Physiology and Genomics, Rennes, France, 2 University of Wuerzburg, Physiological Chemistry, Biocenter, Wuerzburg, Germany, 3 UMR BDR, INRA, ENVA, Université Paris Saclay, Jouy en Josas, France, 4 University of Melbourne, Department of Anatomy & Neuroscience, Parkville, Victoria, Australia, 5 University of Heidelberg, Centre for Organismal Studies (COS), Department of Developmental Biology, Heidelberg, Germany, 6 Comprehensive Cancer Center Mainfranken, University Hospital, Wuerzburg, Germany, 7 Hagler Institute for Advanced Study and Department of Biology, Texas A&M University, College Station, Texas, United States of America

¶Authorship confirmed by corresponding author.

* amaury.herpin@inra.fr

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

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

OPEN ACCESS

**Citation:** Herpin A, Schmidt C, Kneitz S, Gobé C, Regensburger M, Le Cam A, et al. (2019) A novel evolutionary conserved mechanism of RNA stability regulates synexpression of primordial germ cell-specific genes prior to the sex-determination stage in medaka. PLoS Biol 17(4): e3000185. https://doi.org/10.1371/journal.pbio.3000185

**Academic Editor:** Yukiko M Yamashita, University of Michigan, UNITED STATES

**Received:** September 28, 2018

**Accepted:** March 4, 2019

**Published:** April 4, 2019

**Copyright:** © 2019 Herpin et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the Deutsche Forschungsgemeinschaft by a grant (SCHA 408/12-1; HE 7135/2-1) to AH and MS as well as Crédits Incitatifs Phase 2015/Emergence to AH. AH was additionally funded by the project AquaCRISPR (ANR-16-COFA-0004-01). The

## Abstract

Dmrt1 is a highly conserved transcription factor, which is critically involved in regulation of gonad development of vertebrates. In medaka, a duplicate of dmrt1—acting as master sex-determining gene—has a tightly timely and spatially controlled gonadal expression pattern. In addition to transcriptional regulation, a sequence motif in the 3' UTR (D3U-box) mediates transcript stability of dmrt1 mRNAs from medaka and other vertebrates. We show here that in medaka, two RNA-binding proteins with antagonizing properties target this D3U-box, promoting either RNA stabilization in germ cells or degradation in the soma. The D3U-box is also conserved in other germ-cell transcripts, making them responsive to the same RNA binding proteins. The evolutionary conservation of the D3U-box motif within *dmrt1* genes of metazoans—together with preserved expression patterns of the targeting RNA binding proteins in subsets of germ cells—suggest that this new mechanism for controlling RNA stability is not restricted to fishes but might also apply to other vertebrates.

## Author summary

The development of the gonads in vertebrates is mainly regulated by *dmrt1*, a master sex-determining gene that has a timely and spatially controlled gonadal expression pattern. In addition to transcriptional regulation, a sequence motif located in the 3' UTR (D3U-box) mediates transcript stability of dmrt1 mRNAs. However, this regulation is complex, and
funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

**Abbreviations:** BAC, bacterial artificial chromosome; BSA, bovine serum albumin; BSF, bicoid stability factor; Carbio, Carolina Biological Supplies; celf2, CUGBP Elav-like family member 2; CRISPR/Cas9, clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9; D3U-box, dmrt1 3' UTR box; dph, days post hatching; EMSA, electrophoretic mobility shift assay; GFP, green fluorescent protein; GO, gene ontology; Lrpprc, leucine rich pentatricopeptide repeat containing; MES-1, medaka embryonic stem cells; miR-224, microRNA 224; OI-BSF, Oryzias latipes Bicoid Stability Factor; OI-CUG-BP, Oryzias latipes CUG-binding protein; ORF, open reading frame; PARN, poly(a)-specific ribonuclease; PBS, phosphate-buffered saline; PFA, paraformaldehyde; PGC, primordial germ cell; WNT4, Wnt family member 4.

the mechanisms are barely understood. Here, we analyse further dmrt1 regulation in the fish model medaka and show that two RNA-binding proteins with antagonizing properties target the D3U-box and promote either RNA stabilization in germ cells or degradation in the soma. The evolutionary conservation of the D3U-box motif within metazoan *dmrt1* genes—together with preserved expression patterns of the targeting RNA binding proteins in subsets of germ cells—suggest that this new mechanism for controlling RNA stability is not restricted to fishes but might also apply to other vertebrates.

## Introduction

The gonads of vertebrates are characterized by the intimate association of germ cells and supporting somatic cells [1–4]. The precursor cells of the soma are derived from the embryonic lateral plate mesoderm, whereas germ cells originate from the germline lineage [5–9]. To carry out their highly specialized biological functions, the somatic gonadal primordium and the germline cells together must establish timely regulated programs of gene expression [1,10,11].

The *mab-3/doublesex/dmrt1* gene orthologs are, among metazoans, the most evolutionary conserved key regulators of the earliest phases of gonad development. They control complex gene regulatory networks specifying male gonadal primordium development as well as gonadal maintenance [12–15]. Remarkably, besides being firmly anchored within the regulatory network at critical nodes, *dmrt1* genes were found to act as upstream male sex determiners in organisms as phylogenetically diverse as flatworm [16], water flea [17], frog [18], flatfish [19], birds [20] and medaka [21,22].

During the last decade, much has been learned about how *dmrt1* as the most versatile sex gene triggers and controls gonad development. In human, it is a critical dosage-sensitive sex-determining gene, such that haploinsufficiency leads to XY male-to-female sex reversal and infertility [23,24]. In mice, it is required for male gonadal differentiation of somatic and germ cells [25–27], although *Dmrt1* appears to be dispensable for primary sex determination [15]. *Dmrt1* also plays the decisive role in maintaining the cellular identity of the adult testis, most obvious from the fact that its malfunction in adult mutant mice leads to transdifferentiation of Sertoli to granulosa-like cells and feminization of a fully developed testis [27]. Consequently, the action range of *Dmrt1* is not restricted to initiation of the male gonadal phenotype during early development but also contributes to the active suppression of the female networks via repression of two ‘anti-testis’ pathways, Foxl2 and *Wnt* family member 4 (*WNT4*/β-catenin (see [28] for review).

In sex determination model fish medaka, male sex determination is implemented by a male-specific primordial germ cell (PGC) mitotic arrest due to the activity of a Y-chromosome-specific duplicate version of *dmrt1*, designated *dmrt1bY* [29]. In *dmrt1* knockout mice, germ cells fail to arrest mitosis [30]. Further work on *dmrt1* has shown it to be a transcriptional gatekeeper controlling mitosis versus meiosis decision in male germ cells [26]. Thus, *dmrt1* in mice and *dmrt1bY* in medaka appear to be regulators of germ cell proliferation.

Despite its well-characterized crucial functions for gonad development in many vertebrates, the mechanisms that regulate the complex temporal and spatial expression pattern and guarantee precise levels of *dmrt1* transcripts are only barely understood. Diverse regulatory mechanisms have been occasionally reported. Indirect transcriptional regulation of *dmrt1* upon steroid treatments has been described in several fish species (see [13] for review). Gonadal dimorphic expression of *dmrt1* has been suggested to be possibly under the control of differential CpG methylation of its promoter in two different flatfish species [31,32]. Similarly, in the
red-eared slider turtle (*Trachemys scripta*), DNA methylation dynamics accounting for *dmrt1* sexual dimorphic expression are tightly correlated with temperature [33]. In vitro transcriptional regulation assays revealed that binding sites for *Sp1*, *Egr1* [34], and *Gata4* [35] factors, which are present in the promoters of many genes, are also involved in transcriptional regulation of the rat *dmrt1* gene. And finally, evidence was presented that microRNA 224 (miR-224) promotes differentiation of mouse spermatogonial stem cells via direct targeting of *dmrt1*, decreasing its expression in testes [36]. Certainly in the context of ‘indirect’ regulation, *dmrt1* is one of the most prominent examples.

In medaka, for which a functional duplicate of the autosomal *dmrt1a* gene on the Y chromosome—*dmrt1bY*—became the master regulator of male sex determination [22,37], transcriptional rewiring was brought about by exaptation of two transposable elements, *Izanagi* and *Rex1*, co-opted to act as silencers. These turn off the somatic and the germ cell-specific expressions of the *dmrt1bY* gene [38, 39]. Thus far, two factors, *dmrt1* itself [38] and *sox5* [39], were identified, which turn off *dmrt1bY* expression after it has fulfilled its function in the early developing gonad [38, 39].

We previously identified a 11-bp sequence motif in the 3' UTR of *dmrt1bY* (D3U-box, for *dmrt1* 3' UTR box). This motif confers stability to the mRNA in the developing embryonic gonad, whereas in other tissues, the transcript is rapidly degraded [40], indicating that a post-transcriptional regulation mechanism could play a role in germline expression of *dmrt1* in medaka.

Here, we show that the dmrt1 11-nucleotide *cis*-regulatory D3U-Box motif is a target for two antagonizing RNA binding proteins, *Oryzias latipes* CUG-binding protein (Ol-cug-bp) and *Oryzias latipes* Bicoid Stability Factor (Ol-bsf)—also known as cugbp Elav-like family member (celf) and leucine rich pentatricopeptide repeat containing (lrpprc), respectively, in mammals. In *Drosophila*, the bicoid stability factor (bsf) has initially been shown to be involved in regulating the stability of *bicoid* transcripts during oogenesis through binding structures within the 3' UTR of transcripts that resemble CUG hairpins [41]. Later, bsf was also reported to have a role in regulation of early zygotic genes by binding a short consensus sequence (CAGGUA) in the 5' UTR of genes expressed in the early zygote [42]. Cug-bp is the human homolog of the *Xenopus* eden-bp, which was shown to bind to mRNAs, such as *c-mos*, that exhibit rapid deadenylation following fertilization of oocytes [43]. Previous studies of cug-bp function have focused mainly on the roles of this protein in regulating alternative splicing [44] and also on its ability to modulate translation of several mRNAs [45]. However, as cug-bp is able to functionally substitute for eden-bp to induce deadenylation in *Xenopus* oocyte extracts [46], it seems likely that it also plays a similar role in regulating poly(A) shortening in mammalian cells. Indeed, it was shown that cug-bp can interact with poly(a)-specific ribonuclease (PARN) deadenylase to promote deadenylation of its substrate RNAs [47].

We find that in medaka, the D3U-box is targeted by these two different RNA binding proteins, with Ol-cug-bp1 leading to dmrt1bY degradation unless Ol-bsf is present in germ cells. Moreover, this new mechanism of dmrt1 RNA stability appears to regulate also the abundance of other transcripts specifically expressed in PGCs.

## Results

### Specific enrichment of a conserved *cis*-regulatory motif (D3U-box) in vertebrate genomes

In an initial analysis of dmrt1 post-transcriptional regulation [40], we found that an 11-bp long *cis*-regulatory motif in the D3U-box confers transcript stability in PGCs (Fig 1A). *In vitro* evidence was obtained that the D3U-box possibly mediates its function through protein
<figure><img src="image_4.png" /><figcaption>A</figcaption></figure>

* D3U-box present in male form and spliced out in the female form.

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

"Fish-specific" D3U-box consensus motif

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

"Full-species" D3U-box consensus motif

<figure><img src="image_7.png" /><figcaption>B</figcaption></figure>

**Fig 1. Conservation of the D3U-box motif from ecdysozoans to mammals.** (A) Conservation of the D3U-box motif from *Drosophila* up to mammals. From the different D3U-box sequences among vertebrates, weight-position matrices were deduced and used for genome scans. (B) Occurrences of the D3U-box motif compared to control motifs within the 3' UTR sequences of human, mouse, medaka, or the green spotted puffer (*Tetraodon nigroviridis*). For comparison, scrambled control motifs were used. Absolute values of the occurrences depend on the depth of annotation for every genome. Underlying data for (B) can be found in S1 Data. D3U-box, dmrt1 3' UTR box.

https://doi.org/10.1371/journal.pbio.3000185.g001
binding [40]. Further on this sequence motif was found to be evolutionary conserved in *dmrt1* genes from flies to mammals (Fig 1A; [40]).

To investigate whether this motif is specific for the *dmrt1* genes or also present in other genes, we performed genome-wide searches for human, mouse, medaka, and the green spotted puffer (Fig 1B). The D3U-box was found in more than 250 3' UTRs of genes from human and mice and is also abundant in medaka and green spotted puffer (Fig 1B). Notably, also, a 'split' version of the D3U-box (see Fig 1B) displayed specific enrichment in 3' UTRs (Fig 1B), suggesting that the D3U-box motif might be articulated around two independent *cis*-regulatory sequences, hence putatively targeted by 2 different RNA binding proteins.

Medaka whole transcriptome scans (3' UTR and coding sequences) using the vertebrate D3U-box consensus motif matrix (Fig 1A: 3' UTR sequences; and S1 Fig: 3' UTR and coding sequences) resulted in several hits, including tra2, sox10, misr2, dead end, and vasa (S1 Fig). Like dmrt1 [14,29,48,49], these RNAs are critically involved in germ cell development and maintenance in medaka and many other organisms ([36,50–55] and [56,57] for review).

Furthermore, bioinformatics analyses and literature searches [45,47] revealed that the D3U-box *cis*-regulatory motif is a putative target for 2 evolutionary conserved RNA binding proteins involved in either mRNA degradation or stabilization. These 2 proteins, Ol-cug-bp (also known as CELF in mammals) and Ol-bsf (also known as LRPPRC in mammals), have been shown to specifically recognize CUG repeats and the CAGGU(AG) motif, respectively, which constitute the D3U-box (see S2 Fig for phylogeny and synteny analysis of *Ol-bsf* and *Ol-cug-bps*).

## Ol-BSF and Ol-CUG-BP1 specifically bind to the D3U-box motif

To confirm our bioinformatics prediction, Ol-bsf and Ol-cug-bp1 and 2 proteins were subjected to electrophoretic mobility shift assay (EMSA) using the D3U-box motif as target and different competitors (Fig 2). The in vitro-translated proteins (Fig 2A and 2F) were assayed for binding with radioactively labelled RNA probes. Using the D3U-box motif, mobility shifts were detected for the 2 proteins tested: Ol-bsf (Fig 2B to 2E) and Ol-cug-bp1 (Fig 2G to 2I), indicating that Ol-bsf as well as Ol-cug-bp1 are, in principle, able to bind the D3U-box in vitro. Binding specificities were confirmed by competition of the medaka D3U-box motif for Ol-bsf or Ol-cug-bp1 interactions with either a scrambled D3U-box-derived motif (Fig 2B and 2I) or a minus CUG repeat motif competitor (Fig 2H). The absence of any significant interference with the D3U-box binding indicated the specificity of the observed interactions (Fig 2B, 2I and 2H). Furthermore, competition experiments between radioactively and nonradioactively labelled D3U-boxes resulted in progressive loss of the apparent shifts (Fig 2C and 2G). Notably, a clearly visible shift was also observed when using the *Drosophila* D3U-box sequence together with the medaka Ol-bsf protein (Fig 2D). Altogether, these experiments suggest that the D3U-box is a preferential target for Ol-bsf and Ol-cug-bp1 binding. Of note, performing the very same set of experiments together with the Ol-cug-bp2 protein did not result in any convincing evidence for specific binding to the D3U-box.

## Ol-bsf and Ol-cug-bps antagonistically regulate the expression and stability of reporter constructs harbouring the D3U-box motif and of dmrt1bY transcripts

To monitor a possible effect of Ol-bsf on regulation of the male sex-determination gene in medaka, we generated a dmrt1bY reporter line by introducing the green fluorescent protein (GFP) open reading frame (ORF) fused to the dmrt1bY 3'UTR (including the D3U-box) into exon 1 of a bacterial artificial chromosome (BAC) clone containing the *dmrt1bY* gene and
<figure><img src="image_8.png" /><figcaption>**Fig 2.** EMSA of in vitro-translated Ol-bsf, Ol-cug-bp1, and Ol-cug-bp2 proteins indicating interactions with the medaka D3U-box target sequence. (A and F) Recombinant protein production of FLAG-tagged versions of Ol-bsf and Ol-cug-bp1 and western blotting detection with an anti-FLAG antibody. (B-E and G-I) EMSA using the recombinant medaka Ol-bsf and Ol-cug-bp1 proteins (in vitro translated) and either the medaka D3U-box sequence (B, C, E, G, H, I), or the *Drosophila* D3U-box sequence (D) as radioactively labelled RNA probes. Increasing amounts of unlabelled probes (scrambled [B, I], medaka [C, G, K], *Drosophila* [D, E], minus CUG [H, L] target sequences) were used as competitors. (B-E and G-I) Apparent shifts are observed for the labelled RNA target probes with the Ol-bsf and Ol-cug-bp1 RNA-binding proteins, likely indicating direct interactions. Ratios of RNA probe to RNA competitor: (B, lane 2: 1/5; C, lane 2: 1/1, lane 3: 1/2, lane 4: 1/4; D, lane 2: 1/1, lane 3: 1/2, lane 4: 1/5, lane 5: 1/10; E, lane 2: 1/1, lane 3: 1/2, lane 4: 1/5, lane 5: 1/10; G, lane 2: 1/1, lane 3: 1/2, lane 4: 1/5; H, lane 2: 1/1, lane 3: 1/2, lane 4: 1/5; I, lane 2: 1/1, lane 3: 1/2, lane 4: 1/5, lane 5: 1/10). D3U-box, dmrt1 3' UTR box; EMSA, electrophoretic mobility shift assay; lrpprc, leucine rich pentatricopeptide repeat containing; Ol-bsf, Oryzias latipes Bicoid Stability Factor; Ol-cug-bp, Oryzias latipes CUG-binding protein.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3000185.g002

flanking regions (Fig 3A and 3B). The recombined BAC was then used for establishing a stable transgenic line in which GFP expression most reliably indicates endogenous *dmrt1bY* expression [48,58,59]. Expression of dmrt1bY is highly dynamic during primordial gonad formation, progressively switching from germ cell expression only between stages 26 (1.25 dpf) and 29 (3.1 dpf) to an exclusive somatic expression from stages 33/34 (5 dpf) up to hatching (9 dpf) when the gonad is formed [39,59]. After injection of the Ol-bsf morpholino (see S3 Fig for validation of the morpholino), we found a significant reduction of GFP expression (Fig 3B) and, after Ol-bsf overexpression, a strong increase of the reporter at both mRNA and protein (fluorescence) levels (Fig 3A and 3B).

Next, to obtain a more physiological readout of the role(s) of Ol-bsf and Ol-cug-bps for RNA stability in vivo, the relative abundances of endogenous dmrt1bY transcripts were monitored after modulation of Ol-bsf/Ol-cug-bp1/Ol-cug-bp2 expression in medaka embryos (Fig 3C and 3D). First, we checked for changes in dmrt1bY transcript levels after overexpression of the 2 medaka Ol-cug-bp ohnologs (Ol-cug-bp1 and Ol-cug-bp2, Fig 3C). This resulted in
<figure><img src="image_9.png" /><figcaption>**Fig 3.** In vivo modulation of expression and mRNA stability of reporters and *dmrt1bY* transcripts via the D3U-box. (A and B) Effect of the modulation of the Ol-bsf expression on a GFP reporter for *dmrt1bY* expression. Injections of either Ol-bsf morpholino or capped mRNA result in the subsequent modulations of expression of a GFP:(*dmrt1bY* 3' UTR) reporter/sensor construct in a transgenic line (BAC *Dmrt1bY:GFP:3'UTR*) as monitored by fluorescence (A) and RNA quantifications (B). Because endogenous GFP is lowly expressed already before experimental down-regulation, the more sensitive qPCR method was employed to monitor the morpholino effect (B). Dataset results of 4 independent reverse transcription reactions, resulting in 4 different batches of eggs obtained from different couples. Statistical significance has been assessed by mean of the Wilcoxon-Mann-Whitney test (N = 4). (C and D) Real-time qPCR determination of *dmrt1Y* abundances after Ol-cug-bp1 and Ol-cug-bp2 (C) or Ol-bsf (D) modulation of expression in embryos at stage 18 after injection of either capped mRNAs with the full Ol-bsf/Ol-cug-bps ORFs or a splice morpholino targeting the exon2-intron2 splice junction of *Ol-bsf* (S3 Fig). Dataset results of the analysis of 4 batches of eggs injected with either Ol-cug-bp1, Ol-cug-bp2, Ol-bsf, or BSF-MO. Statistical significance has been assessed by mean of the Wilcoxon-Mann-Whitney test (N = 4). \*p ≤ 0.05, \*\*p ≤ 0.01. (E and F) Real-time qPCR determination of the kinetics of RNA stability for D3U-box-containing reporters in presence or in absence of either Ol-cug-bp1 or Ol-bsf coding mRNAs. Datasets are results of 3 independent reverse transcription reactions, resulting from three different batches of injected eggs. Underlying data for (B, C, D, and F) can be found in S1 Data. BAC, bacterial artificial chromosome; D3U-box, *dmrt1* 3' UTR box; GFP, green fluorescent protein; MO, Morpholino; ns, nonsignificant; Ol-BSF, *Oryzias latipes* Bicoid Stability Factor; Ol-CUG-BP, *Oryzias latipes* CUG-binding protein; ORF, open reading frame; qPCR, quantitative PCR; W.T., wild type.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3000185.g003

decreased *dmrt1bY* mRNA (Fig 3C). Second, the relative abundances of *dmrt1bY* transcripts were recorded after either overexpression or morpholino knockdown of *Ol-bsf* in medaka embryos (Fig 3D). It revealed that higher *Ol-bsf* expression correlates with an increased abundance, while lowering *Ol-bsf* expression resulted in a reduction of *dmrt1bY* transcripts
(Fig 3D). In vivo, D3U-box-induced modulation of RNA stability was further investigated in embryos injected with either control RNAs or RNAs harbouring the D3U-box (Fig 3E). Ratios between control and D3U-box-containing mRNAs were then quantified in absence or in presence of either l-cug-bp1 or Ol-bsf mRNAs (Fig 3F). It revealed that, over time, overexpression of Ol-bsf correlates with an increased stabilization of the D3U-box-containing mRNAs, while overexpression of Ol-cug-bp1 correlates with a decreased stabilization of the D3U-box-containing mRNAs (Fig 3F).

## Expression of Ol-cug-bp1, Ol-cug-bp2, and Ol-bsf during embryonic development and in adult tissues

During the embryonic developmental period, both Ol-cug-bp ohnologs display complementary patterns of expression. Ol-cug-bp1 is expressed at early stages and Ol-cug-bp2 only later when dmrt1bY transcripts appear in the germ cells around stage 25 (S4A and S4C Fig). In adult tissues, Ol-cug-bp1 and Ol-cug-bp2 are expressed at high levels in brain and gonads (S4B and S4D Fig), while Ol-cug-bp2 is additionally expressed in eyes, muscles, and skin (S4D Fig). Of note, both ohnologs are always more expressed in testes compared to ovaries (S4B and S4D Fig). Expression profiling of Ol-bsf revealed that in adult medaka, it is ubiquitously expressed in all adult tissues, with particular high levels in gonads of both sexes (S4E and S4F Fig).

## Ol-bsf is specifically expressed in the germ cells during early gonad primordium formation with correlated levels of expression between Ol-bsf and dmrt1bY at hatching stage

For bioimaging analyses of protein localization over time, we used expression reporter lines for vasa [11] and Ol-bsf (see Materials and methods), respectively. During embryonic development, we noted a distinct spatially and temporarily restricted expression pattern (Fig 4). From fertilization up to stages 16/17, Ol-bsf is expressed throughout the embryo (Fig 4A). Of note, cell transfection of a tagged version of OL-bsf and subsequent immunohistochemistry revealed that bsf protein is localized in the cytoplasm (insert in Fig 4A). From stage 25 onwards—when germ cells line up on both sides of the embryo within the lateral plate mesoderm—progressively, Ol-bsf expression becomes restricted to the PGCs (Fig 4B to 4G) where it is co-expressed with Ol-vas, a specific germ cell marker in medaka [60]. During the following developmental stages (stages 33/34), expression heterogeneity for Ol-bsf between germ cells became obvious (Figs 4H to 4J and 5C and 5D). This heterogeneity was particularly apparent between 4 to 10 days post hatching (dph), when Ol-bsf is higher expressed at the tips of the forming gonads (curly brackets in Fig 4K to 4P and square brackets in Fig 5C and 5D). In summary, Ol-bsf has a highly dynamic expression pattern, switching from an early somatic to a progressively restricted germ cell expression. Within the germ cell pool, the levels of expression show a significant heterogeneity.

Furthermore, dmrt1bY expression shows heterogeneity between individual germ cells (Fig 4Q and 4R). Being also higher expressed at the tip of the primordial gonads, the expression levels of Ol-bsf and dmrt1bY show a clear positive correlation (Fig 4Q and 4R).

## Expression of Ol-bsf in the germ-line stem cells of adult gonads

Given the high abundance of Ol-bsf transcripts in adult gonads detected by qPCR (S4E and S4F Fig), we next monitored expression of Ol-bsf in fully mature gonads of both sexes at cellular resolution (Fig 6). In adult testes, Ol-bsf fluorescence is restricted to two distinct
<figure><img src="image_10.png" /><figcaption>Fig 4. Embryonic expression of medaka Ol-bsf and correlated expression dynamics together with dmrt1bY during gonadal primordium formation. (A to R) Medaka Ol-bsf expression investigated by fluorescence using a transgenic reporter line for which a 732-bp *Ol-bsf* promoter fragment (up to the next upstream gene) drives the expression of the mCherry. Medaka Ol-vas and dmrt1bY expressions were investigated by fluorescence using previously described transgenic reporter lines ([11] and [61], respectively). (A) During neurulation (stages 16/17), Ol-bsf is ubiquitously expressed. (A insert) Cytoplasmic localisation of Ol-bsf after transient transfection of a 3-times FLAG version of Ol-bsf in</figcaption></figure>
a medaka embryonic stem cell line and immunohistochemistry. (B to P) From stage 25 on, Ol-bsf expression is restricted to the germ cells and colocalize with the germ cell marker Ol-vas. (K to P) Four-dph expression of Ol-bsf is apparently heterogeneous within the population of PGCs while dmrt1bY is now expressed in the surrounding somatic cells of the primordial gonad (K to M). (Q and R) By stage 29, Ol-bsf and dmrt1bY are both expressed in the germ cells. Variations within the respective levels of Ol-bsf and dmrt1bY are clearly observable amongst different germ cells although always correlated between the 2 fluorescences (curly brackets compared to square brackets). dph, days post hatching; MES-1, medaka embryonic stem cells; Ol-BSF, Oryzias latipes Bicoid Stability Factor; olf1, *Oryzias latipes* fibroblast-1; PGC, primordial germ cell.

https://doi.org/10.1371/journal.pbio.3000185.g004

subpopulations of germ cells, which are also positive for Ol-vas (Fig 6A to 6C). Diagnosed by condensed nuclear morphology and size [62] and localization [39], these first Ol-bsf-positive cells represent the earliest step of germ cell differentiation while another subpopulation of more mature germ cells is also observed (Fig 6C and 6F). In ovaries, Ol-bsf fluorescence is restricted to the germinal cradle [63] located in the interwoven threadlike ovarian cords at the periphery of the ovary (Fig 6G to 6N). These Ol-bsf-positive cells, representing the smaller-size subpopulation of Ol-vas fluorescent cells, are assigned to germline stem cells and early dividing germ cell lineage [63] (Fig 6O to 6Q).

Furthermore, in mice both bsf (lrpprc) and cug-bp1 (celf1) are expressed in the germ cells within the testis cords and germ cells of the ovary (S5A and S5L Fig). Immunohistochemistry revealed that bsf/lrpprc is expressed only in a subpopulation of germ cells in mice (S5M to S5R Fig).

## Ol-bsf mutant fish display gonadal phenotypes

To delineate the physiological role of Ol-bsf during gonad formation and maintenance, we generated medaka *Ol-bsf* knockout lines after genome editing using the clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) technology (S6 Fig). Homozygous *Ol-bsf* knockout larvae display reduced swimming (S7 Fig) and die within the first 2 weeks after hatching. Heterozygote mutant fish develop normally and produce mature gametes. Histological analysis of gonads of heterozygous mutants of both sexes revealed, however, that ovaries had an accumulation of small-sized oocytes compared to wild type (Fig 7A to 7D, S8 and S9 Figs), whereas testes of heterozygous mutants exhibited reduced amounts of spermatogonial stem cells, with germ cells in advanced stages of spermatogenesis located close to the periphery of the organ (Fig 7E and 7F and S10 Fig). In both sexes, *Ol-bsf* heterozygote mutants present an increase of germ cell committing to gametogenesis. While adult mutant males do not display any observable bias of fertility, adult mutant females have reduced egg production together with lower fertilization rates (S11 Fig).

## Ol-bsf and Ol-cug-bps antagonistically orchestrate expression levels of several germ cell transcripts selectively harbouring the D3U-box

Our matrix scan bioinformatic analysis had revealed the presence of the D3U-box either in 5', 3', or coding regions of several germ cell transcripts (S1 Fig). To find out whether Ol-bsf and Ol-cug-bp might regulate stability of these RNAs in a similar way, like for dmrt1bY, during gonadal development the relative abundances of these transcripts harbouring the box motif (tra2, sox10, misr2, dead end, and vasa; see S1 Fig) were first monitored after either overexpression or morpholino knockdown of Ol-bsf in medaka embryos (Fig 5A1 to 5A10). It revealed that for the majority of these transcripts (sox10, misr2, tra2, and vasa), higher Ol-bsf expression correlated with an increased abundance (Fig 5A1 to 5A5 odd numbers and 5A19), while reduced Ol-bsf expression resulted in lower levels of most of these transcripts (Fig 5A2 to 5A6 even numbers). As exceptions, dead end (Fig 5A7 and 5A8) transcript abundance
Modulation of
Ol-BSF / LRPPRC

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

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

Modulation of
Ol-CUG-BP1 & Ol-CUG-BP2

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

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

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

Fig 5. In vivo transcriptional regulation of putative gene candidates harbouring the D3U-box motif in their transcripts after modulation of the Ol-bsf or Ol-cug-bp expressions and correlation between Ol-bsf and Ol-vas expression in germ cells. (A1 to A10) Real-time qPCR determination of expression of putative Ol-BSF-regulated transcripts in embryos at stage 18 after injection of either capped mRNAs with the full Ol-bsf ORF (left odd-numbered panels) or a splice morpholino targeting the exon2-intron2 splice junction (S3 Fig) of the *Ol-bsf* gene (right even-numbered panels) at one cell stage. Ol-bsf overexpression results in an apparent increase of most RNA target candidates, while reduction of the Ol-bsf expression correlates with a reduction of expression. Dataset results of the analysis of 4 batches of eggs injected with either Ol-bsf or Ol-bsf-MO. Statistical significance has been assessed by mean of the Wilcoxon-Mann-Whitney test (N = 4). (B1 to B10) Real-
time qPCR determination of expression of putative Ol-cug-bp-regulated transcripts in embryos at stage 18 after overexpression of either Ol-cug-bp1 or Ol-cug-bp2. Results are presented as normalized expressions compared to wild type using 3 different housekeeping genes. Dataset results of the analysis of 2 batches of eggs injected with either Ol-cug-bp1 or Ol-cug-bp2. Statistical significance has been assessed by mean of the Wilcoxon-Mann-Whitney test (N = 4). (C and D) Germ cell expression of Ol-bsf in comparison to Ol-vas in double transgenic fluorescent reporter lines. After hatching expression levels of Ol-bsf and Ol-vas become more heterogeneous among germ cells (brackets at the tips of the primordial gonad) although tightly correlated between each other within individual germ cells (arrows; ‘*’ = autofluorescent pigment cells). (E) In vivo modulation of PGC number after Ol-bsf morpholino injection. An apparent increase in PGC number is observed after negative regulation of Ol-bsf expression in early embryos (stage 18). Statistical significance has been assessed by means of the Wilcoxon-Mann-Whitney test (N = 34 and 28 for wild-type and MO-bsf-injected embryos, respectively). (F) Model for D3U-box-mediated mRNA regulation. Overall and in addition to a cytoplasmic localization of Ol-bsf (Fig 4A), Ol-bsf and O-cug-bps might mutually antagonize toward the access to the D3U-box, resulting in either stabilisation (more Ol-bsf binding) or destabilisation (more Ol-cug-bp binding) of the transcripts harbouring the D3U-box. \*p ≤ 0.05; \*\*p ≤ 0.01. Underlying data for (A, B, and E) can be found in S1 Data. D3U-box, *dmrt1 3' UTR box*; MO, Morpholino; ns, nonsignificant; Ol-BSF, Oryzias latipes Bicoid Stability Factor; Ol-CUG-BP, Oryzias latipes CUG-binding protein; PGC, primordial germ cell; qPCR, quantitative PCR.

https://doi.org/10.1371/journal.pbio.3000185.g005

showed either no or an opposite pattern following Ol-bsf expression modulation, while vasa (Fig 5A10) was unaffected in the Ol-bsf-morpholino-treated embryos.

We next checked for the relative abundances of the very same transcripts harbouring the D3U-box (tra2, sox10, misr2, dead end, and vasa) after overexpression of the 2 medaka Ol-cug-bp ohnologs (Ol-cug-bp1 and Ol-cug-bp2, Fig 5B1 to 5B10). With the exception of vasa (Fig 5B9 and 5B10), the majority of the transcripts analysed had lower abundance, whereas medaka Ol-cug-bp ohnologs were overexpressed (Fig 5B1 to 5B8).

In vivo, an apparent correlation between Ol-bsf levels and Ol-vas expression could be visualized in the germ cells of the forming gonadal primordium using fluorescent reporter lines (Fig 5C and 5D). Reduced levels of Ol-bsf expression after morpholino injection led to a significant increase in PGC number at stage 18, (Fig 5E). In line with this observation, microarray data comparing *Ol-bsf*⁻(+)⁻deficient testes (displaying reduced levels of Ol-bsf; S12A Fig) to wild type revealed a general up-regulation of genes involved in germ cell proliferation or differentiation. A significant proportion (10.1%) of the down-regulated genes codes for proteins localized in the mitochondria. Finally, gene ontology (GO) term analysis revealed that, in mutant testes partially depleted for the *Ol-bsf* gene, rRNA processing is particularly affected (S12A Fig).

## Discussion

The expression of most genes is dynamically tightly regulated, temporally and spatially. Such regulations occur at multiple steps, including transcription, splicing, mRNA transport, mRNA stability, translation, protein stability, and post-transcriptional modifications [65,66]. While the importance of complex post-transcriptional regulations—like in the case of nanos, oskar, or bicoid, e.g. [61,67–69]—has been mainly identified through genetic approaches for the development of the germline or oocyte, respectively, such approaches have proven to be much less valuable for finding the expected regulatory proteins that bind specifically to these mRNAs. In medaka fish, expression of the master sex determiner *dmrt1bY* mRNA is very dynamic, occurring first in the PGCs prior to morphological somatic sex differentiation and then quickly switches to an exclusive Sertoli cell localisation [59,70]. Importantly, *dmrt1bY* is expressed in PGCs of male embryos much before its expression in the pre-Sertoli cells at the sex determination stage [59]. This early PGC expression is necessary for the later onset of *dmrt1bY* expression in the pre-Sertoli cells at the sex-determination stage of male development [59]. There, the level of *dmrt1bY* mRNA needs to reach a certain threshold to exert the sex-determining function [21]. This suggested that medaka germ cells exhibit sexually different
<figure><img src="image_16.png" /><figcaption>**Fig 6.** Expression of Ol-bsf in adult gonads. (A to F) In adult testes, Ol-bsf fluorescence (same transgenic lines as in Fig 3) is restricted to the cells located within the lobules and colocalize with a subpopulation of Ol-vas-positive cells. (G to Q) Within the ovary, Ol-bsf fluorescence is restricted to the germinal cradle located in the interwoven threadlike ovarian cords at the periphery of the ovary. The Ol-bsf-positive cells represent the smaller-size subpopulation of Ol-vas fluorescent cells (G to N) as well as early dividing germ cells (O to Q). Arrowheads indicate the germinal cradles (G to N), and brackets indicate early dividing mitotic or meiotic germ cell (O to Q). Ol-BSF, *Oryzias latipes* Bicoid Stability Factor.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3000185.g006
<figure><img src="image_17.png" /><figcaption>**Fig 7.** Ovarian and testicular phenotypes of the *Ol-bsf* mutant fish. Heterozygote mutant fish develop normally. (A to D) In-depth morphological inspection of heterozygote mutant ovaries discloses a significant accumulation of small-sized oocytes compared to wild type (dotted lines in [C] and red arrows in [B] and S8 and S9 Figs). (D) Size distribution of the oocytes in 2 wild-type and 2 *Ol-bsf*(+/-) adult ovaries (see also S8 Fig for details). (E and F) Heterozygote mutant testes exhibit a decreased number of spermatogonia with accumulation of type 2 spermatocytes, spermatids, and sperm within the most external layers of the seminiferous epithelium (see also S10 Fig). The different stages of spermatogenesis were determined according to [64]. Each gonad was sectioned through the mid-sagittal plane. Underlying data for (D) can be found in S1 Data. Ol-BSF, *Oryzias latipes* Bicoid Stability Factor; Sc1, type 1 spermatocytes; Sc2, type 2 spermatocytes; Sg, spermatogonia; Sm, sperm; St, spermatids.</figcaption></figure>

https://doi.org/10.1371/journal.pbio.3000185.g007
characters before the formation of the somatic gonadal primordium depending on dynamic
and tightly timely regulated mechanisms of post-transcriptional regulations [59,70].

We have identified in medaka a *cis*-regulatory 11-bp motif in the 3' UTR of *dmrt1bY* called D3U-box. This motif confers stability to the *dmrt1bY* mRNA in germ cells of the developing embryonic gonad, whereas in other tissues, the transcript is rapidly degraded [40]. The D3U-box motif was found to be highly conserved in the dmrt1 3' UTR in the fish lineage (*O. latipes*, *O. curvinotus*, *Takifugu rubripes*, *Tetraodon nigroviridis*, *Epinephelus coioides*, and *Danio rerio*), as well as in other vertebrates, including *Mus musculus*, *Pan troglodytes*, *Macaca mulatta*, and *Homo sapiens*, and even in the ecdysozoan clade (*Anopheles gambiae* and *Bacto-cera oleae*).

Both phylogenetic conservation and presence of the D3U-box in several germ cell tran-
scripts implied the existence of similarly conserved *trans*-acting factor(s) involved in the
synexpression of those transcripts. To identify such factor(s), we undertook an unbiased
approach centred on the D3U-box sequence and based on the evolutionary conservation of
the 'split' motifs of the D3U-box, implying evolutionary conserved *trans*-acting factors. Fur-
ther bioinformatics analyses and literature searches revealed that the D3U-box motif is a puta-
tive target for 2 RNA-binding proteins, namely cug-bp [45,47] and bsf (also known as lrpprc
in mammals [41,42,71]).

EMSAs indicated that Ol-bsf and Ol-cug-bp1, but not Ol-cug-bp2, specifically target and interact with the different parts of the D3U-box, the 3' and the 5' parts, respectively. Additionally, our results suggested that the observed regulation of *dmrt1bY* transcript abundance is likely to be the result of a differential binding of the 2 RNA-binding proteins (Ol-bsf and Ol-cug-bp1) with antagonistic properties, *trans*-regulating RNA stability via the D3U-box.

Being—like its mammalian counterparts—ubiquitously expressed, Ol-cug-bps are highly expressed in the gonads of both sexes in medaka (S4C to S4F Fig) and mice ([72] and S5 Fig). Interestingly, Ol-bsf is specifically expressed in the germ cells during medaka early gonadal primordium formation as well as in the adult ovary and is cytoplasmically localized. This expression pattern resembles the subcellular localization of *Drosophila* bsf, which is present in cytoplasmic particles in oocytes and surrounding nurse cells [41] and in the cytoplasms and nuclei in early embryos [42].

In vivo, using fluorescent reporter medaka fish lines, we could show that, besides obvious
correlated expression levels between Ol-bsf, *dmrt1bY*, and Ol-vas in a subpopulation of germ
cells of the forming gonadal primordium, medaka Ol-bsf is also preferentially expressed in
adult germline stem cells. This restricted and up to now unreported expression pattern might
reveal another so far underappreciated role for bsf/lrpprc in germ cell physiology. Similarly,
immunohistochemistry localization of mouse bsf/lrpprc revealed its presence in a subpopula-
tion of germ cells in mice (S5 Fig), suggesting an evolutionary conserved function that is not
restricted to exclusive regulation of dmrt1.

Reduced expression of bsf/lrpprc in mammalian cell lines [73,74] or flies [75] or because of
a missense mutation in human (French Canadian Leigh Syndrome [76]) resulted in decreased
levels of mitochondrial mRNAs. This led to respiratory chain dysfunction and increased lac-
tate levels in flies and humans [75,76]. Similarly, in medaka, many transcripts with mitochon-
drial function displayed decreased steady-state levels when Ol-bsf expression was reduced
(S12B Fig). Although knockdown of *bsf* in flies affects climbing ability, fecundity, and life
span [75], mutant medaka hatchlings for *Ol-bsf* comparably displayed a significantly reduced
swimming ability (S7 Fig), suggesting—with all the necessary notes of caution, such as in *Dro-
sophila*—possible muscle weakness in relation to mitochondrial dysfunction and energy
metabolism failure [75].
Ol-bsf-heterozygous mutant fish develop normally. Although adult mutant males do not seem to display any observable bias of fertility, adult mutant females have reduced egg production together with significantly lower fertilization rates. In-depth morphological inspection of heterozygote mutant gonads of both sexes revealed that ovaries display an important accumulation of small-sized oocytes compared to wild type. Mutant testes exhibited a discontinuous spermatogenetic flux, likely reflecting uncontrolled spermatogenesis, independent of the seminiferous epithelial cycle. Such a testicular phenotype can be interpreted in the light of dmrt1 loss of expression as a result of low Ol-bsf expression. Similarly, in mice, loss of dmrt1 in germ cells uncouples meiotic initiation from the seminiferous epithelial cycle, resulting in uncontrolled spermatogenesis, too [26]. The mouse celf1 (also known as cugbp1) is predominantly expressed in testis ([72] and S5E and S5F Fig). There, it was demonstrated that celf1 post-transcriptionally represses cyp19a1 (*aromatase*) mRNA translation, by direct binding, to maintain high concentrations of testosterone compatible with spermiogenesis [77]. This situation is reminiscent of the situation we observed in medaka with a strong repression of aromatase expression when both Ol-cug-bps are overexpressed (S12C1 and S12C2 Fig). Additionally, as observed in medaka, mouse celf/cug-bps and bsf/lrpprc are expressed at relevant levels in germ cells, possibly implying a functional conservation across vertebrates.

Whole transcriptome scans using either the medaka or *Drosophila* D3U-box sequences as query resulted in hits enriched for genes specifically expressed in the germ cells, including tra2, sox10, misr2, dead end, and vasa. We provide evidence that Ol- and Ol-cug-bps antagonistically regulate the expression of germ cell transcripts harbouring the D3U-box motif. Certainly most of these regulations occur via an mRNA decay versus stabilization equilibrium after Ol-cug-bps and Ol-bsf targeted the D3U-box motif.

Lastly, to explain the observed differences in the degree of regulation of D3U-box containing transcripts after either Ol-bsf or Ol-cug-bp1/2 modulations, we consider that this follows the degree of conservation of the D3U-box (S1 Fig). It is also intuitive that the location of the D3U-box (5', 3' or coding regions) within the transcripts is of relevance (S1 Fig). For efficient regulation, reasonably high (tra2) or moderate (sox10) conservation and location within the UTRs appears to be more effective than a moderately conserved motif located in the coding sequences (vasa and dead end). However, multiple, highly conserved D3U-boxes nested within the coding region seem to be efficient as well (misr2). Dead end transcripts, for which the D3U-box was identified within the 5' UTR, however, has only strict conservation for the 3' part (CUGCUG) and is only regulated by Ol-cug-bp1 and Ol-cug-bp2 (Fig 5B7 and 5B8) while it expectedly escapes Ol-bsf regulation (Fig 5A7 and 5A8). Altogether, using complementary approaches, our data suggest that the D3U-box motif is—depending on the cellular context—targeted by 2 antagonizing RNA binding proteins, promoting either RNA stabilization in germ cells or degradation in the soma. This new mechanism of dmrt1 RNA stability appears to also regulate the abundance of other transcripts specifically expressed in PGCs, depending of the preservation of the D3U-box motif.

## Materials and methods

### Bioinformatic analyses

Gene, transcript, and UTR annotation, coordinates, and sequence for human, mouse, medaka, and *Tetraodon* were retrieved from EnsEMBL using the EnsEMBL API (version 54). UTR regions spread across several exons were stitched together per transcript, and the resulting sequence was scanned for the presence of the D3U-box consensus matrix and the other motifs.
## Fish maintenance and breeding

All medaka fish used in this study were taken from closed breeding stocks of the wild-type *Carolina Biological Supplies* (Carbio) strain or transgenic lines produced on a wild-type Car-bio background and were kept under standard conditions [5]. Medaka embryos were staged according to Iwamatsu [78]. The fish used in this study from aquaria housed stocks were kept and sampled in accordance with the applicable EU and national German legislation governing animal experimentation. We hold an authorization (568/300-1870/13) of the Veterinary Office of the District Government of Lower Franconia, Germany, in accordance with the German Animal Protection Law (TierSchG).

## Cell culture, transient cell transfection, and immunohistochemistry

Medaka embryonic stem cells (MES-1) were cultured as described [79]. For transfection, cells were grown to 80% confluence in 6-well plates and transfected with 5 µg expression vector using FuGene (Roche, Germany) reagent as described by the manufacturer. After pCS2:: OL-LRPPRC:3XFLAG transfection for 48 hours, cells were fixed in 4% paraformaldehyde (PFA) for 15 minutes, washed with phosphate-buffered saline (PBS) buffer, and then permea-bilized with 0.1% Triton X-100 in PBS for 10 minutes. After blocking in 1% Bovine Serum Albumin (BSA) for 20 minutes, cells were incubated overnight at 4°C in blocking buffer (1% BSA) together with the primary antibody (3-times FLAG, monoclonal anti-FLAG M2, cate-gory number F1804; Sigma-Aldrich). After extensive washes in PBS, cells were then incubated with Alexa 488 conjugated secondary antibody in 1% BSA for 1 hour. Cell nuclei were stained with Hoechst 33343 (Invitrogen) for 5 minutes (1 µg/mL final concentration) and subse-quently mounted using Mowiol 4-88 (Roth). Confocal images were acquired using a Nikon Eclipse C1 laser-scanning microscope (Nikon) and were fitted with a 60× Nikon objective (PL APO, 1.4 NA) and Nikon image software.

## Phylogenetic analysis

The lrpprc phylogenetic tree was built using the online phylogeny.fr automatic pipeline [80]. *lrpprc* sequences were retrieved from public database sequences in the following species and *lrpprc* homologs were retrieved by tblastn searches on the PhyloFish [81] database (http://phylofish.sigenae.org/) using medaka protein (O. latipes, XP_011482612.1) as bait. Sequences were aligned with MUSCLE (version 3.8.31) configured with default settings. After alignment, ambiguous regions were removed with Gblocks (version 0.91b) using the following parame-ters: minimum length of a block after gap cleaning of 10, no gap positions allowed in the final alignment, rejection of all segments with contiguous nonconserved positions bigger than 4, and a minimum number of sequences for a flank position of 85%. The phylogenetic tree was reconstructed using the neighbour joining method implemented in the BioNJ program [82] with *N* = 100 bootstrapping. The resulting phylogenetic tree was exported as a Newick file and edited in Evolview [83]. The public database for *lrpprc* sequences is as follows: *Cyprinodon variegatus* (XP_015233698.1), *Haplochromis burtoni* (XP_005930623.2), *Fundulus heteroclitus* (XP_012720452.1), *Xiphophorus maculatus* (XP_005804231.1), *Poecilia reticulata* (XP_008426945.1), *Lates calcarifer* (XP_018552533.1), *Salmo salar* (XP_014035445.1), *Kryptolebias marmoratus* (XP_017287636.1), *Larimichthys crocea* (KKF31900.1), *Pygocentrus nattereri* (XP_017546072.1), *Ictalurus punctatus* (XP_017319936.1), *Austrofundulus limnaeus* (XP_013886428.1), *Xenopus tropicalis* (NP_001039203.1), *Cynoglossus semilaevis* (XP_008319858.1), *Callorhinchus milii* (XP_007895038.1), *Astyanax mexicanus* (XP_007255829.1), *Chrysemys picta bellii* (XP_005296166.1), *Scleropages formosus* (KPP63655.1), *Bos taurus* (XP_005212770.1), *Cyprinus carpio* (KTG42350.1 and
XP_018937975.1), *Equus caballus* (XP_005600080.1), *H. sapiens* (BAF82705.1), *Rattus norvegicus* (NP_001008519.1), *M. musculus* (AAH59862.1), *Coturnix japonica* (XP_015712995.1), *Serinus canaria* (XP_009094627.1), *Lepisosteus oculatus* (XP_015218688.1), *Amia calva* (AAC LOC100694568.1), *Esox lucius* (XP_012989242.1), *Latimeria chalumnae* (XP_005999623.1), *Sinocyclocheilus anshuiensis* (XP_016313114.1 and XP_016355337.1), *Sinocyclocheilus rhinocerous* (XP_016417572.1 and XP_016400889.1), *Oreochromis niloticus* (XP_003438484.3), *D. rerio* (NP_001136064.1), *Stegastes partitus* (XP_008296301.1), *T. rubripes* (XP_011601434.1), *Nothobranchius furzeri* (XP_015800103.1), and *Gallus* (XP_001234903.3). The PhyloFish species are as follows:

* (DAA_LPPRC.1).2_Anguilla_anguilla;

* Y_Up_12_k31_Locus_13642_Umbra pygmaea;

* W2_Dp_10_k25_Locus_11398_Dallia pectoralis; S2_Pb_10_k31_Locus_2910_Pantodon bulchozi;

* W_Ha_12_k43_Locus_3235_Hiodon alosoides;

* D2_Om_14_k31_Locus_6093_Oncorhynchus mykiss;

* N_St_1_k25_Locus_84_Salmo trutta;

* O_Ot_10_k25_Locus_4385_Oncorhynchus tshawytscha;

* Q_Sf_1_k25_Locus_1875_Salvelinus fontinalis;

* P_Cl_2_k25_Locus_10853_Coregonus lavaretus;

* I_Tt_3_k25_Locus_2938_Thymallus thymallus;

* E_Aa_2_k25_Locus_4938_Alosa alosa;

* Z_Pa_1_k25_Locus_2290_Plecoglossus altivelis;

* U_Cc_1_k25_Locus_265_Coregonus clupeaformis;

* R_Ok_10_k31_Locus_25108_Oncorhynchus keta;

* G_Ph_4_k37_Locus_379_Pangasianodon hypophthalmus;

* K2_Cb_1_k25_Locus_1206_clarias batrachus;

* M_Pf_1_k25_Locus_1736_Perca fluviatilis;

* G2_Ps_1_k65_Locus_195_Polypterus senegalensis;

* K_Gm_1_k43_Locus_6389_Gadus morhua;

* C_Ob_2_k31_Locus_7388_Osteoglossum bicirrhosum;

* B2_Gp_10_k25_Locus_1874_Gnathonemus petersii;

* V_Sc_10_CL11520Contig1_Scyliorhinus caniculata

The celf2 phylogenetic tree was built using the online phylogeny.fr automatic pipeline [80]. Celf2 sequences were retrieved from Ensembl. Sequences were aligned with MUSCLE (version 3.8.31) configured with default settings. After alignment, ambiguous regions were removed with Gblocks (version 0.91b) using the following parameters: minimum length of a block after gap cleaning of 5, gaps were allowed in the final alignment if they were within an appropriate block, all segments with contiguous nonconserved positions bigger than 8 were rejected, and a
minimum number of sequences for a flank position of 55%. The phylogenetic tree was recon-
structed using the neighbour joining method implemented in the BioNJ program [82] with
N = 100 bootstrapping. The resulting phylogenetic tree was exported as a Newick file and
edited in Evolview [83]. The genomic context around celf1 and celf2 was analysed using the
Genomicus website [84]. A few additional genomes were analysed manually by blasting
(tblastn) some fish NCBI genomes with the spotted gar Celf protein and by manually extract-
ing the NCBI gene annotation around these corresponding celf loci.

Locomotor activity

Locomotor activity was determined at hatching stage (stage 39, 9 dpf) either under dark condi-
tions or with 10% light (approximately equal to 100 lux) by measuring the total distance swum
during a period of 5 minutes (described in [85]) for 12 hatchlings in each condition. Briefly,
the larvae were transferred into 12-well plates containing 1 mL Danieau's solution, and the
plate was placed in a Zebrabox equipped with a video camera, infrared light, and filter (View-
Point Life Sciences, Lyon, France) and the ZebraLab Videotrack software (ViewPoint Life Sci-
ences) for tracking. Following a 5-minute-long habituation period with the same light
conditions as for the experimental period, the locomotor activity was recorded. The detection
threshold was 11, the inactive/small threshold was 0.5 cm/s, and the small/large threshold was
1.0 cm/s. The total distance swum for each individual is equal to the sum of distances reached
during inactivity, small and large movements. *p*-Values were calculated using two-tailed
unpaired *t* test with Excel 14.4.8 (Microsoft), and *p* < 0.05 was considered significant. Bars
and error bars indicate mean ± standard deviation.

EMSA

Binding assays were carried out as previously described [38,40]. In detail, 5'-UGGUUCACGU
CUGCUGCAGGUUCUCUGACUCU-3' for the native D3U-box, 5'-UGGUUCACGUCUGCU
GCAGCUCUCUGACUCU-3' for the Olive fruit fly box (Off-box), 5'-UGGUUCACGUUCU
UCAAGACGCUCUGACUCU-3' for the D3U-scrambled box (S-D3U-box), and 5'-UGGUU
CACGUUCUUCACAGGUUCUCUGACUCU-3' for the D3U-minus CUG box (CUGminus-
D3U-box) were synthesized and end-labelled. For radioactive labelling, 50 pmol of the duplex
5' termini were used, together with 50 pmol of gamma-[³²P]-ATP and 20 units of T4 PNK in
1x-adjusted T4 PNK buffer, and were incubated for 20 minutes at 37 °C. Unincorporated
nucleotides were removed through a Sephadex G-50 spin column. For producing Ol-lrpprc,
Ol-cug-bp1, and Ol-cug-bp2 proteins, pCS2::OL-LRPPRC, pCS2::OL-CUG-BP1 or pCS2::
OL-CUG-BP2 plasmids were linearized using KpnI and then in vitro transcribed using mMes-
sage mMachine kit (Ambion). Finally, Ol-bsf, Ol-cug-bp1, and Ol-cug-bp2 proteins were in
vitro translated using Ambion's Retic Lysate Kit from the previously in vitro-transcribed Ol-
lrpprc, Ol-cug-bp1, and Ol-cug-bp2 capped RNAs. DNA binding reaction contained 10
mMTris-HCl (pH 7.9), 100 mM KCl, 10% glycerol, 5 mM MgCl₂, 1 µg torula rRNA, 0.075%
Triton X-100, 1 mM DTT, 1 µg BSA, 0.5 ng radiolabelled probe, and 2 or 4 µL in vitro transla-
tion mix in a total volume of 20 µL. The amount of 1/10 volume heparin (50 mg/mL) was
added just before loading the binding reaction. For control, reticulocyte lysate alone together
with radiolabelled duplex probe was used and did not result in any shift. Binding reactions
were performed on ice for 10 minutes, and complexes were resolved on a 5% native acrylamide
(37.5:1) gel in 0.5 X TBE and then directly subjected to autoradiography. Of note, due to differ-
ent exposure times—which can be appreciated through the different intensities of the free
probes in the respective figures—autoradiographies should not be compared between each
other.
## Establishment of transgenic fluorescent reporter lines and imaging

For a dynamic and in vivo visualization of endogenous Ol-bsf expression, a transgenic fluorescent reporter line was created. The *Ol-bsf* upstream promoter region of the *Ol-bsf* gene (732 bp up to the next upstream gene) was cloned (BamH1 sites) in front of the mCherry ORF of a meganuclease plasmid (BamH1-BSFp-Fw: AAAGGATCCAGTGTGAGTTCTATCAAGCC TGG; BamH1-BSFp-Rv: AAAGGATCCTTCTGTAGCTGCGTAGAGGAAGATC). For the generation of a stable transgenic line, the meganuclease protocol was used [86]. Briefly, approximately 10 to 15 pg of total vector DNA in a volume of 500 pL injection solution containing *I-SceI* meganuclease was injected into the cytoplasm of 1-cell-staged medaka embryos (Carbio strain). Adult F0 fish were mated to each other, and the offspring were tested for the presence of the transgene by fluorescence check. Siblings from positive F1 fish were raised to adulthood and tested again for fluorescence. Tg[vasa:GFP] and Tg[dmrt1bY:GFP] transgenic lines were described earlier [11,29,48,59].

## Visualization of PGCs

For PGC visualization and counting, the GFP-nos1 3' UTR construct that includes the mmGFP5 ORF cloned upstream of the 3' UTR of the zebrafish *nanos1* gene [5,87] was injected at 1 cell stage. N = 34 and 28 for wild-type and MO-bsf-injected embryos, respectively. For imaging, embryos, hatchlings or tissues were mounted with 1.2% low melting temperature agarose.

## In vivo expression regulation analyses and real-time PCR

For producing Ol-bsf, Ol-cug-bp1, and Ol-cug-bp2 capped mRNAs, pCS2::OL-BSF, pCS2:: OL-CUG-BP1, or pCS2::OL-CUG-BP2 plasmids were linearized using *KpnI* and then in vitro transcribed using mMessage mMachine kit (Ambion). For overexpression, 1 nanolitre was injected into the cytoplasm of 1-cell-stage medaka embryos. For mRNA stability assays, equimolar amounts of control and D3U-box-containing RNAs were injected.

Total RNAs were extracted from fish tissues or embryos using the TRIZOL reagent (Invitrogen) according to the supplier's recommendation. After DNase treatment, reverse transcription was performed with 2 µg total RNA using RevertAid First Strand Synthesis kit (Fermentas) and random primers. Real-time quantitative PCR was carried out with SYBR Green reagents, and amplifications were detected with an i-Cycler (Biorad). All results are averages of at least 2 independent reverse transcription reactions. Error bars represent the standard deviation of the mean. Relative expression levels (according to the equation 2-DeltaCT) were calculated after correction of expression of elongation factor 1 alpha (ef1alpha).

## Morpholino injections

For Ol-bsf knockdown experiments, embryos were injected with a splice morpholino: 5'-TTGATGACTGGCCTGCCAACCTGTC-3' targeting the 3' splice junction of *Ol-bsf* exon 2 (see S3 Fig). The most efficient dose (4 mg/mL) was experimentally determined and the specificity of the targeting confirmed in control experiments (see S3 Fig).

## Microarray

Total RNAs were extracted from adult medaka gonads using the Tri-reagent (Sigma-Aldrich) according to the supplier's recommendation. The total RNA yield was estimated using a NanoDrop ND-1000 spectrophotometer (Labtech, Palaiseau, France), and RNA integrity was checked by means of an Agilent Bioanalyzer (Agilent Technologies, Massy, France). Medaka
gene expression profiling was conducted using an Agilent 8x60K high-density oligonucleotide microarray (GEO platform GPL24100). Labelling and hybridization steps were performed following the 'One-Color Microarray-Based Gene Expression Analysis (Low Input Quick Amp labelling)' Agilent protocol. Briefly, for each sample, 150 ng of total RNA was amplified and labelled using Cy3-CTP. Yield (>825 ng cRNA) and specific activity (>6 pmol of Cy3 per µg of cRNA) of Cy3-cRNA produced were checked with the Nanodrop. The amount of 600 ng of Cy3-cRNA was fragmented and hybridized on a sub-array. Hybridization was carried out for 17 hours at 65 °C in a rotating hybridization oven prior to washing and scanning with an Agilent Scanner (Agilent DNA Microarray Scanner, Agilent Technologies, Massy, France) using the standard parameters for a gene expression 8x60K oligo-array (3 µm and 20 bits). Data were then obtained with the Agilent Feature Extraction software (10.7.3.1) according to the appropriate GE protocol (GE1_107_Sep09). The arrays were normalized (scale normalization) and log-transformed using Genespring Software (version 14.5). A *t* test analysis (*p* < 0.05) was employed to determine the genes that were the most differentially expressed between the 2 conditions. Microarray data sets have been deposited to the GEO-NCBI with the accession number GSE 104726. GO was performed using the panther program (http://geneontology.org/).

## CRISPR/Cas9 genome editing

Identification of CRISPR/Cas9 target sites and design of oligonucleotides of were performed by the use of the ZiFiT software (http://zifit.partners.org/ZiFiT/Disclaimer.aspx). For preparation of sgRNAs, the DR274 plasmid (Addgene number 42250) was first linearized with *Bsa*1, electrophoresed in a 2% agarose gel, and purified. Pairs of complementary oligonucleotides were annealed (40 mM Tris-HCl [pH 8.0], 20 mM MgCl₂, and 50 mM NaCl buffer) by heating at 95 °C for 2 minutes and then cooled down slowly to 25 °C within 1 hour. The double-stranded oligonucleotides were then ligated into the linearized pDR274 vector. Different sgRNAs were designed to target several sites within the Ol-BSF gene (see S6 Fig) in order to create deletions. After linearization with *Dra*1 and *Not*1, respectively, pDR274 and pCS2-*nCas*9n plasmids were used for generating either sgRNAs or *Cas*9 RNAs. The synthesized RNAs were then injected into 1-cell-staged embryos at the following concentrations: 25 ng/µL for each sgRNAs and 100 ng/µL for the *Cas*9 mRNA. CRISPR-positive fish were then screened for mutations using PCR primers flanking the site of deletion (see S6 Fig). The inferred mutant protein is presented in S6C Fig. Mutant fish used in this study have been outcrossed for at least 5 generations.

## Histology

Gonads were fixed for 48 hours in Bouin-Holland fluid and then dehydrated serially in aqueous 70% and 95% ethanol, ethanol/butanol (5:95), and butanol. Tissues were embedded in paraffin, and 5-µm mid-sagittal gonad sections were stained with Regaud's haematoxylin and haematoxylin-eosin-safran.

## Supporting information

S1 Fig. Presence of the D3U-box in different UTRs and coding sequence of gonadal genes. Presence of the D3U-box (matrix-scan) was evaluated within the transcriptome of the medaka fish. D3U-box, *dmrt1* 3' UTR box. (DOCX)
**S2 Fig. Phylogeny and synteny analysis of Ol-bsf and Ol-cug-bp ohnologs.** (A) Circular cladogram representation of the phylogenetic tree of Lrpprc proteins in jawed vertebrates (gnatosthomes). This phylogeny shows that *lrpprc* genes were retained as single copies in most jawed vertebrates even following whole genome duplications (red stars), e.g., the teleost-specific duplication or the salmonid-specific duplication. *Lrpprc* is, however, present in duplicated copies in the Cyprininae (tree branches in red). Bootstraps (*N* = 100) values are indicated in each tree node when judged significant (>0.75). Tree branches are depicted in blue for lobefinal vertebrates and cartilaginous fish and in black for teleosts with the exception of Cyprininae in red. (B) Gene evolution of *celf2* genes in some teleosts. The phylogeny on the left is a dendogram representation of *celf2* gene phylogeny in teleosts given as an indication as only a few nodes are supported by good bootstraps' values (*N* = 100, mentioned in each tree nodes when judged significant, i.e., >0.7). The teleost fish whole genomic duplication (3R) is indicated by a red star. The left part of the figure is a representation of the evolution of the genomic context around the *celf2* gene. After the 3R whole genome duplication, *celf2*—which is a single copy gene on the Chr 8 of the spotted gar genome—was duplicated in two 3R ohnologs, *celf2a* and *celf2b*, that were not retained as 2 copies in all teleost fish. The genomic context of the *celf2a* and *celf2b* paralogous regions clearly indicates a partition of the ancestral region found in spotted gar. The *celf2a* gene was retained in all species investigated, but the *celf2b* gene seems to have been lost in Otophysi or at least in *D. rerio* (Cypriniformes), *Astyanax mexicanus* (Characiformes), and *Ictalurus punctatus* (Siluriformes). *celf2*, CUGBP Elav-like family member 2; Chr 8, Chromosome 8; Lrpprc, leucine rich pentatricopeptide repeat containing; Ol-BSF, *Oryzias latipes* Bicoid Stability Factor; Ol-CUG-BP, *Oryzias latipes* CUG-binding protein. (TIF)

**S3 Fig. Analysis of morpholino efficiency and level of Ol-bsf down-regulation.** For in vivo transient down-regulation of Ol-bsf, a splice morpholino was designed to encompass the splice junction between exon 2 and intron 2 of the *Ol-bsf* gene in order to induce aberrant splicing and frame shit of the ORF. To show to what extend the splicing/activity of *Ol-bsf* was impacted, RT-PCR using exons 1, 2, and 3 spanning primers together with cDNAs from different stages of morpholino-injected embryos was achieved. E2, exon 2; i2, intron 2; Ol-BSF, *Oryzias latipes* Bicoid Stability Factor; RT-PCR, Reverse Transcription-Polymerase Chain Reaction.

(TIF)

**S4 Fig. Real-time PCR quantification of Ol-cug-bp1, Ol-cug-bp2, and Ol-bsf during embryogenesis and in adult tissues.** (A and C) During embryonic development, both Ol-cug-bp ohnologs are expressed in a complementary manner. Being likely maternally deposited the expression of Ol-cug-bp1 rapidly decreases after mid-blastula transition (stage 10) to remain virtually off up to hatching stage. On the other hand, low expression of Ol-cug-bp2 is detected until stage 25 and rapidly increases by stage 33. (B and D) In adult tissues, both Ol-cug-bp ohnologs are expressed in brain, muscles, and gonads; ol-cug-bp2 is additionally expressed in eyes and skin. Both ohnologs are higher expressed in male gonads than in female gonads. (E and F) In adult tissues, Ol-bsf is ubiquitously present although higher expression is observed in gonads of both sexes. Underlying data for (A to F) can be found in S1 Data.

(TIF)

**S5 Fig. Lrrprc and celf1, but not celf2, are expressed in mouse embryonic gonads.** (A to H) ISHs on sagittal sections of 14.5 dpc mouse embryos showed expression of Lrrprc (A to D) and celf1 (E to H) most likely in germ cells within testis cords (B and F) and germ cells within the ovary (D and H). In contrast, no celf2 expression was detected in developing gonads (I–L).
However, celf2 expression was detected in other tissues, such as part of the brain and dorsal root ganglia. Scale bars: 1 mm for A, C, E, G, I, and K; 10 mm for B, D, F, H, J, and L. (M–R) Immunofluorescent detection of LRPPRC (M, N, P, Q) and DDX4/VASA (O, R) in adult mouse testes (M–O) and ovaries (P–R). In adult testes, lrpprc is expressed in one subpopulation of germ cells; compared lrpprc staining on (M) and (N) with vasa staining on (O) where most of the germ cells (except some spermatogonia) remain stained by vasa. According to the position of lrpprc-positive cells (arrowheads in M and N) in the seminiferous tubule (not basal and below round spermatids) and to the fact that lrpprc-positive germ cells are those with the largest nucleus, lrpprc-positive cells seem to be spermatocytes at the pachytene stage. In adult ovaries (P–R), lrpprc is mainly expressed into the oocytes of primordial, primary and young secondary follicles (see arrows on [P] and [Q]). LRPPRC staining disappears from the oocyte of secondary follicles that are clearly stained for vasa in (R) (compared stars in [Q] and [R]). Scale bars: 200 µm for M to R. dph, days post hatching; ISH, in situ hybridisation; LRPPRC, leucine rich pentatricopeptide repeat containing.

(TIF)

**S6 Fig. Generation *Ol-bsf* knockout medaka fish after genome editing by CRISPR/Cas9 method.** (A) Several guide RNA were designed in order to target different locations on the *Ol bsf* gene (targets 1, 2, and 6). (B) After injection of different combinations of guide RNAs together with the Cas9 mRNA, putative edited fish were subjected to RT-PCR using primer sets flanking the cutting sites. Lines displaying deletions within the *Ol bsf* gene (red stars) were kept for further investigations. (C) Deletions obtained within the *Ol bsf* gene (left panel) result in a truncated translated Ol-bsf protein (right panel). CRISPR-Cas9, clustered regularly inter-spaced short palindromic repeats/CRISPR-associated protein 9; Ol-BSF, *Oryzias latipes* Bicoid Stability Factor; RT-PCR, Reverse Transcription- Polymerase Chain Reaction.

(TIF)

**S7 Fig. Locomotor activity of *Ol-bsf* mutant fish versus wild type.** Locomotor activity (*Ol bsf* mutants versus wild type) was determined at hatching stage (stage 39, 9 dpf) either under dark conditions (A) or with 10% light (approximately equal to 100 lux [panel B]) by measuring the total distance swum during a period of 5 minutes. (C) The total distance swum for each individual is equal to the sum of distances reached during inactivity, small and large movements. Bars and error bars indicate mean ± standard deviation. N = 12 for each condition. Underlying data for (C) can be found in S1 Data. dpf, days post fertilization; Ol-BSF, *Oryzias latipes* Bicoid Stability Factor.

(TIF)

**S8 Fig. Ovarian phenotypes of the *Ol-bsf* mutant fish.** Morphological inspection of heterozygote mutant ovaries discloses a significant accumulation of small sized-oocytes compared to wild type (A) and (A1 to A5 for details and statistical analyses). (A) Overall size distribution of the oocytes in 9 wild-type and 9 Ol-BSF(+/-) adult ovaries. Each gonad (testes or ovaries) was sectioned through the mid-sagittal plan (see also Materials and methods). Underlying data for (A) can be found in S1 Data. Ol-BSF, *Oryzias latipes* Bicoid Stability Factor.

(TIF)

**S9 Fig. Ovarian phenotypes of the *Ol-BSF* mutant fish.** (A to R) Mid-sagittal sections of the ovaries utilized for counting the oocytes in S8 Fig. Each gonad (testes or ovaries) was sectioned through the mid-sagittal plan (see also Materials and methods). Underlying data for (A) can be found in S1 Data. Ol-BSF, *Oryzias latipes* Bicoid Stability Factor.

(TIF)
**S10 Fig. Testicular phenotypes of the Ol-bsf mutant fish.** (A–J compared to K–T) Heterozygote mutant testes (A–J) exhibit a decreased number of spermatogonia with accumulation of type 2 spermatocytes, spermatids, and sperm within the most external layers of the seminiferous epithelium (arrowheads) compared to wild-type testes (K–T). Either 10 different wild-type (A–J) or *Ol-bsf*-deficient (K–T) testes were analysed. Mid-sagittal gonad sections were stained with haematoxylin–eosin–safran. Each gonad (testes or ovaries) was sectioned through the mid-sagittal plan (see also *Material and methods*). Ol-BSF, *Oryzias latipes* Bicoid Stability Factor.
(TIF)

**S11 Fig. Fertility test.** Egg numbers and fertilization rates were recorded over a period of 9 days for the following crosses: (A) male *Ol-bsf*(-/+) × female *Ol-bsf*(-/+); (B) male *Ol-bsf*(-/+) × female wild type; (C) male wild type × female *Ol-bsf*(-/+). Underlying data for (A to C) can be found in S1 Data. BSF, bicoid stability factor.
(TIF)

**S12 Fig. Microarray data and mitochondrial gene quantification.** (A) Microarray. Adult testes of either *bsf*⁺/⁻ or wild-type animals were subjected to microarray (see *Materials and methods*). GO term analysis reveals that in mutant testes partially depleted for the *bsf* gene, rRNA processing is particularly affected. Ol-bsf and Ol-cug-bp2 are down- and up-regulated, respectively, in mutant animals compared to wild type. Of note, and in accordance with the literature, a significant proportion (10.1%) of the down-regulated genes code for proteins localized in the mitochondria. Finally, supporting our observations that lowering *ol-bsf* transcription (morpholino injection in Fig 7E) resulted in up-modulation of germ cell number and that mutant gonads presented an increase of germ cells committing to gametogenesis (Fig 6), our microarray analysis reveals a general up-regulation of genes involved in germ cell proliferation or differentiation. (B) RNA levels of different mitochondrial genes (Cox1, Cox2, ND1, ND5, and CytB) were quantified by real-time PCR after BSF-morpholino injections and compared to wild type. Most of the mitochondrial genes are down-regulated when the level of Ol-bsf decreases. (C) Modulation of RNA levels of the cyp19a1 (aromatase) gene after overexpression of Ol-cugbp1 or Ol-cug-bp2. Underlying data for (B and C) can be found in S1 Data. GO, gene ontology.
(TIFF)

**S1 Table. Gene evolution of *cugbp elav-like family member 2* genes in some teleosts. *celf2, cugbp elav-like family member 2*.
(DOCX)**

**S1 Data. Underlying data.**
(XLSX)

**Acknowledgments**

Acknowledgements

The authors thank the fish facilities at the University of Würzburg and at the INRA Rennes (C. Duret, F. Borel, and A. Patinote).

**Author Contributions**

**Conceptualization:** Amaury Herpin, Manfred Schartl.

**Formal analysis:** Amaury Herpin, Clara Gobé, Aurélie Le Cam, Yann Guiguen.

**Funding acquisition:** Amaury Herpin.
**Investigation:** Amaury Herpin, Cornelia Schmidt, Clara Gobé, Martina Regensburger, Aurélie Le Cam, Jérome Montfort, Mateus C. Adolfi, Christina Lillesaar, Dagmar Wilhelm, Michael Kraeusling, Brigitte Mourot, Béatrice Porcon, Laurence Ettwiller, Dirk Dolle, Yann Guiguen, Manfred Schartl.

**Methodology:** Amaury Herpin, Christina Lillesaar, Maëlle Pannetier, Eric Pailhoux, Laurence Ettwiller, Dirk Dolle.

**Project administration:** Amaury Herpin.

**Software:** Susanne Kneitz, Laurence Ettwiller, Dirk Dolle.

**Supervision:** Amaury Herpin, Manfred Schartl.

**Validation:** Amaury Herpin.

**Writing – original draft:** Amaury Herpin, Manfred Schartl.

**Writing – review & editing:** Amaury Herpin, Manfred Schartl.

## References

1. Brennan J, Capel B. One tissue, two fates: molecular genetic events that underlie testis versus ovary development. Nat Rev Genet. 2004; 5: 509–521. https://doi.org/10.1038/nrg1381 PMID: 15211353

2. Devlin RH, Nagahama Y. Sex determination and sex differentiation in fish: an overview of genetic, physiological, and environmental influences. Aquaculture. 2002; 208: 191–364. https://doi.org/10.1016/S0044-8486(02)00057-1

3. Lin Y-T, Capel B. Cell fate commitment during mammalian sex determination. Curr Opin Genet Dev. 2015; 32: 144–152. https://doi.org/10.1016/j.gde.2015.03.003 PMID: 25841206

4. Rios-Rojas C, Bowles J, Koopman P. On the role of germ cells in mammalian gonad development: quiet passengers or back-seat drivers? Reprod Camb Engl. 2015; 149: R181–191. https://doi.org/10.1530/REP-14-0663 PMID: 25628441

5. Herpin A, Fischer P, Liedtke D, Kluever N, Neuner C, Raz E, et al. Sequential SDF1α and b-induced mobility guides Medaka PGC migration. Dev Biol. 2008; 320: 319–327. https://doi.org/10.1016/j.ydbio.2008.03.030 PMID: 18440502

6. Kurokawa H, Saito D, Nakamura S, Katoh-Fukui Y, Ohta K, Baba T, et al. Germ cells are essential for sexual dimorphism in the medaka gonad. Proc Natl Acad Sci U S A. 2007; 104: 16958–16963. https://doi.org/10.1073/pnas.0609932104 PMID: 17940041

7. Marlow F. Primordial Germ Cell Specification and Migration. F1000Research. 2015; 4. https://doi.org/10.12688/f1000research.6995.1 PMID: 26918157

8. Nakamura S, Kobayashi D, Aoki Y, Yokoi H, Ebe Y, Wittbrodt J, et al. Identification and lineage tracing of two populations of somatic gonadal precursors in medaka embryos. Dev Biol. 2006; 295: 678–688. https://doi.org/10.1016/j.ydbio.2006.03.052 PMID: 16682019

9. Paksa A, Raz E. Zebrafish germ cells: motility and guided migration. Curr Opin Cell Biol. 2015; 36: 80–85. https://doi.org/10.1016/j.ceb.2015.07.007 PMID: 26232877

10. Johnson AD, Richardson E, Bachvarova RF, Crother BI. Evolution of the germ line-soma relationship in vertebrate embryos. Reprod Camb Engl. 2011; 141: 291–300. https://doi.org/10.1530/REP-10-0474 PMID: 21228047

11. Tanaka M, Kinoshita M, Kobayashi D, Nagahama Y. Establishment of medaka (Oryzias latipes) transgenic lines with the expression of green fluorescent protein fluorescence exclusively in germ cells: a useful model to monitor germ cells in a live vertebrate. Proc Natl Acad Sci U S A. 2001; 98: 2544–2549. https://doi.org/10.1073/pnas.041315498 PMID: 11226275

12. Herpin A, Schartl M. Dtmt1 genes at the crossroads: a widespread and central class of sexual development factors in fish. FEBS J. 2011; 278: 1010–1019. https://doi.org/10.1111/j.1742-4658.2011.08030.x PMID: 21281449

13. Herpin A, Schartl M. Dtmt1 genes at the crossroads: a widespread and central class of sexual development factors in fish. FEBS J. 2011; 278: 1010–1019. https://doi.org/10.1111/j.1742-4658.2011.08030.x PMID: 21281449

14. Matson CK, Zarkower D. Sex and the singular DM domain: insights into sexual regulation, evolution and plasticity. Nat Rev Genet. 2012; 13: 163–174. https://doi.org/10.1038/nrg3161 PMID: 22310892
15. Raymond CS, Murphy MW, O'Sullivan MG, Bardwell VJ, Zarkower D. Dmrt1, a gene related to worm and fly sexual regulators, is required for mammalian testis differentiation. *Genes Dev.* 2000; 14: 2587–2595. PMID: 11040213

16. Chong T, Collins JJ, Brubacher JL, Zarkower D, Newmark PA. A sex-specific transcription factor controls male identity in a simultaneous hermaphrodite. *Nat Commun.* 2013; 4: 1814. https://doi.org/10.1038/ncomms2811 PMID: 23652002

17. Kato Y, Kobayashi K, Watanabe H, Iguchi T. Environmental sex determination in the branchiopod crustacean Daphnia magna: deep conservation of a Doublesex gene in the sex-determining pathway. *PLoS Genet.* 2011; 7: e1001345. https://doi.org/10.1371/journal.pgen.1001345 PMID: 21455482

18. Yoshimoto S, Okada E, Umemoto H, Tamura K, Uno Y, Nishida-Umehara C, et al. A W-linked DM-domain gene, DM-W, participates in primary ovary development in Xenopus laevis. *Proc Natl Acad Sci U S A.* 2008; 105: 2469–2474. https://doi.org/10.1073/pnas.0712244105 PMID: 18268317

19. Chen S, Zhang G, Shao C, Huang Q, Liu G, Zhang P, et al. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. *Nat Genet.* 2014; 46: 253–260. https://doi.org/10.1038/ng.2890 PMID: 24487278

20. Smith CA, Roeszler KN, Ohnesorg T, Cummins DM, Farlie PG, Doran TJ, et al. The avian Z-linked gene DMRT1 is required for male sex determination in the chicken. *Nature*. 2009; 461: 267–271. https://doi.org/10.1038/nature08298 PMID: 19710650

21. Matsuda M, Nagahama Y, Shinomiya A, Sato T, Matsuda C, Kobayashi T, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish. *Nature*. 2002; 417: 559–563. https://doi.org/10.1038/nature751 PMID: 12037570

22. Nanda I, Kondo M, Hornung U, Asakawa S, Winkler C, Shimizu A, et al. A duplicated copy of DMRT1 in the sex-determining region of the Y chromosome of the medaka, Oryzias latipes. *Proc Natl Acad Sci U S A.* 2002; 99: 11778–11783. https://doi.org/10.1073/pnas.182314699 PMID: 12193652

23. Ledig S, Hiort O, Wünsch L, Wieacker P. Partial deletion of DMRT1 causes 46,XY ovotesticular disorder of sexual development. *Eur J Endocrinol.* 2012; 167: 119–124. https://doi.org/10.1530/EJE-12-0136 PMID: 22573722

24. Raymond CS, Parker ED, Kettlewell JR, Brown LG, Page DC, Kusz K, et al. A region of human chromosome 9p required for testis development contains two genes related to known sexual regulators. *Hum Mol Genet.* 1999; 8: 989–996. PMID: 10332030

25. Kim S, Bardwell VJ, Zarkower D. Cell type-autonomous and non-autonomous requirements for Dmrt1 in postnatal testis differentiation. *Dev Biol.* 2007; 307: 314–327. https://doi.org/10.1016/j.ydbio.2007.04.046 PMID: 17540358

26. Matson CK, Murphy MW, Griswold MD, Yoshida S, Bardwell VJ, Zarkower D. The mammalian double-sex homolog DMRT1 is a transcriptional gatekeeper that controls the mitosis versus meiosis decision in male germ cells. *Dev Cell.* 2010; 19: 612–624. https://doi.org/10.1016/j.devcel.2010.09.010 PMID: 20951351

27. Matson CK, Murphy MW, Sarver AL, Griswold MD, Bardwell VJ, Zarkower D. DMRT1 prevents female reprogramming in the postnatal mammalian testis. *Nature*. 2011; 476: 101–104. https://doi.org/10.1038/nature10239 PMID: 21775990

28. Herpin A, Schartl M. Plasticity of gene-regulatory networks controlling sex determination: of masters, slaves, usual suspects, newcomers, and usurpators. *EMBO Rep.* 2015; 16: 1260–1274. https://doi.org/10.15252/embr.201540667 PMID: 26358957

29. Herpin A, Schindler D, Kraiss A, Hornung U, Winkler C, Schartl M. Inhibition of primordial germ cell proliferation by the medaka male determining gene Dmrt I bY. *BMC Dev Biol.* 2007; 7: 99. https://doi.org/10.1186/1471-213X-7-99 PMID: 17760954

30. Krentz AD, Murphy MW, Kim S, Cook MS, Capel B, Zhu R, et al. The DM domain protein DMRT1 is a dose-sensitive regulator of fetal germ cell proliferation and pluripotency. *Proc Natl Acad Sci U S A.* 2009; 106: 22323–22328. https://doi.org/10.1073/pnas.0905431106 PMID: 20007774

31. Shao C, Li Q, Chen S, Zhang P, Lian J, Hu Q, et al. Epigenetic modification and inheritance in sexual reversal of fish. *Genome Res.* 2014; 24: 604–615. https://doi.org/10.1101/gr.162172.113 PMID: 24487721

32. Wen AY, You F, Sun P, Li J, Xu DD, Wu ZH, et al. CpG methylation of dmrt1 and cyp19a promoters in relation to their sexual dimorphic expression in the Japanese flounder Paralichthys olivaceus. *J Fish Biol.* 2014; 84: 193–205. https://doi.org/10.1111/jfb.12277 PMID: 24372528

33. Ge C, Ye J, Zhang H, Zhang Y, Sun W, Sang Y, et al. Dmrt1 induces the male pathway in a turtle species with temperature-dependent sex determination. *Dev Camb Engl.* 2017; 144: 2222–2233. https://doi.org/10.1242/dev.152033 PMID: 28506988
34. Lei N, Heckert LL. Sp1 and Egr1 regulate transcription of the Dmrt1 gene in Sertoli cells. Biol Reprod. 2002; 66: 675-684. PMID: 11870074

35. Lei N, Heckert LL. Gata4 regulates testis expression of Dmrt1. Mol Cell Biol. 2004; 24: 377-388. https://doi.org/10.1128/MCB.24.1.377-388.2004 PMID: 14673170

36. Cui J, Shen X, Zhao H, Nagahama Y. Genome-wide analysis of Sox genes in Medaka (Oryzias latipes) and their expression pattern in embryonic development. Cytogenet Genome Res. 2011; 134: 283-294. https://doi.org/10.1159/000329480 PMID: 21720154

37. Matsuda M, Shinomiya A, Kinoshita M, Suzuki A, Kobayashi T, Paul-Prasanth B, et al. DMY gene induces male development in genetically female (XX) medaka fish. Proc Natl Acad Sci U S A. 2007; 104: 3865-3870. https://doi.org/10.1073/pnas.0611707104 PMID: 17360444

38. Herpin A, Braasch I, Kraeussling M, Schmidt C, Thoma EC, Nakamura S, et al. Transcriptional rewiring of the sex determining dmrt1 gene duplicate by transposable elements. PLoS Genet. 2010; 6: e1000844. https://doi.org/10.1371/journal.pgen.1000844 PMID: 20169179

39. Schartl M, Schories S, Wakamatsu Y, Nagao Y, Hashimoto H, Bertin C, et al. Sox5 is involved in germ-cell regulation and sex determination in medaka following co-option of nested transposable elements. BMC Biol. 2018; 16: 16. https://doi.org/10.1186/s12915-018-0485-8 PMID: 29378592

40. Herpin A, Nakamura S, Wagner TU, Tanaka M, Schartl M. A highly conserved cis-regulatory motif directs differential gonadal synexpression of Dmrt1 transcripts during gonad development. Nucleic Acids Res. 2009; 37: 1510-1520. https://doi.org/10.1093/nar/gkn1065 PMID: 19139075

41. Mancebo R, Zhou X, Shillinglaw W, Henzel W, Macdonald PM. BSF binds specifically to the bicoid mRNA 3' untranslated region and contributes to stabilization of bicoid mRNA. Mol Cell Biol. 2001; 21: 3462-3471. https://doi.org/10.1128/MCB.21.10.3462-3471.2001 PMID: 11313472

42. De Renzis S, Elemento O, Tavazoie S, Wieschaus EF. Unmasking activation of the zygotic genome using chromosomal deletions in the Drosophila embryo. PLoS Biol. 2007; 5: e117. https://doi.org/10.1371/journal.pbio.0050117 PMID: 17456005

43. Paillard L, Omilli F, Legagneux V, Bassez T, Maniey D, Osborne HB. EDEN and EDEN-BP, a cis element and an associated factor that mediate sequence-specific mRNA deadenylation in Xenopus embryos. EMBO J. 1998; 17: 278-287. https://doi.org/10.1093/emboj/17.1.278 PMID: 9427761

44. Philips AV, Timchenko LT, Cooper TA. Disruption of splicing regulated by a CUG-binding protein in myotonic dystrophy. Science. 1998; 280: 737-741. PMID: 9563950

45. Timchenko LT, Miller JW, Timchenko NA, DeVore DR, Datar KV, Lin L, et al. Identification of a (CUG)n triplet repeat RNA-binding protein and its expression in myotonic dystrophy. Nucleic Acids Res. 1996; 24: 4407-4414. PMID: 8948631

46. Paillard L, Legagneux V, Maniey D, Osborne HB. c-Jun ARE targets mRNA deadenylation by an EDEN-BP (embryo deadenylation element-binding protein)-dependent pathway. J Biol Chem. 2002; 277: 3232-3235. https://doi.org/10.1074/jbc.M109362200 PMID: 11707455

47. Moraes KCM, Wilusz CJ, Wilusz J. CUG-BP binds to RNA substrates and recruits PARN deadenylase. RNA N Y N. 2006; 12: 1084-1091. https://doi.org/10.1261/rna.59606 PMID: 16601207

48. Hornung U, Herpin A, Schartl M. Expression of the male determining gene dmrt1bY and its autosomal coorthologue dmrt1a in medaka. Sex Dev Genet Mol Biol Evol Endocrinol Embryol Pathol Sex Determ Differ. 2007; 1: 197-206. https://doi.org/10.1159/000102108 PMID: 18391530

49. Webster KA, Schach U, Ordaz A, Steinfeld JS, Draper BW, Siegfried KR. Dmrt1 is necessary for male sexual development in zebrafish. Dev Biol. 2017; 422: 33-46. https://doi.org/10.1016/j.ydbio.2016.12.008 PMID: 27940159

50. Shiraishi E, Imazato H, Yamamoto T, Yokoi H, Abe S, Kitano T. Identification of two teleost homologs of the Drosophila sex determination factor, transformer-2 in medaka (Oryzias latipes). Mech Dev. 2004; 121: 991-996. https://doi.org/10.1016/j.mod.2004.04.013 PMID: 15210204

51. Amrein H, Gorman M, Nöthiger R. The sex-determining gene tra-2 of Drosophila encodes a putative RNA binding protein. Cell. 1988; 55: 1025-1035. PMID: 3144434

52. Shiraishi E, Yoshinaga N, Miura T, Yokoi H, Wakamatsu Y, Abe S-I, et al. Mullerian inhibiting substance is required for germ cell proliferation during early gonadal differentiation in medaka (Oryzias latipes). Endocrinology. 2008; 149: 1813-1819. https://doi.org/10.1210/en.2007-1535 PMID: 18162517

53. Weidinger G, Stebler J, Slanchev K, Dumstrei K, Wise C, Lovell-Badge R, et al. dead end, a novel vertebrate germ plasm component, is required for zebrafish primordial germ cell migration and survival. Curr Biol CB. 2003; 13: 1429-1434. PMID: 12932328

54. Gross-Thebing T, Yigit S, Pfeiffer J, Reichman-Fried M, Bandemer J, Ruckert C, et al. The Vertebrate Protein Dead End Maintains Primordial Germ Cell Fate by Inhibiting Somatic Differentiation. Dev Cell. 2017; 43: 704-715.e5. https://doi.org/10.1016/j.devcel.2017.11.019 PMID: 29257950
55. Hartung O, Forbes MM, Marlow FL. Zebrafish vasa is required for germ-cell differentiation and maintenance. Mol Reprod Dev. 2014; 81: 946-961. https://doi.org/10.1002/mrd.22414 PMID: 25257909

56. Siegfried KR, Nüsslein-Volhard C. Germ line control of female sex determination in zebrafish. Dev Biol. 2008; 324: 277-287. https://doi.org/10.1016/j.ydbio.2008.09.025 PMID: 18930041

57. Siegfried KR. In search of determinants: gene expression during gonadal sex differentiation. J Fish Biol. 2010; 76: 1879-1902. https://doi.org/10.1111/j.1095-8649.2010.02594.x PMID: 20557645

58. Nakamura S, Watakabe I, Nishimura T, Toyoda A, Taniguchi Y, Tanaka M. Analysis of medaka sox9 orthologue reveals a conserved role in germ cell maintenance. PLoS ONE. 2012; 7: e29982. https://doi.org/10.1371/journal.pone.0029982 PMID: 22253846

59. Nishimura T, Herpin A, Kimura T, Hara I, Kawasaki T, Nakamura S, et al. Analysis of a novel gene, Sdgc, reveals sex chromosome-dependent differences of medaka germ cells prior to gonad formation. Dev Camb Engl. 2014; 141: 3363-3369. https://doi.org/10.1242/dev.106864 PMID: 25078651

60. Herpin A, Rohr S, Riedel D, Kluever N, Raz E, Schartl M. Specification of primordial germ cells in medaka (Oryzias latipes). BMC Dev Biol. 2007; 7: 3. https://doi.org/10.1186/1471-213X-7-3 PMID: 17217535

61. Herpin A, Rohr S, Riedel D, Kluever N, Raz E, Schartl M. Specification of primordial germ cells in medaka (Oryzias latipes). BMC Dev Biol. 2007; 7: 3. https://doi.org/10.1186/1471-213X-7-3 PMID: 17217535

62. Herpin A, Adolfi MC, Nicol B, Hinzmann M, Schmidt C, Klughammer J, et al. Divergent expression regulation of gonad development genes in medaka shows incomplete conservation of the downstream regulatory network of vertebrate sex determination. Mol Biol Evol. 2013; 30: 2328-2346. https://doi.org/10.1093/molbev/mst130 PMID: 23883523

63. Nakamura S, Kobayashi K, Nishimura T, Higashijima S, Tanaka M. Identification of germline stem cells in the ovary of the teleost medaka. Science. 2010; 328: 1561-1563. https://doi.org/10.1126/science.1185473 PMID: 20488987

64. Uribe MC, Grier HJ, Mejía-Roa V. Comparative testicular structure and spermatogenesis in bony fishes. Spermatogenesis. 2015; 4. https://doi.org/10.4161/21565562.2014.983400 PMID: 26413405

65. Keene JD. RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007; 8: 533-543. https://doi.org/10.1038/nrg2111 PMID: 17572691

66. Moore MJ. From birth to death: the complex lives of eukaryotic mRNAs. Science. 2005; 309: 1514-1518. https://doi.org/10.1126/science.1111443 PMID: 16141059

67. Arvola RM, Weidmann CA, Tanaka Hall TM, Goldstrohm AC. Combinatorial control of messenger RNAs by Pumilio, Nanos and Brain Tumor Proteins. RNA Biol. 2017; 1-12. https://doi.org/10.1080/15476286.2017.1306168 PMID: 28318367

68. Lasko P. Posttranscriptional regulation in Drosophila oocytes and early embryos. Wiley Interdiscip Rev RNA. 2011; 2: 408-416. https://doi.org/10.1002/wrna.70 PMID: 21957026

69. Lehmann R. Germ Plasm Biogenesis—An Oskar-Centric Perspective. Curr Top Dev Biol. 2016; 116: 679-707. https://doi.org/10.1016/bs.ctdb.2015.11.024 PMID: 26970648

70. Herpin A, Adolfi MC, Nicol B, Hinzmann M, Schmidt C, Klughammer J, et al. Divergent expression regulation of gonad development genes in medaka shows incomplete conservation of the downstream regulatory network of vertebrate sex determination. Mol Biol Evol. 2013; 30: 2328-2346. https://doi.org/10.1093/molbev/mst130 PMID: 23883523

71. Sterky FH, Ruzzenente B, Gustafsson CM, Samuelsson T, Larsson N-G. LRPPRC is a mitochondrial matrix protein that is conserved in metazoans. Biochem Biophys Res Commun. 2010; 398: 759-764. https://doi.org/10.1016/j.bbrc.2010.07.019 PMID: 20633537

72. Kress C, Gautier-Courteille C, Osborne HB, Babinet C, Paillard L. Inactivation of CUG-BP1/CELF1 causes growth, viability, and spermatogenesis defects in mice. Mol Cell Biol. 2007; 27: 1146-1157. https://doi.org/10.1128/MCB.01009-06 PMID: 17130239

73. Gohil VM, Nilsson R, Belcher-Timme CA, Luo B, Root DE, Mootha VK. Mitochondrial and nuclear genomic responses to loss of LRPPRC expression. J Biol Chem. 2010; 285: 13742-13747. https://doi.org/10.1074/jbc.M109.098400 PMID: 20220140

74. Sasarman F, Brunel-Guitton C, Antonicka H, Wai T, Shoubridge EA, LSFC Consortium. LRPPRC and SLIRP interact in a ribonucleoprotein complex that regulates posttranscriptional gene expression in mitochondria. Mol Biol Cell. 2010; 21: 1315-1323. https://doi.org/10.1091/mbc.E10-01-0047 PMID: 20200222

75. Bratic A, Wredenberg A, Grönke S, Stewart JB, Mourier A, Ruzzenente B, et al. The bicoid stability factor controls polyadenylation and expression of specific mitochondrial mRNAs in Drosophila melanogaster. PLoS Genet. 2011; 7: e1002324. https://doi.org/10.1371/journal.pgen.1002324 PMID: 22022283
76. Sasarman F, Nishimura T, Antonicka H, Weraarpachai W, Shoubridge EA, LSFC Consortium. Tissue-specific responses to the LRPPRC founder mutation in French Canadian Leigh Syndrome. Hum Mol Genet. 2015; 24: 480-491. https://doi.org/10.1093/hmg/ddu468 PMID: 25214534

77. Boulanger G, Cibois M, Viet J, Fostier A, Deschamps S, Pastezeur S, et al. Hypogonadism Associated with Cyp19a1 (Aromatase) Posttranscriptional Upregulation in Celf1 Knockout Mice. Mol Cell Biol. 2015; 35: 3244-3253. https://doi.org/10.1128/MCB.00074-15 PMID: 26169831

78. Iwamatsu T. Stages of normal development in the medaka Oryzias latipes. Mech Dev. 2004; 121: 605-618. https://doi.org/10.1016/j.mod.2004.03.012 PMID: 15210170

79. Hong Y, Winkler C, Schartl M. Production of medakafish chimeras from a stable embryonic stem cell line. Proc Natl Acad Sci U S A. 1998; 95: 3679-3684. PMID: 9520425

80. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008; 36: W465-W469. https://doi.org/10.1093/nar/gkn180 PMID: 18424797

81. Pasquier J, Cabau C, Nguyen T, Jouanno E, Severac D, Braasch I, et al. Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database. BMC Genomics. 2016; 17: 368. https://doi.org/10.1186/s12864-016-2709-z PMID: 27189481

82. Gascuel O. BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997; 14: 685-695. https://doi.org/10.1093/oxfordjournals.molbev.a025808 PMID: 9254330

83. Zhang H, Gao S, Lercher MJ, Hu S, Chen W-H. EvolView, an online tool for visualizing, annotating and managing phylogenetic trees. Nucleic Acids Res. 2012; 40: W569-W572. https://doi.org/10.1093/nar/gks576 PMID: 22695796

84. Louis A, Muffato M, Roest Crollius H. Genomicus: five genome browsers for comparative genomics in eukaryota. Nucleic Acids Res. 2013; 41: D700-D705. https://doi.org/10.1093/nar/gks1156 PMID: 23193262

85. Lange M, Norton W, Coolen M, Chaminade M, Merker S, Proft F, et al. The ADHD-linked gene Lphn3.1 controls locomotor activity and impulsivity in zebrafish. Mol Psychiatry. 2012; 17: 855. https://doi.org/10.1038/mp.2012.119 PMID: 22918194

86. Grabher C, Wittbrodt J. Meganuclease and transposon mediated transgenesis in medaka. Genome Biol. 2007; 8 Suppl 1: S10. https://doi.org/10.1186/gb-2007-8-s1-s10 PMID: 18047687

87. Köprunner M, Thisse C, Thisse B, Raz E. A zebrafish nanos-related gene is essential for the development of primordial germ cells. Genes Dev. 2001; 15: 2877-2885. https://doi.org/10.1101/gad.212401 PMID: 11691838