Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals

Phylogenomic analysis reveals the evolutionary history of Paleartic needle-leaved junipers

Needle-leaved junipers (Juniperus sect. Juniperus, Cupressaceae) are coniferous trees and shrubs with red or blue fleshy cones. They are distributed across Asia, Macaronesia and the Mediterranean Basin, with one species (J. communis) having a circumboreal distribution. Here we aim to resolve the phylogeny of this clade to infer its intricate evolutionary history. To do so, we built a comprehensive, time-calibrated phylogeny using genotyping-by-sequencing (GBS) and combine it with species occurrence using phylogeographic tools. Our results provide solid phylogenetic resolution to propose a new taxonomic classification and a biogeographical history of the section. Specifically, we confirm the monophyly of two groups within J. sect. Juniperus: the Asian (blue-cone) species including the circumboreal J. communis, and the Mediterranean-Macaronesian (red-cone) species. In addition, we provide strong phylogenetic evidence for three distinct species (J. badia, J. conferta, J. lutchuensis) previously considered subspecies or varieties, as well as for the differentiation between the eastern and western Mediterranean lineages of J. macrocarpa. Our findings suggest that the Mediterranean basin was the primary center of diversification for Juniperus sect. Juniperus, followed by an East Asian-Tethyan disjunction resulting from uplifts of the Qinghai-Tibetan Plateau and climatic shifts. The colonization history of Macaronesia by red-cone junipers from the western Mediterranean appears to have taken place independently in two different geological periods: the Miocene (Azores) and the Pliocene (Madeira-Canary Islands). Overall, genomic data and phylogenetic analysis are key to consider a new taxonomic proposal and reconstruct the biogeographical history of the iconic needle-leaved junipers across the Paleartic.

Key innovations and niche variation promoted rapid diversification of the widespread Juniperus (Cupressaceae)

The processes of forming lineages undergoing widespread radiations remain a knowledge gap that is fundamental to our understanding of the geographic distributions of species. Although early studies emphasized the importance of dispersal ability and historical migration events, key innovations that promote rapid diversification and/or adaptation to new habitats may also strongly influence distribution ranges. Juniperus is the second largest genus of conifers and is widely distributed throughout the Northern Hemisphere. Here, we used phylogenetic, phenotypic, and climatic data to investigate the contributions of these processes to the wide distribution and rapid diversification of Juniperus. Combining a time-scaled phylogeny and macroevolutionary theory, we show that the key innovations of berry-like seed cones and dioecy promoted the rapid diversification of Juniperus and that increased dispersal ability promoted allopatric speciation. Ecological niches had significant divergence among different clades of Juniperus. Biogeographic results supported multiple long-distance dispersal events and niche variation that contributed to the modern range of Juniperus, while both phenotypic adaptation and ecological opportunity probably drove its distribution range. Our findings suggest that the current widespread distribution is likely the result of significant divergence driven by niche variation in which ecological opportunities from key innovation and phenotypic divergence.

Comparative analysis of chloroplast genomes of seven Juniperus species from Kazakhstan

Juniperus species are shrubs or trees in the family Cupressaceae that play an important role in forest ecosystems. In this study, we report the complete sequences of the plastid (pt) genomes of five Juniperus species collected in Kazakhstan (J. communis, J. sibirica, J. pseudosabina, J. semiglobosa, and J. davurica). The sequences of the pt genomes of the five species were annotated in addition to two full pt genome sequences from J. sabina and J. seravschanica, which we have previously reported. The pt genome sequences of these seven species were compared to the pt genomes of Juniperus species available in the public NCBI database. The total length of the pt genomes of Juniperus species, including previously published pt genome data, ranged from 127,469 bp (J. semiglobosa) to 128,097 bp (J. communis). Each Juniperus plastome consisted of 119 genes, including 82 protein-coding genes, 33 transfer RNA and 4 ribosomal RNA genes. Among the identified genes, 16 contained one or two introns, and 2 tRNA genes were duplicated. A comparative assessment of pt genome sequences suggested the identification of 1145 simple sequence repeat markers. A phylogenetic tree of 26 Juniperus species based on the 82 protein-coding genes separated the Juniperus samples into two major clades, corresponding to the Juniperus and Sabina sections. The analysis of pt genome sequences indicated that accD and ycf2 were the two most polymorphic genes. The phylogenetic evaluation of 26 Juniperus species using these two genes confirmed that they can be efficiently used as DNA barcodes for phylogenetic analyses in the genus. The sequenced plastomes of these Juniperus species have provided a large amount of genetic data that will be valuable for future genomic studies of this genus.

Past climate cooling and orogenesis of the Hengduan Mountains have influenced the evolution of Impatiens sect. Impatiens (Balsaminaceae) in the Northern Hemisphere

BACKGROUND

Impatiens sect. Impatiens is distributed across the Northern Hemisphere and has diversified considerably, particularly within the Hengduan Mountains (HDM) in southwest China. Yet, the infra-sectional phylogenetic relationships are not well resolved, largely due to limited taxon sampling and an insufficient number of molecular markers. The evolutionary history of its diversification is also poorly understood. In this study, plastome data and the most complete sampling to date were used to reconstruct a robust phylogenetic framework for this section. The phylogeny was then used to investigate its biogeographical history and diversification patterns, specifically with the aim of understanding the role played by the HDM and past climatic changes in its diversification.

RESULTS

A stable phylogeny was reconstructed that strongly supported both the monophyly of the section and its division into seven major clades (Clades I-VII). Molecular dating and ancestral area reconstruction suggest that sect. Impatiens originated in the HDM and Southeast China around 11.76 Ma, after which different lineages dispersed to Northwest China, temperate Eurasia, and North America, mainly during the Pliocene and Pleistocene. An intercontinental dispersal event from East Asia to western North America may have occurred via the Bering Land Bridge or Aleutian Islands. The diversification rate was high during its early history, especially with the HDM, but gradually decreased over time both within and outside the HDM. Multiple linear regression analysis showed that the distribution pattern of species richness was strongly associated with elevation range, elevation, and mean annual temperature. Finally, ancestral niche analysis indicated that sect. Impatiens originated in a relatively cool, middle-elevation area.

CONCLUSIONS

We inferred the evolutionary history of sect. Impatiens based on a solid phylogenetic framework. The HDM was the primary source or pump of its diversity in the Northern Hemisphere. Orogeny and climate change may have also shaped its diversification rates, as a steady decrease in the diversification rate coincided with the uplift of the HDM and climate cooling. These findings provide insights into the distribution pattern of sect. Impatiens and other plants in the Northern Hemisphere.

Biogeographic diversification of Actaea (Ranunculaceae): Insights into the historical assembly of deciduous broad-leaved forests in the Northern Hemisphere

The deciduous broad-leaved forests (DBLFs) cover large temperate and subtropical high-altitude regions in the Northern Hemisphere. They are home to rich biodiversity, especially to numerous endemic and relict species. However, we know little about how this vegetation in the Northern Hemisphere has developed through time. Here, we used Actaea (Ranunculaceae), an herbaceous genus almost exclusively growing in the understory of the Northern Hemisphere DBLFs, to shed light on the historical assembly of this biome in the Northern Hemisphere. We present a complete species-level phylogenetic analysis of Actaea based on five plastid and nuclear loci. Using the phylogenetic framework, we estimated divergence times, ancestral ranges, and diversification rates. Phylogenetic analyses strongly support Actaea as monophyletic. Sections Podocarpae and Oligocarpae compose a clade, sister to all other Actaea. The sister relationship between sections Chloranthae and Souliea is strongly supported. Section Dichanthera is not monophyletic unless section Cimicifuga is included. Actaea originated in East Asia, likely the Qinghai-Tibet Plateau, in the late Paleocene (c. 57 Ma), and subsequently dispersed into North America in the middle Eocene (c. 43 Ma) via the Thulean bridge. Actaea reached Europe twice, Japan twice, and Taiwan once, and all these five colonization events occurred in the late Miocene-early Pliocene, a period when sea level dropped. Actaea began to diversify at c. 43 Ma. The section-level diversification took place at c. 27-37 Ma and the species-level diversification experienced accelerations twice, which occurred at c. 15 Ma and c. 5 Ma, respectively. Our findings suggest that the Northern Hemisphere DBLFs might have risen in the middle Eocene and further diversified in the late Eocene-Oligocene, middle Miocene and early Pliocene, in association with climatic deterioration during these four periods.

