Multiple sex-associated regions and a putative sex chromosome in zebrafish revealed by RAD mapping and population genomics

Rbpms2 promotes female fate upstream of the nutrient sensing Gator2 complex component Mios

Reproductive success relies on proper establishment and maintenance of biological sex. In many animals, including mammals, the primary gonad is initially ovary biased. We previously showed the RNA binding protein (RNAbp), Rbpms2, is required for ovary fate in zebrafish. Here, we identified Rbpms2 targets in oocytes (Rbpms2-bound oocyte RNAs; rboRNAs). We identify Rbpms2 as a translational regulator of rboRNAs, which include testis factors and ribosome biogenesis factors. Further, genetic analyses indicate that Rbpms2 promotes nucleolar amplification via the mTorc1 signaling pathway, specifically through the mTorc1-activating Gap activity towards Rags 2 (Gator2) component, Missing oocyte (Mios). Cumulatively, our findings indicate that early gonocytes are in a dual poised, bipotential state in which Rbpms2 acts as a binary fate-switch. Specifically, Rbpms2 represses testis factors and promotes oocyte factors to promote oocyte progression through an essential Gator2-mediated checkpoint, thereby integrating regulation of sexual differentiation factors and nutritional availability pathways in zebrafish oogenesis.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Gonadal transcriptome analysis of sex-biased gene and genome-wide investigation of dmrt gene family in Acanthogobius ommaturus

Acanthogobius ommaturus is one of the largest goby fish, and widely distributed in the Northwestern Pacific as an annual benthic fish. This study aims to report the gonadal transcriptome of A. ommaturus and identify differentially expressed genes (DEGs) between sexes. A total of 5460 (27.94 %) DEGs were detected from genome, with 3301 (16.89 %) biased towards males and 2159 (11.05 %) towards females. Analysis of 76 known vertebrate sex-related genes revealed multiple key genes, including the male-biased genes dmrt1 (Doublesex and Mab-3 related transcription factor 1) and amh (Anti-Mullerian Hormone), and the female-biased genes foxl2 (Forkhead Box L2) and cyp19a1a (Cytochrome P450 Aromatase 19 Subfamily A1). Furthermore, a genome-wide gene family analysis focused on the most significantly differentially expressed male-biased gene, dmrt1, was conducted using the chromosomal-level genome. Six Aodmrt genes were identified and subjected to phylogenetic and protein interaction network analyses. To validate the expression pattern, quantitative real-time PCR (qRT-PCR) was performed and compared with gonadal transcriptome data. The results showed that only dmrt1 exhibited significant male-bias, while the expression levels and sex differences of other dmrt genes in the gonads were inconclusive. Interestingly, the other dmrt genes displayed higher expression levels in other tissues, suggesting currently unknown functions. In conclusion, this study provides valuable genetic information contributing to the understanding of the sex determination mechanism of A. ommaturus and bony fish.

The emerging H3K9me3 chromatin landscape during zebrafish embryogenesis

The structural organization of eukaryotic genomes is contingent upon the fractionation of DNA into transcriptionally permissive euchromatin and repressive heterochromatin. However, we have a limited understanding of how these distinct states are first established during animal embryogenesis. Histone 3 lysine 9 trimethylation (H3K9me3) is critical to heterochromatin formation and bulk establishment of this mark is thought to help drive large-scale remodeling of an initially naive chromatin state during animal embryogenesis. However, a detailed understanding of this process is lacking. Here, we leverage CUT&RUN to define the emerging H3K9me3 landscape of the zebrafish embryo with high sensitivity and temporal resolution. Despite the prevalence of DNA transposons in the zebrafish genome, we found that LTR transposons are preferentially targeted for embryonic H3K9me3 deposition, with different families exhibiting distinct establishment timelines. High signal-to-noise ratios afforded by CUT&RUN revealed new, emerging sites of low-amplitude H3K9me3 that initiated before the major wave of zygotic genome activation (ZGA). Early sites of establishment predominated at specific subsets of transposons and were particularly enriched for transposon sequences with maternal piRNAs and pericentromeric localization. Notably, the number of H3K9me3 enriched sites increased linearly across blastula development, while quantitative comparison revealed a >10-fold genome-wide increase in H3K9me3 signal at established sites over just 30 minutes at the onset of ZGA. Continued maturation of the H3K9me3 landscape was observed beyond the initial wave of bulk establishment.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Macrophage activation drives ovarian failure and masculinization in zebrafish

