Plastid phylogenomics resolves infrafamilial relationships of the Styracaceae and sheds light on the backbone relationships of the Ericales

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

Identification of closely related species in Aspergillus through Analysis of Whole-Genome

The challenge of discriminating closely related species persists, notably within clinical diagnostic laboratories for invasive aspergillosis (IA)-related species and food contamination microorganisms with toxin-producing potential. We employed Analysis of the whole-GEnome (AGE) to address the challenges of closely related species within the genus Aspergillus and developed a rapid detection method. First, reliable whole genome data for 77 Aspergillus species were downloaded from the database, and through bioinformatic analysis, specific targets for each species were identified. Subsequently, sequencing was employed to validate these specific targets. Additionally, we developed an on-site detection method targeting a specific target using a genome editing system. Our results indicate that AGE has successfully achieved reliable identification of all IA-related species (Aspergillus fumigatus, Aspergillus niger, Aspergillus nidulans, Aspergillus flavus, and Aspergillus terreus) and three well-known species (A. flavus, Aspergillus parasiticus, and Aspergillus oryzae) within the Aspergillus section. Flavi and AGE have provided species-level-specific targets for 77 species within the genus Aspergillus. Based on these reference targets, the sequencing results targeting specific targets substantiate the efficacy of distinguishing the focal species from its closely related species. Notably, the amalgamation of room-temperature amplification and genome editing techniques demonstrates the capacity for rapid and accurate identification of genomic DNA samples at a concentration as low as 0.1 ng/μl within a concise 30-min timeframe. Importantly, this methodology circumvents the reliance on large specialized instrumentation by presenting a singular tube operational modality and allowing for visualized result assessment. These advancements aptly meet the exigencies of on-site detection requirements for the specified species, facilitating prompt diagnosis and food quality monitoring. Moreover, as an identification method based on species-specific genomic sequences, AGE shows promising potential as an effective tool for epidemiological research and species classification.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes

BACKGROUND

Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation.

RESULTS

We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > accD-psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47-56 SSRs and 11-14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support.

CONCLUSION

This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

Comparative and phylogenetic analysis of the complete chloroplast genomes of six Polygonatum species (Asparagaceae)

Polygonatum Miller belongs to the tribe Polygonateae of Asparagaceae. The horizontal creeping fleshy roots of several species in this genus serve as traditional Chinese medicine. Previous studies have mainly reported the size and gene contents of the plastomes, with little information on the comparative analysis of the plastid genomes of this genus. Additionally, there are still some species whose chloroplast genome information has not been reported. In this study, the complete plastomes of six Polygonatum were sequenced and assembled, among them, the chloroplast genome of P. campanulatum was reported for the first time. Comparative and phylogenetic analyses were then conducted with the published plastomes of three related species. Results indicated that the whole plastome length of the Polygonatum species ranged from 154,564 bp (P. multiflorum) to 156,028 bp (P. stenophyllum) having a quadripartite structure of LSC and SSC separated by two IR regions. A total of 113 unique genes were detected in each of the species. Comparative analysis revealed that gene content and total GC content in these species were highly identical. No significant contraction or expansion was observed in the IR boundaries among all the species except P. sibiricum1, in which the rps19 gene was pseudogenized owing to incomplete duplication. Abundant long dispersed repeats and SSRs were detected in each genome. There were five remarkably variable regions and 14 positively selected genes were identified among Polygonatum and Heteropolygonatum. Phylogenetic results based on chloroplast genome strongly supported the placement of P. campanulatum with alternate leaves in sect. Verticillata, a group characterized by whorled leaves. Moreover, P. verticillatum and P. cyrtonema were displayed as paraphyletic. This study revealed that the characters of plastomes in Polygonatum and Heteropolygonatum maintained a high degree of similarity. Five highly variable regions were found to be potential specific DNA barcodes in Polygonatum. Phylogenetic results suggested that leaf arrangement was not suitable as a basis for delimitation of subgeneric groups in Polygonatum and the definitions of P. cyrtonema and P. verticillatum require further study.

A phylogenomic study of Iridaceae Juss. based on complete plastid genome sequences

