The role of diet and temperature in shaping cranial diversification of South American human populations: an approach based on spatial regression and divergence rate tests

Aim Understanding the importance of ecological factors in the origin and maintenance of patterns of phenotypic variation among populations, in an explicit geographical context, is one of the main goals of human biology, ecology and evolutionary biology. Here we study the ecological factors responsible for craniofacial variation among human populations from South America. Location South America. Methods We studied a dataset of 718 males from 40 South American populations, coming from groups that inhabited different geographical and ecological regions. Cranial size and shape variation were studied using 30 cranial measurements. We first used spatial correlograms and interpolated maps to address spatial patterns. We then regressed the shape (principal component scores) and size variables against ecology (mean annual temperature and diet) using multiple and multivariate spatial regression. Finally, the expected magnitudes of shape and size divergence under the influence of genetic drift and mutations alone were evaluated using neutral expectation for the divergence rate. Results The spatial correlograms showed a cline affecting the entire South American distribution. Interpolated maps showed that size and allometric shape vary from south‐east to north‐west. Multiple and multivariate regression analyses suggested that diet has the largest and most significant effect on this pattern of size and allometric shape variation. Finally, the results of the divergence rate test suggested that random processes alone cannot account for the morphological divergence exhibited by cranial size and allometric shape scores among southernmost populations. Main conclusions Correlograms, spatial regression and divergence rate analyses showed that although local factors (neutral processes or local environmental conditions) are important to explain spatial interpopulation differentiation in cranial characteristics among these populations, there is significant correlation of cranial size and allometric shape variation with diet. Gene flow among human populations, or local environmental conditions, could explain spatial variation mainly at smaller spatial scales, whereas the large‐scale pattern of the South American dataset is mainly related to the high proportion of carbohydrates and low proportion of proteins consumed.     

Ecological and evolutionary factors in dental morphological diversification among modern human populations from southern South America

The knowledge of processes involved in morphological variation requires the integrated analysis of evolutionary and ecological factors. Here, we investigate the factors responsible for dental variation among human populations from southern South America. The aim of this work is to test the correspondence of dental size and shape variation with geographical, molecular (i.e. mtDNA) and ecological (i.e. climate, diet and food preparation) variables employing comparative phylogenetic methods, which have not previously been extensively applied at a within-species level. The results of the Procrustes analysis show a significant association of shape variables with molecular distance and geography, whereas dental size is not associated with molecular or geographical distances among groups. Phylogenetic generalized least-squares analysis, which takes into account the evolutionary autocorrelation among populations, shows a significant relationship between dental size variation and diet, while temperature and pottery do not correspond with dental size or shape. Specifically, groups with diets rich in carbohydrates, as well as the maritime hunter-gatherers, have the smallest teeth. In summary, our results support ecological factors as the dominant factor on dental size diversification in this region, while evolutionary relationships account for variation in dental shape.     

Effect of bite force and diet composition on craniofacial diversification of Southern South American human populations

Ecological factors can be important to shape the patterns of morphological variation among human populations. Particularly, diet plays a fundamental role in craniofacial variation due to both the effect of the nutritional status—mostly dependent on the type and amount of nutrients consumed—on skeletal growth and the localized effects of masticatory forces. We examine these two dimensions of diet and evaluate their influence on morphological diversification of human populations from southern South America during the late Holocene. Cranial morphology was measured as 3D coordinates defining the face, base and vault. Size, form, and shape variables were obtained for 474 adult individuals coming from 12 samples. Diet composition was inferred from carious lesions and δ¹³C data, whereas bite forces were estimated using traits of main jaw muscles. The spatial structure of the morphological and ecological variables was measured using correlograms. The influence of diet composition and bite force on morphometric variation was estimated by a spatial regression model. Cranial variation and diet composition display a geographical structure, while no geographical pattern was observed in bite forces. Cranial variation in size and form is significantly associated with diet composition, suggesting a strong effect of systemic factors on cranial growth. Conversely, bite forces do not contribute significantly to the pattern of morphological variation among the samples analyzed. Overall, these results show that an association between diet composition and hardness cannot be assumed, and highlight the complex relationship between morphological diversification and diet in human populations. Am J Phys Anthropol 155:114–127, 2014. © 2014 Wiley Periodicals, Inc.     

