Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Heritability and quantitative genetic divergence of serotiny,a fire-persistence plant trait

Background and Aims

Although it is well known that fire acts as a selective pressure shaping plant phenotypes, there are no quantitative estimates of the heritability of any trait related to plant persistence under recurrent fires, such as serotiny. In this study, the heritability of serotiny in Pinus halepensis is calculated, and an evaluation is made as to whether fire has left a selection signature on the level of serotiny among populations by comparing the genetic divergence of serotiny with the expected divergence of neutral molecular markers (Q_ST–F_ST comparison).

Methods

A common garden of P. halepensis was used, located in inland Spain and composed of 145 open-pollinated families from 29 provenances covering the entire natural range of P. halepensis in the Iberian Peninsula and Balearic Islands. Narrow-sense heritability (h²) and quantitative genetic differentiation among populations for serotiny (Q_ST) were estimated by means of an ‘animal model’ fitted by Bayesian inference. In order to determine whether genetic differentiation for serotiny is the result of differential natural selection, Q_ST estimates for serotiny were compared with F_ST estimates obtained from allozyme data. Finally, a test was made of whether levels of serotiny in the different provenances were related to different fire regimes, using summer rainfall as a proxy for fire regime in each provenance.

Key Results

Serotiny showed a significant narrow-sense heritability (h²) of 0·20 (credible interval 0·09–0·40). Quantitative genetic differentiation among provenances for serotiny (Q_ST = 0·44) was significantly higher than expected under a neutral process (F_ST = 0·12), suggesting adaptive differentiation. A significant negative relationship was found between the serotiny level of trees in the common garden and summer rainfall of their provenance sites.

Conclusions

Serotiny is a heritable trait in P. halepensis, and selection acts on it, giving rise to contrasting serotiny levels among populations depending on the fire regime, and supporting the role of fire in generating genetic divergence for adaptive traits.     

Contrasting Levels of Variation in Neutral and Quantitative Genetic Loci on Island Populations of Moor Frogs (<Emphasis Type="Italic">Rana arvalis</Emphasis>)

Reduced levels of genetic variability and a prominent differentiation in both neutral marker genes and phenotypic traits are typical for many island populations as compared to their mainland conspecifics. However, whether genetic diversity in neutral marker genes reflects genetic variability in quantitative traits, and thus, their evolutionary potential, remains typically unclear. Moreover, the phenotypic differentiation on islands could be attributable to phenotypic plasticity, selection or drift; something which seldom has been tested. Using eight polymorphic microsatellite loci and quantitative genetic breeding experiments we conducted a detailed comparison on genetic variability and differentiation between Nordic islands (viz. Gotland, Öland and Læsø) and neighbouring mainland populations of moor frogs (Rana arvalis). As expected, the neutral variation was generally lower in island than in mainland populations. But as opposed to this, higher levels of additive genetic variation (V _A) in body size and tibia length were found on the island of Gotland as compared to the mainland population. When comparing the differentiation seen in neutral marker genes (F _ST) with the differentiation in genes coding quantitative traits (Q _ST) two different evolutionary scenarios were found: while selection might explain a smaller size of moor frogs on Gotland, the differentiation seen in tibia length could be explained by genetic drift. These results highlight the limited utility of microsatellite loci alone in inferring the causes behind an observed phenotypic differentiation, or in predicting the amount of genetic variation in ecologically important quantitative traits.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

Biased morph ratios and skewed mating success contribute to loss of genetic diversity in the distylous Pulmonaria officinalis

Background and Aims

In heterostylous plant species, skewed morph ratios are not uncommon and may arise from a range of factors. Despite the recognized importance of skewed morph ratios on overall reproductive success within populations, little is known about the impact of skewed morph ratios on population genetic diversity and differentiation in heterostylous species. This study specifically aimed to clarify the effect of population size and morph bias on population genetic diversity and differentiation in the temperate forest herb Pulmonaria officinalis. This species is characterized by a distylous breeding system and shows morph-specific differences in reproductive success.

