Natural selection on PHYE by latitude in the Japanese archipelago: insight from locus specific phylogeographic structure in Arcterica nana (Ericaceae)

Decades-long phylogeographic issues: complex historical processes and ecological factors on genetic structure of alpine plants in the Japanese Archipelago

Mountain regions are important places for biodiversity, where organisms could persist throughout prolonged periods and accumulate genetic divergence as well as promote speciation. Roles of mountains for biodiversity have been exclusively discussed in regions that have specifically diverse species or covered with ice-sheets during the Pleistocene glacial periods, whereas the importance of mountainous regions in East Asia has been less disputed. High mountains in the Japanese Archipelago, located at the eastern edge of the Eurasia continent, have one of southernmost populations of alpine and arctic-alpine plants that are also distributed in the northern Pacific and/or the circumarctic regions. Phylogeographic studies on the Japanese alpine plants have excluded their possible ephemeral occurrence during the current warm period, and rather, suggest persistence of alpine plants throughout several cycles of climate changes in the Pleistocene on high mountains in central Honshu, the main island of the Japanese Archipelago. In this review, I look through decade long phylogeographic studies and show complicated patterns of range dynamics of Japanese alpine plants. In addition, I note recent findings of genetic relationships of Japanese populations of alpine and/or arctic-alpine plants with those in northern regions and their possible ecological divergence in the Japanese Archipelago. Taken together, I provide several issues for understanding historical processes that established distribution of alpine plants following climate changes as well as their diversification and propose importance of Japanese populations of alpine plants on biodiversity in alpine communities across broader range, especially in the northern Pacific region.

Divergence in red light responses associated with thermal reversion of phytochrome B between high- and low-latitude species

Phytochromes play a central role in mediating adaptive responses to light and temperature throughout plant life cycles. Despite evidence for adaptive importance of natural variation in phytochromes, little information is known about molecular mechanisms that modulate physiological responses of phytochromes in nature. We show evolutionary divergence in physiological responses relevant to thermal stability of a physiologically active form of phytochrome (Pfr) between two sister species of Brassicaceae growing at different latitudes. The higher latitude species (Cardamine bellidifolia; Cb) responded more strongly to light-limited conditions compared with its lower latitude sister (C. nipponica; Cn). Moreover, CbPHYB conferred stronger responses to both light-limited and warm conditions in the phyB-deficient mutant of Arabidopsis thaliana than CnPHYB: that is Pfr CbphyB was more stable in nuclei than CnphyB. Our findings suggest that fine tuning Pfr stability is a fundamental mechanism for plants to optimise phytochrome-related traits in their evolution and adapt to spatially varying environments, and open a new avenue to understand molecular mechanisms that fine tune phytochrome responses in nature.

ddRAD-seq based phylogeographic study of Sargassum thunbergii (Phaeophyceae, Heterokonta) around Japanese coast

Using genome-wide SNP data obtained from high-throughput techniques based on double digest restriction site-associated DNA sequencing (ddRAD-seq), we elucidated the migration history and genetic diversity of the Japanese population of the ecologically important brown seaweed Sargassum thunbergii (Mertens ex Roth) Kuntze. STRUCTURE and NeighborNet analyses showed a clear genetic differentiation among populations of four geographic regions: Kyushu (POP1); Sea of Japan (POP2); Hokkaido and Tohoku (POP3); and Pacific coast from Kyushu to Kanto (POP4). Approximate Bayesian Computation (ABC) analysis indicated that POP4 diverged first, followed by the separation between POP2 (the largest effective population size) and POP3; POP1 was the last to form, shaped by the mixture of POP2 (73%) and POP4 (27%). High genetic diversity was detected in POP1 and POP2, whereas low genetic diversity was detected in POP3 and POP4. These results indicated that S. thunbergii populations of Kyushu and the Sea of Japan might have been maintained as large and stable populations gathered different lineages from China, Korea and/or Japan.

Late Pleistocene origin of the entire circumarctic range of the arctic-alpine plant Kalmia procumbens

The circumarctic ranges of arctic-alpine plants are thought to have been established in the late Pliocene/early Pleistocene, when the modern arctic tundra was formed in response to climate cooling. Previous findings of range-wide genetic structure in arctic-alpine plants have been thought to support this hypothesis, but few studies have explicitly addressed the temporal framework of the genetic structure. Here, we estimated the demographic history of the genetic structure in the circumarctic Kalmia procumbens using sequences of multiple nuclear loci and examined whether its genetic structure reflects prolonged isolation throughout the Pleistocene. Both Bayesian clustering and phylogenetic analyses revealed genetic distinction between alpine and arctic regions, whereas detailed groupings were somewhat discordant between the analyses. By assuming a population grouping based on the phylogenetic analyses, which likely reflects a deeper intraspecific divergence, we conducted model-based analyses and demonstrated that the intraspecific genetic divergence in K. procumbens likely originated during the last glacial period. Thus, there is no need to postulate range separation throughout the Pleistocene to explain the current genetic structure in this species. This study demonstrates that range-wide genetic structure in arctic-alpine plants does not necessarily result from the late Pliocene/early Pleistocene origin of their circumarctic ranges and emphasizes the importance of a temporal framework of the current genetic structure for understanding the biogeographic history of the arctic flora.

