Search

Background: The blunt snout bream (Megalobrama amblycephala) is an important economic freshwater fish in China with tender flesh and high nutritional value. With the cultivation of superior new varieties and the expansion of breeding scale, it becomes imperative to employ sex-control technology to cultivate monosexual populations of M. amblycephala, thereby preventing the deterioration of desirable traits. The development of specific markers capable of accurately identifying the sex of M. amblycephala would facilitate the determination of the genetic sex of the breeding population before gonad maturation, thereby expediting the processes of sex-controlled breeding of M. amblycephala. Results: A whole-genome re-sequencing was performed for 116 females and 141 males M. amblycephala collected from nine populations. Seven candidate male-specific sequences were identified through comparative analysis of male and female genomes, which were further compared with the sequencing data of 257 individuals, and finally three male-specific sequences were generated. These three sequences were further validated by PCR amplification in 32 males and 32 females to confirm their potential as male-specific molecular markers for M. amblycephala. One of these markers showed potential applicability in M. pellegrini as well, enabling males to be identified using this specific molecular marker. Conclusions: The study provides a high-efficiency and cost-effective approach for the genetic sex identification in two species of Megalobrama. The developed markers in this study have great potential in facilitating sex-controlled breeding of M. amblycephala and M. pellegrini, while also contributing valuable insights into the underlying mechanisms of fish sex determination. [ABSTRACT FROM AUTHOR]

The Sox gene family, a collection of transcription factors widely distributed throughout the animal kingdom, plays a crucial role in numerous developmental processes. Echinoderms occupy a pivotal position in many research fields, such as neuroscience, sex determination and differentiation, and embryonic development. However, to date, no comprehensive study has been conducted to characterize and analyze Sox genes in echinoderms. In the present study, the evolution and expression of Sox family genes across 11 echinoderms were analyzed using bioinformatics methods. The results revealed a total of 70 Sox genes, with counts ranging from 5 to 8 across different echinoderms. Phylogenetic analysis revealed that the identified Sox genes could be categorized into seven distinct classes: the SoxB1 class, SoxB2 class, SoxC class, SoxD class, SoxE class, SoxF class and SoxH class. Notably, the SoxB1, SoxB2, and SoxF genes were ubiquitously present in all the echinoderms studied, which suggests that these genes may be conserved in echinoderms. The spatiotemporal expression patterns observed for Sox genes in the three echinoderms indicated that various Sox members perform distinct functional roles. Notably, SoxB1 is likely involved in echinoderm ovary development, while SoxH may play a crucial role in testis development in starfish and sea cucumber. In general, the present investigation provides a molecular foundation for exploring the Sox gene in echinoderms, providing a valuable resource for future phylogenetic and genomic studies. [ABSTRACT FROM AUTHOR]

Background: Expansion of genomic resources for the Pacific white shrimp (Litopenaeus vannamei), such as the construction of dense genetic linkage maps, is crucial for the application of genomic tools in order to improve economically relevant traits. Sexual dimorphism exists in Pacific white shrimp, and the mapping of the sex-determination region in this species may help in future reproductive applications. We have constructed male, female, and sex-averaged high-density genetic maps using a 50 K single-nucleotide polymorphism (SNP) array, followed by a genome-wide association study (GWAS) to identify genomic regions associated with sex in white shrimp. Results: The genetic map yielded 15,256 SNPs assigned to 44 linkage groups (LG). The lengths of the male, female, and sex-averaged maps were 5,741.36, 5,461.20 and 5,525.26 cM, respectively. LG18 was found to be the largest for both sexes, whereas LG44 was the shortest for males and LG31 for females. A sex-determining region was found in LG31 with 21 statistically significant SNPs. The most important SNP was previously identified as a sex-linked marker and was able to identify 99% of the males and 88% of the females. Although other significant markers had a lower ability to determine sex, putative genes were intercepted or close to them. The oplophorus-luciferin 2-monooxygenase, serine/arginine repetitive matrix protein and spermine oxidase genes were identified as candidates with possible participation in important processes of sexual differentiation in shrimp. Conclusions: Our results provide novel genomic resources for shrimp, including a high-density linkage map and new insights into the sex-determining region in L. vannamei, which may be usefulfor future genetics and reproduction applications. [ABSTRACT FROM AUTHOR]

Abstract Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

Background: The Asian yellow pond turtle (Mauremys mutica) is an important commercial freshwater aquaculture species in China. This species is a highly sexually dimorphic species, with males growing at a faster rate than females and exhibits temperature-dependent sex determination (TSD), in which the incubation temperature during embryonic development determines the sexual fate. However, the mechanisms of the sex determination or sex differentiation in the Asian yellow pond turtle are remain a mystery. Results: Temperature-specific gonadal transcriptomics of the Asian yellow pond turtle were performed during the thermosensitive period (stage 15) using RNA-seq technology to identify candidate genes that initiate gonadal differentiation. We uncovered candidates that were the first to respond to temperature. These candidates were sexually dimorphic in expression, reflecting differences in gonadal (Cirbp, Runx1) and germline differentiation (Vasa, Nanos1, Piwil2), gametogenesis (Hmgb3, Zar1, Ovoinhibitor-like, Kif4), steroid hormone biosynthesis (Hsd17b5, Hsd17b6), heat shock (Dnajb6, Hsp90b1, Hsp90aa1) and transient receptor potential channel genes (Trpm1, Trpm4, Trpm6, Trpv1). Conclusions: Our work will provide important genetic information to elucidate the mechanisms of sex control in the Asian yellow pond turtles, and will contribute important genetic resources for further studies of temperature-dependent sex determination in turtles. [ABSTRACT FROM AUTHOR]

Abstract 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.

Abstract Background Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish. Results In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways. Conclusions Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish.

In recent years, some common themes in the development of sex-specific traits in different animal lineages have started to emerge since the discovery of the Dmrt (doublesex-mab3-related transcription factor gene) genes. Bivalves are characterized by a diversity of sexual systems, including simultaneous hermaphroditism, sequential hermaphroditism, and strict gonochorism. However, to date, no research has focused on the genome-wide characterization and analysis of Dmrt genes in bivalves. In this study, the identification and analysis of Dmrt genes in 15 bivalves were performed using bioinformatics methods. A total of 55 Dmrt genes were retrieved in the studied bivalve genomes. The number of Dmrt genes in different species ranged from 3 to 5. The phylogenetic tree showed that Dmrt genes in bivalves can be subdivided into 5 classes: the Dmrt2-like class, Dmrt3-like class, Dmrt4/5-like class, Dsx-like class, and scallop-specific Dmrt class. The Ka/Ks ratios suggested that all Dmrt classes underwent purifying selection pressure. Furthermore, the spatiotemporal expression of Dmrt genes in four bivalve species suggested that different Dmrt genes may have different functions, and scallop-specific Dmrt genes may play a key role in sex determination/differentiation. In general, this study provides a molecular basis for in-depth examination of the functions of Dmrt genes and phylogenomic analyses in bivalves. [ABSTRACT FROM AUTHOR]