The plastid genome has proven to be an effective tool for examining deep correlations in plant phylogenetics, owing to its highly conserved structure, uniparental inheritance, and limited variation in evolutionary rates. Iridaceae, comprising more than 2,000 species, includes numerous economically significant taxa that are frequently utilized in food industries and medicines and for ornamental and horticulture purposes. Molecular studies on chloroplast DNA have confirmed the position of this family in the order Asparagales with non-asparagoids. The current subfamilial classification of Iridaceae recognizes seven subfamilies-Isophysioideae, Nivenioideae, Iridoideae, Crocoideae, Geosiridaceae, Aristeoideae, and Patersonioideae-which are supported by limited plastid DNA regions. To date, no comparative phylogenomic studies have been conducted on the family Iridaceae. We assembled and annotated (de novo) the plastid genomes of 24 taxa together with seven published species representing all the seven subfamilies of Iridaceae and performed comparative genomics using the Illumina MiSeq platform. The plastomes of the autotrophic Iridaceae represent 79 protein-coding, 30 tRNA, and four rRNA genes, with lengths ranging from 150,062 to 164,622 bp. The phylogenetic analysis of the plastome sequences based on maximum parsimony, maximum likelihood, and Bayesian inference analyses suggested that Watsonia and Gladiolus were closely related, supported by strong support values, which differed considerably from recent phylogenetic studies. In addition, we identified genomic events, such as sequence inversions, deletions, mutations, and pseudogenization, in some species. Furthermore, the largest nucleotide variability was found in the seven plastome regions, which can be used in future phylogenetic studies. Notably, three subfamilies-Crocoideae, Nivenioideae, and Aristeoideae-shared a common ycf2 gene locus deletion. Our study is a preliminary report of a comparative study of the complete plastid genomes of 7/7 subfamilies and 9/10 tribes, elucidating the structural characteristics and shedding light on plastome evolution and phylogenetic relationships within Iridaceae. Additionally, further research is required to update the relative position of Watsonia within the tribal classification of the subfamily Crocoideae.

Unravelling the palaeobiogeographical history of the living fossil genus Rehderodendron (Styracaceae) with fossil and extant pollen and fruit data

BACKGROUND

The relict genus Rehderodendron (Styracaceae), the species of which are restricted to mostly warm temperate to tropical climate in East Asia today, is known from fossil fruits and pollen in Europe during warmer periods from the lower Eocene to Pliocene. To infer which extant species are most closely related to the fossils, new data of pollen and fruit morphologiesy of six extant species, and additional new data of fossil pollen and previously described fossil fruits of Rehderodendron, are compared.

RESULTS

Both fossil pollen and fruits resemble a morphological mixture of the extant species R. indochinense, R. kwantungense, R. macrocarpum, and R. microcarpum, thus implying that these extant taxa and the fossil European taxa represent an old Eurasian lineage, whereas the pollen and fruit morphology of the extant R. kweichowense and R. truongsonense differ considerably from the fossils and other extant species investigated, and are considered to have evolved independently.

CONCLUSIONS

The palaeobiogeographical history of Rehderodendron reveals that its fossil members of the European lineage were most prominent during climatic optima such as the Palaeocene-Eocene Thermal Maximum (PETM), Early Eocene Climate Optimum (EECO) and Middle Miocene Thermal Maximum (MMTM). However, when during the Pliocene the climate changed to colder and less humid conditions, the genus went extinct in Europe but migrated eastwards, most likely in two dispersal events along the Tethys Sea prior to extinction. One of the former most westerly stepping stones is suggested by the refugial occurrence of R. microcarpum in the southeastern Himalaya, whereas R. macrocarpum and R. kwangtungense, the taxa distributed more to the east, might have migrated eastwards already before the Miocene.

The Complete Chloroplast Genome Sequence of Laportea bulbifera (Sieb. et Zucc.) Wedd. and Comparative Analysis with Its Congeneric Species

Laportea bulbifera (L. bulbifera) is an important medicinal plant of Chinese ethnic minorities, with high economic and medicinal value. However, the medicinal materials of the genus Laportea are prone to be misidentified due to the similar morphological characteristics of the original plants. Thus, it is crucial to discover their molecular marker points and to precisely identify these species for their exploitation and conservation. Here, this study reports detailed information on the complete chloroplast (cp) of L. bulbifera. The result indicates that the cp genome of L. bulbifera of 150,005 bp contains 126 genes, among them, 37 tRNA genes and 81 protein-coding genes. The analysis of repetition demonstrated that palindromic repeats are more frequent. In the meantime, 39 SSRs were also identified, the majority of which were mononucleotides Adenine-Thymine (A-T). Furthermore, we compared L. bulbifera with eight published Laportea plastomes, to explore highly polymorphic molecular markers. The analysis identified four hypervariable regions, including rps16, ycf1, trnC-GCA and trnG-GCC. According to the phylogenetic analysis, L. bulbifera was most closely related to Laportea canadensis (L. canadensis), and the molecular clock analysis speculated that the species originated from 1.8216 Mya. Overall, this study provides a more comprehensive analysis of the evolution of L. bulbifera from the perspective of phylogenetic and intrageneric molecular variation in the genus Laportea, which is useful for providing a scientific basis for further identification, taxonomic, and evolutionary studies of the genus.

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.

Phylogenomics, plastome structure and species identification in Mahonia (Berberidaceae)

BACKGROUND

Elucidating the phylogenetic relationships within species-rich genera is essential but challenging, especially when lineages are assumed to have been going through radiation events. Mahonia Nutt. (Berberidaceae) is a genus with cosmopolitan distribution, comprising approximately 100 species, two of which are known as Caulis Mahoniae (M. bealei and M. fortunei) with crucial pharmacological significance in Chinese herbal medicine. Mahonia is a taxonomically challenging genus, and intrageneric phylogenetic relationships still need to be explored using genome data. Universal DNA barcodes and floral morphological attributes have limited discriminatory power in Mahonia.

RESULTS

