The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens

Decoding Codon Usage Patterns in High-Risk Human Papillomavirus Genomes: A Comprehensive Analysis

Human Papillomavirus (HPV) is a major contributor to various human cancers, particularly cervical cancer. Despite its significant impact, the codon usage bias in high-risk HPV types has not been extensively studied. Understanding this bias, however, could provide valuable insights into the virus itself and inform the optimization of vaccine design. This study explores codon usage bias within the genomes of 17 high-risk HPV types (HPV-16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, and 82) through comparative analysis. While overall codon usage preference across these genotypes is not highly significant, a notable trend emerges in the preference for codons ending in A or U, with 24 out of 26 favored codons (Relative Synonymous Codon Usage > 1) ending in A or U. Moreover, no common optimal codons are shared among the 17 genomes. The study also identifies the underrepresentation of CpG and ApA dinucleotides, alongside the overrepresentation of CpA and UpG, which likely contribute to codon usage preferences that may influence viral replication and immune evasion strategies. Integrated analysis further suggests that natural selection is the primary force driving codon usage bias in these high-risk HPV genomes. Additionally, these HPVs exhibit a limited set of favored codons shared with humans, potentially minimizing competition for translation resources. This study offers new insights into codon usage bias in high-risk HPVs and underscores the importance of this understanding for optimizing vaccine design.

Phylogenomics and plastome evolution of Lithospermeae (Boraginaceae)

BACKGROUND

Lithospermeae is the largest tribe within Boraginaceae. The tribe has been the focus of multiple phylogenetic studies over the last 15 years, with most focused on one genus or a few genera. In the present study, we newly sequenced 69 species of Lithospermeae and relatives to analyze the phylogenomic relationships among its members as well as the evolution of the plastid genome.

RESULTS

The phylogeny of Lithospermeae resolved from the plastid genome and nrDNA cistron is generally congruent with prior studies, but is better resolved and supported. Increasing character sampling across the plastid genome results in gradually more similar trees to that from the entire plastid genome. Overall, plastid genome structure was quite consistent across Lithospermeae. Codon Usage Bias (CUB) analyses demonstrate that across Lithospermeae plastid genomes were rich in AT and poor in GC. Mutation may play a greater role than selection across the plastid genome of Lithospermeae. The present study is the first to highlight the CUB characteristics of Lithospermeae species, which can help elucidate the mechanisms underlying patterns of molecular evolution and improve the expression levels of exogenous genes by codon optimization.

CONCLUSIONS

This study provides a comprehensive phylogenomic analysis of Lithospermeae, significantly enhancing our understanding of the phylogenetic relationships and plastid genome evolution within this largest tribe of Boraginaceae. By utilizing an expanded genomic sampling approach, we have achieved increased resolution and support among the evolutionary relationships of the tribe, in line with but improving upon previous studies. The analyses of plastid genome structure revealed consistency across Lithospermeae, with a notable CUB. This study marks the first investigation into the CUB of Lithospermeae species and sets the stage for further research on the molecular evolution of plastid genomes across Boraginaceae.

Codon Bias of the DDR1 Gene and Transcription Factor EHF in Multiple Species

Milk production is an essential economic trait in cattle, and understanding the genetic regulation of this trait can enhance breeding strategies. The discoidin domain receptor 1 (DDR1) gene has been identified as a key candidate gene that influences milk production, and ETS homologous factor (EHF) is recognized as a critical transcription factor that regulates DDR1 expression. Codon usage bias, which affects gene expression and protein function, has not been fully explored in cattle. This study aims to examine the codon usage bias of DDR1 and EHF transcription factors to understand their roles in dairy production traits. Data from 24 species revealed that both DDR1 and EHF predominantly used G/C-ending codons, with the GC3 content averaging 75.49% for DDR1 and 61.72% for EHF. Synonymous codon usage analysis identified high-frequency codons for both DDR1 and EHF, with 17 codons common to both genes. Correlation analysis indicated a negative relationship between the effective number of codons and codon adaptation index for both DDR1 and EHF. Phylogenetic and clustering analyses revealed similar codon usage patterns among closely related species. These findings suggest that EHF plays a crucial role in regulating DDR1 expression, offering new insights into genetically regulating milk production in cattle.

Identification, molecular evolution, codon bias, and expansion analysis of NLP transcription factor family in foxtail millet (Setaria italica L.) and closely related crops

The NODULE-INCEPTION-like protein (NLP) family is a plant-specific transcription factor (TF) family involved in nitrate transport and assimilation in plants, which are essential for improving plant nitrogen use efficiency. Currently, the molecular nature and evolutionary trajectory of NLP genes in the C4 model crop foxtail millet are unknown. Therefore, we performed a comprehensive analysis of NLP and molecular evolution in foxtail millet by scanning the genomes of foxtail millet and representative species of the plant kingdom. We identified seven NLP genes in the foxtail millet genome, all of which are individually and separately distributed on different chromosomes. They were not structurally identical to each other and were mainly expressed on root tissues. We unearthed two key genes (Si5G004100.1 and Si6G248300.1) with a variety of excellent characteristics. Regarding its molecular evolution, we found that NLP genes in Gramineae mainly underwent dispersed duplication, but maize NLP genes were mainly generated via WGD events. Other factors such as base mutations and natural selection have combined to promote the evolution of NLP genes. Intriguingly, the family in plants showed a gradual expansion during evolution with more duplications than losses, contrary to most gene families. In conclusion, this study advances the use of NLP genetic resources and the understanding of molecular evolution in cereals.

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.

Combination of hydrophobicity and codon usage bias determines sorting of model K+ channel protein to either mitochondria or endoplasmic reticulum

