Comparison between Wild and Hatchery Populations of Korean Pen Shell (Atrina pectinata) Using Microsatellite DNA Markers

Pen shell (Atrina pectinata) is a popular food source with a high commercial value in a number of Asian Pacific areas. The natural A. pectinata population has been declining continuously over the past several decades. Microsatellite DNA markers are a useful DNA-based tool for monitoring the genetic variation of pen shell populations. In this study, 20 polymorphic microsatellite (MS) DNA markers were identified from a partial genomic pen shell DNA library enriched in CA repeats, and used to compare allelic variation between wild and hatchery pen shell populations in Korea. A total of 438 alleles were detected at the 20 MS loci in the two populations. All loci were easily amplified and demonstrated allelic variability, with the number of alleles ranging from 5 to 35 in the wild population and from 5 to 22 in the farmed population. The average observed and expected heterozygosities were 0.69 and 0.82, respectively, in the hatchery samples and 0.69 and 0.83, respectively, in the wild samples. Statistical analysis of fixation index (F(ST)) and analysis of molecular variance (AMOVA) showed minor, but significant, genetic differences between the wild and hatchery populations (F(ST) = 0.0106, CI(95%) = 0.003-0.017). These microsatellite loci may be valuable for future aquaculture and population genetic studies for developing conservation and management plans. Further studies with additional pen shell samples are needed to conclusively determine the genetic diversity between the wild and hatchery populations.     

Isolation and Characterization of Microsatellite Loci in the Chinese Cobra Naja atra (Elapidae)

We characterize thirteen polymorphic microsatellite loci isolated from Naja atra genomic libraries, which were enriched for AC-motif microsatellites. The thirteen loci were screened on a group of 48 individuals from two populations, one in Yong'an and the other in Ganzhou. These markers revealed a relatively high degree of genetic diversity (4-12 alleles per locus) and heterozygosity (Ho ranged from 0.213-0.854 and He ranged from 0.301-0.838). Tests for departure from Hardy-Weinberg equilibrium and for linkage disequilibrium were conducted for each of the two populations separately. After sequential Bonferroni correction, none of the 13 loci showed significant departures from Hardy-Weinberg equilibrium. Hierarchical analysis of molecular variance indicated that a small but significant (P < 0.001) proportion (16.0%) of the total variation in the microsatellite DNA data were attributable to differences among populations, indicating geographical structuring and restricted gene flow. It could be attributable to the Wuyi mountains in the area having a sufficiently isolating effect to significantly reduce gene flow. Our microsatellite data also showed a low N(m) (1.31) value in the two populations from mainland China. Thus, the Yong'an and Ganzhou populations could be treated as distinct evolutionarily significant units (ESUs). The high level of polymorphism revealed by these microsatellite markers will be useful for the study of gene flow, population structure and evolutionary history of N. atra.     

Phylogeography and Genetic Differentiation among Populations of the Moon Turban Snail Lunella granulata Gmelin, 1791 (Gastropoda: Turbinidae)

We examined the genetic variation and phylogeographic relationships among 10 populations of Lunella granulata from mainland China, Penghu Archipelago, Taiwan Island, and Japan using mitochondrial COI and 16S markers. A total of 45 haplotypes were obtained in 112 specimens, and relatively high levels of haplotype diversity (h = 0.903) and low levels of nucleotide diversity (π = 0.0046) were detected. Four major phylogenetic lineage clusters were revealed and were concordant with their geographic distribution, agreeing with the haplotype network. These results suggested that geographic barrier isolating effects were occurring among the populations. This hypothesis was also supported by a significant genetic differentiation index (F_ST = 0.709) and by a spatial analysis of molecular variance (SAMOVA) analysis. A mismatch distribution analysis, neutrality tests and Bayesian skyline plots found a single significant population expansion. This expansion occurred on the coast of mainland China before 20–17 ka. Consequently, although the dispersal ability of the planktonic stage and the circulation of ocean currents generally promote genetic exchanges among populations, L. granulata has tended to maintain distinct genetic groups that reflect the respective geographic origins of the constituent lineages. Although the circulation of ocean currents, in principle, may still play a role in determining the genetic composition of populations, long-distance migration between regions is difficult even at the planktonic stage.     

