Skip to main content

Codon usage pattern of the ancestor of green plants revealed through Rhodophyta

Abstract

Rhodophyta are among the closest known relatives of green plants. Studying the codons of their genomes can help us understand the codon usage pattern and characteristics of the ancestor of green plants. By studying the codon usage pattern of all available red algae, it was found that although there are some differences among species, high-bias genes in most red algae prefer codons ending with GC. Correlation analysis, Nc-GC3s plots, parity rule 2 plots, neutrality plot analysis, differential protein region analysis and comparison of the nucleotide content of introns and flanking sequences showed that the bias phenomenon is likely to be influenced by local mutation pressure and natural selection, the latter of which is the dominant factor in terms of translation accuracy and efficiency. It is worth noting that selection on translation accuracy could even be detected in the low-bias genes of individual species. In addition, we identified 15 common optimal codons in seven red algae except for G. sulphuraria for the first time, most of which were found to be complementary and bound to the tRNA genes with the highest copy number. Interestingly, tRNA modification was found for the highly degenerate amino acids of all multicellular red algae and individual unicellular red algae, which indicates that highly biased genes tend to use modified tRNA in translation. Our research not only lays a foundation for exploring the characteristics of codon usage of the red algae as green plant ancestors, but will also facilitate the design and performance of transgenic work in some economic red algae in the future.

Peer Review reports

Introduction

Plants in a broad sense, known as Archaeplastida, are a large and diverse group, including green plants (Viridiplantae), red algae (Rhodophyta) and glaucophytes (Glaucophyta), all of which can carry out photosynthesis in their chloroplasts evolved from cyanobacteria through the primary endosymbiosis process [1,2,3,4,5,6]. Some fossil [7,8,9,10], molecular clock [11,12,13,14,15] and phylogenetic studies [16,17,18,19] have shown that Rhodophyta appeared approximately 1.5 billion years ago and may be the closest ancestors of Viridiplantae (Fig. 1).

Fig. 1
figure 1

Simplified phylogeny of plant evolution. Mya represents millions of years ago. Approximate divergence dates are taken from Parfrey, et al. [14]. Cyanidioschyzon merolae and Chondrus crispus are representatives of unicellular and multicellular red algae respectively

Codons are the bridge between genes and proteins. In a genome, the phenomenon of unequal usage of synonymous codons is called codon bias, which exists widely in different species [20]. The level of codon usage bias in a gene is commonly determined using the value of Nc (effective number of codons) ranging from 20 to 61, which reflects the degree of preference for codon usage [21]. A large number of researchers have proposed that natural selection, mutation pressure and even genetic drift can affect codon usage patterns in species [22,23,24]. In some species, natural selection is the main factor driving codon usage, and if they have a set of optimal codons for translation, it will provide them with a selective advantage over species with other codons under selection pressure [24, 25]. Because the frequency of codon usage matched the cellular tRNA population, Ikemura [25] suggested that the optimal codon matched the most abundant isoaccepting tRNA and the optimal codons were generally frequently used in highly expressed genes [25, 26].

Suzuki, et al. [27] studied the codon adaptation of plastid genomes of 12 red algae, and found that the relative strength of selection varied greatly in the Rhodophyta, with some showing strong bias, and others showing weaker bias. Li, et al. [28] analyzed the codon usage pattern of the chloroplast genome of Porphyra umbilicalis and found that it was mainly affected by natural selection, mutation pressure and nucleotide composition. Lee, et al. [29] studied expressed sequence tags (ESTs) to analyze codon usage in Griffithsia okiensis, Chondrus crispus and Porphyra yezoensis and found that they had higher GC-content bias and that most of the optimal codons of Griffithsia okiensis as well as Chondrus crispus ended with C, but that of Porphyra yezoensis ended with G. However, there have been few studies on the codon patterns related to the nuclear genome of red algae.

In this study, we performed a detailed analysis of codon usage patterns, influencing factors, optimal codons and corresponding tRNA genes of all nuclear coding genes in four unicellular Rhodophyta and four multicellular Rhodophyta, which might not only lay a foundation for exploring the characteristics of codon usage of red algae as green plant ancestors, but also will aid in the design and performance of transgenic work in some economic red algae to maximize corresponding protein yield in the future.

Results

Direction and strength of codon usage bias in the rhodophyta

The size of the genome and the number of CDSs extracted from eight red algae (C. merolae, G. sulphuraria, C. yangmingshanensis, C. crispus, G. chorda, G. domingensis, P. purpureum and P. umbilicalis) are presented in Table 1, showing that they are 16–23 Mb and 4803–9898 in unicellular red algae and 78–105 Mb and 9603–13360 in multicellular species, respectively. The effective number of codons (Nc) can be used to evaluate the degree of codon usage bias. According to Table 1, the Nc values ranged from 40 to 56 in these species, in which P. umbilicalis shows the strongest codon usage bias, followed by G. sulphuraria, P. purpureum and C. crispus. In addition, the codon usage biases of the other four species were similar.

Table 1 Whole genome and codon usage statistics in the coding sequence of eight red algae

The G or C nucleotide content of the third synonymous position (GC3s) can indicate the bias direction of the synonymous codon of a species or a gene. According to Table 1, G. sulphuraria is the only species that exhibits poor GC bias. The GC3s content of the CDSs was only approximately 0.313 ± 0.049(mean ± SD). The GC3s contents of G. chorda and G. domingensis were medium at approximately 0.5, specifically 0.504 ± 0.087 and 0.505 ± 0.099, respectively. Contrary to G. sulphuraria, there was strong GC bias across the other five red algae, in which P. umbilicalis had the highest GC3s (0.816 ± 0.084), followed by P. purpureum (0.635 ± 0.073).

In short, P. purpureum and P. umbilicalis show a strong preference for codons ending with GC. C. yangmingshanensis, C. merolae and C. crispus show a preference for codons ending with GC, but this preference is not strong. G. chorda and G. domingensis almost show a neutral preference for four codons. G. sulphuraria is different from the other species and strongly prefers codons ending with AT. Therefore, there are different codon usage biases among different red algae.

