The role of sex chromosomes in sexual dimorphism: discordance between molecular and phenotypic data

Genome-Wide Identification and Characterization of gh/prl/sl Family in Cynoglossus semilaevis

The Chinese tongue sole (Cynoglossus semilaevis) is a marine flatfish of significant economic value, characterized by pronounced female-biased sexual size dimorphism (SSD). Sexual differences of cell number and gene expression within the PIT-1 lineage of the pituitary gland may be crucial for interpreting the female-biased SSD of C. semilaevis. Among hormones secreted by PIT-1 cell lineage, growth hormone (gh), prolactin (prl), prolactin 2 (prl2), and somatolactin (sl) comprise a gene family within the extensive superfamily of class-1 helical cytokines. To better understand the function of the gh/prl/sl in teleost SSD, we firstly identified five genes of the gh/prl/sl family (gh, sl, prl, prl2a, and prl2b) and their receptors (ghra, ghrb, prlra, prlrb, and prlr-like) from C. semilaevis at the genome-wide level. Phylogenetic analyses revealed that the gh/prl/sl family and their receptors were each clustered into five distinct groups. More microsatellites were revealed in the intron 2 of gh gene of female rather than the male and pseudo-male individuals, which is positively correlated with its sexual expression pattern. Interaction network prediction indicated that gh, prl, and sl may collectively contribute to individual growth and development. A FRET experiment showed that ghra can act as a receptor for sl. Additionally, the transcripts of the gh/prl/sl family and their receptors exhibited varying abundances in the pituitary, brain, gonad, and liver of both female and male C. semilaevis, with most ligands showing the highest abundance in the female pituitary. Furthermore, gh and sl were found to be maternally expressed. The knock-down of gh, prl, and sl in the pituitary cells could lead to the expression change of igf1, c-fos, and sos2. This study provided a foundation for further functional characterization of the gh/prl/sl gene family, contributing to a deeper understanding of the growth and reproductive mechanisms in C. semilaevis.

Enrichment of hard sweeps on the X chromosome compared to autosomes in six Drosophila species

The X chromosome, being hemizygous in males, is exposed one-third of the time increasing the visibility of new mutations to natural selection, potentially leading to different evolutionary dynamics than autosomes. Recently, we found an enrichment of hard selective sweeps over soft selective sweeps on the X chromosome relative to the autosomes in a North American population of Drosophila melanogaster. To understand whether this enrichment is a universal feature of evolution on the X chromosome, we analyze diversity patterns across 6 commonly studied Drosophila species. We find an increased proportion of regions with steep reductions in diversity and elevated homozygosity on the X chromosome compared to autosomes. To assess if these signatures are consistent with positive selection, we simulate a wide variety of evolutionary scenarios spanning variations in demography, mutation rate, recombination rate, background selection, hard sweeps, and soft sweeps and find that the diversity patterns observed on the X are most consistent with hard sweeps. Our findings highlight the importance of sex chromosomes in driving evolutionary processes and suggest that hard sweeps have played a significant role in shaping diversity patterns on the X chromosome across multiple Drosophila species.

Evolution and expression patterns of the neo-sex chromosomes of the crested ibis

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

Enrichment of hard sweeps on the X chromosome compared to autosomes in six Drosophila species

The X chromosome, being hemizygous in males, is exposed one third of the time increasing the visibility of new mutations to natural selection, potentially leading to different evolutionary dynamics than autosomes. Recently, we found an enrichment of hard selective sweeps over soft selective sweeps on the X chromosome relative to the autosomes in a North American population of Drosophila melanogaster. To understand whether this enrichment is a universal feature of evolution on the X chromosome, we analyze diversity patterns across six commonly studied Drosophila species. We find an increased proportion of regions with steep reductions in diversity and elevated homozygosity on the X chromosome compared to autosomes. To assess if these signatures are consistent with positive selection, we simulate a wide variety of evolutionary scenarios spanning variations in demography, mutation rate, recombination rate, background selection, hard sweeps, and soft sweeps, and find that the diversity patterns observed on the X are most consistent with hard sweeps. Our findings highlight the importance of sex chromosomes in driving evolutionary processes and suggest that hard sweeps have played a significant role in shaping diversity patterns on the X chromosome across multiple Drosophila species.

Biometric indices, growth pattern, and physiological status of captive-reared indigenous Yellowtail brood catfish, Pangasius pangasius (Hamilton, 1822)

The current study was conducted to comprehend the variation of the growth pattern and the physiological status of Pangasius pangasius (Hamilton, 1822) reared in the pond from July 2021 to September 2021. A total of 90 brood individuals were collected from the Meghna River and studied for the present experiment. The overall growth pattern for P. pangasius was isometric (b = 3.00) in the Meghna River, while the male fish showed positive (b > 3.00) and the female showed a negative allometric growth pattern (b < 3.00). The Fulton's condition factor (KF) value was greater than 1, indicating a healthy population and a habitat with abundant food sources. Moreover, the KF value was found to be substantially related to total body mass. On the other hand, both sexes of P. pangasius had relative weight average values of more than 100, simply indicates the fish is naturally obese and has enough stored energy to maintain its physiology. The calculated form factors suggested an elongated body shape, which is typical of many riverine fishes. In addition, a small number of morphological traits varied significantly in the current study. Likewise, in the case of morphometric features, the principal component analysis showed a considerable degree of connectedness between the male and female individuals. The blood values showed no significant difference between the sexes. This might result from feeding fish the same kind of food and maintaining them in the same environment. However, the higher temperature might have caused slight blood abnormalities in both sexes. The current findings strongly support raising this fish in captivity and would offer useful information for fish farm owners, entrepreneurs, stakeholders, and other concerns in Bangladesh and surrounding nations.

DNA methylation is associated with codon degeneracy in a species of bumblebee

