Phylogeny, historical biogeography, and diversification of angiosperm order Ericales suggest ancient Neotropical and East Asian connections

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense

BACKGROUND

Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun.

RESULTS

A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM).

CONCLUSION

Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

Phylogeographical and population genetics of Polyspora sweet in China provides insights into its phylogenetic evolution and subtropical dispersal

BACKGROUND

Geological movements and climatic fluctuations stand as pivotal catalysts driving speciation and phylogenetic evolution. The genus Polyspora Sweet (Theaceae), prominently found across the Malay Archipelagos and Indochina Peninsula in tropical Asia, exhibits its northernmost distribution in China. In this study, we investigated the evolutionary and biogeographical history of the genus Polyspora in China, shedding light on the mechanisms by which these species respond to ancient geological and climatic fluctuations.

METHODS

Phylogenetic relationships of 32 representative species of Theaceae were reconstructed based on the chloroplast genome and ribosome 18-26 S rRNA datasets. Species divergence time was estimated using molecular clock and five fossil calibration. The phylogeography and population genetics in 379 individuals from 32 populations of eight species were analyzed using chloroplast gene sequences (trnH-psbA, rpoB-trnC and petN-psbM), revealing the glacial refugia of each species, and exploring the causes of the phylogeographic patterns.

RESULTS

We found that Chinese Polyspora species diverged in the middle Miocene, showing a tropical-subtropical divergence order. A total of 52 haplotypes were identified by the combined chloroplast sequences. Chinese Polyspora exhibited a distinct phylogeographical structure, which could be divided into two clades and eight genealogical subdivisions. The divergence between the two clades occurred approximately 20.67 Ma. Analysis of molecular variance revealed that the genetic variation mainly occurred between species (77.91%). At the species level, Polyspora axillaris consists of three lineages, while P. speciosa had two lineages. The major lineages of Chinese Polyspora diverged between 12 and 15 Ma during the middle to late Miocene. The peak period of haplotype differentiation in each species occurred around the transition from the last interglacial to the last glacial period, approximately 6 Ma ago.

CONCLUSION

The primary geographical distribution pattern of Chinese Polyspora was established prior to the last glacial maximum, and the population historical dynamics were relatively stable. The geological and climatic turbulence during the Quaternary glacial period had minimal impact on the distribution pattern of the genus. The genus coped with Quaternary climate turbulence by glacial in situ survival in multiple refuges. The Sino-Vietnam border and Nanling corridor might be the genetic mixing center of Polyspora.

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea)

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore 'traditional' taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 "barcode" sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

Decoupling speciation and extinction reveals both abiotic and biotic drivers shaped 250 million years of diversity in crocodile-line archosaurs

Whereas living representatives of Pseudosuchia, crocodylians, number fewer than 30 species, more than 700 pseudosuchian species are known from their 250-million-year fossil record, displaying far greater ecomorphological diversity than their extant counterparts. With a new time-calibrated tree of >500 species, we use a phylogenetic framework to reveal that pseudosuchian evolutionary history and diversification dynamics were directly shaped by the interplay of abiotic and biotic processes over hundreds of millions of years, supported by information theory analyses. Speciation, but not extinction, is correlated with higher temperatures in terrestrial and marine lineages, with high sea level associated with heightened extinction in non-marine taxa. Low lineage diversity and increased speciation in non-marine species is consistent with opportunities for niche-filling, whereas increased competition may have led to elevated extinction rates. In marine lineages, competition via increased lineage diversity appears to have driven both speciation and extinction. Decoupling speciation and extinction, in combination with ecological partitioning, reveals a more complex picture of pseudosuchian evolution than previously understood. As the number of species threatened with extinction by anthropogenic climate change continues to rise, the fossil record provides a unique window into the drivers that led to clade success and those that may ultimately lead to extinction.

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

The evolutionary responses of life-history strategies to climatic variability in flowering plants

The evolution of annual or perennial strategies in flowering plants likely depends on a broad array of temperature and precipitation variables. Previous documented climate life-history correlations in explicit phylogenetic frameworks have been limited to certain clades and geographic regions. To gain insights which generalize to multiple lineages we employ a multi-clade approach analyzing 32 groups of angiosperms across eight climatic variables. We utilize a recently developed method that accounts for the joint evolution of continuous and discrete traits to evaluate two hypotheses: annuals tend to evolve in highly seasonal regions prone to extreme heat and drought; and annuals tend to have faster rates of climatic niche evolution than perennials. We find that temperature, particularly highest temperature of the warmest month, is the most consistent climatic factor influencing the evolution of annual strategy in flowering plants. Unexpectedly, we do not find significant differences in rates of climatic niche evolution between perennial and annual lineages. We propose that annuals are consistently favored in areas prone to extreme heat due to their ability to escape heat stress as seeds, but they tend to be outcompeted by perennials in regions where extreme heat is uncommon or nonexistent.

The evolution of Ericaceae flowers and their pollination syndromes at a global scale

PREMISE

Floral evolution in large clades is difficult to study not only because of the number of species involved, but also because they often are geographically widespread and include a diversity of outcrossing pollination systems. The cosmopolitan blueberry family (Ericaceae) is one such example, most notably pollinated by bees and multiple clades of nectarivorous birds.

METHODS

We combined data on floral traits, pollination ecology, and geography with a comprehensive phylogeny to examine the structuring of floral diversity across pollination systems and continents. We focused on ornithophilous systems to test the hypothesis that some Old World Ericaceae were pollinated by now-extinct hummingbirds.

RESULTS

Despite some support for floral differentiation at a continental scale, we found a large amount of variability within and among landmasses, due to both phylogenetic conservatism and parallel evolution. We found support for floral differentiation in anther and corolla traits across pollination systems, including among different ornithophilous systems. Corolla traits show inconclusive evidence that some Old World Ericaceae were pollinated by hummingbirds, while anther traits show stronger evidence. Some major shifts in floral traits are associated with changes in pollination system, but shifts within bee systems are likely also important.

CONCLUSIONS

Studying the floral evolution of large, morphologically diverse, and widespread clades is feasible. We demonstrate that continent-specific radiations have led to widespread parallel evolution of floral morphology. We show that traits outside of the perianth may hold important clues to the ecological history of lineages.

Profile of a flower: How rates of morphological evolution drive floral diversification in Ericales and angiosperms

PREMISE

Recent studies of floral disparity in the asterid order Ericales have shown that flowers vary strongly among families and that disparity is unequally distributed between the three flower modules (perianth, androecium, gynoecium). However, it remains unknown whether these patterns are driven by heterogeneous rates of morphological evolution or other factors.

METHODS

Here, we compiled a data set of 33 floral characters scored for 414 species of Ericales sampled from 346 genera and all 22 families. We conducted ancestral state reconstructions using an equal-rates Markov model for each character. We estimated rates of morphological evolution for Ericales and for a separate angiosperm-wide data set of 19 characters and 792 species, creating "rate profiles" for Ericales, angiosperms, and major angiosperm subclades. We compared morphological rates among flower modules within each data set separately and between data sets, and we compared rates among angiosperm subclades using the angiosperm data set.

RESULTS

