Transcriptome-wide mining, characterization, and development of microsatellite markers in Lychnis kiusiana (Caryophyllaceae).

BACKGROUND: Lychnis kiusiana Makino is an endangered perennial herb native to wetland areas in Korea and Japan. Despite its conservational and evolutionary significance, population genetic resources are lacking for this species. Next-generation sequencing has been accepted as a rapid and cost-effective solution for the identification of microsatellite markers in nonmodel plants. RESULTS: Using Illumina HiSeq 2000 sequencing technology, we assembled 67,498,600 reads into 91,900 contigs and identified 11,403 microsatellite repeat motifs in 9563 contigs. A total of 4510 microsatellite-containing transcripts had Gene Ontology (GO) annotations, and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified 124 pathways with significant scores. Many microsatellites in the L. kiusiana leaf transcriptome were linked to genes involved in the plant response to light intensity, salt stress, temperature stimulus, and nutrient and water deprivation. A total of 12,486 single-nucleotide polymorphisms (SNPs) were identified on transcripts harboring microsatellites. The analysis of nucleotide substitution rates for 2389 unigenes indicated that 39 genes were under strong positive selection. The primers of 6911 microsatellites were designed, and 40 of 50 selected primer pairs were consistently and successfully amplified from 51 individuals. Twenty-five of these were polymorphic, and the average number of alleles per SSR locus was 6.96, with a range from 2 to 15. The observed and expected heterozygosities ranged from 0.137 to 0.902 and 0.131 to 0.827, respectively, and locus-specific FIS estimates ranged from - 0.116 to 0.290. Eleven of the 25 primer pairs were successfully amplified in three additional species of Lychnis: 56% in L. wilfordii, 64% in L. cognata and 80% in L. fulgens. CONCLUSIONS: The transcriptomic SSR markers of Lychnis kiusiana provide a valuable resource for understanding the population genetics, evolutionary history, and effective conservation management of this species. Furthermore, the identified microsatellite loci linked to the annotated genes should be useful for developing functional markers of L. kiusiana. The developed markers represent a potentially valuable source of transcriptomic SSR markers for population genetic analyses with moderate levels of cross-taxon portability.

Genetic diversity and population structure of Lychnis wilfordii (Caryophyllaceae) with newly developed 17 microsatellite markers.

Lychnis wilfordii (Regel) Maxim. is a perennial plant designated as an endangered species by the Korean government because of rapid reduction in its population size. Thus, a population genetic study of this species is needed to establish the strategy for management and conservation based on scientific evidences. The goals of this study were to develop useful microsatellite markers for L. wilfordii and to understand current genetic status of L. wilfordii in Korean peninsula. Seventeen microsatellite markers were identified using next-generation sequencing and bioinformatic analysis and then analyzed genetic diversity in one hundred forty-five individuals from Korea (KI1, KI2, and KP), China (CX, CF) and Russia (RP). Analysis of molecular variance (AMOVA), principal coordinates analysis (PCoA) and STRUCTURE results consistently showed discontinuity among L. wilfordii populations. AMOVA showed that the percentage of variation among populations was 53%, which was higher than the variation within populations (19%). PCoA showed that the populations were divided into three genetic clusters, (1) Chinese (CX, CF), (2) Russian (RP) populations and Korean populations (KI1, KI2) excluding KP, and (3) the KP population. In particular, KP, the most southern population on the Korean peninsula, showed significantly lower observed and expected heterozygosity, number of effective alleles, and Shannon index (I) than those of KI1 and KI2. L. wilfordii showed high differentiation between populations with low genetic diversity within populations. Among Korean populations, KP is likely to be affected by genetic drift due to small population size, low genetic diversity and limited gene flow.

The complete chloroplast genome sequences of Lychnis wilfordii and Silene capitata and comparative analyses with other Caryophyllaceae genomes.

