Genetic diversity and population structure of six South African Acacia mearnsii breeding populations based on SSR markers

Black wattle (Acacia mearnsii) has great economic value as a commercial source of tannins, timber and a source of firewood for local and international markets. It has been suggested that to maximize the genetic gain of A. mearnsii plantations in South Africa, the gene pool that exist within ICFR needs to be broadened via introduction of new genotypes with diverse traits. In this work, 282 A. mearnsii samples sourced from the ICFR breeding program were genotyped using 11 cross-species SSR markers. Our results showed low to moderate genetic differentiation (F_ST) among the six breeding subpopulations, with positive inbreeding (F_IS) values that could be attributed to an historical inbreeding event. Low levels of relatedness could however indicate some mechanism of inbreeding avoidance. The effects from a recent supplementation of genetic material from two native Australian populations were observed through genetic structuring analyses. Analysis of molecular variance (AMOVA) revealed that significant genetic variation was mainly distributed within populations (75%) and among individuals (23%). The results provide significant information on A. mearnsii population genetic diversity and structure, which can be used for conservation of the current subpopulations and future tree improvement programs.

     

Congruence between environmental parameters, morphology and genetic structure in Australia’s most widely distributed eucalypt, Eucalyptus camaldulensis

Eucalyptus camaldulensis is one of the most widely utilised eucalypts. It is also the only eucalypt that occurs across the Australian continent, playing a key ecological role as fauna habitat and in riverbank stabilisation. Despite its ecological and economic importance, uncertainty remains regarding the delineation of genetic and morphological variants. Nine hundred and ninety trees from 97 populations, representing the species’ geographic range were genotyped using 15 microsatellite loci and patterns of diversity compared with restriction fragment length polymorphisms in 29 of these populations. Both markers showed that despite having a riverine distribution, downstream seed dispersal has had less influence than geographic distance on dispersal patterns. Spatial patterns in the distribution of microsatellite genotypes were compared with environmental parameters and boundaries defined by river systems, drainage basins and proposed subspecies. Significant genetic differences among populations within river systems indicated that rivers should not be treated as a single genetic entity in conservation or breeding programmes. Strong geographic trends were evident with 40% of variation in genetic diversity explained by latitude and moisture index. Isolation by distance and significant correlations between genetic distance and environmental parameters for most loci suggest historical factors have had more influence than selection on current patterns of distribution of genetic diversity. Geographic structuring of molecular variation, together with congruence between genetic and morphological variation indicate that E. camaldulensis should be treated as a number of subspecies rather than a single variable taxon. High levels of genetic diversity and geographic trends in the distribution of variation provide a firm basis for further exploration of the species’ genetic resources.     

Potential genetic consequences of a recent bottleneck in the Amur tiger of

The Amur tiger, Panthera tigris altaica, is a highly endangered felid whose range and population size has been severely reduced in recent times. At present, the wild population is estimated at 490 individuals, having rebounded from the 20–30 tigers remaining following a severe bottleneck in the 1940's. The current study presents preliminary data on the patterns and levels of genetic variation in the mitochondrial DNA control region using DNA extracted from non-invasively sampled faecal material, collected throughout the entire range of P. t. altaica in the Russian Far East. Analysis of 82 scat samples representing at least 27 individuals revealed extremely low levels of CR haplotype diversity, characterized by a single widespread haplotype (96.4%) and two rare variants, each differing by a single step within the hypervariable I (2.4%) and central conserved regions (1.2%), respectively. A comparison with previous data on cytochrome bvariation in 14 captive individuals revealed a potentially greater amount of genetic variation represented in captivity relative to that found in the wild population. The extremely low levels of mitochondrial DNA variation in the wild population is discussed in light of the demographic processes that might have shaped these patterns as well as the potential bias introduced through analysis of fecal samples. These results highlight the continuing need to assess levels of genetic variation even in recovering populations that are increasing in number and underscore the important role that captive breeding programs may play in preserving remnant genetic diversity of endangered species.     

秦岭山区美容杜鹃五个野生种群遗传多样性的ISSR分析