The androecium exhibits the highest evolutionary rates across most characters, whereas most perianth and gynoecium characters evolve more slowly in both Ericales and angiosperms. Both high and low rates of morphological evolution can result in high floral disparity in Ericales. Analyses of an angiosperm-wide floral data set reveal that this pattern appears to be conserved across most major angiosperm clades.

CONCLUSIONS

Elevated rates of morphological evolution in the androecium of Ericales may explain the higher disparity reported for this floral module. Comparing rates of morphological evolution through rate profiles proves to be a powerful tool in understanding floral evolution.

The phylogeny and global biogeography of Primulaceae based on high-throughput DNA sequence data

The angiosperm family Primulaceae is morphologically diverse and distributed nearly worldwide. However, phylogenetic uncertainty has obstructed the identification of major morphological and biogeographic transitions within the clade. We used target capture sequencing with the Angiosperms353 probes, taxon-sampling encompassing nearly all genera of the family, tree-based sequence curation, and multiple phylogenetic approaches to investigate the major clades of Primulaceae and their relationship to other Ericales. We generated dated phylogenetic trees and conducted broad-scale biogeographic analyses as well as stochastic character mapping of growth habit. We show that Ardisia, a pantropical genus and the largest in the family, is not monophyletic, with at least 19 smaller genera nested within it. Neotropical members of Ardisia and several smaller genera form a clade, an ancestor of which arrived in the Neotropics and began diversifying about 20 Ma. This Neotropical clade is most closely related to Elingamita and Tapeinosperma, which are most diverse on islands of the Pacific. Both Androsace and Primula are non-monophyletic by the inclusion of smaller genera. Ancestral state reconstructions revealed that there have either been parallel transitions to an herbaceous habit in Primuloideae, Samolus, and at least three lineages of Myrsinoideae, or a common ancestor of nearly all Primulaceae was herbaceous. Our results provide a robust estimate of phylogenetic relationships across Primulaceae and show that a revised classification of Myrsinoideae and several other clades within the family is necessary to render all genera monophyletic.

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.

The lineage-specific evolution of the oleosin family in Theaceae

Oleosins play essential roles in stabilization of lipid droplets (LDs) and seed oil production. However, evolution of this gene family has not been reported in Theaceae, a large plant family that contains many important tea and oil tea species. In this study, a total of 65 oleosin genes were identified in nine genome-sequenced Theaceae species. Among these genomes, the gene number of oleosin showed significant difference, with Camellia sinensis var. sinensis cv. Shuchazao and Camellia lanceoleosa displayed more oleosin numbers than other species. Phylogenetic analyses revealed that Theaceae oleosin genes were classified into three clades (U, SL, SH) respectively. Proteins within the same clade had similar gene structure and motif composition. Segmental duplication was the primary driving force for the evolution of oleosin genes in Shuchazao (SCZ), Huangdan (HD), C.lanceoleosa (Cla), and wild tea (DASZ). Synteny analysis showed that most oleosin genes displayed inter-species synteny among tea and oil tea species. Expression analysis demonstrated that oleosin genes were specifically expressed in seed and kernel of Huangdan (HD) and C.lanceoleosa. Moreover, expression divergence was observed in paralogous pairs and ∼1-2 oleosin genes in each clade have become activate. This study leads to a comprehensive understanding of evolution of oleosin family in Theaceae, and provides a rich resource to further address the functions of oleosin in tea and oil tea species.

Genomic-based microsatellite development for Ternstroemia (Pentaphylacaceae) and transferability to other Ericales

BACKGROUND

The genus Ternstroemia is associated with the vulnerable tropical montane cloud forest in Mexico and with other relevant vegetation types worldwide. It contains threatened and pharmacologically important species and has taxonomic issues regarding its species limits. This study describes 38 microsatellite markers generated using a genomic-based approach.

METHODS AND RESULTS

We tested 23 of these markers in a natural population of Ternstroemia lineata. These markers are highly polymorphic (all loci polymorphic with 3-14 alleles per locus and expected heterozygosity between 0.202 and 0.908), most of them (19 out of 23) are in Hardy-Weinberg Equilibrium and free of null alleles (18 out of 23). Also we found no evidence of linkage among them. Finally, we tested the transferability to six other American species of Ternstroemia, two other Pentaphylacaceae species, and four species from different families within the order Ericales.

CONCLUSIONS

These molecular resources are promising tools to investigate genetic diversity loss and as barcodes for ethnopharmacological applications and species delimitation in the family Pentaphylacaceae and some Ericales, among other applications.

Knockout of phytoene desaturase gene using CRISPR/Cas9 in highbush blueberry

Among the New Plant Breeding Techniques (NPBTs), the CRISPR/Cas9 system represents a useful tool for target gene editing, improving the traits of the plants rapidly. This technology allows targeting one or more sequences simultaneously, as well as introducing new genetic variations by homology-directed recombination. However, the technology of CRISPR/Cas9 remains a challenge for some polyploid woody species, since all the different alleles for which the mutation is required must be simultaneously targeted. In this work we describe improved protocols adapting the CRISPR/Cas9 system to highbush blueberry (Vaccinium corymbosum L.), using Agrobacterium-mediated transformation. As a proof of concept, we targeted the gene encoding for phytoene desaturase, whose mutation disrupts chlorophyll biosynthesis allowing for the visual assessment of knockout efficiency. Leaf explants of in vitro-cultured blueberry cv. Berkeley has been transformed with a CRISPR/Cas9 construct containing two guide RNAs (gRNA1 and gRNA2) targeting two conserved gene regions of pds and subsequently maintained on a selection medium enriched with kanamycin. After 4 weeks in culture on the selection medium, the kanamycin-resistant lines were isolated, and the genotyping of these lines through Sanger sequencing revealed successful gene editing. Some of mutant shoot lines included albino phenotypes, even if the editing efficiencies were quite low for both gRNAs, ranging between 2.1 and 9.6% for gRNA1 and 3.0 and 23.8 for gRNA2. Here we showed a very effective adventitious shoot regeneration protocol for the commercial cultivar of highbush blueberry "Berkeley", and a further improvement in the use of CRISPR/Cas9 system in Vaccinium corymbosum L., opening the way to the breeding mediated by biotechnological approaches.

Comparative analyses of three complete Primula mitogenomes with insights into mitogenome size variation in Ericales

BACKGROUND

Although knowledge of the sizes, contents, and forms of plant mitochondrial genomes (mitogenomes) is increasing, little is known about the mechanisms underlying their structural diversity. Evolutionary information on the mitogenomes of Primula, an important ornamental taxon, is more limited than the information on their nuclear and plastid counterparts, which has hindered the comprehensive understanding of Primula mitogenomic diversity and evolution. The present study reported and compared three Primula mitogenomes and discussed the size expansion of mitogenomes in Ericales.

RESULTS

Mitogenome master circles were sequenced and successfully assembled for three Primula taxa and were compared with publicly available Ericales mitogenomes. The three mitogenomes contained similar gene contents and varied primarily in their structures. The Primula mitogenomes possessed relatively high nucleotide diversity among all examined plant lineages. In addition, high nucleotide diversity was found among Primula species between the Mediterranean and Himalaya-Hengduan Mountains. Most predicted RNA editing sites appeared in the second amino acid codon, increasing the hydrophobic character of the protein. An early stop in atp6 caused by RNA editing was conserved across all examined Ericales species. The interfamilial relationships within Ericales and interspecific relationships within Primula could be well resolved based on mitochondrial data. Transfer of the two longest mitochondrial plastid sequences (MTPTs) occurred before the divergence of Primula and its close relatives, and multiple independent transfers could also occur in a single MTPT sequence. Foreign sequence [MTPTs and mitochondrial nuclear DNA sequences (NUMTs)] uptake and repeats were to some extent associated with changes in Ericales mitogenome size, although none of these relationships were significant overall.

