Genome sequences and population genomics reveal climatic adaptation and genomic divergence between two closely related sweetgum species.

Understanding the genetic basis of population divergence and adaptation is an important goal in population genetics and evolutionary biology. However, the relative roles of demographic history, gene flow, and/or selective regime in driving genomic divergence, climatic adaptation, and speciation in non-model tree species are not yet fully understood. To address this issue, we generated whole-genome resequencing data of Liquidambar formosana and L. acalycina, which are broadly sympatric but altitudinally segregated in the Tertiary relict forests of subtropical China. We integrated genomic and environmental data to investigate the demographic history, genomic divergence, and climatic adaptation of these two sister species. We inferred a scenario of allopatric species divergence during the late Miocene, followed by secondary contact during the Holocene. We identified multiple genomic islands of elevated divergence that mainly evolved through divergence hitchhiking and recombination rate variation, likely fostered by long-term refugial isolation and recent differential introgression in low-recombination genomic regions. We also found some candidate genes with divergent selection signatures potentially involved in climatic adaptation and reproductive isolation. Our results contribute to a better understanding of how late Tertiary/Quaternary climatic change influenced speciation, genomic divergence, climatic adaptation, and introgressive hybridization in East Asia's Tertiary relict flora. In addition, they should facilitate future evolutionary, conservation genomics, and molecular breeding studies in Liquidambar, a genus of important medicinal and ornamental values.

Phylogenomics and diversification drivers of the Eastern Asian - Eastern North American disjunct Podophylloideae.

Evolutionary and biogeographic processes determine species richness patterns of vascular plants between Eastern Asia (EA) and Eastern North America (ENA). However, the strikingly higher species richness of EA relative to ENA remains poorly understood from this perspective. Here, we studied the relative importance of biogeographical, evolutionary and ecological factors underlying differences in species richness between EA and ENA in Podophylloideae (Berberidaceae, Ranunculales; in total 10 spp. in EA vs. 2 spp. in ENA). Based on large-scale transcriptome data, our phylogenomic analyses strongly supported Podophylloideae and its two multi-species genera, i.e. Dysosma (EA) and Diphylleia (EA/ENA), as monophyletic groups. Sinopodophyllum hexandrum (EA) was identified as sister to the remainder of Podophylloideae. Dysosma (7 spp.) was recovered as sister to Podophyllum peltatum (ENA), forming an EA-ENA disjunct pair with a strong bias of species diversity in the EA counterpart. Our biogeographic analyses support the 'out-of-Tibet' hypothesis, suggesting that Podophylloideae started to diversify in the Himalaya-Hengduan Mountains (Mid-Miocene) and migrated eastward (since the Late Miocene) into Central-eastern China, Japan, and ENA (only P. peltatum and Diphylleia cymosa). Overall, we conclude that the striking species diversity anomaly between EA and ENA in Podophylloideae may be explained by a combination of (1) a longer period of time available to accumulate species in EA; and (2) a greater diversification rate in EA, which might have been promoted by greater physiographic and environmental heterogeneity in this region.

Genome sequencing and transcriptome analyses provide insights into the origin and domestication of water caltrop (Trapa spp., Lythraceae).

Humans have domesticated diverse species from across the plant kingdom; however, our current understanding of plant domestication is largely founded on major cereal crops. Here, we examine the evolutionary processes and genetic basis underlying the domestication of water caltrop (Trapa spp., Lythraceae), a traditional, yet presently underutilized non-cereal crop that sustained early Chinese agriculturalists. We generated a chromosome-level genome assembly of tetraploid T. natans, and then divided the allotetraploid genome into two subgenomes. Based on resequencing data from 57 accessions, representing cultivated diploid T. natans, wild T. natans (2x and 4x) and diploid T. incisa, we showed that water caltrop was likely first domesticated in the Yangtze River Valley as early as 6300 yr BP, and experienced a second improvement c. 800 years ago. We also provided strong support for an allotetraploid origin of T. natans within the past 230 000-310 000 years. By integrating selective sweep and transcriptome profiling analyses, we identified a number of genes potentially selected and/or differentially expressed during domestication, some of which likely contributed not only to larger fruit sizes but also to a more vigorous root system, facilitating nutrient uptake, environmental stress response and underwater photosynthesis. Our results shed light on the evolutionary and domestication history of water caltrop, one of the earliest domesticated crops in China. This study has implications for genomic-assisted breeding of this presently underutilized aquatic plant, and improves our general understanding of plant domestication.