The role of mutation pressure in codon usage

To confirm the influence of mutation pressure on the codon usage pattern in plant ancestor species, we calculated the GC3s of coding sequences, as well as the GC of introns and flanking DNA in high-bias, medium-biased and low-bias gene categories from eight red algae genomes. The results are shown in Fig. 2 and Supplementary Table S1. High-bias genes in the seven red algae preferred codons ending with GC, except for those in G. sulphuraria. In addition, the GC3s of the coding sequences and GC of flanking DNA and introns decreased with the decreasing of preference level in the seven species except for G. sulphuraria. In other words, both GC3s and GC were the highest in the high-bias genes, while they were the lowest in the low-bias genes. However, the GC patter of introns was relatively complex. For four multicellular species, the variation trend of the GC content of introns was consistent with that of the codon bias. However, it was different from that in the four unicellular species. For example, the GC content of the introns in P. purpureum was the highest among the high-bias genes, but the lowest in G. sulphuraria and remained relatively stable in C. yangmingshanensis. By pairwise comparison in three different bias categories (high bias and mid bias, mid bias and low bias) for eight species, it was found that the difference was very significant in most comparisons (Supplementary Table S1). In particular, there was no intron in the high-bias genes and only one intron in the low-bias genes in C. merolae. Therefore, it was impossible to carry out t tests between different categories in C. merolae (Supplementary Table S1). Additionally, except for G. domingensis, the GC3s in the mid-bias genes was similar that in the low bias genes, which was the same as the content of GC in the noncoding regions among the three bias categories (Supplementary Table S1). However, in the high-bias gene category, GC3s was markedly different from the content of GC, which indicates that natural selection plays a dominant role in codon usage across the high-bias gene category. Although the GC3s of CDSs was not related to the GC content of introns in unicellular red algae, a strong correlation was found in multicellular red algae (Supplementary Figure S1 and Supplementary Table S2). Additionally, the stop codons showed a preference for UAA over both UAG and UGA with G in highly biased genes of all seven species (Supplementary Table S3), which proves that the GC bias was not caused mainly by mutation pressure.

Fig. 2
figure 2

Comparison of the average GC3s for coding sequences and GC content for noncoding DNA in red algae. The line chart represents the standard deviation of each category value. Red dots represent the value of GC3s shown on the left y-axis (from left to right: high-bias genes, mid-bias genes and low-bias genes). The GC of noncoding is displayed on the right y-axis (flanking DNA: green dots; intron DNA: blue dots. From left to right: high-bias genes, mid-bias genes and low-bias genes). A: unicellular red algal group; B: multicellular red algal group

Overall, mutation pressure exerts a limited influence, so it is supposed that natural selection strongly affects the codon usage pattern of red algae by driving the high GC content of highly biased genes, especially in unicellular species.

Nc-GC3s plot

The GC3s content of each gene in the genome of red algae was taken as the abscissa, and Nc was taken as the ordinate (Fig. 3). In the plot, only some of the actual observed values of each red algal gene are very close to the expected values, indicating that the codon use of these genes is almost entirely caused by mutation, and most of the gene loci are located at a position deviating from the expected curve, indicating that the codon preference of red algal genes is more affected by selection.

Fig. 3
figure 3

The relationship between the effective number of codons (Nc) and GC content at the third synonymous codon position (GC3s) in eight red algae. In each plot, the curved line represents the expected position of genes under a neutral mutation model. A: unicellular red algal group; B: multicellular red algal group

Neutrality plot analysis

To determine the main influencing factors of codon usage preference of the whole genome of red algae, neutrality plot analysis of the coding sequences of eight red algae was performed, as shown in Fig. 4. The slope of the plots shows the proportion of mutation pressure influencing the codon usage pattern. According to Fig. 4, mutation pressure accounts for less than 50% across all red algae, except G. domingensis(54.1%). For the four unicellular red algae, the slope is not more than 5% (C. merolae: 2.97%, C. yangmingshanensis: 2.75%, G. sulphuraria: 4.27% and P. purpureum: 1.27%) which indicates that natural selection and other factors play a dominant role in the formation of codon usage patterns. Similarly, in the four multicellular red algae (C. crispus, G. chorda, P. umbilicalis and G. domingensis), natural selection and other factors accounted for 74.2%, 84.7%, 86.9% and 45.9%, respectively. Moreover, as shown in Supplementary Table S7, GC12 and GC3 showed weak significant correlations in the six red algae (all r < 0.4, P < 0.01), strong significant correlations in G. domingensis (r = 0.730, P < 0.01) and general correlations in P. umbilicalis (r = 0.569, P < 0.01). In summary, mutation pressure may drive the codon usage bias of G. domingensis more strongly, while natural selection plays a leading role in the formation of codon usage patterns among the other seven red algae.

Fig. 4
figure 4

Neutrality plot analysis of eight red algae between GC12 (the mean GC content at the first and second positions) and GC3 (GC content at the third codon position). The black solid line represents the correlation line, and the equation of the correlation line is shown below the plot. The red dots represent the high-bias genes, the blue dots represent the medium-bias genes, and the green dots represent the low-bias genes. A: unicellular red algal group; B: multicellular red algal group

Parity rule 2 plot

The third base of the codon is closely related to the formation of codon usage preference. Parity rule analysis (PR2) of all genes across the 8 red algae is shown in Fig. 5, revealing that the points representing genes are unevenly distributed in four regions. Most of the points for the four unicellular red algae are scattered in the fourth quadrant (the ratio of A3/(A3 + T3) < 0.5 and G3/(G3 + C3) > 0.5) while the fewest are found in the first quadrant, which shows that at the third base of the codon, T is used more frequently than A, and G is used more frequently than C. In three multicellular red algae, most of the points are scattered in the third quadrant (the ratio of A3/(A3 + T3) < 0.5 and G3/(G3 + C3) < 0.5), illustrating that T and G are used more frequently for the third base. However, for G. domingensis, the number of points in the fourth quadrant was slightly greater than that in the third quadrant. In short, the third base of the codons shows a preference for G, C and T across all eight species.

