RNA sequencing reveals sexually dimorphic gene expression before gonadal differentiation in chicken and allows comprehensive annotation of the W-chromosome

Overview of chicken embryo genes related to sex differentiation

Sex determination in chickens at an early embryonic stage has been a longstanding challenge in poultry production due to the unique ZZ:ZW sex chromosome system and various influencing factors. This review has summarized the genes related to the sex differentiation of chicken early embryos (mainly Dmrt1, Sox9, Amh, Cyp19a1, Foxl2, Tle4z1, Jun, Hintw, Ube2i, Spin1z, Hmgcs1, Foxd1, Tox3, Ddx4, cHemgn and Serpinb11 in this article), and has found that these contributions enhance our understanding of the genetic basis of sex determination in chickens, while identifying potential gene targets for future research. This knowledge may inform and guide the development of sex screening technologies for hatching eggs and support advancements in gene-editing approaches for chicken embryos. Moreover, these insights offer hope for enhancing animal welfare and promoting conservation efforts in poultry production.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

Extraordinary variability in gene activation and repression programs during gonadal sex differentiation across vertebrates

Genes involved in gonadal sex differentiation have been traditionally thought to be fairly conserved across vertebrates, but this has been lately questioned. Here, we performed the first comparative analysis of gonadal transcriptomes across vertebrates, from fish to mammals. Our results unambiguously show an extraordinary overall variability in gene activation and repression programs without a phylogenetic pattern. During sex differentiation, genes such as dmrt1, sox9, amh, cyp19a and foxl2 were consistently either male- or female-enriched across species while many genes with the greatest expression change within each sex were not. We also found that downregulation in the opposite sex, which had only been quantified in the mouse model, was also prominent in the rest of vertebrates. Finally, we report 16 novel conserved markers (e.g., fshr and dazl) and 11 signaling pathways. We propose viewing vertebrate gonadal sex differentiation as a hierarchical network, with conserved hub genes such as sox9 and amh alongside less connected and less conserved nodes. This proposed framework implies that evolutionary pressures may impact genes based on their level of connectivity.

Sex-specific transcriptome of the chicken chorioallantoic membrane

Dimorphism between male and female embryos has been demonstrated in many animal species, including chicken species. Likewise, extraembryonic membranes such as the chorioallantoic membrane (CAM) are likely to exhibit a sex-specific profile. Analysis of the previously published RNA-seq data of the chicken CAM sampled at two incubation times, revealed 783 differentially expressed genes between the CAM of male and female embryos. The expression of some of these genes is sex-dependant only at one or other stage of development, while 415 genes are sex-dependant at both developmental stages. These genes include well-known sex-determining and sex-differentiation genes (DMRT1, HEGM, etc.), and are mainly located on sex chromosomes. This study provides evidence that gene expression of extra-embryonic membranes is differentially regulated between male and female embryos. As such, a better characterisation of associated mechanisms should facilitate the identification of new sex-specific biomarkers.

Gene Expression Profiling Reveals Potential Players of Sex Determination and Asymmetrical Development in Chicken Embryo Gonads

Despite the notable progress made in recent years, the understanding of the genetic control of gonadal sex differentiation and asymmetrical ovariogenesis in chicken during embryonic development remains incomplete. This study aimed to identify potential key genes and speculate about the mechanisms associated with ovary and testis development via an analysis of the results of PacBio and Illumina transcriptome sequencing of embryonic chicken gonads at the initiation of sexual differentiation (E4.5, E5.5, and E6.5). PacBio sequencing detected 328 and 233 significantly up-regulated transcript isoforms in females and males at E4.5, respectively. Illumina sequencing detected 95, 296 and 445 DEGs at E4.5, E5.5, and E6.5, respectively. Moreover, both sexes showed asymmetrical expression in gonads, and more DEGs were detected on the left side. There were 12 DEGs involved in cell proliferation shared between males and females in the left gonads. GO analysis suggested that coagulation pathways may be involved in the degradation of the right gonad in females and that blood oxygen transport pathways may be involved in preventing the degradation of the right gonad in males. These results provide a comprehensive expression profile of chicken embryo gonads at the initiation of sexual differentiation, which can serve as a theoretical basis for further understanding the mechanism of bird sex determination and its evolutionary process.

Preliminary Study on Expression and Function of the Chicken W Chromosome Gene MIER3 in Embryonic Gonads

The sex chromosomes of birds are designated Z and W. The male is homogamous (ZZ), and the female is heterogamous (ZW). The chicken W chromosome is a degenerate version of the Z chromosome and harbors only 28 protein-coding genes. We studied the expression pattern of the W chromosome gene MIER3 (showing differential expression during gonadogenesis) in chicken embryonic gonads and its potential role in gonadal development. The W copy of MIER3 (MIER3-W) shows a gonad-biased expression in chicken embryonic tissues which was different from its Z copy. The overall expression of MIER3-W and MIER3-Z mRNA and protein is correlated with the gonadal phenotype being higher in female gonads than in male gonads or female-to-male sex-reversed gonads. Chicken MIER3 protein is highly expressed in the nucleus, with relatively lower expression in the cytoplasm. Overexpression of MIER3-W in male gonad cells suggested its effect on the GnRH signaling pathway, cell proliferation, and cell apoptosis. MIER3 expression is associated with the gonadal phenotype. MIER3 may promote female gonadal development by regulating EGR1 and αGSU genes. These findings enrich our knowledge of chicken W chromosome genes and support a more systematic and in-depth understanding of gonadal development in chickens.

Comparative Transcriptome Analyses of Leg Muscle during Early Growth between Geese (Anser cygnoides) Breeds Differing in Body Size Characteristics

Goose is an important poultry commonly raised for meat. The early growth performance of geese significantly influences their market weight and slaughter weight, affecting the poultry industry's economic benefits. To identify the growth surge between the Shitou goose and the Wuzong goose, we collected the early growth body traits from 0 to 12 weeks. In addition, we investigated the transcriptomic changes in leg muscles at the high growth speed period to reveal the difference between the two geese breeds. We also estimated the growth curve parameters under three models, including the logistic, von Bertalanffy, and Gompertz models. The results showed that except for body length and keel length, the best-fitting model between the body weight and body size of the Shitou and Wuzong was the logistic model. The growth turning points of Shitou and Wuzong were 5.954 and 4.944 weeks, respectively, and the turning point of their body weight was 1459.01 g and 478.54 g, respectively. Growth surge occurred at 2-9 weeks in Shitou goose and at 1-7 weeks in Wuzong goose. The body size traits of the Shitou goose and Wuzong goose showed a trend of rapid growth in the early stage and slow growth in the later stage, and the Shitou goose growth was higher than the Wuzong goose. For transcriptome sequencing, a total of 87 differentially expressed genes (DEGs) were identified with a fold change ≥ 2 and a false discovery rate < 0.05. Many DEGs have a potential function for growth, such as CXCL12, SSTR4, FABP5, SLC2A1, MYLK4, and EIF4E3. KEGG pathway analysis identified that some DEGs were significantly enriched in the calcium signaling pathway, which may promote muscle growth. The gene-gene interaction network of DEGs was mainly related to the transmission of cell signals and substances, hematological system development, and functions. This study can provide theoretical guidance for the production and breeding management of the Shitou goose and Wuzong goose and help reveal the genetic mechanisms underlying diverse body sizes between two goose breeds.

