Variation and selection on codon usage bias across an entire subphylum

The evolution and convergence of mutation spectra across mammals

Despite the key role genetic mutations play in shaping phenotypic differences between species, little is currently known about the evolution of germline mutation spectra across mammals. Domesticated species are likely particularly interesting case studies because of their high mutation rates and complex evolutionary histories, which can span multiple founding events and genetic bottlenecks. Here we have developed a new reusable workflow, nSPECTRa, that can undertake the key steps in characterising mutation spectra, from determining ancestral alleles to characterising multiple forms of variation. We apply nSPECTRa to seven species, including several that have undergone domestication, and highlight how nSPECTRa can provide important insights into mutation rate evolution. While mutation spectra most often show marked differences between species and even breeds, certain mutation types have risen to a high frequency in subpopulations of different species, indicative of convergent evolution in mutation rates. This includes the previously characterized TCC- > TTC change enriched among European humans, which is also enriched among East Asian cattle. We show Indicine cattle are particularly interesting examples of how different mutation spectra segregate within a population and subsequently spread across the globe. Together, this work has important implications for understanding the control and evolution of mammalian mutation rates.

The complete mitochondrial genome of Sinojackia microcarpa: evolutionary insights and gene transfer

BACKGROUND

As a dicotyledonous plant within the Styracaceae family, Sinojackia microcarpa (S. microcarpa) is notable for its library-shaped fruit and sparse distribution, serving as a model system for studying the entire tree family. However, the scarcity of genomic data, particularly concerning the mitochondrial and nuclear sequences of S. microcarpa, has substantially impeded our understanding of its evolutionary traits and fundamental biological mechanisms.

RESULTS

This study presents the first complete mitochondrial genome sequence of S. microcarpa and conducts a comparative analysis of its protein-encoding genes across eight plant species. Our analysis revealed that the mitochondrial genome of S. microcarpa spans 687,378 base pairs and contains a total of 59 genes, which include 37 protein-coding genes (PCGs), 20 transfer RNA (tRNA) genes, and 2 ribosomal RNA (rRNA) genes. Sixteen plastid-derived fragments strongly linked with mitochondrial genes, including one intact plastid-related gene (rps7), were identified. Additionally, Ka/Ks ratio analysis revealed that most mitochondrial genes are under purifying selection, with a few genes, such as nad9 and ccmB, showing signs of relaxed or adaptive evolution. An analysis of twenty-nine protein-coding genes from twenty-four plant species reveals that S. microcarpa exhibits a closer evolutionary relationship with species belonging to the genus Camellia. The findings of this study provide new genomic data that enhance our understanding of S. microcarpa, and reveal its mitochondrial genome's evolutionary proximity to other dicotyledonous species.

CONCLUSIONS

Overall, this research enhances our understanding of the evolutionary and comparative genomics of S. microcarpa and other plants in the Styracaceae family and lays the foundation for future genetic studies and evolutionary analyses in the Styracaceae family.

Macroevolutionary changes in natural selection on codon usage reflects evolution of the tRNA pool across a budding yeast subphylum

Across taxonomical domains, synonymous codons of an amino acid are found to be used at unequal frequencies within genes. This codon usage bias (CUB) is highly variable across species. Genome-wide CUB reflects a balance between adaptive and non-adaptive microevolutionary processes within a species. Variation in microevolutionary processes results in across-species variation in CUB. As CUB is tightly linked to important molecular and biophysical processes, it is critical to understand how changes to these processes are linked to changes in microevolutionary processes. We employed a population genetics model to quantify natural selection and mutation biases on a per-codon basis across the Saccharomycotina budding yeast subphylum. We found that the strength of natural selection and mutation biases varied significantly between closely related yeasts. Across-species variation in natural selection reflected the evolution of tRNA gene copy number. Additionally, we found evidence that changes to tRNA modification expression can contribute to changes in natural selection across species independent of changes to tGCN. Both lines of evidence support the link between the evolution of the tRNA pool and natural selection in codon usage through changes in the translation efficiency of a codon. Furthermore, we found that changes to tGCN often reflected changes to genome-wide GC%, suggesting changes to the tRNA pool reflect changes to mutation bias. Our work establishes how changes in microevolutionary processes impact changes in molecular mechanisms, ultimately shaping the macroevolutionary variation of a trait.

Significance statement

Codon usage bias (CUB) - the non-uniform usage of synonymous codons - is a feature of all genomes and varies across closely related species. Differences in CUB imply differences in the underlying microevolutionary processes (natural selection, mutation bias) driving CUB. CUB is hypothesized to be tightly linked to key molecular processes, particularly mRNA translation. We used a population genetics model to quantify natural selection and mutation bias on a per-codon basis across 327 budding yeasts. We found high variability in the microevolution of CUB and showed that changes in natural selection were correlated with the evolution of the tRNA pool. Our work establishes how variation in molecular mechanisms relates to variation in microevolution, shaping variation in a trait across species.

Lost in translation: conserved amino acid usage despite extreme codon bias in foraminifera

Analyses of codon usage in eukaryotes suggest that amino acid usage responds to GC pressure so AT-biased substitutions drive higher usage of amino acids with AT-ending codons. Here, we combine single-cell transcriptomics and phylogenomics to explore codon usage patterns in foraminifera, a diverse and ancient clade of predominantly uncultivable microeukaryotes. We curate data from 1,044 gene families in 49 individuals representing 28 genera, generating perhaps the largest existing dataset of data from a predominantly uncultivable clade of protists, to analyze compositional bias and codon usage. We find extreme variation in composition, with a median GC content at fourfold degenerate silent sites below 3% in some species and above 75% in others. The most AT-biased species are distributed among diverse non-monophyletic lineages. Surprisingly, despite the extreme variation in compositional bias, amino acid usage is highly conserved across all foraminifera. By analyzing nucleotide, codon, and amino acid composition within this diverse clade of amoeboid eukaryotes, we expand our knowledge of patterns of genome evolution across the eukaryotic tree of life.IMPORTANCEPatterns of molecular evolution in protein-coding genes reflect trade-offs between substitution biases and selection on both codon and amino acid usage. Most analyses of these factors in microbial eukaryotes focus on model species such as Acanthamoeba, Plasmodium, and yeast, where substitution bias is a primary contributor to patterns of amino acid usage. Foraminifera, an ancient clade of single-celled eukaryotes, present a conundrum, as we find highly conserved amino acid usage underlain by divergent nucleotide composition, including extreme AT-bias at silent sites among multiple non-sister lineages. We speculate that these paradoxical patterns are enabled by the dynamic genome structure of foraminifera, whose life cycles can include genome endoreplication and chromatin extrusion.

Oleaginous Yeast Biology Elucidated With Comparative Transcriptomics

Extremophilic yeasts have favorable metabolic and tolerance traits for biomanufacturing- like lipid biosynthesis, flavinogenesis, and halotolerance - yet the connection between these favorable phenotypes and strain genotype is not well understood. To this end, this study compares the phenotypes and gene expression patterns of biotechnologically relevant yeasts Yarrowia lipolytica, Debaryomyces hansenii, and Debaryomyces subglobosus grown under nitrogen starvation, iron starvation, and salt stress. To analyze the large data set across species and conditions, two approaches were used: a "network-first" approach where a generalized metabolic network serves as a scaffold for mapping genes and a "cluster-first" approach where unsupervised machine learning co-expression analysis clusters genes. Both approaches provide insight into strain behavior. The network-first approach corroborates that Yarrowia upregulates lipid biosynthesis during nitrogen starvation and provides new evidence that riboflavin overproduction in Debaryomyces yeasts is overflow metabolism that is routed to flavin cofactor production under salt stress. The cluster-first approach does not rely on annotation; therefore, the coexpression analysis can identify known and novel genes involved in stress responses, mainly transcription factors and transporters. Therefore, this work links the genotype to the phenotype of biotechnologically relevant yeasts and demonstrates the utility of complementary computational approaches to gain insight from transcriptomics data across species and conditions.

Complete mitogenome characteristics and phylogenetic analysis of traditional Chinese medicinal plant Tinospora sagittata (Oliv.) Gagnep. from the Menispermaceae family

