Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants

The Diversity of Melia azedarach L. from China Based on Transcriptome-Developed SSR Marker

Melia azedarach L. is a native tree species that can be used in a comprehensive way and is widely distributed in all provinces south of the Yellow River in China. Genetic diversity analysis of different M. azedarach germplasm sources is an important basic work for the selection, evaluation, and genetic improvement of M. azedarach germplasm resources. In this study, 100 pairs of SSR primers were designed and synthesized based on M. azedarach transcriptome data, and 16 pairs of reliable SSR primers were finally selected. The developed primers were used to analyze the genetic diversity of M. azedarach from 15 sources in 10 provinces in East, Central, and South China. The results showed that the frequency of the M. azedarach transcriptome SSR loci was high, and the distribution density was high. There were 15 sources of M. azedarach genetic diversity at a moderate level, and genetic variation was mainly present within the sources. The present study further enriches the existing SSR marker database of the M. azedarach family and can provide a reference for genetic diversity analysis and molecularly assisted breeding of M. azedarach plants at the genomic level.

Genetic Structure of Native Blue Honeysuckle Populations in the Western and Eastern Eurasian Ranges

Blue honeysuckle (Lonicera caerulea L.) is a promising berry crop producing edible early-ripening berries with a valuable chemical composition. We evaluated the genetic diversity of native L. caerulea populations from the western (Baltic states) and eastern (the Russian Far East and Japan) edges of the Eurasian range using inter-simple sequence repeat (ISSR) and chloroplast DNA (psbA-trnH and trnL-trnF) markers. The genetic relationships of populations and genotypes were analyzed using principal coordinate and cluster analyses (neighbor joining and Bayesian clustering). Sampling was carried out in two disjunct areas of this circumpolar species and the analyses showed clustering of individuals and populations according to geographic origin. The analysis of genetic structure based on ISSR markers showed that the studied populations of L. caerulea were highly differentiated. However, sequence analysis of two chloroplast DNA (cpDNA) regions revealed no phylogeographic structure among the populations. We also found that the eastern populations of blue honeysuckle had significantly greater genetic diversity parameters than the populations from the Baltic region. This finding correlates with the endangered status of blue honeysuckle in the Baltic states.

Impacts of Chronic Habitat Fragmentation on Genetic Diversity of Natural Populations of Prunus persica in China

Wild peach is an important resource for improving existing peach varieties. However, the extant populations of wild peach show fragmented distribution due to human disturbance and geographic isolation. In this study, we used natural populations (or wild populations) of Prunus persica (Rosaceae) to assess the genetic effects of habitat fragmentation. A total of 368 individuals sampled from 16 natural populations were analyzed using 23 polymorphic simple sequence repeat (SSR) markers. Prunus persica maintained low within-population genetic variation and high level of genetic differentiation. Two genetic clusters were revealed based on three different methods (UPGMA, PCoA, and STRUCTURE). All populations showed a significant heterozygosity deficiency and most extant populations experienced recent reduction in population size. A significant isolation by distance (IBD) was observed with Mantel's test. Compared to historical gene flow, contemporary gene flow was restricted among the studied populations, suggesting a decrease in gene flow due to habitat fragmentation. Habitat fragmentation has impacted population genetic variation and genetic structure of P. persica. For breeding and conservation purpose, collecting as many individuals as possible from multiple populations to maximize genetic diversity was recommended during the process of germplasm collection. In addition, populations from central China had higher genetic diversity, suggesting these populations should be given priority for conservation and germplasm collection.

Adaptive Divergence without Distinct Species Relationships Indicate Early Stage Ecological Speciation in Species of the Rhododendronpseudochrysanthum Complex Endemic to Taiwan

The testing association of environmental variables with genetic and epigenetic variation could be crucial to deciphering the effects of environmental factors playing roles as selective drivers in ecological speciation. Although ecological speciation may occur in closely related species, species boundaries may not be established over a short evolutionary timescale. Here, we investigated the genetic and epigenetic variations using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP), respectively, and tested their associations with environmental variables in populations of four closely related species in the R. pseudochrysanthum complex. No distinctive species relationships were found using genetic clustering analyses, neighbor-joining tree, and neighbor-net tree based on the total AFLP variation, which is suggestive of the incomplete lineage sorting of ancestral variation. Nonetheless, strong isolation-by-environment and adaptive divergence were revealed, despite the significant isolation-by-distance. Annual mean temperature, elevation, normalized difference vegetation index, and annual total potential evapotranspiration were found to be the most important environmental variables explaining outlier genetic and epigenetic variations. Our results suggest that the four closely related species of the R. pseudochrysanthum complex share the polymorphism of their ancestor, but reproductive isolation due to ecological speciation can occur if local environmental divergence persists over time.

Molecular and metabolic traits of some Egyptian species of Cassia L. and Senna Mill (Fabaceae-Caesalpinioideae)

The genus Cassia and Senna have been classified under subfamily Caesalpinioideae of family Fabaceae (Leguminosae) of order Fabales. There is a scarce taxonomical studies of the genus Cassia and Senna inhabiting Egyptian environments, thus, the main objective of the current was to revise and authenticate the phylogenetic relationship between studied taxa of the species of the genera Cassia and Senna in Egypt using the recent tools of ITS barcoding, RAPD analysis and metabolic profiling, in comparing to the traditional taxonomical features. From the cluster analysis of the traditional 27 morphological characters, the studied taxa were categorized into two major clades with an average taxonomic distance of 4.3. The clade I include Cassia fistula, C. renigera, C. javanica L subsp. nodosa and C. roughiia that belongs to series Obolospermae, and C. grandis that belongs to series Grandes. The clade (II) includes Senna surattensis and S. alata at taxonomic level 3.6. The taxonomical description of the studied taxa was confirmed from the molecular analysis of ITS sequences and RAPD analysis. The ITS sequences of the tested plants species C. fistula L, C. grandis MD4, C. javanica subsp. nodosa MD7, C. roxburghii MD5, C. renigera MD5 were deposited at genbank with accession numbers MW367973, MZ960447, MW386305, MW326753 and MW32685, respectively. While, the ITS sequences of the S. surrattensis and S. alata were deposited into genbank accession # MD14 MW367670 and MD20 MW412635, respectively. Thus, from the molecular analysis, two clades were clearly separated into Clade I of Cassia and Clade II of Senna. The cluster I represented by C. fistula, C. renigera, C. roxburghii, and C. javanica sub nodosa, and the cluster II represented by S. alata and S. surattensis. From the PCA of RAPD, a clearly discrimination between the two Taxa was observed revealing the characteristic grouping of Cassia and Senna. The species Senna alata and Senna surattensis were grouped together, but the species of C. renigera, C. javanica, C. roxburghii and C. grandis was grouped on a distinct group. The separation of Cassia and Senna species into two clusters verify the segregation of the genus Cassia L. senso lato into two distinct genera namely Senna P. and Cassia L. The morphological, molecular traits of the studied plants were authenticated from the metabolic profiling by GC-MS analysis. Among the 23 identified metabolites, four compounds namely hexadecanoic acid, methyl ester, 9-Octadecenoic acid (Z)-ethyl ester and Vitamin E were detected with fluctuated concentrations, among C. fistula, C. grandis, C. javanica subsp. nodosa and C. roxburghii. Conclusively, the traditional morphological features, molecular barcoding using ITS sequences, RAPD analysis and metabolic traits by GC-MS analysis, authenticates the taxonomical diversity of the genus Cassia and Senna.

Effects of climate change on the distribution of wild Akebia trifoliata

Understanding the impacts and constraints of climate change on the geographical distribution of wild Akebia trifoliata is crucial for its sustainable management and economic development as a medicinal material or fruit. In this study, according to the first‐hand information obtained from field investigation, the distribution and response to climate change of A. trifoliata were studied by the MaxEnt model and ArcGIS. The genetic diversity and population structure of 21 natural populations of A. trifoliata were studied by simple sequence repeat (SSR) markers. The results showed that the most important bioclimatic variable limiting the distribution of A. trifoliata was the Mean Temperature of Coldest Quarter (bio11). Under the scenarios SSP1‐2.6 and SSP2‐4.5, the suitable area of A. trifoliata in the world will remain stable, and the suitable area will increase significantly under the scenarios of SSP3‐7.0 and SSP5‐8.5. Under the current climate scenario, the suitable growth regions of A. trifoliata in China were 79.9–122.7°E and 21.5–37.5°N. Under the four emission scenarios in the future, the geometric center of the suitable distribution regions of Akebia trifoliata in China will move to the north. The clustering results of 21 populations of A. trifoliata analyzed by SSR markers showed that they had a trend of evolution from south to north.

