Conflicting phylogenetic signals in plastomes of the tribe Laureae (Lauraceae)

Genomic insights and phylogenetics of the mitochondrial genome of Cryptocarya

The tropical genus Cryptocarya is known for its valuable timber and the constituents in these trees show potential for medicinal properties. However, the phylogenetic relationships among species in Asia remain unclear. Here, we report the first mitochondrial genome for Cryptocarya kwangtungensis, consisting of 758,020 bp, including 43 protein-coding genes, 23 tRNA genes and three rRNA genes, with 234 simple sequence repeats, and 1275 dispersed repeats, 35 homologous DNA fragments between the mitogenome and the plastome. Comparative genomic analysis indicated frequent recombination events among the sequences of five magnoliids mitogenomes and only five conserved clustered genes. Further phylogenetic analyses based on 91 mitochondrial regions and nuclear ribosomal cistron sequences of 21 species compound three well-resolved congruent groups for the Cryptocarya species in Asia, both of which support the genus divide into three clades, suggesting that the mitogenome sequence can provide strongly supported relationships within the genus in the family Lauraceae.

Sixteen New Complete Plastid Genomes in the Tribe Loteae (Leguminosae): Structure and Phylogenetic Analysis

The tribe Loteae (Papilioniodeae-Leguminosae), according to plastid data, belongs to the Robinioid clade, which also includes the tribes Robinieae and Sesbanieae. The tribe Loteae contains 16 genera and about two hundred seventy-five species, of which the plastid genomes of five species have been studied to date. The main objectives of our study were to obtain new information on the plastid genome structure of the Loteae representatives in order to assess plastid genome variability and reconstruct phylogenetic relationships within the tribe Loteae. We performed sequencing, assembly, structural and phylogenetic analyses of the Loteae plastid genomes. All assembled Loteae plastomes showed a quadripartite structure with an overall length ranging from 150,069 to 152,206 bp and showed relative stability of inverted repeat borders. The Loteae plastomes demonstrated full collinearity; the most variable sites of the studied plastomes were found in petN-trnC and rps16-accD spacers from the LSC region and in the ycf1 gene within the SSC. All inferred relationships attained maximal support with the Hippocrepis lineage separated first, followed by Coronilla and Anthyllis; Lotus is a sister group to the clade Acmispon + Ornithopus. In this study, completely resolved relationships representing a backbone of plastid phylogeny were produced. The obtained results demonstrated that plastid genomes in the tribe Loteae are structurally conservative in contrast to the closely related tribes Robinieae and Sesbanieae.

Rampant intraspecific variation of plastid genomes in Gentiana section Chondrophyllae

Exploring the level of intraspecific diversity in taxa experienced radiation is helpful to understanding speciation and biodiversity assembly. Gentiana section Chondrophyllae sensu lato encompasses more than 180 species and occupies more a half of species in the genus. In this study, we collected samples across the range of three species (Gentiana aristata, G. crassuloides and G. haynaldii) in section Chondrophyllae s.l., and recovered the intra-species variation by comparing with closely related taxon. Using 25 newly sequenced plastid genomes together with previously published data, we compared structural differences, quantified the variations in plastome size, and measured nucleotide diversity in various regions. Our results showed that the plastome size variation in the three Chondrophyllae species ranged from 285 to 628 bp, and the size variation in LSC, IR and SSC ranged from 236 to 898 bp, 52 to 393 bp and 135 to 356 bp, respectively. Nucleotide diversity of plastome or any of the four regions was much higher than the control species. The average nucleotide diversity in plastomes of the three species ranged from 0.0010 to 0.0023 in protein coding genes, and from 0.0023 to 0.0061 in intergenic regions. More repeat sequence variations were detected within the three Chondrophyllae species than the control species. Various plastid sequence matrixes resulted in different backbone topology in two target species, showed uncertainty in phylogenetic relationship based inference. In conclusion, our results recovered that species of G. section Chondrophyllae s.l. has high intraspecific plastome variation, and provided insights into the radiation in this speciose lineage.