BACKGROUND

Tinospora sagittata, a member belongs to the genus Tinospora of Menispermaceae family. Its tuberous roots have been used as traditional Chinese medicine (TCM) for pharmacological properties and are commonly known name as "Jin Guo Lan". Although its plastome and nuclear genome had been sequenced, its mitochondrial genome has not been explored, which significantly hampers conservation efforts and further research for this species. In addition, previous efforts based on multiple molecular markers providing profound insights into an intergeneric phylogenetic framework for Burasaieae and sampled species of T. sagittata are placed in a superclades, species delimitation of T. sagittata still need to be comprehensively evaluated.

RESULTS

Flow cytometry revealed that Tinospora sagittata has two cytotypes and a wide range in genome sizes. We further sequenced and assembled the organelle genomes of T. sagittata, including the mitogenome (513,210-513,215 bp) and plastome (163,621-164,006 bp). The plastomes were highly similar in gene content and exhibited a typical quadripartite structure, but a translocation as well as two inversions were detected in mitogenomes. The repeats patterns in both organelles are generally similar, but significant difference in the codon bias of the genes of Tinospora organelle genomes. Interesting, both organelle genomes had shown that inter-gene spacer regions could be used as effective molecular markers for further phylogenetic analyses and species identification. Comprehensive analysis of protein coding genes of organelle genomes showed that significant difference in Ka, Ks, and Ka/Ks values among the organelle genomes. Phylogenetic analysis identified a tree that was basically consistent with the phylogeny of Ranunculales described in the APG IV system.

CONCLUSIONS

We provided a high-quality and well-annotated organelle genome for Tinospora sagittata. The study present here advances our understanding of the intricate interplay between plastome and mitogenome. Moreover, our results also laid the foundation for further studying the course, tempo and mode of organelle genome evolution of Menispermaceae.

Complete mitochondrial genome assembly and comparative analysis of Colocasia esculenta

Colocasia esculenta ranks as the fifth most important tuber crop and is known for its high nutritional and medicinal value. However, there is no research on its mitochondrial genome, hindering in-depth exploration of its genomic resources and genetic relationships. Using second- and third-generation sequencing technologies, we assembled and annotated the mitogenome of C. esculenta. Its mitogenome mainly consists of five circular DNA molecules, with a total length of 594,811 bp and a GC content of 46.25%. A total of 55 genes, 157 simple sequence repeats, 29 tandem repeat sequences, 202 dispersed repeat sequences, and 625 RNA editing sites were detected. Most protein-coding genes use ATG as the start codon, and the third position of the codon tends to be A or T (U). GAA, AUU, and UUU are the most common codons in C. esculenta mitochondria. Finally, based on 28 representative plant species, a phylogenetic tree was constructed, revealing a close relationship between C. esculenta and Araceae. This study provides comprehensive information on C. esculenta, laying a foundation for crop genetics and molecular breeding.

Comparative genomics of thermosynechococcaceae and thermostichaceae: insights into codon usage bias

Members of the families Thermosynechococcaceae and Thermostichaceae are well-known unicellular thermophilic cyanobacteria and a non-thermophilic genus Pseudocalidococcus was newly classified into the former. Analysis of the codon usage bias (CUB) of cyanobacterial species inhabiting different thermal and non-thermal niches will benefit the understanding of their genetic and evolutionary characteristics. Herein, the CUB and codon context patterns of protein-coding genes were systematically analyzed and compared between members of the two families. Overall, the nucleotide composition and CUB indices were found to differ between thermophiles and non-thermophiles. The thermophiles showed a higher G/C content in the codon base composition and tended to end with G/C compared to the non-thermophiles. Correlation analysis indicated significant associations between codon base composition and CUB indices. The results of the effective number of codons, parity-rule 2, neutral and correspondence analyses indicated that mutational pressure and natural selection primarily account for CUB in these cyanobacterial species, but the primary driving forces exhibit variation among genera. Moreover, the optimal codons identified based on relative synonymous codon usage values were found to differ among genera and even within genera. In addition, codon context pattern analysis revealed the specificity of the sequence context of start and stop codons among genera. Intriguingly, the clustering of codon context patterns appeared to be more related to thermotolerance than to phylogenomic relationships. In conclusion, this study facilitates the understanding of the characteristics and sources of variation of CUB and the evolution of the surveyed cyanobacterial clades with different thermotolerance and provides insights into their adaptation to different environments.

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. It has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. We characterized codon usage across 1,154 strains from 1,051 species from the fungal subphylum Saccharomycotina to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns. We found a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is distinct between the 12 orders to such a degree that yeasts can be classified with an accuracy >90% using a machine learning algorithm. We also characterized the degree to which codon usage bias is impacted by translational selection. We found it was influenced by a combination of features, including the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that avoidance of the CGN codons is associated with a decline in arginine tRNA function. Consistent with previous findings, codon usage bias within the Saccharomycotina is shaped by genomic features and GC bias. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Selection on synonymous codon usage in soybean (Glycine max) WRKY genes

The WRKY transcription factor gene family in soybean [Glycine max (L.) Merr.] (GmWRKY) is critical for the plant's development and stress responses. This study examines the evolutionary dynamics of the GmWRKY gene family, focusing on its synonymous codon usage bias (CUB) in a comprehensive set of 179 coding sequences. CUB was analyzed using various indices, revealing a preference for A/T-ending codons and relatively low codon bias. Codon adaptation index (CAI) analysis suggested that these genes are optimized for efficient translation despite relatively low bias, reflecting a balance between codon diversity and translation efficiency. Neutrality and NC plots indicated that selective forces dominate over mutational forces in shaping codon usage, while selection signature analysis showed purifying selection being prevalent across the gene family. However, episodic positive selection was also detected in certain clades, highlighting potential adaptive diversification in response to environmental stress. Additionally, promoter binding site analysis uncovered correlations between codon usage and transcriptional regulation, indicating a context-dependent relationship between CUB and gene expression. Phylogenetic analysis identified 11 well-supported clades in the modern GmWRKY gene family and ancestral sequence reconstruction revealed more relaxed codon preferences and reduced selection constraints in modern GmWRKY genes, potentially linked to neofunctionalization and adaptation to environmental changes. These findings provide a framework for optimizing gene expression in transgenic soybean crops with resilience. Further functional validation of positively selected genes is recommended to elucidate their role in stress responses.

Genome-wide patterns of selection-drift variation strongly associate with organismal traits across the green plant lineage

There are many gaps in our knowledge of how life cycle variation and organismal body architecture associate with molecular evolution. Using the diverse range of green algal body architectures and life cycle types as a test case, we hypothesize that increases in cytomorphological complexity are likely to be associated with a decrease in the effective population size, because larger-bodied organisms typically have smaller populations, resulting in increased drift. For life cycles, we expect haploid-dominant lineages to evolve under stronger selection intensity relative to diploid-dominant life cycles owing to masking of deleterious alleles in heterozygotes. We use a genome-scale data set spanning the phylogenetic diversity of green algae and phylogenetic comparative approaches to measure the relative selection intensity across different trait categories. We show stronger signatures of drift in lineages with more complex body architectures compared with unicellular lineages, which we consider to be a consequence of smaller effective population sizes of the more complex algae. Significantly higher rates of synonymous as well as nonsynonymous substitutions relative to other algal body architectures highlight that siphonous and siphonocladous body architectures, characteristic of many green seaweeds, form an interesting test case to study the potential impacts of genome redundancy on molecular evolution. Contrary to expectations, we show that levels of selection efficacy do not show a strong association with life cycle types in green algae. Taken together, our results underline the prominent impact of body architecture on the molecular evolution of green algal genomes.

Complete mitochondrial genome of Agropyron cristatum reveals gene transfer and RNA editing events

BACKGROUND

As an important forage in arid and semi-arid regions, Agropyron cristatum provides livestock with exceptionally high nutritional value. Additionally, A. cristatum exhibits outstanding genetic characteristics to endure drought and disease. Therefore, rich genetic diversity serves as a cornerstone for the improvement of major food crops. The purposes of this study were to systematically describe mitogenome of A.cristatum and preliminarily analyze its internal variations.