Limited genetic diversity and high differentiation in Angelica dahurica resulted from domestication: insights to breeding and conservation

BACKGROUND

Angelica dahurica belongs to the Apiaceae family, whose dry root is a famous traditional Chinese medicine named as "Bai zhi". There are two cultivars (A. dahurica cv. 'Hangbaizhi' and A. dahurica cv. 'Qibaizhi'), which have been domesticated for thousands of years. Long term artificial selection has led to great changes in root phenotypes of the two cultivars, and also decreased their adaptability to environment. We proposed hypothesis that the cultivars may have lost some of the genetic diversity found in the wild species and may be highly differentiated from the latter during the domestication process. However, few studies have been carried out on how domestication affected the genetic variation of this species. Here, we accessed the levels of genetic variation and differentiation within and between wild A. dahurica populations and two cultivars using 12 microsatellite markers.

RESULTS

The results revealed that the genetic diversity of the cultivars was much lower than that of wild A. dahurica, and A. dahurica cv. 'Qibaizhi' had lower genetic diversity compared to A. dahurica cv. 'Hangbaizhi'. AMOVA analysis showed significant genetic differentiation between the wild and cultivated A. dahurica populations, and between A. dahurica cv. 'Hangbaizhi' and A. dahurica cv. 'Qibaizhi'. Results from Bayesian, UPGMA, NJ and PcoA clustering analysis indicated that all 15 populations were assigned to two genetic clusters corresponding to the wild and cultivated populations. Bayesian clustering analysis further divided the cultivated populations into two sub-clusters corresponding to the two cultivars.

CONCLUSIONS

Our study suggests that the domestication process is likely the major factor resulting in the loss of genetic diversity in cultivated A. dahurica populations and in significant genetic differentiation from the wild populations due to founder effect and/or artificially directional selections. This large-scale analysis of population genetics could provide valuable information for genetic resources conservation and breeding programs of Angelica dahurica.

Molecular phylogeography and species distribution modelling evidence of 'oceanic' adaptation for Actinidia eriantha with a refugium along the oceanic-continental gradient in a biodiversity hotspot

BACKGROUND

Refugia is considered to be critical for maintaining biodiversity; while discerning the type and pattern of refugia is pivotal for our understanding of evolutionary processes in the context of conservation. Interglacial and glacial refugia have been studied throughout subtropical China. However, studies on refugia along the oceanic-continental gradient have largely been ignored. We used a liana Actinidia eriantha, which occurs across the eastern moist evergreen broad-leaved forests of subtropical China, as a case study to test hypotheses of refugia along the oceanic-continental gradient and 'oceanic' adaptation.

RESULTS

The phylogeographic pattern of A. eriantha was explored using a combination of three cpDNA markers and 38 nuclear microsatellite loci, Species distribution modelling and dispersal corridors analysis. Our data showed intermediate levels of genetic diversity [haplotype diversity (h_T) = 0.498; unbiased expected heterozygosity (UH_E) = 0.510] both at the species and population level. Microsatellite loci revealed five clusters largely corresponding to geographic regions. Coalescent time of cpDNA lineages was dated to the middle Pliocene (ca. 4.03 Ma). Both geographic distance and climate difference have important roles for intraspecific divergence of the species. The Zhejiang-Fujian Hilly Region was demonstrated to be a refugium along the oceanic-continental gradient of the species and fit the 'refugia in refugia' pattern. Species distribution modelling analysis indicated that Precipitation of Coldest Quarter (importance of 44%), Temperature Seasonality (29%) and Mean Temperature of Wettest Quarter (25%) contributed the most to model development. By checking the isolines in the three climate layers, we found that A. eriantha prefer higher precipitation during the coldest quarter, lower seasonal temperature difference and lower mean temperature during the wettest quarter, which correspond to 'oceanic' adaptation. Actinidia eriantha expanded to its western distribution range along the dispersal corridor repeatedly during the glacial periods.

CONCLUSIONS

Overall, our results provide integrated evidence demonstrating that the Zhejiang-Fujian Hilly Region is a refugium along the oceanic-continental gradient of Actinidia eriantha in subtropical China and that speciation is attributed to 'oceanic' adaptation. This study gives a deeper understanding of the refugia in subtropical China and will contribute to the conservation and utilization of kiwifruit wild resources in the context of climate change.

A Change in Conservation Status of Pachyphytum caesium (Crassulaceae), a Threatened Species from Central Mexico Based on Genetic Studies

Simple Summary

Species decline has accelerated in recent decades, with rare species often being the first to go extinct, mainly due to low population sizes. This study worked with Pachyphytum cesium (Crassulaceae), an endemic species from central Mexico threatened by human activities and climate change. To increase our knowledge of the species, we analyze the diversity and genetic structure of all known populations of P. caesium to change their current genetic status and propose conservation strategies. The results indicate that this species presents low to moderate levels of genetic diversity and incipiently structured populations due to low genetic flow. We found that genetic parameters are essential to the conservation strategies and status vulnerability reclassification. Based on our results, we propose reclassifying the conservation status of the species in danger of extinction; hence a series of conservation strategies are provided to prevent its disappearance.

Abstract

Genetic studies in rare species are important to determine their patterns of genetic diversity among populations and implement conservation plans aimed to reduce extinction risks. Pachyphytum caesium is an endemic species with extremely small populations in central Mexico. This work analyzes the diversity and genetic structure of Pachyphytum cesium (Crassulaceae) to change the conservation status and propose conservation strategies. Twelve dominant ISSR markers were used to describe the genetic diversity of all known populations. Additionally, we carried out two simulations to explore how the loss of individuals or the local populations extinction affect the genetics parameters of the species. The average results indicated moderate levels of genetic diversity (number of alleles = 89.7 ± 3.9, expected heterozygosity = 0.212 ± 0.0, and percentage of polymorphic loci = 56.1 ± 3.9), parameters that decreased significantly when simulations were performed in P. caesium. Additionally, a genetic structure of the populations was found with low gene flow (Nm = 0.92). Genetic parameters are negatively affected by the loss of individuals and the local extinction of populations. Based on our results, we propose to reclassify the conservation status of the species in danger of extinction, hence a series of conservation strategies are provided to prevent its disappearance.

Population Genetic Structure Analysis Reveals Significant Genetic Differentiation of the Endemic Species Camellia chekiangoleosa Hu. with a Narrow Geographic Range

In order to protect and utilize the germplasm resource better, it is highly necessary to carry out a study on the genetic diversity of Camellia chekiangoleosa Hu. However, systematic research on population genetics analysis of the species is comparatively rare. Herein, 16 highly variable simple sequence repeat (SSR) markers were used for genetic structure assessment in 12 natural C. chekiangoleosa populations. The genetic diversity of C. chekiangoleosa was low (h = 0.596), within which, central populations (such as Damaoshan (DMS), Sanqingshan (SQS), and Gutianshan (GTS)) at the junction of four main mountain ranges presented high diversity and represented the center of the C. chekiangoleosa diversity distribution; the Hengshan (HS) population in the west showed the lowest diversity, and the diversity of the eastern and coastal populations was intermediate. C. chekiangoleosa exhibited a high level of genetic differentiation, and the variation among populations accounted for approximately 24% of the total variation. The major reasons for this situation are the small population scale and bottleneck effects in some populations (HS and Lingshan (LS)), coupled with inbreeding within the population and low gene flow among populations (Nm = 0.796). To scientifically protect the genetic diversity of C. chekiangoleosa, in situ conservation measures should be implemented for high-diversity populations, while low-diversity populations should be restored by reintroduction.

Dispersal into the Qinghai-Tibet plateau: evidence from the genetic structure and demography of the alpine plant Triosteum pinnatifidum