美容杜鹃是杜鹃花科杜鹃花属常绿观花植物,也是秦岭地区的特有属,在维持当地生态平衡方面发挥着重要作用,但目前资源正在遭受严重的破坏。该文通过ISSR分子标记技术对秦岭地区5个美容杜鹃野生种群进行遗传多样性分析,了解美容杜鹃不同种群的遗传分化,从而为美容杜鹃野生种质资源保护策略的制定提供理论依据。用8个ISSR引物对5个天然种群的90个单株进行扩增,共扩增出78条带(平均每条引物产生9.75条带),其中65条带是多态的,多态位点占总位点的百分率为83%。种群总的Nei’s基因多样性指数为0.3386,Shannon信息多态性指数为0.4972,表明美容杜鹃总的遗传多样性水平较高。种群多态位点百分率为82.71%~90.25%,Shannon信息多态性指数为0.4161~0.5867,Nei’s基因多样性指数为0.3044~0.4122,表明美容杜鹃不同种群遗传多样性水平差异较大。其中镇安木王种群和柞水牛背梁种群的遗传多样性水平较高。AMOVA分析表明种群内的遗传变异(91.22%)大于种群间的遗传变异(8.78%)。UPGMA聚类分析表明5个种群的遗传分化程度与地理距离没有相关性。因此建议尽可能地保护美容杜鹃所有的天然种群原生境条件,以最大限度保护其遗传变异。由于镇安和柞水种群具有较高的遗传多样性,因此建议优先对这2个种群实施就地保护和迁地保护。     

基于SSR标记的紫花风铃木群体遗传多样性分析

为揭示紫花风铃木(Handroanthus impetiginosus)的群体遗传变异特征,对广东省6市12个群体72份种质材料进行遗传多样性和群体遗传结构分析。结果表明,9对引物共扩增出123个等位基因位点,引物的平均多态信息量为0.754,具有较高的多态性。12个群体均具有较高的遗传多样性,群体间平均有效等位基因数为3.272个,平均Shannon指数为1.159。AMOVA分析表明群体间遗传分化程度相对较低,群体内遗传分化程度较高。群体的总体遗传分化系数为0.077,处于中等程度。基于Structure分析、主坐标分析和NJ聚类分析均可将12个群体分为2大类群,分组结果具有一定相似性,表明供试紫花风铃木群体遗传结构较为简单。这为紫花风铃木优良种质资源的利用、遗传变异和科学育种提供了理论参考。     

Population structure and genetic diversity of wild <Emphasis Type="Italic">Helianthus</Emphasis> species from Mozambique

The production of sunflower suffered a major decline in Mozambique after its independence in 1975. Civil war, human activities and environmental damage subjected the species to an ecological stress contributing to reduce the number and size of wild populations. As this reduction is often related to a loss of genetic variation we estimated the genetic diversity within and among populations of wild Helianthus from five districts of Mozambique using RAPD markers. The 44 accessions studied grouped into four major clusters exhibiting structured variability with regard to geographic origin. A high level of genetic diversity (He = 0.350 and I = 0.527) was retained at the population level. The genetic variation among populations was high (59.7%), which is consistent with low gene flow (Nm = 0.338). The proportion of total genetic diversity residing among these populations should be kept in mind to devise different conservation strategies in order to preserve these populations. Currently wild Helianthus genetic resources present in Maputo and Sofala are on the edge of extinction mainly due to excessive urbanization. Therefore, conservation of what remains of this plant genetic diversity is essential for sustainable utilization and can be useful for breeding programs.     

陕西唐棣遗传多样性和遗传分化的AFLP分析

采用扩增片段长度多态性分子标记技术对陕西省分布的6个野生唐棣居群的96个个体进行了遗传多样性分析, 以明确野生唐棣资源的亲缘关系,为唐棣资源的保护、良种选育和开发利用提供理论依据。结果显示：(1)从64对引物组合中筛选出8对扩增条带清晰、多态性高的引物组合,共扩增出277条清晰条带,其中多态性条带116条,多态性位点百分率为42.86%。(2)UPGMA聚类、主坐标分析(PCoA)和遗传结构分析结果相似,将6个陕西野生唐棣居群分成2大支,秦岭南北居群间遗传分化明显,且群体间存在一定基因流。(3)分子方差分析(AMOVA)结果显示遗传变异主要存在于居群内(63%),居群间遗传变异为37%。Mantel检验表明陕西唐棣居群地理距离与遗传距离之间无明显相关性(r = 0.192,P = 0.220)。研究表明,AFLP分子标记可以准确、有效地用于唐棣遗传多样性分析;唐棣遗传变异主要来源于居群内,居群间的基因交流有限;陕西野生唐棣遗传多样性水平较低,但部分居群的遗传多样性较高。该研究结果可为野生唐棣资源的保护、良种选育和开发利用提供理论依据。     

