Codon usage patterns and evolution of HSP60 in birds

Sturnidae sensu lato Mitogenomics: Novel Insights into Codon Aversion, Selection, and Phylogeny

The Sturnidae family comprises 123 recognized species in 35 genera. The taxa Mimidae and Buphagidae were formerly treated as subfamilies within Sturnidae. The phylogenetic relationships among the Sturnidae and related taxa (Sturnidae sensu lato) remain unresolved due to high rates of morphological change and concomitant morphological homoplasy. This study presents five new mitogenomes of Sturnidae sensu lato and comprehensive mitogenomic analyses. The investigated mitogenomes exhibit an identical gene composition of 37 genes-including 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes-and one control region (CR). The most important finding of this study is drawn from CAM analyses. The surprisingly unique motifs for each species provide a new direction for the molecular species identification of avian. Furthermore, the pervasiveness of the natural selection of PCGs is found in all examined species when analyzing their nucleotide composition and codon usage. We also determine the structures of mt-tRNA, mt-rRNA, and CR structures of Sturnidae sensu lato. Lastly, our phylogenetic analyses not only well support the monophyly of Sturnidae, Mimidae, and Buphagidae, but also define nine stable subclades. Taken together, our findings will enable the further elucidation of the evolutionary relationships within Sturnidae sensu lato.

Analysis of the Mitochondrial COI Gene and Genetic Diversity of Endangered Goose Breeds

The mitochondrial cytochrome c oxidase subunit I (COI) genes of six endangered goose breeds (Xupu, Yangjiang, Yan, Wuzong, Baizi, and Lingxian) were sequenced and compared to assess the genetic diversity of endangered goose breeds. By constructing phylogenetic trees and evolutionary maps of genetic relationships, the affinities and degrees of genetic variations among the six different breeds were revealed. A total of 92 polymorphic sites were detected in the 741 bp sequence of the mtDNA COI gene after shear correction, and the GC content of the processed sequence (51.11%) was higher than that of the AT content (48.89%). The polymorphic loci within the populations of five of the six breeds (Xupu, Yangjiang, Yan, Baizi, and Lingxian) were more than 10, the haplotype diversity > 0.5, and the nucleotide diversity (Pi) > 0.005, with the Baizi geese being the exception. A total of 35 haplotypes were detected based on nucleotide variation among sequences, and the goose breed haplotypes showed a central star-shaped dispersion; the FST values were -0.03781 to 0.02645, The greatest genetic differentiation (FST = 0.02645) was observed in Yan and Wuzong breeds. The most frequent genetic exchange (Nm > 15.00) was between the Wuzong and Yangjiang geese. An analysis of molecular variance showed that the population genetic variation mainly came from within the population; the base mismatch differential distribution analysis of the goose breeds and the Tajima's D and Fu's Fs neutral detection of the historical occurrence dynamics of their populations were negative (p > 0.10). The distribution curve of the base mismatches showed a multimodal peak, which indicated that the population tended to be stabilised. These results provide important genetic information for the conservation and management of endangered goose breeds and a scientific basis for the development of effective conservation strategies.

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.

Exploring Plastomic Resources in Sempervivum (Crassulaceae): Implications for Phylogenetics

The plastid organelle is vital for photosynthesis and energy production. Advances in sequencing technology have enabled the exploration of plastomic resources, offering insights into plant evolution, diversity, and conservation. As an important group of horticultural ornamentals in the Crassulaceae family, Sempervivum plants are known for their unique rosette-like structures and reproduction through offsets. Despite their popularity, the classification status of Sempervivum remains uncertain, with only a single plastome sequence currently available. Furthermore, codon usage bias (CUB) is a widespread phenomenon of the unbalanced usage of synonymous codons in the coding sequence (CDS). However, due to the limited available plastid data, there has been no research that focused on the CUB analysis among Sempervivum until now. To address these gaps, we sequenced and released the plastomes of seven species and one subspecies from Sempervivum, revealing several consistent patterns. These included a shared 110 bp extension of the rps19 gene, 14 hypervariable regions (HVRs) with distinct nucleotide diversity (π: 0.01173 to 0.02702), and evidence of selective pressures shaping codon usage. Notably, phylogenetic analysis robustly divided the monophyletic clade into two sections: Jovibarba and Sempervivum. In conclusion, this comprehensive plastomic resource provides valuable insights into Sempervivum evolution and offers potential molecular markers for DNA barcoding.

A new perspective on codon usage, selective pressure, and phylogenetic implications of the plastomes in the Telephium clade (Crassulaceae)

