Signatures of demographic history and natural selection in the human major histocompatibility complex Loci

Many lines of evidence show that several HLA loci have experienced balancing selection. However, distinguishing among demographic and selective explanations for patterns of variation observed with HLA genes remains a challenge. In this study we address this issue using data from a diverse set of human populations at six classical HLA loci and, employing a comparative genomics approach, contrast results for HLA loci to those for non-HLA markers. Using a variety of analytic methods, we confirm and extend evidence for selection acting on several HLA loci. We find that allele frequency distributions for four of the six HLA loci deviate from neutral expectations and show that this is unlikely to be explained solely by demographic factors. Other features of HLA variation are explained in part by demographic history, including decreased heterozygosity and increased LD for populations at greater distances from Africa and a similar apportionment of genetic variation for HLA loci compared to putatively neutral non-HLA loci. On the basis of contrasts among different HLA loci and between HLA and non-HLA loci, we conclude that HLA loci bear detectable signatures of both natural selection and demographic history.     

Evolutionary genetics of Metridium senile. II. Geographic patterns of allozyme variation

Electrophoretic surveys have demonstrated that populations of the sea anemone Metridium senile along the northeast coast of the United States are polymorphic at four enzyme loci. Phosphoglucose isomerase (PGI) has two alleles in most populations, phosphoglucomutase (PGM) has three alleles, and two leucine aminopeptidase loci have two common alleles each. Phosphoglucose isomerase displays clinal variation and an apparent association with environmental temperature. Phosphoglucomutase shows clinal variation north of Cape Cod for two of the three alleles, while the two leucine aminopeptidase loci are not clinal. All loci show a great deal of variation in populations on Cape Cod, but there is no apparent systematic pattern to this variation. Temperature may be a selective agent in the maintenance of the PGI and PGM clines, although other possibilities cannot presently be completely excluded.Supported by Grant T-4 from the Health Research and Services Foundation, NSF DEB77-14442, NIH GM25809, and NIH GM28024.     

Extensive ribosomal DNA amplification during Andean common bean (Phaseolus vulgaris L.) evolution

The extent of 5S and 45S ribosomal DNA (rDNA) variation was investigated in wild and domesticated common beans (Phaseolus vulgaris) chosen to represent the known genetic diversity of the species. 5S and 45S rDNA probes were localized on mitotic chromosomes of 37 accessions by fluorescent in situ hybridization (FISH). The two 5S rDNA loci were largely conserved within the species, whereas a high variation in the number of 45S rDNA loci and changes in position of loci and number of repeats per locus were observed. Domesticated accessions from the Mesoamerican gene pool frequently had three 45S rDNA loci per haploid genome, and rarely four. Domesticated accessions from Andean gene pool, particularly from the race Peru, showed six, seven, eight or nine loci, but seven loci were found in all three races of this gene pool. Between three and eight loci were observed in accessions resulting from crosses between Andean and Mesoamerican genotypes. The presence of two to eight 45S rDNA loci in wild common beans from different geographic locations indicates that the 45S rDNA amplification observed in the Andean lineage took place before domestication. Our data suggest that ectopic recombination between terminal chromosomal regions might be the mechanism responsible for this variation.     

THE ECOLOGICAL GENETICS OF INTRODUCED POPULATIONS OF THE GIANT TOAD,BUFO MARINUS. III. GEOGRAPHICAL PATTERNS OF VARIATION

The allele frequencies at ten polymorphic loci are described from 31 Bufo marinus populations in the Moreton Bay region in southeastern Queensland, Australia and the variation of these is found to be non-random in all cases. The pattern of non-randomness varies among loci, being clinal in two instances. The allele frequencies at the same ten loci are also described for 12 populations sampled from throughout B. marinus' Australian range. The frequency variation on this larger geographical scale is non-random at all but two loci (Mpi and Hbdh) and also varies among loci, in this case being clinal in four instances. In both cases, the patterns of variation are most reasonably explained as having resulted from genetic drift occurring during the recent range expansion which B. marinus is known to have experienced in Australia. It seems that natural selection has played little, if any, role in generating the observed gene frequency patterns. These results emphasize the need for caution in interpreting geographical patterns of variation. They show that even when clinal patterns exist at some loci but not at others, one cannot conclude that the patterns result from natural selection, unless the demographic histories of the studied populations are known and are inconsistent with the alternative hypothesis that the patterns result from genetic drift.     

Genetic evidence for population expansion in Hydrotaea irritans (Fallèn) (Diptera: Muscidae)

Allozyme variation in the sheep headfly Hydrotaea irritans was studied on two spatial scales. Geographic variation among seven Danish and one Dutch population revealed significant but rather low genetic differentiation with F _ST = 0.01 over all loci. The Dutch population was on average not more different from the Danish populations than the Danish populations from each other. Allele frequencies were very skewed with the most common allele always exceeding 0.85 and usually 0.9 in frequency, but with many rare alleles at some loci. Tests for neutrality of the variation at the nine polymorphic loci revealed highly significant deviations from expected homozygosity in this species, which was not found in a comparative analysis of allozyme variation at similar loci of seven other Hydrotaea species. To explain the peculiar observed pattern of allozyme variation in H. irritans , it is suggested that this species has successfully expanded its range and spread through northern and central Europe in the recent past. Alternatively, H. irritans may have recently invaded a new niche, resulting in increased abundance of the species and subsequent dispersal to former areas of the species distribution.     

Genetic variation at allozyme and RAPD loci in sessile oak Quercus petraea (Matt.) Liebl.: the role of history and geography

The nuclear genetic variation within and among 21 populations of sessile oak was estimated at 31 RAPD loci in conjunction with previous estimates of variation at eight allozyme loci. The aim of the study was to assess the relative role of isolation-by-distance and postglacial history on patterns of nuclear variation. Because of its small effective population size and maternal transmission, the chloroplast genome is a good marker of population history. Both kinds of nuclear variation (RAPD and allozyme) were therefore compared, first, to the geographical distances among populations and, secondly, to chloroplast DNA restriction polymorphism in the same populations. Multiple Mantel tests were used for this purpose. Although RAPDs revealed less genetic diversity than allozymes, levels of genetic differentiation ( G _ST) were identical. The standard genetic distance calculated at all RAPD loci was correlated with geographical distances but not with the genetic distance calculated from chloroplast DNA data. Conversely, allozyme variation was correlated with chloroplast DNA variation, but not with geography. Possibly, divergent selection at two allozyme loci during the glacial period could explain this pattern. Because of its greater number of loci assayed, RAPDs probably provided a less biased picture of the relative role of geography and history.     

Distribution and biomass of microphytes measured by benthic chlorophyll a in a tropical lagoon (New Caledonia,South Pacific)

Starch-gel electrophoretic techniques were applied to the investigation of molecular genetic variation in populations of the rough periwinkle Littorina saxatilis. The investigation comprised two phases: a) technique development to resolve as many as possible of the allozyme loci reported in the literature as having been screened in the genus Littorina and in Melarhaphe neritoides; b) the use of these loci to assess levels of genetic variation in and patterns of genetic differentiation among populations of L. saxatilis from a relatively isolated group of populations from Galway Bay, Ireland. More than 43 allozyme loci (of which four were screened for the first time here in this species), coding for 37 enzymes, were investigated and thirteen of these loci (including two loci screened for the first time here) were found to be variable and reliably scorable. Samples from five pairs of transects were collected from Inismór, Aran Islands, from sites with known exposure levels; one transect within each pair was collected from an exposed site and the other from a nearby, but relatively sheltered site. UPGMA for eleven loci, (ARK and PGDH were excluded from cluster and FST analysis as they were unscorable in a few samples), showed that the samples cluster mostly by pair, reflecting their geographic origin and is indicative of little gene flow between populations. Levels of population differentiation were high among samples from the top of the shore, but unusually so at AAT-1 which showed nearly three times the mean FST value for the eleven loci. There was also a significant regression of frequency of AAT-1100 against level of exposure. In addition, among midshore samples, there was a consistently higher frequency of AAT-1100 in sheltered habitats. These results support the findings of others, indicating that this locus may be subject to natural selection.     

Enzyme variation in the brook lamprey,Lampetra planeri (Bloch), a member of the vertebrate group Agnatha

A natural population of the brook lamprey, Lampetra planeri, was assayed for electrophoretically detectable variation at 30 enzyme loci. The mean heterozygosity per locus of this primitive vertebrate, a member of the Agnatha, was found to be 0.076 ± 0.031, a value similar to those recorded for other vertebrates. The high chromosome numbers recorded for this and related species have been attributed to polyploidy, but our studies do not indicate the existence of large numbers of duplicated loci. Indeed, several enzymes that are encoded by duplicate loci in other vertebrate species appear to be encoded by single loci in the lamprey. It is suggested that studies on the biology and taxonomy of lampreys will benefit greatly from an electrophoretic approach.     

Population genetic consequences of feeding habits in some forest lepidoptera

By surveying variation at allozyme loci in several phytophagous lepidopteran species (Geometridae), we have tested two hypotheses about the relationship of genetic variation to environmental heterogeneity: (1) that allozyme polymorphisms may exist because of associations between genotypes and "niches" (different host plants, in this instance), and (2) that the overall genetic variation of a species is correlated with environmental heterogeneity (or breadth of the species' overall ecological niche).—Genetic differentiation among samples of oligophagous or polyphagous species taken from different host species was observed in one of three species, at only one of seven polymorphic loci. The data thus provide no evidence for pronounced genetic substructuring, or "host race" formation in these sexually reproducing species, although host plant-genotype associations in a parthenogenetic moth give evidence of the potential for diversifying selection.—In a comparison of allozyme variation in polyphagous ("generalized") and oligophagous ("specialized") species, heterozygosity appeared to be higher in specialized species, at all polymorphic loci but one. It is possible that this unexpected result arises from a functional relation between breadth of diet and genetic variation.     

A sustained dietary change increases epigenetic variation in isogenic mice

Epigenetic changes can be induced by adverse environmental exposures, such as nutritional imbalance, but little is known about the nature or extent of these changes. Here we have explored the epigenomic effects of a sustained nutritional change, excess dietary methyl donors, by assessing genomic CpG methylation patterns in isogenic mice exposed for one or six generations. We find stochastic variation in methylation levels at many loci; exposure to methyl donors increases the magnitude of this variation and the number of variable loci. Several gene ontology categories are significantly overrepresented in genes proximal to these methylation-variable loci, suggesting that certain pathways are susceptible to environmental influence on their epigenetic states. Long-term exposure to the diet (six generations) results in a larger number of loci exhibiting epigenetic variability, suggesting that some of the induced changes are heritable. This finding presents the possibility that epigenetic variation within populations can be induced by environmental change, providing a vehicle for disease predisposition and possibly a substrate for natural selection.     

A Reanalysis of Protein Polymorphism in Drosophila Melanogaster, D. Simulans, D. Sechellia and D. Mauritiana: Effects of Population Size and Selection

Comparison of synonymous and nonsynonymous variation/substitution within and between species at individual genes has become a widely used general approach to detect the effect of selection versus drift. The sibling species group comprised of two cosmopolitan (Drosophila melanogaster and Drosophila simulans) and two island (Drosophila mauritiana and Drosophila sechellia) species has become a model system for such studies. In the present study we reanalyzed the pattern of protein variation in these species, and the results were compared against the patterns of nucleotide variation obtained from the literature, mostly available for melanogaster and simulans. We have mainly focused on the contrasting patterns of variation between the cosmopolitan pair. The results can be summarized as follows: (1) As expected the island species D. mauritiana and D. sechellia showed much less variation than the cosmopolitan species D. melanogaster and D. simulans. (2) The chromosome 2 showed significantly less variation than chromosome 3 and X in all four species which may indicate effects of past selective sweeps. (3) In contrast to its overall low variation, D. mauritiana showed highest variation for X-linked loci which may indicate introgression from its sibling, D. simulans. (4) An average population of D. simulans was as heterozygous as that of D. melanogaster (14.4% v.s. 13.9%) but the difference was large and significant when considering only polymorphic loci (37.2% v.s. 26.1%). (5) The species-wise pooled populations of these two species showed similar results (all loci = 18.3% v.s. 20.0%, polymorphic loci = 47.2% v.s. 37.6%). (6) An average population of D. simulans had more low-frequency alleles than D. melanogaster, and the D. simulans alleles were found widely distributed in all populations whereas the D. melanogaster alleles were limited to local populations. As a results of this, pooled populations of D. melanogaster showed more polymorphic loci than those of D. simulans (48.0% v.s. 32.0%) but the difference was reduced when the comparison was made on the basis of an average population (29.1% v.s. 21.4%). (7) While the allele frequency distributions within populations were nonsignificant in both D. melanogaster and D. simulans, melanogaster had fewer than simulans, but more than expected from the neutral theory, low frequency alleles. (8) Diallelic loci with the second allele with a frequency less than 20% had similar frequencies in all four species but those with the second allele with a frequency higher than 20% were limited to only melanogaster the latter group of loci have clinal (latitudinal) patterns of variation indicative of balancing selection. (9) The comparison of D. simulans/D. melanogaster protein variation gave a ratio of 1.04 for all loci and 1.42 for polymorphic loci, against a ratio of approximately 2-fold difference for silent nucleotide sites. This suggests that the species ratios of protein and silent nucleotide polymorphism are too close to call for selective difference between silent and allozyme variation in D. simulans. In conclusion, the contrasting levels of allozyme polymorphism, distribution of rare alleles, number of diallelic loci and the patterns of geographic differentiation between the two species suggest the role of natural selection in D. melanogaster, and of possibly ancient population structure and recent worldwide migration in D. simulans. Population size differences alone are insufficient as an explanation for the patterns of variation between these two species.     

Generating single-copy nuclear gene data for a recent adaptive radiation

Recent adaptive radiations provide an exceptional opportunity to understand the processes of speciation and adaptation. However, reconstructing the phylogenetic history of recent and rapidly evolving clades often requires the use of multiple, independent gene genealogies. Nuclear introns are an obvious source of the necessary data but their use is often limited because degenerate primers can amplify paralogous loci. To identify PCR primers for a large number of loci in an especially rapid adaptive radiation, that of the flowering plant genus Aquilegia (Ranunculaceae), we developed an efficient method for amplifying multiple single-copy nuclear loci by sequencing a modest number of clones from a cDNA library and designing PCR primers; with one primer anchored in the 3' untranslated region (3'-UTR) and one primer in the coding region of each gene. Variation between paralogous loci evolves more quickly in 3'-UTR regions compared to adjacent exons, and therefore we achieved high specificity for isolating orthologous loci. Furthermore, we were able to identify genes containing large introns by amplifying genes from genomic DNA and comparing the PCR product size to that predicted from their cDNA sequence. In Aquilegia eight out of eleven loci were isolated with this method and six of these loci had introns. Among four genes sequenced for samples spanning the phylogenetic breadth of the genus, we found sequence variation at levels similar to that observed in ITS, further supporting the recent and rapid radiation in Aquilegia. We assessed the orthology of amplification products by phylogenetic congruence among loci, the presence of two well established phylogenetic relationships, and similarity among loci for levels of sequence variation. Higher levels of variation among samples for one locus suggest possible paralogy. Overall, this method provides an efficient means of isolating predominantly single-copy loci from both low and high-copy gene families, providing ample nuclear variation for reconstructing species-level phylogenies in non-model taxa.     

DNA polymorphism in Cab locus of tomato induced by tissue culture.

Plants were regenerated from callus induced from leaf disc explants of a tomato F1 hybrid heterozygous for three marker loci anthocyaninless (a), without anthocyanin (aw), and hairless (hl). Regenerants were studied for somaclonal variation at the phenotypic level by scoring for variation in the marker loci, and at the DNA level by probing geomic DNA blots with a chlorophyll a/b binding protein (Cab-3C) cDNA sequence. While no variation was observed at the phenotypic level in over 950 somaclones studied, DNA polymorphism for the Cab locus could be detected in two out of 17 somaclones tested. Tissue culture induced variation at the phenotypic level for specific loci is very low (less than 0.001 for a, aw or hl) but DNA sequence changes are induced at much greater frequency (approximately 0.1 for a multicopy gene family such as Cab).     

Statistical properties of the variation at linked microsatellite loci: implications for the history of human Y chromosomes [published erratum appears in Mol Biol Evol 1997 Mar;14(3):354]

It has recently been suggested that observed levels of variation at microsatellite loci can be used to infer patterns of selection in genomes and to assess demographic history. In order to evaluate the feasibility of these suggestions it is necessary to know something about how levels of variation at microsatellite loci are expected to fluctuate due simply to stochasticity in the processes of mutation and inheritance (genetic sampling). Here we use recently derived properties of the stepwise mutation model to place confidence intervals around the variance in repeat score that is expected at mutation-drift equilibrium and outline a statistical test for whether an observed value differs significantly from expectation. We also develop confidence intervals for the time course of the buildup of variation following a complete elimination of variation, such as might be caused by a selective sweep or an extreme population bottleneck. We apply these methods to the variation observed at human Y-specific microsatellites. Although a number of authors have suggested the possibility of a very recent sweep, our analyses suggest that a sweep or extreme bottleneck is unlikely to have occurred anytime during the last approximately 74,000 years. To generate this result we use a recently estimated mutation rate for microsatellite loci of 5.6 x 10(-4) along with the variation observed at autosomal microsatellite loci to estimate the human effective population size. This estimate is 18,000, implying an effective number of 4,500 Y chromosomes. One important general conclusion to emerge from this study is that in order to reject mutation-drift equilibrium at a set of linked microsatellite loci it is necessary to have an unreasonably large number of loci unless the observed variance is far below that expected at mutation-drift equilibrium.      

