Inducible degradation of dosage compensation protein DPY-27 facilitates isolation of Caenorhabditis elegans males for molecular and biochemical analyses

Biological sex affects numerous aspects of biology, yet how sex influences different biological processes have not been extensively studied at the molecular level. Caenorhabditis elegans, with both hermaphrodites (functionally females as adults) and males, is an excellent system to uncover how sex influences physiology. Here, we describe a method to isolate large quantities of C. elegans males by conditionally degrading DPY-27, a component of the dosage compensation complex essential for hermaphrodite, but not male, development. We show that germ cells from males isolated following DPY-27 degradation undergo meiosis and spermiogenesis like wild type and these males are competent to mate and sire viable offspring. We further demonstrate the efficacy of this system by analyzing gene expression and performing affinity pull-downs from male worm extracts.

Infertility due to defective sperm flagella caused by an intronic deletion in DNAH17 that perturbs splicing

Artificial insemination in pig (Sus scrofa domesticus) breeding involves the evaluation of the semen quality of breeding boars. Ejaculates that fulfill predefined quality requirements are processed, diluted and used for inseminations. Within short time, eight Swiss Large White boars producing immotile sperm that had multiple morphological abnormalities of the sperm flagella were noticed at a semen collection center. The eight boars were inbred on a common ancestor suggesting that the novel sperm flagella defect is a recessive trait. Transmission electron microscopy cross-sections revealed that the immotile sperm had disorganized flagellar axonemes. Haplotype-based association testing involving microarray-derived genotypes at 41,094 SNPs of six affected and 100 fertile boars yielded strong association (P = 4.22 × 10-15) at chromosome 12. Autozygosity mapping enabled us to pinpoint the causal mutation on a 1.11 Mb haplotype located between 3,473,632 and 4,587,759 bp. The haplotype carries an intronic 13-bp deletion (Chr12:3,556,401-3,556,414 bp) that is compatible with recessive inheritance. The 13-bp deletion excises the polypyrimidine tract upstream exon 56 of DNAH17 (XM_021066525.1: c.8510-17_8510-5del) encoding dynein axonemal heavy chain 17. Transcriptome analysis of the testis of two affected boars revealed that the loss of the polypyrimidine tract causes exon skipping which results in the in-frame loss of 89 amino acids from DNAH17. Disruption of DNAH17 impairs the assembly of the flagellar axoneme and manifests in multiple morphological abnormalities of the sperm flagella. Direct gene testing may now be implemented to monitor the defective allele in the Swiss Large White population and prevent the frequent manifestation of a sterilizing sperm tail disorder in breeding boars.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

A novel c.1759T>G variant in follicle-stimulating hormone-receptor gene is concordant with male determination in the flathead grey mullet (Mugil cephalus)

Various master key regulators (MKRs) that control a binary switch of sex determination (SD) have been found in fish; these provide an excellent model for the study of vertebrate genetic SD. The SD region in flathead grey mullet has been previously mapped to a 1 Mbp region harboring 27 genes, of which one is follicle-stimulating hormone receptor (fshr). Although this gene is involved in gonad differentiation and function, it has not been considered as an MKR of SD. We systematically investigated polymorphism in mullet fshr using DNA shotgun sequences, and compared them between males and females. Capable of encoding nonconservative amino acid substitutions, c.1732G>A and c.1759T>G exhibited association with sex on a population level (N = 83; P ≤ 6.7 × 10-19). Hence, 1732 A and 1759 G represent a male-specific haplotype of the gene, designated as "fshry." Additional flanking SNPs showed a weaker degree of association with sex, delimiting the SD critical region to 143 nucleotides on exon 14. Lack of homozygotes for fshry, and the resulting divergence from Hardy-Weinberg equilibrium (N = 170; P ≤ 3.9 × 10-5), were compatible with a male heterogametic model (XY/XX). Capable of replacing a phenylalanine with valine, c.1759T>G alters a conserved position across the sixth transmembrane domain of vertebrate FSHRs. Amino acid substitutions in this position in vertebrates are frequently associated with constant receptor activation and consequently with FSH/FSHR signaling alteration; thus, indicating a potential role of fshr as an MKR of SD.

Meiotic Double-Strand Break Processing and Crossover Patterning Are Regulated in a Sex-Specific Manner by BRCA1-BARD1 in Caenorhabditis elegans

Meiosis is regulated in a sex-specific manner to produce two distinct gametes, sperm and oocytes, for sexual reproduction. To determine how meiotic recombination is regulated in spermatogenesis, we analyzed the meiotic phenotypes of mutants in the tumor suppressor E3 ubiquitin ligase BRC-1-BRD-1 complex in Caenorhabditis elegans male meiosis. Unlike in mammals, this complex is not required for meiotic sex chromosome inactivation, the process whereby hemizygous sex chromosomes are transcriptionally silenced. Interestingly, brc-1 and brd-1 mutants show meiotic recombination phenotypes that are largely opposing to those previously reported for female meiosis. Fewer meiotic recombination intermediates marked by the recombinase RAD-51 were observed in brc-1 and brd-1 mutants, and the reduction in RAD-51 foci could be suppressed by mutation of nonhomologous-end-joining proteins. Analysis of GFP::RPA-1 revealed fewer foci in the brc-1brd-1 mutant and concentration of BRC-1-BRD-1 to sites of meiotic recombination was dependent on DNA end resection, suggesting that the complex regulates the processing of meiotic double-strand breaks to promote repair by homologous recombination. Further, BRC-1-BRD-1 is important to promote progeny viability when male meiosis is perturbed by mutations that block the pairing and synapsis of different chromosome pairs, although the complex is not required to stabilize the RAD-51 filament as in female meiosis under the same conditions. Analyses of crossover designation and formation revealed that BRC-1-BRD-1 inhibits supernumerary COs when meiosis is perturbed. Together, our findings suggest that BRC-1-BRD-1 regulates different aspects of meiotic recombination in male and female meiosis.

Evidence That Runt Acts as a Counter-Repressor of Groucho During Drosophila melanogaster Primary Sex Determination

Runx proteins are bifunctional transcription factors that both repress and activate transcription in animal cells. Typically, Runx proteins work in concert with other transcriptional regulators, including co-activators and co-repressors to mediate their biological effects. In Drosophila melanogaster the archetypal Runx protein, Runt, functions in numerous processes including segmentation, neurogenesis and sex determination. During primary sex determination Runt acts as one of four X-linked signal element (XSE) proteins that direct female-specific activation of the establishment promoter (Pe) of the master regulatory gene Sex-lethal (Sxl). Successful activation of SxlPe requires that the XSE proteins overcome the repressive effects of maternally deposited Groucho (Gro), a potent co-repressor of the Gro/TLE family. Runx proteins, including Runt, contain a C-terminal peptide, VWRPY, known to bind to Gro/TLE proteins to mediate transcriptional repression. We show that Runt’s VWRPY co-repressor-interaction domain is needed for Runt to activate SxlPe. Deletion of the Gro-interaction domain eliminates Runt-ability to activate SxlPe, whereas replacement with a higher affinity, VWRPW, sequence promotes Runt-mediated transcription. This suggests that Runt may activate SxlPe by antagonizing Gro function, a conclusion consistent with earlier findings that Runt is needed for Sxl expression only in embryonic regions with high Gro activity. Surprisingly we found that Runt is not required for the initial activation of SxlPe. Instead, Runt is needed to keep SxlPe active during the subsequent period of high-level Sxl transcription suggesting that Runt helps amplify the difference between female and male XSE signals by counter-repressing Gro in female, but not in male, embryos.