Background: Hong Kong catfish (Clarias fuscus) is an ecologically and economically important species that is widely distributed in freshwater regions of southern China. Hong Kong catfish has significant sexual growth dimorphism. The genome assembly of the Hong Kong catfish would facilitate study of the sex determination and evolution mechanism of the species. Results: The first high-quality chromosome-level genome of the Hong Kong catfish was constructed. The total genome was 933.4 Mb, with 416 contigs and a contig N50 length of 8.52 Mb. Using high-throughput chromosome conformation capture (Hi-C) data, the genome assembly was divided into 28 chromosomes with a scaffold N50 length of 36.68 Mb. A total of 23,345 protein-coding genes were predicted in the genome, and 94.28% of the genes were functionally annotated in public databases. Phylogenetic analysis indicated that C. fuscus and Clarias magur diverged approximately 63.7 million years ago. The comparative genome results showed that a total of 60 unique, 353 expanded and 851 contracted gene families were identified in Hong Kong catfish. A sex-linked quantitative trait locus identified in a previous study was located in a sex-determining region of 30.26 Mb (0.02 to 30.28 Mb) on chromosome 13 (Chr13), the predicted Y chromosome. This QTL region contained 785 genes, of which 18 were identified as sex-related genes. Conclusions: This study is the first to report the chromosome-level genome assembly of Hong Kong catfish. The study provides an excellent genetic resource that will facilitate future studies of sex determination mechanisms and evolution in fish. [ABSTRACT FROM AUTHOR]

Background: Red-tail catfish (Hemibagrus wyckioides) is an important commercially farmed catfish in southern China. Males of red-tail catfish grow faster than females, suggesting that all-male catfish will produce more significant economic benefits in aquaculture practice. However, little research has been reported on sex determination and gonadal development in red-tail catfish. Results: In this study, we performed the first transcriptomic analysis of male and female gonads at four developmental stages at 10, 18, 30, and 48 days post hatching (dph) using RNA-seq technology. A total of 23,588 genes were screened in 24 sequenced samples, of which 28, 213, 636, and 1381 differentially expressed genes (DEGs) were detected at four developmental stages, respectively. Seven candidate genes of sex determination and differentiation were further identified. Real-time quantitative PCR (RT-qPCR) further confirmed that anti-Mullerian hormone (amh), growth differentiation factor 6a (gdf6a), testis-specific gene antigen 10 (tsga10), and cytochrome P450 family 17 subfamily A (cyp17a) were highly expressed mainly in the male, while cytochrome P450 family 19 subfamily A polypeptide 1b (cyp19a1b), forkhead box L2 (foxl2), and hydroxysteroid 17-beta dehydrogenase 1 (hsd17b1) were highly expressed in the female. The KEGG pathway enrichment data showed that these identified DEGs were mainly involved in steroid hormone biosynthesis and TGF-β signaling pathways. Conclusions: Based on RNA-seq data of gonads at the early developmental stages, seven DEGs shared by the four developmental stages were identified, among which amh and gdf6a may be the male-biased expression genes, while foxl2, cyp19a1b and hsd17b1 may be the female-biased expression genes in red-tail catfish. Our study will provide crucial genetic information for the research on sex control in red-tail catfish, as well as for exploring the evolutionary processes of sex determination mechanisms in fish. [ABSTRACT FROM AUTHOR]

Abstract Background Sex determination mechanisms are complicated and diverse across taxonomic categories. Sea cucumber Apostichopus japonicus is a benthic echinoderm, which is the closest group of invertebrates to chordate, and important economic and ecologically aquaculture species in China. A. japonicus is dioecious, and no phenotypic differences between males and females can be detected before sexual maturation. Identification of sex determination locus will broaden knowledge about sex-determination mechanism in echinoderms, which allows for the identification of sex-linked markers and increases the efficiency of sea cucumber breeding industry. Results Here, we integrated assembly of a novel chromosome-level genome and resequencing of female and male populations to investigate the sex determination mechanisms of A. japonicus. We built a chromosome-level genome assembly AJH1.0 using Hi-C technology. The assembly AJH1.0 consists of 23 chromosomes ranging from 22.4 to 60.4 Mb. To identify the sex-determination locus of A. japonicus, we conducted genome-wide association study (GWAS) and analyses of distribution characteristics of sex-specific SNPs and fixation index FST. The GWAS analysis showed that multiple sex-associated loci were located on several chromosomes, including chromosome 4 (24.8%), followed by chromosome 9 (10.7%), chromosome 17 (10.4%), and chromosome 18 (14.1%). Furthermore, analyzing the homozygous and heterozygous genotypes of plenty of sex-specific SNPs in females and males confirmed that A. japonicus might have a XX/XY sex determination system. As a physical region of 10 Mb on chromosome 4 included the highest number of sex-specific SNPs and higher FST values, this region was considered as the candidate sex determination region (SDR) in A. japonicus. Conclusions In the present study, we integrated genome-wide association study and analyses of sex-specific variations to investigate sex determination mechanisms. This will bring novel insights into gene regulation during primitive gonadogenesis and differentiation and identification of master sex determination gene in sea cucumber. In the sea cucumber industry, investigation of molecular mechanisms of sex determination will be helpful for artificial fertilization and precise breeding.

Abstract Background In some vertebrate species, gene-environment interactions can determine sex, driving bipotential gonads to differentiate into either ovaries or testes. In the central bearded dragon (Pogona vitticeps), the genetic influence of sex chromosomes (ZZ/ZW) can be overridden by high incubation temperatures, causing ZZ male to female sex reversal. Previous research showed ovotestes, a rare gonadal phenotype with traits of both sexes, develop during sex reversal, leading to the hypothesis that sex reversal relies on high temperature feminisation to outcompete the male genetic cue. To test this, we conducted temperature switching experiments at key developmental stages, and analysed the effect on gonadal phenotypes using histology and transcriptomics. Results We found sexual fate is more strongly influenced by the ZZ genotype than temperature. Any exposure to low temperatures (28 °C) caused testes differentiation, whereas sex reversal required longer exposure to high temperatures. We revealed ovotestes exist along a spectrum of femaleness to male-ness at the transcriptional level. We found inter-individual variation in gene expression changes following temperature switches, suggesting both genetic sensitivity to, and the timing and duration of the temperature cue influences sex reversal. Conclusions These findings bring new insights to the mechanisms underlying sex reversal, improving our understanding of thermosensitive sex systems in vertebrates.

Abstract Background The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers. Results The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females). Conclusions The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

