Comparative analysis of the complete chloroplast genome sequences of six species of Pulsatilla Miller, Ranunculaceae

Comparative analysis of chloroplast genomes of Pulsatilla species reveals evolutionary and taxonomic status of newly discovered endangered species Pulsatilla saxatilis

BACKGROUND

Pulsatilla saxatilis, a new species of the genus Pulsatilla has been discovered. The morphological information of this species has been well described, but its chloroplast genome characteristics and comparison with species of the same genus remain to be reported.

RESULTS

Our results showed that the total length of chloroplast (cp.) genome of P. saxatilis is 162,659 bp, with a GC content of 37.5%. The cp. genome contains 134 genes, including 90 known protein-coding genes, 36 tRNA genes, and 8 rRNA genes. P. saxatilis demonstrated similar characteristics to other species of genus Pulsatilla. Herein, we compared cp. genomes of 10 species, including P. saxatilis, and found that the cp. genomes of the genus Pulsatilla are extremely similar, with a length of 162,322-163,851 bp. Furthermore, The SSRs of Pulsatilla ranged from 10 to 22 bp in length. Among the four structural regions of the cp. genome, most long repeats and SSRs were detected in the LSC region, followed by that in the SSC region, and least in IRA/ IRB regions. The most common types of long repeats were forward and palindromic repeats, followed by reverse repeats, and only a few complementary repeats were found in 10 cp. genomes. We also analyzed nucleotide diversity and identified ccsA_ndhD, rps16_trnK-UUU, ccsA, and rbcL, which could be used as potential molecular markers for identification of Pulsatilla species. The results of the phylogenetic tree constructed by connecting the sequences of high variation regions were consistent with those of the cp. gene phylogenetic tree, and the species more closely related to P. saxatilis was identified as the P. campanella.

CONCLUSION

It was determined that the closest species to P. saxatilis is P. campanella, which is the same as the conclusion based on pollen grain characteristics, but different from the P. chinensis determined based on morphological characteristics. By revealing information on the chloroplast characteristics, development, and evolution of the cp. genome and the potential molecular markers, this study provides effective molecular data regarding the evolution, genetic diversity, and species identification of the genus Pulsatilla.

Chloroplast genome analyses of Caragana arborescens and Caragana opulens

BACKGROUND

Numerous species within the genus Caragana have high ecological and medicinal value. However, species identification based on morphological characteristics is quite complicated in the genus. To address this issue, we analyzed complete plastid genome data for the genus.

RESULTS

We obtained chloroplast genomes of two species, Caragana arborescens and Caragana opulens, using Illumina sequencing technology, with lengths of 129,473 bp and 132,815 bp, respectively. The absence of inverted repeat sequences in the two species indicated that they could be assigned to the inverted repeat-lacking clade (IRLC). The genomes included 111 distinct genes (4 rRNA genes, 31 tRNA genes, and 76 protein-coding genes). In addition, 16 genes containing introns were identified in the two genomes, the majority of which contained a single intron. Repeat analyses revealed 129 and 229 repeats in C. arborescens and C. opulens, respectively. C. arborescens and C. opulens genomes contained 277 and 265 simple sequence repeats, respectively. The two Caragana species exhibited similar codon usage patterns. rpl20-clpP, rps19-rpl2, and rpl23-ycf2 showed the highest nucleotide diversity (pi). In an analysis of sequence divergence, certain intergenic regions (matK-rbcL, psbM-petN, atpA-psbI, petA-psbL, psbE-petL, and rps7-rps12) were highly variable. A phylogenetic analysis showed that C. arborescens and C. opulens were related and clustered together with four other Caragana species. The genera Astragalus and Caragana were relatively closely related.

CONCLUSIONS

The present study provides valuable information about the chloroplast genomes of C. arborescens and C. opulens and lays a foundation for future phylogenetic research and molecular marker development.

The complete chloroplast genome of Pulsatilla chinensis f. alba D. K. Zang (Ranunculaceae, Pulsatilla Miller)

Pulsatilla chinensis f. alba D. K. Zang 1993 is a forma of Pulsatilla chinensis (Bge.) Regel, the root of P. chinensis is traditional Chinese medicine called Pulsatillae radix. The biggest difference between P. chinensis f. alba and P. chinensis is that P. chinensis f. alba sepals is white. The complete chloroplast genome of P. chinensis f. alba was sequenced using the Illumina NovaSeq platform for the first time. The lengths of the genome, large single-copy (LSC), small single-copy (SSC), two inverted repeats (IRs), and GC content were 163,654 bp, 82,355 bp, 19,069 bp, 31,115 bp, and 37.2%, respectively. It had 134 genes, including 90 protein-coding genes, 36 tRNA genes, and eight rRNA genes. The maximum-likelihood tree indicated that P. chinensis f. alba had a closer relationship with P. chinensis. This study would provide a theoretical basis for the further study of Pulsatilla plants genetics phylogenetic research.