Gene expression of male pathway genes sox9 and amh during early sex differentiation in a reptile departs from the classical amniote model

BACKGROUND

Sex determination is the process whereby the bipotential embryonic gonads become committed to differentiate into testes or ovaries. In genetic sex determination (GSD), the sex determining trigger is encoded by a gene on the sex chromosomes, which activates a network of downstream genes; in mammals these include SOX9, AMH and DMRT1 in the male pathway, and FOXL2 in the female pathway. Although mammalian and avian GSD systems have been well studied, few data are available for reptilian GSD systems.

RESULTS

We conducted an unbiased transcriptome-wide analysis of gonad development throughout differentiation in central bearded dragon (Pogona vitticeps) embryos with GSD. We found that sex differentiation of transcriptomic profiles occurs at a very early stage, before the gonad consolidates as a body distinct from the gonad-kidney complex. The male pathway genes dmrt1 and amh and the female pathway gene foxl2 play a key role in early sex differentiation in P. vitticeps, but the central player of the mammalian male trajectory, sox9, is not differentially expressed in P. vitticeps at the bipotential stage. The most striking difference from GSD systems of other amniotes is the high expression of the male pathway genes amh and sox9 in female gonads during development. We propose that a default male trajectory progresses if not repressed by a W-linked dominant gene that tips the balance of gene expression towards the female trajectory. Further, weighted gene expression correlation network analysis revealed novel candidates for male and female sex differentiation.

CONCLUSION

Our data reveal that interpretation of putative mechanisms of GSD in reptiles cannot solely depend on lessons drawn from mammals.

Overview of Avian Sex Reversal

Sex determination and differentiation are processes by which a bipotential gonad adopts either a testicular or ovarian cell fate, and secondary sexual characteristics adopt either male or female developmental patterns. In birds, although genetic factors control the sex determination program, sex differentiation is sensitive to hormones, which can induce sex reversal when disturbed. Although these sex-reversed birds can form phenotypes opposite to their genotypes, none can experience complete sex reversal or produce offspring under natural conditions. Promising evidence indicates that the incomplete sex reversal is associated with cell autonomous sex identity (CASI) of avian cells, which is controlled by genetic factors. However, studies cannot clearly describe the regulatory mechanism of avian CASI and sex development at present, and these factors require further exploration. In spite of this, the abundant findings of avian sex research have provided theoretical bases for the progress of gender control technologies, which are being improved through interdisciplinary co-operation and will ultimately be employed in poultry production. In this review, we provide an overview of avian sex determination and differentiation and comprehensively summarize the research progress on sex reversal in birds, especially chickens. Importantly, we describe key issues faced by applying gender control systems in poultry production and chronologically summarize the development of avian sex control methods. In conclusion, this review provides unique perspectives for avian sex studies and helps scientists develop more advanced systems for sex regulation in birds.

Avian Sex Determination: A Chicken and Egg Conundrum

BACKGROUND

Primary sex determination is the developmental process that results in the sexual differentiation of the gonads. Vertebrate sex determination is generally considered to follow the model based on the mammalian system, where a sex-specific master regulatory gene activates one of the two different gene networks that underlie testis and ovary differentiation.

SUMMARY

It is now known that, while many of the molecular components of these pathways are conserved across different vertebrates, a wide variety of different trigger factors are utilized to initiate primary sex determination. In birds, the male is the homogametic sex (ZZ), and significant differences exist between the avian system of sex determination and that of mammals. For example, DMRT1, FOXL2, and estrogen are key factors in gonadogenesis in birds, but none are essential for primary sex determination in mammals.

KEY MESSAGE

Gonadal sex determination in birds is thought to depend on a dosage-based mechanism involving expression of the Z-linked DMRT1 gene, and it may be that this "mechanism" is simply an extension of the cell autonomous sex identity associated with avian tissues, with no sex-specific trigger required.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

Blood transcriptome comparison between sexes and their function in 4-week Rhode Island red chickens

Sex is a major biological factor in the development and physiology of a sexual reproductive organism, and its role in the growing process is needed to be investigated in various species. We compare blood transcriptome between 5 males and 5 females in 4-week-old Rhode Island Red chickens and perform functional annotation of differentially expressed genes (DEGs). The results are as follows. 141 and 109 DEGs were located in autosomes and sex chromosomes, respectively. The gene ontology (GO) terms are significantly (p < 0.05) enriched, which were limb development, inner ear development, positive regulation of dendrite development, the KEGG pathway the TGF-beta signaling pathway, and melanogenesis (p < 0.05). These pathways are related to morphological maintenance and growth of the tissues. In addition, the SMAD2W and the BMP5 were involved in the TGF-beta signaling pathway, and both play an important role in maintaining tissue development. The major DEGs related to the development of neurons and synapses include the up-regulated NRN1, GDF10, SLC1A1, BMP5, NBEA, and NRXN1. Also, 7 DEGs were validated using RT-qPCR with high correlation (r ² = 0.74). In conclusion, the differential expression of blood tissue in the early growing chicken was enriched in TGF-beta signaling and related to the development of neurons and synapses including SMAD2W and BMP5. These results suggest that blood in the early growing stage is differentially affected in tissue development, nervous system, and pigmentation by sex. For future research, experimental characterization of DEGs and a holistic investigation of various tissues and growth stages will be required.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Prediction of sex-determination mechanisms in avian primordial germ cells using RNA-seq analysis

In birds, sex is determined through cell-autonomous mechanisms and various factors, such as the dosage of DMRT1. While the sex-determination mechanism in gonads is well known, the mechanism in germ cells remains unclear. In this study, we explored the gene expression profiles of male and female primordial germ cells (PGCs) during embryogenesis in chickens to predict the mechanism underlying sex determination. Male and female PGCs were isolated from blood and gonads with a purity > 96% using flow cytometry and analyzed using RNA-seq. Prior to settlement in the gonads, female circulating PGCs (cPGCs) obtained from blood displayed sex-biased expression. Gonadal PGCs (gPGCs) also exhibited sex-biased expression, and the number of female-biased genes detected was higher than that of male-biased genes. The female-biased genes in gPGCs were enriched in some metabolic processes. To reveal the mechanisms underlying the transcriptional regulation of female-biased genes in gPGCs, we performed stimulation tests. Retinoic acid stimulation of cultured gPGCs derived from male embryos resulted in the upregulation of several female-biased genes. Overall, our results suggest that sex determination in avian PGCs involves aspects of both cell-autonomous and somatic-cell regulation. Moreover, it appears that sex determination occurs earlier in females than in males.

Sex-limited chromosomes and non-reproductive traits

Sex chromosomes are typically viewed as having originated from a pair of autosomes, and differentiated as the sex-limited chromosome (e.g. Y) has degenerated by losing most genes through cessation of recombination. While often thought that degenerated sex-limited chromosomes primarily affect traits involved in sex determination and sex cell production, accumulating evidence suggests they also influence traits not sex-limited or directly involved in reproduction. Here, we provide an overview of the effects of sex-limited chromosomes on non-reproductive traits in XY, ZW or UV sex determination systems, and discuss evolutionary processes maintaining variation at sex-limited chromosomes and molecular mechanisms affecting non-reproductive traits.

