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

Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20-30 bp and 10-20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Repetitive Sequences, Codon Usage Bias and Phylogenetic Analysis of the Plastome of Miliusa glochidioides

Annonaceae is the largest family in Magnoliales, exhibiting the greatest diversity among and within genera. In this study, we conducted an analysis of repetitive sequences and codon usage bias in the previously acquired plastome of Miliusa glochidioides. Using a concatenated dataset of shared genes, we constructed the phylogenetic relationships among 27 Annonaceae species. The results showed that the size of the plastomes in the Annonaceae ranged from 159 to 202 kb, with the size of the inverted repeat region ranging from 40 to 65 kb. Within the plastome of M. glochidioides, we identified 42 SSRs, 36 tandem repeats, and 9 dispersed repeats. These SSRs consist of three nucleotide types and eight motif types, with a preference for A/T bases, primarily located in the large single-copy regions and intergenic spacers. Tandem and dispersed repeat sequences were predominantly detected in the IR region. Through codon usage bias analysis, we identified 30 high-frequency codons and 11 optimal codons. The plastome of M. glochidioides demonstrated relatively weak codon usage bias, favoring codons with A/T endings, primarily influenced by natural selection. Phylogenetic analysis revealed that all four subfamilies formed monophyletic groups, with Cananga odorata (Ambavioideae) and Anaxagorea javanica (Anaxagoreoideae) successively nested outside Annonoideae + Malmeoideae. These findings improve our understanding of the plastome of M. glochidioides and provide additional insights for studying plastome evolution in Annonaceae.

Analysis of chloroplast genome characteristics and codon usage bias in 14 species of Annonaceae

For the study of species evolution, chloroplast gene expression, and transformation, the chloroplast genome is an invaluable resource. Codon usage bias (CUB) analysis is a tool that is utilized to improve gene expression and investigate evolutionary connections in genetic transformation. In this study, we analysed chloroplast genome differences, codon usage patterns and the sources of variation on CUB in 14 Annonaceae species using bioinformatics tools. The study showed that there was a significant variation in both gene sizes and numbers between the 14 species, but conservation was still maintained. It's worth noting that there were noticeable differences in the IR/SC sector boundary and the types of SSRs among the 14 species. The mono-nucleotide repeat type was the most common, with A/T repeats being more prevalent than G/C repeats. Among the different types of repeats, forward and palindromic repeats were the most abundant, followed by reverse repeats, and complement repeats were relatively rare. Codon composition analysis revealed that all 14 species had a frequency of GC lower than 50%. Additionally, it was observed that the proteins in-coding sequences of chloroplast genes tend to end with A/T at the third codon position. Among these species, 21 codons exhibited bias (RSCU > 1), and there were 8 high-frequency (HF) codons and 5 optimal codons that were identical across the species. According to the ENC-plot and Neutrality plot analysis, natural selection had less impact on the CUB of A. muricate and A. reticulata. Based on the PR2-plot, it was evident that base G had a higher frequency than C, and T had a higher frequency A. The correspondence analysis (COA) revealed that codon usage patterns different in Annonaceae.

Comparative analysis of codon usage bias in chloroplast genomes of ten medicinal species of Rutaceae

Rutaceae family comprises economically important plants due to their extensive applications in spices, food, oil, medicine, etc. The Rutaceae plants is able to better utilization through biotechnology. Modern biotechnological approaches primarily rely on the heterologous expression of functional proteins in different vectors. However, several proteins are difficult to express outside their native environment. The expression potential of functional genes in heterologous systems can be maximized by replacing the rare synonymous codons in the vector with preferred optimal codons of functional genes. Codon usage bias plays a critical role in biogenetic engineering-based research and development. In the current study, 727 coding sequences (CDSs) obtained from the chloroplast genomes of ten Rutaceae plant family members were analyzed for codon usage bias. The nucleotide composition analysis of codons showed that these codons were rich in A/T(U) bases and preferred A/T(U) endings. Analyses of neutrality plots, effective number of codons (ENC) plots, and correlations between ENC and codon adaptation index (CAI) were conducted, which revealed that natural selection is a major driving force for the Rutaceae plant family's codon usage bias, followed by base mutation. In the ENC vs. CAI plot, codon usage bias in the Rutaceae family had a negligible relationship with gene expression level. For each sample, we screened 12 codons as preferred and high-frequency codons simultaneously, of which GCU encoding Ala, UUA encoding Leu, and AGA encoding Arg were the most preferred codons. Taken together, our study unraveled the synonymous codon usage pattern in the Rutaceae family, providing valuable information for the genetic engineering of Rutaceae plant species in the future.