Genomic Structure, Evolutionary Origins, and Reproductive Function of a Large Amplified Intrinsically Disordered Protein-Coding Gene on the X Chromosome (Laidx) in Mice

Mouse sex chromosomes are enriched for co-amplified gene families, present in tens to hundreds of copies. Co-amplification of Slx/Slxl1 on the X chromosome and Sly on the Y chromosome are involved in dose-dependent meiotic drive, however the role of other co-amplified genes remains poorly understood. Here we demonstrate that the co-amplified gene family on the X chromosome, Srsx, along with two additional partial gene annotations, is actually part of a larger transcription unit, which we name Laidx. Laidx is harbored in a 229 kb amplicon that represents the ancestral state as compared to a 525 kb Y-amplicon containing the rearranged Laidy. Laidx contains a 25,011 nucleotide open reading frame, predominantly expressed in round spermatids, predicted to encode an 871 kD protein. Laidx has orthologous copies with the rat and also the 825-MY diverged parasitic Chinese liver fluke, Clonorchis sinensis, the likely result of a horizontal gene transfer of rodent Laidx to an ancestor of the liver fluke. To assess the male reproductive functions of Laidx, we generated mice carrying a multi-megabase deletion of the Laidx-ampliconic region. Laidx-deficient male mice do not show detectable reproductive defects in fertility, fecundity, testis histology, and offspring sex ratio. We speculate that Laidx and Laidy represent a now inactive X vs. Y chromosome conflict that occurred in an ancestor of present day mice.

Sex Chromosome Evolution in Muscid Flies

Sex chromosomes and sex determining genes can evolve fast, with the sex-linked chromosomes often differing between closely related species. Population genetics theory has been developed and tested to explain the rapid evolution of sex chromosomes and sex determination. However, we do not know why the sex chromosomes are divergent in some taxa and conserved in others. Addressing this question requires comparing closely related taxa with conserved and divergent sex chromosomes to identify biological features that could explain these differences. Cytological karyotypes suggest that muscid flies (e.g., house fly) and blow flies are such a taxonomic pair. The sex chromosomes appear to differ across muscid species, whereas they are conserved across blow flies. Despite the cytological evidence, we do not know the extent to which muscid sex chromosomes are independently derived along different evolutionary lineages. To address that question, we used genomic and transcriptomic sequence data to identify young sex chromosomes in two closely related muscid species, horn fly (Haematobia irritans) and stable fly (Stomoxys calcitrans). We provide evidence that the nascent sex chromosomes of horn fly and stable fly were derived independently from each other and from the young sex chromosomes of the closely related house fly (Musca domestica). We present three different scenarios that could have given rise to the sex chromosomes of horn fly and stable fly, and we describe how the scenarios could be distinguished. Distinguishing between these scenarios in future work could identify features of muscid genomes that promote sex chromosome divergence.

The Drosophila Post-mating Response: Gene Expression and Behavioral Changes Reveal Perdurance and Variation in Cross-Tissue Interactions

Examining cross-tissue interactions is important for understanding physiology and homeostasis. In animals, the female gonad produces signaling molecules that act distally. We examine gene expression in Drosophila melanogaster female head tissues in 1) virgins without a germline compared to virgins with a germline, 2) post-mated females with and without a germline compared to virgins, and 3) post-mated females mated to males with and without a germline compared to virgins. In virgins, the absence of a female germline results in expression changes in genes with known roles in nutrient homeostasis. At one- and three-day(s) post-mating, genes that change expression are enriched with those that function in metabolic pathways, in all conditions. We systematically examine female post-mating impacts on sleep, food preference and re-mating, in the strains and time points used for gene expression analyses and compare to published studies. We show that post-mating, gene expression changes vary by strain, prompting us to examine variation in female re-mating. We perform a genome-wide association study that identifies several DNA polymorphisms, including four in/near Wnt signaling pathway genes. Together, these data reveal how gene expression and behavior in females are influenced by cross-tissue interactions, by examining the impact of mating, fertility, and genotype.

A ZZ/ZW Sex Chromosome System in Cephalochordate Amphioxus

Sex determination is remarkably variable among animals with examples of environmental sex determination, male heterogametic (XX/XY) and female heterogametic (ZZ/ZW) chromosomal sex determination, and other genetic mechanisms. The cephalochordate amphioxus occupies a key phylogenetic position as a basal chordate and outgroup to vertebrates, but its sex determination mechanism is unknown. During the course of generating Nodal mutants with transcription activator-like effector nucleases (TALENs) in amphioxus Branchiostoma floridae, serendipitously, we generated three mutant strains that reveal the sex determination mechanism of this animal. In one mutant strain, all heterozygous mutant offspring over three generations were female and all wild-type descendants were male. This pattern suggests the Nodal allele targeted is on a female-specific W chromosome. A second mutant showed the same W-linked inheritance pattern, with a female heterozygote passing the mutation only to daughters. In a third mutant strain, both male and female offspring could be heterozygous, but a female heterozygote passed the mutation only to sons. This pattern is consistent with the targeted allele being on a Z chromosome. We found an indel polymorphism linked to a Nodal allele present in most females, but no males in our cultured population. Together, these results indicate that Nodal is sex chromosome-linked in B. floridae, and that B. floridae has a ZZ/ZW sex chromosome system.

Protection Against XY Gonadal Sex Reversal by a Variant Region on Mouse Chromosome 13

XY C57BL/6J (B6) mice harboring a Mus musculus domesticus-type Y chromosome (Y POS ), known as B6.Y POS mice, commonly undergo gonadal sex reversal and develop as phenotypic females. In a minority of cases, B6.Y POS males are identified and a proportion of these are fertile. This phenotypic variability on a congenic B6 background has puzzled geneticists for decades. Recently, a B6.Y POS colony was shown to carry a non-B6-derived region of chromosome 11 that protected against B6.Y POS sex reversal. Here. we show that a B6.Y POS colony bred and archived at the MRC Harwell Institute lacks the chromosome 11 modifier but instead harbors an ∼37 Mb region containing non-B6-derived segments on chromosome 13. This region, which we call Mod13, protects against B6.Y POS sex reversal in a proportion of heterozygous animals through its positive and negative effects on gene expression during primary sex determination. We discuss Mod13's influence on the testis determination process and its possible origin in light of sequence similarities to that region in other mouse genomes. Our data reveal that the B6.Y POS sex reversal phenomenon is genetically complex and the explanation of observed phenotypic variability is likely dependent on the breeding history of any local colony.

The Role of Anti-Müllerian Hormone in Testis Differentiation Reveals the Significance of the TGF-β Pathway in Reptilian Sex Determination

Anti-Müllerian hormone (Amh, or Müllerian-inhibiting substance, Mis), a member of TGF-β superfamily, has been well documented in some vertebrates as initiator or key regulator in sexual development, and particularly in fish. However, its functional role has not yet been identified in reptiles. Here, we characterized the Amh gene in the Chinese soft-shelled turtle Pelodiscus sinensis, a typical reptilian species exhibiting ZZ/ZW sex chromosomes. The messenger RNA of Amh was initially expressed in male embryonic gonads by stage 15, preceding gonadal sex differentiation, and exhibited a male-specific expression pattern throughout embryogenesis. Moreover, Amh was rapidly upregulated during female-to-male sex reversal induced by aromatase inhibitor letrozole. Most importantly, Amh loss of function by RNA interference led to complete feminization of genetic male (ZZ) gonads, suppression of the testicular marker Sox9, and upregulation of the ovarian regulator Cyp19a1 Conversely, overexpression of Amh in ZW embryos resulted in female-to-male sex reversal, characterized by the formation of a testis structure, ectopic activation of Sox9, and a remarkable decline in Cyp19a1 Collectively, these findings provide the first solid evidence that Amh is both necessary and sufficient to drive testicular development in a reptilian species, P. sinensis, highlighting the significance of the TGF-β pathway in reptilian sex determination.