Fig. 5
figure 5

Parity rule 2 plot analysis of eight red algae (GC bias on the x-axis and AT bias on that y-axis). The center point at 0.5 represents A = T and G = C, which means that there is no codon usage deviation between the two DNA strands. The red dots represent the high-bias genes, the blue dots represent the medium-bias genes, and the green dots represent the low-bias genes. A: unicellular red algal group; B: multicellular red algal group

Identification of optimal codons and major trna genes

Optimal codons can also be used to determine the influence strength of natural selection on codon bias. The results of the optimal codons for all 8 red algae are shown in Table 2 and Supplementary Table S3, from which it is seen the optimal codons of G. sulphuraria is end in A and U, except for one (UUG) ending in G. For the other seven species, 163 of 165 optimal codons end in G or C, except for CGU and GGU ending in uracil for C. yangmingshanensis. Furthermore, CGU and GGU are used with the lowest frequency in these optimal codons of arginine and glycine respectively (Supplementary Table S3). Significantly, there was no optimal codon ending in adenine for the seven species. Among the 18 amino acids including synonymous codons, phenylalanine is the only amino acid without an optimal codon in C. yangmingshanensis. Nine twofold degenerate amino acids have the same optimal codons, and seven of the nine threefold, fourfold and sixfold degenerate amino acids (except for valine and threonine) share at least one optimal codon, which indicates that the characteristics of optimal codons are significantly similar across all seven species.

Table 2 Optimal codons designated for the eight species of red algae

tRNA genes were identified to compare their optimal codons with those of major tRNA genes for all species. Major tRNA genes were defined as the most abundant tRNA gene for one amino acid. The tRNA genes of all eight species are shown as Supplementary Table S4, indicating that the copy number of tRNA genes was higher in multicellular species (97–165) than in unicellular red algae (22–26), except for G. sulphuraria (82). However, a few tRNA genes, such as tRNA Ile, appeared only in the unicellular species C. merolae and C. yangmingshanensis.

Except for methionine and tryptophan, 18 common amino acids encoded by 59 synonymous codons included nine twofold degenerate amino acids, one threefold degenerate amino acid, five fourfold degenerate amino acids and three sixfold degenerate amino acids. For the nine twofold degenerate amino acids, all major tRNA genes matched their optimal codons in G. sulphuraria (including 6 G-U matches) and most major tRNA genes corresponding to optimal codons (47/62) were also identified in the other 7 species (Table 2 and Supplementary Table S6). However, for threefold, fourfold and sixfold degenerate amino acids, there were only thirteen matches between optimal codons and the major tRNA genes in the seven species other than G. sulphuraria, in which the common match of glycine (GGC/GCC) simultaneously appeared in the six species other than G. sulphuraria and P. purpureum (Table 2). However, in most major tRNA genes in multicellular red algae (26/32), there were many bases of adenine at the anticodon position (Supplementary Tables S4, S5 and S6). It is inferred that these adenines are deaminated into inosine to pair with the cytosine of the optimal codon by the wobble effect.

Because adenine deamination appears rarely in unicellular red algae, the following analysis was performed on multicellular red algae. Among the six optimal codons perfectly matching major tRNA genes, CCG was complementary to the major tRNA genes of proline (GGC) in G. domingensis. AGC and CGC were complementary to the major tRNA genes of serine (GCU) and arginine (GCG) in C. crispus. GUG, ACG and UCG were complementary to the major tRNA genes of valine (CAC), threonine (CGU) and serine (CGA) in P. umbilicalis (Table 2, Supplementary Tables S4, S5 and S6). Except for these codons, for threefold to sixfold-degenerate amino acids, it was found that other optimal codons ended with cytosine. The matching major tRNA genes had adenine at the wobble site (Table 2 and Supplementary Tables S4, S5 and S6), which indicates that in the process of tRNA gene evolution, adenosines at the sites were deaminated to inosines to match the cytosine at the degenerate position. Overall, the major tRNA genes of most amino acids always matched their most frequently optimal codons in these species.

Surprisingly, using tRNAscan-SE software, the tRNASeC gene, the tRNA gene for selenocysteine (anticodon: TCA), was found to pair with the UGA codon in all four unicellular red algae but not in the four multicellular red algae (Supplementary Tables S4 and S6).

Correlation analysis

The expression of genes can be roughly assessed through the codon bias parameters, i.e., Fop (frequency of optimal codons), the CAI (codon adaptation index) and the CBI (codon bias index) [30,31,32,33]. Correlation analysis of all codon bias indicators of the protein-coding genes in all eight red algae is shown in Supplementary Table S7, showing that Nc was significantly correlated with Fop, the CAI and the CBI in all species. The correlation coefficient r was lower than 0.55(most r = 0.4–0.55, few r < 0.4) which indicates that the gene expression level may have some influence on the codon usage pattern of 8 species.

Protein domain codons show stronger codon preferences than nondomain codons

After filtering, the number of genes actually used this analysis of each species is shown in Supplementary Table S8. The codon bias of different regions of all genes in three categories is shown in Fig. 6 and Supplementary Table S8. For the high-bias genes category, Fop of the domain was markedly higher than that of the nondomain in 6 species (C. merolae, P. purpureum, C. crispus, G. chorda, G. domingensis and P. umbilicalis) and for the medium-bias gene category, Fop of the domain was markedly higher than that of the nondomain in seven species(except for G. sulphuraria), but for the low-bias gene category, only 2 species showed Fop values of in the domain that were markedly higher than those of the nondomain(C. crispus and P. umbilicalis). This phenomenon is consistent with the accuracy of translation selection, which is an important driving factor of codon usage in red algae, and it is more apparent in the high-bias category and the medium-bias category and even found in the low-bias category of individual species.

Fig. 6
figure 6

Average Fop value in domain codons and nondomain sequence codons. The red dots represent the average values of domain codons, and the blue-green dots represent the average values of nondomain codons. The line graph shows the standard deviation of each bias category value (*** P < 0.001; **P < 0.01; *P < 0.05). From left to right: high-bias genes, medium-bias genes and low-bias genes. A: unicellular red algal group; B: multicellular red algal group

