Phylotranscriptomic analyses reveal deep gene tree discordance in Camellia (Theaceae).

Gene tree discordance is a significant legacy of biological evolution. Multiple factors can result in incongruence among genes, such as introgression, incomplete lineage sorting (ILS), gene duplication or loss. Resolving the background of gene tree discordance is a critical way to uncover the process of species diversification. Camellia, the largest genus in Theaceae, has controversial taxonomy and systematics due in part to a complex evolutionary history. We used 60 transcriptomes of 55 species, which represented 15 sections of Camellia to investigate its phylogeny and the possible causes of gene tree discordance. We conducted gene tree discordance analysis based on 1,617 orthologous low-copy nuclear genes, primarily using coalescent species trees and polytomy tests to distinguish hard and soft conflict. A selective pressure analysis was also performed to assess the impact of selection on phylogenetic topology reconstruction. Our results detected different levels of gene tree discordance in the backbone of Camellia, and recovered rapid diversification as one of the possible causes of gene tree discordance. Furthermore, we confirmed that none of the currently proposed sections of Camellia was monophyletic. Comparisons among datasets partitioned under different selective pressure regimes showed that integrating all orthologous genes provided the best phylogenetic resolution of the species tree of Camellia. The findings of this study reveal rapid diversification as a major source of gene tree discordance in Camellia and will facilitate future investigation of reticulate relationships at the species level in this important plant genus.

Phylogenomic insights into the origin and evolutionary history of evergreen broadleaved forests in East Asia under Cenozoic climate change.

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

Species discrimination in Schima (Theaceae): Next-generation super-barcodes meet evolutionary complexity.

Plastid genome and nuclear ribosomal DNA (nrDNA) arrays, proposed recently as "super-barcodes," might provide additional discriminatory power and overcome the limitations of traditional barcoding loci, yet super-barcodes need to be tested for their effectiveness in more plant groups. Morphological homoplasy among Schima species makes the genus a model for testing the efficacy of super-barcodes. In this study, we generated multiple data sets comprising standard DNA barcodes (matK, rbcL, trnH-psbA, nrITS) and super-barcodes (plastid genome, nrDNA arrays) across 58 individuals from 12 out of 13 species of Schima from China. No samples were correctly assigned to species using standard DNA barcodes and nrDNA arrays, while only 27.27% of species with multiple accessions were distinguished using the plastid genome and its partitioned data sets-the lowest estimated rate of super-barcode success in the literature so far. For Schima and other taxa with similarly recently divergence and low levels of genetic variation, incomplete lineage sorting, hybridization or taxonomic oversplitting are all possible causes of the failure. Taken together, our study suggests that by no means are super-barcodes immune to the challenges imposed by evolutionary complexity. We therefore call for developing multilocus nuclear markers for species discrimination in plant groups.

Opening a door to the spatiotemporal history of plants from the tropical Indochina Peninsula to subtropical China.

The complexity of global biodiversity in the tropical Indochina Peninsula and subtropical China bioregions has fascinated biologists for decades, but little is known about the spatiotemporal patterns in these regions. Accordingly, the aims of present study were to investigate the evolutionary and distribution patterns of Engelhardia in these regions and establish a model for examining biogeographic patterns and geological events throughout the tropical Indochina Peninsula and subtropical China. The effects of geological events occurring in the area between the Indochina Peninsula and subtropical China bioregions on the two trees species (i.e., E. roxburghiana and E. fenzelii) were evaluated. A robust phylogenetic framework of 884 individuals from 79 populations was used to generate time-calibrated cytoplasmic and nuclear phylogenetic frameworks based on cpDNA, nrDNA, and nSSR data, respectively. When considered along with ancestral area reconstructions, the genetic data were also used to assess and reconstruct the species' population genetic structure and diversity. These analyses yielded important information about the (1) historical distribution relationships between the tropical and subtropical flora of China; (2) effects of the East Asian summer monsoon (EASM) on the evolutionary history of Asia's plants; and (3) importance of biogeography in conservation planning. Although cytoplasmic-nuclear discordance indicated cpDNA and nrDNA were subject to distinct evolutionary mechanisms that reflected respective evolutionary histories of the plastid and nuclear genomes of prior demographic and biogeographic events. The tropical elements of Engelhardia occupied the Indochina Peninsula during the early Eocene, whereas the subtropical elements were transformed from the tropical elements during Miocene cooling and the onset of the EASM at the Oligocene-Miocene boundary, intensified during the late Miocene and Pliocene, facilitating the transformation of Engelhardia from the tropical Indochina Peninsula to subtropical China. Demographic history provided insights into prominent planning frameworks in conservation biology, namely that subtropical China functioned as a refugium during past climate oscillations and will continue to serve in this capacity in the future.

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.