Comparative analysis of medicinal plants Scutellaria baicalensis and common adulterants based on chloroplast genome sequencing

BACKGROUND

Scutellaria baicalensis Georgi has been extensively used as a medicinal herb in China for over 2000 years. They may be intentionally or inadvertently substituted or blended with comparable species in the local market, threatening clinical medication safety. Molecular markers are effective tools to prevent misidentification and eliminate doping and falsification among Scutellaria plants. This study screened four highly variable regions to identify Scutellaria and its adulterants. In addition, a phylogenetic analysis was performed using the complete cp genome combined with published Scutellaria species samples. Moreover, a comparative analysis of the cp genomes was conducted to investigate the cp genome evolution of S. baicalensis.

RESULTS

The complete cp genome of five species of Scutellaria was sequenced for the first time, and four previously published Scutellaria species were re-sequenced. They all exhibited a conserved quadripartite structure in their cp genomes, including two distinct regions, namely a small and large single copy region, respectively, and two inverted repeats encompassing the majority of ribosomal RNA genes. Furthermore, the nine species exhibited high conservation from aspects of the genome structure, codon usage, repeat sequences, and gene content. Four highly variable regions (matK-rps16, ndhC-trnV-UAC, psbE-petL, and rps16-trnQ-UUG) may function as potential molecular markers for differentiating S. baicalensis from its adulterants. Additionally, the monophyly of Scutellaria was ascertained and could be reclassified into two subgenera, subgenus Anaspis and subgenus Scutellaria, as evidenced by the phylogenetic analyses on sequences of cp genome and shared protein-coding sequences. According to the molecular clock analysis, it has been inferred that the divergence of Scutellaria occurred at approximately 4.0 Mya during the Pliocene Epoch.

CONCLUSION

Our study provides an invaluable theoretical basis for further Scutellaria species identification, phylogenetics, and evolution analysis.

The complete chloroplast genome of the halophyte flowering plant Suaeda monoica from Jeddah, Saudi Arabia

The complete chloroplast genome (plastome) of the annual flowering halophyte herb Suaeda monoica Forssk. ex J. F. Gmel. family (Amaranthaceae) that grows in Jeddah, Saudi Arabia, was identified for the first time in this study. Suaeda monoica is a medicinal plant species whose taxonomic classification remains controversial. Further, studying the species is useful for current conservation and management efforts. In the current study, the full chloroplast genome S. monoica was reassembled using whole-genome next-generation sequencing and compared with the previously published chloroplast genomes of Suaeda species. The chloroplast genome size of Suaeda monoica was 151,789 bp, with a single large copy of 83,404 bp, a small single copy of 18,007 bp and two inverted repeats regions of 25,189 bp. GC content in the whole genome was 36.4%. The cp genome included 87 genes that coded for proteins, 37 genes coding for tRNA, 8 genes coding for rRNA and one non-coding pseudogene. Five chloroplast genome features were compared between S. monoica and S. japonica, S. glauca, S. salsa, S. malacosperma and S. physophora. Among Suaeda genus and equal to most angiosperms chloroplast genomes, the RSCU values were conservative. Two pseudogenes (accD and ycf1), rpl16 intron and ndhF-rpl32 intergenic spacer, were highlighted as suitable DNA barcodes for different Suaeda species. Phylogenetic analyses show Suaeda cluster into three main groups; one in which S. monoica was closer to S. salsa. The obtained result provided valuable information on the characteristics of the S. monoica chloroplast genome and the phylogenetic relationships.

The economic burden of coronary heart disease in mainland China

OBJECTIVES

The aim of this study was to systematically evaluate the current economic burden of coronary heart disease (CHD) in mainland China and provide a reference for the formulation of policies to reduce the economic burden of CHD.

STUDY DESIGN

A systematic literature review was conducted of empirical studies on the economic burden of CHD over the past 20 years.

METHODS

PubMed, Web of Science, Embase, China Knowledge Resource Integrated Database and the WANFANG database were comprehensively searched for relevant articles published between 1 January 2000 and 22 December 2021. Content analysis was used to extract the data, and Stata 17.0 software was used for analysis. The median values were used to describe trends.