Morphological and phylogenetic features of the Crimean population of Juniperus deltoides R.P. Adams

Juniperus deltoides is a relict species from the Tertiary Period. It is a typical representative of the Mediterranean group of the section Juniperus. It is included in the Red Books of the Republic of Crimea and the city of Sevastopol. Until recently, it was believed that a population of J. oxycedrus grew in Crimea. Currently, J. deltoides is described as a cryptic species, morphologically difficult to distinguish from J. oxycedrus. As a result, it became necessary to conduct a series of detailed studies to determine the morphological and phylogenetic features of the Crimean cryptic population in order to identify it as being one of the species of the cryptic pair. The studies were carried out in two stages: at the first stage, the morphological features of the vegetative and generative organs and their difference from J. oxycedrus were determined; the second stage included genetic research. The length of the needles of the Crimean population is 12.94 ± 0.19 mm, which corresponds to the Eastern Italian population of J. deltoides. At the same time, the width of the needles is 1.39 ± 0.02 mm, which is typical of the Portuguese population of J. oxycedrus. The dimensions of the cones are d1 (conditional height) = 7.54 ± 0.14 mm, and d2 (conditional width) = 9.11 ± 0.09 mm, which is more in line with J. deltoides. The shapes of the cones are very diverse. Some individuals have cones, the covering scales of which are visually indistinguishable, and their tops are completely fused. A similar phenomenon is characteristic of the Western Mediterranean populations of J. oxycedrus. Morphological analysis of the vegetative and generative organs of J. deltoides showed that when these two traits are combined, it is not possible to reliably distinguish between J. deltoides and J. oxycedrus individuals. Nuclear (ITS internal transcribed spacer) and chloroplast (petN-psbM, trnS-trnG) non-coding regions of the genome were used for genetic analysis. Studies have shown that the nuclear regions of genes have greater variability than chloroplast regions. The sequences obtained in this work formed a clade with J. deltoides samples 9430 and 9431 (BAYLU) growing in Turkey, which makes it possible to assign the samples studied to J. deltoides.

Time-calibrated phylogenies reveal mediterranean and pre-mediterranean origin of the thermophilous vegetation of the Canary Islands

BACKGROUND AND AIMS

The Canary Islands have strong floristic affinities with the Mediterranean Basin. One of the most characteristic and diverse vegetation belts of the archipelago is the thermophilous woodland (between 200 and 900 m.a.s.l.). This thermophilous plant community consists of many non-endemic species shared with the Mediterranean Floristic Region together with Canarian endemic species. Consequently, phytogeographic studies have historically proposed the hypothesis of an origin of the Canarian thermophilous species following the establishment of the summer-dry mediterranean climate in the Mediterranean Basin around 2.8 million years ago.

METHODS

Time-calibrated phylogenies for 39 plant groups including Canarian thermophilous species were primarily analysed to infer colonization times. In particular, we used 26 previously published phylogenies together with 13 new time-calibrated phylogenies (including newly generated plastid and nuclear DNA sequence data) to assess whether the time interval between stem and crown ages of Canarian thermophilous lineages postdates 2.8 Ma. For lineages postdating this time threshold, we additionally conducted ancestral area reconstructions to infer the potential source area for colonization.

KEY RESULTS

A total of 43 Canarian thermophilous lineages were identified from 39 plant groups. Both mediterranean (16) and pre-mediterranean (9) plant lineages were found. However, we failed to determine the temporal origin for 18 lineages because a stem-crown time interval overlaps with the 2.8-Ma threshold. The spatial origin of thermophilous lineages was also heterogeneous, including ancestral areas from the Mediterranean Basin (nine) and other regions (six).

CONCLUSIONS

Our findings reveal an unexpectedly heterogeneous origin of the Canarian thermophilous species in terms of colonization times and mainland source areas. A substantial proportion of the lineages arrived in the Canaries before the summer-dry climate was established in the Mediterranean Basin. The complex temporal and geographic origin of Canarian thermophilous species challenges the view of the Canary Islands (and Madeira) as a subregion within the Mediterranean Floristic Region.

Historical biogeography and diversification of ringless Amanita (section Vaginatae) support an African origin and suggest niche conservatism in the Americas

Ectomycorrhizal fungi (ECM) sustain nutrient recycling in most terrestrial ecosystems, yet we know little about what major biogeographical events gave rise to present-day diversity and distribution patterns. Given the strict relationship between some ECM lineages and their hosts, geographically well-sampled phylogenies are central to understanding major evolutionary processes of fungal biodiversity patterns. Here, we focus on Amanita sect. Vaginatae to address global diversity and distribution patterns. Ancestral-state-reconstruction based on a 4-gene timetree with over 200 species supports an African origin between the late Paleocene and the early Eocene (ca. 56 Ma). Major biogeographic "out-of-Africa" events include multiple dispersal events to Southeast Asia (ca. 45-21 Ma), Madagascar (ca. 18 Ma), and the current Amazonian basin (ca. 45-36 Ma), the last two likely trans-oceanic. Later events originating in Southeast Asia involve Nearctic dispersal to North America (ca. 20-5 Ma), Oceania (Australia and New Zealand; ca. 15 Ma), and Europe (ca. 10-5 Ma). Subsequent dispersals were also inferred from Southeast Asia to East Asia (ca. 4 Ma); from North America to East Asia (ca. 11-8 Ma), Southeast Asia (ca. 19-2 Ma), Northern Andes (ca. 15 Ma), and Europe (ca. 15-2 Ma), respectively; and from the Amazon to the Caribbean region (ca. 25-20 Ma). Finally, we detected a significant increase in the net diversification rates in the branch leading to most northern temperate species in addition to higher state-dependent diversification rates in temperate lineages, consistent with previous findings. These results suggest that species of sect. Vaginatae likely have higher dispersal ability and higher adaptability to new environments, in particular compared to those of its sister clade, sect. Caesareae. Overall, the much wider distribution of A. sect. Vaginatae, from pan-tropical to pan-arctic, provides a unique window to understanding niche conservatism across a species-rich clade of ECM fungi.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae)

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

Origin and evolutionary history of Populus (Salicaceae): Further insights based on time divergence and biogeographic analysis

INTRODUCTION

Populus (Salicaceae) species harbour rich biodiversity and are widely distributed throughout the Northern Hemisphere. However, the origin and biogeography of Populus remain poorly understood.

METHODS

We infer the divergence times and the historical biogeography of the genus Populus through phylogenetic analysis of 34 chloroplast fragments based on a large sample.

RESULTS AND DISCUSSION