Abstract Background Tambaqui (Colossoma macropomum, Cuvier, 1818) is the most economically important native freshwater fish species in Brazil. It can reach a total length of over 1 m and a weight of over 40 kg. The species displays a clear sex dimorphism in growth performance, with females reaching larger sizes at harvest. In aquaculture, the production of monosex populations in selective breeding programmes has been therefore identified as a key priority. Results In the present study, a genetic linkage map was generated by double digest restriction-site associated DNA (ddRAD) sequencing from 248 individuals sampled from two F1 families. The map was constructed using 14,805 informative SNPs and spanned 27 linkage groups. From this, the tambaqui draft genome was improved, by ordering the scaffolds into chromosomes, and sex-linked markers were identified. A total of 235 markers on linkage group 26 showed a significant association with the phenotypic sex, supporting an XX/XY sex determination system in the species. The four most informative sex-linked markers were validated on another 206 sexed individuals, demonstrating an accuracy in predicting sex ranging from 90.0 to 96.7%. Conclusions The genetic mapping and novel sex-linked DNA markers identified and validated offer new tools for rapid progeny sexing, thus supporting the development of monosex female production in the industry while also supporting breeding programmes of the species.

Abstract Background Sex chromosomes are in some species largely undifferentiated (homomorphic) with restricted sex determination regions. Homomorphic but different sex chromosomes are found in the closely related genera Populus and Salix indicating flexible sex determination systems, ideal for studies of processes involved in sex chromosome evolution. We have performed genome-wide association studies of sex and analysed sex chromosomes in a population of 265 wild collected Salix viminalis accessions and studied the sex determining locus. Results A total of 19,592 markers were used in association analyses using both Fisher’s exact tests and a single-marker mixed linear model, which resulted in 48 and 41 sex-associated (SA) markers respectively. Across all 48 SA markers, females were much more often heterozygous than males, which is expected if females were the heterogametic sex. The majority of the SA markers were, based on positions in the S. purpurea genome, located on chromosome 15, previously demonstrated to be the sex chromosome. Interestingly, when mapping the genotyping-by-sequencing sequence tag harbouring the two SA markers with the highest significance to the S. viminalis genomic scaffolds, five regions of very high similarity were found: three on a scaffold that represents a part of chromosome 15, one on a scaffold that represents a part of chromosome 9 and one on a scaffold not anchored to the genome. Based on segregation differences of the alleles at the two marker positions and on differences in PCR amplification between females and males we conclude that females had multiple copies of this DNA fragment (chromosome 9 and 15), whereas males only had one (chromosome 9). We therefore postulate that the female specific sequences have been copied from chromosome 9 and inserted on chromosome 15, subsequently developing into a hemizygous W chromosome linked region. Conclusions Our results support that sex determination in S. viminalis is controlled by one locus on chromosome 15. The segregation patterns observed at the SA markers furthermore confirm that S. viminalis females are the heterogametic sex. We also identified a translocation from chromosome 9 to the W chromosome.

Background: In the life history of Daphnia, the reproductive mode of parthenogenesis and sexual reproduction alternate in aquatic ecosystem, which are often affected by environmental and genetic factors. Recently, the sex-biased genes are of great significance for clarifying the origin and evolution of reproductive transformation and the molecular regulation mechanism of sex determination in Daphnia. Although some genes on reproductive transition of Daphnia had been researched, molecular mechanism on the maintenance of sexually dimorphic phenotypes of Daphnia are still not well known, including differentially expressed genes in different life-history stages. Results: In this study, four life-history stages of Daphnia sinensis, juvenile female (JF), parthenogenetic female (PF), sexual female (SF) and male (M), were performed for transcriptome, and male-biased genes were screened. A total of 110437 transcripts were obtained and assembled into 22996 unigenes. In the four life-history stages (JF, PF, SF and M), the number of unique unigenes is respectively 2863, 445, 437 and 586, and the number of common unigenes is 9708. The differentially expressed genes (DEGs) between male and other three female stages (M vs JF, M vs PF and M vs SF) were 4570, 4358 and 2855, respectively. GO gene enrichment analysis showed that the up-regulated genes in male were mainly enriched in hydrolase activity and peptidase activity. Thirty-six genes in male were significantly higher expression than in the three female stages, including one Doublesex (Dsx) gene, one laminin gene, five trypsin genes and one serine protease genes, and one chitin synthase gene and two chitinase genes. Conclusions: Our results showed that thirty-six candidate genes may be as the male-biased genes involving in the maintenance of sexually dimorphic phenotypes. This work will provide a reference for further exploring the functional genes related to sex differentiation in Daphnia species. Moreover, according to previous investigations, we thought that the expression level of functional genes may be related to the life-history stages of organisms, and may be also affected by different Daphnia species. [ABSTRACT FROM AUTHOR]

Abstract Background In the animal kingdom, mollusca is an important phylum of the Lophotrochozoa. However, few studies have investigated the molecular cascade of sex determination/early gonadal differentiation within this phylum. The oyster Crassostrea gigas is a sequential irregular hermaphrodite mollusc of economic, physiological and phylogenetic importance. Although some studies identified genes of its sex-determining/−differentiating pathway, this particular topic remains to be further deepened, in particular with regard to the expression patterns. Indeed, these patterns need to cover the entire period of sex lability and have to be associated to future sex phenotypes, usually impossible to establish in this sequential hermaphrodite. This is why we performed a gonadal RNA-Seq analysis of diploid male and female oysters that have not changed sex for 4 years, sampled during the entire time-window of sex determination/early sex differentiation (stages 0 and 3 of the gametogenetic cycle). This individual long-term monitoring gave us the opportunity to explain the molecular expression patterns in the light of the most statistically likely future sex of each oyster. Results The differential gene expression analysis of gonadal transcriptomes revealed that 9723 genes were differentially expressed between gametogenetic stages, and 141 between sexes (98 and 43 genes highly expressed in females and males, respectively). Eighty-four genes were both stage- and sex-specific, 57 of them being highly expressed at the time of sex determination/early sex differentiation. These 4 novel genes including Trophoblast glycoprotein-like, Protein PML-like, Protein singed-like and PREDICTED: paramyosin, while being supported by RT-qPCR, displayed sexually dimorphic gene expression patterns. Conclusions This gonadal transcriptome analysis, the first one associated with sex phenotypes in C. gigas, revealed 57 genes highly expressed in stage 0 or 3 of gametogenesis and which could be linked to the future sex of the individuals. While further study will be needed to suggest a role for these factors, some could certainly be original potential actors involved in sex determination/early sex differentiation, like paramyosin and could be used to predict the future sex of oysters.