RESULT

The A. cristatum mitogenome was a single-ring molecular structure of 381,065 bp that comprised 52 genes, including 35 protein-coding, 3 rRNA and 14 tRNA genes. Among these, two pseudoprotein-coding genes and multiple copies of tRNA genes were observed. A total of 320 repetitive sequences was found to cover more than 10% of the mitogenome (105 simple sequences, 185 dispersed and 30 tandem repeats), which led to a large number of fragment rearrangements in the mitogenome of A. cristatum. Leucine was the most frequent amino acid (n = 1087,10.8%) in the protein-coding genes of A. cristatum mitogenome, and the highest usage codon was ATG (initiation codon). The number of A/T changes at the third base of the codon was much higher than that of G/C. Among 23 PCGs, the range of Pi values is from 0.0021 to 0.0539, with an average of 0.013. Additionally, 81 RNA editing sites were predicted, which were considerably fewer than those reported in other plant mitogenomes. Most of the RNA editing site base positions were concentrated at the first and second codon bases, which were C to T transitions. Moreover, we identified 95 sequence fragments (total length of 34, 343 bp) that were transferred from the chloroplast to mitochondria genes, introns, and intergenic regions. The stability of the tRNA genes was maintained during this process. Selection pressure analysis of 23 protein-coding genes shared by 15 Poaceae plants, showed that most genes were subjected to purifying selection during evolution, whereas rps4, cob, mttB, and ccmB underwent positive selection in different plants. Finally, a phylogenetic tree was constructed based on 22 plant mitogenomes, which showed that Agropyron plants have a high degree of independent heritability in Triticeae.

CONCLUSION

The findings of this study provide new data for a better understanding of A. cristatum genes, and demonstrate that mitogenomes are suitable for the study of plant classifications, such as those of Agropyron. Moreover, it provides a reference for further exploration of the phylogenetic relationships within Agropyron species, and establishes a theoretical basis for the subsequent development and utilization of A. cristatum plant germplasm resources.

Plastome evolution of Engelhardia facilitates phylogeny of Juglandaceae

BACKGROUND

Engelhardia (Juglandaceae) is a genus of significant ecological and economic importance, prevalent in the tropics and subtropics of East Asia. Although previous efforts based on multiple molecular markers providing profound insights into species delimitation and phylogeography of Engelhardia, the maternal genome evolution and phylogeny of Engelhardia in Juglandaceae still need to be comprehensively evaluated. In this study, we sequenced plastomes from 14 samples of eight Engelhardia species and the outgroup Rhoiptelea chiliantha, and incorporated published data from 36 Juglandaceae and six outgroup species to test phylogenetic resolution. Moreover, comparative analyses of the plastomes were conducted to investigate the plastomes evolution of Engelhardia and the whole Juglandaceae family.

RESULTS

The 13 Engelhardia plastomes were highly similar in genome size, gene content, and order. They exhibited a typical quadripartite structure, with lengths from 161,069 bp to 162,336 bp. Three mutation hotspot regions (TrnK-rps16, ndhF-rpl32, and ycf1) could be used as effective molecular markers for further phylogenetic analyses and species identification. Insertion and deletion (InDels) may be an important driving factor for the evolution of plastomes in Juglandoideae and Engelhardioideae. A total of ten codons were identified as the optimal codons in Juglandaceae. The mutation pressure mostly contributed to shaping codon usage. Seventy-eight protein-coding genes in Juglandaceae experienced relaxed purifying selection, only rpl22 and psaI genes showed positive selection (Ka/Ks > 1). Phylogenetic results fully supported Engelhardia as a monophyletic group including two sects and the division of Juglandaceae into three subfamilies. The Engelhardia originated in the Late Cretaceous and diversified in the Late Eocene, and Juglandaceae originated in the Early Cretaceous and differentiated in Middle Cretaceous. The phylogeny and divergence times didn't support rapid radiation occurred in the evolution history of Engelhardia.

CONCLUSION

Our study fully supported the taxonomic treatment of at the section for Engelhardia species and three subfamilies for Juglandaceae and confirmed the power of phylogenetic resolution using plastome sequences. Moreover, our results also laid the foundation for further studying the course, tempo and mode of plastome evolution of Engelhardia and the whole Juglandaceae family.

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. The biased use of synonymous codons has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. The Saccharomycotina, the fungal subphylum containing the yeasts Saccharomyces cerevisiae and Candida albicans , has been a model system for studying codon usage. We characterized codon usage across 1,154 strains from 1,051 species to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns across the subphylum. We found evidence of a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is also distinct between the 12 orders within the subphylum to such a degree that yeasts can be classified into orders with an accuracy greater than 90% using a machine learning algorithm trained on codon usage. We also characterized the degree to which codon usage bias is impacted by translational selection. Interestingly, the degree of translational selection was influenced by a combination of genome features and assembly metrics that included the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs. The order contains 24 species, and 23 are computationally predicted to lack tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that extreme avoidance of the CGN codons is associated with a decline in arginine tRNA function. Codon usage bias within the Saccharomycotina is generally consistent with previous investigations in fungi, which show a role for both genomic features and GC bias in shaping codon usage. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

The Patterns of Codon Usage between Chordates and Arthropods are Different but Co-evolving with Mutational Biases

Different frequencies amongst codons that encode the same amino acid (i.e. synonymous codons) have been observed in multiple species. Studies focused on uncovering the forces that drive such codon usage showed that a combined effect of mutational biases and translational selection works to produce different frequencies of synonymous codons. However, only few have been able to measure and distinguish between these forces that may leave similar traces on the coding regions. Here, we have developed a codon model that allows the disentangling of mutation, selection on amino acids and synonymous codons, and GC-biased gene conversion (gBGC) which we employed on an extensive dataset of 415 chordates and 191 arthropods. We found that chordates need 15 more synonymous codon categories than arthropods to explain the empirical codon frequencies, which suggests that the extent of codon usage can vary greatly between animal phyla. Moreover, methylation at CpG sites seems to partially explain these patterns of codon usage in chordates but not in arthropods. Despite the differences between the two phyla, our findings demonstrate that in both, GC-rich codons are disfavored when mutations are GC-biased, and the opposite is true when mutations are AT-biased. This indicates that selection on the genomic coding regions might act primarily to stabilize its GC/AT content on a genome-wide level. Our study shows that the degree of synonymous codon usage varies considerably among animals, but is likely governed by a common underlying dynamic.

Complete chloroplast genomes of 13 species of sect. Tuberculata Chang (Camellia L.): genomic features, comparative analysis, and phylogenetic relationships

BACKGROUND

Sect. Tuberculata belongs to Camellia, and its members are characterized by a wrinkled pericarp and united filaments. All the plants in this group, which are endemic to China, are highly valuable for exploring the evolution of Camellia and have great potential for use as an oil source. However, due to the complex and diverse phenotypes of these species and the difficulty of investigating them in the field, their complex evolutionary history and interspecific definitions have remained largely unelucidated.

RESULTS

Therefore, we newly sequenced and annotated 12 chloroplast (cp) genomes and retrieved the published cp genome of Camellia anlungensis Chang in sect. Tuberculata. In this study, comparative analysis of the cp genomes of the thirteen sect. Tuberculata species revealed a typical quadripartite structure characterized by a total sequence length ranging from 156,587 bp to 157,068 bp. The cp.genome arrangement is highly conserved and moderately differentiated. A total of 130 to 136 genes specific to the three types were identified by annotation, including protein-coding genes (coding sequences (CDSs)) (87-91), tRNA genes (35-37), and rRNA genes (8). The total observed frequency ranged from 23,045 (C. lipingensis) to 26,557 (C. anlungensis). IR region boundaries were analyzed to show that the ycf1 gene of C. anlungensis is located in the IRb region, while the remaining species are present only in the IRa region. Sequence variation in the SSC region is greater than that in the IR region, and most protein-coding genes have high codon preferences. Comparative analyses revealed six hotspot regions (tRNA-Thr(GGT)-psbD, psbE-petL, ycf15-tRNA-Leu(CAA), ndhF-rpl32, ndhD, and trnL(CAA)-ycf15) in the cp genomes that could serve as potential molecular markers. In addition, the results of phylogenetic tree construction based on the cp genomes showed that the thirteen sect. Tuberculata species formed a monophyletic group and were divided into two evolutionarily independent clades, confirming the independence of the section.