Population demographic history of a temperate shrub, Rhododendron weyrichii (Ericaceae), on continental islands of Japan and South Korea

Continental islands provide opportunities for testing the effects of isolation and migration on genetic variation in plant populations. In characteristic of continental islands is that the geographic connections between these islands, which are currently distinguished by seaways, have experienced fluctuations caused by sea‐level changes due to climate oscillations during the Quaternary. Plant populations on the islands have migrated between these islands via the exposed seafloors or been isolated. Here, we examined the demographic history of a temperate shrub, Rhododendron weyrichii, which is distributed in the southwestern parts of the Japanese archipelago and on an island of South Korea, using statistical phylogeographic approaches based on the DNA sequences of two chloroplast and eight nuclear loci in samples analyzed from 18 populations on eight continental islands, and palaeodistribution modeling. Time estimates for four island populations indicate that the durations of vicariance history are different between these populations, and these events have continued since the last glacial or may have predated the last glacial. The constancy or expansion of population sizes on the Japanese islands, and in contrast a bottleneck in population size on the Korean island Jeju, suggests that these islands may have provided different conditions for sustaining populations. The result of palaeodistribution modeling indicates that the longitudinal range of the species as a whole has not changed greatly since the last glacial maximum. These results indicate that exposed seafloors during the glacial period formed both effective and ineffective migration corridors. These findings may shed light on the effects of seafloor exposure on the migration of plants distributed across continental islands.

Importance of demographic history for phylogeographic inference on the arctic-alpine plant Phyllodoce caerulea in East Asia

Arctic-alpine plants have enormous ranges in the Northern Hemisphere. Phylogeographic studies have provided insights into their glacial survival as well as their postglacial colonization history. However, our understanding of the population dynamics of disjunct alpine populations in temperate regions remains limited. During Pleistocene cold periods, alpine populations of arctic-alpine species in East Asia were either connected to an ice-free Beringia refugium or they persisted with prolonged isolation after their establishment. To estimate which of these scenarios is more likely, we elucidated the genetic structure of Phyllodoce caerulea (Ericaceae) in Beringia and northern Japan, East Asia. Sequence variation in multiple nuclear loci revealed that P. caerulea can be distinguished into northern and southern groups. A demographic analysis demonstrated that the north-south divergence did not predate the last glacial period and detected introgression from Phyllodoce aleutica, relative widely distributed in East Asia, exclusively into the southern group. Therefore, although there has been genetic divergence between northern Japan and Beringia in P. caerulea, the divergence is unlikely to have resulted from their prolonged geographic separation throughout several cycles of glacial and interglacial periods. Instead, our study suggests that the introgression contributed to the genetic divergence of P. caerulea and that the range of P. caerulea was plausibly connected between northern Japan and Beringia during the last glacial period. Overall, our study not only provides a biogeographic insight into alpine populations of arctic-alpine plants in East Asia but also emphasizes the importance of careful interpretation of genetic structure for inferring phylogeographic history.

Adaptive Molecular Evolution of PHYE in Primulina, a Karst Cave Plant

Limestone Karst areas possess high levels of biodiversity and endemism. Primulina is a typical component of Karst endemic floras. The high species richness and wide distribution in various Karst microenvironments make the genus an idea model for studying speciation and local adaptation. In this study, we obtained 10 full-length sequences of the phytochrome PHYE from available transcriptome resources of Primulina and amplified partial sequences of PHYE from the genomic DNA of 74 Primulina species. Then, we used maximum-likelihood approaches to explore molecular evolution of PHYE in this Karst cave plant. The results showed that PHYE was dominated by purifying selection in both data sets, and two sites were identified as potentially under positive selection. Furthermore, the ω ratio varies greatly among different functional domains of PHYE and among different species lineages. These results suggest that potential positive selection in PHYE might have played an important role in the adaption of Primulina to heterogeneous light environments in Karst regions, and different species lineages might have been subjected to different selective pressures.

Persistent history of the bird-dispersed arctic-alpine plant Vaccinium vitis-idaea L. (Ericaceae) in Japan