When the K+ channel-like protein Kesv from Ectocarpus siliculosus virus 1 is heterologously expressed in mammalian cells, it is sorted to the mitochondria. This targeting can be redirected to the endoplasmic reticulum (ER) by altering the codon usage in distinct regions of the gene or by inserting a triplet of hydrophobic amino acids (AAs) into the protein's C-terminal transmembrane domain (ct-TMD). Systematic variations in the flavor of the inserted AAs and/or its codon usage show that a positive charge in the inserted AA triplet alone serves as strong signal for mitochondria sorting. In cases of neutral AA triplets, mitochondria sorting are favored by a combination of hydrophilic AAs and rarely used codons; sorting to the ER exhibits the inverse dependency. This propensity for ER sorting is particularly high when a common codon follows a rarer one in the AA triplet; mitochondria sorting in contrast is supported by codon uniformity. Since parameters like positive charge, hydrophobic AAs, and common codons are known to facilitate elongation of nascent proteins in the ribosome the data suggest a mechanism in which local changes in elongation velocity and co-translational folding in the ct-TMD influence intracellular protein sorting.

Intragenic codon usage in proteobacteria: Translational selection, IS expansion and genomic shrinkage

This manuscript presents a method to systematically study intragenic variations in codon usage using correspondence analysis and the effective number of codons. The method is applied to >1100 proteobacteria. Codon usage biases (measured as inertia) increases with genome size, the same is true for the percentage of inertia explained by the first axis. It is shown that there is often a relaxed or more uniform codon usage near the gene termini. Ithis is not seen n small genomes, notably those of intracellular organisms like Buchnera aphidicola or Rickettsia prowazekii where translational selection plays less of a role. When genes from E. coli, for which translational selection is well described, are split into low, intermediate and high expression, respectively, it is shown that the intragenic codon usage pattern with more uniform usage at termini exist across all three expression groups. Furthermore, the correspondence analysis reveals a unique pattern in Bordetella pertussis due to IS expansion. This study thus shows that translational selection, genome shrinkage and IS expansion result in characteristic patterns in intragenic codon usage.

Plastome Characterization and Phylogenomic Analysis Yield New Insights into the Evolutionary Relationships among the Species of the Subgenus Bryocles (Hosta; Asparagaceae) in East Asia

The genus Hosta, which has a native distribution in temperate East Asia and a number of species ranging from 23 to 40, represents a taxonomically important and ornamentally popular plant. Despite its taxonomic and horticultural importance, the genus Hosta has remained taxonomically challenging owing to insufficient diagnostic features, continuous morphological variation, and the process of hybridization and introgression, making species circumscription and phylogenetic inference difficult. In this study, we sequenced 11 accessions of Hosta plastomes, including members of three geographically defined subgenera, Hosta, Bryocles, and Giboshi, determined the characteristics of plastomes, and inferred their phylogenetic relationships. We found highly conserved plastomes among the three subgenera, identified several mutation hotspots that can be used as barcodes, and revealed the patterns of codon usage bias and RNA editing sites. Five positively selected plastome genes (rbcL, rpoB, rpoC2, rpl16, and rpl20) were identified. Phylogenetic analysis suggested (1) the earliest divergence of subg. Hosta, (2) non-monophyly of subg. Bryocles and its two sections (Lamellatae and Stoloniferae), (3) a sister relationship between H. sieboldiana (subg. Giboshi) and H. ventricosa (subg. Bryocles), and (4) reciprocally monophyletic and divergent lineages of H. capitata in Korea and Japan, requiring further studies of their taxonomic distinction.

Analysis of codon usage of Horseshoe Bat Hepatitis B virus and its host

In the present analysis, codon usage strategies and base distribution of Horseshoe bat hepatitis B virus (HBHBV) were analyzed and compared with its host Rhinolophus sinicus, as no work was yet reported. The magnitude of synonymous codon usage bias (CUB) in the virus and its host was low with higher proportion of the base C. Notably, 21 more frequently used codons, 19 less frequently used codons and 3 underrepresented codons (TCG, ACG and GCG) were found to be similar in both virus and its host coding sequences. Neutrality plot analysis reported greater role of natural selection in HBHBV (67.84%) and R. sinicus (76.90%) over mutation pressure. Base skewness and protein properties also influenced the CUB of genes. Further, codon usage analysis depicted, HBHBV and R. sinicus had many similarities in codon usage patterns that might reflect viral adaptation to its host.

Synonymous but Not Silent: The Codon Usage Code for Gene Expression and Protein Folding

Codon usage bias, the preference for certain synonymous codons, is found in all genomes. Although synonymous mutations were previously thought to be silent, a large body of evidence has demonstrated that codon usage can play major roles in determining gene expression levels and protein structures. Codon usage influences translation elongation speed and regulates translation efficiency and accuracy. Adaptation of codon usage to tRNA expression determines the proteome landscape. In addition, codon usage biases result in nonuniform ribosome decoding rates on mRNAs, which in turn influence the cotranslational protein folding process that is critical for protein function in diverse biological processes. Conserved genome-wide correlations have also been found between codon usage and protein structures. Furthermore, codon usage is a major determinant of mRNA levels through translation-dependent effects on mRNA decay and translation-independent effects on transcriptional and posttranscriptional processes. Here, we discuss the multifaceted roles and mechanisms of codon usage in different gene regulatory processes.

Characterization and comparative analysis among plastome sequences of eight endemic Rubus (Rosaceae) species in Taiwan

Genus Rubus represents the second largest genus of the family Rosaceae in Taiwan, with 41 currently recognized species across three subgenera (Chamaebatus, Idaoeobatus, and Malochobatus). Despite previous morphological and cytological studies, little is known regarding the overall phylogenetic relationships among the Rubus species in Taiwan, and their relationships to congeneric species in continental China. We characterized eight complete plastomes of Taiwan endemic Rubus species: subg. Idaeobatus (R. glandulosopunctatus, R. incanus, R. parviaraliifolius, R rubroangustifolius, R. taitoensis, and R. taiwanicolus) and subg. Malachobatus (R. kawakamii and R. laciniastostipulatus) to determine their phylogenetic relationships. The plastomes were highly conserved and the size of the complete plastome sequences ranged from 155,566 to 156,236 bp. The overall GC content ranged from 37.0 to 37.3%. The frequency of codon usage showed similar patterns among species, and 29 of the 73 common protein-coding genes were positively selected. The comparative phylogenomic analysis identified four highly variable intergenic regions (rps16/trnQ, petA/psbJ, rpl32/trnL-UAG, and trnT-UGU/trnL-UAA). Phylogenetic analysis of 31 representative complete plastomes within the family Rosaceae revealed three major lineages within Rubus in Taiwan. However, overall phylogenetic relationships among endemic species require broader taxon sampling to gain new insights into infrageneric relationships and their plastome evolution.