Wild and hatchery populations of Korean starry flounder (Platichthys stellatus) compared using microsatellite DNA markers

Starry flounder (Platichthys stellatus) is an important sport and food fish found around the margins of the North Pacific. Aquaculture production of this species in Korea has increased because of its commercial value. Microsatellite DNA markers are a useful DNA-based tool for monitoring the genetic variation of starry flounder populations. In this study, 12 polymorphic microsatellite DNA markers were identified from a partial genomic starry flounder DNA library enriched in CA repeats, and used to compare allelic variation between wild and hatchery starry flounder populations in Korea. All loci were readily amplified and demonstrated high allelic diversity, with the number of alleles ranging from 6 to 18 in the wild population and from 2 to 12 in the farmed population. A total of 136 alleles were detected at the 12 microsatellite loci in the two populations. The mean observed and expected heterozygosities were 0.62 and 0.68, respectively, in the hatchery samples and 0.67 and 0.75, respectively, in the wild samples. These results indicate lower genetic variability in the hatchery population as compared to the wild population. Significant shifts in allelic frequencies were detected at eight loci, which resulted in a small but significant genetic differences between the wild and hatchery populations (F(ST) = 0.043, P < 0.05). Further studies with additional starry flounder sample collections are needed for comprehensive determinations of the genetic varieties between the wild and hatchery populations. These microsatellite loci may be valuable for future population genetic studies, monitoring the genetic variation for successful aquaculture management and the preservation of aquatic biodiversity.     

Application of Novel Polymorphic Microsatellite Loci Identified in the Korean Pacific Abalone (Haliotis diversicolor supertexta (Haliotidae)) in the Genetic Characterization of Wild and Released Populations

The small abalone, Haliotis diversicolor supertexta, of the family Haliotidae, is one of the most important species of marine shellfish in eastern Asia. Over the past few decades, this species has drastically declined in Korea. Thus, hatchery-bred seeds have been released into natural coastal areas to compensate for the reduced fishery resources. However, information on the genetic background of the small abalone is scarce. In this study, 20 polymorphic microsatellite DNA markers were identified using next-generation sequencing techniques and used to compare allelic variation between wild and released abalone populations in Korea. Using high-throughput genomic sequencing, a total of 1516 (2.26%; average length of 385 bp) reads containing simple sequence repeats were obtained from 86,011 raw reads. Among the 99 loci screened, 28 amplified successfully, and 20 were polymorphic. When comparing allelic variation between wild and released abalone populations, a total of 243 different alleles were observed, with 18.7 alleles per locus. High genetic diversity (mean heterozygosity = 0.81; mean allelic number = 15.5) was observed in both populations. A statistical analysis of the fixation index (F_ST) and analysis of molecular variance (AMOVA) indicated limited genetic differences between the two populations (F_ST = 0.002, p > 0.05). Although no significant reductions in the genetic diversity were found in the released population compared with the wild population (p > 0.05), the genetic diversity parameters revealed that the seeds released for stock abundance had a different genetic composition. These differences are likely a result of hatchery selection and inbreeding. Additionally, all the primer pair sets were effectively amplified in another congeneric species, H. diversicolor diversicolor, indicating that these primers are useful for both abalone species. These microsatellite loci may be valuable for future aquaculture and population genetic studies aimed at developing conservation and management plans for these two abalone species.     

Integrative Transcriptome,Genome and Quantitative Trait Loci Resources Identify Single Nucleotide Polymorphisms in Candidate Genes for Growth Traits in Turbot