Social insects display extreme phenotypic differences between sexes and castes even though the underlying genome can be almost identical. Epigenetic processes have been proposed as a possible mechanism for mediating these phenotypic differences. Using whole genome bisulfite sequencing of queens, males, and reproductive female workers we have characterised the sex- and caste-specific methylome of the bumblebee Bombus terrestris. We have identified a potential role for DNA methylation in histone modification processes which may influence sex and caste phenotypic differences. We also find differentially methylated genes generally show low levels of DNA methylation which may suggest a separate function for lowly methylated genes in mediating transcriptional plasticity, unlike highly methylated genes which are usually involved in housekeeping functions. We also examined the relationship between the underlying genome and the methylome using whole genome re-sequencing of the same queens and males. We find DNA methylation is enriched at zero-fold degenerate sites. We suggest DNA methylation may be acting as a targeted mutagen at these sites, providing substrate for selection via non-synonymous changes in the underlying genome. However, we did not see any relationship between DNA methylation and rates of positive selection in our samples. In order to fully assess a possible role for DNA methylation in adaptive processes a specifically designed study using natural population data is needed.

Gene-poor Y-chromosomes substantially impact male trait heritabilities and may help shape sexually dimorphic evolution

How natural selection facilitates sexually dimorphic evolution despite a shared genome is unclear. The patrilineal inheritance of Y-chromosomes makes them an appealing solution. However, they have largely been dismissed due to their gene-poor, heterochromatic nature and because the additive genetic variation necessary for adaptive evolution is theoretically difficult to maintain. Further, previous empirical work has revealed mostly Y-linked sign epistatic variance segregating within populations, which can often impede adaptive evolution. To assess the evolutionary impact of Y-linked variation, we established replicate populations in Drosophila simulans containing multiple Y-chromosomes (YN populations) or a single Y-chromosome variant (Y1 populations) drawn from a single population. We estimated male and female heritabilities for several traits known to be influenced by Y-chromosomes, including the number of sternopleural bristles, abdominal bristles, sex comb teeth, and tibia length. A decrease in YN heritabilities compared with Y1 would be consistent with Y-chromosome variation being sign epistatic. A decrease in Y1 heritabilities would be consistent with Y-chromosome variation being additive, though additive-by-additive epistatic variation cannot be entirely dismissed. Female heritability estimates served as controls and were not expected to differ. We found male Y1 populations exhibited lower heritabilities for all traits except tibia length; consistent with Y-linked additivity (on average YN trait heritabilities were 25% greater than Y1). Female estimates showed no difference. These data suggest Y-chromosomes should play an important role in male trait evolution and may even influence sexually dimorphic evolution by shaping traits shared by both sexes.

Transcriptome analysis of differentially expressed circular RNAs in the testis and ovary of golden pompano (Trachinotus blochii)

The artificial breeding of golden pompano (Trachinotus blochii) has expanded greatly in recent years, and after long-term breeding efforts, clear sexual dimorphisms have been observed in T. blochii growth traits, with females growing faster. As sponges of microRNA (miRNAs), circular RNAs (CircRNAs) can alleviate miRNA inhibition of target mRNA. However, few studies have examined sex-related CircRNAs and none of those have looked at T. blochii. To further understand the role of CircRNAs in sex differentiation and sexual size dimorphism in T. blochii, six CircRNA libraries were constructed from the testes and ovaries of T. blochii. A total of 1522 CircRNAs were found distributed over all 24 chromosomes of T. blochii. 135 differentially expressed CircRNAs (DECs) were identified by screening, These DECs were then subjected to GO enrichment, which found 47 enriched pathways. A number of CircRNAs were enriched in cellular processes and metabolic processes. According to the KEGG pathway analysis, a series of sex differentiation pathways were enriched, including the GnRH, calcium, and MAPK signaling pathways. Furthermore, we selected two CircRNAs from the DECs named circ-cacna1b and circ-octc. We found that the cacna1b gene is regulated by 7 miRNAs, 3 of which were regulated by circ-cacna1b, i.e., mmu-miR-138-5p, fru-miR-138, and pma-miR-138b. In addition, the miRNA named pma-miR-138b can regulate sex-related genes, such as sox9 and dmrt1, among others. The co-expression network of CircRNA-miRNA-mRNA showed circ-cacna1b may play a crucial role in T. blochii sex differentiation by regulating pma-miR-138b to affect the expression of sex differentiation genes. The circ-octc may be one of the largest contributors to sexual size dimorphism during growth through its effect on lipid metabolism. These findings could broaden our understanding of CircRNAs and provide new insight into their function in sex differentiation and growth.

Enrichment of Hard Sweeps on the X Chromosome in Drosophila melanogaster

The characteristic properties of the X chromosome, such as male hemizygosity and its unique inheritance pattern, expose it to natural selection in a way that can be different from the autosomes. Here, we investigate the differences in the tempo and mode of adaptation on the X chromosome and autosomes in a population of Drosophila melanogaster. Specifically, we test the hypothesis that due to hemizygosity and a lower effective population size on the X, the relative proportion of hard sweeps, which are expected when adaptation is gradual, compared with soft sweeps, which are expected when adaptation is rapid, is greater on the X than on the autosomes. We quantify the incidence of hard versus soft sweeps in North American D. melanogaster population genomic data with haplotype homozygosity statistics and find an enrichment of the proportion of hard versus soft sweeps on the X chromosome compared with the autosomes, confirming predictions we make from simulations. Understanding these differences may enable a deeper understanding of how important phenotypes arise as well as the impact of fundamental evolutionary parameters on adaptation, such as dominance, sex-specific selection, and sex-biased demography.

Sex hormones, intestinal inflammation, and the gut microbiome: Major influencers of the sexual dimorphisms in obesity

Obesity is defined as the excessive accumulation of body fat and is associated with an increased risk of developing major health problems such as cardiovascular disease, diabetes and stroke. There are clear sexual dimorphisms in the epidemiology, pathophysiology and sequelae of obesity and its accompanying metabolic disorders, with females often better protected compared to males. This protection has predominantly been attributed to the female sex hormone estrogen and differences in fat distribution. More recently, the sexual dimorphisms of obesity have also been attributed to the differences in the composition and function of the gut microbiota, and the intestinal immune system. This review will comprehensively summarize the pre-clinical and clinical evidence for these sexual dimorphisms and discuss the interplay between sex hormones, intestinal inflammation and the gut microbiome in obesity. Major gaps and limitations of this rapidly growing area of research will also be highlighted in this review.

