Gentiana lutea L. (G. lutea L.) is an endangered plant, patchily distributed along the mountains of Central and Southern Europe. In this study, inter-simple sequence repeat (ISSR) markers were used to investigate the genetic variation in this species within and among populations of G. lutea L. var. aurantiaca of the Cantabrian Mountains (Northwest Iberian Peninsula). Samples of G. lutea L. collected at different locations of the Pyrenees and samples of G. lutea L. subsp. vardjanii of the Dolomites Alps were also analyzed for comparison. Using nine ISSR primers, 106 bands were generated, and 89.6% of those were polymorphic. The populations from the Northwest Iberian Peninsula were clustered in three different groups, with a significant correlation between genetic and geographic distances. Gentiana lutea L. var. aurantiaca showed 19.8% private loci and demonstrated a remarkable level of genetic variation, both among populations and within populations; those populations with the highest level of isolation show the lowest genetic variation within populations. The low number of individuals, as well as the observed genetic structure of the analyzed populations makes it necessary to protect them to ensure their survival before they are too small to persist naturally.
alpine plants; aurantiaca; Cantabrian Mountains; genetic diversity; geographical isolation; Gentiana lutea L.; small population; habitat fragmentation; over-exploitation
The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure.
The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France.
The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.
Cytonuclear incongruence; Gene flow; Phylogeography; Secondary contact; Speciation; Vicariance
Repeated Quaternary glaciations have significantly shaped the present distribution and diversity of several European species in aquatic and terrestrial habitats. To study the phylogeography of freshwater invertebrates, patterns of intraspecific variation have been examined primarily using mitochondrial DNA markers that may yield results unrepresentative of the true species history. Here, population genetic parameters were inferred for a montane aquatic caddisfly, Thremma gallicum, by sequencing a 658-bp fragment of the mitochondrial CO1 gene, and 12,514 nuclear RAD loci. T. gallicum has a highly disjunct distribution in southern and central Europe, with known populations in the Cantabrian Mountains, Pyrenees, Massif Central, and Black Forest. Both datasets represented rangewide sampling of T. gallicum. For the CO1 dataset, this included 352 specimens from 26 populations, and for the RAD dataset, 17 specimens from eight populations. We tested 20 competing phylogeographic scenarios using approximate Bayesian computation (ABC) and estimated genetic diversity patterns. Support for phylogeographic scenarios and diversity estimates differed between datasets with the RAD data favouring a southern origin of extant populations and indicating the Cantabrian Mountains and Massif Central populations to represent highly diverse populations as compared with the Pyrenees and Black Forest populations. The CO1 data supported a vicariance scenario (north–south) and yielded inconsistent diversity estimates. Permutation tests suggest that a few hundred polymorphic RAD SNPs are necessary for reliable parameter estimates. Our results highlight the potential of RAD and ABC-based hypothesis testing to complement phylogeographic studies on non-model species.
Biogeography; freshwater ecology; next-generation sequencing; phylogeography; Pleistocene glaciations
Europe is characterised by several high mountain systems dominating major parts of its area, and these structures have strongly influenced the evolution of taxa. For species now restricted to these high mountain systems, characteristic biogeographical patterns of differentiation exist. (i) Many local endemics are found in most of the European high mountain systems especially in the Alps and the more geographically peripheral regions of Europe. Populations isolated in these peripheral mountain ranges often have strongly differentiated endemic genetic lineages, which survived and evolved in the vicinity of these mountain areas over long time periods. (ii) Populations of taxa with wide distributions in the Alps often have two or more genetic lineages, which in some cases even have the status of cryptic species. In many cases, these lineages are the results of several centres of glacial survival in the perialpine areas. Similar patterns also apply to the other geographically extended European high mountain systems, especially the Pyrenees and Carpathians. (iii) Populations from adjoining high mountain systems often show similar genetic lineages, a phenomenon best explained by postglacial retreat to these mountains from one single differentiation centre between them. (iv) The populations of a number of species show gradients of genetic diversity from a genetically richer East to a poorer West. This might indicate better glacial survival conditions for this biogeographical group of species in the more eastern parts of Europe.
Infectious keratoconjunctivitis (IKC) is a clinical condition affecting eyes of domestic and wild Caprinae worldwide, and Mycoplasma conjunctivae is considered the primary causative agent of IKC in sheep, goats and wild Caprinae. Domestic ruminants from high mountain habitats share grazing areas with wild mountain ungulates, such as chamois (Rupicapra spp.), Alpine ibex (Capra ibex) and European mouflon (Ovis aries musimon), and domestic sheep seem to act as M. conjunctivae reservoir. In this study, the presence of M. conjunctivae in domestic sheep and goats from the two main mountain ranges of Northern Spain, the Pyrenees and the Cantabrian Mountains, has been investigated.
Eye swabs were obtained from 439 domestic small ruminants selected from flocks that seasonally graze in alpine meadows during three consecutive years (2011-2012-2013). Seventy-nine out of the 378 domestic sheep (20.9%) tested positive to a M. conjunctivae specific real time-PCR (rt-PCR) in at least one eye, whereas all the 61 sampled domestic goats were negative. Statistically significant higher prevalence and higher proportion of infected flocks (P < 0.001) was observed in the Pyrenees (25.7%; 12 flocks out of 13), where M. conjunctivae is widespread and probably endemic in domestic sheep, than in the Cantabrian Mountains (7.8%; one flock out of six). Twenty-five sheep (three from the Pyrenees and 22 from the Cantabrian Mountains) which showed clinical signs consistent with infectious keratoconjunctivitis (IKC) were negative by rt-PCR. In contrast, 62 out of the 71 (87.3%) M. conjunctivae-positive sheep from the Pyrenees and the eight positive sheep from the Cantabrian Mountains were asymptomatic.
This study provides rt-PCR-based evidences of M. conjunctivae maintenance in domestic sheep, as well as a relationship between prevalence in domestic sheep and previously reported M. conjunctivae and IKC in wild ruminants. Domestic goats do not seem to play an important role in the epidemiology of M. conjunctivae in alpine habitats from Northern Spain.
Goat; Sheep; Infectious keratoconjunctivitis; Mycoplasma conjunctivae; Spain; Pyrenees; Cantabrian mountains
Genetic variation within plant species is determined by a number of factors such as reproductive mode, breeding system, life history traits and climatic events. In alpine regions, plants experience heterogenic abiotic conditions that influence the population's genetic structure. The aim of this study was to investigate the genetic structure and phylogeographic history of the subalpine shrub Rhododendron ferrugineum across the Pyrenees and the links between the populations in the Pyrenees, the Alps and Jura Mountains. We used 27 microsatellite markers to genotype 645 samples from 29 Pyrenean populations, three from the Alps and one from the Jura Mountains. These data were used to estimate population genetics statistics such as allelic richness, observed heterozygosity, expected heterozygosity, fixation index, inbreeding coefficient and number of migrants. Genetic diversity was found to be higher in the Alps than in the Pyrenees suggesting colonization waves from the Alps to the Pyrenees. Two separate genetic lineages were found in both the Alps and Pyrenees, with a substructure of five genetic clusters in the Pyrenees where a loss of genetic diversity was noted. The strong differentiation among clusters is maintained by low gene flow across populations. Moreover, some populations showed higher genetic diversity than others and presented rare alleles that may indicate the presence of alpine refugia. Two lineages of R. ferrugineum have colonized the Pyrenees from the Alps. Then, during glaciation events R. ferrugineum survived in the Pyrenees in different refugia such as lowland refugia at the eastern part of the chain and nunataks at high elevations leading to a clustered genetic pattern.
This study was focused on conservation genetics of threatened Hippocampus guttulatus on the Atlantic coast of NW Iberian Peninsula. Information about spatial structure and temporal stability of wild populations was obtained based on microsatellite markers, and used for monitoring a captive breeding program firstly initiated in this zone at the facilities of the Institute of Marine Research (Vigo, Spain). No significant major genetic structure was observed regarding the biogeographical barrier of Cape Finisterre. However, two management units under continuous gene flow are proposed based on the allelic differentiation between South-Atlantic and Cantabrian subpopulations, with small to moderate contemporary effective size based on single-sample methods. Temporal stability was observed in South-Atlantic population samples of H. guttulatus for the six-year period studied, suggesting large enough effective population size to buffer the effects of genetic drift within the time frame of three generations. Genetic analysis of wild breeders and offspring in captivity since 2009 allowed us to monitor the breeding program founded in 2006 in NW Spain for this species. Similar genetic diversity in the renewed and founder broodstock, regarding the wild population of origin, supports suitable renewal and rearing processes to maintain genetic variation in captivity. Genetic parentage proved single-brood monogamy in the wild and in captivity, but flexible short- and long-term mating system under captive conditions, from strict monogamy to polygamy within and/or among breeding seasons. Family analysis showed high reproductive success in captivity under genetic management assisted by molecular relatedness estimates to avoid inbreeding. This study provides genetic information about H. guttulatus in the wild and captivity within an uncovered geographical range for this data deficient species, to be taken into account for management and conservation purposes.
Knowledge of the natural genetic variation and structure in a species is important for developing appropriate conservation strategies. We aimed to examine (i) the patterns and levels of morphological and genetic variability within/among populations of Plumbago zeylanica and ascertain whether these variations are dependent on geographical conditions; and (ii) evaluated genetic differentiation and population structure within the species. ISSR and RAPD markers were used to study the genetic structure and variation within and among populations collected from 13 widespread regions in India. High levels of genetic diversity and significantly high genetic differentiation were revealed by both the markers among all studied populations.
Knowledge of the natural genetic variation and structure in a species is important for developing appropriate conservation strategies. As genetic diversity analysis among and within populations of Plumbago
zeylanica remains unknown, we aimed (i) to examine the patterns and levels of morphological and genetic variability within/among populations and ascertain whether these variations are dependent on geographical conditions; and (ii) to evaluate genetic differentiation and population structure within the species. A total of 130 individuals from 13 populations of P. zeylanica were collected, covering the entire distribution area of species across India. The genetic structure and variation within and among populations were evaluated using inter-simple sequence repeat (ISSR) and randomly amplified DNA polymorphism (RAPD) markers. High levels of genetic diversity and significantly high genetic differentiation were revealed by both the markers among all studied populations. High values of among-population genetic diversity were found, which accounted for 60 % of the total genetic variance. The estimators of genetic diversity were higher in northern and eastern populations than in southern and western populations indicating the possible loss of genetic diversity during the spread of this species to Southern India. Bayesian analysis, unweighted pair group method with arithmetic average cluster analysis and principal coordinates analysis all showed similar results. A significant isolation-by-distance pattern was revealed in P. zeylanica by ISSR (r = 0.413, P = 0.05) and RAPD (r = 0.279, P = 0.05) analysis. The results obtained suggest an urgent need for conservation of existing natural populations along with extensive domestication of this species for commercial purpose.
Genetic diversity; molecular markers; Plumbago zeylanica; population structure
Background and Aims
The widely accepted paradigm that the modern genetic structure of plant species in the northern hemisphere has been largely determined by recolonization from refugia after the last glacial maximum fails to explain the presence of cold-tolerant species at intermediate latitudes. Another generally accepted paradigm is that mountain ridges act as important barriers causing genetic isolation of species, but this too has been challenged in recent studies. The aims of the work reported here were to determine the genetic diversity and distribution patterns of extant natural populations of an endangered cool temperate species, Faxinus mandshurica, and to examine whether these two paradigms are appropriate when applied to this species over a wide geographical scale.
1435 adult individuals were sampled from 30 natural populations across the main and central range of the species, covering major mountain ranges across North-east China (NEC). Genetic variation was estimated based on nine polymorphic nuclear microsatellite loci. Phylogeographical analyses were employed using various approaches, including Bayesian clustering, spatial analysis of molecular variance, Monmonier's algorithm, neighbor-joining trees, principal co-ordinate analysis and isolation by distance.
Genetic diversity within populations was relatively high, and no significant recent bottlenecks were detected in any of the populations. A significant negative correlation between intra-population genetic diversity and latitude was identified. In contrast, genetic differentiation among all the populations examined was extremely low and no clear geographic genetic structure was identified, with the exception of one distinct population.
The modern genetic structure in this species can be explained by extensive gene flow, an absence of mountains acting as barriers, and the presence of a wide refuge across NEC rather than multiple small refugia. Intra-population genetic variation along latitudes is probably associated with the systematically northward shifts of forest biomes in eastern China during the mid-Holocene. To determine important genetic patterns and identify resources for conservation, however, it will be necessary to examine differentially inherited genetic markers exposed to selection pressures (e.g. chloroplast DNA) and to investigate different generations.
Fraxinus mandshurica; nuclear microsatellites; latitude variation; historical migration; fossil pollen; spatial genetic structure; genetic barriers
Opisthopappus Shih is endemic to the Taihang Mountains, China. It grows in the crevice of cliffs and is in fragmented distribution. This genus consists of two species, namely, O. taihangensis (Ling) Shih and O. longilobus Shih, which are both endangered plants in China. This study adopted intersimple sequence repeat markers (ISSR) to analyze the genetic diversity and genetic structure from different levels (genus, species, and population) in this genus. A total of 253 loci were obtained from 27 primers, 230 of which were polymorphic loci with a proportion of polymorphic bands (PPB) of up to 90.91% at genus level. At species level, both O. taihangensis (PPB = 90.12%, H = 0.1842, and I = 0.289) and O. longilobus (PPB = 95.21%, H = 0.2226, and I = 0.3542) have high genetic diversity. Their respective genetic variation mostly existed within the population. And genetic variation in O. longilobus (84.95%) was higher than that in O. taihangensis (80.45%). A certain genetic differentiation among populations in O. taihangensis was found (Gst = 0.2740, Φst = 0.196) and genetic differentiation in O. longilobus was very small (Gst = 0.1034, Φst = 0.151). Gene flow in different degrees (Nm = 1.325 and 4.336, resp.) and mating system can form the existing genetic structures of these two species. Furthermore, genetic differentiation coefficient (Gst = 0.0453) between species and the clustering result based on the genetic distance showed that interspecific differentiation between O. taihangensis and O. longilobus was not significant and could occur lately.
The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples.
Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene.
For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula.
In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period.
The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.
A molecular phylogeny of the species from the Trechus brucki clade (previously Trechus uhagoni group)based on fragments of four mitochondrial genes and one nuclear gene is given. We describe Trechus (Trechus) bouilloni
sp. n. from the western pre–Pyrenees: Sierras de Urbasa–Andía, Navarra, Spain. The species was collected in mesovoid shallow substratum (mss), a subterranean environment. Molecular as well as morphological evidences demonstrate that the new species belongs to the Trechus brucki clade. A narrow endemic species of high altitude in western French Pyrenees merged with Trechus brucki Fairmaire, 1862a, Trechus bruckoides
sp. n., is described. A lectotype is designated for Trechus brucki and Trechus planiusculus Fairmaire, 1862b (junior synonym of Trechus brucki). The species group is redefined based on molecular and morphological characters, and renamed as the brucki group, as Trechus brucki was the first described species of the clade. A unique synapomorphy of the male genitalia, a characteristic secondary sclerotization of the sperm duct, which is shared by all the species of the brucki group sensu novo, is described and illustrated. The Trechus brucki group sensu novo is composed of Trechus beusti (Schaufuss, 1863), Trechus bouilloni
sp. n., Trechus brucki, Trechus bruckoides
sp. n., Trechus grenieri Pandellé, 1867, T. uhagoni uhagoni Crotch, 1869, T. uhagoni ruteri Colas, 1935 and Trechus pieltaini Jeannel, 1920. We discuss the taxonomy of the group and provide illustrations of structures showing the differences between the species, along with distribution data and biogeographical comments.
Carabidae; Trechini; Trechus brucki group
; new species; molecular phylogeny; subterranean environment; Pyrenees; France; Spain
Samia cynthia ricini (Lepidoptera:Saturniidae), the Indian eri silkworm, contributes significantly to the production of commercial silk and is widely distributed in the Brahmaputra river valley in North-Eastern India. Due to over exploitation coupled with rapid deforestation, most of the natural populations of S. cynthia ricini are dwindling rapidly and its preservation has become an important goal. Assessment of the genetic structure of each population is a prerequisite for a sustainable conservation program. DNA fingerprinting to detect genetic variation has been used in different insect species not only between populations, but also between individuals within a population. Since, information on the genetic basis of phenotypic variability and genetic diversity within the S. cynthia ricini populations is scanty, inter simple sequence repeat (ISSR) system was used to assess genetic diversity and differentiation among six commercially exploited S. cynthia ricini populations. Twenty ISSR primers produced 87% of inter population variability among the six populations. Genetic distance was lowest between the populations Khanapara (E5) and Mendipathar (E6) (0.0654) and highest between Dhanubhanga (E4) and Titabar (E3) (0.3811). Within population, heterozygosity was higher in Borduar (E2) (0.1093) and lowest in Titabar (E3) (0.0510). Highest gene flow (0.9035) was between E5 and E6 and the lowest (0.2172) was between E3 and E5. Regression analysis showed positive correlation between genetic distance and geographic distance among the populations. The high GST value (0.657) among the populations combined with low gene flow contributes significantly to the genetic differentiation among the S. cynthia ricini populations. Based on genetic diversity, these populations can be considered as different ecotypes and in situ conservation of them is recommended.
Eri phenotype; geographic isolation; gene flow; heterozygosity
Amphibians in general are poor dispersers and highly philopatric, and landscape features often have important impacts on their population genetic structure and dispersal patterns. Numerous studies have suggested that genetic differentiation among amphibian populations are particularly pronounced for populations separated by mountain ridges. The Tsinling Mountain range of northern China is a major mountain chain that forms the boundary between the Oriental and Palearctic zoogeographic realms. We studied the population structure of the Chinese wood frog (Rana chensinensis) to test whether the Tsinling Mountains and the nearby Daba Mountains impose major barriers to gene flow.
Using 13 polymorphic microsatellite DNA loci, 523 individuals from 12 breeding sites with geographical distances ranging from 2.6 to 422.8 kilometers were examined. Substantial genetic diversity was detected at all sites with an average of 25.5 alleles per locus and an expected heterozygosity ranging from 0.504 to 0.855, and two peripheral populations revealed significantly lower genetic diversity than the central populations. In addition, the genetic differentiation among the central populations was statistically significant, with pairwise FST values ranging from 0.0175 to 0.1625 with an average of 0.0878. Furthermore, hierarchical AMOVA analysis attributed most genetic variation to the within-population component, and the between-population variation can largely be explained by isolation-by-distance. None of the putative barriers detected from genetic data coincided with the location of the Tsinling Mountains.
The Tsinling and Daba Mountains revealed no significant impact on the population genetic structure of R. chensinensis. High population connectivity and extensive juvenile dispersal may account for the significant, but moderate differentiation between populations. Chinese wood frogs are able to use streams as breeding sites at high elevations, which may significantly contribute to the diminishing barrier effect of mountain ridges. Additionally, a significant decrease in genetic diversity in the peripheral populations supports Mayr's central-peripheral population hypothesis.
Oncomelania hupensis is the unique intermediate host of Schistosoma japonicum, which causes schistosomiasis endemic in the Far East, and especially in mainland China. O. hupensis largely determines the parasite's geographical range. How O. hupensis's genetic diversity is distributed geographically in mainland China has never been well examined with DNA sequence data.
In this study we investigate the genetic variation among O. hupensis from different geographical origins using the combined complete internal transcribed spacer 1 (ITS1) and ITS2 regions of nuclear ribosomal DNA. 165 O. hupensis isolates were obtained in 29 localities from 7 provinces across mainland China: lake/marshland and hill regions in Anhui, Hubei, Hunan, Jiangxi and Jiangsu provinces, located along the middle and lower reaches of Yangtze River, and mountainous regions in Sichuan and Yunnan provinces. Phylogenetic and haplotype network analyses showed distinct genetic diversity and no shared haplotypes between populations from lake/marshland regions of the middle and lower reaches of the Yangtze River and populations from mountainous regions of Sichuan and Yunnan provinces. The genetic distance between these two groups is up to 0.81 based on Fst, and branch time was estimated as 2–6 Ma. As revealed in the phylogenetic tree, snails from Sichuan and Yunnan provinces were also clustered separately. Geographical separation appears to be an important factor accounting for the diversification of the two groups of O. hupensis in mainland China, and probably for the separate clades between snails from Sichuan and Yunnan provinces. In lake/marshland and hill regions along the middle and lower reaches of the Yangtze River, three clades were identified in the phylogenetic tree, but without any obvious clustering of snails from different provinces.
O. hupensis in mainland China may have considerable genetic diversity, and a more complex population structure than expected. It will be of significant importance to consider the genetic diversity of O. hupensis when assessing co-evolutionary interactions with S. japonicum.
The intermediate host of Schistosoma japonicum in Asia is the snail Oncomelania hupensis, which can be separated phenotypically into ribbed- and smooth-shelled morphotypes. In China, the typical morphotype is ribbed-shelled, with its distribution restricted to mainland China. Smooth-shelled snails with varix are also distributed in China, which are considered to belong to the same subspecies as the ribbed-shelled snails. In this study we investigate the genetic variation among O. hupensis from different geographical origins using combined complete ITS1 and ITS2 regions of nuclear ribosomal DNA. Snails including ribbed-shelled and smooth-shelled (but with varix on the shell) from the lake/marshland region of the middle and lower reaches of the Yangtze River, and smooth-shelled snails from mountainous regions of Sichuan and Yunnan provinces, were genetically distinct with no shared haplotypes detected. Furtheremore, the snails from Sichuan and Yunnan provinces were clustered in separate clades in the phylogenetic tree, and three clades were observed for snails from the middle and lower reaches of the Yangtze River. The population diversity of O. hupensis in China is thus considered large, and evolutionary relationships in the host-parasite system of O. hupensis-S. japonicum may be of interest for further research.
Geographic barriers and Quaternary climate changes are two major forces driving the evolution, speciation, and genetic structuring of extant organisms. In this study, we used Pinus armandii and eleven other Asian white pines (subsection Strobus, subgenus Pinus) to explore the influences of geographic factors and Pleistocene climatic oscillations on species in South China, a region known to be centers of plant endemism and biodiversity hotspots. Range-wide patterns of genetic variation were investigated using chloroplast and mitochondrial DNA markers, with extensive sampling throughout the entire range of P. armandii. Both cpDNA and mtDNA revealed that P. armandii exhibits high levels of genetic diversity and significant population differentiation. Three geographically distinct subdivisions corresponding to the Qinling-Daba Mountains (QDM), Himalaya-Hengduan Mountains (HHM) and Yungui Plateau (YGP) were revealed in mainland China by cpDNA. Their break zone was located in the southeastern margin of the Qinghai-Tibetan Plateau (QTP). A series of massive mountains, induced by the QTP uplift, imposed significant geographic barriers to genetic exchange. The disjunct distribution patterns of ancestral haplotypes suggest that a large continuous population of the white pines may have existed from southwest to subtropical China. Repeated range shifts in response to the Pleistocene glaciations led to the isolation and diversification of the subtropical species. The two Taiwanese white pines share a common ancestor with the species in mainland China and obtain their chloroplasts via long-distance pollen dispersal from North Asian pines. Distinct genetic patterns were detected in populations from the Qinling-Daba Mountains, Yungui Plateau, Himalaya-Hengduan Mountains, and subtropical China, indicating significant contributions of geographic factors to the genetic differentiation in white pines. Our study depicts a clear picture of the evolutionary history of Chinese white pines and highlights the heterogeneous contributions of geography and Pleistocene climatic fluctuations to the extremely high plant species diversity and endemism in South China.
The origin and colonisation history after the Quaternary ice ages remain largely unresolved for many plant lineages, mainly owing to a lack of fine-scale studies. Here, we present a molecular phylogeny and a phylogeographic analysis of Antirrhinum, an important model system in plant biology, in the Pyrenees range. Our goal was to reconstruct the evolutionary and colonisation history of four taxa endemic to this region (A. majus subsp. majus, A. majus. subsp. striatum, A. molle, and A. sempervirens) by using a dense sampling strategy, with a total of 452 individuals from 99 populations whose collective distribution spans nearly the entirety of the Pyrenees and adjacent mountains.
Phylogenetic and phylogeographic analyses of the sequences of two plastid (trnS-trnG and trnK-matK) regions revealed the following: (i) historical relationship between the Pyrenees and Iberia (but not with the Alps); (ii) the long persistence of populations in the Pyrenees, at least since the Late Pleistocene; (iii) three different colonisation histories for populations from the Western, Central, and Eastern Pyrenees; (iv) the deep phylogeographic separation of the eastern and western populations; and (v) the colonisation of southern France from the Eastern Pyrenees.
The present study underlines the enormous influence of the glacial history of the mountain ranges on the current configuration of intra- and inter-specific genetic diversity in Antirrhinum, as well as the importance of periglacial areas for the survival of species during glacial periods of the Quaternary.
Antirrhinum; Phylogeny; Phylogeography; Pyrenees; Quaternary
Both Rheum palmatum and R. tanguticum are important but endangered medicinal plants endemic to China. In this study, we aimed to (i) investigate the level and pattern of genetic variability within/among populations of those species; (ii) evaluate genetic differentiation between both species and its relationships and ascertain whether both species are consistent with their current taxonomical treatment as separate species; and (iii) discuss the implications for the effective conservation of two species.
Total 574 individuals from 30 populations of R. palmatum and R. tanguticum were collected, covering the entire distribution range of two species in China. The genetic variation within and among 30 populations was evaluated using inter-simple sequence repeat (ISSR) markers.
Twelve selected ISSR primers generated a total of 175 fragments, 173 (98.86%) of which were polymorphic. The Nei's gene diversity (H) and Shannon's index (I) of both species were high at species level (H = 0.3107, I = 0.4677 for R. palmatum; H = 0.2848, I = 0.4333 for R. tanguticum). But for both species, the genetic diversity was low at population level, and average within-population diversity of R. palmatum was H = 0.1438, I = 0.2151, and that of R. tanguticum was H = 0.1415, I = 0.2126. The hierarchical AMOVA revealed high levels of among-population genetic differentiation in both species, in line with the gene differentiation coefficient and the limited among-population gene flow (R. palmatum: Φst = 0.592, Gst = 0.537, Nm = 0.432; R. tanguticum: Φst = 0.567, Gst = 0.497, Nm = 0.507). By contrast, only 6.52% of the total genetic variance was partitioned between R. palmatum and R. tanguticum. Bayesian analysis, UPGMA cluster analysis, and PCoA analysis all demonstrated the similar results. A significant isolation-by-distance pattern was revealed in R. palmatum (r = 0.547, P = 0.010), but not in R. tanguticum (r = 0.241, P = 0.100). Based on these results, effective conservation strategies were proposed for these two species. The small molecular variance between R. palmatum and R. tanguticum revealed that they had a common ancestor, and we considered that these two species might not be good species.
Knowledge about the population genetic variation of the endangered orchid, Cypripedium japonicum, is conducive to the development of conservation strategies. Here, we examined the levels and partitioning of inter-simple sequence repeat (ISSR) diversity (109 loci) in five populations of this orchid to gain insight into its genetic variation and population structure in Eastern and Central China. It harbored considerably lower levels of genetic diversity both at the population (percentage of polymorphic loci (PPL) = 11.19%, Nei’s gene diversity (H) = 0.0416 and Shannon’s information index (I) = 0.0613) and species level (PPL = 38.53%, H = 0.1273 and I = 0.1928) and a significantly higher degree of differentiation among populations (the proportion of the total variance among populations (Φpt) = 0.698) than those typical of ISSR-based studies in other orchid species. Furthermore, the Nei’s genetic distances between populations were independent of the corresponding geographical distances. Two main clusters are shown in an arithmetic average (UPGMA) dendrogram, which is in agreement with the results of principal coordinate analysis (PCoA) analysis and the STRUCTURE program. In addition, individuals within a population were more similar to each other than to those in other populations. Based on the genetic data and our field survey, the development of conservation management for this threatened orchid should include habitat protection, artificial gene flow and ex situ measures.
Cypripedium japonicum; genetic diversity; ISSR; Orchidaceae
Background and Aims
High mountain ranges of the Mediterranean Basin harbour a large number of narrowly endemic plants. In this study an investigation is made of the levels and partitioning of genetic diversity in Narcissus longispathus, a narrow endemic of south-eastern Spanish mountains characterized by a naturally fragmented distribution due to extreme specialization on a rare habitat type. By using dense sampling of populations across the species' whole geographical range, genetic structuring at different geographical scales is also examined.
Using horizontal starch-gel electrophoresis, allozyme variability was screened at 19 loci for a total of 858 individuals from 27 populations. The data were analysed by means of standard statistical approaches in order to estimate gene diversity and the genetic structure of the populations.
Narcissus longispathus displayed high levels of genetic diversity and extensive diversification among populations. At the species level, the percentage of polymorphic loci was 68 %, with average values of 2·1, 0·11 and 0·14 for the number of alleles per locus, observed heterozygosity and expected heterozygosity, respectively. Southern and more isolated populations tended to have less genetic variability than northern and less-isolated populations. A strong spatial patterning of genetic diversity was found at the various spatial scales. Gene flow/drift equilibrium occurred over distances <4 km. Beyond that distance divergence was relatively more influenced by drift. The populations studied seem to derive from three panmictic units or ‘gene pools’, with levels of admixture being greatest in the central and south-eastern portions of the species' range.
In addition to documenting a case of high genetic diversity in a narrow endemic plant with naturally fragmented populations, the results emphasize the need for dense population sampling and examination of different geographical scales for understanding population genetic structure in habitat specialists restricted to ecological islands.
Allozymes; genetic diversity; geographical scale; habitat isolation; Narcissus longispathus; Mediterranean endemism; mountain range; natural fragmented distribution
The genetic diversity and population structure of Salvia lachnostachys Benth were assessed. Inter Simple Sequence Repeat (ISSR) molecular markers were used to investigate the restricted distribution of S. lachnostachys in Parana State, Brazil. Leaves of 73 individuals representing three populations were collected. DNA was extracted and submitted to PCR-ISSR amplification with nine tested primers. Genetic diversity parameters were evaluated. Our analysis indicated 95.6% polymorphic loci (stress value 0.02) with a 0.79 average Simpson’s index. The Nei-Li distance dendrogram and principal component analysis largely recovered the geographical origin of each sample. Four major clusters were recognized representing each collected population. Nei’s gene diversity and Shannon’s information index were 0.25 and 0.40 respectively. As is typical for outcrossing herbs, the majority of genetic variation occurred at the population level (81.76%). A high gene flow (Nm = 2.48) was observed with a correspondingly low fixation index. These values were generally similar to previous studies on congeneric species. The results of principal coordinate analysis (PCA) and of arithmetic average (UPGMA) were consistent and all three populations appear distinct as in STRUCTURE analysis. In addition, this analysis indicated a majority intrapopulation genetic variation. Despite the human pressure on natural populations our study found high levels of genetic diversity for S. lachnostachys. This was the first molecular assessment for this endemic species with medicinal proprieties and the results can guide for subsequent bioprospection, breeding programs or conservation actions.
intraspecific diversity; ISSR; Labiatae; Southern Brazil
Tree peonies are great ornamental plants associated with a rich ethnobotanical history in Chinese culture and have recently been used as an evolutionary model. The Qinling Mountains represent a significant geographic barrier in Asia, dividing mainland China into northern (temperate) and southern (semi–tropical) regions; however, their flora has not been well analyzed. In this study, the genetic differentiation and genetic structure of Paeonia rockii and the role of the Qinling Mountains as a barrier that has driven intraspecific fragmentation were evaluated using 14 microsatellite markers.
Twenty wild populations were sampled from the distributional range of P. rockii. Significant population differentiation was suggested (FST value of 0.302). Moderate genetic diversity at the population level (HS of 0.516) and high population diversity at the species level (HT of 0.749) were detected. Significant excess homozygosity (FIS of 0.076) and recent population bottlenecks were detected in three populations. Bayesian clusters, population genetic trees and principal coordinate analysis all classified the P. rockii populations into three genetic groups and one admixed Wenxian population. An isolation-by-distance model for P. rockii was suggested by Mantel tests (r = 0.6074, P<0.001) and supported by AMOVA (P<0.001), revealing a significant molecular variance among the groups (11.32%) and their populations (21.22%). These data support the five geographic boundaries surrounding the Qinling Mountains and adjacent areas that were detected with Monmonier's maximum-difference algorithm.
Our data suggest that the current genetic structure of P. rockii has resulted from the fragmentation of a formerly continuously distributed large population following the restriction of gene flow between populations of this species by the Qinling Mountains. This study provides a fundamental genetic profile for the conservation and responsible exploitation of the extant germplasm of this species and for improving the genetic basis for breeding its cultivars.
Despite its wide distribution across the entire Himalayan range, the current status of Podophyllum hexandrum, a highly important anti-cancerous herb, remains endangered. Genetic diversity characterization of 24 populations comprising of 209 individuals representing the whole of the Indian Himalayas revealed that regardless of geographic location, all of the populations are intermixed and are composed broadly of two types of genetic populations. Our findings also suggested that these populations have evolved well in response to the environment. This study will help in the formulation of conservation programs for P. hexandrum populations in this region.
Intraspecific genetic variation in natural populations governs their potential to overcome challenging ecological and environmental conditions. In addition, knowledge of this variation is critical for the conservation and management of endangered plant taxa. Found in the Himalayas, Podophyllum hexandrum is an endangered high-elevation plant species that has great medicinal importance. Here we report on the genetic diversity analysis of 24 P. hexandrum populations (209 individuals), representing the whole of the Indian Himalayas. In the present study, seven amplified fragment length polymorphism (AFLP) primer pairs generated 1677 fragments, of which 866 were found to be polymorphic. Neighbour joining clustering, principal coordinate analysis and STRUCTURE analysis clustered 209 individuals from 24 populations of the Indian Himalayan mountains into two major groups with a significant amount of gene flow (Nm = 2.13) and moderate genetic differentiation Fst(0.196), G′st(0.20). This suggests that, regardless of geographical location, all of the populations from the Indian Himalayas are intermixed and are composed broadly of two types of genetic populations. High variance partitioned within populations (80 %) suggests that most of the diversity is restricted to the within-population level. These results suggest two possibilities about the ancient population structure of P. hexandrum: either all of the populations in the geographical region of the Indian Himalayas are remnants of a once-widespread ancient population, or they originated from two types of genetic populations, which coexisted a long time ago, but subsequently separated as a result of long-distance dispersal and natural selection. High variance partitioned within the populations indicates that these populations have evolved in response to their respective environments over time, but low levels of heterozygosity suggest the presence of historical population bottlenecks.
AMOVA; amplified fragment length polymorphism (AFLP); Baker's rule; genetic structure; Indian Himalayas; Podophyllum hexandrum; self-pollination.
Michelia coriacea, a critically endangered tree, has a restricted and fragmented distribution in Southeast Yunnan Province, China. The genetic diversity, genetic structure and gene flow in the three extant populations of this species were detected by 10 inter-simple sequence repeat (ISSR) markers and 11 simple sequence repeat (SSR) markers. Examination of genetic diversity revealed that the species maintained a relatively high level of genetic diversity at the species level (percentage of polymorphic bands) PPB = 96.36% from ISSRs; PPL (percentage of polymorphic loci) = 95.56% from SSRs, despite several fragmental populations. Low levels of genetic differentiation among the populations of M. coriacea were detected by Nei’s Gst = 0.187 for ISSR and Wright’s Fst = 0.090 for SSR markers, which is further confirmed by Bayesian model-based STRUCTURE and PCoA analysis that could not reveal a clear separation between populations, although YKP was differentiated to other two populations by ISSR markers. Meanwhile, AMOVA analysis also indicated that 22.84% and 13.90% of genetic variation existed among populations for ISSRs and SSRs, respectively. The high level of genetic diversity, low genetic differentiation, and the population, structure imply that the fragmented habitat and the isolated population of M. coriacea may be due to recent over-exploitation. Conservation and management of M. coriacea should concentrate on maintaining the high level of genetic variability through both in and ex-situ conservation actions.
Michelia coriacea; genetic diversity; critically endangered plant; ISSR markers; SSR markers
Sky islands, formed by the highest reaches of mountain tracts physically isolated from one another, represent one of the biodiversity-rich regions of the world. Comparative studies of geographically isolated populations on such islands can provide valuable insights into the biogeography and evolution of species on these islands. The Western Ghats mountains of southern India form a sky island system, where the relationship between the island structure and the evolution of its species remains virtually unknown despite a few population genetic studies.
Methods and Principal Findings
We investigated how ancient geographic gaps and glacial cycles have partitioned genetic variation in modern populations of a threatened endemic bird, the White-bellied Shortwing Brachypteryx major, across the montane Shola forests on these islands and also inferred its evolutionary history. We used Bayesian and maximum likelihood-based phylogenetic and population-genetic analyses on data from three mitochondrial markers and one nuclear marker (totally 2594 bp) obtained from 33 White-bellied Shortwing individuals across five islands. Genetic differentiation between populations of the species correlated with the locations of deep valleys in the Western Ghats but not with geographical distance between these populations. All populations revealed demographic histories consistent with population founding and expansion during the Last Glacial Maximum. Given the level of genetic differentiation north and south of the Palghat Gap, we suggest that these populations be considered two different taxonomic species.
Conclusions and Significance
Our results show that the physiography and paleo-climate of this region historically resulted in multiple glacial refugia that may have subsequently driven the evolutionary history and current population structure of this bird. The first avian genetic study from this biodiversity hotspot, our results provide insights into processes that may have impacted the speciation and evolution of the endemic fauna of this region.