BMP15 is a conserved regulator of ovarian development and maintenance in vertebrates. In humans, premature ovarian insufficiency is caused by autoimmunity and genetic factors, including mutation of BMP15. The cellular mechanisms underlying ovarian failure caused by BMP15 mutation and immune contributions are not understood. Using zebrafish, we established a causal link between macrophage activation and ovarian failure, which, in zebrafish, causes sex reversal. We define a germline-soma signaling axis that activates macrophages and drives ovarian failure and female-to-male sex reversal. Germline loss of zebrafish Bmp15 impairs oogenesis and initiates this cascade. Single-cell RNA sequencing and genetic analyses implicate ovarian somatic cells that express conserved macrophage-activating ligands as mediators of ovarian failure and sex reversal. Genetic ablation of macrophages or elimination of Csf1Rb ligands, Il34 or Csf1a, delays or blocks premature oocyte loss and sex reversal. The axis identified here provides insight into the cells and pathways governing oocyte and ovary maintenance and potential therapeutic targets to preserve female fertility.

Age-Dependent Growth-Related QTL Variations in Pacific Abalone, Haliotis discus hannai

Pacific abalone is a high-value, commercially important marine invertebrate. It shows low growth as well as individual and yearly growth variation in aquaculture. Marker-assisted selection breeding could potentially resolve the problem of low and variable growth and increase genetic gain. Expression of quantitative trait loci (QTLs) for growth-related traits, viz., body weight, shell length, and shell width were analyzed at the first, second, and third year of age using an F1 cross population. A total of 37 chromosome-wide QTLs were identified in linkage groups 01, 02, 03, 04, 06, 07, 08, 10, 11, 12, and 13 at different ages. None of the QTLs detected at any one age were expressed in all three age groups. This result suggests that growth-related traits at different ages are influenced by different QTLs in each year. However, multiple-trait QTLs (where one QTL affects all three traits) were detected each year that are also age-specific. Eleven multiple-trait QTLs were detected at different ages: two QTLs in the first year; two QTLs in the second year; and seven QTLs in the third year. As abalone hatcheries use three-year-old abalone for breeding, QTL-linked markers that were detected at the third year of age could potentially be used in marker-assisted selection breeding programs.

Limited evidence for extensive genetic differentiation between X and Y chromosomes in Hybognathus amarus (Cypriniformes: Leuciscidae)

Sex determination systems and genetic sex differentiation across fishes are highly diverse but are unknown for most Cypriniformes, including Rio Grande silvery minnow (Hybognathus amarus). In this study, we aimed to detect and validate sex-linked markers to infer sex determination system and to demonstrate the utility of combining several methods for sex-linked marker detection in nonmodel organisms. To identify potential sex-linked markers, Nextera-tagmented reductively amplified DNA (nextRAD) libraries were generated from 66 females, 64 males, and 60 larvae of unknown sex. These data were combined with female and male de novo genomes from Nanopore long-read sequences. We identified five potential unique male nextRAD-tags and one potential unique male contig, suggesting an XY sex determination system. We also identified two single-nucleotide polymorphisms (SNPs) in the same contig with values of FST, allele frequencies, and heterozygosity conforming with expectations of an XY system. Through PCR we validated the marker containing the sex-linked SNPs and a single nextRAD-tag sex-associated marker but it was not male specific. Instead, more copies of this locus in the male genome were suggested by enhanced amplification in males. Results are consistent with an XY system with low differentiation between sex-determining regions. Further research is needed to confirm the level of differentiation between the sex chromosomes. Nonetheless, this study highlighted the power of combining reduced representation and whole-genome sequencing for identifying sex-linked markers, especially when reduced representation sequencing does not include extensive variation between sexes, either because such variation is not present or not captured.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

No evidence for sex chromosomes in natural populations of the cichlid fish Astatotilapia burtoni

Sex determination (SD) is not conserved among teleost fishes and can even differ between populations of the same species. Across the outstandingly species-rich fish family Cichlidae, more and more SD systems are being discovered. Still, the picture of SD evolution in this group is far from being complete. Lake Tanganyika and its affluent rivers are home to Astatotilapia burtoni, which belongs to the extremely successful East African cichlid lineage Haplochromini. Previously, in different families of an A. burtoni laboratory strain, an XYW system and an XY system have been described. The latter was also found in a second laboratory strain. In a laboratory-reared family descending from a population of the species' southern distribution, a second XY system was discovered. Yet, an analysis of sex chromosomes for the whole species distribution is missing. Here, we examined the genomes of 11 natural populations of A. burtoni, encompassing a wide range of its distribution, for sex-linked regions. We did not detect signs of differentiated sex chromosomes and also not the previously described sex chromosomal systems present in laboratory lines, suggesting different SD systems in the same species under natural and (long-term) artificial conditions. We suggest that SD in A. burtoni is more labile than previously assumed and consists of a combination of non-genetic, polygenic, or poorly differentiated sex chromosomes.

Metzincin metalloproteases in PGC migration and gonadal sex conversion