Speciation along the elevation gradient: Divergence of Roscoea species within the south slope of the Himalayas.

The Himalayas with dramatic elevation gradient is one of the global biodiversity hotspots. Although origin of biodiversity of the Himalayas is of great concern, the speciation process within the Himalayas is poorly known. Roscoea within the Himalayas serve as a good model system to test the speciation process along an elevation gradient. 32,375 unlinked SNPs were used to reconstruct phylogenetic relationships and introgression analyses in D-statistics and Fastsimicoal2. Species distribution modeling (SDM) was used to simulate habitat shift of Roscoea species during climate changes. Phylogeny suggested that the speciation order, except R. capitata, was from highland to lowland. D-statistics analyses suggested significant bidirectional ancient introgression between elevation-neighboring clades but no introgression between R. capitata and othern clades and no introgression among extant species. Fastsimicoal2 suggested interspecific introgressions were asymmetric. SDM predicted that habitats of Roscoea shifted to low elevation during cooling age. These results suggested that the sudden uplift of the Himalayas likely promoted speciation by vicariance, and climate cooling drove species divergence towards lower elevation. This study provides explanations for the origin of biodiversity within the Himalayas, and an insight to understand speciation along elevation in the mountainous regions.

The complete chloroplast genome sequence of Laplacea alpestris and its phylogenetic position.

Laplacea alpestris is a member of the genus Laplacea, which distributes in Central and South America. Genetic information of L. alpestris would provide guidance for the phylogenetic position of this species. Here, we reported and characterized its complete chloroplast (cp) genome using Illumina pair-end sequencing data. The total chloroplast genome size of this species was 157,211 bp, including inverted repeats (IRs) of 26,103 bp, separated by a large single copy (LSC) and a small single copy (SSC) of 86,749 and 18,256 bp, respectively. A total of 132 genes, including 37 tRNA, 8 rRNA, and 87 protein-coding genes were identified. Phylogenetic analysis showed that L. alpestris formed a monophyletic clade with Laplacea fruticosa, and then grouped with Apterosperma oblata. The systematic position of Southeast Asian Laplacea species needs further studies.

Insights into the historical assembly of East Asian subtropical evergreen broadleaved forests revealed by the temporal history of the tea family.

Subtropical evergreen broadleaved forests (EBLFs) inhabit large areas of East Asia. Although paleovegetation reconstructions have revealed that the subtropical EBLFs existed in Southwest China during the Miocene, the historical construction of these forests remains poorly known. Here, we used the tea family (Theaceae), a characteristic component of the subtropical EBLFs, to gain new insights into the assembly of this important biome. Using a robust phylogenetic framework of Theaceae based on plastome and nuclear ribosomal DNA sequence data, the temporal history of the family was reconstructed. Data from other characteristic components of subtropical EBLFs, including Fagaceae, Lauraceae and Magnoliaceae, were also integrated. Most of the essential elements of the subtropical EBLFs appear to have originated around the Oligocene-Miocene (O-M) boundary. However, small woody lineages (e.g. Camellia, Hartia) from Theaceae were dated to the late Miocene. Accelerated net diversification rates within Theaceae were also detected near the O-M transition period and the late Miocene. Our results suggest that two independent intensifications of the East Asian summer monsoon (EASM) around the O-M boundary and the late Miocene may have facilitated the historical assembly of the subtropical EBLFs in East Asia.