Triosteum pinnatifidum Maxim., an alpine plant, is traditionally used for several medicinal purposes. Here, both chloroplast DNA sequences and nuclear low copy sequence markers were used to investigate the genetic diversity and population structure of T. pinnatifidum. Materials were collected from thirteen localities in the northeast Qinghai-Tibet Plateau (QTP) and adjacent highlands and advanced analytical toolkits were used to access their origin and range shifts. The results revealed a higher level of population differentiation based on chloroplast DNA (cpDNA) concatenated sequences compared with the nuclear DNA sequences (F _ST = 0.654 for cpDNA, F _ST = 0.398 for AT103), indicating that pollen flow was still extensive in T. pinnatifidum. A decline in haplotype variation was observed from the plateau edge and adjoining highlands toward the platform of the QTP. The hypothesis "dispersal into the QTP," proposing that T. pinnatifidum experienced migration from the plateau edge and adjacent highlands to the platform, was supported. These results were in line with the hypothesis that multiple refugia exist on the plateau edge and adjacent highlands rather than on the plateau platform. Our unimodal mismatch distribution, star-like network supported a recent expansion in T. pinnatifidum.

The Genetic Differences and Structure of Selected Important Populations of the Endangered Taxus baccata in the Czech Republic

Taxus baccata L. (common yew) is an endangered tree species in the Czech Republic. However, its natural occurrence has not been adequately protected in all areas of the country. The aim of this study is to determine whether the yew population in the newly established Mařeničky seed orchard (TS_L) enables mixing with other Czech yew populations. Using a set of nuclear microsatellites, the genetic diversity in the Lužické Mountains (TS_L) and in selected wild-provenance populations from the Czech Republic (Jílovské yews, TS_J; Březinské yews, TS_B, and yews from Moravský Karst, TS_M) was studied, as they could be donor sources for potential translocation activities. We observed that the level of genetic diversity within the four Czech yew units that were investigated was high. An analysis of the molecular variance (AMOVA) showed 7% variation among populations, and the genetic differentiation values were low to moderate (FST = 0.042–0.108). According to a STRUCTURE analysis, high genetic similarity was observed between the TS_L and TS_B units. Our results provide important genetic suggestions on how conservation management can be designed to maximize its success.

Geographical context outweighs habitat disturbance effects in explaining mistletoe population genetic differentiation at a regional scale

Genetic differentiation depends on ecological and evolutionary processes that operate at different spatial and temporal scales. While the geographical context is likely to determine large-scale genetic variation patterns, habitat disturbance events will probably influence small-scale genetic diversity and gene flow patterns. Therefore, the genetic diversity patterns that we observe today result from the combination of both processes, but they are rarely assessed simultaneously. We determined the population structure and genetic diversity of a hemiparasitic mistletoe (Tristerix corymbosus) from the temperate rainforests of southern Chile to determine the effects of geographical context and habitat disturbance at a regional scale and if it is affected by the abundance and occurrence of its seed disperser mutualist (the arboreal marsupial Dromiciops gliroides). We genotyped 359 individuals from 12 populations using single nucleotide polymorphisms, across three different geographical contexts and four disturbance conditions. We also used camera traps to estimate the abundance and occurrence of the seed disperser. Our results suggest that genetic differences among populations are related more to geographical context than to habitat disturbance. However, as disturbance increased, D. gliroides abundance and occurrence decreased, and mistletoe inbreeding index (F_IS ) increased. We also found highly uneven gene flow among study sites. Despite the high levels of disturbance that these temperate rainforests are facing, our results suggest that mistletoe genetic differentiation at a regional scale was more influenced by historical events. However, habitat disturbance can indirectly affect mistletoe population genetic differentiation via the seed dispersal process, which may increase levels of inbreeding.

Do we have infraspecific taxa of Salvia macrosiphon Boiss. (Lamiaceae) in Iran?

BACKGROUND

Salvia macrosiphon is an aromatic perennial species of Lamiaceae family that grows naturally in different parts of Iran. This herb is widely used in folk and modern medicine. Although in Flora Iranica and Flora of Iran, no infraspecific taxonomic rank has been detected for S. macrosiphon, some infraspecific taxonomic ranks have been reported. In the current study, we evaluated the genetic diversity and structure of 11 populations of this species to detect inter and intrapopulation genetic diversity and to survey the possibility of infraspecific taxonomic ranks in this species.

METHODS AND RESULTS

We utilized the modified C-TAB protocol for DNA extraction and amplified the genomes using several SCoT molecular markers. We calculated of genetic diversity and polymorphism parameters using GenAlex 6.4, Geno-Dive ver.2, PopGene, PopART and Structure 2.3.4. The parameters of genetic polymorphism differed between the populations. Moreover, a low rate of gene flow supported a moderate level of population's genetic diversity and differentiation. According to haplotypes network (TCS) analysis, a high level of genetic mutation has occurred among the individuals of some populations leading to high intrapopulation diversity. On the basis of structure analysis and Nei's genetic distance, the examined populations were classified into four genotypes.

CONCLUSIONS

The clustering pattern of the populations in each group was not related to geographical distance or phytogeography. It seems that the wide geographic distribution, a small gene flow rate and the occurrence of a high level of genetic mutation lead to infraspecific genetic differentiation in the species and we suppose some infraspecific ranks exist for it.

Chloroplast and Nuclear Genetic Diversity Explain the Limited Distribution of Endangered and Endemic Thuja sutchuenensis in China

Narrow-ranged species face challenges from natural disasters and human activities, and to address why species distributes only in a limited region is of great significance. Here we investigated the genetic diversity, gene flow, and genetic differentiation in six wild and three cultivated populations of Thuja sutchuenensis, a species that survive only in the Daba mountain chain, using chloroplast simple sequence repeats (cpSSR) and nuclear restriction site-associated DNA sequencing (nRAD-seq). Wild T. sutchuenensis populations were from a common ancestral population at 203 ka, indicating they reached the Daba mountain chain before the start of population contraction at the Last Interglacial (LIG, ∼120-140 ka). T. sutchuenensis populations showed relatively high chloroplast but low nuclear genetic diversity. The genetic differentiation of nRAD-seq in any pairwise comparisons were low, while the cpSSR genetic differentiation values varied with pairwise comparisons of populations. High gene flow and low genetic differentiation resulted in a weak isolation-by-distance effect. The genetic diversity and differentiation of T. sutchuenensis explained its survival in the Daba mountain chain, while its narrow ecological niche from the relatively isolated and unique environment in the "refugia" limited its distribution.

Assessment of Genetic Diversity and Population Structure of the Endangered Astragalus exscapus subsp. transsilvanicus through DNA-Based Molecular Markers

Astragalus exscapus L. subsp. transsilvanicus (Schur) Nyár. (Fabaceae) is a rare plant endemic to the Transylvanian Plateau, represented by 24 identified populations. Limited information on the genetic variation and population structure is available, which obstructs efficient measures for conservation strategy. The present study aimed to analyze the genetic diversity and population structure of eight populations of A. exscapus subsp. transsilvanicus revealed by sequence-related amplified polymorphism (SRAP) markers. A total of 164 bands were amplified, 96.7% of which (159) were polymorphic. Nei’s gene diversity index (He) was estimated to be 0.228 at the population level and 0.272 at the subspecies level. The genetic differentiation among populations (Gst) was 0.165, which indicated a low proportion of total genetic diversity. The analysis of molecular variance (AMOVA) indicated that 17% of the total variation of A. exscapus subsp. transsilvanicus is found among the populations, while 83% was found within the populations. A UPGMA dendrogram, principal coordinate analysis, and the STRUCTURE software grouped the populations into two clusters uncorrelated with the provenience of the 125 individuals, which might be attributed to fragmentation processes, insect pollination, population size, and specific environmental conditions of the habitats.

Phylogeography of Iris loczyi (Iridaceae) in Qinghai-Tibet Plateau revealed by chloroplast DNA and microsatellite markers