CONCLUSIONS

In summary, we obtained the cp genomes of thirteen sect. Tuberculata plants and performed the first comparative analysis of this group. These results will help us better characterize the plants in this section, deepen our understanding of their genetic characteristics and phylogenetic relationships, and lay the theoretical foundation for their accurate classification, elucidation of their evolutionary changes, and rational development and utilization of this section in the future.

Multi-Chromosomal mitochondrial genome of medicinal plant Acorus tatarinowii (Acoraceae): Firstly reported from Acorales Order

Acorus tatarinowii Schott (A. tatarinowii), a well-known traditional Chinese medicinal plant renowned for its high medicinal value, but its mitochondrial genome (mitogenome) is still unexplored. In this study, we meticulously assembled the complete mitochondrial genome of A. tatarinowii using a combination of Illumina short reads and Oxford Nanopore long reads. Our findings revealed that A. tatarinowii possesses a complex chromosomal structural mitogenome, comprising two linear chromosomes and seven circular chromosomes. This mitogenome spans 1.81 Mb in length with a GC content of 38.29 %. Notably, it contained 24 unique mitochondrial core genes, seven unique variable genes, 17 tRNA genes, and three rRNA genes. Analyses of codon usage, most protein-coding genes (PCGs) exhibited a common codon usage preference, with RSCU values greater than 1, and the codon with the highest RSCU value was UAA(End, 1.90). We conducted a thorough analysis of repeat sequences, the distribution of repetitive sequences in nine mitochondrial chromosomes showed distinct patterns. Moreover, we identified 82 and 12 homologous fragments by comparing the sequences of chloroplast and nuclear genomes to the A. tatarinowii mitogenome, respectively. Lastly, We predicted a total of 234 potential RNA editing sites in 28 unique PCGs and discovered that the nad4 gene has been edited the most often, at 26 times. Our results contribute to the enrichment of mitochondrial genome resources for Acoraceae, and the mitogenome also can be used as a reference for other species.

Natural trait variation across Saccharomycotina species

Among molecular biologists, the group of fungi called Saccharomycotina is famous for its yeasts. These yeasts in turn are famous for what they have in common-genetic, biochemical, and cell-biological characteristics that serve as models for plants and animals. But behind the apparent homogeneity of Saccharomycotina species lie a wealth of differences. In this review, we discuss traits that vary across the Saccharomycotina subphylum. We describe cases of bright pigmentation; a zoo of cell shapes; metabolic specialties; and species with unique rules of gene regulation. We discuss the genetics of this diversity and why it matters, including insights into basic evolutionary principles with relevance across Eukarya.

Systematic analysis of the codon usage patterns of African swine fever virus genome coding sequences reveals its host adaptation phenotype

African swine fever (ASF) is a severe haemorrhagic disease caused by the African swine fever virus (ASFV), transmitted by ticks, resulting in high mortality among domestic pigs and wild boars. The global spread of ASFV poses significant economic threats to the swine industry. This study employs diverse analytical methods to explore ASFV's evolution and host adaptation, focusing on codon usage patterns and associated factors. Utilizing phylogenetic analysis methods including neighbour-joining and maximum-likelihood, 64 ASFV strains were categorized into four clades. Codon usage bias (CUB) is modest in ASFV coding sequences. This research identifies multiple factors - such as nucleotide composition, mutational pressures, natural selection and geographical diversity - contributing to the formation of CUB in ASFV. Analysis of relative synonymous codon usage reveals CUB variations within clades and among ASFVs and their hosts. Both Codon Adaptation Index and Similarity Index analyses confirm that ASFV strains are highly adapted to soft ticks (Ornithodoros moubata) but less so to domestic pigs, which could be a result of the long-term co-evolution of ASFV with ticks. This study sheds light on the factors influencing ASFV's codon usage and fitness dynamics, enriching our understanding of its evolution, adaptation and host interactions.

Degradation of polypropylene by fungi Coniochaeta hoffmannii and Pleurostoma richardsiae

The urgent need for better disposal and recycling of plastics has motivated a search for microbes with the ability to degrade synthetic polymers. While microbes capable of metabolizing polyurethane and polyethylene terephthalate have been discovered and even leveraged in enzymatic recycling approaches, microbial degradation of additive-free polypropylene (PP) remains elusive. Here we report the isolation and characterization of two fungal strains with the potential to degrade pure PP. Twenty-seven fungal strains, many isolated from hydrocarbon contaminated sites, were screened for degradation of commercially used textile plastic. Of the candidate strains, two identified as Coniochaeta hoffmannii and Pleurostoma richardsiae were found to colonize the plastic fibers using scanning electron microscopy (SEM). Further experiments probing degradation of pure PP films were performed using C. hoffmannii and P. richardsiae and analyzed using SEM, Raman spectroscopy and Fourier transform infrared spectroscopy with attenuated total reflectance (FTIR-ATR). The results showed that the selected fungi were active against pure PP, with distinct differences in the bonds targeted and the degree to which each was altered. Whole genome and transcriptome sequencing was conducted for both strains and the abundance of carbohydrate active enzymes, GC content, and codon usage bias were analyzed in predicted proteomes for each. Enzymatic assays were conducted to assess each strain's ability to degrade naturally occurring compounds as well as synthetic polymers. These investigations revealed potential adaptations to hydrocarbon-rich environments and provide a foundation for further investigation of PP degrading activity in C. hoffmannii and P. richardsiae.

Informatic Capabilities of Translation and Its Implications for the Origins of Life

The ability to encode and convert heritable information into molecular function is a defining feature of life as we know it. The conversion of information into molecular function is performed by the translation process, in which triplets of nucleotides in a nucleic acid polymer (mRNA) encode specific amino acids in a protein polymer that folds into a three-dimensional structure. The folded protein then performs one or more molecular activities, often as one part of a complex and coordinated physiological network. Prebiotic systems, lacking the ability to explicitly translate information between genotype and phenotype, would have depended upon either chemosynthetic pathways to generate its components-constraining its complexity and evolvability- or on the ambivalence of RNA as both carrier of information and of catalytic functions-a possibility which is still supported by a very limited set of catalytic RNAs. Thus, the emergence of translation during early evolutionary history may have allowed life to unmoor from the setting of its origin. The origin of translation machinery also represents an entirely novel and distinct threshold of behavior for which there is no abiotic counterpart-it could be the only known example of computing that emerged naturally at the chemical level. Here we describe translation machinery's decoding system as the basis of cellular translation's information-processing capabilities, and the four operation types that find parallels in computer systems engineering that this biological machinery exhibits.

Natural variation in codon bias and mRNA folding strength interact synergistically to modify protein expression in Saccharomyces cerevisiae

Codon bias and mRNA folding strength (mF) are hypothesized molecular mechanisms by which polymorphisms in genes modify protein expression. Natural patterns of codon bias and mF across genes as well as effects of altering codon bias and mF suggest that the influence of these 2 mechanisms may vary depending on the specific location of polymorphisms within a transcript. Despite the central role codon bias and mF may play in natural trait variation within populations, systematic studies of how polymorphic codon bias and mF relate to protein expression variation are lacking. To address this need, we analyzed genomic, transcriptomic, and proteomic data for 22 Saccharomyces cerevisiae isolates, estimated protein accumulation for each allele of 1,620 genes as the log of protein molecules per RNA molecule (logPPR), and built linear mixed-effects models associating allelic variation in codon bias and mF with allelic variation in logPPR. We found that codon bias and mF interact synergistically in a positive association with logPPR, and this interaction explains almost all the effects of codon bias and mF. We examined how the locations of polymorphisms within transcripts influence their effects and found that codon bias primarily acts through polymorphisms in domain-encoding and 3' coding sequences, while mF acts most significantly through coding sequences with weaker effects from untranslated regions. Our results present the most comprehensive characterization to date of how polymorphisms in transcripts influence protein expression.