Mapping of the Sex Determining Region on Linkage Group 12 of Guppy (Poecilia reticulata)

Poecilia reticulata is one of the most popular ornamental fish species with a higher demand for males due to their large colorful fins. The objectives of this study were mapping of the sex-determining (SD) region on linkage group 12 of guppy, and identification of a sex specific marker. We generated eight polymorphic microsatellite markers distributed along the distal 5.4 Mbp sequence of the previously identified SD region on linkage group (LG) 12. The markers were tested for association with sex in 156 individuals of the Red Blonde and Flame strains, and 126 progeny of four full-sibs Red Blonde families. A male-specific allele was found for microsatellite gu1066 at position of 25.3 Mbp on LG12 for both strains, and gu832 at position 24.4 Mbp for the Flame strain. Thus, a region of 0.9 Mbp between these markers, harboring 27 annotated genes, was selected for analysis. Based on association of copy number variation and sex determination we mapped a duplicated region between LGs 9 and 12, of 1.26 Mbp, containing 59 genes on LG12. The common interval between the segment bounded by gu1066 and gu832, and the duplicated region of 0.43 Mbp on LG12 harbors 11 genes of which 6 have multiple copies (54%). Growth arrest and DNA damage inducible gamma-like (GADD45G-like) is a plausible positional and functional candidate gene for its role in male fertility. We characterized the genomic structure of the gene in guppy, and found two isoforms; but no sex-biased differences were evident in genomic sequence and gene expression.

Fine Mapping Using Whole-Genome Sequencing Confirms Anti-Müllerian Hormone as a Major Gene for Sex Determination in Farmed Nile Tilapia (Oreochromis niloticus L.)

Nile tilapia (Oreochromis niloticus) is one of the most cultivated and economically important species in world aquaculture. Intensive production promotes the use of monosex animals, due to an important dimorphism that favors male growth. Currently, the main mechanism to obtain all-male populations is the use of hormones in feeding during larval and fry phases. Identifying genomic regions associated with sex determination in Nile tilapia is a research topic of great interest. The objective of this study was to identify genomic variants associated with sex determination in three commercial populations of Nile tilapia. Whole-genome sequencing of 326 individuals was performed, and a total of 2.4 million high-quality bi-allelic single nucleotide polymorphisms (SNPs) were identified after quality control. A genome-wide association study (GWAS) was conducted to identify markers associated with the binary sex trait (males = 1; females = 0). A mixed logistic regression GWAS model was fitted and a genome-wide significant signal comprising 36 SNPs, spanning a genomic region of 536 kb in chromosome 23 was identified. Ten out of these 36 genetic variants intercept the anti-Müllerian (Amh) hormone gene. Other significant SNPs were located in the neighboring Amh gene region. This gene has been strongly associated with sex determination in several vertebrate species, playing an essential role in the differentiation of male and female reproductive tissue in early stages of development. This finding provides useful information to better understand the genetic mechanisms underlying sex determination in Nile tilapia.

Genomic Data Reveal Conserved Female Heterogamety in Giant Salamanders with Gigantic Nuclear Genomes

Systems of genetic sex determination and the homology of sex chromosomes in different taxa vary greatly across vertebrates. Much progress remains to be made in understanding systems of genetic sex determination in non-model organisms, especially those with homomorphic sex chromosomes and/or large genomes. We used reduced representation genome sequencing to investigate genetic sex determination systems in the salamander family Cryptobranchidae (genera Cryptobranchus and Andrias), which typifies both of these inherent difficulties. We tested hypotheses of male- or female-heterogamety by sequencing hundreds of thousands of anonymous genomic regions in a panel of known-sex cryptobranchids and characterized patterns of presence/absence, inferred zygosity, and depth of coverage to identify sex-linked regions of these 56 gigabase genomes. Our results strongly support the hypothesis that all cryptobranchid species possess homologous systems of female heterogamety, despite maintenance of homomorphic sex chromosomes over nearly 60 million years. Additionally, we report a robust, non-invasive genetic assay for sex diagnosis in Cryptobranchus and Andrias which may have great utility for conservation efforts with these endangered salamanders. Co-amplification of these W-linked markers in both cryptobranchid genera provides evidence for long-term sex chromosome stasis in one of the most divergent salamander lineages. These findings inform hypotheses about the ancestral mode of sex determination in salamanders, but suggest that comparative data from other salamander families are needed. Our results further demonstrate that massive genomes are not necessarily a barrier to effective genome-wide sequencing and that the resulting data can be highly informative about sex determination systems in taxa with homomorphic sex chromosomes.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

Cyfip1 Haploinsufficiency Increases Compulsive-Like Behavior and Modulates Palatable Food Intake in Mice: Dependence on Cyfip2 Genetic Background, Parent-of Origin, and Sex

Binge eating (BE) is a heritable trait associated with eating disorders and involves episodes of rapid, large amounts of food consumption. We previously identified cytoplasmic FMR1-interacting protein 2 (Cyfip2) as a genetic factor underlying compulsive-like BE in mice. CYFIP2 is a homolog of CYFIP1 which is one of four paternally-deleted genes in patients with Type I Prader-Willi Syndrome (PWS), a neurodevelopmental disorder whereby 70% of cases involve paternal 15q11-q13 deletion. PWS symptoms include hyperphagia, obesity (if untreated), cognitive deficits, and obsessive-compulsive behaviors. We tested whether Cyfip1 haploinsufficiency (+/−) would enhance compulsive-like behavior and palatable food (PF) intake in a parental origin- and sex-dependent manner on two Cyfip2 genetic backgrounds, including the BE-prone C57BL/6N (Cyfip2^N/N) background and the BE-resistant C57BL/6J (Cyfip2^J/J) background. Cyfip1^+/− mice showed increased compulsive-like behavior on both backgrounds and increased PF intake on the Cyfip2^N/N background. In contrast, maternal Cyfip1 haploinsufficiency on the BE-resistant Cyfip2^J/J background induced a robust escalation in PF intake in wild-type Cyfip1^J/J males while having no effect in Cyfip1^J/- males. Notably, induction of behavioral phenotypes in wild-type males following maternal Fmr1^+/− has previously been reported. In the hypothalamus, there was a paternally-enhanced reduction in CYFIP1 protein whereas in the nucleus accumbens, there was a maternally-enhanced reduction in CYFIP1 protein. Nochange in FMR1 protein (FMRP) was observed in Cyfip1^+/− mice, regardless of parental origin. To summarize, Cyfip1 haploinsufficiency increased compulsive-like behavior and induced genetic background-dependent, sex-dependent, and parent-of-origin-dependent effects on PF consumption and CYFIP1 expression that could have relevance for neurodevelopmental and neuropsychiatric disorders.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

Little Evidence of Antagonistic Selection in the Evolutionary Strata of Fungal Mating-Type Chromosomes (Microbotryum lychnidis-dioicae)