Genetic structuring in wild populations of two important medicinal plant species as inferred from AFLP markers

Systematic analysis of germplasm diversity and genetic relationship among cultivars is critical for development of appropriate conservation and breeding strategies. This approach has been applied to gain an insight about the molecular variance that exists in wild population of two important medicinal plant species of India that have a long history of therapeutic usage in herbal medicine. Adhatoda vasica and Andrographis paniculata, members of the family Acanthaceae, have wide geographical and climatic distribution across India suggesting a large amount of genetic diversity available for resource management and breeding programs. In this study we have assessed the genetic diversity of both these species distributed in five varied geo-environmental regions, using Amplified Fragment Length Polymorphism (AFLP) fingerprinting with selected primer combinations and statistical analysis. Cluster analysis and analysis of molecular variance also suggested a very high genetic variation. Detailed analyses of the results predict that the genetic variation found in A. vasica was more discrete that reflected strongly in the populations studied, whereas the genetic variation in A. paniculata was relatively uniform. Considering significantly large sample size and distinctive characteristics of the selected populations, this work contributes valuable insights that can be used to engineer conservation and utilization strategies for these species.     

Isozyme variation inFaidherbia albida (Leguminosae, Mimosoideae)

Genetic variation has been assessed in 30 populations (931 families) ofFaidherbia albida (Leguminosae, Mimosoideae) from across its entire African range, using six isozyme loci identified by five enzyme systems. Among the populations studied a null allele was proposed to explain the absence ofLap-1 activity in populations from southern and eastern Africa. The mean percentage of polymorphic loci per population, the mean number of alleles per locus and the mean genetic diversity within populations were 31.7%, 1.6 and 0.127 respectively. Genetic diversity was greatest in populations from West Africa and lowest in populations from eastern/southern Africa, with Ethiopian/Sudanese populations intermediate. The overall degree of genetic differentiation between populations (G_ST) indicated that approximately 56% of the enzyme variation resided within populations. Clustering of Nei's unbaised genetic distances calculated between all populations produced a dendrogram that generally followed the geographic distribution of the populations. Two major groups were identified that may be considered the eastern/southern African and the Ethiopian/West African clusters. Within the Ethiopian/West African cluster two subclusters could be recognised, one broadly corresponding to those populations from Ethiopia/Sudan and the other to those populations from West Africa. The implications of these results for theories regarding the origin ofF. albida in Africa are discussed.     

Assessment of genetic diversity and population structure of Chinese wild almond, Amygdalus nana, using EST- and genomic SSRs

Amygdalus nana L., commonly known as wild almond, is an endangered wild relative of cultivated almond, which has great potential in almond crop breeding. In this study, we used microsatellite (SSR) loci derived from both expressed sequence tag (EST) and anonymous genomic sequence to explore the genetic diversity and population structure of A. nana in Xinjiang of China. Seven natural populations were collected across the whole distribution of A. nana in China, including populations from both inside (four populations) and outside (three populations) the established protected areas. A total of 22 and 19 alleles were detected from the seven pairs of EST and genomic SSR loci, respectively. Generally, the genomic SSRs showed lower levels of variation than EST-SSRs, which may partially due to the higher cross-species transferability in EST-SSRs than in genomic SSRs. The population-level genetic diversity (A = 1.84, P = 50.00%, H_o = 0.3491, H_E = 0.2271) was lower than cultivated almond and several wild fruit species with similar breeding system. Most of the genetic variation (82.16%) was partitioned within populations. In particular, the population collected from Tacheng County (outside the protected areas) had the highest levels of genetic diversity and had significantly different genetic constitution from other populations.     

