Recombination, balancing selection and phylogenies in MHC and self-incompatibility genes

The distribution of self-incompatibility systems in angiosperms: the relationship between mating system diversity, life span, growth habit and latitude in a changing global environment

BACKGROUND AND AIMS

There is ample theoretical and experimental evidence that angiosperms harbouring self-incompatibility (SI) systems are likely to respond to global changes in unique ways relative to taxa with other mating systems. In this paper, we present an updated database on the prevalence of SI systems across angiosperms and examine the relationship between the presence of SI and latitude, biomes, life-history traits and management conditions to evaluate the potential vulnerability of SI taxa to climate change and habitat disturbance.

METHODS

We performed literature searches to identify studies that employed controlled crosses, microscopic analyses and/or genetic data to classify taxa as having SI, self-compatibility (SC), partial self-compatibility (PSC) or self-sterility (SS). Where described, the site of the SI reaction and the presence of dimorphic versus monomorphic flowers were also recorded. We then combined this database on the distribution of mating systems with information about the life span, growth habit, management conditions and geographic distribution of taxa. Information about the geographic distribution of taxa was obtained from a manually curated version of the Global Biodiversity Information Facility database, and from vegetation surveys encompassing nine biomes. We employed multinomial logit regression to assess the relationship between mating system and life-history traits, management condition, latitude and latitude-squared using self-compatible taxa as the baseline. Additionally, we employed LOESS regression to examine the relationship between the probability of SI and latitude. Finally, by summarizing information at the family level, we plotted the distribution of SI systems across angiosperms, including information about the presence of SI or dioecy and the inferred reaction site of the SI system when known, as well as the proportion of taxa in a family for which information is available.

KEY RESULTS

We obtained information about the SI status of 5686 hermaphroditic taxa, of which 55% exhibit SC and the remaining 45% harbour SI, SS or PSC. Highlights of the multinomial logit regression include that taxa with PSC have a greater odds of being short-lived (OR = 1.3) or long-lived (OR = 1.57) perennials relative to SC ones, and that SS/SI taxa (pooled) are less likely to be annuals (OR = 0.64) and more likely to be long-lived perennials (OR = 1.32). SS/SI taxa had a greater odds of being succulent (OR = 2.4) or a tree (OR = 2.05), and were less likely to be weeds (OR = 0.34). Further, we find a quadratic relationship between the probability of being self-incompatible with latitude: SI taxa were more common in the tropics, a finding that was further supported by the vegetation surveys, which showed fewer species with SS/SI in temperate and northern latitudes compared with Mediterranean and tropical biomes.

CONCLUSIONS

We conclude that in the short-term habitat fragmentation, pollinator loss and temperature increases may negatively impact plants with SI systems, particularly long-lived perennial and woody species dominant in tropical forests. In the longer term, these and other global changes are likely to select for self-compatible or partially self-compatible taxa, which, due to the apparent importance of SI as a driver of plant diversification across the angiosperm tree of life, may globally influence plant species richness.

Genomic implications of the repeated shift to self-fertilization across a species' geographic distribution

The ability to self-fertilize often varies among closely related hermaphroditic plant species, though, variation can also exist within species. In the North American Arabidopsis lyrata, the shift from self-incompatibility (SI) to selfing established in multiple regions independently, mostly since recent postglacial range expansion. This has made the species an ideal model for the investigation of the genomic basis of the breakdown of SI and its population genetic consequences. By comparing nearby selfing and outcrossing populations across the entire species' geographic distribution, we investigated variation at the self-incompatibility (S-)locus and across the genome. Furthermore, a diallel crossing experiment on one mixed-mating population was performed to gain insight into the inheritance of mating system variation. We confirmed that the breakdown of SI had evolved in several S-locus backgrounds. The diallel suggested the involvement of biparental contributions with dominance relations. Though, the population-level genome-wide association study did not single out clear-cut candidate genes but several regions with one near the S-locus. On the implication side, selfing as compared to outcrossing populations had less than half of the genomic diversity, while the number and length of runs of homozygosity (ROHs) scaled with the degree of inbreeding. Selfing populations with a history of long expansion had the longest ROHs. The results highlight that mating system shift to selfing, its genetic underpinning and the likely negative genomic consequences for evolutionary potential can be strongly interlinked with past range dynamics.

The genetic architecture of the load linked to dominant and recessive self-incompatibility alleles in Arabidopsis halleri and Arabidopsis lyrata

The long-term balancing selection acting on mating types or sex-determining genes is expected to lead to the accumulation of deleterious mutations in the tightly linked chromosomal segments that are locally 'sheltered' from purifying selection. However, the factors determining the extent of this accumulation are poorly understood. Here, we took advantage of variations in the intensity of balancing selection along a dominance hierarchy formed by alleles at the sporophytic self-incompatibility system of the Brassicaceae to compare the pace at which linked deleterious mutations accumulate among them. We first experimentally measured the phenotypic manifestation of the linked load at three different levels of the dominance hierarchy. We then sequenced and phased polymorphisms in the chromosomal regions linked to 126 distinct copies of S-alleles in two populations of Arabidopsis halleri and three populations of Arabidopsis lyrata. We find that linkage to the S-locus locally distorts phylogenies over about 10-30 kb along the chromosome. The more intense balancing selection on dominant S-alleles results in greater fixation of linked deleterious mutations, while recessive S-alleles accumulate more linked deleterious mutations that are segregating. Hence, the structure rather than the overall magnitude of the linked genetic load differs between dominant and recessive S-alleles. Our results have consequences for the long-term evolution of new S-alleles, the evolution of dominance modifiers between them, and raise the question of why the non-recombining regions of some sex and mating type chromosomes expand over evolutionary times while others, such as the S-locus of the Brassicaceae, remain restricted to small chromosomal regions.

Facilitated introgression from domestic goat into Alpine ibex at immune loci

Hybridization can result in the transfer of adaptive genetic material from one species to another, known as adaptive introgression. Bottlenecked (and hence genetically depleted) species are expected to be particularly receptive to adaptive introgression, since introgression can introduce new or previously lost adaptive genetic variation. The Alpine ibex (Capra ibex), which recently recovered from near extinction, is known to hybridize with the domestic goat (Capra aegagrus hircus), and signals of introgression previously found at the major histocompatibility complex were suggested to potentially be adaptive. Here, we combine two ancient whole genomes of Alpine ibex with 29 modern Alpine ibex genomes and 31 genomes representing six related Capra species to investigate the genome-wide patterns of introgression and confirm the potential relevance of immune loci. We identified low rates of admixture in modern Alpine ibex through various F statistics and screening for putative introgressed tracts. Further results based on demographic modelling were consistent with introgression to have occurred during the last 300 years, coinciding with the known species bottleneck, and that in each generation, 1-2 out of 100 Alpine ibex had a domestic goat parent. The putatively introgressed haplotypes were enriched at immune-related genes, where the adaptive value of alternative alleles may give individuals with otherwise depleted genetic diversity a selective advantage. While interbreeding with domestic species is a prevalent issue in species conservation, in this specific case, it resulted in putative adaptive introgression. Our findings highlight the complex interplay between hybridization, adaptive evolution, and the potential risks and benefits associated with anthropogenic influences on wild species.