Eurasia is the likely location of the early divergences of Salicaceae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by recurrent spread to the remainder of the Old World and the New World beginning in the Eocene; the extant Populus species began to diversity during the early Oligocene (approximately 27.24 Ma), climate changes during the Oligocene may have facilitated the diversification of modern poplar species; three separate lineages of Populus from Eurasia colonized North America in the Cenozoic via the Bering Land Bridges (BLB); We hypothesize that the present day disjunction in Populus can be explained by two scenarios: (i) Populus likely originated in Eurasia and subsequently colonized other regions, including North America; and (ii) the fact that the ancestor of the genus Populus that was once widely distributed in the Northern Hemisphere and eventually wiped out due to the higher extinction rates in North America, similar to the African Rand flora. We hypothesize that disparities in extinction across the evolutionary history of Populus in different regions shape the modern biogeography of Populus. Further studies with dense sampling and more evidence are required to test these hypotheses. Our research underscores the significance of combining phylogenetic analyses with biogeographic interpretations to enhance our knowledge of the origin, divergence, and distribution of biodiversity in temperate plant floras.

Miocene Cupressinoxylon from Gökçeada (Imbros), Turkey with Protophytobia cambium mining and the study of ecological signals of wood anatomy

Premise

The recognition of the Miocene Climate Optimum (MCO) in terrestrial palaeoenvironments of the Eastern Mediterranean is restricted to Lesbos and Lemnos Islands, Greece. This area is significant for its wood microfossils. A recently-discovered fossil wood assemblage from Gökçeada (Imbros) Island, Turkey, including tree species similar to the Greek findings, is thought to have an early Miocene age. Here, we revise the age of the latter plant fossiliferous locality, re-evaluate the area for the study of MCO for the terrestrial palaeoecosystems of the Eastern Mediterranean and the nomenclature errors referring to the occurrence of fossil wood. We present the plant-insect-environment interactions using detailed anatomical descriptions, of an extinct conifer and its extinct cambium miner feeding traces observed in its secondary xylem.

Methods

Three thin sections were prepared with standard palaeoxylotomical techniques from a small section of the silicified wood; the sections were observed under a light microscope. The anatomy of the conifer and its damage patterns were compared with those of extant and fossil Cupressaceae and Agromyzidae, respectively.

Pivotal results

The common anatomical features of the studied wood specimen and Hesperocyparis macrocarpa (Hartw.) Bartel and a shared characteristic (the number of the cross-field pits - a feature we consider of diagnostic value) with Xanthocyparis vietnamensis Farjon & T.H. Nguyên led to its assignment to the Hesperocyparis-Xanthocyparis-Callitropsis clade. The detailed study of the wound scars and anatomical abnormalities, the anatomical-environmental associations, and structural-functional reactions follow the identification of the wood's anatomy sensu Carlquist providing decisive results.

Conclusions

Based on the distinctive characteristics presented, we identify our macrofossil as Cupressinoxylon matromnense Grambast, a stem or an extinct lineage of the Hesperocyparis-Xanthocyparis vietnamensis-Callitropsis nootkatensis clade with feeding traces of the fossil cambium miner of the genus Protophytobia Süss (Diptera: Agromyzidae), and anatomical damage and reaction tissue on adventitious shoots. The use of Protopinaceae and Pinoxylon F. H. Knowlton from the eastern Mediterranean are re-evaluated and corrections are provided. The age of the studied plant fossiliferous locality in Gökçeada is revised as middle Miocene, allowing the proposal of an eastern Mediterranean MCO hotspot, including Lesbos, Lemnos, and Gökçeada (Imbros) Islands.

The stepwise Indian-Eurasian collision and uplift of the Himalayan-Tibetan plateau drove the diversification of high-elevation Scytodes spiders

The Cenozoic Indian-Eurasian collision and uplift of the Himalayan-Tibetan Plateau (HTP) are among the most important geological events in the world. They have affected the diversification of regional biota of many taxonomic groups on and around the HTP. However, the exact timing and model of the collision and uplift events and speciation on and around the HTP are still in debate. The Himalayas group of Scytodes spitting spiders (Araneae: Scytodidae) are distributed at high elevations of the HTP and northern Indochina. Here, we reconstruct a dated molecular phylogeny for pan-Himalayan Scytodes spiders, including the Himalayas group, with full geographical sampling of the species from the HTP and Indochina. We test a hypothesis to explain that the rich montane biodiversity of the region is uplift-driven diversification-that orogeny drives rapid in situ speciation of the resident Scytodes lineages. Our findings revealed that the separation of the Himalayas clade from the Myanmar clade took place during the middle Oligocene, reflecting the final collision of India with Eurasia. The deep divergences among three clades (the Himalayas, the Myanmar and the Indochina clades) occurred from the middle Eocene to the middle Oligocene, corresponding to two early uplift events of the HTP. The evolutionary split between the Himalayas + Myanmar and Indochina clades were simultaneous with the rapid lateral extrusion of Indochina by the initial Himalayan uplift around the Eocene. This study highlights the importance of the diversification of dispersal-limited, high-elevation invertebrates as independent lines of evidence to reflect key tectonic events in the Himalayan-Indochina region, supporting the stepwise model for the Indian-Eurasian collision and uplift of HTP.

Molecular Identification of African Nymphaea Species (Water Lily) Based on ITS, trnT-trnF and rpl16

The genus Nymphaea L. (water lily) is the most diverse genus in the family Nymphaeaceae, with more than 50 species worldwide, including 11 species distributed in Africa. The complex and variable morphology of Nymphaea makes it extremely difficult to accurately identify species based on morphological characteristics alone. DNA barcoding has the potential to identify species accurately. In this study, 158 Nymphaea populations from seven African countries were collected for species identification by ITS, trnT-trnF and rpl16. Additionally, the three candidate DNA barcodes were evaluated for genetic distance and barcoding gap. Based on the comprehensive analysis of sequence similarity, genetic distance method and phylogenetic tree, a total of 137 populations of seven Nymphaea species from African were well-identified, including N. lotus, N. petersiana, N. zenkeri, N. nouchali var. caerulea, N. micrantha and N. guineensis. ITS has more obvious advantages over trnT-trnF, rpl16 and trnT-trnF+rpl16 in the intraspecific and interspecific variation differences and barcoding gap and can identify most species. trnT-trnF and rpl16 can identify some species that cannot be identified by ITS. The results showed that it is more appropriate to apply the combination of ITS and trnT-trnF (or rpl16) as the DNA barcoding of Nymphaea. Additionally, this study further enriches the DNA barcoding database of Nymphaea and provides a reference basis for studying taxonomy, phylogenetics and evolutionary origin of Nymphaea.

Mummified Wood of Juniperus (Cupressaceae) from the Late Miocene of Taman Peninsula, South Russia

Juniperus L. is the second-largest genus of conifers, having the widest distribution of all conifer genera. Its phytogeographic history is, however, obscure due to its very poor fossil record. We described a wood of Juniperus sp. from the lower Maeotian sediments of the Popov Kamen section, Taman Peninsula, South Russia, in order to clarify its taxonomic position shedding light on the phytogeographic history of the genus. This fossil wood was well-preserved by mummification, which allowed for it to be studied by the same methods as used for the anatomical examination of modern woods. The wood from the Popov Kamen section shows the greatest similarity to the extant Mediterranean species J. excelsa, belonging to the section Sabina. This is the first reliable macrofossil evidence of the sect. Sabina from Eurasia convincingly dated to the Miocene. The age of the mummified wood from the Popov Kamen section is consistent with molecular dating of diversification of the lineage comprising juniper species of the sect. Sabina from Europe, Asia and eastern Africa. The wood of Juniperus sp. has not been buried in situ, as it was found in the relatively deep-water marine sediments. The available coeval pollen series and macrofossils of Cupressaceae from the surrounding regions suggest that this wood was likely transferred by sea current from the northwestern side of the Black Sea, which was a part of the Eastern Paratethys.

Incipient speciation, high genetic diversity, and ecological divergence in the alligator bark juniper suggest complex demographic changes during the Pleistocene