CONCLUSIONS

The present study revealed relatively conserved gene contents, gene clusters, RNA editing, and MTPTs but considerable structural variation in Primula mitogenomes. Relatively high nucleotide diversity was found in the Primula mitogenomes. In addition, mitogenomic genes, collinear gene clusters, and locally collinear blocks (LCBs) all showed phylogenetic signals. The evolutionary history of MTPTs in Primula was complicated, even in a single MTPT sequence. Various reasons for the size variation observed in Ericales mitogenomes were found.

Potential allopolyploid origin of Ericales revealed with gene-tree reconciliation

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales.

Phylogenomics Resolves the Phylogeny of Theaceae by Using Low-Copy and Multi-Copy Nuclear Gene Makers and Uncovers a Fast Radiation Event Contributing to Tea Plants Diversity

Simple Summary

The Theaceae includes more than 300 species of great morphological diversity and has immense economic, cultural, and ornamental values. However, the evolutionary history of this family remains elusive. We integrated 91 genomes and transcriptome datasets of Theaceae and successfully resolved the phylogeny of Theaceae including relatives of cultivated tea plants from both extensive low-copy and multi-copy nuclear gene markers. Bayes-based molecular dating revealed that the ancestor of the tea family originated slightly earlier than the K-Pg boundary (Mass extinction events including the extinction of dinosaurs) with early diversification of three tribes associated with the Early Eocene Climatic Optimum. Further speciation analysis suggested a sole significant diversification shift rate in the common ancestor of Camellia associated with the Mid-Miocene Climatic Optimum. Collectively, polyploidy events, and key morphological innovation characters, such as pericarp with seed coat hardening, could possibly contribute to the Theaceae species diversity.

Abstract

Tea is one of the three most popular nonalcoholic beverages globally and has extremely high economic and cultural value. Currently, the classification, taxonomy, and evolutionary history of the tea family are largely elusive, including phylogeny, divergence, speciation, and diversity. For understanding the evolutionary history and dynamics of species diversity in Theaceae, a robust phylogenetic framework based on 1785 low-copy and 79,103 multi-copy nuclear genes from 91 tea plant genomes and transcriptome datasets had been reconstructed. Our results maximumly supported that the tribes Stewartieae and Gordonieae are successive sister groups to the tribe Theeae from both coalescent and super matrix ML tree analyses. Moreover, in the most evolved tribe, Theeae, the monophyletic genera Pyrenaria, Apterosperma, and Polyspora are the successive sister groups of Camellia. We also yield a well-resolved relationship of Camellia, which contains the vast majority of Theaceae species richness. Molecular dating suggests that Theaceae originated in the late L-Cretaceous, with subsequent early radiation under the Early Eocene Climatic Optimal (EECO) for the three tribes. A diversification rate shift was detected in the common ancestors of Camellia with subsequent acceleration in speciation rate under the climate optimum in the early Miocene. These results provide a phylogenetic framework and new insights into factors that likely have contributed to the survival of Theaceae, especially a successful radiation event of genus Camellia members to subtropic/tropic regions. These novel findings will facilitate the efficient conservation and utilization of germplasm resources for breeding cultivated tea and oil-tea. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Theaceae.

The symmetry spectrum in a hybridising, tropical group of rhododendrons

Many diverse plant clades possess bilaterally symmetrical flowers and specialised pollination syndromes, suggesting that these traits may promote diversification. We examined the evolution of diverse floral morphologies in a species-rich tropical radiation of Rhododendron. We used restriction-site associated DNA sequencing on 114 taxa from Rhododendron sect. Schistanthe to reconstruct phylogenetic relationships and examine hybridisation. We then captured and quantified floral variation using geometric morphometric analyses, which we interpreted in a phylogenetic context. We uncovered phylogenetic conflict and uncertainty caused by introgression within and between clades. Morphometric analyses revealed flower symmetry to be a morphological continuum without clear transitions between radial and bilateral symmetry. Tropical Rhododendron species that began diversifying into New Guinea c. 6 million years ago expanded into novel floral morphological space. Our results showed that the evolution of tropical Rhododendron is characterised by recent speciation, recurrent hybridisation and the origin of floral novelty. Floral variation evolved via changes to multiple components of the corolla that are only recognised in geometric morphometrics with both front and side views of flowers.

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.

Comparative genomics elucidates the origin of a supergene controlling floral heteromorphism

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?

Plastome Phylogenomics of Aucuba (Garryaceae)

Aucuba (Garryaceae), which includes approximately ten evergreen woody species, is a genus endemic to East Asia. Their striking morphological features give Aucuba species remarkable ornamental value. Owing to high levels of morphological divergence and plasticity, species definitions of Aucuba remain perplexing and problematic. Here, we sequenced and characterized the complete plastid genomes (plastomes) of three Aucuba species: Aucuba chlorascens, Aucuba eriobotryifolia, and Aucuba japonica. Incorporating Aucuba plastomes available in GenBank, a total of seven Aucuba plastomes, representing six out of ten species of Aucuba, were used for comparative plastome analysis, phylogenetic analysis and divergence time estimation in this study. Comparative analyses revealed that plastomes of Aucuba are highly conserved in size, structure, gene content, and organization, and exhibit high levels of sequence similarity. Phylogenetic reconstruction based on 68 plastid protein-coding genes strongly supported the monophyly of Garryales, Garryaceae and Aucuba. Aucuba eriobotryifolia was sister to the other Aucuba species examined, consistent with its unique fused anther locule. The divergence time of Aucuba was estimated to be approximately late Miocene. Extant Aucuba species derived from recent divergence events associated with the establishment of monsoonal climates in East Asia and climatic fluctuations.

Molecular phylogeny of Asian Ardisia (Myrsinoideae, Primulaceae) and their leaf-nodulated endosymbionts, Burkholderia s.l. (Burkholderiaceae)

The genus Ardisia (Myrsinoideae, Primulaceae) has 16 subgenera and over 700 accepted names, mainly distributed in tropical Asia and America. The circumscription of Ardisia is not well-defined and sometimes confounded with the separation of some small genera. A taxonomic revision focusing on Ardisia and allies is necessary. In the Ardisia subgenus Crispardisia, symbiotic association with leaf-nodule bacteria is a unique character within the genus. The endosymbionts are vertically transmitted, highly specific and highly dependent on the hosts, suggesting strict cospeciation may have occurred in the evolutionary history. In the present study, we aimed to establish a phylogenetic framework for further taxonomic revision. We also aimed to test the cospeciation hypothesis of the leaf-nodulate Ardisia and their endosymbiotic bacteria. Nuclear ITS and two chloroplast intergenic spaces were used to reconstruct the phylogeny of Asian Ardisia and relatives in Myrsinoideae, Primulaceae. The 16S-23S rRNA were used to reconstruct the bacterial symbionts’ phylogeny. To understand the evolutionary association of the Ardisia and symbionts, topology tests and cophylogenetic analyses were conducted. The molecular phylogeny suggested Ardisia is not monophyletic, unless Sardiria, Hymenandra, Badula and Oncostemum are included. The results suggest the generic limit within Myrsinoideae (Primulaceae) needs to be further revised. The subgenera Crispardisia, Pimelandra, and Stylardisia were supported as monophyly, while the subgenus Bladhia was separated into two distant clades. We proposed to divide the subgenus Bladhia into subgenus Bladhia s.str. and subgenus Odontophylla. Both of the cophylogenetic analyses and topology tests rejected strict cospeciation hypothesis between Ardisia hosts and symbiotic Burkholderia. Cophylogenetic analyses showed general phylogenetic concordance of Ardisia and Burkholderia, and cospeciation events, host-switching events and loss events were all inferred.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

