Phylogenomics resolves evolutionary relationships and provides insights into floral evolution in the tribe Shoreeae (Dipterocarpaceae)

From sequences to sustainability: Exploring dipterocarp genomes for oleoresin production, timber quality, and conservation

Dipterocarp species dominate tropical forest ecosystems and provide key ecological and economic value through their use of aromatic resins, medicinal chemicals, and high-quality timber. However, habitat loss and unsustainable logging have endangered many Dipterocarpaceae species. Genomic strategies provide new opportunities for both elucidating the molecular pathways underlying these desirable traits and informing conservation efforts for at-risk taxa. This review summarizes the progress in dipterocarp genomics analysis and applications. We describe 16 recently published Dipterocarpaceae genome sequences, representing crucial genetic blueprints. Phylogenetic comparisons delineate evolutionary relationships among species and provide frameworks for pinpointing functional changes underlying specialized metabolism and wood development patterns. We also discuss connections revealed thus far between specific gene families and both oleoresin biosynthesis and wood quality traits-including the identification of key terpenoid synthases and cellulose synthases likely governing pathway flux. Moreover, the characterization of adaptive genomic markers offers vital resources for supporting conservation practices prioritizing resilient genotypes displaying valuable oleoresin and timber traits. Overall, progress in dipterocarp functional and comparative genomics provides key tools for addressing the intertwined challenges of preserving biodiversity in endangered tropical forest ecosystems while sustainably deriving aromatic chemicals and quality lumber that support diverse human activities.

Pervasive Introgression During Rapid Diversification of the European Mountain Genus Soldanella (L.) (Primulaceae)

Hybridization is a key mechanism involved in lineage diversification and speciation, especially in ecosystems that experienced repeated environmental oscillations. Recently radiated plant groups, which have evolved in mountain ecosystems impacted by historical climate change provide an excellent model system for studying the impact of gene flow on speciation. We combined organellar (whole-plastome) and nuclear genomic data (RAD-seq) with a cytogenetic approach (rDNA FISH) to investigate the effects of hybridization and introgression on evolution and speciation in the genus Soldanella (snowbells, Primulaceae). Pervasive introgression has already occurred among ancestral lineages of snowbells and has persisted throughout the entire evolutionary history of the genus, regardless of the ecology, cytotype, or distribution range size of the affected species. The highest extent of introgression has been detected in the Carpathian species, which is also reflected in their extensive karyotype variation. Introgression occurred even between species with dysploid and euploid cytotypes, which were considered to be reproductively isolated. The magnitude of introgression detected in snowbells is unprecedented in other mountain genera of the European Alpine System investigated hitherto. Our study stresses the prominent evolutionary role of hybridization in facilitating speciation and diversification on the one hand, but also enriching previously isolated genetic pools. [chloroplast capture; diversification; dysploidy; European Alpine system; introgression; nuclear-cytoplasmic discordance; ribosomal DNA.].

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Down, then up: non-parallel genome size changes and a descending chromosome series in a recent radiation of the Australian allotetraploid plant species, Nicotiana section Suaveolentes (Solanaceae)

BACKGROUND AND AIMS

The extent to which genome size and chromosome numbers evolve in concert is little understood, particularly after polyploidy (whole-genome duplication), when a genome returns to a diploid-like condition (diploidization). We study this phenomenon in 46 species of allotetraploid Nicotiana section Suaveolentes (Solanaceae), which formed <6 million years ago and radiated in the arid centre of Australia.

METHODS

We analysed newly assessed genome sizes and chromosome numbers within the context of a restriction site-associated nuclear DNA (RADseq) phylogenetic framework.

KEY RESULTS

RADseq generated a well-supported phylogenetic tree, in which multiple accessions from each species formed unique genetic clusters. Chromosome numbers and genome sizes vary from n = 2x = 15 to 24 and 2.7 to 5.8 pg/1C nucleus, respectively. Decreases in both genome size and chromosome number occur, although neither consistently nor in parallel. Species with the lowest chromosome numbers (n = 15-18) do not possess the smallest genome sizes and, although N. heterantha has retained the ancestral chromosome complement, n = 2x = 24, it nonetheless has the smallest genome size, even smaller than that of the modern representatives of ancestral diploids.

CONCLUSIONS