The most recent glacial cycles of the Pleistocene affected the distribution, population sizes, and levels of genetic structure of temperate-forest species in the main Mexican mountain systems. Our objective was to investigate the effects these cycles had on the genetic structure and distribution of a dominant species of the "mexical" vegetation across North and Central America. We studied the genetic diversity of Juniperus deppeana, a conifer distributed from the Southwestern United States to the highlands of Central America. We combined information of one plastid marker and two nuclear markers to infer phylogeographic structure, genetic diversity and demographic changes. We also characterized the climatic niche for each variety to infer the plausible area of suitability during past climatic conditions and to evaluate climatic niche discontinuities along with the species distribution. We found a marked phylogeographic structure separating the populations North and South of the Isthmus of Tehuantepec, with populations to the South of this barrier forming a distinct genetic cluster corresponding to Juniperus deppeana var. gamboana. We also found signals of population expansion in the Northern genetic cluster. Ecological niche modeling results confirmed climatic niche differences and discontinuities among J. deppeana varieties and heterogeneous responses to climatic oscillations. Overall, J. deppeana's genetic diversity has been marked by distribution shifts, population growth and secondary contact the North, and in situ permanence in the South since the last interglacial to the present. High genetic variation suggests a wide and climatically diverse distribution during climatic oscillations. We detected the existence of two main genetic clusters, supporting previous proposals that Juniperus deppeana and Juniperus gamboana may be considered two separate species.

Nomenclatural notes of Sabinaconvallium var. microsperma (Cupressaceae)

The name Sabinaconvalliumvar.microsperma W.C.Cheng & W.T.Wang was not validly published when it was first described in 1975, but was validated in 1978 at the same time as the specific combination Sabinaconvallium (Rehder & E.H.Wilson) W.C.Cheng & W.T.Wang was validly published in Flora Reipublicae Popularis Sinicae. Under Art. 41.6 of the Shenzhen Code, other names based onSabinaconvalliumvar.microsperma were valid, including Sabinamicrosperma (W.C.Cheng & L.K.Fu) W.C.Cheng & L.K.Fu, Juniperusconvalliumvar.microsperma (W.C.Cheng & L.K.Fu) Silba, Juniperusmicrosperma (W.C.Cheng & L.K.Fu) R.P.Adams, despite the reference of the basionym being erroneously cited when these authors made the new combinations.

Recent advances on phylogenomics of gymnosperms and a new classification

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.

Phylogeographical Pattern and Population Evolution History of Indigenous Elymus sibiricus L. on Qinghai-Tibetan Plateau

Elymus sibiricus L. is a perennial allotetraploid belonging to Triticeae of Poaceae, Elymus L., as the type species of genus Elymus L. The existing geographical distribution pattern and genetic spatial structure of E. sibiricus on Qinghai-Tibetan Plateau (QTP) are not yet clear. In this study, population genetic structure and demography history of 216 individuals from 44 E. sibiricus populations on QTP were studied used specific-locus amplified fragment sequencing (SLAF-seq). The result of genetic diversity showed that there was no single genetic diversity center was observed across all E. sibiricus populations. The results of genetic variation showed that 44 populations were clearly divided into the following three groups: Qinghai Plateau (Group I), South Tibet (Group II), and Hengduan Mountains (Group III). From the three analyses of AMOVA, Mantel test and Treemix, strong genetic differentiation across all populations and low genetic differentiation among populations within three groups. Molecular dating indicated that E. sibiricus diverged at 16.08 Ma (during the early Miocene) can be linked to the Himalayan Motion stage of QTP uplift. It is speculated that the reasons affecting the current phylogeographical pattern are as follows: (1) The environmental changes due to the uplift of the QTP; (2) The geographic distance between the populations (Groups I and III are close in geographic distance, and gene flow are frequent); (3) Geographical barriers (the Tanggula and Bayangela Mountains between Groups I and II). This study provides new evidence and historical perspective to the future exploration of the evolution and geographic distribution pattern of Elymus L.

Incomplete lineage sorting and local extinction shaped the complex evolutionary history of the Paleogene relict conifer genus, Chamaecyparis (Cupressaceae)

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the "boreotropical" flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Essential Oils and Extracts of Juniperus macrocarpa Sm. and Juniperus oxycedrus L.: Comparative Phytochemical Composition and Anti-Proliferative and Antioxidant Activities

In this work, we conducted a comparative phytochemical, chemotaxonomic, and biological study of essential oils (EOs) and extracts (ethyl acetate and methanol) obtained from the leaves of Juniperusmacrocarpa and J. oxycedrus. The dominant compounds of J. macrocarpa EO, analysed by gas chromatography (GC) and gas chromatography-mass spectrometry (GC-MS), are α-pinene, sabinene, manoyl oxide, and germacrene D, whereas α-pinene, limonene, (Z,E)-farnesol, β-pinene, and γ-cadinene are the most representative volatiles of J. oxycedrus EOs. A multivariate analysis of EOs, included a selection of literature data comparing our samples to samples of J. oxycedrus/macrocarpa/deltoides from the Mediterranean area, was performed. As evident by high-performance liquid chromatography (HPLC) analyses, apigenin, (−)-epicatechin, and luteolin were abundant in J. oxycedrus extracts, while gallic acid, kaempferol-3-O-glucoside, and protocatechuic acid were the dominant constituents of J. macrocarpa extracts. EOs and extracts have been investigated for their potential antioxidant properties and anti-proliferative activity against lung adenocarcinoma (A549), breast cancer (MCF-7 and MDA-MB-231), and lung large cell carcinoma (COR-L23) human cell lines. The methanol and ethyl acetate extracts of J. oxycedrus exerted the most valuable antioxidant activity and exhibited the most promising activity against the COR-L23 cell line with an IC₅₀ of 26.0 and 39.1 μg/mL, respectively, lower than that obtained with the positive control (IC₅₀ of 45.5 μg/mL). To the best of our knowledge, this is the first report highlighting the anti-proliferative activity of J. oxycedrus and J. macrocarpa extracts against this lung cancer cell line. Our results indicate that J. oxycedrus may be considered a source of natural compounds with antioxidant and anti-proliferative effects that could be suitable for future applications.

Evaluation of the Antimicrobial and Antivirulent Potential of Essential Oils Isolated from Juniperus oxycedrus L. ssp. macrocarpa Aerial Parts

As a consequence of the worsening situation with multidrug-resistant (MDR) pathogens and a disparity in the commercialization of novel antimicrobial agents, scientists have been prompted to seek out new compounds with antimicrobial activity from a wide range of sources, including medicinal plants. In the present study, the antibacterial, antifungal, anti-virulence, and resistance-modulating properties of the essential oil from the Sardinian endemic Juniperus oxycedrus L. ssp. macrocarpa aerial parts were evaluated. The GC/MS analysis showed that the main compounds in the oil were α-pinene (56.63 ± 0.24%), limonene (14.66 ± 0.11%), and β-pinene (13.42 ± 0.09%). The essential oil showed potent antibacterial activity against Gram-positive bacteria (0.25-2 v/v%) and Salmonella spp. (4 v/v%). The strongest fungicidal activity was recorded against Candida auris sessile cells (median FICI was 0.088) but not against C. albicans biofilms (median FICI was 1). The oil showed potent efflux pump inhibitory properties in the case of Staphylococcus aureus and Escherichia coli. The therapeutic potential of Juniperus may be promising for future more extensive research and in vivo tests to develop new drugs against antibiotic and antifungal resistance.

Productivity and Seasonality Drive Total Soil Respiration in Semi-Arid Juniper Woodlands (Juniperus thurifera L., Southern Spain)

