Ecotype variation in the endemic tree Callicarpa subpubescens on small oceanic islands: genetic, phenotypic, and environmental insights.

Callicarpa subpubescens, endemic to the Ogasawara Islands, is suggested to have multiple ecotypes in the Hahajima Islands, specifically in the central part of the Ogasawara Islands. In this study, associations between genetic groups and spatial distribution, habitat, leaf morphology, size structure, and flowering time of each genetic group were investigated on Hahajima and the satellite Imoutojima Islands. Genetic groups were identified using EST-SSR markers, revealing four ecotypes named based on morphological features: Dwarf (D), Glabrescent (G), Tall (T), and Middle (M), with M being a result of the hybridization of G and T. Ecotype D, adapted to dry environments, is characterized by small tree size, dense thick leaves with abundant hairs, and is distributed in dry scrub. Ecotype G, adapted to understory of mesic forests, lacks leaf hairs. Ecotype T, adapted to the canopy of mesic forests, has hairy leaves and is tall in tree height. Ecotype M, adapted to the canopy of mesic scrub or edges of mesic forests, has hairy leaves but with a shorter tree height than ecotype T. Flowering peaks differed among all ecotype pairs except G and M, but the flowering times more or less overlapped among all ecotypes, suggesting that pre-mating isolation among ecotypes is not perfect. Post-mating isolation is considered absent, as there were no differences in the results, germination, and survival rates of one-year seedlings among inter- and intra-ecotype crossings. The existence of such ecotypes provides valuable insights into the ongoing speciation processes adapting to the oceanic island environments.

Genetic differentiation and evolution of broad-leaved evergreen shrub and tree varieties of Daphniphyllum macropodum (Daphniphyllaceae).

Tree form evolution is an important ecological specialization for woody species, but its evolutionary process with adaptation is poorly understood, especially on the microevolutionary scale. Daphniphyllum macropodum comprises two varieties: a tree variety growing in a warm temperate climate with light snowfall and a shrub variety growing in a cool temperate climate with heavy snowfall in Japan. Chloroplast DNA variations and genome-wide single-nucleotide polymorphisms across D. macropodum populations and D. teijsmannii as an outgroup were used to reveal the evolutionary process of the shrub variety. Population genetic analysis indicated that the two varieties diverged but were weakly differentiated. Approximate Bayesian computation analysis supported a scenario that assumed migration between the tree variety and the southern populations of the shrub variety. We found migration between the two varieties where the distributions of the two varieties are in contact, and it is concordant with higher tree height in the southern populations of the shrub variety than the northern populations. The genetic divergence between the two varieties was associated with snowfall. The heavy snowfall climate is considered to have developed since the middle Quaternary in this region. The estimated divergence time between the two varieties suggests that the evolution of the two varieties may be concordant with such paleoclimatic change.

Genetic and demographic signatures accompanying the evolution of the selfing syndrome in Daphne kiusiana, an evergreen shrub.

BACKGROUND AND AIMS: The evolution of mating systems from outcrossing to self-fertilization is a common transition in flowering plants. This shift is often associated with the 'selfing syndrome', which is characterized by less visible flowers with functional changes to control outcrossing. In most cases, the evolutionary history and demographic dynamics underlying the evolution of the selfing syndrome remain poorly understood. METHODS: Here, we characterize differences in the demographic genetic consequences and associated floral-specific traits between two distinct geographical groups of a wild shrub, Daphne kiusiana, endemic to East Asia; plants in the eastern region (southeastern Korea and Kyushu, Japan) exhibit smaller and fewer flowers compared to those of plants in the western region (southwestern Korea). Genetic analyses were conducted using nuclear microsatellites and chloroplast DNA (multiplexed phylogenetic marker sequencing) datasets. KEY RESULTS: A high selfing rate with significantly increased homozygosity characterized the eastern lineage, associated with lower levels of visibility and herkogamy in the floral traits. The two lineages harboured independent phylogeographical histories. In contrast to the western lineage, the eastern lineage showed a gradual reduction in the effective population size with no signs of a severe bottleneck despite its extreme range contraction during the last glacial period. CONCLUSIONS: Our results suggest that the selfing-associated morphological changes in D. kiusiana are of relatively old origin (at least 100 000 years ago) and were driven by directional selection for efficient self-pollination. We provide evidence that the evolution of the selfing syndrome in D. kiusiana is not strongly associated with a severe population bottleneck.

