Selection versus random drift: long-term polymorphism persistence in small populations (evidence and modelling)

A Numerical Model Supports the Evolutionary Advantage of Recombination Plasticity in Shifting Environments

AbstractNumerous empirical studies have witnessed an increase in meiotic recombination rate in response to physiological stress imposed by unfavorable environmental conditions. Thus, inherited plasticity in recombination rate is hypothesized to be evolutionarily advantageous in changing environments. Previous theoretical models proceeded from the assumption that organisms increase their recombination rate when the environment becomes more stressful and demonstrated the evolutionary advantage of such a form of plasticity. Here, we numerically explore a complementary scenario-when the plastic increase in recombination rate is triggered by the environmental shifts. Specifically, we assume increased recombination in individuals developing in a different environment than their parents and, optionally, also in offspring of such individuals. We show that such shift-inducible recombination is always superior when the optimal constant recombination implies an intermediate rate. Moreover, under certain conditions, plastic recombination may also appear beneficial when the optimal constant recombination is either zero or free. The advantage of plastic recombination was better predicted by the range of the population's mean fitness over the period of environmental fluctuations, compared with the geometric mean fitness. These results hold for both panmixia and partial selfing, with faster dynamics of recombination modifier alleles under selfing. We think that recombination plasticity can be acquired under the control of environmentally responsive mechanisms, such as chromatin epigenetics remodeling.

Heterozygote advantage and pleiotropy contribute to intraspecific color trait variability

The persistence of intrapopulation phenotypic variation typically requires some form of balancing selection because drift and directional selection eventually erode genetic variation. Heterozygote advantage remains a classic explanation for the maintenance of genetic variation in the face of selection. However, examples of heterozygote advantage, other than those associated with disease resistance, are rather uncommon. Across most of its distribution, males of the aposematic moth Arctia plantaginis have two hindwing phenotypes determined by a heritable one locus-two allele polymorphism (genotypes: WW/Wy = white morph, yy = yellow morph). Using genotyped moths, we show that the presence of one or two copies of the yellow allele affects several life-history traits. Reproductive output of both males and females and female mating success are negatively affected by two copies of the yellow allele. Females carrying one yellow allele (i.e., Wy) have higher fertility, hatching success, and offspring survival than either homozygote, thus leading to strong heterozygote advantage. Our results indicate strong female contribution especially at the postcopulatory stage in maintaining the color polymorphism. The interplay between heterozygote advantage, yellow allele pleiotropic effect, and morph-specific predation pressure may exert balancing selection on the color locus, suggesting that color polymorphism may be maintained through complex interactions between natural and sexual selection.

Spatial genetic diversity in the Cape mole-rat, Georychus capensis: Extreme isolation of populations in a subterranean environment

The subterranean niche harbours animals with extreme adaptations. These adaptations decrease the vagility of taxa and, along with other behavioural adaptations, often result in isolated populations characterized by small effective population sizes, high inbreeding, population bottlenecks, genetic drift and consequently, high spatial genetic structure. Although information is available for some species, estimates of genetic diversity and whether this variation is spatially structured, is lacking for the Cape mole-rat (Georychus capensis). By adopting a range-wide sampling regime and employing two variable mitochondrial markers (cytochrome b and control region), we report on the effects that life-history, population demography and geographic barriers had in shaping genetic variation and population genetic patterns in G. capensis. We also compare our results to information available for the sister taxon of the study species, Bathyergus suillus. Our results show that Georychus capensis exhibits low genetic diversity relative to the concomitantly distributed B. suillus, most likely due to differences in habitat specificity, habitat fragmentation and historical population declines. In addition, the isolated nature of G. capensis populations and low levels of population connectivity has led to small effective population sizes and genetic differentiation, possibly aided by genetic drift. Not surprisingly therefore, G. capensis exhibits pronounced spatial structure across its range in South Africa. Along with geographic distance and demography, other factors shaping the genetic structure of G. capensis include the historical and contemporary impacts of mountains, rivers, sea-level fluctuations and elevation. Given the isolation and differentiation among G. capensis populations, the monotypic genus Georychus may represent a species complex.

Local and regional scale genetic variation in the Cape dune mole-rat, Bathyergus suillus

