Multilocus patterns of genetic variation across Cryptosporidium species suggest balancing selection at the gp60 locus

Cryptosporidium is an apicomplexan protozoan that lives in most vertebrates, including humans. Its gp60 gene is functionally involved in its attachment to host cells, and its high level of genetic variation has made it the reference marker for sample typing in epidemiological studies. To understand the origin of such high diversity and to determine the extent to which this classification applies to the rest of the genome, we analysed the patterns of variation at gp60 and nine other nuclear loci in isolates of three Cryptosporidium species. Most loci showed low genetic polymorphism (π_S <1%) and similar levels of between‐species divergence. Contrastingly, gp60 exhibited very different characteristics: (i) it was nearly ten times more variable than the other loci; (ii) it displayed a significant excess of polymorphisms relative to between‐species differences in a maximum‐likelihood Hudson–Kreitman–Aguadé test; (iii) gp60 subtypes turned out to be much older than the species they were found in; and (iv) showed a significant excess of polymorphic variants shared across species from random expectations. These observations suggest that this locus evolves under balancing selection and specifically under negative frequency‐dependent selection (FDS). Interestingly, genetic variation at the other loci clusters very well within the groups of isolates defined by gp60 subtypes, which may provide new tools to understand the genome‐wide patterns of genetic variation of the parasite in the wild. These results suggest that gp60 plays an active and essential role in the life cycle of the parasite and that genetic variation at this locus might be essential for the parasite's long‐term success.     

Selection outweighs drift at a fine scale: Lack of MHC differentiation within a family living lizard across geographically close but disconnected rocky outcrops

The highly polymorphic genes of the major histocompatibility complex (MHC) are involved in disease resistance, mate choice and kin recognition. Therefore, they are widely used markers for investigating adaptive variation. Although selection is the key driver, gene flow and genetic drift also influence adaptive genetic variation, sometimes in opposing ways and with consequences for adaptive potential. To further understand the processes that generate MHC variation, it is helpful to compare variation at the MHC with that at neutral genetic loci. Differences in MHC and neutral genetic variation are useful for inferring the relative influence of selection, gene flow and drift on MHC variation. To date, such investigations have usually been undertaken at a broad spatial scale. Yet, evolutionary and ecological processes can occur at a fine spatial scale, particularly in small or fragmented populations. We investigated spatial patterns of MHC variation among three geographically close, naturally discrete, sampling sites of Egernia stokesii, an Australian lizard. The MHC of E. stokesii has recently been characterized, and there is evidence for historical selection on the MHC. We found E. stokesii MHC weakly differentiated among sites compared to microsatellites, suggesting selection, acting similarly at each site, has outweighed any effects of low gene flow or of genetic drift on E. stokesii MHC variation. Our findings demonstrate the strength of selection in shaping patterns of MHC variation or consistency at a fine spatial scale.     

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.     

Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Balancing selection and genetic drift create unusual patterns of MHCIIβ variation in Galápagos mockingbirds

The extracellular subunit of the major histocompatibility complex MHCIIβ plays an important role in the recognition of pathogens and the initiation of the adaptive immune response of vertebrates. It is widely accepted that pathogen‐mediated selection in combination with neutral micro‐evolutionary forces (e.g. genetic drift) shape the diversity of MHCIIβ, but it has proved difficult to determine the relative effects of these forces. We evaluated the effect of genetic drift and balancing selection on MHCIIβ diversity in 12 small populations of Galápagos mockingbirds belonging to four different species, and one larger population of the Northern mockingbird from the continental USA. After genotyping MHCIIβ loci by high‐throughput sequencing, we applied a correlational approach to explore the relationships between MHCIIβ diversity and population size by proxy of island size. As expected when drift predominates, we found a positive effect of population size on the number of MHCIIβ alleles present in a population. However, the number of MHCIIβ alleles per individual and number of supertypes were not correlated with population size. This discrepancy points to an interesting feature of MHCIIβ diversity dynamics: some levels of diversity might be shaped by genetic drift while others are independent and possibly maintained by balancing selection.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.     

Contrasting patterns of polymorphism and selection in bacterial‐sensing toll‐like receptor 4 in two house mouse subspecies