Growth traits represent a main goal in aquaculture breeding programs and may be related to adaptive variation in wild fisheries. Integrating quantitative trait loci (QTL) mapping and next generation sequencing can greatly help to identify variation in candidate genes, which can result in marker-assisted selection and better genetic structure information. Turbot is a commercially important flatfish in Europe and China, with available genomic information on QTLs and genome mapping. Muscle and liver RNA-seq from 18 individuals was carried out to obtain gene sequences and markers functionally related to growth, resulting in a total of 20,447 genes and 85,344 single nucleotide polymorphisms (SNPs). Many growth-related genes and SNPs were identified and placed in the turbot genome and genetic map to explore their co-localization with growth-QTL markers. Forty-five SNPs on growth-related genes were selected based on QTL co-localization and relevant function for growth traits. Forty-three SNPs were technically feasible and validated in a wild Atlantic population, where 91% were polymorphic. The integration of functional and structural genomic resources in turbot provides a practical approach for QTL mining in this species. Validated SNPs represent a useful set of growth-related gene markers for future association, functional and population studies in this flatfish species.     

Tetra-Repeat Microsatellite Markers for the Masu Salmon (Oncorhynchus masou masou) and Its Application in Cross-Subspecies Amplification

We developed tetranucleotide-repeat microsatellite markers for the masu salmon (Oncorhynchus masou) complex. 454 pyrosequencing was used to discover repeat motifs, and seven polymorphic microsatellite-primer sets were identified. The number of alleles detected at each locus ranged from four to 24 and the expected heterozygosity varied from 0.57 to 0.92. Cross-subspecies amplification for O. m. masou, O. m. ishikawae and O. m. subsp. was successful. These microsatellites can be utilized in studies of genetic structure, genetic diversity, and intra- and inter-subspecific hybridization, making a contribution to conservation and management of the Oncorhynchus masou complex.     

Polymorphic Microsatellite Loci Isolated from the Squalidus argentatus Using PCR-Based Isolation of Microsatellite Arrays (PIMA)

Squalidus argentatus (Sauvage and Dabry de Thiersant 1874) is a small-sized freshwater fish which is distributed in Mainland China, Hainan Island and Taiwan. The populations of S. argentatus have dropped sharply probably due to overharvesting and water pollution recently. Eleven polymorphic microsatellite markers were developed for the cyprinid fish S. argentatus. These new markers were tested on 43 individuals collected from Yangtze River and Qiantang River. The number of alleles, observed and expected heterozygosity per locus, in two populations ranged from 3 to 14, from 0.333 to 0.954 and from 0.480 to 0.928, respectively. Only two loci are significantly deviated from Hardy-Weinberg expectations due to the heterozygote deficiency. No significant linkage disequilibrium was detected between the pairwise comparisons of these loci. These polymorphic microsatellite loci will enable us to study the genetic variation, population structure, and conservation genetics of this species in the future.     

Introgression Between Cultivars and Wild Populations of Momordica charantia L. (Cucurbitaceae) in Taiwan

The landrace strains of Momordica charantia are widely cultivated vegetables throughout the tropics and subtropics, but not in Taiwan, a continental island in Southeast Asia, until a few hundred years ago. In contrast, the related wild populations with smaller fruit sizes are native to Taiwan. Because of the introduction of cultivars for agricultural purposes, these two accessions currently exhibit a sympatric or parapatric distribution in Taiwan. In this study, the cultivars and wild samples from Taiwan, India, and Korea were collected for testing of their hybridization and evolutionary patterns. The cpDNA marker showed a clear distinction between accessions of cultivars and wild populations of Taiwan and a long divergence time. In contrast, an analysis of eight selectively neutral nuclear microsatellite loci did not reveal a difference between the genetic structures of these two accessions. A relatively short divergence time and frequent but asymmetric gene flows were estimated based on the isolation-with-migration model. Historical and current introgression from cultivars to wild populations of Taiwan was also inferred using MIGRATE-n and BayesAss analyses. Our results showed that these two accessions shared abundant common ancestral polymorphisms, and the timing of the divergence and colonization of the Taiwanese wild populations is consistent with the geohistory of the Taiwan Strait land bridge of the Last Glacial Maximum (LGM). Long-term and recurrent introgression between accessions indicated the asymmetric capacity to receive foreign genes from other accessions. The modern introduction of cultivars of M. charantia during the colonization of Taiwan by the Han Chinese ethnic group enhanced the rate of gene replacement in the native populations and resulted in the loss of native genes.     

Phylogeographical Analysis on Squalidus argentatus Recapitulates Historical Landscapes and Drainage Evolution on the Island of Taiwan and Mainland China