We sequenced 17 representative plastomes and integrated three published plastome data together to conduct comparative and phylogenetic analyses. We found that Mahonia and Berberis share a large IR expansion (~ 12 kb), which is recognized as a typical character of Berberideae. Repeated sequences are revealed in the species of Mahonia, which are valuable for further population genetic studies. Using a comparative plastome analysis, we determined eight hypervariable regions whose discriminative power is comparable to that of the whole plastid genomes. The incongruence of the ITS and the plastome tree topologies may be ascribed to ancestral hybridization events and/or to incomplete lineage sorting. In addition, we suggest that leaf epidermal characters could help to distinguish closely related species in Mahonia.

CONCLUSIONS

We propose an integrative approach combining special barcodes and micromorphological traits to circumscribe Mahonia species. The results cast a new light on the development of an integrative method for accurate species circumscription and provide abundant genetic resources for further research on Mahonia.

Comparative plastome analysis of Musaceae and new insights into phylogenetic relationships

Background

Musaceae is an economically important family consisting of 70-80 species. Elucidation of the interspecific relationships of this family is essential for a more efficient conservation and utilization of genetic resources for banana improvement. However, the scarcity of herbarium specimens and quality molecular markers have limited our understanding of the phylogenetic relationships in wild species of Musaceae. Aiming at improving the phylogenetic resolution of Musaceae, we analyzed a comprehensive set of 49 plastomes for 48 species/subspecies representing all three genera of this family.

Results

Musaceae plastomes have a relatively well-conserved genomic size and gene content, with a full length ranging from 166,782 bp to 172,514 bp. Variations in the IR borders were found to show phylogenetic signals to a certain extent in Musa. Codon usage bias analysis showed different preferences for the same codon between species and three genera and a common preference for A/T-ending codons. Among the two genes detected under positive selection (dN/dS > 1), ycf2 was indicated under an intensive positive selection. The divergent hotspot analysis allowed the identification of four regions (ndhF-trnL, ndhF, matK-rps16, and accD) as specific DNA barcodes for Musaceae species.

Bayesian and maximum likelihood phylogenetic analyses using full plastome resulted in nearly identical tree topologies with highly supported relationships between species. The monospecies genus Musella is sister to Ensete, and the genus Musa was divided into two large clades, which corresponded well to the basic number of n = x = 11 and n = x =10/9/7, respectively. Four subclades were divided within the genus Musa. A dating analysis covering the whole Zingiberales indicated that the divergence of Musaceae family originated in the Palaeocene (59.19 Ma), and the genus Musa diverged into two clades in the Eocene (50.70 Ma) and then started to diversify from the late Oligocene (29.92 Ma) to the late Miocene. Two lineages (Rhodochlamys and Australimusa) radiated recently in the Pliocene /Pleistocene periods.

Conclusions

The plastome sequences performed well in resolving the phylogenetic relationships of Musaceae and generated new insights into its evolution. Plastome sequences provided valuable resources for population genetics and phylogenetics at lower taxon.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08454-3.

Comparative genomics and phylogenetic relationships of two endemic and endangered species (Handeliodendron bodinieri and Eurycorymbus cavaleriei) of two monotypic genera within Sapindales

BACKGROUND

Handeliodendron Rehder and Eurycorymbus Hand.-Mazz. are the monotypic genera in the Sapindaceae family. The phylogenetic relationship of these endangered species Handeliodendron bodinieri (Lévl.) Rehd. and Eurycorymbus cavaleriei (Lévl.) Rehd. et Hand.-Mazz. with other members of Sapindaceae s.l. is not well resolved. A previous study concluded that the genus Aesculus might be paraphyletic because Handeliodendron was nested within it based on small DNA fragments. Thus, their chloroplast genomic information and comparative genomic analysis with other Sapindaceae species are necessary and crucial to understand the circumscription and plastome evolution of this family.

RESULTS

The chloroplast genome sizes of Handeliodendron bodinieri and Eurycorymbus cavaleriei are 151,271 and 158,690 bp, respectively. Results showed that a total of 114 unique genes were annotated in H. bodinieri and E. cavaleriei, and the ycf1 gene contained abundant SSRs in both genomes. Comparative analysis revealed that gene content, PCGs, and total GC content were remarkably similar or identical within 13 genera from Sapindaceae, and the chloroplast genome size of four genera was generally smaller within the family, including Acer, Dipteronia, Aesculus, and Handeliodendron. IR boundaries of the H. bodinieri showed a significant contraction, whereas it presented a notable expansion in E. cavaleriei cp genome. Ycf1, ndhC-trnV-UAC, and rpl32-trnL-UAG-ccsA were remarkably divergent regions in the Sapindaceae species. Analysis of selection pressure showed that there are a few positively selected genes. Phylogenetic analysis based on different datasets, including whole chloroplast genome sequences, coding sequences, large single-copy, small single-copy, and inverted repeat regions, consistently demonstrated that H. bodinieri was sister to the clade consisting of Aesculus chinensis and A. wangii and strongly support Eurycorymbus cavaleriei as sister to Dodonaea viscosa.

CONCLUSION