The Telephium clade of the Crassulaceae family contains many medicinal, ornamental, and ecologically restorative plants. However, the phylogenetic relationships within the clade remain debated, and comprehensive analyses of codon usage and selection pressure in Telephium plastomes are limited. In this study, we assembled and annotated four plastomes and performed extensive analyses. The plastomes exhibited a typical quadripartite structure and high conservation. The lengths ranged from 151,357 bp to 151,641 bp with 134 genes identified. The GC content was the highest within IR, followed by LSC, and lowest in the SSC region. Meanwhile, a unique inversion was observed within the LSC region of Meterostachys sikokianus. Polymorphisms analysis revealed minimum nucleotide diversity in the IR regions, with over ten highly polymorphic regions identified. Phylogenetically, two subclades formed within the monophyletic Telephium clade, with Umbilicus as the sister group to the remaining Hylotelephium subclade members. Notably, no significant positive selection was found among the 79 plastid genes, which showed varying evolutionary patterns. However, 19 genes contained codons under positive selection. The specific functions of these sites require further investigation. Synonymous codon usage was biased and conserved across the tested plastomes, shaped by natural selection, mutations and other factors of varying influence. We also identified 34 taxon-specific codon aversion motifs from 49 plastid genes. Our plastomic analyses elucidate phylogenetic relationships and evolutionary patterns in this medicinal clade, providing a foundation for further research on these ecologically and pharmaceutically important plants.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

Chromosome-Level Analysis of the Pelochelys cantorii Genome Provides Insights to Its Immunity, Growth and Longevity

The Asian giant soft-shelled turtle, Pelochelys cantorii (Trionychidae), is one of the largest aquatic turtles in China and was designated as a First-Grade Protected Animal in China in 1989. Previous investigation based on a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies acquired a high-quality chromosome-level genome of Pc. cantorii. In this study, comparative genomic analysis between Pc. cantorii and 16 other vertebrate genomes indicated that turtles separated from the ancestor of archosaurians approximately 256.6 (95% highest posterior density interval, 263.6-251.9) million years ago (Mya) (Upper Permian to Triassic) and that Pc. cantorii separated from the ancestor of Pd. sinensis and R. swinhoei approximately 59.3 (95% highest posterior density interval, 64.3-54.3) Mya. Moreover, several candidate genes, such as VWA5A, ABCG2, A2M and IGSF1, associated with tumor suppression, growth and age were expanded, implicating their potential roles in the exceptional longevity of turtles. This new chromosome-level assembly has important scientific value in the study of conservation of Pc. cantorii and also enriches the evolutionary investigation of turtle species.

Genome-wide survey of catalase genes in Brassica rapa, Brassica oleracea, and Brassica napus: identification, characterization, molecular evolution, and expression profiling of BnCATs in response to salt and cadmium stress

Catalase (CAT, EC 1.11.1.6), one of the most important antioxidant enzymes, can control excess levels of H₂O₂ produced under oxidative stress in plants. In this study, 16, 8, and 7 CAT genes in the genome of Brassica napus, B. rapa, and B. oleracea were identified, respectively. Phylogenetic studies showed that CATs could be divided into two main groups, each containing specific monocotyledon and dicotyledon subgroups. Motifs, gene structure, and intron phase of CATs in B. napus, Brassica rapa, and Brassica oleracea are highly conserved. Analysis of codon usage bias showed the mutation pressure and natural selection of the codon usage of CATs. Segmental duplication and polyploid were major factors in the expansion of this gene family in B. napus, and genes have experienced  negative selection during evolution. Existence of hormones and stress-responsive cis-elements and identifying miRNA molecules affecting CATs showed that these genes are complexly regulated at the transcriptional and posttranscriptional levels. Based on RNA-seq data, CATs are divided into two groups; the first group has moderate and specific expression in flowers, leaves, stems, and roots, while the second group shows expression in most tissues. qRT-PCR analysis showed that the expression of these genes is dynamic and has a specific expression consistent with other CAT genes in response to salinity and cadmium (Cd) stresses. These results provide information for further investigation of the function of CAT genes in response to stresses and the development of tolerant  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.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.

Plastomes of Bletilla (Orchidaceae) and Phylogenetic Implications

The genus Bletilla is a small genus of only five species distributed across Asia, including B. chartacea, B. foliosa, B. formosana, B. ochracea and B. striata, which is of great medicinal importance. Furthermore, this genus is a member of the key tribe Arethuseae (Orchidaceae), harboring an extremely complicated taxonomic history. Recently, the monophyletic status of Bletilla has been challenged, and the phylogenetic relationships within this genus are still unclear. The plastome, which is rich in both sequence and structural variation, has emerged as a powerful tool for understanding plant evolution. Along with four new plastomes, this work is committed to exploring plastomic markers to elucidate the phylogeny of Bletilla. Our results reveal considerable plastomic differences between B. sinensis and the other three taxa in many aspects. Most importantly, the specific features of the IR junction patterns, novel pttRNA structures and codon aversion motifs can serve as useful molecular markers for Bletilla phylogeny. Moreover, based on maximum likelihood and Bayesian inference methods, our phylogenetic analyses based on two datasets of Arethuseae strongly imply that Bletilla is non-monophyletic. Accordingly, our findings from this study provide novel potential markers for species identification, and shed light on the evolution of Bletilla and Arethuseae.