Development of a functional gonad includes migration of primordial germ cells (PGCs), differentiations of somatic and germ cells, formation of primary follicles or spermatogenic cysts with somatic gonadal cells, development and maturation of gametes, and subsequent releasing of mature germ cells. These processes require extensive cellular and tissue remodeling, as well as broad alterations of the surrounding extracellular matrix (ECM). Metalloproteases, including MMPs (matrix metalloproteases), ADAMs (a disintegrin and metalloproteinases), and ADAMTS (a disintegrin and metalloproteinase with thrombospondin motifs), are suggested to have critical roles in the remodeling of the ECM during gonad development. However, few research articles and reviews are available on the functions and mechanisms of metalloproteases in remodeling gonadal ECM, gonadal development, or gonadal differentiation. Moreover, most studies focused on the roles of transcription and growth factors in early gonad development and primary sex determination, leaving a significant knowledge gap on how differentially expressed metalloproteases exert effects on the ECM, cell migration, development, and survival of germ cells during the development and differentiation of ovaries or testes. We will review gonad development with focus on the evidence of metalloprotease involvements, and with an emphasis on zebrafish as a model for studying gonadal sex differentiation and metalloprotease functions.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

Environmentally-induced sex reversal in fish with chromosomal vs. polygenic sex determination

Sex ratio depends on sex determination mechanisms and is a key demographic parameter determining population viability and resilience to natural and anthropogenic stressors. There is increasing evidence that the environment can alter sex ratio even in genetically sex-determined species (GSD), as elevated temperature can cause female-to-male sex reversal (neomales). Alarmingly, neomales are being discovered in natural populations of several fish, amphibian and reptile species worldwide. Understanding the basis of neomale development is important for conservation biology. Among GSD species, it is unknown whether those with chromosomal sex determination (CSD), the most common system, will better resist the influence of high temperature than those with polygenic sex determination (PSD). Here, we compared the effects of elevated temperature in two wild zebrafish strains, Nadia (NA) and Ekkwill (EKW), which have CSD with a ZZ/ZW system, against the AB laboratory strain, which has PSD. First, we uncovered novel sex genotypes and the results showed that, at control temperature, the masculinization rate roughly doubled with the addition of each Z chromosome, while some ZW and WW fish of the wild strains became neomales. Surprisingly, we found that at elevated temperatures WW fish were just as likely as ZW fish to become neomales and that all strains were equally susceptible to masculinization. These results demonstrate that the Z chromosome is not essential for male development and that the dose of W buffers masculinization at the control temperature but not at elevated temperature. Furthermore, at the elevated temperature the testes of neomales, but not of normal males, contained more spermatozoa than at the control temperature. Our results show in an unprecedented way that, in a global warming scenario, CSD species may not necessarily be better protected against the masculinizing effect of elevated temperature than PSD species, and reveal genotype-by-temperature interactions in male sex determination and spermatogenesis.

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas

Background

The fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation.

Methods

Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950.

Results