Analysis of codon usage bias and evolution in the chloroplast genome of Mesona chinensis Benth

Mesona chinensis Benth (MCB) is one of the main economic crops in tropical and subtropical areas. To understand the codon usage bias (CUB) in M. chinensis Benth, chloroplast genome is essential to study its genetic law, molecular phylogenetic relationships, and exogenous gene expression. Results showed that the GC content of 53 CDS sequences was 37.95%, and GC1, GC2, and GC3 content were 46.02%, 38.26%, and 29.85%, respectively. The general GC content order was GC1>GC2>GC3. Moreover, the majority of genes had an effective number of codon (ENC) value greater than 40, except ndhE, rps8, and rps18. Correlation analysis results revealed that the GC content was significantly correlated with GC1, GC2, GC3, and ENC. Neutrality plot analysis, ENC-plot analysis, and PR2-plot analysis presented that the CUB of M. chinensis Benth chloroplast genome was mainly affected by mutation and selection. In addition, GGG, GCA, and TCC were found to be the optimal codons. Furthermore, results of cluster analysis and evolutionary tree showed that M. chinensis Benth was closely related to Ocimum basilicum, indicating that there was a certain correlation between the CUB of the chloroplast gene and the genetic relationship of plant species. Overall, the study on the CUB of chloroplast genome laid a basis for genetic modification and phylogenetic research of M. chinensis Benth chloroplast genome.

CUBAP: an interactive web portal for analyzing codon usage biases across populations

Synonymous codon usage significantly impacts translational and transcriptional efficiency, gene expression, the secondary structure of both mRNA and proteins, and has been implicated in various diseases. However, population-specific differences in codon usage biases remain largely unexplored. Here, we present a web server, https://cubap.byu.edu, to facilitate analyses of codon usage biases across populations (CUBAP). Using the 1000 Genomes Project, we calculated and visually depict population-specific differences in codon frequencies, codon aversion, identical codon pairing, co-tRNA codon pairing, ramp sequences, and nucleotide composition in 17,634 genes. We found that codon pairing significantly differs between populations in 35.8% of genes, allowing us to successfully predict the place of origin for African and East Asian individuals with 98.8% and 100% accuracy, respectively. We also used CUBAP to identify a significant bias toward decreased CTG pairing in the immunity related GTPase M (IRGM) gene in East Asian and African populations, which may contribute to the decreased association of rs10065172 with Crohn's disease in those populations. CUBAP facilitates in-depth gene-specific and codon-specific visualization that will aid in analyzing candidate genes identified in genome-wide association studies, identifying functional implications of synonymous variants, predicting population-specific impacts of synonymous variants and categorizing genetic biases unique to certain populations.

A code within the genetic code: codon usage regulates co-translational protein folding

The genetic code is degenerate, and most amino acids are encoded by two to six synonymous codons. Codon usage bias, the preference for certain synonymous codons, is a universal feature of all genomes examined. Synonymous codon mutations were previously thought to be silent; however, a growing body evidence now shows that codon usage regulates protein structure and gene expression through effects on co-translational protein folding, translation efficiency and accuracy, mRNA stability, and transcription. Codon usage regulates the speed of translation elongation, resulting in non-uniform ribosome decoding rates on mRNAs during translation that is adapted to co-translational protein folding process. Biochemical and genetic evidence demonstrate that codon usage plays an important role in regulating protein folding and function in both prokaryotic and eukaryotic organisms. Certain protein structural types are more sensitive than others to the effects of codon usage on protein folding, and predicted intrinsically disordered domains are more prone to misfolding caused by codon usage changes than other domain types. Bioinformatic analyses revealed that gene codon usage correlates with different protein structures in diverse organisms, indicating the existence of a codon usage code for co-translational protein folding. This review focuses on recent literature on the role and mechanism of codon usage in regulating translation kinetics and co-translational protein folding.

Comprehensive study of the genes involved in chlorophyll synthesis and degradation pathways in some monocot and dicot plant species

Chlorophyll (Chl) biosynthesis is one of the most important cellular processes essential for plant photosynthesis. Chl degradation pathway is also important catabolic process occurs during leaf senescence, fruit ripening and under biotic or abiotic stress conditions. Here we have systematically investigated the molecular evolution, gene structure, compositional analysis along with ENc plot, correspondence analysis and codon usage bias of the proteins and encoded genes involved in Chl metabolism from monocots and dicots. The gene and species specific phylogenetic trees using amino acid sequences showed clear clustering formation of the selected species based on monocots and dicots but not supported by 18S rRNA. Nucleotide composition of the encoding genes showed that average GC%, GC1%, GC2% and GC3% were higher in monocots. RSCU analysis depicts that genes from monocots for both pathways and genes for synthesis pathway from dicots only biased to G/C-ending synonymous codons but in degradation pathway most optimal codons (except UUG) in dicots biased to A/U-ending synonymous codons. We found strong evidence of episodic diversifying selection at several amino acid sites in all genes investigated. Conserved domain and gene structures were observed for the genes with varying lengths of introns and exons, involved in Chl metabolism along with some intronless genes within synthesis pathway. ENc and correspondence analyses suggested the mutational or selection constraint on the genes to shape the codon usage. These comprehensive studies may be helpful in further research in molecular phylogenetics and genomics and to better understand the evolutionary dynamics of Chl metabolic pathway.Communicated by Ramaswamy H. Sarma.

Comparison of the Whole-Plastome Sequence between the Bonin Islands Endemic Rubus boninensis and Its Close Relative, Rubus trifidus (Rosaceae), in the Southern Korean Peninsula

