Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

Unraveling Selection in the Mitochondrial Genome of Drosophila

We examine mitochondrial DNA variation at the cytochrome b locus within and between three species of Drosophila to determine whether patterns of variation conform to the predictions of neutral molecular evolution. The entire 1137-bp cytochrome b locus was sequenced in 16 lines of Drosophila melanogaster, 18 lines of Drosophila simulans and 13 lines of Drosophila yakuba. Patterns of variation depart from neutrality by several test criteria. Analysis of the evolutionary clock hypothesis shows unequal rates of change along D. simulans lineages. A comparison within and between species of the ratio of amino acid replacement change to synonymous change reveals a relative excess of amino acid replacement polymorphism compared to the neutral prediction, suggestive of slightly deleterious or diversifying selection. There is evidence for excess homozygosity in our world wide sample of D. melanogaster and D. simulans alleles, as well as a reduction in the number of segregating sites in D. simulans, indicative of selective sweeps. Furthermore, a test of neutrality for codon usage shows the direction of mutations at third positions differs among different topological regions of the gene tree. The analyses indicate that molecular variation and evolution of mtDNA are governed by many of the same selective forces that have been shown to govern nuclear genome evolution and suggest caution be taken in the use of mtDNA as a ``neutral' molecular marker.     

Genetic localization and sequential electrophoresis of glucose-6-phosphate dehydrogenase in Drosophila melanogaster

Studies were undertaken to investigate two critical aspects of the glucose-6-phosphate dehydrogenase polymorphism in Drosophila melanogaster. The first investigation unequivocally maps the genetic site of the G6PD locus to the X chromosome. The second study subjects a set of isochromosomal lines to sequential electrophoresis in an attempt to uncover common molecular heterogeneity within the global polymorphism, assuming that this variation may have gone undetected under conventional electrophoretic conditions. The genetic site was mapped following the segregation of the two common electrophoretic alleles, a so-called null allele, and two rare electrophoretic variants. From the pooled results, the Zw locus mapped to 62.9 on the X chromosome relative to the flanking markers car (at 62.5) and sw (at 64.7). A set of 126 iso-X chromosomal lines of diverse geographic origin was subjected to sequential electrophoresis under three different acrylamide conditions in addition to the conventional starch electrophoretic system. No additional variation beyond the common diallele polymorphism was seen.     

Spatial and environmental correlates of intraspecific morphological variation in three species of passerine birds from the Purus–Madeira interfluvium,Central Amazonia

Biogeographic studies in Amazonia typically describe biodiversity across interfluvia, rarely within them, where geographic variability in morphological traits might be observed. We tested for intraspecific phenotypic variation in three bird species within the Purus–Madeira interfluvium (Central Amazon) and whether phenotypes were correlated with environmental heterogeneity or geographic distance among sites. We compared coloration indexes derived from reflectance spectra and morphometrics of up to five adult individuals of each sex among 11 sites within the interfluvium and contrasted them with proxies for geographic distance and environmental variation (tree basal area and bird community). Environmental heterogeneity was minimally spatially autocorrelated, and there were no obvious geographical barriers to dispersal in the study region. The null hypothesis was that we would see either no phenotypic variation or random variation that was not explained by the tested variables. Half of the cases analyzed showed intraspecific morphological variation. Coloration varied more frequently than morphometrics, and color was better explained by environmental heterogeneity, particularly in males, whereas brightness also varied with geographic distance. Geographic distance explained the only case of variation in morphometrics. Our results indicate that coloration, particularly plumage brightness, is more labile than morphometric traits and that plumage color might be under stronger effects of local adaptation than brightness, which also seems to be under effects of neutral drift and gene flow among populations. Higher frequencies of association between male coloration and the environment suggest a role of non-arbitrary mechanisms of sexual selection on the expression of male phenotypes, whereas arbitrary intersexual selection might explain the randomly distributed variation that is not explained by environmental heterogeneity or geographic distance. We revealed intraspecific phenotypic variation in a spatial extent usually not considered in biogeographic studies in the Amazon and demonstrate that both local adaptation and neutral drift are important to explain intraspecific trait diversification at this geographical scale.     

Neutral and Nonneutral Mitochondrial Genetic Variation in Deep-Sea Clams from the Family Vesicomyidae