Plastome evolution of Aeonium and Monanthes (Crassulaceae): insights into the variation of plastomic tRNAs, and the patterns of codon usage and aversion

This study reported 13 new plastomes from Aeonium and Monanthes, and observed new markers for phylogeny and DNA barcoding, such as novel tRNA structures and codon usage bias and aversion. The Macaronesian clade of Crassulaceae consists of three genera: Aichryson, with about 15 species; Monanthes, with about 10 species; Aeonium, with about 40 species. Within this clade, Aeonium, known as "the botanical equivalent of Darwin's finches", is regarded as an excellent model plant for researching adaptive evolution. Differing from the well-resolved relationships among three genera of the Macaronesian clade, the internal branching patterns within the genus Aeonium are largely unclear. In this study, we first reported 13 new plastomes from genus Aeonium and the closely related genus Monanthes. We further performed comprehensive analyses of the plastomes, with focuses on the secondary structures of pttRNAs and the patterns of codon usage and aversion. With a typical circular and quadripartite structure, the 13 plastomes ranged from 149,900 to 151,030 bp in size, and the unique pattern in IR junctions might become a family-specific marker for Crassulaceae species. Surprisingly, the π values of plastomes from Monanthes were almost twice those from Aeonium. Most importantly, we strongly recommend that highly polymorphic regions, novel putative pttRNA structures, patterns of codon usage bias and aversion derived from plastomes might have phylogenetic implications, and could act as new markers for DNA barcoding of plants. The results of phylogenetic analyses strongly supported a clear internal branching pattern in Macaronesian clade (represented by Aeonium and Monanthes), with higher nodal support values. The findings reported here will provide new insights into the variation of pttRNAs, and the patterns of codon usage and aversion of the family Crassulaceae.

Codon usage divergence of important functional genes in Mycobacterium tuberculosis

Sequence characteristics are usually used to explain the adaptive ability to hosts, metabolism, genetic diversity, drug resistance, and infectivity of Mycobacterium tuberculosis. Exploring the codon usage pattern of coding sequences in Mycobacterium tuberculosis is of great significance. In the present study, two hundred random complete genomes of Mycobacterium tuberculosis were downloaded from the National Center for Biotechnology Information database. The important codon usage pattern, such as the codon bias index, the effective number of codons, the relative synonymous codon usage as well as the base component, of twenty one specific functional genes were counted or calculated. The differences of the relative synonymous codon usage values among those functional genes, and the summation of the standard deviations of codon usage parameters were used to evaluate the divergence degree of the concerned genes. The results show that among the concerned genes, 1) all genes are high GC sequences, the codon usage frequency corresponding to each amino acid of these functional genes had a significant bias; 2) the genes of those with high effective number of codons, such as the coding sequences of Myco-bacterial membrane protein large family, usually have higher divergences; and 3) genes with lower divergences, such as the ag85A and the sigH, are usually highly conserved and are often used as drug target genes. The findings of the present work would improve new understandings on the evolution of Mycobacterium tuberculosis and on the measures to prevent and control tuberculosis from the gene engineering.

The Codon Usage Bias Analysis of Free-Living Ciliates' Macronuclear Genomes and Clustered Regularly Interspaced Short Palindromic Repeats/Cas9 Vector Construction of Stylonychia lemnae

Ciliates represent higher unicellular animals, and several species are also important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome in ciliates can not only promote a better understanding of the genetic mode and evolution history of these organisms but also help optimize codons to improve the gene editing efficiency of model ciliates. In this study, macronuclear genome sequences of nine free-living ciliates were analyzed with CodonW software to calculate the following indices: the guanine-cytosine content (GC); the frequency of the nucleotides U, C, A, and G at the third position of codons (U3s, C3s, A3s, G3s); the effective number of codons (ENC); the correlation between GC at the first and second positions (GC12); the frequency of the nucleotides G + C at the third position of synonymous codons (GC3s); the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, neutrality plot analysis, and correlation analysis were performed to explore the factors that influence codon preference. The results showed that the GC contents in nine ciliates' MAC genomes were lower than 50% and appeared AT-rich. The base compositions of GC12 and GC3s are markedly distinct and the codon usage pattern and evolution of ciliates are affected by genetic mutation and natural selection. According to the synonymous codon analysis, the codons of most ciliates ended with A or U and eight codons were the general optimal codons of nine ciliates. A clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9) expression vector of Stylonychia lemnae was constructed by optimizing the macronuclear genome codon and was successfully used to knock out the Adss gene. This is the first such extensive investigation of the MAC genome CUB of ciliates and the initial successful application of the CRISPR/Cas9 technique in free-living ciliates.