Recombination suppression on sex chromosomes often extends in a stepwise manner, generating evolutionary strata of differentiation between sex chromosomes. Sexual antagonism is a widely accepted explanation for evolutionary strata, postulating that sets of genes beneficial in only one sex are successively linked to the sex-determining locus. The anther-smut fungus Microbotryum lychnidis-dioicae has mating-type chromosomes with evolutionary strata, only some of which link mating-type genes. Male and female roles are non-existent in this fungus, but mating-type antagonistic selection can also generate evolutionary strata, although the life cycle of the fungus suggests it should be restricted to few traits. Here, we tested the hypothesis that mating-type antagonism may have triggered recombination suppression beyond mating-type genes in M. lychnidis-dioicae by searching for footprints of antagonistic selection in evolutionary strata not linking mating-type loci. We found that these evolutionary strata (i) were not enriched in genes upregulated in the haploid phase, where cells are of alternative mating types, (ii) carried no gene differentially expressed between mating types, and (iii) carried no genes displaying footprints of specialization in terms of protein sequences (d_N/d_S) between mating types after recommended filtering. Without filtering, eleven genes showed signs of positive selection in the strata not linking mating-type genes, which constituted an enrichment compared to autosomes, but their functions were not obviously involved in antagonistic selection. Thus, we found no strong evidence that antagonistic selection has contributed to extending recombination suppression beyond mating-type genes. Alternative hypotheses should therefore be explored to improve our understanding of the sex-related chromosome evolution.

Chromosome-Wide Evolution and Sex Determination in the Three-Sexed Nematode Auanema rhodensis

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis, males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM (doublesex/mab-3) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis, instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

Identification of rfk-1, a Meiotic Driver Undergoing RNA Editing in Neurospora

Sk-2 is a meiotic drive element that was discovered in wild populations of Neurospora fungi over 40 years ago. While early studies quickly determined that Sk-2 transmits itself through sexual reproduction in a biased manner via spore killing, the genetic factors responsible for this phenomenon have remained mostly unknown. Here, we identify and characterize rfk-1, a gene required for Sk-2-based spore killing. The rfk-1 gene contains four exons, three introns, and two stop codons, the first of which undergoes RNA editing to a tryptophan codon during sexual development. Translation of an unedited rfk-1 transcript in vegetative tissue is expected to produce a 102-amino acid protein, whereas translation of an edited rfk-1 transcript in sexual tissue is expected to produce a protein with 130 amino acids. These findings indicate that unedited and edited rfk-1 transcripts exist and that these transcripts could have different roles with respect to the mechanism of meiotic drive by spore killing. Regardless of RNA editing, spore killing only succeeds if rfk-1 transcripts avoid silencing caused by a genome defense process called meiotic silencing by unpaired DNA (MSUD). We show that rfk-1's MSUD avoidance mechanism is linked to the genomic landscape surrounding the rfk-1 gene, which is located near the Sk-2 border on the right arm of chromosome III. In addition to demonstrating that the location of rfk-1 is critical to spore-killing success, our results add to accumulating evidence that MSUD helps protect Neurospora genomes from complex meiotic drive elements.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

The mir-35 Family Links Maternal Germline Sex to Embryonic Viability in Caenorhabditis elegans

The germline sex determination pathway in C. elegans determines whether germ cells develop as oocytes or sperm, with no previously known effect on viability. The mir-35 family of microRNAs are expressed in the C. elegans germline and embryo and are essential for both viability and normal hermaphroditic sex determination, preventing aberrant male gene expression in XX hermaphrodite embryos. Here we show that combining feminizing mutations with partial loss of function of the mir-35 family results in enhanced penetrance embryonic lethality that preferentially kills XO animals. This lethal phenotype is due to altered signaling through the germline sex determination pathway, and maternal germline feminization is sufficient to induce enhanced lethality. These findings reveal a surprising pleiotropy of sperm-fate promoting pathways on organismal viability. Overall, our results demonstrate an unexpectedly strong link between sex determination and embryonic viability, and suggest that in wild type animals, mir-35 family members buffer against misregulation of pathways outside the sex determination program, allowing for clean sex reversal rather than deleterious effects of perturbing sex determination genes.

The Role of miRNAs in Drosophila melanogaster Male Courtship Behavior

Drosophila melanogaster courtship, although stereotypical, continually changes based on cues received from the courtship subject. Such adaptive responses are mediated via rapid and widespread transcriptomic reprogramming, a characteristic now widely attributed to microRNAs (miRNAs), along with other players. Here, we conducted a large-scale miRNA knockout screen to identify miRNAs that affect various parameters of male courtship behavior. Apart from identifying miRNAs that impact male-female courtship, we observed that miR-957 mutants performed significantly increased male-male courtship and "chaining" behavior, whereby groups of males court one another. We tested the effect of miR-957 reduction in specific neuronal cell clusters, identifying miR-957 activity in Doublesex (DSX)-expressing and mushroom body clusters as an important regulator of male-male courtship interactions. We further characterized the behavior of miR-957 mutants and found that these males court male subjects vigorously, but do not elicit courtship. Moreover, they fail to lower courtship efforts toward females with higher levels of antiaphrodisiac pheromones. At the level of individual pheromones, miR-957 males show a reduced inhibitory response to both 7-Tricosene (7-T) and cis-vaccenyl acetate, with the effect being more pronounced in the case of 7-T. Overall, our results indicate that a single miRNA can contribute to the regulation of complex behaviors, including detection or processing of chemicals that control important survival strategies such as chemical mate-guarding, and the maintenance of sex- and species-specific courtship barriers.

Female Sex Development and Reproductive Duct Formation Depend on Wnt4a in Zebrafish

Wnt4 is a key regulator of ovary development in mammals, but its role in other vertebrates is unknown. Here, Kossack et al. show that zebrafish wnt4a is the ortholog of mammalian Wnt4 and is expressed....

In laboratory strains of zebrafish, sex determination occurs in the absence of a typical sex chromosome and it is not known what regulates the proportion of animals that develop as males or females. Many sex determination and gonad differentiation genes that act downstream of a sex chromosome are well conserved among vertebrates, but studies that test their contribution to this process have mostly been limited to mammalian models. In mammals, WNT4 is a signaling ligand that is essential for ovary and Müllerian duct development, where it antagonizes the male-promoting FGF9 signal. Wnt4 is well conserved across all vertebrates, but it is not known if Wnt4 plays a role in sex determination and/or the differentiation of sex organs in nonmammalian vertebrates. This question is especially interesting in teleosts, such as zebrafish, because they lack an Fgf9 ortholog. Here we show that wnt4a is the ortholog of mammalian Wnt4, and that wnt4b was present in the last common ancestor of humans and zebrafish, but was lost in mammals. We show that wnt4a loss-of-function mutants develop predominantly as males and conclude that wnt4a activity promotes female sex determination and/or differentiation in zebrafish. Additionally, both male and female wnt4a mutants are sterile due to defects in reproductive duct development. Together these results strongly argue that Wnt4a is a conserved regulator of female sex determination and reproductive duct development in mammalian and nonmammalian vertebrates.

Heterochromatin-Enriched Assemblies Reveal the Sequence and Organization of the Drosophila melanogaster Y Chromosome

Heterochromatic regions of the genome are repeat-rich and poor in protein coding genes, and are therefore underrepresented in even the best genome assemblies. One of the most difficult regions of the genome to assemble are sex-limited chromosomes. The Drosophila melanogaster Y chromosome is entirely heterochromatic, yet has wide-ranging effects on male fertility, fitness, and genome-wide gene expression. The genetic basis of this phenotypic variation is difficult to study, in part because we do not know the detailed organization of the Y chromosome. To study Y chromosome organization in D. melanogaster, we develop an assembly strategy involving the in silico enrichment of heterochromatic long single-molecule reads and use these reads to create targeted de novo assemblies of heterochromatic sequences. We assigned contigs to the Y chromosome using Illumina reads to identify male-specific sequences. Our pipeline extends the D. melanogaster reference genome by 11.9 Mb, closes 43.8% of the gaps, and improves overall contiguity. The addition of 10.6 MB of Y-linked sequence permitted us to study the organization of repeats and genes along the Y chromosome. We detected a high rate of duplication to the pericentric regions of the Y chromosome from other regions in the genome. Most of these duplicated genes exist in multiple copies. We detail the evolutionary history of one sex-linked gene family, crystal-Stellate While the Y chromosome does not undergo crossing over, we observed high gene conversion rates within and between members of the crystal-Stellate gene family, Su(Ste), and PCKR, compared to genome-wide estimates. Our results suggest that gene conversion and gene duplication play an important role in the evolution of Y-linked genes.

