- Research article
- Open Access
Identification of candidate gonadal sex differentiation genes in the chicken embryo using RNA-seq
BMC Genomics volume 16, Article number: 704 (2015)
- The Erratum to this article has been published in BMC Genomics 2016 17:169
Despite some advances in recent years, the genetic control of gonadal sex differentiation during embryogenesis is still not completely understood. To identify new candidate genes involved in ovary and testis development, RNA-seq was used to define the transcriptome of embryonic chicken gonads at the onset of sexual differentiation (day 6.0/stage 29).
RNA-seq revealed more than 1000 genes that were transcribed in a sex-biased manner at this early stage of gonadal differentiation. Comparison with undifferentiated gonads revealed that sex biased expression was derived primarily from autosomal rather than sex-linked genes. Gene ontology and pathway analysis indicated that many of these genes encoded proteins involved in extracellular matrix function and cytoskeletal remodelling, as well as tubulogenesis. Several of these genes are novel candidate regulators of gonadal sex differentiation, based on sex-biased expression profiles that are altered following experimental sex reversal. We further characterised three female-biased (ovarian) genes; calpain-5 (CAPN5), G-protein coupled receptor 56 (GPR56), and FGFR3 (fibroblast growth factor receptor 3). Protein expression of these candidates in the developing ovaries suggests that they play an important role in this tissue.
This study provides insight into the earliest steps of vertebrate gonad sex differentiation, and identifies novel candidate genes for ovarian and testicular development.
Gonadal sex differentiation in vertebrate embryos involves sexually dimorphic gene expression, leading to ovarian or testicular development. While a number of key gonadal differentiation genes have been identified and validated in a variety of organisms, large gaps in our understanding still exist [1, 2]. This likely underlies the significant number of human Disorders of Sex Development (DSD’s) that remain unexplained at the genetic level . Identification of novel regulators of embryonic ovarian and testis differentiation is therefore required, and will greatly improve our understanding of normal and aberrant gonad sexual differentiation. The chicken embryo provides a unique model system for studying gonad development because development occurs in ovo (in the egg), allowing gonadal development to be directly manipulated . In chickens and other birds, as in mammals, sex is determined by sex chromosome inheritance. However, birds have a ZZ male: ZW female sex chromosome system, which is the opposite of the XY male: XX female system found in mammals. Genes located on the sex chromosomes drive differentiation of the “bipotential” gonads into ovaries (in ZW females) or testes (in ZZ males). Several key genes involved in mammalian gonadal sex differentiation are conserved in chicken, including male up-regulated DMRT1, SOX9 and AMH, and female up-regulated Aromatase, FOXL2 and WNT4/RSPO1 [2–12]. However, the Y-chromosome linked SRY gene, which initiates testis development in mammals, is absent in birds. The best candidate master sex switch in birds is DMRT1, which is present on the Z sex chromosome and is expressed more highly in male gonads compared to females prior to and during sexual differentiation. Knockdown and over-expression of DMRT1 causes feminisation and masculinisation of the gonads respectively [3, 7]. Although a W-linked female (ovarian) determinant may yet exist, our extensive analysis has so far not produced a convincing candidate gene [13, 14].
In chicken embryos, the gonads form on the medioventral surface of the mesonephric kidneys at embryonic day 3 (E3 - Hamburger and Hamilton stage 18) . At this stage, they are undifferentiated or “bipotential”. At E6 (HH29), gonads begin morphological differentiation into testes in ZZ embryos or unilateral ovary in ZW embryos. Sex-specific gonadal morphology emerges at this time, and a small number of genes showing sexually dimorphic expression have been identified, notably the ovarian determinants FOXL2 and aromatase in females [11, 16, 17], and DMRT1, HEMOGEN and SOX9 in males [1, 4, 8, 10, 18–21]. However, in chicken, the exact functional relationships among these genes are unclear. Furthermore, human DSD’s imply the existence of other as yet unidentified gonad-determining genes in vertebrates generally.
Previous screens in the mouse embryo have employed various methods to isolate novel genes regulating gonadal sex differentiation, such as high throughput whole mount in situ hybridization (WISH) , differential display [23, 24], cDNA subtractive hybridization [25–28] cDNA microarrays [19, 29, 30] and RNA-seq on sorted germ cells . However, technical limitations often mean that these approaches assay only a small fraction of the transcriptome, and some studies have sampled stages when the gonads are larger but when differentiation has already occurred. Nevertheless, several of these large-scale studies have successfully identified hundreds of genes with sexually dimorphic expression in embryonic mouse gonads. A microarray study on sorted gonadal somatic cells (using the Sf1-GFP reporter) from 10.5 and 11.5dpc found numerous sex-biased genes specific to the early differentiating Sertoli and granulosa cells . Another study used microarrays on gonadal somatic cells to identify 2306 genes expressed in a sex-specific manner prior to, during and after gonadal sex differentiation (E10.5 – E13.5) . An independent microarray screen compared sorted supporting cells (Sry or Sox9 – eGFP), interstitial or stromal cells (Mafb-eGFP) and germ cells (Oct4-eGFP) from XY and XX mouse gonads at E11.5, E12.5 and E13.5 . These studies as well as the others mentioned above identified novel candidates such as Aard, Cst9, Mmd2, Cbln4, some of which have been further characterised in gonad development [24, 35–38]. Whilst these studies represent extensive transcriptome profiling of embryonic gonads, deep RNA sequencing (RNA-seq) is considerably more sensitive . In addition, RNA-seq reveals expression for all mRNAs, not only genes covered by microarray, allowing an unbiased analysis.
In the current study, RNA-seq is used to assess sex-biased gene expression at the onset gonadal sex differentiation in embryonic chicken gonads. Our rationale is based on the assumption that genes important for ovary versus testis development will show sex-biased expression. We previously performed RNA sequencing on E4.5 gonads, at the bi-potential stage prior to sexual differentiation . By comparing gene expression from this early stage to E6 (the onset of morphological differentiation), we find novel genes and pathways that are activated sex-specifically at the time of sexual differentiation; more than 1000 genes are transcribed in a sex-biased manner at this stage, and the majority of these become biased between E4.5 and E6. This is primarily due to an increased expression of autosomal, rather than sex-linked, genes between these time points. Gene ontology and pathway analysis indicate that these genes are involved in extracellular signalling as well as driving morphological changes such as ECM, cytoskeletal remodelling and tubulogenesis. We also find significant enrichment of genes reported in the embryonic mouse gonadal screens in the genes identified in our chicken screen, highlighting conservation of molecular mechanisms underlying gonad differentiation. Several novel candidate gonadal sex differentiation genes are identified and validated in vivo. Three female-biased (ovarian) genes are further characterised; calpain-5 (CAPN5), G-protein coupled receptor 56 (GPR56), and FGFR3 (fibroblast growth factor receptor 3). These genes may play an important role during ovarian development. This study provides insight into the earliest steps of gonad sex differentiation, and identifies novel candidate genes for ovarian and testis development.
Significant sex-biased gene expression in early differentiating gonads (E6.0/stage 29)
To assess gene expression in the early differentiating testes and ovaries, paired gonads from E6 (HH29) were dissected and mRNA isolated as previously described . For each sex, 2 pools of 18 gonad pairs (from 18 individuals) were collected. Messenger RNA was then processed and sequenced (Methods), resulting in an average of 46 000 000 read pairs per sample. We have previously performed the same analysis on E4.5 paired gonads . Expression for each gene and differential expression (DE) analysis were then calculated (Methods). Gene expression is calculated in Fragments Per Kilobase of Exon Per Million Fragments Mapped (FPKM). Significance is calculated using False Discovery Rate (FDR) (see Methods).
RNA-seq analysis detected expression of 17108 genes. As the gonads begin to differentiate at E6, a total of 1003 genes showed significant DE between the sexes (FDR <0.001). Of these, 347 were female-biased and 656 were male-biased (Fig. 1a). In contrast, at the earlier time point of E4.5, when the gonads are morphologically undifferentiated, only 267 genes were significantly DE (193 male biased genes and 74 female biased genes; FDR <0.001) (Fig. 1a, b). Thus, between E4.5 and E6 as the gonads begin to differentiate, a four-fold increase in sexually dimorphic transcript number was observed (from 1.56 to 5.86 % of all expressed genes). Substantial overlap was observed between the genes that were sexually DE at E4.5 and those at E6; 165 genes in males and 62 in females (85.5 and 83.8 % of the E4.5 DE genes were also DE at E6.5 in males and females respectively) (Fig. 1a). Indeed, of the 267 DE genes at E4.5, only 40 (15 %) lost this bias at E6 (Fig. 1a) and the majority of these (all but 5) reside on the sex chromosomes [W, Z or suspected W chromosome (un-random)] . This may reflect changes in dosage compensation of these sex chromosome genes or changes in the expression of W linked genes [14, 40–42].
Hierarchical clustering was performed on the 400 most significantly DE genes (those with the lowest P-values) from E4.5 and E6 gonads. The clustering of genes on the basis of expression level is represented in a heat map (Fig. 1b). W-linked genes clustered together and showed high differential expression in both female samples at both time-points (Fig. 1b), consistent with our previous observations . Similarly, Z-linked genes clustered together. Importantly, autosomal genes that were DE were primarily present in E6 and not E4.5 gonads, with female DE genes clustering together for each sample and male DE genes clustering. The overall number of DE sex-linked genes did not greatly change between the time-points (E4.5 = 224, E6 = 304). This indicates that the increase in biased genes expression between E4.5 and E6 was largely due to transcriptional changes in autosomal genes, which increased from 43 DE genes at E4.5 (16 %) to 699 (70 %) at E6 (Fig. 1c). Taken together, these data show that between E4.5 and E6, there was a significant increase in the number of autosomal genes showing sexually dimorphic expression. These genes are novel candidates for contributing to gonad sexual differentiation.
Known sex-differentiation genes show dimorphic expression
The expression profiles of known or putative regulators of gonadal sexual differentiation were assessed at E4.5 and E6 (Table 1). Of the known sexual regulators based on previous studies in chicken or mammals, only one autosomal gene was significantly DE at E4.5 in females; FOXL2 (FDR <0.001). At E6, FOXL2 continued to be differentially expressed, and in addition CYP19A1/Aromatase, HSD17B1 and NROB1 (DAX-1) transcripts became female enriched (Table 1). Four genes showed significantly male-biased expression at both time-points; HSD17B4, DMRT1, HEMOGEN, and AMH (FDR <0.001). Inhibin A and B trended towards sexual dimorphism (FDR 0.0027, 0.0012 respectively) at E4.5, and by E6 both were significantly male-enriched (Table 1). SOX9 showed male-biased expression at E6 only, confirming previous studies suggesting that it lies downstream of Z-linked genes DMRT1 and HEMOGEN [3, 20]. These data imply that major components of the sex determination pathways are activated between E4.5 and E6 (developmental stages HH 25/26 to HH 29/30). Several other sex-development genes show gonadal expression but no sexual dimorphism (Table 1). Indeed, many known regulators of bipotential gonad formation such as LHX9, EMX2, GATA4 and FOG2 were very highly expressed in both female and male E4.5 gonads, while expression decreased significantly in E6 gonads (Table 1). This confirms that the E4.5 gonads are still undifferentiated, whereas at E6 the differentiation process has commenced, reflected by down-regulation of genes required for the formation of the undifferentiated gonad and up-regulation of sex-differentiation genes.
While no RNA-seq analysis has yet been reported for whole embryonic mouse gonads, large-scale transcriptional screens have been carried out . To assess conservation of sexually dimorphic genes in other organisms, chicken DE genes at E6 (P-value <0.05) were compared with those from three separate large-scale mouse microarray screens [32, 33, 43]. The first of these studies assayed somatic cells from XX and XY embryonic mouse gonads at E10.5, E11.5, E12.5 and E13.5, with E11.5 – 12.5 being equivalent to chicken E6 (stage 29/30) . The second study used isolated somatic cells from XX and XY gonads at 10.5 and 11.5dpc and identified genes showing male up-regulation or down-regulation specific to 11.5dpc, as well as female up-regulation at 11.5dpc . In the third study, isolated supporting cells, stromal/interstitials cells and germ cells were profiled from XX and XY mice . All three of these studies used different methods of assessing differential expression, but all used a P-value of 0.05 to assess DE genes at each time point. To compare our data to these studies, we used a P-value of 0.05 to find E6 DE genes in the chicken. This results in a list of 2123 genes. This gene list was then compared with those identified as sex-biased at 11.5 and 12.5dpc from Nef et al.  and Jameson et al. , and with those identified in the other study . Fisher’s exact tests were also carried out to test for statistically significant enrichment of each data set in our set. When we compared each data set to our DE gene list we found significant enrichment (P-value <0.001) for each study (Table 2). Our RNA-seq study was carried out on whole gonads, and will hence encompass expression from various populations of cells include supporting cells and germ cells. By comparing to each cell lineage from the Jameson et al. study, we could also see that DE genes associated with supporting cell (443), and interstitial cells (210) are represented in our data set at a significant level (Table 2 P-value <0.001). Germ cells were enriched with a P-value of 0.00123 (163 genes). Overall, we found 76 genes in common between the three mouse studies alones (specifically in the supporting cell populations) (Table 2). Of these, 20 genes were also DE in our study (Table 2). Thus significant overlap between differentiating mouse and chickens gonads was observed, suggesting that similar genetic pathways are engaged between mammals and avians, validating the use of chicken embryonic gonads to screen for novel sex determination candidates.
20 genes were common between all studies of somatic cells and our study (Table 2). These genes were Col14a1, Cyp11a1, Ogt, Slc20a1, Lect1, Nedd9, Prss12, Fst, Hpse, Gpm6b, Tesc, Atp1a1, Dapk1, Adm, Scarb2, Mef2c, Tyro3, Gdnf, Bace2 and Nac1. Some of these genes have been implicated in sex development previously. CYP11A1 mutations cause congenital adrenal insufficiency and with partial or complete 46,XY sex reversal in humans. GDNF or glial cell line derived neurotrophic factor exerts a proliferative effect on developing mouse Sertoli cells . Finally, Tescalin (Tesc) has been previously described as an early testis gene [45, 46], although no gonadal phenotypes have been described in a null mutant mouse . We also analysed the overlap between these studies and our dataset using a P-value cut off of 0.001 to determine DE genes. This dataset has 1003 genes, half of what we have using the less stringent P-value of 0.05. Consistently, we found reduced numbers of overlapping genes, but enrichment was still significant between all groups (Additional file 1: Table S3).
Transcriptional changes in the differentiating gonads reflect morphological differentiation
To elucidate novel pathways engaged during sexual differentiation of ovaries and testes, pathway analysis was performed on the dataset. To this end, genes were filtered that were only sexually dimorphic at E6 (FDR <0.001), but not in the undifferentiated gonads (E4.5). Genes previously implicated in sex determination or differentiation were removed, to reveal novel activated pathways. Gene ontology (GO) analysis was carried out on this list of 774 genes (Additional file 2: Table S1) using DAVID [48, 49]. We assessed three ontological categories: Biological process (BP), Cellular component (CC) and Molecular function (MF) (Table 3). The top 3 BP categories in DE genes were tube development, cell proliferation and respiratory tube development. The CCs most represented were extracellular and cell surface. The MFs were lipid binding, growth factor activity and ATPase activity (Table 3).
To assess the probable functions of the genes in this data set, we used DAVID functional analysis of KEGG pathways . The top five KEGG pathways that were significantly engaged (based on calculated P-value <0.05) were focal adhesion, dorsal-ventral axis formation, regulation of actin cytoskeleton, ECM-receptor interaction, and cytokine-cytokine receptor interaction (Table 3). Numerous DE genes were implicated in each of these pathways (i.e. Additional file 3: Figure S1A, B). The high number of genes involved in extracellular function and in the engagement of pathways involved in actin cytoskeletal remodelling and ECM function is consistent with the significant morphological changes that are taking place as the gonads differentiate into ovaries or testis. Indeed, the embryonic ovary and testis differ significantly in their morphology; the ovary develops a thickened cortex to house germ cells, while the medulla becomes fragmented, and the testis organizes into basement membrane bound seminiferous cords that enclose germ cells and supporting Sertoli cells. These morphological events begin at E6. Vasculature and seminiferous cord development are important in the testis , which was reflected here in the number of genes implicated in tube formation and development. Panther pathway analysis also revealed that the TGF-β pathway was significantly represented (P-value = 0.05), including BMP2, 3, 4 and 7 and BMPR2 as well as JUN and Myostatin. Indeed, BMP signalling has been implicated in ovarian and follicle development [50–53], although the exact role of these BMPs at this early stage needs to be elucidated.
Novel candidates for testis and ovarian development
Known sex determination genes were effectively identified from the RNA-seq data set, many of which had sexually dimorphic expression only at E6. Based on this analysis, the sequencing data was interrogated to identify novel genes involved in chicken gonadal sex differentiation. For potential testis-determining genes, the data was filtered for genes showing sexually dimorphic expression at E6 (FDR < 0.001) which were male biased (Log2 fold change/Fem:male ratio of less than 1) and with a male gonadal expression of more than 10 FPKM. FPKM of 10 was chosen as a cut-off based on the expression of known gonad genes and on previous studies in our lab that suggest in vivo methods such as WISH do not detect the expression of genes with less than this value (data not shown). These testis genes were ranked according to P-value and the top 45 genes were tabulated (Additional file 4: Table S2). This list includes previously mentioned known sex-regulators such as DMRT1, HEMOGEN, AMH and SOX9. From this list, four candidate genes were chosen for further study, based on gene ontology and expression. Specifically we chose genes with a variety of expression levels/fold changes, and which are thought to be involved in transcriptional regulation or morphogenesis/cell differentiation. These were laminin1 (LAMA1), BMP receptor 2 (BMPR2), zinc-finger 385B (ZNF385b) and a novel, uncharacterised gene on the Z chromosome that we designated NZP (Table 4). To identify genes involved in ovarian development, genes were filtered for significant sexual dimorphism at E6 (FDR <0.001), female bias (Log2 fold change > than 1) and with a female expression greater than FPKM 10. All W-linked genes were removed as these have been analyzed previously , and the remaining genes were ranked according to P-value (Additional file 4: Table S2), and the top 45 taken. From this list, three novel candidate genes for ovarian development were chosen based on their expression profiles and on their gene ontology. These were calpain-5 (CAPN5), G-protein coupled receptor 56 (GPR56) and fibroblast growth factor receptor 3 (FGFR3) (Table 4).
Whole-mount in situ hybridisation (WISH) of chicken embryonic urogenital systems (UGS) at E6 was used to confirm mRNA expression of the seven candidate genes (Figs. 2 and 3). By this method, genes could be assessed for their degree of gonad specificity, and their expression could be localised to a gonadal cell type. WISH confirmed male-biased expression in the E6 gonads for all four of the testis candidate genes; LAMA1, BMPR2, ZNF385B, and NZP (Fig. 2a-d versus i-l). Over-staining and sectioning showed that all four of these genes were expressed in the seminiferous cords of the developing testis (Fig. 2e-h). BMPR2, ZNF385b and NZP had only low level, diffuse expression in the ovary (Fig. 2i-o), while LAMA1 had expression in the juxta-cortical medulla of the ovary (Fig. 2l, p). None of these genes were notably expressed at other sites within the urogenital system.
All three candidate ovarian genes showed strong female-biased expression at E6.0 (stage 29/30) (Fig. 3a-c versus g-i) in WISH. CAPN5 and GPR56 were strongly expressed in female gonads but only weakly expressed in males (Fig. 3a, b and g, h). GPR56 and FGFR3 had no expression in other sites within the urogenital system, while CAPN5 showed expression in the adrenal glands of both males and females (Fig. 3a, g arrows). Sectioning revealed that CAPN5 was expressed in a subset of medullary cells close to the cortex – the juxta-cortical medulla (Fig. 3d). GPR56 and FGFR3 were expressed primarily in the cortex of the ovary (Fig. 3e, f). Taken together these results indicate that the gonad RNA-seq data was robustly validated by in situ hybridisation.
Novel sex-differentiation candidates show altered expression following experimental sex-reversal of gonads
RNA-seq and WISH analysis revealed the expression of three female- and four male-biased gonadal genes. To test the involvement of these genes in the gonadal differentiation pathways, experimental sex reversal was carried out using the aromatase inhibitor, fadrozole. Fadrozole inhibits the activity of the aromatase enzyme, blocking oestrogen production, and resulting in female to male sex reversal of the embryonic gonads [54, 55]. Fadrozole was injected into fertile eggs at E3 (stage 19) and gonads were harvested from embryos at E9 (stage 33). WISH was used to compare gene expression to that of stage-matched controls.
The male-biased genes, BMPR2 and ZNF385b (Fig. 4a-c and d-f) were up-regulated in sex-reversed females (Fig. 4c, f) compared to control females (Fig. 4b, e), suggesting that the expression of these is a consequence of male sex differentiation. While NZP also showed up-regulation in the treated female (Fig. 4i), expression of this gene is low at this time-point in both sexes (Fig. 4g, h). Similarly, LAMA1 expression was weaker in the male at this time point, and only a weak activation was observed in the female treated with fadrozole (Fig. 4j-l). Further functional analysis of these three candidate testis genes is currently underway. Antibodies have been raised against ZNF385b and knockdown is being conducted using our established virally delivered RNAi approach . The female-biased genes CAPN5 (Fig. 4m-o), GPR56 (Fig. 4p-r) and FGFR3 (Fig. 4s-u) all showed down-regulation in fadrozole treated female gonads (Fig. 4o, r, u) compared to the untreated female. In GPR56, only the left gonad of females showed robust expression at this point (Fig. 4q). This was lost in fadrozole treated gonads (Fig. 4r). These data show that these genes are expressed in female-specific cell types as a consequence of the early steps of ovarian differentiation. Their expression is subsequently lost when ovaries are sex reversed. These data indicate that the RNA-seq screen was used to successfully identify candidate genes for ovarian or testis development. As our interest lies in identifying novel genes involved in ovary differentiation, we further focused on CAPN5, GPR56 and FGFR3.
Candidate ovarian genes are translated during gonad development
Only a small number of genes have so far been shown to play a deterministic role in ovarian differentiation and development . Screens such as the one described here provide additional candidate genes, which will fill gaps in our understanding of ovarian morphogenesis. Therefore we focused on the female-biased candidates FGFR3, CAPN5 and GPR56. FGFR3 is a receptor for fibroblast growth factors (FGFs) which have been shown to play a central role in gonad differentiation and germ cell development [56–58]. A role for FGFR3 has not been described in the ovary during embryogenesis. Calpain-5 (CAPN5) encodes an intracellular calcium-dependent cysteine protease, and has similarity to the Caenorhabditis elegans sex determination gene tra-3. CAPN5 has not been implicated in gonad development in vertebrates. GPR56 is an adhesion-associated G-protein coupled receptor that mediates intracellular signal transduction. GPR56 null mice have impaired testis formation and infertility . No phenotype has been reported in female mice.
To study these proteins, a commercial antibody was used for FGFR3 and in-house antibodies were raised against CAPN5 and GPR56 proteins. These new antibodies were tested for specificity (Additional file 5: Figure S2). To this end, open reading frames of chicken CAPN5 and GPR56 were cloned and expressed them in a chicken fibroblastic cell line (Additional file 5: Figure S2). Two polyclonal antibodies were tested for each protein, at various concentrations (data not shown). For GPR56, while low background staining was observed in untransfected DF1 cells (Additional file 5: Figure S2A), strong expression was observed when transfection with the GPR56 plasmid was carried out (Additional file 5: Figure S2B). This staining appeared to be cytoplasmic, excluded from the nucleus (Additional file 5: Figure S2Bi, Bii). CAPN5 specific staining is observed only in the transfected cells and not in the control cells (Additional file 5: Figure S2C, D), and it appears punctate throughout the cytoplasm (Di, Dii). These antibodies were then used to assay expression in embryonic gonads. CAPN5 was detected in female gonads at E6 in a subset of cells in the medulla. By E8.5 (HH stage 32) this expression had increased and was evident in a subset of juxta-cortical medullary cells (Fig. 5b). CAPN5 appeared to be cytoplasmic in these cells (Fig. 5c). This protein expression pattern was consistent with the mRNA expression pattern. Some low level punctate staining was observed in the males from E8.5 (Fig. 5a). GPR56 was also expressed from E6.5 (HH 30) in the female, with only background staining in the male. Expression of GPR56 was highest in the cortex of the ovary where it appeared to be cytoplasmic and at the plasma membrane (Fig. 5e, f). (See also Additional file 5: Figure S2E, F for additional images.) The FGFR3 antibody detected protein in the germ cells in both sexes (Fig. 5g, h), where it is cytoplasmic or at the plasma membrane (Fig. 5Gi). Staining was stronger in the ovary (Fig. 5h), especially in the cortex (Fig. 5i), which may reflect its expression in somatic cortical cells as well as germ cells.
The female-biased expression of these proteins, especially CAPN5 and GPR56, indicates that they might play a role in ovarian development from the time of sexual differentiation, and further validates RNA-seq as a useful tool for the identification of novel candidate gonad differentiation genes. Further studies including functional analysis are currently underway in our laboratory to assess the precise role of these genes.
Sexual differentiation of the gonads during embryonic development involves the sex-biased expression and action of several well-characterised regulators of ovarian or testicular development. However, only a small number of genes have been identified and many cases of unexplained human DSDs point to the existence of uncharacterised genetic regulators of embryonic gonadal development. RNA-seq was employed here to identify novel candidate sex genes, using embryonic chicken gonads as a model system. This study is the first to report RNA-seq in entire gonads at the time of sexual differentiation in a vertebrate species. We have combined previous analysis of female and male gonads at E4.5, prior to sexual differentiation, with new data from E6 gonads. By comparing these two time-points, a detailed view of the gonad transcriptome during sexual differentiation has been revealed.
While the expression of several known or putative gonad sex determinants has previously been assessed in embryonic chicken gonads [10, 16, 19, 60–63], this has often been based on assessment of suspected genes and pathways, employing whole mount in situ, qRT-PCR, or protein expression. These screens have been successful in identifying candidates, some of which are also found here (e.g., CLDN11, INHA, InhibinA and B and BMPR2). . In addition, some of these studies have identified transcriptional differences between the left and right ovary, the latter of which regresses in chicken [64, 65]. We did not separate left and right ovary in this study. Nevertheless, these approaches are biased towards known pathways inferred from other vertebrates. They lack both the sensitivity and detail of RNA-seq. Indeed, RNA-seq can reveal early changes in gene expression undetectable by other methods.
Firstly, the RNA-seq data has provided information on the expression of known or putative sex genes. These results are interesting when considering the inferred molecular hierarchy of genes and their signalling in the gonads, both in chicken and other species. For example, in females FOXL2 shows sexual dimorphism from E4.5, preceding that of aromatase, which does not become significantly DE until E6. This is consistent with previous studies that have suggested that FOXL2 acts upstream of other female genes, at least at the transcriptional level [11, 16, 66]. Indeed, in other species ranging from fishes to mammals, FOXL2 is a putative or known activator of aromatase gene expression [67, 68]. Given that Aromatase expression and oestrogen synthesis is essential for ovary formation in birds, FOXL2 likely occupies a key position in the female pathway. In addition, a recent study has suggested that the main role of FOXL2 is as a repressor of male genes , and this may explain the importance of early FOXL2 expression where it would drive the gonadal primordium towards an ovarian fate by repressing male genes and activating female genes.
In the avian embryo, and indeed other vertebrates, the activator of FOXL2 has not been clearly defined, and while a W-chromosome linked gene may fulfil this role in the chicken, no convincing candidate exists. Alternatively, FOXL2 may be activated by a factor common to both sexes that is antagonized by elevated expression of Z-linked DMRT1. DMRT1 can antagonise FOXL2 expression, as evidenced in DMRT1 knockdown gonads, which exhibit elevated FOXL2 transcription and feminization . Both Z-linked DMRT1 and HEMOGEN show significant male-bias at E4.5, as does AMH, preceding SOX9, which becomes male up-regulated at E6. This indicates that DMRT1/Hemogen activate SOX9 transcription, which is supported by recent over-expression studies in ovo [3, 10, 20].
As the gonads enter the phase of sexual differentiation, they up-regulate sexually dimorphic gene expression, most of which is autosomal. Given that some known sex determining genes are dimorphic prior to this activation (e.g., DMRT1 and FOXL2), these activated genes at E6 (stage 29/30) are likely downstream regulators of ovarian versus testis development. Gene ontology and KEGG pathway analysis suggests that many of these genes are involved in the restructuring of the gonad into the morphologically distinct testis or ovary. There was an over-representation of genes involved in ECM and cytoskeleton formation and regulation in the datasets. Indeed, a study that used network approaches to analyse sex-biased gene expression in chickens suggested that genes of the same sex bias tended to be more connected to each other than expected, and that when genes were grouped into sex-biased modules these showed functional homogeneity .
Our RNA-seq was carried out on whole chicken gonads and therefore encompasses expression from both supporting cell types (presumptive granulosa and Sertoli) as well as from interstitial cells and from germ cells. To assess the relative importance of each cell type and the potential overlap with mouse gonadal expression profiles we compared our screen and three high throughput microarray screens carried out on embryonic mouse gonads/cells at a similar stage of differentiation [32–34]. We found significant enrichment (P-value <0.001) between the genes identified in each study and our own. In addition, these analyses indicate that we have found DE genes from not just supporting cell types, but also from the germ cell and stromal/interstitial cell pools, albeit to a lesser extent. Interestingly, we found only 76 genes were common between the three mouse studies (specifically in the supporting cells only). This probably reflects the subtle differences between each of these screens such as the transgenic mice used, the sorting methods, the microarray platforms and the analyses. Interestingly though, of these 76 genes 20 were also DE in our data set. The presence of several already characterised gonadal genes in this list means that these represent excellent candidates for further investigation.
The RNA-seq analysis also provided novel candidate genes for sexual differentiation in female and male chicken embryos. By comparing expression in females and males at E6 numerous significantly sexually dimorphic genes not previously implicated in gonad development were identified. To validate these lists, several male and female biased genes of interest were chosen, based on their gene ontology and expression. These were validated in vivo using WISH to assess their expression in the UGS and in sex reversed gonads. For all genes assayed, expression profiles in WISH mimicked those in the RNA-seq data, which we have also found previously with W genes . For male-biased genes, expression was present in the testis cords, and was up-regulated when female gonads were experimentally sex-reversed using the aromatase inhibitor fadrozole. The female genes assayed were either present in the medulla (CAPN5) or the cortex (FGFR3 and GPR56) of the developing ovary. In the presence of fadrozole, these female-biased genes had reduced expression, which was similar to the level observed in ZZ males. These results suggest that RNA-seq has successfully identified candidate genes for gonadal sexual differentiation in the avian system.
Only a small number of genes implicated in ovarian development have been described. Many of these have derived from analysis of human or animal disorders of sex development, i.e. WNT4 , R-spondin1 [72, 73], FOXL2  reviewed in . There are still gaps in our understanding of both the initial triggers and the downstream modulators of ovarian development. This study has identified numerous female-biased genes, which may play a role in ovarian differentiation. Indeed, we focused our studies on three of these genes, using specific antibodies against the chicken proteins. Protein expression was similar to that revealed by the RNA-seq and WISH analyses. Closer investigation into expression provided insight into the potential role of these genes.
One candidate gene was FGFR3, an FGF receptor. FGF9 plays a key role in testis differentiation, and in the regulation of germ cell development [56–58]. Several FGF receptors exist in vertebrate species. Loss of FGFR2 in mice testis leads to partial sex reversal , but intriguingly both FGFR2 and FGFR3 show significant female-biased expression in our E6 data. FGFR3 has been shown to interact with FGF1 and FGF9  and is thought to be important in the post-natal testis . A role for FGFR3 has not yet been described in the ovary during embryogenesis. Interestingly, FGF9 itself was very weakly expressed in E6 gonads, and was not sex-biased, and so the role of this FGF may not be conserved outside mammals. (Although a role could be played by a related FGF). Furthermore, we did not identify any FGF ligand showing male-biased expression at E6. In contrast, FGF16 showed robust female-biased expression. This FGF is also female-biased in Nile tilapia gonads , and may play a role in early oocyte development. It is interesting to speculate that FGF16 signalling though FGFR2/3 in the female may contribute to ovary differentiation or female germ cell maturation in chickens.
CAPN5 was expressed in the adrenal glands of both sexes, and the gonads of the female only, where it was found in a subset of cells present in the juxta-cortical medulla. The 3βHSD enzyme, which is involved in the synthesis of sex steroids testosterone and oestrogen is also expressed in a similar subset of cells as well as the adrenal glands (Lambeth et al. Submitted). This indicates that different cell lineages exist in the medulla and it appears that CAPN5 is expressed in a subset of these. These data suggest that this novel candidate may be involved in the differentiation or function of steroid producing cells in the ovary as well as in the adrenal glands of both sexes.
In contrast, our third ovarian candidate, GPR56, was expressed in cortex of the differentiating ovary as it proliferates during ovarian development. WNT4 and RSPO1 are expressed in the cortex of the developing ovary, and their activated transcription factor β-catenin appears to be stabilized there . Recently, another G-protein coupled receptor, Lgr5 has been shown to be expressed in cortical granulosa precursor cells, and plays a role in germ cell differentiation and survival . Indeed, LGR5 binds R-spondin proteins and enhances RSPO/WNT/β-catenin signalling [80–83]. It will be interesting to determine if GPR56 also interacts with this pathway in the cortex of the chicken ovary.
These genes are all interesting candidates for ovarian development and functional studies including overexpression and knockdown are now underway in our laboratory to establish their exact roles. Most recently, we have developed the method of in ovo electroporation to deliver genes or antisense RNAs directly into developing gonads [3, 84]. This methodology is being used to manipulate the novel candidate genes identified in the current study.
Several caveats to this type of analysis that focuses on sexually dimorphic genes is that not all key genes involved in ovarian versus testicular development will be expressed in a sex-biased fashion. Some genes may be expressed in both sexes but only have a key role on one sex. In mouse embryos, Gadd45g and Map3k4 are expressed in the early gonad equally in both sexes, yet mutations in these genes lead only to sex reversal in XY males [85–87] because they act on the SRY gene (present on the Y chromosome) through the p38 MAPK pathway. Typically these genes have been identified in human DSD cases or mouse mutants. In addition, this analysis was based on genes correctly annotated in the Galgal4 chicken genome, thus missing genes that are unannotated or lie in regions that have typically been difficult to sequence. De novo transcript assembly of the mRNAs expressed at E6 may reveal additional novel candidates for sex determination.
This study represents the first RNA-seq analysis of embryonic gonads at the time of sexual differentiation. These data have provided a detailed view of the transcriptional changes that occur when sexual differentiation begins. Over 1000 genes were significantly sexually DE at this stage, most of which are autosomal genes. Many of these likely contribute to the morphological changes taking place, such as ECM remodelling and tubulogenesis. These data provide several new candidate ovary or testis differentiation genes, some of which we have validated in vivo. In particular, we show that CAPN5 and GPR56 are excellent candidates for ovarian development. Functional studies on these novel candidates are now underway in our laboratory, with a view to clarifying their roles in gonadal development and differentiation.
Embryo incubation and tissue collection
All experimental work adhered to the Australian code of practice for the care and use of animals for scientific purposes (7th edition, 2004) and to in-house animal handling guidelines (Murdoch Childrens Research Institute and CSIRO Australian Animal Health Laboratory). Ethics was from MCRI - Chicken Ethics: Institute Animal Ethics number A627. IBC (Institute Biosafety Committee): 103ext-2007, and 194–2013 EXEMPT. All experiments were carried out on a single line of Specific Pathogen Free embryos (SPAFAS) from the White Leghorn strain of chick (Lohman-LSL) and incubated at 37.8 °C with rocking in a humidified incubator. For RNA-seq, SPF eggs were incubated until stage 29 (embryonic day 6) or E4.5 see . Paired gonads were removed and stored at −80 °C. Embryos were sexed by PCR as described previously , briefly, a small piece of limb tissue was digested in PCR compatible buffer containing proteinase K (200 μg/ml at 50 °C for at least 30 min), followed by rapid PCR amplification of the sex-linked, female-specific Xho1 sequence. Amplification of 18S ribosomal RNA genomic sequence was used as the internal control in a duplex reaction.
RNA extraction and sequencing
Tissues were pooled according to sex. Eighteen paired gonads (i.e. left and right gonads from 18 individuals) were pooled for each replicate (two male replicates and two female replicates). Pooling has been documented as an appropriate way to prepare samples for expression analysis . Total RNA was extracted using the RNeasy micro kit (QIAGEN). On column DNAsing was performed to remove contaminating genomic DNA. The resulting 4 RNA replicates were poly A-selected, reverse transcribed, fragmented, bar-coded and sequenced using the Illumina HiSeq2000 at Australian Genome Research Facility (AGRF) in Melbourne. We sequenced 100 base pair, paired ends reads.
Differential expression (DE) analysis on the raw reads was carried out by testing the female read counts against male read counts using edgeR , using count data from E6 and earlier time-points from . P-values were adjusted for multiple testing using the Benjamini and Hochberg false discovery rate (FDR) method . Sequenced read pairs (from both E6 and E4.5 time-points) were mapped to the chicken genome (galGal4) using the TopHat 2.0.6 software . We used Ensembl version 73, which has 17108 genes. The Ensembl gene annotation was provided to TopHat to aid alignment, with other settings default. Read pairs were summarized to Ensembl gene-level counts using featureCounts . Default edgeR settings were used including normalisation using TMM. Homology of the read-pairs with Ensembl genes was then calculated, and expression for each gene calculated in Fragments Per Kilobase Of Exon Per Million Fragments Mapped (FPKM).
Genes were sorted based on their minimum P-value (irrespective of time-point) and the top 400 DE genes were hierarchical clustered. Clustering was performed on the proportion of counts in each sample with respect to the sum across samples, after normalising for library size. The clustering and heatmap were created using the hclust and heatmap.2 functions in R respectively (using default settings - agglomeration clustering using complete linkage). Gene ontology and KEGG pathway analysis were carried out using DAVID according to the default actions.
Gene lists (converted to mouse universal gene symbol) were compiled for each significantly sexually dimorphic gene (generally P-value <0.05, FC of 1.5 or 2) described in each study using DAVID software (from either Ensembl ID, microarray ID or Unigene ID). Gene lists were then compared using R (version 3.2.0), and one-sided Fisher exact test was used to obtain a P-value for enrichment.
Embryos for in vivo studies
Fertile eggs (SPF) were obtained from CSIRO, Werribee and incubated under humid conditions at 37.8 °C. Embryos were harvested at various days throughout development and staged according to the morphological criteria of Hamburger and Hamilton . The UGS including the mesonephric kidneys, gonads and Müllerian duct were dissected out in PBS and fixed for either WISH or immunofluorescence.
Whole-mount in situ hybridisation
Plasmids containing fragments of mRNA for each gene were ordered from the BBSRC ChickEST Database. These were verified by PCR and sequencing. The clones were as follows; ZNF385b - ChEST470h3, BMPR2 - ChEST10191l3, NZP – ChEST203m10, Laminin - ChEST791a8 CAPN5 – ChEST992n12, FGFR3 – ChEST779m22. For the GPR56 and the ZNF385b the ORFs were also used to generate probes in pGEM T-easy. Probes were made using the either T7, Sp6 or T3 RNA polymerase and labelled with digoxygenin (DIG). In situ hybridization was performed as described previously [9, 94]. Briefly, UGS were dissected from embryos and fixed overnight at 4 °C in 4 % paraformaldehyde. At least three embryos were used for each sex and time point or condition. They were then processed for in situ hybridization with the DIG-labelled antisense riboprobe. For negative controls, a DIG-labelled sense riboprobe was generated for each gene. Alkaline phosphatase conjugated Anti-DIG antibodies were used, and the chromagen was BCIP/NBT. Following WISH, specimens were photographed again under bright field and were then over-stained before being imbedded in OCT for cryosectioning (sections were cut at 10, 14 and 18 μm). Staining was analysed by photography under bright field and photos taken of sections using 5x and 20x lenses. Both sense and antisense probes were initially tested for all genes. Staining was only observed with antisense probes, and not in the sense controls (data not shown).
Urogenital tissues were dissected from chicken embryos and briefly fixed for 15 min in 4 % paraformaldehyde/PBS at room temperature. Tissues were then cryoprotected by immersion in 30 % sucrose/PBS (overnight at 4 °C), infiltrated with OCT embedding compound. Ten micron frozen sections on slides were treated as previously described . Secondary antibodies were AlexFluor 488 donkey anti-rabbit IgG, 1:1000, and AlexaFluor 594 donkey anti-mouse IgG, at 1:1500 from Invitrogen). Sections were mounted in Fluorosave (Calbiochem). Images were taken on a Zeiss Axiovision M1.
Antibodies were raised against the chicken GPR56 and CAPN5 proteins at GenScript. Two affinity purified polyclonal rabbit antibodies were raised for each, against the peptides VSEPIDLTEGDYTTC or CHLQDRGNRRSNDLP for CAPN5 and CAFASPKEENREVQG or CWREDGTASSGNWDS, for GPR56. The FGFR3 antibody was from Bioss (bs-0165R) and used at 1:100.
To test the GPR56 antibody in vitro, the GPR56 ORF was cloned from chicken gonad cDNA using the primers GPR56_Fow GAACCCAGGCTGAGAGCCAG, GPR56_Rev GCAAACACCCTCTGCCAGGC. Capn5_Fow GAAAGGATGTTTTCCTCA, Capn5_Rev TACACTTGCAGGTGGTAGA. GPR56 was cloned into the Tol2 plasmid under the CMV promoter. Chicken fibroblasts were transfected according to the Lipofectamine 2000 protocol (Life Technologies) with either empty or GPR56 Tol2 vector. Anti-GPR56 antibodies 1 and 2 were tested at various concentrations (1 :500, 1: 2000, 1:10,000) and assessed for specific staining using immunofluorescence as described above. CAPN5 ORF was cloned into the RCAS.BP.A avian retroviral vector. This was transfected into cells, which were passaged at least once to allow the spread of virus. These cells and control/uninfected cells were then fixed and stained with the two antibodies as above. Anti-P27 was used to check viral infection (not shown).
Availability of supporting data
RNA-seq data for both time points has been uploaded to the DDBJ database as part of the Chicken BioProject (ref number SRP014719).
Disorders of sex development
Hamburger and Hamilton stage
Green fluorescent protein
Fragments Per Kilobase of Exon Per Million Fragments Mapped
False discovery rate
Kyoto Encyclopedia of Genes and Genomes
Specific Pathogen Free embryos
Eggers S, Sinclair A. Mammalian sex determination—insights from humans and mice. Chromosome Res. 2012;20:215–38.
Cutting A, Chue J, Smith CA. Just how conserved is vertebrate sex determination? Dev Dyn 2013;242(4):380-387.
Lambeth LS, Raymond CS, Roeszler KN, Kuroiwa A, Nakata T, Zarkower D, et al. Over-expression of DMRT1 induces the male pathway in embryonic chicken gonads. Dev Biol 2014;389(2):160-172.
Smith CA, Sinclair AH. Sex determination: insights from the chicken. Bioessays. 2004;26:120–32.
Ayers KL, Sinclair AH, Smith CA. The Molecular Genetics of Ovarian Differentiation in the Avian Model. Sex Dev. 2012;7:80–94.
Sinclair A, Smith C. Females battle to suppress their inner male. Cell. 2009;139:1051–3.
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–71.
Smith CA, Smith MJ, Sinclair AH. Gene expression during gonadogenesis in the chicken embryo. Gene. 1999;234:395–402.
Smith CA, McClive PJ, Western PS, Reed KJ, Sinclair AH. Conservation of a sex-determining gene. Nature. 1999;402:601–2.
Vaillant S, Magre S, Dorizzi M, Pieau C, Richard-Mercier N. Expression of AMH, SF1, and SOX9 in gonads of genetic female chickens during sex reversal induced by an aromatase inhibitor. Dev Dyn. 2001;222:228–37.
Chue J, Smith CA. Sex determination and sexual differentiation in the avian model. FEBS J. 2011;278:1027–34.
Lambeth LS, Cummins DM, Doran TJ, Sinclair AH, Smith CA. Overexpression of Aromatase Alone is Sufficient for Ovarian Development in Genetically Male Chicken Embryos. PLoS ONE. 2013;8:e68362.
Smith CA, Roeszler KN, Sinclair AH. Genetic evidence against a role for W-linked histidine triad nucleotide binding protein (HINTW) in avian sex determination. Int J Dev Biol. 2009;53:59–67.
Ayers KL, Davidson NM, Demiyah D, Roeszler KN, Grützner F, Sinclair AH, et al. RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken embryos and allows comprehensive annotation of W-chromosome genes. Genome Biol. 2013;14:R26.
Hamburger V, Hamilton HL. A series of normal stages in the development of the chick embryo. J Morphol. 1951;88:49–92.
Govoroun MS, Pannetier M, Pailhoux E, Cocquet J, Brillard J-P, Couty I, et al. Isolation of chicken homolog of the FOXL2 gene and comparison of its expression patterns with those of aromatase during ovarian development. Dev Dyn. 2004;231:859–70.
Smith CA, Andrews JE, Sinclair AH. Gonadal sex differentiation in chicken embryos: expression of estrogen receptor and aromatase genes. J Steroid Biochem Mol Biol. 1997;60:295–302.
Shan Z, Nanda I, Wang Y, Schmid M, Vortkamp A, Haaf T. Sex-specific expression of an evolutionarily conserved male regulatory gene, DMRT1, in birds. Cytogenet Cell Genet. 2000;89:252–7.
Smith CA, Katz M, Sinclair AH. DMRT1 is upregulated in the gonads during female-to-male sex reversal in ZW chicken embryos. Biol Reprod. 2003;68:560–70.
Nakata T, Ishiguro M, Aduma N, Izumi H, Kuroiwa A. Chicken hemogen homolog is involved in the chicken-specific sex-determining mechanism. Proc Natl Acad Sci. 2013;110:3417–22.
Ayers KL, Smith CA, Lambeth LS. The molecular genetics of avian sex determination and its manipulation. Genesis. 2013;51(5):1-12.
Wertz K, Herrmann BG. Large-scale screen for genes involved in gonad development. Mech Dev. 2000;98:51–70.
Nordqvist K, Töhönen V. An mRNA differential display strategy for cloning genes expressed during mouse gonad development. Int J Dev Biol. 1997;41:627–38.
Töhönen V, Osterlund C, Nordqvist K. Testatin: a cystatin-related gene expressed during early testis development. Proc Natl Acad Sci U S A. 1998;95:14208–13.
McClive PJ, Sinclair AH. Rapid DNA extraction and PCR-sexing of mouse embryos. Mol Reprod Dev. 2001;60:225–6.
Menke DB, Page DC. Sexually dimorphic gene expression in the developing mouse gonad. Gene Expr Patterns. 2002;2:359–67.
Bowles J, Bullejos M, Koopman P. A subtractive gene expression screen suggests a role for vanin-1 in testis development in mice. Genesis. 2000;27:124–35.
Boyer A, Lussier JG, Sinclair AH, McClive PJ, Silversides DW. Pre-sertoli specific gene expression profiling reveals differential expression of Ppt1 and Brd3 genes within the mouse genital ridge at the time of sex determination. Biol Reprod. 2004;71:820–7.
Grimmond S, Van Hateren N, Siggers P, Arkell R, Larder R, Soares MB, et al. Sexually dimorphic expression of protease nexin-1 and vanin-1 in the developing mouse gonad prior to overt differentiation suggests a role in mammalian sexual development. Hum Mol Genet. 2000;9:1553–60.
Small CL, Shima JE, Uzumcu M, Skinner MK, Griswold MD. Profiling gene expression during the differentiation and development of the murine embryonic gonad. Biol Reprod. 2005;72:492–501.
Lesch BJ, Dokshin GA, Young RA, McCarrey JR, Page DC. A set of genes critical to development is epigenetically poised in mouse germ cells from fetal stages through completion of meiosis. Proc Natl Acad Sci. 2013;110:16061–6.
Beverdam A, Koopman P. Expression profiling of purified mouse gonadal somatic cells during the critical time window of sex determination reveals novel candidate genes for human sexual dysgenesis syndromes. Hum Mol Genet. 2006;15:417–31.
Nef S, Schaad O, Stallings NR, Cederroth CR, Pitetti JL, Schaer G, et al. Gene expression during sex determination reveals a robust female genetic program at the onset of ovarian development. Dev Biol. 2005;287:361–77.
Jameson SA, Natarajan A, Cool J, Defalco T, Maatouk DM, Mork L, et al. Temporal transcriptional profiling of somatic and germ cells reveals biased lineage priming of sexual fate in the fetal mouse gonad. PLoS Genet. 2012;8:e1002575.
Svingen T, Beverdam A, Verma P, Wilhelm D, Koopman P. Aard is specifically up-regulated in Sertoli cells during mouse testis differentiation. Int J Dev Biol. 2007;51:255–8.
Bradford ST, Hiramatsu R, Maddugoda MP, Bernard P, Chaboissier MC, SINCLAIR A, et al. The cerebellin 4 precursor gene is a direct target of SRY and SOX9 in mice. Biol Reprod. 2009;80:1178–88.
Hasegawa K, Chuma S, Tada T, Sakurai T, Tamura M, Suemori H, et al. Testatin transgenic and knockout mice exhibit normal sex-differentiation. Biochem Biophys Res Commun. 2006;341:369–75.
Töhönen V, Frygelius J, Mohammadieh M, Kvist U, Pelliniemi LJ, O’Brien K, et al. Normal sexual development and fertility in testatin knockout mice. Mol Cell Biol. 2005;25:4892–902.
Oshlack A, Robinson MD, Young MD. From RNA-seq reads to differential expression results. Genome Biol. 2010;11:220.
Zhang SO, Mathur S, Hattem G, Tassy O, Pourquié O. Sex-dimorphic gene expression and ineffective dosage compensation of Z-linked genes in gastrulating chicken embryos. BMC Genomics. 2010;11:13.
Itoh Y, Melamed E, Yang X, Kampf K, Wang S, Yehya N, et al. Dosage compensation is less effective in birds than in mammals. J Biol. 2007;6:2.
Moghadam HK, Pointer MA, Wright AE, Berlin S, Mank JE. W chromosome expression responds to female-specific selection. Proceedings of the National Academy of Sciences 2012.
Coveney D, Ross AJ, Slone JD, Capel B. A microarray analysis of the XX Wnt4 mutant gonad targeted at the identification of genes involved in testis vascular differentiation. Gene Expr Patterns. 2007;7:82–92.
Wu Z, Templeman JL, Smith RA, Mackay S. Effects of glial cell line-derived neurotrophic factor on isolated developing mouse Sertoli cells in vitro. J Anat. 2005;206:175–84.
Perera EM, Martin H, Seeherunvong T, Kos L, Hughes IA, Hawkins JR, et al. Tescalcin, a novel gene encoding a putative EF-hand Ca(2+)-binding protein, Col9a3, and renin are expressed in the mouse testis during the early stages of gonadal differentiation. Endocrinology. 2001;142:455–63.
Bao Y, Hudson QJ, Perera EM, Akan L, Tobet SA, Smith CA, et al. Expression and evolutionary conservation of the tescalcin gene during development. Gene Expr Patterns. 2009;9:273–81.
Ukarapong S, Bao Y, Perera EM, Berkovitz GD. Megakaryocyte development is normal in mice with targeted disruption of Tescalcin. Exp Cell Res. 2012;318:662–9.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35(Web Server issue):W169–75.
Onagbesan OM, Bruggeman V, Van As P, Tona K, Williams J, Decuypere E. BMPs and BMPRs in chicken ovary and effects of BMP-4 and −7 on granulosa cell proliferation and progesterone production in vitro. Am J Physiol Endocrinol Metab. 2003;285:E973–83.
Kashimada K, Pelosi E, Chen H, Schlessinger D, Wilhelm D, Koopman P. FOXL2 and BMP2 act cooperatively to regulate follistatin gene expression during ovarian development. Endocrinology. 2011;152:272–80.
Hoshino A, Koide M, Ono T, Yasugi S. Sex-specific and left-right asymmetric expression pattern of Bmp7 in the gonad of normal and sex-reversed chicken embryos. Dev Growth Differ. 2005;47:65–74.
Guioli S, Lovell-Badge R. PITX2 controls asymmetric gonadal development in both sexes of the chick and can rescue the degeneration of the right ovary. Development. 2007;134:4199–208.
Elbrecht A, Smith RG. Aromatase enzyme activity and sex determination in chickens. Science. 1992;255:467–70.
Vaillant S, Dorizzi M, Pieau C, Richard-Mercier N. Sex reversal and aromatase in chicken. J Exp Zool. 2001;290:727–40.
Kim Y, Kobayashi A, Sekido R, DiNapoli L, Brennan J, Chaboissier M-C, et al. Fgf9 and Wnt4 act as antagonistic signals to regulate mammalian sex determination. PLoS Biol. 2006;4:e187.
Bowles J, Feng C-W, Spiller C, Davidson T-L, Jackson A, Koopman P. FGF9 suppresses meiosis and promotes male germ cell fate in mice. Dev Cell. 2010;19:440–9.
DiNapoli L, Batchvarov J, Capel B. FGF9 promotes survival of germ cells in the fetal testis. Development. 2006;133:1519–27.
Chen G, Yang L, Begum S, Xu L. GPR56 is essential for testis development and male fertility in mice. Dev Dyn. 2010;239:3358–67.
Nomura O, Nakabayashi O, Nishimori K, Yasue H, Mizuno S. Expression of five steroidogenic genes including aromatase gene at early developmental stages of chicken male and female embryos. J Steroid Biochem Mol Biol. 1999;71:103–9.
Yamamoto I, Tsukada A, Saito N, Shimada K. Profiles of mRNA expression of genes related to sex differentiation of the gonads in the chicken embryo. Poult Sci. 2003;82:1462–7.
Oreal E, Pieau C, Mattei MG, Josso N, Picard JY, Carre-Eusebe D, et al. Early expression of AMH in chicken embryonic gonads precedes testicular SOX9 expression. Dev Dyn. 1998;212:522–32.
Smith CA, Shoemaker CM, Roeszler KN, Queen J, Crews D, Sinclair AH. Cloning and expression of R-Spondin1 in different vertebrates suggests a conserved role in ovarian development. BMC Dev Biol. 2008;8:72.
Carré G-A, Couty I, Hennequet-Antier C, Govoroun MS. Gene expression profiling reveals new potential players of gonad differentiation in the chicken embryo. PLoS ONE. 2011;6:e23959.
Intarapat S, Stern CD. Sexually dimorphic and sex-independent left-right asymmetries in chicken embryonic gonads. PLoS ONE. 2013;8:e69893.
Hudson QJ, Smith CA, Sinclair AH. Aromatase inhibition reduces expression of FOXL2 in the embryonic chicken ovary. Dev Dyn. 2005;233:1052–5.
Pannetier M, Fabre S, Batista F, Kocer A, Renault L, Jolivet G, et al. FOXL2 activates P450 aromatase gene transcription: towards a better characterization of the early steps of mammalian ovarian development. J Mol Endocrinol. 2006;36:399–413.
Wang DS, Kobayashi T, Zhou LY, Paul-Prasanth B, Ijiri S, Sakai F, et al. Foxl2 up-regulates aromatase gene transcription in a female-specific manner by binding to the promoter as well as interacting with ad4 binding protein/steroidogenic factor 1. Mol Endocrinol. 2007;21:712–25.
Elzaiat M, Jouneau L, Thépot D, Klopp C, Allais-Bonnet A, Cabau C, et al. High-Throughput Sequencing Analyses of XX Genital Ridges Lacking FOXL2 Reveal DMRT1 Up-Regulation Before SOX9 Expression During the Sex-Reversal Process in Goats. Biol Reprod. 2014;91:153.
Frings O, Mank JE, Alexeyenko A, Sonnhammer ELL. Network analysis of functional genomics data: application to avian sex-biased gene expression. ScientificWorldJournal. 2012;2012:130491.
Vainio S, Heikkilä M, Kispert A, Chin N, McMahon AP. Female development in mammals is regulated by Wnt-4 signalling. Nature. 1999;397:405–9.
Tomaselli S, Megiorni F, De Bernardo C, Felici A, Marrocco G, Maggiulli G, et al. Syndromic true hermaphroditism due to an R-spondin1 (RSPO1) homozygous mutation. Hum Mutat. 2007;29:220–6.
Parma P, Radi O, Vidal V, Chaboissier M-C, Dellambra E, Valentini S, et al. R-spondin1 is essential in sex determination, skin differentiation and malignancy. Nat Genet. 2006;38:1304–9.
Boulanger L, Pannetier M, Gall L, Allais-Bonnet A, Elzaiat M, Le Bourhis D, et al. FOXL2 is a female sex-determining gene in the goat. Curr Biol. 2014;24:404–8.
Bagheri-Fam S, Sim H, Bernard P, Jayakody I, Taketo MM, Scherer G, et al. Loss of Fgfr2 leads to partial XY sex reversal. Dev Biol. 2008;314:71–83.
Chellaiah A, Yuan W, Chellaiah M, Ornitz DM. Mapping ligand binding domains in chimeric fibroblast growth factor receptor molecules. Multiple regions determine ligand binding specificity. J Biol Chem. 1999;274:34785–94.
Goriely A, Hansen RMS, Taylor IB, Olesen IA, Jacobsen GK, McGowan SJ, et al. Activating mutations in FGFR3 and HRAS reveal a shared genetic origin for congenital disorders and testicular tumors. Nat Genet. 2009;41:1247–52.
Sun Y-L, Zeng S, Ye K, Yang C, Li M-H, Huang B-F, et al. Involvement of FGF9/16/20 subfamily in female germ cell development of the Nile tilapia, Oreochromis niloticus. Fish Physiol Biochem. 2012;38:1427–39.
Rastetter RH, Bernard P, Palmer JS, Chassot A-A, Chen H, Western PS, et al. Marker genes identify three somatic cell types in the fetal mouse ovary. Dev Biol. 2014;394:242–52.
Carmon KS, Gong X, Lin Q, Thomas A, Liu Q. R-spondins function as ligands of the orphan receptors LGR4 and LGR5 to regulate Wnt/beta-catenin signaling. Proc Natl Acad Sci. 2011;108:11452–7.
de Lau W, Barker N, Low TY, Koo B-K, Li VSW, Teunissen H, et al. Lgr5 homologues associate with Wnt receptors and mediate R-spondin signalling. Nature 2011;476:293-7.
Ruffner H, Sprunger J, Charlat O, Leighton-Davies J, Grosshans B, Salathe A, et al. R-Spondin potentiates Wnt/β-catenin signaling through orphan receptors LGR4 and LGR5. PLoS ONE. 2012;7:e40976.
Chen P-H, Chen X, Lin Z, Fang D, He X. The structural basis of R-spondin recognition by LGR5 and RNF43. Genes Dev. 2013;27:1345–50.
Ayers KL, Cutting AD, Roeszler KN, Sinclair AH, Smith CA. DMRT1 is required for Müllerian duct formation in the chicken embryo. Dev Biol. 2015.
Gierl MS, Gruhn WH, von Seggern A, Maltry N, Niehrs C. GADD45G functions in male sex determination by promoting p38 signaling and Sry expression. Dev Cell. 2012;23:1032–42.
Warr N, Carré G-A, Siggers P, Faleato JV, Brixey R, Pope M, et al. Gadd45γ and Map3k4 interactions regulate mouse testis determination via p38 MAPK-mediated control of Sry expression. Dev Cell. 2012;23:1020–31.
Bogani D, Siggers P, Brixey R, Warr N, Beddow S, Edwards J, et al. Loss of mitogen-activated protein kinase kinase kinase 4 (MAP3K4) reveals a requirement for MAPK signalling in mouse sex determination. PLoS Biol. 2009;7:e1000196.
Clinton M, Haines L, Belloir B, McBride D. Sexing chick embryos: a rapid and simple protocol. Br Poult Sci. 2001;42:134–8.
Kendziorski C, Irizarry RA, Chen K-S, Haag JD, Gould MN. On the utility of pooling biological samples in microarray experiments. Proc Natl Acad Sci U S A. 2005;102:4252–7.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26:139–40.
Benjamini YAHY. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.
Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.
Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Andrews JE, Smith CA, Sinclair AH. Sites of estrogen receptor and aromatase expression in the chicken embryo. Gen Comp Endocrinol. 1997;108:182–90.
We thank the Australian Genomic Research Institute (AGRF) in Melbourne. AO is supported by NHMRC Career Development Fellowship (1051481). CS is funded through an Australian Research Council (ARC) Future Fellowship. KA and AS are supported by an NHMRC Program grant (1074258). MCRI is supported by the Victorian Government’s Operational Infrastructure Support Program.
The authors declare they have no competing interests.
KA participated in study design, dissected and extracted RNA, carried out data interpretation, performed in vivo validations and drafted the manuscript. ND participated in study design and carried out bioinformatic analysis on RNA-seq. LL helped with dissections and in vivo validations. AS participated in study design and provided critical analysis of data. AO participated in study design and in bioinformatic analysis. CS conceived of the study and participated in its design, dissected and extracted RNA and aided with in vivo validation. All authors read and approved the final manuscript.
Mouse screen comparison using chicken DE genes (p-value <0.001). Results from one-sided Fisher exact test using chicken DE genes with p-value 0.001 as a cutoff. (XLSX 34 kb)
Genes used for DAVID analysis of gene ontology and KEGG pathways. Total genes were filtered for those showing only DE expression (P-value <0.001) at E6 and not E4.5, to isolate those genes that become dimorphically expressed as the gonad begin to differentiate. We removed well-characterized sex genes such as FOXL2 and SOX9 to allow identification of additionally activated pathways. This resulted in a list of 774 genes, which were subject to DAVID analysis. In the table expression (shown as FPKM) is shown for female and males at E4.5 and E6. P-values for the differences between females and males are shown for each time point (BH_PValue 4.5 or 6). P-values are adjusted using Benjamini & Hochberg method. (XLSX 193 kb)
Example of top KEGG pathways showing sex-biased expression in E6 gonads. DAVID KEGG pathway analysis was carried out on all genes showing a significant sexually dimorphic expression (P-value <0.001) at E6, but excluding those at E4.5 and known sex genes (Additional file 2: Table S1). This revealed several additional pathways that were significantly represented. (A) Focal adhesion pathway, which was the most highly represented KEGG pathway in DAVID analysis of all genes. Red stars indicate those genes that were in our DE list. (B) Pathway for the regulation of actin cytoskeleton. Red starts indicate genes found in our list. (TIFF 11721 kb)
Top 45 female and top 45 male biased genes. For male genes we filtered for sexually dimorphic genes at E6 (p < 0.001) which were male biased (Fem:male ratio of less than 0), with a male gonadal expression of more than FPKM 10 at E6. These were ranked according to P-value. We took the top 45 genes. For female genes we filtered for sexually dimorphic genes at E6 (p < 0.001) which were female biased (Fem:male ratio of greater than 0), with a female gonadal expression of more than FPKM 10 at E6. In the table expression (shown as FPKM) is shown for female and males at E4.5 and E6. P-values for the differences between females and males are shown for each time point (BH_PValue 4.5 or 6). P-values are adjusted using Benjamini & Hochberg method. (XLSX 68 kb)
Validation of newly raised antibodies for chicken GPR56 and CAPN5. Antibodies were tested at various concentrations (not shown) in chicken DF1 fibroblasts either transfected with a control plasmid or a plasmid expression protein of choice. (A) Only low background staining is shown for GPR56 in control cells, whereas when transfected with GPR56 expressing plasmid staining was observed (B). This staining is mostly cytoplasmid (see Bi, Bii). (C) Some weak background staining was seen in control cells, compared to those over-expressing CAPN5 (D). This staining is punctate and cytoplasmic (Di, Dii), similar to what is seen in vivo. Additional examples of in vivo staining for GPR56 in male E8.5 gonads (E) and female (F). (TIFF 11721 kb)
About this article
Cite this article
Ayers, K.L., Lambeth, L.S., Davidson, N.M. et al. Identification of candidate gonadal sex differentiation genes in the chicken embryo using RNA-seq. BMC Genomics 16, 704 (2015). https://doi.org/10.1186/s12864-015-1886-5
- Sex determination
- Sexual differentiation
- Chicken embryo