Methods

Genetic diversity was determined for 27 P. officinalis populations in northern Belgium by using eight recently developed microsatellite markers. Multiple regressions were used to assess the relationship between genetic diversity, morph bias and population size, and F_ST-values were calculated for short- and long-styled morphs separately to study genetic differentiation as a function of morph type.

Key Results

For all genetic measures used, morph bias was more important in explaining patterns of genetic diversity than population size, and in all cases patterns of population genetic diversity followed a quadratic function, which showed a symmetrical decrease in genetic diversity with increasing morph bias. However, probably due to the reproductive advantage of L-morphs relative to S-morphs, maximum genetic diversity was found in populations showing an excess of L-morphs (60·7 % L-morph). On the other hand, no significant difference in pairwise genetic distances between populations was observed between L- (0·107) and S-morphs (0·106).

Conclusions

Our results indicate that significant deviations from equal morph ratios not only affect plant reproductive success but also population genetic diversity of heterostylous plant species. Hence, when defining conservation measures for populations of heterostylous plant species, morph ratios should be considered as an important trait affecting their long-term population viability.     

Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.     

Morphological and microsatellite diversity associated with ecological factors in natural populations of Medicago laciniata Mill. (Fabaceae)

Genetic variability in 10 natural Tunisian populations of Medicago laciniata were analysed using 19 quantitative traits and 12 polymorphic microsatellite loci. A large degree of genetic variability within-populations and among-populations was detected for both quantitative characters and molecular markers. High genetic differentiation among populations for quantitative traits was seen, with Q _ST = 0.47, and F _ST = 0.47 for microsatellite markers. Several quantitative traits displayed no statistical difference in the levels of Q _ST and F _ST . Further, significant correlations between quantitative traits and eco-geographical factors suggest that divergence in the traits among populations may track environmental differences. There was no significant correlation between genetic variability at quantitative traits and microsatellite markers within populations. The site-of-origin of eco-geographical factors explain between 18.13% and 23.40% of genetic variance among populations at quantitative traits and microsatellite markers, respectively. The environmental factors that most influence variation in measured traits among populations are assimilated phosphorus (P₂0₅) and mean annual rainfall, followed by climate and soil texture, altitude and organic matter. Significant associations between eco-geographical factors and gene diversity, H _e , were established in five-microsatellite loci suggesting that these simple sequence repeats (SSRs) are not necessarily biologically neutral.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

RAD sequencing of common whelk,Buccinum undatum,reveals fine‐scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F _ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation‐by‐distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome‐wide association studies.     

Analysis of Genetic Variation in Natural Populations of Medicago truncatula of Southern Tunisian Ecological Areas, Using Morphological Traits and SSR Markers

We used 19 quantitative traits and 14 microsatellite markers (SSRs) to analyze the genetic variation in four natural populations of the model legume Medicago truncatula sampled in southern Tunisia. The greatest genetic variation of quantitative traits and molecular markers occurred within populations (>71%). In contrast to quantitative population differentiation (Q _ST ?=?0.09), a high level of molecular differentiation (F _ST ?=?0.23) was found among populations. The majority of quantitative traits exhibited Q _ST values significantly less than F _ST values, suggesting that selection may be acting to suppress differentiation for these traits. There was no significant correlation between genetic variation of quantitative traits and molecular markers within populations. On the other hand, significant correlations were found between measured quantitative characters and the site-of-origin environmental factors. The eco-geographical factors with the greatest influence on the variation of measured traits among populations were altitude, followed by soil texture, assimilated phosphorus (P₂O₅) and organic matter. Nevertheless, there were no consistent patterns of associations between gene diversity (He) and eco-geographical factors.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Divergence at neutral and non-neutral loci in Drosophila buzzatii populations and their hybrids