The complete chloroplast genomes of Lychnis wilfordii and Silene capitata were determined and compared with ten previously reported Caryophyllaceae chloroplast genomes. The chloroplast genome sequences of L. wilfordii and S. capitata contain 152,320 bp and 150,224 bp, respectively. The gene contents and orders among 12 Caryophyllaceae species are consistent, but several microstructural changes have occurred. Expansion of the inverted repeat (IR) regions at the large single copy (LSC)/IRb and small single copy (SSC)/IR boundaries led to partial or entire gene duplications. Additionally, rearrangements of the LSC region were caused by gene inversions and/or transpositions. The 18 kb inversions, which occurred three times in different lineages of tribe Sileneae, were thought to be facilitated by the intermolecular duplicated sequences. Sequence analyses of the L. wilfordii and S. capitata genomes revealed 39 and 43 repeats, respectively, including forward, palindromic, and reverse repeats. In addition, a total of 67 and 56 simple sequence repeats were discovered in the L. wilfordii and S. capitata chloroplast genomes, respectively. Finally, we constructed phylogenetic trees of the 12 Caryophyllaceae species and two Amaranthaceae species based on 73 protein-coding genes using both maximum parsimony and likelihood methods.

Colonization and diversification in the African 'sky islands': insights from fossil-calibrated molecular dating of Lychnis (Caryophyllaceae).

The flora on the isolated high African mountains or 'sky islands' is remarkable for its peculiar adaptations, local endemism and striking biogeographical connections to remote parts of the world. Ages of the plant lineages and the timing of their radiations have frequently been debated but remain contentious as there are few estimates based on explicit models and fossil-calibrated molecular clocks. We used the plastid region maturaseK (matK) and a Caryophylloflora paleogenica fossil to infer the age of the genus Lychnis, and constructed a data set of three plastid (matK; a ribosomal protein S16 (rps16); and an intergenic spacer (psbE-petL)) and two nuclear (internal transcribed spacer (ITS) and a region spanning exon 18-24 in the second largest subunit of RNA polymerase II (RPB2)) loci for joint estimation of the species tree and divergence time of the African representatives. The time of divergence of the African high-altitude Lychnis was placed in the late Miocene to early Pliocene. A single speciation event was inferred in the early Pliocene; subsequent speciation took place sporadically from the late Pliocene to the middle Pleistocene. We provide further support for a Eurasian origin of the African 'sky islands' floral elements, which seem to have been recruited via dispersals at different times: some old, as in Lychnis, and others very recent. We show that dispersal and diversification within Africa play an important role in shaping these isolated plant communities.

Genetic variation and selfing rate in Lychnis flos-cuculi along an industrial pollution gradient.

We studied nine populations of a meadow mixed-mating plant Lychnis flos-cuculi growing in a gradient of copper smelter emissions. We hypothesize that metal tolerant populations in the polluted areas have experienced a loss of genetic variation and are more selfing than the populations from the unpolluted areas. One hundred and thirty-five parental plants and 1059 offspring were genotyped with six microsatellite markers. Selfing rates were assessed manually, with Rmes, Mltr and Colony2. Soil toxicity, population density and pollinators' activity were estimated in the studied areas. Populations from the heavily polluted area have experienced a strong founder effect. However, at present, they are characterized by high densities. A recent genetic explosion was registered for the population from the most polluted site, probably due to forest thinning under pollution effects. Selfing rates estimated with different approaches agreed well only for populations with high genetic variation; they comprised 0-0.23 and were similar between polluted and clean areas. Self-fertilization in L. flos-cuculi hardly represents a mechanism for the fixation of advantageous alleles and a barrier for gene flow from non-tolerant populations. The employment of different methods of selfing rate estimation in populations with low genetic variation appears to be necessary, though not a guarantee of reliable conclusions.

The structural and functional connectivity of the grassland plant Lychnis flos-cuculi.

Understanding the relationship between structural and functional connectivity is essential for successful restoration and conservation management, particularly in intensely managed agricultural landscapes. We evaluated the relationship between structural and functional connectivity of the wetland plant Lychnis flos-cuculi in a fragmented agricultural landscape using landscape genetic and network approaches. First, we studied the effect of structural connectivity, such as geographic distance and various landscape elements (forest, agricultural land, settlements and ditch verges), on gene flow among populations as a measurement of functional connectivity. Second, we examined the effect of structural graph-theoretic connectivity measures on gene flow among populations and on genetic diversity within populations of L. flos-cuculi. Among landscape elements, forests hindered gene flow in L. flos-cuculi, whereas gene flow was independent of geographic distance. Among the structural graph-theoretic connectivity variables, only intrapopulation connectivity, which was based on population size, had a significant positive effect on gene flow, that is, more gene flow took place among larger populations. Unexpectedly, interpopulation connectivity of populations, which takes into account the spatial location and distance among populations, did not influence gene flow in L. flos-cuculi. However, higher observed heterozygosity and lower inbreeding was observed in populations characterised by higher structural interpopulation connectivity. This finding shows that a spatially coherent network of populations is significant for maintaining the genetic diversity of populations. Nevertheless, lack of significant relationships between gene flow and most of the structural connectivity measures suggests that structural connectivity does not necessarily correspond to functional connectivity.

