GSE160494
Transcriptomic and epigenetic preservation of genetic sex identity in estrogen-feminized male chicken embryonic gonads
Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within one year. To identify the feminization-resistant “memory” of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (e.g., DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex. Overall design: Transcriptomes and DNA methylomes of day-19 chicken embryonic gonads with or without exposure to 17alpha-ethynylestradiol in ovo were examined by deep sequencing methods RNA-seq and MBD-seq, respectively.