Comprehensive Analysis of Codon Usage Bias in Seven Epichloë Species and Their Peramine-Coding Genes

Codon usage bias in chloroplast genes implicate adaptive evolution of four ginger species

Codon usage bias (CUB) refers to different codons exhibiting varying frequencies of usage in the genome. Studying CUB is crucial for understanding genome structure, function, and evolutionary processes. Herein, we investigated the codon usage patterns and influencing factors of protein-coding genes in the chloroplast genomes of four sister genera (monophyletic Roscoea and Cautleya, and monophyletic Pommereschea and Rhynchanthus) from the Zingiberaceae family with contrasting habitats in southwestern China. These genera exhibit distinct habitats, providing a unique opportunity to explore the adaptive evolution of codon usage. We conducted a comprehensive analysis of nucleotide composition and codon usage on protein-coding genes in the chloroplast genomes. The study focused on understanding the relationship between codon usage and environmental adaptation, with a particular emphasis on genes associated with photosynthesis. Nucleotide composition analysis revealed that the overall G/C content of the coding genes was ˂ 48%, indicating an enrichment of A/T bases. Additionally, synonymous and optimal codons were biased toward ending with A/U bases. Natural selection is the primary factor influencing CUB characteristics, particularly photosynthesis-associated genes. We observed differential gene expressions related to light adaptation among sister genera inhabiting different environments. Certain codons were favored under specific conditions, possibly contributing to gene expression regulation in particular environments. This study provides insights into the adaptive evolution of these sister genera by analyzing CUB and offers theoretical assistance for understanding gene expression and regulation. In addition, the data support the relationship between RNA editing and CUB, and the findings shed light on potential research directions for investigating adaptive evolution.

Codon Usage Provide Insights into the Adaptation of Rice Genes under Stress Condition

Plants experience different stresses, i.e., abiotic, or biotic, and to combat them, plants re-program the expression of growth-, metabolism-, and resistance-related genes. These genes differ in their synonymous codon usage frequency and show codon usage bias. Here, we investigated the correlation among codon usage bias, gene expression, and underlying mechanisms in rice under abiotic and biotic stress conditions. The results indicated that genes with higher expression (up- or downregulated) levels had high GC content (≥60%), a low effective number of codon usage (≤40), and exhibited strong biases towards the codons with C/G at the third nucleotide position, irrespective of stress received. TTC, ATC, and CTC were the most preferred codons, while TAC, CAC, AAC, GAC, and TGC were moderately preferred under any stress (abiotic or biotic) condition. Additionally, downregulated genes are under mutational pressure (R² ≥ 0.5) while upregulated genes are under natural selection pressure (R² ≤ 0.5). Based on these results, we also identified the possible target codons that can be used to design an optimized set of genes with specific codons to develop climate-resilient varieties. Conclusively, under stress, rice has a bias towards codon usage which is correlated with GC content, gene expression level, and gene length.

An investigation of codon usage pattern analysis in pancreatitis associated genes

BACKGROUND

Pancreatitis is an inflammatory disorder resulting from the autoactivation of trypsinogen in the pancreas. The genetic basis of the disease is an old phenomenon, and evidence is accumulating for the involvement of synonymous/non-synonymous codon variants in disease initiation and progression.

RESULTS

The present study envisaged a panel of 26 genes involved in pancreatitis for their codon choices, compositional analysis, relative dinucleotide frequency, nucleotide disproportion, protein physical properties, gene expression, codon bias, and interrelated of all these factors. In this set of genes, gene length was positively correlated with nucleotide skews and codon usage bias. Codon usage of any gene is dependent upon its AT and GC component; however, AGG, CGT, and CGA encoding for Arg, TCG for Ser, GTC for Val, and CCA for Pro were independent of nucleotide compositions. In addition, Codon GTC showed a correlation with protein properties, isoelectric point, instability index, and frequency of basic amino acids. We also investigated the effect of various evolutionary forces in shaping the codon usage choices of genes.

CONCLUSIONS

This study will enable us to gain insight into the molecular signatures associated with the disease that might help identify more potential genes contributing to enhanced risk for pancreatitis. All the genes associated with pancreatitis are generally associated with physiological function, and mutations causing loss of function, over or under expression leads to an ailment. Therefore, the present study attempts to envisage the molecular signature in a group of genes that lead to pancreatitis in case of malfunction.

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)₂Cys₆) 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)₂Cys₆ 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)₂Cys₆ genes. The perseverance of GC-rich codons of Zn(II)₂Cys₆ 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.

Codon Usage is Influenced by Compositional Constraints in Genes Associated with Dementia

