Sexual dimorphism results from sex-biased gene expression, which evolves when selection acts differently on males and females. While there is an intimate connection between sex-biased gene expression and sex-specific selection, few empirical studies have studied this relationship directly. Here we compare the two on a genome-wide scale in humans and flies. We find a distinctive “Twin Peaks” pattern in humans that relates the strength of sex-specific selection, quantified by genetic divergence between male and female adults at autosomal loci, to the degree of sex-biased expression. Genes with intermediate degrees of sex-biased expression show evidence of ongoing sex-specific selection, while genes with either little or completely sex-biased expression do not. This pattern apparently results from differential viability selection in males and females acting in the current generation. The Twin Peaks pattern is also found in Drosophila using a different measure of sex-specific selection acting on fertility. We develop a simple model that successfully recapitulates the Twin Peaks. Our results suggest that many genes with intermediate sex-biased expression experience ongoing sex-specific selection in humans and flies.

Previous efforts to uncover the genetic underpinnings of ongoing ecological speciation of the M and S forms of the African malaria vector Anopheles gambiae revealed two centromere-proximal islands of genetic divergence on X and chromosome 2. Under the assumption of considerable ongoing gene flow between M and S, these persistently divergent genomic islands were widely considered to be “speciation islands”. In the course of microarray-based divergence mapping, we discovered a third centromere-associated island of divergence on chromosome 3, which was validated by targeted re-sequencing. To test for genetic association between the divergence islands on all three chromosomes, SNP-based assays were applied in four natural populations of M and S spanning West, Central and East Africa. Genotyping of 517 female M and S mosquitoes revealed nearly complete linkage disequilibrium between the centromeres of the three independently assorting chromosomes. These results suggest that despite the potential for inter-form gene flow through hybridization, actual (realized) gene flow between M and S may be substantially less than commonly assumed, and may not explain most shared variation. Moreover, the possibility of very low gene flow calls into question whether diverged pericentromeric regions-- characterized by reduced levels of variation and recombination-- are in fact instrumental rather than merely incidental to the speciation process.

The association between fitness-related phenotypic traits and an environmental gradient offers one of the best opportunities to study the interplay between natural selection and migration. In cases in which specific genetic variants also show such clinal patterns, it may be possible to uncover the mutations responsible for local adaptation. The malaria vector, Anopheles gambiae, is associated with a latitudinal cline in aridity in Cameroon; a large inversion on chromosome 2L of this mosquito shows large differences in frequency along this cline, with high frequencies of the inverted karyotype present in northern, more arid populations and an almost complete absence of the inverted arrangement in southern populations. Here we use a genome resequencing approach to investigate patterns of population divergence along the cline. By sequencing pools of individuals from both ends of the cline as well as in the center of the cline-where the inversion is present in intermediate frequency-we demonstrate almost complete panmixia across collinear parts of the genome and high levels of differentiation in inverted parts of the genome. Sequencing of separate pools of each inversion arrangement in the center of the cline reveals large amounts of gene flux (i.e., gene conversion and double crossovers) even within inverted regions, especially away from the inversion breakpoints. The interplay between natural selection, migration, and gene flux allows us to identify several candidate genes responsible for the match between inversion frequency and environmental variables. These results, coupled with similar conclusions from studies of clinal variation in Drosophila, point to a number of important biological functions associated with local environmental adaptation.

Y chromosomes control essential male functions in many species, including sex determination and fertility. However, because of obstacles posed by repeat-rich heterochromatin, knowledge of Y chromosome sequences is limited to a handful of model organisms, constraining our understanding of Y biology across the tree of life. Here, we leverage long single-molecule sequencing to determine the content and structure of the nonrecombining Y chromosome of the primary African malaria mosquito, Anopheles gambiae. We find that the An. gambiae Y consists almost entirely of a few massively amplified, tandemly arrayed repeats, some of which can recombine with similar repeats on the X chromosome. Sex-specific genome resequencing in a recent species radiation, the An. gambiae complex, revealed rapid sequence turnover within An. gambiae and among species. Exploiting 52 sex-specific An. gambiae RNA-Seq datasets representing all developmental stages, we identified a small repertoire of Y-linked genes that lack X gametologs and are not Y-linked in any other species except An. gambiae, with the notable exception of YG2, a candidate male-determining gene. YG2 is the only gene conserved and exclusive to the Y in all species examined, yet sequence similarity to YG2 is not detectable in the genome of a more distant mosquito relative, suggesting rapid evolution of Y chromosome genes in this highly dynamic genus of malaria vectors. The extensive characterization of the An. gambiae Y provides a long-awaited foundation for studying male mosquito biology, and will inform novel mosquito control strategies based on the manipulation of Y chromosomes.

The African malaria mosquito Anopheles gambiae is diversifying into ecotypes known as M and S forms. This process is thought to be promoted by adaptation to different larval habitats, but its genetic underpinnings remain elusive. To identify candidate targets of divergent natural selection in M and S, we performed genomewide scanning in paired population samples from Mali, followed by resequencing and genotyping from five locations in West, Central, and East Africa. Genome scans revealed a significant peak of M-S divergence on chromosome 3L, overlapping five known or suspected immune response genes. Resequencing implicated a selective target at or near the TEP1 gene, whose complement C3-like product has antiparasitic and antibacterial activity. Sequencing and allele-specific genotyping showed that an allelic variant of TEP1 has been swept to fixation in M samples from Mali and Burkina Faso and is spreading into neighboring Ghana, but is absent from M sampled in Cameroon, and from all sampled S populations. Sequence comparison demonstrates that this allele is related to, but distinct from, TEP1 alleles of known resistance phenotype. Experimental parasite infections of advanced mosquito intercrosses demonstrated a strong association between this TEP1 variant and resistance to both rodent malaria and the native human malaria parasite Plasmodium falciparum. Although malaria parasites may not be direct agents of pathogen-mediated selection at TEP1 in nature-where larvae may be the more vulnerable life stage-the process of adaptive divergence between M and S has potential consequences for malaria transmission.ecological speciation | malaria vector | population genomics | thioester immune gene