The evolutionary history of rice azaleas (Rhododendron tschonoskii alliance) involved niche evolution to a montane environment.

PREMISE: The formation of isolated montane geography on islands promotes evolution, speciation, and then radiation if there are ecological changes. Thus, investigating evolutionary histories of montane species and associated ecological changes may help efforts to understand how endemism formed in islands' montane floras. To explore this process, we investigated the evolutionary history of the Rhododendron tschonoskii alliance, which grows in montane environments of the Japanese archipelago and the Korean Peninsula. METHODS: We studied the five species in the R. tschonoskii alliance and 30 outgroup species, using genome-wide single-nucleotide polymorphisms and cpDNA sequences, in association with environmental analyses. RESULTS: The monophyletic R. tschonoskii alliance diverged since the late Miocene. Species in the alliance currently inhabit a cold climatic niche that is largely different from that of the outgroup species. We observed clear genetic and niche differentiations between the taxa of the alliance. CONCLUSIONS: The association of the alliance's evolution with the formation of cooler climates on mountains indicates that it was driven by global cooling since the mid-Miocene and by rapid uplift of mountains since the Pliocene. The combination of geographic and climatic isolation promoted high genetic differentiation between taxa, which has been maintained by climatic oscillations since the Quaternary.

Different population size change and migration histories created genetic diversity of three oaks in Tokai region, central Japan.

To understand genetic diversity in focal species, it is important to consider the possibility of speciation with gene flow, especially in species with porous genomes such as oaks. We studied genetic diversity and structure in three oaks, Quercus mongolica var. mongolicoides (QM), Q. mongolica var. crispula (QC) and Q. serrata (QS), growing in the Tokai region, central Japan. QM is semi-endemic to the region while the others are common taxa. We also conducted demographic modeling to infer their population size change and migration histories using an approximate Bayesian computation (ABC) approach. The three taxa showed distinct genetic structures but there was genetic admixture among the taxa, especially between QM and QC. ABC analysis of population size change revealed that the population size of QM was stable during and after the last glacial period, while QC and QS showed population expansion after the last glacial maximum. ABC analysis of population divergence and migration revealed that continuous gene flow between QM and QC after their divergence was supported, while between QM and QS, and between QC and QS, secondary contact after sufficient isolation was supported. These historical migration patterns among the three taxa indicate that QM and QC are currently in the early stage or gray zone of speciation, whereas speciation of the other two taxon pairs is considered to have almost been established. Observed gene flow patterns and strength between QM and QC, and between QM and QS, were explained by both flowering patterns and historical distributions, but those between QC and QS were not.

Comparative Mitogenomic Analysis Reveals Gene and Intron Dynamics in Rubiaceae and Intra-Specific Diversification in Damnacanthus indicus.

The dynamic evolution of mitochondrial gene and intron content has been reported across the angiosperms. However, a reference mitochondrial genome (mitogenome) is not available in Rubiaceae. The phylogenetic utility of mitogenome data at a species level is rarely assessed. Here, we assembled mitogenomes of six Damnacanthus indicus (Rubiaceae, Rubioideae) representing two varieties (var. indicus and var. microphyllus). The gene and intron content of D. indicus was compared with mitogenomes from representative angiosperm species and mitochondrial contigs from the other Rubiaceae species. Mitogenome structural rearrangement and sequence divergence in D. indicus were analyzed in six individuals. The size of the mitogenome in D. indicus varied from 417,661 to 419,435 bp. Comparing the number of intact mitochondrial protein-coding genes in other Gentianales taxa (38), D. indicus included 32 genes representing several losses. The intron analysis revealed a shift from cis to trans splicing of a nad1 intron (nad1i728) in D. indicus and it is a shared character with the other four Rubioideae taxa. Two distinct mitogenome structures (type A and B) were identified. Two-step direct repeat-mediated recombination was proposed to explain structural changes between type A and B mitogenomes. The five individuals from two varieties in D. indicus diverged well in the whole mitogenome-level comparison with one exception. Collectively, our study elucidated the mitogenome evolution in Rubiaceae along with D. indicus and showed the reliable phylogenetic utility of the whole mitogenome data at a species-level evolution.

Genetic diversity, structure, and demography of Pandanus boninensis (Pandanaceae) with sea drifted seeds, endemic to the Ogasawara Islands of Japan: Comparison between young and old islands.