Quaternary climate oscillations and complex topography have tremendous effects on current distribution and genetic structure of species, and hence the Qinghai-Tibet Plateau (QTP), the largest plateau in the world, has become a hotspot for many phylogeographic studies. However, little is known about the phylogeographic pattern of herbaceous plants in QTP. Here, we investigate the genetic diversity, population structure and historical dynamics of Iris loczyi, using five chloroplast DNA (cpDNA) fragments and seven microsatellite markers. A total of 15 populations, and 149 individuals were sampled throughout the QTP. High genetic diversity was detected both in cpDNA (H _d = 0.820) and SSR (H _o = 0.689, H _e = 0.699). Ten cpDNA haplotypes and 163 alleles were identified. AMOVA and clustering analyses revealed obvious differentiation between regions. The N _st, G _st and Mantel test showed significant phylogeographic structure of I. loczyi. The neutrality test and mismatch distribution analyses indicated that I. loczyi could not have undergone a historical population expansion, but population XS from the Qilian Mountain area could have experienced a local expansion. Bottleneck analyses indicated that I. loczyi had not experienced bottleneck recently. Based on cpDNA and SSR results, the Qilian Mountain area was inferred as a potential glacial refuge, and the southern Tibet valley was considered as a 'microrefugia' for I. loczyi. These findings provided new insights into the location of glacial refuges for the species distributed in QTP, and supplemented more plant species data for the response of QTP species to the Quaternary climate.

Genetic diversity and structure of Rhododendron meddianum, a plant species with extremely small populations

Rhododendron meddianum is a critically endangered species with important ornamental value and is also a plant species with extremely small populations. In this study, we used double digest restriction-site-associated DNA sequencing (ddRAD) technology to assess the genetic diversity, genetic structure and demographic history of the three extant populations of R. meddianum. Analysis of SNPs indicated that R. meddianum populations have a high genetic diversity (π = 0.0772 ± 0.0024, H _E  = 0.0742 ± 0.002). Both F _ST values (0.1582-0.2388) and AMOVA showed a moderate genetic differentiation among the R. meddianum populations. Meanwhile, STRUCTURE, PCoA and NJ trees indicated that the R. meddianum samples were clustered into three distinct genetic groups. Using the stairway plot, we found that R. meddianum underwent a population bottleneck about 70,000 years ago. Furthermore, demographic models of R. meddianum and its relative, Rhododendron cyanocarpum, revealed that these species diverged about 3.05 (2.21-5.03) million years ago. This divergence may have been caused by environmental changes that occurred after the late Pliocene, e.g., the Asian winter monsoon intensified, leading to a drier climate. Based on these findings, we recommend that R. meddianum be conserved through in situ, ex situ approaches and that its seeds be collected for germplasm.

Life in the desert: The impact of geographic and environmental gradients on genetic diversity and population structure of Ivesia webberi

For range-restricted species with disjunct populations, it is critical to characterize population genetic structure, gene flow, and factors that influence functional connectivity among populations in order to design effective conservation programs. In this study, we genotyped 314 individuals from 16 extant populations of Ivesia webberi, a United States federally threatened Great Basin Desert using six microsatellite loci. We assessed the effects of Euclidean distance, landscape features, and ecological dissimilarity on the pairwise genetic distance of the sampled populations, while also testing for a potential relationship between I. webberi genetic diversity and diversity in the vegetative communities. The results show low levels of genetic diversity overall (H e = 0.200-0.441; H o = 0.192-0.605) and high genetic differentiation among populations. Genetic diversity was structured along a geographic gradient, congruent with patterns of isolation by distance. Populations near the species' range core have relatively high genetic diversity, supporting in part a central-marginal pattern, while also showing some evidence for a metapopulation dynamic. Peripheral populations have lower genetic diversity, significantly higher genetic distances, and higher relatedness. Genotype cluster admixture results suggest a complex dispersal pattern among populations with dispersal direction and distance varying on the landscape. Pairwise genetic distance strongly correlates with elevation, actual evapotranspiration, and summer seasonal precipitation, indicating a role for isolation by environment, which the observed phenological mismatches among the populations also support. The significant correlation between pairwise genetic distance and floristic dissimilarity in the germinated soil seed bank suggests that annual regeneration in the plant communities contribute to the maintenance of genetic diversity in I. webberi.

Genetic diversity of wild and cultivated Coffea canephora in northeastern DR Congo and the implications for conservation

PREMISE

Many cultivated coffee varieties descend from Coffea canephora, commonly known as Robusta coffee. The Congo Basin has a century-long history of Robusta coffee cultivation and breeding, and is hypothesized to be the region of origin of many of the cultivated Robusta varieties. Since little is known about the genetic composition of C. canephora in this region, we assessed the genetic diversity of wild and cultivated C. canephora shrubs in the Democratic Republic of the Congo.

METHODS

Using 18 microsatellite markers, we studied the genetic composition of wild and backyard-grown C. canephora shrubs in the Tshopo and Ituri provinces and multiple accessions from the INERA Yangambi Coffee Collection. We assessed genetic clustering patterns, genetic diversity, and genetic differentiation between populations.

RESULTS

Genetic differentiation was relatively strong between wild and cultivated C. canephora shrubs, and both gene pools harbored multiple unique alleles. Strong genetic differentiation was also observed between wild populations. The level of genetic diversity in wild populations was similar to that of the INERA Yangambi Coffee Collection, but local wild genotypes were mostly missing from that collection. Shrubs grown in the backyards were genetically similar to the breeding material from INERA Yangambi.

CONCLUSIONS

Most C. canephora that is grown in local backyards originated from INERA breeding programs, while a few shrubs were obtained directly from surrounding forests. The INERA Yangambi Coffee Collection could benefit from an enrichment with local wild genotypes to increase the genetic resources available for breeding purposes and to support ex situ conservation.

SSR-Sequencing Reveals the Inter- and Intraspecific Genetic Variation and Phylogenetic Relationships among an Extensive Collection of Radish (Raphanus) Germplasm Resources

Simple Summary

Raphanus is an important genus of Brassicaceae and has undergone a lengthy evolutionary process. However, the inter- and intraspecific phylogenetic relationships and genetic diversity are not well understood. To elucidate these issues, we SSR-sequenced 939 wild, semi-wild and cultivated accessions, and discovered that Europe was the origin center of radishes with diverse European wild radishes, and Europe, South Asia and East Asia might be three independent domestication centers. There was considerable genetic differentiation within European cultivated radishes. European primitive cultivated radish exhibited gene flow with black radish/oil radish and rat-tail radish. Among Asian cultivated radishes, rat-tail radish was a sister to the clade of Chines big radish (including Japanese wild radish), suggesting that they may share the most recent common ancestry. Japanese wild radish had strong gene exchange with Japanese/Korea big radish, oil radish and rat-tail radish. American wild radish developed from natural hybridization between European wild radish and European small radish. All these demonstrated that European primitive cultivated type, American wild radish and Japanese wild radish might have played indispensable roles in radish evolution. Our study provides new perspectives into the origin, evolution and genetic diversity of Raphanus and facilitates the conservation and exploitation of radish germplasm resources.

Abstract

Raphanus has undergone a lengthy evolutionary process and has rich diversity. However, the inter- and intraspecific phylogenetic relationships and genetic diversity of this genus are not well understood. Through SSR-sequencing and multi-analysis of 939 wild, semi-wild and cultivated accessions, we discovered that the European wild radish (EWR) population is separated from cultivated radishes and has a higher genetic diversity. Frequent intraspecific genetic exchanges occurred in the whole cultivated radish (WCR) population; there was considerable genetic differentiation within the European cultivated radish (ECR) population, which could drive radish diversity formation. Among the ECR subpopulations, European primitive cultivated radishes (EPCRs) with higher genetic diversity are most closely related to the EWR population and exhibit a gene flow with rat-tail radishes (RTRs) and black radishes (BRs)/oil radishes (ORs). Among Asian cultivated radishes (ACRs), Chinese big radishes (CBRs) with a relatively high diversity are furthest from the EWR population, and most Japanese/Korean big radishes (JKBRs) are close to CBR accessions, except for a few old Japanese landraces that are closer to the EPCR. The CBR and JKBR accessions are independent of RTR accessions; however, phylogenetic analysis indicates that the RTR is sister to the clade of CBR (including JWR), which suggests that the RTR may share the most recent common ancestry with CBRs and JWRs. In addition, Japanese wild radishes (JWRs), (namely, R. sativus forma raphanistroides) are mainly scattered between CBRs and EPCRs in PCoA analysis. Moreover, JWRs have a strong gene exchange with the JKBR, OR and RTR subpopulations. American wild radishes (AWRs) are closely related to European wild and cultivated radishes, and have a gene flow with European small radishes (ESRs), suggesting that the AWR developed from natural hybridization between the EWR and the ESR. Overall, this demonstrates that Europe was the origin center of the radish, and that Europe, South Asia and East Asia appear to have been three independent domestication centers. The EPCR, AWR and JWR, as semi-wild populations, might have played indispensable transitional roles in radish evolution. Our study provides new perspectives into the origin, evolution and genetic diversity of Raphanus and facilitates the conservation and exploitation of radish germplasm resources.