Estimation of translation selection strength in eight red algae

Analyses of all eight genomes showed that natural selection plays an important role in translation accuracy. Researchers have designed many methods to determine the selection strength of codon usage. Here two different methods were used to estimate the strength of natural selection in Rhodophyta.

Assuming that the codon bias is in a mutation-drift-selection equilibrium, the odds ratio was used to estimate the translational selection strength in enterobacteria [34]. The odds ratios of nine twofold degenerate amino acids from eight red algae are shown in Supplementary Table S9, which reveals that except for Phe in C. yangmingshanensis, the odds ratios ranged from 1 to 11, which were similar to those in enterobacteria [34], holozoan protists [35], Drosophila melanogaster [36] and Arabidopsis thaliana [36].

Based on the population genetic model developed by Bulmer [24], widespread research on the codon usage of 80 bacterial genomes was performed to gain the most important statistical parameter, S (the log of the odds ratio) [22], which was used to estimate the translation selection strength. dos Reis, et al. [36] reported that Ŝ, the weighted average value of S in each genome can be used to accurately estimate the codon selection strength in eukaryotes. The S and Ŝ of eight red algae are shown in Supplementary Table S9 and Table 1, in which S ranged from 0 to 3 and Ŝ ranged from 0.6 to 1.5. Except for Phe, the S of other amino acids were greater than 1 in P. purpureum, C. crispus and P. umbilicalis. For the nine twofold degenerate amino acids, the S of 6 amino acids was more than 1 in G. sulphuraria, but S was less than 1 in the remaining species. Corresponding to S, Ŝ was greater than 1 in P. purpureum, C. crispus and P. umbilicalis, and the Ŝ for P. purpureum was similar to that for Monosiga brevicollis, and those for P. umbilicalis and C. crispus were slightly higher than those for Monosiga brevicollis [35]. The Ŝ for G. sulphuraria was calculated to be 0.97, which is similar to that for Drosophila melanogaster [36]. The Ŝ values of the other four species are comparable to those of some diatoms [14]. Previous studies have shown that there is selectivity in the codon usage of Monosiga brevicollis, Drosophila melanogaster and diatoms [14, 35, 36]. Therefore, there is also selection pressure on the codon usage of eight red algal genomes.

Discussion

Our analysis revealed similarities and differences in CUB in red algae. Most red algae prefer codons ending with CG, and especially in P. umbilicalis, the codons exhibit very strong GC bias (Table 1). The value is similar to the GC3 values of some species(0.70–0.86) in Chlorophyta reported by Li, et al. [37]. P. purpureum also has a strong bias for GC3, with values similar to the GC3 of Volvox carteri(0.697) in Chlorophyta and some species(0.6–0.63) in Monocotyledons [37]. C. yangmingshanensis, C. merolae and C. crispus showed a marginal preference for codons ending with GC, similar to the GC3 of Selaginella moellendorffii (0.593) in Pteridophyta and most species(0.57–0.63) in Monocotyledons [37]. G. chorda and G. domingensis show an almost neutral preference for four codons, with values similar to the GC3 of Physcomitrella patens (0.5) in Bryophyta, Musa acuminate in Monocotyledons (0.528), Linum usitatissimum(0.501) and Eucalyptus grandis(0.522) in Dicotyledons [37]. G. sulphuraria is different from other species and showed a strong preference for codons ending with AT, which is similar to the findings in Picea abies(0.438) [37], Pinus taeda(0.428) [37] as well as Taxus contorta (0.401) [38] of Gymnosperms, and most species (0.355–0.486) of Dicotyledons. Overall, in the eight red algae, as the ancestors of green plants, show a large degree of variation in the GC3 content of genes of the nuclear genomes during the evolution, in which P. umbilicalis is similar to Chlorophyta, G. chorda and G. domingensis are similar to Bryophyta, C. yangmingshanensis, C. merolae and C. crispus are similar to Pteridophyta, G. sulphuraria is similar to gymnosperms, P. purpureum is similar to Monocotyledons, and G. sulphuraria is similar to Eudicotyledons.

In addition, among the 7 species of red algae, the frequency of using codons ending with CG is significantly higher than that of AU, which is similar to some species in Chlorophyta, Pteridophyta and Monocotyledons [37]. The frequency of using codons ending with AU is significantly higher in G. sulphuraria, which is similar to observations in some species of Dicotyledons, Bryophyta and Gymnosperms [37]. Schonknecht, et al. [39] speculate that G. sulphuraria might have obtained a large number of AT-rich genes from extreme thermophilic bacteria and archaea through gene transfer.

Codon usage bias is gradually formed in various organisms mainly through selection, mutation and drift [20, 24]. Comparing plastid genomes [27, 28], the effects of mutations and drift on codon bias are limited in the nuclear genomes of red algae, and selection plays a dominant role in the formation of codon preference, especially in high-bias genes, which affects CUB by improving translation efficiency and translation accuracy. First, in our research, the Fop, CAI and CBI values of high-bias genes were higher than those of medium- and low-bias genes (Supplementary Table S10), consistent with findings in Escherichia coli [40], which indicates that red algae improve translation efficiency by using more optimal codons in high-bias genes. Second, in high- and medium-bias genes and even in individual low-bias genes of most species, the Fop values in the structural domain were greater than those in the nondomain (Fig. 6), in accordance with the results for Holozoan protists [35] and Escherichia coli [40], which indicates that red algae improve translation accuracy by using more optimal codons in the protein domain (conserved sites).

Interestingly, many tRNA modifications were found in the highly degenerate amino acids of all multicellular red algae and individual unicellular red algae in our research, which appear in Holozoan protists [35], most plants [41] and animals [42,43,44] by tRNA deamination under ADAT [41, 45,46,47].