LC-MS/DIA-based strategy for comprehensive flavonoid profiling: an Ocotea spp. applicability case

We introduce a liquid chromatography - mass spectrometry with data-independent acquisition (LC-MS/DIA)-based strategy, specifically tailored to achieve comprehensive and reliable glycosylated flavonoid profiling. This approach facilitates in-depth and simultaneous exploration of all detected precursors and fragments during data processing, employing the widely-used open-source MZmine 3 software. It was applied to a dataset of six Ocotea plant species. This framework suggested 49 flavonoids potentially newly described for these plant species, alongside 45 known features within the genus. Flavonols kaempferol and quercetin, both exhibiting O-glycosylation patterns, were particularly prevalent. Gas-phase fragmentation reactions further supported these findings. For the first time, the apigenin flavone backbone was also annotated in most of the examined Ocotea species. Apigenin derivatives were found mainly in the C-glycoside form, with O. porosa displaying the highest flavone : flavonol ratio. The approach also allowed an unprecedented detection of kaempferol and quercetin in O. porosa species, and it has underscored the untapped potential of LC-MS/DIA data for broad and reliable flavonoid profiling. Our study annotated more than 50 flavonoid backbones in each species, surpassing the current literature.

Intraspecific differentiation of Lindera obtusiloba as revealed by comparative plastomic and evolutionary analyses

Lindera obtusiloba Blume is the northernmost tree species in the family Lauraceae, and it is a key species in understanding the evolutionary history of this family. The species of L. obtusiloba in East Asia has diverged into the Northern and Southern populations, which are geographically separated by an arid belt. Though the morphological differences between populations have been observed and well documented, intraspecific variations at the plastomic level have not been systematically investigated to date. Here, ten chloroplast genomes of L. obtusiloba individuals were sequenced and analyzed along with three publicly available plastomes. Comparative plastomic analysis suggests that both the Northern and the Southern populations share similar overall structure, gene order, and GC content in their plastomes although the size of the plasome and the level of intraspecific variability do vary between the two populations. The Northern have relatively larger plastomes while the Southern population possesses higher intraspecific variability, which could be attributed to the complexity of the geological environments in the South. Phylogenomic analyses also support the split of the Northern and Southern clades among L. obtusiloba individuals. However, there is no obvious species boundary between var. obtusiloba and var. heterophylla in the Southern population, indicating that gene flow could still occur between these two varieties, and this could be used as a good example of reticulate evolution. It is also found that a few photosynthesis-related genes are under positive selection, which is mainly related to the geological and environmental differences between the Northern and the Southern regions. Our results provide a reference for phylogenetic analysis within species and suggest that phylogenomic analyses with a sufficient number of nuclear and chloroplast genomic target loci from widely distributed individuals could provide a deeper understanding of the population evolution of the widespread species.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

New insights into the plastome evolution of Lauraceae using herbariomics

BACKGROUND

The family Lauraceae possesses ca. 50 genera and 2,500-3,000 species that are distributed in the pantropics. Only half of the genera of the family were represented in previously published plastome phylogenies because of the difficulty of obtaining research materials. Plastomes of Hypodaphnideae and the Mezilaurus group, two lineages with unusual phylogenetic positions, have not been previously reported and thus limit our full understanding on the plastome evolution of the family. Herbariomics, promoted by next generation sequencing technology, can make full use of herbarium specimens, and provides opportunities to fill the sampling gap.

RESULTS