Comparative chloroplast genomes of eleven Schima (Theaceae) species: Insights into DNA barcoding and phylogeny.

Schima is an ecologically and economically important woody genus in tea family (Theaceae). Unresolved species delimitations and phylogenetic relationships within Schima limit our understanding of the genus and hinder utilization of the genus for economic purposes. In the present study, we conducted comparative analysis among the complete chloroplast (cp) genomes of 11 Schima species. Our results indicate that Schima cp genomes possess a typical quadripartite structure, with conserved genomic structure and gene order. The size of the Schima cp genome is about 157 kilo base pairs (kb). They consistently encode 114 unique genes, including 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, with 17 duplicated in the inverted repeat (IR). These cp genomes are highly conserved and do not show obvious expansion or contraction of the IR region. The percent variability of the 68 coding and 93 noncoding (>150 bp) fragments is consistently less than 3%. The seven most widely touted DNA barcode regions as well as one promising barcode candidate showed low sequence divergence. Eight mutational hotspots were identified from the 11 cp genomes. These hotspots may potentially be useful as specific DNA barcodes for species identification of Schima. The 58 cpSSR loci reported here are complementary to the microsatellite markers identified from the nuclear genome, and will be leveraged for further population-level studies. Phylogenetic relationships among the 11 Schima species were resolved with strong support based on the cp genome data set, which corresponds well with the species distribution pattern. The data presented here will serve as a foundation to facilitate species identification, DNA barcoding and phylogenetic reconstructions for future exploration of Schima.

Phylogeny and historical biogeography of Isodon (Lamiaceae): rapid radiation in south-west China and Miocene overland dispersal into Africa.

Rapid organismal radiations occurring on the Qinghai-Tibetan Plateau (QTP) and the mechanisms underlying Asia-Africa intercontinental disjunctions have both attracted much attention from evolutionary biologists. Here we use the genus Isodon (Lamiaceae), a primarily East Asian lineage with disjunct species in central and southern Africa, as a case study to shed light upon these processes. The molecular phylogeny and biogeographic history of Isodon were reconstructed using sequences of three plastid markers, the nuclear ribosomal internal transcribed spacer (nrITS), and a low-copy nuclear gene (LEAFY intron II). The evolution of chromosome numbers in this genus was also investigated using probabilistic models. Our results support a monophyletic Isodon that includes the two disjunct African species, both of which likely formed through allopolyploidy. An overland migration from Asia to Africa through Arabia during the early Miocene is proposed as the most likely explanation for the present disjunct distribution of Isodon. The opening of the Red Sea in the middle Miocene may appear to have had a major role in disrupting floristic exchange between Asia and Africa. In addition, a rapid radiation of Isodon was suggested to occur in the late Miocene. It corresponds with one of the major uplifts of the QTP and subsequent aridification events. Our results support the hypothesis that geological and climatic events play important roles in driving biological diversification of organisms distributed in the QTP area.

Isolation and characterization of microsatellite markers for a worldwide invasive weed, Chromolaena odorata (Asteraceae).

PREMISE OF THE STUDY: Microsatellite loci were isolated and characterized for a worldwide invasive weed Chromolaena odorata (Asteraceae) to elucidate the population genetic structure and invasive history. • METHODS AND RESULTS: A total of 14 microsatellite primer pairs were developed using the Fast Isolation by AFLP of Sequences Containing repeats (FIASCO) protocol, and their polymorphism was assessed in two natural populations of C. odorata from Mexico and Trinidad and Tobago. Eleven loci showed polymorphism and eight of these loci were successfully amplified in Ageratina adenophora, another invasive weed related to C. odorata. • CONCLUSIONS: These microsatellite markers are useful for investigating the population genetic structure and the history of range expansion of these invasive species.