Exploring the diversity of andean berries from northern Peru based on molecular analyses

More than 12,000 species have been listed under the category of berries, and most of them belong to the orders Ericales and Rosales. Recent phylogenetic studies using molecular data have revealed disagreements with morphological approaches mainly due to diverse floral arrangements, which has proven to be a problem when recognizing species. Therefore, the use of multilocus sequence data is essential to establish robust species boundaries. Although berries are common in Andean cloud forests, diversity of these taxa has not been extensively evaluated in the current context of DNA-based techniques. In this regard, this study characterized morphologically and constructed multilocus phylogenies using four molecular markers, two chloroplast markers (matK and rbcL) and two nuclear markers (ITS and GBSSI-2). Specimens did not show diagnostic features to delimit species of berries. A total of 125 DNA-barcodes of andean berries were newly generated for the four molecular markers. The multilocus phylogenies constructed from these markers allowed the identification of 24 species grouped into the order Ericales (Cavendishia = 1, Clethra = 2, Disterigma = 2, Gaultheria = 4, Thibaudia = 4, Vaccinium = 3) and Rosales (Rubus = 8), incorporating into the Peruvian flora four new records (Disterigma ecuadorense, Disterigma synanthum, Vaccinium meridionale and Rubus glabratus) and revealing the genus Rubus as the most diverse group of berries in the Amazonas region. The results of this study showed congruence in all the multilocus phylogenies, with internal transcribed spacer (ITS) showing the best resolution to distinguish the species. These species were found in coniferous forests, dry and humid forests, rocky slopes, and grasslands at 2,506–3,019 masl from Amazonas region. The integration of morphological and DNA-based methods is recommended to understand the diversity of berries along the Peruvian Andean cloud forest.

Abstract in Quechua language

Qhawarqan astawan chunka iskayniyuq waranqa especiekuna bayasmanta huch’uy mit’a maypichus hatun rak’i chayaqi ordenkunata Ericaleswan Rosaleswan. Chayraqpi Khuski filogeneticamanta rurachiy allincharqan chanikuna molecularkuna willarqan ayñi rikunawanta morfologicokunamanta, qaylla llapan rantichay t’ika tiktutaywan ñawray, ima kay kaqta qhawacgirqan kay huk champay pachaman riqsiypa especiekunamanta. Hina kaqtintaq, chanikuna qatikipaykunamanta multilocus hat’alliy tiksipmi takyachiypaq saywakuna sinchikuna especiekunamanta. Pana bayaskuna kanku allatinkuna sach’a-sach’api phuyusqa anti runap, ñawran manan karqan achka kamaykuy kunan pacha allwiyaraykupi takyasqakuna ADN. Chayrayku, Noqanchispa taqwi allincharqan huk filogenia multilocus, rarachikupúnmi tawa molecular marcadorkuna, caspa iskay markadorkunawan cloroplastomanta (matK, rbcL) iskay markadorkunawan nuclearkunamanta (ITS, GBSSI-2). Kaykunawan filogeniamanta huniqamuran kikinchay iskay chunka tawayoq especies ima tantaqamuran q'anchis generospi (Cavendishia=1, Clethra=2, Disterigma=2, Gaultheria=4, Thibaudia=4, Vaccinium=3, Rubus=8), kaykunata huñuyqamuranta piruwanu llacha kamay tawa musuq quillqakamachikuta (Disterigma ecuadorense, Disterigma synanthum, Vaccinium meridionale, Rubus glabratus). Nocaykuq lluqsisqan kuwirinti rikuchirurqan llapankuna filogeniaspi multilocusmanta, kaspa espaciador transcrito interno (ITS) pi rikuchina kutuwi mihur rantichay riqsiypaq especiekunata.

Abstract in Awajun language

Dekanauwai juú weantug 12000 sag nagkaikiut, júna nejég tente ainawai nuintushkam kuashtai Ericales nuigtu Rosales weantui. Molecularesjai takasmaug juki filogeneticos augtus yamá dekai antugnaiñasmauwa nuna Morfologicosjai disa umikmaug, juka waignawai kuashag yagkunum, juwai dekaata tamanum kuashat utugchata ama nunuka. Nunui asamtai multilocus takasmauwa nujai dekanui wajukut ainawa pipish tumaig aidaush. Tujashkam kuashtai tentee nejég ainaug ikam naig yujagkim amuamua nunuig, wajupá kuashtakit tusajig ashi dekapasjig ADNjain dischamui. Nuni tamaugmak, ii augtusag duka takasé filogenia multilocus dekamua nujai, takasji ipák usumat marcadores molecularesjai, jimag marcadores cloroplastosjai (matK nuigtu rbcL) nuigtu jimag marcadores nuclearesjai (ITS nuigtu GBSSI-2). Juu filogenias dekaji 24 sag nagkaikiut tuwaka 7 generosnug tuwaka awa nunu (Cavendishia=1, Clethra=2, Disterigma=2, Gaultheria=4, Thibaudia=4, Vaccinium=3, Rubus=8), juui dekanai yamajam ipák usumat ajag perunum awanunu (Disterigma ecuadorense, Disterigma synanthum, Vaccinium meridionale nuigtu Rubus glabratus).

Admixture may be extensive among hyperdominant Amazon rainforest tree species

Admixture is a mechanism by which species of long-lived plants may acquire novel alleles. However, the potential role of admixture in the origin and maintenance of tropical plant diversity is unclear. We ask whether admixture occurs in an ecologically important clade of Eschweilera (Parvifolia clade, Lecythidaceae), which includes some of the most widespread and abundant tree species in Amazonian forests. Using target capture sequencing, we conducted a detailed phylogenomic investigation of 33 species in the Parvifolia clade and investigated specific hypotheses of admixture within a robust phylogenetic framework. We found strong evidence of admixture among three ecologically dominant species, E. coriacea, E. wachenheimii and E. parviflora, but a lack of evidence for admixture among other lineages. Accepted species were largely distinguishable from one another, as was geographic structure within species. We show that hybridization may play a role in the evolution of the most widespread and ecologically variable Amazonian tree species. While admixture occurs among some species of Eschweilera, it has not led to widespread erosion of most species' genetic or morphological identities. Therefore, current morphological based species circumscriptions appear to provide a useful characterization of the clade's lineage diversity.

Phylogeographic pattern of a cryptoviviparous mangrove, Aegiceras corniculatum, in the Indo-West Pacific, provides insights for conservation actions