The distribution of genetic variation is determined through the interaction of life history, morphology and habitat specificity of a species in conjunction with landscape structure. While numerous studies have investigated this interplay of factors in species inhabiting aquatic, riverine, terrestrial, arboreal and saxicolous systems, the fossorial system has remained largely unexplored. In this study we attempt to elucidate the impacts of a subterranean lifestyle coupled with a heterogeneous landscape on genetic partitioning by using a subterranean mammal species, the Cape dune mole-rat (Bathyergus suillus), as our model. Bathyergus suillus is one of a few mammal species endemic to the Cape Floristic Region (CFR) of the Western Cape of South Africa. Its distribution is fragmented by rivers and mountains; both geographic phenomena that may act as geographical barriers to gene-flow. Using two mitochondrial fragments (cytochrome b and control region) as well as nine microsatellite loci, we determined the phylogeographic structure and gene-flow patterns at two different spatial scales (local and regional). Furthermore, we investigated genetic differentiation between populations and applied Bayesian clustering and assignment approaches to our data. Nearly every population formed a genetically unique entity with significant genetic structure evident across geographic barriers such as rivers (Berg, Verlorenvlei, Breede and Gourits Rivers), mountains (Piketberg and Hottentots Holland Mountains) and with geographic distance at both spatial scales. Surprisingly, B. suillus was found to be paraphyletic with respect to its sister species, B. janetta-a result largely overlooked by previous studies on these taxa. A systematic revision of the genus Bathyergus is therefore necessary. This study provides a valuable insight into how the biology, life-history and habitat specificity of animals inhabiting a fossorial system may act in concert with the structure of the surrounding landscape to influence genetic distinctiveness and ultimately speciation.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

Microsatellite Diversity in Natural Populations of Ascomycetous Fungus, Emericella Nidulans, from Different Climatic-Edaphic Conditions in Israel

The genetic divergence of Israeli populations of the soil ascomycetous fungusEmericella nidulanswas studied on regional and local scales using fifteen microsatellite (SSR) markers. The study was performed in the framework of the "Evolution Canyon" research program at the Institute of Evolution, University of Haifa, in three "Evolution Canyons" (ECs): EC I (Mt. Carmel), EC II (western Upper Galilee), and EC III (the southern Negev desert). The first two canyons (EC I and EC II) are located in the northern part of Israel at a distance of 38 km apart; EC III is located southward at a distance of nearly 350 km from the northern ECs. In each canyon,E. nidulansstrains were isolated from opposite slopes and, in EC III, from the valley bottom. All three EC populations ofE. nidulanswere found to be genetically distinct. The estimated genetic divergences correspond to geographical distances and ecological differences between the three studied canyons. On a regional scale, simple sequence repeat (SSR) polymorphism tends to increase with severity of ecological conditions. In general, both environmental parameters (soil moisture and temperature) and genetic factors (predicted number of repeats in SSR markers, distance from marker to centromere, codon evolutionary chronologies, and hydrophobic vs. hydrophilic character of encoded amino acid) influenced genetic diversity ofE. nidulanspopulations.

Lineage-specific differences in the amino acid substitution process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size. To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.

Karyotypic and molecular polymorphisms in Ctenomys torquatus (Rodentia: Ctenomyidae): taxonomic considerations

The rodent genus Ctenomys (tuco-tucos) comprises more than 60 described species, and shows extraordinary inter- and intraspecific karyotypic variation. The most widely distributed species of Ctenomys in Brazil is C. torquatus. Although several cytogenetic studies have been done, the karyotypic variability of this species is still poorly known. In this paper we report two new diploid numbers for C. torquatus: 2n = 40 and 2n = 42, both showing AN = 72. The new distribution limits of C. torquatus here reported include localities in the southern, central and western parts of Rio Grande do Sul (RS) State in southern Brazil. The phylogenetic relationship between C. torquatus from Alegrete, RS, and Ctenomys sp. from Corrientes, Argentina, is described by means of mtDNA cytochrome b analysis. Although both entities share similar karyotypes and sperm morphology, these two species are not phylogenetically close.

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus)

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.

The importance of immune gene variability (MHC) in evolutionary ecology and conservation