Background: Sex determination mechanisms are complicated and diverse across taxonomic categories. Sea cucumber Apostichopus japonicus is a benthic echinoderm, which is the closest group of invertebrates to chordate, and important economic and ecologically aquaculture species in China. A. japonicus is dioecious, and no phenotypic differences between males and females can be detected before sexual maturation. Identification of sex determination locus will broaden knowledge about sex-determination mechanism in echinoderms, which allows for the identification of sex-linked markers and increases the efficiency of sea cucumber breeding industry. Results: Here, we integrated assembly of a novel chromosome-level genome and resequencing of female and male populations to investigate the sex determination mechanisms of A. japonicus. We built a chromosome-level genome assembly AJH1.0 using Hi-C technology. The assembly AJH1.0 consists of 23 chromosomes ranging from 22.4 to 60.4 Mb. To identify the sex-determination locus of A. japonicus, we conducted genome-wide association study (GWAS) and analyses of distribution characteristics of sex-specific SNPs and fixation index FST. The GWAS analysis showed that multiple sex-associated loci were located on several chromosomes, including chromosome 4 (24.8%), followed by chromosome 9 (10.7%), chromosome 17 (10.4%), and chromosome 18 (14.1%). Furthermore, analyzing the homozygous and heterozygous genotypes of plenty of sex-specific SNPs in females and males confirmed that A. japonicus might have a XX/XY sex determination system. As a physical region of 10 Mb on chromosome 4 included the highest number of sex-specific SNPs and higher FST values, this region was considered as the candidate sex determination region (SDR) in A. japonicus. Conclusions: In the present study, we integrated genome-wide association study and analyses of sex-specific variations to investigate sex determination mechanisms. This will bring novel insights into gene regulation during primitive gonadogenesis and differentiation and identification of master sex determination gene in sea cucumber. In the sea cucumber industry, investigation of molecular mechanisms of sex determination will be helpful for artificial fertilization and precise breeding. [ABSTRACT FROM AUTHOR]

Background: Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development. Results: Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination. Conclusions: We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology. [ABSTRACT FROM AUTHOR]

Background: Aedes aegypti is a medically-important mosquito vector that transmits arboviruses including yellow fever, dengue, chikungunya, and Zika viruses to humans. The mosquito exhibits typical sexually dimorphic behaviors such as courtship, mating, host seeking, bloodfeeding, and oviposition. All these behaviors are mainly regulated by the brain; however, little is known about the function and neuron composition of the mosquito brain. In this study, we generated an initial atlas of the adult male and female brain of Ae. aegypti using 10xGenomics based single-nucleus RNA sequencing. Results: We identified 35 brain cell clusters in male and female brains, and 15 of those clusters were assigned to known cell types. Identified cell types include glia (astrocytes), Kenyon cells, (ventral) projection neurons, monoaminergic neurons, medulla neurons, and proximal medulla neurons. In addition, the cell type compositions of male and female brains were compared to each other showing that they were quantitatively distinct, as 17 out of 35 cell clusters varied significantly in their cell type proportions. Overall, the transcriptomes from each cell cluster looked very similar between the male and female brain as only up to 25 genes were differentially expressed in these clusters. The sex determination factor Nix was highly expressed in neurons and glia of the male brain, whereas doublesex (dsx) was expressed in all neuron and glia cell clusters of the male and female brain. Conclusions: An initial cell atlas of the brain of the mosquito Ae. aegypti has been generated showing that the cellular compositions of the male and female brains of this hematophagous insect differ significantly from each other. Although some of the rare brain cell types have not been detected in our single biological replicate, this study provides an important basis for the further development of a complete brain cell atlas as well as a better understanding of the neurobiology of the brains of male and female mosquitoes and their sexually dimorphic behaviors. [ABSTRACT FROM AUTHOR]

Abstract Background A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits. Results By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map’ length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation. Conclusions We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

Abstract Background Müllerian ducts are paired embryonic tubes that give rise to the female reproductive tract in vertebrates. Many disorders of female reproduction can be attributed to anomalies of Müllerian duct development. However, the molecular genetics of Müllerian duct formation is poorly understood and most disorders of duct development have unknown etiology. In this study, we describe for the first time the transcriptional landscape of the embryonic Müllerian duct, using the chicken embryo as a model system. RNA sequencing was conducted at 1 day intervals during duct formation to identify developmentally-regulated genes, validated by in situ hybridization. Results This analysis detected hundreds of genes specifically up-regulated during duct morphogenesis. Gene ontology and pathway analysis revealed enrichment for developmental pathways associated with cell adhesion, cell migration and proliferation, ERK and WNT signaling, and, interestingly, axonal guidance. The latter included factors linked to neuronal cell migration or axonal outgrowth, such as Ephrin B2, netrin receptor, SLIT1 and class A semaphorins. A number of transcriptional modules were identified that centred around key hub genes specifying matrix-associated signaling factors; SPOCK1, HTRA3 and ADGRD1. Several novel regulators of the WNT and TFG-β signaling pathway were identified in Müllerian ducts, including APCDD1 and DKK1, BMP3 and TGFBI. A number of novel transcription factors were also identified, including OSR1, FOXE1, PRICKLE1, TSHZ3 and SMARCA2. In addition, over 100 long non-coding RNAs (lncRNAs) were expressed during duct formation. Conclusions This study provides a rich resource of new candidate genes for Müllerian duct development and its disorders. It also sheds light on the molecular pathways engaged during tubulogenesis, a fundamental process in embryonic development.

Abstract Background Goldfish is an important model for various areas of research, including neural development and behavior and a species of significant importance in aquaculture, especially as an ornamental species. It has a male heterogametic (XX/XY) sex determination system that relies on both genetic and environmental factors, with high temperatures being able to produce female-to-male sex reversal. Little, however, is currently known on the molecular basis of genetic sex determination in this important cyprinid model. Here we used sequencing approaches to better characterize sex determination and sex-chromosomes in an experimental strain of goldfish. Results Our results confirmed that sex determination in goldfish is a mix of environmental and genetic factors and that its sex determination system is male heterogametic (XX/XY). Using reduced representation (RAD-seq) and whole genome (pool-seq) approaches, we characterized sex-linked polymorphisms and developed male specific genetic markers. These male specific markers were used to distinguish sex-reversed XX neomales from XY males and to demonstrate that XX female-to-male sex reversal could even occur at a relatively low rearing temperature (18 °C), for which sex reversal has been previously shown to be close to zero. We also characterized a relatively large non-recombining region (~ 11.7 Mb) on goldfish linkage group 22 (LG22) that contained a high-density of male-biased genetic polymorphisms. This large LG22 region harbors 373 genes, including a single candidate as a potential master sex gene, i.e., the anti-Mullerian hormone gene (amh). However, no sex-linked polymorphisms were detected in the coding DNA sequence of the goldfish amh gene. Conclusions These results show that our goldfish strain has a relatively large sex locus on LG22, which is likely the Y chromosome of this experimental population. The presence of a few XX males even at low temperature also suggests that other environmental factors in addition to temperature could trigger female-to-male sex reversal. Finally, we also developed sex-linked genetic markers, which will be important tools for future research on sex determination in our experimental goldfish population. However, additional work would be needed to explore whether this sex locus is conserved in other populations of goldfish.

Background: The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers. Results: The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females). Conclusions: The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations. [ABSTRACT FROM AUTHOR]