Morphological differentiation of aboriginal human populations from Tierra del Fuego (Patagonia): implications for South American peopling

This study aims to integrate the craniofacial morphological variation of southern South American populations with the results of mtDNA haplogroup variation, to discuss the South America peopling. Because the causes of morphological differentiation of Fueguian populations are still a controversial subject, the comparison with neutral variation could contribute to elucidate them. Samples of human remains from South America regions were used to analyze the evolutionary relationships. Several craniofacial traits observed in frontal and lateral view were analyzed by means of geometric morphometrics techniques, and the evolutionary relationships based on morphological and molecular data were established in base to ordination analyses. The results from the facial skeleton agree with those obtained from mtDNA haplogroup frequencies, with La Pampa/Chaco samples detached from the Patagonian samples. Hence, the same mechanism that accounts for the pattern of frequency of haplogroups could explain the variation found in facial skeleton among the samples. It is suggested that such geographic pattern of craniofacial and molecular diversity may reflect the effect of genetic drift that occurred in the small founding populations isolated by distance or geographic barriers. Conversely, the results obtained using the traits from the lateral view slightly differ from the molecular results, showing differences between southernmost Patagonian and the other samples. Therefore, mechanisms other than genetic drift (e.g., natural selection) could have acted to shape the pattern observed in some craniofacial structures present in the lateral view, characterized by the fact that the southernmost Patagonian samples display the most robust and dolichocephalic crania.     

Origins of clonal diversity in the hypervariable asexual ostracode Cypridopsis vidua

Prior allozyme studies have indicated that populations of the asexual ostracode, Cypridopsis vidua (Müller), show extraordinary clonal diversity. Based on a joint examination of allozyme variation and sequence divergence at the mitochondrial cytochrome c oxidase I (COI) gene, the present analysis provides new insights concerning the origins of this variation. The results establish that populations of C. vidua in one recently deglaciated region of North America are not only allozymically diverse, but also include several divergent mitochondrial DNA (mtDNA) lineages. The extent of sequence divergence among these lineages is so large as to suggest their diversification over the past 7–8 million years. The patterning of genetic divergence among co-occurring clones makes it apparent that much of the mtDNA and allozyme diversity in local populations owes its origins to recurrent colonization events. However, in situ mutational diversification also appears to explain some variation. The mechanisms enabling the sustained coexistence of such a large array of closely allied genotypes remain unclear, but there is an apparent difference in equilibrium diversity between benthic and planktonic asexual organisms.     

East-West cranial differentiation in pre-Columbian populations from Central and North America

In a recent study we found that crania from South Amerindian populations on each side of the Andes differ significantly in terms of craniofacial shape. Western populations formed one morphological group, distributed continuously over 14,000km from the Fuegian archipelago (southern Chile) to the Zulia region (northwestern Venezuela). Easterners formed another group, distributed from the Atlantic Coast up to the eastern foothills of the Andes. This differentiation is further supported by several genetic studies, and indirectly by ecological and archaeological studies. Some authors suggest that this dual biological pattern is consistent with differential rates of gene flow and genetic drift operating on both sides of the Cordillera due to historical reasons. Here we show that such East-West patterning is also observable in North America. We suggest that the "ecological zones model" proposed by Dixon, explaining the spread of the early Americans along a Pacific dispersal corridor, combined with the evolution of different population dynamics in both regions, is the most parsimonious mechanism to explain the observed patterns of within- and between-group craniofacial variability.     

Glacial refugia in a maritime temperate climate: Cicada (Kikihia subalpina) mtDNA phylogeography in New Zealand