Evolution of the vls Antigenic Variability Locus of the Lyme Disease Pathogen and Development of Recombinant Monoclonal Antibodies Targeting Conserved VlsE Epitopes

VlsE (variable major protein-like sequence, expressed) is an outer surface protein of the Lyme disease pathogen (Borreliella species) responsible for its within-host antigenic variation and a key diagnostic biomarker of Lyme disease. However, the high sequence variability of VlsE poses a challenge to the development of consistent VlsE-based diagnostics and therapeutics. In addition, the standard diagnostic protocols detect immunoglobins elicited by the Lyme pathogen, not the presence of the pathogen or its derived antigens. Here, we described the development of recombinant monoclonal antibodies (rMAbs) that bound specifically to conserved epitopes on VlsE. We first quantified amino-acid sequence variability encoded by the vls genes from 13 B. burgdorferi genomes by evolutionary analyses. We showed broad inconsistencies of the sequence phylogeny with the genome phylogeny, indicating rapid gene duplications, losses, and recombination at the vls locus. To identify conserved epitopes, we synthesized peptides representing five long conserved invariant regions (IRs) on VlsE. We tested the antigenicity of these five IR peptides using sera from three mammalian host species including human patients, the natural reservoir white-footed mouse (Peromyscus leucopus), and VlsE-immunized New Zealand rabbits (Oryctolagus cuniculus). The IR4 and IR6 peptides emerged as the most antigenic and reacted strongly with both the human and rabbit sera, while all IR peptides reacted poorly with sera from natural hosts. Four rMAbs binding specifically to the IR4 and IR6 peptides were identified, cloned, and purified. Given their specific recognition of the conserved epitopes on VlsE, these IR-specific rMAbs are potential novel diagnostic and research agents for direct detection of Lyme disease pathogens regardless of strain heterogeneity. IMPORTANCE Current diagnostic protocols of Lyme disease indirectly detect the presence of antibodies produced by the patient upon infection by the bacterial pathogen, not the pathogen itself. These diagnostic tests tend to underestimate early-stage bacterial infections before the patients develop robust immune responses. Further, the indirect tests do not distinguish between active or past infections by the Lyme disease bacteria in a patient sample. Here, we described novel monoclonal antibodies that have the potential to become the basis of direct and definitive diagnostic detection of the Lyme disease pathogen, regardless of its genetic heterogeneity.

Genetic diversity, evolution and selection in the major histocompatibility complex DRB and DQB loci in the family Equidae

Background

The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family.

Results

Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3.

Conclusions

The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae.

Evolution of self-incompatibility in the Brassicaceae: Lessons from a textbook example of natural selection

Self-incompatibility (SI) is a self-recognition genetic system enforcing outcrossing in hermaphroditic flowering plants and results in one of the arguably best understood forms of natural (balancing) selection maintaining genetic variation over long evolutionary times. A rich theoretical and empirical population genetics literature has considerably clarified how the distribution of SI phenotypes translates into fitness differences among individuals by a combination of inbreeding avoidance and rare-allele advantage. At the same time, the molecular mechanisms by which self-pollen is specifically recognized and rejected have been described in exquisite details in several model organisms, such that the genotype-to-phenotype map is also pretty well understood, notably in the Brassicaceae. Here, we review recent advances in these two fronts and illustrate how the joint availability of detailed characterization of genotype-to-phenotype and phenotype-to-fitness maps on a single genetic system (plant self-incompatibility) provides the opportunity to understand the evolutionary process in a unique perspective, bringing novel insight on general questions about the emergence, maintenance, and diversification of a complex genetic system.

Selection, drift, and introgression shape MHC polymorphism in lizards

The major histocompatibility complex (MHC) has long served as a model for the evolution of adaptive genetic diversity in wild populations. Pathogen-mediated selection is thought to be a main driver of MHC diversity, but it remains elusive to what degree selection shapes MHC diversity in complex biogeographical scenarios where other evolutionary processes (e.g. genetic drift and introgression) may also be acting. Here we focus on two closely related green lizard species, Lacerta trilineata and L. viridis, to address the evolutionary forces acting on MHC diversity in populations with different biogeographic structure. We characterized MHC class I exon 2 and exon 3, and neutral diversity (microsatellites), to study the relative importance of selection, drift, and introgression in shaping MHC diversity. As expected, positive selection was a significant force shaping the high diversity of MHC genes in both species. Moreover, introgression significantly increased MHC diversity in mainland populations, with a primary direction of gene flow from L. viridis to L. trilineata. Finally, we found significantly fewer MHC alleles in island populations, but maintained MHC sequence and functional diversity, suggesting that positive selection counteracted the effect of drift. Overall, our data support that different evolutionary processes govern MHC diversity in different biogeographical scenarios: positive selection occurs broadly while introgression acts in sympatry and drift when the population sizes decrease.

A hypervariable immunoglobulin superfamily member from Crassostrea gigas functions as pattern recognition receptor with opsonic activity