The results indicate that decreases in genome size and chromosome number occur in parallel down to a chromosome number threshold, n = 20, below which genome size increases, a phenomenon potentially explained by decreasing rates of recombination over fewer chromosomes. We hypothesize that, more generally in plants, major decreases in genome size post-polyploidization take place while chromosome numbers are still high because in these stages elimination of retrotransposons and other repetitive elements is more efficient. Once such major genome size change has been accomplished, then dysploid chromosome reductions take place to reorganize these smaller genomes, producing species with small genomes and low chromosome numbers such as those observed in many annual angiosperms, including Arabidopsis.

Thirteen Dipterocarpoideae genomes provide insights into their evolution and borneol biosynthesis

Dipterocarpoideae, the largest subfamily of the Dipterocarpaceae, is a dominant component of Southeast Asian rainforests and is widely used as a source of wood, damar resin, medicine, and essential oil. However, many Dipterocarpoideae species are currently on the IUCN Red List owing to severe degradation of their habitats under global climate change and human disturbance. Genetic information regarding these taxa has only recently been reported with the sequencing of four Dipterocarp genomes, providing clues to the function and evolution of these species. Here, we report on 13 high-quality Dipterocarpoideae genome assemblies, ranging in size from 302.6 to 494.8 Mb and representing the five most species-rich genera in Dipterocarpoideae. Molecular dating analyses support the Western Gondwanaland origin of Dipterocarpaceae. Based on evolutionary analysis, we propose a three-step chromosome evolution scenario to describe the karyotypic evolution from an ancestor with six chromosomes to present-day species with 11 and 7 chromosomes. We discovered an expansion of genes encoding cellulose synthase (CesA), which is essential for cellulose biosynthesis and secondary cell-wall formation. We functionally identified five bornyl diphosphate synthase (BPPS) genes, which specifically catalyze the biosynthesis of borneol, a natural medicinal compound extracted from damar resin and oils, thus providing a basis for large-scale production of natural borneol in vitro.

A chromosome-scale genome and transcriptomic analysis of the endangered tropical tree Vatica mangachapoi (Dipterocarpaceae)

Vatica mangachapoi is a tropical tree species native to Southeast Asia. It has long been valued as a timber species because the wood resists decay, but it is now considered vulnerable to extinction due to habitat loss and overexploitation. Here, we present the first chromosome-level genome assembly of V. mangachapoi that we created by combining data from PacBio long read sequencing with Hi-C proximity ligation and Illumina short-read sequencing. The assembled genome was 456.21 Mb, containing 11 chromosome and a BUSCO score of 93.4%. From the newly assembled genome, 46,811 protein-coding genes were predicted. Repetitive DNA accounted for 53% of the genome. Phylogenomic and gene family analyses showed that V. mangachapoi diverged from a common ancestor of Gossypium raimondii 70 million years ago. Transcriptome analyses found 227 genes that were differentially expressed in the leaves of plants grown in normal soil relative to plants grown in dry, coastal, sandy soil. For these genes, we identified three significantly enriched with GO terms: responses to organonitrogen compounds, chitin-triggered immunity, and wound response. This genome provides an important comparative benchmark not only for future conservation work on V. mangachapoi but also for phylogenomics work on Dipterocarpaceae.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

Rhododendron kuomeianum (Ericaceae), a new species from northeastern Yunnan (China), based on morphological and genomic data

Rhododendron kuomeianum Y.H. Chang, J. Nielsen & Y.P. Ma, a new species of Rhododendron (Ericaceae) within subsect. Maddenia in sect. Rhododendron from Yiliang County, NE Yunnan, China, is described and illustrated. The new species is similar to R. valentinianum, but it can be easily distinguished by its sparse scales on the abaxial surface of the leaf blade, fewer flowers per inflorescence and white corolla with pale red margins. There are also differences in the widths of calyx lobes, leaf blade shape and indumentum characteristics of the petiole between the new species and Rhododendron linearilobum. We confirmed that R. kuomeianum is a new species closely related to R. valentinianum and R. changii with phylogenomic studies of 10 species within this subsection based on restriction site-associated DNA sequencing (RAD-seq) data. These phylogenomic analyses also clarified additional taxonomic problems in this subsection previously raised by morphological analysis. Our findings make a strong case for using next-generation sequencing to explore phylogenetic relationships and identify new species, especially in plants groups with complicated taxonomic problems.

Spatial and Ecological Drivers of Genetic Structure in Greek Populations of Alkanna tinctoria (Boraginaceae), a Polyploid Medicinal Herb