Understanding the biological significance of Pleistocene glaciations requires knowledge of the nature and extent of habitat refugia during glacial maxima. An opportunity to examine evidence of glacial forest refugia in a maritime, Southern Hemisphere setting is found in New Zealand, where the extent of Pleistocene forests remains controversial. We used the mitochondrial phylogeography of a forest-edge cicada ( Kikihia subalpina ) to test the hypothesis that populations of this species survived throughout South Island during the Last Glacial Maximum. We also compared mitochondrial DNA phylogeographic patterns with male song patterns that suggest allopatric divergence across Cook Strait. Cytochrome oxidase I and II sequences were analyzed using network analysis, maximum-likelihood phylogenetic estimation, Bayesian dating and Bayesian skyline plots. K. subalpina haplotypes from North Island and South Island form monophyletic clades that are concordant with song patterns. Song divergence corresponds to approximately 2% genetic divergence, and Bayesian dating suggests that the North Island and South Island population-lineages became isolated around 761 000 years bp . Almost all South Island genetic variation is found in the north of the island, consistent with refugia in Marlborough Sounds, central Nelson and northwest Nelson. All central and southern South Island and Stewart Island haplotypes are extremely similar to northern South Island haplotypes, a 'northern richness/southern purity' pattern that mirrors genetic patterns observed in many Northern Hemisphere taxa. Proposed southern South Island forest habitat fragments may have been too small to sustain populations of K. subalpina , and/or they may have harboured ecological communities with no modern-day analogues.     

Discordant patterns of evolutionary differentiation in two Neotropical treefrogs

Comparative studies of codistributed taxa test the degree to which historical processes have shaped contemporary population structure. Discordant patterns of lineage divergence among taxa indicate that species differ in their response to common historical processes. The complex geologic landscape of the Isthmus of Central America provides an ideal setting to test the effects of vicariance and other biogeographic factors on population history. We compared divergence patterns between two codistributed Neotropical frogs ( Dendropsophus ebraccatus and Agalychnis callidryas ) that exhibit colour pattern polymorphisms among populations, and found significant differences between them in phenotypic and genetic divergence among populations. Colour pattern in D. ebraccatus did not vary with genetic or geographic distance, while colour pattern co-varied with patterns of gene flow in A. callidryas . In addition, we detected significant species differences in the phylogenetic history of populations, gene flow among them, and the extent to which historical diversification and recent gene flow have been restricted by five biogeographic barriers in Costa Rica and Panama. We inferred that alternate microevolutionary processes explain the unique patterns of diversification in each taxon. Our study underscores how differences in selective regimes and species-typical ecological and life-history traits maintain spatial patterns of diversification.     

Genetic divergence disclosing a rapid prehistorical dispersion of native americans in central and South america

An accurate estimate of the divergence time between Native Americans is important for understanding the initial entry and early dispersion of human beings in the New World. Current methods for estimating the genetic divergence time of populations could seriously depart from a linear relationship with the true divergence for multiple populations of a different population size and significant population expansion. Here, to address this problem, we propose a novel measure to estimate the genetic divergence time of populations. Computer simulation revealed that the new measure maintained an excellent linear correlation with the population divergence time in complicated multi-population scenarios with population expansion. Utilizing the new measure and microsatellite data of 21 Native American populations, we investigated the genetic divergences of the Native American populations. The results indicated that genetic divergences between North American populations are greater than that between Central and South American populations. None of the divergences, however, were large enough to constitute convincing evidence supporting the two-wave or multi-wave migration model for the initial entry of human beings into America. The genetic affinity of the Native American populations was further explored using Neighbor-Net and the genetic divergences suggested that these populations could be categorized into four genetic groups living in four different ecologic zones. The divergence of the population groups suggests that the early dispersion of human beings in America was a multi-step procedure. Further, the divergences suggest the rapid dispersion of Native Americans in Central and South Americas after a long standstill period in North America.     

Variable extent of parallelism in respiratory,circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non‐parallelism in gill size divergence related to differences in genetic background. There was also a non‐significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non‐parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.     

Biogeography of the túngara frog, Physalaemus pustulosus: a molecular perspective