Pandanus boninensis, endemic to the Ogasawara Islands, Japan, is distributed on both the older Bonin and younger Volcano Islands. In this study, we conducted population genetic analyses of P. boninensis on these islands to examine the population diversity and structure across old and young islands, to assess potential differences in population demography with island age, and to collect any evidence of migration between old and young islands. We found that the genetic diversity of expressed sequence tag (EST)-based microsatellite (SSR) markers, the nucleotide diversity of nuclear DNA sequences, and the haplotype diversity of chloroplast DNA on young islands were lower than those on old islands. Clustering analyses of EST-SSR indicated that populations on old islands were strongly diverged from those on young islands. Approximate Bayesian computation analysis of EST-SSR suggested that population expansion occurred on old islands while population reduction occurred on young islands. We also found evidence of migration among old islands (mostly from south to north), while it appears that there have been very few migration events between old and young islands. These differences could be due to the fact that young islands tend to be geographically isolated and support smaller populations that began a shorter time ago from limited founders. The P. boninensis populations on the Volcano Islands are interesting from an evolutionary perspective as they constitute a classic example of the early stages of progressive colonization on oceanic islands with small effective population sizes and low genetic diversity.

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.

Prognostic index consisting of early post-transplant variables <2 weeks in adult living-donor liver transplantation.

AIM: Six-month recipient mortality after adult-to-adult living-donor liver transplantation (LDLT) remains high. Early and accurate prediction of recipient outcome and continuous monitoring of recipient severity after surgery are both essential for guiding appropriate care. This study was designed to identify early post-transplant parameters associated with 6-month mortality, and thereby to construct a discriminatory prognostic index (PI). METHODS: We retrospectively analyzed 400 consecutive primary adult-to-adult LDLTs in our center (2006-2017). Perioperative variables were comprehensively analyzed for their accuracy in predicting recipient mortality by comparing the area under the receiver operating characteristic (AUROC) of each factor. RESULTS: The AUROCs of preoperative predictive factors, for example, Model for End-stage Liver Disease (MELD) score and donor age, were 0.56 and 0.64, respectively, whereas those of post-transplant platelet count (PLT), total bilirubin (T-BIL), and prothrombin time - international normalized ratio (INR) on postoperative day (POD)-7-14 were 0.71/0.84, 0.68/0.82, and 0.71/0.78, respectively. Logistic regression analysis provided a formula: PIPOD-14 = 3.39 + 0.12 × PLTPOD-14 - 0.09 × T-BILPOD-14 - 1.23 × INRPOD-14 , indicating a high AUROC of 0.87. Recipient 6-month survival with PIPOD-14 < 2.38 (n = 173) was 71.7%, whereas that with PIPOD-14 ≥ 2.38 (n = 222) was 97.7% (P < 0.001). The AUROCs of PIPOD-7 were as high as 0.8 in the subgroups with younger donors (<50 years of age), right lobe grafts, ABO-identical/compatible combinations, or low MELD score (<20), indicating usefulness of PI to identify unexpectedly complicated cases within the first week. CONCLUSIONS: A novel, post-transplant survival estimator, PI, accurately predicts recipient 6-month mortality within 1-2 weeks after adult LDLT. Daily monitoring of PI could facilitate early interventions including retransplantation in critically ill patients.

Correction to: Population genetic structure and demography of Magnolia kobus: variety borealis is not supported genetically.

The original article was updated.

[A Case Report of a Difficult to Diagnose, Extremely Well-Differentiated Stomach Adenocarcinoma Diagnosed by Biopsy].

A Japanese man in his 60s presented with complaints of epigastric pain and weight loss. A gastrointestinal endoscopy revealed multiple gastric ulcers and an irregular mound located on the wall of the lower gastric body along the greater curvature, which was suspected to be cancerous. A biopsy revealed that it was a Group 2 tumor even though the biopsy was repeated 4 times. He was referred to our hospital and 3 biopsies were performed. The final result classified the tumor as Group 4. The patient underwent surgery and the pathological examination revealed an extremely well-differentiated adenocarcinoma( EWDA). An EWDA is characterized by a well-formed mucinous gland with little or no nuclear atypia, which makes preoperative biopsy diagnosis difficult.

Population genetic structure and demography of Magnolia kobus: variety borealis is not supported genetically.

