Neutral behavior of shared polymorphism

Trans-Specific Polymorphisms Between Cryptic Daphnia Species Affect Fitness and Behavior

Shared polymorphisms, loci with identical alleles across species, are of unique interest in evolutionary biology as they may represent cases of selection maintaining ancient genetic variation post-speciation, or contemporary selection promoting convergent evolution. In this study, we investigate the abundance of shared polymorphism between two members of the Daphnia pulex species complex. We test whether the presence of shared mutations is consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting or convergent evolution. We analyzed over 2,000 genomes from six taxa in the D. pulex species group and examined the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that North American and European D. pulex diverged over 10 million years ago, yet retained tens of thousands of shared polymorphisms. We suggest that the number of shared polymorphisms between North American and European D. pulex cannot be fully explained by hybridization or incomplete lineage sorting alone. We show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting convergent and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild.

Hybridization has localized effect on genetic variation in closely related pine species

BACKGROUND

Hybridization is a known phenomenon in nature but its genetic impact on populations of parental species remains less understood. We investigated the evolutionary consequences of the interspecific gene flow in several contact zones of closely related pine species. Using a set of genetic markers from both nuclear and organellar genomes, we analyzed four hybrid zones (384 individuals) and a large panel of reference allopatric populations of parental taxa (2104 individuals from 96 stands).

RESULTS

We observed reduced genetic diversity in maternally transmitted mitochondrial genomes of pure pine species and hybrids from contact zones compared to reference allopatric populations. The distribution of mtDNA haplotypes followed geographic rather than species boundaries. Additionally, no new haplotypes emerged in the contact zones, instead these zones contained the most common local variants. However, species diverged significantly at nuclear genomes and populations in contact zones exhibited similar or higher genetic diversity compared to the reference stands. There were no signs of admixture in any allopatric population, while clear admixture was evident in the contact zones, indicating that hybridization has a geographically localized effect on the genetic variation of the analyzed pine species.

CONCLUSIONS

Our results suggest that hybrid zones act as sinks rather than melting pots of genetic diversity. Hybridization influences sympatric populations but is confined to contact zones. The spectrum of parental species ancestry in hybrids reflects the old evolutionary history of the sympatric populations. These findings also imply that introgression may play a crucial role in the adaptation of hybrids to specific environments.

Balancing selection and the functional effects of shared polymorphism in cryptic Daphnia species

The patterns of genetic variation within and between related taxa represent the genetic history of a species. Shared polymorphisms, loci with identical alleles across species, are of unique interest as they may represent cases of ancient selection maintaining functional variation post-speciation. In this study, we investigate the abundance of shared polymorphism in the Daphnia pulex species complex. We test whether shared mutations are consistent with the action of balancing selection or alternative hypotheses such as hybridization, incomplete lineage sorting, or convergent evolution. We analyzed over 2,000 genomes from North American and European D. pulex and several outgroup species to examine the prevalence and distribution of shared alleles between the focal species pair, North American and European D. pulex. We show that while North American and European D. pulex diverged over ten million years ago, they retained tens of thousands of shared alleles. We found that the number of shared polymorphisms between North American and European D. pulex cannot be explained by hybridization or incomplete lineage sorting alone. Instead, we show that most shared polymorphisms could be the product of convergent evolution, that a limited number appear to be old trans-specific polymorphisms, and that balancing selection is affecting young and ancient mutations alike. Finally, we provide evidence that a blue wavelength opsin gene with trans-specific polymorphisms has functional effects on behavior and fitness in the wild. Ultimately, our findings provide insights into the genetic basis of adaptation and the maintenance of genetic diversity between species.

Botany and geogenomics: Constraining geological hypotheses in the neotropics with large-scale genetic data derived from plants

Decades of empirical research have revealed how the geological history of our planet shaped plant evolution by establishing well-known patterns (e.g., how mountain uplift resulted in high rates of diversification and replicate radiations in montane plant taxa). This follows a traditional approach where botanical data are interpreted in light of geological events. In this synthesis, I instead describe how by integrating natural history, phylogenetics, and population genetics, botanical research can be applied alongside geology and paleontology to inform our understanding of past geological and climatic processes. This conceptual shift aligns with the goals of the emerging field of geogenomics. In the neotropics, plant geogenomics is a powerful tool for the reciprocal exploration of two long standing questions in biology and geology: how the dynamic landscape of the region came to be and how it shaped the evolution of the richest flora. Current challenges that are specific to analytical approaches for plant geogenomics are discussed. I describe the scale at which various geological questions can be addressed from biological data and what makes some groups of plants excellent model systems for geogenomics research. Although plant geogenomics is discussed with reference to the neotropics, the recommendations given here for approaches to plant geogenomics can and should be expanded to exploring long-standing questions on how the earth evolved with the use of plant DNA.

Inferring the Demographic History and Inheritance Mode of Tetraploid Species Using ABC

Genomic patterns of diversity and divergence are impacted by certain life history traits, reproductive systems, and demographic history. The latter is characterized by fluctuations in population sizes over time, as well as by temporal patterns of introgression. For a given organism, identifying a demographic history that deviates from the standard neutral model allows a better understanding of its evolution but also helps to reduce the risk of false positives when screening for molecular targets of natural selection. Tetraploid organisms and beyond have demographic histories that are complicated by the mode of polyploidization, the mode of inheritance, and different scenarios of gene flow between sub-genomes and diploid parental species. Here we provide guidelines for experimenters wishing to address these issues through a flexible statistical framework: approximate Bayesian computation (ABC). The emphasis is on the general philosophy of the approach to encourage future users to exploit the enormous flexibility of ABC beyond the limitations imposed by generalist data analysis pipelines.

The roles of balancing selection and recombination in the evolution of rattlesnake venom