Immunoglobulin superfamily (IgSF), an extensive collection of proteins possessing at least one immunoglobulin-like (Ig-like) domain, performs a wide range of functions in recognition, binding or adhesion process of cells. In the present study, a cysteine-rich motif associated immunoglobulin domain containing protein (designated CgCAICP-1) was identified in Pacific oyster Crassostrea gigas. The deduced protein sequence of CgCAICP-1 contained 534 amino acidresidues, with three Ig domains which were designated as IG1, IG2 and IG3, and a cysteine-rich motif between the first and second Ig domain. The mRNA transcripts of CgCAICP-1 were highly expressed in hemocytes and up-regulated significantly (p < 0.05) after the stimulation of lipopolysaccharides (LPS), but not peptidoglycan (PGN). The recombinant CgCAICP-1 protein (rCgCAICP-1) exhibited binding activity to various pathogen-associated molecular patterns (PAMPs) including LPS, PGN, mannose (Man) and D-galactose (D-Gal), and microorganisms including Vibrio splendidus, Escherichia coli, Staphylococcus aureus, Micrococcus luteus and Pichia pastoris. The phagocytic rates of oyster hemocytes towards Gram-negative bacteria V. splendidus and Gram-positive bacteria M. luteus were significantly enhanced (p < 0.05) after pre-incubation of microbes with rCgCAICP-1. Furthermore, the transcripts of CgCAICP-1 exhibited high level of polymorphism among individuals. The ratio of nonsynonymous and synonymous distances (dN/dS) for AA'BCC'D strands of IG1 (the possible binding sites 1, pbs1) across all allelic variants was 2.09 (p < 0.05), while the ratio for the non-pbs regions was less than 1.0. The 1248 bp fragment amplified from the 5' end of CgCAICP-1 open reading frame (ORF) from 24 transcript variants could be divided artificially into seven regions of 50 elements, and all of the allelic variants might be derived from these elements by point mutation and recombination processes. These results collectively suggested that CgCAICP-1 might function as an important pattern recognition receptor (PRR) to recognize various PAMPs and facilitated the phagocytosis of oyster hemocytes towards both Gram-positive and Gram-negative bacteria. Diverse isoforms of CgCAICP-1 were generated through point mutation and recombination processes and maintained by balancing selection, which would provide a broader spectrum of interaction surface and be associated with immune resistance of oysters to infectious pathogens.

Global population structure and adaptive evolution of aflatoxin-producing fungi

Aflatoxins produced by several species in Aspergillus section Flavi are a significant problem in agriculture and a continuous threat to human health. To provide insights into the biology and global population structure of species in section Flavi, a total of 1,304 isolates were sampled across six species (A. flavus, A. parasiticus, A. nomius, A. caelatus, A. tamarii, and A. alliaceus) from single fields in major peanut‐growing regions in Georgia (USA), Australia, Argentina, India, and Benin (Africa). We inferred maximum‐likelihood phylogenies for six loci, both combined and separately, including two aflatoxin cluster regions (aflM/alfN and aflW/aflX) and four noncluster regions (amdS, trpC, mfs and MAT), to examine population structure and history. We also employed principal component and STRUCTURE analysis to identify genetic clusters and their associations with six different categories (geography, species, precipitation, temperature, aflatoxin chemotype profile, and mating type). Overall, seven distinct genetic clusters were inferred, some of which were more strongly structured by G chemotype diversity than geography. Populations of A. flavus S in Benin were genetically distinct from all other section Flavi species for the loci examined, which suggests genetic isolation. Evidence of trans‐speciation within two noncluster regions, whereby A. flavus S_BG strains from Australia share haplotypes with either A. flavus or A. parasiticus, was observed. Finally, while clay soil and precipitation may influence species richness in Aspergillus section Flavi, other region‐specific environmental and genetic parameters must also be considered.

A genome-wide scan for genes under balancing selection in Drosophila melanogaster

BACKGROUND

In the history of population genetics balancing selection has been considered as an important evolutionary force, yet until today little is known about its abundance and its effect on patterns of genetic diversity. Several well-known examples of balancing selection have been reported from humans, mice, plants, and parasites. However, only very few systematic studies have been carried out to detect genes under balancing selection. We performed a genome scan in Drosophila melanogaster to find signatures of balancing selection in a derived (European) and an ancestral (African) population. We screened a total of 34 genomes searching for regions of high genetic diversity and an excess of SNPs with intermediate frequency.

RESULTS

In total, we found 183 candidate genes: 141 in the European population and 45 in the African one, with only three genes shared between both populations. Most differences between both populations were observed on the X chromosome, though this might be partly due to false positives. Functionally, we find an overrepresentation of genes involved in neuronal development and circadian rhythm. Furthermore, some of the top genes we identified are involved in innate immunity.

CONCLUSION

Our results revealed evidence of genes under balancing selection in European and African populations. More candidate genes have been found in the European population. They are involved in several different functions.

Disentangling host, pathogen, and environmental determinants of a recently emerged wildlife disease: lessons from the first 15 years of amphibian chytridiomycosis research

The amphibian fungal disease chytridiomycosis, which affects species across all continents, recently emerged as one of the greatest threats to biodiversity. Yet, many aspects of the basic biology and epidemiology of the pathogen, Batrachochytrium dendrobatidis (Bd), are still unknown, such as when and from where did Bd emerge and what is its true ecological niche? Here, we review the ecology and evolution of Bd in the Americas and highlight controversies that make this disease so enigmatic. We explore factors associated with variance in severity of epizootics focusing on the disease triangle of host susceptibility, pathogen virulence, and environment. Reevaluating the causes of the panzootic is timely given the wealth of data on Bd prevalence across hosts and communities and the recent discoveries suggesting co-evolutionary potential of hosts and Bd. We generate a new species distribution model for Bd in the Americas based on over 30,000 records and suggest a novel future research agenda. Instead of focusing on pathogen "hot spots," we need to identify pathogen "cold spots" so that we can better understand what limits the pathogen's distribution. Finally, we introduce the concept of "the Ghost of Epizootics Past" to discuss expected patterns in postepizootic host communities.

The effects of recombination, mutation and selection on the evolution of the Rp1 resistance genes in grasses

Plant immune genes, or resistance genes, are involved in a co-evolutionary arms race with a diverse range of pathogens. In agronomically important grasses, such R genes have been extensively studied because of their role in pathogen resistance and in the breeding of resistant cultivars. In this study, we evaluate the importance of recombination, mutation and selection on the evolution of the R gene complex Rp1 of Sorghum, Triticum, Brachypodium, Oryza and Zea. Analyses show that recombination is widespread, and we detected 73 independent instances of sequence exchange, involving on average 1567 of 4692 nucleotides analysed (33.4%). We were able to date 24 interspecific recombination events and found that four occurred postspeciation, which suggests that genetic introgression took place between different grass species. Other interspecific events seemed to have been maintained over long evolutionary time, suggesting the presence of balancing selection. Significant positive selection (i.e. a relative excess of nonsynonymous substitutions (dN /dS >1)) was detected in 17-95 codons (0.42-2.02%). Recombination was significantly associated with areas with high levels of polymorphism but not with an elevated dN /dS ratio. Finally, phylogenetic analyses show that recombination results in a general overestimation of the divergence time (mean = 14.3%) and an alteration of the gene tree topology if the tree is not calibrated. Given that the statistical power to detect recombination is determined by the level of polymorphism of the amplicon as well as the number of sequences analysed, it is likely that many studies have underestimated the importance of recombination relative to the mutation rate.