Background: Tambaqui (Colossoma macropomum, Cuvier, 1818) is the most economically important native freshwater fish species in Brazil. It can reach a total length of over 1 m and a weight of over 40 kg. The species displays a clear sex dimorphism in growth performance, with females reaching larger sizes at harvest. In aquaculture, the production of monosex populations in selective breeding programmes has been therefore identified as a key priority. Results: In the present study, a genetic linkage map was generated by double digest restriction-site associated DNA (ddRAD) sequencing from 248 individuals sampled from two F1 families. The map was constructed using 14,805 informative SNPs and spanned 27 linkage groups. From this, the tambaqui draft genome was improved, by ordering the scaffolds into chromosomes, and sex-linked markers were identified. A total of 235 markers on linkage group 26 showed a significant association with the phenotypic sex, supporting an XX/XY sex determination system in the species. The four most informative sex-linked markers were validated on another 206 sexed individuals, demonstrating an accuracy in predicting sex ranging from 90.0 to 96.7%. Conclusions: The genetic mapping and novel sex-linked DNA markers identified and validated offer new tools for rapid progeny sexing, thus supporting the development of monosex female production in the industry while also supporting breeding programmes of the species. [ABSTRACT FROM AUTHOR]

Background: Sex chromosomes are in some species largely undifferentiated (homomorphic) with restricted sex determination regions. Homomorphic but different sex chromosomes are found in the closely related genera Populus and Salix indicating flexible sex determination systems, ideal for studies of processes involved in sex chromosome evolution. We have performed genome-wide association studies of sex and analysed sex chromosomes in a population of 265 wild collected Salix viminalis accessions and studied the sex determining locus. Results: A total of 19,592 markers were used in association analyses using both Fisher's exact tests and a single-marker mixed linear model, which resulted in 48 and 41 sex-associated (SA) markers respectively. Across all 48 SA markers, females were much more often heterozygous than males, which is expected if females were the heterogametic sex. The majority of the SA markers were, based on positions in the S. purpurea genome, located on chromosome 15, previously demonstrated to be the sex chromosome. Interestingly, when mapping the genotyping-by-sequencing sequence tag harbouring the two SA markers with the highest significance to the S. viminalis genomic scaffolds, five regions of very high similarity were found: three on a scaffold that represents a part of chromosome 15, one on a scaffold that represents a part of chromosome 9 and one on a scaffold not anchored to the genome. Based on segregation differences of the alleles at the two marker positions and on differences in PCR amplification between females and males we conclude that females had multiple copies of this DNA fragment (chromosome 9 and 15), whereas males only had one (chromosome 9). We therefore postulate that the female specific sequences have been copied from chromosome 9 and inserted on chromosome 15, subsequently developing into a hemizygous W chromosome linked region. Conclusions: Our results support that sex determination in S. viminalis is controlled by one locus on chromosome 15. The segregation patterns observed at the SA markers furthermore confirm that S. viminalis females are the heterogametic sex. We also identified a translocation from chromosome 9 to the W chromosome. [ABSTRACT FROM AUTHOR]

Background: The lower Dipteran fungus fly, Sciara coprophila, has many unique biological features that challenge the rule of genome DNA constancy. For example, Sciara undergoes paternal chromosome elimination and maternal X chromosome nondisjunction during spermatogenesis, paternal X elimination during embryogenesis, intrachromosomal DNA amplification of DNA puff loci during larval development, and germline-limited chromosome elimination from all somatic cells. Paternal chromosome elimination in Sciara was the first observation of imprinting, though the mechanism remains a mystery. Here, we present the first draft genome sequence for Sciara coprophila to take a large step forward in addressing these features. Results: We assembled the Sciara genome using PacBio, Nanopore, and Illumina sequencing. To find an optimal assembly using these datasets, we generated 44 short-read and 50 long-read assemblies. We ranked assemblies using 27 metrics assessing contiguity, gene content, and dataset concordance. The highest-ranking assemblies were scaffolded using BioNano optical maps. RNA-seq datasets from multiple life stages and both sexes facilitated genome annotation. A set of 66 metrics was used to select the first draft assembly for Sciara. Nearly half of the Sciara genome sequence was anchored into chromosomes, and all scaffolds were classified as X-linked or autosomal by coverage. Conclusions: We determined that X-linked genes in Sciara males undergo dosage compensation. An entire bacterial genome from the Rickettsia genus, a group known to be endosymbionts in insects, was co-assembled with the Sciara genome, opening the possibility that Rickettsia may function in sex determination in Sciara. Finally, the signal level of the PacBio and Nanopore data support the presence of cytosine and adenine modifications in the Sciara genome, consistent with a possible role in imprinting. [ABSTRACT FROM AUTHOR]

Background: In the animal kingdom, mollusca is an important phylum of the Lophotrochozoa. However, few studies have investigated the molecular cascade of sex determination/early gonadal differentiation within this phylum. The oyster Crassostrea gigas is a sequential irregular hermaphrodite mollusc of economic, physiological and phylogenetic importance. Although some studies identified genes of its sex-determining/−differentiating pathway, this particular topic remains to be further deepened, in particular with regard to the expression patterns. Indeed, these patterns need to cover the entire period of sex lability and have to be associated to future sex phenotypes, usually impossible to establish in this sequential hermaphrodite. This is why we performed a gonadal RNA-Seq analysis of diploid male and female oysters that have not changed sex for 4 years, sampled during the entire time-window of sex determination/early sex differentiation (stages 0 and 3 of the gametogenetic cycle). This individual long-term monitoring gave us the opportunity to explain the molecular expression patterns in the light of the most statistically likely future sex of each oyster. Results: The differential gene expression analysis of gonadal transcriptomes revealed that 9723 genes were differentially expressed between gametogenetic stages, and 141 between sexes (98 and 43 genes highly expressed in females and males, respectively). Eighty-four genes were both stage- and sex-specific, 57 of them being highly expressed at the time of sex determination/early sex differentiation. These 4 novel genes including Trophoblast glycoprotein-like, Protein PML-like, Protein singed-like and PREDICTED: paramyosin, while being supported by RT-qPCR, displayed sexually dimorphic gene expression patterns. Conclusions: This gonadal transcriptome analysis, the first one associated with sex phenotypes in C. gigas, revealed 57 genes highly expressed in stage 0 or 3 of gametogenesis and which could be linked to the future sex of the individuals. While further study will be needed to suggest a role for these factors, some could certainly be original potential actors involved in sex determination/early sex differentiation, like paramyosin and could be used to predict the future sex of oysters. [ABSTRACT FROM AUTHOR]

Background: The clearhead icefish, Protosalanx hyalocranius, is an economically important fishery species in China. Since 1980s, P. hyalocranius was widely introduced into lakes and reservoirs of northern China for aquaculture. However, the lack of a rapid and cost-effective sex identification method based on sex specific genetic markers has hindered study on sex determination mechanisms and breeding applications. Results: Female-specific genomic regions were discovered by comparing whole genome re-sequencing data of both males and females. Two female-specific genomic regions larger than 50 bp were identified, and one (598 bp) contained a putative FOXI gene, which was paralogous to another FOXI gene with sex-associated SNPs. The two FOXI sequences displayed significant length difference with nine deletions of total length of 230 bp. This deletion-type structural variation could be easily and efficiently detected by traditional PCR and agarose gel electrophoresis with one 569 bp band for males and two bands (569 and 339 bp) for females, which were validated in 50 females and 40 males with known phenotypic sexes. Conclusions: The results provided structural genomic evidence for the ZZ/ZW sex determination system in P. hyalocranius discovered in our previous study with association analysis of SNPs. Moreover, the female-specific markers and rapid and cost-effective PCR-based genetic sex identification method should have applications in further studies of sex determination mechanism for this species. [ABSTRACT FROM AUTHOR]