RESULTS

A total of 35 studies were included in this review. The annual median per-capita hospitalisation expense and the average expense per hospitalisation were $3544.40 ($891.64-$18,371.46) and $5407.34 ($1139.93-$8277.55), respectively. The median ratio on medical consumables expenses, drug expenses, medical examination expenses and treatment expenses were 41.59% (12.40%-63.73%), 26.90% (7.30%-60.00%), 9.45% (1.65%-33.40%) and 10.10% (2.36%-66.00%), respectively. The median per-capita hospitalisation expense in the eastern, central and western regions were $9374.45 ($2056.13-$18,371.46), $4751.5 ($2951.95-$8768.93) and $3251.25 ($891.64-$13,986.38), respectively. The median average expense per hospitalisation in the eastern and central regions were $6177.15 ($1679.15-$8277.55) and $1285.49 ($1239.93-$2197.36), respectively. The median average length of stay in the eastern, central and western regions were 9.3 days, 15.2 days and 16.1 days, respectively.

CONCLUSIONS

The economic burden of CHD is more severe in mainland China than in developed countries, especially in terms of the direct economic burden. In terms of the types of direct medical expenses, a proportion of medical examination expenses, treatment expenses and drug expenses were lowest in the eastern region, but medical consumables expenses were the highest in this region. This study provides guidance for the formulation of policies to reduce the economic burden of CHD in mainland China.

Elucidation of the mechanism of Yiqi Tongluo Granule against cerebral ischemia/reperfusion injury based on a combined strategy of network pharmacology, multi-omics and molecular biology

BACKGROUND

Ischemic stroke is caused by local lesions of the central nervous system and is a severe cerebrovascular disease. A traditional Chinese medicine, Yiqi Tongluo Granule (YQTL), shows valuable therapeutic effects. However, the substances and mechanisms remain unclear.

PURPOSE

We combined network pharmacology, multi-omics, and molecular biology to elucidate the mechanisms by which YQTL protects against CIRI.

STUDY DESIGN

We innovatively created a combined strategy of network pharmacology, transcriptomics, proteomics and molecular biology to study the active ingredients and mechanisms of YQTL. We performed a network pharmacology study of active ingredients absorbed by the brain to explore the targets, biological processes and pathways of YQTL against CIRI. We also conducted further mechanistic analyses at the gene and protein levels using transcriptomics, proteomics, and molecular biology techniques.

RESULTS

YQTL significantly decreased the infarction volume percentage and improved the neurological function of mice with CIRI, inhibited hippocampal neuronal death, and suppressed apoptosis. Fifteen active ingredients of YQTL were detected in the brains of rats. Network pharmacology combined with multi-omics revealed that the 15 ingredients regulated 19 pathways via 82 targets. Further analysis suggested that YQTL protected against CIRI via the PI3K-Akt signaling pathway, MAPK signaling pathway, and cAMP signaling pathway.

CONCLUSION

We confirmed that YQTL protected against CIRI by inhibiting nerve cell apoptosis enhanced by the PI3K-Akt signaling pathway.

Complete Chloroplast Genome Determination of Ranunculus sceleratus from Republic of Korea (Ranunculaceae) and Comparative Chloroplast Genomes of the Members of the Ranunculus Genus

Ranunculus sceleratus (family: Ranunculaceae) is a medicinally and economically important plant; however, gaps in taxonomic and species identification limit its practical applicability. This study aimed to sequence the chloroplast genome of R. sceleratus from Republic of Korea. Chloroplast sequences were compared and analyzed among Ranunculus species. The chloroplast genome was assembled from Illumina HiSeq 2500 sequencing raw data. The genome was 156,329 bp and had a typical quadripartite structure comprising a small single-copy region, a large single-copy region, and two inverted repeats. Fifty-three simple sequence repeats were identified in the four quadrant structural regions. The region between the ndhC and trnV-UAC genes could be useful as a genetic marker to distinguish between R. sceleratus populations from Republic of Korea and China. The Ranunculus species formed a single lineage. To differentiate between Ranunculus species, we identified 16 hotspot regions and confirmed their potential using specific barcodes based on phylogenetic tree and BLAST-based analyses. The ndhE, ndhF, rpl23, atpF, rps4, and rpoA genes had a high posterior probability of codon sites in positive selection, while the amino acid site varied between Ranunculus species and other genera. Comparison of the Ranunculus genomes provides useful information regarding species identification and evolution that could guide future phylogenetic analyses.