In this study, we sequenced five new plastomes (including four genera which are reported for the first time, viz. Chlorocardium, Hypodaphnis, Licaria and Sextonia) from herbarium specimens using genome skimming to conduct a comprehensive analysis of plastome evolution of Lauraceae as a means of sampling representatives of all major clades of the family. We identified and recognized six types of plastomes and revealed that at least two independent loss events at the IR-LSC boundary and an independent expansion of SSC occurred in the plastome evolution of the family. Hypodaphnis possesses the ancestral type of Lauraceae with trnI-CAU, rpl23 and rpl2 duplicated in the IR regions (Type-I). The Mezilaurus group shares the same plastome structure with the core Lauraceae group in the loss of trnI-CAU, rpl23 and rpl2 in the IRa region (Type-III). Two new types were identified in the Ocotea group: (1) the insertion of trnI-CAU between trnL-UAG and ccsA in the SSC region of Licaria capitata and Ocotea bracteosa (Type-IV), and (2) trnI-CAU and pseudogenizated rpl23 inserted in the same region of Nectandra angustifolia (Type-V). Our phylogeny suggests that Lauraceae are divided into nine major clades largely in accordance with the plastome types. The Hypodaphnideae are the earliest diverged lineage supported by both robust phylogeny and the ancestral plastome type. The monophyletic Mezilaurus group is sister to the core Lauraceae.

CONCLUSIONS

By using herbariomics, we built a more complete picture of plastome evolution and phylogeny of the family, thus providing a convincing case for further use of herbariomics in phylogenetic studies of the Lauraceae.

Plastome evolution and phylogenomic insights into the evolution of Lysimachia (Primulaceae: Myrsinoideae)

BACKGROUND

Lysimachia L., the second largest genus within the subfamily Myrsinoideae of Primulaceae, comprises approximately 250 species worldwide. China is the species diversity center of Lysimachia, containing approximately 150 species. Despite advances in the backbone phylogeny of Lysimachia, species-level relationships remain poorly understood due to limited genomic information. This study analyzed 50 complete plastomes for 46 Lysimachia species. We aimed to identify the plastome structure features and hypervariable loci of Lysimachia. Additionally, the phylogenetic relationships and phylogenetic conflict signals in Lysimachia were examined.

RESULTS

These fifty plastomes within Lysimachia had the typical quadripartite structure, with lengths varying from 152,691 to 155,784 bp. Plastome size was positively correlated with IR and intron length. Thirteen highly variable regions in Lysimachia plastomes were identified. Additionally, ndhB, petB and ycf2 were found to be under positive selection. Plastid ML trees and species tree strongly supported that L. maritima as sister to subg. Palladia + subg. Lysimachia (Christinae clade), while the nrDNA ML tree clearly placed L. maritima and subg. Palladia as a sister group.

CONCLUSIONS

The structures of these plastomes of Lysimachia were generally conserved, but potential plastid markers and signatures of positive selection were detected. These genomic data provided new insights into the interspecific relationships of Lysimachia, including the cytonuclear discordance of the position of L. maritima, which may be the result of ghost introgression in the past. Our findings have established a basis for further exploration of the taxonomy, phylogeny and evolutionary history within Lysimachia.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

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.

Untying the Gordian knot of plastid phylogenomic conflict: A case from ferns

Phylogenomic studies based on plastid genome have resolved recalcitrant relationships among various plants, yet the phylogeny of Dennstaedtiaceae at the level of family and genera remains unresolved due to conflicting plastid genes, limited molecular data and incomplete taxon sampling of previous studies. The present study generated 30 new plastid genomes of Dennstaedtiaceae (9 genera, 29 species), which were combined with 42 publicly available plastid genomes (including 24 families, 27 genera, 42 species) to explore the evolution of Dennstaedtiaceae. In order to minimize the impact of systematic errors on the resolution of phylogenetic inference, we applied six strategies to generate 30 datasets based on CDS, intergenic spacers, and whole plastome, and two tree inference methods (maximum-likelihood, ML; and multispecies coalescent, MSC) to comprehensively analyze the plastome-scale data. Besides, the phylogenetic signal among all loci was quantified for controversial nodes using ML framework, and different topologies hypotheses among all datasets were tested. The species trees based on different datasets and methods revealed obvious conflicts at the base of the polypody ferns. The topology of the "CDS-codon-align-rm3" (CDS with the removal of the third codon) matrix was selected as the primary reference or summary tree. The final phylogenetic tree supported Dennstaedtiaceae as the sister group to eupolypods, and Dennstaedtioideae was divided into four clades with full support. This robust reconstructed phylogenetic backbone establishes a framework for future studies on Dennstaedtiaceae classification, evolution and diversification. The present study suggests considering plastid phylogenomic conflict when using plastid genomes. From our results, reducing saturated genes or sites can effectively mitigate tree conflicts for distantly related taxa. Moreover, phylogenetic trees based on amino acid sequences can be used as a comparison to verify the confidence of nucleotide-based trees.