This study revealed that the cp genome size of the Hippocastanoideae was generally smaller compared to the other subfamilies within Sapindaceae, and three highly divergent regions could be used as the specific DNA barcodes within Sapindaceae. Phylogenetic results strongly support that the subdivision of four subfamilies within Sapindaceae, and Handeliodendron is not nested within the genus Aesculus.

Comparative plastomes of Carya species provide new insights into the plastomes evolution and maternal phylogeny of the genus

Carya, in the Juglandiodeae subfamily, is to a typical temperate-subtropical forest-tree genus for studying the phylogenetic evolution and intercontinental disjunction between eastern Asia (EA) and North America (NA). Species of the genus have high economic values worldwide for their high-quality wood and the rich healthy factors of their nuts. Although previous efforts based on multiple molecular markers or genome-wide SNPs supported the monophyly of Carya and its two EA and NA major subclades, the maternal phylogeny of Carya still need to be comprehensively evaluated. The variation of Carya plastome has never been thoroughly characterized. Here, we novelly present 19 newly generated plastomes of congeneric Carya species, including the recently rediscovered critically endangered C. poilanei. The overall assessment of plastomes revealed highly conservative in the general structures. Our results indicated that remarkable differences in several plastome features are highly consistent with the EA-NA disjunction and showed the relatively diverse matrilineal sources among EA Carya compared to NA Carya. The maternal phylogenies were conducted with different plastome regions and full-length plastome datasets from 30 plastomes, representing 26 species in six genera of Juglandoideae and Myrica rubra (as root). Six out of seven phylogenetic topologies strongly supported the previously reported relationships among genera of Juglandoideae and the two subclades of EA and NA Carya, but displayed significant incongruencies between species within the EA and NA subclades. The phylogenetic tree generated from full-length plastomes demonstrated the optimal topology and revealed significant geographical maternal relationships among Carya species, especially for EA Carya within overlapping distribution areas. The full-length plastome-based phylogenetic topology also strongly supported the taxonomic status of five controversial species as separate species of Carya. Historical and recent introgressive hybridization and plastid captures might contribute to plastome geographic patterns and inconsistencies between topologies built from different datasets, while incomplete lineage sorting could account for the discordance between maternal topology and the previous nuclear genome data-based phylogeny. Our findings highlight full-length plastomes as an ideal tool for exploring maternal relationships among the subclades of Carya, and potentially in other outcrossing perennial woody plants, for resolving plastome phylogenetic relationships.

Deep Insights Into the Plastome Evolution and Phylogenetic Relationships of the Tribe Urticeae (Family Urticaceae)

Urticeae s.l., a tribe of Urticaceae well-known for their stinging trichomes, consists of more than 10 genera and approximately 220 species. Relationships within this tribe remain poorly known due to the limited molecular and taxonomic sampling in previous studies, and chloroplast genome (CP genome/plastome) evolution is still largely unaddressed. To address these concerns, we used genome skimming data-CP genome and nuclear ribosomal DNA (18S-ITS1-5.8S-ITS2-26S); 106 accessions-for the very first time to attempt resolving the recalcitrant relationships and to explore chloroplast structural evolution across the group. Furthermore, we assembled a taxon rich two-locus dataset of trnL-F spacer and ITS sequences across 291 accessions to complement our genome skimming dataset. We found that Urticeae plastomes exhibit the tetrad structure typical of angiosperms, with sizes ranging from 145 to 161 kb and encoding a set of 110-112 unique genes. The studied plastomes have also undergone several structural variations, including inverted repeat (IR) expansions and contractions, inversion of the trnN-GUU gene, losses of the rps19 gene, and the rpl2 intron, and the proliferation of multiple repeat types; 11 hypervariable regions were also identified. Our phylogenomic analyses largely resolved major relationships across tribe Urticeae, supporting the monophyly of the tribe and most of its genera except for Laportea, Urera, and Urtica, which were recovered as polyphyletic with strong support. Our analyses also resolved with strong support several previously contentious branches: (1) Girardinia as a sister to the Dendrocnide-Discocnide-Laportea-Nanocnide-Zhengyia-Urtica-Hesperocnide clade and (2) Poikilospermum as sister to the recently transcribed Urera sensu stricto. Analyses of the taxon-rich, two-locus dataset showed lower support but was largely congruent with results from the CP genome and nuclear ribosomal DNA dataset. Collectively, our study highlights the power of genome skimming data to ameliorate phylogenetic resolution and provides new insights into phylogenetic relationships and chloroplast structural evolution in Urticeae.

Chloroplast genome sequencing based on genome skimming for identification of Eriobotryae Folium

Background

Whole chloroplast genome (cpDNA) sequence is becoming widely used in the phylogenetic studies of plant and species identification, but in most cases the cpDNA were acquired from silica gel dried fresh leaves. So far few reports have been available to describe cpDNA acquisition from crude drugs derived from plant materials, the DNA of which usually was seriously damaged during their processing. In this study, we retrieved cpDNA from the commonly used crude drug Eriobotryae Folium (Pipaye in Chinese, which is the dried leaves of Eriobotrya japonica, PPY) using genome skimming technique.