Phylogeographical analyses on Squalidus argentatus samples from thirteen localities within mainland China and Taiwan were conducted for biogeographic studies, as their dispersal strictly depends on geological evolution of the landmasses. A total of 95 haplotypes were genotyped for mtDNA cyt b gene in 160 specimens from nine river systems. Relatively high levels of haplotype diversity (h = 0.984) and low levels of nucleotide diversity (π = 0.020) were detected in S. argentatus. Two major phylogenetic haplotype groups, A and B, were revealed via phylogenetic analysis. The degree of intergroup divergence (3.96%) indicates that these groups diverged about 4.55 myr (million years) ago. Haplotype network and population analyses indicated significant genetic structure (F(ST) = 0.775), largely concordant with the geographical location of the populations. According to SAMOVA analysis, we divided these populations into four units: Yangtze-Pearl, Qiantang-Minjiang, Jiulong-Beijiang and Taiwan groups. Mismatch distribution analysis, neutrality tests and Bayesian skyline plots indicated a significant population expansion for lineage A and B, approximately dated 0.35 and 0.04 myr ago, respectively. We found strong geographical organization of the haplotype clades across different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow.     

Molecular Characterization of a New Wheat-Thinopyrum intermedium Translocation Line with Resistance to Powdery Mildew and Stripe Rust

A new wheat-Thinopyrum translocation line CH13-21 was selected from the progenies derived from a cross between wheat-Th. intermedium partial amphiploid TAI7047 and wheat line Mianyang11. CH13-21 was characterized by using genomic in situ hybridization (GISH), multicolor-GISH (mc-GISH), multicolor-fluorescence in situ hybridization (mc-FISH) and chromosome-specific molecular markers. When inoculated with stripe rust and powdery mildew isolates, CH13-21 displayed novel resistance to powdery mildew and stripe rust which inherited from its Thinopyrum parent. The chromosomal counting analyses indicated that CH13-21 has 42 chromosomes, with normal bivalent pairing at metaphase I of meiosis. GISH probed by Th. intermedium genomic DNA showed that CH13-21 contained a pair of wheat-Th. intermedium translocated chromosomes. Sequential mc-FISH analyses probed by pSc119.2 and pAs1 clearly revealed that chromosome arm 6BS of CH13-21 was replaced by Thinopyrum chromatin in the translocation chromosome. The molecular markers analysis further confirmed that the introduced Th. intermedium chromatin in CH13-21 belonged to the long arm of homoeologous group 6 chromosome. Therefore, CH13-21 was a new T6BS.6Ai#1L compensating Robertsonian translocation line. It concludes that CH13-21 is a new genetic resource for wheat breeding programs providing novel variation for disease resistances.     

Development of New Microsatellite DNA Markers from Apostichopus japonicus and Their Cross-Species Application in Parastichopus parvimensis and Pathallus mollis

Twenty microsatellite DNA markers were developed for sea cucumber and used to investigate polymorphisms of 60 wild Apostichopus japonicus individuals collected from China. It revealed that all the markers were polymorphic. A total of 164 alleles were detected at 20 loci. The number of alleles per locus varied from 3 to 17 with an average of 8.2, and the expected heterozygosities of each locus ranged from 0.03 to 0.89 with an average of 0.64. Cross-species amplification was also conducted in Parastichopus parvimensis collected from the United States and Pathallus mollis collected from Peru. The result showed that 17 loci amplified Parastichopus parvimensis DNAs while only 4 loci could amplify Pathallus mollis DNAs. All of the polymorphic markers would be useful for future genetic breeding and the assessment of genetic variation within sea cucumbers.     

Homogeneous Nature of Malaysian Marine Fish Epinephelus fuscoguttatus (Perciformes; Serranidae): Evidence Based on Molecular Markers,Morphology and Fourier Transform Infrared Analysis