Species delimitations by morphological and by genetic markers are not always congruent. Magnolia kobus consists of two morphologically different varieties, kobus and borealis. The latter variety is characterized by larger leaves than the former. For the conservation of M. kobus genetic resources in natural forests, the relationships between morphological and genetic variation should be clarified. We investigated variations in nuclear microsatellites, chloroplast DNA (cpDNA) sequences and leaf morphological traits in 23 populations of M. kobus over the range of species. Two genetically divergent lineages, northern and southern were detected and their geographical boundary was estimated to be at 39°N. The northern lineage consisted of two genetic clusters and a single cpDNA haplotype, while the southern one had multiple genetic clusters and cpDNA haplotypes. The northern lineage showed significantly lower genetic diversity than the southern. Approximate Bayesian computation indicated that the northern and southern lineages had experienced, respectively, population expansion and long-term stable population size. The divergence time between the two lineages was estimated to be 565,000 years ago and no signature of migration between the two lineages after divergence was detected. Ecological niche modeling showed that the potential distribution area in northern Japan at the last glacial maximum was very small. It is thus considered that the two lineages have experienced different population histories over several glacial-inter-glacial cycles. Individuals of populations in the central to northern part of Honshu on the Sea of Japan side and in Hokkaido had large leaf width and area. These leaf characteristics corresponded with those of variety borealis. However, the delimitation of the northern and southern lineages detected by genetic markers (39°N) was not congruent with that detected by leaf morphologies (36°N). It is therefore suggested that variety borealis is not supported genetically and the northern and southern lineages should be considered separately when identifying conservation units based not on morphology but on genetic markers.

Impact of Donor Age on Recipient Survival in Adult-to-Adult Living-donor Liver Transplantation.

OBJECTIVE: To investigate the influence of donor age on recipient outcome after living-donor partial liver transplantation (LDLT). BACKGROUND: Donor age is a well-known prognostic factor in deceased donor liver transplantation; however, its role in LDLT remains unclear. METHODS: We retrospectively analyzed 315 consecutive cases of primary adult-to-adult LDLT in our center between April 2006 and March 2014. Recipients were divided into 5 groups according to the donor age: D-20s (n = 60); D-30s (n = 72); D-40s (n = 57); D-50s (n = 94); and D-60s (n = 32). The recipient survival and the association with various clinical factors were investigated. RESULTS: Recipient survival proportions were significantly higher in D-20s compared with all the other groups (P = 0.008, < 0.001, < 0.001, and = 0.006, vs D-30s, -40s, -50s, and -60s, respectively), whereas there was no association between recipient survival and their own age. There are 3 typical relationships between donors and recipients in adult-to-adult LDLT: from child-to-parent, between spouses/siblings, and from parent-to-child. The overall survival in child-to-parent was significantly higher than in spouses/siblings (P = 0.002) and in parent-to-child (P = 0.005), despite significantly higher recipient age in child-to-parent [59 (42-69) years, P < 0.001]. Contrastingly, parent-to-child exhibited the lowest survival, despite the youngest recipient age [26 (20-43) years, P < 0.001]. In addition, younger donor age exhibited significantly better recipient survival both in hepatitis C virus-related and in non-hepatitis C virus diseases. Univariate and multivariate analyses both demonstrated that donor age and graft-type (right-sided livers) are independent prognostic factors for recipient survival. CONCLUSIONS: Donor age is an independent, strong prognostic factor in adult-to-adult LDLT.

Hydrogen Flush After Cold Storage as a New End-Ischemic Ex Vivo Treatment for Liver Grafts Against Ischemia/Reperfusion Injury.