Results

We successfully recovered cpDNA sequences and rDNA sequences from the crude drug PPY, and bioinformatics analysis showed a high overall consistency between the cpDNA obtained from the crude drugs and fresh samples. In the ML tree, each species formed distinct monophyletic clades based on cpDNA sequence data, while the phylogenetic relationships between Eriobotrya species were poorly resolved based on ITS and ITS2.

Conclusion

Our results demonstrate that both cpDNA and ITS/ITS2 are effective for identifying PPY and its counterfeits derived from distantly related species (i.e. Dillenia turbinata and Magnolia grandiflora), but cpDNA is more effective for distinguishing the counterfeits derived from the close relatives of Eriobotrya japonica, suggesting the potential of genome skimming for retrieving cpDNA from crude drugs used in Traditional Chinese Medicine for their identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00728-0.

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.

The chloroplast genomes of four Bupleurum (Apiaceae) species endemic to Southwestern China, a diversity center of the genus, as well as their evolutionary implications and phylogenetic inferences

Background

As one of the largest genera in Apiaceae, Bupleurum L. is well known for its high medicinal value. The genus has frequently attracted the attention of evolutionary biologist and taxonomist for its distinctive characteristics in the Apiaceae family. Although some chloroplast genomes data have been now available, the changes in the structure of chloroplast genomes and selective pressure in the genus have not been fully understood. In addition, few of the species are endemic to Southwest China, a distribution and diversity center of Chinese Bupleurum. Endemic species are key components of biodiversity and ecosystems, and investigation of the chloroplast genomes features of endemic species in Bupleurum will be helpful to develop a better understanding of evolutionary process and phylogeny of the genus. In this study, we analyzed the sequences of whole chloroplast genomes of 4 Southwest China endemic Bupleurum species in comparison with the published data of 17 Bupleurum species to determine the evolutionary characteristics of the genus and the phylogenetic relationships of Asian Bupleurum.

Results

The complete chloroplast genome sequences of the 4 endemic Bupleurum species are 155,025 bp to 155,323 bp in length including a SSC and a LSC region separated by a pair of IRs. Comparative analysis revealed an identical chloroplast gene content across the 21 Bupleurum species, including a total of 114 unique genes (30 tRNA genes, 4 rRNA genes and 80 protein-coding genes). Chloroplast genomes of the 21 Bupleurum species showed no rearrangements and a high sequence identity (96.4–99.2%). They also shared a similar tendency of SDRs and SSRs, but differed in number (59–83). In spite of their high conservation, they contained some mutational hotspots, which can be potentially exploited as high-resolution DNA barcodes for species discrimination. Selective pressure analysis showed that four genes were under positive selection. Phylogenetic analysis revealed that the 21 Bupleurum formed two major clades, which are likely to correspond to their geographical distribution.

Conclusions

The chloroplast genome data of the four endemic Bupleurum species provide important insights into the characteristics and evolution of chloroplast genomes of this genu, and the phylogeny of Bupleurum.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08008-z.

Plastome structure and phylogenetic relationships of Styracaceae (Ericales)

BACKGROUND

The Styracaceae are a woody, dicotyledonous family containing 12 genera and an estimated 160 species. Recent studies have shown that Styrax and Sinojackia are monophyletic, Alniphyllum and Bruinsmia cluster into a clade with an approximately 20-kb inversion in the Large Single-Copy (LSC) region. Halesia and Pterostyrax are not supported as monophyletic, while Melliodendron and Changiostyrax always form sister clades. Perkinsiodendron and Changiostyrax are newly established genera of Styracaceae. However, the phylogenetic relationship of Styracaceae at the generic level needs further research.

RESULTS

We collected 28 complete plastomes of Styracaceae, including 12 sequences newly reported here and 16 publicly available sequences, comprising 11 of the 12 genera of Styracaceae. All species possessed the typical quadripartite structure of angiosperm plastomes, with sequence differences being minor, except for a large 20-kb (14 genes) inversion found in Alniphyllum and Bruinsmia. Seven coding sequences (rps4, rpl23, accD, rpoC1, psaA, rpoA and ndhH) were identified to possess positively selected sites. Phylogenetic reconstructions based on seven data sets (i.e., LSC, SSC, IR, Coding, Non-coding, combination of LSC + SSC and concatenation of LSC + SSC + one IR) produced similar topologies. In our analyses, all genera were strongly supported as monophyletic. Styrax was sister to the remaining genera. Alniphyllum and Bruinsmia form a clade. Halesia diptera does not cluster with Perkinsiodendron, while Perkinsiodendron and Rehderodendron form a clade. Changiostyrax is sister to a clade of Pterostyrax and Sinojackia.

CONCLUSION

Overall, our results demonstrate the power of plastid phylogenomics in improving estimates of phylogenetic relationships among genera. This study also provides insight into plastome evolution across Styracaceae.

Complete Chloroplast Genome of Clethra fargesii Franch., an Original Sympetalous Plant from Central China: Comparative Analysis, Adaptive Evolution, and Phylogenetic Relationships

