Chloroplast Genomes of Genus Tilia: Comparative Genomics and Molecular Evolution

Comparison of key floral components between two Tilia species and among different processing technologies as revealed by widely targeted metabolomic profiling

BACKGROUND

Tilia has a long history of cultivation and holds high ornamental and economic value. The volatile aroma compounds of Tilia flowers have significant characteristics that contribute to their ornamental appeal, and affect the flavor of floral tea. Here, widely targeted metabolomic analyses were conducted to investigate the aroma active compounds in the fresh inflorescences of Tilia cordata Mill. (Tc) and Tilia miqueliana Maxim. (Tm), and in samples prepared by freeze drying, air drying, and oven drying.

RESULTS

We identified 442 volatile organic compounds by headspace solid-phase microextraction gas chromatography-mass spectrometry. Terpenoids were the most abundant and diverse group, while heterocyclic compounds were the main contributors to the aroma profile. Notably, 46 aroma compounds were identified as primary contributors to the characteristic aroma of Tilia, including abhexon, 2-isobutyl-3-methoxypyrazine, (Z)-6-nonenal, methyl benzoate, (E)-2-hexenal, 1-hexanol, 2-thiophenemethanethiol, p-cymene, furaneol, and (Z)-4-heptenal. The concentration of volatile organic compounds was higher in Tc than in Tm, indicating a more pronounced aroma character of Tc. For both Tc and Tm, the aroma compounds were better preserved and present at higher concentrations in freeze-dried samples than in air-dried and oven-dried samples.

CONCLUSION

These results provide a foundation for further research on the molecular mechanisms of aroma formation in Tilia flowers and on aroma as a cue for insect pollination. Furthermore, the results provide a basis for the development and commercialization of Tilia floral teas and other related products. © 2025 Society of Chemical Industry.

Comparative analyses of chloroplast genomes of Theobroma cacao from northern Peru

Theobroma cacao is the most economically important species within the genus Theobroma. Despite its importance, the intraspecific relationships of this species has not been fully elucidated due to insufficient molecular information. To facilitate a better understanding of the intraspecific evolutionary relationships of T. cacao, Sequencing technology has been to decode the plastid genomes, with the objective of identify potential DNA barcode genetic markers, explore intraspecific relationships, and infer divergence times. The plastid genome of the seven cocoa genotypes analyzed in this study, exhibited a typical angiosperm genomic structure. However, the structure of each plastid genome reflects notable changes in each genotype; for example, the infA gene was present in all the analyzed samples, unlike in previously published cocoa plastid genomes, while the complete ycf1 gene sequence has potential for use as DNA Barcoding in T. cacao. The estimated age of the node connecting T. cacao and T. grandiflorum, which was 10.11 Ma, supports this indication. It can be inferred that T. cacao diverged at approximately 7.55 Ma, and it is highly likely that T. cacao populations diversified during the Pliocene or Miocene. Therefore, it is crucial to perform mitochondrial and nuclear-based analyses on a broader spectrum of cocoa samples to validate these evolutionary mechanisms, including genetic estimates and divergence. This approach enables a deeper understanding of the evolutionary relationships among cocoa.

The complete chloroplast genomes and comparative study of the two tung trees of Vernicia (Euphorbiaceae)

BACKGROUND

Vernicia montana and V. fordii are economically important woody oil species in the Euphorbiaceae that have great industrial oil and ornamental greening properties, however, the wild resources of Vernicia trees have been reduced because of their habitat destruction. Considering the diverse economic and ecological importance of Vernicia species, it is important to collect more molecular data to determine the genetic differences between V. montana and V. fordii.

RESULTS

We sequenced, assembled, and annotated the complete chloroplast (CP) genome of two tung trees based on the genome skimming approach. The whole CP genomes of V. montana and V. fordii were 163,518 bp and 161,495 bp in length, both including a pair of inverted repeats separated by a large single-copy and a small single-copy region. We detected a total number of 311 tandem repeats, 100 dispersed repeats, and 255 simple repeats from V. montana and V. fordii CP genomes. The mean value of nucleotide diversity between the two species was 0.0122, and the average Ka/Ks ratio across all coding genes was 0.3483. Comparative chloroplast genome analysis showed that the coding regions were more conserved than the non-coding regions. The phylogenetic relationships yielded by the complete genome sequences showed that V. montana was closely related to V. fordii and is considered as a sister group.

CONCLUSIONS

We sequenced, assembled, annotated, and analyzed the CP genome of two tung trees, which will be useful in investigating the conservation genetics and potential breeding applications of this oil shrub.