Rubus boninensis is a rare endemic species found on the Bonin Islands with a very restricted distribution. It is morphologically most closely related to Rubus trifidus, occurring widely in the southern Korean peninsula and Japan. This species pair provides a good example of anagenetic speciation on an oceanic island in the northwestern Pacific Ocean-R. trifidus as a continental progenitor and R. boninensis as an insular derivative species. In this study, we firstly characterized the complete plastome of R. boninensis and R. trifidus and compared this species pair to another anagenetically derived species pair (R. takesimensis-R. crataegifolius). The complete plastome of R. trifidus was 155,823 base pairs (bp) long, slightly longer (16 bp) than that of R. boninensis (155,807 bp). No structural or content rearrangements were found between the species pair. Eleven hotspot regions, including trnH/psbA, were identified between R. trifidus and R. boninensis. Phylogenetic analysis of 19 representative plastomes within the family Rosaceae suggested sister relationships between R. trifidus and R. boninensis, and between R. crataegifolius and R. takesimensis. The plastome resources generated by the present study will help elucidate plastome evolution and resolve phylogenetic relationships within highly complex and reticulated lineages of the genus Rubus.

In Silico Analyses of Burial Codon Bias Among the Species of Dipterocarpaceae Through Molecular and Phylogenetic Data

Introduction:

DNA barcode, a molecular marker, is used to distinguish among the closely related species, and it can be applied across a broad range of taxa to understand ecology and evolution. MaturaseK gene (matK) and rubisco bisphosphate carboxylase/oxygenase form I gene (rbcL) of the chloroplast are highly conserved in a plant system, which are used as core barcode. This present endeavor entails the comprehensive examination of the under threat plant species based on success of discrimination on DNA barcode under selection pressure.

Result:

The family Dipterocarpaceae comprising of 15 genera is under threat due to some factors, namely, deforestation, habitat alteration, poor seed, pollen dispersal, etc. Species of this family was grouped into 6 clusters for matK and 5 clusters and 2 sub-clusters for rbcL in the phylogenetic tree by using neighbor-joining method. Cluster I to cluster VI of matK and cluster I to cluster V of rbcL genes were analyzed by various codon and substitution bias tools. Mutational pressure guided the codon bias which was favored by the avoidance of higher GC content and significant negative correlation between GC12 and GC3 (in sub-cluster I of cluster I [0.03 < P], cluster I [0.00001 < P], and cluster II [0.01 < P] of rbcL, and cluster IV [0.013 < P] of matK). After refining the results, it could be speculated that the lower null expectation values (R = 0.5 or <0.5) were less divergent from the evolutionary perspective. Apart from that, the higher null expectation values (R = >0.85) also showed the same result, which possibly could be due to the negative impact of very high and low transition rate than transversion.

Conclusion:

Through the analysis of inter-generic, inter/intra-specific variation and phylogenetic data, it was found that both selection and mutation played an important role in synonymous codon choice in these genes, but they acted inconsistently on the genes, both matK and rbcL. In vitro stable proteins of both matK and rbcL were selected through natural selection rather than mutational selection. matK gene had higher individual discrimination and barcode success compared with rbcL. These discriminatory approaches may describe the problem related to the extinction of plant species. Hence, it becomes very imperative to identify and detect the under threat plant species in advance.

Factors affecting synonymous codon usage of housekeeping genes in Drosophila melanogaster

Housekeeping genes (HK genes) are required for cell survival and the maintenance of basic cellular functions. The investigation of factors affecting codon usage patterns in HK genes of insects can help in understanding the molecular evolution of insects and aid the development of insect pest management strategies. In this study, we employed bioinformatics approaches to analyze the codon usage bias (CUB) of HK genes in the insect model organism, Drosophila melanogaster. A comparison of CUB between 1107 HK genes and 1084 high tissue specificity genes suggested that HK genes have higher CUB in D. melanogaster. In addition, we found that CUB inversely correlates with the non-synonymous substitution rate of HK genes. Therefore, we attempted to identify the factors that potentially influence the codon usage pattern of HK genes. Our results suggest that mutation pressure and natural selection highly correlate with CUB in the HK genes of D. melanogaster and that two topological properties of HK proteins (proportion of protein interacting length and protein connectivity) also correlate with CUB in the HK genes of D. melanogaster. This study provides insight into CUB in the HK genes of D. melanogaster, and the results can support future investigations of potential applications in agricultural and biomedical field.

New insights into the codon usage patterns of the bactericidal/permeability-increasing (BPI) gene across nine species

Bactericidal/permeability-increasing (BPI) protein is a member of a new generation of proteins known as super-antibiotics that are implicated as endotoxin neutralising agents. Non-uniform usage of synonymous codons for a specific amino acid during translation of a protein is known as codon usage bias (CUB). Analysis of CUB and compositional dynamics of coding sequences could contribute to a better understanding of the molecular mechanism and the evolution of a particular gene. In this study, we performed CUB analysis of the complete coding sequences of the BPI gene from nine different species. The codon usage patterns of BPI across different species were found to be influenced by GC bias, particularly GC3s, with a moderate bias in the codon usage of BPI. We found significant similarities in the codon usage patterns in BPI gene among closely related species, such as Sus_scrofa and Bos_taurus. Moreover, we observed evolutionary conservation of the most over-represented codon CUG for the amino acid leucine in the BPI gene across all species. In conclusion, our analysis provides a novel insight into the codon usage patterns of BPI. This information facilitates an improved understanding of the structural, functional and evolutionary significance of BPI gene among species, and provides a theoretical reference for developing antiseptic drug proteins with high efficiency across species.

Genome-Wide Analysis of the Synonymous Codon Usage Patterns in Riemerella anatipestifer

Riemerella anatipestifer (RA) belongs to the Flavobacteriaceae family and can cause a septicemia disease in poultry. The synonymous codon usage patterns of bacteria reflect a series of evolutionary changes that enable bacteria to improve tolerance of the various environments. We detailed the codon usage patterns of RA isolates from the available 12 sequenced genomes by multiple codon and statistical analysis. Nucleotide compositions and relative synonymous codon usage (RSCU) analysis revealed that A or U ending codons are predominant in RA. Neutrality analysis found no significant correlation between GC12 and GC₃ (p > 0.05). Correspondence analysis and ENc-plot results showed that natural selection dominated over mutation in the codon usage bias. The tree of cluster analysis based on RSCU was concordant with dendrogram based on genomic BLAST by neighbor-joining method. By comparative analysis, about 50 highly expressed genes that were orthologs across all 12 strains were found in the top 5% of high CAI value. Based on these CAI values, we infer that RA contains a number of predicted highly expressed coding sequences, involved in transcriptional regulation and metabolism, reflecting their requirement for dealing with diverse environmental conditions. These results provide some useful information on the mechanisms that contribute to codon usage bias and evolution of RA.