Long-Term Balancing Selection in LAD1 Maintains a Missense Trans-Species Polymorphism in Humans, Chimpanzees, and Bonobos

Balancing selection maintains advantageous genetic and phenotypic diversity in populations. When selection acts for long evolutionary periods selected polymorphisms may survive species splits and segregate in present-day populations of different species. Here, we investigate the role of long-term balancing selection in the evolution of protein-coding sequences in the Homo-Pan clade. We sequenced the exome of 20 humans, 20 chimpanzees, and 20 bonobos and detected eight coding trans-species polymorphisms (trSNPs) that are shared among the three species and have segregated for approximately 14 My of independent evolution. Although the majority of these trSNPs were found in three genes of the major histocompatibility locus cluster, we also uncovered one coding trSNP (rs12088790) in the gene LAD1. All these trSNPs show clustering of sequences by allele rather than by species and also exhibit other signatures of long-term balancing selection, such as segregating at intermediate frequency and lying in a locus with high genetic diversity. Here, we focus on the trSNP in LAD1, a gene that encodes for Ladinin-1, a collagenous anchoring filament protein of basement membrane that is responsible for maintaining cohesion at the dermal-epidermal junction; the gene is also an autoantigen responsible for linear IgA disease. This trSNP results in a missense change (Leucine257Proline) and, besides altering the protein sequence, is associated with changes in gene expression of LAD1.

Positive selection drives the evolution of a major histocompatibility complex gene in an endangered Mexican salamander species complex

Immune gene evolution can be critical to species survival in the face of infectious disease. In particular, polymorphism in the genes of the major histocompatibility complex (MHC) helps vertebrates combat novel and diverse pathogens by increasing the number of pathogen-derived proteins that can initiate the host's acquired immune response. In this study, we used a combination of presumably adaptive and neutral markers to investigate MHC evolution in populations of five salamander species within the Ambystoma velasci complex, a group consisting of 15 recently diverged species, several of which are endangered. We isolated 31 unique MHC class II β alleles from 75 total individuals from five species in this complex. MHC heterozygosity was significantly lower than expected for all five species, and we found no clear relationship between number of MHC alleles and species range, life history, or level of heterozygosity. We inferred a phylogeny representing the evolutionary history of Ambystoma MHC, with which we found signatures of positive selection on the overall gene, putative peptide-binding residues, and allelic lineages. We identified several instances of trans-species polymorphism, a hallmark of balancing selection observed in other groups of closely related species. In contrast, we did not detect comparable allelic diversity or signatures of selection on neutral loci. Additionally, we identified 17 supertypes among the 44 unique Ambystoma alleles, indicating that these sequences may encode functionally distinct MHC variants. We therefore have strong evidence that positive selection is a major evolutionary force driving patterns of MHC polymorphism in this recently radiated species complex.

Lineage sorting in apes

Recombination allows different parts of the genome to have different genealogical histories. When a species splits in two, allelic lineages sort into the two descendant species, and this lineage sorting varies along the genome. If speciation events are close in time, the lineage sorting process may be incomplete at the second speciation event and lead to gene genealogies that do not match the species phylogeny. We review different recent approaches to model lineage sorting along the genome and show how it is possible to learn about population sizes, natural selection, and recombination rates in ancestral species from application of these models to genome alignments of great ape species.

Advantageous diversity maintained by balancing selection in humans

Most human polymorphisms are neutral or slightly deleterious, but some genetic variation is advantageous and maintained in populations by balancing selection. Considered a rarity and overlooked for years, balanced polymorphisms have recently received renewed attention with several lines of evidence showing their relevance in human evolution. From theoretical work on its role in adaptation to empirical studies that identify its targets, recent developments have showed that balancing selection is more prevalent than previously thought. Here we review these developments and discuss their implications in our understanding of the influence of balancing selection in human evolution. We also review existing evidence on the biological functions that benefit most from advantageous diversity, and the functional consequences of these variants. Overall, we argue that balancing selection must be considered an important selective force in human evolution.

Gradual molecular evolution of a sex determination switch through incomplete penetrance of femaleness

Some genes regulate phenotypes that are either present or absent. They are often important regulators of developmental switches and are involved in morphological evolution. We have little understanding of the molecular mechanisms by which these absence/presence gene functions have evolved, because the phenotype and fitness of molecular intermediate forms are unknown. Here, we studied the sex-determining switch of 14 natural sequence variants of the csd gene among 76 genotypes of the honeybee (Apis mellifera). Heterozygous genotypes (different specificities) of the csd gene determine femaleness, while hemizygous genotypes (single specificity) determine maleness. Homozygous genotypes of the csd gene (same specificity) are lethal. We found that at least five amino acid differences and length variation between Csd specificities in the specifying domain (PSD) were sufficient to regularly induce femaleness. We estimated that, on average, six pairwise amino acid differences evolved under positive selection. We also identified a natural evolutionary intermediate that showed only three amino acid length differences in the PSD relative to its parental allele. This genotype showed an intermediate fitness because it implemented lethality regularly and induced femaleness infrequently (i.e., incomplete penetrance). We suggest incomplete penetrance as a mechanism through which new molecular switches can gradually and adaptively evolve.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation

BACKGROUND

Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1).

RESULTS

Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases.

CONCLUSIONS

Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

The effect of fluctuating selection on the genealogy at a linked site

The genealogical consequences of temporally fluctuating selection at linked neutrally-evolving sites are studied using coalescent processes structured by genetic backgrounds. Surprisingly, although between-generation fluctuating selection and within-generation fecundity variance polymorphism lead to indistinguishable allele frequency dynamics at the selected site, I show that these two scenarios affect the genealogical structure of the population in distinctive ways. Whereas coalescence times are elevated at the selected site when heterozygotes have lower within-generation fecundity variance than either homozygote, fluctuating selection typically depresses the depth of the genealogy at tightly-linked sites. More importantly, these results indicate that fluctuating selection will have a characteristic signature combining an excess of selected variation at the affected site and a slight reduction in neutral variation at tightly linked sites. This observation suggests that it may be possible to distinguish balanced polymorphisms maintained by environmental variation from those maintained by heterozygote advantage.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

Recent and ancient signature of balancing selection around the S-locus in Arabidopsis halleri and A. lyrata

