Use of morphological traits and microsatellite markers to characterize the Tunisian cultivated and wild figs (Ficus carica L.)

We used 8 morphological traits and 17 simple sequence repeats loci to characterize 71 cultivated and wild Tunisian fig trees (Ficus carica L.). Significant morphological differences were inferred from leaf traits. The statistical analysis showed two major fig groups that indicated a common morphological basis. A total of 74 SSR alleles was revealed, defining 63 unique multilocus genotypes indicating a substantial genetic diversity. Based on multilocus SSR genotypes an identification key was established using MFC30, MFC3, MFC11 and MFC19 loci to identify figs. Analysis of variance components of linkage disequilibrium shown that among 136 pairs of loci, 32 present a significant gametic disequilibrium. The parameter D′²_IS (0.1284) was greater than D′²_ST (0.0079), a pointer of close to zero variance in total simple, and consequently the more pronounced independence of the 17 SSR loci. The majority of Ohta's variance components of linkage disequilibrium followed a pattern caused by genetic drift or a non-systematic disequilibrium profiles and natural selection occurs only for LMFC24-MFC8 pair loci in cultivated figs. Our results suggest that the morphological and SSR markers are suitable to characterize figs and should be recommended in conservation management strategy.     

Analysis of multilocus zygotic associations

While nonrandom associations between zygotes at different loci (zygotic associations) frequently occur in Hardy-Weinberg disequilibrium populations, statistical analysis of such associations has received little attention. In this article, we describe the joint distributions of zygotes at multiple loci, which are completely characterized by heterozygosities at individual loci and various multilocus zygotic associations. These zygotic associations are defined in the same fashion as the usual multilocus linkage (gametic) disequilibria on the basis of gametic and allelic frequencies. The estimation and test procedures are described with details being given for three loci. The sampling properties of the estimates are examined through Monte Carlo simulation. The estimates of three-locus associations are not free of bias due to the presence of two-locus associations and vice versa. The power of detecting the zygotic associations is small unless different loci are strongly associated and/or sample sizes are large (>100). The analysis of zygotic associations not only offers an effective means of packaging numerous genic disequilibria required for a complete characterization of multilocus structure, but also provides opportunities for making inference about evolutionary and demographic processes through a comparative assessment of zygotic association vs. gametic disequilibrium for the same set of loci in nonequilibrium populations.     

Comparative assessment of genetic diversity in wild and primitive cultivated barley in a center of diversity

Wild barley (Hordeum spontaneum K.) and indigenous primitive varieties of cultivated barley (Hordeum vulgare L.), collected from 43 locations in four eastern Mediterranean countries, Jordan, Syria, Turkey and Greece, were electrophoretically assayed for genetic diversity at 16 isozyme loci. Contrary to a common impression, cultivated barley populations were found to maintain a level of diversity similar to that in its wild progenitor species. Apportionment of overall diversity in the region showed that in cultivated barley within-populations diversity was of higher magnitude than the between-populations component. Neighboring populations of wild and cultivated barleys showed high degree of genetic identity. Groups of 3 or 4 isozyme loci were analyzed to detect associations among loci. Multilocus associations of varying order were detected for all three groups chosen for the analysis. Some of the association terms differed between the two species in the region. Although there was no clear evidence for decrease in diversity attributable to the domestication of barley in the region, there was an indication of different multilocus organizations in the two closely related species.     

Worldwide pattern of multilocus structure in barley determined by discrete log-linear multivariate analyses