Characterization of the complete mitochondrial genomes of five hard ticks and phylogenetic implications

Ticks are blood-sucking ectoparasites with significant medical and veterinary importance, capable of transmitting bacteria, protozoa, fungi, and viruses that cause a variety of human and animal diseases worldwide. In the present study, we sequenced the complete mitochondrial (mt) genomes of five hard tick species and analyzed features of their gene contents and genome organizations. The complete mt genomes of Haemaphysalis verticalis, H. flava, H. longicornis, Rhipicephalus sanguineus and Hyalomma asiaticum were 14855 bp, 14689 bp, 14693 bp, 14715 bp and 14722 bp in size, respectively. Their gene contents and arrangements are the same as those of most species of metastriate Ixodida, but distinct from species of genus Ixodes. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes with two different computational algorithms (Bayesian inference and maximum likelihood) revealed the monophylies of the genera Rhipicephalus, Ixodes and Amblyomma, however, rejected the monophyly of the genus Haemaphysalis. To our knowledge, this is the first report of the complete mt genome of H. verticalis. These datasets provide useful mtDNA markers for further studies of the identification and classification of hard ticks.

Codon Optimization Improves the Prediction of Xylose Metabolism from Gene Content in Budding Yeasts

Xylose is the second most abundant monomeric sugar in plant biomass. Consequently, xylose catabolism is an ecologically important trait for saprotrophic organisms, as well as a fundamentally important trait for industries that hope to convert plant mass to renewable fuels and other bioproducts using microbial metabolism. Although common across fungi, xylose catabolism is rare within Saccharomycotina, the subphylum that contains most industrially relevant fermentative yeast species. The genomes of several yeasts unable to consume xylose have been previously reported to contain the full set of genes in the XYL pathway, suggesting the absence of a gene-trait correlation for xylose metabolism. Here, we measured growth on xylose and systematically identified XYL pathway orthologs across the genomes of 332 budding yeast species. Although the XYL pathway coevolved with xylose metabolism, we found that pathway presence only predicted xylose catabolism about half of the time, demonstrating that a complete XYL pathway is necessary, but not sufficient, for xylose catabolism. We also found that XYL1 copy number was positively correlated, after phylogenetic correction, with xylose utilization. We then quantified codon usage bias of XYL genes and found that XYL3 codon optimization was significantly higher, after phylogenetic correction, in species able to consume xylose. Finally, we showed that codon optimization of XYL2 was positively correlated, after phylogenetic correction, with growth rates in xylose medium. We conclude that gene content alone is a weak predictor of xylose metabolism and that using codon optimization enhances the prediction of xylose metabolism from yeast genome sequence data.

Analysis of synonymous codon usage patterns in mitochondrial genomes of nine Amanita species

Introduction

Codon basis is a common and complex natural phenomenon observed in many kinds of organisms.

Methods

In the present study, we analyzed the base bias of 12 mitochondrial core protein-coding genes (PCGs) shared by nine Amanita species.

Results

The results showed that the codons of all Amanita species tended to end in A/T, demonstrating the preference of mitochondrial codons of Amanita species for a preference for this codon. In addition, we detected the correlation between codon base composition and the codon adaptation index (CAI), codon bias index (CBI), and frequency of optimal codons (FOP) indices, indicating the influence of base composition on codon bias. The average effective number of codons (ENC) of mitochondrial core PCGs of Amanita is 30.81, which is <35, demonstrating the strong codon preference of mitochondrial core PCGs of Amanita. The neutrality plot analysis and PR2-Bias plot analysis further demonstrated that natural selection plays an important role in Amanita codon bias. In addition, we obtained 5–10 optimal codons (ΔRSCU > 0.08 and RSCU > 1) in nine Amanita species, and GCA and AUU were the most widely used optimal codons. Based on the combined mitochondrial sequence and RSCU value, we deduced the genetic relationship between different Amanita species and found large variations between them.

Discussion

This study promoted the understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Genes enriched in A/T-ending codons are co-regulated and conserved across mammals

Codon usage influences gene expression distinctly depending on the cell context. Yet, the importance of codon bias in the simultaneous turnover of specific groups of protein-coding genes remains to be investigated. Here, we find that genes enriched in A/T-ending codons are expressed more coordinately in general and across tissues and development than those enriched in G/C-ending codons. tRNA abundance measurements indicate that this coordination is linked to the expression changes of tRNA isoacceptors reading A/T-ending codons. Genes with similar codon composition are more likely to be part of the same protein complex, especially for genes with A/T-ending codons. The codon preferences of genes with A/T-ending codons are conserved among mammals and other vertebrates. We suggest that this orchestration contributes to tissue-specific and ontogenetic-specific expression, which can facilitate, for instance, timely protein complex formation.

Molecular characterizations of genes in chloroplast genomes of the genus Arachis L. (Fabaceae) based on the codon usage divergence

Studies on the molecular characteristics of chloroplast genome are generally important for clarifying the evolutionary processes of plant species. The base composition, the effective number of codons, the relative synonymous codon usage, the codon bias index, and their correlation coefficients of a total of 41 genes in 21 chloroplast genomes of the genus Arachis were investigated to further perform the correspondence and clustering analyses, revealing significantly higher variations in genomes of wild species than those of the cultivated taxa. The codon usage patterns of all 41 genes in the genus Arachis were AT-rich, suggesting that the natural selection was the main factor affecting the evolutionary history of these genomes. Five genes (i.e., ndhC, petD, atpF, rpl14, and rps11) and five genes (i.e., atpE, psbD, psaB, ycf2, and rps12) showed higher and lower base usage divergences, respectively. This study provided novel insights into our understanding of the molecular evolution of chloroplast genomes in the genus Arachis.

Comprehensive analysis of codon bias in 13 Ganoderma mitochondrial genomes

Introduction

Codon usage bias is a prevalent phenomenon observed across various species and genes. However, the specific attributes of codon usage in the mitochondrial genome of Ganoderma species remain unknown.

Methods

In this study, we investigated the codon bias of 12 mitochondrial core protein-coding genes (PCGs) in 9 Ganoderma species, including 13 Ganoderma strains.

Results

The codons of all Ganoderma strains showed a preference for ending in A/T. Additionally, correlations between codon base composition and the codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) were identified, demonstrating the impact of base composition on codon bias. Various base bias indicators were found to vary between or within Ganoderma strains, including GC3s, the CAI, the CBI, and the FOP. The results also revealed that the mitochondrial core PCGs of Ganoderma have an average effective number of codons (ENC) lower than 35, indicating strong bias toward certain codons. Evidence from neutrality plot and PR2-bias plot analysis indicates that natural selection is a major factor affecting codon bias in Ganoderma. Additionally, 11 to 22 optimal codons (ΔRSCU>0.08 and RSCU>1) were identified in 13 Ganoderma strains, with GCA, AUC, and UUC being the most widely used optimal codons in Ganoderma. By analyzing the combined mitochondrial sequences and relative synonymous codon usage (RSCU) values, the genetic relationships between or within Ganoderma strains were determined, indicating variations between them. Nevertheless, RSCU-based analysis illustrated the intra- and interspecies relationships of certain Ganoderma species.

Discussion

This study deepens our insight into the synonymous codon usage characteristics, genetics, and evolution of this important fungal group.

Selection for Translational Efficiency in Genes Associated with Alphaproteobacterial Gene Transfer Agents