RNA Sequencing of the Pituitary Gland and Association Analyses Reveal PRKG2 as a Candidate Gene for Growth and Carcass Traits in Chinese Ningdu Yellow Chickens

Growth and carcass traits are of great economic importance to the chicken industry. The candidate genes and mutations associated with growth and carcass traits can be utilized to improve chicken growth. Therefore, the identification of these genes and mutations is greatly importance. In this study, a total of 17 traits related to growth and carcass were measured in 399 Chinese Ningdu yellow chickens. RNA sequencing (RNA-seq) was performed to detect candidate genes using 12 pituitary gland samples (six per group), which exhibited extreme growth and carcass phenotypes: either a high live weight and carcass weight (H group) or a low live weight and carcass weight (L group). A differential expression analysis, utilizing RNA-seq, between the H and L groups identified 428 differentially expressed genes (DEGs), including 110 up-regulated genes and 318 down-regulated genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the identified genes showed a significant enrichment of 158 GO terms and two KEGG pathways, including response to stimulus and neuroactive ligand-receptor interaction, respectively. Furthermore, RNA-seq data, qRT–PCR, and quantitative trait transcript (QTT) analysis results suggest that the PRKG2 gene is an important candidate gene for growth and carcass traits of Chinese Ningdu yellow chickens. More specifically, association analyses of a single nucleotide polymorphism (SNP) in PRKG2 and growth and carcass traits showed that the SNP rs16400745 was significantly associated with 12 growth and carcass traits (P < 0.05), such as carcass weight (P = 9.68E-06), eviscerated weight (P = 3.04E-05), and semi-eviscerated weight (P = 2.14E-04). Collectively, these results provide novel insights into the genetic basis of growth in Chinese Ningdu yellow chickens and the SNP rs16400745 reported here could be incorporated into the selection programs involving this breed.

The m6A methylation regulates gonadal sex differentiation in chicken embryo

Background

As a ubiquitous reversible epigenetic RNA modification, N6-methyladenosine (m6A) plays crucial regulatory roles in multiple biological pathways. However, its functional mechanisms in sex determination and differentiation during gonadal development of chicken embryos are not clear. Therefore, we established a transcriptome-wide m6A map in the female and male chicken left gonads of embryonic day 7 (E7) by methylated RNA immunoprecipitation sequencing (MeRIP-seq) to offer insight into the landscape of m6A methylation and investigate the post-transcriptional modification underlying gonadal differentiation.

Results

The chicken embryonic gonadal transcriptome was extensively methylated. We found 15,191 and 16,111 m6A peaks in the female and male left gonads, respectively, which were mainly enriched in the coding sequence (CDS) and stop codon. Among these m6A peaks, we identified that 1013 and 751 were hypermethylated in females and males, respectively. These differential peaks covered 281 and 327 genes, such as BMP2, SMAD2, SOX9 and CYP19A1, which were primarily associated with development, morphogenesis and sex differentiation by functional enrichment. Further analysis revealed that the m6A methylation level was positively correlated with gene expression abundance. Furthermore, we found that YTHDC2 could regulate the expression of sex-related genes, especially HEMGN and SOX9, in male mesonephros/gonad mingle cells, which was verified by in vitro experiments, suggesting a regulatory role of m6A methylation in chicken gonad differentiation.

Conclusions

This work provided a comprehensive m6A methylation profile of chicken embryonic gonads and revealed YTHDC2 as a key regulator responsible for sex differentiation. Our results contribute to a better understanding of epigenetic factors involved in chicken sex determination and differentiation and to promoting the future development of sex manipulation in poultry industry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00710-6.

RNA-Seq Analysis Reveals Expression Regulatory Divergence of W-Linked Genes between Two Contrasting Chicken Breeds

Simple Summary

Understanding the mode of gene expression and regulation is essential for understanding the evolutionary process. Many previous studies tried to explain regulatory changes at the autosomal level, but little research has extended these explorations to the field of sex chromosomes due to their complex sex-limit features. Here, we first adopted an innovative method of identifying regulatory divergence of W-linked genes. Compared with cis-regulatory divergence, trans acting genes were more extensive in the W chromosome. We also found that divergent sex specific selection cannot strongly affect the expression evolution of the W chromosome. This insensitivity to selection may be one of the reasons why regulatory divergence is so small between autosomal and sex chromosomes.

Abstract

The regulation of gene expression is a complex process involving organism function and phenotypic diversity, and is caused by cis- and trans- regulation. While prior studies identified the regulatory pattern of the autosome rewiring in hybrids, the role of gene regulation in W sex chromosomes is not clear due to their degradation and sex-limit expression. Here, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibited broad differences in gender-related traits, and assessed the expression of the genes on the W chromosome to disentangle the contribution of cis- and trans-factors to expression divergence. We found that female-specific selection does not have a significant effect on W chromosome gene-expression patterns. For different tissues, there were most parental divergence expression genes in muscle, and also more heterosis compared with two other tissues. Notably, a broader pattern of trans regulation in the W chromosome was observed, which is consistent with autosomes. Taken together, this work describes the regulatory divergence of W-linked genes between two contrasting breeds and indicates sex chromosomes have a unique regulation and expression mechanism.

H3K27ac chromatin acetylation and gene expression analysis reveal sex- and situs-related differences in developing chicken gonads

BACKGROUND

Birds exhibit a unique asymmetry in terms of gonad development. The female left gonad generates a functional ovary, whereas the right gonad regresses. In males, both left and right gonads would develop into testes. How is this left/right asymmetry established only in females but not in males remains unknown. The epigenetic regulation of gonadal developmental genes may contribute to this sex disparity. The modification of histone tails such as H3K27ac is tightly coupled to chromatin activation and gene expression. To explore whether H3K27ac marked chromatin activation is involved in the asymmetric development of avian gonads, we probed genome-wide H3K27ac occupancy in left and right gonads from both sexes and related chromatin activity profile to the expression of gonadal genes. Furthermore, we validated the effect of chromatin activity on asymmetric gonadal development by manipulating the chromatin histone acetylation levels.

METHODS

The undifferentiated gonads from both sides of each sex were collected and subjected to RNA-Seq and H3K27ac ChIP-Seq experiments. Integrated analysis of gene expression and active chromatin regions were performed to identify the sex- and situs-specific regulation and expression of gonadal genes. The histone deacetylase inhibitor trichostatin A (TSA) was applied to the undifferentiated female right gonads to assess the effect of chromatin activation on gonadal gene expression and cell proliferation.

RESULTS

Even before sex differentiation, the gonads already show divergent gene expression between different sexes and between left/right sides in females. The sex-specific H3K27ac chromatin distributions coincide with the higher expression of male/female specification genes in each sex. Unexpectedly, the H3K27ac marked chromatin activation show a dramatic difference between left and right gonads in both sexes, although the left/right asymmetric gonadal development was observed only in females but not in males. In females, the side-specific H3K27ac occupancy instructs the differential expression of developmental genes between the pair of gonads and contributes to the development of left but not right gonad. However, in males, the left/right discrepancy of H3K27ac chromatin distribution does not drive the side-biased gene expression or gonad development. The TSA-induced retention of chromatin acetylation causes up-regulation of ovarian developmental genes and increases cell proliferation in the female right gonad.

CONCLUSIONS