The origin of snake venom involved duplication and recruitment of non-venom genes into venom systems. Several studies have predicted that directional positive selection has governed this process. Venom composition varies substantially across snake species and venom phenotypes are locally adapted to prey, leading to coevolutionary interactions between predator and prey. Venom origins and contemporary snake venom evolution may therefore be driven by fundamentally different selection regimes, yet investigations of population-level patterns of selection have been limited. Here, we use whole-genome data from 68 rattlesnakes to test hypotheses about the factors that drive genomic diversity and differentiation in major venom gene regions. We show that selection has resulted in long-term maintenance of genetic diversity within and between species in multiple venom gene families. Our findings are inconsistent with a dominant role of directional positive selection and instead support a role of long-term balancing selection in shaping venom evolution. We also detect rapid decay of linkage disequilibrium due to high recombination rates in venom regions, suggesting that venom genes have reduced selective interference with nearby loci, including other venom paralogues. Our results provide an example of long-term balancing selection that drives trans-species polymorphism and help to explain how snake venom keeps pace with prey resistance.

Ecological Load and Balancing Selection in Circumboreal Barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography.

Spectral Diversification and Trans-Species Allelic Polymorphism during the Land-to-Sea Transition in Snakes

Snakes are descended from highly visual lizards [1] but have limited (probably dichromatic) color vision attributed to a dim-light lifestyle of early snakes [2-4]. The living species of front-fanged elapids, however, are ecologically very diverse, with ∼300 terrestrial species (cobras, taipans, etc.) and ∼60 fully marine sea snakes, plus eight independently marine, amphibious sea kraits [1]. Here, we investigate the evolution of spectral sensitivity in elapids by analyzing their opsin genes (which are responsible for sensitivity to UV and visible light), retinal photoreceptors, and ocular lenses. We found that sea snakes underwent rapid adaptive diversification of their visual pigments when compared with their terrestrial and amphibious relatives. The three opsins present in snakes (SWS1, LWS, and RH1) have evolved under positive selection in elapids, and in sea snakes they have undergone multiple shifts in spectral sensitivity toward the longer wavelengths that dominate below the sea surface. Several relatively distantly related Hydrophis sea snakes are polymorphic for shortwave sensitive visual pigment encoded by alleles of SWS1. This spectral site polymorphism is expected to confer expanded "UV-blue" spectral sensitivity and is estimated to have persisted twice as long as the predicted survival time for selectively neutral nuclear alleles. We suggest that this polymorphism is adaptively maintained across Hydrophis species via balancing selection, similarly to the LWS polymorphism that confers allelic trichromacy in some primates. Diving sea snakes thus appear to share parallel mechanisms of color vision diversification with fruit-eating primates.

Sex chromosome inversions enforce reproductive isolation across an avian hybrid zone

Across hybrid zones, the sex chromosomes are often more strongly differentiated than the autosomes. This is regularly attributed to the greater frequency of reproductive incompatibilities accumulating on sex chromosomes and their exposure in the heterogametic sex. Working within an avian hybrid zone, we explore the possibility that chromosome inversions differentially accumulate on the Z chromosome compared to the autosomes and thereby contribute to Z chromosome differentiation. We analyse the northern Australian hybrid zone between two subspecies of the long-tailed finch (Poephila acuticauda), first described based on differences in bill colour, using reduced-representation genomic sequencing for 293 individuals over a 1,530-km transect. Autosomal differentiation between subspecies is minimal. In contrast, 75% of the Z chromosome is highly differentiated and shows a steep genomic cline, which is displaced 350 km to the west of the cline in bill colour. Differentiation is associated with two or more putative chromosomal inversions, each predominating in one subspecies. If inversions reduce recombination between hybrid incompatibilities, they are selectively favoured and should therefore accumulate in hybrid zones. We argue that this predisposes inversions to differentially accumulate on the Z chromosome. One genomic region affecting bill colour is on the Z, but the main candidates are on chromosome 8. This and the displacement of the bill colour and Z chromosome cline centres suggest that bill colour has not strongly contributed to inversion accumulation. Based on cline width, however, the Z chromosome and bill colour both contribute to reproductive isolation established between this pair of subspecies.

Lima A, Nakamura AM, Fernandes F, Sobrinho I, de Brito RA. Evidence for Introgression Among Three Species of the Anastrepha fraterculus Group, a Radiating Species Complex of Fruit Flies

Introgression should no longer be considered as rare a phenomenon as once thought, since several studies have recently documented gene flow between closely related and radiating species. Here, we investigated evolutionary relationships among three closely related species of fruit flies of the Anastrepha fraterculus group (Anastrepha fraterculus, A. obliqua and A. sororcula). We sequenced a set of 20 genes and implemented a combined populational and phylogenetic inference with a model selection approach by an ABC framework in order to elucidate the demographic history of these species. The phylogenetic histories inferred from most genes showed a great deal of discordance and substantial shared polymorphic variation. The analysis of several population and speciation models reveal that this shared variation is better explained by introgression rather than convergence by parallel mutation or incomplete lineage sorting. Our results consistently showed these species evolving under an isolation with migration model experiencing a continuous and asymmetrical pattern of gene flow involving all species pairs, even though still showed a more closely related relationship between A. fraterculus and A. sororcula when compared with A. obliqua. This suggests that these species have been exchanging genes since they split from their common ancestor ∼2.6 MYA ago. We also found strong evidence for recent population expansion that appears to be consequence of anthropic activities affecting host crops of fruit flies. These findings point that the introgression here found may have been driven by genetic drift and not necessary by selection, which has implications for tracking and managing fruit flies.

Searching for ancient balanced polymorphisms shared between Neanderthals and Modern Humans

Hominin evolution is characterized by adaptive solutions often rooted in behavioral and cognitive changes. If balancing selection had an important and long-lasting impact on the evolution of these traits, it can be hypothesized that genes associated with them should carry an excess of shared polymorphisms (trans- SNPs) across recent Homo species. In this study, we investigate the role of balancing selection in human evolution using available exomes from modern (Homo sapiens) and archaic humans (H. neanderthalensis and Denisovan) for an excess of trans-SNP in two gene sets: one associated with the immune system (IMMS) and another one with behavioral system (BEHS). We identified a significant excess of trans-SNPs in IMMS (N=547), of which six of these located within genes previously associated with schizophrenia. No excess of trans-SNPs was found in BEHS, but five genes in this system harbor potential signals for balancing selection and are associated with psychiatric or neurodevelopmental disorders. Our approach evidenced recent Homo trans-SNPs that have been previously implicated in psychiatric diseases such as schizophrenia, suggesting that a genetic repertoire common to the immune and behavioral systems could have been maintained by balancing selection starting before the split between archaic and modern humans.