Summary Data from the electrophoretic assay for seven enzyme loci of 1,032 accessions of cultivated barley, Hordeum vulgare L., from the USDA world barley collection were analyzed for multilocus structure using discrete log-linear multivariate techniques. Three major steps were involved in the analysis: (i) identification and elimination of terms that have inconsequential effects in multilocus association; (ii) construction of a log-linear model that best describes the complete multilocus structure of the genetic system; and (iii) evaluation of each of the association terms included in the model. The results of analyses of two subsets of loci show that the multilocus genetic system of cultivated barley, including loci located on different chromosomes, is organized into hierarchically structured complexes of loci. Multilocus structure differs in various geographical regions of the world. The structure of barleys from Southwest Asia, the putative center of origin for cultivated barley, is intermediate for both subsets of loci. Differences increased progressively across the Eurasian-African landmasses in each direction with increasing distance from Southwest Asia, with the consequence that the barleys from West Europe, East Asia, and Ethiopia are maximally different from those of Southwest Asia and Middle South Asia.     

The effects of multilocus balancing selection on neutral variability

We studied the effect of multilocus balancing selection on neutral nucleotide variability at linked sites by simulating a model where diallelic polymorphisms are maintained at an arbitrary number of selected loci by means of symmetric overdominance. Different combinations of alleles define different genetic backgrounds that subdivide the population and strongly affect variability. Several multilocus fitness regimes with different degrees of epistasis and gametic disequilibrium are allowed. Analytical results based on a multilocus extension of the structured coalescent predict that the expected linked neutral diversity increases exponentially with the number of selected loci and can become extremely large. Our simulation results show that although variability increases with the number of genetic backgrounds that are maintained in the population, it is reduced by random fluctuations in the frequencies of those backgrounds and does not reach high levels even in very large populations. We also show that previous results on balancing selection in single-locus systems do not extend to the multilocus scenario in a straightforward way. Different patterns of linkage disequilibrium and of the frequency spectrum of neutral mutations are expected under different degrees of epistasis. Interestingly, the power to detect balancing selection using deviations from a neutral distribution of allele frequencies seems to be diminished under the fitness regime that leads to the largest increase of variability over the neutral case. This and other results are discussed in the light of data from the Mhc.     

Genome and population dynamics under selection and neutrality: an example of S-allele diversity in wild cherry (Prunus avium L.)

Self-incompatibility, a common attribute of plant development, forms a classical paradigm of balancing selection in natural populations, in particular negative frequency-dependent selection. Under negative frequency-dependent selection population genetics theory predicts that the S-locus, being in command of self-incompatibility, keeps numerous alleles in equal frequencies demonstrating a wide allelic range. Moreover, while natural populations exhibit a higher within population genetic diversity, a reduction of population differentiation and increase of effective migration rate is expected in comparison to neutral loci. Allelic frequencies were investigated in terms of distribution and genetic structure at the gametophytic self-incompatibility locus in five wild cherry (Prunus avium L.) populations. Comparisons were also made between the differentiation at the S-locus and at the SSR loci. Theoretical expectations under balancing selection were congruent to the results observed. The S-locus showed broad multiplicity (16 S-alleles), high genetic diversity, and allelic isoplethy. Genetic structure at the self-incompatibility locus was almost four times lower than at 11 nSSR loci. Analysis of molecular variance revealed that only 5?% of the total genetic variation concerns differentiation among populations. In conclusion, the wealth of S-allelic diversity found in natural wild cherry populations in Greece is useful not only in advancing basic population genetics research of self-incompatibility systems in wild cherry but also in the development of breeding programs.     

Genome-Wide Association Analysis of Aluminum Tolerance in Cultivated and Tibetan Wild Barley

Tibetan wild barley (Hordeum vulgare L. ssp. spontaneum), originated and grown in harsh enviroment in Tibet, is well-known for its rich germpalsm with high tolerance to abiotic stresses. However, the genetic variation and genes involved in Al tolerance are not totally known for the wild barley. In this study, a genome-wide association analysis (GWAS) was performed by using four root parameters related with Al tolerance and 469 DArT markers on 7 chromosomes within or across 110 Tibetan wild accessions and 56 cultivated cultivars. Population structure and cluster analysis revealed that a wide genetic diversity was present in Tibetan wild barley. Linkage disequilibrium (LD) decayed more rapidly in Tibetan wild barley (9.30 cM) than cultivated barley (11.52 cM), indicating that GWAS may provide higher resolution in the Tibetan group. Two novel Tibetan group-specific loci, bpb-9458 and bpb-8524 were identified, which were associated with relative longest root growth (RLRG), located at 2H and 7H on barely genome, and could explain 12.9% and 9.7% of the phenotypic variation, respectively. Moreover, a common locus bpb-6949, localized 0.8 cM away from a candidate gene HvMATE, was detected in both wild and cultivated barleys, and showed significant association with total root growth (TRG). The present study highlights that Tibetan wild barley could provide elite germplasm novel genes for barley Al-tolerant improvement.     