Dementia is a clinical syndrome characterized by progressive cognitive decline, and the symptoms could be gradual, persistent, and progressive. In the present study, we investigated 47 genes that have been linked to dementia. Compositional, selectional, and mutational forces were seen to be involved. Nucleotide components that influenced A- and GC-affected codon usages bias at all three codon positions. The influence of these two compositional constraints on codon usage bias (CUB) was positive for nucleotide A and negative for GC. Nucleotide A also experienced the highest mutational force, and GC-ending codons were preferred over AT-ending codons. A high bias toward GC-ending codons enhances the gene expression level, evidenced by the positive association between CAI- and GC-ending codons. Unusual behavior of the TTG codon showing an inverse relationship with the GC-ending codon and negative influence of gene expression, behavior contrary to all other GC-ending codons, shows an operative selectional force. Furthermore, parity analysis, higher translational selection value, preference of GC-ending codons over AT-ending codons, and association of gene length with gene expression refer to the dominant role of selection pressure with compositional constraint and mutational force-shaping codon usage.

Codon Usage Bias Correlates With Gene Length in Neurodegeneration Associated Genes

Codon usage analysis is a crucial part of molecular characterization and is used to determine the factors affecting the evolution of a gene. The length of a gene is an important parameter that affects the characteristics of the gene, such as codon usage, compositional parameters, and sometimes, its functions. In the present study, we investigated the association of various parameters related to codon usage with the length of genes. Gene expression is affected by nucleotide disproportion. In sixty genes related to neurodegenerative disorders, the G nucleotide was the most abundant and the T nucleotide was the least. The nucleotide T exhibited a significant association with the length of the gene at both the overall compositional level and the first and second codon positions. Codon usage bias (CUB) of these genes was affected by pyrimidine and keto skews. Gene length was found to be significantly correlated with codon bias in neurodegeneration associated genes. In gene segments with lengths below 1,200 bp and above 2,400 bp, CUB was positively associated with length. Relative synonymous CUB, which is another measure of CUB, showed that codons TTA, GTT, GTC, TCA, GGT, and GGA exhibited a positive association with length, whereas codons GTA, AGC, CGT, CGA, and GGG showed a negative association. GC-ending codons were preferred over AT-ending codons. Overall analysis indicated that the association between CUB and length varies depending on the segment size; however, CUB of 1,200–2,000 bp gene segments appeared not affected by gene length. In synopsis, analysis suggests that length of the genes correlates with various imperative molecular signatures including A/T nucleotide disproportion and codon choices. In the present study we additionally evaluated various molecular features and their correlation with different indices of codon usage, like the Codon Adaptation Index (CAI) and Relative Dynonymous Codon Usage (RSCU) of codons. We also considered the impact of gene fragment size on different molecular features in genes related to neurodegeneration. This analysis will aid our understanding of and in potentially modulating gene expression in cases of defective gene functioning in clinical settings.

Calculating and comparing codon usage values in rare disease genes highlights codon clustering with disease-and tissue- specific hierarchy

We designed a novel strategy to define codon usage bias (CUB) in 6 specific small cohorts of human genes. We calculated codon usage (CU) values in 29 non-disease-causing (NDC) and 31 disease-causing (DC) human genes which are highly expressed in 3 distinct tissues, kidney, muscle, and skin. We applied our strategy to the same selected genes annotated in 15 mammalian species. We obtained CUB hierarchical clusters for each gene cohort which showed tissue-specific and disease-specific CUB fingerprints. We showed that DC genes (especially those expressed in muscle) display a low CUB, well recognizable in codon hierarchical clustering. We defined the extremely biased codons as "zero codons" and found that their number is significantly higher in all DC genes, all tissues, and that this trend is conserved across mammals. Based on this calculation in different gene cohorts, we identified 5 codons which are more differentially used across genes and mammals, underlining that some genes have favorite synonymous codons in use. Since of the muscle genes clear clusters, and, among these, dystrophin gene surprisingly does not show any "zero codon" we adopted a novel approach to study CUB, we called "mapping-on-codons". We positioned 2828 dystrophin missense and nonsense pathogenic variations on their respective codon, highlighting that its frequency and occurrence is not dependent on the CU values. We conclude our strategy consents to identify a hierarchical clustering of CU values in a gene cohort-specific fingerprints, with recognizable trend across mammals. In DC muscle genes also a disease-related fingerprint can be observed, allowing discrimination between DC and NDC genes. We propose that using our strategy which studies CU in specific gene cohorts, as rare disease genes, and tissue specific genes, may provide novel information about the CUB role in human and medical genetics, with implications on synonymous variations interpretation and codon optimization algorithms.

Codon usage bias predicts the functional MYB10 gene in Populus

