Comparative Analysis of Codon Bias in the Chloroplast Genomes of Theaceae Species

Theaceae species are dicotyledonous angiosperms with extremely high ornamental and economic value. The chloroplast genome is traditionally used to study species evolution, expression of chloroplast genes and chloroplast transformation. Codon usage bias (CUB) analysis is beneficial for investigations of evolutionary relationships and can be used to improve gene expression efficiency in genetic transformation research. However, there are relatively few systematic studies of the CUB in the chloroplast genomes of Theaceae species. In this study, CUB and nucleotide compositions parameters were determined by the scripts written in the Perl language, CodonW 1.4.2, CU.Win2000, RStudio and SPSS 23.0. The chloroplast genome data of 40 Theaceae species were obtained to analyse the codon usage (CU) characteristics of the coding regions and the influence of the source of variation on CUB. To explore the relationship between the CUB and gene expression levels in these 40 Theaceae plastomes, the synonymous codon usage order (SCUO) and measure independent of length and composition (MILC) values were determined. Finally, phylogenetic analysis revealed the genetic evolutionary relationships among these Theaceae species. Our results showed that based on the chloroplast genomes of these 40 Theaceae species, the CUB was for codons containing A/T bases and those that ended with A/T bases. Moreover, there was great commonality in the CUB of the Theaceae species according to comparative analysis of relative synonymous codon usage (RSCU) and relative frequency of synonymous codon (RFSC): these species had 29 identical codons with bias (RSCU > 1), and there were 19 identical high-frequency codons. The CUB of Theaceae species is mainly affected by natural selection. The SCUO value of the 40 Theaceae species was 0.23 or 0.24, and the chloroplast gene expression level was moderate, according to MILC values. Additionally, we observed a positive correlation between the SCUO and MILC values, which indicated that CUB might affect gene expression. Furthermore, the phylogenetic analysis showed that the evolutionary relationships in these 40 Theaceae species were relatively conserved. A systematic study on the CUB and expression of Theaceae species provides further evidence for their evolution and phylogeny.

Introduction

Theaceae, a characteristic family of subtropical evergreen broad-leaved forests, is widely distributed in subtropical regions (Zhang et al., 2021). Different Camellia (Camellia sinensis L. Kuntze.) varieties are popular in many countries around the world, including India, Vietnam and Japan (Teixeira and Sousa, 2021; Wang et al., 2021). Camellia originated in Southeast Asia, and 80% of the species of this genus come from China. Theaceae species have high ornamental value and high economic value, such as medicine, food and drink (Yuan et al., 2021). Tea is water soaked with Camellia leaves, known as “oriental magic water”, which is one of three most famous drinks in the world (Balentine et al., 1997). In addition, the seeds of Camellia are rich in oil that can be extracted and used for various purposes, such as edible and industrial oils.

The chloroplast genome usually consists of double-stranded closed circular DNA molecules, but the chloroplast DNA of a few plants is linear. The chloroplast genome is relatively small, generally ranging from 107 to 218 kb in length, accounting for approximately 10%–20% of total plant DNA (Yao et al., 2015). The chloroplast not only has a relatively stable evolutionary rate but also has a relatively slow rate of nucleotide substitution. Thus, the investigation of the chloroplast genome is of great significance for revealing the evolutionary rules, genetic relationships, expression of chloroplast genes and chloroplast transformation among plants (Olmstead and Palmer, 1994).

Codons are the link between amino acids, proteins, and genetic materials in organisms. Meanwhile, codons play an irreplaceable role in transmitting genetic information. There are 64 codons that collectively code for 20 amino acids. Except for methionine (Met, coded by AUG) and tryptophan (Trp, coded by UGG), which have unique codons, the remaining amino acids generally have 2–6 corresponding codons; this phenomenon is known as codon degeneracy (Ikemura, 1985). Codons coding the same amino acid are called synonymous codons, and their frequencies of use in the translation process are different. This phenomenon of unbalanced codon use is called codon usage bias (CUB; Sharp and Li., 1987; Wang et al., 2020). The study of CUB is beneficial for exploring genetic evolution, understanding gene expression characteristics and providing guidance for molecular breeding.

In recent years, the rapid development of high-throughput sequencing technology has led to increasingly enriched chloroplast genome databases, laying a solid data foundation for plant genome structure and phylogenetic analyses. Thus, the chloroplast genome has been studied more extensively and in-depth. Moreover, related research on CUB in the chloroplast genome has been carried out in a variety of plants, including Amaranthus hybridus (Chaney et al., 2016), Camellia oleifera (Zhang et al., 2019), Medicago truncatula (Guo et al., 2020), Nicotiana tabacum (Arévalo-Gallegos et al., 2020); meanwhile, Magnoliaceae (Chang et al., 2006), Bambusoideae (Zhang et al., 2011), Euphorbiaceae plants (Wang et al., 2020) and Oryza species (Chakraborty et al., 2020). There are numerous reports on the chloroplast genome of the Theaceae species, but relatively rare on its CUB.

From the perspective of the coding sequence, we systematically analysed the codon usage (CU) characteristics of the chloroplast genomes of Theaceae species, the source of variation, the high-frequency (HF) codons, and the relationship between the CUB and gene expression levels based on measures independent of length and composition (MILC) and synonymous codon usage order (SCUO) parameters, and phylogenetic analysis was used to investigate the genetic relationships between the plants. This study revealed the relationships between the CUB and gene expression frequency, and provided a reference direction for further genetic evolution research and transgenic research in these 40 Theaceae species.

Materials and Methods

Source of Material