Nucleotide sequences at two mitochondrial genes from 57 individuals representing eight species of deep-sea clams (Vesicomyidae) were examined for variation consistent with the neutral model of molecular evolution. One gene, cytochrome oxidase subunit I (COI), deviated from the expectations of neutrality by containing an excess of intraspecific nonsynonymous polymorphism. Additionally, one species, Calyptogena kilmeri, showed a significant excess of rare polymorphism specifically at the COI locus. In contrast, a second mitochondrial gene, the large-subunit 16S ribosomal RNA gene (16S), showed little deviation from neutrality either between or within species. Together, COI and 16S show no deviation from neutral expectations by the HKA test, produce congruent phylogenetic relationships between species, and show correlated numbers of fixed differences between species and polymorphism within species. These patterns of both neutral and nonneutral evolution within the mitochondrial genome are most consistent with a model where intraspecific nonsynonymous polymorphism at COI is near neutrality. In addition to examining the forces of molecular evolution, we extend hypotheses about interspecific relationships within this family for geographical locations previously unexamined by molecular methods including habitats near the Middle Atlantic, the Aleutian Trench, and Costa Rica. Received: 10 March 1999 / Accepted: 13 September 1999     

Molecular Population Genetics of Sex Determination Genes: The Transformer Gene of Drosophila Melanogaster

The transformer locus (tra) produces an RNA processing protein that alternatively splices the doublesex pre-mRNA in the sex determination hierarchy of Drosophila melanogaster. Comparisons of the tra coding region among Drosophila species have revealed an unusually high degree of divergence in synonymous and nonsynonymous sites. In this study, we tested the hypothesis that the tra gene will be polymorphic in synonymous and nonsynonymous sites within species by investigating nucleotide sequence variation in eleven tra alleles within D. melanogaster. Of the 1063 nucleotides examined, two synonymous sites were polymorphic and no amino acid variation was detected. Three statistical tests were used to detect departures from an equilibrium neutral model. Two tests failed to reject a neutral model of molecular evolution because of low statisitical power associated with low levels of genetic variation (Tajima/Fu and Li). The Hudson, Kreitman, and Aguade test rejected a neutral model when the tra region was compared to the 5'-flanking region of alcohol dehydrogenase (Adh). The lack of variability in the tra gene is consistent with a recent selective sweep of a beneficial allele in or near the tra locus.     

Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids

Phenotypically diverse Lake Malawi cichlids exhibit similar genomes. The extensive sharing of genetic polymorphism among forms has both intrigued and frustrated biologists trying to understand the nature of diversity in this and other rapidly evolving systems. Shared polymorphism might result from hybridization and/or the retention of ancestrally polymorphic alleles. To examine these alternatives, we used new genomic tools to characterize genetic differentiation in widespread, geographically structured populations of Labeotropheus fuelleborni and Metriaclima zebra. These phenotypically distinct species share mitochondrial DNA (mtDNA) haplotypes and show greater mtDNA differentiation among localities than between species. However, Bayesian analysis of nuclear single nucleotide polymorphism (SNP) data revealed two distinct genetic clusters corresponding perfectly to morphologically diagnosed L. fuelleborni and M. zebra. This result is a function of the resolving power of the multi‐locus dataset, not a conflict between nuclear and mitochondrial partitions. Locus‐by‐locus analysis showed that mtDNA differentiation between species (F_CT) was nearly identical to the median single‐locus SNP F_CT. Finally, we asked whether there is evidence for gene flow at sites of co‐occurrence. We used simulations to generate a null distribution for the level of differentiation between co‐occurring populations of L. fuelleborni and M. zebra expected if there was no hybridization. The null hypothesis was rejected for the SNP data; populations that co‐occur at rock reef sites were slightly more similar than expected by chance, suggesting recent gene flow. The coupling of numerous independent markers with extensive geographic sampling and simulations utilized here provides a framework for assessing the prevalence of gene flow in recently diverged species.     

Neutral and Non-Neutral Evolution of Drosophila Mitochondrial DNA