Genetic studies have typically inferred the effects of human impact by documenting patterns of genetic differentiation and levels of genetic diversity among potentially isolated populations using selective neutral markers such as mitochondrial control region sequences, microsatellites or single nucleotide polymorphism (SNPs). However, evolutionary relevant and adaptive processes within and between populations can only be reflected by coding genes. In vertebrates, growing evidence suggests that genetic diversity is particularly important at the level of the major histocompatibility complex (MHC). MHC variants influence many important biological traits, including immune recognition, susceptibility to infectious and autoimmune diseases, individual odours, mating preferences, kin recognition, cooperation and pregnancy outcome. These diverse functions and characteristics place genes of the MHC among the best candidates for studies of mechanisms and significance of molecular adaptation in vertebrates. MHC variability is believed to be maintained by pathogen-driven selection, mediated either through heterozygote advantage or frequency-dependent selection. Up to now, most of our knowledge has derived from studies in humans or from model organisms under experimental, laboratory conditions. Empirical support for selective mechanisms in free-ranging animal populations in their natural environment is rare. In this review, I first introduce general information about the structure and function of MHC genes, as well as current hypotheses and concepts concerning the role of selection in the maintenance of MHC polymorphism. The evolutionary forces acting on the genetic diversity in coding and non-coding markers are compared. Then, I summarise empirical support for the functional importance of MHC variability in parasite resistance with emphasis on the evidence derived from free-ranging animal populations investigated in their natural habitat. Finally, I discuss the importance of adaptive genetic variability with respect to human impact and conservation, and implications for future studies.

Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae)

The Chatham Island black robin, Petroica traversi, is a highly inbred, endangered passerine with extremely low levels of variation at hypervariable neutral DNA markers. In this study we investigated variation in major histocompatibility complex (MHC) class II genes in both the black robin and its nonendangered relative, the South Island robin Petroica australis australis. Previous studies have shown that Petroica have at least four expressed class II B MHC genes. In this study, the sequences of introns flanking exon 2 of these loci were characterized to design primers for peptide-binding region (PBR) sequence analysis. Intron sequences were comprised of varying numbers of repeated units, with highly conserved regions immediately flanking exon 2. Polymerase chain reaction primers designed to this region amplified three or four sequences per black robin individual, and eight to 14 sequences per South Island robin individual. MHC genes are fitness-related genes thought to be under balancing selection, so they may be more likely to retain variation in bottlenecked populations. To test this, we compared MHC variation in the black robin with artificially bottlenecked populations of South Island robin, and with their respective source populations, using restriction fragment length polymorphism analyses and DNA sequencing of the PBR. Our results indicate that the black robin is monomorphic at class II B MHC loci, while both source and bottlenecked populations of South Island robin have retained moderate levels of variation. Comparison of MHC variation with minisatellite DNA variation indicates that genetic drift outweighs balancing selection in determining MHC diversity in the bottlenecked populations. However, balancing selection appears to influence MHC diversity over evolutionary timescales, and the effects of gene conversion are evident.

Quantum speciation in Aegilops: molecular cytogenetic evidence from rDNA cluster variability in natural populations

Data are presented on quantum speciation in the Sitopsis section of the genus Aegilops (Poaceae, Monocotyledones). Two small, peripheral, isolated, wild populations of annual cross-pollinated Ae. speltoides and annual self-pollinated Ae. sharonensis are located 30 m apart on different soil types. Despite the close proximity of the two species and their close relatedness, no mixed groups are known. Comparative molecular cytogenetic analysis based on the intrapopulation variability of rRNA-encoding DNA (rDNA) chromosomal patterns of individual Ae. speltoides geno-types revealed an ongoing dynamic process of permanent chromosomal rearrangements. Chromosomal mutations can arise de novo and can be eliminated. Analysis of the progeny of the investigated genotypes testifies that inheritance of de novo rDNA sites happens frequently. Heterologous recombination and/or transposable elements-mediated rDNA transfer seem to be the mechanisms for observed chromosomal repatterning. Consequently, several modified genomic forms, intermediate between Ae. speltoides and Ae. sharonensis, permanently arise in the studied wild population of Ae. speltoides, which make it possible to recognize Ae. sharonensis as a derivative species of Ae. speltoides, as well as to propose rapidness and canalization of quantum speciation in Sitopsis species.

Effects of habitat fragmentation and changes of dispersal behaviour after a recent population decline on the genetic variability of noncoding and coding DNA of a monogamous Malagasy rodent