Female-limited X chromosome evolution reveals that lifespan is mainly modulated by interlocus rather than intralocus sexual conflict

Abstract

Sexual dimorphism in somatic investment may be shaped by two distinct forms of sexual conflict; under intralocus sexual conflict (IASC), males and females have different optimal levels of somatic investment but are constrained from reaching their respective optima by their shared genome, while under interlocus sexual conflict (IRSC), males and females have different optimal sexual strategies, which could have direct or indirect effects on levels of somatic investment. We investigated effects of IASC and IRSC on two aspects of somatic investment, immune defence strategies and longevity, using previously established female-limited experimental evolution lines in Drosophila melanogaster. We found little evidence for any effect of either type of sexual conflict on investment in the immune defence resistance or tolerance. Nor did we find convincing evidence that longevity is subject to IASC in this species. However, we did find evidence that increased female control over mating rate had important and opposite effects on longevity between the sexes. Specifically, females that had adapted to high levels of female control over mating had a longer lifespan when kept in mixed-sex groups, while males had shorter longevity, perhaps due to increased investment in post-copulatory sexual selection. These novel results show that female control over mating rates may have important and unexpected effects on patterns of somatic investment.

Significance statement

Sexual conflict occurs between the two sexes over numerous life history traits, and it is complex to disentangle how these traits interact and affect each other. Here we use a long-term evolution experiment to investigate sexual dimorphism in somatic maintenance. We found no effect of feminising the X chromosome on female immune defence. However, we did find that increased female control over mating rate resulted in longer female lifespan, but reduced male lifespan, and that these effects were dependent on social context (isolated or in mixed-sex groups). Unlike previous studies on the effect of sexual conflict on longevity, our experiment did not manipulate environmental conditions nor the adult sex ratio, which is likely to reduce both pre- and post-copulatory sexual selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00265-022-03231-4.

X chromosome agents of sexual differentiation

Understanding sex differences in physiology and disease requires the identification of the molecular agents that cause phenotypic sex differences. Two groups of such agents are genes located on the sex chromosomes, and gonadal hormones. The former have coherent linkage to chromosomes that form differently in the two sexes under the influence of genomic forces that are not related to reproductive function, whereas the latter have a direct or indirect relationship to reproduction. Evidence published in the past 5 years supports the identification of several agents of sexual differentiation encoded by the X chromosome in mice, including Kdm5c, Kdm6a, Ogt and Xist. These X chromosome agents have wide pleiotropic effects, potentially influencing sex differences in many different tissues, a characteristic shared with the gonadal hormones. The identification of X chromosome agents of sexual differentiation will facilitate understanding of complex intersecting gene pathways underlying sex differences in disease.

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.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Direct Full-Length RNA Sequencing Reveals an Important Role of Epigenetics During Sexual Reversal in Chinese Soft-Shelled Turtle

Sex dimorphism is a key feature of Chinese soft-shelled turtle (Pelodiscus sinensis). The males (M) have higher econosmic value than females (F) due to wider calipash and faster growth. Exogenous hormones like estradiol and methyltestosterone can induce sexual reversal to form new phenotypes (pseudo-female, PF; pseudo-male, PM) without changing the genotype. The possibility of inducing sexual reversal is particularly important in aquaculture breeding, but the underlying biological mechanisms remain unclear. Here we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in P. sinensis. Nanopore sequencing of the four gender types (M, F, PF, and PM) showed that the distribution of read length and gene expression was more similar between same-sex phenotypes than same-sex genotypes. Compared to turtles with an M phenotype, alternative splicing was more pronounced in F turtles, especially at alternative 3′ splice sites, alternative 5′ splice sites, and alternative first exons. Furthermore, the two RNA methylation modifications m5C and m6A were differentially distributed across gender phenotypes, with the M type having more modification sites in coding sequence regions, but fewer modification sites in 3′UTR regions. Quantitative analysis of enriched m6A RNAs revealed that the N6-methylated levels of Odf2, Pacs2, and Ak1 were significantly higher in M phenotype individuals, while the N6-methylated levels of Ube2o were reduced after sexual reversal from both M and F phenotypes. Taken together, these findings reveal an important role of epigenetics during sexual reversal in Chinese soft-shelled turtles.

Experimental evolution reveals sex-specific dominance for surviving bacterial infection in laboratory populations of Drosophila melanogaster

Males and females are subjected to distinct kinds of selection pressures, often leading to the evolution of sex‐specific genetic architecture, an example being sex‐specific dominance. Sex‐specific dominance reversals (SSDRs), where alleles at sexually antagonistic loci are at least partially dominant in the sex they benefit, have been documented in Atlantic salmon, rainbow trout, and seed beetles. Another interesting feature of many sexually reproducing organisms is the asymmetric inheritance pattern of X chromosomes, which often leads to distinct evolutionary outcomes on X chromosomes compared to autosomes. Examples include the higher efficacy of sexually concordant selection on X chromosomes, and X chromosomes being more conducive to the maintenance of sexually antagonistic polymorphisms under certain conditions. Immunocompetence is a trait that has been extensively investigated for sexual dimorphism with growing evidence for sex‐specific or sexually antagonistic variation. X chromosomes have been shown to harbor substantial immunity‐related genetic variation in the fruit fly, Drosophila melanogaster. Here, using interpopulation crosses and cytogenetic cloning, we investigated sex‐specific dominance and the role of the X chromosome in improved postinfection survivorship of laboratory populations of D. melanogaster selected against pathogenic challenge by Pseudomonas entomophila. We could not detect any contribution of the X chromosome to the evolved immunocompetence of our selected populations, as well as to within‐population variation in immunocompetence. However, we found strong evidence of sex‐specific dominance related to surviving bacterial infection. Our results indicate that alleles that confer a survival advantage to the selected populations are, on average, partially dominant in females but partially recessive in males. This could also imply an SSDR for overall fitness, given the putative evidence for sexually antagonistic selection affecting immunocompetence in Drosophila melanogaster. We also highlight sex‐specific dominance as a potential mechanism of sex differences in immunocompetence, with population‐level sex differences primarily driven by sex differences in heterozygotes.