A Dense Linkage Map of Lake Victoria Cichlids Improved the Pundamilia Genome Assembly and Revealed a Major QTL for Sex-Determination

Genetic linkage maps are essential for comparative genomics, high quality genome sequence assembly and fine scale quantitative trait locus (QTL) mapping. In the present study we identified and genotyped markers via restriction-site associated DNA (RAD) sequencing and constructed a genetic linkage map based on 1,597 SNP markers of an interspecific F2 cross of two closely related Lake Victoria cichlids (Pundamilia pundamilia and P sp. 'red head'). The SNP markers were distributed on 22 linkage groups and the total map size was 1,594 cM with an average marker distance of 1.01 cM. This high-resolution genetic linkage map was used to anchor the scaffolds of the Pundamilia genome and estimate recombination rates along the genome. Via QTL mapping we identified a major QTL for sex in a ∼1.9 Mb region on Pun-LG10, which is homologous to Oreochromis niloticus LG 23 (Ore-LG23) and includes a well-known vertebrate sex-determination gene (amh).

Sex Determination in Ceratopteris richardii Is Accompanied by Transcriptome Changes That Drive Epigenetic Reprogramming of the Young Gametophyte

The fern Ceratopteris richardii is an important model for studies of sex determination and gamete differentiation in homosporous plants. Here we use RNA-seq to de novo assemble a transcriptome and identify genes differentially expressed in young gametophytes as their sex is determined by the presence or absence of the male-inducing pheromone called antheridiogen. Of the 1,163 consensus differentially expressed genes identified, the vast majority (1,030) are up-regulated in gametophytes treated with antheridiogen. GO term enrichment analyses of these DEGs reveals that a large number of genes involved in epigenetic reprogramming of the gametophyte genome are up-regulated by the pheromone. Additional hormone response and development genes are also up-regulated by the pheromone. This C. richardii gametophyte transcriptome and gene expression dataset will prove useful for studies focusing on sex determination and differentiation in plants.

An Association Mapping Framework To Account for Potential Sex Difference in Genetic Architectures

Over the past few years, genome-wide association studies have identified many trait-associated loci that have different effects on females and males, which increased attention to the genetic architecture differences between the sexes. The between-sex differences in genetic architectures can cause a variety of phenomena such as differences in the effect sizes at trait-associated loci, differences in the magnitudes of polygenic background effects, and differences in the phenotypic variances. However, current association testing approaches for dealing with sex, such as including sex as a covariate, cannot fully account for these phenomena and can be suboptimal in statistical power. We present a novel association mapping framework, MetaSex, that can comprehensively account for the genetic architecture differences between the sexes. Through simulations and applications to real data, we show that our framework has superior performance than previous approaches in association mapping.

Identification of General Patterns of Sex-Biased Expression in Daphnia, a Genus with Environmental Sex Determination

Daphnia reproduce by cyclic-parthenogenesis, where phases of asexual reproduction are intermitted by sexual production of diapause stages. This life cycle, together with environmental sex determination, allow the comparison of gene expression between genetically identical males and females. We investigated gene expression differences between males and females in four genotypes of Daphnia magna and compared the results with published data on sex-biased gene expression in two other Daphnia species, each representing one of the major phylogenetic clades within the genus. We found that 42% of all annotated genes showed sex-biased expression in D. magna This proportion is similar both to estimates from other Daphnia species as well as from species with genetic sex determination, suggesting that sex-biased expression is not reduced under environmental sex determination. Among 7453 single copy, one-to-one orthologs in the three Daphnia species, 707 consistently showed sex-biased expression and 675 were biased in the same direction in all three species. Hence these genes represent a core-set of genes with consistent sex-differential expression in the genus. A functional analysis identified that several of them are involved in known sex determination pathways. Moreover, 75% were overexpressed in females rather than males, a pattern that appears to be a general feature of sex-biased gene expression in Daphnia.

Dosage-Dependent Expression Variation Suppressed on the Drosophila Male X Chromosome

DNA copy number variation is associated with many high phenotypic heterogeneity disorders. We systematically examined the impact of Drosophila melanogaster deletions on gene expression profiles to ask whether increased expression variability owing to reduced gene dose might underlie this phenotypic heterogeneity. Indeed, we found that one-dose genes have higher gene expression variability relative to two-dose genes. We then asked whether this increase in variability could be explained by intrinsic noise within cells due to stochastic biochemical events, or whether expression variability is due to extrinsic noise arising from more complex interactions. Our modeling showed that intrinsic gene expression noise averages at the organism level and thus cannot explain increased variation in one-dose gene expression. Interestingly, expression variability was related to the magnitude of expression compensation, suggesting that regulation, induced by gene dose reduction, is noisy. In a remarkable exception to this rule, the single X chromosome of males showed reduced expression variability, even compared with two-dose genes. Analysis of sex-transformed flies indicates that X expression variability is independent of the male differentiation program. Instead, we uncovered a correlation between occupancy of the chromatin-modifying protein encoded by males absent on the first (mof) and expression variability, linking noise suppression to the specialized X chromosome dosage compensation system. MOF occupancy on autosomes in both sexes also lowered transcriptional noise. Our results demonstrate that gene dose reduction can lead to heterogeneous responses, which are often noisy. This has implications for understanding gene network regulatory interactions and phenotypic heterogeneity. Additionally, chromatin modification appears to play a role in dampening transcriptional noise.

Long-Term Fragility of Y Chromosomes Is Dominated by Short-Term Resolution of Sexual Antagonism

The evolution of heteromorphic sex chromosomes has fascinated biologists, inspiring theoretical models, experimental studies, and studies of genome structure. This work has produced a clear model, in which heteromorphic sex chromosomes result from repeated fixations of inversions (or other recombination suppression mechanisms) that tether sexually antagonistic alleles to sex-determining regions, followed by the degeneration of these regions induced by the lack of sex chromosome recombination in the heterogametic sex. However, current models do not predict if inversions are expected to preferentially accumulate on one sex-chromosome or another, and do not address if inversions can accumulate even when they cause difficulties in pairing between heteromorphic chromosomes in the heterogametic sex increasing aneuploidy or meiotic arrest. To address these questions, we developed a population genetic model in which the sex chromosome aneuploidy rate is elevated when males carry an inversion on either the X or Y chromosome. We show that inversions fix more easily when male-beneficial alleles are dominant, and that inversions on the Y chromosome fix with lower selection coefficients than comparable X chromosome inversions. We further show that sex-chromosome inversions can often invade and fix despite causing a substantial increase in the risk of aneuploidy. As sexual antagonism can lead to the fixation of inversions that increase sex chromosomes aneuploidy (which underlies genetic diseases including Klinefelter and Turner syndrome in humans) selection could subsequently favor diverse mechanisms to reduce aneuploidy-including alternative meiotic mechanisms, translocations to, and fusions with, the sex chromosomes, and sex chromosome turnover.

Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