Complete chloroplast genome and phylogenetic analysis of Anemone shikokiana

BACKGROUND

Anemone shikokiana (Makino) Makino, disjunctly distributed in Shandong Peninsula of China and Shikoku Island of Japan, is a rare and endangered species. To provide genetic information and understand its phylogeny, we conducted research on the chloroplast (cp) genome of A. shikokiana.

METHODS AND RESULTS

The complete cp genome sequence of A. shikokiana was constructed in this study. The results showed that the cp genome of A. shikokiana has a typical quadripartite cyclic with a total length of 159,286 bp. In total, 111 unique genes were identified, including 78 protein-coding genes, 29 tRNA-coding genes and 4 rRNA-coding genes. A total of 37 long repeat sequences and 67 microsatellites were found in this cp genome. The cp genome of A. shikokiana was compared with eleven other Anemone cp genomes available from the Genbank database. We found some variations among the different genomes, especially in the LSC and SSC regions, and identified some regions as potential molecular markers such as ycf1, ndhE, ndhD, ndhF-trnL, ndhA and ndhF. The results of phylogenetic analysis suggested that A. narcissiflora was the closest relative of A. shikokiana.

CONCLUSIONS

The results filled the gap of cp genome sequence information of A. shikokiana, laying the foundation to explore the evolutionary relationships of A. shikokiana in future studies. It provided a valuable genetic resource for the molecular identification and phylogenetic study of Anemone.

Are the Organellar Genomes Useful for Fine Scale Population Structure Analysis of Endangered Plants?-A Case Study of Pulsatilla patens (L.) Mill

Pulsatilla patens is a rare and endangered species in Europe and its population resources have significantly decreased over the past decades. Previous genetic studies of this species made it possible to estimate the genetic diversity of the European population and to describe the structure of chloroplast and mitochondrial genomes. The main aim of these studies was to characterize the variability of chloroplast and mitochondrial genomes in more detail at the intra-population and inter-population levels. Our study presents new organelle genome reference sequences that allow the design of novel markers that can be the starting point for testing hypotheses, past and modern biogeography of rare and endangered species P. patens, and adaptive responses of this species to changing environments. The study included sixteen individuals from five populations located in Northeastern Poland. Comparative analysis of 16 P. patens plastomes from 5 populations enabled us to identify 160 point mutations, including 64 substitutions and 96 InDels. The most numerous detected SNPs and Indels (75%) were accumulated in three intergenic spacers: ndhD-ccsA, rps4-rps16, and trnL(UAG)-ndhF. The mitogenome dataset, which was more than twice as large as the plastome (331 kbp vs. 151 kbp), revealed eight times fewer SNPs (8 vs. 64) and six times fewer InDels (16 vs. 96). Both chloroplast and mitochondrial genome identified the same number of haplotypes-11 out of 16 individuals, but both organellar genomes slightly differ in haplotype clustering. Despite the much lower variation, mitogenomic data provide additional resolution in the haplotype detection of P. patens, enabling molecular identification of individuals, which were unrecognizable based on the plastome dataset.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

The complete chloroplast genome of Dictamnus dasycarpus Turcz. (Rutaceae)

Dictamnus dasycarpus Turcz. 1842 is a medicinal plant of China. Its dry root bark is called BAIXIANPI, which is a common traditional Chinese medicine. Here, we report the complete chloroplast genome of D. dasycarpus. The length of the genome, large single-copy (LSC), small single-copy (SSC), inverted repeat (IR), and GC content was 157,056 bp, 84,497 bp, 18,487 bp, 27,036 bp, and 38.5%, respectively. A total of 132 genes were annotated, including 87 protein coding, eight rRNA, and 37 tRNA genes. Interestingly, 15 genes contained single intron while two others contained two introns. The phylogenetic tree showed the two D. dasycarpus (D. albus) clustered in a clade, which was sister to clade formed by the species of Melicope, Tetradium, Phellodendron, and Zanthoxylum.

The Complete Chloroplast Genome of Endangered Species Stemona parviflora: Insight into the Phylogenetic Relationship and Conservation Implications