Pangenome Evidence for Higher Codon Usage Bias and Stronger Translational Selection in Core Genes of Escherichia coli

Codon usage bias, as a combined interplay from mutation and selection, has been intensively studied in Escherichia coli. However, codon usage analysis in an E. coli pangenome remains unexplored and the relative importance of mutation and selection acting on core genes and strain-specific genes is unknown. Here we perform comprehensive codon usage analyses based on a collection of multiple complete genome sequences of E. coli. Our results show that core genes that are present in all strains have higher codon usage bias than strain-specific genes that are unique to single strains. We further explore the forces in influencing codon usage and investigate the difference of the major force between core and strain-specific genes. Our results demonstrate that although mutation may exert genome-wide influences on codon usage acting similarly in different gene sets, selection dominates as an important force to shape biased codon usage as genes are present in an increased number of strains. Together, our results provide important insights for better understanding genome plasticity and complexity as well as evolutionary mechanisms behind codon usage bias.

Analysis of synonymous codon usage in spike protein gene of infectious bronchitis virus

Infectious bronchitis virus (IBV) is responsible for causing respiratory, renal, and urogenital diseases in poultry. IBV infection in poultry leads to high mortality rates in affected flocks and to severe economic losses due to a drop in egg production and a reduced gain in live weight of the broiler birds. IBV-encoded spike protein (S) is the major protective immunogen for the host. Although the functions of the S protein have been well studied, the factors shaping synonymous codon usage bias and nucleotide composition in the S gene have not been reported yet. In the present study, we analyzed the relative synonymous codon usage and effective number of codons (Nc) using the 53 IBV S genes. The major trend in codon usage variation was studied using correspondence analysis. The plot of Nc values against GC3 as well as the correlation between base composition and codon usage bias suggest that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in the S gene. Interestingly, no association of aromaticity, degree of hydrophobicity, and aliphatic index was observed with the codon usage variation in IBV S genes. The study represents a comprehensive analysis of IBV S gene codon usage patterns and provides a basic understanding of the codon usage bias.

Compositional Constraint Is the Key Force in Shaping Codon Usage Bias in Hemagglutinin Gene in H1N1 Subtype of Influenza A Virus

It is vital to unravel the codon usage bias in order to gain insights into the evolutionary forces dictating the viral evolution process. Influenza A virus has attracted attention of many investigators over the years due to high mutation rate and being cross-specific shift operational in the viral genome. Several authors have reported that the codon usage bias is low in influenza A viruses, citing mutational pressure as the decisive force shaping up the codon usage in these viruses. In this study, complete coding sequences of hemagglutinin genes for H1N1 subtype of influenza A virus have been explored for the possible codon usage bias acting upon these genes. The results indicate overall low bias with peaking ENC values. The GC content is found to be substantially low as against AT content in the silent codon sites. Significant correlations were observed in between the compositional parameters versus AT₃, implying the possible role of the latter in shaping codon usage profile in the viral hemagglutinin. The data showed conspicuously that the sequences were A redundant with most codons preferring nucleotide A over others in the third synonymous codon site. The results indicated the pivotal role of compositional pressure affecting codon usage in this virus.

Analysis of codon use features of stearoyl-acyl carrier protein desaturase gene in Camellia sinensis

The stearoyl-acyl carrier protein desaturase (SAD) gene widely exists in all kinds of plants. In this paper, the Camellia sinensis SAD gene (CsSAD) sequence was firstly analyzed by Codon W, CHIPS, and CUSP programs online, and then compared with genomes of the tea plant, other species and SAD genes from 11 plant species. The results show that the CsSAD gene and the selected 73 of C. sinensis genes have similar codon usage bias. The CsSAD gene has a bias toward the synonymous codons with A and T at the third codon position, the same as the 73 of C. sinensis genes. Compared with monocotyledons such as Triticum aestivum and Zea mays, the differences in codon usage frequency between the CsSAD gene and dicotyledons such as Arabidopsis thaliana and Nicotiana tobacum are less. Therefore, A. thaliana and N. tobacum expression systems may be more suitable for the expression of the CsSAD gene. The analysis result of SAD genes from 12 plant species also shows that most of the SAD genes are biased toward the synonymous codons with G and C at the third codon position. We believe that the codon usage bias analysis presented in this study will be essential for providing a theoretical basis for discussing the structure and function of the CsSAD gene.

Targeted sequence capture and resequencing implies a predominant role of regulatory regions in the divergence of a sympatric lake whitefish species pair (Coregonus clupeaformis)

Latest technological developments in evolutionary biology bring new challenges in documenting the intricate genetic architecture of species in the process of divergence. Sympatric populations of lake whitefish represent one of the key systems to investigate this issue. Despite the value of random genotype-by-sequencing methods and decreasing cost of sequencing technologies, it remains challenging to investigate variation in coding regions, especially in the case of recently duplicated genomes as in salmonids, as this greatly complicates whole genome resequencing. We thus designed a sequence capture array targeting 2773 annotated genes to document the nature and the extent of genomic divergence between sympatric dwarf and normal whitefish. Among the 2728 genes successfully captured, a total of 2182 coding and 10,415 noncoding putative single-nucleotide polymorphisms (SNPs) were identified after applying a first set of basic filters. A genome scan with a quality-refined selection of 2203 SNPs identified 267 outlier SNPs in 210 candidate genes located in genomic regions potentially involved in whitefish divergence and reproductive isolation. We found highly heterogeneous FST estimates among SNP loci. There was an overall low level of coding polymorphism, with a predominance of noncoding mutations among outliers. The heterogeneous patterns of divergence among loci confirm the porous nature of genomes during speciation with gene flow. Considering that few protein-coding mutations were identified as highly divergent, our results, along with previous transcriptomic studies, imply that changes in regulatory regions most likely had a greater role in the process of whitefish population divergence than protein-coding mutations. This study is the first to demonstrate the efficiency of large-scale targeted resequencing for a nonmodel species with such a large and unsequenced genome.