Comparative and phylogenetic analysis of the chloroplast genomes of four commonly used medicinal cultivars of Chrysanthemums morifolium

'Boju' and 'Huaiju' are cultivars of the Chrysanthemum (Chrysanthemum morifolium Ramat.) in the family Asteraceae, valued for their medicinal, tea, and ornamental properties, and valued by individuals. However, the yield and quality of medicinal chrysanthemums are limited by the characteristics of the germplasm resources, including the identification at the varieties and cultivation levels. Currently, research characterizing the chloroplast genomes of medicinal Chrysanthemum flowers is relatively limited. This study conducted chloroplast whole-genome sequencing on two cultivars of Chrysanthemum, 'Boju' and 'Huaiju', and compared them with the previously published chloroplast genomes of 'Hangbaiju' and 'Gongju'. The study analyzed the chloroplast genome structures of these four medicinal Chrysanthemums, identifying mutation hotspots and clarifying their phylogenetic relationships. The chloroplast genome sizes of four medicinal Chrysanthemum cultivation products ranged from 151,057 to 151,109 bp, with GC content ranging from 37.45% to 37.76%. A total of 134 genes were identified, including 89 protein-coding genes, 37 ribosomal RNA genes, and 8 transfer RNA genes. Comparative genomic analysis revealed 159 large repeat sequences, 276 simple sequence repeats, 1 gene, and 8 intergenic regions identified as highly variable regions. Nucleotide diversity (Pi) values were high (≥ 0.004) for the petN-psbM, trnR-UCU-trnT-GGU, trnT-GGU-psbD, ndhC-trnV-UCA, ycf1, ndhI-ndhG, trnL-UGA-rpl32, rpl32-ndhF, and ndhF-ycf1 fragments, aiding in variety identification. Phylogenetic analysis revealed consistent results between maximum likelihood and Bayesian inference trees, showing that the four medicinal Chrysanthemum cultivars, along with their wild counterparts and related species, evolved as a monophyletic group, forming a sister clade to Artemisia and Ajania. Among the six Chrysanthemum species, the wild Chrysanthemum diverged first (Posterior probability = 1, bootstrap = 1,000), followed by Ajania, while C. indicum and C. morifolium clustered together (Bootstrap = 100), indicating their closest genetic relationship. The chloroplast whole-genome data and characteristic information provided in this study can be used for variety identification, genetic conservation, and phylogenetic analysis within the family Asteraceae.

Insights into Aquilaria phylogenetics through comparative plastomic resources

The plastid is an essential organelle for its role in photosynthesis and energy production and its genomic information is always employed as important evolutionary markers to explore the relationship among species. Agarwood (Aquilaria), prized for its aromatic blend, finds extensive use in various cultures as incense and perfume. Despite its high economic importance, the phylogenetic status among Aquilaria based on plastomes remains ambiguous due to the lack of available plastomic resources. To bridge this knowledge gap, 22 Aquilaria plastomes were newly sequenced, similar variation patterns in this genus were determined, including a shared 16 bp extension of the rps19 gene and seven highly variable regions. The analysis highlighted the highest prevalence of the A/T motif among simple sequence repeats in these plastomes. Further phylogenetic analysis revealed Aquilaria's phylogenetic implications with an expanded dataset. This comprehensive plastomic resource not only enhances our understanding of Aquilaria evolution but also presents potential molecular markers for DNA barcoding.

Identification of key genes involved in lignin and flavonoid accumulation during Tilia tuan seed maturation

KEY MESSAGE

Transcriptomics and phenotypic data analysis identified 24 transcription factors (TFs) that play key roles in regulating the competitive accumulation of lignin and flavonoids. Tilia tuan Szyszyl. (T. tuan) is a timber tree species with important ecological and commercial value. However, its highly lignified pericarp results in a low seed germination rate and a long dormancy period. In addition, it is unknown whether there is an interaction between the biosynthesis of flavonoids and lignin as products of the phenylpropanoid pathway during seed development. To explore the molecular regulatory mechanism of lignin and flavonoid biosynthesis, T. tuan seeds were harvested at five stages (30, 60, 90, 120, and 150 days after pollination) for lignin and flavonoid analyses. The results showed that lignin accumulated rapidly in the early and middle stages (S1, S3, and S4), and rapid accumulation of flavonoids during the early and late stages (S1 and S5). High-throughput RNA sequencing analysis of developing seeds identified 50,553 transcripts, including 223 phenylpropanoid biosynthetic pathway genes involved in lignin accumulation grouped into 3 clusters, and 106 flavonoid biosynthetic pathway genes (FBPGs) grouped into 2 clusters. Subsequent WGCNA and time-ordered gene co-expression network (TO-GCN) analysis revealed that 24 TFs (e.g., TtARF2 and TtWRKY15) were involved in flavonoids and lignin biosynthesis regulation. The transcriptome data were validated by qRT-PCR to analyze the expression profiles of key enzyme-coding genes. This study revealed that there existed a competitive relationship between flavonoid and lignin biosynthesis pathway during the development of T. tuan seeds, that provide a foundation for the further exploration of molecular mechanisms underlying lignin and flavonoid accumulation in T. tuan seeds.