Taxonomic confusion exists within the genus Epinephelus due to the lack of morphological specializations and the overwhelming number of species reported in several studies. The homogenous nature of the morphology has created confusion in the Malaysian Marine fish species Epinephelus fuscoguttatus and Epinephelus hexagonatus. In this study, the partial DNA sequence of the 16S gene and mitochondrial nucleotide sequences of two gene regions, Cytochrome Oxidase Subunit I and III were used to investigate the phylogenetic relationship between them. In the phylogenetic trees, E. fuscoguttatus was monophyletic with E. hexagonatus species and morphology examination shows that no significant differences were found in the morphometric features between these two taxa. This suggests that E. fuscoguttatus is not distinguishable from E. hexagonatus species, and that E. fuscoguttatus have been identified to be E. hexagonatus species is likely attributed to differences in environment and ability to camouflage themselves under certain conditions. Interestingly, this finding was also supported by Principal Component Analysis on Attenuated Total Reflectance–Fourier-transform Infrared (ATR-FTIR) data analysis. Molecular, morphological and meristic characteristics were combined with ATR-FTIR analysis used in this study offer new perspectives in fish species identification. To our knowledge, this is the first report of an extensive genetic population study of E. fuscoguttatus in Malaysia and this understanding will play an important role in informing genetic stock-specific strategies for the management and conservation of this highly valued fish.     

Construction and Analysis of Siberian Tiger Bacterial Artificial Chromosome Library with Approximately 6.5-Fold Genome Equivalent Coverage

Bacterial artificial chromosome (BAC) libraries are extremely valuable for the genome-wide genetic dissection of complex organisms. The Siberian tiger, one of the most well-known wild primitive carnivores in China, is an endangered animal. In order to promote research on its genome, a high-redundancy BAC library of the Siberian tiger was constructed and characterized. The library is divided into two sub-libraries prepared from blood cells and two sub-libraries prepared from fibroblasts. This BAC library contains 153,600 individually archived clones; for PCR-based screening of the library, BACs were placed into 40 superpools of 10 × 384-deep well microplates. The average insert size of BAC clones was estimated to be 116.5 kb, representing approximately 6.46 genome equivalents of the haploid genome and affording a 98.86% statistical probability of obtaining at least one clone containing a unique DNA sequence. Screening the library with 19 microsatellite markers and a SRY sequence revealed that each of these markers were present in the library; the average number of positive clones per marker was 6.74 (range 2 to 12), consistent with 6.46 coverage of the tiger genome. Additionally, we identified 72 microsatellite markers that could potentially be used as genetic markers. This BAC library will serve as a valuable resource for physical mapping, comparative genomic study and large-scale genome sequencing in the tiger.     

Quantitative Trait Loci Mapping for Bacterial Blight Resistance in Rice Using Bulked Segregant Analysis

Oryza meyeriana is highly resistant to rice bacterial blight (BB) and this resistance trait has been transferred to cultivated rice (O. sativa) using asymmetric somatic hybridization. However, no resistance genes have yet been cloned. In the present study, a progeny of the somatic hybridization with high BB resistance was crossed with a rice cultivar with high BB susceptibility to develop an F₂ population. Using bulked segregant analysis (BSA), 17 polymorphic markers that were linked to rice BB resistance were obtained through scanning a total of 186 simple sequence repeats (SSR) and sequence-tagged site (STS) markers, evenly distributed on 12 chromosomes. A genetic linkage map was then constructed based on the 17 linkage markers and the F₂ segregating population, which was followed by mapping for quantitative trait loci (QTLs) for BB resistance. Three QTLs were identified on chromosomes 1, 3 and 5, respectively, and the alleles of the resistant parent at any of the QTLs increased BB resistance. All of the three QTLs had a strong effect on resistance, explaining about 21.5%, 12.3% and 39.2% of the resistance variance, respectively. These QTLs were different from the loci of the BB resistance genes that have been identified in previous studies. The QTLs mapped in this work will facilitate the isolation of novel BB resistance genes and their utilization in rice resistance breeding.     

Development and characterization of 20 microsatellite mrkers for Chinese back sleeper, Bostrychus sinensis