Investigation of CaMV-host co-evolution through synonymous codon pattern

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Decoding Evolution of Rubioideae: Plastomes Reveal Sweet Secrets of Codon Usage, Diagnostides, and Superbarcoding

Galium genus belongs to the Rubiaceae family, which consists of approximately 14,000 species. In comparison to its well-known relatives, the plastomes of the Galium genus have not been explored so far. The plastomes of this genus have a typical, quadripartite structure, but differ in gene content, since the infA gene is missing in Galium palustre and Galium trfidum. An evaluation of the effectiveness of using entire chloroplast genome sequences as superbarcodes for accurate plant species identification revealed the high potential of this method for molecular delimitation within the genus and tribe. The trnE-UUC-psbD region showed the biggest number of diagnostides (diagnostic nucleotides) which might be new potential barcodes, not only in Galium, but also in other closely related genera. Relative synonymous codon usage (RSCU) appeared to be connected with the phylogeny of the Rubiaceae family, showing that during evolution, plants started preferring specific codons over others.

Twelve newly assembled jasmine chloroplast genomes: unveiling genomic diversity, phylogenetic relationships and evolutionary patterns among Oleaceae and Jasminum species

BACKGROUND

Jasmine (Jasminum), renowned for its ornamental value and captivating fragrance, has given rise to numerous species and accessions. However, limited knowledge exists regarding the evolutionary relationships among various Jasminum species.

RESULTS

In the present study, we sequenced seven distinct Jasminum species, resulting in the assembly of twelve high-quality complete chloroplast (cp) genomes. Our findings revealed that the size of the 12 cp genomes ranged from 159 to 165 kb and encoded 134-135 genes, including 86-88 protein-coding genes, 38-40 tRNA genes, and 8 rRNA genes. J. nudiflorum exhibited a larger genome size compared to other species, mainly attributed to the elevated number of forward repeats (FRs). Despite the typically conservative nature of chloroplasts, variations in the presence or absence of accD have been observed within J. sambac. The calculation of nucleotide diversity (Pi) values for 19 cp genomes indicated that potential mutation hotspots were more likely to be located in LSC regions than in other regions, particularly in genes ycf2, rbcL, atpE, ndhK, and ndhC (Pi > 0.2). Ka/Ks values revealed strong selection pressure on the genes rps2, atpA, rpoA, rpoC1, and rpl33 when comparing J. sambac with the three most closely related species (J. auriculatum, J. multiflorum, and J. dichotomum). Additionally, SNP identification, along with the results of Structure, PCA, and phylogenetic tree analyses, divided the Jasminum cp genomes into six groups. Notably, J. polyanthum showed gene flow signals from both the G5 group (J. nudiflorum) and the G3 group (J. tortuosum and J. fluminense). Phylogenetic tree analysis reflected that most species from the same genus clustered together with robust support in Oleaceae, strongly supporting the monophyletic nature of cp genomes within the genus Jasminum.

CONCLUSION

Overall, this study provides comprehensive insights into the genomic composition, variation, and phylogenetic relationships among various Jasminum species. These findings enhance our understanding of the genetic diversity and evolutionary history of Jasminum.

Intraspecific and Intrageneric Genomic Variation across Three Sedum Species (Crassulaceae): A Plastomic Perspective