Clethra fargesii, an essential ecological and endemic woody plant of the genus Clethra in Clethraceae, is widely distributed in Central China. So far, there have been a paucity of studies on its chloroplast genome. In the present study, we sequenced and assembled the complete chloroplast genome of C. fargesii. We also analyzed the chloroplast genome features and compared them to Clethra delavayi and other closely related species in Ericales. The complete chloroplast genome is 157,486 bp in length, including a large single-copy (LSC) region of 87,034 bp and a small single-copy (SSC) region of 18,492 bp, separated by a pair of inverted repeat (IR) regions of 25,980 bp. The GC content of the whole genome is 37.3%, while those in LSC, SSC, and IR regions are 35.4%, 30.7%, and 43.0%, respectively. The chloroplast genome of C. fargesii encodes 132 genes in total, including 87 protein-coding genes (PCGs), 37 tRNA genes, and eight rRNA genes. A total of 26,407 codons and 73 SSRs were identified in C. fargesii chloroplast genome. Additionally, we postulated and demonstrated that the structure of the chloroplast genome in Clethra species may present evolutionary conservation based on the comparative analysis of genome features and genome alignment among eight Ericales species. The low Pi values revealed evolutionary conservation based on the nucleotide diversity analysis of chloroplast genome in two Clethra species. The low selection pressure was shown by a few positively selected genes by adaptive evolution analysis using 80 coding sequences (CDSs) of the chloroplast genomes of two Clethra species. The phylogenetic tree showed that Clethraceae and Ericaceae are sister clades, which reconfirm the previous hypothesis that Clethra is highly conserved in the chloroplast genome using 75 CDSs of chloroplast genome among 40 species. The genome information and analysis results presented in this study are valuable for further study on the intraspecies identification, biogeographic analysis, and phylogenetic relationship in Clethraceae.

The large single-copy (LSC) region functions as a highly effective and efficient molecular marker for accurate authentication of medicinal Dendrobium species

Having great medicinal values, Dendrobium species of "Fengdou" (DSFs) are a taxonomically complex group in Dendrobium genus including many closely related and recently diverged species. Traditionally used DNA markers have been proved to be insufficient in authenticating many species of this group. Here, we investigated 101 complete plastomes from 23 DSFs, comprising 72 newly sequenced and 29 documented, which all exhibited well-conserved genomic organization and gene order. Plastome-wide comparison showed the co-occurrence of single nucleotide polymorphisms (SNPs) and insertions/deletions (indels), which can be explained by both the repeat-associated and indel-associated mutation hypotheses. Moreover, guanine-cytosine (GC) content was found to be negatively correlated with the three divergence variables (SNPs, indels and repeats), indicating that GC content may reflect the level of the local sequence divergence. Our species authentication analyses revealed that the relaxed filtering strategies of sequence alignment had no negative impact on species identification. By assessing the maximum likelihood (ML) trees inferred from different datasets, we found that the complete plastome and large single-copy (LSC) datasets both successfully identified all 23 DSFs with the maximum bootstrap values. However, owing to the high efficiency of LSC in species identification, we recommend using LSC for accurate authentication of DSFs.

Plastome phylogenomic study of Gentianeae (Gentianaceae): widespread gene tree discordance and its association with evolutionary rate heterogeneity of plastid genes

BACKGROUND

Plastome-scale data have been prevalent in reconstructing the plant Tree of Life. However, phylogenomic studies currently based on plastomes rely primarily on maximum likelihood inference of concatenated alignments of plastid genes, and thus phylogenetic discordance produced by individual plastid genes has generally been ignored. Moreover, structural and functional characteristics of plastomes indicate that plastid genes may not evolve as a single locus and are experiencing different evolutionary forces, yet the genetic characteristics of plastid genes within a lineage remain poorly studied.

RESULTS

We sequenced and annotated 10 plastome sequences of Gentianeae. Phylogenomic analyses yielded robust relationships among genera within Gentianeae. We detected great variation of gene tree topologies and revealed that more than half of the genes, including one (atpB) of the three widely used plastid markers (rbcL, atpB and matK) in phylogenetic inference of Gentianeae, are likely contributing to phylogenetic ambiguity of Gentianeae. Estimation of nucleotide substitution rates showed extensive rate heterogeneity among different plastid genes and among different functional groups of genes. Comparative analysis suggested that the ribosomal protein (RPL and RPS) genes and the RNA polymerase (RPO) genes have higher substitution rates and genetic variations among plastid genes in Gentianeae. Our study revealed that just one (matK) of the three (matK, ndhB and rbcL) widely used markers show high phylogenetic informativeness (PI) value. Due to the high PI and lowest gene-tree discordance, rpoC2 is advocated as a promising plastid DNA barcode for taxonomic studies of Gentianeae. Furthermore, our analyses revealed a positive correlation of evolutionary rates with genetic variation of plastid genes, but a negative correlation with gene-tree discordance under purifying selection.

CONCLUSIONS