Dynamics of Correlated Genetic Systems. V. Rates of Decay of Linkage Disequilibria in Experimental Populations of DROSOPHILA MELANOGASTER

The dynamic behavior of four-locus gametic frequency distributions was studied in five replicate cage populations of Drosophila melanogaster for up to 50 generations. The joint frequency distributions were resolved into gene frequencies and various disequilibrium measures. In addition, F statistics for marginal single-locus genotypic frequency distributions were followed through time. The gene frequency, disequilibrium and F statistics were obtained for four chromosome 3 enzyme marker loci [isocitrate dehydrogenase (3–27.1), esterase-6 (3–36.8), phosphoglucomutase (3–43.4) and esterase-C (3–49.0)]. The initial structure of the experimental populations featured random mating proportions, and two complementary gametic types with respect to the marker loci, thus assuring complete pairwise linkage disequilibrium among the markers.——The experimental results indicate: (1) the between-replicate variance in gene frequency varied substantially among loci, with isocitrate dehydrogenase showing the greatest between-replicate variance, and esterase-C the least. (2) The F statistics initially were strongly negative but decayed to the neighborhood of zero for all marker loci except esterase-C. The rate at which the F statistics approached zero varied among the marker loci, indicating substantial differences in the distribution of selective effects along the chromosome. The centromeric region, marked by esterase-C, shows the strongest selective effects. (3) The rate of decay of linkage disequilibrium was much faster than expected for pairs of neutral loci, averaging 1.82 times the neutral rate over all replicates and pairs of loci. This acceleration, which was observed for all six pairwise combinations of loci, was interpreted as resulting from the interaction between selection and recombination. Our experimental results are consistent with many investigations of linkage disequilibrium in natural populations of Drosophila melanogaster that show little or no disequilibrium among enzyme loci. (4) A fortuitous contamination of two cages revealed an apparent regulatory interaction between the migrant and nonmigrant chromosomes at the esterase-C locus. The migrant chromosomes were very rapidly absorbed into the recipient populations, despite this interaction. This result suggests that the dynamics of migration in populations may be phenomenologically richer than anticipated by simple theory.     

Genetic diversity of six isozyme loci in cultivated barley of Tibet

Summary A random sample of 463 accessions of cultivated barley from the Tibet Hordeum germplasm collection was assayed electorphoretically for genetic diversity at six isozyme loci. Two loci (Acp-1 and Got-1) were found to be monomorphic and extensive variation was detected at the remaining four loci (Est-1, Est-2, Est-3 and Est-4). The allelic composition of Tibetan barley appeared to be distinct as compared to the results of previous studies of barleys from other parts of the world. Partitioning of genetic diversity showed that approximately 96% of the total variation was maintained at the within-subregion level and only about 4% was accounted for by differentiation among the eight subregions. Analysis of multilocus genotypes revealed non-random association of the alleles at the four loci, both in the entire sample and in all the subregions, although the four major multilocus genotypes did not show significant departure from the expectation based on complete random association. The possible causes for the establishment of these multilocus associations were discussed.     

Development and Characterization of Polymorphic EST-SSR and Genomic SSR Markers for Tibetan Annual Wild Barley