We revealed that left/right asymmetry in H3K27ac marked chromatin activation exists in both sexes, but this discrepancy gives rise to asymmetric gonadal development only in females. Other mechanisms overriding the chromatin activation would control the symmetric development of male gonads in chicken.

Protein expression reveals a molecular sexual identity of avian primordial germ cells at pre-gonadal stages

In poultry, in vitro propagated primordial germ cells (PGCs) represent an important tool for the cryopreservation of avian genetic resources. However, several studies have highlighted sexual differences exhibited by PGCs during in vitro propagation, which may compromise their reproductive capacities. To understand this phenomenon, we compared the proteome of pregonadal migratory male (ZZ) and female (ZW) chicken PGCs propagated in vitro by quantitative proteomic analysis using a GeLC-MS/MS strategy. Many proteins were found to be differentially abundant in chicken male and female PGCs indicating their early sexual identity. Many of the proteins more highly expressed in male PGCs were encoded by genes localised to the Z sex chromosome. This suggests that the known lack of dosage compensation of the transcription of Z-linked genes between sexes persists at the protein level in PGCs, and that this may be a key factor of their autonomous sex differentiation. We also found that globally, protein differences do not closely correlate with transcript differences indicating a selective translational mechanism in PGCs. Male and female PGC expressed protein sets were associated with differential biological processes and contained proteins known to be biologically relevant for male and female germ cell development, respectively. We also discovered that female PGCs have a higher capacity to uptake proteins from the cell culture medium than male PGCs. This study presents the first evidence of an early predetermined sex specific cell fate of chicken PGCs and their sexual molecular specificities which will enable the development of more precise sex-specific in vitro culture conditions for the preservation of avian genetic resources.

Genetic Regulation of Avian Testis Development

As in other vertebrates, avian testes are the site of spermatogenesis and androgen production. The paired testes of birds differentiate during embryogenesis, first marked by the development of pre-Sertoli cells in the gonadal primordium and their condensation into seminiferous cords. Germ cells become enclosed in these cords and enter mitotic arrest, while steroidogenic Leydig cells subsequently differentiate around the cords. This review describes our current understanding of avian testis development at the cell biology and genetic levels. Most of this knowledge has come from studies on the chicken embryo, though other species are increasingly being examined. In chicken, testis development is governed by the Z-chromosome-linked DMRT1 gene, which directly or indirectly activates the male factors, HEMGN, SOX9 and AMH. Recent single cell RNA-seq has defined cell lineage specification during chicken testis development, while comparative studies point to deep conservation of avian testis formation. Lastly, we identify areas of future research on the genetics of avian testis development.

Temperature-Induced Sex Reversal in Reptiles: Prevalence, Discovery, and Evolutionary Implications

Sex reversal is the process by which an individual develops a phenotypic sex that is discordant with its chromosomal or genotypic sex. It occurs in many lineages of ectothermic vertebrates, such as fish, amphibians, and at least one agamid and one scincid reptile species. Sex reversal is usually triggered by an environmental cue that alters the genetically determined process of sexual differentiation, but it can also be caused by exposure to exogenous chemicals, hormones, or pollutants. Despite the occurrence of both temperature-dependent sex determination (TSD) and genetic sex determination (GSD) broadly among reptiles, only 2 species of squamates have thus far been demonstrated to possess sex reversal in nature (GSD with overriding thermal influence). The lack of species with unambiguously identified sex reversal is not necessarily a reflection of a low incidence of this trait among reptiles. Indeed, sex reversal may be relatively common in reptiles, but little is known of its prevalence, the mechanisms by which it occurs, or the consequences of sex reversal for species in the wild under a changing climate. In this review, we present a roadmap to the discovery of sex reversal in reptiles, outlining the various techniques that allow new occurrences of sex reversal to be identified, the molecular mechanisms that may be involved in sex reversal and how to identify them, and approaches for assessing the impacts of sex reversal in wild populations. We discuss the evolutionary implications of sex reversal and use the central bearded dragon (Pogona vitticeps) and the eastern three-lined skink (Bassiana duperreyi) as examples of how species with opposing patterns of sex reversal may be impacted differently by our rapidly changing climate. Ultimately, this review serves to highlight the importance of understanding sex reversal both in the laboratory and in wild populations and proposes practical solutions to foster future research.

Insights into Gonadal Sex Differentiation Provided by Single-Cell Transcriptomics in the Chicken Embryo

Although the genetic triggers for gonadal sex differentiation vary across species, the cell biology of gonadal development was long thought to be largely conserved. Here, we present a comprehensive analysis of gonadal sex differentiation, using single-cell sequencing in the embryonic chicken gonad during sexual differentiation. The data show that chicken embryonic-supporting cells do not derive from the coelomic epithelium, in contrast to other vertebrates studied. Instead, they derive from a DMRT1⁺/PAX2⁺/WNT4⁺/OSR1⁺ mesenchymal cell population. We find a greater complexity of gonadal cell types than previously thought, including the identification of two distinct sub-populations of Sertoli cells in developing testes and derivation of embryonic steroidogenic cells from a differentiated supporting-cell lineage. Altogether, these results indicate that, just as the genetic trigger for sex differs across vertebrate groups, cell lineage specification in the gonad may also vary substantially.

The Curious Case of Avian Sex Determination

Ioannidis and colleagues show that the gene DMRT1 is the master regulator of testis development in the chicken. Yet, remarkably, when this gene is deleted in genetic males and gonads form ovaries, the body remains male. This debunks the notion that somatic sex is driven primarily by hormones in birds.

Sexual Development of the Hermaphroditic Scallop Argopecten irradians Revealed by Morphological, Endocrine and Molecular Analysis

Simultaneous or functional hermaphrodites possessing both ovary and testis at the same time are good materials for studying sexual development. However, previous research on sex determination and differentiation was mainly conducted in gonochoristic species and studies on simultaneous hermaphrodites are still limited. In this study, we conducted a combined morphological, endocrine and molecular study on the gonadal development of a hermaphroditic scallop Argopecten irradians aged 2–10 month old. Morphological analysis showed that sex differentiation occurred at 6 months of age. By examining the dynamic changes of progesterone, testosterone and estradiol, we found testosterone and estradiol were significantly different between the ovaries and testes almost throughout the whole process, suggesting the two hormones may be involved in scallop sex differentiation. In addition, we identified two critical sex-related genes FoxL2 and Dmrt1L, and investigated their spatiotemporal expression patterns. Results showed that FoxL2 and Dmrt1L were female- and male-biased, respectively, and mainly localized in the germ cells and follicular cells, indicating their feasibility as molecular markers for early identification of sex. Further analysis on the changes of FoxL2 and Dmrt1L expression in juveniles showed that significant sexual dimorphic expression of FoxL2 occurred at 2 months of age, earlier than that of Dmrt1L. Moreover, FoxL2 expression was significantly correlated with estradiol/testosterone ratio (E₂/T). All these results indicated that molecular sex differentiation occurs earlier than morphological sex differentiation, and FoxL2 may be a key driver that functions through regulating sex steroid hormones in the scallop. This study will deepen our understanding of the molecular mechanism underlying sex differentiation and development in spiralians.

Sex identification in embryos and adults of Darwin's finches