Spermatogenesis is a fundamental process in male reproductive biology and depends on precise balance between self-renewal and differentiation of male germ cells. However, the regulative factors for controlling the balance are poorly understood. In this study, we examined the roles of amh and dmrt1 in male germ cell development by generating their mutants with Crispr/Cas9 technology in zebrafish. Amh mutant zebrafish displayed a female-biased sex ratio, and both male and female amh mutants developed hypertrophic gonads due to uncontrolled proliferation and impaired differentiation of germ cells. A large number of proliferating spermatogonium-like cells were observed within testicular lobules of the amh-mutated testes, and they were demonstrated to be both Vasa- and PH3-positive. Moreover, the average number of Sycp3- and Vasa-positive cells in the amh mutants was significantly lower than in wild-type testes, suggesting a severely impaired differentiation of male germ cells. Conversely, all the dmrt1-mutated testes displayed severe testicular developmental defects and gradual loss of all Vasa-positive germ cells by inhibiting their self-renewal and inducing apoptosis. In addition, several germ cell and Sertoli cell marker genes were significantly downregulated, whereas a prominent increase of Insl3-positive Leydig cells was revealed by immunohistochemical analysis in the disorganized dmrt1-mutated testes. Our data suggest that amh might act as a guardian to control the balance between proliferation and differentiation of male germ cells, whereas dmrt1 might be required for the maintenance, self-renewal, and differentiation of male germ cells. Significantly, this study unravels novel functions of amh gene in fish.

Distinct and Cooperative Roles of amh and dmrt1 in Self-Renewal and Differentiation of Male Germ Cells in Zebrafish

Spermatogenesis is a fundamental process in male reproductive biology and depends on precise balance between self-renewal and differentiation of male germ cells. However, the regulative factors for controlling the balance are poorly understood. In this study, we examined the roles of amh and dmrt1 in male germ cell development by generating their mutants with Crispr/Cas9 technology in zebrafish. Amh mutant zebrafish displayed a female-biased sex ratio, and both male and female amh mutants developed hypertrophic gonads due to uncontrolled proliferation and impaired differentiation of germ cells. A large number of proliferating spermatogonium-like cells were observed within testicular lobules of the amh-mutated testes, and they were demonstrated to be both Vasa- and PH3-positive. Moreover, the average number of Sycp3- and Vasa-positive cells in the amh mutants was significantly lower than in wild-type testes, suggesting a severely impaired differentiation of male germ cells. Conversely, all the dmrt1-mutated testes displayed severe testicular developmental defects and gradual loss of all Vasa-positive germ cells by inhibiting their self-renewal and inducing apoptosis. In addition, several germ cell and Sertoli cell marker genes were significantly downregulated, whereas a prominent increase of Insl3-positive Leydig cells was revealed by immunohistochemical analysis in the disorganized dmrt1-mutated testes. Our data suggest that amh might act as a guardian to control the balance between proliferation and differentiation of male germ cells, whereas dmrt1 might be required for the maintenance, self-renewal, and differentiation of male germ cells. Significantly, this study unravels novel functions of amh gene in fish.

Hill-Robertson Interference Reduces Genetic Diversity on a Young Plant Y-Chromosome

X and Y chromosomes differ in effective population size (N_e ), rates of recombination, and exposure to natural selection, all of which can affect patterns of genetic diversity. On Y chromosomes with suppressed recombination, natural selection is expected to eliminate linked neutral variation, and lower the N_e of Y compared to X chromosomes or autosomes. However, female-biased sex ratios and high variance in male reproductive success can also reduce Y-linked N_e , making it difficult to infer the causes of low Y-diversity. Here, we investigate the factors affecting levels of polymorphism during sex chromosome evolution in the dioecious plant Rumexhastatulus (Polygonaceae). Strikingly, we find that neutral diversity for genes on the Y chromosome is, on average, 2.1% of the value for their X-linked homologs, corresponding to a chromosome-wide reduction of 93% compared to the standard neutral expectation. We demonstrate that the magnitude of this diversity loss is inconsistent with reduced male N_e caused by neutral processes. Instead, using forward simulations and estimates of the distribution of deleterious fitness effects, we show that Y chromosome diversity loss can be explained by purifying selection acting in aggregate over a large number of genetically linked sites. Simulations also suggest that our observed level of Y-diversity is consistent with the joint action of purifying and positive selection, but only for models in which there were fewer constrained sites than we empirically estimated. Given the relatively recent origin of R. hastatulus sex chromosomes, our results imply that Y-chromosome degeneration in the early stages may be largely driven by selective interference rather than by neutral genetic drift of silenced Y-linked genes.

Sex Chromosome Evolution, Heterochiasmy, and Physiological QTL in the Salmonid Brook Charr Salvelinus fontinalis

Whole-genome duplication (WGD) can have large impacts on genome evolution, and much remains unknown about these impacts. This includes the mechanisms of coping with a duplicated sex determination system and whether this has an impact on increasing the diversity of sex determination mechanisms. Other impacts include sexual conflict, where alleles having different optimums in each sex can result in sequestration of genes into nonrecombining sex chromosomes. Sex chromosome development itself may involve sex-specific recombination rate (i.e., heterochiasmy), which is also poorly understood. The family Salmonidae is a model system for these phenomena, having undergone autotetraploidization and subsequent rediploidization in most of the genome at the base of the lineage. The salmonid master sex determining gene is known, and many species have nonhomologous sex chromosomes, putatively due to transposition of this gene. In this study, we identify the sex chromosome of Brook Charr Salvelinus fontinalis and compare sex chromosome identities across the lineage (eight species and four genera). Although nonhomology is frequent, homologous sex chromosomes and other consistencies are present in distantly related species, indicating probable convergence on specific sex and neo-sex chromosomes. We also characterize strong heterochiasmy with 2.7-fold more crossovers in maternal than paternal haplotypes with paternal crossovers biased to chromosome ends. When considering only rediploidized chromosomes, the overall heterochiasmy trend remains, although with only 1.9-fold more recombination in the female than the male. Y chromosome crossovers are restricted to a single end of the chromosome, and this chromosome contains a large interspecific inversion, although its status between males and females remains unknown. Finally, we identify quantitative trait loci (QTL) for 21 unique growth, reproductive, and stress-related phenotypes to improve knowledge of the genetic architecture of these traits important to aquaculture and evolution.

Rapid Evolution of Ovarian-Biased Genes in the Yellow Fever Mosquito (Aedes aegypti)

Males and females exhibit highly dimorphic phenotypes, particularly in their gonads, which is believed to be driven largely by differential gene expression. Typically, the protein sequences of genes upregulated in males, or male-biased genes, evolve rapidly as compared to female-biased and unbiased genes. To date, the specific study of gonad-biased genes remains uncommon in metazoans. Here, we identified and studied a total of 2927, 2013, and 4449 coding sequences (CDS) with ovary-biased, testis-biased, and unbiased expression, respectively, in the yellow fever mosquito Aedes aegypti The results showed that ovary-biased and unbiased CDS had higher nonsynonymous to synonymous substitution rates (dN/dS) and lower optimal codon usage (those codons that promote efficient translation) than testis-biased genes. Further, we observed higher dN/dS in ovary-biased genes than in testis-biased genes, even for genes coexpressed in nonsexual (embryo) tissues. Ovary-specific genes evolved exceptionally fast, as compared to testis- or embryo-specific genes, and exhibited higher frequency of positive selection. Genes with ovary expression were preferentially involved in olfactory binding and reception. We hypothesize that at least two potential mechanisms could explain rapid evolution of ovary-biased genes in this mosquito: (1) the evolutionary rate of ovary-biased genes may be accelerated by sexual selection (including female-female competition or male-mate choice) affecting olfactory genes during female swarming by males, and/or by adaptive evolution of olfactory signaling within the female reproductive system (e.g., sperm-ovary signaling); and/or (2) testis-biased genes may exhibit decelerated evolutionary rates due to the formation of mating plugs in the female after copulation, which limits male-male sperm competition.