Some studies have shown that transgenic design using optimal codons can improve the level of heterologous gene expression [48,49,50]. Some of the red algae studied here have important ecological, edible, medicinal, ecological, industrial and other value [51, 52], and red algae are also important materials for studying the origin and evolution of green plants. Our data provide important reference information for the successful expression of artificial transgenes. In transgenic engineering, using the optimal codon of a species can maximize the production of its encoded protein, increasing it by three orders of magnitude [53]. Table 2 shows the optimal codons and optimal termination codons of the eight species, which will make it possible to successfully apply transgenic design for each red alga in the future.

Conclusion

Overall, by studying the codon usage patterns of C. merolae, G. sulphuraria, C. yangmingshanensis, C. crispus, G. chorda, G. domingensis, P. purpureum and P. umbilicalis, it was found that high-bias genes in green plant ancestors show a preference for codons ending with GC. Correlation analysis indicated that the codon bias pattern was significantly correlated with Fop, the CAI and the CBI. Nc-GC3s plot analysis showed that the codon bias pattern of most genes is affected by selection because of deviation from the expected curve. Neutrality plot analysis indicates that natural selection plays a dominant role in the formation of codon usage patterns, which was also proven by comparing protein domain and nondomain codons. The analysis comparing the nucleotide content of intron and flanking sequences showed that natural selection strongly affects the codon usage pattern by driving the high GC content of high-bias genes, especially in unicellular species. Finally, we have obtained 15 common optimal codons in the seven red algae except for G. sulphuraria and the characteristics of related tRNA genes. Our results lay a foundation studying transgenic expression of economic red algae and the evolutionary study of codon usage characteristics of green plant ancestors.

Materials and methods

Data source

To analyze the codon usage pattern of green plant ancestors, the genome and annotated files of all red algae were downloaded from the NCBI Genome database, including four unicellular Rhodophyta (Cyanidioschyzon merolae, Cyanidiococcus yangmingshanensis, Galdieria sulphuraria and Porphyridium purpureum) and four multicellular Rhodophyta (Chondrus crispus, Gracilariopsis chorda, Gracilaria domingensis and Porphyra umbilicalis). The versions of the genome and its annotation file from the eight species are shown in Supplementary Table S11. By using TBtools software, we obtained the coding sequences of these species [54]. To reduce sampling error, we screened the sequences based on the following considerations: (1) the length of each CDS was a multiple of three; (2) ATG was the start codon, and UAA, UAG or UGA was the stop codon, (3) there was no termination codon in the middle of the CDS (Excluding the coding sequences with a premature termination); (4) the length of the sequence was greater than 300 bp; (5) repetitive sequences were removed, and the longest transcript was selected; (6) pseudogenes were eliminated; and (7) sequences with poor sequencing quality were excluded.

Codon usage analysis

Codon W 1.4.2 [55] was used to calculate the GC, GC1, GC2, GC3, Nc, CAI, CBI, FOP and RSCU values of each species and SPSS Statistics 27 was used to perform the correlation analysis among these parameters. According to the order of Nc values, the 5–10% of genes with the lowest values were chosen as the high-bias gene category (high-expression gene category), the 5–10% of genes with the highest values were chosen as the low-bias gene category (low-expression gene category), and the rest with the middle values were chosen as the medium-bias gene category. The detailed percentages for each species were as follows: 5% (C. merolae), 7% (C. yangmingshanensis), 5% (G. sulphuraria), 5% (P. purpureum), 8% (C. crispus), 9% (G. domingensis), 10% (G. chorda) and 5% (P. umbilicalis). The codons with an the RSCU difference between the high-expression gene category and low-expression gene category greater than 0.08 was defined as the optimal codons [56].

Defining the GC content of noncoding regions

The introns and flanking DNA of genes were extracted and screened from the genome and its annotation files according to the following rule. If conditions permit, we can extract 200 bp flanking DNA from each 5’ and 3’ end of each gene. When the intergenic region was < 200 bp, or a gene was located at the end of a contig or a scaffold, we extracted the maximum possible length of flanking DNA [35]. We extracted and connected all introns of each gene to calculate their GC content through Codon W, which is similar to the flanking DNA. Finally, the average value and standard deviation of each category were calculated.

Nc-GC3s plot analysis

Nc-GC3s plot analysis was used to explore the relationship between GC3s and Nc values [21]. In the rectangular coordinate system, the GC3s values of each gene were defined as the abscissa, and the Nc values were defined as the ordinate. After that, we drew an expected gene curve, which represented the expected position of genes only under neutral mutation pressure. Finally, we confirmed the effects of neutral mutations and natural selection on the codon bias of species through a comparison of expected gene values and actual values [21].

Neutrality plot analysis

Neutrality plot analysis was used to assess the effects of different drivers on codon bias. In the plot, the GC12(average of GC1 and GC2) of each gene was regarded as the x-axis and GC3 as the y-axis, establishing a rectangular coordinate system for regression analysis. The slope of the regression curve reflected the degree of influence of different drivers, mutation pressure or natural selection. When the slope was close to 0, most of the influence was natural selection [57].

Parity rule 2 plot

Using MEGA X software, we calculated the values of A3, T3, C3 and G3 for all genes of each species, regarding G3/(G3 + C3) and A3/(A3 + T3) as the abscissa and ordinate, respectively. The influence of mutation on codon usage bias can be inferred from the proportion of four bases in the figure. If preference was affected only by mutation, then the frequency of bases A/T and C/G at codon 3 would be similar, that is, A≈T or G≈C; otherwise it may be affected by other factors such as natural selection [58, 59].

Prediction and screening of tRNA genes

The downloaded genome sequence files of 8 red algae were scanned using the default parameters of tRNAscan-SE v.2.0.9 software to predict tRNA genes [60].

Frequency of optimal codon usage for domain and nondomain codons

To determine whether codon bias affects different regions of genes in red algae, the frequency of optimal codons (Fop) was measured and compared between the domain and nondomain regions of the three bias categories. First, each gene was divided into putative functional domain and nondomain codons according to the NCBI's Preserved Domain Database (CDD). However, according to the method of Southworth [35], the following genes were excluded: (1) genes that did not code annotation regions were excluded from the analysis; (2) when the functional domain spanned all codons of the whole gene, the gene was excluded; and (3) when the functional domain spanned all codons except for the start codon and/or the stop codon, the gene was excluded. Second, the values of Fop for the domain and nondomain regions were determined by Codon W. Because some domain and nondomain regions may be too short to calculate the value of Nc, Fop was used to determine the codon use bias of the regions [35].

