A transcriptome-based resolution for a key taxonomic controversy in Cupressaceae

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.

Efficient purging of deleterious mutations contributes to the survival of a rare conifer

Cupressaceae is a conifer family rich in plants of horticultural importance, including Cupressus, Chamaecyparis, Juniperus, and Thuja, yet genomic surveys are lacking for this family. Cupressus gigantea, one of the many rare conifers that are threatened by climate change and anthropogenic habitat fragmentation, plays an ever-increasing role in ecotourism in Tibet. To infer how past climate change has shaped the population evolution of this species, we generated a de novo chromosome-scale genome (10.92 Gb) and compared the species' population history and genetic load with that of a widespread close relative, C. duclouxiana. Our demographic analyses, based on 83 re-sequenced individuals from multiple populations of the two species, revealed a sharp decline of population sizes during the first part of the Quaternary. However, populations of C. duclouxiana then started to recover, while C. gigantea populations continued to decrease until recently. The total genomic diversity of C. gigantea is smaller than that of C. duclouxiana, but contrary to expectations, C. gigantea has fewer highly and mildly deleterious mutations than C. duclouxiana, and simulations and statistical tests support purifying selection during prolonged inbreeding as the explanation. Our results highlight the evolutionary consequences of decreased population size on the genetic burden of a long-lived endangered conifer with large genome size and suggest that genetic purging deserves more attention in conservation management.

Phylogeny and evolution of Cupressaceae: Updates on intergeneric relationships and new insights on ancient intergeneric hybridization

After the merger of the former Taxodiaceae and Cupressaceae s.s., currently the conifer family Cupressaceae (sensu lato) comprises seven subfamilies and 32 genera, most of which are important components of temperate and mountainous forests. With the exception of a recently published genus-level phylogeny of gymnosperms inferred from sequence analysis of 790 orthologs, previous phylogenetic studies of Cupressaceae were based mainly on morphological characters or a few molecular markers, and did not completely resolve the intergeneric relationships. In this study, we reconstructed a robust and well-resolved phylogeny of Cupressaceae represented by all 32 genera, using 1944 genes (Orthogroups) generated from transcriptome sequencing. Reticulate evolution analyses detected a possible ancient hybridization that occurred between ancestors of two subclades of Cupressoideae, including Microbiota-Platycladus-Tetraclinis (MPT) and Juniperus-Cupressus-Hesperocyparis-Callitropsis-Xanthocyparis (JCHCX), although both concatenation and coalescent trees are highly supported. Moreover, divergence time estimation and ancestral area reconstruction indicate that Cupressaceae very likely originated in Asia in the Triassic, and geographic isolation caused by continental separation drove the vicariant evolution of the two subfamilies Cupressoideae and Callitroideae in the northern and southern hemispheres, respectively. Evolutionary analyses of some morphological characters suggest that helically arranged linear-acicular leaves and imbricate bract-scale complexes represent ancestral states, and the shift from linear-acicular leaves to scale-like leaves was associated with the shift from helical to decussate arrangement. Our study sheds new light on phylogeny and evolutionary history of Cupressaceae, and strongly suggests that both dichotomous phylogenetic and reticulate evolution analyses be conducted in phylogenomic studies.

Incomplete lineage sorting and reticulate evolution mask species relationships in Brunelliaceae, an Andean family with rapid, recent diversification

PREMISE

To date, phylogenetic relationships within the monogeneric Brunelliaceae have been based on morphological evidence, which does not provide sufficient phylogenetic resolution. Here we use target-enriched nuclear data to improve our understanding of phylogenetic relationships in the family.

METHODS

We used the Angiosperms353 toolkit for targeted recovery of exonic regions and supercontigs (exons + introns) from low copy nuclear genes from 53 of 70 species in Brunellia, and several outgroup taxa. We removed loci that indicated biased inference of relationships and applied concatenated and coalescent methods to infer Brunellia phylogeny. We identified conflicts among gene trees that may reflect hybridization or incomplete lineage sorting events and assessed their impact on phylogenetic inference. Finally, we performed ancestral-state reconstructions of morphological traits and assessed the homology of character states used to define sections and subsections in Brunellia.

RESULTS

Brunellia comprises two major clades and several subclades. Most of these clades/subclades do not correspond to previous infrageneric taxa. There is high topological incongruence among the subclades across analyses.

CONCLUSIONS

Phylogenetic reconstructions point to rapid species diversification in Brunelliaceae, reflected in very short branches between successive species splits. The removal of putatively biased loci slightly improves phylogenetic support for individual clades. Reticulate evolution due to hybridization and/or incomplete lineage sorting likely both contribute to gene-tree discordance. Morphological characters used to define taxa in current classification schemes are homoplastic in the ancestral character-state reconstructions. While target enrichment data allows us to broaden our understanding of diversification in Brunellia, the relationships among subclades remain incompletely understood.

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.

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.

Phylotranscriptomics of Theaceae: generic-level relationships, reticulation and whole-genome duplication

BACKGROUND AND AIMS

Theaceae, with three tribes, nine genera and more than 200 species, are of great economic and ecological importance. Recent phylogenetic analyses based on plastomic data resolved the relationships among the three tribes and the intergeneric relationships within two of those tribes. However, generic-level relationships within the largest tribe, Theeae, were not fully resolved. The role of putative whole-genome duplication (WGD) events in the family and possible hybridization events among genera within Theeae also remain to be tested further.