The complete chloroplast genome of Durio zibethinus L. cultivar Ri6 (Helicteroideae, Malvaceae)

Durian, a member of the Malvaceae family, is famous for its delicious fruits, which have strong scents and are rich in nutrients. In this study, we sequenced and characterized the complete chloroplast genome of Durio zibethinus L. 1774 cultivar Ri6, a popular durian cultivar in Vietnam, using the Illumina Hiseq platform. The results showed a circular chloroplast genome composed of a large single copy of 96,115 bp, a small single copy of 20,819 bp, and two inverted repeat regions of 24,185 bp. This genome consisted of 79 protein-coding genes, 30 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of this genome was 35.7%. Phylogenetic analysis inferred from 78 protein-coding regions revealed monophyly of Durio species and a close relationship between D. zibethinus cultivar Ri6 and cultivar Mongthong. This study provides essential information for further studies examining genetic population, breedings, and species identification among Durio taxa and cultivars.

Plastome diversity and evolution in mosses: Insights from structural characterization, comparative genomics, and phylogenetic analysis

Mosses play a significant role in ecology, evolution, and the economy. They belong to the nonvascular plant kingdom and are considered the closest living relatives of the first terrestrial plants. The circular chloroplast DNA molecules (plastomes) of mosses contain all the genetic information essential for chloroplast functions and represent the source of the evolutionary history of these organisms. This study comprehensively analyzed the plastomes of 47 moss species belonging to 14 orders, focusing on their size, GC content, gene loss, gene content, synteny, and evolution. The findings revealed great differences among plastome sizes, with Takakia lepidozioides (Takakiopsida) and Funaria hygrometrica (Funariales) having the largest and smallest plastomes, respectively. Moss plastomes included 69 to 89 protein-coding genes, 8 rRNA genes, and 34 to 42 tRNA genes, resulting in the total number of genes in a plastome ranging between 115 and 138. Various genes have been lost from the plastomes of different moss species, with Atrichum angustatum lacking the highest number of genes. This study also examined plastome synteny and moss evolution using comparative genomics and repeat sequence analysis. The results demonstrated that synteny and similarity levels varied across the 47 moss examined species, with some exhibiting structure similarity and others displaying structural inversions. Maximum likelihood and Bayesian approaches were used to construct a phylogenetic tree using 36 concatenated protein-coding genes, and the results revealed that the genera Sphagnum and Takakia are sister groups to the other mosses. Additionally, it was found that Tetraphidales, Polytrichales, Buxbaumiales, and Diphysciales are closely related. This research describes the evolutionary diversity of mosses and offers guidelines for future studies in this field. The findings also highlight the need for more investigations into the factors regulating plastome size variation in these plants.

Comparative analyses of complete chloroplast genomes reveal interspecific difference and intraspecific variation of Tripterygium genus

The genus Tripterygium was of great medicinal value and attracted much attention on the taxonomic study using morphological and molecular methods. In this study, we assembled 12 chloroplast genomes of Tripterygium to reveal interspecific difference and intraspecific variation. The sequence length (156,692-157,061 bp) and structure of Tripterygium were conserved. Comparative analyses presented abundant variable regions for further study. Meanwhile, we determined the ndhB gene under positive selection through adaptive evolution analysis. And the phylogenetic analyses based on 15 chloroplast genomes supported the monophyly of Tripterygium hypoglaucum and the potential sister relationship between Tripterygium wilfordii and Tripterygium regelii. Molecular dating analysis indicated that the divergence time within Tripterygium was approximately 5.99 Ma (95% HPD = 3.11-8.68 Ma). The results in our study provided new insights into the taxonomy, evolution process, and phylogenetic construction of Tripterygium using complete plastid genomes.