[Investigation on reproductive characteristics of Polygonatum cyrtonema].

The study is aimed to investigate the reproductive biology characteristics of Polygonatum cyrtonema, especially including phenology, flower bud differentiation, flowering timing, floral traits, pollen vigor and stigma receptivity. The results showed that P. cyrtonema forms inflorescence before the leaves spread. In the wild, P. cyrtonema is mainly pollinated by insects such as bumblebees, with a seed setting rate of 65.12%. The seed setting rate of indoor single plant isolation or self-pollination enclosed by parchment paper bag is 0, indicating that it is self-incompatible. In Lin'an city, seedlings begin to emerge from mid-March to early April(the temperature is higher than 7.5 ℃), buds begin to emerge from the end of March to mid-April, and then undergo the full bloom stage from mid-to-late April, and the final flowering stage from the end of April to mid-May. The whole flowering period lasts 36 to 45 days. There are obvious differences in the phenology of different provenances. The flowers come into bloom from the base to the top along the aboveground main axis, which usually contain 4-22 inflorescences with(2-) 4-10(-21) flowers per inflorescence. The flowering pe-riod for a single plant is 26-38 days. The single flower lasts about 20-25 days from budding to opening and withers 2 days after pollination, and then the ovary will gradually expand. If unpollinated, it will continue to bloom for 3-5 days and then wither. Flower development period is significantly related to pollen vigor and stigma remittance. The pollen viability is the highest when the flower is fully opened with anthers gathering on the stigma, and the receptivity is the strongest when the stigma protrudes out of the perianth and secretes mucus. The fruits and seeds ripen in October, and proper shading can ensure the smooth development and maturity of the seeds. This study provides a basis for the hybrid breeding and seed production of P. cyrtonema.

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas.

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

Genomic insights into historical population dynamics, local adaptation, and climate change vulnerability of the East Asian Tertiary relict Euptelea (Eupteleaceae).

The warm-temperate and subtropical climate zones of East Asia are a hotspot of plant species richness and endemism, including a noticeable number of species-poor Tertiary relict tree genera. However, little is understood about when East Asian Tertiary relict plants diversified, how they responded demographically to past environmental change, and to what extent their current genomic composition (and adaptive capacity) might mitigate the effects of global warming. Here, we obtained genomic (RAD-SNP) data for 171 samples from two extant species of Euptelea in China (24 E. pleiosperma populations) and Japan (11 E. polyandra populations) to elucidate their divergence and demographic histories, genome-wide associations with current environmental variables, and genomic vulnerability to future climate change. Our results indicate that Late Miocene changes in climate and/or sea level promoted species divergence, whereas Late Pliocene uplifting in southwest China likely fostered lineage divergence within E. pleiosperma. Its subsequent range expansion into central/east (CE) China bears genomic signatures of climate-driven selection, yet extant CE populations are predicted to be most vulnerable to future climate change. For E. polyandra, geography was the only significant predictor of genomic variation. Our findings indicate a profound impact of Late Neogene geological and climate change on the evolutionary history of Euptelea, with much stronger signals of local adaptation left in China than in Japan. This study deepens our understanding of the complex evolutionary forces that influence the distribution of genetic variation of Tertiary relict trees, and provides insights into their susceptibility to global change and potential for adaptive responses. Our results lay the groundwork for future conservation and restoration programs for Euptelea.

Phylogenomic insights into the temporal-spatial divergence history, evolution of leaf habit and hybridization in Stachyurus (Stachyuraceae).