Physalaemus pustulosus, a small leptodactylid frog with South American affinities, ranges across northern South America through Middle America to southern Mexico. To investigate its geographic variation and evolutionary origins, we analysed the presumptive gene products of 14 allozyme loci and sequenced a portion of the mitochondrial COI gene from individuals sampled throughout the distribution. Generally, allozyme dissimilarities and sequence divergences are correlated with each other and with geographic proximity. The greatest discontinuity in genetic variation was found between populations in Middle America vs. South America + Panama. Based on two Bayesian MCMC (Markov chain Monte Carlo) divergence time estimates involving two independent temporal constraints, the timing of the separation of northern and southern túngara frog lineages is significantly older than the time since completion of the current Panama land bridge. P. pustulosus first invaded Middle America from South America about 6-10 million years ago giving rise to the northern lineage. The southern lineage then invaded Panama independently after land bridge completion. Despite millions of years of independent evolution, the multilocus allozyme data revealed that western Panama populations represent a contact zone containing individuals with alleles from both groups present.     

Genetic diversification without obvious genitalic morphological divergence in harvestmen (Opiliones, Laniatores, Sclerobunus robustus) from montane sky islands of western North America

The southern Rocky Mountains and adjacent Intermontane Plateau Highlands region of western North America is a geographically diverse area with an active geologic history. Given the topological complexity and extensive geologic activity, organisms inhabiting this region are expected to show some degree of morphological and genetic divergence, especially populations found on the southern montane 'sky islands' of this region. Here we examine the phylogeographic history and diversification of a montane forest inhabiting harvestmen, Sclerobunus robustus, using a combination of genetic and morphological data. Divergence time estimates indicate that much of the diversification within and between major groups S. robustus predate the Pleistocene glacial cycles. The most widespread subspecies, Sclerobunus robustus robustus, is recovered as six genetically distinct, geographically cohesive mitochondrial phylogroups. Gene tree data for a single nuclear gene reveals congruent, albeit slightly more conservative, patterns of genetic divergence. Despite high levels of genetic divergence throughout their distribution, phylogroups show extreme conservation in somatic and reproductive morphology. This uncoupling of morphological and genetic differentiation may be due to morphological conservatism associated with a conserved microhabitat preference. Based on these data, it is obvious that S. robustus has undergone some level of cryptic diversification.     

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery‐reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996–2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make‐up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.     

中国板栗居群间等位酶基因频率的空间分布

中国板栗21个自然群间等位酶遗传变异的空间自相关分析及F－统计分析结果表明：其多数等位基因频率在居群间呈随机分布模式,缺乏一定的空间结构;而部分等位基因表现为渐变或双向渐变的非随机分布模式,又具有特定空间结构。中国板栗遗传变异空间结构模式的形成可能是长距离基因流、自然气候、人类活动、地理距离隔离等诸因素综合作用的结果。文中还就居群等位基因分布格局的成因进行了讨论：在第四纪冰川后,中国板栗以长江流域中下游的孑遗中心为起点,等位基因分别沿着向北和向南的不同方向迁移形成现在的居群结构;季风气候和人类活动干扰是削弱居群分化的主要因素,而基于环境梯度的选择,是形成由北向南渐变分布的原因。     

High levels of variation despite genetic fragmentation in populations of the endangered mountain pygmy-possum, Burramys parvus, in alpine Australia

In endangered mammals, levels of genetic variation are often low and this is accompanied by genetic divergence among populations. The mountain pygmy-possum (Burramys parvus) is an endangered marsupial restricted to the alpine region of Victoria and New South Wales, Australia. By scoring variation at eight microsatellite loci, we found that B. parvus populations exhibit high levels of genetic divergence and fall into three distinct groups from the northern, central and southern areas of the distribution of this species, consistent with previous assessments of mitochondrial DNA variation. F(ST) values between populations from these regions ranged from 0.19 to 0.54. Within the central area, there was further genetic fragmentation, and a linear association between genetic and geographical distance. This pattern is likely to reflect limited dispersal across barriers despite the fact that individual B. parvus can move several kilometres. Levels of genetic variation within populations were high with the exception of a southern population where there was evidence of inbreeding. From a conservation perspective, all three areas where B. parvus are found should be considered as separate gene pools; management of populations within these areas needs to take into account the low gene flow between populations, as well as threats posed by roads, resorts and other developments in the alpine region. The low genetic variability and inbreeding in the southern population is of particular concern given the high levels of variability in other B. parvus populations.     