While interactions among demography, behaviour and genetic structure are well-documented for neutral genetic markers, the role of these parameters and the effects of genetic drift and selection are considerably less well understood in functional genes, such as the major histocompatibility complex (MHC). In this study, the consequences of habitat fragmentation and the effects of a current population decline on noncoding (mitochondrial DNA) and two coding MHC loci (DQA, DRB) with different functional importance were investigated in the small remnant subdivided population of the endangered Malagasy giant jumping rat (Hypogeomys antimena). Both neutral and selective markers revealed a significant genetic differentiation between the two remnant subpopulations. The FST values were much lower in the MHC DQA and DRB genes than in the mitochondrial data. The MHC DRB loci display the effects of both balancing selection (high sequence diversity, four times higher nonsynonymous than synonymous substitutions in the functionally important antigen-binding site positions, twice the average heterozygosity of individual amino acids at the positions identified as part of the antigen-binding site (ABS) than those outside the ABS and nonselective forces including genetic drift. Simultaneously with a current population decline offspring reduced their dispersal distances. No substantial effects were detected within the first 6 years of reduced gene flow in either mitochondrial or MHC markers.

Evolution of genome-phenome diversity under environmental stress

The genomic era revolutionized evolutionary biology. The enigma of genotypic-phenotypic diversity and biodiversity evolution of genes, genomes, phenomes, and biomes, reviewed here, was central in the research program of the Institute of Evolution, University of Haifa, since 1975. We explored the following questions. (i) How much of the genomic and phenomic diversity in nature is adaptive and processed by natural selection? (ii) What is the origin and evolution of adaptation and speciation processes under spatiotemporal variables and stressful macrogeographic and microgeographic environments? We advanced ecological genetics into ecological genomics and analyzed globally ecological, demographic, and life history variables in 1,200 diverse species across life, thousands of populations, and tens of thousands of individuals tested mostly for allozyme and partly for DNA diversity. Likewise, we tested thermal, chemical, climatic, and biotic stresses in several model organisms. Recently, we introduced genetic maps and quantitative trait loci to elucidate the genetic basis of adaptation and speciation. The genome-phenome holistic model was deciphered by the global regressive, progressive, and convergent evolution of subterranean mammals. Our results indicate abundant genotypic and phenotypic diversity in nature. The organization and evolution of molecular and organismal diversity in nature at global, regional, and local scales are nonrandom and structured; display regularities across life; and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection, including diversifying, balancing, cyclical, and purifying selective regimes, interacting with, but ultimately overriding, the effects of mutation, migration, and stochasticity.

Variation on islands: major histocompatibility complex (Mhc) polymorphism in populations of the Australian bush rat

Loss of genetic variation in small, isolated populations is commonly observed at neutral or nearly neutral loci. In this study, the loss of genetic variation was assessed in island populations for a locus of major histocompatibility complex (Mhc), a locus shown to be under the influence of balancing selection. A total of 36 alleles was found at the second exon of RT1.Ba in 14 island and two mainland populations of Rattus fuscipes greyii. Despite this high overall diversity, a substantial lack of variation was observed in the small island populations, with 13 islands supporting only one to two alleles. Two populations, Waldegrave and Williams Islands, showed moderately high levels of heterozygosity (52-56%) which were greater than expected under neutrality, suggesting the action of balancing selection. However, congruence between the level of variation at this Mhc locus and in previous allozyme electrophoresis and mitochondrial DNA studies highlights the dominant influence of genetic drift and population factors, such as bottlenecks and structuring in the founding population, in the loss of genetic variation in these small, isolated populations.

Good genes, oxidative stress and condition-dependent sexual signals

The immune and the detoxication systems of animals are characterized by allelic polymorphisms, which underlie individual differences in ability to combat assaults from pathogens and toxic compounds. Previous studies have shown that females may improve offspring survival by selecting mates on the basis of sexual ornaments and signals that honestly reveal health. In many cases the expression of these ornaments appears to be particularly sensitive to oxidative stress. Activated immune and detoxication systems often generate oxidative stress by an extensive production of reactive metabolites and free radicals. Given that tolerance or resistance to toxic compounds and pathogens can be inherited, female choice should promote the evolution of male ornaments that reliably reveal the status of the bearers' level of oxidative stress. Hence, oxidative stress may be one important agent linking the expression of sexual ornaments to genetic variation in fitness-related traits, thus promoting the evolution of female mate choice and male sexual ornamentation, a controversial issue in evolutionary biology ever since Darwin.