A comparative phytochemical study of nine Lauraceae species by using chemometric data analysis

The diversity of secondary metabolites of individual plants results from multiple enzymatic processes in planta and various environmental factors, such as temperature, moisture, and soil conditions. Chemical composition analysis of plants can lead to a new method to understand relationship among comparable plants along with biological classification such as genetic and anatomical method. In this study, the chemical diversity of nine different Lauraceae species was investigated, and the plant samples were chemically analyzed and classified. Multivariate analysis methods, such as PLS-DA, were used to select important metabolites distinguishing the nine Lauraceae species. The selected metabolites were identified through preparative LC-MS or MS/MS fragment pattern analysis. In addition, the chemical dendrogram for the nine Lauraceae species was interpreted through molecular network analysis and compared with the genetic dendrogram. This approach enabled us to compare the complete chemical compositions of multiple plant samples to identify relationships among plants.

Plastome structure, phylogenomics, and divergence times of tribe Cinnamomeae (Lauraceae)

BACKGROUND

Tribe Cinnamomeae is a species-rich and ecologically important group in tropical and subtropical forests. Previous studies explored its phylogenetic relationships and historical biogeography using limited loci, which might result in biased molecular dating due to insufficient parsimony-informative sites. Thus, 15 plastomes were newly sequenced and combined with published plastomes to study plastome structural variations, gene evolution, phylogenetic relationships, and divergence times of this tribe.

RESULTS

Among the 15 newly generated plastomes, 14 ranged from 152,551 bp to 152,847 bp, and the remaining one (Cinnamomum chartophyllum XTBGLQM0164) was 158,657 bp. The inverted repeat (IR) regions of XTBGLQM0164 contained complete ycf2, trnI^CAU, rpl32, and rpl2. Four hypervariable plastid loci (ycf1, ycf2, ndhF-rpl32-trnL^UAG, and petA-psbJ) were identified as candidate DNA barcodes. Divergence times based on a few loci were primarily determined by prior age constraints rather than by DNA data. In contrast, molecular dating using complete plastid protein-coding genes (PCGs) was determined by DNA data rather than by prior age constraints. Dating analyses using PCGs showed that Cinnamomum sect. Camphora diverged from C. sect. Cinnamomum in the late Oligocene (27.47 Ma).

CONCLUSIONS

This study reports the first case of drastic IR expansion in tribe Cinnamomeae, and indicates that plastomes have sufficient parsimony-informative sites for molecular dating. Besides, the dating analyses provide preliminary insights into the divergence time within tribe Cinnamomeae and can facilitate future studies on its historical biogeography.

Complete plastid genome of Litsea honghoensis Liou 1933, an endangered evergreen species in China

Litsea honghoensis Liou is an endangered tree endemic to south China. In this study, the first complete plastid genome of L. honghoensis was presented, which had a length of 152,605 base pairs (bp) with a GC content of 39.20%. The genome consisted of a large single-copy (LSC) region of 93,560 bp, a small single-copy (SSC) region of 18,905 bp, and two inverted repeat regions (IRa and IRb) of 20,070 bp. There were 125 genes in the plastid genome, including 81 protein-coding genes, 36 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. Phylogenomic analysis based on 52 complete plastomes of Laureae in the family Lauraceae supports the close relationships among L. honghoensis, Lindera communis, Lindera nacusua, Lindera angustifolia, and Lindera glauca.

Fingerprints of climatic changes through the late Cenozoic in southern Asian flora: Magnolia section Michelia (Magnoliaceae)

BACKGROUND AND AIMS