Phylogeny and diversification of Valerianaceae (Dipsacales) in the southern Andes

The southern Andean clade of Valeriana provides an excellent model for the study of biogeography. Here we provide new data to help clarify phylogenetic relationships among the South American valerians, with special focus on taxa found in the southern Andes. We found that the southern Andean taxa formed a clade in maximum likelihood and maximum parsimony analyses, and used a Bayesian relaxed clock method to estimate divergence times within Valerianaceae. Our temporal results were similar to other studies, but we found greater variance in our estimates, suggesting that the species of Valeriana have been on the South American continent for some time, and have been successful at exploiting new niche opportunities that reflects the contemporary radiation. Regardless of the time frame for the radiation of the clade, the uptick in the rate of diversification in Valerianaceae appears correlated with a dispersal event from Central to South America. The appearance of Valeriana in the southern Andes (13.7 Ma) corresponds with the transition from closed forest on the western side of the Andes in central Chile to a more open Mediterranean woodland environment. This would suggest that the high species richness of Valerianaceae in South America is the result of multiple, smaller radiations such as the one in the southern Andes, that may or may not be geographically isolated. These smaller radiations may also be driven by species moving into new biomes (migration from a temperate to a more Mediterranean-type climate and into alpine). The degree to which different ecological and geological factors interact to drive diversification is difficult to ascertain. Likewise, without a better-resolved phylogeny it is impossible to determine the directionality of dispersal in this group; did they colonize the southern Andes first, then move northward as the central Andean alpine habitat became more widely available or vice versa?     

HISTORICAL BIOGEOGRAPHY AND CONTEMPORARY PATTERNS OF MITOCHONDRIAL DNA VARIATION IN WHITE-TAILED DEER FROM THE SOUTHEASTERN UNITED STATES

Mitochondrial DNA (mtDNA) was used to characterize patterns of geographic variation among white-tailed deer (Odocoileus virginianus) populations in the southeastern United States. Fifteen restriction enzymes were employed to survey and map 99 restriction sites in 142 deer from 18 localities in five southeastern states. Phylogenetic analysis revealed three primary groups of haplotypes: (1) southern Florida and the Florida Keys, (2) the remainder of peninsular Florida northward to South Carolina, and (3) the Florida panhandle westward to Mississippi. Geographical heterogeneity in haplotype frequencies suggests that stochastic lineage sorting or isolation by distance are not important determinates of mtDNA differentiation among deer populations. The pattern of mtDNA variation in white-tailed deer is concordant spatially with those observed in unrelated taxa suggesting the common influence of historical biogeographic events. The data (1) support previous hypotheses that relate contemporary patterns of intraspecific phylogeography in northern Florida to the physiogeographic history of the region; and (2) suggest that genetic differentiation in southern Florida may be attributable to episodes of Pleistocene dispersal. Despite potentially high vagility and human intervention, ecological and demographic characteristics of deer have effectively preserved the historical pattern of intraspecific mtDNA differentiation.     

Speciation along a latitudinal gradient: The origin of the Neotropical cycad sister pair Dioon sonorense–D. vovidesii (Zamiaceae)