We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (F_ST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

Development of a High-Density 665 K SNP Array for Rainbow Trout Genome-Wide Genotyping

Single nucleotide polymorphism (SNP) arrays, also named « SNP chips », enable very large numbers of individuals to be genotyped at a targeted set of thousands of genome-wide identified markers. We used preexisting variant datasets from USDA, a French commercial line and 30X-coverage whole genome sequencing of INRAE isogenic lines to develop an Affymetrix 665 K SNP array (HD chip) for rainbow trout. In total, we identified 32,372,492 SNPs that were polymorphic in the USDA or INRAE databases. A subset of identified SNPs were selected for inclusion on the chip, prioritizing SNPs whose flanking sequence uniquely aligned to the Swanson reference genome, with homogenous repartition over the genome and the highest Minimum Allele Frequency in both USDA and French databases. Of the 664,531 SNPs which passed the Affymetrix quality filters and were manufactured on the HD chip, 65.3% and 60.9% passed filtering metrics and were polymorphic in two other distinct French commercial populations in which, respectively, 288 and 175 sampled fish were genotyped. Only 576,118 SNPs mapped uniquely on both Swanson and Arlee reference genomes, and 12,071 SNPs did not map at all on the Arlee reference genome. Among those 576,118 SNPs, 38,948 SNPs were kept from the commercially available medium-density 57 K SNP chip. We demonstrate the utility of the HD chip by describing the high rates of linkage disequilibrium at 2–10 kb in the rainbow trout genome in comparison to the linkage disequilibrium observed at 50–100 kb which are usual distances between markers of the medium-density chip.

Sex-specific differences in zebrafish brains

In this systematic review, we highlight the differences between the male and female zebrafish brains to understand their differentiation and their use in studying sex-specific neurological diseases. Male and female brains display subtle differences at the cellular level which may be important in driving sex-specific signaling. Sex differences in the brain have been observed in humans as well as in non-human species. However, the molecular mechanisms of brain sex differentiation remain unclear. The classical model of brain sex differentiation suggests that the steroid hormones derived from the gonads are the primary determinants in establishing male and female neural networks. Recent studies indicate that the developing brain shows sex-specific differences in gene expression prior to gonadal hormone action. Hence, genetic differences may also be responsible for differentiating the brain into male and female types. Understanding the signaling mechanisms involved in brain sex differentiation could help further elucidate the sex-specific incidences of certain neurological diseases. The zebrafish model could be appropriate for enhancing our understanding of brain sex differentiation and the signaling involved in neurological diseases. Zebrafish brains show sex-specific differences at the hormonal level, and recent advances in RNA sequencing have highlighted critical sex-specific differences at the transcript level. The differences are also evident at the cellular and metabolite levels, which could be important in organizing sex-specific neuronal signaling. Furthermore, in addition to having one ortholog for 70% of the human gene, zebrafish also shares brain structural similarities with other higher eukaryotes, including mammals. Hence, deciphering brain sex differentiation in zebrafish will help further enhance the diagnostic and pharmacological intervention of neurological diseases.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Control of meiotic chromosomal bouquet and germ cell morphogenesis by the zygotene cilium

A hallmark of meiosis is chromosomal pairing, which requires telomere tethering and rotation on the nuclear envelope via microtubules, driving chromosome homology searches. Telomere pulling toward the centrosome forms the "zygotene chromosomal bouquet". Here, we identified the "zygotene cilium" in oocytes. This cilium provides a cable system for the bouquet machinery, extending throughout the germline cyst. Using zebrafish mutants and live manipulations, we demonstrate that the cilium anchors the centrosome to counterbalance telomere pulling. The cilium is essential for bouquet and synaptonemal complex formation, oogenesis, ovarian development, and fertility. Thus, a cilium represents a conserved player in zebrafish and mouse meiosis, which sheds light on reproductive aspects in ciliopathies, and suggests that cilia can control chromosomal dynamics.

Zebrafish Establish Female Germ Cell Identity by Advancing Cell Proliferation and Meiosis

Zebrafish is a popular research model; but its mechanism of sex determination is unclear and the sex of juvenile fish cannot be distinguished. To obtain fish with defined sex, we crossed domesticated zebrafish with the Nadia strain that has a female-dominant W segment. These fish were placed on a ziwi:GFP background to facilitate sorting of fluorescent germ cells for transcriptomic analysis. We analyzed the transcriptomes of germ cells at 10–14 days postfertilization (dpf), when sex dimorphic changes started to appear. Gene ontology showed that genes upregulated in the 10-dpf presumptive females are involved in cell cycles. This correlates with our detection of increased germ cell numbers and proliferation. We also detected upregulation of meiotic genes in the presumptive females at 14 dpf. Disruption of a meiotic gene, sycp3, resulted in sex reversal to infertile males. The germ cells of sycp3 mutants could not reach diplotene and underwent apoptosis. Preventing apoptosis by disrupting tp53 restored female characteristics in sycp3 mutants, demonstrating that adequate germ cells are required for female development. Thus, our transcriptome and gene mutation demonstrate that initial germ cell proliferation followed by meiosis is the hallmark of female differentiation in zebrafish.

Bloom syndrome helicase contributes to germ line development and longevity in zebrafish

RecQ helicases—also known as the “guardians of the genome”—play crucial roles in genome integrity maintenance through their involvement in various DNA metabolic pathways. Aside from being conserved from bacteria to vertebrates, their importance is also reflected in the fact that in humans impaired function of multiple RecQ helicase orthologs are known to cause severe sets of problems, including Bloom, Werner, or Rothmund-Thomson syndromes. Our aim was to create and characterize a zebrafish (Danio rerio) disease model for Bloom syndrome, a recessive autosomal disorder. In humans, this syndrome is characterized by short stature, skin rashes, reduced fertility, increased risk of carcinogenesis, and shortened life expectancy brought on by genomic instability. We show that zebrafish blm mutants recapitulate major hallmarks of the human disease, such as shortened lifespan and reduced fertility. Moreover, similarly to other factors involved in DNA repair, some functions of zebrafish Blm bear additional importance in germ line development, and consequently in sex differentiation. Unlike fanc genes and rad51, however, blm appears to affect its function independent of tp53. Therefore, our model will be a valuable tool for further understanding the developmental and molecular attributes of this rare disease, along with providing novel insights into the role of genome maintenance proteins in somatic DNA repair and fertility.

17β-Trenbolone binds to androgen receptor, decreases number of primordial germ cells, modulates expression of genes related to sexual differentiation, and affects sexual differentiation in zebrafish (Danio rerio)

Exposure to 17β-trenbolone caused a skewed sex ratio in fish. However, the molecular initiating event and key molecular event(s) remain unknown. In this study, zebrafish were exposed to 17β-trenbolone at nominal concentrations of 2 ng/L, 20 ng/L, 200 ng/L, and 2000 ng/L from fertilization to 60 days post fertilization (dpf). First, the sex ratio at 60 dpf was calculated to evaluate adverse outcomes on sexual differentiation. 17β-Trenbolone caused a skewed sex ratio toward males, with intersex individuals observed in the 20 ng/L group and all-male populations found in the 200 ng/L and 2000 ng/L groups. Then, the distribution and number of primordial germ cells, the expression of sex differentiation-related genes, and plasma vitellogenin concentrations were detected in wild-type zebrafish and the EGFP-nanos-3'UTR transgenic line using whole-mount in situ hybridization, real-time PCR, EGFP fluorescence quantification, and enzyme-linked immunosorbent assay. The results indicated that 17β-trenbolone exposure decreased the number of primordial germ cells at 1 dpf and 3 dpf, decreased expression of ovarian differentiation-related genes foxl2 and cyp19a1a at 60 dpf, increased expression of testis differentiation-related genes dmrt1, sox9a, and amh at 60 dpf, and decreased plasma vitellogenin levels at 60 dpf, revealing the key molecular events at different time points involved in affected sexual differentiation by 17β-trenbolone. Finally, molecular docking showed that 17β-trenbolone docked into ligand-binding domain of zebrafish androgen receptor with high binding energy (-3.72 kcal/mol), suggesting that binding to androgen receptor is the molecular initiating event affecting sexual differentiation by 17β-trenbolone. We found that 17β-trenbolone can bind to the zebrafish androgen receptor, decrease the number of primordial germ cells during the early embryonic stage, modulate the expression of genes related to sexual differentiation during gonadal differentiation, and eventually cause a skewed sex ratio toward males in adults.

Sexual determination in zebrafish

Zebrafish have emerged as a major model organism to study vertebrate reproduction due to their high fecundity and external development of eggs and embryos. The mechanisms through which zebrafish determine their sex have come under extensive investigation, as they lack a definite sex-determining chromosome and appear to have a highly complex method of sex determination. Single-gene mutagenesis has been employed to isolate the function of genes that determine zebrafish sex and regulate sex-specific differentiation, and to explore the interactions of genes that promote female or male sexual fate. In this review, we focus on recent advances in understanding of the mechanisms, including genetic and environmental factors, governing zebrafish sex development with comparisons to gene functions in other species to highlight conserved and potentially species-specific mechanisms for specifying and maintaining sexual fate.

Labile sex chromosomes in the Australian freshwater fish family Percichthyidae

Sex-specific ecology has management implications, but rapid sex-chromosome turnover in fishes hinders sex-marker development for monomorphic species. We used annotated genomes and reduced-representation sequencing data for two Australian percichthyids, Macquarie perch Macquaria australasica and golden perch M. ambigua, and whole genome resequencing for 50 Macquarie perch of each sex, to identify sex-linked loci and develop an affordable sexing assay. In silico pool-seq tests of 1,492,004 Macquarie perch SNPs revealed that a 275-kb scaffold was enriched for gametologous loci. Within this scaffold, 22 loci were sex-linked in a predominantly XY system, with females being homozygous for the X-linked allele at all 22, and males having the Y-linked allele at >7. Seven XY-gametologous loci (all males, but no females, are heterozygous or homozygous for the male-specific allele) were within a 146-bp region. A PCR-RFLP sexing assay targeting one Y-linked SNP, tested in 66 known-sex Macquarie perch and two of each sex of three confamilial species, plus amplicon sequencing of 400 bp encompassing the 146-bp region, revealed that the few sex-linked positions differ between species and between Macquarie perch populations. This indicates sex-chromosome lability in Percichthyidae, supported by nonhomologous scaffolds containing sex-linked loci for Macquarie- and golden perches. The present resources facilitate genomic research in Percichthyidae, including formulation of hypotheses about candidate genes of interest such as transcription factor SOX1b that occurs in the 275-kb scaffold ~38 kb downstream of the 146-bp region containing seven XY-gametologous loci. Sex-linked markers will be useful for determining genetic sex in some populations and studying sex chromosome turnover.

Meiotic Chromosome Dynamics in Zebrafish

Recent studies in zebrafish have revealed key features of meiotic chromosome dynamics, including clustering of telomeres in the bouquet configuration, biogenesis of chromosome axis structures, and the assembly and disassembly of the synaptonemal complex that aligns homologs end-to-end. The telomere bouquet stage is especially pronounced in zebrafish meiosis and sub-telomeric regions play key roles in mediating pairing and homologous recombination. In this review, we discuss the temporal progression of these events in meiosis prophase I and highlight the roles of proteins associated with meiotic chromosome architecture in homologous recombination. Finally, we discuss the interplay between meiotic mutants and gonadal sex differentiation and future research directions to study meiosis in living cells, including cell culture.

Identification of chromatin states during zebrafish gastrulation using CUT&RUN and CUT&Tag

BACKGROUND

Cell fate decisions are governed by interactions between sequence-specific transcription factors and a dynamic chromatin landscape. Zebrafish offer a powerful system for probing the mechanisms that drive these cell fate choices, especially in the context of early embryogenesis. However, technical challenges associated with conventional methods for chromatin profiling have slowed progress toward understanding the exact relationships between chromatin changes, transcription factor binding, and cellular differentiation during zebrafish embryogenesis.

RESULTS

To overcome these challenges, we adapted the chromatin profiling methods Cleavage Under Targets and Release Using Nuclease (CUT&RUN) and CUT&Tag for use in zebrafish and applied these methods to generate high-resolution enrichment maps for H3K4me3, H3K27me3, H3K9me3, RNA polymerase II, and the histone variant H2A.Z using tissue isolated from whole, mid-gastrula stage embryos. Using this data, we identify a subset of genes that may be bivalently regulated during both zebrafish and mouse gastrulation, provide evidence for an evolving H2A.Z landscape during embryo development, and demonstrate the effectiveness of CUT&RUN for detecting H3K9me3 enrichment at repetitive sequences.

CONCLUSIONS

Our results demonstrate the power of combining CUT&RUN and CUT&Tag methods with the strengths of the zebrafish system to define emerging chromatin landscapes in the context of vertebrate embryogenesis.

Diversity of Modes of Reproduction and Sex Determination Systems in Invertebrates, and the Putative Contribution of Genetic Conflict

About eight million animal species are estimated to live on Earth, and all except those belonging to one subphylum are invertebrates. Invertebrates are incredibly diverse in their morphologies, life histories, and in the range of the ecological niches that they occupy. A great variety of modes of reproduction and sex determination systems is also observed among them, and their mosaic-distribution across the phylogeny shows that transitions between them occur frequently and rapidly. Genetic conflict in its various forms is a long-standing theory to explain what drives those evolutionary transitions. Here, we review (1) the different modes of reproduction among invertebrate species, highlighting sexual reproduction as the probable ancestral state; (2) the paradoxical diversity of sex determination systems; (3) the different types of genetic conflicts that could drive the evolution of such different systems.

Female-specific genomic regions and molecular sex identification of the clearhead icefish (Protosalanx hyalocranius)

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07830-9.

Masculinization of Zebrafish Through Partial Depletion of Primordial Germ Cells by Injecting Diluted Morpholino Oligonucleotides into Embryos

The regulation of reproduction in zebrafish, the prime model of fish research, is not fully understood. An efficient tool to gain a better understanding of this complicated process is utilization of severely sex-biased families or groups. Here, we describe a method for partial depletion of primordial germ cells (PGCs) that leads to eventual masculinization of zebrafish. The technique is based on injecting early embryos with diluted morpholino oligonucleotides that temporarily interfere with the production of Dead end (Dnd), an RNA-binding protein essential for PGC survival. In addition, we also propose the use of eviscerated trunk, as a suitable alternative for examining gonadal expression in juvenile zebrafish.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of “personalized medicine.” Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Using RAD-seq to develop sex-linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double-digest restriction site-associated DNA sequencing (ddRAD-seq) to identify sex-linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD-seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non-native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female-linked and 19 putatively male-linked sequences. Four female- and eight male-linked markers amplified in all three life cycle stages. Using one female- and one male-linked marker that were sex-specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non-native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD-seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex-linked markers in other haplodiplontic macroalgae for which genomes are lacking.

Knockout of myoc Provides Evidence for the Role of Myocilin in Zebrafish Sex Determination Associated with Wnt Signalling Downregulation

Myocilin is a secreted glycoprotein with a poorly understood biological function and it is mainly known as the first glaucoma gene. To explore the normal role of this protein in vivo we developed a myoc knockout (KO) zebrafish line using CRISPR/Cas9 genome editing. This line carries a homozygous variant (c.236_239delinsAAAGGGGAAGGGGA) that is predicted to result in a loss-of-function of the protein because of a premature termination codon p.(V75EfsX60) that resulted in a significant reduction of myoc mRNA levels. Immunohistochemistry showed the presence of myocilin in wild-type embryonic (96 h post-fertilization) anterior segment eye structures and caudal muscles. The protein was also detected in different adult ocular and non-ocular tissues. No gross macroscopic or microscopic alterations were identified in the KO zebrafish, but, remarkably, we observed absence of females among the adult KO animals and apoptosis in the immature juvenile gonad (28 dpf) of these animals, which is characteristic of male development. Transcriptomic analysis showed that adult KO males overexpressed key genes involved in male sex determination and presented differentially expressed Wnt signalling genes. These results show that myocilin is required for ovary differentiation in zebrafish and provides in vivo support for the role of myocilin as a Wnt signalling pathway modulator. In summary, this myoc KO zebrafish line can be useful to investigate the elusive function of this protein, and it provides evidence for the unexpected function of myocilin as a key factor in zebrafish sex determination.

Network architecture and sex chromosome turnovers: Do epistatic interactions shape patterns of sex chromosome replacement?

Recent studies have revealed an astonishing diversity of sex chromosomes in many vertebrate lineages, prompting questions about the mechanisms of sex chromosome turnover. While there is considerable population genetic theory about the evolutionary forces promoting sex chromosome replacement, this theory has not yet been integrated with our understanding of the molecular and developmental genetics of sex determination. Here, we review recent data to examine four questions about how the structure of gene networks influences the evolution of sex determination. We argue that patterns of epistasis, arising from the structure of genetic networks, may play an important role in regulating the rates and patterns of sex chromosome replacement.

A map of cis-regulatory elements and 3D genome structures in zebrafish

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species¹. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

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.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

Loss of dmrt1 restores zebrafish female fates in the absence of cyp19a1a but not rbpms2a/b

Sex determination and differentiation is a complex process regulated by multiple factors, including factors from the germline or surrounding somatic tissue. In zebrafish, sex-determination involves establishment of a bipotential ovary that undergoes sex-specific differentiation and maintenance to form the functional adult gonad. However, the relationships among these factors are not fully understood. Here, we identify potential Rbpms2 targets and apply genetic epistasis experiments to decipher the genetic hierarchy of regulators of sex-specific differentiation. We provide genetic evidence that the crucial female factor rbpms2 is epistatic to the male factor dmrt1 in terms of adult sex. Moreover, the role of Rbpms2 in promoting female fates extends beyond repression of Dmrt1, as Rbpms2 is essential for female differentiation even in the absence of Dmrt1. In contrast, female fates can be restored in mutants lacking both cyp19a1a and dmrt1, and prolonged in bmp15 mutants in the absence of dmrt1. Taken together, this work indicates that cyp19a1a-mediated suppression of dmrt1 establishes a bipotential ovary and initiates female fate acquisition. Then, after female fate specification, Cyp19a1a regulates subsequent oocyte maturation and sustains female fates independently of Dmrt1 repression.

Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758)

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.