Twenty microsatellite markers were isolated and characterized from the Chinese black sleeper, Bostrychus sinensis. Loci were screened in 30 individuals from Taiwan. For each locus, the number of alleles varied from 4 to 22 with mean expected and observed heterozygosity of 0.79 and 0.66, respectively. One locus significantly deviated from Hardy-Weinberg equilibrium after Bonferroni correction and no significant linkage disequilibrium was detected. This set of microsatellites will provide a suitable tool for population genetic studies of Chinese black sleeper.     

Development of Simple Sequence Repeats (SSR) Markers in Setaria italica (Poaceae) and Cross-Amplification in Related Species

Foxtail millet is one of the world's oldest cultivated crops. It has been adopted as a model organism for providing a deeper understanding of plant biology. In this study, 45 simple sequence repeats (SSR) markers of Setaria italica were developed. These markers showing polymorphism were screened in 223 samples from 12 foxtail millet populations around Taiwan. The most common dinucleotide and trinucleotide repeat motifs are AC/TG (84.21%) and CAT (46.15%). The average number of alleles (N(a)), the average heterozygosities observed (H(o)) and expected (H(e)) are 3.73, 0.714, 0.587, respectively. In addition, 24 SSR markers had shown transferability to six related Poaceae species. These new markers provide tools for examining genetic relatedness among foxtail millet populations and other related species. It is suitable for germplasm management and protection in Poaceae.     

High Genetic Diversity and Low Differentiation of Michelia coriacea (Magnoliaceae), a Critically Endangered Endemic in Southeast Yunnan, China

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 G(st) = 0.187 for ISSR and Wright's F(st) = 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.     

Genomic Resources of Three Pulsatilla Species Reveal Evolutionary Hotspots,Species-Specific Sites and Variable Plastid Structure in the Family Ranunculaceae

Background: The European continent is presently colonized by nine species of the genus Pulsatilla, five of which are encountered only in mountainous regions of southwest and south-central Europe. The remaining four species inhabit lowlands in the north-central and eastern parts of the continent. Most plants of the genus Pulsatilla are rare and endangered, which is why most research efforts focused on their biology, ecology and hybridization. The objective of this study was to develop genomic resources, including complete plastid genomes and nuclear rRNA clusters, for three sympatric Pulsatilla species that are most commonly found in Central Europe. The results will supply valuable information about genetic variation, which can be used in the process of designing primers for population studies and conservation genetics research. The complete plastid genomes together with the nuclear rRNA cluster can serve as a useful tool in hybridization studies. Methodology/principal findings: Six complete plastid genomes and nuclear rRNA clusters were sequenced from three species of Pulsatilla using the Illumina sequencing technology. Four junctions between single copy regions and inverted repeats and junctions between the identified locally-collinear blocks (LCB) were confirmed by Sanger sequencing. Pulsatilla genomes of 120 unique genes had a total length of approximately 161–162 kb, and 21 were duplicated in the inverted repeats (IR) region. Comparative plastid genomes of newly-sequenced Pulsatilla and the previously-identified plastomes of Aconitum and Ranunculus species belonging to the family Ranunculaceae revealed several variations in the structure of the genome, but the gene content remained constant. The nuclear rRNA cluster (18S-ITS1-5.8S-ITS2-26S) of studied Pulsatilla species is 5795 bp long. Among five analyzed regions of the rRNA cluster, only Internal Transcribed Spacer 2 (ITS2) enabled the molecular delimitation of closely-related Pulsatilla patens and Pulsatilla vernalis. Conclusions/significance: The determination of complete plastid genome and nuclear rRNA cluster sequences in three species of the genus Pulsatilla is an important contribution to our knowledge of the evolution and phylogeography of those endangered taxa. The resulting data can be used to identify regions that are particularly useful for barcoding, phylogenetic and phylogeographic studies. The investigated taxa can be identified at each stage of development based on their species-specific SNPs. The nuclear and plastid genomic resources enable advanced studies on hybridization, including identification of parent species, including their roles in that process. The identified nonsynonymous mutations could play an important role in adaptations to changing environments. The results of the study will also provide valuable information about the evolution of the plastome structure in the family Ranunculaceae.     

Genetic Variation of the Endangered Gentiana lutea L. var. aurantiaca (Gentianaceae) in Populations from the Northwest Iberian Peninsula

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.