Gonadal sex differentiation and reproduction are the keys to the perpetuation of favorable gene combinations and positively selected traits. In vertebrates, several gonad development features that differentiate tetrapods and ﬁshes are likely to be, at least in part, related to the water-to-land transition. The collection of information from basal sarcopterygians, coelacanths, and lungﬁshes, is crucial to improve our understanding of the molecular evolution of pathways involved in reproductive functions, since these organisms are generally regarded as “living fossils” and as the direct ancestors of tetrapods. Here, we report for the ﬁrst time the characterization of >50 genes related to sex differentiation and gametogenesis in Latimeria menadoensis and Protopterus annectens. Although the expression proﬁles of most genes is consistent with the intermediate position of basal sarcopterygians between actinopterygian ﬁsh and tetrapods, their phylogenetic placement and presence/absence patterns often reveal a closer afﬁnity to the tetrapod orthologs. On the other hand, particular genes, for example, the male gonad factor gsdf (Gonadal Soma-Derived Factor), provide examples of ancestral traits shared with actinopterygians, which disappeared in the tetrapod lineage. Key words: sex differentiation, gametogenesis, ovary, testis, evolution, Protopterus annectens, Latimeria menadoensis. Introduction sex determination modes, which either have a genetic basis or Sexual dimorphism is one of the most pervasive and diverse depend on other factors, for example, temperature popula- features of animal morphology, physiology, and behavior. The tion density, and visual cues. This divergence matches the process of gonad development involves complex and ﬁnely great variety of the initial triggers at the top of the gene net- tuned mechanisms that regulate gametogenesis, molecular work that direct the development toward a male or female sex determination, and gonad differentiation. Despite the phenotype, as well as the diversity of downstream acting generality of the phenomenon itself, sex determination, and components that implement and translate the initial signal differentiation mechanisms largely vary across metazoan taxa, into the development and differentiation of either testis or as the result of repeated, independent lineage-speciﬁc evolu- ovary from the bipotential embryonal gonad anlage (Herpin tion, and quick modiﬁcation of the underlying molecular et al. 2013). This situation is particularly evident in teleost pathways. Within vertebrates, mammals and birds share a ﬁshes, the largest group of extant chordates, where closely common genetic sex determining mechanism, which also dis- related species often display notably different sex determina- plays a remarkable degree of conservation at the molecular tion mechanisms (Kikuchi and Hamaguchi 2013). Despite this level. Others animal groups are characterized by highly diverse puzzling variability among species, some common principles The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact email@example.com 1430 Genome Biol. Evol. 10(6):1430–1444. doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE are emerging from the analysis of the growing informa- Materials and Methods tion available about metazoan gonad sex development Male and female P. annectens specimens were obtained from molecular networks, for example, the recurrent involve- a local ﬁsh importer. Two testis, three ovaries, and one brain ment of certain molecules and pathways as primary male and liver sample each from a male and a female specimen sex determinants. Interestingly, some of these common- were used for total RNA isolation and sequencing as described alities mark teleosts as clearly different from those of previously (Biscotti et al. 2016). other vertebrate classes and, in some instances, of all Pieces of the gonads used in the transcriptome analyzes other tetrapods. These peculiarities might even include were ﬁxed in 4% formaldehyde solution for 24 h at 4 C. the presence or absence of genes that are critical for go- Subsequently, the tissues were dehydrated, embedded in par- nadsex formationand gametogenesisinother branches afﬁn, and then serially sectioned at 5 mm thickness. The sec- of the vertebrate phylogenetic tree. The study of sex de- tions were counterstained with hematoxylin and eosin (HE). termination and differentiation mechanisms in basal sar- The sex and the gonadal maturation stage of the animals copterygians represents a key to unveil whether the used for the generation of transcriptomes were conﬁrmed by unique features displayed by teleosts exemplify an ances- histological analysis (ﬁg. 1). The male gonad of P. annectens tral situation or they correspond to teleost-speciﬁc special- was classiﬁed as anastomosing tubular testis type, where izations. Coelacanths and lungﬁshes represent the closest tubules do not terminate at the testis periphery (Parenti and extant relatives to the tetrapod lineage, even though the Grier 2004). This type is common in the more basal bony latter have been shown to be phylogenetically closer to ﬁshes, for example, salmonids, cyprinids, and lepisosteids the tetrapod stem (Biscotti et al. 2016; Irisarri and Meyer (Uribe et al. 2014). Both male gonadal tissues were at the 2016). Moreover, the morphological similarity between prereproductive stage (Schulz et al. 2010), with the MG1 the extant sarcopterygian ﬁsh and their fossil ancestors containing only few spermatogonias surrounded by Sertoli make them of outstanding scientiﬁc interest, considering cells. In MG6, the formation of individual tubules and their that they may also retain, to some extent, components of lumina are observed, and the germinal epithelium contains the molecular machinery of the last common ancestor of higher amounts of spermatogonia. The oogenesis in lung- tetrapods and ﬁsh. Due to the unique phylogenetic posi- ﬁshes is similar to the teleost (Chaves 1992; Lubzens et al. tion of lungﬁsh, their molecular studies are pivotal in un- 2010), and all three female gonads presented previtellogenic derstanding the evolution of vertebrates and the changes and vitellogenic oocytes. No postovulatory follicles were that accompanied the water-to-land transition. However, observed. molecular studies in lungﬁsh and coelacanths in particular Sequences corresponding to key sex development-related are severely hampered by the difﬁcult availability and genes were obtained by tBLASTn sequence homology maintenance of specimens in laboratory conditions. searches performed against the P. annectens and Latimeria Taking into account the impossibility of performing func- menadoensis transcriptomes (Canapa et al. 2012; Pallavicini tional studies in these species, the analysis of the gene et al. 2013; Biscotti et al. 2016) and the genomes of L. cha- expression proﬁles of sex differentiation and gametogen- lumnae (Amemiya et al. 2013), D. rerio (Howe et al. 2013), esis genes are the ﬁrst choice for obtaining insights about and O. latipes (Kasahara et al. 2007). The sequence of sex gonad determination and differentiation. differentiation and gametogenesis-related genes previously Due to the paucity of biological material available from characterized in other species were used as queries and sim- coelacanths, the only study published so far on this topic ilarity thresholds were set on a case-by-case basis. Based on reported the expression of a few genes in the testis of the expression pattern and the so far functional analyzes of Latimeria menadoensis (Forconi et al. 2013). Whole genome those genes, we classiﬁed them in male- or female-related studies in Dipnoi have been long hindered by the enormous genes. The retrieved transcript sequences were translated us- size and high repeat content of their genomes. However, the ing the Sequence Translation (https://www.ebi.ac.uk/Tools/st/), publication of a reference transcriptome for the West African permitting to identify UTRs and CDSs. The orthology of each lungﬁsh, Protopterus annectens (Biscotti et al. 2016)helped to transcript was conﬁrmed using NCBI BLAST (http://blast.ncbi. overcome this major obstacle, ﬁnally allowing RNA-Seq-based nlm.nih.gov/Blast.cgi). Each transcript was used as query gene expression studies in this species. against the database<all nonredundant GenBank CDS The aim of this study was to expand the knowledge on translationsþ PDBþ SwissProtþ PIRþ PRF> excluding envi- basal sarcopterygian sex differentiation and gametogenesis by ronmental samples from WGS projects. characterizing 51 genes related to this process and assessing Orthologous sequences used in the phylogenetic analyses their expression levels in gonads. We ﬁnd that the intermedi- were retrieved from GenBank (http://www.ncbi.nlm.nih.gov/) ate phylogenetic placement of lungﬁsh, between acanthop- or ENSEMBL (http://www.ensembl.org/index.html) (for acces- terygian ﬁsh and tetrapods, is matched by the patterns of sion numbers see supplementary table 1, Supplementary molecular evolution of the sex differentiation and gameto- Material online). Multiple alignments of the amino acid genesis gene network. Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1431 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE FIG.1.—Morphology of the gonads of Protopterus annectens.(A) Testes (MG1) revealing spermatogonias (asterisks) surrounded by Sertoli cells (arrowheads). (B) Testes (MG6) with completely formed tubules and lumina (Lu), and germinal epithelium containing spermatogonia (asterisks) and Sertoli cells (arrowheads). (C–E) Ovarian tissue (FG 2, 3, 4) containing pre–vitellogenic (PVO) and vitellogenic oocytes (VO). sequences were obtained with ClustalW (http://www.ge- to 0.75 and 0.98, respectively. Data from zebraﬁsh and me- nome.jp/tools/clustalw/) using default parameters. daka were retrieved from the PhyloFish project database Phylogenetic analyses were performed by Bayesian infer- (Pasquier et al. 2016), BioProject PRJNA255848 and ence using MrBayes (version 3.2) (Huelsenbeck et al. 2001). PRJNA255889. In medaka, the sequence of the fgf24 tran- Substitution models as deﬁned by posterior probabilities, sta- script, encoded by a gene present in a poorly assembled ge- tionarity, generations, sampling, burnin, and speciﬁc tree nomic region, was retrieved from the de novo assembled building parameters and rooting details are speciﬁed in the transcriptome. Raw reads counts were then converted to ﬁgure legend description of each tree. TPM values according to the formula suggested by Wagner The x rate (nonsynonymous/synonymous mutations) was et al. (2012) (supplementary ﬁg. 1, Supplementary Material obtained with KaKs_calculator (Zhang et al. 2006) according online). Considering that absolute TPM expression values, in to Goldman and Yang (1994) (supplementary table 2, absence of biological replicates, might be quite variable be- Supplementary Material online). Poorly aligned positions and tween individuals and are therefore inappropriate for inter- divergent regions of alignments were eliminated by Gblocks species comparisons, we assigned expression ranges with the (Castresana 2000; Talavera and Castresana 2007). following categories for each gene: absent< 10 The expression levels of the selected genes in the different TPM,þ between 10 and 200,þþ between 200 and tissues and species were determined as described by Biscotti 700,þþþ between 700 and 2,000,þþþþ over 2,000. et al. (2016). Brieﬂy, to ensure comparability among species, a Real-time quantitative PCRs of ﬁve genes (amh, ctnnb1, set of conserved orthologous genes present in L. menadoen- dmrt1, fgf9,and rspo1, for primers sequences see supplemen- sis, P. annectens,and D. rerio was selected as a reference for tary table 3, Supplementary Material online) were carried out the calculation of Transcripts Per Million (TPM). RNA-seq data in order to validate the transcriptome data. Reverse transcrip- from lungﬁsh (Biscotti et al. 2016) (male and female gonads, tion was done using RevertAid First Strand Synthesis kit brain, and liver), coelacanth (Pallavicini et al. 2013) (male go- (Fermentas) and random primers. The reactions were per- nad, liver, and skeletal muscle), zebraﬁsh (brain, liver, muscle, formed with SYBR Green reagent and ampliﬁcations were male and female gonads), and medaka (male and female detected with a Mastercycler ep realplex (Eppendorf). All gonads) were mapped to the reference nucleotide sequences results are averages of three independent PCR reactions with the CLC Genomics Workbench v.9 (Qiagen, Hilden, from cDNA preparations of three females and two males. Germany), setting length and similarity fraction parameters Transcript levels of target genes were normalized against 1432 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE Table 1 Identiﬁed Genes in Lungﬁsh and Coelacanth Growth Factors amh amhr2 fgf9 fgf16 fgf20 fgf24 gsdf rspo1 wnt4 ptch2 hhip dhh pdgfa pdgfb pdgfra pdgfrb sdf1 cxcr4 Steroidogenic Enzymes srd5a1 srd5a2 srd5a3 cyp11b cyp19a1 Transcription Factors dmrt1 dmrt3 dmrt6 foxl2 foxl3 sf1 dax1 sox3 sox5 sox8 sox9 sox10 wt1 gata4 lrpprc pax2 Meiosis Regulation Genes aldh1a1 aldh1a2 aldh1a3 cyp26a1 cyp26b1 cyp26c1 stra8 Follistatin and b Catenin fst ctnnb1 Steroid Hormone Receptor Genes ar esra esrb Underlined genes were not found in lungﬁsh. the lungﬁsh eukaryotic translation initiation factor 4 gamma 1 growth factor-b (TGF-b) superfamily, which causes regression (eif4g1)gene. of the Mu ¨ llerian ducts during testicular differentiation in mammals. The high expression values of amh observed in male gonads of several vertebrates have been linked to the Results and Discussion proliferation and differentiation of spermatogonia (Skaar et al. In this paper, the presence of 51 genes known to be involved 2011; Hu et al. 2015). amh is also expressed in females, where in gonad sex differentiation and gametogenesis development it is restricted to the granulosa cells of ovarian follicles (Piprek, in other vertebrates was assessed in the transcriptomes of the Pecio, Laskowska-Kaszub, Kubiak, et al. 2013). Amh mainly lungﬁsh P. annectens. This analysis allowed the identiﬁcation acts through the Amhr2 signaling system in the development of 44 transcripts (table 1), all showing a complete CDS (sup- of gonadal tissue. Orthologous genes have been identiﬁed in plementary table 4, Supplementary Material online). The ﬁve several teleosts (Klu ¨ ver et al. 2007; Li et al. 2015; Pfennig et al. genes used for validation of the transcriptome by qRT-PCR 2015), despite the absence of Mu ¨ llerian ducts in these organ- showed mRNA levels consistent with the transcriptome data isms. In addition, Amh or its type II receptor act as Y-linked of the respective genes (supplementary table 5, male sex determining genes in some teleost species (Hattori Supplementary Material online). For comparison’s sake we et al. 2013; Kikuchi and Hamaguchi 2013; Herpin and Schartl also studied the sarcopterygian ﬁsh Latimeria, identifying 17 2015). On the other hand, the Mu ¨ llerian duct develops in new genes (table 1 and supplementary table 6, female coelacanths and lungﬁsh (Hildebrand and Goslow Supplementary Material online) in addition to those investi- 2001) and in a basal group of actinopterygian ﬁshes, the gated in our previous work (Forconi et al. 2013). The expres- sturgeons (Wrobel 2003). In lungﬁsh and in the teleosts, sion levels (absent, low, intermediate, high) of the identiﬁed amh is preferentially expressed in male gonads and a robust genes were compared between the two basal sarcopterygians expression of this gene could be also observed in the testis of andwiththe teleosts Danio rerio and Oryzias latipes to provide Latimeria (table 2). the basis for a better understanding of the gene regulatory Like coelacanth and most other ﬁsh, lungﬁsh possess a network governing gametogenesis and sex differentiation bona ﬁde ortholog of amhr2 (supplementary ﬁg. 2, and their evolution. Supplementary Material online). Unexpectedly, this gene is missing in the zebraﬁsh genome. This appears to be species-speciﬁc feature, as amhr2 orthologs are present in Growth Factors and Their Receptors other teleost genomes. In medaka, amhr2 is highly expressed Anti-Mu ¨ llerian hormone (Amh), also known as Mu¨llerian in testis (supplementary ﬁg. 1, Supplementary Material on- duct-inhibitory substance, is a member of the transforming line), and mutation of this gene leads to overproliferation of Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1433 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE Table 2 Growth Factors and Receptors Generally Related to Male Sex Determination, Differentiation, and Reproductive Functions Lungﬁsh Coelacanth Zebraﬁsh Medaka Testis Ovary Testis Ovary Testis Ovary Testis Ovary amh þ þ N.R. þþþþ þþ þ amhr2 þ þ N.R. Absent þþ þ gsdf þ þ N.R. þþþþ þ þþþþ þþ fgf9 þ N.R. Absent Absent fgf16 þ N.R. fgf20a Absent N.R. fgf20b þ fgf24 Absent N.R. þþþ sdf1a þþþþ þ þ N.R. þþþ þ þ sdf1b þ þþ þ cxcr4a þþþþ þ þþ N.R. þþþþ cxcr4b þþ þþ a a pdgfaa þ þ þ N.R. þ pdgfab þ pdgfra þ þ N.R. þþþ pdgfb þ þ N.R. þþ pdgfrb þþ N.R. þþþ Approximately 3-fold higher in ovary than testis. the germ cells in both testis and ovary, showing an important gene was already present in the common ancestor of role of Amh signaling in germ cell proliferation in adult gonad Chondrichthyes and Osteichthyes. (Morinaga et al. 2007). This gene is expressed in the testis of The high expression of gsdf in P. annectens testis (table 2)is Latimeria, whereas in the lungﬁsh its expression is only found in line with data previously observed in coelacanth (Forconi in female gonads (table 2). This ﬁnding supports the presence et al. 2013) conﬁrming the important role of this gene for of female expression bias of amhr2 in lungﬁsh, like in tetra- male gonad development also in basal sarcopterygians. pods (Fagerberg et al. 2014; Jansson et al. 2016), but not in Furthermore, this pattern is in accordance with the well- Latimeria,where amhr2 is exclusively expressed in testis, like in described dimorphic expression of this gene in other ﬁsh spe- teleost ﬁsh, where a reverse organ expression pattern has cies (Kaneko 2015). gsdf can even act as master sex deter- been observed (Lubieniecki et al. 2015). mining gene (Myosho et al. 2012; Zhu et al. 2016)or as a Gonadal Soma-Derived Factor (GSDF) is another member male sex initiator, being activated by Dmrt1/Dmy (Zhang et al. of the TGF-b superfamily. This gene is particularly interesting 2016). because of its crucial function in ﬁsh testis development In summary, our results point out that a gsdf gene with (Zhang et al. 2016). Prior to our identiﬁcation of a homolo- testis-speciﬁc function was retained in the sarcopterygian gous sequence in Latimeria (Forconi et al. 2013), gsdf was ancestors of tetrapods, only to be lost during the transition hypothesized to be a teleost innovation. Here, we report a from water to land. transcript encoding for gsdf also in P. annectens,whose R-spondin1 (Rspo1) is a regulator of the WNT/b-catenin orthology assessment is strongly supported by phylogenetic signaling pathway. This gene has been characterized as a fe- inference (ﬁg. 2A). Moreover, as a distinctive feature, Gsdfs male sexual regulator in mammals, as it suppresses the male lack a glycine residue in the CxGxC motif, different from Inha axis and maintains oocyte survival by positively regulating and Amh proteins which have this amino acid (ﬁg. 2B). The WNT4 signaling. In tetrapods and several teleosts, rspo1 presence of glycine in this motif is essential for heterodimer shows a sexually dimorphic expression pattern with higher formation of glycoprotein hormones (Kinoshita et al. 2006). values in the female gonad (Zhou et al. 2016), suggesting The functional meaning of this amino acid change for Gsdf is that Rspo1/Wnt/b-catenin signaling is an ancient pathway unknown, due to the lack of information about the signaling of ovarian differentiation in vertebrates. We found a quite pathway of this gene. The identiﬁcation of a gsdf sequence in high expression for rspo1 in lungﬁsh testis (table 3), and P. annectens implies its subsequent loss in the tetrapod line- much lower levels in ovaries. In zebraﬁsh, rspo1 is equally age. Moreover, the additional identiﬁcation of a gsdf se- expressed in ovary and testis, while in medaka no expression quence in the cartilaginous ﬁsh C. milii evidences that this was detected. rspo1 is also expressed in the ovary and testis of 1434 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE FIG.2.—Gsdf phylogeny. (A) Phylogenetic tree of Gsdf, Amh, and Inha amino acid sequences. Bayesian inference: 2,000,000 generations, sampling every 100, Jones substitution model, stationarity deﬁned as when the average SD of split frequencies approaching 0.007, burn-in set to 5,000, midpoint rooting. Numbers close to nodes indicate posterior probability values. Only values >95 were reported. The sequences in bold were obtained in this work. (B) Multiple alignment of Gsdf, Amh, and Inha amino acid sequences showing the CxGxC conserved motif and the lack of the glycine residue, which is diagnostic for Gsdf proteins. Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1435 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE Table 3 Growth Factors and Receptors Generally Related to Female Sex Determination, Differentiation, and Reproductive Functions Lungﬁsh Coelacanth Zebraﬁsh Medaka Testis Ovary Testis Ovary Testis Ovary Testis Ovary dhh þþ N.R. þ þþ ptch2 N.R. þ þ hhip þ N.R. þþ þ þ þ ctnnb1 þþþ þþþþ þþþ N.R. þþþþ þþþþ þþþ þþ fst þ N.R. þ þ a a wnt4a þ N.R. þ þ a a wnt4b þ þ þ rspo1 þþ þ þ N.R. þþ Approximately 2- to 3-fold higher in testis than ovary. other teleosts, suggesting that it might be involved in both of fgf9 in C. milii dates back the origin of this gene to the ovarian and testicular development and spermatogenesis (Wu common ancestor of Chondrichthyes and Osteichthyes, im- et al. 2016). In other vertebrates, RSPO1 is widely expressed plying its loss in the lineage of Actinopterygii. Neither fgf9 nor (Han et al. 2015) and its transcript has been also detected in the other members of the FGF9/16/20 subfamily are adult human testis (Parma et al. 2006). Our ﬁndings suggest expressed in testis of lungﬁsh and coelacanth (table 2). In that theroleof rspo1 in lungﬁsh male gonads is not con- lungﬁsh, we veriﬁed the absence of fgf20 both at the tran- nected to the action of wnt4, asthisgene is poorly expressed scriptome level. Moreover, fgf9 is strongly expressed in lung- in lungﬁsh and coelacanth, while both wnt4 co-orthologs are ﬁsh ovary (table 2). This may indicate that the important well expressed in the gonads of zebraﬁsh. Medaka, on the function of fgf9 for testis development is an innovation ac- other hand, presents relatively low expression of wnt4b in quired after the transition from water to land. testis. To evaluate the possible involvement of another member ß-catenin (Ctnnb1) is the critical molecule that transmits of the FGF gene family as a functional homolog of the tetra- the activity of the female determining pathway signals of pod Fgf9, as recently proposed by Leerberg et al. (2017),we Wnt4 and R-spondin to the nucleus to activate effector genes. took a closer look at fgf24, a member of the FGF8/17/18/24 High expression values were recorded for ctnnb1 in all four family (ﬁg. 4). The isolation of this gene in the elephant shark species (table 3). While no strong sex bias in expression is supported its presence in the common ancestor of gna- present in zebraﬁsh and medaka, the lungﬁsh ovary shows thostomes (Jovelin et al. 2010). Within the sarcopterygian considerably higher levels than testis, in contrast with the ex- lineage, fgf24 could be only identiﬁed in Latimeria pression pattern of wnt4 and rspo1 (table 3). (Amemiya et al. 2013) but not in tetrapods (Jovelin et al. Follistatin (Fst) encodes a TGFß superfamily binding protein 2010). In addition, no orthologous sequence to fgf24 that acts downstream of WNT4 to maintain germ cell survival could be retrieved from the P. annectens transcriptome, in the cortical domain of the mammalian ovary (Yao et al. which might either indicate the lack of expression of this 2004). This gene is present in both basal sarcopterygians. The gene, or its loss. If conﬁrmed at the genome level, this lack expression of fst in lungﬁsh, in line with the absence of ﬁnding might indicate that fgf24 was lost in the common detectable expression for wnt4 (table 3), may be connected to ancestor of lungﬁsh and tetrapods. This gene is expressed the prereproductive stage of the female gonads. in the early somatic and in adult gonads of zebraﬁsh Fibroblast Growth Factors (Fgf) are a family of small pro- Leerberg et al. (2017) (table 2) and was suggested as a teins that have multiple functions in regulating cellular devel- plausible candidate to substitute for fgf9 function in tele- opment and growth. Numerous studies in mice have shown osts. In addition, fgf24 is expressed only in testis of adult that FGF9 plays an important role in testis development, cre- medaka. ating a positive feedback loop with Sox9 and antagonizing Stromal-cell-derived factor 1 (Sdf1, also known as Cxcl12) Wnt4, an ovarian factor (Jameson et al. 2012). No orthologs is a chemokine involved in a wide variety of body functions. It of fgf9 have been identiﬁedinactinopterygian ﬁshsofar (Itoh binds to CXC receptor 4 (Cxcr4), a G-coupled transmembrane and Konishi 2007). The identiﬁcation of fgf9 in Latimeria has glycoprotein located on the surface membranes of primordial previously suggested that this gene was present in the com- germ cells. The interaction between Sdf1 and its receptor mon ancestor of Sarcopterygii (Forconi et al. 2013), a hypoth- regulates migration and proliferation of various cell types, esis which is now further supported by the retrieval of this which importantly comprise also teleosts primordial germ sequence also in P. annectens (ﬁg. 3). Moreover, the presence cells. Previous studies assigned a central role to Sdf1 and 1436 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE FIG.3.—Fgf9/16/20 phylogeny. Phylogenetic analysis of Fgf9, Fgf16, Fgf20. Bayesian inference: 1,000,000 generations, sampling every 100, Jones substitution model, stationarity deﬁned as when the average SD of split frequencies approaching 0.0055, burn-in set to 2,500, midpoint rooting. Numbers close to nodes indicate posterior probability values. Only values >95 were reported. The sequences in bold were obtained in this work. CXCR4 in the development of embryonic and postnatal go- Steroidogenic Enzymes nad (Yang et al. 2013; Westernstro ¨ er et al. 2014; Fernandez Steroid-5-Alpha-Reductase Alpha Polypeptides 1, 2, 3 et al. 2015). In human adults, the expression of SDF1 was (Srd5a1, Srd5a2, Srd5a3) are proteins involved in the conver- observed in Sertoli cells. The expression of Cxcr4 in adult sion of testosterone to dihydrotestosterone (DHT), a very po- mouse testis (Yang et al. 2013) further suggests a role for tent androgen in teleosts ﬁshes and tetrapods. Unexpectedly, the ligand/receptor in spermatogenesis. Both genes are pref- we detected higher expression of srd5a1 in lungﬁsh ovary, in erentially expressed in lungﬁsh male gonads (table 2), at par- contrast to its absence in zebraﬁsh and the low expression in ticularly high levels in the more immature testis sample. In female medaka (supplementary table 7, Supplementary addition, the expression of these genes has been also recently Material online). Srd5a1 is predominantly expressed in the reported in both male and female gonads of the teleost female gonads of frogs but also in rat (Lephart et al. 1992) Oreochromis mossambicus (Amat-Fernandez et al. 2017). and human (Haning et al. 1996; Mahendroo and Russell Platelet-derived growth factors (Pdgf) and their receptors 1999), consistently with our observations in lungﬁsh (supple- are involved in regulating the development of various tissues mentary table 7, Supplementary Material online). This ﬁnding (Miura et al. 2002). In particular, they play a fundamental role suggests that Srd5a1 might confer the ability of metabolizing in the development of gonadal tissue in mammalian fetal and testosterone into DHT in the sarcopterygian lineage, while a adult testis (Li et al. 2016). In zebraﬁsh and medaka, this gene different isozyme is likely used in teleost ﬁsh. is preferentially expressed in testis (table 2). However, in both Cytochrome P450 19A1 (Cyp19a1 or aromatase) cata- sarcopterygians, the expression of pdgfs and their receptors in lyzes the ﬁnal step in the biosynthesis of estrogens, which testis were quite low. The lungﬁsh pdgfa and pdgfrb ortho- have a key role in ovarian differentiation and maintenance logs unexpectedly showed the highest expression levels in in ﬁsh, amphibians, reptiles, and birds (Bruggeman et al. ovary (table 2). Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1437 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE FIG.4.—Fgf8/17/18/24 phylogeny. Phylogenetic analysis of Fgf8, Fgf17, Fgf18, Fgf24. Bayesian inference: 2,000,000 generations, sampling every 100, Jones substitution model, stationarity deﬁned as when the average SD of split frequencies approaching 0.0062, burn-in set to 5,000, midpoint rooting. Numbers close to nodes indicate posterior probability values. Only values >95 were reported. The sequences in bold were obtained in this work. 2002; Dumond et al. 2008; Lubieniecki et al. 2015). The detect 11-ketotestosterone in Australian lungﬁsh (Joss et al. expression patterns observed in the present study is con- 1996). sistent with a conserved function for this gene also in basal sarcopterygians (supplementary table 8, Steroid Hormone Receptors Supplementary Material online). Cytochrome P450 11c1 (Cyp11c1). In most teleost ﬁsh the The expression of steroidogenic enzymes is usually connected predominant androgen is 11-ketotestosterone, which is pro- to the production of biologically active estrogens and andro- duced from testosterone by the action of the cyp11c1 gene. gens. The action of these steroid hormones requires the pres- However, tetrapod cyp11B homologs function predominantly ence of the cognate receptors, whose transcription is often in adrenal steroid biosynthesis and are not expressed in testis further up-regulated by the presence of hormones in a posi- (Schiffer et al. 2015). A gene pertaining to the CYP11 sub- tive feedback loop. We found a very unexpected expression family is present in lungﬁsh as well as in the coelacanth (sup- pattern for these receptors in lungﬁsh. Indeed, while the an- plementary ﬁg. 3, Supplementary Material online). The drogen receptor (ar) was expressed in the mature testis of absence of expression in the gonads of these two species is coelacanth and teleost, no expression could be detected in in line with the results of a previous study, which failed to lungﬁsh male gonad nor in any other tissue (supplementary 1438 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE Table 4 Transcription Factors Generally Related to Male Sex Determination, Differentiation, and Reproductive Functions Lungﬁsh Coelacanth Zebraﬁsh Medaka Testis Ovary Testis Ovary Testis Ovary Testis Ovary dmrt1 þþ þ þþ N.R. þþþþ þ þþ þ dmrt3 N.R. þ þ dmrt6 þþ þ þþþ N.R. Absent Absent sox3 N.R. þ þ sox5 þ N.R. þþþ þ sox8a þ N.R. þþþ sox8b þþþþ þ a a sox9a þ þ þþ N.R. þ þ þ b b sox9b þ þ þ sox10 þ N.R. þ gata4 þþ N.R. þ þ þ lrpprc þþ þþþ þþ N.R. þþþ þþþ þ þþ sf1 þþþ þ þ N.R. þþþþ þþþþ þþ þ wt1a þþ N.R. þþ þ þ þ wt1b þþ pax2a þ N.R. þ pax2b þ In ovary barley above background, but robust expression in testis. Approximately 3-fold higher in ovary than testis. table 7, Supplementary Material online). Whether this ab- Although no dmrt6 gene has been annotated yet in any of the sence is linked to the nonreproductive stage of the gonads fully sequenced teleost genomes, a recent publication or to the use of another nuclear receptor for sensing andro- reported dmrt6-like sequences in spotted gar, cichlids, and gen levels in lungﬁsh remains an interesting question to be siluriforms. This dmrt6-like gene was expressed in spermato- approached in the future. Teleost esrs were also preferentially cytes and its function was shown to be involved in spermato- expressedintestis(supplementary table 8, Supplementary genesis in tilapia (Zhang, Wang, et al. 2014). The loss of this Material online). On the other hand, mammalian ESRsare gene in some teleost lineages could represent another aspect produced at very low levels in testis (Cavaco et al. 2009). of the extreme variability of the genetic network involved in The exceptionally high levels identiﬁed in lungﬁsh male sex determination and differentiation in ﬁsh. Overall, the pres- gonads may be related to the premature status of the tissues ence and expression of dmrt6 in basal sarcopterygians sup- available for this study (supplementary table 7, Supplementary ports an ancient and conserved function for this gene in testis Material online). development and function. The Steroidogenic Factor 1 (Sf1) was mainly expressed in lungﬁsh and coelacanth testis, while in zebraﬁsh it showed a Transcription Factors broader tissue distribution (table 4). Sf1 is expressed in both sexes in mice fetal gonads, but in later developmental stages it Doublesex- and mab-3-related transcriptional factor 1 continues to be expressed only in testes, where it is activated (Dmrt1), together with Amh and Sox9, is important for ver- and maintained by Wilms tumor 1 (WT1). On the other hand, tebrate testicular differentiation. Its expression of dmrt1 in Sf1 is down-regulated by Forkhead box L2 (FOXL2) in ovary lungﬁsh (table 4) is in agreement with the dimorphic expres- (Takasawa et al. 2014). The low expression of sf1 observed in sion pattern observed from ﬁsh to mammals (Websteretal. lungﬁsh ovary, opposed to the high expression detected in 2017) and suggest a major role in adult testis function. zebraﬁsh female gonad, is in line with the mammalian In contrast to other vertebrates (Bratus and Słota 2009), no pattern. trace of dmrt3 expression was detected in the two basal sar- SRY-related HMG box (SOX) family proteins are transcrip- copterygians (table 4). On the other hand, the expression of tion factors involved in many developmental processes. In dmrt6 in lungﬁsh and coelacanth testis is intriguing (table 4). particular, sox3, sox5, sox8, sox9,and sox10 have a function In mice, this gene is likely involved in directing the transition in gonad development. SOX3 belongs to subfamily B and from spermatogonia to spermatocytes (Zhang, Murphy, et al. is the ancestor of SRY (Marshall-Graves 2008). Although 2014). This study allows to extend the taxonomic range dis- sox3 is not required for sex determination in mammals tribution of dmrt6 from tetrapods to basal sarcopterygians. Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1439 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE (Weiss et al. 2003), it is the master male sex determining gene in the ﬁsh Oryzias dancena (Takehana 2014). sox5, pertaining to the subfamily D, is a regulator of male germ cells in medaka (Schartl et al. 2018). It is expressed in mice postnatal testis and, together with Sox9 and Sox13, it may play a role in spermatogenesis (Daigle et al. 2015). The other three genes belong to subfamily E: SOX9 is activated by SRY in mammals and is required for sex determination in tetrapods (Jakob and Lovell-Badge 2011); whereas Sox8 is activated by SOX9 and acts in a positive feedback loop on Sox9 in testis differentia- tion, as shown in mice. All ﬁve sox genes are present in the lungﬁsh transcriptome. In Latimeria, in addition to the sox8/9/ 10 sequences identiﬁed in our previous work (Forconi et al. 2013), we can report the identiﬁcation of sox3 (from the ge- nome) and sox5 (from transcriptome data). The orthology status of the isolated sequences was assessed by phylogenetic analysis (supplementary ﬁg. 4, Supplementary Material online) FIG.5.—Stra8 phylogeny. Phylogenetic analysis of Stra8. Bayesian that assigned each sequence to the corresponding sox gene inference: 1,000,000 generations, sampling every 100, Jones substitution clade. model, stationarity deﬁned as when the average SD of split frequencies The tissue distribution of sox3 transcripts conﬁrms that the approaching 0.0033, burn-in set to 2,500, midpoint rooting. Numbers almost exclusive role of this gene in brain is also retained in close to nodes indicate posterior probability values. Only values >95 basal sarcopterygians (table 4). Together with SOX1 and were reported. The sequences in bold were obtained in this work. SOX2, the other two paralogs of the B1 group of the SOX family, SOX3 plays a key role in central nervous system devel- coelacanth, lungﬁsh, and zebraﬁsh. Expression of foxl3 in opment in mouse (Rogers et al. 2013). sox8 was neither the gonad was observed only in medaka, indicating a diver- expressed in Latimeria (Forconi et al. 2013) nor in lungﬁsh gent role of this gene between different vertebrates. This testis (table 4). In mammals, the action of SOX8 in the gene has not been identiﬁed in placental mammals and in FGF9/SOX9 interaction loop is important for testis mainte- amphibians (Bertho et al. 2016) and, except from the study in nance (Barrionuevo et al. 2009). The lack of expression of medaka, its gonadal function remains unclear. both sox8 and fgf9 hints that this interaction loop, which is crucial in mammalian sex determination, is not yet in place in Meiosis Regulation Genes basal sarcopterygians. Both sox9 and wt1 are expressed in lungﬁsh and coelacanth testis, in line with their known pat- Meiosis is one of the most important processes in sexual re- tern of distribution in tetrapods. production, as it produces haploid gametes. However, regu- Forkhead transcription factors (Fox) are involved in many lation of the meiotic program is still insufﬁciently understood biological processes, including the regulation of cell prolifera- in most species. In mammals, the morphogen retinoic acid tion and differentiation, tissue development, and various met- triggers meiotic entry through up-regulation of the abolic processes. Foxl2 plays a key role in the early Stimulated by retinoic acid gene 8 (Stra8). This gene encodes development of vertebrate female gonads and, in adults, a retinoic acid-responsive protein that plays a critical role in has an ovarian function. In particular FOXL2 suppresses the regulation of meiotic initiation during spermatogenesis Dmrt1 for maintaining female fate (Matson et al. 2011) and oogenesis (Feng et al. 2014). Commonly present in tetra- and, together with WNT4, cooperates in regulating Fst ex- pods, this gene is absent in some teleosts species, including pression during ovarian development (Garcia-Ortiz et al. zebraﬁsh (Rodrıguez-Marı et al. 2013). A more comprehen- 2009). It regulates estrogen synthesis via up-regulation of sive analysis of the evolution of stra8 suggests that this gene Cyp19A1a (Pannetier et al. 2006). In medaka, the paralogous was lost in Acanthomorpha and, independently, in the gene foxl3 encodes a protein that may play a role in sexual Cypriniform lineage within Actinopterygii (Pasquier et al. fate decision in germ cells by suppressing spermatogenesis 2016). We can however report the presence of this gene in (Nishimura et al. 2015). While both these genes are present the two basal sarcopterygian species. Phylogenetic analysis in the genome of L. chalumnae, only foxl2 was found in the P. positioned the lungﬁsh and coelacanth sequences at a basal annectens transcriptome. The sequences identiﬁed were used position compared with tetrapods (ﬁg. 5). for phylogenic analyzes (supplementary ﬁg. 5, Supplementary The level of the retinoic acid as a meiosis signal is deter- Material online). foxl2 was expressed in lungﬁsh female mined by the antagonistic action of enzymes involved in its gonads (supplementary table 9, Supplementary Material on- synthesis and degradation: the enzyme ALDH1a2, a member line). No expression of foxl3 was however detected in of the retinaldehyde dehydrogenase family, is responsible for 1440 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE FIG.6.—Summarizing scheme of sex-related genes lost in the main vertebrate lineages. The scheme shows the sex genes lost in tetrapods, lungﬁsh, coelacanth, and teleosts. The colours are referred to the functional groups, respectively. its synthesis, while CYP26B1 encodes a major retinoic acid- (supplementary table 2, Supplementary Material online). All degrading enzyme. In mice, down-regulation of Cyp26b1 in genes appear to evolve under purifying selection (Ka/Ks<1) femalesleadstohighlevels of retinoic acid, whichin turn with the exception of the x rate obtained between the two induce the expression of the premeiotic marker Stra8, thereby gsdf co-orthologs in D. rerio. allowing female germ cells to enter meiosis. Cyp26b1 acts as a Our analysis of lungﬁsh ovary and testis transcriptomes and meiosis-inhibiting factor in mammals, birds, and amphibians the comparison with similar data sets from both zebraﬁsh and (Bowles and Koopman 2010; Piprek, Pecio, Laskowska- medaka, and with the testis transcriptome of coelacanth Kaszub, Kloc, et al. 2013; Yu et al. 2013), while in teleosts revealed a situation that reﬂects the peculiar phylogenetic this role is exerted by cyp26a1 (Adolﬁ et al. 2016). position of the Dipnoi. We identiﬁed some genes whose ex- We found high expression of cyp26b1 and stra8 in pression pattern in lungﬁsh and, as far as testis is concerned, lungﬁsh and coelacanth male gonads suggesting that also in Latimeria is more similar to teleosts than tetrapods. As basal sarcopterygians follow the tetrapod pattern (sup- an example, gsdf, which is an important male determining plementary table 10, Supplementary Material online). gene in teleosts, but which is missing in the tetrapod lineage, However, signiﬁcant differences were observed in the was present and highly expressed in testis of both lungﬁsh expression of aldh1as: aldh1a1 showed higher expres- and coelacanth. sion in lungﬁsh male gonads, and is therefore the only On the other hand, a number of other genes present in gene responsible for retinoic acid synthesis in coelacanth lungﬁsh display a behavior more similar to tetrapods than testis (supplementary table 10, Supplementary Material teleost ﬁshes, for example, srd5a1, which was highly online). Such lineage speciﬁc subfunctionalization of ret- expressed in lungﬁsh ovary, like in tetrapods, or b-catenin, inoic acid metabolizing enzymes has also been reported the transcription factor readout of the female gonad spe- in a comparative study among different teleost ﬁsh ciﬁc tetrapod signaling cascade, which was also highly (Canestro et al. 2009). cyp26c1 was not found in the expressed in lungﬁsh ovary. Intriguingly, the upstream ef- lungﬁsh transcriptome, marking a possible lineage- fector R-spondin1 does not follow this pattern. At the speciﬁc loss. same time wnt4, which is an ovary speciﬁc gene in tetra- pods and an upstream regulator of ß-catenin, is highly expressed in teleost testis and did not show signiﬁcant Conclusions expression in lungﬁsh and coelacanth gonad samples. Our analysis of sex differentiation and gametogenesis genes Sf1 is a male speciﬁc gene in mammals and lungﬁsh go- in lungﬁsh and coelacanth revealed a complex picture. None nad, which was also highly expressed in Latimeria testis. of the gene-encoded protein products showed any unusual stra8 is absent from many major teleost lineages, but it is structural feature, like missing or additional domains or unex- present and highly expressed in the testis of basal sarcop- pectedly high sequence divergence. We calculated the x rates terygians, following the tetrapod proﬁle. Dmrt6, a gene for genes that were of special interest due to their appear- which has an important function in mammalian spermato- ance/disappearance at the split of sarcopterygians and actino- genesis and which is only seldom found in teleost pterygians (fgf9, foxl3, gsdf,and stra8), in lungﬁsh, genomes, showed high expression in lungﬁsh and coela- coelacanth and other representative vertebrates canth testis, indicating functional conservation. The Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1441 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE reasons underpinning the possible loss of such an impor- the conception of the study and the experimental design and tant gene in most teleosts—similar but independent from have given ﬁnal approval for the version to be published. stra8—are totally unclear. Finally, some lungﬁsh and coelacanth genes follow nei- Literature Cited ther the teleost nor the tetrapod scheme. Fgf9 is an im- Adolﬁ MC, Herpin A, Regensburger M, Sacquegno J, Waxman JS. 2016. portant male sex determining gene in tetrapods, which is Retinoic acid and meiosis induction in adult versus embryonic gonads absent in teleosts. Our analysis demonstrates that fgf9 has of medaka. Sci Rep. 6:34281. been speciﬁcally lost in the teleost lineage and the ab- Amat-Fernandez J, et al. 2017. Molecular characterization of sdf1 and cxcr4 in the Mozambique tilapia, Oreochromis mossambicus.Anim sence of expression in the lungﬁsh and coelacanth testis Reprod Sci. 176:51–63. may indicate that its role as a male gene was acquired Amemiya CT, et al. 2013. The African coelacanth genome provides later along the evolution of the tetrapod lineage. Sox8 is insights into tetrapod evolution. Nature 496(7445):311–316. a crucial component of the Sox9-driven regulatory loop Barrionuevo F, et al. 2009. Testis cord differentiation after the sex deter- that determines testis development and Sertoli cell main- mination stage is independent of Sox9 but fails in the combined ab- sence of Sox9 and Sox8. Dev Biol. 327(2):301–312. tenance in mice. Despite the high expression of sox8 in Bertho S, et al. 2016. Foxl2 and its relatives are evolutionary conserved zebraﬁsh testis, the lack of expression in coelacanth, lung- players in gonadal sex differentiation. Sex Dev. 10(3):111–129. ﬁsh, and medaka testis suggests that this regulatory loop Biscotti MA, et al. 2016. The lungﬁsh transcriptome: a glimpse into mo- is a tetrapod innovation. lecular evolution events at the transition from water to land. Sci Rep. Our phylogenetic analyses (ﬁg. 6) pointed out that several 6:21571. Bowles J, Koopman P. 2010. Sex determination in mammalian germ cells: genes that are important for sex differentiation and gameto- extrinsic versus intrinsic factors. Reproduction 139(6):943–958. genesis in one branch of the vertebrate tree of life can be Bratus A, Słota E. 2009. Comparative cytogenetic and molecular studies of totally absent in another branch, further suggesting that other DM domain genes in pig and cattle. Cytogenet Genome Res. 126(1– related genes may eventually take over their function, as ex- 2):180–185. empliﬁed by the case of fgf9 and fgf24. In all cases however, Bruggeman V, Van AP, Decuypere E. 2002. Developmental endocrinology of the reproductive axis in the chicken embryo. Comp Biochem Physiol the presence/absence of many sex differentiation and game- A Mol Integr Physiol. 131(4):839–846. togenesis genes across different animal groups appears to be Canapa A, et al. 2012. Composition and phylogenetic analysis of vitello- mostly due to the lineage-speciﬁc loss of ancestral genes genin coding sequences in the Indonesian coelacanth Latimeria mena- rather than to the emergence of novel genes. This implies doensis. J Exp Zool B Mol Dev Evol. 318(5):404–416. that the majority of sex differentiation and gametogenesis ~ Canestro C, Catchen JM, Rodrıguez-MarıA, Yokoi H, Postlethwait JH. 2009. Consequences of lineage-speciﬁc gene loss on functional evo- genes were already present that at the base of vertebrate lution of surviving paralogs: aLDH1A and retinoic acid signaling in evolution, after the 1R and 2R whole genome duplication vertebrate genomes. PLoS Genet. 5(5):e1000496. events. Castresana J. 2000. Selection of conserved blocks from multiple align- ments for their use in phylogenetic analysis. Mol Biol Evol. 17(4):540–552. Supplementary Material Cavaco JE, Laurentino SS, Barros A, Sousa M, Socorro S. 2009. Estrogen Supplementary data areavailableat Genome Biology and receptors alpha and beta in human testis: both isoforms are expressed. Syst Biol Reprod Med. 55(4):137–144. Evolution online. Daigle M, Roumaud P, Martin LJ. 2015. Expressions of Sox9, Sox5, and Sox13 transcription factors in mice testis during postnatal develop- ment. Mol Cell Biochem. 407(1–2):209–221. Acknowledgments Dumond H, et al. 2008. Sexual development of the urodele amphibian Pleurodeles waltl. Sex Dev. 2(2):104–114. This work was supported by grants supplied by the Deutsche Fagerberg L, et al. 2014. Analysis of the human tissue-speciﬁc expression Forschungsgemeinschaft (SCHA 408/10-1; 12-1) to M.S. This by genome-wide integration of transcriptomics and antibody-based publication was funded by the German Research Foundation proteomics. Mol Cell Proteomics 13(2):397–406. (DFG) and the University of Wuerzburg in the funding pro- Feng CW, Bowles J, Koopman P. 2014. Control of mammalian germ cell entry into meiosis. Mol Cell Endocrinol. 382(1):488–497. gramme Open Access Publishing. We thank the PhyloSex con- Fernandez JA, et al. 2015. Primordial germ cell migration in the yellowtail sortium for making the zebraﬁsh and medaka RNA-seq data kingﬁsh (Seriola lalandi) and identiﬁcation of stromal cell-derived factor available. M.A.B. and A.C. were involved in the identiﬁcation 1. Gen Comp Endocrinol. 213:16–23. and characterization of sex genes, phylogenetic analyses, Forconi M, et al. 2013. Characterization of sex determination and sex interpreted the data and wrote the paper. M.F. and M.B. differentiation genes in Latimeria. PLoS One 8(4):e56006. Garcia-Ortiz JE, et al. 2009. Foxl2 functions in sex determination and his- did the dN/dS analyses and contributed to interpretation of togenesis throughout mouse ovary development. BMC Dev Biol. 9:36. the data. A.P. and M.G. subjected in expression analyses. Goldman N, Yang Z. 1994. A codon-based model of nucleotide substitu- M.C.A. did the histological analyses and contributed to the tion for protein-coding DNA sequences. Mol Biol Evol. 11(5):725–736. data interpretation. M.S. conceived the study, interpreted the Han YQ, et al. 2015. Molecular cloning and tissue distribution proﬁles of data and wrote the paper. All authors made contributions to the chicken R-spondin1 gene. Genet Mol Res. 14(2):3090–3097. 1442 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Sex Genes in Lungﬁsh and Coelacanth GBE Haning RV Jr, Tantravahi U, Zhao Q, Hackett RJ, Canick JA. 1996. 5alpha- Li M, et al. 2015. A tandem duplicate of anti-Mullerian hormone with a reductase 1 and 2 expression and activity in human ovarian follicles, missense SNP on the Y chromosome is essential for male sex deter- stroma and corpus luteum as compared to neonatal foreskin. J Steroid mination in nile tilapia, Oreochromis niloticus. PLoS Genet. Biochem Mol Biol. 59(2):199–204. 11(11):e1005678. Hattori RS, Stru ¨ ssmann CA, Fernandino JI, Somoza GM. 2013. Genotypic Li X, et al. 2016. Regulation of seminiferous tubule-associated stem Leydig sex determination in teleosts: insights from the testis-determining cells in adult rat testes. Proc Natl Acad Sci U S A. 113(10):2666–2671. amhy gene. Gen Comp Endocrinol. 192:55–59. Lubieniecki KP, et al. 2015. Expression analysis of sex-determining path- Herpin A, et al. 2013. Divergent expression regulation of gonad develop- way genes during development in male and female Atlantic salmon ment genes in medaka shows incomplete conservation of the down- (Salmo salar). Physiol Genomics 47(12):581–587. stream regulatory network of vertebrate sex determination. Mol Biol Lubzens E, Young G, Bobe J, Cerda J. 2010. Oogenesis in teleosts: how eggs are formed. Gen Comp Endocrinol. 165(3):367–189. Evol. 30(10):2328–2346. Herpin A, Schartl M. 2015. Plasticity of gene-regulatory networks control- Mahendroo MS, Russell DW. 1999. Male and female isoenzymes of ste- ling sex determination: of masters, slaves, usual suspects, newcomers, roid 5alpha-reductase. Rev Reprod. 4(3):179–183. and usurpators. EMBO Rep. 16(10):1260–1274. Marshall-Graves JA. 2008. Weird animal genomes and the evolution of Hildebrand M, Goslow GE Jr. 2001. Analysis of vertebrate structure. New vertebrate sex and sex chromosomes. Annu Rev Genet. York: John Wiley and Sons. 42(1):565–586. Howe K, et al. 2013. The zebraﬁsh reference genome sequence and its Matson CK, et al. 2011. DMRT1 prevents female reprogramming in the relationship to the human genome. Nature 496(7446):498–503. postnatal mammalian testis. Nature 476(7358):101–104. Hu Q, Guo W, Gao Y, Tang R, Li D. 2015. Molecular cloning and charac- Miura T, Miura C, Konda Y, Yamauchi K. 2002. Spermatogenesis-prevent- terization of amh and dax1 genes and their expression during sex ing substance in Japanese eel. Development 129(11):2689–2697. inversion in rice-ﬁeld eel Monopterus albus. Sci Rep. 5:16667. Morinaga C, et al. 2007. The hotei mutation of medaka in the anti- Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP. 2001. Bayesian infer- Mullerian hormone receptor causes the dysregulation of germ cell ence of phylogeny and its impact on evolutionary biology. Science and sexual development. Proc Natl Acad Sci U S A. 294(5550):2310–2314. 104(23):9691–9696. Irisarri I, Meyer A. 2016. The identiﬁcation of the closest living relative(s) of Myosho T, et al. 2012. Tracing the emergence of a novel sex-determining tetrapods: phylogenomic lessons for resolving short ancient interno- gene in medaka, Oryzias luzonensis. Genetics 191(1):163–170. des. Syst Biol. 65(6):1057–1075. Nishimura T, et al. 2015. Sex determination. Foxl3 is a germ cell-intrinsic Itoh N, Konishi M. 2007. The zebraﬁsh FGF family. Zebraﬁsh factor involved in sperm-egg fate decision in medaka. Science 4(3):179–186. 349(6245):328–331. Jakob S, Lovell-Badge R. 2011. Sex determination and the control of Sox9 Pallavicini A, et al. 2013. Analysis of the transcriptome of the Indonesian expression in mammals. FEBS J. 278(7):1002–1009. coelacanth Latimeria menadoensis. BMC Genomics 14:538. Jameson SA, Lin YT, Capel B. 2012. Testis development requires the re- Pannetier M, et al. 2006. FOXL2 activates P450 aromatase gene transcrip- pression of Wnt4 by Fgf signaling. Dev Biol. 370(1):24–32. tion: towards a better characterization of the early steps of mamma- Jansson E, Mattsson A, Goldstone J, Berg C. 2016. Sex-dependent lian ovarian development. J Mol Endocrinol. 36(3):399–413. expression of anti-Mu ¨ llerian hormone (amh) and amh receptor 2 dur- Parenti LR, Grier HJ. 2004. Evolution and phylogeny of gonad morphology ing sex organ differentiation and characterization of the Mu¨llerian in bony ﬁshes. Integr Comp Biol. 44(5):333–348. duct development in Xenopus tropicalis. Gen Comp Endocrinol. Parma P, et al. 2006. R-spondin1 is essential in sex determination, skin 229:132–144. differentiation and malignancy. Nat Genet. 38(11):1304–1309. Joss JM, Edwards A, Kime DE. 1996. In vitro biosynthesis of androgens in Pasquier J, et al. 2016. Gene evolution and gene expression after whole the Australian lungﬁsh, Neoceratodus forsteri. Gen Comp Endocrinol. genome duplication in ﬁsh: the PhyloFish database. BMC Genomics 101(3):256–263. 17:368. Jovelin R, et al. 2010. Evolution of developmental regulation in the verte- Pfennig F, Standke A, Gutzeit HO. 2015. The role of Amh signaling in brate FgfD subfamily. J Exp Zool B Mol Dev Evol. 314(1):33–56. teleost ﬁsh–Multiple functions not restricted to the gonads. Gen Kaneko H, et al. 2015. Gonadal soma-derived factor (gsdf), a TGF-beta Comp Endocrinol. 223:87–107. superfamily gene, induces testis differentiation in the teleost ﬁsh Piprek RP, Pecio A, Laskowska-Kaszub K, Kubiak JZ, Szymura JM. 2013. Oreochromis niloticus. Mol Cell Endocrinol. 415:87–99., Sexual dimorphism of AMH, DMRT1 and RSPO1 localization in the Kasahara M, et al. 2007. The medaka draft genome and insights into developing gonads of six anuran species. Int J Dev Biol. 57(11– vertebrate genome evolution. Nature 447(7145):714–719. 12):891–895. Kikuchi K, Hamaguchi S. 2013. Novel sex-determining genes in ﬁsh and Piprek RP, Pecio A, Laskowska-Kaszub K, Kloc M, et al. 2013. Retinoic acid sex chromosome evolution. Dev Dyn. 242(4):339–353. homeostasis regulates meiotic entry in developing anuran gonads and Kinoshita K, Kusunoki M, Miyai K. 2006. Analysis of the three dimensional in Bidder’s organ through Raldh2 and Cyp26b1 proteins. Mech Dev. structure of the CXGXC motif in the CMGCC and CAGYC regions of 130(11–12):613–627. alpha-and beta-subunits of human chorionic gonadotropin: impor- Rodr ıguez-Mar ı A, et al. 2013. Retinoic acid metabolic genes, meiosis, and tance of glycine residue (G) in the motif. Endocr J. 53(1):51–58. gonadal sex differentiation in zebraﬁsh. PLoS One 8(9):e73951. Klu ¨ ver N, et al. 2007. Differential expression of anti-mu ¨ llerian hormone Rogers N, et al. 2013. Expression of the murine transcription factor SOX3 (amh) and anti-mu ¨ llerian hormone receptor type II (amhrII) in the tel- during embryonic and adult neurogenesis. Gene Expr Patterns eost medaka. Dev Dyn. 236(1):271–281. 13(7):240–248. Leerberg DM, Sano K, Draper BW. 2017. Fibroblast growth factor signal- Schartl M, et al. 2018. Sox5 is involved in germ-cell regulation and sex ing is required for early somatic gonad development in zebraﬁsh. PLoS determination in medaka following co-option of nested transposable Genet. 13(9):e1006993. elements. BMC Biol. 16(1):16. Lephart ED, Doody KJ, McPhaul MJ, Simpson ER. 1992. Inverse relationship Schiffer L, Anderko S, Hannemann F, Eiden-Plach A, Bernhardt R. 2015. between ovarian aromatase cytochrome P450 and 5 alpha-reductase The CYP11B subfamily. J Steroid Biochem Mol Biol. 151:38–51. enzyme activities and mRNA levels during the estrous cycle in the rat. J Schulz RW, et al. 2010. Spermatogenesis in ﬁsh. Gen Comp Endocrinol. Steroid Biochem Mol Biol. 42(5):439–447. 165(3):390–411. Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 1443 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018 Biscotti et al. GBE Skaar KS, et al. 2011. Proteolytically activated, recombinant anti-mullerian teleosts, Nile tilapia (Oreochromis niloticus). Gen Comp Endocrinol. hormone inhibits androgen secretion, proliferation, and differentiation 230-231:177–185. of spermatogonia in adult zebraﬁsh testis organ cultures. Yang QE, Kim D, Kaucher A, Oatley MJ, Oatley JM. 2013. CXCL12-CXCR4 Endocrinology 152(9):3527–3540. signaling is required for the maintenance of mouse spermatogonial Takasawa K, et al. 2014. FOXL2 transcriptionally represses Sf1 expression stem cells. J Cell Sci. 126(Pt 4):1009–1020. by antagonizing WT1 during ovarian development in mice. FASEB J. Yao HH, et al. 2004. Follistatin operates downstream of Wnt4 in mam- 28(5):2020–2028. malian ovary organogenesis. Dev Dyn. 230(2):210–215. Takehana Y, et al. 2014. Co-option of Sox3 as the male-determining factor Yu M, et al. 2013. RALDH2, the enzyme for retinoic acid synthesis, medi- on the Y chromosome in the ﬁsh Oryzias dancena. Nat Commun. ates meiotic initiation in germ cells of the female embryonic chickens. 5:4157. Amino Acids 44(2):405–412. Talavera G, Castresana J. 2007. Improvement of phylogenies after remov- Zhang T, Murphy MW, Gearhart MD, Bardwell VJ, Zarkower D. 2014. The ing divergent and ambiguously aligned blocks from protein sequence mammalian Doublesex homolog DMRT6 coordinates the transition alignments. System Biol. 56(4):564–577. between mitotic and meiotic developmental programs during sperma- Uribe MC, Grier HJ, Mej ıa-Roa V. 2014. Comparative testicular structure togenesis. Development 141(19):3662–3671. and spermatogenesis in bony ﬁshes. Spermatogenesis 4(3):e983400. Zhang X, Wang H, et al. 2014. Isolation of doublesex- and mab-3-related Wagner GP, Kin K, Lynch VJ. 2012. Measurement of mRNA abundance transcription factor 6 and its involvement in spermatogenesis in tilapia. using RNA-seq data: rPKM measure is inconsistent among samples. Biol Reprod. 91(6):136. Theory Biosci. 131(4):281–285. Zhang X, et al. 2016. Autosomal gsdf acts as a male sex initiator in the ﬁsh Webster KA, et al. 2017. Dmrt1 is necessary for male sexual development medaka. Sci Rep. 6:19738. in zebraﬁsh. Dev Biol. 422(1):33–46. Zhang Z, et al. 2006. KaKs calculator: calculating Ka and Ks through model Weiss J, et al. 2003. Sox3 is required for gonadal function, but not sex selection and model averaging. Genomics Proteomics Bioinformatics determination, in males and females. Mol Cell Biol. 4(4):259–263. 23(22):8084–8091. Zhou L, et al. 2016. Rspo1-activated signalling molecules are sufﬁcient to Westernstro ¨ er B, et al. 2014. Proﬁling of Cxcl12 receptors, Cxcr4 and induce ovarian differentiation in XY medaka (Oryzias latipes). Sci Rep. Cxcr7 in murine testis development and a spermatogenic depletion 6:19543. model indicates a role for Cxcr7 in controlling Cxcl12 activity. PLoS Zhu Y, Wang C, Chen X, Guan G. 2016. Identiﬁcation of gonadal soma- One 9(12):e112598. derived factor involvement in Monopterus albus (protogynous rice Wrobel KH. 2003. The genus Acipenser as a model for vertebrate urogen- ﬁeld eel) sex change. Mol Biol Rep. 43(7):629–637. ital development: the muellerian duct. Anat Embryol. 206(4):255–271. Wu L, Yang P, Luo F, Wang D, Zhou L. 2016. R-spondin1 signaling path- way is required for both the ovarian and testicular development in a Associate editor: Maria Costantini 1444 Genome Biol. Evol. 10(6):1430–1444 doi:10.1093/gbe/evy101 Advance Access publication May 29, 2018 Downloaded from https://academic.oup.com/gbe/article-abstract/10/6/1430/5020732 by Ed 'DeepDyve' Gillespie user on 20 June 2018
Genome Biology and Evolution – Oxford University Press
Published: May 29, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
All the latest content is available, no embargo periods.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud