The Genetic Structure of a Tribal Population, the Yanomama Indians. Xv. Patterns Inferred by Autocorrelation Analysis

Fifteen allele frequencies have previously been determined for 50 villages of the Yanomama, an Amerindian tribe from southern Venezuela and northern Brazil. These frequencies were subjected to spatial autocorrelation analysis to investigate their population structure. There are significant spatial patterns for most allele frequencies. Clinical patterns, investigated by one-dimensional and directional spatial correlograms, were relatively few in number and were moderate in strength. Overall, however, there is a marked decline in genetic similarity with geographic distance. The results are compatible with a hierarchic population structure superimposed on the geography, and generated by a stochastic fission-fusion model of village propagation, followed by localized gene flow. Strong temporal autocorrelations of allele frequencies based on linguistic-historical distances representing time since divergence were also found. There appears to be a stronger relation between geography and linguistic-historical hierarchic subdivisions than between either feature and genetic distances. These findings confirm by different approaches the results of earlier analyses concerning the important roles of both stochastic and social factors in determining village allele frequencies and the occurrence within this tribe of some allele frequency clines most likely due to the operation of chance historical processes.     

A powerful regression-based method for admixture mapping of isolation across the genome of hybrids

We propose a novel method that uses natural admixture between divergent lineages (hybridization) to investigate the genetic architecture of reproductive isolation and adaptive introgression. Our method employs multinomial regression to estimate genomic clines and to quantify introgression for individual loci relative to the genomic background (clines in genotype frequency along a genomic admixture gradient). Loci with patterns of introgression that deviate significantly from null expectations based on the remainder of the genome are potentially subject to selection and thus of interest to understanding adaptation and the evolution of reproductive isolation. Using simulations, we show that different forms of selection modify these genomic clines in predictable ways and that our method has good power to detect moderate to strong selection for multiple forms of selection. Using individual-based simulations, we demonstrate that our method generally has a low false positive rate, except when genetic drift is particularly pronounced (e.g. low population size, low migration rates from parental populations, and substantial time since initial admixture). Additional individual-based simulations reveal that moderate selection against heterozygotes can be detected as much as 50 c m away from the focal locus directly experiencing selection, but is not detected at unlinked loci. Finally, we apply our analytical method to previously published data sets from a mouse ( Mus musculus and M. domesticus ) and two sunflower ( Helianthus petiolaris and H. annuus ) hybrid zones. This method should be applicable to numerous species that are currently the focus of research in evolution and ecology and should help bring about new insights regarding the processes underlying the origin and maintenance of biological diversity.     

The effect of the Neolithic expansion on European molecular diversity

We performed extensive and realistic simulations of the colonization process of Europe by Neolithic farmers, as well as their potential admixture and competition with local Palaeolithic hunter-gatherers. We find that minute amounts of gene flow between Palaeolithic and Neolithic populations should lead to a massive Palaeolithic contribution to the current gene pool of Europeans. This large Palaeolithic contribution is not expected under the demic diffusion (DD) model, which postulates that agriculture diffused over Europe by a massive migration of individuals from the Near East. However, genetic evidence in favour of this model mainly consisted in the observation of allele frequency clines over Europe, which are shown here to be equally probable under a pure DD or a pure acculturation model. The examination of the consequence of range expansions on single nucleotide polymorphism (SNP) diversity reveals that an ascertainment bias consisting of selecting SNPs with high frequencies will promote the observation of genetic clines (which are not expected for random SNPs) and will lead to multimodal mismatch distributions. We conclude that the different patterns of molecular diversity observed for Y chromosome and mitochondrial DNA can be at least partly owing to an ascertainment bias when selecting Y chromosome SNPs for studying European populations.     

THE COLONIZATION OF DROSOPHILA SUBOBSCURA IN CHILE. II. CLINES IN THE CHROMOSOMAL ARRANGEMENTS

Drosophila subobscura is a Palearctic species that was first detected in the New World in Puerto Montt (Chile) in February 1978. Since that time, it has spread over a broad area and increased in population density. The South American populations exhibit a high level of chromosomal polymorphism: 20 different arrangements exist, distributed among five chromosomes. Chromosomal arrangement heterozygosity varies from 0.55 to 0.61 in the nine populations examined. Incipient clines in the frequencies of the arrangements are appearing; these clines follow the same latitudinal direction as in the Old World. Wing length significantly decreases with latitude, as it does in Europe. The colonization of South America by D. subobscura appears to be a major natural experiment with outcomes that duplicate the distributional patterns—in chromosomal polymorphism and in wing length—observed in the Old World, thereby strongly supporting the adaptive significance of these patterns. The data show a very rapid effect of natural selection promoting genetic differentiation among natural populations.     

Secondary Contact and Admixture between Independently Invading Populations of the Western Corn Rootworm,Diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a Bayesian analysis of the population structure and in an approximate Bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed.     

Is floral divergence sufficient to maintain species boundaries upon secondary contact in Mediterranean food-deceptive orchids?

Analyzing the processes that determine whether species boundaries are maintained on secondary contact may shed light on the early phase of speciation. In Anacamptis morio and Anacamptis longicornu, two Mediterranean orchid sister-species, we used molecular and morphological analyses, together with estimates of pollination success and experimental crosses, to assess whether floral isolation can shelter the species' genomes from genetic admixture on secondary contact. We found substantial genetic and morphological homogenization in sympatric populations in combination with an apparent lack of postmating isolation. We further detected asymmetric introgression in the sympatric populations and an imbalance in cytotype representation, which may be due either to a difference in flowering phenology or else be a consequence of cytonuclear incompatibilities. Estimates of genetic clines for markers across sympatric zones revealed markers that significantly deviated from neutral expectations. We observed a significant correlation between spur length and reproductive success in sympatric populations, which may suggest that directional selection is the main cause of morphological differentiation in this species pair. Our results suggest that allopatric divergence has not led to the evolution of sufficient reproductive isolation to prevent genomic admixture on secondary contact in this orchid species pair.     

Adaptive potential of northernmost tree populations to climate change, with emphasis on Scots pine (Pinus sylvestris L.)

The adaptive potential of the northernmost Pinus sylvestris L. (and other northern tree) populations is considered by examining first the current patterns of quantitative genetic adaptive traits, which show high population differentiation and clines. We then consider the postglacial history of the populations using both paleobiological and genetic data. The current patterns of diversity at nuclear genes suggest that the traces of admixture are mostly visible in mitochondrial DNA variation patterns. There is little evidence of increased diversity due to admixture between an eastern and western colonization lineage, but no signal of reduced diversity (due to sequential bottlenecks) either. Quantitative trait variation in the north is not associated with the colonizing lineages. The current clines arose rapidly and may be based on standing genetic variation. The initial phenotypic response of Scots pine in the north is predicted to be increased survival and growth. The genetic responses are examined based on quantitative genetic predictions of sustained selection response and compared with earlier simulation results that have aimed at more ecological realism. The phenotypic responses of increased growth and survival reduce the opportunity for selection and delay the evolutionary responses. The lengthening of the thermal growing period also causes selection on the critical photoperiod in the different populations. Future studies should aim at including multiple ecological and genetic factors in evaluating potential responses.     

Different cell size and cell number contribution in two newly established and one ancient body size cline of Drosophila subobscura

Latitudinal genetic clines in body size occur in many ectotherms including Drosophila species. In the wing of D. melanogaster, these clines are generally based on latitudinal variation in cell number. In contrast, differences in wing area that evolve by thermal selection in the laboratory are in general based on cell size. To investigate possible reasons for the different cellular bases of these two types of evolutionary response, we compared the newly established North and South American wing size clines of Drosophila subobscura. The new clines are based on latitudinal variation in cell area in North America and cell number in South America. The ancestral European cline is also based on latitudinal variation in cell number. The difference in the cellular basis of wing size variation in the American clines, which are roughly the same age, together with the similar cellular basis of the new South American cline and the ancient European one, suggest that the antiquity of a cline does not explain its cellular basis. Furthermore, the results indicate that wing size as a whole, rather than its cellular basis, is under selection. The different cellular bases of different size clines are most likely explained either entirely by chance or by different patterns of genetic variance--or its expression--in founding populations.     

Catawba Indians: morphology, genetics, and history

The Catawba Indians of South Carolina appear to be descendants of the once flourishing Siouan tribe now admixed primarily with Whites of the area. Study of physical traits and blood types was made in 1962 at the time they terminated reservation status. The data permit comparisions with other populations and an estimation of the genetic contribution of presumed parental populations to the hybrid. Skeletal remains of Siouan Indians from archeological sites in the area provide a basis for morphological comparisons. Similar ties of the Catawba to the Indian Knoll people of Kentucky and to the Occaneechi people of early historic times are noted. Results of blood typings of these ancient bones also enable some tentative serological comparisons with the present day Catawba Indians. Certain unusual factors in the history of the present population have affected their isolation and admixture, and this interaction of culture and biology is explored.     

Multivariate Analysis of Gametic Disequilibrium in the Yanomama

The gametic disequilibria between all possible pairs of loci were examined for a set of eight codominant loci in each of fifty Yanomama villages, using a multivariate correlation analysis which reduces the results to a single measure of departure from multiple-locus-gametic equilibrium. Thirty-two of the fifty villages departed significantly from multiple-locus gametic equilibrium. The largest contributions to the departure from multiple-locus equilibrium were due to the disequilibria between MN and Ss and between Rh(Cc) and Rh(Ee), indicating the effects of tight linkage. After removing the effects of these obvious sources of disequilibrium, sixteen of the fifty villages still remained significantly out of equilibrium. The disequilibrium between any particular pair of loci was highly erratic from village to village, and (with the exception of the MN-Ss and Cc-Ee disequilibria) averaged out very close to zero overall, suggesting a lack of systematic forces (epistatic selection). The departure from equilibrium in any one village is in excess of that expected from random sampling alone, and is attributed primarily to the fission-fusion mode of village formation operative in the Yanomama and the fact that a single village consists of a few extended lineages. Village allele frequencies are highly correlated across loci, and most of the non-independence is accounted for by large correlations in the average allelic frequencies of different loci for related villages. It is suggested that these correlations also are due to territorial expansion and population growth. For the tribe as a whole, all but the tightly linked markers of the MNSs and Rh complexes are approximately uncorrelated, and large departures from multiple-locus Hardy-Weinberg expectation are primarily due to substantial Wahlund variance within the tribe. There is no need to postulate a role for selection in these disequilibria.     

An experimental study of evolution in progress: Clines for quantitative traits in colonizing and Palearctic populations ofDrosophila

Summary Drosophila subobscura has recently colonized the American continent and is an excellent model for studying evolution in action. Previous analyses have shown that these colonizing populations have significant latitudinal clines for the frequencies of some chromosomal arrangements that parallel those clines found in the Old World. These results strongly suggest that this polymorphism is adaptive. In the present study, significant latitudinal clines for continuous morphometric variables (flies are larger in the north) have been detected in Old World populations ofD. subobscura. The adaptive nature of these clines is reinforced by the fact that parallel latitudinal clines for body size have also been detected inDrosophila obscura, a closely related sympatric species, as well as previously in otherDrosophila. On the other hand, no significant latitudinal clines for continuous morphometric traits, not even when using an overall size index, have been detected in colonizing populations ofD. subobscura. This is a rather surprising result given the number of generations that have elapsed since the species was detected in America and given that significant clines in chromosomal inversions are already established. Thus, the adaptive response of quantitative morphometric variables is not as rapid as that found for chromosomal inversions. Nevertheless, canonical correlation analysis suggests that significant latitudinal clines for body size might soon be detected in the American continent. The results obtained here are included in a projected time series with the aim of documenting size evolution in action.     

Biochemical genetics of Fundulus heteroclitus (L.). I. Temporal and spatial variation in gene frequencies of Ldh-B,Mdh-A,Gpi-B,and Pgm-A

Natural populations of Fundulus heteroclitus show extensive spatial variation in gene frequencies at four unlinked polymorphic loci. Large clinal changes in gene frequencies were found for Ldh-B, Mdh-A, and Gpi-B, whereas the spatial variation for the Pgm-B locus was small. Since the geographical area over which these clines are found is characterized by a steep thermal gradient, the clines in gene frequency are correlated with a directional change in mean water temperature. Maximum gene diversity of these four loci was correlated with annual fluctuations in water temperature. Temporal stability of the allelic frequencies was established for a 2–4 year period.This work was supported by NSF Grants GB37548 and DEB 76-19877 and by a grant from the National Geographic Society. A. R. P. is an NIH trainee supported by a training grant (No. HD00139) to the Department of Biology.This is contribution No. 960 of the Department of Biology, Johns Hopkins University.     

Population biology of the crustose red alga Lithophyllum incrustans Phil. 3. The effects of local environmental variables

Populations of the crustose red alga Lithophyllum incrustans were analysed for their age structure and reproductive characteristics. Two different types of age structure were observed, and shown to be associated with different surface area to volume ratios of the rock pools from which the populations were sampled. Geographic clines and surface area to volume ratios of the pools influenced the reproductive effort and the proportion of sexual plants in the populations. The relationship between temperature, salinity, surface area to volume ratio of the pools and the life history of the alga is discussed.     

Rapid evolution and range expansion of an invasive plant are driven by provenance–environment interactions

To improve our ability to prevent and manage biological invasions, we must understand their ecological and evolutionary drivers. We are often able to explain invasions after they happen, but our predictive ability is limited. Here, we show that range expansions of introduced Pinus taeda result from an interaction between genetic provenance and climate and that temperature and precipitation clines predict the invasive performance of particular provenances. Furthermore, we show that genotypes can occupy climate niche spaces different from those observed in their native ranges and, at least in our case, that admixture is not a main driver of invasion. Genotypes respond to climate in distinct ways, and these interactions affect the ability of populations to expand their ranges. While rapid evolution in introduced ranges is a mechanism at later stages of the invasion process, the introduction of adapted genotypes is a key driver of naturalisation of populations of introduced species.     

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well‐known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well‐known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Bo?i?evi? et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.     