We analyzed the relationship between forest productivity (joint effect of forest maturity and soil quality) and total soil respiration (µmol CO2 m−2 s−1) in semi-arid juniper woodlands (young woodlands growing in abandoned farmlands with deeper soils, and mature woodlands in lithic soils), and investigated the seasonal variation in soil CO2 efflux as a function of soil temperature and the soil water content. We measured the soil CO2 efflux from twelve cylinders in the soil over a three-year period using LI-6400 equipment. The results show that, in the more productive site (young woodland), soil CO2 efflux was higher due to greater respiration, mainly in the driest periods. Soil respiration followed a seasonal trend, being higher in spring and decreasing in cold periods. In both juniper woodlands and especially in the older forest, the CO2 efflux rates were low (<2.5 for Q10), typical of slow-growing species. Soil respiration was controlled by soil temperature without drought and in the temperate-warm season, whereas respiration showed sensitivity to soil water content in periods when edaphic humidity was low (but only in the more productive, young forest, which seemed to show better adaptation to drought), and under high soil moisture (soil water > 25%) for both woodlands, coinciding with warm temperatures in the spring. This period also corresponded to the highest CO2 efflux recorded in both woodlands. The accumulation of organic C seems to also be important to maintain elevated soil respiration in summer, especially in young woodlands. Thus, apart from microclimatic conditions, factors related to productivity regulate respiratory activity.

Disjunction and Vicariance Between East and West Asia: A Case Study on Euonymus sect. Uniloculares Based on Plastid Genome Analysis

Scientists have long been captivated by biogeographic disjunctions, and disjunctions between East Asia and North America have been particularly well-studied at the genus and family levels. By contrast, disjunctions between eastern and western Asia have received less attention. Euonymus L. is taxonomically divided into two sections based on the number of cells in anthers as follows: E. sect. Uniloculares has one-celled anthers and occurs mainly in Asia, whereas E. sect. Biloculares has two-celled anthers and is distributed globally. We used Illumina sequencing to investigate the genomes of four species in sect. Uniloculares. The chloroplast (cp) genomes are highly conserved (157,290-158,094 bp). Pseudogenisation of ndhF and intron loss in rps16 was detected. Based on the cp genomes of the four species of E. sect. Uniloculares, we propose a novel hypothesis of disjunction between eastern and western Asia. Biogeographic reconstruction and molecular dating revealed that sect. Uniloculares separated from its sect. Biloculares forebears 4.0 Mya during the Pliocene era. The radial diversification of sect. Uniloculares from East Asia and the establishment of the western Asian clade during the Pleistocene era (1.9 Mya) were the results of both dispersal and vicariance, making the section the youngest diverged clade conforming to age estimation. The centre of origin of sect. Uniloculares was determined to be in East Asia. Disjunctions and diversification between eastern and western Asia in sect. Uniloculares are thought to have been caused by changes in monsoon patterns, temperature variations, and the emergence of the Gobi Desert.

Phenotypic Comparison of Three Populations of Juniperus turbinata Guss. in North-Eastern Morocco

Juniperus turbinata Guss. is a native species of Morocco; however, an exhaustive taxonomic description based on phenotypical characterization of north-eastern Moroccan population species is lacking, which might lead to taxonomic confusion. In order to expound the phenotypic description of J. turbinata of the north-eastern Moroccan population and to examine the taxonomic differences within it; a comparative analysis of cones, leaves, and seeds was performed between three populations. A total of 280 samples were compared on the basis of nine measured and eight calculated traits. The results reveal significant interpopulation changes in the studied characteristics of cones, leaves, and seeds. The most discriminating traits were associated with the proportion between cone diameter and number of seeds. We detected the lowest number of seeds in coastal population when compared to other localities, but at the same time, the seeds from the littoral were the longest and the widest. In addition, the semi-continental population had the highest quantity of seeds, and leaves had intermediate values for the majority of the assessed traits. The phenotypical difference between populations demonstrates a certain adaptability of the species in a biogeographical pattern. This study is a contribution to completing the description of patterns of phenotypical differences of the Phoenician juniper in the Mediterranean region, and confirms its evolutionary plasticity linked to adaptation to local environmental conditions.

Species Diversification of the Coniferous Pathogenic Fungal Genus Coniferiporia (Hymenochaetales, Basidiomycota) in Association with Its Biogeography and Host Plants

Coniferiporia, belonging to Hymenochaetaceae and now segregated from Phellinidium, is a wood-inhabiting fungal genus with three species, each having a specific geographic distribution and a strong host specificity as a forest pathogen of coniferous trees. In this study, the species diversity of Coniferiporia is further clarified with the aid of a wider sampling and multilocus-based phylogenetic analysis, which reveals a new species Coniferiporia uzbekistanensis. The molecular clock and ancestral geographic origin analyses indicate that the ancestor of Coniferiporia emerged in one of the Pinaceae and Cupressaceae, then jumped to the other plant family originated in eastern Eurasia 17.01 million years ago (Mya; 95% highest posterior density: 9.46 to 25.86 Mya), and later extended its distribution to western North America, Central Asia, and eastern Europe. Coniferiporia sulphurascens speciated on Pinaceae in eastern Eurasia 8.78 Mya (9.46 to 25.86 Mya) and then extended its distribution to western North America and eastern Europe. Coniferiporia qilianensis and C. uzbekistanensis speciated on Juniperus przewalskii in eastern Eurasia 3.67 Mya (0.36 to 8.02 Mya) and on Juniperus polycarpos in Central Asia 4.35 Mya (0.94 to 8.37 Mya), respectively. The speciation event of Coniferiporia weirii occurred 4.45 Mya (0.77 to 9.33 Mya) right after the emergence of its host, the endemic Cupressaceae species Thuja plicata, and soon after, this fungus evolved to also inhabit another endemic Cupressaceae species Calocedrus decurrens. In summary, this study for the first time unambiguously clarified and timed the adaptive evolutionary event of Coniferiporia in association with its biogeography and host plants.

The number of tree species on Earth

One of the most fundamental questions in ecology is how many species inhabit the Earth. However, due to massive logistical and financial challenges and taxonomic difficulties connected to the species concept definition, the global numbers of species, including those of important and well-studied life forms such as trees, still remain largely unknown. Here, based on global ground-sourced data, we estimate the total tree species richness at global, continental, and biome levels. Our results indicate that there are ∼73,000 tree species globally, among which ∼9,000 tree species are yet to be discovered. Roughly 40% of undiscovered tree species are in South America. Moreover, almost one-third of all tree species to be discovered may be rare, with very low populations and limited spatial distribution (likely in remote tropical lowlands and mountains). These findings highlight the vulnerability of global forest biodiversity to anthropogenic changes in land use and climate, which disproportionately threaten rare species and thus, global tree richness.

Innate Immunomodulatory Activity of Cedrol, a Component of Essential Oils Isolated from Juniperus Species

Little is known about the immunomodulatory activity of essential oils isolated from Juniperus species. Thus, we isolated essential oils from the cones and leaves of eight juniper species found in Montana and in Kazakhstan, including J. horizontalis, J. scopolorum, J. communis, J. seravschanica, J. sabina, J. pseudosabina, J. pseudosabina subsp. turkestanica, and J. sibirica. We report here the chemical composition and innate immunomodulatory activity of these essential oils. Compositional analysis of the 16 samples of Juniper essential oils revealed similarities and differences between our analyses and those previously reported for essential oils from this species. Our studies represent the first analysis of essential oils isolated from the cones of four of these Juniper species. Several essential oil samples contained high levels of cedrol, which was fairly unique to three Juniper species from Kazakhstan. We found that these essential oils and pure (+)-cedrol induced intracellular Ca²⁺ mobilization in human neutrophils. Furthermore, pretreatment of human neutrophils and N-formyl peptide receptor 1 and 2 (FPR1 and FPR2) transfected HL60 cells with these essential oils or (+)-cedrol inhibited agonist-induced Ca²⁺ mobilization, suggesting these responses were desensitized by this pretreatment. In support of this conclusion, pretreatment with essential oils from J. seravschanica cones (containing 16.8% cedrol) or pure (+)-cedrol inhibited human neutrophil chemotaxis to N-formyl peptide. Finally, reverse pharmacophore mapping predicted several potential kinase targets for cedrol. Thus, our studies have identified cedrol as a novel neutrophil agonist that can desensitize cells to subsequent stimulation by N-formyl peptide.