One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark

With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.

Sex Chromosome Evolution: So Many Exceptions to the Rules

Genomic analysis of many nonmodel species has uncovered an incredible diversity of sex chromosome systems, making it possible to empirically test the rich body of evolutionary theory that describes each stage of sex chromosome evolution. Classic theory predicts that sex chromosomes originate from a pair of homologous autosomes and recombination between them is suppressed via inversions to resolve sexual conflict. The resulting degradation of the Y chromosome gene content creates the need for dosage compensation in the heterogametic sex. Sex chromosome theory also implies a linear process, starting from sex chromosome origin and progressing to heteromorphism. Despite many convergent genomic patterns exhibited by independently evolved sex chromosome systems, and many case studies supporting these theoretical predictions, emerging data provide numerous interesting exceptions to these long-standing theories, and suggest that the remarkable diversity of sex chromosomes is matched by a similar diversity in their evolution. For example, it is clear that sex chromosome pairs are not always derived from homologous autosomes. In addition, both the cause and the mechanism of recombination suppression between sex chromosome pairs remain unclear, and it may be that the spread of recombination suppression is a more gradual process than previously thought. It is also clear that dosage compensation can be achieved in many ways, and displays a range of efficacy in different systems. Finally, the remarkable turnover of sex chromosomes in many systems, as well as variation in the rate of sex chromosome divergence, suggest that assumptions about the inevitable linearity of sex chromosome evolution are not always empirically supported, and the drivers of the birth-death cycle of sex chromosome evolution remain to be elucidated. Here, we concentrate on how the diversity in sex chromosomes across taxa highlights an equal diversity in each stage of sex chromosome evolution.