Sedum is the largest succulent genus in Crassulaceae. Because of predominant maternal inheritance, little recombination, and slow evolution, plastomes can serve as powerful super barcodes for inter- or intra-species phylogenetic analyses. While previous research has focused on plastomes between Sedum species, intra-species studies are scarce. Here, we sequenced plastomes from three Sedum species (Sedum alfredii, Sedum plumbizincicola, and Sedum japonicum) to understand their evolutionary relationships and plastome structural evolution. Our analyses revealed minimal size and GC content variation across species. However, gene distribution at IR boundaries, repeat structures, and codon usage patterns showed diversity at both inter-specific and intra-specific levels. Notably, an rps19 gene expansion and a bias toward A/T-ending codons were observed. Codon aversion motifs also varied, potentially serving as markers for future studies. Phylogenetic analyses confirmed the non-monophyly of Sedum and divided the Acre clade into two groups. Individuals from the same species clustered together, with strong support for the relationships between S. alfredii, S. tricarpum, and S. plumbizincicola. Additionally, S. japonicum clearly affiliates with the Acre clade. This study provides valuable insights into both intra-specific and intra-generic plastome variation in Sedum, as well as overall plastome evolution within the genus.

Plastid genome data provide new insights into the dynamic evolution of the tribe Ampelopsideae (Vitaceae)

BACKGROUND

Ampelopsideae J. Wen & Z.L. Nie is a small-sized tribe of Vitaceae Juss., including ca. 47 species from four genera showing a disjunct distribution worldwide across all the continents except Antarctica. There are numerous species from the tribe that are commonly used as medicinal plants with immune-modulating, antimicrobial, and anti-hypertensive properties. The tribe is usually recognized into three clades, i.e., Ampelopsis Michx., Nekemias Raf., and the Southern Hemisphere clade. However, the relationships of the three clades differ greatly between the nuclear and the plastid topologies. There has been limited exploration of the chloroplast phylogenetic relationships within Ampelopsideae, and studies on the chloroplast genome structure of this tribe are only available for a few individuals. In this study, we aimed to investigate the evolutionary characteristics of plastid genomes of the tribe, including their genome structure and evolutionary insights.

RESULTS

We sequenced, assembled, and annotated plastid genomes of 36 species from the tribe and related taxa in the family. Three main clades were recognized within Ampelopsideae, corresponding to Ampelopsis, Nekemias, and the Southern Hemisphere lineage, respectively, and all with 100% bootstrap supports. The genome sequences and content of the tribe are highly conserved. However, comparative analyses suggested that the plastomes of Nekemias demonstrate a contraction in the large single copy region and an expansion in the inverted repeat region, and possess a high number of forward and palindromic repeat sequences distinct from both Ampelopsis and the Southern Hemisphere taxa.

CONCLUSIONS

Our results highlighted plastome variations in genome length, expansion or contraction of the inverted repeat region, codon usage bias, and repeat sequences, are corresponding to the three lineages of the tribe, which probably faced with different environmental selection pressures and evolutionary history. This study provides valuable insights into understanding the evolutionary patterns of plastid genomes within the Ampelopsideae of Vitaceae.

Codon usage bias and phylogenetic analysis of chloroplast genome in 36 gracilariaceae species

Gracilariaceae is a group of marine large red algae and main source of agar with important economic and ecological value. The codon usage patterns of chloroplast genomes in 36 species from Graciliaceae show that GC range from 0.284 to 0.335, the average GC3 range from 0.135 to 0.243 and the value of ENC range from 35.098 to 42.327, which indicates these genomes are rich in AT and prefer to use codons ending with AT in these species. Nc plot, PR2 plot, neutrality plot analyses and correlation analysis indicate that these biases may be caused by multiple factors, such as natural selection and mutation pressure, but prolonged natural selection is the main driving force influencing codon usage preference. The cluster analysis and phylogenetic analysis show that the differentiation relationship of them is different and indicate that codons with weak or unbiased preferences may also play an irreplaceable role in these species' evolution. In addition, we identified 26 common high-frequency codons and 8-18 optimal codons all ending in A/U in these 36 species. Our results will not only contribute to carrying out transgenic work in Gracilariaceae species to maximize the protein yield in the future, but also lay a theoretical foundation for further exploring systematic classification of them.

Genome-wide identification and characterization of TCP gene family in Dendrobium nobile and their role in perianth development