Stachyuraceae, an East Asian endemic family of shrubs or small trees, comprises a single genus, Stachyurus (c. 11 spp.). Despite previous extensive studies based on both morphology and molecular data, species relationships within Stachyurus are still unresolved. Here, we employed transcriptome data aiming to elucidate the evolutionary history of Stachyurus and investigate possible causes of phylogenetic incongruence among individual gene trees in this genus. Our transcriptome phylogeny strongly supports four major clades of Stachyurus, with S. praecox from Japan being resolved as sister to the remainder of the genus on the Asian mainland. The deciduous S. praecox in Japan appears to have originated in the late Miocene, while the remainder diversified and expanded on the mainland over late Miocene to Pliocene/early Pleistocene times. These latter episodes of diversification and expansion were likely promoted by changes in paleoclimate and orogeny (e.g., late Miocene uplift of the Hengduan Mts. and/or enforcement of the East Asian summer monsoon). Species of this genus evolved from a deciduous ancestor, followed by multiple and independent transitions in leaf habit, possibly reflecting climate-related adaptations. Phylogenetic incongruence observed among individual gene trees may be attributable to incomplete lineage sorting following ancient rapid diversification and frequent interspecific gene flow caused by hybridization events. In sum, this study demonstrates the potential usefulness of genome-wide phylogenetic incongruence and network analyses for reconstructing complex evolutionary histories in rapidly diversifying and naturally hybridizing species groups.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region.

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Contrasting patterns of genetic structure and phylogeography in the marine agarophytes Gelidiophycus divaricatus and G. freshwateri (Gelidiales, Rhodophyta) from East Asia.

The evolutionary and population demographic history of marine red algae in East Asia is poorly understood. Here, we reconstructed the phylogeographies of two upper intertidal species endemic to East Asia, Gelidiophycus divaricatus and G. freshwateri. Phylogenetic and phylogeographic inferences of 393 mitochondrial cox1, 128 plastid rbcL, and 342 nuclear ITS2 sequences were complemented with ecological niche models. Gelidiophycus divaricatus, a southern species adapted to warm water, is characterized by a high genetic diversity and a strong geographical population structure, characteristic of stable population sizes and sudden reduction to recent expansion. In contrast, G. freshwateri, a northern species adapted to cold temperate conditions, is genetically relatively homogeneous with a shallow population structure resulting from steady population growth and recent equilibrium. The overlap zone of the two species roughly matches summer and winter isotherms, indicating that surface seawater temperature is a key feature influencing species range. Unidirectional genetic introgression was detected at two sites on Jeju Island where G. divaricatus was rare while G. freshwateri was common, suggesting the occurrence of asymmetric natural hybrids, a rarely reported event for rhodophytes. Our results illustrate that Quaternary climate oscillations have left strong imprints on the current day genetic structure and highlight the importance of seawater temperature and sea level change in driving speciation in upper intertidal seaweed species.

Phylogenomics recovers monophyly and early Tertiary diversification of Dipteronia (Sapindaceae).

Dipteronia (Sapindaceae) is an ancient relict woody genus, and contains just two extant species endemic to Southwestern and Central China. As sharing numerous morphological characters, Dipteronia and Acer have long been considered as sister groups forming the traditional family Aceraceae. However, molecular phylogenetics has generally not resolved the phylogenetic placement of Dipteronia, especially not in its expected position as sister to Acer. In this study, we present large-scale, phylogenomic data sets, incorporating complete chloroplast (cp) genome sequences and transcriptome data from 13 Sapindaceae species, to resolve the phylogenetic relationship between Dipteronia and Acer. Moreover, the impact of long-branch attraction (LBA) artefacts and robustness of inferred topologies are assessed by long-branch excluding and coalescent-based methods. Corroborating classical morphology-based classifications, both cp genome and nuclear datasets (2466 co-orthologous genes and 273 co-SCNGs) recover Dipteronia and Acer as mutually monophyletic groups. In addition, our fossil-calibrated molecular phylogenies date the crown of the two extant Dipteronia species to the Paleocene/Eocene boundary, implying that these morphologically highly similar taxa are amongst the oldest 'living fossils' of the East Asian Flora.

Plastome organization, genome-based phylogeny and evolution of plastid genes in Podophylloideae (Berberidaceae).

Species of Podophylloideae (Berberidaceae, Ranunculales) are of great pharmacogenetic importance and represent the classic biogeographic disjunction between eastern Asia (EA; 10 ssp.) and eastern North America (ENA; 2 ssp.). However, previous molecular studies of this group suffered from low phylogenetic resolution and/or insufficient marker variability. This study is the first to report whole-plastome sequence data for all 12 species of Podophylloideae (14 individuals) and a close relative, Achlys triphylla. These 15 plastomes proved highly similar in overall size (156,240-157,370 bp), structure, gene order and content, also when compared to other Ranunculales, but also revealed some structural variations caused by the expansion or contraction of the inverted repeats (IRs) into or out of adjacent single-copy regions. Our phylogenomic analysis, based on 63 plastome-derived protein-coding genes (CDS), supported the monophyly of Podophylloideae and its two major genera (EA: Dysosma, EA/ENA: Diphylleia), with Podophyllum peltatum L. (ENA) being more closely related to Diphylleia than to the group's earliest diverging species, Sinopodophyllum hexandrum (EA). Furthermore, within this subfamily/dataset, matK was identified as the fastest evolving gene, which proved to be under positive selection especially in more recently derived, lower-elevation lineages of Dysosma, possibly reflecting an adaptive response to novel environmental (i.e. subtropical compared to higher-elevation/alpine) conditions. Finally, several highly variable noncoding regions were identified in the plastomes of Podophylloideae and Ranunculales. These highly variable loci should be the best choices for future phylogenetic, phylogeographic, and population-level genetic studies. Overall, our results demonstrate the power of plastid phylogenomics to improve phylogenetic resolution, and contribute to a better understanding of plastid gene evolution in Podophylloideae.

