Small islands represent a common feature in the Mediterranean and host a significant fraction of its biodiversity. However, the distribution of plant species richness across spatial scales—from local communities (alpha) to whole islands (gamma)—is largely unknown, and so is the influence of environmental, geographical, and topographical factors. By building upon classic biogeographic theory, we used the species–area relationship and about 4500 vegetation plots in 54 Central Mediterranean small islands to identify hotspots of plant species richness and the underlying spatial determinants across scales. To do so, we fitted and averaged eight species–area models on gamma and alpha richness against island area and plot size, respectively. Based on positive deviations from the fitted curves, we identified 12 islands as cross-scale hotspots. These islands encompassed around 70% of species and habitat richness, as well as almost 50% of the rarest species in the data set, while occupying less than 40% of the total island surface. By fitting generalized linear mixed models, we found that gamma richness was mainly explained by island area and was weakly related to mean annual temperature (positively) and annual precipitation (negatively). As for alpha richness, after accounting for the idiosyncratic effect of habitats and islands, plot size and gamma richness remained the only significant predictors, showing a positive relationship. This work contributes to the understanding of the patterns and drivers of plant diversity in Central Mediterranean small islands and outlines a useful methodology for the prioritization of conservation efforts.

China has the most numerous Danxia and Karst landscapes, which serve as special terrestrial islands harboring ample endemic species, though how did these endemic species spread among those isolated sites is still an unresolved issue. To address this question, we explored the phylogeographical structure and demographic history of Firmiana danxiaensis, a tree species endemic to Danxia and Karst landscapes. We collected 295 samples (28 populations) of F. danxiaensis. Plastid genomes were assembled for 25 representative samples. Sanger sequencing of four plastid regions and restriction-site-associated DNA sequencing were performed on the 28 populations. The phylogenetic tree constructed from plastid genomes and restriction site-associated DNA sequencing (RAD-seq) data supported that F. danxiaensis originated from Mount Danxia and Nanxiong Basin, spread to Karst landscapes near Yingde City, and then back to Danxia Mountain and the Nanxiong Basin. In the Nanxiong Basin, the latter arrivals captured the plastid of the former. Population analyses revealed strong population structure among and within Danxia and Karst landscapes, possibly due to low seed and pollen dispersal abilities of the species. The demographic and ecological niche modeling approaches suggested that F. danxiaensis have widely occurred in the southeast of China during the last glacial period, and later retreated to the cliffs of Danxia and Karst landscapes due to temperature rising and competition failure. The declining of the effective population size of the species throughout the postglacial period suggested that global warming, agriculture, and industrial civilizations could have affected the survival of this species, and more measures should be taken to conserve these species.

Parrotia C. A. Meyer (Hamamelidaceae) is a relictual genus with only two extant species disjunctly distributed in the subtropical forests of East China and temperate forests of North Iran. Fossil records suggest that Parrotia was widespread in Europe and Asia during the Miocene, but its fossils are predominantly based on pollen and leaves. In this paper we describe the first fossil flower of Parrotia based on an exceptionally well-preserved amber inclusion from the middle Miocene of Zhangpu, Southeast China. The fossil flower was investigated with light microscope and microcomputed tomography techniques. Parrotia zhiyanii sp. nov. is a small apetalous staminate flower subtended by a pair of prominent bracts. The androecium consists of 12 stamens, and each stamen consists of a short, slender filament and a prominent, elongated anther. The anthers with short simple trichomes on the outer surface and a prominent apical connective extension are opened by longitudinal slits. Unlike its living relatives, the new Parrotia from Zhangpu grew in a Miocene seasonal tropical rainforest.