Darwin’s finches are an iconic example of adaptive radiation and evolution under natural selection. Comparative genetic studies using embryos of Darwin’s finches have shed light on the possible evolutionary processes underlying the speciation of this clade. Molecular identification of the sex of embryonic samples is important for such studies, where this information often cannot be inferred otherwise. We tested a fast and simple chicken embryo protocol to extract DNA from Darwin’s finch embryos. In addition, we applied minor modifications to two of the previously reported PCR primer sets for CHD1, a gene used for sexing adult passerine birds. The sex of all 29 tested embryos of six species of Darwin’s finches was determined successfully by PCR, using both primer sets. Next to embryos, hatchlings and fledglings are also impossible to distinguish visually. This extends to juveniles of sexually dimorphic species which are yet to moult in adult-like plumage and beak colouration. Furthermore, four species of Darwin’s finches are monomorphic, males and females looking alike. Therefore, sex assessment in the field can be a source of error, especially with respect to juveniles and mature monomorphic birds outside of the mating season. We caught 567 juveniles and adults belonging to six species of Darwin’s finches and only 44% had unambiguous sex-specific morphology. We sexed 363 birds by PCR: individuals sexed based on marginal sex specific morphological traits; and birds which were impossible to classify in the field. PCR revealed that for birds with marginal sex specific traits, sexing in the field produced a 13% error rate. This demonstrates that PCR based sexing can improve field studies on Darwin’s finches, especially when individuals with unclear sex-related morphology are involved. The protocols used here provide an easy and reliable way to sex Darwin’s finches throughout ontogeny, from embryos to adults.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

Expression profiling of sexually dimorphic genes in the Japanese quail, Coturnix japonica

Research on avian sex determination has focused on the chicken. In this study, we established the utility of another widely used animal model, the Japanese quail (Coturnix japonica), for clarifying the molecular mechanisms underlying gonadal sex differentiation. In particular, we performed comprehensive gene expression profiling of embryonic gonads at three stages (HH27, HH31 and HH38) by mRNA-seq. We classified the expression patterns of 4,815 genes into nine clusters according to the extent of change between stages. Cluster 2 (characterized by an initial increase and steady levels thereafter), including 495 and 310 genes expressed in males and females, respectively, contained five key genes involved in gonadal sex differentiation. A GO analysis showed that genes in this cluster are related to developmental processes including reproductive structure development and developmental processes involved in reproduction were significant, suggesting that expression profiling is an effective approach to identify novel candidate genes. Based on RNA-seq data and in situ hybridization, the expression patterns and localization of most key genes for gonadal sex differentiation corresponded well to those of the chicken. Our results support the effectiveness of the Japanese quail as a model for studies gonadal sex differentiation in birds.

Transcriptome Sequencing and Comparative Analysis of Amphoteric ESCs and PGCs in Chicken (Gallus gallus)

Chicken (Gallus gallus) pluripotent embryonic stem cells (ESCs) and primordial germ cells (PGCs) can be broadly applied in the research of developmental and embryonic biology, but the difference between amphoteric ESCs and PGCs is still elusive. This study determined the sex of collected samples by identifying specific sex markers via polymerase chain reaction (PCR) and fluorescence activated cell sorter (FACS). RNA-seq was utilized to investigate the transcriptomic profile of amphoteric ESCs and PGCs in chicken. The results showed no significant differentially expressed genes (DEGs) in amphoteric ESCs and 227 DEGs exhibited in amphoteric PGCs. Moreover, those 227 DEGs were mainly enriched in 17 gene ontology (GO) terms and 27 pathways according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Furthermore, qRT-PCR was performed to verify RNA-seq results, and the results demonstrated that Notch1 was highly expressed in male PGCs. In summary, our results provided a knowledge base of chicken amphoteric ESCs and PGCs, which is helpful for future research in relevant biological processes.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Identification of Y chromosome markers in the eastern three-lined skink (Bassiana duperreyi) using in silico whole genome subtraction

Background

Homologous sex chromosomes can differentiate over time because recombination is suppressed in the region of the sex determining locus, leading to the accumulation of repeats, progressive loss of genes that lack differential influence on the sexes and sequence divergence on the hemizygous homolog. Divergence in the non-recombining regions leads to the accumulation of Y or W specific sequence useful for developing sex-linked markers. Here we use in silico whole-genome subtraction to identify putative sex-linked sequences in the scincid lizard Bassiana duperreyi which has heteromorphic XY sex chromosomes.

Results

We generated 96.7 × 10⁹ 150 bp paired-end genomic sequence reads from a XY male and 81.4 × 10⁹ paired-end reads from an XX female for in silico whole genome subtraction to yield Y enriched contigs. We identified 7 reliable markers which were validated as Y chromosome specific by polymerase chain reaction (PCR) against a panel of 20 males and 20 females.

Conclusions

The sex of B. duperreyi can be reversed by low temperatures (XX genotype reversed to a male phenotype). We have developed sex-specific markers to identify the underlying genotypic sex and its concordance or discordance with phenotypic sex in wild populations of B. duperreyi. Our pipeline can be applied to isolate Y or W chromosome-specific sequences of any organism and is not restricted to sequence residing within single-copy genes. This study greatly improves our knowledge of the Y chromosome in B. duperreyi and will enhance future studies of reptile sex determination and sex chromosome evolution.

Transcriptome analysis of genes related to gonad differentiation and development in Muscovy ducks

BACKGROUND

Sex-related genes play a crucial role in gonadal differentiation into testes or ovaries. However, the genetic control of gonadal differentiation in Muscovy ducks remains unknown. Therefore, the objective of our study was to screen new candidate genes associated with ovarian and testicular development.

RESULTS

In this study, 24 males before gonadal differentiation (MB), 24 females before gonadal differentiation (FB), 24 males after gonadal differentiation (MA) and 24 females after gonadal differentiation (FA) were selected from Putian Muscovy ducks, forming 4 groups. RNA-Seq revealed 101.76 Gb of clean reads and 2800 differentially expressed genes (DEGs), including 46 in MB vs FB, 609 in MA vs FA, 1027 in FA vs FB, and 1118 in MA vs MB. A total of 146 signalling pathways were enriched by KEGG analysis, among which 20, 108, 108 and 116 signalling pathways were obtained in MB vs FB, MA vs MB, MA vs FA and FA vs FB, respectively. In further GO and KEGG analyses, a total of 21 candidate genes related to gonad differentiation and development in Muscovy ducks were screened. Among these, 9 genes were involved in the differentiation and development of the testes, and 12 genes were involved in the differentiation and development of the ovaries. In addition, RNA-Seq data revealed 2744 novel genes.

CONCLUSIONS

RNA-Seq data revealed 21 genes related to gonadal differentiation and development in Muscovy ducks. We further identified 12 genes, namely, WNT5B, HTRA3, RSPO3, BMP3, HNRNPK, NIPBL, CREB3L4, DKK3, UBE2R2, UBPL3KCMF1, ANXA2, and OSR1, involved in the differentiation and development of ovaries. Moreover, 9 genes, namely, TTN, ATP5A1, DMRT1, DMRT3, AMH, MAP3K1, PIK3R1, AGT and ADAMTSL1, were related to the differentiation and development of testes. Moreover, after gonadal differentiation, DMRT3, AMH, PIK3R1, ADAMTSL1, AGT and TTN were specifically highly expressed in males. WNT5B, ANXA2 and OSR1 were specifically highly expressed in females. These results provide valuable information for studies on the sex control of Muscovy ducks and reveal novel candidate genes for the differentiation and development of testes and ovaries.