Inferring spatial patterns and drivers of population divergence of Neolitsea sericea (Lauraceae), based on molecular phylogeography and landscape genomics.

The relative roles of geography, climate and ecology in driving population divergence and (incipient) speciation has so far been largely neglected in studies addressing the evolution of East Asia's island flora. Here, we employed chloroplast and ribosomal DNA sequences and restriction site-associated DNA sequencing (RADseq) loci to investigate the phylogeography and drivers of population divergence of Neolitsea sericea. These data sets support the subdivision of N. sericea populations into the Southern and Northern lineages across the 'Tokara gap'. Two distinct sublineages were further identified for the Northern lineage of N. sericea from the RADseq data. RADseq was also used along with approximate Bayesian computation to show that the current distribution and differentiation of N. sericea populations resulted from a combination of relatively ancient migration and successive vicariant events that likely occurred during the mid to late Pleistocene. Landscape genomic analyses showed that, apart from geographic barriers, barrier, potentially local adaptation to different climatic conditions appears to be one of the major drivers for lineage diversification of N. sericea.

Development and Application of Genomic Resources in an Endangered Palaeoendemic Tree, Parrotia subaequalis (Hamamelidaceae) From Eastern China.

Parrotia subaequalis is an endangered palaeoendemic tree from disjunct montane sites in eastern China. Due to the lack of effective genomic resources, the genetic diversity and population structure of this endangered species are not clearly understood. In this study, we conducted paired-end shotgun sequencing (2 × 125 bp) of genomic DNA for two individuals of P. subaequalis on the Illumina HiSeq platform. Based on the resulting sequences, we have successfully assembled the complete chloroplast genome of P. subaequalis, as well as identified the polymorphic chloroplast microsatellites (cpSSRs), nuclear microsatellites (nSSRs) and mutational hotspots of chloroplast. Ten polymorphic cpSSR loci and 12 polymorphic nSSR loci were used to genotype 96 individuals of P. subaequalis from six populations to estimate genetic diversity and population structure. Our results revealed that P. subaequalis exhibited abundant genetic diversity (e.g., cpSSRs: Hcp = 0.862; nSSRs: HT = 0.559) and high genetic differentiation (e.g., cpSSRs: RST = 0.652; nSSRs: RST = 0.331), and characterized by a low pollen-to-seed migration ratio (r ≈ 1.78). These genetic patterns are attributable to its long evolutionary histories and low levels of contemporary inter-population gene flow by pollen and seed. In addition, lack of isolation-by-distance pattern and strong population genetic structuring in both marker systems, suggests that long-term isolation and/or habitat fragmentation as well as genetic drift may have also contributed to the geographic differentiation of P. subaequalis. Therefore, long-term habitat protection is the most important methods to prevent further loss of genetic variation and a decrease in effective population size. Furthermore, both cpSSRs and nSSRs revealed that P. subaequalis populations consisted of three genetic clusters, which should be considered as separated conservation units.

Molecular phylogenetics and biogeography of the mint tribe Elsholtzieae (Nepetoideae, Lamiaceae), with an emphasis on its diversification in East Asia.