Genetic diversity in the US ex situ populations of the endangered Puerto Rican Boa,Chilabothrus inornatus

The Puerto Rican Boa (Chilabothrus inornatus) was placed on the US Endangered Species List in 1970. Progress has been made since to clarify the recovery status of this species, though the design of a new recovery plan must include information regarding genetic variation within and among populations of this species. While measures of genetic diversity in wild populations of this species are finally becoming available, relative genetic diversity represented in ex situ populations is unknown, which hampers efforts to develop an ex situ species management plan. Here, we provide an analysis of genetic diversity in US public and private collections (zoos and breeders) using mitochondrial sequence data and five highly polymorphic nuclear microsatellite loci. We analyzed 50 boas from the US ex situ population and determined overall genetic diversity and relatedness among these individuals. We then compared these data to mitochondrial and microsatellite data obtained from 176 individuals from wild populations across the native range of the species. We found little inbreeding and a large amount of retained genetic diversity in the US ex situ population of C. inornatus relative to wild populations. Genetic diversity in the ex situ population is similar to that found in wild populations. Ours is only the second explicit attempt to characterize genetic diversity at the molecular level in ex situ populations of boid snakes. We anticipate that these results will inform current breeding strategies as well as offer additional information that will facilitate the continuation of ex situ conservation breeding or management in boas.     

Analysis of genetic diversity in the endangered Hucho taimen from China

Hucho taimen are listed as endangered in China. The population size has declined recently, prompting an increase in the level of listing from grade three in 2002 to grade five in 2006. We analyzed the genetic diversity of wild populations using 17 microsatellite markers to establish a scientific basis for conservation of this species. We collected tissue samples from four populations in the Heilongjiang River basin: Huma River (HM), Hutou (HT), Haiqing (HQ), and Zhuaji (ZJ). A total of 21 loci were amplified, 18 of which were polymorphic. The number of alleles per locus ranged from 2 to 9 (mean: 4.1905). There were 13 highly polymorphic loci and 5 moderately polymorphic loci. Analysis of five genetic diversity parameters (Na, Ne, Ho, He, and PIC) suggested moderate levels of diversity within the populations. The populations were ranked HT > HQ > ZJ > HM, but the differences in diversity were not statistically significant (P > 0.05). A comparison of variation among all four populations suggested Hardy–Weinberg disequilibrium at 20% of the loci. Genetic differentiation (Fst) was 0.0644 and the gene flow among populations was estimated at 3.36 individuals per generation. The majority of diversity (93.88%) occurred among individuals within a population. In contrast, relatively little (6.12%) of the genetic diversity was distributed between the populations. An analysis of genetic differentiation and genetic distance between pairs of populations revealed that both parameters were higher in comparisons of the HM population to the HT, HQ, and ZJ populations than among the three latter populations. This suggests that the HM population has a distinct genetic structure. We hypothesize that habitat degradation and excessive fishing, not low genetic diversity, has caused the decline in H. taimen populations. However, this species should be protected from further declines in genetic diversity.     

广东省猴耳环遗传多样性研究

为了解猴耳环(Archidendron clypearia)种质资源的遗传多样性,以广东省12个野生猴耳环群体的146份种质资源为材料,采用SSR分子标记技术对其遗传多样性和亲缘关系进行分析。结果表明,21对SSR引物共检测到249个等位基因,平均每对SSR引物检测的等位基因数(Na)为11.857,有效等位基因数(Ne)为3.500,期望杂合度(He)为0.718,多态信息含量(PIC)为0.676;12个群体中博罗群体的Shannon多样性指数(I=0.528)和有效等位基因数(Ne=0.716)均最大,是遗传多样性最丰富的群体;群体间的遗传分化系数为0.071,AMOVA分析表明,猴耳环的遗传变异主要在群体内(97%),群体内的遗传分化大于群体间。聚类分析表明,遗传系数在0.16时,可将12个群体分为6大类,与主坐标分析的结果大致相同。这为发掘、利用与保护猴耳环群体种质资源,开展猴耳环优良品种的遗传育种提供重要的理论依据。     