Background and Aims: Quantifying genetic variation is fundamental to understand a species' demographic trajectory and its ability to adapt to future changes. In comparison with diploids, however, genetic variation and factors fostering genetic divergence remain poorly studied in polyploids due to analytical challenges. Here, by employing a ploidy-aware framework, we investigated the genetic structure and its determinants in polyploid Alkanna tinctoria (Boraginaceae), an ancient medicinal herb that is the source of bioactive compounds known as alkannin and shikonin (A/S). From a practical perspective, such investigation can inform biodiversity management strategies. Methods: We collected 14 populations of A. tinctoria within its main distribution range in Greece and genotyped them using restriction site-associated DNA sequencing. In addition, we included two populations of A. sieberi. By using a ploidy-aware genotype calling based on likelihoods, we generated a dataset of 16,107 high-quality SNPs. Classical and model-based analysis was done to characterize the genetic structure within and between the sampled populations, complemented by genome size measurements and chromosomal counts. Finally, to reveal the drivers of genetic structure, we searched for associations between allele frequencies and spatial and climatic variables. Key Results: We found support for a marked regional structure in A. tinctoria along a latitudinal gradient in line with phytogeographic divisions. Several analyses identified interspecific admixture affecting both mainland and island populations. Modeling of spatial and climatic variables further demonstrated a larger contribution of neutral processes and a lesser albeit significant role of selection in shaping the observed genetic structure in A. tinctoria. Conclusion: Current findings provide evidence of strong genetic structure in A. tinctoria mainly driven by neutral processes. The revealed natural genomic variation in Greek Alkanna can be used to further predict variation in A/S production, whereas our bioinformatics approach should prove useful for the study of other non-model polyploid species.

Phylogenomics resolves deep subfamilial relationships in Malvaceae s.l

Malvaceae s.l., the most diverse family within Malvales, includes well-known species of great economic importance like cotton, cacao, and durian. Despite numerous phylogenetic analyses employing multiple markers, relationships between several of its nine subfamilies, particularly within the largest lineage /Malvadendrina, remain unclear. In this study, we attempted to resolve the relationships within the major clades of Malvaceae s.l. using plastid genomes of 48 accessions representing all subfamilies. Maximum likelihood and Bayesian analyses recovered a fully resolved and well-supported topology confirming the split of the family into /Byttneriina (/Grewioideae +/Byttnerioideae) and /Malvadendrina. Within /Malvadendrina, /Helicteroideae occupied the earliest branching position, followed by /Sterculioideae, /Brownlowioideae, /Tiliodeae, and /Dombeyoideae formed a clade sister to /Malvatheca (/Malvoideae +/Bombacoideae), a grouping morphologically supported by the lack of androgynophore. Results from dating analyses suggest that all subfamilies originated during hot or warm phases in the Late Cretaceous to Paleocene. This study presents a well-supported phylogenetic framework for Malvaceae s.l. that will aid downstream revisions and evolutionary studies of this economically important plant family.

Pollination syndromes in the 21st century: where do we stand and where may we go?

Pollination syndromes, recurring suites of floral traits appearing in connection with specific functional pollinator groups, have served for decades to organise floral diversity under a functional-ecological perspective. Some potential caveats, such as over-simplification of complex plant-animal interactions or lack of empirical observations, have been identified and discussed in recent years. Which of these caveats do indeed cause problems, which have been solved and where do future possibilities lie? I address these questions in a review of the pollination-syndrome literature of 2010 to 2019. I show that the majority of studies was based on detailed empirical pollinator observations and could reliably predict pollinators based on a few floral traits such as colour, shape or reward. Some traits (i.e. colour) were less reliable in predicting pollinators than others (i.e. reward, corolla width), however. I stress that future studies should consider floral traits beyond those traditionally recorded to expand our understanding of mechanisms of floral evolution. I discuss statistical methods suitable for objectively analysing the interplay of system-specific evolutionary constraints, pollinator-mediated selection and adaptive trade-offs at microecological and macroecological scales. I exemplify my arguments on an empirical dataset of floral traits of a neotropical plant radiation in the family Melastomataceae.

Phylogenomic Relationships of Diploids and the Origins of Allotetraploids in Dactylorhiza (Orchidaceae)