Analytical Biases Associated with GC-Content in Molecular Evolution

Molecular evolution is being revolutionized by high-throughput sequencing allowing an increased amount of genome-wide data available for multiple species. While base composition summarized by GC-content is one of the first metrics measured in genomes, its genomic distribution is a frequently neglected feature in downstream analyses based on DNA sequence comparisons. Here, we show how base composition heterogeneity among loci and taxa can bias common molecular evolution analyses such as phylogenetic tree reconstruction, detection of natural selection and estimation of codon usage. We then discuss the biological, technical and methodological causes of these GC-associated biases and suggest approaches to overcome them.

Mistaken Identity: Another Bias in the Use of Relative Genetic Divergence Measures for Detecting Interspecies Introgression

Measures of genetic divergence have long been used to identify evolutionary processes operating within and between species. However, recent reviews have described a bias in the use of relative divergence measures towards incorrectly identifying genomic regions that are seemingly immune to introgression. Here, we present a novel and opposite bias of relative divergence measures: misidentifying regions of introgression between sister species. We examine two distinct haplotypes of intermediate frequency within Drosophila pseudoobscura at the DPSX009 locus. One of these haplotypes had lower relative divergence than another to sister species D. persimilis. Although we and others initially presumed one haplotype have spread via introgression between D. pseudoobscura and D. persimilis, absolute divergence measures and individual sequence analysis suggest that haplotype structuring occurred as the result of within-species processes. The potential for this type of misinference may occur with any haplotype that recently spread within a species. We conclude that absolute measures of genetic divergence are necessary for confirming putative regions of introgression.

High regional genetic diversity and lack of host-specificity in Ostrinia nubilalis (Lepidoptera: Crambidae) as revealed by mtDNA variation

The European corn borer (Ostrinia nubilalis) infests a wide array of host plants and is considered one of the most serious pests of maize in Europe. Recent studies suggest that individuals feeding on maize in Europe should be referred to O. nubilalis (sensu nov.), while those infesting dicots as Ostrinia scapulalis (sensu nov.). We test if the clear genetic distinctiveness among individuals of O. nubilalis living on maize vs. dicots is tracked by mitochondrial DNA (mtDNA). We used fragments of COI and COII genes of 32 individuals traditionally recognized as O. nubilalis collected on three host plants, maize, mugwort and hop, growing in different parts of Poland. In addition, we reconstructed the mtDNA phylogeny of Ostrinia species based on our data and sequences retrieved from GenBank to assess host and/or biogeographic patterns. We also compared haplotype variation found in Poland (east-central Europe) with other regions (Anatolia, Eastern Europe, Balkans, Far East, North America). Our study showed high mtDNA diversity of O. nubilalis in Poland in comparison with other regions and revealed rare haplotypes likely of Asian origin. We did not find distinct mtDNA haplotypes in larvae feeding on maize vs. dicotyledonous plants. Phylogenetic analyses showed an apparent lack of mtDNA divergence among putatively distinct lineages belonging to the O. nubilalis group as identical haplotypes are shared by Asian and European individuals. We argue that human-mediated dispersal, hybridization and sporadic host jumps are likely responsible for the lack of a geographic pattern in mtDNA variation.

Purifying selection shapes the coincident SNP distribution of primate coding sequences

Genome-wide analysis has observed an excess of coincident single nucleotide polymorphisms (coSNPs) at human-chimpanzee orthologous positions, and suggested that this is due to cryptic variation in the mutation rate. While this phenomenon primarily corresponds with non-coding coSNPs, the situation in coding sequences remains unclear. Here we calculate the observed-to-expected ratio of coSNPs (coSNP_O/E) to estimate the prevalence of human-chimpanzee coSNPs, and show that the excess of coSNPs is also present in coding regions. Intriguingly, coSNP_O/E is much higher at zero-fold than at nonzero-fold degenerate sites; such a difference is due to an elevation of coSNP_O/E at zero-fold degenerate sites, rather than a reduction at nonzero-fold degenerate ones. These trends are independent of chimpanzee subpopulation, population size, or sequencing techniques; and hold in broad generality across primates. We find that this discrepancy cannot fully explained by sequence contexts, shared ancestral polymorphisms, SNP density, and recombination rate, and that coSNP_O/E in coding sequences is significantly influenced by purifying selection. We also show that selection and mutation rate affect coSNP_O/E independently, and coSNPs tend to be less damaging and more correlated with human diseases than non-coSNPs. These suggest that coSNPs may represent a “signature” during primate protein evolution.

A genomic perspective on hybridization and speciation

Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.

Selection for population-specific adaptation shaped patterns of variation in the photoperiod pathway genes in Arabidopsis lyrata during post-glacial colonization