Comprehensive analysis of miRNAs, lncRNAs, and mRNAs reveals potential players of sexually dimorphic and left-right asymmetry in chicken gonad during gonadal differentiation

Despite thousands of sex-biased genes being found in chickens, the genetic control of sexually dimorphic and left-right asymmetry during gonadal differentiation is not yet completely understood. This study aimed to identify microRNAs (miRNAs), long noncoding RNAs (lncRNAs), messenger RNAs (mRNAs), and signaling pathways during gonadal differentiation in chick embryos (day 6/stage 29). The left and right gonads were collected for RNA sequencing. Sex-biased, side-biased miRNAs, lncRNAs, mRNAs, and shared differentially expressed miRNAs (DEmiRNA)–differentially expressed mRNAs (DEmRNA)–differentially expressed lncRNAs (DElncRNA) interaction networks were performed. A total of 8 DEmiRNAs, 183 DElncRNAs, and 123 DEmRNAs were identified for the sex-biased genes, and 7 DEmiRNAs, 189 DElncRNAs, and 183 DEmRNAs for the side-biased genes. The results of quantitative real-time PCR were generally consistent with the RNA-sequencing results. The study suggested that miRNAs and lncRNAs regulation were novel gene-specific dosage compensation mechanism and they could contribute to left-right asymmetry of chicken, but sex-biased and side-biased miRNAs, lncRNAs, and mRNAs were independent of each other. The competing endogenous RNA (ceRNA) networks showed that 17 target pairs including miR-7b (CYP19A1, FSHR, GREB1, STK31, CORIN, and TDRD9), miR-211 (FSHR, GREB1, STK31, CORIN, and TDRD9), miR-204 (FSHR, GREB1, CORIN, and TDRD9), and miR-302b-5p (CYP19A1 and TDRD9) may play crucial roles in ovarian development. These analyses provide new clues to uncover molecular mechanisms and signaling networks of ovarian development.

Expression Profile of Chicken Sex Chromosome Gene BTF3 is Linked to Gonadal Phenotype

In birds, the female is heterogametic (ZW) and the male homogametic (ZZ). The small W chromosome comprises only 28 protein coding genes (homologues to Z chromosome counterparts) and a number of repeat regions. Here, we report our analysis of one of these genes, BTF3 (basic transcription factor 3), which exhibits differential expression during gonadogenesis. We measured RNA levels of both Z and W homologues and BTF3 protein levels in male and female gonads during development of the chicken embryo. In addition, BTF3 RNA and protein levels were compared in female gonads (ovary) and in female gonads following treatment to induce sex reversal (testis). Combined BTF3 RNA levels were higher in female gonads than male gonads, while BTF3-Z was expressed at similar levels in males and females. Surprisingly, BTF3 protein levels were higher in male gonads than female gonads at embryonic day 6 (E6), suggesting translational rather than transcriptional regulation. BTF3 protein was expressed in both somatic and germ cells and was restricted to the medulla of the developing ovary in females and the sex cords of the developing testis in males. In addition, in gonadal sex-reversed females, RNA and protein levels of BTF3 were similar to those normally found in male gonads, suggesting that BTF3 expression correlated with the gonadal phenotype.

Comparative Transcriptomics Analyses across Species, Organs, and Developmental Stages Reveal Functionally Constrained lncRNAs

The functionality of long noncoding RNAs (lncRNAs) is disputed. In general, lncRNAs are under weak selective pressures, suggesting that the majority of lncRNAs may be nonfunctional. However, although some surveys showed negligible phenotypic effects upon lncRNA perturbation, key biological roles were demonstrated for individual lncRNAs. Most lncRNAs with proven functions were implicated in gene expression regulation, in pathways related to cellular pluripotency, differentiation, and organ morphogenesis, suggesting that functional lncRNAs may be more abundant in embryonic development, rather than in adult organs. To test this hypothesis, we perform a multidimensional comparative transcriptomics analysis, across five developmental time points (two embryonic stages, newborn, adult, and aged individuals), four organs (brain, kidney, liver, and testes), and three species (mouse, rat, and chicken). We find that, overwhelmingly, lncRNAs are preferentially expressed in adult and aged testes, consistent with the presence of permissive transcription during spermatogenesis. LncRNAs are often differentially expressed among developmental stages and are less abundant in embryos and newborns compared with adult individuals, in agreement with a requirement for tighter expression control and less tolerance for noisy transcription early in development. For differentially expressed lncRNAs, we find that the patterns of expression variation among developmental stages are generally conserved between mouse and rat. Moreover, lncRNAs expressed above noise levels in somatic organs and during development show higher evolutionary conservation, in particular, at their promoter regions. Thus, we show that functionally constrained lncRNA loci are enriched in developing organs, and we suggest that many of these loci may function in an RNA-independent manner.

Regional epigenetic differentiation of the Z Chromosome between sexes in a female heterogametic system

In male heterogametic systems, the X Chromosome is epigenetically differentiated between males and females, to facilitate dosage compensation. For example, the X Chromosome in female mammals is largely inactivated. Relative to well-studied male heterogametic systems, the extent of epigenetic differentiation between male and female Z Chromosomes in female heterogametic species, which often lack complete dosage compensation, is poorly understood. Here, we examined the chromosomal DNA methylation landscapes of male and female Z Chromosomes in two distantly related avian species, namely chicken and white-throated sparrow. We show that, in contrast to the pattern in mammals, male and female Z Chromosomes in these species exhibit highly similar patterns of DNA methylation, which is consistent with weak or absent dosage compensation. We further demonstrate that the epigenetic differences between male and female chicken Z Chromosomes are localized to a few regions, including a previously identified male hypermethylated region 1 (MHM1; CGNC: 80601). We discovered a novel region with elevated male-to-female methylation ratios on the chicken Z Chromosome (male hypermethylated region 2 [MHM2]; CGNC: 80602). The MHM1 and MHM2, despite little sequence similarity between them, bear similar molecular features that are likely associated with their functions. We present evidence consistent with female hypomethylation of MHMs and up-regulation of nearby genes. Therefore, despite little methylation differentiation between sexes, extremely localized DNA methylation differences between male and female chicken Z Chromosomes have evolved and affect expression of nearby regions. Our findings offer new insights into epigenetic regulation of gene expression between sexes in female heterogametic systems.

Miniscule differences between sex chromosomes in the giant genome of a salamander

In the Mexican axolotl (Ambystoma mexicanum), sex is determined by a single Mendelian factor, yet its sex chromosomes do not exhibit morphological differentiation typical of many vertebrate taxa that possess a single sex-determining locus. As sex chromosomes are theorized to differentiate rapidly, species with undifferentiated sex chromosomes provide the opportunity to reconstruct early events in sex chromosome evolution. Whole genome sequencing of 48 salamanders, targeted chromosome sequencing and in situ hybridization were used to identify the homomorphic sex chromosome that carries an A. mexicanum sex-determining factor and sequences that are present only on the W chromosome. Altogether, these sequences cover ~300 kb of validated female-specific (W chromosome) sequence, representing ~1/100,000^th of the 32 Gb genome. Notably, a recent duplication of ATRX, a gene associated with mammalian sex-determining pathways, is one of few functional (non-repetitive) genes identified among these W-specific sequences. This duplicated gene (ATRW) was used to develop highly predictive markers for diagnosing sex and represents a strong candidate for a recently-acquired sex determining locus (or sexually antagonistic gene) in A. mexicanum.