Mixing and Matching Chromosomes during Female Meiosis

Meiosis is a key event in the manufacturing of an oocyte. During this process, the oocyte creates a set of unique chromosomes by recombining paternal and maternal copies of homologous chromosomes, and by eliminating one set of chromosomes to become haploid. While meiosis is conserved among sexually reproducing eukaryotes, there is a bewildering diversity of strategies among species, and sometimes within sexes of the same species, to achieve proper segregation of chromosomes. Here, we review the very first steps of meiosis in females, when the maternal and paternal copies of each homologous chromosomes have to move, find each other and pair. We explore the similarities and differences observed in C. elegans, Drosophila, zebrafish and mouse females.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Transcriptome profiling analysis of sex-based differentially expressed mRNAs and lncRNAs in the brains of mature zebrafish (Danio rerio)

Background

Similar to humans, the zebrafish brain plays a central role in regulating sexual reproduction, maturation and sexual behavior. However, systematic studies of the dimorphic patterns of gene expression in the brain of male and female zebrafish are lacking.

Results

In this study, the mRNA and lncRNA expression profiles were obtained from the brain tissue samples of the three male and three female zebrafish by high-throughput transcriptome sequencing. We identified a total of 108 mRNAs and 50 lncRNAs with sex-based differential expression. We randomly selected four differentially expressed genes for RT-qPCR verification and the results certified that the expression pattern showed a similar trend between RNA-seq and RT-qPCR results. Protein-protein interaction network analysis, Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to obtain the biological significance of differentially expressed mRNA in the brain dimorphism of zebrafish. Finally, a Pearson correlation analysis was performed to construct the co-expression network of the mRNAs and lncRNAs.