TCP is a widely distributed, essential plant transcription factor that regulates plant growth and development. An in-depth study of TCP genes in Dendrobium nobile, a crucial parent in genetic breeding and an excellent model material to explore perianth development in Dendrobium, has not been conducted. We identified 23 DnTCP genes unevenly distributed across 19 chromosomes and classified them as Class I PCF (12 members), Class II: CIN (10 members), and CYC/TB1 (1 member) based on the conserved domain and phylogenetic analysis. Most DnTCPs in the same subclade had similar gene and motif structures. Segmental duplication was the predominant duplication event for TCP genes, and no tandem duplication was observed. Seven genes in the CIN subclade had potential miR319 and -159 target sites. Cis-acting element analysis showed that most DnTCP genes contained many developmental stress-, light-, and phytohormone-responsive elements in their promoter regions. Distinct expression patterns were observed among the 23 DnTCP genes, suggesting that these genes have diverse regulatory roles at different stages of perianth development or in different organs. For instance, DnTCP4 and DnTCP18 play a role in early perianth development, and DnTCP5 and DnTCP10 are significantly expressed during late perianth development. DnTCP17, 20, 21, and 22 are the most likely to be involved in perianth and leaf development. DnTCP11 was significantly expressed in the gynandrium. Specially, MADS-specific binding sites were present in most DnTCP genes putative promoters, and two Class I DnTCPs were in the nucleus and interacted with each other or with the MADS-box. The interactions between TCP and the MADS-box have been described for the first time in orchids, which broadens our understanding of the regulatory network of TCP involved in perianth development in orchids.

Evolutionary analysis and biological characterization of a novel alphabaculovirus isolated from Mythimna separata

Baculoviruses are insect-specific pathogens. Novel baculovirus isolates provide new options for the biological control of pests. Therefore, research into the biological characteristics of newly isolated baculoviruses, including accurate classification and nomenclature, is important. In this study, a baculovirus was isolated from Mythimna separata and its complete genome sequence was determined by next-generation sequencing. The double-stranded DNA genome was 153 882 bp in length, encoding 163 open reading frames. The virus was identified as a variant of Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV) and designated Mamestra brassicae multiple nucleopolyhedrovirus CHN1 (MbMNPV-CHN1) according to ultrastructural analysis, genome comparison and phylogenetic analysis. Phylogenetic inference placed MbMNPV-CHN1 in a clade containing isolates of MacoNPV-A, MacoNPV-B and MbMNPV, which we have designated the Mb-McNPV group. The genomes of isolates in the Mb-McNPV group exhibited a high degree of collinearity with relatively minor differences in the content of annotated open reading frames. The development of codon usage bias in the Mb-McNPV group was affected mainly by natural selection. MbMNPV-CHN1 shows high infectivity against seven species of Lepidoptera. The yield of MbMNPV-CHN1 in the fourth- and fifth-instar M. separata larvae was 6.25×109-1.23×1010 OBs/cadaver. Our data provide insights into the classification, host range and virulence differences among baculoviruses of the Mb-McNPV group, as well as a promising potential new baculoviral insecticide.

Comparative analysis of codon usage patterns in the chloroplast genomes of nine forage legumes

Leguminosae is one of the three largest families of angiosperms after Compositae and Orchidaceae. It is widely distributed and grows in a variety of environments, including plains, mountains, deserts, forests, grasslands, and even waters where almost all legumes can be found. It is one of the most important sources of starch, protein and oil in the food of mankind and also an important source of high-quality forage material for animals, which has important economic significance. In our study, the codon usage patterns and variation sources of the chloroplast genome of nine important forage legumes were systematically analyzed. Meanwhile, we also constructed a phylogenetic tree based on the whole chloroplast genomes and protein coding sequences of these nine forage legumes. Our results showed that the chloroplast genomes of nine forage legumes end with A/T bases, and seven identical high-frequency (HF) codons were detected among the nine forage legumes. ENC-GC3s mapping, PR2 analysis, and neutral analysis showed that the codon bias of nine forage legumes was influenced by many factors, among which natural selection was the main influencing factor. The codon usage frequency showed that the Nicotiana tabacum and Saccharomyces cerevisiae can be considered as receptors for the exogenous expression of chloroplast genes of these nine forage legumes. The phylogenetic relationships of the chloroplast genomes and protein coding genes were highly similar, and the nine forage legumes were divided into three major clades. Among the clades Melilotus officinalis was more closely related to Medicago sativa, and Galega officinalis was more closely related to Galega orientalis. This study provides a scientific basis for the molecular markers research, species identification and phylogenetic studies of forage legumes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01421-0.