Tibetan annual wild barley is rich in genetic variation. This study was aimed at the exploitation of new SSRs for the genetic diversity and phylogenetic analysis of wild barley by data mining. We developed 49 novel EST-SSRs and confirmed 20 genomic SSRs for 80 Tibetan annual wild barley and 16 cultivated barley accessions. A total of 213 alleles were generated from 69 loci with an average of 3.14 alleles per locus. The trimeric repeats were the most abundant motifs (40.82%) among the EST-SSRs, while the majority of the genomic SSRs were di-nuleotide repeats. The polymorphic information content (PIC) ranged from 0.08 to 0.75 with a mean of 0.46. Besides this, the expected heterozygosity (He) ranged from 0.0854 to 0.7842 with an average of 0.5279. Overall, the polymorphism of genomic SSRs was higher than that of EST-SSRs. Furthermore, the number of alleles and the PIC of wild barley were both higher than that of cultivated barley, being 3.12 vs 2.59 and 0.44 vs 0.37. Indicating more polymorphism existed in the Tibetan wild barley than in cultivated barley. The 96 accessions were divided into eight subpopulations based on 69 SSR markers, and the cultivated genotypes can be clearly separated from wild barleys. A total of 47 SSR-containing EST unigenes showed significant similarities to the known genes. These EST-SSR markers have potential for application in germplasm appraisal, genetic diversity and population structure analysis, facilitating marker-assisted breeding and crop improvement in barley.     

Genetic diversity and multilocus associations in Cunninghamia lanceolata (Lamb.) Hook from The People's Republic of China

Open-pollinated seeds were assayed for allozyme polymorphisms at 24 loci to assess genetic diversity and multilocous associations in 16 populations of Cunninghamia lanceolata (Lamb.) Hook in the People's Republic of China. On average, the percentage of polymorphic loci was 88.0, the number of alleles per locus was 3.0, and the expected heterozygosity was 0.394. The distribution of genetic diversity was not correlated with the geographic and climatic variables of the populations. However, allele frequencies correlated linearly with the mean annual temperature of the populations at Mdh-1, Mdh-2, Mnr-2, Pgi-1, and Skdh-1 and with the altitude of the populations at Aph-4 and 6Pg-2. Of the total gene diversity 6% was attributed to among-population differentiation; 94% resided within populations. Two-locus gametic disequilibria were found in 15 of the 16 populations, and higher-order gametic disquilibria were significant in most populations. The gametic disequilibria did not correlate with geographic and climatic variables. The results suggest that population subdivision, founder effect, occurrence across diverse environments, a mating system dominated by inbreeding, and historical events from 2000 years of cultivation are contributing factors in the generation and maintenance of the multilocus genetic structure in this conifer.     

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F(ST) = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13% of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.     

On the origin and domestication history of Barley (Hordeum vulgare)

Remains of barley (Hordeum vulgare) grains found at archaeological sites in the Fertile Crescent indicate that about 10,000 years ago the crop was domesticated there from its wild relative Hordeum spontaneum. The domestication history of barley is revisited based on the assumptions that DNA markers effectively measure genetic distances and that wild populations are genetically different and they have not undergone significant change since domestication. The monophyletic nature of barley domestication is demonstrated based on allelic frequencies at 400 AFLP polymorphic loci studied in 317 wild and 57 cultivated lines. The wild populations from Israel-Jordan are molecularly more similar than are any others to the cultivated gene pool. The results provided support for the hypothesis that the Israel-Jordan area is the region in which barley was brought into culture. Moreover, the diagnostic allele I of the homeobox gene BKn-3, rarely but almost exclusively found in Israel H. spontaneum, is pervasive in western landraces and modern cultivated varieties. In landraces from the Himalayas and India, the BKn-3 allele IIIa prevails, indicating that an allelic substitution has taken place during the migration of barley from the Near East to South Asia. Thus, the Himalayas can be considered a region of domesticated barley diversification.     

Linked vs unlinked markers: multilocus microsatellite haplotype-sharing as a tool to estimate gene flow and introgression