Conclusions

We found that 12 new lncRNAs not only have significant gender specificity in the brain of zebrafish, and this finding may provide a clue to further study of the functional difference between male and female zebrafish brain.

Let's talk about sex: A rigorous statistical framework to assign the sex of individuals from reduced-representation sequencing data

Molecular markers have been used to identify the sex of sampled individuals for several decades, but the time-consuming development phase prevented their application in many systems. Recently, a growing number of papers have applied reduced-representation sequencing (RRS) protocols to the identification of sex-specific markers without the use of test crosses or prior genomic information. While such an approach has great advantages in terms of versatility and ease of use, the "shotgun sequencing" nature of RRS data sets leads to a high amount of missing data, which results in statistical challenges to the confident assignment of sex to individuals. In this issue of Molecular Ecology Resources, Stovall et al. (Molecular Ecology Resources, 18, 2018) provide a statistical framework to answer two questions: (1) how many individuals of one sex only must possess a genotype for this locus to be considered significantly sex-specific? and (2) How many sex-specific loci must an individual of unknown sex possess (in a given data set) to be confidently assigned a sex? The statistical pipeline introduced, and applied to samples of New Zealand fur seal (Arctocephalus forsteri) to identify 90 sex-specific loci, should be broadly applicable to a large number of species and constitutes a nice addition to the molecular ecology toolkit in the genomics era.