An in silico analysis of troponin I mutations in hypertrophic cardiomyopathy of Indian origin

Hypertrophic Cardiomyopathy (HCM) is an autosomal dominant disorder of the myocardium which is hypertrophied resulting in arrhythmias and heart failure leading to sudden cardiac death (SCD). Several sarcomeric proteins and modifier genes have been implicated in this disease. Troponin I, being a part of the Troponin complex (troponin I, troponin C, troponin T), is an important gene for sarcomeric function. Four mutations (1 novel) were identified in Indian HCM cases, namely, Pro82Ser, Arg98Gln, Arg141Gln and Arg162Gln in Troponin I protein, which are in functionally significant domains. In order to analyse the effect of the mutations on protein stability and protein-protein interactions within the Troponin complex, an in silico study was carried out. The freely available X-ray crystal structure (PDB ID: 1JIE) was used as the template to model the protein followed by loop generation and development of troponin complex for both the troponin I wild type and four mutants (NCBI ID: PRJNA194382). The structural study was carried out to determine the effect of mutation on the structural stability and protein-protein interactions between three subunits in the complex. These mutations, especially the arginine to glutamine substitutions were found to result in local perturbations within the troponin complex by creating/removing inter/intra molecular hydrogen bonds with troponin T and troponin C. This has led to a decrease in the protein stability and loss of important interactions between the three subunits. It could have a significant impact on the disease progression when coupled with allelic heterogeneity which was observed in the cases carrying these mutations. However, this can be further confirmed by functional studies on protein levels in the identified cases.

Characterization of Synonymous Codon Usage Bias in the Riemerella Anatipestifer OmpA/MotB Gene

The analysis on codon usage bias of OmpA/MotB gene of Riemerella anatipestifer (RA) may provide a basis for understanding the evolution and pathogenesis of RA and for selecting appropriate host expression systems to improve the expression of target genes in vivo and in vitro. In our study, a comparative analysis of the codon usage bias in the newly discovered RA OmpA/MotB gene and the OmpA/MotB gene of 20 reference flavobacteriaceae was performed. The results of the codon adaptation indes (CAI), effective number of codon (ENC), and GC3s values indicated that synonymous codon usage bias in the OmpA/MotB gene of flavobacteriaceae. The results showed that codon usage bias in the RA OmpA/MotB gene was strong bias towards the synonymous codons with A and T at the third codon position. A high level of diversity in codon usage bias existed, and the effective number of codons used in a gene plot revealed that the G+C compositional constraint is the main factor that determines the codon usage bias in OmpA/MotB gene of flavobacteriaceae. Comparison of the codon usage in the OmpA/MotB gene of different organisms revealed that there were 31 codons showing distinct codon usage differences between the RA and E. coli, 41 between the RA and humans, but 29 between the RA and yeast. Therefore the yeast expression system may be more suitable for the expression of RA OmpA/MotB gene. These results may improve our understanding of the evolution, pathogenesis and functional studies of RA, as well as contribute significantly to the area of flavobacteriaceae research.

Coding nucleic acids are chaperons for protein folding: a novel theory of protein folding

The arguments for nucleic acid chaperons are reviewed and three new lines of evidence are added. (1) It was found that amino acids encoded by codons in short nucleic acid loops frequently form turns and helices in the corresponding protein structures. (2) The amino acids encoded by partially complementary (1st and 3rd nucleotides) codons are more frequently co-located in the encoded proteins than expected by chance. (3) There are significant correlations between thermodynamic changes (ddG) caused by codon mutations in nucleic acids and the thermodynamic changes caused by the corresponding amino acid mutations in the encoded proteins. We conclude that the concept of the Proteomic Code and nucleic acid chaperons seems correct from the bioinformatics point of view, and we expect to see direct biochemical experiments and evidence in the near future.

[Analysis of codon use features of CBF gene in Camellia sinensis]

CBF (C-repeat-binding factor) transcription factor exists widely in all kinds of plants. It is an important regulative factor in the process of plant resistance adversity. In this paper, Camellia sinensis CBF1 gene sequence was analyzed by Codon W, CHIPS, and CUSP programs online, and then compared with C. sinensis genes, genomes in other species, and CBF genes from 39 plant species. It is important to identify the codon usage of CsCBF1 gene and select appropriate expression systems. The results showed that CsCBF1 gene and selected 70 C. sinensis genes had distinct usage differences. CsCBF1 gene was bias toward the synonymous codons with G and C at the third codon position, but 70 C. sinensis genes were bias toward the synonymous codons with A and T. The differences in codon usage frequency between CsCBF1 gene and dicotyledons such as Arabidopsis thaliana and Nicotiana tobacum were less than monocotyledons such as wheat (Triticum aestivum) and corn (Zea mays). Therefore, A. thaliana and N. tobacum expression systems may be more suitable for the expression of CsCBF1 gene. The analysis results of CBF genes from 40 plant species also showed that most of the CBF genes were bias toward the synonymous codons with G and C at the third codon position. The reason of this phenomenon is possible due to special functions of these genes.

PXR and CAR single nucleotide polymorphisms influence plasma efavirenz levels in South African HIV/AIDS patients

Background

This study investigated variation in NR1I2 and NR1I3 and its effect on plasma efavirenz levels in HIV/AIDS patients. Variability in plasma drug levels has largely led research on identifying causative variants in drug metabolising enzyme (DME) genes, with little focus on the nuclear receptor genes NR1I2 and NR1I3, coding for PXR and CAR, respectively, that are involved in regulating DMEs.

Methods

464 Bantu-speaking South Africans comprising of HIV/AIDS patients on efavirenz-based treatment (n=301) and 163 healthy subjects were genotyped for 6 SNPs in NR1I2 and NR1I3. 32 of the 301 patients had their DNA binding domains (DBDs) in NR1I2 and NR1I3 sequenced.

Results