We have explored the use of multilocus microsatellite haplotypes to study introgression from cultivated (Malus domestica) into wild apple (Malus sylvestris), and to study gene flow among remnant populations of M. sylvestris. A haplotype consisted of alleles at microsatellite loci along one chromosome. As destruction of haplotypes through recombination occurs much faster than loss of alleles due to genetic drift, the lifespan of a multilocus haplotype is much shorter than that of the underlying alleles. When different populations share the same haplotype, this may indicate recent gene flow between populations. Similarly, haplotypes shared between two species would be a strong signal for introgression. As the expected lifespan of a haplotype depends on the strength of the linkage, the length [in centiMorgans (cM)] of the haplotype shared contains information on the number of generations passed. This application of shared haplotypes is distinct from using haplotype-sharing to detect association between markers and a certain trait. We inferred haplotypes for four to eight microsatellite loci on Linkage Group 10 of apple from genotype data using the program phase, and then identified those haplotypes shared between populations and species. Compared with a Bayesian analysis of unlinked microsatellite loci using the program structure, haplotype-sharing detected a partially different set of putative hybrids. Cultivated haplotypes present in M. sylvestris were short (< 1.5 cM), indicating that introgression had taken place many generations ago, except for two Belgian plants that contained a haplotype of 47.1 cM, indicating recent introgression. In the estimation of gene flow, F(ST) based on unlinked loci indicated small (0.032-0.058) but statistically significant differentiation between some populations only. However, various M. sylvestris haplotypes were shared in nearly all pairwise comparisons of populations, and their length indicated recent gene flow. Hence, all Dutch populations should be considered as one conservation unit. The added value of using sharing of multilocus microsatellite haplotypes as a source of population genetic information is discussed.     

Importance of genetic characterization and conservation of plant genetic resources: The breeding system and genetic diversity of wild soybean (Glycine soja)

Abstract Plant genetic resources play an important role in the improvement of cultivated plants. To characterize and evaluate the ecological and reproductive features of wild soybean ( Glycine soja Sieb. et Zucc.), which is the most probable ancestor of cultivated soybean ( G. max (L) Merr.), the breeding system and genetic diversity of G. soja were investigated. The extent of natural cross-pollination of G. soja was estimated in four populations along the Omono River in Akita Prefecture, Japan by examining allozyme variation. Although it has been previously believed that G. soja is autogamous, as is cultivated soybean, the mean multilocus outcrossing rate ( t m) estimate was 13%. These values are much higher than the outcrossing rate previously reported for both G. soja and G. max . Frequent visits by honeybees and carpenter bees to flowers were also observed, which supported this conjecture. Furthermore, to evaluate the genetic variation of G. soja as a genetic resource, the genetic structure of 447 populations over Japan were analyzed. Wild soybean populations had a higher degree of variation of isozyme loci. The G ST coefficient of gene differentian values among the sites within the district were particularly high, revealing that the isozyme genotype was greatly different among site populations and homogeneous within the sites. The genetic differentiation among nine districts was observed in the allele frequencies of a few loci, indicating that geographic isolation in the wild soybean population was effectively created through the distance between the districts. The difference in the allele frequency among the districts may be produced under genetic drift. Finally, the importance of the preservation of natural plant populations and the habitats of wild progenitors (i.e. the in situ conservation of plant genetic resources) was emphasized.     

Inferring ancient Agave cultivation practices from contemporary genetic patterns