METHODS

Transcriptomes or low-depth whole-genome sequencing of 57 species of Theaceae, as well as additional plastome sequence data, were generated. Using a dataset of low-copy nuclear genes, we reconstructed phylogenetic relationships using concatenated, species tree and phylogenetic network approaches. We further conducted molecular dating analyses and inferred possible WGD events by examining the distribution of the number of synonymous substitutions per synonymous site (Ks) for paralogues in each species. For plastid protein-coding sequences , phylogenies were reconstructed for comparison with the results obtained from analysis of the nuclear dataset.

RESULTS

Based on the 610 low-copy nuclear genes (858 606 bp in length) investigated, Stewartieae was resolved as sister to the other two tribes. Within Theeae, the Apterosperma-Laplacea clade grouped with Pyrenaria, leaving Camellia and Polyspora as sister. The estimated ages within Theaceae were largely consistent with previous studies based mainly on plastome data. Two reticulation events within Camellia and one between the common ancestor of Gordonia and Schima were found. All members of the tea family shared two WGD events, an older At-γ and a recent Ad-β; both events were also shared with the outgroups (Diapensiaceae, Pentaphylacaceae, Styracaceae and Symplocaceae).

CONCLUSIONS

Our analyses using low-copy nuclear genes improved understanding of phylogenetic relationships at the tribal and generic levels previously proposed based on plastome data, but the phylogenetic position of the Apterosperma-Laplacea clade needs more attention. There is no evidence for extensive intergeneric hybridization within Theeae or for a Theaceae-specific WGD event. Land bridges (e.g. the Bering land bridge) during the Late Oligocene may have permitted the intercontinental plant movements that facilitated the putative ancient introgression between the common ancestor of Gordonia and Schima.

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.

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.

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

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 phylogenetic analyses reveal a deep dichotomy in the conifer-inhabiting genus Trisetacus (Eriophyoidea: Nalepellidae), with the two lineages differing in their female genital morphology and host associations

We analyzed the phylogenetic relationships of the genus Trisetacus using two genes [cytochrome c oxidase subunit I (COI) and D1-D2 region of 28S rDNA (D1-D2 28S)], a representive taxon sampling (nearly 40% of known diversity), and a large set of close and distant outgroups. Our analyses suggest the presence of a dichotomy between Trisetacus associated with Cupressaceae and Pinaceae. The following smaller molecular clades were found: Pin-1 (bud mites, twig sheath mites, bark gall mites, and endoparasitic mites from pinaceans), Pin-2 (needle sheath mites from pines), Pin-2a (putative Nearctic group of needle sheath mites), Pin-2b (putative Palearctic group of needle sheath mites), Cup-1 and 2 (bud, cone, seed mites and mites living under bark scales from cupressaceans). The monophyly of the recently proposed subgenus Brevithecus nested within clade Cup-2 was confirmed. Ancestral character reconstruction analyses recovered: (1) Pinaceae as the ancestral hosts of Nalepellidae and Trisetacus, (2) repetitive reductions of the spermathecal tube independently occurred in two lineages of Trisetacus from Cupressaceae, and (3) several mite habitats on host (galls, cones, twig sheaths, seeds, inside leaves, and under scales) are evolutionarily derived states, whereas living in buds or needle sheaths are ancestral states for Trisetacus clades Cup and Pin. Using confocal microscopy, we identified six basic types of the female internal genitalia of Trisetacus based on shapes of the spermatheca and spermathecal tube. These genitalic types are strongly correlated with lineages recovered by molecular phylogenetic analyses, suggesting that the female genital morphology is both evolutionarily conserved and is a factor influencing macroevolutionary patterns in this group of mites.

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.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

In-depth transcriptome characterization uncovers distinct gene family expansions for Cupressus gigantea important to this long-lived species' adaptability to environmental cues

Background

Cupressus gigantea, a rare and endangered tree species with remarkable medicinal value, is endemic to the Tibetan Plateau. Yet, little is known about the underlying genetics of the unique ecological adaptability of this extremely long-lived conifer with a large genome size. Here, we present its first de novo and multi-tissue transcriptome in-depth characterization.

Results

We performed Illumina paired-end sequencing and RNA libraries assembly derived from terminal buds, male and female strobili, biennial leaves, and cambium tissues taken from adult C. gigantea. In total, large-scale high-quality reads were assembled into 101,092 unigenes, with an average sequence length of 1029 bp, and 6848 unigenes (6.77%) were mapped against the KEGG databases to identify 292 pathways. A core set of 41,373 genes belonging to 2412 orthologous gene families shared between C. gigantea and nine other plants was revealed. In addition, we identified 2515 small to larger-size gene families containing in total 9223 genes specific to C. gigantea, and enriched for gene ontologies relating to biotic interactions. We identified an important terpene synthases gene family expansion with its 121 putative members.

Conclusions

This study presents the first comprehensive transcriptome characterization of C. gigantea. Our results will facilitate functional genomic studies to support genetic improvement and conservation programs for this endangered conifer.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5584-6) contains supplementary material, which is available to authorized users.