Evidence for genetic hitchhiking effect associated with insecticide resistance in Aedes aegypti.

Information on genetic variation within and between populations is critical for understanding the evolutionary history of mosquito populations and disease epidemiology. Previous studies with Drosophila suggest that genetic variation of selectively neutral loci in a large fraction of genome may be constrained by fixation of advantageous mutations associated with hitchhiking effect. This study examined restriction fragment length polymorphisms of four natural Aedes aegypti mosquito populations from Trinidad and Tobago, at 16 loci. These populations have been subjected to organophosphate (OP) insecticide treatments for more than two decades, while dichlor-diphenyltrichlor (DDT) was the insecticide of choice prior to this period. We predicted that genes closely linked to the OP target loci would exhibit reduced genetic variation as a result of the hitchhiking effect associated with intensive OP insecticide selection. We also predicted that genetic variability of the genes conferring resistance to DDT and loci near the target site would be similar to other unlinked loci. As predicted, reduced genetic variation was found for loci in the general chromosomal region of a putative OP target site, and these loci generally exhibited larger F(ST) values than other random loci. In contrast, the gene conferring resistance to DDT and its linked loci show polymorphisms and genetic differentiation similar to other random loci. The reduced genetic variability and apparent gene deletion in some regions of chromosome 1 likely reflect the hitchhiking effect associated with OP insecticide selection.     

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.     