Overall, our results demonstrate the heterogeneity of nucleotide substitution rates and genetic characteristics among plastid genes providing new insights into plastome evolution, while highlighting the necessity of considering gene-tree discordance into phylogenomic studies based on plastome-scale data.

Plastome Evolution in Dolomiaea (Asteraceae, Cardueae) Using Phylogenomic and Comparative Analyses

Dolomiaea is a medicinally important genus of Asteraceae endemic to alpine habitats of the Qinghai-Tibet Plateau (QTP) and adjacent areas. Despite significant medicinal value, genomic resources of Dolomiaea are still lacking, impeding our understanding of its evolutionary history. Here, we sequenced and annotated plastomes of four Dolomiaea species. All analyzed plastomes share the gene content and structure of most Asteraceae plastomes, indicating the conservation of plastome evolutionary history of Dolomiaea. Eight highly divergent regions (rps16-trnQ, trnC-petN, trnE-rpoB, trnT-trnL-trnF, psbE-petL, ndhF-rpl32-trnL, rps15-ycf1, and ycf1), along with a total of 51-61 simple sequence repeats (SSRs) were identified as valuable molecular markers for further species delimitation and population genetic studies. Phylogenetic analyses confirmed the evolutionary position of Dolomiaea as a clade within the subtribe Saussureinae, while revealing the discordance between the molecular phylogeny and morphological treatment. Our analysis also revealed that the plastid genes, rpoC2 and ycf1, which are rarely used in Asteraceae phylogenetic inference, exhibit great phylogenetic informativeness and promise in further phylogenetic studies of tribe Cardueae. Analysis for signatures of selection identified four genes that contain sites undergoing positive selection (atpA, ndhF, rbcL, and ycf4). These genes may play important roles in the adaptation of Dolomiaea to alpine environments. Our study constitutes the first investigation on the sequence and structural variation, phylogenetic utility and positive selection of plastomes of Dolomiaea, which will facilitate further studies of its taxonomy, evolution and conservation.

Structural variation of the complete chloroplast genome and plastid phylogenomics of the genus Asteropyrum (Ranunculaceae)

Two complete chloroplast genome sequences of Asteropyrum, as well as those of 25 other species from Ranunculaceae, were assembled using both Illumina and Sanger sequencing methods to address the structural variation of the cp genome and the controversial systematic position of the genus. Synteny and plastome structure were compared across the family. The cp genomes of the only two subspecies of Asteropyrum were found to be differentiated with marked sequence variation and different inverted repeat-single copy (IR-SC) borders. The plastomes of both subspecies contains 112 genes. However, the IR region of subspecies peltatum carries 27 genes, whereas that of subspecies cavaleriei has only 25 genes. Gene inversions, transpositions, and IR expansion-contraction were very commonly detected in Ranunculaceae. The plastome of Asteropyrum has the longest IR regions in the family, but has no gene inversions or transpositions. Non-coding regions of the cp genome were not ideal markers for inferring the generic relationships of the family, but they may be applied to interpret species relationship within the genus. Plastid phylogenomic analysis using complete cp genome with Bayesian method and partitioned modeling obtained a fully resolved phylogenetic framework for Ranunculaceae. Asteropyrum was detected to be sister to Caltha, and diverged early from subfamily Ranunculoideae.

Plastid phylogenomic insights into the evolution of the Caprifoliaceae s.l. (Dipsacales)

The family Caprifoliaceae s.l. is an asterid angiosperm clade of ca. 960 species, most of which are distributed in temperate regions of the northern hemisphere. Recent studies show that the family comprises seven major clades: Linnaeoideae, Zabelia, Morinoideae, Dipsacoideae, Valerianoideae, Caprifolioideae, and Diervilloideae. However, its phylogeny at the subfamily or genus level remains controversial, and the backbone relationships among subfamilies are incompletely resolved. In this study, we utilized complete plastome sequencing to resolve the relationships among the subfamilies of the Caprifoliaceae s.l. and clarify several long-standing controversies. We generated and analyzed plastomes of 48 accessions of Caprifoliaceae s.l., representing 44 species, six subfamilies and one genus. Combined with available Caprifoliaceae s.l. plastomes on GenBank and 12 outgroups, we analyzed a final dataset of 68 accessions. Genome structure was strongly conserved in general, although the boundaries between the Inverted Repeat were found to have contracted across Caprifoliaceae s.l. to exclude rpl2, rps19, and ycf1, all or parts of which are typically present in the IR of most angiosperms. The ndhF gene was found to have been inverted in all plastomes of Adoxaceae. Phylogenomic analyses of 68 complete plastomes yielded a highly supported topology that strongly supported the monophyly of Zabelia and its sister relationship to Morinoideae. Moreover, a clade of Valerianoideae + Dipsacoideae was recovered as sister to a clade of Linnaeoideae + Zabelia + Morinoideae clade, and Heptacodium was sister to remaining Caprifolioideae. The Diervilloideae and Caprifolioideae were successively sister to all other Caprifoliaceae s.l. Major lineages of Caprifoliaceae s.l. were estimated to have diverged from the Upper Cretaceous to the Eocene (50-100 Ma), whereas within-genus diversification was dated to the Oligocene and later, concomitant with global cooling and drying. Our results demonstrate the power of plastid phylogenomics in improving estimates of phylogeny among genera and subfamilies, and provide new insights into plastome evolution across Caprifoliaceae s.l.