The coding regions (coding sequences, CDSs) of the chloroplast genomes of 40 Theaceae species, namely, 33 species of Camellia, one species of Hartia (Hartia laotica) and six species of Stewartia, were downloaded and collected from the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov/) on May 17, 2021. Detailed information on the chloroplast genomes of the 40 Theaceae species is provided in Table 1. To ensure the accuracy of the experimental results, the scripts written in the Perl language (Wang et al., 2020) were used to process the original CDSs of the 40 chloroplast genomes, and the selected sequences were required to meet the following conditions: 1) the CDS is composed of only A, T, G, and C bases; 2) the sequence length is a multiple of three; 3) the sequence starts with the start codon ATG and ends with the stop codon TAA, TAG or TGA; 4) there is no stop codon in the middle of the sequence (means to exclude prematurely terminated coding sequences); 5) the length of the CDS is ≥ 300 bp (He et al., 2016; Li et al., 2016; Wang et al., 2020).

TABLE 1

TABLE 1. Chloroplast genome information of 40 Theaceae species.

Coding Usage Index

GC1, GC2, and GC3 respectively refer to the GC content of the first, second, and third base positions, respectively, of the chloroplast codons in the 40 Theaceae species. The GC1, GC2, and GC3 contents and the average GC content among the three sites of the 40 Theaceae species were calculated by using the scripts written in the Perl language (Wang et al., 2020). The CU index of the chloroplast CDSs was analysed by the CodonW 1.4.2 program, including T3s, C3s, A3s, and G3s, representing the frequency of T/C/A/G usage at the third base of the codon. The GC content was calculated as a percentage of the total base content, and the RSCU was determined.

Analysis of Relative Synonymous Codon Usage (RSCU), Relative Frequency of Synonymous Codon (RFSC) and High-Frequency (HF)

RSCU is a statistical index used to weigh the relative frequency of each synonymous codon, which describes the ratio of the actual observed usage frequency of a certain codon to the expected usage frequency of the codon. Codons with an RSCU value > 1 have a positive bias, as they are used more frequently than any other synonymous codons; while an RSCU value < 1 represents negative bias, indicating that these codons are used less frequently than any other synonymous codons; an RSCU value = 1 indicates no bias, which means that synonymous codons are used equally and that CU is random (Sharp and Li, 1986). For the calculation formula of RSCU, please refer to the following formula (Eq. 1) by Sharp and Li. (1986):

$RSC{U}_{ij}=\frac{x_{ij}}{{{\displaystyle \sum }}_j^{n_i}{x}_{ij}}{n}_i(1)$

A Toolkit for Biologists integrating various biological data-handling tools (TBtools) software was used to create a heatmap of the average RSCU values (Chen et al., 2020).

RFSC refers to the ratio of the number of codons observed in the plants to the total number of synonymous codons. The calculation formula for the frequency of use under the influence of gene length and amino acid abundance is detailed in the study of Zhou et al. (2007a) and Zhou et al. (2007b). The RFSC value can be calculated according to the actual number of codons using the CodonW 1.4.2 program (Sharp and Li, 1986). RFSC can be calculated using the following Eq. 2 (Sharp and Li, 1986):

$RFSC=\frac{x_{ij}}{{{\displaystyle \sum }}_j^{n_i}{x}_{ij}}(2)$

When the RFSC value of a codon is greater than 60% or more than 1.5 times the average frequency of its synonymous codons, this codon is called an HF codon (Zhou et al., 2007a; Zhou et al., 2007b).

Source Analysis of Codon Variation

Correspondence analysis (COA) multivariate statistical methods are often used to study the variance and gene distribution of RSCU, and each CDS was represented as a 59-dimensional vector (excluding the codon ATG encoding Met, the codon TGG encoding Trp, and the three stop codons TAA, TAG, and TGA) (Sharp and Li, 1986). The COA in CodonW 1.4.2 is generally used to study the CU variation in the genome. To reduce the influence of amino acid composition on codon utilization, each dimension corresponds to the RSCU value of each relevant codon, eliminating the influence of the stop codons TAA, TAG, and TGA.

The effective number of codons (ENc) evaluates codon use bias at the genomic level; the lower the ENc value is, the more significant the codon use preference is (Wright, 1990; Zhang et al., 2020).

$ENc=2+S+\frac{29}{S^2+{(1-S)}^2}(3)$

In the ENc-GC3s plot, GC3s and ENc values were used as horizontal and vertical coordinates, respectively, to analyse the relationship between the base composition and CUB (Wang and Hickey, 2007).

Parity rule 2 (PR2) plot analysis is widely used to explain whether there is an impact on nucleotide composition due to mutation pressure and selection pressure in DNA double strands, which includes G3/(G3+C3) and A3/(A3+T3) as the horizontal and vertical coordinates, respectively. The positions of A = T and G = C (PR2) are the centres, and the coordinates are (0.5, 0.5; Meyer, 2021). The centre position of the graph indicates that there is no deviation of mutation and selectivity (replacement rate) between two complementary DNA strands. When the proportions of G and C (or A and T) are similar, the gene CUB is affected only by mutation pressure; if the composition ratio is too large, it indicates that the CUB originates from natural selection and other factors (Sueoka, 1999).

Neutral analysis is a method for quantitatively analyzing the influence of directed mutation pressure and natural selection on CUB. We examined the neutral plot and calculated the GC12 and GC3 values of the genome using the scripts written in the Perl language (Wang et al., 2020). When the slope of the regression curve goes to 0 and there is no significant correlation between GC12 and GC3, it is fully influenced by natural selection. In contrast, when the slope is close to or equal to 1 and the correlation is significant, mutation pressure is speculated to have a major effect on genes (Sueoka, 1988).

SCUO Analysis

To better illustrate the relationship between the CUB and gene expression, we used the CU.Win2000 program to calculate the SCUO of the chloroplast sequences of the 40 Theaceae species (Wan et al., 2006). The SCUO value represents the CUB of synonymous codons relative to the whole sequence. SCUO values range between 0 and 1. The larger the value is, the stronger the CUB is. The calculation formula is explained in detail in the study by Wan et al. (2006).

MILC Analysis