Enzyme polymorphism in the freshwater mussel species, Aspatharia pfeifferiana (Mollusca, Lamellibranchiata) from southern Africa

Gene products of 31 protein coding loci in Aspatharia pfeifferiana examined by horizontal starch gel electrophoresis revealed genetic variation at 10 (32.26%) of the loci studied. The mean number of alleles per locus was 1±39 (±0±11), and the average heterozygosity value was 0±074 (±0±026). Not only is this the first account of electrophoretic variants of this species, but also the results show that the species has less genetic variation than most other molluscs. Genetic variation compares favourably with values obtained for other invertebrate species and for fish from the same geographical region. A. pfeifferiana has sufficient variation to adapt to environmental changes or to be used in selection programmes. The results are useful for conservation stocking and it will undoubtedly aid in improving the systematics and phylogeny of freshwater clams in general.     

The Quantitative Genetics of Phenotypic Robustness

Phenotypic robustness, or canalization, has been extensively investigated both experimentally and theoretically. However, it remains unknown to what extent robustness varies between individuals, and whether factors buffering environmental variation also buffer genetic variation. Here we introduce a quantitative genetic approach to these issues, and apply this approach to data from three species. In mice, we find suggestive evidence that for hundreds of gene expression traits, robustness is polymorphic and can be genetically mapped to discrete genomic loci. Moreover, we find that the polymorphisms buffering genetic variation are distinct from those buffering environmental variation. In fact, these two classes have quite distinct mechanistic bases: environmental buffers of gene expression are predominantly sex-specific and trans-acting, whereas genetic buffers are not sex-specific and often cis-acting. Data from studies of morphological and life-history traits in plants and yeast support the distinction between polymorphisms buffering genetic and environmental variation, and further suggest that loci buffering different types of environmental variation do overlap with one another. These preliminary results suggest that naturally occurring polymorphisms affecting phenotypic robustness could be abundant, and that these polymorphisms may generally buffer either genetic or environmental variation, but not both.     

小麦体细胞无性系SSR位点的遗传变异特性分析

研究结果表明:(1)小麦体细胞无性系SSR位点变异类型有:扩增片段迁移率的变大或变小、扩增片段缺失以及新的扩增片段;(2)变异特点为:变异频率与基因型有关,不同染色体组上的SSR位点变异频率不同,而不同无性系后代的SSR位点变异频率也不同;(3)同一SSR位点的变异类型在同一基因型的无性系后代中变异表现一致,在不同基因型无性系后代中的变异表现不同,有的SSR位点在无性系后代中表现出一致的变异,而有的则不一致。     

PERSPECTIVE: HIGHLY VARIABLE LOCI AND THEIR INTERPRETATION IN EVOLUTION AND CONSERVATION

Although highly variable loci, such as microsatellite loci, are revolutionizing both evolutionary and conservation biology, data from these loci need to be carefully evaluated. First, because these loci often have very high within-population heterozygosity, the magnitude of differentiation measures may be quite small. For example, maximum G_ST values for populations with no common alleles at highly variable loci may be small and are at maximum less than the average within-population homozygosity. As a result, measures that are variation independent are recommended for highly variable loci. Second, bottlenecks or a reduction in population size can generate large genetic distances in a short time for these loci. In this case, the genetic distance may be corrected for low variation in a population and tests to detect bottlenecks are advised. Third, statistically significant differences may not reflect biologically meaningful differences both because the patterns of adaptive loci may not be correlated with highly variable loci and statistical power with these markers is so high. As an example of this latter effect, the statistical power to detect a one-generation bottleneck of different sizes for different numbers of highly variable loci is discussed. All of these concerns need to be incorporated in the utilization and interpretation of patterns of highly variable loci for both evolutionary and conservation biology.