Cold storage (CS) remains the gold standard for organ preservation worldwide, although it is inevitably associated with ischemia/reperfusion injury (IRI). Molecular hydrogen (H2 ) is well known to have antioxidative properties. However, its unfavorable features, ie, inflammability, low solubility, and high tissue/substance permeability, have hampered its clinical application. To overcome such obstacles, we developed a novel reconditioning method for donor organs named hydrogen flush after cold storage (HyFACS), which is just an end-ischemic H2 flush directly to donor organs ex vivo, and, herein, we report its therapeutic impact against hepatic IRI. Whole liver grafts were retrieved from Wistar rats. After 24-hour CS in UW solution, livers were cold-flushed with H2 solution (1.0 ppm) via the portal vein (PV), the hepatic artery (HA), or both (PV + HA). Functional integrity and morphological damages were then evaluated by 2-hour oxygenated reperfusion at 37°C. HyFACS significantly lowered portal venous pressure, transaminase, and high mobility group box protein 1 release compared with vehicle-treated controls (P < 0.01). Hyaluronic acid clearance was significantly higher in the HyFACS-PV and -PV + HA groups when compared with the others (P < 0.01), demonstrating the efficacy of the PV route to maintain the sinusoidal endothelia. In contrast, bile production and lactate dehydrogenase leakage therein were both significantly improved in HyFACS-HA and -PV + HA (P < 0.01), representing the superiority of the arterial route to attenuate biliary damage. Electron microscopy consistently revealed that sinusoidal ultrastructures were well maintained by portal HyFACS, while microvilli in bile canaliculi were well preserved by arterial flush. As an underlying mechanism, HyFACS significantly lowered oxidative damages, thus improving the glutathione/glutathione disulfide ratio in liver tissue. In conclusion, HyFACS significantly protected liver grafts from IRI by ameliorating oxidative damage upon reperfusion in the characteristic manner with its route of administration. Given its safety, simplicity, and cost-effectiveness, end-ischemic HyFACS may be a novel pretransplant conditioning for cold-stored donor organs.

The origin of wild populations of Toxicodendron succedaneum on mainland Japan revealed by genetic variation in chloroplast and nuclear DNA.

Toxicodendron succedaneum: (L.) Kuntze is a tree cultivated for the production of sumac wax, which is extracted from the mesocarp. There are several hypotheses regarding the origin of T. succedaneum on mainland Japan. In this study, the geographical distribution of genetic variation in 13 wild populations on Honshu, Shikoku, Kyushu, and Ryukyu Islands, Japan was investigated and compared with that of individuals from continental Asia. Seven chloroplast DNA haplotypes of T. succedaneum were observed in Japan and could be divided into three lineages based on relatedness between haplotypes. One of these lineages was also observed in continental Asia, and the others were genetically distant from the haplotypes that originated on the continent, with one considered to have originated on the Ryukyu Islands, and the other on mainland Japan. The genetic diversity of both chloroplast and nuclear DNA was lower in populations from Ryukyu Islands than in populations from mainland Japan. Bayesian clustering based on nuclear genotypes showed a clear difference between the groups from Ryukyu Islands and mainland Japan. Based on approximate Bayesian computation analysis of polymorphic data for both genomes, it was inferred that wild populations of T. succedaneum on mainland Japan consist of both lineages with natural distribution on mainland Japan and those introduced from Ryukyu Islands and continental Asia.

Reduced incompatibility in the production of second generation hybrids between two Magnolia species revealed by Bayesian gene dispersal modeling.

PREMISE OF THE STUDY: Hybrid zones are areas where gene flow between related species is currently occurring, so information on the compatibility between related species and their hybrids is essential for predicting the dynamics of such zones generated by introgressive hybridization. In this study, we quantified the compatibility among Magnolia stellata, M. salicifolia, and their hybrids in a hybrid zone using gene dispersal modeling. METHODS: After determining the genealogical classes of adult trees in the hybrid zone, the paternity of 574 open-pollinated seeds from 37 known maternal trees was analyzed with microsatellite markers. A neighborhood-based Bayesian gene dispersal model developed by us for estimating compatibility was then applied to the paternity data. KEY RESULTS: When M. stellata or M. salicifolia were mothers, interspecific mating to produce F1 hybrids yielded significant incompatibility, but backcrossing with F1 hybrids did not. Furthermore, when F1 hybrids became mothers, no significant incompatibility resulted from backcrossing to parental species or intra-F1 mating to produce F2 hybrids. The estimated proportion of F1 hybrids in the outcrossed seeds (1.7%) in the hybrid zone was much lower than that in the adult trees (14.0%). CONCLUSIONS: While it is difficult to obtain F1 hybrids, their low incompatibility makes it easy to produce advanced generation hybrids, once they have been successfully obtained. Although the production of F1 seeds is rare, heterosis and/or weak selection pressure in an empty niche between the parental species' niches may have contributed to the increased proportion of adult F1 hybrids in the hybrid zone.

Genome sequencing of Metrosideros polymorpha (Myrtaceae), a dominant species in various habitats in the Hawaiian Islands with remarkable phenotypic variations.

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

Adaptive radiation of the Callicarpa genus in the Bonin Islands revealed through double-digest restriction site-associated DNA sequencing analysis.