Disentangling Verbal Arguments: Intralocus Sexual Conflict in Haplodiploids

AbstractIn haplodiploids, (1) alleles spend twice as many generations in females as in males, (2) males are never heterozygous and therefore express recessive alleles, and (3) males sire daughters but not sons. Intralocus sexual conflict therefore operates differently in haplodiploids than in diploids and shares strong similarities with loci on X (or Z) chromosomes. The common co-occurrence of all three features makes it difficult to pinpoint their respective roles. However, they do not always co-occur in nature, and missing cases can be additionally studied with hypothetical life cycles. We model sexually antagonistic alleles in eight different sex determination systems and find that arguments 1 and 2 promote invasion and fixation of female-beneficial and male-beneficial alleles, respectively; argument 2 also improves prospects for polymorphism. Argument 3 harms the invasion prospects of sexually antagonistic alleles (irrespective of which sex benefits) but promotes fixation should invasion nevertheless occur. Disentangling the features helps to evaluate the validity of previous verbal arguments and yields better-informed predictions about intralocus sexual conflict under different sex determination systems, including hitherto undiscovered ones.

MicroRNAs May Play an Important Role in Sexual Reversal Process of Chinese Soft-Shelled Turtle, Pelodiscus sinensis

The Chinese soft-shelled (Pelodiscus sinensis) turtle exhibits obvious sex dimorphism, which leads to the higher economic and nutritional value of male individuals. Exogenous hormones can cause the transformation from male to female phenotype during gonadal differentiation. However, the molecular mechanism related to the sexual reversal process is unclear. In this study, we compared the difference between the small RNAs of male, female, and pseudo-female turtles by small RNA-seq to understand the sexual reversal process of Chinese soft-shelled turtles. A certain dose of estrogen can cause the transformation of Chinese soft-shelled turtles from male to female, which are called pseudo-female individuals. The result of small RNA-seq has revealed that the characteristics of pseudo-females are very similar to females, but are strikingly different from males. The number of the microRNAs (miRNAs) of male individuals was significantly less than the number of female individuals or pseudo-female individuals, while the expression level of miRNAs of male individuals were significantly higher than the other two types. Furthermore, we found 533 differentially expressed miRNAs, including 173 up-regulated miRNAs and 360 down-regulated miRNAs, in the process of transformation from male to female phenotype. Cluster analysis of the total 602 differential miRNAs among females, males, and pseudo-females showed that miRNAs played a crucial role during the sexual differentiation. Among these differential miRNAs, we found 12 miRNAs related to gonadal development and verified their expression by qPCR. The TR-qPCR results confirmed the differential expression of 6 of the 12 miRNAs: miR-26a-5p, miR-212-5p, miR-202-5p, miR-301a, miR-181b-3p and miR-96-5p were involved in sexual reversal process, which was consistent with the results of omics. Using these six miRNAs and some of their target genes, we constructed a network diagram related to gonadal development. We suggest that these miRNAs may play an important role in the process of effective sex reversal, which would contribute to the breeding of all male strains of Chinese soft-shelled turtles.

Rapid evolution of sexual size dimorphism facilitated by Y-linked genetic variance

Sexual dimorphism is ubiquitous in nature but its evolution is puzzling given that the mostly shared genome constrains independent evolution in the sexes. Sex differences should result from asymmetries between the sexes in selection or genetic variation but studies investigating both simultaneously are lacking. Here, we combine a quantitative genetic analysis of body size variation, partitioned into autosomal and sex chromosome contributions and ten generations of experimental evolution to dissect the evolution of sexual body size dimorphism in seed beetles (Callosobruchus maculatus) subjected to sexually antagonistic or sex-limited selection. Female additive genetic variance (V_A) was primarily linked to autosomes, exhibiting a strong intersexual genetic correlation with males ([Formula: see text] = 0.926), while X- and Y-linked genes further contributed to the male V_A and X-linked genes contributed to female dominance variance. Consistent with these estimates, sexual body size dimorphism did not evolve in response to female-limited selection but evolved by 30-50% under male-limited and sexually antagonistic selection. Remarkably, Y-linked variance alone could change dimorphism by 30%, despite the C. maculatus Y chromosome being small and heterochromatic. Our results demonstrate how the potential for sexual dimorphism to evolve depends on both its underlying genetic basis and the nature of sex-specific selection.

Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish

Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.

Differential gene expression patterns during gametophyte development provide insights into sex differentiation in the dioicous kelp Saccharina japonica

BACKGROUND

In brown algae, dioicy is the prevalent sexual system, and phenotypic differences between male and female gametophytes have been found in many dioicous species. Saccharina japonica show remarkable sexual dimorphism in gametophytes before gametogenesis. A higher level of phenotypic differentiation was also found in female and male gametes after gametogenesis. However, the patterns of differential gene expression throughout gametophyte development and how these changes might relate to sex-specific fitness at the gamete stage in S. japonica are not well known.

RESULTS

In this study, differences in gene expression between male and female gametophytes in different developmental stages were investigated using comparative transcriptome analysis. Among the 20,151 genes expressed in the haploid gametophyte generation, 37.53% were sex-biased. The abundance of sex-biased genes in mature gametophytes was much higher than that in immature gametophytes, and more male-biased than female-biased genes were observed in the mature stage. The predicted functions of most sex-biased genes were closely related to the sex-specific characteristics of gametes, including cell wall biosynthesis, sperm motility, and sperm and egg recognition. In addition, 51 genes were specifically expressed in males in both stages, showing great potential as candidate male sex-determining region (SDR) genes.

CONCLUSIONS

This study describes a thorough investigation into differential gene expression between male and female gametophytes in the dioicous kelp S. japonica. A large number of sex-biased genes in mature gametophytes may be associated with the divergence of phenotypic traits and physiological functions between female gametes (eggs) and male gametes (sperm) during sexual differentiation. These genes may mainly come from new sex-biased genes that have recently evolved in the S. japonica lineage. The duplication of sex-biased genes was detected, which may increase the number of sex-biased genes after gametogenesis in S. japonica to some extent. The excess of male-biased genes over female-biased genes in the mature stage may reflect the different levels of sexual selection across sexes. This study deepens our understanding of the regulation of sex development and differentiation in the dioicous kelp S. japonica.