Balancing selection can maintain different alleles over long evolutionary times. Beyond this direct effect on the molecular targets of selection, balancing selection is also expected to increase neutral polymorphism in linked genome regions, in inverse proportion to their genetic map distances from the selected sites. The genes controlling plant self-incompatibility are subject to one of the strongest forms of balancing selection, and they show clear signatures of balancing selection. The genome region containing those genes (the S-locus) is generally described as nonrecombining, and the physical size of the region with low recombination has recently been established in a few species. However, the size of the region showing the indirect footprints of selection due to linkage to the S-locus is only roughly known. Here, we improved estimates of this region by surveying synonymous polymorphism and estimating recombination rates at 12 flanking region loci at known physical distances from the S-locus region boundary, in two closely related self-incompatible plants Arabidopsis halleri and A. lyrata. In addition to studying more loci than previous studies and using known physical distances, we simulated an explicit demographic scenario for the divergence between the two species, to evaluate the extent of the genomic region whose diversity departs significantly from neutral expectations. At the closest flanking loci, we detected signatures of both recent and ancient indirect effects of selection on the S-locus flanking genes, finding ancestral polymorphisms shared by both species, as well as an excess of derived mutations private to either species. However, these effects are detected only in a physically small region, suggesting that recombination in the flanking regions is sufficient to quickly break up linkage disequilibrium with the S-locus. Our approach may be useful for distinguishing cases of ancient versus recently evolved balancing selection in other systems.

A sequence-based approach demonstrates that balancing selection in classical human leukocyte antigen (HLA) loci is asymmetric

Balancing selection has maintained human leukocyte antigen (HLA) allele diversity, but it is unclear whether this selection is symmetric (all heterozygotes are comparable and all homozygotes are comparable in terms of fitness) or asymmetric (distinct heterozygote genotypes display greater fitness than others). We tested the hypothesis that HLA is under asymmetric balancing selection in populations by estimating allelic branch lengths from genetic sequence data encoding peptide-binding domains. Significant deviations indicated changes in the ratio of terminal to internal branch lengths. Such deviations could arise even if no individual alleles present a strikingly altered branch length (e.g. if there is an overall distortion, with all or many terminal branches being longer than expected). DQ and DP loci were also analyzed as haplotypes. Using allele frequencies for 419 distinct populations in 10 geographical regions, we examined population differentiation in alleles within and between regions, and the relationship between allelic branch length and frequency. The strongest evidence for asymmetrical balancing selection was observed for HLA-DRB1, HLA-B and HLA-DPA1, with significant deviation (P ≤ 1.1 × 10(-4)) in about half of the populations. There were significant results at all loci except HLA-DQB1/DQA1. We observed moderate genetic variation within and between geographic regions, similar to the rest of the genome. Branch length was not correlated with allele frequency. In conclusion, sequence data suggest that balancing selection in HLA is asymmetric (some heterozygotes enjoy greater fitness than others). Because HLA polymorphism is crucial for pathogen resistance, this may manifest as a frequency-dependent selection with fluctuation in the fitness of specific heterozygotes over time.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this genomic region.

Creation and maintenance of variation in allorecognition Loci: molecular analysis in various model systems

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated conspecifics. Effective allorecognition systems are critical to the survival of organisms; they prevent inbreeding and facilitate fusions between close relatives. Where the loci governing allorecognition outcomes have been identified, the corresponding proteins often exhibit exceptional polymorphism. Two important questions about this polymorphism remain unresolved: how is it created, and how is it maintained. Because the genetic bases of several allorecognition systems have now been identified, including alr1 and alr2 in Hydractinia, fusion histocompatibility in Botryllus, the het (vic) loci in fungi, tgrB1 and tgrC1 in Dictyostelium, and self-incompatibility (SI) loci in several plant families, we are now poised to achieve a clearer understanding of how these loci evolve. In this review, we summarize what is currently known about the evolution of allorecognition loci, highlight open questions, and suggest future directions.

Evolutionary genetics of an S-like polymorphism in Papaveraceae with putative function in self-incompatibility

BACKGROUND

Papaver rhoeas possesses a gametophytic self-incompatibility (SI) system not homologous to any other SI mechanism characterized at the molecular level. Four previously published full length stigmatic S-alleles from the genus Papaver exhibited remarkable sequence divergence, but these studies failed to amplify additional S-alleles despite crossing evidence for more than 60 S-alleles in Papaver rhoeas alone.

METHODOLOGY/PRINCIPAL FINDINGS

Using RT-PCR we identified 87 unique putative stigmatic S-allele sequences from the Papaveraceae Argemone munita, Papaver mcconnellii, P. nudicuale, Platystemon californicus and Romneya coulteri. Hand pollinations among two full-sib families of both A. munita and P. californicus indicate a strong correlation between the putative S-genotype and observed incompatibility phenotype. However, we also found more than two S-like sequences in some individuals of A. munita and P. californicus, with two products co-segregating in both full-sib families of P. californicus. Pairwise sequence divergence estimates within and among taxa show Papaver stigmatic S-alleles to be the most variable with lower divergence among putative S-alleles from other Papaveraceae. Genealogical analysis indicates little shared ancestral polymorphism among S-like sequences from different genera. Lack of shared ancestral polymorphism could be due to long divergence times among genera studied, reduced levels of balancing selection if some or all S-like sequences do not function in incompatibility, population bottlenecks, or different levels of recombination among taxa. Preliminary estimates of positive selection find many sites under selective constraint with a few undergoing positive selection, suggesting that self-recognition may depend on amino acid substitutions at only a few sites.

CONCLUSIONS/SIGNIFICANCE

Because of the strong correlation between genotype and SI phenotype, sequences reported here represent either functional stylar S-alleles, tightly linked paralogs of the S-locus or a combination of both. The considerable complexity revealed in this study shows we have much to learn about the evolutionary dynamics of self-incompatibility systems.

Origin and diversification dynamics of self-incompatibility haplotypes

Self-incompatibility (SI) is a genetic system found in some hermaphrodite plants. Recognition of pollen by pistils expressing cognate specificities at two linked genes leads to rejection of self pollen and pollen from close relatives, i.e., to avoidance of self-fertilization and inbred matings, and thus increased outcrossing. These genes generally have many alleles, yet the conditions allowing the evolution of new alleles remain mysterious. Evolutionary changes are clearly necessary in both genes, since any mutation affecting only one of them would result in a nonfunctional self-compatible haplotype. Here, we study diversification at the S-locus (i.e., a stable increase in the total number of SI haplotypes in the population, through the incorporation of new SI haplotypes), both deterministically (by investigating analytically the fate of mutations in an infinite population) and by simulations of finite populations. We show that the conditions allowing diversification are far less stringent in finite populations with recurrent mutations of the pollen and pistil genes, suggesting that diversification is possible in a panmictic population. We find that new SI haplotypes emerge fastest in populations with few SI haplotypes, and we discuss some implications for empirical data on S-alleles. However, allele numbers in our simulations never reach values as high as observed in plants whose SI systems have been studied, and we suggest extensions of our models that may reconcile the theory and data.