This study identified the historical geoclimatic factors which caused low genetic diversity and strong phylogeographic structure in a cryptoviviparous mangrove. The phylogeographic pattern was used to suggest conservation actions. Phylogeographic studies are used to understand the spatial distribution and evolution of genetic diversity, and have major conservation implications, especially for threatened taxa like the mangroves. This study aimed to assess the phylogeographic pattern of Aegiceras corniculatum, a cryptoviviparous mangrove, across its distribution range in the Indo-West Pacific (IWP) region. We genotyped 398 samples, collected from 37 populations, at four chloroplast DNA (cpDNA) loci, and identified the influence of historical processes on the contemporary population structure of the species. Low genetic diversity at the population level was observed. The evolutionary relationship between 12 cpDNA haplotypes suggested a strong phylogeographic structure, which was further validated by the clustering algorithms and proportioning of maximum variation among hierarchical population groups. The magnitude and direction of historical gene flow indicated that the species attained its wide distribution from its likely ancestral area of the Malay Archipelago. The divergence time estimates of the haplotypes indicated that the geoclimatic changes during the Pleistocene, especially the glacial sea-level changes and emergence of landmasses, hindered genetic exchange and created genetic differentiation between the phylogenetic groups. The species overwintered the last glacial maxima in multiple refugia in the IWP, as identified by the environmental niche modelling. Overall, our findings indicated that ancient glacial vicariance had influenced the present genetic composition of A. corniculatum, which was maintained by the current demographic features of this region. We discussed how these findings can be used to prioritize areas for conservation actions, restore disturbed habitats and prevent further genetic erosion.

Around the world in 40 million years: Phylogeny and biogeography of Tecomeae (Bignoniaceae)

Intercontinental disjunct distributions can arise from vicariance, long distance dispersal, or both. Tecomeae (Bignoniaceae) are a nearly cosmopolitan clade of flowering plants providing us with an excellent opportunity to investigate global distribution patterns. While the tribe contains only about 57 species, it has achieved a distribution that is not only pantropical, but also extends into the temperate zones in both the Northern and Southern hemispheres. This distribution is similar to the distribution of its sister group, a clade of about 750 spp. that includes most remaining taxa in Bignoniaceae. To infer temporal and spatial patterns of dispersal, we generated a phylogeny of Tecomeae by gathering sequence data from chloroplast and nuclear markers for 41 taxa. Fossil calibrations were used to determine divergence times, and ancestral states were reconstructed to infer its biogeographic history. We found support for a South American origin and a crown age of the tribe estimated at ca. 40 Ma. Two dispersal events seem to have happened during the Eocene-Oligocene, one from South America to the Old World, and another from South America to North America. Furthermore, two other dispersal events seem to have taken place during the Miocene, one from North America to Asia, and another from Australia to South America. We suggest that intercontinental dispersal via land bridges and island hopping, as well as sweepstakes of long distance dispersal from the Eocene to the present explain the global distribution of Tecomeae.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

Evolution of Dispersal, Habit, and Pollination in Africa Pushed Apocynaceae Diversification After the Eocene-Oligocene Climate Transition

Apocynaceae (the dogbane and milkweed family) is one of the ten largest flowering plant families, with approximately 5,350 species and diverse morphology and ecology, ranging from large trees and lianas that are emblematic of tropical rainforests, to herbs in temperate grasslands, to succulents in dry, open landscapes, and to vines in a wide variety of habitats. Despite a specialized and conservative basic floral architecture, Apocynaceae are hyperdiverse in flower size, corolla shape, and especially derived floral morphological features. These are mainly associated with the development of corolline and/or staminal coronas and a spectrum of integration of floral structures culminating with the formation of a gynostegium and pollinaria—specialized pollen dispersal units. To date, no detailed analysis has been conducted to estimate the origin and diversification of this lineage in space and time. Here, we use the most comprehensive time-calibrated phylogeny of Apocynaceae, which includes approximately 20% of the species covering all major lineages, and information on species number and distributions obtained from the most up-to-date monograph of the family to investigate the biogeographical history of the lineage and its diversification dynamics. South America, Africa, and Southeast Asia (potentially including Oceania), were recovered as the most likely ancestral area of extant Apocynaceae diversity; this tropical climatic belt in the equatorial region retained the oldest extant lineages and these three tropical regions likely represent museums of the family. Africa was confirmed as the cradle of pollinia-bearing lineages and the main source of Apocynaceae intercontinental dispersals. We detected 12 shifts toward accelerated species diversification, of which 11 were in the APSA clade (apocynoids, Periplocoideae, Secamonoideae, and Asclepiadoideae), eight of these in the pollinia-bearing lineages and six within Asclepiadoideae. Wind-dispersed comose seeds, climbing growth form, and pollinia appeared sequentially within the APSA clade and probably work synergistically in the occupation of drier and cooler habitats. Overall, we hypothesize that temporal patterns in diversification of Apocynaceae was mainly shaped by a sequence of morphological innovations that conferred higher capacity to disperse and establish in seasonal, unstable, and open habitats, which have expanded since the Eocene-Oligocene climate transition.

Gene Duplication and Differential Expression of Flower Symmetry Genes in Rhododendron (Ericaceae)

Bilaterally symmetric flowers have evolved over a hundred times in angiosperms, yet orthologs of the transcription factors CYCLOIDEA (CYC), RADIALIS (RAD), and DIVARICATA (DIV) are repeatedly implicated in floral symmetry changes. We examined these candidate genes to elucidate the genetic underpinnings of floral symmetry changes in florally diverse Rhododendron, reconstructing gene trees and comparing gene expression across floral organs in representative species with radial and bilateral flower symmetries. Radially symmetric R. taxifolium Merr. and bilaterally symmetric R. beyerinckianum Koord. had four and five CYC orthologs, respectively, from shared tandem duplications. CYC orthologs were expressed in the longer dorsal petals and stamens and highly expressed in R. beyerinckianum pistils, whereas they were either ubiquitously expressed, lost from the genome, or weakly expressed in R. taxifolium. Both species had two RAD and DIV orthologs uniformly expressed across all floral organs. Differences in gene structure and expression of Rhododendron RAD compared to other asterids suggest that these genes may not be regulated by CYC orthologs. Our evidence supports CYC orthologs as the primary regulators of differential organ growth in Rhododendron flowers, while also suggesting certain deviations from the typical asterid gene regulatory network for flower symmetry.

Phytogeographic history of the Tea family inferred through high-resolution phylogeny and fossils

The tea family (Theaceae) has a highly unusual amphi-Pacific disjunct distribution: most extant species in the family are restricted to subtropical evergreen broadleaf forests in East Asia, while a handful of species occur exclusively in the subtropical and tropical Americas. Here we used an approach that integrates the rich fossil evidence of this group with phylogenies in biogeographic analysis to study the processes behind this distribution pattern. We first combined genome-skimming sequencing with existing molecular data to build a robust species-level phylogeny for c.140 Theaceae species, resolving most important unclarified relationships. We then developed an empirical Bayesian method to incorporate distribution evidence from fossil specimens into historical biogeographic analyses and used this method to account for the spatiotemporal history of Theaceae fossils. We compared our method with an alternative Bayesian approach and show that it provides consistent results while significantly reduces computational demands which allows analyses of much larger datasets. Our analyses revealed a circumboreal distribution of the family from the early Cenozoic to the Miocene and inferred repeated expansions and retractions of the modelled distribution in the Northern Hemisphere, suggesting that the current Theaceae distribution could be the remnant of a larger continuous distribution associated with the boreotropical forest that has been hypothesized to occupy most of the northern latitudes in the early Cenozoic. These results contradict with studies that only considered current species distributions and showcase the necessity of integrating fossil and molecular data in phylogeny-based parametric biogeographic models to improve the reliability of inferred biogeographical events.