Latitude is correlated with environmental components that determine the distribution of biodiversity. In combination with geographic factors, latitude‐associated environmental variables are expected to influence speciation, but empirical evidence on how those factors interplay is scarce. We evaluated the genetic and environmental variation among populations in the pair of sister species Dioon sonorense–D. vovidesii, two cycads distributed along a latitudinal environmental gradient in northwestern Mexico, to reveal their demographic histories and the environmental factors involved in their divergence. Using genome‐wide loci data, we determined the species delimitation, estimated the gene flow, and compared multiple demographic scenarios of divergence. Also, we estimated the variation of climatic variables among populations and used ecological niche models to test niche overlap between species. The effect of geographic and environmental variables on the genetic variation among populations was evaluated using linear models. Our results suggest the existence of a widespread ancestral population that split into the two species ~829 ky ago. The geographic delimitation along the environmental gradient occurs in the absence of major geographic barriers, near the 28th parallel north, where a zonation of environmental seasonality exists. The northern species, D. vovidesii, occurs in more seasonal environments but retains the same niche of the southern species, D. sonorense. The genetic variation throughout populations cannot be solely explained by stochastic processes; the latitudinal‐associated seasonality has been an additive factor that strengthened the species divergence. This study represents an example of how speciation can be achieved by the effect of the latitude‐associated factors on the genetic divergence among populations.     

Population genetics of Eldana saccharina Walker (Lepidoptera: Pyralidae) and the implications for management using biocontrol

The African sugarcane stalk borer, Eldona saccharina Walker (Lepidoptera: Pyralidae), is widely distributed throughout sub-Saharan Africa and is an important insect pest of maize and sugarcane. The insect shows significant variation in behaviour, host plant and natural enemy guild in different regions. Several attempts to redistribute the natural enemies of E. saccharina from West Africa to South Africa were unsuccessful. The significant behavioural, host plant and natural enemy variations as well as failures of biocontrol attempts evoked a hypothesis of genetic diversification. To evaluate this hypothesis a molecular analysis was conducted on geographically isolated populations of E. saccharina from East, North, South and West Africa, using the cytochrome c oxidase subunit I (COI) region of the mitochondrial genome. The results revealed that E. saccharina populations are separated into four major units corresponding to the West Africa, Rift Valley, South/East Africa and southern African populations. Mitochondrial DNA divergence among the four populations ranged from 1% to 4.98%. To examine the impact of the observed genetic variation on the fertility of inter-population crosses, a mating experiment was conducted between the Rift valley and South African population to produce an F1 generation, and these were backcrossed with the South African parent population. Fertility of eggs produced by the F1/parent population cross was significantly reduced when compared to fertility of the "true" South African line, and the F1/F1 cross. The contributions of the observed genetic differences and inter population incompatibility for the failure of previous biocontrol attempts are discussed and recommendations on future biocontrol practices are given.     

WEEDY ADAPTATION IN SETARIA SPP. I. ISOZYME ANALYSIS OF GENETIC DIVERSITY AND POPULATION GENETIC STRUCTURE IN SETARIA VIRIDIS

Setaria viridis is an important self-pollinating, cosmopolitan weed of temperate regions worldwide. Allozyme markers were used to investigate genetic diversity and structure in 168 accessions (including four S. italica) collected mainly from North America and Eurasia. Genetic diversity in green foxtail, and its population genetic structure, provided important clues about this weed's evolutionary history. Genetic diversity was low, with marked population differentiation: the percentage of polymorphic loci was 25% (0.95 criterion); mean number of alleles per locus was 1.86; mean panmictic heterozygosity was 0.07; and the coefficient of population genetic differentiation was 0.65. A common genotype occurred in 25 accessions distributed in six countries from both the Old World and New World, in a wide variety of ecological situations. Relatively little genetic divergence occurred between Eurasia and North America, with Nei's unbiased genetic identity between the two regions equaling 1.0. Populations from these two continents also had equivalent genetic diversity. Within North America, regional differentiation was indicated by northern and southern groups separated at 43.5° N latitude. No geographic pattern in genetic diversity was found within Eurasia. The size of the geographic range from which populations were sampled was not an accurate indicator of the extent of genetic diversity found among populations from that region. These results suggest that present patterning among green foxtail populations in North America is the consequence of multiple introductions into the New World followed by local adaptation and regional differentiation. Finally, S. italica and several green foxtail varieties did not differ isozymatically from typical forms of green foxtail. This supports the view that S. italica and S. viridis are conspecific, that the former (foxtail millet) is a domesticated form of the latter, and also questions the taxonomic validity of formally recognizing morphological varieties within green foxtail.