Characterization of the complete chloroplast genome of the ornamental plant Osmanthus cooperi

Osmanthus cooperi is an evergreen ornamental plant belonging to the olive family. In this study, its complete chloroplast genome was assembled from the whole genome Illumina sequencing data. The circular genome is 155,262 bp long, and comprises a pair of inverted repeat regions (IRs, 25,685 bp each), a large single-copy region (LSC, 86,525 bp), and a small single-copy region (SSC, 17,367 bp). It encodes 132 genes, including 8 rRNA genes, 36 tRNAs genes, and 88 protein-coding genes. The GC content of O. cooperi cp genome is 37.8%. Phylogenetic analysis indicates that O. cooperi is close to O. fragrans in Oleaceae.

The first complete chloroplast genome sequence from Polemonium chinense (Polemoniaceae)

In this study, the complete chloroplast genome of Polemonium chinense was determined by next-generation sequencing technology. The plastome is 155,578 bp in length, with a large single-copy region (LSC) of 85,864 bp and a small single-copy region (SSC) of 18,246 bp, separated by a pair of inverted repeat regions (IRs) of 25,734 bp. It contains 135 genes, including 89 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 1 pseudogene. The overall GC content is 37.1%, while the corresponding values in the LSC, SSC and IR region are 35.1, 30.3, and 42.9%, respectively. The phylogenetic tree shows that Polemoniaceae is sister to Fouquieriaceae.

Plastome phylogenomics of Saussurea (Asteraceae: Cardueae)

BACKGROUND

Saussurea DC. is one of the largest and most morphologically heterogeneous genera in Asteraceae. The relationships within Saussurea have been poorly resolved, probably due an early, rapid radiation. To examine plastome evolution and resolve backbone relationships within Saussurea, we sequenced the complete plastomes of 17 species representing all four subgenera.

RESULTS

All Saussurea plastomes shared the gene content and structure of most Asteraceae plastomes. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection. Phylogenomic analyses of 20 Saussurea plastomes that alternatively included nucleotide or amino acid sequences of all protein-coding genes, vs. the nucleotide sequence of the entire plastome, supported the monophyly of Saussurea and identified three clades within it. Three of the four traditional subgenera were recovered as paraphyletic. Seven plastome regions were identified as containing the highest nucleotide variability.

CONCLUSIONS

Our analyses reveal both the structural conservatism and power of the plastome for resolving relationships in congeneric taxa. It is very likely that differences in topology among data sets is due primarily to differences in numbers of parsimony-informative characters. Our study demonstrates that the current taxonomy of Saussurea is likely based at least partly on convergent morphological character states. Greater taxon sampling will be necessary to explore character evolution and biogeography in the genus. Our results here provide helpful insight into which loci will provide the most phylogenetic signal in Saussurea and Cardueae.

Expected and unexpected evolution of plant RNA editing factors CLB19, CRR28 and RARE1: retention of CLB19 despite a phylogenetically deep loss of its two known editing targets in Poaceae

BACKGROUND

C-to-U RNA editing in mitochondria and chloroplasts and the nuclear-encoded, RNA-binding PPR proteins acting as editing factors present a wide field of co-evolution between the different genetic systems in a plant cell. Recent studies on chloroplast editing factors RARE1 and CRR28 addressing one or two chloroplast editing sites, respectively, found them strictly conserved among 65 flowering plants as long as one of their RNA editing targets remained present.

RESULTS

Extending the earlier sampling to 117 angiosperms with high-quality genome or transcriptome data, we find more evidence confirming previous conclusions but now also identify cases for expected evolutionary transition states such as retention of RARE1 despite loss of its editing target or the degeneration of CRR28 truncating its carboxyterminal DYW domain. The extended angiosperm set was now used to explore CLB19, an "E+"-type PPR editing factor targeting two chloroplast editing sites, rpoAeU200SF and clpPeU559HY, in Arabidopsis thaliana. We found CLB19 consistently conserved if one of the two targets was retained and three independent losses of CLB19 after elimination of both targets. The Ericales show independent regains of the ancestrally lost clpPeU559HY editing, further explaining why multiple-target editing factors are lost much more rarely than single target factors like RARE1. The retention of CLB19 despite loss of both editing targets in some Ericaceae, Apocynaceae and in Camptotheca (Nyssaceae) likely represents evolutionary transitions. However, the retention of CLB19 after a phylogenetic deep loss in the Poaceae rather suggests a yet unrecognized further editing target, for which we suggest editing event ndhAeU473SL.

CONCLUSION

Extending the scope of studies on plant organelle RNA editing to further taxa and additional nuclear cofactors reveals expected evolutionary transitions, strikingly different evolutionary dynamics for multiple-target editing factors like CLB19 and CRR28 and suggests additional functions for editing factor CLB19 among the Poaceae.