Genetic variation in the endangered <Emphasis Type="Italic">Anisodus tanguticus</Emphasis> (Solanaceae), an alpine perennial endemic to the Qinghai-Tibetan Plateau

We used random amplified polymorphic DNA markers (RAPDs) to assess genetic variation between- and within-populations of Anisodus tanguticus (Solanaceae), an endangered perennial endemic to the Qinghai-Tibetan Plateau with important medicinal value. We recorded a total of 92 amplified bands, using 12 RAPD primers, 76 of which (P = 82.61%) were polymorphic, and calculated values of H_t and H_sp of 0.3015 and 0.4459, respectively, suggesting a remarkably high rate of genetic variation at the species level. The average within-population diversity also appeared to be high, with P, H_e and H_pop values of 55.11%, 0.1948 and 0.2918, respectively. Analyses of molecular variance (AMOVA) showed that among- and between-population genetic variation accounted for 67.02% and 32.98% of the total genetic variation, respectively. In addition, Nei’s coefficient of differentiation (G_ST) was found to be high (0.35), confirming the relatively high level of genetic differentiation among the populations. These differentiation coefficients are higher than mean corresponding coefficients for outbreeding species, but lower than reported coefficients for some rare species from this region. The genetic structure of A. tanguticus has probably been shaped by its breeding attributes, biogeographic history and human impact due to collection for medicinal purposes. The observed genetic variations suggest that as many populations as possible should be considered in any planned in situ or ex situ conservation programs for this species.     

An AFLP analysis of genetic diversity and structure of Caragana microphylla populations in Inner Mongolia steppe, China

Genetic diversity and structure of five natural populations of Caragana microphylla from the Inner Mongolia steppe were estimated using AFLP markers. Five pairs of primers generated a total of 312 bands among 90 individuals, with percentage of polymorphic bands (PPB) being 63% at the population level and 76% at the species level, respectively. The genetic diversity within populations was correlated significantly with the soil N:P ratio. AMOVA analysis revealed high genetic variations within populations (95.5%). The estimated number of migrants per generation (Nm) was 10.72, indicating a high level of gene flow among populations. There was no significant correlation (r = 0.36) between genetic distance and geographical distance. These results were discussed in terms of eco-geographical variations among populations, together with the life history traits and breeding system of the species. The knowledge obtained may have important implications for better conservation and wise use of the vegetation dominate by C. microphylla.     

Hierarchical genetic structure in fragmented populations of the Little Pocket Mouse (Perognathus longimembris) in Southern California

The geographic genetic structure, based onsequence variation of an 810 base pair fragmentof the mitochondrial cytochrome b gene,is described for populations of five subspeciesof the Little Pocket Mouse, Perognathuslongimembris, from Southern California. Oneof these, P. l. pacificus (Pacific PocketMouse), is listed as Endangered by the U.S.Federal Government. Sixty-two uniquehaplotypes were recovered from 99 individualssampled. Phylogenetic analyses of thesevariants do not identify regionallyreciprocally monophyletic lineages concordantwith the current subspecies designations, butmost haplotypes group by subspecies in networksgenerated by either statistical parsimony ormolecular variance parsimony. Moreover, asubstantial proportion of the total pool ofhaplotype variation is attributed to thesesubspecies, or to local populations withingeographic segments of each, indicating theirrelative evolutionary independence. The pooledextant populations of the endangered PacificPocket Mouse exhibit the same levels ofnucleotide and haplotype diversity as other,presumptively less-impacted populations ofadjacent subspecies, although the sample fromDana Point, Orange County, has markedly lowhaplotype diversity in comparison to allothers. These populations also show a geneticsignature of population expansion rather thanone of decline. Both pieces of evidence are atodds with current empirical populationestimates, which reinforces the fact thatpresent-day patterns of genetic diversity arethe product of coalescent history and will notnecessarily reflect recent anthropogenic, orother, perturbations. Comparison of haplotypevariation within and among extant populationsof the Pacific Pocket Mouse with those obtainedfrom museum samples collected more than 70years ago suggests that the pattern ofpopulation differentiation and diversity was inplace before the post-World War II exponentialurbanization of Southern California.     