Ongoing global warming is a challenge for humankind. A series of drastic climatic changes have been proven to have occurred throughout the Cenozoic based on a variety of geological evidence, which helps to better understand our planet's future climate. Notably, extant biomes have recorded drastic environmental shifts. The climate in southern Asia, which hosts high biodiversity, is deeply impacted by the Asian monsoon. The origins and evolutionary dynamics of biomes occurring between the tropics and sub-tropics in southern Asia have probably been deeply impacted by climatic changes; however, these aspects remain poorly studied. We tested whether the evolutionary dynamics of the above biomes have recorded the drastic, late Cenozoic environmental shifts, by focusing on Magnolia section Michelia of the family Magnoliaceae.

METHODS

We established a fine time-calibrated phylogeny of M. section Michelia based on complete plastid genomes and inferred its ancestral ranges. Finally, we estimated the evolutionary dynamics of this section through time, determining its diversification rate and the dispersal events that occurred between tropical and sub-tropical areas.

KEY RESULTS

The tropical origin of M. section Michelia was dated to the late Oligocene; however, the diversification of its core group (i.e. M. section Michelia subsection Michelia) has occurred mainly from the late Miocene onward. Two key evolutionary shifts (dated approx. 8 and approx. 3 million years ago, respectively) were identified, each of them probably in response to drastic climatic changes.

CONCLUSION

Here, we inferred the underlying evolutionary dynamics of biomes in southern Asia, which probably reflect late Cenozoic climatic changes. The occurrence of modern Asian monsoons was probably fundamental for the origin of M. section Michelia; moreover, the occurrence of asymmetric dispersal events between the tropics and sub-tropics hint at an adaptation strategy of M. section Michelia to global cooling, in agreement with the tropical conservatism hypothesis.

Plastid genome evolution of a monophyletic group in the subtribe Lauriineae (Laureae, Lauraceae)

Litsea, a non-monophyletic group of the tribe Laureae (Lauraceae), plays important roles in the tropical and subtropical forests of Asia, Australia, Central and North America, and the islands of the Pacific. However, intergeneric relationships between Litsea and Laurus, Lindera, Parasassafras and Sinosassafras of the tribe Laureae remain unresolved. In this study, we present phylogenetic analyses of seven newly sequenced Litsea plastomes, together with 47 Laureae plastomes obtained from public databases, representing six genera of the Laureae. Our results highlight two highly supported monophyletic groups of Litsea taxa. One is composed of 16 Litsea taxa and two Lindera taxa. The 18 plastomes of these taxa were further compared for their gene structure, codon usage, contraction and expansion of inverted repeats, sequence repeats, divergence hotspots, and gene evolution. The complete plastome size of newly sequenced taxa varied between 152,377 bp (Litsea auriculata) and 154,117 bp (Litsea pierrei). Seven of the 16 Litsea plastomes have a pair of insertions in the IRa (trnL-trnH) and IRb (ycf2) regions. The 18 plastomes of Litsea and Lindera taxa exhibit similar gene features, codon usage, oligonucleotide repeats, and inverted repeat dynamics. The codons with the highest frequency among these taxa favored A/T endings and each of these plastomes had nine divergence hotspots, which are located in the same regions. We also identified six protein coding genes (accD, ndhJ, rbcL, rpoC2, ycf1 and ycf2) under positive selection in Litsea; these genes may play important roles in adaptation of Litsea species to various environments.

Molecular evolution and phylogenetic relationships of Ligusticum (Apiaceae) inferred from the whole plastome sequences

Background

The genus Ligusticum belongs to Apiaceae, and its taxonomy has long been a major difficulty. A robust phylogenetic tree is the basis of accurate taxonomic classification of Ligusticum. We herein used 26 (including 14 newly sequenced) plastome-scale data to generate reliable phylogenetic trees to explore the phylogenetic relationships of Chinese Ligusticum.

Results

We found that these plastid genomes exhibited diverse plastome characteristics across all four currently identified clades in China, while the plastid protein-coding genes were conserved. The phylogenetic analyses by the concatenation and coalescent methods obtained a more robust molecular phylogeny than prior studies and showed the non-monophyly of Chinese Ligusticum. In the concatenation-based phylogeny analyses, the two datasets yielded slightly different topologies that may be primarily due to the discrepancy in the number of variable sites.