Genetic variability and population structure of Ethiopian chickpea (Cicer arietinum L.) germplasm

Evaluation of the genetic diversity and an understanding of the genetic structure and relationships of chickpea genotypes are valuable to design efficient germplasm conservation strategies and crop breeding programs. Information is limited, in these regards, for Ethiopian chickpea germplasms. Therefore, the present study was carried out to estimate the genetic diversity, population structure, and relationships of 152 chickpea genotypes using simple sequence repeats (SSR) markers. Twenty three SSR markers exhibited polymorphism producing a total of 133 alleles, with a mean of 5.8 alleles per locus. Analyses utilizing various genetic-based statistics included pairwise population Nei’s genetic distance, heterozygosity, Shannon’s information index, polymorphic information content, and percent polymorphism. These analyses exemplified the existence of high genetic variation within and among chickpea genotypes. The 152 genotypes were divided into two major clusters based on Nei’s genetic distances. The exotic genotypes were grouped in one cluster exclusively showing that these genotypes are distinct to Ethiopian genotypes, while the patterns of clustering of Ethiopian chickpea genotypes based on their geographic region were not consistent because of the seed exchange across regions. Model-based population structure clustering identified two discrete populations. These finding provides useful insight for chickpea collections and ex-situ conservation and national breeding programs for widening the genetic base of chickpea.

Distinct Genetic Structure Reflects Ploidy Level Differentiation in Newly Discovered, Extremely Small Populations of Xanthocyparis vietnamensis from Southwestern China

Population genetic assessment is crucial for the conservation and management of threatened species. Xanthocyparis vietnamensis is an endangered species that is currently restricted to karst mountains in southwestern China and Vietnam. This rare conifer was first recorded in 2002 from northern Vietnam and then in 2013 from Guangxi, China, yet nothing is known about its genetic diversity nor ploidy level variation, although previous cytological study suggest that Vietnamese populations are tetraploids. There have been about 45 individuals found to date in Guangxi, China. Here, we genotyped 33 X. vietnamensis individuals using 20 newly developed, polymorphic microsatellite loci, to assess the genetic variability of its extremely small populations. The genetic diversity of X. vietnamensis (H_E = 0.511) was lower than that of two other heliophile species, Calocedrus macrolepis and Fokienia hodginsii, which have similar distribution ranges. This is consistent with the signature of a genetic bottleneck detected in X. vietnamensis. Although the population genetic differentiation coefficient across loci is moderate (F_ST = 0.125), STRUCTURE analysis revealed two distinct genetic clusters, namely the northern and southern population groups; DAPC analysis grouped the southern populations together in one cluster separate from the northern populations; AMOVA analysis detected a significant genetic differentiation between the two population groups (F_RT = 0.089, p < 0.05), and BARRIER analysis detected a genetic barrier between them. Moreover, we detected differentiation in ploidy level between northern and southern populations, sampled individuals from the former and the later are all diploid and tetraploid cytotypes with mean genome sizes of 26.08 and 48.02 pg/2C, respectively. We deduced that heterogeneous geomorphology and historical events (e.g., human deforestation, Quaternary climate oscillations) may have contributed to population fragmentation and small population size in X. vietnamensis. Considering both genetic and ploidy level differentiation, we propose that two different management units (northern and southern) should be considered and a combination of in situ and ex situ conservation measures should be employed to preserve populations of this endangered species in southwestern China in the light of our findings.

Patterns of genetic diversity, spatial genetic structure and gene flow in Campomanesia xanthocarpa: insights from SSR markers of different genomic origins

Campomanesia xanthocarpa (Mart.) O. Berg is a South American fruit tree species with important ecological and medicinal properties, which remnants are currently found mainly in isolated forest fragments. In this study, SSR markers from three different genomic origins (gene-linked, nuclear neutral, and organellar) were used to evaluate the patterns of genetic diversity, fine-scale spatial genetic structure and historical gene flow in fragmented forest formations of C. xanthocarpa from the Atlantic Forest in southern Brazil. Our results show that the forest fragments present moderate to high levels of genetic diversity in comparison to species presenting similar life traits, although a trend opposite to expected was observed concerning gene-linked and neutral SSR markers. The fine-scale spatial genetic structure revealed different patterns in short and large distance classes, with a distinct influence of gene-linked and neutral markers in driving the genetic structure in each distance class. The presence of an isolation-by-adaptation pattern implies the need for maintenance of the current remnants to assure the conservation of the private alleles. Finally, as the genetic diversity is found predominantly within forest fragments, programs of seed collection and/or genetic rescue should prioritize a larger number of individuals within each fragment, to increase the sampled diversity.

Genetic diversity in North American Cercis Canadensis reveals an ancient population bottleneck that originated after the last glacial maximum

Understanding of the present-day genetic diversity, population structure, and evolutionary history of tree species can inform resource management and conservation activities, including response to pressures presented by a changing climate. Cercis canadensis (Eastern Redbud) is an economically valuable understory tree species native to the United States (U.S.) that is also important for forest ecosystem and wildlife health. Here, we document and explain the population genetics and evolutionary history of this deciduous tree species across its distributed range. In this study, we used twelve microsatellite markers to investigate 691 wild-type trees sampled at 74 collection sites from 23 Eastern U.S. states. High genetic diversity and limited gene flow were revealed in wild, natural stands of C. canadensis with populations that are explained by two major genetic clusters. These findings indicate that an ancient population bottleneck occurred coinciding with the last glacial maximum (LGM) in North America. The structure in current populations likely originated from an ancient population in the eastern U.S. that survived LGM and then later diverged into two contemporary clusters. Data suggests that populations have expanded since the last glaciation event from one into several post-glacial refugia that now occupy this species’ current geographic range. Our enhanced understanding benchmarks the genetic variation preserved within this species and can direct future efforts in conservation, and resource utilization of adaptively resilient populations that present the greatest genetic and structural diversity.

How to achieve a higher selection plateau in forest tree breeding? Fostering heterozygote × homozygote relationships in optimal contribution selection in the case study of Populus nigra

In breeding, optimal contribution selection (OCS) is one of the most effective strategies to balance short- and long-term genetic responses, by maximizing genetic gain and minimizing global coancestry. Considering genetic diversity in the selection dynamic-through coancestry-is undoubtedly the reason for the success of OCS, as it avoids preliminary loss of favorable alleles. Originally formulated with the pedigree relationship matrix, global coancestry can nowadays be assessed with one of the possible formulations of the realized genomic relationship matrix. Most formulations were optimized for genomic evaluation, but few for the management of coancestry. We introduce here an alternative formulation specifically developed for genomic OCS (GOCS), intended to better control heterozygous loci, and thus better account for Mendelian sampling. We simulated a multigeneration breeding program with mate allocation and under GOCS for twenty generations, solved with quadratic programming. With the case study of Populus nigra, we have shown that, although the dynamic was mainly determined by the trade-off between genetic gain and genetic diversity, better formulations of the genomic relationship matrix, especially those fostering individuals carrying multiple heterozygous loci, can lead to better short-term genetic gain and a higher selection plateau.

Deciphering allelic variability and population structure in buckwheat: An analogy between the efficiency of ISSR and SSR markers