Gene transfer agents (GTAs) are virus-like elements that are encoded by some bacterial and archaeal genomes. The production of GTAs can be induced by carbon depletion and results in host lysis and the release of virus-like particles that contain mostly random fragments of the host DNA. The remaining members of a GTA-producing population act as GTA recipients by producing proteins needed for GTA-mediated DNA acquisition. Here, we detected a codon usage bias toward codons with more readily available tRNAs in the RcGTA-like GTA genes of alphaproteobacterial genomes. Such bias likely improves the translational efficacy during GTA gene expression. While the strength of codon usage bias fluctuates substantially among individual GTA genes and across taxonomic groups, it is especially pronounced in Sphingomonadales, whose members are known to inhabit nutrient-depleted environments. By screening genomes for gene families with trends in codon usage biases similar to those in GTA genes, we found a gene that likely encodes head completion protein in some GTAs where it appeared missing, and 13 genes previously not implicated in the GTA life cycle. The latter genes are involved in various molecular processes, including the homologous recombination and transport of scarce organic matter. Our findings provide insights into the role of selection for translational efficiency in the evolution of GTA genes and outline genes that are potentially involved in the previously hypothesized integration of GTA-delivered DNA into the host genome. IMPORTANCE Horizontal gene transfer (HGT) is a fundamental process that drives evolution of microorganisms. HGT can result in a rapid dissemination of beneficial genes within and among microbial communities and can be achieved via multiple mechanisms. One peculiar HGT mechanism involves viruses "domesticated" by some bacteria and archaea (their hosts). These so-called gene transfer agents (GTAs) are encoded in hosts' genomes, produced under starvation conditions, and cannot propagate themselves as viruses. We show that GTA genes are under selection to improve the efficiency of their translation when the host activates GTA production. The selection is especially pronounced in bacteria that occupy nutrient-depleted environments. Intriguingly, several genes involved in incorporation of DNA into a genome are under similar selection pressure, suggesting that they may facilitate the integration of GTA-delivered DNA into the host genome. Our findings underscore the potential importance of GTAs as a mechanism of HGT under nutrient-limited conditions, which are widespread in microbial habitats.

Chloroplast gene expression level is negatively correlated with evolutionary rates and selective pressure while positively with codon usage bias in Ophioglossum vulgatum L

BACKGROUND

Characterization of the key factors determining gene expression level has been of significant interest. Previous studies on the relationship among evolutionary rates, codon usage bias, and expression level mostly focused on either nuclear genes or unicellular/multicellular organisms but few in chloroplast (cp) genes. Ophioglossum vulgatum is a unique fern and has important scientific and medicinal values. In this study, we sequenced its cp genome and transcriptome to estimate the evolutionary rates (dN and dS), selective pressure (dN/dS), gene expression level, codon usage bias, and their correlations.

RESULTS

The correlation coefficients between dN, dS, and dN/dS, and Transcripts Per Million (TPM) average values were -0.278 (P = 0.027 < 0.05), -0.331 (P = 0.008 < 0.05), and -0.311 (P = 0.013 < 0.05), respectively. The codon adaptation index (CAI) and tRNA adaptation index (tAI) were significantly positively correlated with TPM average values (P < 0.05).

CONCLUSIONS

Our results indicated that when the gene expression level was higher, the evolutionary rates and selective pressure were lower, but the codon usage bias was stronger. We provided evidence from cp gene data which supported the E-R (E stands for gene expression level and R stands for evolutionary rate) anti-correlation.

Evolution and codon usage bias of mitochondrial and nuclear genomes in Aspergillus section Flavi

The fungal genus Aspergillus contains a diversity of species divided into taxonomic sections of closely related species. Section Flavi contains 33 species, many of industrial, agricultural, or medical relevance. Here, we analyze the mitochondrial genomes (mitogenomes) of 20 Flavi species-including 18 newly assembled mitogenomes-and compare their evolutionary history and codon usage bias patterns to their nuclear counterparts. Codon usage bias refers to variable frequencies of synonymous codons in coding DNA and is shaped by a balance of neutral processes and natural selection. All mitogenomes were circular DNA molecules with highly conserved gene content and order. As expected, genomic content, including GC content, and genome size differed greatly between mitochondrial and nuclear genomes. Phylogenetic analysis based on 14 concatenated mitochondrial genes predicted evolutionary relationships largely consistent with those predicted by a phylogeny constructed from 2,422 nuclear genes. Comparing similarities in interspecies patterns of codon usage bias between mitochondrial and nuclear genomes showed that species grouped differently by patterns of codon usage bias depending on whether analyses were performed using mitochondrial or nuclear relative synonymous usage values. We found that patterns of codon usage bias at gene level are more similar between mitogenomes of different species than the mitogenome and nuclear genome of the same species. Finally, we inferred that, although most genes-both nuclear and mitochondrial-deviated from the neutral expectation for codon usage, mitogenomes were not under translational selection while nuclear genomes were under moderate translational selection. These results contribute to the study of mitochondrial genome evolution in filamentous fungi.

Composition and Codon Usage Pattern Results in Divergence of the Zinc Binuclear Cluster (Zn(II)2Cys6) Sequences among Ascomycetes Plant Pathogenic Fungi

Zinc binuclear cluster proteins (ZBC; Zn(II)2Cys6) are unique to the fungi kingdom and associated with a series of functions, viz., the utilization of macromolecules, stress tolerance, and most importantly, host-pathogen interactions by imparting virulence to the pathogen. Codon usage bias (CUB) is the phenomenon of using synonymous codons in a non-uniform fashion during the translation event, which has arisen because of interactions among evolutionary forces. The Zn(II)2Cys6 coding sequences from nine Ascomycetes plant pathogenic species and model system yeast were analysed for compositional and codon usage bias patterns. The clustering analysis diverged the Ascomycetes fungi into two clusters. The nucleotide compositional and relative synonymous codon usage (RSCU) analysis indicated GC biasness toward Ascomycetes fungi compared with the model system S. cerevisiae, which tends to be AT-rich. Further, plant pathogenic Ascomycetes fungi belonging to cluster-2 showed a higher number of GC-rich high-frequency codons than cluster-1 and was exclusively AT-rich in S. cerevisiae. The current investigation also showed the mutual effect of the two evolutionary forces, viz. natural selection and compositional constraints, on the CUB of Zn(II)2Cys6 genes. The perseverance of GC-rich codons of Zn(II)2Cys6 in Ascomycetes could facilitate the invasion process. The findings of the current investigation show the role of CUB and nucleotide composition in the evolutionary divergence of Ascomycetes plant pathogens and paves the way to target specific codons and sequences to modulate host-pathogen interactions through genome editing and functional genomics tools.

The complete mitochondrial genome of Talpa martinorum (Mammalia: Talpidae), a mole species endemic to Thrace: genome content and phylogenetic considerations

The complete mitogenome sequence of Talpa martinorum, a recently described Balkan endemic mole, was assembled from next generation sequence data. The mitogenome is similar to that of the three other Talpa species sequenced to date, being 16,835 bp in length, and containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, an origin of L-strand replication, and a control region or D-loop. Compared to other Talpa mitogenomes sequenced to date, that of T. martinorum differs in the length of D-loop and stop codon usage. TAG and T-- are the stop codons for the ND1 and ATP8 genes, respectively, in T. martinorum, whilst TAA acts as a stop codon for both ND1 and ATP8 in the other three Talpa species sequenced. Phylogeny reconstructions based on Maximum Likelihood and Bayesian inference analyses yielded phylogenies with similar topologies, demonstrating that T. martinorum nests within the western lineage of the genus, being closely related to T. aquitania and T. occidentalis.

Strong Selectional Forces Fine-Tune CpG Content in Genes Involved in Neurological Disorders as Revealed by Codon Usage Patterns

Neurodegenerative disorders cause irreversible damage to the neurons and adversely affect the quality of life. Protein misfolding and their aggregation in specific parts of the brain, mitochondrial dysfunction, calcium load, proteolytic stress, and oxidative stress are among the causes of neurodegenerative disorders. In addition, altered metabolism has been associated with neurodegeneration as evidenced by reductions in glutamine and alanine in transient global amnesia patients, higher homocysteine-cysteine disulfide, and lower methionine decline in serum urea have been observed in Alzheimer's disease patients. Neurodegeneration thus appears to be a culmination of altered metabolism. The study's objective is to analyze various attributes like composition, physical properties of the protein, and factors like selectional and mutational forces, influencing codon usage preferences in a panel of genes involved directly or indirectly in metabolism and contributing to neurodegeneration. Various parameters, including gene composition, dinucleotide analysis, Relative synonymous codon usage (RSCU), Codon adaptation index (CAI), neutrality and parity plots, and different protein indices, were computed and analyzed to determine the codon usage pattern and factors affecting it. The correlation of intrinsic protein properties such as the grand average of hydropathicity index (GRAVY), isoelectric point, hydrophobicity, and acidic, basic, and neutral amino acid content has been found to influence codon usage. In genes up to 800 amino acids long, the GC3 content was highly variable, while GC12 content was relatively constant. An optimum CpG content is present in genes to maintain a high expression level as required for genes involved in metabolism. Also observed was a low codon usage bias with a higher protein expression level. Compositional parameters and nucleotides at the second position of codons played essential roles in explaining the extent of bias. Overall analysis indicated that the dominance of selection pressure and compositional constraints and mutational forces shape codon usage.