Codon usage characterization and phylogenetic analysis of the mitochondrial genome in Hemerocallis citrina

BACKGROUND

Hemerocallis citrina Baroni is a traditional vegetable crop widely cultivated in eastern Asia for its high edible, medicinal, and ornamental value. The phenomenon of codon usage bias (CUB) is prevalent in various genomes and provides excellent clues for gaining insight into organism evolution and phylogeny. Comprehensive analysis of the CUB of mitochondrial (mt) genes can provide rich genetic information for improving the expression efficiency of exogenous genes and optimizing molecular-assisted breeding programmes in H. citrina.

RESULTS

Here, the CUB patterns in the mt genome of H. citrina were systematically analyzed, and the possible factors shaping CUB were further evaluated. Composition analysis of codons revealed that the overall GC (GCall) and GC at the third codon position (GC3) contents of mt genes were lower than 50%, presenting a preference for A/T-rich nucleotides and A/T-ending codons in H. citrina. The high values of the effective number of codons (ENC) are indicative of fairly weak CUB. Significant correlations of ENC with the GC3 and codon counts were observed, suggesting that not only compositional constraints but also gene length contributed greatly to CUB. Combined ENC-plot, neutrality plot, and Parity rule 2 (PR2)-plot analyses augmented the inference that the CUB patterns of the H. citrina mitogenome can be attributed to multiple factors. Natural selection, mutation pressure, and other factors might play a major role in shaping the CUB of mt genes, although natural selection is the decisive factor. Moreover, we identified a total of 29 high-frequency codons and 22 optimal codons, which exhibited a consistent preference for ending in A/T. Subsequent relative synonymous codon usage (RSCU)-based cluster and mt protein coding gene (PCG)-based phylogenetic analyses suggested that H. citrina is close to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in evolutionary terms, reflecting a certain correlation between CUB and evolutionary relationships.

CONCLUSIONS

There is weak CUB in the H. citrina mitogenome that is subject to the combined effects of multiple factors, especially natural selection. H. citrina was found to be closely related to Asparagus officinalis, Chlorophytum comosum, Allium cepa, and Allium fistulosum in terms of their evolutionary relationships as well as the CUB patterns of their mitogenomes. Our findings provide a fundamental reference for further studies on genetic modification and phylogenetic evolution in H. citrina.

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

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

Comparative and phylogenetic analysis of the complete chloroplast genomes of Uncaria (Rubiaceae) species

The genus Uncaria is famous for its high medicinal value. However, the high morphological similarities and unclear interspecific genetic relationships have posed challenges to the classification and identification of Uncaria species. Here, we newly sequenced six chloroplast genomes of Uncaria species: U. hirsuta, U. rhynchophylla, U. rhynchophylloides, U. homomalla, U. sinensis, and U. lancifolia. Comparisons among the chloroplast genomes of Uncaria species showed their conservation in structure, gene content, and order. Ten highly variable loci could be potentially used as specific molecular markers in the identification of Uncaria species. The third position of codons tended to use A/U base, and natural selection contributed more to the formation of codon usage bias in comparison to mutation pressure. Four genes (rbcL, ndhF, rps8, and ycf2) were detected to be subjected to positive selection. Phylogenetic analysis showed that the genus Uncaria was a monophyletic group, belonging to the tribe Naucleeae. Moreover, U. sinensis was not a variant of U. rhynchophylla. U. rhynchophylloides and U. rhynchophylla were not the same species. The results of the comparative and phylogenetic analysis provide valuable references for further research studies of classification, identification, breeding improvement, and phylogenetic relationships in Uncaria species.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

Genome-wide identification and molecular evolution of NAC gene family in Dendrobium nobile