Abstract Background A complete understanding of the genetic basis for sexual determination and differentiation is necessary in order to implement efficient breeding schemes at early stages of development. Atlantic salmon belongs to the family Salmonidae of fishes and represents a species of great commercial value. Although the species is assumed to be male heterogametic with XY sex determination, the precise genetic basis of sexual development remains unclear. The complexity is likely associated to the relatively recent salmonid specific whole genome duplication that may be responsible for certain genome instability. This instability together with the capacity of the sex-determining gene to move across the genome as reported by previous studies, may explain that sexual development genes are not circumscribed to the same chromosomes in all members of the species. In this study, we have used a 220 K SNP panel developed for Atlantic salmon to identify the chromosomes explaining the highest proportion of the genetic variance for sex as well as candidate regions and genes associated to sexual development in this species. Results Results from regional heritability analysis showed that the chromosomes explaining the highest proportion of variance in these populations were Ssa02 (heritability = 0.42, SE = 0.12) and Ssa21 (heritability = 0.26, SE = 0.11). After pruning by linkage disequilibrium, genome-wide association analyses revealed 114 SNPs that were significantly associated with sex, being Ssa02 the chromosome containing a greatest number of regions. Close examination of the candidate regions evidenced important genes related to sex in other species of Class Actinopterygii, including SDY, genes from family SOX, RSPO1, ESR1, U2AF2A, LMO7, GNRH-R, DND and FIGLA. Conclusions The combined results from regional heritability analysis and genome-wide association have provided new advances in the knowledge of the genetic regulation of sex determination in Atlantic salmon, supporting that Ssa02 is the candidate chromosome for sex in this species and suggesting an alternative population lineage in Spanish wild populations according to the results from Ssa21.

Abstract Background Phlebotomine sand flies (Diptera, Nematocera) are important vectors of several pathogens, including Leishmania parasites, causing serious diseases of humans and dogs. Despite their importance as disease vectors, most aspects of sand fly biology remain unknown including the molecular basis of their reproduction and sex determination, aspects also relevant for the development of novel vector control strategies. Results Using comparative genomics/transcriptomics data mining and transcriptional profiling, we identified the sex determining genes in phlebotomine sand flies and proposed the first model for the sex determination cascade of these insects. For all the genes identified, we produced manually curated gene models, developmental gene expression profile and performed evolutionary molecular analysis. We identified and characterized, for the first time in a Nematocera species, the transformer (tra) homolog which exhibits both conserved and novel features. The analysis of the tra locus in sand flies and its expression pattern suggest that this gene is able to autoregulate its own splicing, as observed in the fruit fly Ceratitis capitata and several other insect species. Conclusions Our results permit to fill the gap about sex determination in sand flies, contribute to a better understanding of this developmental pathway in Nematocera and open the way for the identification of sex determining orthologs in other species of this important Diptera sub-order. Furthermore, the sex determination genes identified in our work also provide the opportunity of future biotechnological applications to control natural population of sand flies, reducing their impact on public health.

Abstract Background Elevated water temperature, as is expected through climate change, leads to masculinization in fish species with sexual plasticity, resulting in changes in population dynamics. These changes are one important ecological consequence, contributing to the risk of extinction in small and inbred fish populations under natural conditions, due to male-biased sex ratio. Here we investigated the effect of elevated water temperature during embryogenesis on sex ratio and sex-biased gene expression profiles between two different tissues, namely gonad and caudal fin of adult zebrafish males and females, to gain new insights into the molecular mechanisms underlying sex determination (SD) and colour patterning related to sexual attractiveness. Results Our study demonstrated sex ratio imbalances with 25.5% more males under high-temperature condition, resulting from gonadal masculinization. The result of transcriptome analysis showed a significantly upregulated expression of male SD genes (e.g. dmrt1, amh, cyp11c1 and sept8b) and downregulation of female SD genes (e.g. zp2.1, vtg1, cyp19a1a and bmp15) in male gonads compared to female gonads. Contrary to expectations, we found highly differential expression of colour pattern (CP) genes in the gonads, suggesting the ‘neofunctionalisation’ of those genes in the zebrafish reproduction system. However, in the caudal fin, no differential expression of CP genes was identified, suggesting the observed differences in colouration between males and females in adult fish may be due to post-transcriptional regulation of key enzymes involved in pigment synthesis and distribution. Conclusions Our study demonstrates male-biased sex ratio under high temperature condition and support a polygenic SD (PSD) system in laboratory zebrafish. We identify a subset of pathways (tight junction, gap junction and apoptosis), enriched for SD and CP genes, which appear to be co-regulated in the same pathway, providing evidence for involvement of those genes in the regulation of phenotypic sexual dimorphism in zebrafish.

Abstract Background Alternative splicing (AS) may play an important role in gonadal sex determination (GSD) in mammals. The present study was designed to identify differentially expressed isoforms and AS modifications accompanying GSD in mice. Results Using deep RNA-sequencing, we performed a transcriptional analysis of XX and XY gonads during sex determination on embryonic days 11 (E11) and 12 (E12). Analysis of differentially expressed genes (DEG) identified hundreds of genes related to GSD and early sex differentiation that may represent good candidates for sex reversal. Expression at time point E11 in males was significantly enriched in RNA splicing and mRNA processing Gene Ontology terms. Differentially expressed isoform analysis identified hundreds of specific isoforms related to GSD, many of which showed no differences in the DEG analysis. Hundreds of AS events were identified as modified at E11 and E12. Female E11 gonads featured sex-biased upregulation of intron retention (in genes related to regulation of transcription, protein phosphorylation, protein transport and mRNA splicing) and exon skipping (in genes related to chromatin repression) suggesting AS as a post-transcription mechanism that controls sex determination of the bipotential fetal gonad. Conclusion Our data suggests an important role of splicing regulatory mechanisms for sex determination in mice.

Background: Müllerian ducts are paired embryonic tubes that give rise to the female reproductive tract in vertebrates. Many disorders of female reproduction can be attributed to anomalies of Müllerian duct development. However, the molecular genetics of Müllerian duct formation is poorly understood and most disorders of duct development have unknown etiology. In this study, we describe for the first time the transcriptional landscape of the embryonic Müllerian duct, using the chicken embryo as a model system. RNA sequencing was conducted at 1 day intervals during duct formation to identify developmentally-regulated genes, validated by in situ hybridization. Results: This analysis detected hundreds of genes specifically up-regulated during duct morphogenesis. Gene ontology and pathway analysis revealed enrichment for developmental pathways associated with cell adhesion, cell migration and proliferation, ERK and WNT signaling, and, interestingly, axonal guidance. The latter included factors linked to neuronal cell migration or axonal outgrowth, such as Ephrin B2, netrin receptor, SLIT1 and class A semaphorins. A number of transcriptional modules were identified that centred around key hub genes specifying matrix-associated signaling factors; SPOCK1, HTRA3 and ADGRD1. Several novel regulators of the WNT and TFG-β signaling pathway were identified in Müllerian ducts, including APCDD1 and DKK1, BMP3 and TGFBI. A number of novel transcription factors were also identified, including OSR1, FOXE1, PRICKLE1, TSHZ3 and SMARCA2. In addition, over 100 long non-coding RNAs (lncRNAs) were expressed during duct formation. Conclusions: This study provides a rich resource of new candidate genes for Müllerian duct development and its disorders. It also sheds light on the molecular pathways engaged during tubulogenesis, a fundamental process in embryonic development. [ABSTRACT FROM AUTHOR]