Elsholtzia and its allied genera such as Collinsonia and Perilla (tribe Elsholtzieae, Lamiaceae) are an ecologically and economically important plant group consisting of ~71 species, with most species distributed in East and Southeast Asia, and several species in North America. Their phylogeny and historical biogeography resulting in a distant intercontinental disjunction are poorly understood. Here we use two nuclear (ETS, ITS) and five chloroplast (rbcL, matK, trnL-F, ycf1, ycf1-rps15) fragments to reconstruct the phylogeny, biogeographic history, and patterns of diversification of Elsholtzieae. The tribe Elsholtzieae is monophyletic and divided into five clades. The woody Elsholtzia species are nested within herbaceous ones and are inferred to have evolved from herbaceous ancestors. Molecular dating shows that the five major clades were established during the Eocene period, but most of the modern diversity did not originate until the Miocene. The divergence between the New World Collinsonia and the Old World Mosla-Keiskea-Perilla clade was dated to the mid-Miocene. Ancestral area reconstructions suggest that the tribe originated in East Asia, and then dispersed to Southeast Asia and North America. Overall, our findings highlight the important roles of the uplifts of the Qinghai-Tibetan Plateau (QTP) and climate changes from Middle Miocene onwards in promoting species diversification of Elsholtzieae.

Comparative Genomics and Phylogenomics of East Asian Tulips (Amana, Liliaceae).

The genus Amana Honda (Liliaceae), when it is treated as separate from Tulipa, comprises six perennial herbaceous species that are restricted to China, Japan and the Korean Peninsula. Although all six Amana species have important medicinal and horticultural uses, studies focused on species identification and molecular phylogenetics are few. Here we report the nucleotide sequences of six complete Amana chloroplast (cp) genomes. The cp genomes of Amana range from 150,613 bp to 151,136 bp in length, all including a pair of inverted repeats (25,629-25,859 bp) separated by the large single-copy (81,482-82,218 bp) and small single-copy (17,366-17,465 bp) regions. Each cp genome equivalently contains 112 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 78 protein coding genes. Gene content, gene order, AT content, and IR/SC boundary structure are nearly identical among all Amana cp genomes. However, the relative contraction and expansion of the IR/SC borders among the six Amana cp genomes results in length variation among them. Simple sequence repeat (SSR) analyses of these Amana cp genomes indicate that the richest SSRs are A/T mononucleotides. The number of repeats among the six Amana species varies from 54 (A. anhuiensis) to 69 (Amana kuocangshanica) with palindromic (28-35) and forward repeats (23-30) as the most common types. Phylogenomic analyses based on these complete cp genomes and 74 common protein-coding genes strongly support the monophyly of the genus, and a sister relationship between Amana and Erythronium, rather than a shared common ancestor with Tulipa. Nine DNA markers (rps15-ycf1, accD-psaI, petA-psbJ, rpl32-trnL, atpH-atpI, petD-rpoA, trnS-trnG, psbM-trnD, and ycf4-cemA) with number of variable sites greater than 0.9% were identified, and these may be useful for future population genetic and phylogeographic studies of Amana species.

Quaternary climate change drives allo-peripatric speciation and refugial divergence in the Dysosma versipellis-pleiantha complex from different forest types in China.

Subtropical China harbours the world's most diverse temperate flora, but little is known about the roles of geographical and eco-climatic factors underlying the region's exceptionally high levels of species diversity and endemism. Here we address this key question by investigating the spatio-temporal and ecological processes of divergence within the Dysosma versipellis-pleiantha species complex, endemic to subtropical China. Our cpDNA phylogeny showed that this monophyletic group of understory herbs is derived from a Late Pliocene ancestor of the Qinghai-Tibetan Plateau (QTP)/Southwest China. Genetic and ENM data in conjunction with niche differentiation analyses support that the early divergence of D. versipellis and D. pleiantha proceeded through allo-peripatric speciation, possibly triggered by Early Pleistocene climate change, while subsequent climate-induced cycles of range contractions/expansions enhanced the eco-geographical isolation of both taxa. Furthermore, modelling of population-genetic data indicated that major lineage divergences within D. versipellis likely resulted from long-term allopatric population isolation in multiple localized refugia over the last glacial/interglacial periods, and which in turn fostered endemic species formation (D. difformis, D. majoensis) from within D. versipellis in Southwest China. These findings point to an overriding role of Quaternary climate change in triggering essentially allopatric (incipient) speciation in this group of forest-restricted plant species in subtropical China.

Ten Years of Landscape Genomics: Challenges and Opportunities.