The Bonin Islands, comprised of the Mukojima, Chichijima, and Hahajima Islands, are known for their isolated and distinctive habitats, hosting a diverse array of endemic flora and fauna. In these islands, adaptive radiation has played a remarkable role in speciation, particularly evident in the Callicarpa genus that is represented by three species: Callicarpa parvifolia and C. glabra exclusive to the Chichijima Islands, and Callicarpa subpubescens, distributed across the entire Bonin Islands. Notably, C. subpubescens exhibits multiple ecotypes, differing in leaf hair density, flowering time, and tree size. In this study, we aimed to investigate species and ecotype diversification patterns, estimate divergence times, and explore cryptic species within Callicarpa in the Bonin Islands using phenotypic and genetic data (double-digest restriction site-associated DNA sequencing). Genetic analysis revealed that C. parvifolia and C. glabra both formed single, distinct genetic groups. Conversely, C. subpubescens consisted of six genetic groups corresponding to different ecotypes and regions, and a hybrid group resulting from the hybridization between two of these genetic groups. Population demography analysis focusing on six Chichijima and Hahajima Islands-based species/ecotypes indicated that all species and ecotypes except one ecotype diverged simultaneously around 73-77 kya. The star-shaped neighbor-net tree also suggests the simultaneous divergence of species and ecotypes. The species and ecotypes that simultaneously diverged adapted to dry environments and understory forests, suggesting that aridification may have contributed to this process of adaptive radiation. Moreover, leaf morphology, flowering time, and genetic analyses suggested the presence of two cryptic species and one hybrid species within C. subpubescens.

Phylogenetic, population structure, and population demographic analyses reveal that Vicia sepium in Japan is native and not introduced.

Vicia sepium (bush vetch) is a perennial legume widely distributed throughout the Eurasian continent. However, its distribution in Japan is limited to Mt. Ibuki and small parts of central and southern Hokkaido. Therefore, each Japanese V. sepium lineage has been considered to have been introduced separately from Europe. Here, we examined whether the species was introduced or not on the basis of cpDNA sequences and genome-wide SNPs from Japanese and overseas samples. Both the cpDNA haplotype network and the nuclear DNA phylogenetic tree showed that Japanese V. sepium is monophyletic. Furthermore, although the nuclear DNA phylogenetic tree also showed that each lineage is clearly monophyletic, genetic admixture of the genetic cluster dominated in the Hokkaido lineage was also detected in the Mt. Ibuki lineage. Population divergence analysis showed that the two lineages diverged during the last glacial period. The Mt. Ibuki lineage showed a sudden population decline 300-400 years ago, indicating that some anthropogenic activity might be involved, while the Hokkaido lineage showed a gradual population decline from 5000 years ago. Consequently, these two lineages show low current genetic diversity compared with overseas lineages. These results show that the Japanese V. sepium is not introduced but is native.

The genetic structure and demographic history of Zabelia tyaihyonii, endemic to Korean limestone karst forests, based on genome-wide SNP markers.

Similar to the global phenomenon, many plant species endemic to Korean limestone karst forests are at risk of extinction due to human intervention. Zabelia tyaihyonii is a familiar shrub, called "Hardy abelia" and "Fragrant abelia" growing in the karst forests of Korea, where it is one of the most threatened species. We investigated the genetic structure and demographic history of Z. tyaihyonii, which allow us to develop appropriate conservation and management strategies. The genetic structure was evaluated using a total of 187 samples from 14 populations, covering the entire distribution of Z. tyaihyonii in South Korea. We utilized 254 and 1753 SNP loci obtained via MIG-seq (Multiplexed ISSR Genotyping by sequencing) for structure and demographic analyses, respectively. The population demographic modeling was performed with site frequency spectrum. To gain further historical insights, we also employed ENM (Ecological Niche Modeling). We found two distinct clusters (CLI and CLII) of ancient origin (ca. 490 ka). Despite CLII experiencing a more severe bottleneck, both clusters showed similar levels of genetic diversity, indicating mutual historical gene flow. Their historical distribution range seems to have changed very little. We proposed a historical distribution scenario for Z. tyaihyonii, taking into account its intrinsic factors, and emphasized a more complex response to Quaternary climate change beyond simple allopatric speciation models. These findings provide valuable insights for conservation and management strategies for Z. tyaihyonii.