Between-population outbreeding affects plant defence.

Between-population crosses may replenish genetic variation of populations, but may also result in outbreeding depression. Apart from direct effects on plant fitness, these outbreeding effects can also alter plant-herbivore interactions by influencing plant tolerance and resistance to herbivory. We investigated effects of experimental within- and between-population outbreeding on herbivore resistance, tolerance and plant fitness using plants from 13 to 19 Lychnis flos-cuculi populations. We found no evidence for outbreeding depression in resistance reflected by the amount of leaf area consumed. However, herbivore performance was greater when fed on plants from between-population compared to within-population crosses. This can reflect outbreeding depression in resistance and/or outbreeding effects on plant quality for the herbivores. The effects of type of cross on the relationship between herbivore damage and plant fitness varied among populations. This demonstrates how between-population outbreeding effects on tolerance range from outbreeding depression to outbreeding benefits among plant populations. Finally, herbivore damage strengthened the observed outbreeding effects on plant fitness in several populations. These results raise novel considerations on the impact of outbreeding on the joint evolution of resistance and tolerance, and on the evolution of multiple defence strategies.

Effects of experimental inbreeding on herbivore resistance and plant fitness: the role of history of inbreeding, herbivory and abiotic factors.

Inbreeding is common in plant populations and can affect plant fitness and resistance against herbivores. These effects are likely to depend on population history. In a greenhouse experiment with plants from 17 populations of Lychnis flos-cuculi, we studied the effects of experimental inbreeding on resistance and plant fitness. Depending on the levels of past herbivory and abiotic factors at the site of plant origin, we found either inbreeding or outbreeding depression in herbivore resistance. Furthermore, when not damaged experimentally by snail herbivores, plants from populations with higher heterozygosity suffered from inbreeding depression and those from populations with lower heterozygosity suffered from outbreeding depression. These effects of inbreeding and outbreeding were not apparent under experimental snail herbivory. We conclude that inbreeding effects on resistance and plant fitness depend on population history. Moreover, herbivory can mask inbreeding effects on plant fitness. Thus, understanding inbreeding effects on plant fitness requires studying multiple populations and considering population history and biotic interactions.

Microsatellite variation and structure of 28 populations of the common wetland plant, Lychnis flos-cuculi L., in a fragmented landscape.

Habitat fragmentation is known to cause genetic differentiation between small populations of rare species and decrease genetic variation within such populations. However, common species with recently fragmented populations have rarely been studied in this context. We investigated genetic variation and its relationship to population size and geographical isolation of populations of the common plant species, Lychnis flos-cuculi L., in fragmented fen grasslands. We analysed 467 plants from 28 L. flos-cuculi populations of different sizes (60 000-54 000 flowering individuals) in northeastern Switzerland using seven polymorphic microsatellite loci. Genetic differentiation between populations is small (F(ST) = 0.022; amova; P < 0.001), suggesting that gene flow among populations is still high or that habitat fragmentation is too recent to result in pronounced differentiation. Observed heterozygosity (H(O) = 0.44) significantly deviates from Hardy-Weinberg equilibrium, and within-population inbreeding coefficient F(IS) is high (0.30-0.59), indicating a mixed mating breeding system with substantial inbreeding in L. flos-cuculi. Gene diversity is the only measure of genetic variation which decreased with decreasing population size (R = 0.42; P < 0.05). While our results do not indicate pronounced effects of habitat fragmentation on genetic variation in the still common L. flos-cuculi, the lower gene diversity of smaller populations suggests that the species is not entirely unaffected.