A role for sex-determination genes in life history evolution? Doublesex mediates sexual size dimorphism in the gazelle dung beetle

An organism's fitness depends strongly on its age and size at maturation. Although the evolutionary forces acting on these critical life history traits have been heavily scrutinized, the developmental mechanisms underpinning intraspecific variation in adult size and development time remain much less well-understood. Using RNA interference, I here show that the highly conserved sex-determination gene doublesex (dsx) mediates sexual size dimorphism (SSD) in the gazelle dung beetle Digitonthophagus gazella. Because doublesex undergoes sex-specific splicing and sex-limited isoforms regulate different target genes, this suggests that dsx contributes to the resolution of intralocus sexual conflict in body size. However, these results contrast with previous studies demonstrating that dsx does not affect body size or SSD in Drosophila. This indicates that intraspecific body size variation is underlain by contrasting developmental mechanisms in different insect lineages. Furthermore, although male D. gazella have a longer development time than females, sexual bimaturism was not affected by dsx expression knockdown. In addition, and in contrast to secondary sexual morphology, dsx did not significantly affect nutritional plasticity in life history. Taken together, these findings indicate that dsx signalling contributes to intraspecific life history variation but that dsx's function in mediating sexual dimorphism in life history differs among traits and species. More generally, these findings suggest that genes ancestrally tasked with sex determination have been co-opted into the developmental regulation of life history traits and may represent an underappreciated mechanism of life history evolution.

Impact of male trait exaggeration on sex-biased gene expression and genome architecture in a water strider

BACKGROUND

Exaggerated secondary sexual traits are widespread in nature and often evolve under strong directional sexual selection. Although heavily studied from both theoretical and empirical viewpoints, we have little understanding of how sexual selection influences sex-biased gene regulation during the development of exaggerated secondary sexual phenotypes, and how these changes are reflected in genomic architecture. This is primarily due to the limited availability of representative genomes and associated tissue and sex transcriptomes to study the development of these traits. Here we present the genome and developmental transcriptomes, focused on the legs, of the water strider Microvelia longipes, a species where males exhibit strikingly long third legs compared to females, which they use as weapons.

RESULTS

We generated a high-quality genome assembly with 90% of the sequence captured in 13 scaffolds. The most exaggerated legs in males were particularly enriched in both sex-biased and leg-biased genes, indicating a specific signature of gene expression in association with trait exaggeration. We also found that male-biased genes showed patterns of fast evolution compared to non-biased and female-biased genes, indicative of directional or relaxed purifying selection. By contrast to male-biased genes, female-biased genes that are expressed in the third legs, but not the other legs, are over-represented in the X chromosome compared to the autosomes. An enrichment analysis for sex-biased genes along the chromosomes revealed also that they arrange in large genomic regions or in small clusters of two to four consecutive genes. The number and expression of these enriched regions were often associated with the exaggerated legs of males, suggesting a pattern of common regulation through genomic proximity in association with trait exaggeration.

CONCLUSION

Our findings indicate how directional sexual selection may drive sex-biased gene expression and genome architecture along the path to trait exaggeration and sexual dimorphism.

A First Insight into the Gonad Transcriptome of Hong Kong Catfish (Clarias fuscus)

Hong Kong catfish (Clarias fuscus) exhibit sexual dimorphism, particularly in body size. Due to the fast growth rate of males, the sexual size dimorphism of Hong Kong catfish has become an economically important trait. However, limited knowledge is known about the molecular mechanisms of sex determination and sex differentiation in this species. In this study, a first de novo transcriptome sequencing analysis of testes and ovaries was performed to identify sex-biased genes in Hong Kong catfish. The results showed that a total of 290,291 circular consensus sequences (CCSs) were obtained, from which 248,408 full-length non-chimeric (FLNC) reads were generated. After non-redundant analysis, a total of 37,305 unigenes were predicted, in which 34,342 unigenes were annotated with multiple public databases. Comparative transcriptomic analysis identified 5750 testis-biased differentially expressed genes (DEGs) and 6991 ovary-biased DEGs. The enrichment analysis showed that DEGs were classified into 783 Gene Ontology (GO) terms and 16 Kyoto Encyclopedia of Gene and Genome (KEGG) pathways. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, androgen secretion, gonadal development and steroid biosynthesis pathways. In addition, the expression levels of 23 unigenes were confirmed to validate the transcriptomic data by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first investigation into the transcriptome of Hong Kong catfish testes and ovaries. This study provides an important molecular basis for the sex determination and sex control breeding of Hong Kong catfish.

The Aphid X Chromosome Is a Dangerous Place for Functionally Important Genes: Diverse Evolution of Hemipteran Genomes Based on Chromosome-Level Assemblies

Different evolutionary forces shape gene content and sequence evolution on autosomes versus sex chromosomes. Location on a sex chromosome can favor male-beneficial or female-beneficial mutations depending on the sex determination system and selective pressure on different sexual morphs. An X0 sex determination can lead to autosomal enrichment of male-biased genes, as observed in some hemipteran insect species. Aphids share X0 sex determination; however, models predict the opposite pattern, due to their unusual life cycles, which alternate between all-female asexual generations and a single sexual generation. Predictions include enrichment of female-biased genes on autosomes and of male-biased genes on the X, in contrast to expectations for obligately sexual species. Robust tests of these models require chromosome-level genome assemblies for aphids and related hemipterans with X0 sex determination and obligate sexual reproduction. In this study, we built the first chromosome-level assembly of a psyllid, an aphid relative with X0 sex determination and obligate sexuality, and compared it with recently resolved chromosome-level assemblies of aphid genomes. Aphid and psyllid X chromosomes differ strikingly. In aphids, female-biased genes are strongly enriched on autosomes and male-biased genes are enriched on the X. In psyllids, male-biased genes are enriched on autosomes. Furthermore, functionally important gene categories of aphids are enriched on autosomes. Aphid X-linked genes and male-biased genes are under relaxed purifying selection, but gene content and order on the X is highly conserved, possibly reflecting constraints imposed by unique chromosomal mechanisms associated with the unusual aphid life cycle.