Arctic-alpine plants have expanded and contracted their ranges in response to the Pleistocene climate oscillations. Today, many arctic-alpine plants have vast distributions in the circumarctic region as well as marginal, isolated occurrences in high mountains at lower latitudes. These marginal populations may represent relict, long-standing populations that have persisted for several cycles of cold and warm climate during the Pleistocene, or recent occurrences that either result from southward step-wise migration during the last glacial period or from recent long-distance dispersal. In light of these hypotheses, we investigated the biogeographic history of the marginal Japanese populations of the widespread arctic-alpine plant Vaccinium vitis-idaea (Ericaceae), which is bird-dispersed, potentially over long distances. We sequenced three nuclear loci and one plastid DNA region in 130 individuals from 65 localities covering its entire geographic range, with a focus on its marginal populations in Japan. We found a homogenous genetic pattern across its enormous range based on the loci analysed, in contrast to the geographically structured variation found in a previous study of amplified fragment length polymorphisms in this species. However, we found several unique haplotypes in the Japanese populations, excluding the possibility that these marginal populations result from recent southward migration. Thus, even though V. vitis-idaea is efficiently dispersed via berries, our study suggests that its isolated populations in Japan have persisted during several cycles of cold and warm climate during the Pleistocene.

Molecular evidence for ancient relicts of arctic-alpine plants in East Asia

Following climate cooling at the end of the Tertiary, arctic-alpine plants attained most of their extant species diversity. Because East Asia was not heavily glaciated, the importance of this region as a location for the long-term persistence of these species and their subsequent endemism during the Pleistocene was proposed in early discussions of phytogeography. However, this hypothesis remains to be fully tested. Here, we address this hypothesis by elucidating the phylogenetic history of Phyllodoce (Ericaceae). A phylogenetic tree based on multiple nuclear loci revealed that Phyllodoce nipponica was not derived from widespread species such as the arctic-alpine Phyllodoce caerulea, but rather represented an independent lineage sister to the clade of widespread relatives. Molecular dating indicated a mid-Pleistocene divergence of P. nipponica. These findings exclude the hypothesis that P. nipponica was derived from an arctic-alpine species that extended its range southwards during recent glacial periods. Instead, our results support the hypothesis that P. nipponica is an ancestral species which persisted in the Japanese archipelago during the mid- and late Pleistocene. Our findings demonstrate support for the early proposal and shed light on the importance of the Japanese archipelago for the evolution and persistence of arctic-alpine species.

Natural selection and neutral evolution jointly drive population divergence between alpine and lowland ecotypes of the allopolyploid plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074-0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (FST  = 0.041-0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.

Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus

BACKGROUND

Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance.

RESULTS

Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of up-regulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations.

CONCLUSIONS

Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.

Phylogenetic utility of the nuclear genes AGAMOUS 1 and PHYTOCHROME B in palms (Arecaceae): an example within Bactridinae

BACKGROUND AND AIMS

Molecular phylogenetic studies of palms (Arecaceae) have not yet provided a fully resolved phylogeny of the family. There is a need to increase the current set of markers to resolve difficult groups such as the Neotropical subtribe Bactridinae (Arecoideae: Cocoseae). We propose the use of two single-copy nuclear genes as valuable tools for palm phylogenetics.

METHODS

New primers were developed for the amplification of the AGAMOUS 1 (AG1) and PHYTOCHROME B (PHYB) genes. For the AGAMOUS gene, the paralogue 1 of Elaeis guineensis (EgAG1) was targeted. The region amplified contained coding sequences between the MIKC K and C MADS-box domains. For the PHYB gene, exon 1 (partial sequence) was first amplified in palm species using published degenerate primers for Poaceae, and then specific palm primers were designed. The two gene portions were sequenced in 22 species of palms representing all genera of Bactridinae, with emphasis on Astrocaryum and Hexopetion, the status of the latter genus still being debated.

KEY RESULTS

The new primers designed allow consistent amplification and high-quality sequencing within the palm family. The two loci studied produced more variability than chloroplast loci and equally or less variability than PRK, RPBII and ITS nuclear markers. The phylogenetic structure obtained with AG1 and PHYB genes provides new insights into intergeneric relationships within the Bactridinae and the intrageneric structure of Astrocaryum. The Hexopetion clade was recovered as monophyletic with both markers and was weakly supported as sister to Astrocaryum sensu stricto in the combined analysis. The rare Astrocaryum minus formed a species complex with Astrocaryum gynacanthum. Moreover, both AG1 and PHYB contain a microsatellite that could have further uses in species delimitation and population genetics.

CONCLUSIONS

AG1 and PHYB provide additional phylogenetic information within the palm family, and should prove useful in combination with other genes to improve the resolution of palm phylogenies.