Natural selection coupled with intragenic recombination shapes diversity patterns in the major histocompatibility complex class II genes of the giant panda

Ample variations of the major histocompatibility complex (MHC) genes are essential for vertebrates to adapt to various environmental conditions. In this study, we investigated the genetic variations and evolutionary patterns of seven functional MHC class II genes (one DRA, two DRB, two DQA, and two DQB) of the giant panda. The results showed the presence of two monomorphic loci (DRA and DQB2) and five polymorphic loci with different numbers of alleles (seven at DRB1, six at DRB3, seven at DQA1, four at DQA2, six at DQB1). The presence of balancing selection in the giant panda was supported by the following pieces of evidence: (1) The observed heterozygosity was higher than expected. (2) Amino acid heterozygosity was significantly higher at antigen-binding sites (ABS) compared with non-ABS sequences. (3) The selection parameter omega (d(N)/d(S)) was significantly higher at ABS compared with non-ABS sequences. (4) Approximately 95.45% of the positively selected codons (P>0.95) were located at or adjacent to an ABS. Furthermore, this study showed that (1) The Qinling subspecies exhibited high omega values across each locus (all >1), supporting its extensive positive selection. (2) The Sichuan subspecies displayed small omega at DRB1 (omega<0.72) and DQA2 (omega<0.48), suggesting that these sites underwent strong purifying selection. (3) Intragenic recombination was detected in DRB1, DQA1, and DQB1. The molecular diversity in classic Aime-MHC class II genes implies that the giant panda had evolved relatively abundant variations in its adaptive immunity along the history of host-pathogen co-evolution. Collectively, these findings indicate that natural selection accompanied by recombination drives the contrasting diversity patterns of the MHC class II genes between the two studied subspecies of giant panda.

Evolutionary patterns at the RNase based gametophytic self - incompatibility system in two divergent Rosaceae groups (Maloideae and Prunus)

Background

Within Rosaceae, the RNase based gametophytic self-incompatibility (GSI) system has been studied at the molecular level in Maloideae and Prunus species that have been diverging for, at least, 32 million years. In order to understand RNase based GSI evolution within this family, comparative studies must be performed, using similar methodologies.

Result

It is here shown that many features are shared between the two species groups such as levels of recombination at the S-RNase (the S-pistil component) gene, and the rate at which new specificities arise. Nevertheless, important differences are found regarding the number of ancestral lineages and the degree of specificity sharing between closely related species. In Maloideae, about 17% of the amino acid positions at the S-RNase protein are found to be positively selected, and they occupy about 30% of the exposed protein surface. Positively selected amino acid sites are shown to be located on either side of the active site cleft, an observation that is compatible with current models of specificity determination. At positively selected amino acid sites, non-conservative changes are almost as frequent as conservative changes. There is no evidence that at these sites the most drastic amino acid changes may be more strongly selected.

Conclusions

Many similarities are found between the GSI system of Prunus and Maloideae that are compatible with the single origin hypothesis for RNase based GSI. The presence of common features such as the location of positively selected amino acid sites and lysine residues that may be important for ubiquitylation, raise a number of issues that, in principle, can be experimentally addressed in Maloideae. Nevertheless, there are also many important differences between the two Rosaceae GSI systems. How such features changed during evolution remains a puzzling issue.

Characterization of phylogenetic networks with NetTest

Background

Typical evolutionary events like recombination, hybridization or gene transfer make necessary the use of phylogenetic networks to properly depict the evolution of DNA and protein sequences. Although several theoretical classes have been proposed to characterize these networks, they make stringent assumptions that will likely not be met by the evolutionary process. We have recently shown that the complexity of simulated networks is a function of the population recombination rate, and that at moderate and large recombination rates the resulting networks cannot be categorized. However, we do not know whether these results extend to networks estimated from real data.

Results

We introduce a web server for the categorization of explicit phylogenetic networks, including the most relevant theoretical classes developed so far. Using this tool, we analyzed statistical parsimony phylogenetic networks estimated from ~5,000 DNA alignments, obtained from the NCBI PopSet and Polymorphix databases. The level of characterization was correlated to nucleotide diversity, and a high proportion of the networks derived from these data sets could be formally characterized.

Conclusions