Southeast Asia (SEA) has seen strong climatic oscillations and fluctuations in sea levels during the Quaternary. The impact of past climate changes on the evolution and distribution of local flora in SEA is still poorly understood. Here we aim to infer how the Quaternary climate change affects the evolutionary process and range shifts in two pine species. We investigated the population genetic structure and diversity using cytoplasmic DNA markers, and performed ecological niche modeling to reconstruct the species past distribution and to project range shift under future climates. We found substantial gene flow across the continuous distribution of the subtropical Pinus yunnanensis. In contrast, the tropical Pinus kesiya showed a strong population structure in accordance with its disjunct distribution across montane islands in Indochina and the Philippines. A broad hybrid zone of the two species occurs in southern Yunnan. Asymmetric introgression from the two species was detected in this zone with dominant mitochondrial gene flow from P. yunnanensis and chloroplast gene flow from P. kesiya. The observed population structure suggests a typical postglaciation expansion in P. yunnanensis, and a glacial expansion and interglacial contraction in P. kesiya. Ecological niche modeling supports the inferred demographic history and predicts a decrease in range size for P. kesiya under future climates. Our results suggest that tropical pine species in SEA have undergone evolutionary trajectories different from high latitude species related to their Quaternary climate histories. We also illustrate the need for urgent conservation actions in this fragmented landscape.

Chinese Tajiks are an Indo-Iranian-speaking population in Xinjiang, northwest China. Although the complex demographic history has been characterized, the ancestral sources and genetic admixture of Indo-Iranian-speaking groups in this region remain poorly understood. We here provide the genome-wide genotyping data for over 700 000 single-nucleotide polymorphisms (SNPs) and mtDNA multiplex sequencing data in 64 Chinese male Tajik individuals from two dialect groups, Wakhi and Selekur. We applied principal component analysis (PCA), ADMIXTURE, f-statistics, treemix, qpWave/qpAdm, Admixture-induced Linkage Disequilibrium for Evolutionary Relationships (ALDER), and Fst analyses to infer a fine-scale population genetic structure and admixture history. Our results reveal that Chinese Tajiks showed the closest affinity and similar genetic admixture pattern with ancient Xinjiang populations, especially Xinjiang samples in the historical era. Chinese Tajiks also have gene flow from European and Neolithic Iran farmers-related populations. We observed a genetic substructure in the two Tajik dialect groups. The Selekur-speaking group who lived in the county had more gene flow from East Asians than Wakhi-speaking people who inhabited the village. These results document the population movements contributed to the influx of diverse ancestries in the Xinjiang region.

Climate change is promoting global declines in plant diversity, which are expected to be more critical in islands or island-like ecosystems due to environmental constraints and isolation. The species' vulnerability to climate change (VUL) depends on their ability to cope with changes or mitigate them. Therefore, we investigate the influence of growth and dispersal strategies of species from the Sugarloaf Rock Complex, Brazil, an island-like ecosystem, on their niche breadth (NB), long-dispersal (LD) capacity, and geographical range (GR). Besides, we evaluate the potential use of these strategies as indicators of species' VUL. We found that rock specialists exhibit narrower NB, lower LD capacity, and a more restricted GR when compared to other species. We also found that 63% of rock specialists are found in conservation red-lists and they are more vulnerable to climate change than woody plants. Conversely, self-dispersed plants are expected to be less vulnerable to climate change when compared to species with other dispersal mechanisms. Species vulnerable to climate change are 14 times more likely to be included in conservation red lists, and it might indicate that the species' VUL might also describe the species' vulnerability to other anthropogenic threats. Still, we suggest conservation attention on some species that are expected to be vulnerable to climate change but were not yet included in conservation red lists. We advocate for more efforts to ensure the conservation aspects of different functional groups in which inselbergs might not only offer isolation but also a refuge opportunity.