To test hypotheses of neutral evolution of mitochondrial DNA (mtDNA), nucleotide sequences were determined for 1515 base pairs of the NADH dehydrogenase subunit 5 (ND5) gene in the mitochondrial DNA of 29 lines of Drosophila melanogaster and 9 lines of its sibling species Drosophila simulans. In contrast to the patterns for nuclear genes, where D. melanogaster generally exhibits much less nucleotide polymorphism, the number of segregating sites was slightly higher in a global sample of nine ND5 sequences in D. melanogaster (s = 8) than in the nine lines of D. simulans (s = 6). When compared to variation at nuclear loci, the mtDNA variation in D. melanogaster does not depart from neutral expectations. The ND5 sequences in D. simulans, however, show fewer than half the number of variable sites expected under neutrality when compared to sequences from the period locus. While this reduction in variation is not significant at the 5% level, HKA tests with published restriction data for mtDNA in D. simulans do show a significant reduction of variation suggesting a selective sweep of variation in the mtDNA in this species. Tests of neutral evolution based on the ratios of synonymous and replacement polymorphism and divergence are generally consistent with neutral expectations, although a significant excess of amino acid polymorphism within both species is localized in one region of the protein. The rate of mtDNA evolution has been faster in D. melanogaster than in D. simulans and the population structure of mtDNA is distinct in these species. The data reveal how different rates of mtDNA evolution between species and different histories of neutral and adaptive evolution within species can compromise historical inferences in population and evolutionary biology.     

Nucleotide polymorphism and natural selection at the pantophysin (Pan I) locus in the Atlantic cod, Gadus morhua (L.)

Molecular studies of nucleotide sequence variation have rarely attempted to test hypotheses related to geographically varying patterns of natural selection. The present study tested the role of spatially varying selection in producing significant linkage disequilibrium and large differences in the frequencies of two common alleles at the pantophysin (Pan I) locus among five populations of the Atlantic cod, Gadus morhua. Nucleotide sequences of 124 Pan I alleles showed strong evidence for an unusual mix of balancing and directional selection but no evidence of stable geographically varying selection. The alleles were highly divergent at both the nucleotide level (differing on average by 19 mutations) and at amino acid level (each having experienced three amino acid substitutions since diverging from a common ancestral allele). All six amino acid substitutions occurred in a 56-residue intravesicular loop (IV1 domain) of the vesicle protein and each involved a radical change. An analysis of molecular variation revealed significant heterogeneity in the frequencies of recently derived mutations segregating within both allelic classes, suggesting that two selective sweeps may be presently occurring among populations. The dynamic nature of the Pan I polymorphism in G. morhua and clear departure from equilibrium conditions invalidate a simple model of spatially varying selection.     

Nucleotide Variation at the Gpdh Locus in the Genus Drosophila

The Gpdh locus was sequenced in a broad range of Drosophila species. In contrast to the extreme evolutionary constraint seen at the amino acid level, the synonymous sites evolve at rates comparable to those of other genes. Gpdh nucleotide sequences were used to infer a phylogenetic tree, and the relationships among the species of the obscura group were examined in detail. A survey of nucleotide polymorphism within D. pseudoobscura revealed no amino acid variation in this species. Applying a modified McDonald-Kreitman test, the amino acid divergence between species in the obscura group does not appear to be excessive, implying that drift is adequate to explain the patterns of amino acid change at this locus. In addition, the level of polymorphism at the Gpdh locus in D. pseudoobscura is comparable to that found at other loci, as determined by a Hudson-Kreitman-Aguade test. Thus, the pattern of nucleotide variation within and between species at the Gpdh locus is consistent with a neutral model.     

Testing the generalized neutrality hypothesis

Various tests of the hypothesis of selective neutrality based on gene frequency are now available. These tests take as null hypothesis the concept of “strict neutrality”: all new mutants are required to be selectively identical to each other. For evolutionary questions, however, (as opposed to those of genetic polymorphism), a wider null hypothesis might be of interest. Since deleterious alleles have essentially no evolutionary importance, one might wish to test the null hypothesis that only neutral or deleterious mutations occur. The principal alternative to this hypothesis is that there exists heterotic selection of some form for some alleles tending to maintain a level of genetic polymorphism higher than that under neutrality. In this paper an assessment is made of the usefulness of a test of strict neutrality first proposed by this author (Ewens, 1972) as a test of null hypothesis of “generalized neutrality,” i.e. that only neutral or deleterious alleles occur. At the same time some remarks will be made about estimation of the fundamental parameter θ defining these processes.     

Reduced nucleotide variability at an androgen-binding protein locus (Abpa) in house mice: evidence for positive natural selection.