Food and nutritional security continue to be the issues of concern in developing countries like ours. Exploring the reservoir of high potential unexplored genetic resources could address the world's food and nutritional insecurity. The availability of diverse data and the population structure of any crop germplasm is a valuable genetic resource for discovering genes that can help achieve food and nutritional stability. We used seven ISSR and seven SSR markers to investigate diversity among 63 buckwheat genotypes, including landraces from India's northwestern Himalayas. Various parameters such as percent polymorphism, PIC, resolving power, and marker index was used to evaluate the inequitable efficacy of these markers. We foundthat both marker systems are effective in detecting polymorphism in buckwheat germplasm. Seven ISSRs produced 55 polymorphic bands, while seven SSRs produced 32bands. When compared to ISSRs, SSRs had a greater average PIC value (0.43) than that of (0.36). ISSRs, on the other hand, had a resolving power of (4.38) compared to (1.42) for SSRs. The hierarchical cluster analysis dendrogram divided genotypes into three major clusters. We found that both marker systems were equally accurate in grouping buckwheat genotypes according to their geographical origins. Using 7 ISSR and 7 SSR markers, the model-based STRUCTURE analysis established a population with two sub-populations that correspond to species-based groupings. Within the population, there was a high level of genetic diversity. These results have consequences for both buckwheat breeding and conservation efforts.

History or demography? Determining the drivers of genetic variation in North American plants

Understanding the impact of historical and demographic processes on genetic variation is essential for devising conservation strategies and predicting responses to climate change. Recolonization after Pleistocene glaciations is expected to leave distinct genetic signatures, characterised by lower genetic diversity in previously glaciated regions. Populations’ positions within species ranges also shape genetic variation, following the central‐marginal paradigm dictating that peripheral populations are depauperate, sparse and isolated. However, the general applicability of these patterns and relative importance of historical and demographic factors remains unknown. Here, we analysed the distribution of genetic variation in 91 native species of North American plants by coupling microsatellite data and species distribution modelling. We tested the contributions of historical climatic shifts and the central‐marginal hypothesis on genetic diversity and structure on the whole data set and across subsets based on taxonomic groups and growth forms. Decreased diversity was found with increased distance from potential glacial refugia, coinciding with the expected make‐up of postglacially colonised localities. At the range periphery, lower genetic diversity, higher inbreeding levels and genetic differentiation were reported, following the assumptions of the central‐marginal hypothesis. History and demography were found to have approximately equal importance in shaping genetic variation.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Clone worth? Genetic diversity in obligate apomictic Miconia albicans (Melastomataceae)

Apomixis is the asexual production of seeds by plants and, in theory, would render low genetic diversity and even clonal lineages. However, recent studies have shown otherwise, although is not always clear where the genetic diversity of obligate apomicts comes from. We evaluated the genetic diversity among sister seedlings of M. albicans, an obligate apomictic species in Cerrado, Neotropical Savannas in Central Brazil. A total of 50 seedlings from five individuals were analysed using ISSR primers. We obtained 107 fragments, all with good resolution, consistently observed and replicable. The percentage of polymorphic loci ranged from 28.04% to 33.64% and Shannon's information index (I) averaged 0.173. The expected heterozygosity (He) averaged 0.117, similar to the observed for populations of M. albicans and other selfed species. Only two seedlings showed the same genotype (possible clones), but most differed at least for five loci. Most of variance was among progenies (62%), but we found that 38% was within progenies. Genetic distances separated the progenies in two groups, and analogous analyses between individuals reconstructed the original progenies clustering. The results confirmed a relatively high genetic diversity among sister seedling of this obligatory apomictic plant and clones were rare. This diversity can be generated during development, probably by restitutional meiosis or other recombination processes. These differences may accumulate into lineages and populations well adapted to heterogenous Cerrado environment.

Invasive populations of Spiraea tomentosa (Rosaceae) are genetically diverse but decline during succession in forest habitats

Population genetic and ecological data may help to control invasive plants, which are considered a major threat to natural habitats. In contrast to expected bottleneck events, genetic diversity of such invasive populations may be high due to extensive propagule pressure or admixture. The ecological impact of invasive species has been broadly evaluated in the field; however, long-term studies on the fate of invasive plants are scarce. We analysed genetic diversity and structure in invasive Spiraea tomentosa populations in eastern Germany and western Poland using Amplified Fragment Length Polymorphism. Potential hybridization between co-occurring diploid Sp. tomentosa and tetraploid Sp. douglasii was investigated using Flow Cytometry. The genetic analyses were complemented by data from a 13-year vegetation study in an area invaded by these Spiraea species. We found no evidence for hybridization between Spiraea species. In populations of Sp. tomentosa both genetic diversity (He = 0.26) and genetic structure (ΦPT = 0.27) were high and comparable to other outcrossing woody plants. Low levels of clonality, presence of seedlings and new patches in sites that had been colonized over the last 13 years imply that populations spread via sexual reproduction. In all habitat types, native species diversity declined following Sp. tomentosa invasion. However, detailed aerial mapping of a forest reserve with ongoing succession revealed that Spiraea spp. populations have declined over a 10-year period. Despite its potential for dispersal and negative effects on native plant communities, invasive Spiraea populations may be controlled by increasing canopy cover in forest habitats.

Morphological, phytochemical and genetic diversity of threatened Polygonatum verticillatum (L.) All. populations of different altitudes and habitat types in Himalayan region

Polygonatum verticillatum (L.) All. is an important medicinal herb that belongs to the family Asparagaceae. The rhizome of the species is used in Chyavanprash preparation and several other ayurvedic formulations. Numerous active constituents like saponins, alkaloids, phytohormones, flavonoids, antioxidants, lysine, serine, aspartic acid, diosgenin, β-sitosterol, etc. have been reported from this species. In this study, morphological, phytochemical, antioxidant and genetic variations of 11 distant populations of P. verticillatum were measured. Considerably (P < 0.05) higher variations were recorded among different populations of P. verticillatum using morphological, phytochemical and genetic diversity parameters. AGFW (above ground fresh weights); flavonols, FRAP (Ferric ion reducing antioxidant power) and NO (Nitric Oxide scavenging activity) were recorded maximum in Kafni population. Similarly, a significantly higher above and below ground dry weight was recorded in Mayawati and Surmoli populations respectively. Maximum phenolic content, tannins, and DPPH (2,2-diphenyl-1-picrylhydrazyl) activity were recorded in Milam population. A total of 165 individuals from 11 populations were assessed for genetic diversity using inter-simple sequence repeats (ISSR) marker. High genetic diversity (He = 0.35) was recorded in Himkhola and Surmoli populations while it was observed minimum (0.28) in the Mayawati population. Altitude showed a significant positive correlation with tannins (r = 0.674; P < 005) and DPPH (r = 0.820; P < 0.01). Phenol content exhibited a considerably positive relationship with He (r = 0.606; P < 0.05) and BGFW (r = 0.620; P < 0.05), flavonol displayed a positive correlation with Pp% (r = 0.606; P < 0.05). The population structure of P. verticillatum, exhibited that the optimal value of the K was 3 for its populations as determined by the ΔK statistic structure. Among populations, the amount of gene flow is higher (Nm = 1.717) among all sites. Hence, it can be concluded that P. verticillatum populations possess considerable variability in the collected populations. Likewise, the populations from Kafni, Satbunga and Himkhola with higher morphological, phytochemicals and genetic variability were prioritized and therefore recommended for cultivation and mass multiplication to meet the industrial demand for target species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01044-9.

Predicting stability of DNA bulge at mononucleotide microsatellite

Mononucleotide microsatellites are clinically and forensically crucial DNA sequences due to their high mutability and abundance in the human genome. As a mutagenic intermediate of an indel in a microsatellite and a consequence of probe hybridization after such mutagenesis, a bulge with structural degeneracy sliding within a microsatellite is formed. Stability of such dynamic bulges, however, is still poorly understood despite their critical role in cancer genomics and neurological disease studies. In this paper, we have built a model that predicts the thermodynamics of a sliding bulge at a microsatellite. We first identified 40 common bulge states that can be assembled into any sliding bulges, and then characterized them with toehold exchange energy measurement and the partition function. Our model, which is the first to predict the free energy of sliding bulges with more than three repeats, can infer the stability penalty of a sliding bulge of any sequence and length with a median prediction error of 0.22 kcal/mol. Patterns from the prediction clearly explain landscapes of microsatellites observed in the literature, such as higher mutation rates of longer microsatellites and C/G microsatellites.