Having a comprehensive understanding of genetic differentiation, responses to environmental change and demographic history is critical for genetic improvement and conservation efforts. Forest trees are an excellent resource for understanding population differentiation and adaptive genetic variation due to their ability to adapt to different climates and environments. Cephalotaxus oliveri is a relict conifer endemic to China. In this study, we generated transcriptome data and identified 17 728 high-quality single-nucleotide polymorphisms (SNPs) from 18 populations. We found significant negative correlations between expression diversity and nucleotide diversity within and among populations, suggesting that gene expression and nucleotide diversity have a reciprocal relationship when the species adapts to the environment. The analyses of population structure showed that C. oliveri displayed a striking genetic structure with four groups. BayeScEnv and RDA methods detected the signatures of local adaptation, and identified that 738 outlier SNPs were associated with precipitation, temperature and soil conditions across heterogeneous environmental conditions. Approximate Bayesian computation analyses showed that the first and second divergence occurred in the late Miocene (c. 10.075 million years ago [Ma]) and the middle Pleistocene transition (c. 0.815 Ma), respectively. Ecological niche modeling of C. oliveri revealed signs of westward expansion after the last glacial maximum, while it was predicted to experience significant range contractions in future climate change scenarios. Geographical factors and environmental factors in southern China have played a critical role in establishing the current genetic diversity and population structure of C. oliveri. This study provides an important reference for forest resource management and conservation for C. oliveri.

Rhododendron is the largest genus in Ericaceae and is well known for its diversity and beauty of flowers present in different species, making it a much-revered lineage of ornamental plants. Many species of Rhododendron are intolerant of high temperatures, which are becoming more common and intense in urban areas under global climate change. Therefore, the discovery and description of genes from heat-tolerant Rhododendron lineages are essential in the development of new climate-resilient cultivars. One such species known to be heat tolerant is Rhododendron × pulchrum Sweet. To better understand the genomics of heat tolerance in this species, we assembled a haplotype-resolved and chromosome-scale genome for R. × pulchrum, which had a genome size of 509 Mb; a scaffold N50 of 37 251 370 bp; and contained 35 610 genes. In addition, based on the same reannotation pipeline, we conducted pan-genomic analyses for all seven available chromosome-scale Rhododendron genomes and found 14 415 gene groups shared across all species and 18 018 gene groups distributed in the other species, including 1879 gene groups found in only a single species. Finally, we analyzed the transcriptomic data from heat-treated and non-heat-treated R. × pulchrum plants to quantify the genes that are most important during heat stress in an effort to inform the development of climate-resilient cultivars. This study provides insight into the genome diversity in Rhododendron and targets several genes related to agronomic traits that may help in further analysis.

Despite representing a fraction of the global terrestrial surface area, oceanic islands are disproportionately diverse in species, resulting from high rates of endemicity. Island plants are thought to share a unique phenotype—referred to as an island syndrome—which is thought to be driven by convergent evolution in response to selection by shared abiotic and biotic factors. One aspect of the island plant syndrome that has received relatively little research focus is that island plants are expected to have converged on conservative resource use associated with slow growth rates and weak competitive abilities. Here we tested whether native, woody Hawaiian plant species are phenotypically distinct—with more resource-conservative leaf traits—compared to a globally distributed sample of continental species. Using an archipelago-wide trait data set, we detected that on average, native Hawaiian species had lower leaf nutrient concentrations overall, and lower nutrient concentrations at high leaf mass per area, but no other phenotypic differences compared with continental plants. There was also considerable overlap in the trait spaces of native Hawaiian species and continental species. Our findings indicate that an island plant syndrome for leaf traits is not present in the Hawaiian flora, and that island species can demonstrate extensive variation in their resource-use strategies, on a scale that is comparable with that of continental species worldwide.