Comparative Genomics of the Ectomycorrhizal Sister Species Rhizopogon vinicolor and Rhizopogon vesiculosus (Basidiomycota: Boletales) Reveals a Divergence of the Mating Type B Locus

Divergence of breeding system plays an important role in fungal speciation. Ectomycorrhizal fungi, however, pose a challenge for the study of reproductive biology because most cannot be mated under laboratory conditions. To overcome this barrier, we sequenced the draft genomes of the ectomycorrhizal sister species Rhizopogon vinicolor Smith and Zeller and R. vesiculosus Smith and Zeller (Basidiomycota, Boletales)-the first genomes available for Basidiomycota truffles-and characterized gene content and organization surrounding their mating type loci. Both species possess a pair of homeodomain transcription factor homologs at the mating type A-locus as well as pheromone receptor and pheromone precursor homologs at the mating type B-locus. Comparison of Rhizopogon genomes with genomes from Boletales, Agaricales, and Polyporales revealed synteny of the A-locus region within Boletales, but several genomic rearrangements across orders. Our findings suggest correlation between gene content at the B-locus region and breeding system in Boletales with tetrapolar species possessing more diverse gene content than bipolar species. Rhizopogon vinicolor possesses a greater number of B-locus pheromone receptor and precursor genes than R. vesiculosus, as well as a pair of isoprenyl cysteine methyltransferase genes flanking the B-locus compared to a single copy in R. vesiculosus Examination of dikaryotic single nucleotide polymorphisms within genomes revealed greater heterozygosity in R. vinicolor, consistent with increased rates of outcrossing. Both species possess the components of a heterothallic breeding system with R. vinicolor possessing a B-locus region structure consistent with tetrapolar Boletales and R. vesiculosus possessing a B-locus region structure intermediate between bipolar and tetrapolar Boletales.

Control of a Novel Spermatocyte-Promoting Factor by the Male Germline Sex Determination Factor PHF7 of Drosophila melanogaster

A key aspect of germ cell development is to establish germline sexual identity and initiate a sex-specific developmental program to promote spermatogenesis or oogenesis. Previously, we have identified the histone reader Plant Homeodomain Finger 7 (PHF7) as an important regulator of male germline identity. To understand how PHF7 directs sexual differentiation of the male germline, we investigated the downstream targets of PHF7 by combining transcriptome analyses, which reveal genes regulated by Phf7, with genomic profiling of histone H3K4me2, the chromatin mark that is bound by PHF7. Through these genomic experiments, we identify a novel spermatocyte factor Receptor Accessory Protein Like 1 (REEPL1) that can promote spermatogenesis and whose expression is kept off by PHF7 in the spermatogonial stage. Loss of Reepl1 significantly rescues the spermatogenesis defects in Phf7 mutants, indicating that regulation of Reepl1 is an essential aspect of PHF7 function. Further, increasing REEPL1 expression facilitates spermatogenic differentiation. These results indicate that PHF7 controls spermatogenesis by regulating the expression patterns of important male germline genes.

Revisiting Suppression of Interspecies Hybrid Male Lethality in Caenorhabditis Nematodes

Within the nematode genus Caenorhabditis, Caenorhabditis briggsae and C. nigoni are among the most closely related species known. They differ in sexual mode, with C. nigoni retaining the ancestral XO male-XX female outcrossing system, while C. briggsae recently evolved self-fertility and an XX-biased sex ratio. Wild-type C. briggsae and C. nigoni can produce fertile hybrid XX female progeny, but XO progeny are either 100% inviable (when C. briggsae is the mother) or viable but sterile (when C. nigoni is the mother). A recent study provided evidence suggesting that loss of the Cbr-him-8 meiotic regulator in C. briggsae hermaphrodites allowed them to produce viable and fertile hybrid XO male progeny when mated to C. nigoni Because such males would be useful for a variety of genetic experiments, we sought to verify this result. Preliminary crosses with wild-type C. briggsae hermaphrodites occasionally produced fertile males, but they could not be confirmed to be interspecies hybrids. Using an RNA interference (RNAi) protocol that eliminates any possibility of self-progeny in Cbr-him-8 hermaphrodites, we found sterile males bearing the C. nigoni X chromosome, but no fertile males bearing the C. briggsae X, as in wild-type crosses. Our results suggest that the apparent rescue of XO hybrid viability and fertility is due to incomplete purging of self-sperm prior to mating.

Advances in Understanding Mating Type Gene Organization in the Mushroom-Forming Fungus Flammulina velutipes

The initiation of sexual development in the important edible and medicinal mushroom Flammulina velutipes is controlled by special genes at two different, independent, mating type (MAT) loci: HD and PR. We expanded our understanding of the F. velutipes mating type system by analyzing the MAT loci from a series of strains. The HD locus of F. velutipes houses homeodomain genes (Hd genes) on two separated locations: sublocus HD-a and HD-b. The HD-b subloci contained strain-specific Hd1/Hd2 gene pairs, and crosses between strains with different HD-b subloci indicated a role in mating. The function of the HD-a sublocus remained undecided. Many, but not all strains contained the same conserved Hd2 gene at the HD-a sublocus. The HD locus usually segregated as a whole, though we did detect one new HD locus with a HD-a sublocus from one parental strain, and a HD-b sublocus from the other. The PR locus of F. velutipes contained pheromone receptor (STE3) and pheromone precursor (Pp) genes at two locations, sublocus PR-a and PR-b. PR-a and PR-b both contained sets of strain-specific STE3 and Pp genes, indicating a role in mating. PR-a and PR-b cosegregated in our experiments. However, the identification of additional strains with identical PR-a, yet different PR-b subloci, demonstrated that PR subloci can recombine within the PR locus. In conclusion, at least three of the four MAT subloci seem to participate in mating, and new HD and PR loci can be generated through intralocus recombination in F. velutipes.

To Break or Not To Break: Sex Chromosome Hemizygosity During Meiosis in Caenorhabditis

Meiotic recombination establishes connections between homologous chromosomes to promote segregation. Hemizygous regions of sex chromosomes have no homologous chromosome to recombine with, yet must be transmitted through meiosis. An extreme case of hemizygosity exists in the genus Caenorhabditis, where males have a single X chromosome that completely lacks a homologous partner. To determine whether similar strategies have evolved to accommodate hemizygosity of the X during male meiosis in Caenorhabditis with distinct modes of sexual reproduction, we examined induction and processing of meiotic double strand breaks (DSBs) in androdioecious (hermaphrodite/male) Caenorhabditis elegans and C. briggsae, and gonochoristic (female/male) C. remanei and C. brenneri Analysis of the recombinase RAD-51 suggests more meiotic DSBs are induced in gonochoristic vs. androdioecious species. However, in late prophase in all species, chromosome pairs are restructured into bivalents around a single axis, suggesting that the holocentric nature of Caenorhabditis chromosomes dictates a single crossover per bivalent regardless of the number of DSBs induced. Interestingly, RAD-51 foci were readily observed on the X chromosome of androdioecious male germ cells, while very few were detected in gonochoristic male germ cells. As in C. elegans, the X chromosome in C. briggsae male germ cells undergoes transient pseudosynapsis and flexibility in DSB repair pathway choice. In contrast, in C. remanei and C. brenneri male germ cells, the X chromosome does not undergo pseudosynapsis and appears refractory to SPO-11-induced breaks. Together our results suggest that distinct strategies have evolved to accommodate sex chromosome hemizygosity during meiosis in closely related Caenorhabditis species.

Sex Differences in Drosophila Somatic Gene Expression: Variation and Regulation by doublesex