Distinct Genetic Structure Reflects Ploidy Level Differentiation in Newly Discovered, Extremely Small Populations of Xanthocyparis vietnamensis from Southwestern China

Population genetic assessment is crucial for the conservation and management of threatened species. Xanthocyparis vietnamensis is an endangered species that is currently restricted to karst mountains in southwestern China and Vietnam. This rare conifer was first recorded in 2002 from northern Vietnam and then in 2013 from Guangxi, China, yet nothing is known about its genetic diversity nor ploidy level variation, although previous cytological study suggest that Vietnamese populations are tetraploids. There have been about 45 individuals found to date in Guangxi, China. Here, we genotyped 33 X. vietnamensis individuals using 20 newly developed, polymorphic microsatellite loci, to assess the genetic variability of its extremely small populations. The genetic diversity of X. vietnamensis (H_E = 0.511) was lower than that of two other heliophile species, Calocedrus macrolepis and Fokienia hodginsii, which have similar distribution ranges. This is consistent with the signature of a genetic bottleneck detected in X. vietnamensis. Although the population genetic differentiation coefficient across loci is moderate (F_ST = 0.125), STRUCTURE analysis revealed two distinct genetic clusters, namely the northern and southern population groups; DAPC analysis grouped the southern populations together in one cluster separate from the northern populations; AMOVA analysis detected a significant genetic differentiation between the two population groups (F_RT = 0.089, p < 0.05), and BARRIER analysis detected a genetic barrier between them. Moreover, we detected differentiation in ploidy level between northern and southern populations, sampled individuals from the former and the later are all diploid and tetraploid cytotypes with mean genome sizes of 26.08 and 48.02 pg/2C, respectively. We deduced that heterogeneous geomorphology and historical events (e.g., human deforestation, Quaternary climate oscillations) may have contributed to population fragmentation and small population size in X. vietnamensis. Considering both genetic and ploidy level differentiation, we propose that two different management units (northern and southern) should be considered and a combination of in situ and ex situ conservation measures should be employed to preserve populations of this endangered species in southwestern China in the light of our findings.

Lectotypification of Chamaecyparishodginsii of the Cupressaceae

Recent phylogenetic studies have suggested that the monotypic Fokienia A.Henry & H.H.Thomas is nested within Chamaecyparis Spach, which is in agreement with separate morphological studies. Here the authors confirm a previous taxonomic treatment that incorporated Fokieniahodginsii (Dunn) A.Henry & H.H.Thomas into Chamaecyparis based on the monophyly requirement of taxonomy, i.e. Chamaecyparishodginsii (Dunn) Rushforth. In addition, the type collection of the basionym Cupressushodginsii Dunn was found to contain three sheets of specimens, one in K including a vegetative branch (K000088294) and a separate ovulate cone (K001090486), a second one in A (A00022477), and a third one in IBSC (IBSC0016081). All three specimens are marked with Hongkong Herbarium No. 3505, but only the two specimens in K and IBSC possess similar handwriting of "Cupressushodginsii Dunn". The two specimens should be considered as syntypes according to the Shenzhen Code. The specimen in K is better preserved but it is a mixture according to the collection label: cones from Foochow (Fuzhou) and foliage from Yenping (Nanping). We lectotypified the name Cupressushodginsii with K000088294 because the specimen is well preserved and has enough characters for identification. Moreover, an ovulate cone (K001090486) is on the same sheet.

Stem CO2 Efflux as an Indicator of Forests’ Productivity in Relict Juniper Woodlands (Juniperus thurifera L.) of Southern Spain

There are considerable uncertainties about the C cycle in semi-arid ecosystems. Hence, studies that have focused on Juniperus in Mediterranean woodlands are non-existent. This study provides a survey of the effect of the juniper woodland type (young and mature woodlands; joint effect of maturity and forest productivity) on stem respiration. We checked the seasonal variation of stem respiration, evaluating the effects of stem temperature on stem CO2 efflux. For this, we measured the stem CO2 efflux (µmol CO2 m−2 s−1) over the four seasons on 16 junipers using LI-6400 equipment. The results showed that in the more productive site (young woodland), the stem CO2 efflux was higher. This variable followed a clear seasonal trend, being higher during the spring and progressively decreasing in cold periods. In both juniper woodlands, and especially in the older forests, the Q10 coefficients were low (<2), typical of cold forests and slow-growing species. The exponential model also confirmed that the Q10 was significantly higher in young juniper trees. Thus, stem CO2 efflux was an indicator of the growth in this juniper woodland that is well adapted to a semi-arid climate.

Phylogenomic and Macroevolutionary Evidence for an Explosive Radiation of a Plant Genus in the Miocene

Mountain systems harbor a substantial fraction of global biodiversity and, thus, provide excellent opportunities to study rapid diversification and to understand the historical processes underlying the assembly of biodiversity hotspots. The rich biodiversity in mountains is widely regarded as having arisen under the influence of geological and climatic processes as well as the complex interactions among them. However, the relative contribution of geology and climate in driving species radiation is seldom explored. Here, we studied the evolutionary radiation of Oreocharis (Gesneriaceae), which has diversified extensively throughout East Asia, especially within the Hengduan Mountains (HDM), using transcriptomic data and a time calibrated phylogeny for 88% (111/126) of all species of the genus. In particular, we applied phylogenetic reconstructions to evaluate the extent of incomplete lineage sorting accompanying the early and rapid radiation in the genus. We then fit macroevolutionary models to explore its spatial and diversification dynamics in Oreocharis and applied explicit birth–death models to investigate the effects of past environmental changes on its diversification. Evidence from 574 orthologous loci suggest that Oreocharis underwent an impressive early burst of speciation starting ca. 12 Ma in the Miocene, followed by a drastic decline in speciation toward the present. Although we found no evidence for a shift in diversification rate across the phylogeny of Oreocharis, we showed a difference in diversification dynamics between the HDM and non-HDM lineages, with higher diversification rates in the HDM. The diversification dynamic of Oreocharis is most likely positively associated with temperature-dependent speciation and dependency on the Asian monsoons. We suggest that the warm and humid climate of the mid-Miocene was probably the primary driver of the rapid diversification in Oreocharis, while mountain building of the HDM might have indirectly affected species diversification of the HDM lineage. This study highlights the importance of past climatic changes, combined with mountain building, in creating strong environmental heterogeneity and driving diversification of mountain plants, and suggests that the biodiversity in the HDM cannot directly be attributed to mountain uplift, contrary to many recent speculations.[East Asian monsoons; environmental heterogeneity; Hengduan Mountains; incomplete lineage sorting; Oreocharis; past climate change; rapid diversification; transcriptome.]

Essential Oil Composition and Bioactivity of Two Juniper Species from Bulgaria and Slovakia