Estimation of the strength of translational selection

The strength of translational selection was estimated by using the following two methods according to Eyre-Walker, et al. [34] and dos Reis, et al. [36]. First, the highly expressed gene category and lowly expressed gene category were determined. Then, the optimal codon and the suboptimal codon were identified for both categories. Finally, the parameters of translation selection strength, odds ratios and S were calculated. The odds ratio can be calculated by calculating the relative difference of 2 degenerate synonymous codons between the highly expressed gene category and the lowly expressed gene category (Formula 1). S was the logarithm of the odds ratio used to estimate strength of translational selection [22]. Using this method, the odds ratio and the S of eight red algae were calculated [34].

$$\text{odds ratio }=\frac{\text{f}1H}{\text{f}2H} *\frac{\text{f}2L}{\text{f}1H}$$
(1)

In the above formula, f represents the codon frequency. The subscripts 1 and 2 represent the optimal codon and suboptimal codon respectively; the subscripts H and L represent the high-bias gene category and the low-bias gene category, respectively.

By the second method from dos Reis, the Ŝ of the eight genomes was determined by calculating the weighted average of S for the number of codons of the nine 2-degenerate amino acids in the high-bias gene category to estimate the codon usage bias strength of the species genomes[22, 34, 36].

Availability of data and materials

The taxonomic information and genomes data that support the findings of this study are openly available in [NCBI] at https://www.ncbi.nlm.nih.gov/. The remaining data generated or analyzed during this study are available within the article and its supplementary materials.