Gonadal and Endocrine Analysis of a Gynandromorphic Chicken

Birds have a ZZ male and ZW female sex chromosome system. The relative roles of genetics and hormones in regulating avian sexual development have been revealed by studies on gynandromorphs. Gynandromorphs are rare bilateral sex chimeras, male on one side of the body and female on the other. We examined a naturally occurring gynandromorphic chicken that was externally male on the right side of the body and female on the left. The bird was diploid but with a mix of ZZ and ZW cells that correlated with the asymmetric sexual phenotype. The male side was 96% ZZ, and the female side was 77% ZZ and 23% ZW. The gonads of this bird at sexual maturity were largely testicular. The right gonad was a testis, with SOX9+ Sertoli cells, DMRT1+ germ cells, and active spermatogenesis. The left gonad was primarily testicular, but with some peripheral aromatase-expressing follicles. The bird had low levels of serum estradiol and high levels of testosterone, as expected for a male. Despite the low percentage of ZW cells on that side, the left side had female sex-linked feathering, smaller muscle mass, smaller leg and spur, and smaller wattle than the male side. This indicates that these sexually dimorphic structures must be at least partly independent of sex steroid effects. Even a small percentage of ZW cells appears sufficient to support female sexual differentiation. Given the lack of chromosome-wide dosage compensation in birds, various sexually dimorphic features may arise due to Z-gene dosage differences between the sexes.

FOXL2 and DMRT1L Are Yin and Yang Genes for Determining Timing of Sex Differentiation in the Bivalve Mollusk Patinopecten yessoensis

Sex determination and differentiation have long been a research hotspot in metazoans. However, little is known about when and how sex differentiation occurs in most mollusks. In this study, we conducted a combined morphological and molecular study on sex differentiation in the Yesso scallop Patinopecten yessoensis. Histological examination on gonads from 5- to 13-month-old juveniles revealed that the morphological sex differentiation occurred at 10 months of age. To determine the onset of molecular sex differentiation, molecular markers were screened for early identification of sex. The gonadal expression profiles of eight candidate genes for sex determination or differentiation showed that only two genes displayed sexually dimorphic expression, with FOXL2 being abundant in ovaries and DMRT1L in testes. In situ hybridization revealed that both of them were detected in germ cells and follicle cells. We therefore developed LOG₁₀(DMRT1L/FOXL2) for scallop sex identification and confirmed its feasibility in differentiated individuals. By tracing its changes in 5- to 13-month-old juveniles, molecular sex differentiation time was determined: some scallops differentiate early in September when they are 7 months old, and some do late in December when they are 10 months old. Two kinds of coexpression patterns were found between FOXL2 and DMRT1L: expected antagonism after differentiation and unexpected coordination before differentiation. Our results revealed that scallop sex differentiation co-occurs with the formation of follicles, and molecular sex differentiation is established prior to morphological sex differentiation. Our study will assist in a better understanding of the molecular mechanism underlying bivalve sex differentiation.

A new approach for comprehensively describing heterogametic sex chromosomes

Notwithstanding the rapid developments in sequencing techniques, Y and W sex chromosomes have still been mostly excluded from whole genome sequencing projects due to their high repetitive DNA content. Therefore, Y and W chromosomes are poorly described in most species despite their biological importance. Several methods were developed for identifying Y or W-linked sequences among unmapped scaffolds. However, it is not enough to discover functional regions from short unmapped scaffolds. Here, we provide a new and simple strategy based on k-mer comparison for comprehensive analysis of the W chromosome in Bombyx mori. Using this novel method, we effectively assembled de novo 1281 W-derived genome contigs (totaling 1.9 Mbp), and identified 156 W-linked transcript RNAs and 345 W-linked small RNAs. This method will help in the elucidation of mechanisms of sexual development and exploration of W chromosome biological functions, and provide insights into the evolution of sex chromosomes. Moreover, we showed this method can be employed in identifying heterogametic sex chromosomes (W and Y chromosomes) in many other species where genomic information is still scarce.

Size and Content of the Sex-Determining Region of the Y Chromosome in Dioecious Mercurialis annua, a Plant with Homomorphic Sex Chromosomes

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua, a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

Chromosomal Context Affects the Molecular Evolution of Sex-linked Genes and Their Autosomal Counterparts in Turtles and Other Vertebrates

Sex chromosomes evolve differently from autosomes because natural selection acts distinctly on them given their reduced recombination and smaller population size. Various studies of sex-linked genes compared with different autosomal genes within species support these predictions. Here, we take a novel alternative approach by comparing the rate of evolution between subsets of genes that are sex-linked in selected reptiles/vertebrates and the same genes located in autosomes in other amniotes. We report for the first time the faster evolution of Z-linked genes in a turtle (the Chinese softshell turtle Pelodiscus sinensis) relative to autosomal orthologs in other taxa, including turtles with temperature-dependent sex determination (TSD). This faster rate was absent in its close relative, the spiny softshell turtle (Apalone spinifera), thus revealing important lineage effects, and was only surpassed by mammalian-X linked genes. In contrast, we found slower evolution of X-linked genes in the musk turtle Staurotypus triporcatus (XX/XY) and homologous Z-linked chicken genes. TSD lineages displayed overall faster sequence evolution than taxa with genotypic sex determination (GSD), ruling out global effects of GSD on molecular evolution beyond those by sex-linkage. Notably, results revealed a putative selective sweep around two turtle genes involved in vertebrate gonadogenesis (Pelodiscus-Z-linked Nf2 and Chrysemys-autosomal Tspan7). Our observations reveal important evolutionary changes at the gene level mediated by chromosomal context in turtles despite their low overall evolutionary rate and illuminate sex chromosome evolution by empirically testing expectations from theoretical models. Genome-wide analyses are warranted to test the generality and prevalence of the observed patterns.

Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation

Sexual dimorphism depends on sex-biased gene expression, but the contributions of microRNAs (miRNAs) have not been globally assessed. We therefore produced an extensive small RNA sequencing data set to analyze male and female miRNA expression profiles in mouse, opossum, and chicken. Our analyses uncovered numerous cases of somatic sex-biased miRNA expression, with the largest proportion found in the mouse heart and liver. Sex-biased expression is explained by miRNA-specific regulation, including sex-biased chromatin accessibility at promoters, rather than piggybacking of intronic miRNAs on sex-biased protein-coding genes. In mouse, but not opossum and chicken, sex bias is coordinated across tissues such that autosomal testis-biased miRNAs tend to be somatically male-biased, whereas autosomal ovary-biased miRNAs are female-biased, possibly due to broad hormonal control. In chicken, which has a Z/W sex chromosome system, expression output of genes on the Z Chromosome is expected to be male-biased, since there is no global dosage compensation mechanism that restores expression in ZW females after almost all genes on the W Chromosome decayed. Nevertheless, we found that the dominant liver miRNA, miR-122-5p, is Z-linked but expressed in an unbiased manner, due to the unusual retention of a W-linked copy. Another Z-linked miRNA, the male-biased miR-2954-3p, shows conserved preference for dosage-sensitive genes on the Z Chromosome, based on computational and experimental data from chicken and zebra finch, and acts to equalize male-to-female expression ratios of its targets. Unexpectedly, our findings thus establish miRNA regulation as a novel gene-specific dosage compensation mechanism.