Clear species boundaries are crucial for plans and actions on biodiversity conservation. However, morphological similarities among allied species can result in taxonomic difficulties, thus impeding conservation efforts. In China, Cinnamomum japonicum Siebold is a well-known endangered plant, yet suffers from longstanding taxonomic issues. Here, we explicitly evaluate whether C. japonicum, C. chenii, and C. chekiangense are the same phylogenetic species on the basis of a multi-individual sampling strategy. We identified three sets of low-copy orthologous genes from 19 Lauraceae taxa for phylogenetic inferences. Both the concatenation and coalescent-based phylogenies supported that C. chenii individuals were embedded in the C. japonicum clade, indicating these two taxa are conspecific. Meanwhile, C. chekiangense accessions formed a monophyly which was not sister to C. japonicum. This result, together with the morphological differences that the leaves of C. japonicum are glabrous with a faveolate pattern of venation while those of C. chekiangense have trichomes and inevident lateral veins, led us to consider both as two distinct species. Based on 17 728 neutral single nucleotide polymorphisms (SNPs), the ADMIXTURE analysis suggested that the Chinese C. japonicum populations in Zhoushan Archipelago (=C. chenii) were genetically differentiated from the Japanese and Korean ones. Furthermore, ecological niche modeling predicted that the present distribution area of Chinese C. japonicum is likely to be unsuitable under global warming scenarios. Together with its limited distribution and genetic uniqueness, we recommend that Chinese C. japonicum deserves conservation priorities.

Polyploids are common in Camellia sect. Paracamellia, which contain many important oil crop species. However, their complex evolutionary history is largely unclear. In this study, 22 transcriptomes and 19 plastomes of related species of Camellia were sequenced and assembled, providing the most completed taxa sampling of Camellia sect. Oleifera and C. sect. Paracamellia. Phylogenetic trees were reconstructed with predicted single-copy nuclear genes and plastomes. Phylogenetic trees with nuclear genes demonstrated that C. sect. Oleifera should be merged into C. sect. Paracamellia. Cytonuclear discordance and network analyses suggested hybridizations among polyploid species and relatives. The divergence of major clades in C. sect. Paracamellia was dated to be during the middle to late Miocene from the ancestral Lingnan region, and a rapid diversification during the Quaternary was found, probably through hybridization and polyploidization. The tetraploid Camellia meiocarpa Hu may have originated from hybridization between closely related diploid species. The hexaploid Camellia oleifera C. Abel probably originated from hybridization between closely related diploid and tetraploid (e.g., C. meiocarpa) species. The octoploid Camellia vietnamensis T. C. Huang ex Hu could have originated from hybridization between hexaploid C. oleifera and the closely related diploid species. Hybridization and polyploidization played an important role in generating the rich variation of important fruit traits, especially increased fruit size in polyploid species.

Tsuga (hemlock) is a small genus of 10 extant species in the Pinaceae, with a disjunct distribution in East Asia and eastern and western North America. Reliable species-level identification of Tsuga fossils depends on the discovery of seed cones with intact bracts, but such cones are rare in the fossil record. Here we describe a new fossil species of hemlock as T. weichangensis sp. nov. based on exquisitely preserved seed cones with nearly complete bracts from the Lower Miocene of Weichang, Hebei Province, North China. This fossil species displays a mosaic of characters between Tsuga and Nothotsuga. The well-developed and slightly exserted bract scales of T. weichangensis are reminiscent of Nothotsuga, but other characters, such as nonleaved peduncles and tongue-shaped bract scales, in addition to monosaccate pollen found at the same fossil locality, suggest an affinity closer to Tsuga. Cladistic analysis based on 15 morphological characters and a molecular backbone constraint supports the assignment of these fossil cones to Tsuga rather than Nothotsuga, and places the fossil species of T. weichangensis in an unresolved polytomy within the genus Tsuga. The occurrence of Tsuga seed cone fossils indicate the paleoclimate in the Miocene of Weichang was warmer and more humid than today's climate, which is consistent with the paleoclimate reconstructed by paleopalynology.