MILC is a measure of the gene expression level used to quantify the difference in CU between a gene and some expected distribution of codons. The codon distribution can either be calculated from the background nucleotide composition or derived from a single gene or a gene group (Supek and Vlahovicek, 2005). For each chloroplast genome sequence, CU analysis in R uses the package coRdon 1.13.0. First, we need to count occurrences of each codon in each sequence. The default parameter values are set during the calculation, and the MILC value of each sequence in the set is calculated for the average CU deviation of the entire sample (https://bioconductor.org/packages/devel/bioc/vignettes/coRdon/). A high MILC value indicates a high gene expression standard, and a low MILC value indicates a low gene expression standard (Dorić and Bilela, 2014).

Phylogenetic Analysis of the 40 Theaceae Species

Two codons (AUG and UGG) and three stop codons (UAA, UAG, and UGA) encoding only one amino acid were removed from the chloroplast genomes of the 40 Theaceae species. The RSCU values of the remaining 59 synonymous codons were used as variables to calculate the distance between genes, namely, the squared Euclidian distance (Hou et al., 2016). Then, to explore the relationships between the CUB of the Theaceae species and their phylogenetic evolution in association with species relatedness, the intergroup linkage method in the SPSS 23.0 software (Pallant, 2013) was used to draw a cluster diagram for cluster analysis. The chloroplast genome sequences of 40 Theaceae species were aligned using the Muitiple Sequence Alignment Program (MAFFT v.7.4.0) (Katoh and Standley, 2013). The Maximum Likelihood phylogenetic tree of 40 Theaceae species was constructed with RAxML v.8.2.12 (Stamatakis, 2014) with 1000 bootstrap replicates, referring to the method of Liu et al. (2021). Other parameters were set by default in the software.

Results

Nucleotide Composition Analysis

The chloroplast coding genes of the 40 Theaceae species had very similar base compositions in the whole coding region (Table 2). The GC1, GC2 and GC3 contents were 35.84%-39.61%, 35.95%-39.88% and 35.61%-39.65%, respectively, and the mean values of the three were 37.14%-38.29%, which indicates an overall preference for codons containing A and T and ending in A/T. Therefore, from the perspective of base composition, the CUB of the chloroplast genomes of these 40 Theaceae species is extremely similar.

TABLE 2

TABLE 2. Summary of the SCUO, MILC, GC, GC1, GC2, and GC3 values of the chloroplast genomes of 40 Theaceae species.

RSCU, RFSC Analysis and HF Codon Screening Reveal Biased Use of A/T

As shown in Figure 1, the number of preferred codons (RSCU > 1) in the 40 chloroplast genomes was 29, among which, 28 (96.55%) ended in A or U (excluding UUG, 12 ended in A and 16 ended in U), the probability of the preferred codon ending in A/T, was much higher than that ending in G/C. In all the Theaceae species, AGA, encoding arginine (Arg), was strongly preferred (RSCU > 1.7), followed by GAU, encoding aspartic acid (Asp). The RFSC value was calculated by using the actual number of codons used by the CodonW 1.4.2 program (Supplementary Table S1). In summary, the codons in the chloroplast genomes of the 40 Theaceae species all showed a preference for ending in A/U. Furthermore, in Figure 1, the colours of the 40 Theaceae species for the same codon were extremely similar, indicating that the RSCU values for a given codon in the chloroplast genomes of the Theaceae species were similar. The HF codons in the chloroplast genomes of the 40 Theaceae species are shown in Table 3. The HF codons encoding 19 amino acids (Ala, Cys, Asp, Glu, Phe, Gly, His, Ile, Lys, Leu, Asn, Pro, Gln, Arg, Ser, Thr, Val, Tyr, and TER) are highly similar, with only a few codons showing some differences (Phe, Gln, Val, TER). There were 18 common HF codons (GCU, UGU, GAU, GAA, UUU, GGA, CAU, AUU, AAA, UUA, AAU, CCU, CAA, AGA, UCU, ACU, GUA, and UAU) in the chloroplast genomes, and only the codons UUC, CAG, GUU and UAA were different.

FIGURE 1

FIGURE 1. The RSCU values of the chloroplast genomes in 40 Theaceae species. A gradient from blue to red indicates that the average RSCU value of the codon is from low to high.

TABLE 3

TABLE 3. Screening of high-frequency codons in the chloroplast genomes of the 40 Theaceae species.

Analysis of Codon Variation Suggests Natural Selection Dominates

In the COA of the 40 Theaceae species (Figure 2), the points representing different genes were separated from each other in all plants of Camellia. In the chloroplast genomes of the 40 Theaceae species, axis 1 and axis 2 accounted for 9.32%-10.73% and 8.05%-9.95% of the synonymous CU variation, respectively. Genes with different GC contents are marked with different colours to reveal the influence of GC content on CUB. The GC content of all the genes in the 40 Theaceae species was less than 45%. This high consistency indicates that the sources of CU variation in different Theaceae species are extremely similar.

FIGURE 2

FIGURE 2. COA analysis of the chloroplast genomes in 40 Theaceae species. In all the Theaceae species, the dots of different genes are separated from each other in the figure. The numbers (A–AN) of the 40 species are shown in Table 1.

In the ENc-GC3s map of the 40 Theaceae species (Figure 3), all the genes were scattered in small clusters on the left side. Based on observation and analysis, the chloroplast genomes of the 40 Theaceae species showed scattered small clusters overall. Most of the genes were distributed far away from the standard curve and others were distributed in the appendix of the standard curve. This shows that the CUB of the chloroplast genomes is affected not only by mutation pressure but also by natural selection and other factors.

FIGURE 3

FIGURE 3. ENc-GC3s plot of the chloroplast genomes in 40 Theaceae species. The chloroplast genomes of the 40 Theaceae species are generally scattered in small clusters, and the genes are distributed on the left side of the figure. The numbers (A–AN) of the 40 species are shown in Table 1.

Natural selection may affect the use of bases A and T together with G and C. In this study, we found that the chloroplast genes of the 40 Theaceae species were distributed unevenly in the four quadrants of the PR2 plot (Figure 4), and they were distributed mainly in the regions of G3/(G3+C3) > 0.5 and A3/(A3+T3) < 0.5, indicating that the use of the third codon was unbalanced. Figure 5 shows an important positive correlation between GC12 and GC3 of the chloroplast genomes (P ≤ 0.01), which indicates that all the Theaceae species have been subject to direct mutation pressure on chloroplast gene CU. In the neutral analysis chart, the distribution ranges of GC12 and GC3 were relatively narrow. GC12 was distributed between 0.32 and 0.57, and GC3 was distributed between 0.14 and 0.36. Most points were distributed in small clusters, indicating that both natural selection and mutation pressure affect CUB.

FIGURE 4

FIGURE 4. PR2 plot of the chloroplast genomes in 40 Theaceae species. Points are randomly distributed in the four quadrants, and are mainly distributed in the regions of G3/(G3+C3) > 0.5 and A3/(A3+T3) < 0.5. The numbers (A–AN) of the 40 species are shown in Table 1.

FIGURE 5

FIGURE 5. Neutral analysis of the chloroplast genomes in 40 Theaceae species. The distribution range of GC12 and GC3 is relatively narrow, GC12 is distributed between 0.32 and 0.57, GC3 is distributed between 0.14 and 0.36, most of the points are distributed in small clusters. The numbers (A–AN) of the 40 species are shown in Table 1.

The neutral analysis chart can further reveal the strength of mutation pressure and natural selection. When the slope of the regression line is close to 0, natural selection plays a leading role, and when the slope of the regression line is close to 1, rapid mutation pressure is assumed to play a main role. Among six Camellia species, namely, Camellia amplexicaulis, Camellia brevistyla, Camellia fascicularis, Camellia gauchowensis, Camellia granthamiana, and Camellia leptophylla, the absolute value of the angle at which the regression line slopes was slightly greater than 0.1, including −0.1006, 0.1218, 0.1112, 0.1484, −0.1024 and 0.1024, respectively, indicating that the mutation pressure accounted for approximately 10% of the variation, which means that natural selection plays an important role; additionally, among the other 34 Theaceae species, natural selection had an effect on CUB as high as 90%.

SCUO Analysis

SCUO was used as an important parameter to analyse CUB in the chloroplast genomes of the 40 Theaceae species. The average SCUO value of each chloroplast genome is shown in Table 2. The average SCUO of the chloroplast genome of 14 species of Theaceae species was 0.23, while the average SCUO of the other 26 species of Theaceae species was 0.24, both of which are lower than the median value of 0.5, which indicates that most genes in the chloroplast genome of Theaceae species have a low CUB. In summary, from the perspective of SCUO, the CUB of the chloroplast genomes of the 40 Theaceae species is very similar.

CUB Affects Gene Expression

The MILC value reflects the level of gene expression. The greater the MILC number is, the higher the gene expression levels is. We calculated the MILC value of chloroplast genes of the 40 Theaceae species (Supplementary Table S2), and the results showed that 31 Theaceae species had a MILC value of 0.55, and the other 9 species had a MILC value of 0.56, which revealed that the expression level of most of the chloroplast genes was moderate. Additionally, we used correlation analysis to explore the relationship between the gene expression level and the CUB. We observed a significant positive correlation between the SCUO and MILC values (Table 4), indicating that the gene expression levels may be affected by the CUB.

TABLE 4

TABLE 4. Correlation between the SCUO and MILC values in the different chloroplast genomes of 40 Theaceae species.

Phylogenetic Analysis of the 40 Theaceae Species

The phylogenetic tree constructed based on the cluster analysis of the RSCU values and the ML analysis of the chloroplast genomes has a highly similar topology. The greater the squared Euclidean distance between plants is, the greater the difference in their CUB. According to the RSCU calculations, the squared Euclidean distance between the chloroplast genes of the 40 Theaceae species was extremely small (all less than 0.12). At the same time, the clustering results divided the plants into two major branches: the first branch contained ten species of Camellia, all belonging to the genus Camellia; the second branch contained 30 of Theaceae, including 23 species of Camellia, one species of Hartia (Hartia laotica) and six species of Stewartia species (Figure 6). On the first branch, Camellia brevistyla, Camellia gauchowensis and Camellia yuhsienensis formed a small group, while the remaining plants belonged to another group: Camellia azalea, Camellia petelotii, Camellia grandibracteata, Camellia pubicosta, Camellia leptophylla, Camellia Nitidissima, and Camellia reticulata. On the second branch, two species of Camellia (Camellia monadelpha and Camelliaserrata), one species of Hartia (Hartia laotica) and six species of Stewartia formed a small group, with the two species of Camellia (Camellia renshanxiae and Camellia sasanqua) forming their own small group. The other 21 species of Camellia belonged to another small group. The results showed that the 21 species of Camellia in the second clade were more closely related, while the 10 species in the first clade were more distantly related.

FIGURE 6

FIGURE 6. Phylogenetic analysis of the chloroplast genomes in 40 Theaceae species. The clustering results and the phylogenetic tree are divided into two branches. The topology of the phylogenetic tree is similar to the clustering results based on the codon RSCU values of the chloroplast genome to a certain extent.

The results of the constructed phylogenetic tree based on the genomes are samely divided into two major branches: the first branch contains 33 species of Camellia and all belong to the genus Camellia, and the second branch contains the remaining seven species of Camellia, including one species of Hartia (Hartia laotica) and six species of Stewartia. However, the six species of Stewartia belong to a small branch in phylogenetic tree according to the RSCU values. Overall, the topological structure of the phylogenetic tree and the clustering results based on the RSCU values of the chloroplast genome are similar.

Discussion

Codons are an important core element that connects amino acids, proteins and genetic materials in living organisms (Sundararajan et al., 2016). Studies of the preferred usage bias of codons provide reliable information for the study of protein expression and its corresponding functions (Wang et al., 2018). Many biological factors have been demonstrated to influence CUB, including gene length, tRNA quantity, gene expression level, chain-specific mutation bias, and GC composition (Wan et al., 2004; Meyer, 2021).

Therefore, this study systematically analysed the CUB characteristics, CUB variation sources and phylogenetic relationships of the chloroplast genomes of 40 Theaceae species. As claimed by the results of RSCU, RFSC analysis and HF codon screening, the codons tended to end with A/T, which is consistent with the assumption that “higher plants tend to use codons ending with A/T” (Campbell and Gowri, 1990; Wang et al., 2020). According to Sharp et al. (1988) proposed that “organisms with close genetic relationships have very similar codon usage bias”, the CUB in the chloroplast genomes of the 40 Theaceae species are extremely similar based on the RSCU values and comparative analysis of CU. In addition, a CUB analysis of 18 rice chloroplast genomes revealed that their CUB was low and AT base preference was strong; different rice varieties have different CU patterns that are mainly subjected to mutation pressure and natural selection (Chakraborty et al., 2020). We also analysed the preferred codons and HF codons of the 40 Theaceae species (Campbell and Gowri, 1990); specifically, there was very high similarity in the use of codons among the 40 Theaceae species. Eighteen codons were both preferred codons and HF codons: GCU, UGU, GAU, GAA, UUU, GGA, CAU, AUU, AAA, UUA, AAU, CCU, CAA, AGA, UCU, ACU, GUA, and UAU. Compared to the use of one type of codon (preferred codons or HF codons), the use of the two types (preferred codons and HF codons) may have a higher reference value for the optimization of codons in terms of heterologous expression.

Regarding the source of codon variation, our results further suggested that both mutation pressure and natural selection affect CUB in the chloroplast genomes of the 40 Theaceae species. ENc-GC3s and PR2 plot analysis showed that the CUB of the 40 Theaceae species was influenced by mutation pressure and natural selection, and the similar results were found in other tea plants (Rawal et al., 2020; Rawal et al., 2021). In addition, a natural selection study of the chloroplast genomes of 29 species of Magnoliaceae systematically compared their CU characteristics and elucidated their evolutionary relationships (Chang et al., 2006). Meanwhile, Zhang et al. (2011) get the similar conclusion on 3 species in the bamboo subfamily. Comparisons of plant chloroplast genomes have revealed the main reason for their CUB. Furthermore, the quantitative results of neutral analysis showed that the CUB of the 40 Theaceae species were mainly affected by natural selection, which was consistent with the influencing factors of CUB in the chloroplast genomes of Magnoliaceae (Chang et al., 2006) and Euphorbiaceae plant ( Wang et al., 2020). In summary, natural selection plays a leading role while mutation pressure plays a secondary role in determining the CU of chloroplast genes in 40 Theaceae species.

In our study, the calculated SCUO values were lower than 0.25, which means that most genes in the chloroplast-encoding genome of Theaceae plastomes have a low bias. This result is consistent with the result of chloroplast genes in Oryza species (Chakraborty et al., 2020). In Camellia amplexicaulis, Camellia fascicularis, Camellia granthamiana, Camellia japonica S288C, Camellia kissii, Camellia ptilophylla, Camellia reticulata, Hartia laotica, Stewartia micrantha, Stewartia monadelpha, Stewartia obovata, Stewartia serrata, Stewartia sinii, and Stewartia villosa, the average value of SCUO was 0.23, which reflects the low codon bias. In the results of 18 Oryza species studies, the SCUO values ranged from 0.24 to 0.27, indicating that the CUB was low and varied in most of the rice species. Likewise, the 40 Theaceae species had low chloroplast CUB. We regarded MILC as a measure of gene expression in our current study and observed high MILC values. The MILC value of 31 of the Theaceae species was 0.55, and the value of the other nine kinds of MILC was 0.56. In the related research on Oryza species, the expression level of each chloroplast gene of 18 rice species based on the MILC value was calculated to be in the range of 0.74 to 0.78, which means that most of the chloroplast genes are expressed at high levels in rice and that there are differences in gene expression levels (Chakraborty et al., 2020). The results showed that, from the perspective of MILC, there may be certain differences in gene expression between dicotyledonous plants and monocotyledonous plants of chloroplast genomes. Karl Pearson’s correlation analysis between SCUO and MILC was used to explore the effect of CUB (Sareen et al., 2018). Based on the gene expression level, we observed a significant positive between SCUO and MILC values (P < 0.01; Table 4), indicating that the gene expression levels may be affected by the CUB. All in all, the higher the CUB is, the higher the gene expression level is (Li et al., 2016). In related reports by Zhao et al. (2016) and Liu et al. (2020), there is a close correlation between gene expression frequency and CUB, that is, the stronger the CUB, the higher the gene expression frequency. The efficiency of gene expression can be improved by codon optimization. For example, enhancing the expression of viral genes in mammalian cell lines by complete resynthesis of the gene (Gustafsson et al., 2004). Based on the common use characteristics of the codons identified in the above analysis, the expression efficiency of some exogenous genes in the Theaceae species can be improved, and the functional genes regulating the growth and development of the Theaceae species can be targeted (Peng et al., 2021).

There are many species of plants in Theaceae, and the interspecific relationship of this economically important plant is still a very controversial issue (Vijayan et al., 2009; Yang et al., 2013). The two classification approaches proposed by Chang (1998) and Ming (2000) diverge in many respects, especially with regard to the division of subgenera, sections and species. So far, it is uncertain which system has most accurately described the phylogenetic relationships within the Theaceae family. Therefore, it is necessary to look for other evidence to reconstruct the establishment of the Theaceae classification system. The squared Euclidean distances between plants within the genera Camellia and Stewartia were smaller (≤ 0.1), indicating that the CUB of the chloroplast genes were very similar (Wang, 2007). Based on the cluster analysis map constructed by RSCU values, most species of the same genus are clustered together, indicating that closely related species tend to have the same codon bias pattern. This is consistent with the conclusion in Liu et al. (2005) research that the closer the genetic relationship between species, the more similar the CU patterns. In addition, the clustering results based on RSCU values were similar to the topological structure of the phylogenetic tree constructed based on chloroplast genome sequences. In the topology of the phylogenetic tree, Camellia chekiangoleosa, Camellia crapnelliana, Camellia cuspidat, Camellia granthamiana, Camellia gymnogyna, Camellia impressinervis, Camellia pitardii, Camellia ptilophylla, Camellia renshanxiangiae, Camellia yunnanensis belongs to the first clade, which has the same opinion with Li et al. (2021). And Camellia pitardii was clustered with Camellia crapnelliana, Camellia cuspidate, Camellia impressinervis, and Camellia yunnanensis, indicating that Camellia pitardii might belong to subgenus Camellia, and it means they’re very closely related, which is similar to Li et al. (2021). In addition, there are some inconsistencies between the phylogenetic tree of this study and others’ Theaceae phylogeny. Li et al. (2021) is mainly focus on the phylogeny of Camellia, Camellia renshanxiangiae and Camellia sasanqua belong to the genus Camellia. However, in cluster analysis based on RSCU values, the two species of Camellia grouped into the Stewartia-Hartia clade. This provides a clearer direction for our subsequent research in the further. Therefore, the clustering method based on RSCU values can be used as a supplement to the molecular taxonomy and phylogenetic studies of Theaceae species. Furthermore, based on the genetic relationships revealed by cluster analysis, it is expected that in the future, transgenic technology can be used to conduct molecular experiments on genes underlying important traits in Theaceae species in order to facilitate in-depth studies on the cultivation of Theaceae species.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author Contributions

ZZ and ZJW designed the research. ZJW, QC, YW, ML, and CW analyzed the data and performed the experiments. ZJW, QC, YW, and ML wrote the paper. ZXW, CJ, CX, and HW conducted manuscript editing. ZZ revised the manuscript. All authors contributed to the article and approved to the published version of the article.

Funding

This work was supported by the National Natural Science Foundation of China (32072624; 32172626), the Open Fund of State Key Laboratory of Tea Plant Biology and Utilization (SKLTOF20200129), the Support program for Outstanding Young Talents in Colleges and universities of Anhui Province (gxyq2020040), the Special Scientific Research Project of Graduate Tutor of Hefei Normal University (DSKY17).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.824610/full#supplementary-material

References

Arévalo-Gallegos, S., Varela-Rodríguez, H., Lugo-Aguilar, H., Siqueiros-Cendón, T. S., Iglesias-Figueroa, B. F., Espinoza-Sánchez, E. A., et al. (2020). Transient Expression of a Green Fluorescent Protein in Tobacco and Maize Chloroplast. Electron. J. Biotechnol. 45, 1–9. doi:10.1016/j.ejbt.2020.01.008

CrossRef Full Text | Google Scholar

Balentine, D. A., Wiseman, S. A., and Bouwens, L. (1997). The Chemistry of Tea Flavonoids. Crit. Rev. Food Sci. 37 (8), 693–704. doi:10.1080/10408399709527797

CrossRef Full Text | Google Scholar

Campbell, W. H., and Gowri, G. (1990). Codon Usage in Higher Plants, Green Algae, and Cyanobacteria. Plant Physiol. 92 (1), 1–11. doi:10.1104/pp.92.1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakraborty, S., Yengkhom, S., and Uddin, A. (2020). Analysis of Codon Usage Bias of Chloroplast Genes in Oryza Species. Planta 252, 67. doi:10.1007/s00425-020-03470-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaney, L., Mangelson, R., Ramaraj, T., Jellen, E. N., and Maughan, P. J. (2016). The Complete Chloroplast Genome Sequences for Four Amaranthus Species (Amaranthaceae). Appl. Plant Sci. 4 (9), 1600063. doi:10.3732/apps.1600063

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, C., Lin, H., Lin, I., Chow, T., Chen, H., Chen, W., et al. (2006). The Chloroplast Genome of Phalaenopsis aphrodite (Orchidaceae): Comparative Analysis of Evolutionary Rate with that of Grasses and Its Phylogenetic Implications. Mol. Biol. Evol. 23, 279–291. doi:10.1093/molbev/msj029

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, H. T. (1998). Flora of Reipublicae Popularis Sinicae. Delectis Florae Republicae Popularis Sinicae, Agendae Academiae Sinicae Edita. Beijing, China: Science Press, Vol. 49 (3).

Google Scholar

Chen, C., Chen, H., Zhang, Y., Thomas, H. R., and Xia, R. (2020). TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 13 (8), 1194–1202. doi:10.1016/j.molp.2020.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dorić, S., and Bilela, L. L. (2014). “Comparison of Codon Usage in Mitochondrial Genomes of Rhinolophid and Hipposiderid Bats,” in The 5th International Symposium on Sustainable Development, 39–45.

Google Scholar

Guo, S., Zhang, X., Bai, Q., Zhao, W., Fang, Y., Zhou, S., et al. (2020). Cloning and Functional Analysis of Dwarf Gene Mini Plant 1 (MNP1) in Medicago truncatula. Ijms 21 (14), 4968. doi:10.3390/ijms21144968

PubMed Abstract | CrossRef Full Text | Google Scholar

Gustafsson, C., Govindarajan, S., and Minshull, J. (2004). Codon Bias and Heterologous Protein Expression. Trends Biotechnol. 22 (7), 346–353. doi:10.1016/j.tibtech.2004.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

He, B., Dong, H., Jiang, C., Cao, F., Tao, S., and Xu, L. (2016). Analysis of Codon Usage Patterns in Ginkgo biloba Reveals Codon Usage Tendency from A/U-Ending to G/C-Ending. Sci. Rep. 6, 35927. doi:10.1038/srep35927

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, H., Wang, S., and Wu, Z. (2016). Cluster Analysis and Codon Usage Bias of psb A Genes from Prokaryotic Microalgae. J. Agric. Sci. 30 (7), 1296–1307. doi:10.11869/j.issn.100-8551

CrossRef Full Text | Google Scholar

Ikemura, T. (1985). Codon Usage and tRNA Content in Unicellular and Multicellular Organisms. Mol. Biol. Evol. 2 (1), 13–34. doi:10.1093/oxfordjournals.molbev.a040335

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., and Standley, D. M. (2013). Mafft Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 30, 772–780. doi:10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Hu, Y., He, M., Zhang, B., Wu, W., Cai, P., et al. (2021). Comparative Chloroplast Genomes: Insights into the Evolution of the Chloroplast Genome of Camellia sinensis and the Phylogeny of Camellia. BMC Genom. 22, 138. doi:10.1186/s12864-021-07427-2

CrossRef Full Text | Google Scholar

Li, N., Sun, M., Jiang, Z., Shu, H., and Zhang, S. (2016). Genome-Wide Analysis of the Synonymous Codon Usage Patterns in Apple. J. Integr. Agric. 15, 983–991. doi:10.1016/s2095-3119(16)61333-3

CrossRef Full Text | Google Scholar

Liu, L., Du, X., Guo, C., and Li, D. (2021). Resolving Robust Phylogenetic Relationships of Core Brassicaceae Using Genome Skimming Data. J. Syst. Evol. 59 (3), 442–453. doi:10.1111/jse.12666

CrossRef Full Text | Google Scholar

Liu, Q., Dou, S., Ji, Z., and Xue, Q. (2005). Synonymous Codon Usage and Gene Function are Strongly Related in Oryza sativa. Biosystems 80 (2), 123–131. doi:10.1016/j.biosystems.2004.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Li, Y., Ji, K., Zhu, J., Ling, P., Zhou, T., et al. (2020). Genome-Wide Codon Usage Pattern Analysis Reveals the Correlation Between Codon Usage Bias and Gene Expression in Cuscuta australis. Genomics 112, 2695–2702. doi:10.1016/j.ygeno.2020.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, M. M. (2021). Revisiting the Relationships Between Genomic G + C Content, RNA Secondary Structures, and Optimal Growth Temperature. J. Mol. Evol. 89, 165–171. doi:10.1007/s00239-020-09974-w

CrossRef Full Text | Google Scholar

Ming, T. (2000). Monograph of the genus ‘Camellia’. Kunming, China: Yunnan Science and Technology Press.

Google Scholar

Olmstead, R. G., and Palmer, J. D. (1994). Chloroplast DNA Systematics: A Review of Methods and Data Analysis. Am. J. Bot. 81, 1205–1224. doi:10.1002/j.1537-2197.1994.tb15615.x

CrossRef Full Text | Google Scholar

Pallant, J. (2013). SPSS Survival Manual: A Step guide to data analisis using SPSS for Windows. Aust. N. Z. J. Public Health 37 (6), 597–598. doi:10.1111/1753-6405.12166

CrossRef Full Text | Google Scholar

Peng, J., Zhao, Y., Dong, M., Liu, S., Hu, Z., Zhong, X., et al. (2021). Exploring the Evolutionary Characteristics Between Cultivated Tea and Its Wild Relatives Using Complete Chloroplast Genomes. BMC Ecol. Evol. 21, 71. doi:10.1186/s12862-021-01800-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Rawal, H. C., Borchetia, S., Bera, B., Soundararajan, S., Ilango, R. V. J., Barooah, A. K., et al. (2021). Comparative Analysis of Chloroplast Genomes Indicated Different Origin for Indian Tea (Camellia assamica cv TV1) as Compared to Chinese Tea. Sci. Rep. 11, 110. doi:10.1038/s41598-020-80431-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Rawal, H. C., Kumar, P. M., Bera, B., Singh, N. K., and Mondal, T. K. (2020). Decoding and Analysis of Organelle Genomes of Indian Tea (Camellia assamica) for Phylogenetic Confirmation. Genomics 112, 659–668. doi:10.1016/j.ygeno.2019.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Sareen, K., Bhalja, B. R., Srivastava, S., Swarnkar, Y., and Maheshwari, R. P. (2018). Islanding Detection Technique Based on Karl Pearson's Coefficient of Correlation for Distribution Network with High Penetration of Distributed Generations. JEEPS 19 (3), 20170232. doi:10.1515/ijeeps-2017-0232

CrossRef Full Text | Google Scholar

Sharp, P. M., Elizabeth, C., Desmond, G. H., Shields, D. C., Kenneth, H. W., and Wright, F. (1988). Codon Usage Patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens; A Review of the Considerable Within-Species Diversity. Nucleic Acids Res. 16 (17), 8207–8011. doi:10.1093/nar/16.17.8207

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharp, P. M., and Li, W. (1986). An Evolutionary Perspective on Synonymous Codon Usage in Unicellular Organisms. J. Mol. Evol. 24, 28–38. doi:10.1007/BF02099948

CrossRef Full Text | Google Scholar

Sharp, P. M., and Li, W. (1987). The Codon Adaptation Index-a Measure of Directional Synonymous Codon Usage Bias, and Its Potential Applications. Nucleic Acids Res. 15 (3), 1281–1295. doi:10.1093/nar/15.3.1281

PubMed Abstract | CrossRef Full Text | Google Scholar

Stamatakis, A. (2014). Raxml Version 8: A Tool for Phylogenetic Analysis and Post-Analysis of Large Phylogenies. Bioinformatics 30, 1312–1313. doi:10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Sueoka, N. (1988). Directional Mutation Pressure and Neutral Molecular Evolution. Proc. Natl. Acad. Sci. USA 85 (8), 2653–2657. doi:10.1073/pnas.85.8.2653

PubMed Abstract | CrossRef Full Text | Google Scholar

Sueoka, N. (1999). Translation-Coupled Violation of Parity Rule 2 in Human Genes is Not the Cause of Heterogeneity of the DNA G+C Content of Third Codon Position. Gene 238 (1), 53–58. doi:10.1016/S0378-1119(99)00320-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundararajan, A., Dukowic-Schulze, S., Kwicklis, M., Engstrom, K., Garcia, N., Oviedo, O. J., et al. (2016). Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays. Front. Plant Sci. 7, 1433. doi:10.3389/fpls.2016.01433

PubMed Abstract | CrossRef Full Text | Google Scholar

Supek, F., and Vlahoviček, K. (2005). Comparison of Codon Usage Measures and Their Applicability in Prediction of Microbial Gene Expressivity. BMC Bioinform. 6, 182. doi:10.1186/1471-2105-6-182

CrossRef Full Text | Google Scholar

Teixeira, A. M., and Sousa, C. (2021). A Review on the Biological Activity of Camellia species. Molecules 26 (8), 2178. doi:10.3390/molecules26082178

PubMed Abstract | CrossRef Full Text | Google Scholar

Vijayan, K., Zhang, W., and Tsou, C. (2009). Molecular Taxonomy of Camellia (Theaceae) Inferred from nrITS Sequences. Am. J. Bot. 96, 1348–1360. doi:10.3732/ajb.0800205

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, X., Xu, D., Kleinhofs, A., and Zhou, J. (2004). Quantitative Relationship Between Synonymous Codon Usage Bias and GC Composition Across Unicellular Genomes. BMC Ecol. Evol. 4, 19. doi:10.1186/1471-2148-4-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, X., Zhou, J., and Xu, D. (2006). CodonO: A New Informatics Method for Measuring Synonymous Codon Usage Bias Within and Across Genomes. Int. J. Gen. Syst. 35, 109–125. doi:10.1080/03081070500502967

CrossRef Full Text | Google Scholar

Wang, H., and Hickey, D. A. (2007). Rapid Divergence of Codon Usage Patterns Within the Rice Genome. BMC Ecol. Evol. 7 (Suppl 1), S6. doi:10.1186/1471-2148-7-S1-S6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Xing, H., Yuan, Y., Wang, X., Saeed, M., Tao, J., et al. (2018). Genome-Wide Analysis of Codon Usage Bias in Four Sequenced Cotton Species. PLoS One 13, e0194372. doi:10.1371/journal.pone.0194372

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W. (2007). “Squared Euclidean Distance Based Convolutive Non-Negative Matrix Factorization with Multiplicative Learning Rules for Audio Pattern Separation,” in 2007 IEEE International Symposium on Signal Processing and Information Technology, 347–352. doi:10.1109/isspit.2007.4458186

CrossRef Full Text | Google Scholar

Wang, Y., Zhuang, H., Shen, Y., Wang, Y., and Wang, Z. (2021). The Dataset of Camellia Cultivars Names in the World. Biodivers. Data J. 9, e61646. doi:10.3897/bdj.9.e61646

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Xu, B., Li, B., Zhou, Q., Wang, G., Jiang, X., et al. (2020). Comparative Analysis of Codon Usage Patterns in Chloroplast Genomes of Six Euphorbiaceae Species. Peer J. 8, e8251. doi:10.7717/peerj.8251

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, F. (1990). The 'Effective Number of Codons' Used in a Gene. Gene 87 (1), 23–29. doi:10.1016/0378-1119(90)90491-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Yang, S., Li, H., Yang, J., and Li, D. (2013). Comparative Chloroplast Genomes of Camellia Species. PLoS One 8, e73053. doi:10.1371/journal.pone.0073053

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, X., Tang, P., Li, Z., Li, D., Liu, Y., and Huang, H. (2015). The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis. PLoS One 10, e0129347. doi:10.1371/journal.pone.0129347

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, B., Zhou, S., Liu, C., Zhang, S., Li, J., and Liu, A. (2021). Status of the Application of Exogenous Enzyme Technology for the Development of Natural Plant Resources. Bioprocess Biosyst. Eng. 44, 429–442. doi:10.1007/s00449-020-02463-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Xu, B., Li, J., Zhao, Z., and Liu, Z. (2020). Transit from Autotrophism to Heterotrophism: Sequence Variation and Evolution of Chloroplast Genomes in Orobanchaceae Species. Front. Genet. 11, 542017. doi:10.3389/fgene.2020.542017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Zhao, Y., Yang, G., Peng, J., Chen, S., and Xu, Z. (2019). Determination of the Evolutionary Pressure on Camellia oleifera on Hainan Island Using the Complete Chloroplast Genome Sequence. Peer J. 7, e7210. doi:10.7717/peerj.7210

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Ma, P., and Li, D. (2011). High-Throughput Sequencing of Six Bamboo Chloroplast Genomes: Phylogenetic Implications for Temperate Woody Bamboos (Poaceae: Bambusoideae). PLoS One 6, e20596. doi:10.1371/journal.pone.0020596

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Yang, Y., Ye, P., Ma, J., and Ye, H. (2021). The Complete Chloroplast Genome Sequence of Camellia rostrata S. X. Yang & S. F. Chai (Theaceae), a Critically Endangered Yellow Camellia from Southwest China. Mitochondrial DNA B Resour. 6, 2427–2429. doi:10.1080/23802359.2021.1955028

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Zheng, H., Xu, A., Yan, D., Jiang, Z., Qi, Q., et al. (2016). Analysis of Codon Usage Bias of Envelope Glycoprotein Genes in Nuclear Polyhedrosis Virus (NPV) and Its Relation to Evolution. BMC Genom. 17, 677. doi:10.1186/s12864-016-3021-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, M., Tong, C., and Shi, J. (2007a). Analysis of Codon Usage Between Different Poplar Species. J. Genet. Genom. 34, 555–561. doi:10.1016/s1673-8527(07)60061-7

CrossRef Full Text | Google Scholar

Zhou, M., Tong, C. F., and Shi, J. S. (2007b). A Preliminary Analysis of Synonymous Codon Usage in Poplar Species. Zhi Wu Sheng Li Yu Fen Zi Sheng Wu Xue Xue Bao 33 (4), 285–293. doi:10.2471/BLT.13.118778

PubMed Abstract | CrossRef Full Text | Google Scholar