Complete chloroplast genome of novel Adrinandra megaphylla Hu species: molecular structure, comparative and phylogenetic analysis

Adrinandra megaphylla Hu is a medicinal plant belonging to the Adrinandra genus, which is well-known for its potential health benefits due to its bioactive compounds. This study aimed to assemble and annotate the chloroplast genome of A. megaphylla as well as compare it with previously published cp genomes within the Adrinandra genus. The chloroplast genome was reconstructed using de novo and reference-based assembly of paired-end reads generated by long-read sequencing of total genomic DNA. The size of the chloroplast genome was 156,298 bp, comprised a large single-copy (LSC) region of 85,688 bp, a small single-copy (SSC) region of 18,424 bp, and a pair of inverted repeats (IRa and IRb) of 26,093 bp each; and a total of 51 SSRs and 48 repeat structures were detected. The chloroplast genome includes a total of 131 functional genes, containing 86 protein-coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes. The A. megaphylla chloroplast genome indicated that gene content and structure are highly conserved. The phylogenetic reconstruction using complete cp sequences, matK and trnL genes from Pentaphylacaceae species exhibited a genetic relationship. Among them, matK sequence is a better candidate for phylogenetic resolution. This study is the first report for the chloroplast genome of the A. megaphylla.

Biodiversity at the global scale: the synthesis continues

Traditionally, the generation and use of biodiversity data and their associated specimen objects have been primarily the purview of individuals and small research groups. While deposition of data and specimens in herbaria and other repositories has long been the norm, throughout most of their history, these resources have been accessible only to a small community of specialists. Through recent concerted efforts, primarily at the level of national and international governmental agencies over the last two decades, the pace of biodiversity data accumulation has accelerated, and a wider array of biodiversity scientists has gained access to this massive accumulation of resources, applying them to an ever-widening compass of research pursuits. We review how these new resources and increasing access to them are affecting the landscape of biodiversity research in plants today, focusing on new applications across evolution, ecology, and other fields that have been enabled specifically by the availability of these data and the global scope that was previously beyond the reach of individual investigators. We give an overview of recent advances organized along three lines: broad-scale analyses of distributional data and spatial information, phylogenetic research circumscribing large clades with comprehensive taxon sampling, and data sets derived from improved accessibility of biodiversity literature. We also review synergies between large data resources and more traditional data collection paradigms, describe shortfalls and how to overcome them, and reflect on the future of plant biodiversity analyses in light of increasing linkages between data types and scientists in our field.

Phylogenomic conflict coincides with rapid morphological innovation

Evolutionary biologists have long been fascinated with the episodes of rapid phenotypic innovation that underlie the emergence of major lineages. Although our understanding of the environmental and ecological contexts of such episodes has steadily increased, it has remained unclear how population processes contribute to emergent macroevolutionary patterns. One insight gleaned from phylogenomics is that gene-tree conflict, frequently caused by population-level processes, is often rampant during the origin of major lineages. With the understanding that phylogenomic conflict is often driven by complex population processes, we hypothesized that there may be a direct correspondence between instances of high conflict and elevated rates of phenotypic innovation if both patterns result from the same processes. We evaluated this hypothesis in six clades spanning vertebrates and plants. We found that the most conflict-rich regions of these six clades also tended to experience the highest rates of phenotypic innovation, suggesting that population processes shaping both phenotypic and genomic evolution may leave signatures at deep timescales. Closer examination of the biological significance of phylogenomic conflict may yield improved connections between micro- and macroevolution and increase our understanding of the processes that shape the origin of major lineages across the Tree of Life.

Striking developmental convergence in angiosperm endoparasites

PREMISE

A subset of parasitic plants bear extremely reduced features and grow nearly entirely within their hosts. Until recently, most of these endoparasites were thought to represent a single clade united by their reduced morphology. Current phylogenetic understanding contradicts this assumption and indicates these plants represent distantly related clades, thus offering an opportunity to examine convergence among plants with this life history.

METHODS

We sampled species from Apodanthaceae, Cytinaceae, Mitrastemonaceae, and Rafflesiaceae spanning a range of developmental stages. To provide a broader comparative framework, Santalaceae mistletoes with a similar lifestyle were also analyzed. Microtomography and microscopy were used to analyze growth patterns and the ontogeny of host-parasite vascular connections.

RESULTS

Apodanthaceae, Cytinaceae, Mitrastemonaceae, and Rafflesiaceae species demonstrated a common development characterized by late cell differentiation. These species were also observed to form direct connections to host vessels and to cause severe alterations of host xylem development. Apodanthaceae and Rafflesiaceae species were additionally observed to form sieve elements, which connect with the host phloem. Endophytic Santalaceae species demonstrated a dramatically different developmental pattern, featuring early cell differentiation and tissue organization, and little effect on host anatomy and cambial activity.

CONCLUSIONS

Our results illuminate two distinct developmental trajectories in endoparasites. One involves the retention of embryonic characteristics and late connection with host vessels, as demonstrated in species of Apodanthaceae, Cytinaceae, Mitrastemonaceae, and Rafflesiaceae. The second involves tissue specialization and early connection with host xylem, as exemplified by Santalaceae species. These differences are hypothesized to be related to the absence/presence of photosynthesis in these plants.

Asterid Phylogenomics/Phylotranscriptomics Uncover Morphological Evolutionary Histories and Support Phylogenetic Placement for Numerous Whole-Genome Duplications

Asterids are one of the most successful angiosperm lineages, exhibiting extensive morphological diversity and including a number of important crops. Despite their biological prominence and value to humans, the deep asterid phylogeny has not been fully resolved, and the evolutionary landscape underlying their radiation remains unknown. To resolve the asterid phylogeny, we sequenced 213 transcriptomes/genomes and combined them with other data sets, representing all accepted orders and nearly all families of asterids. We show fully supported monophyly of asterids, Berberidopsidales as sister to asterids, monophyly of all orders except Icacinales, Aquifoliales, and Bruniales, and monophyly of all families except Icacinaceae and Ehretiaceae. Novel taxon placements benefited from the expanded sampling with living collections from botanical gardens, resolving hitherto uncertain relationships. The remaining ambiguous placements here are likely due to limited sampling and could be addressed in the future with relevant additional taxa. Using our well-resolved phylogeny as reference, divergence time estimates support an Aptian (Early Cretaceous) origin of asterids and the origin of all orders before the Cretaceous-Paleogene boundary. Ancestral state reconstruction at the family level suggests that the asterid ancestor was a woody terrestrial plant with simple leaves, bisexual, and actinomorphic flowers with free petals and free anthers, a superior ovary with a style, and drupaceous fruits. Whole-genome duplication (WGD) analyses provide strong evidence for 33 WGDs in asterids and one in Berberidopsidales, including four suprafamilial and seven familial/subfamilial WGDs. Our results advance the understanding of asterid phylogeny and provide numerous novel evolutionary insights into their diversification and morphological evolution.