Disentangling phylogenetic relationships proves challenging for groups that have evolved recently, especially if there is ongoing reticulation. Although they are in most cases immediately isolated from diploid relatives, sets of sibling allopolyploids often hybridize with each other, thereby increasing the complexity of an already challenging situation. Dactylorhiza (Orchidaceae: Orchidinae) is a genus much affected by allopolyploid speciation and reticulate phylogenetic relationships. Here, we use genetic variation at tens of thousands of genomic positions to unravel the convoluted evolutionary history of Dactylorhiza. We first investigate circumscription and relationships of diploid species in the genus using coalescent and maximum likelihood methods, and then group 16 allotetraploids by maximum affiliation to their putative parental diploids, implementing a method based on genotype likelihoods. The direction of hybrid crosses is inferred for each allotetraploid using information from maternally inherited plastid RADseq loci. Starting from age estimates of parental taxa, the relative ages of these allotetraploid entities are inferred by quantifying their genetic similarity to the diploids and numbers of private alleles compared with sibling allotetraploids. Whereas northwestern Europe is dominated by young allotetraploids of postglacial origins, comparatively older allotetraploids are distributed further south, where climatic conditions remained relatively stable during the Pleistocene glaciations. Our bioinformatics approach should prove effective for the study of other naturally occurring, nonmodel, polyploid plant complexes.

Integrating DNA Barcoding and Traditional Taxonomy for the Identification of Dipterocarps in Remnant Lowland Forests of Sumatra

DNA barcoding has been used as a universal tool for phylogenetic inferences and diversity assessments, especially in poorly studied species and regions. The aim of this study was to contrast morphological taxonomy and DNA barcoding, using the three frequently used markers matK, rbcL, and trnL-F, to assess the efficiency of DNA barcoding in the identification of dipterocarps in Sumatra, Indonesia. The chloroplast gene matK was the most polymorphic among these three markers with an average interspecific genetic distance of 0.020. The results of the molecular data were mostly in agreement with the morphological identification for the clades of Anthoshorea, Hopea, Richetia, Parashorea, and Anisoptera, nonetheless these markers were inefficient to resolve the relationships within the Rubroshorea group. The maximum likelihood and Bayesian inference phylogenies identified Shorea as a paraphyletic genus, Anthoshorea appeared as sister to Hopea, and Richetia was sister to Parashorea. A better discriminatory power among dipterocarp species provided by matK and observed in our study suggests that this marker has a higher evolutionary rate than the other two markers tested. However, a combination of several different barcoding markers is essential for reliable identification of the species at a lower taxonomic level.

Iterative allogamy-autogamy transitions drive actual and incipient speciation during the ongoing evolutionary radiation within the orchid genus Epipactis (Orchidaceae)

BACKGROUND AND AIMS

The terrestrial orchid genus Epipactis has become a model system for the study of speciation via transitions from allogamy to autogamy, but close phylogenetic relationships have proven difficult to resolve through Sanger sequencing.

METHODS

We analysed with restriction site-associated sequencing (RAD-seq) 108 plants representing 29 named taxa that together span the genus, focusing on section Epipactis. Our filtered matrix of 12 543 single nucleotide polymorphisms was used to generate an unrooted network and a rooted, well-supported likelihood tree. We further inferred genetic structure through a co-ancestry heat map and admixture analysis, and estimated inbreeding coefficients per sample.

KEY RESULTS

The 27 named taxa of the ingroup were resolved as 11 genuine, geographically widespread species: four dominantly allogamous and seven dominantly autogamous. A single comparatively allogamous species, E. helleborine, is the direct ancestor of most of the remaining species, though one of the derived autogams has generated one further autogamous species. An assessment of shared ancestry suggested only sporadic hybridization between the re-circumscribed species. Taxa with the greatest inclination towards autogamy show less, if any, admixture, whereas the gene pools of more allogamous species contain a mixture alleles found in the autogams.

CONCLUSIONS

This clade is presently undergoing an evolutionary radiation driven by a wide spectrum of genotypic, phenotypic and environmental factors. Epipactis helleborine has also frequently generated many local variants showing inclinations toward autogamy (and occasionally cleistogamy), best viewed as incipient speciation from within the genetic background provided by E. helleborine, which thus becomes an example of a convincingly paraphyletic species. Autogams are often as widespread and ecologically successful as allogams.

Exploring evolution and diversity of Chinese Dipterocarpaceae using next-generation sequencing