Conclusions

Our plastid phylogenomics analyses emphasized that the current circumscription of the Chinese Ligusticum should be reduced, and the taxonomy of Ligusticum urgently needs revision. Wider taxon sampling including the related species of Ligusticum will be necessary to explore the phylogenetic relationships of this genus. Overall, our study provided new insights into the taxonomic classification of Ligusticum and would serve as a framework for future studies on taxonomy and delimitation of Ligusticum from the perspective of the plastid genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02010-z.

A comparative analysis of complete chloroplast genomes of seven Ocotea species (Lauraceae) confirms low sequence divergence within the Ocotea complex

The genus Ocotea (Lauraceae) includes about 450 species, of which about 90% are Neotropical, while the rest is from Macaronesia, Africa and Madagascar. In this study we present the first complete chloroplast genome sequences of seven Ocotea species, six Neotropical and one from Macaronesia. Genome sizes range from 152,630 (O. porosa) to 152,685 bp (O. aciphylla). All seven plastomes contain a total of 131 (114 unique) genes, among which 87 (80 unique) encode proteins. The order of genes (if present) is the same in all Lauraceae examined so far. Two hypervariable loci were found in the LSC region (psbA-trnH, ycf2), three in the SSC region (ycf1, ndhH, trnL(UAG)-ndhF). The pairwise cp genomic alignment between the taxa showed that the LSC and SSC regions are more variable compared to the IR regions. The protein coding regions comprise 25,503-25,520 codons in the Ocotea plastomes examined. The most frequent amino acids encoded in the plastomes were leucine, isoleucine, and serine. SSRs were found to be more frequent in the two dioecious Neotropical Ocotea species than in the four bisexual species and the gynodioecious species examined (87 vs. 75-84 SSRs). A preliminary phylogenetic analysis based on 69 complete plastomes of Lauraceae species shows the seven Ocotea species as sister group to Cinnamomum sensu lato. Sequence divergence among the Ocotea species appears to be much lower than among species of the most closely related, likewise species-rich genera Cinnamomum, Lindera and Litsea.

Plastid phylogenomics of tribe Perseeae (Lauraceae) yields insights into the evolution of East Asian subtropical evergreen broad-leaved forests

BACKGROUND

The East Asian subtropical evergreen broad-leaved forests (EBLFs) harbor remarkable biodiversity. However, their historical assembly remains unclear. To gain new insights into the assembly of this biome, we generated a molecular phylogeny of one of its essential plant groups, the tribe Perseeae (Lauraceae).

RESULTS

Our plastid tree topologies were robust to analyses based on different plastid regions and different strategies for data partitioning, nucleotide substitution saturation, and gap handling. We found that tribe Perseeae comprised six major clades and began to colonize the subtropical EBLFs of East Asia in the early Miocene. The diversification rates of tribe Perseeae accelerated twice in the late Miocene.

CONCLUSIONS

Our findings suggest that the intensified precipitation in East Asia in the early Miocene may have facilitated range expansions of the subtropical EBLFs and establishment of tribe Perseeae within this biome. By the late Miocene, species assembly and diversification within the EBLFs had become rapid.

Leaf epidermal micromorphology defining the clades in Cinnamomum (Lauraceae)

In this study, we sampled 48 species of Asian Cinnamomum covering the species groups that were identified in recent phylogenetic studies and conducted leaf micromorphological observations using both light microscopy (LM) and scanning electron microscopy (SEM). Synapomorphies were determined by means of mapping micromorphological characters on a phylogenetic tree. The results indicate that Cinnamomum exhibits two different types of leaf upper epidermis: Type I has smooth/non-reticulate periclinal walls whereas Type II has reticulate periclinal walls and is unusual in the family Lauraceae. We found that the two types of micromorphological characters are clade-specific, sect. Camphora s.s. possesses Type I leaf upper epidermis, and sect. Cinnamomum s.l. has Type II leaf upper epidermis. Our study also reveals that C.saxatile, a member of sect. Camphora s.l. in the traditional classification, actually has Type II leaf upper epidermis, thus reinforcing the result of a recent molecular phylogeny that has this species in a clade consisting mainly of species of sect. Cinnamomum.