Global patterns and a latitudinal gradient of flower disparity: perspectives from the angiosperm order Ericales

Morphological diversity (disparity) is an essential but often neglected aspect of biodiversity. Hence, it seems timely and promising to re-emphasize morphology in modern evolutionary studies. Disparity is a good proxy for the diversity of functions and interactions with the environment of a group of taxa. In addition, geographical and ecological patterns of disparity are crucial to understand organismal evolution and to guide biodiversity conservation efforts. Here, we analyse floral disparity across latitudinal intervals, growth forms, climate types, types of habitats, and regions for a large and representative sample of the angiosperm order Ericales. We find a latitudinal gradient of floral disparity and a decoupling of disparity from species richness. Other factors investigated are intercorrelated, and we find the highest disparity for tropical trees growing in African and South American forests. Explanations for the latitudinal gradient of floral disparity may involve the release of abiotic constraints and the increase of biotic interactions towards tropical latitudes, allowing tropical lineages to explore a broader area of the floral morphospace. Our study confirms the relevance of biodiversity parameters other than species richness and is consistent with the importance of species interactions in the tropics, in particular with respect to angiosperm flowers and their pollinators.

The complete mitochondrial and plastid genomes of Rhododendron simsii, an important parent of widely cultivated azaleas

The genus Rhododendron of the heather family (Ericaceae) is well known and widely cultivated for their highly ornamental value. The most widely cultivated Rhododendron species is Rhododendron simsii (Indoor azalea). In this study, we assembled the complete linear mitochondrial genome (GenBank accession number MW030508) and quadripartite plastid genome (GenBank accession number MW030509). The mitochondrial genome is 802,707 bp in length with containing 53 unique genes (33 protein-coding, 17 tRNA, and 3 rRNA genes), while the 152,214 bp long plastid genome is smaller and containing 105 unique genes (4 rRNA, 26 tRNA, and 75 protein-coding genes). Phylogenetic analysis showed that the same species relationship with APG system as well as the low supports of branches which is the common characteristic of resolved Ericales phylogenetics.

New genetic markers for Sapotaceae phylogenomics: More than 600 nuclear genes applicable from family to population levels

Some tropical plant families, such as the Sapotaceae, have a complex taxonomy, which can be resolved using Next Generation Sequencing (NGS). For most groups however, methodological protocols are still missing. Here we identified 531 monocopy genes and 227 Short Tandem Repeats (STR) markers and tested them on Sapotaceae using target capture and NGS. The probes were designed using two genome skimming samples from Capurodendron delphinense and Bemangidia lowryi, both from the Tseboneae tribe, as well as the published Manilkara zapota transcriptome from the Sapotoideae tribe. We combined our probes with 261 additional ones previously published and designed for the entire angiosperm group. On a total of 792 low-copy genes, 638 showed no signs of paralogy and were used to build a phylogeny of the family with 231 individuals from all main lineages. A highly supported topology was obtained at high taxonomic ranks but also at the species level. This phylogeny revealed the existence of more than 20 putative new species. Single nucleotide polymorphisms (SNPs) extracted from the 638 genes were able to distinguish lineages within a species complex and to highlight geographical structuration. STR were recovered efficiently for the species used as reference (C. delphinense) but the recovery rate decreased dramatically with the phylogenetic distance to the focal species. Altogether, the new loci will help reaching a sound taxonomic understanding of the family Sapotaceae for which many circumscriptions and relationships are still debated, at the species, genus and tribe levels.

Angiosperms at the edge: Extremity, diversity, and phylogeny

A hallmark of flowering plants is their ability to invade some of the most extreme and dynamic habitats, including cold and dry biomes, to a far greater extent than other land plants. Recent work has provided insight to the phylogenetic distribution and evolutionary mechanisms which have enabled this success, yet needed is a synthesis of evolutionary perspectives with plant physiological traits, morphology, and genomic diversity. Linking these disparate components will not only lead to better understand the evolutionary parallelism and diversification of plants with these two strategies, but also to provide the framework needed for directing future research. We summarize the primary physiological and structural traits involved in response to cold- and drought stress, outline the phylogenetic distribution of these adaptations, and describe the recurring association of these changes with rapid diversification events that occurred in multiple lineages over the past 15 million years. Across these threefold facets of dry-cold correlation (traits, phylogeny, and time) we stress the contrast between (a) the amazing diversity of solutions flowering plants have developed in the face of extreme environments and (b) a broad correlation between cold and dry adaptations that in some cases may hint at deep common origins.

Complete Chloroplast Genome of Argania spinosa: Structural Organization and Phylogenetic Relationships in Sapotaceae

Argania spinosa (Sapotaceae), an important endemic Moroccan oil tree, is a primary source of argan oil, which has numerous dietary and medicinal proprieties. The plant species occupies the mid-western part of Morocco and provides great environmental and socioeconomic benefits. The complete chloroplast (cp) genome of A. spinosa was sequenced, assembled, and analyzed in comparison with those of two Sapotaceae members. The A. spinosa cp genome is 158,848 bp long, with an average GC content of 36.8%. The cp genome exhibits a typical quadripartite and circular structure consisting of a pair of inverted regions (IR) of 25,945 bp in length separating small single-copy (SSC) and large single-copy (LSC) regions of 18,591 and 88,367 bp, respectively. The annotation of A. spinosa cp genome predicted 130 genes, including 85 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. A total of 44 long repeats and 88 simple sequence repeats (SSR) divided into mononucleotides (76), dinucleotides (7), trinucleotides (3), tetranucleotides (1), and hexanucleotides (1) were identified in the A. spinosa cp genome. Phylogenetic analyses using the maximum likelihood (ML) method were performed based on 69 protein-coding genes from 11 species of Ericales. The results confirmed the close position of A. spinosa to the Sideroxylon genus, supporting the revisiting of its taxonomic status. The complete chloroplast genome sequence will be valuable for further studies on the conservation and breeding of this medicinally and culinary important species and also contribute to clarifying the phylogenetic position of the species within Sapotaceae.

Phylogenetic analysis of fossil flowers using an angiosperm-wide data set: proof-of-concept and challenges ahead

PREMISE

Significant paleobotanical discoveries in recent decades have considerably improved our understanding of the early evolution of angiosperms and their flowers. However, our ability to test the systematic placement of fossil flowers on the basis of phylogenetic analyses has remained limited, mainly due to the lack of an adequate, angiosperm-wide morphological data set for extant taxa. Earlier attempts to place fossil flowers phylogenetically were, therefore, forced to make prior qualitative assessments of the potential systematic position of fossils and to restrict phylogenetic analyses to selected angiosperm subgroups.

METHODS

We conduct angiosperm-wide molecular backbone analyses of 10 fossil flower taxa selected from the Cretaceous record. Our analyses make use of a floral trait data set built within the framework of the eFLOWER initiative. We provide an updated version of this data set containing data for 28 floral and two pollen traits for 792 extant species representing 372 angiosperm families.

RESULTS

We find that some fossils are placed congruently with earlier hypotheses while others are found in positions that had not been suggested previously. A few take up equivocal positions, including the stem branches of large clades.

CONCLUSIONS