Juniperus excelsa M. Bieb and J. sabina L. contain essential oil (EO), while J. sabina also contains podophyllotoxin, which is used as a precursor for anti-cancer drugs. Two studies were conducted. The first assessed the variability in the EO profile and podophyllotoxin concentration of the two junipers, depending on the location and tree gender. The main EO constituents of J. excelsa were α-cedrol, α-limonene and α-pinene, while the constituents in J. sabina were sabinene, terpinen-4-ol, myrtenyl acetate and α-cadinol. The podophyllotoxin yield of 18 J. sabina accessions was 0.07–0.32% (w/w), but this was not found in any of the J. excelsa accessions. The second study assessed the effect of hydrodistillation (Clevenger apparatus) and steam distillation (in a semi-commercial apparatus) on the EO profile and bioactivity. The extraction type did not significantly alter the EO composition. The EO profiles of the two junipers and their accessions were different and may be of interest to the industry utilizing juniper leaf EO. Breeding and selection programs could be developed with the two junipers (protected species) in order to identify chemotypes with (1) a high EO content and desirable composition, and (2) a high concentration of podophyllotoxin in J. sabina. Such chemotypes could be established as agricultural crops for the commercial production of podophyllotoxin and EO.

Alpine speciation and morphological innovations: revelations from a species-rich genus in the northern hemisphere

A large number of studies have attempted to determine the mechanisms driving plant diversity and distribution on a global scale, but the diverse and endemic alpine herbs found in harsh environments, showing adaptive evolution, require more studies. Here, we selected 466 species from the genus Saussurea, one of the northern hemisphere's highest-altitude plant genera with high species richness and striking morphological traits, to explore the mechanisms driving speciation and adaptative evolution. We conducted phylogenetic signals analysis and ancestral character estimation to explore the phylogenetic significance of ecological factors. Moreover, we used spatial simultaneous autoregressive (SAR) error models, modified t-tests and partial regression models to quantify the relative effects of ecological factors and morphological diversity upon diversity and endemism of Saussurea. Phylogenetic analyses reveal that geological influences and climate stability exhibit significant phylogenetic signals and that Saussurea originated at a relatively high elevation. Regression models indicate that geological influences and climatic stability significantly affect the diversity and endemism patterns of Saussurea and its morphological innovations. Moreover, morphological innovations in an area show significant contributions to the local diversity and endemism of Saussurea. We conclude that geological influences (mean altitude and topographic heterogeneity), glacial-interglacial climate stability and phylogenetic conservatism have together promoted the speciation and adaptive evolution of the genus Saussurea. In addition, adaptively morphological innovations of alpine species also promote diversification in local regions. Our findings improve the understanding of the distribution pattern of diversity/endemism and adaptive evolution of alpine specie in the whole northern hemisphere.

Past, present, and future geographic range of the relict Mediterranean and Macaronesian Juniperus phoenicea complex

AIM

The aim of this study is to model the past, current, and future distribution of J. phoenicea s.s., J. turbinata, and J. canariensis, based on bioclimatic variables using a maximum entropy model (Maxent) in the Mediterranean and Macaronesian regions.

LOCATION

Mediterranean and Macaronesian.

TAXON

Cupressaceae, Juniperus.

METHODS

Data on the occurrence of the J. phoenicea complex were obtained from the Global Biodiversity Information Facility (GBIF.org), the literature, herbaria, and the authors' field notes. Bioclimatic variables were obtained from the WorldClim database and Paleoclim. The climate data related to species localities were used for predictions of niches by implementation of Maxent, and the model was evaluated with ENMeval.

RESULTS

The potential niches of Juniperus phoenicea during the Last Interglacial period (LIG), Last Glacial Maximum climate (LGM), and Mid-Holocene (MH) covered 30%, 10%, and almost 100%, respectively, of the current potential niche. Climate warming may reduce potential niches by 30% in RCP2.6 and by 90% in RCP8.5. The potential niches of Juniperus turbinata had a broad circum-Mediterranean and Canarian distribution during the LIG and the MH; its distribution extended during the LGM when it was found in more areas than at present. The predicted warming in scenarios RCP2.6 and RCP8.5 could reduce the current potential niche by 30% and 50%, respectively. The model did not find suitable niches for J. canariensis during the LIG and the LGM, but during the MH its potential niche was 30% larger than at present. The climate warming scenario RCP2.6 indicates a reduction in the potential niche by 30%, while RCP8.5 so indicates a reduction of almost 60%.

MAIN CONCLUSIONS

This research can provide information for increasing the protection of the juniper forest and for counteracting the phenomenon of local extinctions caused by anthropic pressure and climate changes.

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next-generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch-Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Polyploidy Expands the Range of Centaurium (Gentianaceae)

The Mediterranean region is one of the most important worldwide hotspots in terms of number of species and endemism, and multiple hypotheses have been proposed to explain how diversification occurred in this area. The contribution of different traits to the diversification process has been evaluated in different groups of plants. In the case of Centaurium (Gentianaceae), a genus with a center of diversity placed in the Mediterranean region, polyploidy seems to have been an important driver of diversification as more than half of species are polyploids. Moreover, ploidy levels are strongly geographically structured across the range of the genus, with tetraploids distributed towards more temperate areas in the north and hexaploids in more arid areas towards the south. We hypothesize that the diversification processes and biodiversity patterns in Centaurium are explained by the coupled formation of polyploid lineages and the colonization of different areas. A MCC tree from BEAST2 based on three nuclear DNA regions of a total of 26 taxa (full sampling, of 18 species and 8 subspecies) was used to perform ancestral area reconstruction analysis in "BioGeoBEARS." Chromosome evolution was analyzed in chromEvol and diversification in BAMM to estimate diversification rates. Our results suggest that two major clades diverged early from the common ancestor, one most likely in the western Mediterranean and the other in a widespread area including west and central Asia (but with high uncertainty in the exact composition of this widespread area). Most ancestral lineages in the western clade remained in or around the western Mediterranean, and dispersal to other areas (mainly northward and eastward), occurred at the tips. Contrarily, most ancestral lineages in the widespread clade had larger ancestral areas. Polyploidization events in the western clade occurred at the tips of the phylogeny (with one exception of a polyploidization event in a very shallow node) and were mainly tetraploid, while polyploidization events occurred in the widespread clade were at the tips and in an ancestral node of the phylogeny, and were mainly hexaploid. We show how ancestral diploid lineages remained in the area of origin, whereas recent and ancestral polyploidization could have facilitated colonization and establishment in other areas.

Molecular phylogeny and ancestral biogeographic reconstruction of Platanthera subgenus Limnorchis (Orchidaceae) using target capture methods

Platanthera is one of the largest genera of temperate orchids in the Holarctic and exemplifies a lineage that has adaptively radiated into diverse habitats within North America, Asia, Europe, North Africa, Borneo, and Sarawak. Major centers of diversity in this genus are North America and eastern Asia. Despite its diversity, a thorough phylogenetic hypothesis for the genus is lacking because no studies have yet sampled taxa exhaustively or developed a robust molecular toolkit. While there is strong evidence that suggests monophyly of subgenus Limnorchis, most taxa in this group have not been included in a phylogenetic analysis. In this study, we developed a new toolkit for Platanthera consisting of genomic information from 617 low-copy nuclear loci. Using a targeted enrichment approach, we collected high-throughput sequence data in 23 accessions of nine of the 12 diploid species of subgenus Limnorchis and outgroup species across Platanthera. A maximum likelihood analysis resolved a strongly supported monophyletic clade for subgenus Limnorchis. Ancestral biogeographic reconstruction indicated that subgenus Limnorchis originated in western North America ca. 3-4.5 Mya from an ancestor that was widespread in western North America and eastern Asia and subsequently diversified in western North America, followed by dispersal of some species to eastern North America. Our results indicate complex biogeographic connections between Asia and North America, and therefore it suggests that Platanthera is a suitable system to test biogeographic hypotheses over time and space in the Holarctic. Our results are also expected to facilitate further study of diversification and biogeographic spread across Platanthera and lay the groundwork for understanding independent origins, biogeography, and morphological diversification of polyploid species within subgenus Limnorchis.