Sex Reversal and Comparative Data Undermine the W Chromosome and Support Z-linked DMRT1 as the Regulator of Gonadal Sex Differentiation in Birds

The exact genetic mechanism regulating avian gonadal sex differentiation has not been completely resolved. The most likely scenario involves a dosage mechanism, whereby the Z-linked DMRT1 gene triggers testis development. However, the possibility still exists that the female-specific W chromosome may harbor an ovarian determining factor. In this study, we provide evidence that the universal gene regulating gonadal sex differentiation in birds is Z-linked DMRT1 and not a W-linked (ovarian) factor. Three candidate W-linked ovarian determinants are HINTW, female-expressed transcript 1 (FET1), and female-associated factor (FAF). To test the association of these genes with ovarian differentiation in the chicken, we examined their expression following experimentally induced female-to-male sex reversal using the aromatase inhibitor fadrozole (FAD). Administration of FAD on day 3 of embryogenesis induced a significant loss of aromatase enzyme activity in female gonads and masculinization. However, expression levels of HINTW, FAF, and FET1 were unaltered after experimental masculinization. Furthermore, comparative analysis showed that FAF and FET1 expression could not be detected in zebra finch gonads. Additionally, an antibody raised against the predicted HINTW protein failed to detect it endogenously. These data do not support a universal role for these genes or for the W sex chromosome in ovarian development in birds. We found that DMRT1 (but not the recently identified Z-linked HEMGN gene) is male upregulated in embryonic zebra finch and emu gonads, as in the chicken. As chicken, zebra finch, and emu exemplify the major evolutionary clades of birds, we propose that Z-linked DMRT1, and not the W sex chromosome, regulates gonadal sex differentiation in birds.

SuperTranscripts: a data driven reference for analysis and visualisation of transcriptomes

Numerous methods have been developed to analyse RNA sequencing (RNA-seq) data, but most rely on the availability of a reference genome, making them unsuitable for non-model organisms. Here we present superTranscripts, a substitute for a reference genome, where each gene with multiple transcripts is represented by a single sequence. The Lace software is provided to construct superTranscripts from any set of transcripts, including de novo assemblies. We demonstrate how superTranscripts enable visualisation, variant detection and differential isoform detection in non-model organisms. We further use Lace to combine reference and assembled transcriptomes for chicken and recover hundreds of gaps in the reference genome.

Gene Expression Profiling Reveals Potential Players of Left-Right Asymmetry in Female Chicken Gonads

Most female birds develop only a left ovary, whereas males develop bilateral testes. The mechanism underlying this process is still not completely understood. Here, we provide a comprehensive transcriptional analysis of female chicken gonads and identify novel candidate side-biased genes. RNA-Seq analysis was carried out on total RNA harvested from the left and right gonads on embryonic day 6 (E6), E12, and post-hatching day 1 (D1). By comparing the gene expression profiles between the left and right gonads, 347 differentially expressed genes (DEGs) were obtained on E6, 3730 were obtained on E12, and 2787 were obtained on D1. Side-specific genes were primarily derived from the autosome rather than the sex chromosome. Gene ontology and pathway analysis showed that the DEGs were most enriched in the Piwi-interactiing RNA (piRNA) metabolic process, germ plasm, chromatoid body, P granule, neuroactive ligand-receptor interaction, microbial metabolism in diverse environments, and methane metabolism. A total of 111 DEGs, five gene ontology (GO) terms, and three pathways were significantly different between the left and right gonads among all the development stages. We also present the gene number and the percentage within eight development-dependent expression patterns of DEGs in the left and right gonads of female chicken.

Characterization of genome-wide segmental duplications reveals a common genomic feature of association with immunity among domestic animals

BACKGROUND

Segmental duplications (SDs) commonly exist in plant and animal genomes, playing crucial roles in genomic rearrangement, gene innovation and the formation of copy number variants. However, they have received little attention in most livestock species.

RESULTS

Aiming at characterizing SDs across the genomes of diverse livestock species, we mapped genome-wide SDs of horse, rabbit, goat, sheep and chicken, and also enhanced the existing SD maps of cattle and pig genomes based on the most updated genome assemblies. We adopted two different detection strategies, whole genome analysis comparison and whole genome shotgun sequence detection, to pursue more convincing findings. Accordingly we identified SDs for each species with the length of from 21.7 Mb to 164.1 Mb, and 807 to 4,560 genes were harboured within the SD regions across different species. More interestingly, many of these SD-related genes were involved in the process of immunity and response to external stimuli. We also found the existence of 59 common genes within SD regions in all studied species except goat. These common genes mainly consisted of both UDP glucuronosyltransferase and Interferon alpha families, implying the connection between SDs and the evolution of these gene families.

CONCLUSIONS

Our findings provide insights into livestock genome evolution and offer rich genomic sources for livestock genomic research.

Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development

Vertebrate sexual fate is decided primarily by the individual's genotype (GSD), by the environmental temperature during development (TSD), or both. Turtles exhibit TSD and GSD, making them ideal to study the evolution of sex determination. Here we analyze temperature-specific gonadal transcriptomes (RNA-sequencing validated by qPCR) of painted turtles (Chrysemys picta TSD) before and during the thermosensitive period, and at equivalent stages in soft-shell turtles (Apalone spinifera-GSD), to test whether TSD's and GSD's transcriptional circuitry is identical but deployed differently between mechanisms. Our data show that most elements of the mammalian urogenital network are active during turtle gonadogenesis, but their transcription is generally more thermoresponsive in TSD than GSD, and concordant with their sex-specific function in mammals [e.g., upregulation of Amh, Ar, Esr1, Fog2, Gata4, Igf1r, Insr, and Lhx9 at male-producing temperature, and of β-catenin, Foxl2, Aromatase (Cyp19a1), Fst, Nf-kb, Crabp2 at female-producing temperature in Chrysemys]. Notably, antagonistic elements in gonadogenesis (e.g., β-catenin and Insr) were thermosensitive only in TSD early-embryos. Cirbp showed warm-temperature upregulation in both turtles disputing its purported key TSD role. Genes that may convert thermal inputs into sex-specific development (e.g., signaling and hormonal pathways, RNA-binding and heat-shock) were differentially regulated. Jak-Stat, Nf-κB, retinoic-acid, Wnt, and Mapk-signaling (not Akt and Ras-signaling) potentially mediate TSD thermosensitivity. Numerous species-specific ncRNAs (including Xist) were differentially-expressed, mostly upregulated at colder temperatures, as were unannotated loci that constitute novel TSD candidates. Cirbp showed warm-temperature upregulation in both turtles. Consistent transcription between turtles and alligator revealed putatively-critical reptilian TSD elements for male (Sf1, Amh, Amhr2) and female (Crabp2 and Hspb1) gonadogenesis. In conclusion, while preliminary, our data helps illuminate the regulation and evolution of vertebrate sex determination, and contribute genomic resources to guide further research into this fundamental biological process.