The adaptive significance of phenotypic differences between females and males can provide insights into sex-specific selection and the evolution of sexual dimorphism. Dioecious plants commonly exhibit sexual dimorphism in height, although its ecological and evolutionary significance have rarely been examined experimentally. Here, we investigate the functional consequences of the temporal reversal of height dimorphism for pollen and seed dispersal in dioecious Rumex hastatulus Baldw., a species in which males are taller than females at flowering and the reverse pattern occurs at fruiting. Populations of this colonizing weed are wind-pollinated and seeds are wind-dispersed. In a glasshouse experiment we manipulated the height of pollen donors and using sex-specific genetic markers compared the paternal success of males of contrasting height and investigated whether seed families showed evidence of sexual dimorphism in early life-history traits. In a second glasshouse experiment using fruiting plants we also examined how female height influenced the distance that seeds were dispersed. We found that taller males had significantly higher siring success than males of equivalent height to flowering females. Similarly, taller females dispersed fruit to greater distances than shorter females. Female seeds were significantly heavier than male seeds and germinated more rapidly, although early seedling growth was greater in males. Our study suggests that the striking sex reversal of height in R. hastatulus likely functions to optimize the contrasting reproductive functions of the sexes by promoting increased pollen and seed dispersal distances. Improved dispersal quality could limit inbreeding and reduce local mate and resource competition within populations.

Habitat stability is important for maintaining biodiversity by preventing species extinction, but this stability is being challenged by climate change. The tropical alpine ecosystem is currently one of the ecosystems most threatened by global warming, and the flora close to the permanent snow line is at high risk of extinction. The tropical alpine ecosystem, found in South and Central America, Malesia and Papuasia, Africa, and Hawaii, is of relatively young evolutionary age, and it has been exposed to changing climates since its origin, particularly during the Pleistocene. Estimating habitat loss and gain between the Last Glacial Maximum (LGM) and the present allows us to relate current biodiversity to past changes in climate and habitat stability. In order to do so, (i) we developed a unifying climate-based delimitation of tropical alpine regions across continents, and (ii) we used this delimitation to assess the degree of habitat stability, that is, the overlap of suitable areas between the LGM and the present, in different tropical alpine regions. Finally, we discuss the link between habitat stability and tropical alpine plant diversity. Our climate-based delimitation approach can be easily applied to other ecosystems using our developed code, facilitating macro-comparative studies of habitat dynamics through time.

The general dynamic model (GDM) of oceanic island biogeography views oceanic islands predominantly as sinks rather than sources of dispersing lineages. To test this, we conducted a biogeographic analysis of a highly successful insular plant taxon, Cyrtandra, and inferred the directionality of dispersal and founder events throughout the four biogeographical units of the Indo-Australian Archipelago (IAA), namely Sunda, Wallacea, Philippines, and Sahul. Sunda was recovered as the major source area, followed by Wallacea, a system of oceanic islands. The relatively high number of events originating from Wallacea is attributed to its central location in the IAA and its complex geological history selecting for increased dispersibility. We also tested if diversification dynamics in Cyrtandra follow predictions of adaptive radiation, which is the dominant process as per the GDM. Diversification dynamics of dispersing lineages of Cyrtandra in the Southeast Asian grade showed early bursts followed by a plateau, which is consistent with adaptive radiation. We did not detect signals of diversity-dependent diversification, and this is attributed to Southeast Asian cyrtandras occupying various niche spaces, evident by their wide morphological range in habit and floral characters. The Pacific clade, which arrived at the immaturity phase of the Pacific Islands, showed diversification dynamics predicted by the island immaturity speciation pulse model (IISP), wherein rates increase exponentially, and their morphological range is controlled by the least action effect favoring woodiness and fleshy fruits. Our study provides a first step toward a framework for investigating diversification dynamics as predicted by the GDM in highly successful insular taxa.