Detailed investigation of variation in genes involved in pathogen recognition is crucial for understanding co‐evolutionary processes between parasites and their hosts. Triggering immediate innate response to invading microbes, Toll‐like receptors (TLRs) belong presently among the best‐studied receptors of vertebrate immunity. TLRs exhibit remarkable interspecific variation and also intraspecific polymorphism is well documented. In humans and laboratory mice, several studies have recently shown that single amino acid substitution may significantly alter receptor function. Unfortunately, data concerning polymorphism in free‐living species are still surprisingly scarce. In this study, we analyzed the polymorphism of Toll‐like receptor 4 (Tlr4) over the Palearctic range of house mouse (Mus musculus). Our results reveal contrasting evolutionary patterns between the two recently (0.5 million years ago) diverged house mouse subspecies: M. m. domesticus (Mmd) and M. m. musculus (Mmm). Comparison with cytochrome b indicates strong directional selection in Mmd Tlr4. Throughout the whole Mmd western Palaearctic region, a single variant of the ligand‐binding region is spread, encoded mainly by one dominant haplotype (71% of Mmd). In contrast, Tlr4 in Mmm is much more polymorphic with several haplotypes at intermediate frequencies. Moreover, we also found clear signals of recombination between two principal haplogroups in Mmm, and we identified eight sites under positive selection in our dataset. Our results suggest that observed differences in Tlr4 diversity may be attributed to contrasting parasite‐mediated selection acting in the two subspecies.     

Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents

Parasite-mediated selection on major histocompatibility complex (MHC) genes has mainly been explored at the intraspecific level, although many molecular studies have revealed trans-species polymorphism. Interspecific patterns of MHC diversity might reveal factors responsible for the long-term evolution of MHC polymorphism. We hypothesize that host taxa harbouring high parasite diversity should exhibit high levels of MHC genetic diversity. We test this assumption using data on rodent species and their helminth parasites compiled from the literature. Controlling for similarity due to common descent, we present evidence indicating that high helminth species richness in rodent species is associated with increased MHC class II polymorphism. Our results are consistent with the idea that parasites sharing a long-term coevolutionary history with their hosts are the agents of selection explaining MHC polymorphism.     

Assessing footprints of selection in commercial Atlantic salmon populations using microsatellite data

Relatively large rates of response to traits of economic importance have been observed in different selection experiments in salmon. Several QTL have been mapped in the salmon genome, explaining unprecedented levels of phenotypic variation. Owing to the relatively large selection intensity, individual loci may be indirectly selected, leaving molecular footprints of selection, together with increased inbreeding, as its likely relatives will share the selected loci. We used population differentiation and levels of linkage disequilibrium in chromosomes known to be harbouring QTL for body weight, infectious pancreatic necrosis resistance and infectious salmon anaemia resistance to assess the recent selection history at the genomic level in Atlantic salmon. The results clearly suggest that the marker SSA0343BSFU on chromosome 3 (body weight QTL) showed strong evidence of directional selection. It is intriguing that this marker is physically mapped to a region near the coding sequence of DVL2 , making it an ideal candidate gene to explain the rapid evolutionary response of this chromosome to selection for growth in Salmo salar. Weak evidence of diversifying selection was observed in the QTL associated with infectious pancreatic necrosis and infectious salmon anaemia resistance. Overall, this study showed that artificial selection has produced important changes in the Atlantic salmon genome, validating QTL in commercial salmon populations used for production purposes according to the recent selection history.     

Ecology and multilevel selection explain aggression in spider colonies

Progress in sociobiology continues to be hindered by abstract debates over methodology and the relative importance of within‐group vs. between‐group selection. We need concrete biological examples to ground discussions in empirical data. Recent work argued that the levels of aggression in social spider colonies are explained by group‐level adaptation. Here, we examine this conclusion using models that incorporate ecological detail while remaining consistent with kin‐ and multilevel selection frameworks. We show that although levels of aggression are driven, in part, by between‐group selection, incorporating universal within‐group competition provides a striking fit to the data that is inconsistent with pure group‐level adaptation. Instead, our analyses suggest that aggression is favoured primarily as a selfish strategy to compete for resources, despite causing lower group foraging efficiency or higher risk of group extinction. We argue that sociobiology will benefit from a pluralistic approach and stronger links between ecologically informed models and data.     

Phenotypic selection on flowering phenology and pollination efficiency traits between Primula populations with different pollinator assemblages

Floral traits have largely been attributed to phenotypic selection in plant–pollinator interactions. However, the strength of this link has rarely been ascertained with real pollinators. We conducted pollinator observations and estimated selection through female fitness on flowering phenology and floral traits between two Primula secundiflora populations. We quantified pollinator‐mediated selection by subtracting estimates of selection gradients of plants receiving supplemental hand pollination from those of plants receiving open pollination. There was net directional selection for an earlier flowering start date at populations where the dominant pollinators were syrphid flies, and flowering phenology was also subjected to stabilized quadratic selection. However, a later flowering start date was significantly selected at populations where the dominant pollinators were legitimate (normal pollination through the corolla tube entrance) and illegitimate bumblebees (abnormal pollination through nectar robbing hole which located at the corolla tube), and flowering phenology was subjected to disruptive quadratic selection. Wider corolla tube entrance diameter was selected at both populations. Furthermore, the strength of net directional selection on flowering start date and corolla tube entrance diameter was stronger at the population where the dominant pollinators were syrphid flies. Pollinator‐mediated selection explained most of the between‐population variations in the net directional selection on flowering phenology and corolla tube entrance diameter. Our results suggested the important influence of pollinator‐mediated selection on floral evolution. Variations in pollinator assemblages not only resulted in variation in the direction of selection but also the strength of selection on floral traits.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Nonlinear frequency‐dependent selection promotes long‐term coexistence between bacteria species