Several Agave species have played an important ethnobotanical role since prehistory in Mesoamerica and semiarid areas to the north, including central Arizona. We examined genetic variation in relict Agave parryi populations northeast of the Mogollon Rim in Arizona, remnants from anthropogenic manipulation over 600 years ago. We used both allozymes and microsatellites to compare genetic variability and structure in anthropogenically manipulated populations with putative wild populations, to assess whether they were actively cultivated or the result of inadvertent manipulation, and to determine probable source locations for anthropogenic populations. Wild populations were more genetically diverse than anthropogenic populations, with greater expected heterozygosity, polymorphic loci, effective number of alleles and allelic richness. Anthropogenic populations exhibited many traits indicative of past active cultivation: fixed heterozygosity for several loci in all populations (nonexistent in wild populations); fewer multilocus genotypes, which differed by fewer alleles; and greater differentiation among populations than was characteristic of wild populations. Furthermore, manipulated populations date from a period when changes in the cultural context may have favoured active cultivation near dwellings. Patterns of genetic similarity among populations suggest a complex anthropogenic history. Anthropogenic populations were not simply derived from the closest wild A. parryi stock; instead they evidently came from more distant, often more diverse, wild populations, perhaps obtained through trade networks in existence at the time of cultivation.     

Haplotype diversity and population structure in cultivated and wild barley evaluated for Fusarium head blight responses

Fusarium head blight (FHB) is a threat to barley (Hordeum vulgare L.) production in many parts of the world. A number of barley accessions with partial resistance have been reported and used in mapping experiments to identify quantitative trait loci (QTL) associated with FHB resistance. Here, we present a set of barley germplasm that exhibits FHB resistance identified through screening a global collection of 23,255 wild (Hordeum vulgare ssp. spontaneum) and cultivated (Hordeum vulgare ssp. vulgare) accessions. Seventy-eight accessions were classified as resistant or moderately resistant. The collection of FHB resistant accessions consists of 5, 27, 46 of winter, wild and spring barley, respectively. The population structure and genetic relationships of the germplasm were investigated with 1,727 Diversity Array Technology (DArT) markers. Multiple clustering analyses suggest the presence of four subpopulations. Within cultivated barley, substructure is largely centered on spike morphology and growth habit. Analysis of molecular variance indicated highly significant genetic variance among clusters and within clusters, suggesting that the FHB resistant sources have broad genetic diversity. The haplotype diversity was characterized with DArT markers associated with the four FHB QTLs on chromosome 2H bin8, 10 and 13 and 6H bin7. In general, the wild barley accessions had distinct haplotypes from those of cultivated barley. The haplotype of the resistant source Chevron was the most prevalent in all four QTL regions, followed by those of the resistant sources Fredrickson and CIho4196. These resistant QTL haplotypes were rare in the susceptible cultivars and accessions grown in the upper Midwest USA. Some two- and six-rowed accessions were identified with high FHB resistance, but contained distinct haplotypes at FHB QTLs from known resistance sources. These germplasm warrant further genetic studies and possible incorporation into barley breeding programs.     

The use of the expectation-maximization (EM) algorithm for maximum likelihood estimation of gametic frequencies of multilocus polymorphic codominant systems based on sampled population data