Tropical forests, a key-category of land ecosystems, are faced with the world's highest levels of habitat conversion and associated biodiversity loss. In tropical Asia, Dipterocarpaceae are one of the economically and ecologically most important tree families, but their genomic diversity and evolution remain understudied, hampered by a lack of available genetic resources. Southern China represents the northern limit for Dipterocarpaceae, and thus changes in habitat ecology, community composition and adaptability to climatic conditions are of particular interest in this group. Phylogenomics is a tool for exploring both biodiversity and evolutionary relationships through space and time using plastome, nuclear and mitochondrial genome. We generated full plastome and Nuclear Ribosomal Cistron (NRC) data for Chinese Dipterocarpaceae species as a first step to improve our understanding of their ecology and evolutionary relationships. We generated the plastome of Dipterocarpus turbinatus, the species with the widest distribution using it as a baseline for comparisons with other taxa. Results showed low level of genomic diversity among analysed range-edge species, and different evolutionary history of the incongruent NRC and plastome data. Genomic resources provided in this study will serve as a starting point for future studies on conservation and sustainable use of these dominant forest taxa, phylogenomics and evolutionary studies.

Multiple auto- and allopolyploidisations marked the Pleistocene history of the widespread Eurasian steppe plant Astragalus onobrychis (Fabaceae)

The Eurasian steppes occupy a significant portion of the worldwide land surface and their biota have been affected by specific past range dynamics driven by ice ages-related climatic fluctuations. The dynamic alterations in conditions during the Pleistocene often triggered reticulate evolution and whole genome duplication events. Employing genomic, genetic and cytogenetic tools as well as morphometry we investigate the intricate evolution of Astragalus onobrychis, a widespread Eurasian steppe plant with diploid, tetraploid and octoploid cytotypes. To analyse the heteroploid RADseq dataset we employ both genotype-based and genotype-free methods that result in highly consistent results, and complement our inference with information from the plastid ycf1 region. We uncover a complex and reticulate evolutionary history, including at least one auto-tetraploidization event and two allo-octoploidization events; one of them involved also genetic contributions from other species, most likely A. goktschaicus. The present genetic structure points to the existence of four main clades within A. onobrychis, which only partly correspond to different ploidies. Time-calibrated diffusion models suggest that diversification within A. onobrychis was associated with ice age-related climatic fluctuations during the last million years. We finally argue for the usefulness of uniparentally inherited plastid markers, even in the genomic era, especially when investigating heteroploid systems.

Whole genomes: the holy grail. A commentary on: 'Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastidgenomes'

This article comments on: Jacqueline Heckenhauer, Ovidiu Paun, Mark W. Chase, Peter S. Ashton, A.S. Kamariah and Rosabelle Samuel. 2019. Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastid genomes. Annals of Botany 123(5): 857–865.

Molecular phylogenomics of the tribe Shoreeae (Dipterocarpaceae) using whole plastid genomes

BACKGROUND AND AIMS

Phylogenetic relationships within tribe Shoreeae, containing the main elements of tropical forests in Southeast Asia, present a long-standing problem in the systematics of Dipterocarpaceae. Sequencing whole plastomes using next-generation sequencing- (NGS) based genome skimming is increasingly employed for investigating phylogenetic relationships of plants. Here, the usefulness of complete plastid genome sequences in resolving phylogenetic relationships within Shoreeae is evaluated.

METHODS

A pipeline to obtain alignments of whole plastid genome sequences across individuals with different amounts of available data is presented. In total, 48 individuals, representing 37 species and four genera of the ecologically and economically important tribe Shoreeae sensu Ashton, were investigated. Phylogenetic trees were reconstructed using maximum parsimony, maximum likelihood and Bayesian inference.

KEY RESULTS

Here, the first fully sequenced plastid genomes for the tribe Shoreeae are presented. Their size, GC content and gene order are comparable with those of other members of Malvales. Phylogenomic analyses demonstrate that whole plastid genomes are useful for inferring phylogenetic relationships among genera and groups of Shorea (Shoreeae) but fail to provide well-supported phylogenetic relationships among some of the most closely related species. Discordance in placement of Parashorea was observed between phylogenetic trees obtained from plastome analyses and those obtained from nuclear single nucleotide polymorphism (SNP) data sets identified in restriction-site associated sequencing (RADseq).

CONCLUSIONS

Phylogenomic analyses of the entire plastid genomes are useful for inferring phylogenetic relationships at lower taxonomic levels, but are not sufficient for detailed phylogenetic reconstructions of closely related species groups in Shoreeae. Discordance in placement of Parashorea was further investigated for evidence of ancient hybridization.