Gene flow across tribal barriers and its effect among the Amazonian Içana river Indians

Demographic information was obtained from 622 individuals of five communities of primarily Baniwa Amerindians living near the Içana river in Brazil. Four of these populations, plus another from the same area, were also studied genetically. The latter investigation included the blood and, in some cases, saliva of 531 subjects, variously tested in relation to 40 genetic systems. Demographically these groups are characterized by young age, high intertribal admixture, low non-Indian admixture, high exogamy but low marital distance and high inbreeding, high fertility but low variance in offspring number, and relatively low mortality. Their gene pool shows a peptidase B variant (PEPB^2BAN1) and “private” polymorphism of carbonic anhydrase₂ (CA₂^BAN1) until now observed only among them. Other distinctive characteristics are the low frequencies of L^NS (0.08), L^Ns (0.09), R^z (0.01), R^O or r (0.02), ACP^A (0.08), GALT^D (0.01), and the relatively high prevalences of Gm^1;11,13,16 (0.05) and Gc¹ (0.82). Tf^Dehi occurs with a low prevalence (0.01). Genetic distance analysis reveals that the one Baniwa sample by history comprised of minimally admixed individuals is quite similar genetically to the Wapishana, another Arawak-speaking tribe some 900 km to the east, and that the genetic distances between the Baniwa communities reflect the amount of historical admixture in a way that indicates which should be excluded from considerations of intertribal genetic distances. Finally, the genetic relation of the Baniwa to the nearby tribes is examined.     

Contrasting plant physiological adaptation to climate in the native and introduced range of Hypericum perforatum

How introduced plants, which may be locally adapted to specific climatic conditions in their native range, cope with the new abiotic conditions that they encounter as exotics is not well understood. In particular, it is unclear what role plasticity versus adaptive evolution plays in enabling exotics to persist under new environmental circumstances in the introduced range. We determined the extent to which native and introduced populations of St. John's Wort (Hypericum perforatum) are genetically differentiated with respect to leaf-level morphological and physiological traits that allow plants to tolerate different climatic conditions. In common gardens in Washington and Spain, and in a greenhouse, we examined clinal variation in percent leaf nitrogen and carbon, leaf delta(13)C values (as an integrative measure of water use efficiency), specific leaf area (SLA), root and shoot biomass, root/shoot ratio, total leaf area, and leaf area ratio (LAR). As well, we determined whether native European H. perforatum experienced directional selection on leaf-level traits in the introduced range and we compared, across gardens, levels of plasticity in these traits. In field gardens in both Washington and Spain, native populations formed latitudinal clines in percent leaf N. In the greenhouse, native populations formed latitudinal clines in root and shoot biomass and total leaf area, and in the Washington garden only, native populations also exhibited latitudinal clines in percent leaf C and leaf delta(13)C. Traits that failed to show consistent latitudinal clines instead exhibited significant phenotypic plasticity. Introduced St. John's Wort populations also formed significant or marginally significant latitudinal clines in percent leaf N in Washington and Spain, percent leaf C in Washington, and in root biomass and total leaf area in the greenhouse. In the Washington common garden, there was strong directional selection among European populations for higher percent leaf N and leaf delta(13)C, but no selection on any other measured trait. The presence of convergent, genetically based latitudinal clines between native and introduced H. perforatum, together with previously published molecular data, suggest that native and exotic genotypes have independently adapted to a broad-scale variation in climate that varies with latitude.     

Chromosome Inversion Polymorphisms in DROSOPHILA MELANOGASTER. I. Latitudinal Clines and Associations between Inversions in Australasian Populations

Nineteen Australasian populations of Drosophila melanogaster have been screened for chromosome inversion polymorphisms. All 15 of the inversion types found are paracentric and autosomal, but only four of these, one on each of the major autosome arms, are common and cosmopolitan. North-south clines occur, with the frequencies of all four of the common cosmopolitan inversions increasing toward the equator. These clines in the Southern Hemisphere mirror north-south clines in the Northern Hemisphere, where the frequencies of all four of the common cosmopolitan inversions again increase towards the equator.-While few of the Australasian populations show significant disequilibrium between linked common cosmopolitan inversions, those that do invariably have excesses of coupling gametes, which is consistent with other reports. We also find nonrandom associations between the two major autosomes, with the northern populations in Australasia (those with high inversion frequencies) tending to be deficient in gametes with common cosmopolitan inversions on both major autosomes, while the southern populations in Australasia (low inversion frequencies) tend to have an excess of this class of gametes.-The clines and the nonrandom associations between the two major autosomes are best interpreted in terms of selection operating to maintain the common cosmopolitan inversion polymorphisms in natural populations of D. melanogaster.