High Genetic Diversity and Low Population Differentiation in Wild Hop (Humulus lupulus L.) from Croatia

Hop (Humulus lupulus L.) is used in the brewing industry as a source of compounds responsible for the bitterness, aroma, and preservative properties of beer. In this study, we used microsatellite markers to investigate genetic diversity and genetic differentiation of wild hop populations sampled in the northwestern part of Croatia. Analysis of 12 microsatellite loci revealed high diversity and weak population differentiation among wild hop populations. A total of 152 alleles were determined with an average of 12.67 alleles per locus. Observed heterozygosity ranged from 0.689 to 0.839 (average 0.767) and expected heterozygosity ranged from 0.725 to 0.789 (average 0.760). A total of 38 private alleles were detected. The data suggest that H. lupulus populations are not affected by recent bottlenecks. The degree of genetic differentiation among populations was low and not significant for most pairwise FST values, except for the pair of geographically most distant populations. The results did not indicate the existence of genetic structure among the sampled populations. The high genetic diversity and low differentiation among populations, combined with the absence of isolation by distance, indicate the existence of substantial gene flow among wild hop populations. Therefore, extensive sampling per population is clearly required to assess the genetic diversity of hop populations. Sampling strategies involving sampling across a large number of localities represented by only a few samples could lead to erroneous conclusions.

Population genetic structure and connectivity of a riparian selfing herb Caulokaempferia coenobialis at a fine-scale geographic level in subtropical monsoon forest

BACKGROUND

Rivers and streams facilitate movement of individuals and their genes across the landscape and are generally recognized as dispersal corridors for riparian plants. Nevertheless, some authors have reported directly contrasting results, which may be attributed to a complex mixture of factors, such as the mating system and dispersal mechanisms of propagules (seed and pollen), that make it difficult to predict the genetic diversity and population structure of riparian species. Here, we investigated a riparian self-fertilizing herb Caulokaempferia coenobialis, which does not use anemochory or zoochory for seed dispersal; such studies could contribute to an improved understanding of the effect of rivers or streams on population genetic diversity and structure in riparian plants. Using polymorphic ISSR and cpDNA loci, we studied the effect at a microgeographic scale of different stream systems (a linear stream, a dendritic stream, and complex transverse hydrological system) in subtropical monsoon forest on the genetic structure and connectivity of C. coenobialis populations across Dinghu Mountain (DH) and Nankun Mountain (NK).

RESULTS

The results indicate that the most recent haplotypes (DH: H7, H8; NK: h6, h7, h11, h12) are not shared among local populations of C. coenobialis within each stream system. Furthermore, downstream local populations do not accumulate genetic diversity, whether in the linear streamside local populations across DH (H: 0.091 vs 0.136) or the dendritic streamside local populations across NK (H: 0.079 vs 0.112, 0.110). Our results show that the connectivity of local C. coenobialis populations across DH and NK can be attributed to historical gene flows, resulting in a lack of spatial genetic structure, despite self-fertilization. Selfing C. coenobialis can maintain high genetic diversity (H = 0.251; I = 0.382) through genetic differentiation (GST = 0.5915; FST = 0.663), which is intensified by local adaptation and neutral mutation and/or genetic drift in local populations at a microgeographic scale.

CONCLUSION

We suggest that streams are not acting as corridors for dispersal of C. coenobialis, and conservation strategies for maintaining genetic diversity of selfing species should be focused on the protection of all habitat types, especially isolated fragments in ecosystem processes.

Assessment of genetic diversity in Amomum tsao-ko Crevost & Lemarié, an important medicine food homologous crop from Southwest China using SRAP and ISSR markers

Amomum tsao-ko Crevost & Lemarié is an important crop that has been widely used in traditional Chinese medicine and daily diets for a long time. In this study, the genetic diversity and relationships of eight cultivated populations of A. tsao-ko grown in Southwest China were examined using sequence-related amplified polymorphism (SRAP) and inter-simple sequence repeat (ISSR) markers. The results showed that 139 (99.29%) of 140 and 185 (99.46%) of 186 bands were polymorphic by SRAP and ISSR primers amplification, respectively. The polymorphic information content of detected bands were 0.270 (SRAP) and 0.232 (ISSR), respectively. The average Nei's gene diversity (H = 0.217) and Shannon's information index (I = 0.348) at the species level generated by SRAP primer were higher than those by ISSR analysis (H = 0.158, I = 0.272). Genetic differentiation coefficients and molecular variance analysis (AMOVA) indicated that the genetic variance of A. tsao-ko mainly occurred within populations rather than among populations. The high genetic identity among populations was revealed by SRAP (0.937) and ISSR (0.963). Using UPGMA cluster analysis, principal coordinate analysis, and population structure analysis, the accessions were categorized into two major groups. Overall, results obtained here will be useful for A. tsao-ko germplasm characterization, conservation, and utilization.

Isolation of Metrosideros ('Ohi'a) Taxa on O'ahu Increases with Elevation and Extreme Environments

Species radiations should be facilitated by short generation times and limited dispersal among discontinuous populations. Hawaii's hyper-diverse, landscape-dominant tree, Metrosideros, is unique among the islands' radiations for its massive populations that occur continuously over space and time within islands, its exceptional capacity for gene flow by both pollen and seed, and its extended life span (ca. >650 years). Metrosideros shows the greatest phenotypic and microsatellite DNA diversity on O'ahu, where taxa occur in tight sympatry or parapatry in mesic and montane wet forest on 2 volcanoes. We document the nonrandom distributions of 12 taxa (including unnamed morphotypes) along elevation gradients, measure phenotypes of ~6-year-old common-garden plants of 8 taxa to verify heritability of phenotypes, and examine genotypes of 476 wild adults at 9 microsatellite loci to compare the strengths of isolation across taxa, volcanoes, and distance. All 8 taxa retained their diagnostic phenotypes in the common garden. Populations were isolated by taxon to a range of degrees (pairwise FST between taxa: 0.004-0.267), and there was no pattern of isolation by distance or by elevation; however, significant isolation between volcanoes was observed within monotypic species, suggesting limited gene flow between volcanoes. Among the infraspecific taxa of Metrosideros polymorpha, genetic diversity and isolation significantly decreased and increased, respectively, with elevation. Overall, 5 of the 6 most isolated taxa were associated with highest elevations or otherwise extreme environments. These findings suggest a principal role for selection in the origin and maintenance of the exceptional diversity that occurs within continuous Metrosideros stands on O'ahu.

Contrasting effects of local environment and grazing pressure on the genetic diversity and structure of Artemisia frigida

Drylands count among the most globally extensive biomes, and while many desert and dry rangeland ecosystems are under threat, genetic structures of dryland species are still rarely studied. Artemisia frigida is one of the most widely distributed plant species in the temperate rangelands of Eurasia and North America, and it also dominates in many habitats of Mongolia due to its tolerance to low temperatures, drought and disturbance. Local environmental conditions and grazing pressure can influence species performance and affect spatial patterns of genetic diversity in contrasting ways, and our study set out to evaluate such effects on the genetic diversity and structure of A. frigida. We first developed new species-specific Simple Sequence Repeats (SSRs) markers using whole genome sequencing. We then analysed 11 populations of A. frigida that had been sampled along a large climatic gradient in Mongolia, which were sub-structured according to three levels of grazing intensity. Estimates of genetic diversity at the population level were high (HO = 0.56, HE = 0.73) and tended to increase with higher precipitation and soil nutrient availability. Grazing had no effect on genetic diversity, however, a high number of grazing-specific indicator alleles was found at grazed sites. Genetic differentiation among populations was extremely low (global GST = 0.034). Analysis of Molecular Variance revealed 5% variance between populations along the climatic gradient, with 3% of the variance being partitioned among different grazing intensity levels. We found no relationship between geographic and genetic distances, and thus no isolation by distance in this widely distributed species. The relatively low genetic structuring suggests that considerable gene flow exists among A. frigida populations across the rangelands of Mongolia, in spite of the pervasive grazing in the region.