Scots pine (Pinus sylvestris L.) is one of the most important tree species of the temperate and boreal zones in Eurasia. Its wide distribution range and current patterns of genetic variation have been influenced by Quaternary climatic oscillations and the demographic processes connected to them. In order to better understand the relationship between evolutionary history and demographic factors in a widespread species with a large genome, we used the single-nucleotide polymorphism (SNP) array to genotype thousands of SNP markers across 62 natural populations (N = 686 trees) of Scots pine in Eurasia. This provides the largest range-wide SNPs' genetic diversity assessment of Scots pine to date. Our findings show evidence of past admixture events between genetic clusters that were retained despite the potential for effective pollen-mediated gene flow across the species' distribution range. We also examined the contemporary population structure of the species and analyzed the range-wide genetic diversity patterns. Phylogenetic analyses and demographic modeling suggest that the observed divergence patterns between genetic lineages likely predate the last glaciation events. Two of the most distinctive groups are represented by trees from the eastern parts of Fennoscandia and Eastern Russia, which have remained separated since the mid-Pleistocene. The patterns of genetic variation also confirm the dual colonization of Fennoscandia and the existence of an admixture zone in Central Europe that was formed during multiple waves of postglacial recolonization. This study provides insights into the genetic relationships of Scots pine populations from Europe and Asia and offers a more comprehensive understanding of the species' history.

Islands in the tropical Pacific Ocean are renowned for high biodiversity and endemism despite having relatively small landmasses. However, our knowledge of how this biodiversity is formed remains limited. The taxon cycle, where well-dispersed, earlier colonizers become displaced from coastal to inland habitats by new waves of colonizers, producing isolated, range-restricted species, has been proposed to explain current biodiversity patterns. Here, we integrate the outcomes of phylogenetic studies in the region to investigate the sources, age, number of colonizations, and diversification of 16 archipelagos in the tropical and subtropical South Pacific. We then evaluate whether the results support the taxon cycle as a plausible mechanism for these observations. We find that most species in the Pacific arrived less than 5 Mya from geographically close sources, suggesting that colonization by new taxa is a frequent and ongoing process. Therefore, our findings are broadly consistent with the theory of the Taxon Cycle, which posits that ongoing colonization results in the gradual displacement of established lineages. Only the oldest archipelagos, New Caledonia and Fiji, do not conform to this trend, having proportionally less recent colonization events, suggesting that the taxon cycle may slow on older islands. This conclusion is further validated by New Caledonia having lower diversification rate estimates than younger islands. We found that diversification rates across archipelagos are negatively correlated with area and age. Therefore, a taxon cycle that slows with island age appears to be a suitable concept for understanding the dynamic nature and biodiversity patterns of the Pacific Islands.

Phylogenomics enhances our understanding of plant radiations in the biodiverse Andes. Our study focuses on Puya, primarily Andean and a part of the Bromeliaceae family. Using a phylogenomic framework based on the Angiosperms353 probe set for 80 species, we explored Puya′s phenotypic evolution and biogeography. Divergence time analyses and ancestral area estimations suggested that Puya originated in Central Coastal Chile around 9 million years ago (Ma). Subsequently, it dispersed to the dry valleys of the Central Andes and Puna regions between 5–8 Ma, leading to the emergence of major lineages. Key events in the last 2–4 million years include the recolonization of Chilean lowlands and dispersal to the northern Andes via Peru's Jalcas, facilitating passage through the Huancabamba depression. This event gave rise to the high-elevation Northern Andes clade. Using phylogenetic comparative methods, we tested the hypothesis that adaptation to the Andes' island-like high-elevation ecosystems was facilitated by unique leaf and floral traits, life history, and inflorescence morphology. Our findings suggest correlations between inflorescence axis compression, protective bract overlap, and high-elevation living, potentially preventing reproductive structure freezing. Semelparity evolved exclusively at high elevations, although its precise adaptive value remains uncertain. Our framework offers insights into Andean evolution, highlighting that lineages adapted to life in dry ecosystems can easily transition to high-elevation biomes. It also underscores how the island-like nature of high-elevation ecosystems influences phenotypic evolution rates. Moreover, it opens avenues to explore genetic mechanisms underlying adaptation to extreme mountain conditions.