References

  1. Adl SM, Simpson AG, Farmer MA, Andersen RA, Anderson OR, Barta JR, Bowser SS, Brugerolle G, Fensome RA, Fredericq S, et al. The new higher level classification of eukaryotes with emphasis on the taxonomy of protists. J Eukaryot Microbiol. 2005;52(5):399–451.

    PubMed  Google Scholar 

  2. Chan CX, Gross J, Yoon HS, Bhattacharya D. Plastid origin and evolution: new models provide insights into old problems. Plant Physiol. 2011;155(4):1552–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Chan CX, Yang EC, Banerjee T, Yoon HS, Martone PT, Estevez JM, Bhattacharya D. Red and Green Algal Monophyly and Extensive Gene Sharing Found in a Rich Repertoire of Red Algal Genes. Curr Biol. 2011;21(4):328–33.

    CAS  PubMed  Google Scholar 

  4. Price DC, Goodenough UW, Roth R, Lee JH, Kariyawasam T, Mutwil M, Ferrari C, Facchinelli F, Ball SG, Cenci U, et al. Analysis of an improved Cyanophora paradoxa genome assembly. DNA Res. 2019;26(4):287–99.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Rodriguez-Ezpeleta N, Brinkmann H, Burey SC, Roure B, Burger G, Loffelhardt W, Bohnert HJ, Philippe H, Lang BF. Monophyly of primary photosynthetic eukaryotes: green plants, red algae, and glaucophytes. Curr Biol. 2005;15(14):1325–30.

    CAS  PubMed  Google Scholar 

  6. Reyes-Prieto A, Weber AP, Bhattacharya D. The origin and establishment of the plastid in algae and plants. Annu Rev Genet. 2007;41:147–68.

    CAS  PubMed  Google Scholar 

  7. Bengtson S, Sallstedt T, Belivanova V, Whitehouse M. Three-dimensional preservation of cellular and subcellular structures suggests 1.6 billion-year-old crown-group red algae. Plos Biol. 2017;15(3):e2000735.

    PubMed  PubMed Central  Google Scholar 

  8. Butterfield NJ. Bangiomorpha pubescens n. gen., n. sp.: implications for the evolution of sex, multicellularity, and the Mesoproterozoic/Neoproterozoic radiation of eukaryotes. Paleobiology. 2000;26(3):386–404.

    Google Scholar 

  9. Xiao S, Knoll AH, Yuan X, Pueschel CM. Phosphatized multicellular algae in the Neoproterozoic Doushantuo Formation, China, and the early evolution of florideophyte red algae. Am J Bot. 2004;91(2):214–27.

    PubMed  Google Scholar 

  10. Xiao S, Yuan X, Steiner M, Knoll AH. Macroscopic Carbonaceous Compressions in a Terminal Proterozoic Shale: A Systematic Reassessment of the Miaohe Biota, South China. J Paleontol. 2002;76(2):347–76.

    Google Scholar 

  11. Becker B. Snow ball earth and the split of Streptophyta and Chlorophyta. Trends Plant Sci. 2013;18(4):180–3.

    CAS  PubMed  Google Scholar 

  12. Ciniglia C, Yoon HS, Pollio A, Pinto G, Bhattacharya D. Hidden biodiversity of the extremophilic Cyanidiales red algae. Mol Ecol. 2004;13(7):1827–38.

    CAS  PubMed  Google Scholar 

  13. Douzery EJ, Snell EA, Bapteste E, Delsuc F, Philippe H. The timing of eukaryotic evolution: does a relaxed molecular clock reconcile proteins and fossils? Proc Natl Acad Sci U S A. 2004;101(43):15386–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Parfrey LW, Lahr DJ, Knoll AH, Katz LA. Estimating the timing of early eukaryotic diversification with multigene molecular clocks. Proc Natl Acad Sci U S A. 2011;108(33):13624–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Schönknecht G, Weber AP, Lercher MJ. Horizontal gene acquisitions by eukaryotes as drivers of adaptive evolution. BioEssays. 2014;36(1):9–20.

    PubMed  Google Scholar 

  16. Initiative OTPT. One thousand plant transcriptomes and the phylogenomics of green plants. Nature. 2019;574(7780):679–85.

    Google Scholar 

  17. Cavalier-Smith T, Chao EE. Phylogeny of choanozoa, apusozoa, and other protozoa and early eukaryote megaevolution. J Mol Evol. 2003;56(5):540–63.

    CAS  PubMed  Google Scholar 

  18. Moustafa A, Beszteri B, Maier UG, Bowler C, Valentin K, Bhattacharya D. Genomic footprints of a cryptic plastid endosymbiosis in diatoms. Science. 2009;324(5935):1724–6.

    CAS  PubMed  Google Scholar 

  19. Sanchez-Baracaldo P, Raven JA, Pisani D, Knoll AH. Early photosynthetic eukaryotes inhabited low-salinity habitats. Proc Natl Acad Sci U S A. 2017;114(37):E7737–45.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Parvathy ST, Udayasuriyan V, Bhadana V. Codon usage bias. Mol Biol Rep. 2022;49(1):539–65.

    CAS  PubMed  Google Scholar 

  21. Wright F. The “effective number of codons” used in a gene. Gene. 1990;87(1):23–9.

    CAS  PubMed  Google Scholar 

  22. Sharp PM, Bailes E, Grocock RJ, Peden JF, Sockett RE. Variation in the strength of selected codon usage bias among bacteria. Nucleic Acids Res. 2005;33(4):1141–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Sharp PM, Li WH. An evolutionary perspective on synonymous codon usage in unicellular organisms. J Mol Evol. 1986;24(1–2):28–38.

    CAS  PubMed  Google Scholar 

  24. Bulmer M. The selection-mutation-drift theory of synonymous codon usage. Genetics. 1991;129(3):897–907.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Ikemura T. Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes: a proposal for a synonymous codon choice that is optimal for the E coli translational system. J Mol Biol. 1981;151(3):389–409.

    CAS  PubMed  Google Scholar 

  26. Lloyd AT, Sharp PM. Synonymous codon usage in Kluyveromyces lactis. Yeast. 1993;9(11):1219–28.

    CAS  PubMed  Google Scholar 

  27. Suzuki H, Morton BR. Codon Adaptation of Plastid Genes. PLoS ONE. 2016;11(5):e0154306.

    PubMed  PubMed Central  Google Scholar 

  28. Li G, Pan Z, Gao S, He Y, Xia Q, Jin Y, Yao H. Analysis of synonymous codon usage of chloroplast genome in Porphyra umbilicalis. Genes Genomics. 2019;41(10):1173–81.

    CAS  PubMed  Google Scholar 

  29. Lee H, Lee HK, An G, Lee YK. Analysis of expressed sequence tags from the red alga Griffithsia okiensis. J Microbiol. 2007;45(6):541–6.

    CAS  PubMed  Google Scholar 

  30. Fox JM, Erill I. Relative codon adaptation: a generic codon bias index for prediction of gene expression. DNA Res. 2010;17(3):185–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Carbone A, Zinovyev A, Képès F. Codon adaptation index as a measure of dominating codon bias. Bioinformatics. 2003;19(16):2005–15.

    CAS  PubMed  Google Scholar 

  32. Ikemura T. Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes. J Mol Biol. 1981;146(1):1–21.

    CAS  PubMed  Google Scholar 

  33. Bennetzen JL, Hall BD. Codon selection in yeast. J Biol Chem. 1982;257(6):3026–31.

    CAS  PubMed  Google Scholar 

  34. Eyre-Walker A, Bulmer M. Synonymous substitution rates in enterobacteria. Genetics. 1995;140(4):1407–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Southworth J, Armitage P, Fallon B, Dawson H, Bryk J, Carr M. Patterns of Ancestral Animal Codon Usage Bias Revealed through Holozoan Protists. Mol Biol Evol. 2018;35(10):2499–511.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. dos Reis M, Wernisch L. Estimating translational selection in eukaryotic genomes. Mol Biol Evol. 2009;26(2):451–61.

    PubMed  Google Scholar 

  37. Li N, Li YY, Zheng CC, Huang JG, Zhang SZ. Genome-wide comparative analysis of the codon usage patterns in plants. Genes Genom. 2016;38(8):723–31.

    CAS  Google Scholar 

  38. Majeed A, Kaur H, Bhardwaj P. Selection constraints determine preference for A/U-ending codons in Taxus contorta. Genome. 2020;63(4):215–24.

    CAS  PubMed  Google Scholar 

  39. Schonknecht G, Chen WH, Ternes CM, Barbier GG, Shrestha RP, Stanke M, Brautigam A, Baker BJ, Banfield JF, Garavito RM, et al. Gene transfer from bacteria and archaea facilitated evolution of an extremophilic eukaryote. Science. 2013;339(6124):1207–10.

    PubMed  Google Scholar 

  40. Yannai A, Katz S, Hershberg R. The Codon Usage of Lowly Expressed Genes Is Subject to Natural Selection. Genome Biol Evol. 2018;10(5):1237–46.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Ramos-Morales E, Bayam E, Del-Pozo-Rodriguez J, Salinas-Giege T, Marek M, Tilly P, Wolff P, Troesch E, Ennifar E, Drouard L, et al. The structure of the mouse ADAT2/ADAT3 complex reveals the molecular basis for mammalian tRNA wobble adenosine-to-inosine deamination. Nucleic Acids Res. 2021;49(11):6529–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Rafels-Ybern A, Torres AG, Grau-Bove X, Ruiz-Trillo I, de Pouplana LR. Codon adaptation to tRNAs with Inosine modification at position 34 is widespread among Eukaryotes and present in two Bacterial phyla. Rna Biol. 2018;15(4–5):500–7.

    PubMed  Google Scholar 

  43. Yu Y, Zhou HX, Kong YM, Pan BH, Chen LX, Wang HB, Hao P, Li X: The Landscape of A-to-I RNA Editome Is Shaped by Both Positive and Purifying Selection. Plos Genetics. 2016;12(7):e1006191.

  44. Karcher D, Bock R. Identification of the chloroplast adenosine-to-inosine tRNA editing enzyme. RNA. 2009;15(7):1251–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Frigole HR, Camacho N, Coma MC, Fernandez-Lozano C, Garcia-Lema J, Rafels-Ybern A, Canals A, Coll M, de Pouplana LR. tRNA deamination by ADAT requires substrate-specific recognition mechanisms and can be inhibited by tRFs. RNA. 2019;25(5):607–19.

    CAS  Google Scholar 

  46. Crick FH. Codon–anticodon pairing: the wobble hypothesis. J Mol Biol. 1966;19(2):548–55.

    CAS  PubMed  Google Scholar 

  47. Tarrant D, von der Haar T. Synonymous codons, ribosome speed, and eucaryotic gene expression regulation. Cell Mol Life Sci. 2014;71(21):4195–206.

  48. Frumkin I, Lajoie MJ, Gregg CJ, Hornung G, Church GM, Pilpel Y. Codon usage of highly expressed genes affects proteome-wide translation efficiency. Proc Natl Acad Sci U S A. 2018;115(21):E4940–9.

    PubMed  PubMed Central  Google Scholar 

  49. Ma L, Cui P, Zhu J, Zhang Z, Zhang Z. Translational selection in human: more pronounced in housekeeping genes. Biol Direct. 2014;9:17.

    PubMed  PubMed Central  Google Scholar 

  50. Ingvarsson PK. Molecular evolution of synonymous codon usage in Populus. BMC Evol Biol. 2008;8:307.

    PubMed  PubMed Central  Google Scholar 

  51. Bixler HJ, Porse H. A decade of change in the seaweed hydrocolloids industry. J Appl Phycol. 2011;23(3):321–35.

    Google Scholar 

  52. Blouin NA, Brodie JA, Grossman AC, Xu P, Brawley SH. Porphyra: a marine crop shaped by stress. Trends Plant Sci. 2011;16(1):29–37.

    CAS  PubMed  Google Scholar 

  53. Gustafsson C, Minshull J, Govindarajan S, Ness J, Villalobos A, Welch M. Engineering genes for predictable protein expression. Protein Expres Purif. 2012;83(1):37–46.

    CAS  Google Scholar 

  54. Chen CJ, Chen H, Zhang Y, Thomas HR, Frank MH, He YH, Xia R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol Plant. 2020;13(8):1194–202.

    CAS  PubMed  Google Scholar 

  55. Peden JF. Analysis of codon usage [PhD dissertation]. University of Nottingham; 1999.

  56. Duret L, Mouchiroud D. Expression pattern and surprisingly, gene length shape codon usage in caenorhabditis, drosophila, and arabidopsis. Proc Natl Acad Sci U S A. 1999;96(8):4482–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Sueoka N. Directional mutation pressure and neutral molecular evolution. Proc Natl Acad Sci U S A. 1988;85(8):2653–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  58. Sueoka N. Near homogeneity of PR2-bias fingerprints in the human genome and their implications in phylogenetic analyses. J Mol Evol. 2001;53(4–5):469–76.

    CAS  PubMed  Google Scholar 

  59. Sueoka N. 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. 1999;238(1):53–8.

    CAS  PubMed  Google Scholar 

  60. Chan PP, Lin BY, Mak AJ, Lowe TM. tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes. Nucleic Acids Res. 2021;49(16):9077–96.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to all those who have supported this work.