We have developed a public web server, NetTest (freely available from the software section at http://darwin.uvigo.es), to formally characterize the complexity of phylogenetic networks. Using NetTest we found that most statistical parsimony networks estimated with the program TCS could be assigned to a known network class. The level of network characterization was correlated to nucleotide diversity and dependent upon the intra/interspecific levels, although no significant differences were detected among genes. More research on the properties of phylogenetic networks is clearly needed.

RNase-based gametophytic self-incompatibility evolution: Questioning the hypothesis of multiple independent recruitments of the S-pollen gene

Multiple independent recruitments of the S-pollen component (always an F-box gene) during RNase-based gametophytic self-incompatibility evolution have recently been suggested. Therefore, different mechanisms could be used to achieve the rejection of incompatible pollen in different plant families. This hypothesis is, however, mainly based on the interpretation of phylogenetic analyses, using a small number of divergent nucleotide sequences. In this work we show, based on a large collection of F-box S-like sequences, that the inferred relationship of F-box S-pollen and F-box S-like sequences is dependent on the sequence alignment software and phylogenetic method used. Thus, at present, it is not possible to address the phylogenetic relationship of F-box S-pollen and S-like sequences from different plant families. In Petunia and Malus/Pyrus the putative S-pollen gene(s) show(s) variability patterns different than expected for an S-pollen gene, raising the question of false identification. Here we show that in Petunia, the unexpected features of the putative S-pollen gene are not incompatible with this gene's being the S-pollen gene. On the other hand, it is very unlikely that the Pyrus SFBB-gamma gene is involved in specificity determination.

Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen-mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d(n)/d(s) ratio and the presence of trans-species polymorphisms suggest that a strong long-term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D-loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000-30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.

Multiple advantageous amino acid variants in the NAT2 gene in human populations

BACKGROUND

Genetic variation at NAT2 has been long recognized as the cause of differential ability to metabolize a wide variety of drugs of therapeutic use. Here, we explore the pattern of genetic variation in 12 human populations that significantly extend the geographic range and resolution of previous surveys, to test the hypothesis that different dietary regimens and lifestyles may explain inter-population differences in NAT2 variation.

METHODOLOGY/PRINCIPAL FINDINGS

The entire coding region was resequenced in 98 subjects and six polymorphic positions were genotyped in 150 additional subjects. A single previously undescribed variant was found (34T>C; 12Y>H). Several aspects of the data do not fit the expectations of a neutral model, as assessed by coalescent simulations. Tajima's D is positive in all populations, indicating an excess of intermediate alleles. The level of between-population differentiation is low, and is mainly accounted for by the proportion of fast vs. slow acetylators. However, haplotype frequencies significantly differ across groups of populations with different subsistence.

CONCLUSIONS/SIGNIFICANCE

Data on the structure of haplotypes and their frequencies are compatible with a model in which slow-causing variants were present in widely dispersed populations before major shifts to pastoralism and/or agriculture. In this model, slow-causing mutations gained a selective advantage in populations shifting from hunting-gathering to pastoralism/agriculture. We suggest the diminished dietary availability of folates resulting from the nutritional shift, as the possible cause of the fitness increase associated to haplotypes carrying mutations that reduce enzymatic activity.

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

Effects of recombination on hitchhiking diversity in the Brassica self-incompatibility locus complex

In self-incompatibility, a number of S haplotypes are maintained by frequency-dependent selection, which results in trans-specific S haplotypes. The region of several kilobases (approximately 40-60 kb) from SP6 to SP2, including self-incompatibility-related genes and some adjacent genes in Brassica rapa, has high nucleotide diversity due to the hitchhiking effect, and therefore we call this region the "S-locus complex." Recombination in the S-locus complex is considered to be suppressed. We sequenced regions of >50 kb of the S-locus complex of three S haplotypes in B. rapa and found higher nucleotide diversity in intergenic regions than in coding regions. Two highly similar regions of >10 kb were found between BrS-8 and BrS-46. Phylogenetic analysis using trans-specific S haplotypes (called interspecific pairs) of B. rapa and B. oleracea suggested that recombination reduced the nucleotide diversity in these two regions and that the genes not involved in self-incompatibility in the S-locus complex and the kinase domain, but not the S domain, of SRK have also experienced recombination. Recombination may reduce hitchhiking diversity in the S-locus complex, whereas the region from the S domain to SP11 would disfavor recombination.

Different positively selected sites at the gametophytic self-incompatibility pistil S-RNase gene in the Solanaceae and Rosaceae (Prunus, Pyrus, and Malus)

In this work we perform a comparative study on the location of positively selected sites (those likely responsible for defining specificity differences) at the S-RNase gene, the pistil component of the gametophytic self-incompatibility system. For Plantaginaceae and Rosaceae (Prunus and Pyrus/Malus) this is the first study of this kind. A clear sign of positive selection was observed for 13, 17, and 27 amino acid sites in Solanaceae, Prunus, and Pyrus/Malus, respectively, using two different methodologies. In Plantaginaceae no clear positively selected sites were identified. Possible reasons for this result are discussed. Indirect experimental evidence suggests that the identified positively selected amino acid sites play a role in specificity determination. The percentage of positively selected sites is similar in Solanaceae and Rosaceae but the location of those sites is different.

Genomic variation in rice: genesis of highly polymorphic linkage blocks during domestication

Genomic regions that are unusually divergent between closely related species or racial groups can be particularly informative about the process of speciation or the operation of natural selection. The two sequenced genomes of cultivated Asian rice, Oryza sativa, reveal that at least 6% of the genomes are unusually divergent. Sequencing of ten unlinked loci from the highly divergent regions consistently identified two highly divergent haplotypes with each locus in nearly complete linkage disequilibrium among 25 O. sativa cultivars and 35 lines from six wild species. The existence of two highly divergent haplotypes in high divergence regions in species from all geographical areas (Africa, Asia, and Oceania) was in contrast to the low polymorphism and low linkage disequilibrium that were observed in other parts of the genome, represented by ten reference loci. While several natural processes are likely to contribute to this pattern of genomic variation, domestication may have greatly exaggerated the trend. In this hypothesis, divergent haplotypes that were adapted to different geographical and ecological environments migrated along with humans during the development of domesticated varieties. If true, these high divergence regions of the genome would be enriched for loci that contribute to the enormous range of phenotypic variation observed among domesticated breeds.

Ancient polymorphism and functional variation in the primate MHC-DQA1 5' cis-regulatory region

Precise regulation of MHC gene expression is critical to vertebrate immune surveillance and response. Polymorphisms in the 5' proximal promoter region of the human class II gene HLA-DQA1 have been shown to influence its transcriptional regulation and may contribute to the pathogenesis of autoimmune diseases. We investigated the evolutionary history of this cis-regulatory region by sequencing the DQA1 5' proximal promoter region in eight nonhuman primate species. We observed unexpectedly high levels of sequence variation and multiple strong signatures of balancing selection in this region. Specifically, the considerable DQA1 promoter region diversity was characterized by abundant shared (or trans-species) polymorphism and a pronounced lack of fixed differences between species. The majority of transcription factor binding sites in the DQA1 promoter region were polymorphic within species, and these binding site polymorphisms were commonly shared among multiple species despite evidence for negative selection eliminating a significant fraction of binding site mutations. We assessed the functional consequences of intraspecific promoter region diversity using a cell line-based reporter assay and detected significant differences among baboon DQA1 promoter haplotypes in their ability to drive transcription in vitro. The functional differentiation of baboon promoter haplotypes, together with the significant deviations from neutral sequence evolution, suggests a role for balancing selection in the evolution of DQA1 transcriptional regulation in primates.

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.

Modified Hudson-Kreitman-Aguade test and two-dimensional evaluation of neutrality tests

There are a number of polymorphism-based statistical tests of neutrality, but most of them focus on either the amount or the pattern of polymorphism. In this article, a new test called the two-dimensional (2D) test is developed. This test evaluates a pair of summary statistics in a two-dimensional field. One statistic should summarize the pattern of polymorphism, while the other could be a measure of the level of polymorphism. For the latter summary statistic, the polymorphism-divergence ratio is used following the idea of the Hudson-Kreitman-Aguadé (HKA) test. To incorporate the HKA test in the 2D test, a summary statistic-based version of the HKA test is developed such that the polymorphism-divergence ratio at a particular region of interest is examined if it is consistent with the average of those in other independent regions.

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.

Variability patterns and positively selected sites at the gametophytic self-incompatibility pollen SFB gene in a wild self-incompatible Prunus spinosa (Rosaceae) population

Current models for the generation of new gametophytic self-incompatibility specificities require that neutral variability segregates within specificity classes. Furthermore, one of the models predicts greater ratios of nonsynonymous to synonymous substitutions in pollen than in pistil specificity genes. All models assume that new specificities arise by mutation only. To test these models, 21 SFB (the pollen S-locus) alleles from a wild Prunus spinosa (Rosaceae) population were obtained. For seven of these, the corresponding S-haplotype was also characterized. The SFB data set was also used to identify positively selected sites. Those sites are likely to be the ones responsible for defining pollen specificities. Of the 23 sites identified as being positively selected, 21 are located in the variable (including a new region described here) and hypervariable regions. Little variability is found within specificity classes. There is no evidence for selective sweeps being more frequent in pollen than in pistil specificity genes. The S-RNase and the SFB genes have only partially correlated evolutionary histories. None of the models is compatible with the variability patterns found in the SFB and the S-haplotype data.

Can a genome change its (hot)spots?

A new study by Jeffreys et al. shows that the rate of recombination in recombination hotspots in humans is not constant through time. This observation adds weight to the idea that hotspots are transient on evolutionary timescales. However, questions remain as to what controls their evolution and how these rapid changes influence broad-scale rates of recombination.

Population genetics of CAPN10 and GPR35: implications for the evolution of type 2 diabetes variants

A positional cloning study of type 2 diabetes in Mexican Americans identified a region, termed "NIDDM1," on chromosome 2q37 with significant linkage evidence. Haplotype combinations at the calpain-10 gene (CAPN10) within this region were shown to increase diabetes risk in several populations. On the basis of the thrifty genotype hypothesis, variants that increase susceptibility to type 2 diabetes under modern lifestyle conditions provided a survival advantage in past environments by increasing the efficiency of energy use and storage. Here, our goal is to make inferences about the evolutionary forces shaping variation in genes in the NIDDM1 region and to investigate the population genetics models that may underlie the thrifty genotype hypothesis. To this end, we surveyed sequence variation in CAPN10 and in an adjacent gene, G-protein-coupled receptor 35 (GPR35), in four population samples from different ethnic groups. These data revealed two distinct deviations from the standard neutral model in CAPN10, whereas GPR35 variation was largely consistent with neutrality. CAPN10 showed a significant deficit of variation in the haplotype class defined by the derived allele at SNP44, a polymorphism that is significantly associated with diabetes in meta-analysis studies. This suggests that this haplotype class was quickly driven to high frequency by positive natural selection. Interestingly, the derived allele at SNP44 is protective against diabetes. CAPN10 also showed a local excess of polymorphism and linkage disequilibrium decay in intron 13. Simulations show that this pattern may be explained by long-standing balancing selection that maintains multiple selected alleles. Alternatively, it is possible that the local mutation and recombination rates changed since the divergence of human and chimpanzee; this scenario does not require the action of natural selection on intron 13 variation.

Inheritance and dominance of self-incompatibility alleles in polyploid Arabidopsis lyrata

Natural populations of diploid Arabidopsis lyrata exhibit the sporophytic type of self-incompatibility system characteristic of Brassicaceae, in which complicated dominance interactions among alleles in the diploid parent determine self-recognition phenotypes of both pollen and stigma. The purpose of this study was to investigate how polyploidy affects this already complex system. One tetraploid population (Arabidopsis lyrata ssp kawasakiana from Japan) showed complete self-compatibility and produced viable selfed progeny for at least three generations subsequent to field collection. In contrast, individuals from a second tetraploid population (A. lyrata ssp petraea from Austria) were strongly self-incompatible (SI). Segregation of SI genotypes in this population followed Mendelian patterns based on a tetrasomic model of inheritance, with two to four alleles per individual, independent segregation of alleles, and little evidence of dosage effects of alleles found in multiple copies. Similar to results from diploids, anomalous compatibility patterns involving particular combinations of individuals occurred at a low frequency in the tetraploids, suggesting altered dominance in certain genetic backgrounds that could be due to the influence of a modifier locus. Overall, dominance relationships among S-alleles in self-incompatible tetraploid families were remarkably similar to those in related diploids, suggesting that this very important and complicated locus has not undergone extensive modification subsequent to polyploidization.

Variable selection pressures across lineages in Trichodesmium and related cyanobacteria based on the heterocyst differentiation protein gene hetR

Due to the irreversible inhibition of nitrogenase by O2, N2 fixation is incompatible with the oxygenic photosynthesis of cyanobacteria. These organisms have therefore evolved various strategies for growing diazotrophically. One group of N2-fixing cyanobacteria has specialized cells, heterocysts, which contain the nitrogenase, lack the oxygenic photosystem II, and are virtually anoxic inside as the result of respiratory activity and a thick glycolipid cell wall. The hetR gene encodes a serine protease which is thought to be involved in the regulation of heterocyst development and in DNA binding. Although hetR is also present in many non-heterocystous N2-fixing cyanobacteria, its function in these organisms is unknown. In this study, hetR sequences of the N2-fixing, non-heterocystous cyanobacterium Trichodesmium spp. and related genera were examined for signatures of selection. In parsimony- or distance-based hetR phylogenies, the filamentous non-heterocystous cyanobacteria Symploca sp. and Leptolyngbya sp. were closest to Trichodesmium sp. However, accommodating molecular attributes of hetR such as nucleotide frequencies and rate heterogeneity in phylogenetic analyses suggested that many other genera could not be excluded as sister taxa of Trichodesmium. Maximum likelihood analysis of the dN/dS ratio (omega) showed that-irrespective of the use of Symploca, Leptolyngbya, or more distant taxa as an outgroup-the lineage between an outgroup and Trichodesmium (omega1=0.02-0.05) and a lineage leading to Trichodesmium erythraeum (omega1=0.02) were under much stronger purifying selection than the other lineages in Trichodesmium (omega0=0.13-0.32). Although the results from the maximum likelihood analyses are most trustworthy because of codon usage bias in Trichodesmium, the results from a simpler tree-based McDonald-Kreitman test were in general agreement. Due to their quite different assumptions, the combination of these two methods of analysis circumvents multiple testing which, in general, is problematic when using branch models. Although the causal selective forces underlying the substitution patterns in hetR have not yet been identified, these findings parallel the variety of physiological, molecular, and behavioral differences in cyanobacteria related to N2 fixation. The heterogeneity of selection pressures in Trichodesmium is more surprising, because multiple adaptation mechanisms have not been described in this genus.