Genetic diversity and structure of Saussurea polylepis (Asteraceae) on continental islands of Korea: Implications for conservation strategies and management

Saussurea polylepis Nakai is an herbaceous perennial endemic to Korea and is highly restricted to several continental islands in the southwestern part of the Korean Peninsula. Given its very narrow geographical distribution, it is more vulnerable to anthropogenic activities and global climate changes than more widely distributed species. Despite the need for comprehensive genetic information for conservation and management, no such population genetic studies of S. polylepis have been conducted. In this study, genetic diversity and population structure were evaluated for 97 individuals from 5 populations (Gwanmaedo, Gageodo, Hongdo, Heusando, and Uido) using 19 polymorphic microsatellites. The populations were separated by a distance of 20–90 km. We found moderate levels of genetic diversity in S. polylepis (Ho = 0.42, He = 0.43). This may be due to long lifespans, outcrossing, and gene flow, despite its narrow range. High levels of gene flow (Nm = 1.76, mean Fst = 0.09), especially from wind-dispersed seeds, would contribute to low levels of genetic differentiation among populations. However, the small population size and reduced number of individuals in the reproductive phase of S. polylepis can be a major threat leading to inbreeding depression and genetic diversity loss. Bayesian cluster analysis revealed three significant structures at K = 3, consistent with DAPC and UPGMA. It is thought that sea level rise after the last glacial maximum may have acted as a geographical barrier, limiting the gene flow that would lead to distinct population structures. We proposed the Heuksando population, which is the largest island inhabited by S. polylepis, as a source population because of its large population size and high genetic diversity. Four management units (Gwanmaedo, Gageodo, Hongdo-Heuksando, and Uido) were suggested for conservation considering population size, genetic diversity, population structure, unique alleles, and geographical location (e.g., proximity).

Growing Out of the Tropical Forests: Gene Flow of Native Mesoamerican Trees Among Forest and Mayan Homegardens

This work aimed to evaluate domestication effects on the genetic structure of two dioecious speciesBrosimum alicastrumSw. (Moraceae) andSpondias purpureaL. (Anacardiaceae), and a heterostylous oneCordia dodecandraA. DC. (Cordiaceae), growing in remnant forests and homegardens within two climatic regions of the Peninsula of Yucatan. The trees ofB. alicastrumandC. dodecandraare propagated by seeds in both population types, while those ofS. purpureaare propagated asexually in the homegardens. ISSRs genetic markers were amplified from foliar tissue of 18 to 21 plants per population type/region combination for each species. Genetic diversity, genetic differentiation, and genetic structure estimators were obtained and compared among species at the regional and population level. We found higher polymorphism (37.5–41), but lower private alleles (4–4.4) and similar heterozygosity (0.1–0.12) in the species with sexual reproduction compared toS. purpurea(34, 8, and 0.11, respectively). Genetic diversity inB. alicastrumpopulations varied with the region; inC. dodecandra, to the population type; and inS. purpurea, to both the population type and the region. Unrestricted gene flow among regions was suggested by low ΦRTinC. dodecandraandS. purpurea(−0.006 and 0.002) but not forB. alicastrum(0.1). Gene flow between populations within the regions for the sexually reproducing species was suggested by lower θII(0.005–0.07 and 0.008–0.1) estimates than those ofS. purpurea(0.09 and 0.13). Even though the lowest paired FST(0.002–0.05) and ΦST(0.002–0.12) values were found between the northeastern forest and homegarden populations for the three species, the dendrogram, Bayesian assignment, and K-Means analyses suggest that the least differentiated populations are southwestern forest and homegarden populations ofB. alicastrumandS. purpurea, and the southwestern forest and northeastern homegarden ofC. dodecandra. The sexual reproduction, biotic interactions, and extensive management ofB. alicastrumandC. dodecandrain the agroforestry and the urban systems may contribute to connectivity between wild and domesticated populations, while inS. purpureathis connectivity is interrupted by the clonal propagation of the species in the homegardens.

Stevia rebaudiana germplasm characterization using microsatellite markers and steviol glycosides quantification by HPLC

Stevia rebaudiana Bertoni, Asteraceae, is an herbaceous perennial plant native to Paraguay. This species is considered since ancient times a medicinal plant with important bioactive compounds and pharmacologic and food properties, namely diterpenes glycosides. The high natural sweetener potential stevioside and rebaudioside A produced by S. rebaudiana plants are suitable sucrose substitutes, and their obtention is influenced by environmental, phytosociological, and genetic factors. The plants’ genetic profile and sweet potential depiction are needed for suitable plant selection for improvement and deployment. Thirty-one S. rebaudiana accessions grown in the same plot where leaves samples were collected in early 2019, were genotyped using six microsatellite markers, including two steviol glycosides biosynthesis functionally involved markers. Additionally, an aqueous extract of each sample was obtained in a water bath and purified by SPE for stevioside and rebaudioside A quantification by normal phase HPLC. Stevioside and rebaudioside A contents varied between 0.53–7.36% (w w⁻¹) and 0.37–3.60% (w w⁻¹), respectively. Two genotypes displayed interesting ratios of rebaudioside A/stevioside (number 3 and 33). The level of genetic similarity between genotypes was tested through a pairwise similarity coefficient, and two groups of individuals had the same fingerprinting. Strong relatedness was found within genotypes, possibly due to cloning, thus, influx of new germplasm ought to be made to prevent mating between relatives, and for further selection and genetic improvement.

Conservation Genetics of Four Critically Endangered Greek Endemic Plants: A Preliminary Assessment

The Mediterranean basin constitutes one of the largest global biodiversity hotspots, hosting more than 11,000 endemic plants, and it is recognised as an area with a high proportion of threatened taxa. Nevertheless, only a tiny fraction of the threatened Mediterranean endemics have their genetic diversity assessed, and we are unaware if and how climate change might impact their conservation status. This is even more pronounced in Eastern Mediterranean countries with a rich endemic flora, such as Greece, which hosts a large portion of the plant taxa assessed at the European level under the IUCN criteria. Using inter simple sequence repeats (ISSR) markers and species distribution models, we analysed the genetic diversity and investigated the impacts of climate change on four critically endangered and extremely narrow and rare Greek island endemic plants, namely Aethionema retsina, Allium iatrouinum, Convolvulus argyrothamnos, and Saponaria jagelii. All four species are facing intense anthropogenic threats and display moderate genetic diversity (uHe: 0.254–0.322), while climate change is expected to have a profound impact on their range size during the coming decades. A combination of in- and ex-situ measures, such as population reinforcement and seed bank conservation, are urgently needed in order to preserve these highly threatened and rare Greek endemics.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Drivers of population divergence and species differentiation in a recent group of indigenous orchids (Vanilla spp.) in Madagascar

With over 25,000 species, orchids are among families with remarkable high rate of diversification. Since Darwin's time, major advances attributed the exceptional diversity of orchids to plant-pollinator interactions. However, unraveling the processes and factors that determine the phenotypic and genotypic variation of natural orchid populations remains a challenge. Here, we assessed genetic population structure and floral differentiation in recently diverged leafless Vanilla species in a world biodiversity hotspot, Madagascar, using seven microsatellite loci and 26 morphometric variables. Additionally, analyses were performed to test for the occurrence of any patterns of isolation by distance, isolation by environment, and isolation by adaptation and to detect possible physical barriers that might have caused genetic discontinuities between populations. Positive inbreeding coefficients detected in 22 populations were probably due to the presence of null alleles, geitonogamy and/or some admixture (sympatric species). In contrast, the only high-altitude population showed an important rate of clonality leading to heterozygote excess. Genetic diversity was maximum in western populations, suggesting a postglacial colonization to the north and south. Clustering analyses identified seven genetic groups characterized by specific floral traits that matched five botanical descriptions in the literature. A contribution of montane refugia and river barriers on population differentiation was detected. We also detected combined effects of IBD/IBE and IBE/IBA on genetic differentiation and suggested this pattern is more likely determined by ecological isolation, although pollinator-mediated divergent selection could not be ruled out for some of the species. Overall, this study provides further insights on speciation in orchids, a group for which Madagascar shows one of the world's highest level of endemism and confirms the importance of the peculiar biogeography of the island in shaping species differentiation.