NAC transcription factors are an important genes that regulate plant growth and development, and can regulate functions such as fruit ripening in plants. Based on genome data of Dendrobium nobile, the NAC gene family was identified and analyzed by bioinformatics methods. In this study, we identified 85 NAC genes in Dendrobium nobile genome, and systematically analyzed the NAC gene family. We found that they were distributed unevenly in the nineteen chromosomes. The amino acid length of D. nobile NAC gene family (DnoNACs) ranged from 80 to 1065, molecular weight ranged from 22.17 to 119.02 kD, and isoelectric point ranged from 4.61~9.26. Its promoter region contains multiple stress responsive elements, including light responsive, gibberellin-responsive, abscisic acid responsiveness, MeJA-responsiveness and drought-inducibility elements. Phylogenetic analysis indicates that the D. nobile NAC gene family is most closely related to Dendrobium catenatum and Dendrobium chrysotoxum. Analysis of SSR loci indicates that the fraction of mononucleotide repeats was the largest, as was the frequency of A/T. Non-coding RNA analysis showed that these 85 NAC genes contain 397 miRNAs. The collinearity analysis shows that 9 collinear locis were found on the chromosomes of D. nobile with Arabidopsis thaliana, and 75 collinear locis with D.chrysotoxum. QRT-PCR experiment under different salt concentration and temperature conditions verified the response mechanism of DnoNAC gene family under stress conditions. Most DnoNAC genes are sensitive to salt stress and temperature stress. The results of this study provide a reference for further understanding the function of NAC gene in D. nobile.

An analysis of codon utilization patterns in the chloroplast genomes of three species of Coffea

BACKGROUND

The chloroplast genome of plants is known for its small size and low mutation and recombination rates, making it a valuable tool in plant phylogeny, molecular evolution, and population genetics studies. Codon usage bias, an important evolutionary feature, provides insights into species evolution, gene function, and the expression of exogenous genes. Coffee, a key crop in the global tropical agricultural economy, trade, and daily life, warrants investigation into its codon usage bias to guide future research, including the selection of efficient heterologous expression systems for coffee genetic transformation.

RESULTS

Analysis of the codon utilization patterns in the chloroplast genomes of three Coffea species revealed a high degree of similarity among them. All three species exhibited similar base compositions, with high A/T content and low G/C content and a preference for A/T-ending codons. Among the 30 high-frequency codons identified, 96.67% had A/T endings. Fourteen codons were identified as ideal. Multiple mechanisms, including natural selection, were found to influence the codon usage patterns in the three coffee species, as indicated by ENc-GC3s mapping, PR2 analysis, and neutral analysis. Nicotiana tabacum and Saccharomyces cerevisiae have potential value as the heterologous expression host for three species of coffee genes.

CONCLUSION

This study highlights the remarkable similarity in codon usage patterns among the three coffee genomes, primarily driven by natural selection. Understanding the gene expression characteristics of coffee and elucidating the laws governing its genetic evolution are facilitated by investigating the codon preferences in these species. The findings can enhance the efficacy of exogenous gene expression and serve as a basis for future studies on coffee evolution.

Comparative Analysis of Codon Usage Patterns in Nuclear and Chloroplast Genome of Dalbergia (Fabaceae)

The Dalbergia plants are widely distributed across more than 130 tropical and subtropical countries and have significant economic and medicinal value. Codon usage bias (CUB) is a critical feature for studying gene function and evolution, which can provide a better understanding of biological gene regulation. In this study, we comprehensively analyzed the CUB patterns of the nuclear genome, chloroplast genome, and gene expression, as well as systematic evolution of Dalbergia species. Our results showed that the synonymous and optimal codons in the coding regions of both nuclear and chloroplast genome of Dalbergia preferred ending with A/U at the third codon base. Natural selection was the primary factor affecting the CUB features. Furthermore, in highly expressed genes of Dalbergia odorifera, we found that genes with stronger CUB exhibited higher expression levels, and these highly expressed genes tended to favor the use of G/C-ending codons. In addition, the branching patterns of the protein-coding sequences and the chloroplast genome sequences were very similar in the systematic tree, and different with the cluster from the CUB of the chloroplast genome. This study highlights the CUB patterns and features of Dalbergia species in different genomes, explores the correlation between CUB preferences and gene expression, and further investigates the systematic evolution of Dalbergia, providing new insights into codon biology and the evolution of Dalbergia plants.