Landscape genomics is a relatively new discipline that aims to reveal the relationship between adaptive genetic imprints in genomes and environmental heterogeneity among natural populations. Although the interest in landscape genomics has increased since this term was coined, studies on this topic remain scarce. Landscape genomics has become a powerful method to scan and determine the genes responsible for the complex adaptive evolution of species at population (mostly) and individual (more rarely) level. This review outlines the sampling strategies, molecular marker types and research categories in 37 articles published during the first 10 years of this field (i.e., 2007-2016). We also address major challenges and future directions for landscape genomics. This review aims to promote interest in conducting additional studies in landscape genomics.

Chloroplast DNA sequencing and detailed microsatellite genotyping of all remnant populations suggests that only a single genet survives of the critically endangered plant Rehmannia japonica.

Rehmannia japonica (Thunb.) Makino ex T. Yamaz. is an endangered perennial herb species in Japan. Although earlier the Japanese considered it a variety of R. glutinosa, recent Japanese taxonomists have consistently regarded it as an independent species. According to the historical literature, Rehmannia japonica seems to have been known in China and Japan in the past. However, Chinese taxonomists do not recognize R. japonica at present. In Japan, only two populations are known, and although these populations flower every year, seed reproduction has not been observed. In this study, we aimed to reveal the phylogenetic relationships and levels of genetic diversity of R. japonica. A haplotype network based on two chloroplast DNA regions (trnL-trnF and rps16) showed that the sequences of R. japonica were distinguishable by three or four sites of indels from the most closely related species, R. chingii, consistent with the separate species status of R. japonica. An analysis of genetic diversity using twelve microsatellite loci showed that all of the ramets of R. japonica collected from two geographically isolated populations had an identical multilocus genotype, including identical heterozygous genotypes at six loci. This result indicated asexual origin of all sampled ramets. This study also suggests that the absence of sexual reproduction of R. japonica is explained by self-incompatibility combined with only a single genet remaining in the R. japonica populations.

Contributions of historical and contemporary geographic and environmental factors to phylogeographic structure in a Tertiary relict species, Emmenopterys henryi (Rubiaceae).

Examining how historical and contemporary geographic and environmental factors contribute to genetic divergence at different evolutionary scales is a central yet largely unexplored question in ecology and evolution. Here, we examine this key question by investigating how environmental and geographic factors across different epochs have driven genetic divergence at deeper (phylogeographic) and shallower (landscape genetic) evolutionary scales in the Chinese Tertiary relict tree Emmenopterys henryi. We found that geography played a predominant role at all levels - phylogeographic clades are broadly geographically structured, the deepest levels of divergence are associated with major geological or pre-Quaternary climatic events, and isolation by distance (IBD) primarily explained population genetic structure. However, environmental factors are clearly also important - climatic fluctuations since the Last Interglacial (LIG) have likely contributed to phylogeographic structure, and the population genetic structure (in our AFLP dataset) was partly explained by isolation by environment (IBE), which may have resulted from natural selection in environments with divergent climates. Thus, historical and contemporary geography and historical and contemporary environments have all shaped patterns of genetic structure in E. henryi, and, in fact, changes in the landscape through time have also been critical factors.

Evolution of East Asia's Arcto-Tertiary relict Euptelea (Eupteleaceae) shaped by Late Neogene vicariance and Quaternary climate change.

BACKGROUND: The evolutionary origin and historical demography of extant Arcto-Tertiary forest species in East Asia is still poorly understood. Here, we reconstructed the evolutionary and population demographic history of the two extant Euptelea species in China (E. pleiosperma) and Japan (E. polyandra). Chloroplast/nuclear DNA sequences and microsatellite loci were obtained from 36 Euptelea populations to explore molecular structure and diversity in relation to past and present distributions based on ecological niche modelling (ENM). Time-calibrated phylogenetic/phylogeographic inferences and niche-identity tests were used to infer the historical process of lineage formation. RESULTS: Euptelea pleiosperma diverged from E. polyandra around the Late Miocene and experienced significant ecological differentiation. A near-simultaneous diversification of six phylogroups occurred during the mid-to-late Pliocene, in response to the abrupt uplift of the eastern Tibetan Plateau and an increasingly cooler and drier climate. Populations of E. pleiosperma seem to have been mostly stationary through the last glacial cycles, while those of E. polyandra reflect more recent climate-induced cycles of range contraction and expansion. CONCLUSIONS: Our results illustrate how Late Neogene climatic/tectonic changes promoted speciation and lineage diversification in East Asia's Tertiary relict flora. They also demonstrate for the first time a greater variation in such species' responses to glacial cycles in Japan when compared to congeners in China.