Previous work has shown that the gene for the alpha subunit of androgen-binding protein, Abpa, may be involved in premating isolation between different subspecies of the house mouse, Mus musculus. We investigated patterns of DNA sequence variation at Abpa within and between species of mice to test several predictions of a model of neutral molecular evolution. Intraspecific variation among 10 Mus musculus domesticus alleles was compared with divergence between M. m. domesticus and M. caroli for Abpa and two X-linked genes, Glra2 and Amg. No variation was observed at Abpa within M. m. domesticus. The ratio of polymorphism to divergence was significantly lower at Abpa than at Glra2 and Amg, despite the fact that all three genes experience similar rates of recombination. Interspecific comparisons among M. m. domesticus, Mus musculus musculus, Mus musculus castaneus, Mus spretus, Mus spicilegus, and Mus caroli revealed that the ratio of nonsynonymous substitutions to synonymous substitutions on a per-site basis (Ka/Ks) was generally greater than one. The combined observations of no variation at Abpa within M. m: domesticus and uniformly high Ka/Ks values between species suggest that positive directional selection has acted recently at this locus.     

Abdominal pigmentation variation in drosophila polymorpha: geographic variation in the trait, and underlying phylogeography

Drosophila polymorpha is a widespread species that exhibits abdominal pigmentation variation throughout its range. To gain insight into this variation we combined phenotypic and genotypic data to test a series of nested hypotheses. First, we tested the null hypothesis that geographic variation in pigmentation is due to neutral factors. We used nested clade analysis to examine the distribution of haplotypes from a nuclear and a mitochrondrial locus. Restricted gene flow via isolation by distance, the primary inference of this phylogeographic analysis, was then used to generate and test the hypothesis of increasing average abdominal pigmentation difference with increasing geographic distance. We found no correlation between geographic distance and phenotypic distance. We then tested the hypothesis that pigmentation is affected by environmental differences among localities. We found a significant effect of habitat type on the average abdominal pigmentation phenotype of different localities. Finally, we tested the hypothesis that pigmentation in D. polymorpha is associated with desiccation resistance. We found that dark individuals of both sexes survived significantly longer in a desiccating environment than light individuals. These patterns combined lead us to hypothesize that abdominal pigmentation variation in D. polymorpha is important in mediating the organism's interactions with local ecological factors.     

Trans-specific evolution of opsin alleles and the maintenance of trichromatic colour vision in Callitrichine primates

Many New World (NW) primates possess a remarkable polymorphism in an X-linked locus, which encodes for the visual pigments (opsins) used for colour vision. Females that are heterozygous for opsin alleles of different spectral sensitivity at this locus have trichromatic colour vision, whereas homozygous females and males are dichromatic, with poor colour discrimination in the red-green range. Here we describe an extensive survey of allelic variation in both exons and introns at this locus within and among species of the Callitrichines (marmosets and tamarins). All five genera of Callitrichines have the X-linked polymorphism, and only the three functional allelic classes described previously (with maximum wavelength sensitivities at about 543 nm, 556 nm and 563 nm) were found among the 16 species and 233 or more X-chromosomes sampled. In spite of the homogenizing effects of gene conversion, phylogenetic analyses provide direct evidence for trans-specific evolution of alleles over time periods of at least 5-6 million years, and up to 14 million years (estimated from independent phylogenies). These conclusions are supported by the distribution of insertions and deletions in introns. The maintenance of polymorphism over these time periods requires an adaptive explanation, which must involve a heterozygote advantage for trichromats. The lack of detection of alleles that are recombinant for spectral sensitivity suggests that such alleles are suboptimal. The two main hypotheses for the selective advantage of trichromacy in primates are frugivory for ripe fruits and folivory for young leaves. The latter can be discounted in Callitrichines, as they are not folivorous.     

Environmental heterogeneity generates opposite gene‐by‐environment interactions for two fitness‐related traits within a population