Analysis of codon usage bias (CUB) in different species can reveal the patterns of genetic information transfer across those species. To better understand the characteristics of MYB10-a key regulator of anthocyanin biosynthesis-and identify the true (functional) MYB10 gene among the two candidates in Populus, we analysed the coding sequences of MYB10 genes in 10 different species using Codon W, CHIPS, CUSP, and CAI. Majority of the optimal amino acid codons of MYB10 genes ended with A/U, and GGA, UCA, GCA, AGA, and CCA were over-represented in all plant species studied. Among the two most promising MYB10 gene candidates in Populus, Potri.17G125700 shared a higher similarity of codon usage with MYB10 genes from other plant species, suggesting that it encodes the functional MYB10 in Populus. We verified this speculation by cloning both candidate MYB10 genes from Populus into vectors to produce transiently transformed seedlings. Colour phenotypes and anthocyanin content of the transiently transformed seedlings indicated that Potri.17G125700 encodes the true MYB10 transcription factor, which positively regulates anthocyanin accumulation in Populus. Furthermore, CUB analysis was used to select the most promising MYB12 candidate in Malus sp. (crabapple). Our results demonstrate the effectiveness of CUB analysis as a promising method to identify the functional gene from a set of candidates in long-living plants with complex genetics.

Comparative study of codon usage profiles of Zingiber officinale and its associated fungal pathogens

Codon usage bias influences the genetic features prevalent in genomes of all the organisms. It also plays a crucial role in establishing the host-pathogen relationship. The present study elucidates the role of codon usage pattern regarding the predilection of fungal pathogens Aspergillus flavus, Aspergillus niger, Fusarium oxysporum and Colletotrichum gloeosporioides towards host plant Zingiber officinale. We found a similar trend of codon usage pattern operative in plant and fungal pathogens. This concurrence might be attributed for the colonization of fungal pathogens in Z. officinale. The transcriptome of both plant and pathogens showed bias towards GC-ending codons. Natural selection and mutational pressure seem to be accountable for shaping the codon usage pattern of host and pathogen. We also identified some distinctive preferred codons in A. flavus, F. oxysporum and Z. officinale that could be regarded as signature codons for the identification of these organisms. Knowledge of favored, avoided and unique codons will help to devise strategies for reducing spice losses due to fungal pathogens.

Evolutionary and codon usage preference insights into spike glycoprotein of SARS-CoV-2

Interaction of SARS-CoV-2 spike glycoprotein with the ACE2 cell receptor is very crucial for virus attachment to human cells. Selected mutations in SARS-CoV-2 S-protein are reported to strengthen its binding affinity to mammalian ACE2. The N501T mutation in SARS-CoV-2-CTD furnishes better support to hotspot 353 in comparison with SARS-CoV and shows higher affinity for receptor binding. Recombination analysis exhibited higher recombination events in SARS-CoV-2 strains, irrespective of their geographical origin or hosts. Investigation further supports a common origin among SARS-CoV-2 and its predecessors, SARS-CoV and bat-SARS-like-CoV. The recombination events suggest a constant exchange of genetic material among the co-infecting viruses in possible reservoirs and human hosts before SARS-CoV-2 emerged. Furthermore, a comprehensive analysis of codon usage bias (CUB) in SARS-CoV-2 revealed significant CUB among the S-genes of different beta-coronaviruses governed majorly by natural selection and mutation pressure. Various indices of codon usage of S-genes helped in quantifying its adaptability in other animal hosts. These findings might help in identifying potential experimental animal models for investigating pathogenicity for drugs and vaccine development experiments.

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.

Codon usage signatures in the genus Cryptococcus: A complex interplay of gene expression, translational selection and compositional bias

The fungal genus Cryptococcus comprises of several diverse species. The pathogens forming Cryptococcus neoformans/ Cryptococcus gatti species complex are of immense clinical significance owing to the high frequency of infections and deaths globally. Three closely related non-pathogenic species namely, Cryptococcus amylolentus, Cryptococcus wingfieldii and Cryptococcus depauperatus are the non-pathogenic ancestral species from which pathogenic lineages have diverged. In the current study, a comprehensive analysis of factors influencing the codon and amino acid usage bias in six pathogenic and three non-pathogenic species was performed. Our results revealed that though compositional bias played a crucial role, translational selection and gene expression were the key determinants of codon usage variations. Analysis of relative dinucleotide abundance and codon context signatures revealed strict avoidance of TpA dinucleotide across genomes. Multivariate statistical analysis based on codon usage data resulted in discrete clustering of pathogens and non-pathogens which correlated with previous reports on their phylogenetic distribution.

The Complete Chloroplast Genomes of Two Lespedeza Species: Insights into Codon Usage Bias, RNA Editing Sites, and Phylogenetic Relationships in Desmodieae (Fabaceae: Papilionoideae)