Codon Usage Analyses Reveal the Evolutionary Patterns among Plastid Genes of Saxifragales at a Larger-Sampling Scale

Saxifragales is a 15-family order of early-divergent Eudicots with a rich morphological diversity and an ancient rapid radiation. Codon usage bias (CUB) analyses have emerged as an essential tool for understanding the evolutionary dynamics in genes. Thus far, the codon utilization patterns had only been reported in four separate genera within Saxifragales. This study provides a comprehensive assessment of the codon manipulation based on 50 plastid genes, covering 11 constituent families at a larger sampling scale. Our results first showed a high preference for AT bases and AT-ending codons. We then used effective number of codons (ENC) to assess a range of codon bias levels in the plastid genes. We also detected high-informative intrafamilial differences of ENC in three families. Subsequently, parity rule 2 (PR2) plot analyses revealed both family-unique and order-shared bias patterns. Most importantly, the ENC plots and neutrality analyses collectively supported the dominant roles of selection in the CUB of Saxifragales plastid genes. Notably, the phylogenetic affinities inferred by both ML and BI methods were consistent with each other, and they all comprised two primary clades and four subclades. These findings significantly enhance our understanding of the evolutionary processes of the Saxifrage order, and could potentially inspire more CUB analyses at higher taxonomic levels.

Comparative analysis of codon usage patterns in chloroplast genomes of ten Epimedium species

BACKGROUND

The Phenomenon of codon usage bias exists in the genomes of prokaryotes and eukaryotes. The codon usage pattern is affected by environmental factors, base mutation, gene flow and gene expression level, among which natural selection and mutation pressure are the main factors. The study of codon preference is an effective method to analyze the source of evolutionary driving forces in organisms. Epimedium species are perennial herbs with ornamental and medicinal value distributed worldwide. The chloroplast genome is self-replicating and maternally inherited which is usually used to study species evolution, gene expression and genetic transformation.

RESULTS

The results suggested that chloroplast genomes of Epimedium species preferred to use codons ending with A/U. 17 common high-frequency codons and 2-6 optimal codons were found in the chloroplast genomes of Epimedium species, respectively. According to the ENc-plot, PR2-plot and neutrality-plot, the formation of codon preference in Epimedium was affected by multiple factors, and natural selection was the dominant factor. By comparing the codon usage frequency with 4 common model organisms, it was found that Arabidopsis thaliana, Populus trichocarpa, and Saccharomyces cerevisiae were suitable exogenous expression receptors.

CONCLUSION

The evolutionary driving force in the chloroplast genomes of 10 Epimedium species probably comes from mutation pressure. Our results provide an important theoretical basis for evolutionary analysis and transgenic research of chloroplast genes.

Comprehensive Analysis of Codon Usage in Quercus Chloroplast Genome and Focus on psbA Gene

Quercus (oak) is an important economic and ecological tree species in the world, and it is the necessary feed for oak silkworm feeding. Chloroplasts play an important role in green plants but the codon usage of oak chloroplast genomes is not fully studied. We examined the codon usage of the oak chloroplast genomes in detail to facilitate the understanding of their biology and evolution. We downloaded all the protein coding genes of 26 non-redundant chloroplast reference genomes, removed short ones and those containing internal stop codons, and finally retained 50 genes shared by all genomes for comparative analyses. The base composition, codon bias, and codon preference are not significantly different between genomes but are significantly different among genes within these genomes. Oak chloroplast genomes prefer T/A-ending codons and avoid C/G-ending codons, and the psbA gene has the same preference except for the codons encoding amino acid Phe. Complex factors such as context-dependent mutations are the major factors affecting codon usage in these genomes, while selection plays an important role on the psbA gene. Our study provided an important understanding of codon usage in the oak chloroplast genomes and found that the psbA gene has nearly the same codon usage preference as other genes in the oak chloroplasts.

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

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