Theory predicts that environmental heterogeneity offers a potential solution to the maintenance of genetic variation within populations, but empirical evidence remains sparse. The live‐bearing fish Xiphophorus variatus exhibits polymorphism at a single locus, with different alleles resulting in up to five distinct melanistic “tailspot” patterns within populations. We investigated the effects of heterogeneity in two ubiquitous environmental variables (temperature and food availability) on two fitness‐related traits (upper thermal limits and body condition) in two different tailspot types (wild‐type and upper cut crescent). We found gene‐by‐environment (G × E) interactions between tailspot type and food level affecting upper thermal limits (UTL), as well as between tailspot type and thermal environment affecting body condition. Exploring mechanistic bases underlying these G × E patterns, we found no differences between tailspot types in hsp70 gene expression despite significant overall increases in expression under both thermal and food stress. Similarly, there was no difference in routine metabolic rates between the tailspot types. The reversal of relative performance of the two tailspot types under different environmental conditions revealed a mechanism by which environmental heterogeneity can balance polymorphism within populations through selection on different fitness‐related traits.     

DNA Polymorphism at the Pgi Locus of a Wild Yam, Dioscorea Tokoro

To study the origin and maintenance mechanisms of the PGI allozyme polymorphism of a wild plant, Dioscorea tokoro, DNA sequences of the entire coding region (1701 bp) and two intronic regions (total 2049 bp) of the Pgi gene as well as a part of the Adh gene (590 bp) were analyzed. Two replacement substitutions were revealed to be responsible for the differentiation of three allozymes alleles (Pgi-a, Pgi-b and Pgi-c) that occur in natural population in intermediate frequencies. Interspecific comparison of DNA sequences identified Pgi-b as the oldest allele, from which two other alleles were derived probably within the last 150,000 years. The level of DNA polymorphism at D. tokoro Pgi locus was low. No elevated level of DNA polymorphism was detected in the close vicinity of the two replacement sites differentiating the three allozymes. Departures from the neutral mutation hypothesis were detected by Fu and Li's and MK tests. The observed patterns of DNA polymorphism are explainable by both (1) the neutral mutation hypothesis with an assumption of small effective size of D. tokoro population, and (2) the positive selection hypothesis that the allele frequencies of Pgi-a and Pgi-c have increased in a short time by their selective advantages.     

Contrasting patterns of variation in MHC loci in the Alpine newt

Major histocompatibility complex (MHC) genes are essential in pathogen recognition and triggering an adaptive immune response. Although they are the most polymorphic genes in vertebrates, very little information on MHC variation and patterns of evolution are available for amphibians, a group known to be declining rapidly worldwide. As infectious diseases are invoked in the declines, information on MHC variation should contribute to devising appropriate conservation strategies. In this study, we examined MHC variation in 149 Alpine newts ( Mesotriton alpestris ) from three allopatric population groups in Poland at the northeastern margin of the distribution of this species. The genetic distinctiveness of the population groups has previously been shown by studies of skin graft rejection, allozymes and microsatellites. Two putative expressed MHC II loci with contrasting levels of variation and clear evidence of gene conversion/recombination between them were detected. The Meal-DAB locus is highly polymorphic (37 alleles), and shows evidence of historical positive selection for amino acid replacements and substantial geographical differentiation in allelic richness. On the contrary, the Meal-DBB locus exhibits low polymorphism (three alleles differing by up to two synonymous substitutions) and a uniform distribution of three alleles among geographical regions. The uniform frequencies of the presumptively neutral Meal-DBB alleles may be explained by linkage to Meal-DAB . We found differences in allelic richness in Meal-DAB between regions, consistent with the hypothesis that genetic drift prevails with increasing distance from glacial refugia. Pseudogene loci appear to have evolved neutrally. The level of DAB variation correlated with variation in microsatellite loci, implying that selection and drift interplayed to produce the pattern of MHC variation observed in marginal populations of the Alpine newt.     

The niche-width variation hypothesis reconfirmed: Validation by genetic diversity in the sessile intertidal cirripedes Chthamalus stellatus and Euraphia depressa (Crustacea,Chthamalidae)

Genetic diversity in a pair of cirripede species, based on electrophoretic analysis of 25 gene loci, is higher for Chthamalus stellatus, the species with a broader biogeographical distribution, and hence, with a broader ecological niche, than that of Euraphia depressa. Comparing the genetic diversity within E. depressa we again report a higher genetic diversity among the specimens exposed to solar radiation (living in the wider ecological niche) versus the specimens confined to the sheltered dark environments in caves or under boulders. The indices of genetic diversity used in this study are polymorphism, heterozygosity, mean number of alleles per locus and gene diversity. We conclude that the niche-width variation hypothesis has been confirmed in this pair of species: genetic diversity is positively correlated with niche breadth.     

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.