Background: Goldfish is an important model for various areas of research, including neural development and behavior and a species of significant importance in aquaculture, especially as an ornamental species. It has a male heterogametic (XX/XY) sex determination system that relies on both genetic and environmental factors, with high temperatures being able to produce female-to-male sex reversal. Little, however, is currently known on the molecular basis of genetic sex determination in this important cyprinid model. Here we used sequencing approaches to better characterize sex determination and sex-chromosomes in an experimental strain of goldfish. Results: Our results confirmed that sex determination in goldfish is a mix of environmental and genetic factors and that its sex determination system is male heterogametic (XX/XY). Using reduced representation (RAD-seq) and whole genome (pool-seq) approaches, we characterized sex-linked polymorphisms and developed male specific genetic markers. These male specific markers were used to distinguish sex-reversed XX neomales from XY males and to demonstrate that XX female-to-male sex reversal could even occur at a relatively low rearing temperature (18 °C), for which sex reversal has been previously shown to be close to zero. We also characterized a relatively large non-recombining region (~ 11.7 Mb) on goldfish linkage group 22 (LG22) that contained a high-density of male-biased genetic polymorphisms. This large LG22 region harbors 373 genes, including a single candidate as a potential master sex gene, i.e., the anti-Mullerian hormone gene (amh). However, no sex-linked polymorphisms were detected in the coding DNA sequence of the goldfish amh gene. Conclusions: These results show that our goldfish strain has a relatively large sex locus on LG22, which is likely the Y chromosome of this experimental population. The presence of a few XX males even at low temperature also suggests that other environmental factors in addition to temperature could trigger female-to-male sex reversal. Finally, we also developed sex-linked genetic markers, which will be important tools for future research on sex determination in our experimental goldfish population. However, additional work would be needed to explore whether this sex locus is conserved in other populations of goldfish. [ABSTRACT FROM AUTHOR]