Sex differences in gene expression have been widely studied in Drosophila melanogaster. Sex differences vary across strains, but many molecular studies focus on only a single strain, or on genes that show sexually dimorphic expression in many strains. How extensive variability is and whether this variability occurs among genes regulated by sex determination hierarchy terminal transcription factors is unknown. To address these questions, we examine differences in sexually dimorphic gene expression between two strains in Drosophila adult head tissues. We also examine gene expression in doublesex (dsx) mutant strains to determine which sex-differentially expressed genes are regulated by DSX, and the mode by which DSX regulates expression. We find substantial variation in sex-differential expression. The sets of genes with sexually dimorphic expression in each strain show little overlap. The prevalence of different DSX regulatory modes also varies between the two strains. Neither the patterns of DSX DNA occupancy, nor mode of DSX regulation explain why some genes show consistent sex-differential expression across strains. We find that the genes identified as regulated by DSX in this study are enriched with known sites of DSX DNA occupancy. Finally, we find that sex-differentially expressed genes and genes regulated by DSX are highly enriched on the fourth chromosome. These results provide insights into a more complete pool of potential DSX targets, as well as revealing the molecular flexibility of DSX regulation.

Sequence of the Gonium pectorale Mating Locus Reveals a Complex and Dynamic History of Changes in Volvocine Algal Mating Haplotypes

Sex-determining regions (SDRs) or mating-type (MT) loci in two sequenced volvocine algal species, Chlamydomonas reinhardtii and Volvox carteri, exhibit major differences in size, structure, gene content, and gametolog differentiation. Understanding the origin of these differences requires investigation of MT loci from related species. Here, we determined the sequences of the minus and plus MT haplotypes of the isogamous 16-celled volvocine alga, Gonium pectorale, which is more closely related to the multicellular V. carteri than to C. reinhardtii. Compared to C. reinhardtii MT, G. pectorale MT is moderately larger in size, and has a less complex structure, with only two major syntenic blocs of collinear gametologs. However, the gametolog content of G. pectorale MT has more overlap with that of V. carteri MT than with C. reinhardtii MT, while the allelic divergence between gametologs in G. pectorale is even lower than that in C. reinhardtii. Three key sex-related genes are conserved in G. pectorale MT: GpMID and GpMTD1 in MT–, and GpFUS1 in MT+. GpFUS1 protein exhibited specific localization at the plus-gametic mating structure, indicating a conserved function in fertilization. Our results suggest that the G. pectorale–V. carteri common ancestral MT experienced at least one major reformation after the split from C. reinhardtii, and that the V. carteri ancestral MT underwent a subsequent expansion and loss of recombination after the divergence from G. pectorale. These data begin to polarize important changes that occurred in volvocine MT loci, and highlight the potential for discontinuous and dynamic evolution in SDRs.

Primary Sex Determination in Drosophila melanogaster Does Not Rely on the Male-Specific Lethal Complex

It has been proposed that the Male Specific Lethal (MSL) complex is active in Drosophila melanogaster embryos of both sexes prior to the maternal-to-zygotic transition. Elevated gene expression from the two X chromosomes of female embryos is proposed to facilitate the stable establishment of Sex-lethal (Sxl) expression, which determines sex and represses further activity of the MSL complex, leaving it active only in males. Important supporting data included female-lethal genetic interactions between the seven msl genes and either Sxl or scute and sisterlessA, two of the X-signal elements (XSE) that regulate early Sxl expression. Here I report contrary findings that there are no female-lethal genetic interactions between the msl genes and Sxl or its XSE regulators. Fly stocks containing the msl3(1) allele were found to exhibit a maternal-effect interaction with Sxl, scute, and sisterlessA mutations, but genetic complementation experiments showed that msl3 is neither necessary nor sufficient for the female-lethal interactions, which appear to be due to an unidentified maternal regulator of Sxl. Published data cited as evidence for an early function of the MSL complex in females, including a maternal effect of msl2, have been reevaluated and found not to support a maternal, or other effect, of the MSL complex in sex determination. These findings suggest that the MSL complex is not involved in primary sex determination or in X chromosome dosage compensation prior to the maternal-to-zygotic transition.

An Examination of Dynamic Gene Expression Changes in the Mouse Brain During Pregnancy and the Postpartum Period. G3 (Bethesda). 2015;6:221-233. [PMID: 26596646 PMCID: PMC4704721 DOI: 10.1534/g3.115.020982]

The developmental transition to motherhood requires gene expression changes that alter the brain to drive the female to perform maternal behaviors. We broadly examined the global transcriptional response in the mouse maternal brain, by examining four brain regions: hypothalamus, hippocampus, neocortex, and cerebellum, in virgin females, two pregnancy time points, and three postpartum time points. We find that overall there are hundreds of differentially expressed genes, but each brain region and time point shows a unique molecular signature, with only 49 genes differentially expressed in all four regions. Interestingly, a set of “early-response genes” is repressed in all brain regions during pregnancy and postpartum stages. Several genes previously implicated in underlying postpartum depression change expression. This study serves as an atlas of gene expression changes in the maternal brain, with the results demonstrating that pregnancy, parturition, and postpartum maternal experience substantially impact diverse brain regions.

A new dataset of spermatogenic vs. oogenic transcriptomes in the nematode Caenorhabditis elegans

The nematode Caenorhabditis elegans is an important model for studies of germ cell biology, including the meiotic cell cycle, gamete specification as sperm or oocyte, and gamete development. Fundamental to those studies is a genome-level knowledge of the germline transcriptome. Here, we use RNA-Seq to identify genes expressed in isolated XX gonads, which are approximately 95% germline and 5% somatic gonadal tissue. We generate data from mutants making either sperm [fem-3(q96)] or oocytes [fog-2(q71)], both grown at 22°. Our dataset identifies a total of 10,754 mRNAs in the polyadenylated transcriptome of XX gonads, with 2748 enriched in spermatogenic gonads, 1732 enriched in oogenic gonads, and the remaining 6274 not enriched in either. These spermatogenic, oogenic, and gender-neutral gene datasets compare well with those of previous studies, but double the number of genes identified. A comparison of the additional genes found in our study with in situ hybridization patterns in the Kohara database suggests that most are expressed in the germline. We also query our RNA-Seq data for differential exon usage and find 351 mRNAs with sex-enriched isoforms. We suggest that this new dataset will prove useful for studies focusing on C. elegans germ cell biology.

Molecular characterization of a heterothallic mating system in Pseudogymnoascus destructans, the Fungus causing white-nose syndrome of bats

White-nose syndrome (WNS) of bats has devastated bat populations in eastern North America since its discovery in 2006. WNS, caused by the fungus Pseudogymnoascus destructans, has spread quickly in North America and has become one of the most severe wildlife epidemics of our time. While P. destructans is spreading rapidly in North America, nothing is known about the sexual capacity of this fungus. To gain insight into the genes involved in sexual reproduction, we characterized the mating-type locus (MAT) of two Pseudogymnoascus spp. that are closely related to P. destructans and homothallic (self-fertile). As with other homothallic Ascomycota, the MAT locus of these two species encodes a conserved α-box protein (MAT1-1-1) as well as two high-mobility group (HMG) box proteins (MAT1-1-3 and MAT1-2-1). Comparisons with the MAT locus of the North American isolate of P. destructans (the ex-type isolate) revealed that this isolate of P. destructans was missing a clear homolog of the conserved HMG box protein (MAT1-2-1). These data prompted the discovery and molecular characterization of a heterothallic mating system in isolates of P. destructans from the Czech Republic. Both mating types of P. destructans were found to coexist within hibernacula, suggesting the presence of mating populations in Europe. Although populations of P. destructans in North America are thought to be clonal and of one mating type, the potential for sexual recombination indicates that continued vigilance is needed regarding introductions of additional isolates of this pathogen.