Spatially varying selection can lead to population-specific adaptation, which is often recognized at the phenotypic level; however, the genetic evidence is weaker in many groups of organisms. In plants, environmental shifts that occur due to colonization of a novel environment may require adaptive changes in the timing of growth and flowering, which are often governed by location-specific environmental cues such as day length. We studied locally varying selection in 19 flowering time loci in nine populations of the perennial herb Arabidopsis lyrata, which has a wide but patchy distribution in temperate and boreal regions of the northern hemisphere. The populations differ in their recent population demographic and colonization histories and current environmental conditions, especially in the growing season length. We searched for population-specific molecular signatures of directional selection by comparing a set of candidate flowering time loci with a genomic reference set within each population using multiple approaches and contrasted the patterns of different populations. The candidate loci possessed approximately 20% of the diversity of the reference loci. On average the flowering time loci had more rare alleles (a smaller Tajima's D) and an excess of highly differentiated sites relative to the reference, suggesting positive selection. The strongest signal of selection was detected in photoperiodic pathway loci in the colonizing populations of Northwestern Europe, whereas no evidence of positive selection was detected in the Central European populations. These findings emphasized the population-specific nature of selection and suggested that photoperiodic adaptation was important during postglacial colonization of the species.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

Identification of Allorecognition Loci in Neurospora crassa by Genomics and Evolutionary Approaches

Understanding the genetic and molecular bases of the ability to distinguish self from nonself (allorecognition) and mechanisms underlying evolution of allorecognition systems is an important endeavor for understanding cases where it becomes dysfunctional, such as in autoimmune disorders. In filamentous fungi, allorecognition can result in vegetative or heterokaryon incompatibility, which is a type of programmed cell death that occurs following fusion of genetically different cells. Allorecognition is genetically controlled by het loci, with coexpression of any combination of incompatible alleles triggering vegetative incompatibility. Herein, we identified, characterized, and inferred the evolutionary history of candidate het loci in the filamentous fungus Neurospora crassa. As characterized het loci encode proteins carrying an HET domain, we annotated HET domain genes in 25 isolates from a natural population along with the N. crassa reference genome using resequencing data. Because allorecognition systems can be affected by frequency-dependent selection favoring rare alleles (i.e., balancing selection), we mined resequencing data for HET domain loci whose alleles displayed elevated levels of variability, excess of intermediate frequency alleles, and deep gene genealogies. From these analyses, 34 HET domain loci were identified as likely to be under balancing selection. Using transformation, incompatibility assays and genetic analyses, we determined that one of these candidates functioned as a het locus (het-e). The het-e locus has three divergent allelic groups that showed signatures of positive selection, intra- and intergroup recombination, and trans-species polymorphism. Our findings represent a compelling case of balancing selection functioning on multiple alleles across multiple loci potentially involved in allorecognition.

Trans-species polymorphism at antimicrobial innate immunity cathelicidin genes of Atlantic cod and related species

Natural selection, the most important force in evolution, comes in three forms. Negative purifying selection removes deleterious variation and maintains adaptations. Positive directional selection fixes beneficial variants, producing new adaptations. Balancing selection maintains variation in a population. Important mechanisms of balancing selection include heterozygote advantage, frequency-dependent advantage of rarity, and local and fluctuating episodic selection. A rare pathogen gains an advantage because host defenses are predominantly effective against prevalent types. Similarly, a rare immune variant gives its host an advantage because the prevalent pathogens cannot escape the host's apostatic defense. Due to the stochastic nature of evolution, neutral variation may accumulate on genealogical branches, but trans-species polymorphisms are rare under neutrality and are strong evidence for balancing selection. Balanced polymorphism maintains diversity at the major histocompatibility complex (MHC) in vertebrates. The Atlantic cod is missing genes for both MHC-II and CD4, vital parts of the adaptive immune system. Nevertheless, cod are healthy in their ecological niche, maintaining large populations that support major commercial fisheries. Innate immunity is of interest from an evolutionary perspective, particularly in taxa lacking adaptive immunity. Here, we analyze extensive amino acid and nucleotide polymorphisms of the cathelicidin gene family in Atlantic cod and closely related taxa. There are three major clusters, Cath1, Cath2, and Cath3, that we consider to be paralogous genes. There is extensive nucleotide and amino acid allelic variation between and within clusters. The major feature of the results is that the variation clusters by alleles and not by species in phylogenetic trees and discriminant analysis of principal components. Variation within the three groups shows trans-species polymorphism that is older than speciation and that is suggestive of balancing selection maintaining the variation. Using Bayesian and likelihood methods positive and negative selection is evident at sites in the conserved part of the genes and, to a larger extent, in the active part which also shows episodic diversifying selection, further supporting the argument for balancing selection.

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.

Footprints of ancient-balanced polymorphisms in genetic variation data from closely related species

When long-lasting, balancing selection can lead to “trans-species” polymorphisms that are shared by two or more species identical by descent. In such cases, the gene genealogy at the selected site clusters by allele instead of by species, and nearby neutral sites also have unusual genealogies because of linkage. While this scenario is expected to leave discernible footprints in genetic variation data, the specific patterns remain poorly characterized. Motivated by recent findings in primates, we focus on the case of a biallelic polymorphism under ancient balancing selection and derive approximations for summaries of the polymorphism data from two species. Specifically, we characterize the length of the segment that carries most of the footprints, the expected number of shared neutral single nucleotide polymorphisms (SNPs), and the patterns of allelic associations among them. We confirm the accuracy of our approximations by coalescent simulations. We further show that for humans and chimpanzees—more generally, for pairs of species with low genetic diversity levels—these patterns are highly unlikely to be generated by neutral recurrent mutations. We discuss the implications for the design and interpretation of genome scans for ancient balanced polymorphisms in primates and other taxa.

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.

Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds

Animals that sustain high levels of aerobic activity under hypoxic conditions (e.g., birds that fly at high altitude) face the physiological challenge of jointly optimizing blood-O2 affinity for O2 loading in the pulmonary circulation and O2 unloading in the systemic circulation. At high altitude, this challenge is especially acute for small endotherms like hummingbirds that have exceedingly high mass-specific metabolic rates. Here we report an experimental analysis of hemoglobin (Hb) function in South American hummingbirds that revealed a positive correlation between Hb-O2 affinity and native elevation. Protein engineering experiments and ancestral-state reconstructions revealed that this correlation is attributable to derived increases in Hb-O2 affinity in highland lineages, as well as derived reductions in Hb-O2 affinity in lowland lineages. Site-directed mutagenesis experiments demonstrated that repeated evolutionary transitions in biochemical phenotype are mainly attributable to repeated amino acid replacements at two epistatically interacting sites that alter the allosteric regulation of Hb-O2 affinity. These results demonstrate that repeated changes in biochemical phenotype involve parallelism at the molecular level, and that mutations with indirect, second-order effects on Hb allostery play key roles in biochemical adaptation.

Linking great apes genome evolution across time scales using polymorphism-aware phylogenetic models

The genomes of related species contain valuable information on the history of the considered taxa. Great apes in particular exhibit variation of evolutionary patterns along their genomes. However, the great ape data also bring new challenges, such as the presence of incomplete lineage sorting and ancestral shared polymorphisms. Previous methods for genome-scale analysis are restricted to very few individuals or cannot disentangle the contribution of mutation rates and fixation biases. This represents a limitation both for the understanding of these forces as well as for the detection of regions affected by selection. Here, we present a new model designed to estimate mutation rates and fixation biases from genetic variation within and between species. We relax the assumption of instantaneous substitutions, modeling substitutions as mutational events followed by a gradual fixation. Hence, we straightforwardly account for shared ancestral polymorphisms and incomplete lineage sorting. We analyze genome-wide synonymous site alignments of human, chimpanzee, and two orangutan species. From each taxon, we include data from several individuals. We estimate mutation rates and GC-biased gene conversion intensity. We find that both mutation rates and biased gene conversion vary with GC content. We also find lineage-specific differences, with weaker fixation biases in orangutan species, suggesting a reduced historical effective population size. Finally, our results are consistent with directional selection acting on coding sequences in relation to exonic splicing enhancers.

Crossing the species barrier: genomic hotspots of introgression between two highly divergent Ciona intestinalis species

Inferring a realistic demographic model from genetic data is an important challenge to gain insights into the historical events during the speciation process and to detect molecular signatures of selection along genomes. Recent advances in divergence population genetics have reported that speciation in face of gene flow occurred more frequently than theoretically expected, but the approaches used did not account for genome-wide heterogeneity (GWH) in introgression rates. Here, we investigate the impact of GWH on the inference of divergence with gene flow between two cryptic species of the marine model Ciona intestinalis by analyzing polymorphism and divergence patterns in 852 protein-coding sequence loci. These morphologically similar entities are highly diverged molecular-wise, but evidence of hybridization has been reported in both laboratory and field studies. We compare various speciation models and test for GWH under the approximate Bayesian computation framework. Our results demonstrate the presence of significant extents of gene flow resulting from a recent secondary contact after >3 My of divergence in isolation. The inferred rates of introgression are relatively low, highly variable across loci and mostly unidirectional, which is consistent with the idea that numerous genetic incompatibilities have accumulated over time throughout the genomes of these highly diverged species. A genomic map of the level of gene flow identified two hotspots of introgression, that is, large genome regions of unidirectional introgression. This study clarifies the history and degree of isolation of two cryptic and partially sympatric model species and provides a methodological framework to investigate GWH at various stages of speciation process.

Molecular hyperdiversity and evolution in very large populations

The genomic density of sequence polymorphisms critically affects the sensitivity of inferences about ongoing sequence evolution, function and demographic history. Most animal and plant genomes have relatively low densities of polymorphisms, but some species are hyperdiverse with neutral nucleotide heterozygosity exceeding 5%. Eukaryotes with extremely large populations, mimicking bacterial and viral populations, present novel opportunities for studying molecular evolution in sexually reproducing taxa with complex development. In particular, hyperdiverse species can help answer controversial questions about the evolution of genome complexity, the limits of natural selection, modes of adaptation and subtleties of the mutation process. However, such systems have some inherent complications and here we identify topics in need of theoretical developments. Close relatives of the model organisms Caenorhabditis elegans and Drosophila melanogaster provide known examples of hyperdiverse eukaryotes, encouraging functional dissection of resulting molecular evolutionary patterns. We recommend how best to exploit hyperdiverse populations for analysis, for example, in quantifying the impact of noncrossover recombination in genomes and for determining the identity and micro-evolutionary selective pressures on noncoding regulatory elements.

Unconstrained evolution in short introns? - an analysis of genome-wide polymorphism and divergence data from Drosophila

An unconstrained reference sequence facilitates the detection of selection. In Drosophila, sequence variation in short introns seems to be least influenced by selection and dominated by mutation and drift. Here, we test this with genome-wide sequences using an African population (Malawi) of D. melanogaster and data from the related outgroup species D. simulans, D. sechellia, D. erecta and D. yakuba. The distribution of mutations deviates from equilibrium, and the content of A and T (AT) nucleotides shows an excess of variance among introns. We explain this by a complex mutational pattern: a shift in mutational bias towards AT, leading to a slight nonequilibrium in base composition and context-dependent mutation rates, with G or C (GC) sites mutating most frequently in AT-rich introns. By comparing the corresponding allele frequency spectra of AT-rich vs. GC-rich introns, we can rule out the influence of directional selection or biased gene conversion on the mutational pattern. Compared with neutral equilibrium expectations, polymorphism spectra show an excess of low frequency and a paucity of intermediate frequency variants, irrespective of the direction of mutation. Combining the information from different outgroups with the polymorphism data and using a generalized linear model, we find evidence for shared ancestral polymorphism between D. melanogaster and D. simulans, D. sechellia, arguing against a bottleneck in D. melanogaster. Generally, we find that short introns can be used as a neutral reference on a genome-wide level, if the spatially and temporally varying mutational pattern is accounted for.

Widespread horizontal genomic exchange does not erode species barriers among sympatric ducks

BACKGROUND

The study of speciation and maintenance of species barriers is at the core of evolutionary biology. During speciation the genome of one population becomes separated from other populations of the same species, which may lead to genomic incompatibility with time. This separation is complete when no fertile offspring is produced from inter-population matings, which is the basis of the biological species concept. Birds, in particular ducks, are recognised as a challenging and illustrative group of higher vertebrates for speciation studies. There are many sympatric and ecologically similar duck species, among which fertile hybrids occur relatively frequently in nature, yet these species remain distinct.

RESULTS

We show that the degree of shared single nucleotide polymorphisms (SNPs) between five species of dabbling ducks (genus Anas) is an order of magnitude higher than that previously reported between any pair of eukaryotic species with comparable evolutionary distances. We demonstrate that hybridisation has led to sustained exchange of genetic material between duck species on an evolutionary time scale without disintegrating species boundaries. Even though behavioural, genetic and ecological factors uphold species boundaries in ducks, we detect opposing forces allowing for viable interspecific hybrids, with long-term evolutionary implications. Based on the superspecies concept we here introduce the novel term "supra-population" to explain the persistence of SNPs identical by descent within the studied ducks despite their history as distinct species dating back millions of years.

CONCLUSIONS

By reviewing evidence from speciation theory, palaeogeography and palaeontology we propose a fundamentally new model of speciation to accommodate our genetic findings in dabbling ducks. This model, we argue, may also shed light on longstanding unresolved general speciation and hybridisation patterns in higher organisms, e.g. in other bird groups with unusually high hybridisation rates. Observed parallels to horizontal gene transfer in bacteria facilitate the understanding of why ducks have been such an evolutionarily successful group of animals. There is large evolutionary potential in the ability to exchange genes among species and the resulting dramatic increase of effective population size to counter selective constraints.

A trans-specific polymorphism in ZC3HAV1 is maintained by long-standing balancing selection and may confer susceptibility to multiple sclerosis

The human ZC3HAV1 gene encodes an antiviral protein. The longest splicing isoform of ZC3HAV1 contains a C-terminal PARP-like domain, which has evolved under positive selection in primates. We analyzed the evolutionary history of this same domain in humans and in Pan troglodytes. We identified two variants that segregate in both humans and chimpanzees; one of them (rs3735007) does not occur at a hypermutable site and accounts for a nonsynonymous substitution (Thr851Ile). The probability that the two trans-specific polymorphisms have occurred independently in the two lineages was estimated to be low (P = 0.0054), suggesting that at least one of them has arisen before speciation and has been maintained by selection. Population genetic analyses in humans indicated that the region surrounding the shared variants displays strong evidences of long-standing balancing selection. Selection signatures were also observed in a chimpanzee population sample. Inspection of 1000 Genomes data confirmed these findings but indicated that search for selection signatures using low-coverage whole-genome data may need masking of repetitive sequences. A case–control study of more than 1,000 individuals from mainland Italy indicated that the Thr851Ile SNP is significantly associated with susceptibility to multiple sclerosis (MS) (odds ratio [OR] = 1.47, 95% confidence intervals [CI]: 1.08–1.99, P = 0.011). This finding was confirmed in a larger sample of 4,416 Sardinians cases/controls (OR = 1.18, 95% CI: 1.037–1.344, P = 0.011), but not in a population from Belgium. We provide one of the first instances of human/chimpanzee trans-specific coding variant located outside the major histocompatibility complex region. The selective pressure is likely to be virus driven; in modern populations, this variant associates with susceptibility to MS, possibly via the interaction with environmental factors.

Nuclear and plastid haplotypes suggest rapid diploid and polyploid speciation in the N Hemisphere Achillea millefolium complex (Asteraceae)

Background

Species complexes or aggregates consist of a set of closely related species often of different ploidy levels, whose relationships are difficult to reconstruct. The N Hemisphere Achillea millefolium aggregate exhibits complex morphological and genetic variation and a broad ecological amplitude. To understand its evolutionary history, we study sequence variation at two nuclear genes and three plastid loci across the natural distribution of this species complex and compare the patterns of such variations to the species tree inferred earlier from AFLP data.

Results

Among the diploid species of A. millefolium agg., gene trees of the two nuclear loci, ncpGS and SBP, and the combined plastid fragments are incongruent with each other and with the AFLP tree likely due to incomplete lineage sorting or secondary introgression. In spite of the large distributional range, no isolation by distance is found. Furthermore, there is evidence for intragenic recombination in the ncpGS gene. An analysis using a probabilistic model for population demographic history indicates large ancestral effective population sizes and short intervals between speciation events. Such a scenario explains the incongruence of the gene trees and species tree we observe. The relationships are particularly complex in the polyploid members of A. millefolium agg.

Conclusions

The present study indicates that the diploid members of A. millefolium agg. share a large part of their molecular genetic variation. The findings of little lineage sorting and lack of isolation by distance is likely due to short intervals between speciation events and close proximity of ancestral populations. While previous AFLP data provide species trees congruent with earlier morphological classification and phylogeographic considerations, the present sequence data are not suited to recover the relationships of diploid species in A. millefolium agg. For the polyploid taxa many hybrid links and introgression from the diploids are suggested.

Drosophila duplication hotspots are associated with late-replicating regions of the genome

Duplications play a significant role in both extremes of the phenotypic spectrum of newly arising mutations: they can have severe deleterious effects (e.g. duplications underlie a variety of diseases) but can also be highly advantageous. The phenotypic potential of newly arisen duplications has stimulated wide interest in both the mutational and selective processes shaping these variants in the genome. Here we take advantage of the Drosophila simulans–Drosophila melanogaster genetic system to further our understanding of both processes. Regarding mutational processes, the study of two closely related species allows investigation of the potential existence of shared duplication hotspots, and the similarities and differences between the two genomes can be used to dissect its underlying causes. Regarding selection, the difference in the effective population size between the two species can be leveraged to ask questions about the strength of selection acting on different classes of duplications. In this study, we conducted a survey of duplication polymorphisms in 14 different lines of D. simulans using tiling microarrays and combined it with an analogous survey for the D. melanogaster genome. By integrating the two datasets, we identified duplication hotspots conserved between the two species. However, unlike the duplication hotspots identified in mammalian genomes, Drosophila duplication hotspots are not associated with sequences of high sequence identity capable of mediating non-allelic homologous recombination. Instead, Drosophila duplication hotspots are associated with late-replicating regions of the genome, suggesting a link between DNA replication and duplication rates. We also found evidence supporting a higher effectiveness of selection on duplications in D. simulans than in D. melanogaster. This is also true for duplications segregating at high frequency, where we find evidence in D. simulans that a sizeable fraction of these mutations is being driven to fixation by positive selection.

DNA duplications are important contributors to the phenotypic differences observed between individuals. These mutations can disrupt the normal functioning of genes and so are often associated with disease. But because they can add genetic information they can also lead to evolutionary change. Understanding how selection and non-random mutation processes shape the distribution of duplications throughout the genome is important to elucidate both the medical and evolutionary impacts of these mutations. Here, we examined the roles of selection and mutation in shaping patterns of duplication polymorphisms across the genomes of the fruit fly Drosophila melanogaster and its sister species, D. simulans. We found that selection is pervasive in both genomes but is more efficient in D. simulans than in D. melanogaster. We also found that these two species have shared duplication hotspots, i.e. orthologous regions experiencing high rates of duplication in the two genomes. After excluding the hypothesis that Drosophila duplication hotspots are associated with regions of the genome rich in segmental duplications (as observed for mammalian genomes), we show that they are associated with late-replicating regions of the genome. Our work therefore proposes a link between DNA replication and rates of duplication across the genome.

Ancestral polymorphisms in Drosophila pseudoobscura and Drosophila miranda

Ancestral polymorphisms are defined as variants that arose by mutation prior to the speciation event that generated the species in which they segregate. Their presence may complicate the interpretation of molecular data and lead to incorrect phylogenetic inferences. They may also be used to identify regions of the genome that are under balancing selection. It is thus important to take into account the contribution of ancestral polymorphisms to variability within species and divergence between species. Here, we extend and improve a method for estimation of the proportion of ancestral polymorphisms within a species, and apply it to a dataset of 33 X-linked and 34 autosomal protein-coding genes for which sequence polymorphism data are available in both Drosophila pseudoobscura and Drosophila miranda, using Drosophila affinis as an outgroup. We show that a substantial proportion of both X-linked and autosomal synonymous variants in these two species are ancestral, and that a small number of additional genes with unusually high sequence diversity seem to have an excess of ancestral polymorphisms, suggestive of balancing selection.

Characterizing the influence of effective population size on the rate of adaptation: Gillespie's Darwin domain

Characterizing the role of effective population size in dictating the rate of adaptive evolution remains a major challenge in evolutionary biology. Depending on the underlying distribution of fitness effects of new mutations, populations of different sizes may differ vastly in their rate of adaptation. Here, we collect polymorphism data at over 100 loci for two closely related Drosophila species with different current effective population sizes (N_e), Drosophila miranda and D. pseudoobscura, to evaluate the prevalence of adaptive evolution versus genetic drift in molecular evolution. Utilizing these large and consistently sampled data sets, we obtain greatly improved estimates of the demographic histories of both species. Specifically, although current N_e differs between these species, their ancestral sizes were much more similar. We find that statistical approaches capturing recent adaptive evolution (using patterns of polymorphisms) detect higher rates of adaptive evolution in the larger D. pseudoobscura population. In contrast, methods aimed at detecting selection over longer time periods (i.e., those relying on divergence data) estimate more similar rates of adaptation between the two species. Thus, our results suggest an important role of effective population size in dictating rates of adaptation and highlight how complicated population histories—as is probably the case for most species—can effect rates of adaptation. Additionally, we also show how different methodologies to detect positive selection can reveal information about different timescales of adaptive evolution.

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.

Mutation rate distribution inferred from coincident SNPs and coincident substitutions

Mutation rate variation has the potential to bias evolutionary inference, particularly when rates become much higher than the mean. We first confirm prior work that inferred the existence of cryptic, site-specific rate variation on the basis of coincident polymorphisms—sites that are segregating in both humans and chimpanzees. Then we extend this observation to a longer evolutionary timescale by identifying sites of coincident substitutions using four species. From these data, we develop analytic theory to infer the variance and skewness of the distribution of mutation rates. Even excluding CpG dinucleotides, we find a relatively large coefficient of variation and positive skew, which suggests that, although most sites in the genome have mutation rates near the mean, the distribution contains a long right-hand tail with a small number of sites having high mutation rates. At least for primates, these quickly mutating sites are few enough that the infinite sites model in population genetics remains appropriate.

Effective population size and the efficacy of selection on the X chromosomes of two closely related Drosophila species

The prevalence of natural selection relative to genetic drift is of central interest in evolutionary biology. Depending on the distribution of fitness effects of new mutations, the importance of these evolutionary forces may differ in species with different effective population sizes. Here, we survey population genetic variation at 105 orthologous X-linked protein coding regions in Drosophila melanogaster and its sister species D. simulans, two closely related species with distinct demographic histories. We observe significantly higher levels of polymorphism and evidence for stronger selection on codon usage bias in D. simulans, consistent with a larger historical effective population size on average for this species. Despite these differences, we estimate that <10% of newly arising nonsynonymous mutations have deleterious fitness effects in the nearly neutral range (i.e., −10 < N_es < 0) in both species. The inferred distributions of fitness effects and demographic models translate into surprisingly high estimates of the fraction of “adaptive” protein divergence in both species (∼85–90%). Despite evidence for different demographic histories, differences in population size have apparently played little role in the dynamics of protein evolution in these two species, and estimates of the adaptive fraction (α) of protein divergence in both species remain high even if we account for recent 10-fold growth. Furthermore, although several recent studies have noted strong signatures of recurrent adaptive protein evolution at genes involved in immunity, reproduction, sexual conflict, and intragenomic conflict, our finding of high levels of adaptive protein divergence at randomly chosen proteins (with respect to function) suggests that many other factors likely contribute to the adaptive protein divergence signature in Drosophila.

Evolution of disease response genes in loblolly pine: insights from candidate genes

Background

Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen.

Methodology/Principal Findings

Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models.

Conclusions/Significance

Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States.

Host-specific differentiation among populations of Venturia inaequalis causing scab on apple, pyracantha and loquat

Patterns of multilocus DNA sequence variation within and between closely related taxa can provide insights into the history of divergence. Here, we report on DNA polymorphism and divergence at six nuclear loci in globally distributed samples of the ascomycete Venturia inaequalis, responsible for scab on apple, loquat, and pyracantha. Isolates from different hosts were differentiated but did not form diagnosable distinct phylogenetic species. Parameters of an Isolation-with-Migration model estimated from the data suggested that the large amount of variation shared among groups more likely resulted from recent splitting than from extensive genetic exchanges. Inferred levels of gene flow among groups were low and more concentrated toward recent times, and we identified two potentially recent one-off shifters from apple and pyracantha to loquat. These findings support a scenario of recent divergence in allopatry followed by introgression through secondary contact, with groups from loquat and pyracantha being the most recently differentiated.

Targets of balancing selection in the human genome

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Historical divergence and gene flow in the genus Zea

Gene flow plays a fundamental role in plant evolutionary history, yet its role in population divergence--and ultimately speciation--remains poorly understood. We investigated gene flow and the modalities of divergence in the domesticated Zea mays ssp. mays and three wild Zea taxa using sequence polymorphism data from 26 nuclear loci. We described diversity across loci and assessed evidence for adaptive and purifying selection at nonsynonymous sites. For each of three divergence events in the history of these taxa, we used approximate Bayesian simulation to estimate population sizes and divergence times and explicitly compare among alternative models of divergence. Our estimates of divergence times are surprisingly consistent with previous data from other markers and suggest rapid diversification of lineages within Zea in the last approximately 150,000 years. We found widespread evidence of historical gene flow, including evidence for divergence in the face of gene flow. We speculate that cultivated maize may serve as a bridge for gene flow among otherwise allopatric wild taxa.

The signature of long-standing balancing selection at the human defensin beta-1 promoter

BACKGROUND

Defensins, small endogenous peptides with antimicrobial activity, are pivotal components of the innate immune response. A large cluster of defensin genes is located on human chromosome 8p; among them the beta defensin 1 (DEFB1) promoterhas been extensively studied since discovery that specific polymorphisms and haplotypes associate with asthma and atopy, susceptibility to severe sepsis, as well as HIV and Candida infection predisposition.

RESULTS

Here, we characterize the sequence variation and haplotype structure of the DEFB1 promoter region in six human populations. In all of them, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, reduced population differentiation and a genealogy with common haplotypes separated by deep branches. Indeed, a significant departure from the expectation of evolutionary neutrality was observed in all populations and the possibility that this is due to demographic history alone was ruled out. Also, we verified that the selection signature is restricted to the promoter region and not due to a linked balanced polymorphism. A phylogeny-based estimation indicated that the two major haplotype clades separated around 4.5 million years ago, approximately the time when the human and chimpanzee lineages split.

CONCLUSION

Altogether, these features represent strong molecular signatures of long-term balancing selection, a process that is thought to be extremely rare outside major histocompatibility complex genes. Our data indicate that the DEFB1 promoter region carries functional variants and support previous hypotheses whereby alleles predisposing to atopic disorders are widespread in modern societies because they conferred resistance to pathogens in ancient settings.

Localization of candidate regions maintaining a common polymorphic inversion (2La) in Anopheles gambiae

Chromosomal inversion polymorphisms are thought to play a role in adaptive divergence, but the genes conferring adaptive benefits remain elusive. Here we study 2La, a common polymorphic inversion in the African malaria vector Anopheles gambiae. The frequency of 2La varies clinally and seasonally in a pattern suggesting response to selection for aridity tolerance. By hybridizing genomic DNA from individual mosquitoes to oligonucleotide microarrays, we obtained a complete map of differentiation across the A. gambiae genome. Comparing mosquitoes homozygous for the 2La gene arrangement or its alternative (2L+^a), divergence was highest at loci within the rearranged region. In the 22 Mb included within alternative arrangements, two ∼1.5 Mb regions near but not adjacent to the breakpoints were identified as being significantly diverged, a conclusion validated by targeted sequencing. The persistent association of both regions with the 2La arrangement is highly unlikely given known recombination rates across the inversion in 2La heterozygotes, thus implicating selection on genes underlying these regions as factors responsible for the maintenance of 2La. Polymorphism and divergence data are consistent with a model in which the inversion is maintained by migration-selection balance between multiple alleles inside these regions, but further experiments will be needed to fully distinguish between the epistasis (coadaptation) and local adaptation models for the maintenance of 2La.

A chromosomal inversion occurs when part of the chromosome breaks, rotates 180 degrees, and rejoins the broken chromosome. The result is a chromosome carrying a segment whose gene order is reversed. Whereas the physical rearrangement itself may have no direct consequences on gene function, recombination between alleles in the rearranged and wild type segments is suppressed. If multiple alleles inside the inverted or original orientations are well adapted to contrasting environmental conditions, suppressed recombination provides a mechanism to keep beneficial allelic combinations from being shuffled between different genetic backgrounds. Working with wild populations of flies, Dobzhansky provided the first evidence that selection was key to maintaining inversion polymorphism. Subsequently, examples of inversion polymorphisms under selection in other organisms have been found, notably in the mosquito that transmits most cases of human malaria, Anopheles gambiae. However, the genes or gene regions conferring fitness advantages have yet to be discovered. In this study, the authors used modern genomics tools to map such regions in an inversion at an unprecedented level of detail, and show that these regions are likely to be responsible for the maintenance of the inversion polymorphism in natural populations. This study lays the groundwork for future efforts to identify the genes themselves and their role in adaptation.

Divergence with gene flow between Ponto-Caspian refugia in an anadromous cyprinid Rutilus frisii revealed by multiple gene phylogeography

The Black and Caspian Seas have experienced alternating periods of isolation and interconnection over many Milankovitch climate oscillations and most recently became separated when the meltwater overflow from the Caspian Sea ceased at the end of the last glaciation. Climate-induced habitat changes have indisputably had profound impacts on distribution and demography of aquatic species, yet uncertainties remain about the relative roles of isolation and dispersal in the response of species shared between the Black and Caspian Sea basins. We examined these issues using phylogeographical analysis of an anadromous cyprinid fish Rutilus frisii. Bayesian coalescence analyses of sequence variation at two nuclear and one mitochondrial genes suggest that the Black and Caspian Seas supported separate populations of R. frisii during the last glaciation. Parameter estimates from the fitted isolation-with-migration model showed that their separation was not complete, however, and that the two populations continued to exchange genes in both directions. These analyses also suggested that majority of migrations occurred during the Pleistocene, showing that the variation shared between the Black and Caspian Seas is the result of ancient dispersal along the temporary natural connections between the basins, rather than of incomplete lineage sorting or recent human-mediated dispersal. Gene flow between the refugial populations was therefore an important source of genetic variation, and we suggest that it facilitated the evolutionary response of the populations to changing climate.