Sex-linked genetic diversity and differentiation in a globally distributed avian species complex

Sex chromosomes often bear distinct patterns of genetic variation due to unique patterns of inheritance and demography. The processes of mutation, recombination, genetic drift and selection also influence rates of evolution on sex chromosomes differently than autosomes. Measuring such differences provides information about how these processes shape genomic variation and their roles in the origin of species. To test hypotheses and predictions about patterns of autosomal and sex-linked genomic diversity and differentiation, we measured population genetic statistics within and between populations and subspecies of the barn swallow (Hirundo rustica) and performed explicit comparisons between autosomal and Z-linked genomic regions. We first tested for evidence of low Z-linked genetic diversity and high Z-linked population differentiation relative to autosomes, then for evidence that the Z chromosome bears greater ancestry information due to faster lineage sorting. Finally, we investigated geographical clines across hybrid zones for evidence that the Z chromosome is resistant to introgression due to selection against hybrids. We found evidence that the barn swallow mating system, demographic history and linked selection each contribute to low Z-linked diversity and high Z-linked differentiation. While incomplete lineage sorting is rampant across the genome, our results indicate faster sorting of ancestral polymorphism on the Z. Finally, hybrid zone analyses indicate barriers to introgression on the Z chromosome, suggesting that sex-linked traits are important in reproductive isolation, especially in migratory divide regions. Our study highlights how selection, gene flow and demography shape sex-linked genetic diversity and underlines the relevance of the Z chromosome in speciation.

X chromosome escapee genes are involved in ischemic sexual dimorphism through epigenetic modification of inflammatory signals

BACKGROUND

Stroke is a sexually dimorphic disease. Previous studies have found that young females are protected against ischemia compared to males, partially due to the protective effect of ovarian hormones, particularly estrogen (E2). However, there are also genetic and epigenetic effects of X chromosome dosage that contribute to stroke sensitivity and neuroinflammation after injury, especially in the aged. Genes that escape from X chromosome inactivation (XCI) contribute to sex-specific phenotypes in many disorders. Kdm5c and kdm6a are X escapee genes that demethylate H3K4me3 and H3K27me3, respectively. We hypothesized that the two demethylases play critical roles in mediating the stroke sensitivity.

METHODS

To identify the X escapee genes involved in stroke, we performed RNA-seq in flow-sorted microglia from aged male and female wild type (WT) mice subjected to middle cerebral artery occlusion (MCAO). The expression of these genes (kdm5c/kdm6a) were confirmed in four core genotypes (FCG) mice and in post-mortem human stroke brains by immunohistochemistry (IHC), Western blot, and RT-PCR. Chromatin immunoprecipitation (ChIP) assays were conducted to detect DNA levels of inflammatory interferon regulatory factor (IRF) 4/5 precipitated by histone H3K4 and H3K27 antibodies. Manipulation of kdm5c/kdm6a expression with siRNA or lentivirus was performed in microglial culture, to determine downstream pathways and examine the regulatory roles in inflammatory cytokine production.

RESULTS

Kdm5c and kdm6a mRNA levels were significantly higher in aged WT female vs. male microglia, and the sex difference also existed in ischemic brains from FCG mice and human stroke patients. The ChIP assay showed the IRF 4/5 had higher binding levels to demethylated H3K4 or H3K27, respectively, in female vs. male ischemic microglia. Knockdown or over expression of kdm5c/kdm6a with siRNA or lentivirus altered the methylation of H3K4 or H3K27 at the IRF4/5 genes, which in turn, impacted the production of inflammatory cytokines.

CONCLUSIONS

The KDM-Histone-IRF pathways are suggested to mediate sex differences in cerebral ischemia. Epigenetic modification of stroke-related genes constitutes an important mechanism underlying the ischemic sexual dimorphism.

Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.

Mapping genomic regions for reproductive traits in beef cattle: Inclusion of the X chromosome

Although the second largest chromosome of the genome, the X chromosome is usually excluded from genome-wide association studies (GWAS). Considering the presence and importance of genes on this chromosome that are involved in reproduction, the aim of this study was to evaluate the effect of its inclusion in GWAS on reproductive traits (scrotal circumference [SC], early pregnancy [P16] and age at first calving [AFC]) in a Nelore herd. Genotype data from 3,263 animals with the above-mentioned phenotypes were used. The results showed an increase in the variances explained by the autosomal markers for all traits when the X chromosome was not included. For SC, there was an increase of more than 10% for the windows on chromosomes 2 and 6. For P16, the effect was increased by almost 20% for windows on chromosome 5. The same pattern was found for AFC, with an increase of more than 10% for the most important windows. The results indicate that the noninclusion of the X chromosome can overestimate the effects of autosomes on SC, P16 and AFC not only because of the additive effect of the X chromosome itself but also because of its epistatic effect on autosomal genes.

Partial Amniote Sex Chromosomal Linkage Homologies Shared on Snake W Sex Chromosomes Support the Ancestral Super-Sex Chromosome Evolution in Amniotes