Genome-wide RAD sequencing resolves the evolutionary history of serrate leaf Juniperus and reveals discordance with chloroplast phylogeny

Juniper (Juniperus) is an ecologically important conifer genus of the Northern Hemisphere, the members of which are often foundational tree species of arid regions. The serrate leaf margin clade is native to topologically variable regions in North America, where hybridization has likely played a prominent role in their diversification. Here we use a reduced-representation sequencing approach (ddRADseq) to generate a phylogenomic data set for 68 accessions representing all 22 species in the serrate leaf margin clade, as well as a number of close and distant relatives, to improve understanding of diversification in this group. Phylogenetic analyses using three methods (SVDquartets, maximum likelihood, and Bayesian) yielded highly congruent and well-resolved topologies. These phylogenies provided improved resolution relative to past analyses based on Sanger sequencing of nuclear and chloroplast DNA, and were largely consistent with taxonomic expectations based on geography and morphology. Calibration of a Bayesian phylogeny with fossil evidence produced divergence time estimates for the clade consistent with a late Oligocene origin in North America, followed by a period of elevated diversification between 12 and 5 Mya. Comparison of the ddRADseq phylogenies with a phylogeny based on Sanger-sequenced chloroplast DNA revealed five instances of pronounced discordance, illustrating the potential for chloroplast introgression, chloroplast transfer, or incomplete lineage sorting to influence organellar phylogeny. Our results improve understanding of the pattern and tempo of diversification in Juniperus, and highlight the utility of reduced-representation sequencing for resolving phylogenetic relationships in non-model organisms with reticulation and recent divergence.

Phylogeny, Taxonomy, and Biogeography of Pterocarya (Juglandaceae)

Relict species play an important role in understanding the biogeography of intercontinental disjunctions. Pterocarya (a relict genus) is the valuable model taxon for studying the biogeography of East Asian versus southern European/West Asian disjunct patterns. This disjunction has not been as well studied as others (e.g., between Eastern Asia and North America). Several phylogenetic studies on Pterocarya have been conducted, but none have provided a satisfactory phylogenetic resolution. Here, we report the first well-resolved phylogeny of Pterocarya using restriction site-associated DNA sequencing data based on the sampling of all taxa across the entire distribution area of the genus. Taxonomic treatments were also clarified by combining morphological traits. Furthermore, fossil-calibrated phylogeny was used to explore the biogeography of Pterocarya. Our results support the existence of two sections in Pterocarya, which is in accordance with morphological taxonomy. Section Platyptera comprises three species: P. rhoifolia, P. macroptera, and P. delavayi. Section Pterocarya also comprises three species: P. fraxinifolia, P. hupehensis, and P. stenoptera. The divergence between the two sections took place during the early Miocene (20.5 Ma). The formation of the Gobi Desert and climate cooling of northern Siberia in the Middle Miocene (15.7 Ma) might have caused the split of the continuous distribution of this genus and the formation of the East Asian versus southern European/West Asian disjunct pattern. Lastly, the divergence between P. hupehensis and P. stenoptera as well as between P. rhoifolia and P. macroptera/P. delavayi (10.0 Ma) supports the late Miocene diversification hypothesis in East Asia.

Molecular Identification Based on Chloroplast Sequences and Anti-complementary Activity Comparison of Juniperus Samples from the Qinghai-Tibet Plateau

Juniperus (Cupressaceae, Pinales) plants are widely distributed in the Qinghai-Tibet Plateau of China. The leaves and twigs of at least 8 Juniperus species (J. pingii, J. pingii var. wilsonii, J. squamata, J. recurva var. coxii, J. saltuaria, J. indica, J. tibetica and J. convallium var. microsperma) have been used as the Tibetan medicine Xuba. At present, it is difficult to distinguish among the original species of Xuba based only on their similar morphological characteristics. However, in our previous studies, 4 Xuba samples from different Juniperus species exhibited significant differences in both anticomplementary activity in vitro and anti-inflammatory effects on acute lung injury in vivo. To identify the effective original species of Xuba reliably, in this study, we developed a sequencing-based DNA molecular technology to distinguish 14 populations of 8 Juniperus species collected from Tibet region, using trnS-G, trnD - T, and petN-psbM genomic regions to build phylogenetic trees. In addition, their anticomplementary activities were evaluated. The results showed that combined sequence of these 3 genomic regions could identify 8 Juniperus species clearly and clustered individuals of one species but from different locations, whichever phylogenetic tree was constructed. Moreover, the anticomplementary activities of the 8 species were clustered into 2 groups. Among them, J. saltuaria and J. recurva var. coxii, which formed an independent branch apart from the other 6 species in phylogenetic trees, were the most potent (CH₅₀: 0.029 - 0.032 mg/mL). Consequently, DNA identification of Juniperus using the combined sequence could provide beneficial guidance for further efficacy evaluation and quality control of Xuba.

Plastome Phylogenomic and Biogeographical Study on Thuja (Cupressaceae)

Investigating the biogeographical disjunction of East Asian and North American flora is key to understanding the formation and dynamics of biodiversity in the Northern Hemisphere. The small Cupressaceae genus Thuja, comprising five species, exhibits a typical disjunct distribution in East Asia and North America. Owing to obscure relationships, the biogeographical history of the genus remains controversial. Here, complete plastomes were employed to investigate the plastome evolution, phylogenetic relationships, and biogeographic history of Thuja. All plastomes of Thuja share the same gene content arranged in the same order. The loss of an IR was evident in all Thuja plastomes, and the B-arrangement as previously recognized was detected. Phylogenomic analyses resolved two sister pairs, T. standishii-T. koraiensis and T. occidentalis-T. sutchuenensis, with T. plicata sister to T. occidentalis-T. sutchuenensis. Molecular dating and biogeographic results suggest the diversification of Thuja occurred in the Middle Miocene, and the ancestral area of extant species was located in northern East Asia. Incorporating the fossil record, we inferred that Thuja likely originated from the high-latitude areas of North America in the Paleocene with a second diversification center in northern East Asia. The current geographical distribution of Thuja was likely shaped by dispersal events attributed to the Bering Land Bridge in the Miocene and subsequent vicariance events accompanying climate cooling. The potential effect of extinction may have profound influence on the biogeographical history of Thuja.

Patterns of genetic diversity in North Africa: Moroccan-Algerian genetic split in Juniperus thurifera subsp. africana

Juniperus thurifera is a key element of the forest communities in arid and semi-arid areas of the western Mediterranean. Previous genetic and morphological investigations suggested that Algerian populations are genetically more similar to European than to Moroccan populations and advocated their recognition at the variety rank. We aimed to investigate the spatial genetic structure in J. thurifera to verify the distinct character of the Algerian population in terms of the genetic breaks reported among several North African taxa. We also modelled species distributions since the Eemian to recognise the impact of past climatic changes on the current pattern of diversity and predict possible changes in species distribution in the future. Species-specific microsatellites were used in the analysis of 11 populations from Algeria, Morocco and Europe. We revealed the significant genetic distinctiveness of the Algerian populations from the Moroccan and European stands that may have important taxonomic and conservation implications. The diversity pattern revealed for J. thurifera reflects the east-west genetic splits reported among some North African plant and animal taxa and suggests an impact of shared historical processes. Additionally, modelling of the distribution allowed us to identify possible glacial refugia and their impact on the modern pattern of differentiation in J. thurifera. Reduction of species occurrence, especially in the European domain, is likely according to the future projections of the species distribution.