Significantly decreased efavirenz plasma concentrations were observed in patients carrying the NR1I3 rs3003596C/C and T/C genotypes (P=0.015 and P=0.010, respectively). Sequencing resulted in the discovery of a further 13 SNPs, 3 of which are novel variants in the DBD of NR1I2. There were significant differences in the distribution of NR1I2 and NR1I3 SNPs between South Africans when compared to Caucasian, Asian and Yoruba population groups.

Conclusion

For the realisation of personalised medicine, PXR and CAR genetic variation should be taken into consideration because of their involvement in the regulation of DMEs.

ABCB1 4036A>G and 1236C>T Polymorphisms Affect Plasma Efavirenz Levels in South African HIV/AIDS Patients

The ABCB1 gene encodes P-glycoprotein, an ATP-dependent drug efflux pump, which is responsible for drug transport across extra- and intra-cellular membranes. The variability in the expression of ABCB1 may contribute to variable plasma efavirenz concentration which results in variability in the levels of suppression of the human immunodeficiency syndrome virus (HIV). The aim of the study was to evaluate the role of polymorphisms in ABCB1 gene on plasma efavirenz levels and treatment response in the form of change in viral load and CD-4 cell count in HIV/AIDS patients receiving efavirenz-containing highly active antiretroviral treatment regimens. Two hundred and eighty-two HIV-infected patients were recruited from Themba Lethu Clinic in Johannesburg and plasma efavirenz drug concentration levels were measured using LC-MS/MS. SNaPshot was used to genotype five known ABCB1 single nucleotide polymorphisms (SNPs). Genotype-phenotype correlations were computed. The ABCB1 4036A/G and 4036G/G genotypes were significantly associated with low plasma efavirenz concentrations (P = 0.0236), while the ABCB1 1236C/T and 1236T/T genotypes were associated with high efavirenz concentrations (P = 0.0282). A haplotype ABCB1 T-G-T-A is reported that is associated with significantly increased plasma efavirenz levels. This is the first report on 61A>G, 2677G>T/A, and 4036A>G SNPs in the South African population. ABCB1 plays a role in determining the plasma concentrations of efavirenz and should be taken into account in future design of assays for genotype-based dosing of efavirenz-containing regimens.

High prevalence of Arginine to Glutamine substitution at 98, 141 and 162 positions in Troponin I (TNNI3) associated with hypertrophic cardiomyopathy among Indians

BACKGROUND

Troponin I (TNNI3) is the inhibitory subunit of the thin filament regulatory complex Troponin, which confers calcium-sensitivity to striated muscle actomyosin ATPase activity. Mutations (2-7%) in this gene had been reported in hypertrophic cardiomyopathy patients (HCM). However, the frequencies of mutations and associated clinical presentation have not been established in cardiomyopathy patients of Indian origin, hence we have undertaken this study.

METHODS

We have sequenced all the exons, including the exon-intron boundaries of TNNI3 gene in 101 hypertrophic cardiomyopathy patients (HCM), along with 160 healthy controls, inhabited in the same geographical region of southern India.

RESULTS

Our study revealed a total of 16 mutations. Interestingly, we have observed Arginine to Glutamine (R to Q) mutation at 3 positions 98, 141 and 162, exclusively in HCM patients with family history of sudden cardiac death. The novel R98Q was observed in a severe hypertrophic obstructive cardiomyopathy patient (HOCM). The R141Q mutation was observed in two familial cases of severe asymmetric septal hypertrophy (ASH++). The R162Q mutation was observed in a ASH++ patient with mean septal thickness of 29 mm, and have also consists of allelic heterogeneity by means of having one more synonymous (E179E) mutation at g.4797: G → A: in the same exon 7, which replaces a very frequent codon (GAG: 85%) with a rare codon (GAA: 14%). Screening for R162Q mutation in all the available family members revealed its presence in 9 individuals, including 7 with allelic heterogeneity (R162Q and E179E) of which 4 were severely affected. We also found 2 novel SNPs, (g.2653; G → A and g.4003 C → T) exclusively in HCM, and in silico analysis of these SNPs have predicted to cause defect in recognition/binding sites for proteins responsible for proper splicing.

CONCLUSION

Our study has provided valuable information regarding the prevalence of TNNI3 mutations in Indian HCM patients and its risk assessment, these will help in genetic counseling and to adopt appropriate treatment strategies.

Insight into pattern of codon biasness and nucleotide base usage in serotonin receptor gene family from different mammalian species

5-HT (5-Hydroxy-tryptamine) or serotonin receptors are found both in central and peripheral nervous system as well as in non-neuronal tissues. In the animal and human nervous system, serotonin produces various functional effects through a variety of membrane bound receptors. In this study, we focus on 5-HT receptor family from different mammals and examined the factors that account for codon and nucleotide usage variation. A total of 110 homologous coding sequences from 11 different mammalian species were analyzed using relative synonymous codon usage (RSCU), correspondence analysis (COA) and hierarchical cluster analysis together with nucleotide base usage frequency of chemically similar amino acid codons. The mean effective number of codon (ENc) value of 37.06 for 5-HT(6) shows very high codon bias within the family and may be due to high selective translational efficiency. The COA and Spearman's rank correlation reveals that the nucleotide compositional mutation bias as the major factors influencing the codon usage in serotonin receptor genes. The hierarchical cluster analysis suggests that gene function is another dominant factor that affects the codon usage bias, while species is a minor factor. Nucleotide base usage was reported using Goldman, Engelman, Stietz (GES) scale reveals the presence of high uracil (>45%) content at functionally important hydrophobic regions. Our in silico approach will certainly help for further investigations on critical inference on evolution, structure, function and gene expression aspects of 5-HT receptors family which are potential antipsychotic drug targets.

Patterns and influencing factor of synonymous codon usage in porcine circovirus

Background

Analysis of codon usage can reveal much about the molecular evolution of the viruses. Nevertheless, little information about synonymous codon usage pattern of porcine circovirus (PCV) genome in the process of its evolution is available. In this study, to give a new understanding on the evolutionary characteristics of PCV and the effects of natural selection from its host on the codon usage pattern of the virus, Patterns and the key determinants of codon usage in PCV were examined.

Methods

We carried out comprehensive analysis on codon usage pattern in the PCV genome, by calculating relative synonymous codon usage (RSCU), effective number of codons (ENC), dinucleotides and nucleic acid content of the PCV genome.