Squamate reptile chromosome 2 (SR2) is thought to be an important remnant of an ancestral amniote super-sex chromosome, but a recent study showed that the Siamese cobra W sex chromosome is also a part of this larger ancestral chromosome. To confirm the existence of an ancestral amniote super-sex chromosome and understand the mechanisms of amniote sex chromosome evolution, chromosome maps of two snake species [Russell's viper: Daboia russelii (DRU) and the common tiger snake: Notechis scutatus (NSC)] were constructed using bacterial artificial chromosomes (BACs) derived from chicken and zebra finch libraries containing amniote sex chromosomal linkages. Sixteen BACs were mapped on the W sex chromosome of DRU and/or NSC, suggesting that these BACs contained a common genomic region shared with the W sex chromosome of these snakes. Two of the sixteen BACs were co-localized to DRU2 and NSC2, corresponding to SR2. Prediction of genomic content from all BACs mapped on snake W sex chromosomes revealed a large proportion of long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) retrotransposons. These results led us to predict that amplification of LINE and SINE may have occurred on snake W chromosomes during evolution. Genome compartmentalization, such as transposon amplification, might be the key factor influencing chromosome structure and differentiation. Multiple sequence alignments of all BACs mapped on snake W sex chromosomes did not reveal common sequences. Our findings indicate that the SR2 and snake W sex chromosomes may have been part of a larger ancestral amniote super-sex chromosome, and support the view of sex chromosome evolution as a colorful myriad of situations and trajectories in which many diverse processes are in action.

A gene's-eye view of sexual antagonism

Females and males may face different selection pressures. Accordingly, alleles that confer a benefit for one sex often incur a cost for the other. Classic evolutionary theory holds that the X chromosome, whose sex-biased transmission sees it spending more time in females, should value females more than males, whereas autosomes, whose transmission is unbiased, should value both sexes equally. However, recent mathematical and empirical studies indicate that male-beneficial alleles may be more favoured by the X chromosome than by autosomes. Here we develop a gene's-eye-view approach that reconciles the classic view with these recent discordant results, by separating a gene's valuation of female versus male fitness from its ability to induce fitness effects in either sex. We use this framework to generate new comparative predictions for sexually antagonistic evolution in relation to dosage compensation, sex-specific mortality and assortative mating, revealing how molecular mechanisms, ecology and demography drive variation in masculinization versus feminization across the genome.

Differences in the contributions of sex linkage and androgen regulation to sex-biased gene expression in juvenile and adult sticklebacks

Different evolutionary interests between males and females can lead to the evolution of sexual dimorphism. However, intersex genetic correlations due to the shared genome can constrain the evolution of sexual dimorphism, resulting in intra-locus sexual conflict. One of the mechanisms resolving this conflict is sex linkage, which allows males and females to carry different alleles on sex chromosomes. Another is a regulatory mutation causing sex-biased gene expression, which is often mediated by gonadal steroids in vertebrates. How do these two mechanisms differ in the contributions to the resolution of intra-locus sexual conflict? The magnitude of sexual conflict often varies between the juvenile and adult stages. Because gonadal steroids change in titre during development, we hypothesized that gonadal steroids play a role in sexual dimorphism expression only at certain developmental stages, whereas sex linkage is more important for sexual dimorphism expressed throughout life. Our brain transcriptome analysis of juvenile and adult threespine sticklebacks showed that the majority of genes that were sex-biased in both stages were sex-linked. The relative contribution of androgen-dependent regulation to the sex-biased transcriptome increased and that of sex linkage declined in adults compared to juveniles. The magnitude of the sex differences was greater in sex-linked genes than androgen-responsive genes, suggesting that sex linkage is more effective than androgen regulation in the production of large sex differences in gene expression. Overall, our data are consistent with the hypothesis that sex linkage is effective in resolving sexual conflict throughout life, whereas androgen-dependent regulation can contribute to temporary resolution of sexual conflict.

The microevolutionary response to male-limited X-chromosome evolution in Drosophila melanogaster reflects macroevolutionary patterns

Due to its hemizygous inheritance and role in sex determination, the X-chromosome is expected to play an important role in the evolution of sexual dimorphism and to be enriched for sexually antagonistic genetic variation. By forcing the X-chromosome to only be expressed in males over >40 generations, we changed the selection pressures on the X to become similar to those experienced by the Y. This releases the X from any constraints arising from selection in females and should lead to specialization for male fitness, which could occur either via direct effects of X-linked loci or trans-regulation of autosomal loci by the X. We found evidence of masculinization via up-regulation of male-benefit sexually antagonistic genes and down-regulation of X-linked female-benefit genes. Potential artefacts of the experimental evolution protocol are discussed and cannot be wholly discounted, leading to several caveats. Interestingly, we could detect evidence of microevolutionary changes consistent with previously documented macroevolutionary patterns, such as changes in expression consistent with previously established patterns of sexual dimorphism, an increase in the expression of metabolic genes related to mito-nuclear conflict and evidence that dosage compensation effects can be rapidly altered. These results confirm the importance of the X in the evolution of sexual dimorphism and as a source for sexually antagonistic genetic variation and demonstrate that experimental evolution can be a fruitful method for testing theories of sex chromosome evolution.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

The genomic footprint of sexual conflict

Genes with sex-biased expression show a number of unique properties and this has been seen as evidence for conflicting selection pressures in males and females, forming a genetic ‘tug-of-war’ between the sexes. However, we lack studies of taxa where an understanding of conflicting phenotypic selection in the sexes has been linked with studies of genomic signatures of sexual conflict. Here, we provide such a link. We used an insect where sexual conflict is unusually well understood, the seed beetle Callosobruchus maculatus, to test for molecular genetic signals of sexual conflict across genes with varying degrees of sex-bias in expression. We sequenced, assembled and annotated its genome and performed population resequencing of three divergent populations. Sex-biased genes showed increased levels of genetic diversity and bore a remarkably clear footprint of relaxed purifying selection. Yet, segregating genetic variation was also affected by balancing selection in weakly female-biased genes, while male-biased genes showed signs of overall purifying selection. Female-biased genes contributed disproportionally to shared polymorphism across populations, while male-biased genes, male seminal fluid protein genes and sex-linked genes did not. Genes showing genomic signatures consistent with sexual conflict generally matched life-history phenotypes known to experience sexually antagonistic selection in this species. Our results highlight metabolic and reproductive processes, confirming the key role of general life-history traits in sexual conflict.

Variation in sexual dimorphism in a wind-pollinated plant: the influence of geographical context and life-cycle dynamics

Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less information is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life-cycle dynamics. Here, we investigated patterns of genetically based sexual dimorphism in vegetative and reproductive traits of a wind-pollinated dioecious plant, Rumex hastatulus, across three life-cycle stages using open-pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1 Y2 ) of the species. The direction and degree of sexual dimorphism was highly variable among populations and life-cycle stages. Sex-specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. Sex-specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life-cycle.

Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb

BACKGROUND AND AIMS

Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.

METHODS

The present study documents differences in gene expression in both above- and below-ground tissues of early pre-reproductive individuals of the wind-pollinated dioecious annual herb, Mercurialis annua, which otherwise shows clear sexual dimorphism only at the adult stage.

KEY RESULTS

Whereas males and females differed in their gene expression at the first leaf stage, sex-biased gene expression peaked just prior to, and after, flowering, as might be expected if sexual dimorphism is partly a response to differential costs of reproduction. Sex-biased genes were over-represented among putative sex-linked genes in M. annua but showed no evidence for more rapid evolution than unbiased genes.

CONCLUSIONS

Sex-biased gene expression in M. annua occurs as early as the first whorl of leaves is produced, is highly dynamic during plant development and varies substantially between vegetative tissues.

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.

Masculinization of gene expression is associated with male quality in Drosophila melanogaster

The signature of sexual selection has been revealed through the study of differences in patterns of genome-wide gene expression, both between the sexes and between alternative reproductive morphs within a single sex. What remains unclear, however, is whether differences in gene expression patterns between individuals of a given sex consistently map to variation in individual quality. Such a pattern, particularly if found in males, would provide unambiguous evidence that the phenotypic response to sexual selection is shaped through sex-specific alterations to the transcriptome. To redress this knowledge gap, we explored whether patterns of sex-biased gene expression are associated with variation in male reproductive quality in Drosophila melanogaster. We measured two male reproductive phenotypes, and their association with sex-biased gene expression, across a selection of inbred lines from the Drosophila Genetic Reference Panel. Genotypes with higher expression of male-biased genes produced males exhibiting shorter latencies to copulation, and higher capacity to inseminate females. Conversely, female-biased genes tended to show negative associations with these male reproductive traits across genotypes. We uncovered similar patterns, by reanalyzing a published dataset from a second D. melanogaster population. Our results reveal the footprint of sexual selection in masculinising the male transcriptome.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long-term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant-specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model-based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population-specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host-plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host-plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.

High Levels of Copy Number Variation of Ampliconic Genes across Major Human Y Haplogroups

Because of its highly repetitive nature, the human male-specific Y chromosome remains understudied. It is important to investigate variation on the Y chromosome to understand its evolution and contribution to phenotypic variation, including infertility. Approximately 20% of the human Y chromosome consists of ampliconic regions which include nine multi-copy gene families. These gene families are expressed exclusively in testes and usually implicated in spermatogenesis. Here, to gain a better understanding of the role of the Y chromosome in human evolution and in determining sexually dimorphic traits, we studied ampliconic gene copy number variation in 100 males representing ten major Y haplogroups world-wide. Copy number was estimated with droplet digital PCR. In contrast to low nucleotide diversity observed on the Y in previous studies, here we show that ampliconic gene copy number diversity is very high. A total of 98 copy-number-based haplotypes were observed among 100 individuals, and haplotypes were sometimes shared by males from very different haplogroups, suggesting homoplasies. The resulting haplotypes did not cluster according to major Y haplogroups. Overall, only two gene families (RBMY and TSPY) showed significant differences in copy number among major Y haplogroups, and the haplogroup of a male could not be predicted based on his ampliconic gene copy numbers. Finally, we did not find significant correlations either between copy number variation and individual's height, or between the former and facial masculinity/femininity. Our results suggest rapid evolution of ampliconic gene copy numbers on the human Y, and we discuss its causes.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

Evolution of sex-specific pace-of-life syndromes: genetic architecture and physiological mechanisms

Sex differences in life history, physiology, and behavior are nearly ubiquitous across taxa, owing to sex-specific selection that arises from different reproductive strategies of the sexes. The pace-of-life syndrome (POLS) hypothesis predicts that most variation in such traits among individuals, populations, and species falls along a slow-fast pace-of-life continuum. As a result of their different reproductive roles and environment, the sexes also commonly differ in pace-of-life, with important consequences for the evolution of POLS. Here, we outline mechanisms for how males and females can evolve differences in POLS traits and in how such traits can covary differently despite constraints resulting from a shared genome. We review the current knowledge of the genetic basis of POLS traits and suggest candidate genes and pathways for future studies. Pleiotropic effects may govern many of the genetic correlations, but little is still known about the mechanisms involved in trade-offs between current and future reproduction and their integration with behavioral variation. We highlight the importance of metabolic and hormonal pathways in mediating sex differences in POLS traits; however, there is still a shortage of studies that test for sex specificity in molecular effects and their evolutionary causes. Considering whether and how sexual dimorphism evolves in POLS traits provides a more holistic framework to understand how behavioral variation is integrated with life histories and physiology, and we call for studies that focus on examining the sex-specific genetic architecture of this integration.

Autosomal and X-Linked Additive Genetic Variation for Lifespan and Aging: Comparisons Within and Between the Sexes in Drosophila melanogaster

Theory makes several predictions concerning differences in genetic variation between the X chromosome and the autosomes due to male X hemizygosity. The X chromosome should: (i) typically show relatively less standing genetic variation than the autosomes, (ii) exhibit more variation in males compared to females because of dosage compensation, and (iii) potentially be enriched with sex-specific genetic variation. Here, we address each of these predictions for lifespan and aging in Drosophila melanogaster. To achieve unbiased estimates of X and autosomal additive genetic variance, we use 80 chromosome substitution lines; 40 for the X chromosome and 40 combining the two major autosomes, which we assay for sex-specific and cross-sex genetic (co)variation. We find significant X and autosomal additive genetic variance for both traits in both sexes (with reservation for X-linked variation of aging in females), but no conclusive evidence for depletion of X-linked variation (measured through females). Males display more X-linked variation for lifespan than females, but it is unclear if this is due to dosage compensation since also autosomal variation is larger in males. Finally, our results suggest that the X chromosome is enriched for sex-specific genetic variation in lifespan but results were less conclusive for aging overall. Collectively, these results suggest that the X chromosome has reduced capacity to respond to sexually concordant selection on lifespan from standing genetic variation, while its ability to respond to sexually antagonistic selection may be augmented.