Genetic Variation in the Endangered Endemic Species Cycas fairylakea (Cycadaceae) in China and Implications for Conservation

Cycas fairylakea is an endangered endemic species in China. Genetic diversity within and among four natural populations of this species in China was investigated using amplified fragments length polymorphism (AFLP). A moderate to low level of intraspecific genetic diversity was detected in this species (at population level: P = 39.57 %, H₀ = 0.244; at species level: P = 60.22%, H₀ = 0.356). The among-population component accounted for, respectively, 25.7 and 31.5% of the genetic variation, according to AMOVA and Shannon’s index, indicating most of the genetic variation was found between individuals within populations. All four populations have opposite pyramid age structure, and few coning individuals, which is still decreasing. Possibly because of habitat degradation and environmental pollution, plant diseases and insect pests in the populations were extremely serious, suggesting that the main factors threatening the survival of C. fairylakea populations were not genetic variation, but human activities and the breeding system of this species.     

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (H_E) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low‐density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.     

Low genetic diversity and significant population structuring in the relict Amentotaxus argotaenia complex (Taxaceae) based on ISSR fingerprinting

Amentotaxus, a genus of the Taxaceae, represents an ancient lineage that has long existed in Eurasia. All Amentotaxus species experienced frequent population expansion and contraction over periodical glaciations in Tertiary and Quaternary. Among them, Amentotaxus argotaenia complex consists of three morphologically alike species, A. argotaenia, Amentotaxus yunnanensis, and Amentotaxus formosana. This complex is distributed in the subtropical region of mainland China and Taiwan where many Pleistocene refugia have been documented. In this study, genetic diversity and population structuring within and between species were investigated based on the inter-simple sequence repeats (ISSR) fingerprinting. Mean genetic diversity within populations was estimated in three ways: (1) the percentage of polymorphic loci out of all loci (P) (2) Neis unbiased expected heterozygosity (He), and (3) Shannons index of phenotypic diversity. For a total of 310 individuals of 15 populations sampled from the three species, low levels of ISSR genetic variation within populations were detected, with P=4.66–16.58%, He=0.0176–0.0645 and Hpop=0.0263–0.0939, agreeing with their seriously threatened status. AMOVA analyses revealed that the differences between species only accounted for 27.38% of the total variation, whereas differences among populations and within populations were 57.70 and 14.92%, respectively, indicating substantial isolation between the patch-like populations. A neighbor-joining tree identified a close affinity between A. yunnanensis and A. formosana. Genetic drift due to small population size, plus limited current gene flow, resulted in significant genetic structuring. Low levels of intrapopulational genetic variation and considerable interpopulational divergence were also attributable to demographic bottlenecks during and/or after the Pleistocene glaciations.     

Correlation between hordatine accumulation, environmental factors and genetic diversity in wild barley (Hordeum spontaneum C. Koch) accessions from the Near East Fertile Crescent

Wild barley shows a large morphological and phenotypic variation, which is associated with ecogeographical factors and correlates with genotypic differences. Diversity of defense related genes and their expression in wild barley has been recognized and has led to attempts to exploit genes from H. spontaneum in breeding programs. The aim of this study was to determine the variation in the accumulation of hordatines, which are Hordeum-specific preformed secondary metabolites with strong and broad antimicrobial activity in vitro, in 50 accessions of H. spontaneum from different habitats in Israel. Differences in the accumulation of hordatines in the seedling stage were significant between different H. spontaneum genotypes from different regional locations and micro-sites. Variation in the hordatine accumulation within genotypes was between 9% and 45%, between genotypes from the same location between 13% and 38%, and between genotypes from different locations up to 121%. Principal component analysis showed that water related factors explain 39%, temperature related factors explain 33% and edaphic factors account for 11% of the observed variation between the populations of H. spontaneum. Genetic analysis of the tested accessions with LP-PCR primers that are specific for genes involved in the biosynthetic pathway of hordatines showed tight correlations between hordatine abundance and genetic diversity of these markers. Multiple regression analyses indicated associations between genetic diversity of genes directly involved in hordatine biosynthesis, ecogeographical factors and the accumulation of hordatines.