Results

PCV genomes have relatively much lower content of GC and codon preference, this result shows that nucleotide constraints have a major impact on its synonymous codon usage. The results of the correspondence analysis indicate codon usage patterns of PCV of various genotypes, various subgenotypes changed greatly, and significant differences in codon usage patterns of Each virus of Circoviridae.There is much comparability between PCV and its host in their synonymous codon usage, suggesting that the natural selection pressure from the host factor also affect the codon usage patterns of PCV. In particular, PCV genotype II is in synonymous codon usage more similar to pig than to PCV genotype I, which may be one of the most important molecular mechanisms of PCV genotype II to cause disease. The calculations results of the relative abundance of dinucleotides indicate that the composition of dinucleotides also plays a key role in the variation found in synonymous codon usage in PCV. Furthermore, geographic factors, the general average hydrophobicity and the aromaticity may be related to the formation of codon usage patterns of PCV.

Conclusion

The results of these studies suggest that synonymous codon usage pattern of PCV genome are the result of interaction between mutation pressure and natural selection from its host. The information from this study may not only have theoretical value in understanding the characteristics of synonymous codon usage in PCV genomes, but also have significant value for the molecular evolution of PCV.

Analysis of synonymous codon usage inShigella flexneri2a strain 301 and otherShigellaandEscherichia colistrains

In this study, we analysed synonymous codon usage in Shigella flexneri 2a strain 301 (Sf301) and performed a comparative analysis of synonymous codon usage patterns in Sf301 and other strains of Shigella and Escherichia coli . Although there was a significant variety in codon usage bias among different Sf301 genes, there was a slight but observable codon usage bias that could primarily be attributable to mutational pressure and translational selection. In addition, the relative abundance of dinucleotides in Sf301 was observed to be independent of the overall base composition but was still caused by differential mutational pressure; this also shaped codon usage. By comparing the relative synonymous codon usage values across different Shigella and E. coli strains, we suggested that the synonymous codon usage pattern in the Shigella genomes was strain specific. This study represents a comprehensive analysis of Shigella codon usage patterns and provides a basic understanding of the mechanisms underlying codon usage bias.

Factors affecting synonymous codon usage bias in chloroplast genome of oncidium gower ramsey

Oncidium Gower Ramsey is a fascinating and important ornamental flower in floral industry. In this research, the complete nucleotide sequence of the chloroplast genome in Oncidium Gower Ramsey was studied, then analyzed using Codonw software. Correspondence analysis and method of effective number of codon as Nc-plot were conducted to analyze synonymous codon usage. According to the corresponding analysis, codon bias in the chloroplast genome of Oncidium Gower Ramsey is related to their gene length, mutation bias, gene hydropathy level of each protein, gene function and selection or gene expression only subtly affect codon usage. This study will provide insights into the molecular evolution study and high-level transgene expression.

Analysis of codon usage in Newcastle disease virus

In this study, the relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents, and dinucleotide were used to investigate codon usage pattern of each protein-coding gene and genome among 31 Newcastle disease virus (NDV) isolates. The result shows that the overall extent of codon usage bias in NDV is low (mean ENC = 56.15 > 40). The good correlation between the (C + G)₁₂% and (G + C)₃% suggests that the mutational pressure, rather than natural selection, is the main factor that determines the codon usage bias and base component in NDV. It is observed that synonymous codon usage pattern in NDV genes is gene function and geography specific, but not host specific. By contrasting synonymous codon usage patterns of different NDV isolates, we suggest that more than one genotype of NDV circulates in waterfowl in USA; and gene length has no significant effect on the variations of synonymous codon usage in these virus genes. CpG under-represented is a characteristic for NDV to fit in its host. These results not only provide an insight into the variation of codon usage pattern among the genomes of NDV, but also may help in understanding the processes governing the evolution of NDV.

Analysis of synonymous codon usage bias in 09H1N1

A novel subtype of influenza A virus 09H1N1 has rapidly spread across the world. Evolutionary analyses of this virus have revealed that 09H1N1 is a triple reassortant of segments from swine, avian and human influenza viruses. In this study, we investigated factors shaping the codon usage bias of 09H1N1 and carried out cluster analysis of 60 strains of influenza A virus from different subtypes based on their codon usage bias. We discovered that more preferentially used codons of 09H1N1 are A-ended or U-ended, and the intra-genomic codon usage bias of 09H1N1 is quite low. Base composition constraint, dinucleotide biases and translational selection are the main factors influencing the codon usage bias of 09H1N1. At the genome level, we find that the codon usage bias of 09H1N1 is similar to H1N1 (A/swine/Kansas/77778/2007H1N1), H9N2 from Asia, H1N2 from Asia and North America and H3N2 from North America. Our results provide insight for understanding the processes governing evolution, regulation of gene expression, and revealing the evolution of 09H1N1.

Codon Usage Patterns in Corynebacterium glutamicum: Mutational Bias, Natural Selection and Amino Acid Conservation

The alternative synonymous codons in Corynebacterium glutamicum, a well-known bacterium used in industry for the production of amino acid, have been investigated by multivariate analysis. As C. glutamicum is a GC-rich organism, G and C are expected to predominate at the third position of codons. Indeed, overall codon usage analyses have indicated that C and/or G ending codons are predominant in this organism. Through multivariate statistical analysis, apart from mutational selection, we identified three other trends of codon usage variation among the genes. Firstly, the majority of highly expressed genes are scattered towards the positive end of the first axis, whereas the majority of lowly expressed genes are clustered towards the other end of the first axis. Furthermore, the distinct difference in the two sets of genes was that the C ending codons are predominate in putatively highly expressed genes, suggesting that the C ending codons are translationally optimal in this organism. Secondly, the majority of the putatively highly expressed genes have a tendency to locate on the leading strand, which indicates that replicational and transciptional selection might be invoked. Thirdly, highly expressed genes are more conserved than lowly expressed genes by synonymous and nonsynonymous substitutions among orthologous genes fromthe genomes of C. glutamicum and C. diphtheriae. We also analyzed other factors such as the length of genes and hydrophobicity that might influence codon usage and found their contributions to be weak.