The impact of intraspecific hybridisation on fitness and morphological traits depends on the history of natural selection and genetic drift, which may have led to differently coadapted gene-complexes in the parental populations. The divergence at neutral and non-neutral loci between populations can be evaluated by estimating F_ST and Q_ST respectively, and hence give an estimate of drift and selection in the populations. Here we investigate (1) whether divergence between populations in quantitative traits (wing size and shape) can be attributed to selection or drift alone, (2) The impact of intraspecific hybridisation on estimators for divergence at neutral (F_ST) and non-neutral loci (Q_ST) in hybrids, (3) If measurement of shape is more informative than size in order to detect divergence in quantitative traits between populations. The aims were addressed by performing two hybridisations between three populations of Drosophila buzzatii, one between populations from Argentina and the Canary Islands (separated for 200 years), and the other between populations from Argentina and Australia (separated for 80 years). We observed the highest divergence at neutral loci between the Argentinean and Canary Island populations, but highest morphological divergence between the Argentinean and Australian populations, indicating that natural selection is acting on the wings. Divergence based on Q_ST measures in the hybrids was sensitive towards increased phenotypic variance (σ²p) within groups and should be used with care when σ²p of populations differ. Our results indicate that measures of shape give a better estimate of divergence at the underlying quantitative traits loci than measures of size.     

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (V_A) and trait differentiation (Q_ST) from small to large population size. Alternatively, small populations might exhibit larger variation in V_A and Q_ST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating V_A and Q_ST using nine fragmented populations of brook trout varying 50‐fold in census size N (179–8416) and 10‐fold in effective number of breeders, N_b (18–135). Across 15 traits, no evidence was found for consistent differences in V_A and Q_ST with population size and almost no evidence for increased variability of V_A or Q_ST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.     

Detecting Signatures of Selection Through Haplotype Differentiation Among Hierarchically Structured Populations

The detection of molecular signatures of selection is one of the major concerns of modern population genetics. A widely used strategy in this context is to compare samples from several populations and to look for genomic regions with outstanding genetic differentiation between these populations. Genetic differentiation is generally based on allele frequency differences between populations, which are measured by F_ST or related statistics. Here we introduce a new statistic, denoted hapFLK, which focuses instead on the differences of haplotype frequencies between populations. In contrast to most existing statistics, hapFLK accounts for the hierarchical structure of the sampled populations. Using computer simulations, we show that each of these two features—the use of haplotype information and of the hierarchical structure of populations—significantly improves the detection power of selected loci and that combining them in the hapFLK statistic provides even greater power. We also show that hapFLK is robust with respect to bottlenecks and migration and improves over existing approaches in many situations. Finally, we apply hapFLK to a set of six sheep breeds from Northern Europe and identify seven regions under selection, which include already reported regions but also several new ones. We propose a method to help identifying the population(s) under selection in a detected region, which reveals that in many of these regions selection most likely occurred in more than one population. Furthermore, several of the detected regions correspond to incomplete sweeps, where the favorable haplotype is only at intermediate frequency in the population(s) under selection.     

Genetic differentiation and species cohesion in two widespread Central American Begonia species

Begonia is one of the ten largest plant genera, with over 1500 species. This high species richness may in part be explained by weak species cohesion, which has allowed speciation by divergence in allopatry. In this study, we investigate species cohesion in the widespread Central American Begonia heracleifolia and Begonia nelumbiifolia, by genotyping populations at microsatellite loci. We then test for post-zygotic reproductive barriers using experimental crosses, and assess whether sterility barriers are related to intraspecific changes in genome size, indicating major genome restructuring between isolated populations. Strong population substructure was found for B. heracleifolia (F_ST=0.364, F′_ST=0.506) and B. nelumbiifolia (F_ST=0.277, F′_ST=0.439), and Bayesian admixture analysis supports the division of most populations into discrete genetic clusters. Moderate levels of inferred selfing (B. heracleifolia s=0.40, B. nelumbiifolia s=0.62) and dispersal limitation are likely to have contributed to significant genetic differentiation (B. heracleifolia Jost''s D=0.274; B. nelumbiifolia D=0.294). Interpopulation crosses involving a divergent B. heracleifolia population with a genome size ∼10% larger than the species mean had a ∼20% reduction in pollen viability compared with other outcrosses, supporting reproductive isolation being polymorphic within the species. The population genetic data suggest that Begonia populations are only weakly connected by gene flow, allowing reproductive barriers to accumulate between the most isolated populations. This supports allopatric divergence in situ being the precursor of speciation in Begonia, and may also be a common speciation mechanism in other tropical herbaceous plant groups.     