Funding

This work was jointly funded by the 2023 Provincial Science and Technology Plan (2023YFH0103) and the 2022 Double Support Project of Sichuan Agricultural University.

Author information

Authors and Affiliations

Authors

Contributions

HPY and TTL identified the research topic and conceived the experiments. YXZ collected and download the data. XYW and LXX participated in the data processing of experiments. TTL, ZM and HPY performed and finished the data analysis for experiments. TTL wrote the manuscript. YC and ZZT participated in partial revision of the article. HPY was responsible for the entire article revision, and finalization. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Huipeng Yao.

Ethics declarations

Ethics approval and consent to participate

There is no any research involving humans or animals in this article.

Consent for publication

Not applicable.

Competing interests

The authors declare that there is no conflict of interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplementary Table S1. GC3s of CDSs, GC content of flanking sequences and introns respectively in the three bias categories from the eight red algae.

Additional file 2:

Supplementary Table S2. Correlation analysis and t test (double tail) for the GC3s of CDSs and GC content of introns among the eight red algae.

Additional file 3:

Supplementary Table S3. Output of optimal codon analysis from Codon W in the eight red algae (S3a-S3h).

Additional file 4:

Supplementary Table S4. Output of tRNAscan-SE from whole genome contigs in the eight red algae (S4a-S4h).

Additional file 5:

Supplementary Table S5. The optimal codon which lack perfectly matching complementary tRNA genes in multicellular red algae.

Additional file 6:

Supplementary Table S6. Number of tRNAs and anticodons of different amino acids in the eight red algae (S6a-S6h).

Additional file 7:

Supplementary Table S7. Correlation analysis and significance test (t test, double tail) for 12 parameters related to codon usage among eight red algae (S7a-S7h).

Additional file 8:

 Supplementary Table S8. The average Fop values for the domain and nondomain codons of the three bias gene categories in the eight red algae.

Additional file 9:

 Supplementary Table S9. Odds ratios and their logarithmic values(Ln) of two degenerate amino acids among the eight red algae.

Additional file 10:

Supplementary Table S10. Mean and standard deviation of CAI, CBI, and Fop for each bias category of the seven red algae (excepting for Galdieria sulphuraria).

Additional file 11:

Supplementary Table S11. Versions of genomes and annotation files and taxonomic information for the eight red algae.

Additional file 12:

Supplementary Figure S1. The plotting for GC3s of CDSs and GC content of introns among eight species. The correlation between two statistical data points is represented by the plotting of the GC3s of CDSs and GC content of introns in the same gene. The black solid line represents the correlation line, and its equation is shown at the top of the plot.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yao, H., Li, T., Ma, Z. et al. Codon usage pattern of the ancestor of green plants revealed through Rhodophyta. BMC Genomics 24, 538 (2023). https://doi.org/10.1186/s12864-023-09586-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09586-w

Keywords