Our study provides an objective approach to test for the phylogenetic position of fossil flowers across angiosperms. Such analyses may provide a complementary tool for paleobotanical studies, allowing for a more comprehensive understanding of fossil phylogenetic relationships in angiosperms. Ongoing work focused on extending the sampling of extant taxa and the number of floral traits will further improve the applicability and accuracy of our approach.

Large-scale survey reveals pervasiveness and potential function of endogenous geminiviral sequences in plants

The family Geminiviridae contains viruses with single-stranded DNA genomes that have been found infecting a wide variety of angiosperm species. The discovery within the last 25 years of endogenous geminivirus-like (EGV) elements within the nuclear genomes of several angiosperms has raised questions relating to the pervasiveness of EGVs and their impacts on host biology. Only a few EGVs have currently been characterized and it remains unclear whether any of these have influenced, or are currently influencing, the evolutionary fitness of their hosts. We therefore undertook a large-scale search for evidence of EGVs within 134 genome and 797 transcriptome sequences of green plant species. We detected homologues of geminivirus replication-associated protein (Rep) genes in forty-two angiosperm species, including two monocots, thirty-nine dicots, and one ANITA-grade basal angiosperm species (Amborella trichopoda). While EGVs were present in the members of many different plant orders, they were particularly common within the large and diverse order, Ericales, with the highest copy numbers of EGVs being found in two varieties of tea plant (Camellia sinensis). Phylogenetic and clustering analyses revealed multiple highly divergent previously unknown geminivirus Rep lineages, two of which occur in C.sinensis alone. We find that some of the Camellia EGVs are likely transcriptionally active, sometimes co-transcribed with the same host genes across several Camellia species. Overall, our analyses expand the known breadths of both geminivirus diversity and geminivirus host ranges, and strengthens support for the hypothesis that EGVs impact the biology of their hosts.

Phylogeny and biogeography of South American marsh pitcher plant genus Heliamphora (Sarraceniaceae) endemic to the Guiana Highlands

Heliamphora is a genus of carnivorous pitcher plants endemic to the Guiana Highlands with fragmented distributions. We present a well resolved, time-calibrated, and comprehensive Heliamphora phylogeny estimated using Bayesian inference and maximum likelihood based on nuclear genes (26S, ITS, and PHYC) and secondary calibration. We used stochastic mapping to infer ancestral states of morphological characters and ecological traits. Our ancestral state estimations revealed that the pitcher drainage structures characteristic of the genus transformed from a hole to a slit in single clade, while other features (scape pubescence and hammock-like growth) have been gained and lost multiple times. Habitat was similarly labile in Heliamphora, with multiple transitions from the ancestral highland habitats into the lowlands. Using a Mantel test, we found closely related species tend to be geographically closely distributed. Placing our phylogeny in a historical context, major clades likely emerged through both vicariance and dispersal during the Miocene with more recent diversification driven by vertical displacement during the Pleistocene glacial-interglacial thermal oscillations. Despite the dynamic climatic history experienced by Heliamphora, the temperature changes brought by global warming pose a significant threat, particularly for those species at the highest elevations.

Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal

Phylogenomic data from a rapidly increasing number of studies provide new evidence for resolving relationships in recently radiated clades, but they also pose new challenges for inferring evolutionary histories. Most existing methods for reconstructing phylogenetic hypotheses rely solely on algorithms that only consider incomplete lineage sorting (ILS) as a cause of intra- or intergenomic discordance. Here, we utilize a variety of methods, including those to infer phylogenetic networks, to account for both ILS and introgression as a cause for nuclear and cytoplasmic-nuclear discordance using phylogenomic data from the recently radiated flowering plant genus Polemonium (Polemoniaceae), an ecologically diverse genus in Western North America with known and suspected gene flow between species. We find evidence for widespread discordance among nuclear loci that can be explained by both ILS and reticulate evolution in the evolutionary history of Polemonium. Furthermore, the histories of organellar genomes show strong discordance with the inferred species tree from the nuclear genome. Discordance between the nuclear and plastid genome is not completely explained by ILS, and only one case of discordance is explained by detected introgression events. Our results suggest that multiple processes have been involved in the evolutionary history of Polemonium and that the plastid genome does not accurately reflect species relationships. We discuss several potential causes for this cytoplasmic-nuclear discordance, which emerging evidence suggests is more widespread across the Tree of Life than previously thought. [Cyto-nuclear discordance, genomic discordance, phylogenetic networks, plastid capture, Polemoniaceae, Polemonium, reticulations.].

Nuclear phylogenomic analyses of asterids conflict with plastome trees and support novel relationships among major lineages

PREMISE

Discordance between nuclear and organellar phylogenies (cytonuclear discordance) is a well-documented phenomenon at shallow evolutionary levels but has been poorly investigated at deep levels of plant phylogeny. Determining the extent of cytonuclear discordance across major plant lineages is essential not only for elucidating evolutionary processes, but also for evaluating the currently used framework of plant phylogeny, which is largely based on the plastid genome.

METHODS

We present a phylogenomic examination of a major angiosperm clade (Asteridae) based on sequence data from the nuclear, plastid, and mitochondrial genomes as a means of evaluating currently accepted relationships inferred from the plastome and exploring potential sources of genomic conflict in this group.

RESULTS

We recovered at least five instances of well-supported cytonuclear discordance concerning the placements of major asterid lineages (i.e., Ericales, Oncothecaceae, Aquifoliales, Cassinopsis, and Icacinaceae). We attribute this conflict to a combination of incomplete lineage sorting and hybridization, the latter supported in part by previously inferred whole-genome duplications.

CONCLUSIONS

Our results challenge several long-standing hypotheses of asterid relationships and have implications for morphological character evolution and for the importance of ancient whole-genome duplications in early asterid evolution. These findings also highlight the value of reevaluating broad-scale angiosperm and green-plant phylogeny with nuclear genomic data.

A consensus phylogenomic approach highlights paleopolyploid and rapid radiation in the history of Ericales

PREMISE

Large genomic data sets offer the promise of resolving historically recalcitrant species relationships. However, different methodologies can yield conflicting results, especially when clades have experienced ancient, rapid diversification. Here, we analyzed the ancient radiation of Ericales and explored sources of uncertainty related to species tree inference, conflicting gene tree signal, and the inferred placement of gene and genome duplications.

METHODS

We used a hierarchical clustering approach, with tree-based homology and orthology detection, to generate six filtered phylogenomic matrices consisting of data from 97 transcriptomes and genomes. Support for species relationships was inferred from multiple lines of evidence including shared gene duplications, gene tree conflict, gene-wise edge-based analyses, concatenation, and coalescent-based methods, and is summarized in a consensus framework.

RESULTS

Our consensus approach supported a topology largely concordant with previous studies, but suggests that the data are not capable of resolving several ancient relationships because of lack of informative characters, sensitivity to methodology, and extensive gene tree conflict correlated with paleopolyploidy. We found evidence of a whole-genome duplication before the radiation of all or most ericalean families, and demonstrate that tree topology and heterogeneous evolutionary rates affect the inferred placement of genome duplications.

CONCLUSIONS

We provide several hypotheses regarding the history of Ericales, and confidently resolve most nodes, but demonstrate that a series of ancient divergences are unresolvable with these data. Whether paleopolyploidy is a major source of the observed phylogenetic conflict warrants further investigation.