Hemimetabolous insects elucidate the origin of sexual development via alternative splicing

Insects are the only known animals in which sexual differentiation is controlled by sex-specific splicing. The doublesex transcription factor produces distinct male and female isoforms, which are both essential for sex-specific development. dsx splicing depends on transformer, which is also alternatively spliced such that functional Tra is only present in females. This pathway has evolved from an ancestral mechanism where dsx was independent of tra and expressed and required only in males. To reconstruct this transition, we examined three basal, hemimetabolous insect orders: Hemiptera, Phthiraptera, and Blattodea. We show that tra and dsx have distinct functions in these insects, reflecting different stages in the changeover from a transcription-based to a splicing-based mode of sexual differentiation. We propose that the canonical insect tra-dsx pathway evolved via merger between expanding dsx function (from males to both sexes) and narrowing tra function (from a general splicing factor to dedicated regulator of dsx).

Scanning of Genetic Variants and Genetic Mapping of Phenotypic Traits in Gilthead Sea Bream Through ddRAD Sequencing

Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.

A Search for Sex-Linked Loci in the Agamid Lizard, Calotes versicolor

The Indian garden lizard, Calotes versicolor, lacks cytologically recognizable sex chromosomes, and its mechanism of sex determination is unclear. We evaluated genotype-to-sex-phenotype association using RAD-seq in wild-caught males and females, 30 of each sex. Of 210,736 unique, 96-nt long RAD-tags, 48% contained polymorphisms, 23% of which were present in at least 40 of 60 individuals. Twenty one RAD-tags neared, but none achieved, the inclusion criteria for sex enrichment, as expected if C. versicolor lacks highly differentiated sex chromosomes. Three RAD-tags with alleles most strongly associated with sex tended to be heterozygous in females and to lack male-specific alleles, suggesting a ZW female/ZZ male system. Putative female alleles, however, were present in some males and lacking from some females, suggesting either recombination between these markers and the sex locus or sex reversal due to environmental or genetic factors. Paired-end, 250-nt reads from 1 male provided a fragmented draft genome assembly. Four sex-associated RAD-tags were identical to portions of 4 unique C. versicolor genomic contigs rather than linked to a single putative sex-linked region. The lack of strongly sex-linked loci coupled with weak evidence for temperature-associated sex determination intensifies the need for further investigation of the puzzling sex determination mechanism in C. versicolor.

Genome and Phylogenetic Analysis of Genes Involved in the Immune System of Solea senegalensis - Potential Applications in Aquaculture

Global aquaculture production continues to increase rapidly. One of the most important species of marine fish currently cultivated in Southern Europe is Solea senegalensis, reaching more than 300 Tn in 2017. In the present work, 14 Bacterial Artificial Chromosome (BAC) clones containing candidate genes involved in the immune system (b2m, il10, tlr3, tap1, tnfα, tlr8, trim25, lysg, irf5, hmgb2, calr, trim16, and mx), were examined and compared with other species using multicolor Fluorescence in situ Hybridization (mFISH), massive sequencing and bioinformatic analysis to determine the genomic surroundings and syntenic chromosomal conservation of the genomic region contained in each BAC clone. The mFISH showed that the groups of genes hmgb2-trim25-irf5-b2m; tlr3-lysg; tnfα-tap1, and il10-mx-trim16 were co-localized on the same chromosomes. Synteny results suggested that the studied BACs are placed in a smaller number of chromosomes in S. senegalensis that in other species. Phylogenetic analyses suggested that the evolutionary rate of immune system genes studied is similar among the taxa studied, given that the clustering obtained was in accordance with the accepted phylogenetic relationships among these species. This study contributes to a better understanding of the structure and function of the immune system of the Senegalese sole, which is essential for the development of new technologies and products to improve fish health and productivity.