Stemona parviflora is an endangered species, narrowly endemic to Hainan and Southwest Guangdong. The taxonomic classification of S. parviflora remains controversial. Moreover, studying endangered species is helpful for current management and conservation. In this study, the first complete chloroplast genome of S. parviflora was assembled and compared with other Stemona species. The chloroplast genome size of S. parviflora was 154,552 bp, consisting of 87 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and one pseudogene. The ψycf1 gene was lost in the cp genome of S. sessilifolia, but it was detected in four other species of Stemona. The inverted repeats (IR) regions have a relatively lower length variation compared with the large single copy (LSC) and small single copy (SSC) regions. Long repeat sequences and simple sequence repeat (SSR) were detected, and most SSR were distributed in the LSC region. Codon usage bias analyses revealed that the RSCU value of the genus Stemona has almost no difference. As with most angiosperm chloroplast genomes, protein-coding regions were more conservative than the inter-gene spacer. Seven genes (atpI, ccsA, cemA, matK, ndhA, petA, and rpoC1) were detected under positive selection in different Stemona species, which may result from adaptive evolution to different habitats. Phylogenetic analyses show the Stemona cluster in two main groups; S. parviflora were closest to S. tuberosa. A highly suitable region of S. parviflora was simulated by Maxent in this study; it is worth noting that the whole territory of Taiwan has changed to a low fitness area and below in the 2050 s, which may not be suitable for the introduction and cultivation of S. parviflora. In addition, limited by the dispersal capacity of S. parviflora, it is necessary to carry out artificial grafts to expand the survival areas of S. parviflora. Our results provide valuable information on characteristics of the chloroplast genome, phylogenetic relationships, and potential distribution range of the endangered species S. parviflora.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Two Korean Endemic Clematis Chloroplast Genomes: Inversion, Reposition, Expansion of the Inverted Repeat Region, Phylogenetic Analysis, and Nucleotide Substitution Rates

Previous studies on the chloroplast genome in Clematis focused on the chloroplast structure within Anemoneae. The chloroplast genomes of Cleamtis were sequenced to provide information for studies on phylogeny and evolution. Two Korean endemic Clematis chloroplast genomes (Clematis brachyura and C. trichotoma) range from 159,170 to 159,532 bp, containing 134 identical genes. Comparing the coding and non-coding regions among 12 Clematis species revealed divergent sites, with carination occurring in the petD-rpoA region. Comparing other Clematis chloroplast genomes suggested that Clematis has two inversions (trnH-rps16 and rps4), reposition (trnL-ndhC), and inverted repeat (IR) region expansion. For phylogenetic analysis, 71 protein-coding genes were aligned from 36 Ranunculaceae chloroplast genomes. Anemoneae (Anemoclema, Pulsatilla, Anemone, and Clematis) clades were monophyletic and well-supported by the bootstrap value (100%). Based on 70 chloroplast protein-coding genes, we compared nonsynonymous (dN) and synonymous (dS) substitution rates among Clematis, Anemoneae (excluding Clematis), and other Ranunculaceae species. The average synonymoussubstitution rates (dS)of large single copy (LSC), small single copy (SSC), and IR genes in Anemoneae and Clematis were significantly higher than those of other Ranunculaceae species, but not the nonsynonymous substitution rates (dN). This study provides fundamental information on plastid genome evolution in the Ranunculaceae.

Characterization of the complete chloroplast genome of Populus deltoides Zhonglin 2025

The complete chloroplast genome of Populus deltoides was characterized by reference-based assembly using whole-genome sequencing data. The total chloroplast genome size of Populus deltoides included a pair of inverted repeat regions (IRs) of 27,649 bp each, a small single-copy region (SSC) of 16,563 bp, and large single-copy region (LSC) of 85,096 bp, which was 156,957 bp in length. A total of 109 genes were predicted from the chloroplast genome, including 83 protein-coding genes, 22 tRNA genes, and four rRNA genes. The GC content of chloroplast genome for Populus deltoides was 36.68%. The phylogenetic analysis based on the reported chloroplast genomes of Populus showed that the chloroplast of the Populus deltoides is most closely related to the Populus fremontii. The complete chloroplast genome of Populus deltoides provides new insights into Populus evolutionary and genomic studies.

The complete chloroplast genome of Ligusticum jeholense (Umbelliferae, Ligusticum L.)

Ligusticum jeholense is an important medicinal plant. Chloroplast genome information is helpful for the development of molecular markers and the study of plant phylogeny. In this study, we report the complete chloroplast genome sequence of L. jeholense. The genome sequence is 148,493 bp in size (GenBank accession number MN652885), with 37.25% GC contents. There are 127 genes in the genome, including 83 known protein-coding genes (PCGs), 36 transfer RNAs (tRNAs), and 8 ribosomal RNAs (rRNAs). The maximum-likelihood method are used to construct phylogenetic tree of 32 species. These data will provide certain theoretical basis for plant genetics research.