Estimating Population-Level Coancestry Coefficients by an Admixture F Model

In this article, we develop an admixture F model (AFM) for the estimation of population-level coancestry coefficients from neutral molecular markers. In contrast to the previously published F model, the AFM enables disentangling small population size and lack of migration as causes of genetic differentiation behind a given level of F_ST. We develop a Bayesian estimation scheme for fitting the AFM to multiallelic data acquired from a number of local populations. We demonstrate the performance of the AFM, using simulated data sets and real data on ninespine sticklebacks (Pungitius pungitius) and common shrews (Sorex araneus). The results show that the parameterization of the AFM conveys more information about the evolutionary history than a simple summary parameter such as F_ST. The methods are implemented in the R package RAFM.     

基于形态及分子标记的濒危植物夏蜡梅自然居群的遗传变异研究

比较了濒危植物夏蜡梅(Sinocalycanthus chinensis)大明山、大雷山、龙须山3个居群的果实与种子的形态变异,采用ISSR分子标记技术分析DNA序列的变异,综合评价夏蜡梅3个居群的遗传变异。结果如下:果实性状中,果柄长、果实长、果实重、每果种子数等指标均以大明山居群最大、大雷山居群次之、龙须山居群最小,它们之间差异显著(P<0.05)。种子重、种子长、种子宽、种子厚等指标均以大雷山居群最大,与大明山居群和龙须山居群差异显著(P<0.05)。基于果实形态特征的表型分化系数(VST)在0.5518—0.9750之间,平均为0.8930,所有果实形态指标的变异大部分存在于居群间。基于种子形态特征的VST在0.1669—0.8678之间,平均为0.6240,除种子长外,其他种子形态指标的变异也大部分存在于居群间。ISSR分析表明,3个居群的多态位点百分率、Shannon信息指数、Nei基因多样性均是大明山居群最高、龙须山居群次之、大雷山居群最低。居群间的遗传分化系数(GST)为0.6050。基于ISSR分子标记数据的聚类结果显示大明山居群先与龙须山居群聚在一起,再与大雷山居群相聚,这与基于种子性状特征的聚类结果相似,而与基于果实性状特征的聚类结果不同。3个居群的种子形态特征变异系数与遗传多样性水平具有显著的正相关,且基于果实形态特征估算的VST要高于GST,而基于种子形态特征估算的VST仅略高于GST,表明夏蜡梅不同居群间的种子形态变异主要由遗传变异所造成的,而环境因子在果实形态变异中起了重要作用。Mantel检验显示3个居群基于分子标记数据的遗传距离矩阵和基于果实及种子形态特征的欧氏距离矩阵之间的相关性不显著,表明果实与种子在居群间出现的表型分化除了受遗传因素影响外,还受到其它环境因子的强烈影响。     

Vulnerability to cavitation,hydraulic efficiency,growth and survival in an insular pine (Pinus canariensis)

Background and Aims

It is widely accepted that hydraulic failure due to xylem embolism is a key factor contributing to drought-induced mortality in trees. In the present study, an attempt is made to disentangle phenotypic plasticity from genetic variation in hydraulic traits across the entire distribution area of a tree species to detect adaptation to local environments.

Methods

A series of traits related to hydraulics (vulnerability to cavitation and hydraulic conductivity in branches), growth performance and leaf mass per area were assessed in eight Pinus canariensis populations growing in two common gardens under contrasting environments. In addition, the neutral genetic variability (F_ST) and the genetic differentiation of phenotypic variation (Q_ST) were compared in order to identify the evolutionary forces acting on these traits.