Negative frequency‐dependent selection (NFDS) is an important mechanism for species coexistence and for the maintenance of genetic polymorphism. Long‐term coexistence nevertheless requires NFDS interactions to be resilient to further evolution of the interacting species or genotypes. For closely related genotypes, NFDS interactions have been shown to be preserved through successive rounds of evolution in coexisting lineages. On the contrary, the evolution of NFDS interactions between distantly related species has received less attention. Here, we tracked the co‐evolution of Escherichia coli and Citrobacter freundii that initially differ in their ecological characteristics. We showed that these two bacterial species engaged in an NFDS interaction particularly resilient to further evolution: despite a very strong asymmetric rate of adaptation, their coexistence was maintained owing to an NFDS pattern where fitness increases steeply as the frequency decreases towards zero. Using a model, we showed how and why such NFDS pattern can emerge. These findings provide a robust explanation for the long‐term maintenance of species at very low frequencies.     

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole‐genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.     

Selective modes determine evolutionary rates,gene compactness and expression patterns in Brassica

It has been well documented that most nuclear protein‐coding genes in organisms can be classified into two categories: positively selected genes (PSGs) and negatively selected genes (NSGs). The characteristics and evolutionary fates of different types of genes, however, have been poorly understood. In this study, the rates of nonsynonymous substitution (K_a) and the rates of synonymous substitution (K_s) were investigated by comparing the orthologs between the two sequenced Brassica species, Brassica rapa and Brassica oleracea, and the evolutionary rates, gene structures, expression patterns, and codon bias were compared between PSGs and NSGs. The resulting data show that PSGs have higher protein evolutionary rates, lower synonymous substitution rates, shorter gene length, fewer exons, higher functional specificity, lower expression level, higher tissue‐specific expression and stronger codon bias than NSGs. Although the quantities and values are different, the relative features of PSGs and NSGs have been largely verified in the model species Arabidopsis. These data suggest that PSGs and NSGs differ not only under selective pressure (K_a/K_s), but also in their evolutionary, structural and functional properties, indicating that selective modes may serve as a determinant factor for measuring evolutionary rates, gene compactness and expression patterns in Brassica.     

Antagonism between sexual and natural selection in experimental populations of Saccharomyces cerevisiae

Trade-offs between life-history components are a central concept of evolution and ecology. Sexual and natural selection seem particularly apt to impose antagonistic selective pressures. When sex is not integrated into reproduction, as in Saccharomyces cerevisiae, natural selection can impair or even eliminate it. In this study, a genetic trade-off between the sexual and asexual phases of the yeast life cycle was suggested by sharp declines in the mating and sporulation abilities of unrelated genotypes that were propagated asexually in minimal growth medium and in mice. When sexual selection was applied to populations that had previously evolved asexually, sexual fitness increased but asexual fitness declined. No such negative correlation was observed when sexual selection was applied to an ancestral strain: sexual and asexual fitness both increased. Thus, evolutionary history affected the evolution of genetic correlations, as fitness increases in a population already well adapted to the environment were more likely to come at the expense of sexual functions.     

Segregation of male‐sterility alleles across a species boundary

Hybrid zones may serve as bridges permitting gene flow between species, including alleles influencing the evolution of breeding systems. Using greenhouse crosses, we assessed the likelihood that a hybrid zone could serve as a conduit for transfer of nuclear male‐sterility alleles between a gynodioecious species and a hermaphroditic species with very rare females in some populations. Segregation patterns in progeny of crosses between rare females of hermaphroditic Schiedea menziesii and hermaphroditic plants of gynodioecious Schiedea salicaria heterozygous at the male‐sterility locus, and between female S. salicaria and hermaphroditic plants from the hybrid zone, were used to determine whether male‐sterility was controlled at the same locus in the parental species and the hybrid zone. Segregations of females and hermaphrodites in approximately equal ratios from many of the crosses indicate that the same nuclear male‐sterility allele occurs in the parent species and the hybrid zone. These rare male‐sterility alleles in S. menziesii may result from gene flow from S. salicaria through the hybrid zone, presumably facilitated by wind pollination in S. salicaria. Alternatively, rare male‐sterility alleles might result from a reversal from gynodioecy to hermaphroditism in S. menziesii, or possibly de novo evolution of male sterility. Phylogenetic analysis indicates that some species of Schiedea have probably evolved separate sexes independently, but not in the lineage containing S. salicaria and S. menziesii. High levels of selfing and expression of strong inbreeding depression in S. menziesii, which together should favour females in populations, argue against a reversal from gynodioecy to hermaphroditism in S. menziesii.