Intragenomic variation in non-adaptive nucleotide biases causes underestimation of selection on synonymous codon usage

Patterns of non-uniform usage of synonymous codons vary across genes in an organism and between species across all domains of life. This codon usage bias (CUB) is due to a combination of non-adaptive (e.g. mutation biases) and adaptive (e.g. natural selection for translation efficiency/accuracy) evolutionary forces. Most models quantify the effects of mutation bias and selection on CUB assuming uniform mutational and other non-adaptive forces across the genome. However, non-adaptive nucleotide biases can vary within a genome due to processes such as biased gene conversion (BGC), potentially obfuscating signals of selection on codon usage. Moreover, genome-wide estimates of non-adaptive nucleotide biases are lacking for non-model organisms. We combine an unsupervised learning method with a population genetics model of synonymous coding sequence evolution to assess the impact of intragenomic variation in non-adaptive nucleotide bias on quantification of natural selection on synonymous codon usage across 49 Saccharomycotina yeasts. We find that in the absence of a priori information, unsupervised learning can be used to identify genes evolving under different non-adaptive nucleotide biases. We find that the impact of intragenomic variation in non-adaptive nucleotide bias varies widely, even among closely-related species. We show that the overall strength and direction of translational selection can be underestimated by failing to account for intragenomic variation in non-adaptive nucleotide biases. Interestingly, genes falling into clusters identified by machine learning are also physically clustered across chromosomes. Our results indicate the need for more nuanced models of sequence evolution that systematically incorporate the effects of variable non-adaptive nucleotide biases on codon frequencies.

Comprehensive analysis of complete chloroplast genome and phylogenetic aspects of ten Ficus species

BACKGROUND

The large genus Ficus comprises approximately 800 species, most of which possess high ornamental and ecological values. However, its evolutionary history remains largely unknown. Plastome (chloroplast genome) analysis had become an essential tool for species identification and for unveiling evolutionary relationships between species, genus and other rank groups. In this work we present the plastomes of ten Ficus species.

RESULTS

The complete chloroplast (CP) genomes of eleven Ficus specimens belonging to ten species were determined and analysed. The full length of the Ficus plastome was nearly 160 kbp with a similar overall GC content, ranging from 35.88 to 36.02%. A total of 114 unique genes, distributed in 80 protein-coding genes, 30 tRNAs, and 4 rRNAs, were annotated in each of the Ficus CP genome. In addition, these CP genomes showed variation in their inverted repeat regions (IR). Tandem repeats and mononucleotide simple sequence repeat (SSR) are widely distributed across the Ficus CP genome. Comparative genome analysis showed low sequence variability. In addition, eight variable regions to be used as potential molecular markers were proposed for future Ficus species identification. According to the phylogenetic analysis, these ten Ficus species were clustered together and further divided into three clades based on different subgenera. Simultaneously, it also showed the relatedness between Ficus and Morus.

CONCLUSION

The chloroplast genome structure of 10 Ficus species was similar to that of other angiosperms, with a typical four-part structure. Chloroplast genome sizes vary slightly due to expansion and contraction of the IR region. And the variation of noncoding regions of the chloroplast genome is larger than that of coding regions. Phylogenetic analysis showed that these eleven sampled CP genomes were divided into three clades, clustered with species from subgenus Urostigma, Sycomorus, and Ficus, respectively. These results support the Berg classification system, in which the subgenus Ficus was further decomposed into the subgenus Sycomorus. In general, the sequencing and analysis of Ficus plastomes, especially the ones of species with no or limited sequences available yet, contribute to the study of genetic diversity and species evolution of Ficus, while providing useful information for taxonomic and phylogenetic studies of Ficus.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

BioKIT: a versatile toolkit for processing and analyzing diverse types of sequence data