Key Results

The variability for hydraulic traits was largely due to phenotypic plasticity. Nevertheless, the vulnerability to cavitation displayed a significant genetic variability (approx. 5 % of the explained variation), and a significant genetic × environment interaction (between 5 and 19 % of the explained variation). The strong correlation between vulnerability to cavitation and survival in the xeric common garden (r = –0·81; P < 0·05) suggests a role for the former in the adaptation to xeric environments. Populations from drier sites and higher temperature seasonality were less vulnerable to cavitation than those growing at mesic sites. No trade-off between xylem safety and efficiency was detected. Q_ST of parameters of the vulnerability curve (0·365 for P₅₀ and the slope of the vulnerability curve and 0·452 for P₈₈) differed substantially from F_ST (0·091), indicating divergent selection. In contrast, genetic drift alone was found to be sufficient to explain patterns of differentiation for xylem efficiency and growth.

Conclusions

The ability of P. canariensis to inhabit a wide range of ecosystems seemed to be associated with high phenotypic plasticity and some degree of local adaptations of xylem and leaf traits. Resistance to cavitation conferred adaptive potential for this species to adapt successfully to xeric conditions.     

Comparison of the extent of genetic variation of Vallisneria natans and its sympatric congener V. spinulosa in lakes of the middle-lower reaches of the Yangtze River

Vallisneria natans and Vallisneria spinulosa are two morphologically very similar and sympatrically dominant submerged macrophytes in lakes of the middle-lower reaches of the Yangtze River. Genetic variation was compared based on a total of 196 individuals from six V. natans populations and 201 individuals from seven V. spinulosa populations. Using eight ISSR primers, a total of 139 and 129 DNA fragments were generated with 121 being polymorphic in V. natans and 99 in V. spinulosa. The two species maintained higher genetic variation both at the species and population levels in comparison with other aquatic macrophytes. A higher level of genetic diversity among populations was found in V. natans than in V. spinulosa: the percentage of polymorphic loci (PPL) in V. natans was 52-62% vs. 38-47% in V. spinulosa; gene diversity (H) was 0.21 in V. natans vs. 0.17 in V. spinulosa.Both an analysis of molecular variance (AMOVA) and F-estimation (F_ST) showed that most of the total genetic variation resided within populations of both species (AMOVA: 85% and 80%; F_ST: 0.132 and 0.202), indicating low genetic differentiation between populations. Principal coordinates analysis (PCA) indicated evident gene flow between populations of both species. The outcrossing reproductive mode and pervasive gene flow might have played important roles in maintaining high genetic diversity and in shaping low population differentiation of the two Vallisneria species, while the extent of clonal growth might account for the different levels of population divergence between them.     

Genetic variation in cultivated Rheum tanguticum populations

To examine whether cultivation reduced genetic variation in the important Chinese medicinal plant Rheum tanguticum, the levels and distribution of genetic variation were investigated using ISSR markers. Fifty-eight R. tanguticum individuals from five cultivated populations were studied. Thirteen primers were used and a total of 320 DNA bands were scored. High levels of genetic diversity were detected in cultivated R. tanguticum (PPB = 82.19, H = 0.2498, H_B = 0.3231, I = 0.3812) and could be explained by the outcrossing system, as well as long-lived and human-mediated seed exchanges. Analysis of molecular variance (AMOVA) showed that more genetic variation was found within populations (76.1%) than among them (23.9%). This was supported by the coefficient of gene differentiation (G_st = 0.2742) and Bayesian analysis (θ_B = 0.1963). The Mantel test revealed no significant correlation between genetic and geographic distances among populations (r = 0.1176, p = 0.3686). UPGMA showed that the five cultivated populations were separated into three clusters, which was in good accordance with the results provided by the Bayesian software STRUCTURE (K = 3). A short domestication history and no artificial selection may be an effective way of maintaining and conserving the gene pools of wild R. tanguticum.