Abstract Background Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex FST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings. Results Intersex standardized variance (FST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with FST values ≥ 0.100 and/or with increased expected heterozygosity (H exp , one-sided T-test, p

Abstract Background The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper. Results A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data. Conclusions The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Abstract Background The mandarin fish (Siniperca chuatsi) is an important and widely cultured fish in China. However, the lack of selective breeding of mandarin fish in previous decades has resulted in a decline in the growth rate of pond-cultured fish, a shortened period of sexual maturity, and reduced disease resistance; these issues seriously affect the quality and safety of the fish products. Therefore, it is necessary to establish a selective breeding program for the mandarin fish to improve the economical traits of the fish and to sustain the development of the mandarin fish industry. Results We constructed a high-density linkage map for it based on double digest restriction site associated DNA sequencing (ddRAD-Sequencing). This map contained 3283 dimorphic single nucleotide polymorphism markers and 24 linkage groups (LGs). The total map-length was 1972.01 cM, with an average interlocus distance of 0.61 cM. One significant quantitative trait locus (QTL) for sex determination trait was detected on LG23, which was supported by five markers, clustered between 60.27 and 68.71 cM. The highest logarithm of odds value (17.73) was located at 60.27 cM, near the marker r1_73194, accounting for 53.3% of the phenotypic variance. Genotypes of all the male fish on r1_33008 were homozygous, whereas those of all females were heterozygous. Thus, LG23 was considered a sex-related linkage group. Eleven significant QTLs, for three growth traits, at two growth stages and the increased values were distributed on four LGs; their contributions to the phenotypic variation were quite low (12.4–17.2%), suggesting that multiple genes affected the growth traits. Conclusion This high-resolution genetic map provides a valuable resource for fine-mapping of important traits and for identification of sex-related markers that should facilitate breeding of all-female mandarin fish for aquaculture and mechanistic studies on sex determination.

Abstract Background Tilapias are the second most farmed fishes in the world and a sustainable source of food. Like many other fish, tilapias are sexually dimorphic and sex is a commercially important trait in these fish. In this study, we developed a significantly improved assembly of the tilapia genome using the latest genome sequencing methods and show how it improves the characterization of two sex determination regions in two tilapia species. Results A homozygous clonal XX female Nile tilapia (Oreochromis niloticus) was sequenced to 44X coverage using Pacific Biosciences (PacBio) SMRT sequencing. Dozens of candidate de novo assemblies were generated and an optimal assembly (contig NG50 of 3.3Mbp) was selected using principal component analysis of likelihood scores calculated from several paired-end sequencing libraries. Comparison of the new assembly to the previous O. niloticus genome assembly reveals that recently duplicated portions of the genome are now well represented. The overall number of genes in the new assembly increased by 27.3%, including a 67% increase in pseudogenes. The new tilapia genome assembly correctly represents two recent vasa gene duplication events that have been verified with BAC sequencing. At total of 146Mbp of additional transposable element sequence are now assembled, a large proportion of which are recent insertions. Large centromeric satellite repeats are assembled and annotated in cichlid fish for the first time. Finally, the new assembly identifies the long-range structure of both a ~9Mbp XY sex determination region on LG1 in O. niloticus, and a ~50Mbp WZ sex determination region on LG3 in the related species O. aureus. Conclusions This study highlights the use of long read sequencing to correctly assemble recent duplications and to characterize repeat-filled regions of the genome. The study serves as an example of the need for high quality genome assemblies and provides a framework for identifying sex determining genes in tilapia and related fish species.

The article discusses a study that examined the key components of the sex determining cascade, identifying the first transformer homolog in a Nematocera species. Topics include the in silico translation of the five cDNAs, the use of Ppesod gene as endogenous positive control, and sand fly life cycle and developmental expression analyses of sex determining genes in Phlebotomus perniciosus.

Background: Alternative splicing (AS) may play an important role in gonadal sex determination (GSD) in mammals. The present study was designed to identify differentially expressed isoforms and AS modifications accompanying GSD in mice. Results: Using deep RNA-sequencing, we performed a transcriptional analysis of XX and XY gonads during sex determination on embryonic days 11 (E11) and 12 (E12). Analysis of differentially expressed genes (DEG) identified hundreds of genes related to GSD and early sex differentiation that may represent good candidates for sex reversal. Expression at time point E11 in males was significantly enriched in RNA splicing and mRNA processing Gene Ontology terms. Differentially expressed isoform analysis identified hundreds of specific isoforms related to GSD, many of which showed no differences in the DEG analysis. Hundreds of AS events were identified as modified at E11 and E12. Female E11 gonads featured sexbiased upregulation of intron retention (in genes related to regulation of transcription, protein phosphorylation, protein transport and mRNA splicing) and exon skipping (in genes related to chromatin repression) suggesting AS as a posttranscription mechanism that controls sex determination of the bipotential fetal gonad. Conclusion: Our data suggests an important role of splicing regulatory mechanisms for sex determination in mice. [ABSTRACT FROM AUTHOR]

Background: The factors determining sex in teleosts are diverse. Great efforts have been made to characterize the underlying genetic network in various species. However, only seven master sex-determining genes have been identified in teleosts. While the function of a few genes involved in sex determination and differentiation has been studied, we are far from fully understanding how genes interact to coordinate in this process. Results: To enable systematic insights into fish sexual differentiation, we generated a dynamic co-expression network from tilapia gonadal transcriptomes at 5, 20, 30, 40, 90, and 180 dah (days after hatching), plus 45 and 90 dat (days after treatment) and linked gene expression profiles to both development and sexual differentiation. Transcriptomic profiles of female and male gonads at 5 and 20 dah exhibited high similarities except for a small number of genes that were involved in sex determination, while drastic changes were observed from 90 to 180 dah, with a group of differently expressed genes which were involved in gonadal differentiation and gametogenesis. Weighted gene correlation network analysis identified changes in the expression of Borealin, Gtsf1, tesk1, Zar1, Cdn15, and Rpl that were correlated with the expression of genes previously known to be involved in sex differentiation, such as Foxl2, Cyp19a1a, Gsdf, Dmrt1, and Amh. Conclusions: Global gonadal gene expression kinetics during sex determination and differentiation have been extensively profiled in tilapia. These findings provide insights into the genetic framework underlying sex determination and sexual differentiation, and expand our current understanding of developmental pathways during teleost sex determination. [ABSTRACT FROM AUTHOR]

Background: The assembly and annotation of a genome is a valuable resource for a species, with applications ranging from conservation genomics to gene discovery. Genomic resource development is especially important for species in culture, such as the California Yellowtail (Seriola dorsalis), the likely candidate for the establishment of commercial offshore aquaculture production in southern California. Genomic resource development for this species will improve the understanding of sex and other phenotypic traits, and allow for rapid increases in genetic improvement for and economic gain in culture production. Results: We describe the assembly and annotation of the S. dorsalis genome, and present resequencing data from 45 male and 45 female wild-caught S. dorsalis used to identify a sex-determining region and marker in this species. The genome assembly captured approximately 93% of the total 685 MB genome with an average coverage depth of 180×. Using the assembled genome, resequencing data from the 90 fish were aligned to place boundaries on the sex-determining region. Sex-specific markers were developed based on a female-specific, 61 nucleotide deletion identified in that region. We hypothesize that Estradiol 17-beta-dehydrogenase is the putative sex-determining gene and propose a plausible genetic mechanism for ZW sex determination in S. dorsalis involving a female-specific deletion of a transcription factor binding motif that may be targeted by Sox3. Conclusions: Understanding the mechanism of sex determination and development of assays to determine sex is critical both for management of wild fisheries and for development of efficient and sustainable aquaculture practices. In addition, this genome assembly for S. dorsalis will be a substantial resource for a variety of future research applications. [ABSTRACT FROM AUTHOR]

Background: Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex FST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings. Results: Intersex standardized variance (FST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with FST values ≤ 0.100 and/or with increased expected heterozygosity (Hexp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex FST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex FST genes proximal to the M locus, whereas SenAae had inter-sex FST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased Hexp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation. Conclusions: The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae. [ABSTRACT FROM AUTHOR]

Background: The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. Results: Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180x. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. Conclusions: The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda. [ABSTRACT FROM AUTHOR]

Background: Tilapias are the second most farmed fishes in the world and a sustainable source of food. Like many other fish, tilapias are sexually dimorphic and sex is a commercially important trait in these fish. In this study, we developed a significantly improved assembly of the tilapia genome using the latest genome sequencing methods and show how it improves the characterization of two sex determination regions in two tilapia species. Results: A homozygous clonal XX female Nile tilapia (Oreochromis niloticus) was sequenced to 44X coverage using Pacific Biosciences (PacBio) SMRT sequencing. Dozens of candidate de novo assemblies were generated and an optimal assembly (contig NG50 of 3.3Mbp) was selected using principal component analysis of likelihood scores calculated from several paired-end sequencing libraries. Comparison of the new assembly to the previous O. niloticus genome assembly reveals that recently duplicated portions of the genome are now well represented. The overall number of genes in the new assembly increased by 27.3%, including a 67% increase in pseudogenes. The new tilapia genome assembly correctly represents two recent vasa gene duplication events that have been verified with BAC sequencing. At total of 146Mbp of additional transposable element sequence are now assembled, a large proportion of which are recent insertions. Large centromeric satellite repeats are assembled and annotated in cichlid fish for the first time. Finally, the new assembly identifies the long-range structure of both a ~9Mbp XY sex determination region on LG1 in O. niloticus, and a ~50Mbp WZ sex determination region on LG3 in the related species O. aureus. Conclusions: This study highlights the use of long read sequencing to correctly assemble recent duplications and to characterize repeat-filled regions of the genome. The study serves as an example of the need for high quality genome assemblies and provides a framework for identifying sex determining genes in tilapia and related fish species. [ABSTRACT FROM AUTHOR]

Background: Biased sex ratios are common among dioecious plant species despite the theoretical prediction of selective advantage of even sex ratios. Albeit the high prevalence of deviations from even sex ratios, the genetic causes to sex biases are rarely known outside of a few model species. Here we present a mechanism underlying the female biased sex ratio in the dioecious willow species Salix viminalis. Results: We compared the segregation pattern of genome-wide single nucleotide polymorphism markers in two contrasting bi-parental pedigree populations, the S3 with even sex ratio and the S5 with a female biased sex ratio. With the segregation analysis and comparison between the two populations, we were able to demonstrate that sex determination and sex ratio distortion are controlled by different genetic mechanisms. We furthermore located the sex ratio distorter locus to a Z/W-gametologous region on chromosome 15, which was in close linkage with the sex determination locus. Interestingly, all males in the population with biased sex ratio have in this sex ratio distorter locus the same genotype, meaning that males with the Z1/Z3-genotype were missing from the population, thereby creating the 2:1 female biased sex ratio. Conclusions: We attribute the absence of Z1/Z3 males to an allelic incompatibility between maternally and paternally inherited alleles in this sex ratio distorter locus. Due to the tight linkage with the sex determination locus only male individuals are purged from the population at an early age, presumably before or during seed development. We showed that such allelic incompatibility could be stably maintained over evolutionary times through a system of overdominant or pseudooverdominant alleles. Thus, it is possible that the same mechanism generates the female biased sex ratio in natural willow populations. [ABSTRACT FROM AUTHOR]