Bioinformatic analysis-such as genome assembly quality assessment, alignment summary statistics, relative synonymous codon usage, file format conversion, and processing and analysis-is integrated into diverse disciplines in the biological sciences. Several command-line pieces of software have been developed to conduct some of these individual analyses, but unified toolkits that conduct all these analyses are lacking. To address this gap, we introduce BioKIT, a versatile command line toolkit that has, upon publication, 42 functions, several of which were community-sourced, that conduct routine and novel processing and analysis of genome assemblies, multiple sequence alignments, coding sequences, sequencing data, and more. To demonstrate the utility of BioKIT, we conducted a comprehensive examination of relative synonymous codon usage across 171 fungal genomes that use alternative genetic codes, showed that the novel metric of gene-wise relative synonymous codon usage can accurately estimate gene-wise codon optimization, evaluated the quality and characteristics of 901 eukaryotic genome assemblies, and calculated alignment summary statistics for 10 phylogenomic data matrices. BioKIT will be helpful in facilitating and streamlining sequence analysis workflows. BioKIT is freely available under the MIT license from GitHub (https://github.com/JLSteenwyk/BioKIT), PyPi (https://pypi.org/project/jlsteenwyk-biokit/), and the Anaconda Cloud (https://anaconda.org/jlsteenwyk/jlsteenwyk-biokit). Documentation, user tutorials, and instructions for requesting new features are available online (https://jlsteenwyk.com/BioKIT).

An orthologous gene coevolution network provides insight into eukaryotic cellular and genomic structure and function

The evolutionary rates of functionally related genes often covary. We present a gene coevolution network inferred from examining nearly 3 million orthologous gene pairs from 332 budding yeast species spanning ~400 million years of evolution. Network modules provide insight into cellular and genomic structure and function. Examination of the phenotypic impact of network perturbation using deletion mutant data from the baker's yeast Saccharomyces cerevisiae, which were obtained from previously published studies, suggests that fitness in diverse environments is affected by orthologous gene neighborhood and connectivity. Mapping the network onto the chromosomes of S. cerevisiae and Candida albicans revealed that coevolving orthologous genes are not physically clustered in either species; rather, they are often located on different chromosomes or far apart on the same chromosome. The coevolution network captures the hierarchy of cellular structure and function, provides a roadmap for genotype-to-phenotype discovery, and portrays the genome as a linked ensemble of genes.

Genome-wide analysis of TALE superfamily in Triticum aestivum reveals TaKNOX11-A is involved in abiotic stress response

BACKGROUND

Three-amino-loop-extension (TALE) superfamily genes are widely present in plants and function directly in plant growth and development and abiotic stress response. Although TALE genes have been studied in many plant species, members of the TALE family have not been identified in wheat.

RESULTS

In this study, we identified 70 wheat TALE protein candidate genes divided into two subfamilies, KNOX (KNOTTED-like homeodomain) and BEL1-like (BLH/BELL homeodomain). Genes in the same subfamily or branch in the phylogenetic tree are similar in structure, and their encoded proteins have similar motifs and conserved structures. Wheat TALE genes are unevenly distributed on 21 chromosomes and expanded on the fourth chromosome. Through gene duplication analysis, 53 pairs of wheat TALE genes were determined to result from segmental duplication events, and five pairs were caused by tandem duplication events. The Ka/Ks between TALE gene pairs indicates a strong purification and selection effect. There are multiple cis-elements in the 2000 bp promoter sequence that respond to hormones and abiotic stress, indicating that most wheat TALE genes are involved in the growth, development, and stress response of wheat. We also studied the expression profiles of wheat TALE genes in different developmental stages and tissues and under different stress treatments. We detected the expression levels of four TALE genes by qRT-PCR, and selected TaKNOX11-A for further downstream analysis. TaKNOX11-A enhanced the drought and salt tolerances of Arabidopsis thaliana. TaKNOX11-A overexpressing plants had decreased malondialdehyde content and increased proline content, allowing for more effective adaptation of plants to unfavorable environments.

CONCLUSIONS

We identified TALE superfamily members in wheat and conducted a comprehensive bioinformatics analysis. The discovery of the potential role of TaKNOX11-A in drought resistance and salt tolerance provides a basis for follow-up studies of wheat TALE family members, and also provides new genetic resources for improving the stress resistance of wheat.

Kingdom-Wide Analysis of Fungal Transcriptomes and tRNAs Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

A detailed comparative analysis of codon usage bias in Alongshan virus

Alongshan virus (ALSV) is an emerging tick-borne pathogen that infects humans, causing febrile disease. ALSV uses Ixodes Persulcatus ticks to infect humans with a wide range of signs, from asymptomatic to encephalitis-like syndrome. There is an increasing public health concern about the ALSV infection. To get insight into the impacts of viral relations with their hosts on viral ability, survival, and evasion from hosts immune systems remain unknown. The codon usage is a driving force in viral genome evolution; therefore, we enrolled 41 ALSV strains in codon usage analysis to elucidate the molecular evolutionary dynamics of ALSV. The results indicate that the overall codon usage among ALSV isolates is relatively similar and slightly biased. Base compositions for the cds were in order of G >A >C >U and in the third position of codons G3 >A3 >C3 >T3. The RSCU values revealed that the more frequently used codons were mostly GC ended. Different codon preferences in ALSV genes in relation to codon usage of H. sapiens and Ixodes Persulcatus genes were found. Neutrality plot was determined to reveal the superiority of natural selection over directional mutation pressure in causing CUB based on GC12 versus GC3 contents. The results of these studies suggest that the emergence of ALSV in China, Russia and Finland may also be reflected in ALSV codon usage. Altogether, the presence of both mutation pressure and natural selection effect in shaping the codon usage patterns of ALSV.

Strategies and Patterns of Codon Bias in Molluscum Contagiosum Virus

Trends associated with codon usage in molluscum contagiosum virus (MCV) and factors governing the evolution of codon usage have not been investigated so far. In this study, attempts were made to decipher the codon usage trends and discover the major evolutionary forces that influence the patterns of codon usage in MCV with special reference to sub-types 1 and 2, MCV-1 and MCV-2, respectively. Three hypotheses were tested: (1) codon usage patterns of MCV-1 and MCV-2 are identical; (2) SCUB (synonymous codon usage bias) patterns of MCV-1 and MCV-2 slightly deviate from that of human host to avoid affecting the fitness of host; and (3) translational selection predominantly shapes the SCUB of MCV-1 and MCV-2. Various codon usage indices viz. relative codon usage value, effective number of codons and codon adaptation index were calculated to infer the nature of codon usage. Correspondence analysis and correlation analysis were performed to assess the relative contribution of silent base contents and significance of codon usage indices in defining bias in codon usage. Among the tested hypotheses, only the second and third hypotheses were accepted.

Codon Usage Bias in Autophagy-Related Gene 13 in Eukaryotes: Uncovering the Genetic Divergence by the Interplay Between Nucleotides and Codon Usages

Synonymous codon usage bias is a universal characteristic of genomes across various organisms. Autophagy-related gene 13 (atg13) is one essential gene for autophagy initiation, yet the evolutionary trends of the atg13 gene at the usages of nucleotide and synonymous codon remains unexplored. According to phylogenetic analyses for the atg13 gene of 226 eukaryotic organisms at the nucleotide and amino acid levels, it is clear that their nucleotide usages exhibit more genetic information than their amino acid usages. Specifically, the overall nucleotide usage bias quantified by information entropy reflected that the usage biases at the first and second codon positions were stronger than those at the third position of the atg13 genes. Furthermore, the bias level of nucleotide ‘G’ usage is highest, while that of nucleotide ‘C’ usage is lowest in the atg13 genes. On top of that, genetic features represented by synonymous codon usage exhibits a species-specific pattern on the evolution of the atg13 genes to some extent. Interestingly, the codon usages of atg13 genes in the ancestor animals (Latimeria chalumnae, Petromyzon marinus, and Rhinatrema bivittatum) are strongly influenced by mutation pressure from nucleotide composition constraint. However, the distributions of nucleotide composition at different codon positions in the atg13 gene display that natural selection still dominates atg13 codon usages during organisms’ evolution.

Analysis of Codon Usage Patterns of Six Sequenced Brachypodium distachyon Lines Reveals a Declining CG Skew of the CDSs from the 5'-ends to the 3'-ends

Brachypodium distachyon, a new monocotyledonous model plant, has received wide attention in biological research due to its small genome and numerous genetic resources. Codon usage bias is an important feature of genes and genomes, and it can be used in transgenic and evolutionary studies. In this study, the nucleotide compositions and patterns of codon usage bias were calculated using Codon W. Additionally, an ENC plot, Parity rule 2 and correspondence analyses were used to explore the major factors influencing codon usage bias patterns. The numbers of hydrogen bonds and skews were used to analyze the GC trend in the 5'-ends of the coding sequences. The results showed that minor differences in the codon usage bias patterns were revealed by the ENC plot, Parity rule 2 and correspondence analyses. The analyses of the CG-skew and the number of hydrogen bonds showed a declining trend in the number of cytosines at the 5'-ends of the CDSs (from the 5'-ends to the 3'-ends), indicating that GC may play a major role in codon usage bias. In addition, our results laid a foundation for the study of codon usage bias patterns in Brachypodium genus and suggested that the GC plays a major role in determining these patterns.

Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.)

Cystatins, as reversible inhibitors of papain-like and legumain proteases, have been identified in several plant species. Although the cystatin family plays crucial roles in plant development and defense responses to various stresses, this family in wheat (Triticum aestivum L.) is still poorly understood. In this study, 55 wheat cystatins (TaCystatins) were identified. All TaCystatins were divided into three groups and both the conserved gene structures and peptide motifs were relatively conserved within each group. Homoeolog analysis suggested that both homoeolog retention percentage and gene duplications contributed to the abundance of the TaCystatin family. Analysis of duplication events confirmed that segmental duplications played an important role in the duplication patterns. The results of codon usage pattern analysis showed that TaCystatins had evident codon usage bias, which was mainly affected by mutation pressure. TaCystatins may be regulated by cis-acting elements, especially abscisic acid and methyl jasmonate responsive elements. In addition, the expression of all selected TaCystatins was significantly changed following viral infection and cold stress, suggesting potential roles in response to biotic and abiotic challenges. Overall, our work provides new insights into TaCystatins during wheat evolution and will help further research to decipher the roles of TaCystatins under diverse stress conditions.

Comparative Analysis of Genomic and Transcriptome Sequences Reveals Divergent Patterns of Codon Bias in Wheat and Its Ancestor Species

The synonymous codons usage shows a characteristic pattern of preference in each organism. This codon usage bias is thought to have evolved for efficient protein synthesis. Synonymous codon usage was studied in genes of the hexaploid wheat Triticum aestivum (AABBDD) and its progenitor species, Triticum urartu (AA), Aegilops tauschii (DD), and Triticum turgidum (AABB). Triticum aestivum exhibited stronger usage bias for G/C-ending codons than did the three progenitor species, and this bias was especially higher compared to T. turgidum and Ae. tauschii. High GC content is a primary factor influencing codon usage in T. aestivum. Neutrality analysis showed a significant positive correlation (p<0.001) between GC12 and GC3 in the four species with regression line slopes near zero (0.16–0.20), suggesting that the effect of mutation on codon usage was only 16–20%. The GC3s values of genes were associated with gene length and distribution density within chromosomes. tRNA abundance data indicated that codon preference corresponded to the relative abundance of isoaccepting tRNAs in the four species. Both mutation and selection have affected synonymous codon usage in hexaploid wheat and its progenitor species. GO enrichment showed that GC biased genes were commonly enriched in physiological processes such as photosynthesis and response to acid chemical. In some certain gene families with important functions, the codon usage of small parts of genes has changed during the evolution process of T. aestivum.