Complete plastid genome sequence of Neolitsea aciculata (Laurales: Lauraceae), an evergreen broad-leaved tree endemic to East Asia

In this study, we sequenced the complete plastid genome (plastome) of Neolitsea aciculata, an evergreen broad-leaved tree endemic to East Asia, a woody component of East Asian warm-temperate and subtropical forests across China, Korea, and Japan. The plastome of N. aciculata is assembled as a single contig (152,722 bp). A large and a small single copy (93,785 and 18,795 bp, respectively) of the genome are separated by a pair of inverted repeats (20,071 bp). The genome consists of 126 genes, including 80 protein-coding, eight ribosomal RNA, and 36 transfer RNA genes. Two genes in the IR region (ycf1 and ycf2) are pseudogenized. Our phylogenetic analysis revealed the phylogenetic position of N. aciculata in a highly supported clade of the genus Neolitsea along with other two congeners, N. pallens and N. sericea.

Defining coalescent genes: Theory meets practice in organelle phylogenomics

The species tree paradigm that dominates current molecular systematic practice infers species trees from collections of sequences under assumptions of the multispecies coalescent (MSC), i.e., that there is free recombination between the sequences and no (or very low) recombination within them. These coalescent genes (c-genes) are thus defined in an historical rather than molecular sense, and can in theory be as large as an entire genome or as small as a single nucleotide. A debate about how to define c-genes centers on the contention that nuclear gene sequences used in many coalescent analyses undergo too much recombination, such that their introns comprise multiple c-genes, violating a key assumption of the MSC. Recently a similar argument has been made for the genes of plastid (e.g., chloroplast) and mitochondrial genomes, which for the last 30 or more years have been considered to represent a single c-gene for the purposes of phylogeny reconstruction because they are non-recombining in a historical sense. Consequently, it has been suggested that these genomes should be analyzed using coalescent methods that treat their genes-over 70 protein-coding genes in the case of most plastid genomes (plastomes)-as independent estimates of species phylogeny, in contrast to the usual practice of concatenation, which is appropriate for generating gene trees. However, although recombination certainly occurs in the plastome, as has been recognized since the 1970's, it is unlikely to be phylogenetically relevant. This is because such historically effective recombination can only occur when plastomes with incongruent histories are brought together in the same plastid. However, plastids sort rapidly into different cell lineages and rarely fuse. Thus, because of plastid biology, the plastome is a more canonical c-gene than is the average multi-intron mammalian nuclear gene. The plastome should thus continue to be treated as a single estimate of the underlying species phylogeny, as should the mitochondrial genome. The implications of this long-held insight of molecular systematics for studies in the phylogenomic era are explored.

Data Types and the Phylogeny of Neoaves

The phylogeny of Neoaves, the largest clade of extant birds, has remained unclear despite intense study. The difficulty associated with resolving the early branches in Neoaves is likely driven by the rapid radiation of this group. However, conflicts among studies may be exacerbated by the data type analyzed. For example, analyses of coding exons typically yield trees that place Strisores (nightjars and allies) sister to the remaining Neoaves, while analyses of non-coding data typically yield trees where Mirandornites (flamingos and grebes) is the sister of the remaining Neoaves. Our understanding of data type effects is hampered by the fact that previous analyses have used different taxa, loci, and types of non-coding data. Herein, we provide strong corroboration of the data type effects hypothesis for Neoaves by comparing trees based on coding and non-coding data derived from the same taxa and gene regions. A simple analytical method known to minimize biases due to base composition (coding nucleotides as purines and pyrimidines) resulted in coding exon data with increased congruence to the non-coding topology using concatenated analyses. These results improve our understanding of the resolution of neoavian phylogeny and point to a challenge—data type effects—that is likely to be an important factor in phylogenetic analyses of birds (and many other taxonomic groups). Using our results, we provide a summary phylogeny that identifies well-corroborated relationships and highlights specific nodes where future efforts should focus.