The genus Lespedeza (tribe: Desmodieae) consists of about 40 species that have high medicinal and economic value. However, in this genus, using morphological characters, the species identification is quite complicated, which can be solved by the analysis of the complete chloroplast genomes. As primary organelle genomes, the complete genome sequences of chloroplasts (cp) provide unique molecular information to study the divergence of species, RNA editing, and phylogeny. Therefore, to the best of our knowledge, for the first time, we sequenced the complete cp genomes of two representative Lespedeza species: Lespedeza davurica and Lespedeza cuneata. The cp genomes of both the species were found to be 149,010 bp in length, exhibiting the typical angiosperm chloroplast structure containing four regions. The Lespedeza cp genomes showed similar conserved gene contents, order, and orientations with a total GC content of 35.0%. A total of 128 genes, including 83 protein-coding genes, 37 tRNAs, and eight rRNAs, were identified from each genome. Unique molecular features of the two Lespedeza cp genome sequences were obtained by performing the analysis of repeats, sequence divergence, codon usage, and predicting the RNA editing sites in addition to phylogenetic analysis with other key genera in tribe Desmodieae. Using the two datasets, the phylogenetic relationship of Lespedeza species among Deasmodieae was discovered, suggesting that whole cp genomes provided useful information for phylogenetic studies of these species.

Comparative Genomic Analysis of Soil Dwelling Bacteria Utilizing a Combinational Codon Usage and Molecular Phylogenetic Approach Accentuating on Key Housekeeping Genes

Soil is a diversified and complex ecological niche, home to a myriad of microorganisms particularly bacteria. The physico-chemical complexities of soil results in a plethora of physiological variations to exist within the different types of soil dwelling bacteria, giving rise to a wide variation in genome structure and complexity. This serves as an attractive proposition to analyze and compare the genome of a large number soil bacteria to comprehend their genome complexity and evolution. In this study a combination of codon usage and molecular phylogenetics of the whole genome and key housekeeping genes like infB (translation initiation factor 2), trpB (tryptophan synthase, beta subunit), atpD (ATP synthase, beta subunit), and rpoB (RNA polymerase, beta subunit) of 92 soil bacterial species spread across the entire eubacterial domain and residing in different soil types was performed. The results indicated the direct relationship of genome size with codon bias and coding frequency in the studied bacteria. The codon usage profile demonstrated by the gene trpB was found to be relatively different from the rest of the housekeeping genes with a large number of bacteria having a greater percentage of genes with Nc values less than the Nc of trpB. The results from the overall codon usage bias profile also depicted that the codon usage bias in the key housekeeping genes of soil bacteria was majorly due to selectional pressure and not mutation. The analysis of hydrophobicity of the gene product encoded by the rpoB coding sequences demonstrated tight clustering across all the soil bacteria suggesting conservation of protein structure for maintenance of form and function. The phylogenetic affinities inferred using 16S rRNA gene and the housekeeping genes demonstrated conflicting signals with trpB gene being the noisiest one. The housekeeping gene atpD was found to depict the least amount of evolutionary change in the soil bacteria considered in this study except in two Clostridium species. The phylogenetic and codon usage analysis of the soil bacteria consistently demonstrated the relatedness of Azotobacter chroococcum with different species of the genus Pseudomonas.

Comparative in silico analysis of ftsZ gene from different bacteria reveals the preference for core set of codons in coding sequence structuring and secondary structural elements determination

The deluge of sequence information in the recent times provide us with an excellent opportunity to compare organisms on a large genomic scale. In this study we have tried to decipher the variation in the gene organization and structuring of a vital bacterial gene called ftsZ which codes for an integral component of the bacterial cell division, the FtsZ protein. FtsZ is homologous to tubulin protein and has been found to be ubiquitous in eubacteria. FtsZ is showing increasing promise as a target for antibacterial drug discovery. Our study of ftsZ protein from 143 different bacterial species spanning a wider range of morphological and physiological type demonstrates that the ftsZ gene of about ninety three percent of the organisms show relatively biased codon usage profile and significant GC deviation from their genomic GC content. Comparative codon usage analysis of ftsZ and a core housekeeping gene rpoB demonstrated that codon usage pattern of ftsZ CDS is shaped by natural selection to a large extent and mimics that of a housekeeping gene. We have also detected a tendency among the different organisms to utilize a core set of codons in structuring the ftsZ coding sequence. We observed that the compositional frequency of the amino acid serine in the FtsZ protein appears to be a indicator of the bacterial lifestyle. Our meticulous analysis of the ftsZ gene linked with the corresponding FtsZ protein show that there is a bias towards the use of specific synonymous codons particularly in the helix and strand regions of the multi-domain FtsZ protein. Overall our findings suggest that in an indispensable and vital protein such as FtsZ, there is an inherent tendency to maintain form for optimized performance in spite of the extrinsic variability in coding features.