Estimation of gametic frequencies in multilocus polymorphic systems based on the numerical distribution of multilocus genotypes in a population sample ("analysis without pedigrees") is difficult because some gametes are not recognized in the data obtained. Even in the case of codominant systems, where all alleles can be recognized by genotypes, so that direct estimation of the frequencies of genes (alleles) is possible ("complete data"), estimation of the frequencies of multilocus gametes based on the data on multilocus genotypes is sometimes impossible, whether population data or even family data are used for studying genotypic segregation or analysis of linkage ("incomplete data"). Such "incomplete data" are analyzed based on the corresponding genetic models using the expectation-maximization (EM) algorithm. In this study, the EM algorithm based on the random-marriage model for a nonsubdivided population was used to estimate gametic frequencies. The EM algorithm used in the study does not set any limitations on the number of loci and the number of alleles of each locus. Locus and alleles are identified by numeration making possible to arrange loops. In each combination of alleles for a given combination of m out of L loci (L is the total number of loci studied), all alleles are assigned value 1, and the remaining alleles are assigned value 0. The sum of zeros and unities for each gamete is its gametic value (h), and the sum of the gametic values of the gametes that form a given genotype is the genotypic value (g) of this genotype. Then, gametes with the same h are united into a single class, which reduces the number of the estimated parameters. In a general case of m loci, this procedure yields m + 1 classes of gametes and 2m + 1 classes of genotypes with genotypic values g = 0, 1, 2, ..., 2m. The unknown frequencies of the m + 1 classes of gametes can be represented as functions of the gametic frequencies whose maximum likelihood estimations (MLEs) have been obtained in all previous EM procedures and the only unknown frequency (Pm(m)) that is to be estimated in the given EM procedure. At the expectation step, the expected frequencies (Fm(g) of the genotypes with genotypic values g are expressed in terms of the products of the frequencies of m + 1 classes of gametes. The data on genotypes are the numbers (ng) of individuals with genotypic values g = 0, 1, 2, 3, ..., 2m. The maximization step is the maximization of the logarithm of the likelihood function (LLF) for ng values. Thus, the EM algorithm is reduced, in each case, to solution of only one equation with one unknown parameter with the use of the ng values, i.e., the numbers of individuals after the corresponding regrouping of the data on the individuals' genotypes. Treatment of the data obtained by Kurbatova on the MNSs and Rhesus systems with alleles C, Cw, c, D, d, E, e with the use of Weir's EM algorithm and the EM algorithm suggested in this study yielded similar results. However, the MLEs of the parameters obtained with the use of either algorithm often converged to a wrong solution: the sum of the frequencies of all gametes (4 and 12 gametes for MNSs and Rhesus, respectively) was not equal to 1.0 even if the global maximum of LLF was reached for each of them (as it was for MNSs with the use of Weir's EM algorithm), with each parameter falling within admissible limits (e.g., [0, min(PN,Ps)] for PNs). The chi 2 function is suggested to be used as a goodness-of-fit function for the distribution of genotypes in a sample in order to select acceptable solutions. However, the minimum of this function only guarantee the acceptability of solutions if all limitations on the parameters are met: the sum of estimations of gametic frequencies is 1.0, each frequency falls within the admissible limits, and the "gametic algebra" is complied with (none of the frequencies is negative).     

Genetic and morphological contrasts between wild and anthropogenic populations of Agave parryi var. huachucensis in south-eastern Arizona

Background and Aims

At least seven species of Agave, including A. parryi, were cultivated prehistorically in Arizona, serving as important sources of food and fibre. Many relict populations from ancient cultivation remain in the modern landscape, offering a unique opportunity to study pre-Columbian plant manipulation practices. This study examined genetic and morphological variation in six A. p. var. huachucensis populations of unknown origin to compare them with previous work on A. parryi populations of known origin, to infer their cultivation history and to determine whether artificial selection is evident in populations potentially managed by early agriculturalists.

Methods

Six A. p. var. huachucensis and 17 A. parryi populations were sampled, and morphometric, allozyme and microsatellite data were used to compare morphology and genetic structure in purportedly anthropogenic and wild populations, as well as in the two taxa. Analysis of molecular variance and Bayesian clustering were performed to partition variation associated with taxonomic identity and hypothesized evolutionary history, to highlight patterns of similarity among populations and to identify potential wild sources for the planting stock.

Key Results A

p. var. huachucensis and A. parryi populations differed significantly both morphologically and genetically. Like A. parryi, wild A. p. var. huachucensis populations were more genetically diverse than the inferred anthropogenic populations, with greater expected heterozygosity, percentage of polymorphic loci and number of alleles. Inferred anthropogenic populations exhibited many traits indicative of past active cultivation: greater morphological uniformity, fixed heterozygosity for several loci (non-existent in wild populations), fewer multilocus genotypes and strong differentiation among populations.

Conclusions

Where archaeological information is lacking, the genetic signature of many Agave populations in Arizona can be used to infer their evolutionary history and to identify potentially fruitful sites for archaeological investigation of ancient settlements and cultivation practices. The same approach can clearly be adopted for other species in similar situations.