Analysis of synonymous codon usage in H5N1 virus and other influenza A viruses

Evolution shapes and conserves genomic signatures in viruses

The genomic signature of an organism captures the characteristics of repeated oligonucleotide patterns in its genome 1, such as oligomer frequencies, GC content, and differences in codon usage. Viruses, however, are obligate intracellular parasites that are dependent on their host cells for replication, and information about genomic signatures in viruses has hitherto been sparse.Here, we investigate the presence and specificity of genomic signatures in 2,768 eukaryotic viral species from 105 viral families, aiming to illuminate dependencies and selective pressures in viral genome evolution. We demonstrate that most viruses have highly specific genomic signatures that often also differ significantly between species within the same family. The species-specificity is most prominent among dsDNA viruses and viruses with large genomes. We also reveal consistent dissimilarities between viral genomic signatures and those of their host cells, although some viruses present slight similarities, which may be explained by genetic adaptation to their native hosts. Our results suggest that significant evolutionary selection pressures act upon viral genomes to shape and preserve their genomic signatures, which may have implications for the field of synthetic biology in the construction of live attenuated vaccines and viral vectors.

Compositional Features and Codon Usage Pattern of Genes Associated with Parkinson's Disease

Codon usage bias (CUB) is the phenomenon of non-uniform usage of synonymous codons in which some codons are more used than others and it helps in understanding the molecular organization of genome. Bioinformatic approach was used to analyze the protein-coding sequences of genes associated with Parkinson's disease (PD) to explore compositional features and codon usage pattern as no details work was reported yet. The average improved effective number of codons (Nc) and Nc prime were 42.74 and 44.26 respectively, indicated that CUB was low in these genes. In most of the genes, the overall GC content was almost 50% and GC content at the 1st codon position was the highest while GC content at the 2nd codon position was lowest. Relative synonymous codon usage (RSCU) analysis elucidated over-represented (p > 1.6) and under-represented codons (p < 0.6). The GTG (Val) is the only codon over-represented in all genes. Over-represented codons except (GTG) were A or T ending while under-represented codons (except ACT) were G or C ending. The codons namely TTA (Leu), CTA (Leu), ATC (Ile), ATA (Ile), AGT (Ser), AAC (Asn), TGT (Cys), TGC (Cys), CGC (Arg), AGA (Arg), and AGG (Arg) were absent in SNCA1 to SNCA8 genes. The codon TCG (Ser) was absent in all genes except UCHL1 and PINK1. Correspondence analysis (COA) revealed that the pattern of codon usage differs among genes associated with PD. Neutrality plot analysis indicated some of the points are diagonal distribution suggested that mutation pressure influenced the CUB in genes associated with PD.

Comprehensive Genomics Investigation of Neboviruses Reveals Distinct Codon Usage Patterns and Host Specificity

Neboviruses (NeVs) from the Caliciviridae family have been linked to enteric diseases in bovines and have been detected worldwide. As viruses rely entirely on the cellular machinery of the host for replication, their ability to thrive in a specific host is greatly impacted by the specific codon usage preferences. Here, we systematically analyzed the codon usage bias in NeVs to explore the genetic and evolutionary patterns. Relative Synonymous Codon Usage and Effective Number of Codon analyses indicated a marginally lower codon usage bias in NeVs, predominantly influenced by the nucleotide compositional constraints. Nonetheless, NeVs showed a higher codon usage bias for codons containing G/C at the third codon position. The neutrality plot analysis revealed natural selection as the primary factor that shaped the codon usage bias in both the VP1 (82%) and VP2 (57%) genes of NeVs. Furthermore, the NeVs showed a highly comparable codon usage pattern to bovines, as reflected through Codon Adaptation Index and Relative Codon Deoptimization Index analyses. Notably, yak NeVs showed considerably different nucleotide compositional constraints and mutational pressure compared to bovine NeVs, which appear to be predominantly host-driven. This study sheds light on the genetic mechanism driving NeVs' adaptability, evolution, and fitness to their host species.

Comprehensive examination on codon usage bias pattern of the Bovine Ephemeral fever virus

Bovine Ephemeral Fever Virus (BEFV) is a non-contagious virus that commonly infects cattle and water buffalo, reduces milk productivity, decreases the quality of beef, and causes an adverse economic impact on the global livestock industry. However, the evolution of BEFV is unclear, and uncertainty exists regarding its global geodynamics. Consequently, this study aims to comprehend the pattern of viral evolution and gene expression in the BEFV genes G, M, N, and P, including synonymous codons. Additionally, we performed recombination analyses, which exclusively detected recombination signals in the G- and P-genes. Subsequently, a phylogenetic tree was constructed to validate and support these findings. The codon usage bias results showed that the BEFV-selected genes were influenced by both natural and mutation pressure. Furthermore, nucleotide A is more abundant in all the selected genes. The eNC values, ranging from 42.99 to 47.10, revealed the presence of moderate codon usage bias, where gene P exhibited the highest and gene G had the lowest codon usage bias. The neutrality and PR-2 plots, specified codon usage patterns of the genes, are also being shaped by strong selectional pressure. This comprehensive analysis of BEFV genes (G, M, N, and P) sheds light on the molecular evolutionary patterns, co-adaptation, and different genes expression in diverse regions, facilitating the development of preventative programs and insights into viral pathogenesis and vaccine design.Communicated by Ramaswamy H. Sarma.

Codon Usage Bias Analysis in Macronuclear Genomes of Ciliated Protozoa

Ciliated protozoa (ciliates) are unicellular eukaryotes, several of which are important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome of ciliates can promote a better understanding of the genetic mode and evolutionary history of these organisms and help optimize codons to improve gene editing efficiency in model ciliates. In this study, the following indices were calculated: the guanine-cytosine (GC) content, the frequency of the nucleotides at the third position of codons (T3, C3, A3, G3), the effective number of codons (ENc), GC content at the 3rd position of synonymous codons (GC3s), and the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, Neutrality plot analysis, ENc plot analysis, and correlation analysis were employed to explore the main influencing factors of CUB. The results showed that the GC content in the MAC genomes of each of 21 ciliate species, the genomes of which were relatively complete, was lower than 50%, and the base compositions of GC and GC3s were markedly distinct. Synonymous codon analysis revealed that the codons in most of the 21 ciliates ended with A or T and four codons were the general putative optimal codons. Collectively, our results indicated that most of the ciliates investigated preferred using the codons with anof AT-ending and that codon usage bias was affected by gene mutation and natural selection.

Insights into the Epidemiology, Phylodynamics, and Evolutionary Changes of Lineage GI-7 Infectious Bronchitis Virus

Infectious bronchitis virus (IBV) is distributed worldwide and causes significant losses in the poultry industry. In recent decades, lineages GI-19 and GI-7 have become the most prevalent IBV strains in China. However, the molecular evolution and phylodynamics of the lineage GI-7 IBV strains remain largely unknown. In this study, we identified 19 IBV strains from clinical samples from January 2021 to June 2022 in China, including 12 strains of GI-19, 3 strains of GI-7, and 1 strain each of GI-1, GI-9, GI-13, and GI-28. These results indicated that lineages GI-19 and GI-7 IBVs are still the most prevalent IBVs in China. Here, we investigated the evolution and transmission dynamics of lineage GI-7 IBVs. Our results revealed that the Taiwan province might be the origin of lineage GI-7 IBVs and that South China plays an important role in the spread of IBV. Furthermore, we found low codon usage bias of the S1 gene in lineage GI-7 IBVs. This allowed IBV to replicate in the host during evolution as a result of reduced competition, mainly driven by natural selection and mutational pressure, where the role of natural selection is more prominent. Collectively, our results reveal the genetic diversity and evolutionary dynamics of lineage GI-7 IBVs, which could assist in the prevention and control of viral infection.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Codon Pattern and Compositional Constraints Determination of Genes Associated with Chronic Periodontitis

Genome-wide association studies showed the relationship of NIN, ABHD12B, WHAMM, AP3B2, and SIGLEC5 with chronic periodontitis. The study’s objective was to investigate different molecular patterns and evolutionary forces acting on the mentioned genes. The investigation of molecular patterns encompasses the study of compositional parameters, expression profile, physical properties of genes, codon preferences, degree of codon bias, determination of the most influential codons, and assessment of actions of evolutionary forces, such as mutations and natural selection. The overall compositional analysis revealed the dominance of A and G nucleotides compared to T and C. A relatively low codon usage bias is observed. The CTG codon is the most overused codon, followed by TCC. The genes, AP3B2 and SIGLEC5, preferred GC-ending codons, while NIN, ABHD12B, and WHAMM preferred AT-ending codons. The presence of directional mutational force and natural selection was found to operate codon usage in genes envisaged, and selective forces were dominant over mutational forces. Apart from mutation and selection forces, compositional constraints also played imperative roles. The study enriched our knowledge of specific molecular patterns associated with the set of genes significantly associated with chronic periodontitis. Further studies are warranted to identify more genetic signatures associated with the disease.

Analysis of the Compositional Features and Codon Usage Pattern of Genes Involved in Human Autophagy

Autophagy plays an intricate role in paradigmatic human pathologies such as cancer, and neurodegenerative, cardiovascular, and autoimmune disorders. Autophagy regulation is performed by a set of autophagy-related (ATG) genes, first recognized in yeast genome and subsequently identified in other species, including humans. Several other genes have been identified to be involved in the process of autophagy either directly or indirectly. Studying the codon usage bias (CUB) of genes is crucial for understanding their genome biology and molecular evolution. Here, we examined the usage pattern of nucleotide and synonymous codons and the influence of evolutionary forces in genes involved in human autophagy. The coding sequences (CDS) of the protein coding human autophagy genes were retrieved from the NCBI nucleotide database and analyzed using various web tools and software to understand their nucleotide composition and codon usage pattern. The effective number of codons (ENC) in all genes involved in human autophagy ranges between 33.26 and 54.6 with a mean value of 45.05, indicating an overall low CUB. The nucleotide composition analysis of the autophagy genes revealed that the genes were marginally rich in GC content that significantly influenced the codon usage pattern. The relative synonymous codon usage (RSCU) revealed 3 over-represented and 10 under-represented codons. Both natural selection and mutational pressure were the key forces influencing the codon usage pattern of the genes involved in human autophagy.

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.

The genetic evolution and codon usage pattern of severe fever with thrombocytopenia syndrome virus

Severe fever with thrombocytopenia syndrome (SFTS) is a newly emerging zoonotic infectious disease caused by the SFTS virus (SFTSV), which has been continuously circulating in Eastern Asia in recent years. Although the evolution of SFTSV has been investigated, the evolutionary changes associated with codon usage have not been reported. Thus, a comprehensive genetic and codon usage bias analysis of SFTSV was conducted to elucidate the genetic diversity and evolutionary relationships in a novel perspective. The study amplified and sequenced fifteen SFTSV strains from a prefecture of Zhejiang Province, Eastern China in 2020, where SFTS cases have been continuously reported in the past decade. Phylogenetic analysis was conducted based on the complete coding sequences of SFTSV segments. It suggested that all SFTSV strains circulating in Zhejiang were clustered with Japanese and Korean strains, which belonged to two different genotypes. Meanwhile, thirty-nine genetic reassortants classified into nineteen different reassortment forms were identified, while 45 recombination events in 41 SFTSV strains were found. Codon usage patterns were further analyzed to understand the evolutionary changes in relation to genotype and host. And it revealed that codon usage bias was mainly driven by natural selection rather than mutation pressure. In addition, the codon adaptation index (CAI) analysis demonstrated the strong adaptability of SFTSV to Gallus gallus and Homo sapiens. Similarity index (SiD) analysis indicated that Haemaphysalis longicornis posed a strong selection pressure to SFTSV. In conclusion, this study revealed that the genetic diversity of SFTSV is gradually increasing. The codon usage analysis suggested that codon usage bias of SFTSV was mainly driven by natural selection, and SFTSV has evolved host-specific codon usage patterns. This contributes to the development of control measures against SFTSV.

Influencing elements of codon usage bias in Birnaviridae and its evolutionary analysis

Birnaviridae is a family of double stranded (ds) RNA virus with non-enveloped virions and 2-segmented genomes. These viruses are known to cause diseases in many hosts. Virus of this family has affected the fish and poultry economy in a wide sector. Unevenness in the use of synonymous codons for a particular amino acid in the coding strand of DNA is known as codon usage bias (CUB). Codons that code the same amino acid are used with variable frequency in a variety of life forms. To understand the pattern of CUB in Birnaviridae, we carried out bioinformatics study to understand the properties of coding sequences of proteins. ENC value of Birnaviridae suggested low CUB. Nucleotide analysis revealed high GC content. Parameters such as RSCU values, nucleotide skewness, translational selection, parity plot and neutrality plot were studied to investigate the pattern of codon use and it was clear that both mutational pressure and natural selection contributed to the designing of CUB in Birnaviridae family. The neutrality plot revealed natural selection to dominate the structuring of CUB and hence remained the major CUB determinant in Birnaviridae. Outcome of our study exemplified the pattern of codon use in the Birnaviridae genomes and contributed the basic primary data for fundamental evolutionary research on them.

Aldahish A, Uddin A. Analysis of Codon Usage of Speech Gene FoxP2 among Animals

The protein-coding gene FoxP2 (fork head box protein P2) plays a major role in communication and evolutionary changes. The present study carried out a comprehensive codon usage bias analysis in the FoxP2 gene among a diverse group of animals including fishes, birds, reptiles, and mammals. We observed that in the genome of fishes for the FoxP2 gene, codons ending with C or G were most frequently used, while in birds, reptiles, and mammals, codons ending with T or A were most frequently used. A higher ENC value was observed for the FoxP2 gene indicating a lower CUB. Parity role two-bias plots suggested that apart from mutation pressure, other factors such as natural selection might have influenced the CUB. The frequency distribution of the ENC observed and ENC expected ratio revealed that mutation pressure plays a key role in the patterns of codon usage of FoxP2. Besides, correspondence analysis exposed the composition of the nucleobase under mutation bias affects the codon usage of the FoxP2 gene. However, neutrality plots revealed the major role of natural selection over mutation pressure in the CUB of FoxP2. In addition, the codon usage patterns for FoxP2 among the selected genomes suggested that nature has favored nearly all the synonymous codons for encoding the corresponding amino acid. The uniform usage of 12 synonymous codons for FoxP2 was observed among the species of birds. The amino acid usage frequency for FoxP2 revealed that the amino acids Leucine, Glutamine, and Serine were predominant over other amino acids among all the species of fishes, birds, reptiles, and mammals.

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.

Insight into Codon Utilization Pattern of Tumor Suppressor Gene EPB41L3 from Different Mammalian Species Indicates Dominant Role of Selection Force

Uneven codon usage within genes as well as among genomes is a usual phenomenon across organisms. It plays a significant role in the translational efficiency and evolution of a particular gene. EPB41L3 is a tumor suppressor protein-coding gene, and in the present study, the pattern of codon usage was envisaged. The full-length sequences of the EPB41L3 gene for the human, brown rat, domesticated cattle, and Sumatran orangutan available at the NCBI were retrieved and utilized to analyze CUB patterns across the selected mammalian species. Compositional properties, dinucleotide abundance, and parity analysis showed the dominance of A and G whilst RSCU analysis indicated the dominance of G/C-ending codons. The neutrality plot plotted between GC12 and GC3 to determine the variation between the mutation pressure and natural selection indicated the dominance of selection pressure (R = 0.926; p < 0.00001) over the three codon positions across the gene. The result is in concordance with the codon adaptation index analysis and the ENc-GC3 plot analysis, as well as the translational selection index (P2). Overall selection pressure is the dominant pressure acting during the evolution of the EPB41L3 gene.

Composition, codon usage pattern, protein properties, and influencing factors in the genomes of members of the family Anelloviridae

The present study was carried out on 62 genome sequences of members of the family Anelloviridae, as there have been no reports of genome analysis of these DNA viruses using a bioinformatics approach. The genes were found to be rich in AC content with low codon usage bias (CUB). Relative synonymous codon usage (RSCU) values identified the preferred codons for each amino acid in the family. The codon AGA was overrepresented, while the codons TCG, TTG, CGG, CGT, ACG, GCG and GAT were underrepresented in all of the genomes. A significant correlation was found between the effective number of codons (ENC) and base constraints, indicating that compositional properties might have influenced the CUB. A highly significant correlation was observed between the overall base content and the base content at the third codon position, indicating that mutations might have affected the CUB. A highly significant positive correlation was observed between GC12 and GC3 (r = 0.904, p < 0.01), which indicated that directional mutation pressure influenced all three codon positions. A neutrality plot revealed that the contribution of mutation and natural selection in determining the CUB was 58.6% and 41.4%, respectively.

Analysis of mitochondrial protein-coding genes of Antheraea assamensis: Muga silkworm of Assam

To understand the synonymous codon usage pattern in mitochondrial genome of Antheraea assamensis, we analyzed the 13 mitochondrial protein-coding genes of this species using a bioinformatic approach as no work was reported yet. The nucleotide composition analysis suggested that the percentages of A, T, G,and C were 33.73, 46.39, 9.7 and 10.17, respectively and the overall GC content was 19.86, that is, lower than 50% and the genes were AT rich. The mean effective number of codons of mitochondrial protein-coding genes was 36.30 and it indicated low codon usage bias (CUB). Relative synonymous codon usage analysis suggested overrepresented and underrepresented codons in each gene and the pattern of codon usage was different among genes. Neutrality plot analysis revealed a narrow range of distribution for GC content at the third codon position and some points were diagonally distributed, suggesting both mutation pressure and natural selection influenced the CUB.

Human Schlafen 11 exploits codon preference discrimination to attenuate viral protein synthesis of prototype foamy virus (PFV)

Recently, the Schlafen (SLFN) proteins have been identified as a novel interferon-stimulated family with antiviral properties. In this study, we reported that SLFN11 inhibited prototype foamy virus (PFV) replication. Over-expression of human SLFN11 reduced viral production, while knockdown of SLFN11 enhanced viral infectivity. In addition, SLFN11 from cattle and African green monkey also suppressed PFV production. Both the ATPase activity and helicase activity of SLFN11 were required for its inhibitory function. Dephosphorylation activated the antiviral activity of SLFN11. More importantly, SLFN11 inhibited the expression of viral protein, which was rescued by viral gene codon optimization. Together, our results demonstrated that SLFN11 impaired PFV viral protein synthesis by exploiting the distinct codon usage between the virus and the host. These findings further broaden our understanding of the antiviral properties of the SLFN family and the molecular mechanism of PFV latent infection.

Genome-wide analysis of codon usage pattern in herpesviruses and its relation to evolution

The preferential use of a specific codon, out of a group of synonymous codons encoding the same amino acid, in a gene transcript results from the bias in codon choice. Various evolutionary forces namely mutation pressure and natural selection influence the pattern of codon usage i.e. distinct for each gene/genome. We investigated the pattern of codon usage of eight human herpesvirus genomes and compared them with two other herpesvirus genomes namely murine herpesvirus 68 and bovine herpesvirus type 1.1 to elucidate its compositional features, pattern of codon usage across the genomes and report the differences of codon usage pattern of human herpesviruses from that of other two other viruses. We also identified the similarity of the codon usage of human herpesviruses with its host (human). The genes were found to be CG rich in HHV2, HHV3, HHV4, HHV6, HHV7 and BH genomes while TA rich in HHV1, HHV5, HHV8 and MH genomes. The codon usage bias (CUB) of genes was low. A highly significant correlation was found among compositional contents depicting the role of mutational pressure along with natural selection in framing CUB. Several more frequently used codons as well as less frequently used codons were identified to be similar between each human virus and its host (human), while murine herpesvirus 68 and bovine herpesvirus type 1.1 genomes did not possess similar adaptation strategy as human herpesviruses to human (host), thus we could conclude that viral CUB might have been shaped as per their host's nature for better surveillance. Neutrality plot revealed mutational pressure mostly influenced the CUB of HHV1, HHV8 and MH viruses, while natural selection had a major impact in the CUB of HHV2, HHV3, HHV4, HHV5, HHV6, HHV7 and BH genomes.

Insights into Genomic Epidemiology, Evolution, and Transmission Dynamics of Genotype VII of Class II Newcastle Disease Virus in China

Newcastle disease virus (NDV) is distributed worldwide and has caused significant losses to the poultry industry. Almost all virulent NDV strains belong to class II, among which genotype VII is the predominant genotype in China. However, the molecular evolution and phylodynamics of class II genotype VII NDV strains in China remained largely unknown. In this study, we identified 13 virulent NDV including 11 genotype VII strains and 2 genotype IX strains, from clinical samples during 1997 to 2019. Combined NDV sequences submitted to GenBank, we investigate evolution, and transmission dynamics of class II NDVs in China, especially genotype VII strains. Our results revealed that East and South China have the most genotypic diversity of class II NDV, and East China might be the origin of genotype VII NDVs in China. In addition, genotype VII NDVs in China are presumably transmitted by chickens, as the virus was most prevalent in chickens. Furthermore, codon usage analysis revealed that the F genes of genotype VII NDVs have stronger adaptation in chickens, and six amino acids in this gene are found under positive selection via selection model analysis. Collectively, our results revealed the genetic diversity and evolutionary dynamics of genotype VII NDVs in China, providing important insights into the epidemiology of these viruses in China.

Analysis of the Codon Usage Pattern of HA and NA Genes of H7N9 Influenza A Virus

Novel H7N9 influenza virus transmitted from birds to human and, since March 2013, it has caused five epidemic waves in China. Although the evolution of H7N9 viruses has been investigated, the evolutionary changes associated with codon usage are still unclear. Herein, the codon usage pattern of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA), was studied to understand the evolutionary changes in relation to host, epidemic wave, and pathogenicity. Both genes displayed a low codon usage bias, with HA higher than NA. The codon usage was driven by mutation pressure and natural selection, although the main contributing factor was natural selection. Additionally, the codon adaptation index (CAI) and deoptimization (RCDI) illustrated the strong adaptability of H7N9 to Gallus gallus. Similarity index (SiD) analysis showed that Homo sapiens posed a stronger selection pressure than Gallus gallus. Thus, we assume that this may be related to the gradual adaptability of the virus to human. In addition, the host strong selection pressure was validated based on CpG dinucleotide content. In conclusion, this study analyzed the usage of codons of two genes of H7N9 and expanded our understanding of H7N9 host specificity. This aids into the development of control measures against H7N9 influenza virus.

Sequence analysis of SARS-CoV-2 genome reveals features important for vaccine design

As the SARS-CoV-2 pandemic is rapidly progressing, the need for the development of an effective vaccine is critical. A promising approach for vaccine development is to generate, through codon pair deoptimization, an attenuated virus. This approach carries the advantage that it only requires limited knowledge specific to the virus in question, other than its genome sequence. Therefore, it is well suited for emerging viruses, for which we may not have extensive data. We performed comprehensive in silico analyses of several features of SARS-CoV-2 genomic sequence (e.g., codon usage, codon pair usage, dinucleotide/junction dinucleotide usage, RNA structure around the frameshift region) in comparison with other members of the coronaviridae family of viruses, the overall human genome, and the transcriptome of specific human tissues such as lung, which are primarily targeted by the virus. Our analysis identified the spike (S) and nucleocapsid (N) proteins as promising targets for deoptimization and suggests a roadmap for SARS-CoV-2 vaccine development, which can be generalizable to other viruses.

Codon usage pattern and evolutionary forces of mitochondrial ND genes among orders of class Amphibia

In this study, we used a bioinformatics approach to analyze the nucleotide composition and pattern of synonymous codon usage in mitochondrial ND genes in three amphibian groups, that is, orders Anura, Caudata, and Gymnophiona to identify the commonality and the differences of codon usage as no research work was reported yet. The high value of the effective number of codons revealed that the codon usage bias (CUB) was low in mitochondrial ND genes among the orders. Nucleotide composition analysis suggested that for each gene, the compositional features differed among Anura, Caudata, and Gymnophiona and the GC content was lower than AT content. Furthermore, a highly significant difference (p < .05) for GC content was found in each gene among the orders. The heat map showed contrasting patterns of codon usage among different ND genes. The regression of GC12 on GC3 suggested a narrow range of GC3 distribution and some points were located in the diagonal, indicating both mutation pressure and natural selection might influence the CUB. Moreover, the slope of the regression line was less than 0.5 in all ND genes among orders, indicating natural selection might have played the dominant role whereas mutation pressure had played a minor role in shaping CUB of ND genes across orders.

Update of the current knowledge on genetics, evolution, immunopathogenesis, and transmission for coronavirus disease 19 (COVID-19)

In December 2019, an acute respiratory disease caused by novel species of coronavirus (SARS-CoV-2), emerged in China and has spread throughout the world. On 11th March 2020, the World Health Organization (WHO) officially declared coronavirus disease 19 (COVID-19) a pandemic, severe coronavirus-mediated human disease. Based on genomic and phylogenetic studies, SARS-CoV-2 might originate from bat coronaviruses and infects humans directly or through intermediate zoonotic hosts. However, the exact origin or the host intermediate remains unknown. Genetically, SARS-CoV-2 is similar to several existing coronaviruses, particularly SARS-CoV, but differs by silent and non-silent mutations. The virus uses different transmission routes and targets cells and tissues with angiotensin-converting enzyme 2 (ACE2) protein, which makes it contagious. COVID-19 shares both the main clinical features and excessive/dysregulated cell responses with the two previous Middle East respiratory syndrome coronavirus (MERS) and severe acute respiratory syndrome coronavirus (SARS) epidemics. In this review, we provide an update of the current knowledge on the COVID-19 pandemic. Gaining a deeper understanding of SARS-CoV-2 structure, transmission routes, and molecular responses, will assist in the prevention and control of COVID-19 outbreaks in the future.

Analysis of compositional properties and codon usage bias of mitochondrial CYB gene in anura, urodela and gymnophiona

We delineated the pattern of synonymous codon usage bias (CUB) and its determinants in mitochondrial CYB gene of respiratory chain across different amphibian groups namely orders anura, urodela and gymnophiona. We observed that CUB was low in CYB gene of amphibia. The gymnophionans had comparatively high bias followed by urodeles and anurans. The codons namely TCA, CCA, CAA, CGA, TGA, AAA and ACA were over-represented in all three orders. The codons such as GCC and TCC were over-represented in anura whereas in urodela, the over-represented codons were TTA, CTA, ATA, GTA, GAA, GGA and GCA. In gymnophiona, GCC, TTA, CTA, ATA, GTA, GAA and GGA codons were over-represented. The regression analysis between effective number of codons (ENC) and nucleobase at the 3rd position revealed that nucleobase A and C influenced CUB positively in order anura, while in urodela and gymnophiona, nucleobase A and T influenced the CUB positively. Mutation pressure and natural selection mutually illustrate the CUB of CYB gene (complex III gene) of amphibia as elucidated by correlation analysis between 3rd nucleotide in a codon and overall nucleotide content of the gene. However, neutrality plot showed that natural selection was the dominant evolutionary factor of CUB.

Codon usage similarity between viral and some host genes suggests a codon-specific translational regulation

The codon usage pattern is a specific characteristic of each species; however, the codon usage of all of the genes in a genome is not uniform. Intriguingly, most viruses have codon usage patterns that are vastly different from the optimal codon usage of their hosts. How viral genes with different codon usage patterns are efficiently expressed during a viral infection is unclear. An analysis of the similarity between viral codon usage and the codon usage of the individual genes of a host genome has never been performed. In this study, we demonstrated that the codon usage of human RNA viruses is similar to that of some human genes, especially those involved in the cell cycle. This finding was substantiated by its concordance with previous reports of an upregulation at the protein level of some of these biological processes. It therefore suggests that some suboptimal viral codon usage patterns may actually be compatible with cellular translational machineries in infected conditions.

From SARS and MERS CoVs to SARS-CoV-2: Moving toward more biased codon usage in viral structural and nonstructural genes

Background

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) is an emerging disease with fatal outcomes. In this study, a fundamental knowledge gap question is to be resolved by evaluating the differences in biological and pathogenic aspects of SARS‐CoV‐2 and the changes in SARS‐CoV‐2 in comparison with the two prior major COV epidemics, SARS and Middle East respiratory syndrome (MERS) coronaviruses.

Methods

The genome composition, nucleotide analysis, codon usage indices, relative synonymous codons usage, and effective number of codons (ENc) were analyzed in the four structural genes; Spike (S), Envelope (E), membrane (M), and Nucleocapsid (N) genes, and two of the most important nonstructural genes comprising RNA‐dependent RNA polymerase and main protease (Mpro) of SARS‐CoV‐2, Beta‐CoV from pangolins, bat SARS, MERS, and SARS CoVs.

Results

SARS‐CoV‐2 prefers pyrimidine rich codons to purines. Most high‐frequency codons were ending with A or T, while the low frequency and rare codons were ending with G or C. SARS‐CoV‐2 structural proteins showed 5 to 20 lower ENc values, compared with SARS, bat SARS, and MERS CoVs. This implies higher codon bias and higher gene expression efficiency of SARS‐CoV‐2 structural proteins. SARS‐CoV‐2 encoded the highest number of over‐biased and negatively biased codons. Pangolin Beta‐CoV showed little differences with SARS‐CoV‐2 ENc values, compared with SARS, bat SARS, and MERS CoV.

Conclusion

Extreme bias and lower ENc values of SARS‐CoV‐2, especially in Spike, Envelope, and Mpro genes, are suggestive for higher gene expression efficiency, compared with SARS, bat SARS, and MERS CoVs.

Analysis of preferred codon usage in the coronavirus N genes and their implications for genome evolution and vaccine design

•

The nucleotide variations among the N genes of 13 different coronaviruses (CoVs) were interpreted.

•

Overall, 18 amino acids observed with varying preferred codons.

•

The effective number of codon values ranged from 40.43 to 53.85, revealing a slight codon bias.

•

A highly significant correlation between GC3s and ENc values was observed in porcine epidemic diarrhea CoV, followed by Middle East respiratory syndrome CoV.

The nucleocapsid (N) protein of a coronavirus plays a crucial role in virus assembly and in its RNA transcription. It is important to characterize a virus at the nucleotide level to discover the virus’s genomic sequence variations and similarities relative to other viruses that could have an impact on the functions of its genes and proteins. This entails a comprehensive and comparative analysis of the viral genomes of interest for preferred nucleotides, codon bias, nucleotide changes at the 3^rd position (NT3s), synonymous codon usage and relative synonymous codon usage. In this study, the variations in the N proteins among 13 different coronaviruses (CoVs) were analysed at the nucleotide and amino acid levels in an attempt to reveal how these viruses adapt to their hosts relative to their preferred codon usage in the N genes. The results revealed that, overall, eighteen amino acids had different preferred codons and eight of these were over-biased. The N genes had a higher AT% over GC% and the values of their effective number of codons ranged from 40.43 to 53.85, indicating a slight codon bias. Neutrality plots and correlation analyses showed a very high level of GC3s/GC correlation in porcine epidemic diarrhea CoV (pedCoV), followed by Middle East respiratory syndrome-CoV (MERS CoV), porcine delta CoV (dCoV), bat CoV (bCoV) and feline CoV (fCoV) with r values 0.81, 0.68, -0.47, 0.98 and 0.58, respectively. These data implied a high rate of evolution of the CoV genomes and a strong influence of mutation on evolutionary selection in the CoV N genes. This type of genetic analysis would be useful for evaluating a virus’s host adaptation, evolution and is thus of value to vaccine design strategies.

Multivariate analyses of codon usage of SARS-CoV-2 and other betacoronaviruses

Coronavirus disease 2019 (COVID-19) is a global health concern as it continues to spread within China and beyond. The causative agent of this disease, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), belongs to the genus Betacoronavirus, which also includes severe acute respiratory syndrome-related coronavirus (SARSr-CoV) and Middle East respiratory syndrome-related coronavirus (MERSr-CoV). Codon usage of viral genes are believed to be subjected to different selection pressures in different host environments. Previous studies on codon usage of influenza A viruses helped identify viral host origins and evolution trends, however, similar studies on coronaviruses are lacking. In this study, we compared the codon usage bias using global correspondence analysis (CA), within-group CA and between-group CA. We found that the bat RaTG13 virus best matched the overall codon usage pattern of SARS-CoV-2 in orf1ab, spike and nucleocapsid genes, while the pangolin P1E virus had a more similar codon usage in membrane gene. The amino acid usage pattern of SARS-CoV-2 was generally found similar to bat and human SARSr-CoVs. However, we found greater synonymous codon usage differences between SARS-CoV-2 and its phylogenetic relatives on spike and membrane genes, suggesting these two genes of SARS-CoV-2 are subjected to different evolutionary pressures.

Multivariate analyses of codon usage of SARS-CoV-2 and other betacoronaviruses

Coronavirus disease 2019 (COVID-19) is a global health concern as it continues to spread within China and beyond. The causative agent of this disease, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), belongs to the genus Betacoronavirus, which also includes severe acute respiratory syndrome-related coronavirus (SARSr-CoV) and Middle East respiratory syndrome-related coronavirus (MERSr-CoV). Codon usage of viral genes are believed to be subjected to different selection pressures in different host environments. Previous studies on codon usage of influenza A viruses helped identify viral host origins and evolution trends, however, similar studies on coronaviruses are lacking. In this study, we compared the codon usage bias using global correspondence analysis (CA), within-group CA and between-group CA. We found that the bat RaTG13 virus best matched the overall codon usage pattern of SARS-CoV-2 in orf1ab, spike and nucleocapsid genes, while the pangolin P1E virus had a more similar codon usage in membrane gene. The amino acid usage pattern of SARS-CoV-2 was generally found similar to bat and human SARSr-CoVs. However, we found greater synonymous codon usage differences between SARS-CoV-2 and its phylogenetic relatives on spike and membrane genes, suggesting these two genes of SARS-CoV-2 are subjected to different evolutionary pressures.

The dynamic proteome of influenza A virus infection identifies M segment splicing as a host range determinant

Pandemic influenza A virus (IAV) outbreaks occur when strains from animal reservoirs acquire the ability to infect and spread among humans. The molecular basis of this species barrier is incompletely understood. Here we combine metabolic pulse labeling and quantitative proteomics to monitor protein synthesis upon infection of human cells with a human- and a bird-adapted IAV strain and observe striking differences in viral protein synthesis. Most importantly, the matrix protein M1 is inefficiently produced by the bird-adapted strain. We show that impaired production of M1 from bird-adapted strains is caused by increased splicing of the M segment RNA to alternative isoforms. Strain-specific M segment splicing is controlled by the 3' splice site and functionally important for permissive infection. In silico and biochemical evidence shows that avian-adapted M segments have evolved different conserved RNA structure features than human-adapted sequences. Thus, we identify M segment RNA splicing as a viral host range determinant.

Comprehensive analysis of synonymous codon usage patterns in orf3 gene of porcine epidemic diarrhea virus in China

The ORF3 protein of porcine epidemic diarrhea virus (PEDV) is found to function as an ion channel which influences virus virulence and production. Taking consideration of the importance of PEDV orf3 gene, we have performed comprehensive analysis to investigate its synonymous codon usage patterns. In this study, the results of base composition analysis showed A/T rich and G/C poor in PEDV orf3 genes, and the most abundant base was nucleotide T. The relative synonymous codon usage value in each codon revealed that codon usage bias existed. The mean ENC value of each gene was 48.75, indicating a low codon usage bias, as well as a relatively instable change in PEDV orf3 genes. The general correlation analysis between base composition and codon usage bias indicated that mutational bias has an impact on the PEDV codon usage bias. Neutral analysis suggested that natural selection pressure takes a more important influence than mutational bias in shaping codon usage bias. Moreover, other factors including hydrophobicity and aromaticity have been also found to influence the codon usage variation among the PEDV orf3 genes. This study not only represents the most systematic analysis of codon usage patterns in PEDV orf3 genes, but also provides a basic shaping mechanism of the codon usage bias.

•

ORF3 genes are A/T rich and C/G poor, and codon usage bias existed.

•

Natural selection pressure plays a more important role than mutational bias in codon usage bias.

•

Hydrophobicity and aromaticity also influence the codon usage variation.

Insights into the nucleotide composition and codon usage pattern of human tumor suppressor genes

Tumor suppressor genes encode different proteins that inhibit the uncontrolled proliferation of cell growth and tumor development. To acquire clues for predicting gene expression level, it is essential to understand the codon usage bias (CUB) of genes to characterize genome which possesses its own compositional characteristics and unique coding sequences. We used bioinformatic tools to analyze the codon usage patterns of 637 human tumor suppressor genes as no work was reported earlier. The mean effective number of codons of these genes was 48, indicating low CUB. Our results exhibited a significant positive correlation among different nucleotide compositions and the codons ending with C base was most frequently used along with the most over-represented codon CTG and GTG codifying leucine and valine amino acid, respectively, in human tumor suppressor genes. The neutrality plot showed a significant positive correlation (Pearson, r = 0. 646; P < .01) suggesting that mutation on GC bias might affect the CUB. However, the linear regression coefficient of GC₁₂ on GC₃ in human tumor suppressor genes suggested that natural selection played a major role while mutation pressure played a minor role in the codon usage patterns of tumor suppressor genes in human. Our study would throw light into the factors that affect CUB and the codon usage patterns in the human tumor suppressor genes.

Analysis of Synonymous Codon Usage Bias in Flaviviridae Virus

Background

Flaviviridae viruses are single-stranded, positive-sense RNA viruses, which threat human constantly mediated by mosquitoes, ticks, and sandflies. Considering the recent increase in the prevalence of the family virus and its risk potential, we investigated the codon usage pattern to understand its evolutionary processes and provide some useful data to develop the medications for most of Flaviviridae viruses.

Results

The overall extent of codon usage bias in 65 Flaviviridae viruses is low with the average value of GC contents being 50.5% and the highest value being 55.9%; the lowest value is 40.2%. ENC values of Flaviviridae virus genes vary from 48.75 to 57.83 with a mean value of 55.56. U- and A-ended codons are preferred in the Flaviviridae virus. Correlation analysis shows that the positive correlation between ENC value and GC content at the third nucleotide positions was significant in this family virus. The result of analysis of ENC, neutrality plot analysis, and correlation analysis revealed that codon usage bias of all the viruses was affected mainly by natural selection. Meanwhile, according to correspondence analysis (CoA) based on RSCU and phylogenetic analysis, the Flaviviridae viruses mainly are made up of two groups, Group I (Yellow fever virus, Apoi virus, Tembusu virus, Dengue virus 1, and others) and Group II (West Nile virus lineage 2, Japanese encephalitis virus, Usutu virus, Kedougou virus, and others).

Conclusions

All in, the bias of codon usage pattern is affected not only by compositional constraints but also by natural selection. Phylogenetic analysis also illustrates that codon usage bias of virus can serve as an effective means of evolutionary classification in Flaviviridae virus.

Compositional properties and codon usage pattern of mitochondrial ATP gene in different classes of Arthropoda

Codon usage bias (CUB) is defined as the usage of synonymous codons unequally for an amino acid in a gene transcript. It is influenced by both mutation pressure and natural selection and is a species-specific property. In our current study, we used bioinformatic methods to investigate the coding sequences of mitochondrial adenosine triphosphate gene (MT-ATP) in different classes of arthropoda to know the codon usage pattern of the gene as no work was described earlier. The analysis of compositional properties suggested that the gene is AT rich. The effective number of codons revealed the CUB of both ATP6 and ATP8 gene was moderate. Heat map showed that the codons ending with AT were negatively associated with GC3 while the codons ending with GC were positively associated with GC3 in all the classes of arthropoda. Correspondence study revealed that the pattern of codon usage of ATP6 and ATP8 genes differed across classes. Neutrality plot suggested the codon usage bias of these two genes in phylum arthropoda was influenced by both mutation pressure and natural selection.

Novel Approaches for The Development of Live Attenuated Influenza Vaccines

Influenza virus still represents a considerable threat to global public health, despite the advances in the development and wide use of influenza vaccines. Vaccination with traditional inactivate influenza vaccines (IIV) or live-attenuated influenza vaccines (LAIV) remains the main strategy in the control of annual seasonal epidemics, but it does not offer protection against new influenza viruses with pandemic potential, those that have shifted. Moreover, the continual antigenic drift of seasonal circulating influenza viruses, causing an antigenic mismatch that requires yearly reformulation of seasonal influenza vaccines, seriously compromises vaccine efficacy. Therefore, the quick optimization of vaccine production for seasonal influenza and the development of new vaccine approaches for pandemic viruses is still a challenge for the prevention of influenza infections. Moreover, recent reports have questioned the effectiveness of the current LAIV because of limited protection, mainly against the influenza A virus (IAV) component of the vaccine. Although the reasons for the poor protection efficacy of the LAIV have not yet been elucidated, researchers are encouraged to develop new vaccination approaches that overcome the limitations that are associated with the current LAIV. The discovery and implementation of plasmid-based reverse genetics has been a key advance in the rapid generation of recombinant attenuated influenza viruses that can be used for the development of new and most effective LAIV. In this review, we provide an update regarding the progress that has been made during the last five years in the development of new LAIV and the innovative ways that are being explored as alternatives to the currently licensed LAIV. The safety, immunogenicity, and protection efficacy profile of these new LAIVs reveal their possible implementation in combating influenza infections. However, efforts by vaccine companies and government agencies will be needed for controlled testing and approving, respectively, these new vaccine methodologies for the control of influenza infections.

Genetic Evolution and Molecular Selection of the HE Gene of Influenza C Virus

Influenza C virus (ICV) was first identified in humans and swine, but recently also in cattle, indicating a wider host range and potential threat to both the livestock industry and public health than was originally anticipated. The ICV hemagglutinin-esterase (HE) glycoprotein has multiple functions in the viral replication cycle and is the major determinant of antigenicity. Here, we developed a comparative approach integrating genetics, molecular selection analysis, and structural biology to identify the codon usage and adaptive evolution of ICV. We show that ICV can be classified into six lineages, consistent with previous studies. The HE gene has a low codon usage bias, which may facilitate ICV replication by reducing competition during evolution. Natural selection, dinucleotide composition, and mutation pressure shape the codon usage patterns of the ICV HE gene, with natural selection being the most important factor. Codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analysis revealed that the greatest adaption of ICV was to humans, followed by cattle and swine. Additionally, similarity index (SiD) analysis revealed that swine exerted a stronger evolutionary pressure on ICV than humans, which is considered the primary reservoir. Furthermore, a similar tendency was also observed in the M gene. Of note, we found HE residues 176, 194, and 198 to be under positive selection, which may be the result of escape from antibody responses. Our study provides useful information on the genetic evolution of ICV from a new perspective that can help devise prevention and control strategies.

Analysis of codon usage patterns and influencing factors in Nipah virus

•

Nipah virus (NiV) genes are AT-rich.

•

Codon usage bias of NiV genes is low.

•

Patterns of codon usage bias differ across the genomes of NiV.

•

Both mutation pressure and natural selection influenced codon usage bias of NiV genes.

Codon usage bias (CUB) is the unequal usage of synonymous codons of an amino acid in which some codons are used more often than others and is widely used in understanding molecular biology, genetics, and functional regulation of gene expression. Nipah virus (NiV) is an emerging zoonotic paramyxovirus that causes fatal disease in both humans and animals. NiV was first identified during an outbreak of a disease in Malaysia in 1998 and then occurred periodically since 2001 in India, Bangladesh, and the Philippines. We used bioinformatics tools to analyze the codon usage patterns in a genome-wide manner among 11 genomes of NiV as no work was reported yet. The compositional properties revealed that the overall GC and AT contents were 41.96 and 58.04%, respectively i.e. Nipah virus genes were AT-rich. Correlation analysis between overall nucleotide composition and its 3^rd codon position suggested that both mutation pressure and natural selection might influence the CUB across Nipah genomes. Neutrality plot revealed natural selection might have played a major role while mutation pressure had a minor role in shaping the codon usage bias in NiV genomes.

A cross-talk on compositional dynamics and codon usage patterns of mitochondrial CYB gene in Echinodermata

Codon usage bias (CUB) refers to a phenomenon in which some synonymous codons are used in mature mRNA at a higher frequency than other members codifying the same amino acid. CUB is mainly determined by mutation pressure and natural selection. We used bioinformatic tools to analyze the protein coding sequences of mitochondrial CYB gene in different classes of Echinodermata to understand the patterns of codon usage. The ENC values of CYB gene in five different classes of Echinodermata were 41.64, 30.33, 43.63, 41.11, and 41.33, which suggested that the CUB of this gene was low. The relative synonymous codon usage (RSCU) values showed that the patterns of over-represented and under-represented codons were different among different classes. Correspondence analysis indicated that the plots of CYB gene were different across classes, suggesting that the pattern of codon usage was also different among five classes under study. Highly significant correlation (p < .01) between overall nucleotide composition and its 3rd codon position indicated that both mutational pressure and natural selection had an influence on the codon usage bias of CYB gene. Furthermore, PR-2 bias plot analysis showed that both mutation pressure and natural selection might have affected the pattern of codon usage in CYB gene of Echinodermata.

DNA compositional dynamics and codon usage patterns of M1 and M2 matrix protein genes in influenza A virus

Influenza A virus subtype H3N2 has been a serious health issue across the globe with approximately 36 thousand annual casualties in the United States of America only. Co-circulation in avian, swine and human hosts has led to frequent mutations in the virus genome, due to which development of successful antivirals against the virus has become a formidable challenge. Recently, focussed research is being carried out targeting the matrix proteins of this strain as vaccine candidates. This study is carried out to unravel the key features of the genes encoding the matrix proteins that manoeuvre the codon usage profile in the H3N2 strains. The findings reveal differential codon choice for both matrix protein 1 and matrix protein 2. The overall codon usage bias is less pronounced in both the datasets which is evident from higher value of effective number of codons (>55). Comparison of the codon usage for both the genes under study with that of humans revealed that the viral codon usage is not fully optimized for the human host conditions. Both the genes enrolled in the study showed variation which was reflected in almost all the indices used for codon usage studies. Neutrality analysis revealed a weak role of mutation pressure while selection was the major contributor towards codon usage.

Gene Birth Contributes to Structural Disorder Encoded by Overlapping Genes

The same nucleotide sequence can encode two protein products in different reading frames. Overlapping gene regions encode higher levels of intrinsic structural disorder (ISD) than nonoverlapping genes (39% vs. 25% in our viral dataset). This might be because of the intrinsic properties of the genetic code, because one member per pair was recently born de novo in a process that favors high ISD, or because high ISD relieves increased evolutionary constraint imposed by dual-coding. Here, we quantify the relative contributions of these three alternative hypotheses. We estimate that the recency of de novo gene birth explains [Formula: see text] or more of the elevation in ISD in overlapping regions of viral genes. While the two reading frames within a same-strand overlapping gene pair have markedly different ISD tendencies that must be controlled for, their effects cancel out to make no net contribution to ISD. The remaining elevation of ISD in the older members of overlapping gene pairs, presumed due to the need to alleviate evolutionary constraint, was already present prior to the origin of the overlap. Same-strand overlapping gene birth events can occur in two different frames, favoring high ISD either in the ancestral gene or in the novel gene; surprisingly, most de novo gene birth events contained completely within the body of an ancestral gene favor high ISD in the ancestral gene (23 phylogenetically independent events vs. 1). This can be explained by mutation bias favoring the frame with more start codons and fewer stop codons.

Comprehensive Analysis of Codon Usage on Rabies Virus and Other Lyssaviruses

Rabies virus (RABV) and other lyssaviruses can cause rabies and rabies-like diseases, which are a persistent public health threat to humans and other mammals. Lyssaviruses exhibit distinct characteristics in terms of geographical distribution and host specificity, indicative of a long-standing diversification to adapt to the environment. However, the evolutionary diversity of lyssaviruses, in terms of codon usage, is still unclear. We found that RABV has the lowest codon usage bias among lyssaviruses strains, evidenced by its high mean effective number of codons (ENC) (53.84 ± 0.35). Moreover, natural selection is the driving force in shaping the codon usage pattern of these strains. In summary, our study sheds light on the codon usage patterns of lyssaviruses, which can aid in the development of control strategies and experimental research.

Codon Usage Pattern of Genes Involved in Central Nervous System

Codon usage bias (CUB) is the non-uniform usage of synonymous codons in which some codons are more preferred to others in the transcript. Analysis of codon usage bias has applications in understanding the basics of molecular biology, genetics, gene expression, and molecular evolution. To understand the patterns of codon usage in genes involved in the central nervous system (CNS), we used bioinformatic approaches to analyze the protein-coding sequences of genes involved in the CNS. The improved effective number of codons (ENC) suggested that the overall codon usage bias was low. The relative synonymous codon usage (RSCU) revealed that the most frequently occurring codons had a G or C at the third codon position. The codons namely TCC, AGC, CTG, CAG, CGC, ATC, ACC, GTG, GCC, GGC, and CGG (average RSCU > 1.6) were over-represented. Both mutation pressure and natural selection might affect the codon usage pattern as evident from correspondence and parity plot analyses. The overall GC content (59.93) was higher than AT content, i.e., genes were GC-rich. The correlation of GC12 with GC3 suggested that mutation pressure might affect the codon usage pattern.

Genetic and evolutionary analysis of emerging H3N2 canine influenza virus

The H3N2 canine influenza virus (CIV) originated from an avian species. Since its emergence, it has circulated in multiple states and has caused pandemics among dog populations; however, no comprehensive studies have explored the causes driving these ongoing cases. The study of the codon usage patterns of viruses can reveal the genetic changes required for the viruses to adapt to new hosts and the external environment. Here we performed a thorough genetic, evolutionary, and codon usage analysis. We identified three evolutionary H3N2 CIV clades from a timescaled phylogenetic tree, namely, Origin, China, and Korea/USA, by principal component analysis (PCA). Additionally, we found a low codon usage bias and that mutation pressure, natural selection, and dinucleotide abundance shape the codon usage bias of H3N2 CIVs, with natural selection being more crucial than the others. Moreover, the human codon adaptation index was similar to that of dogs (the natural host) and cats. In addition, the H3N2 CIV similarity index values were higher than those of the avian influenza virus (AIV), suggesting viral adaptation to the host. Therefore, H3N2 CIVs may pose a potential risk to public health in the future, and further epidemiologic, evolutionary, and pathogenetic studies are required.

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.

Genetic and codon usage bias analyses of polymerase genes of equine influenza virus and its relation to evolution

BACKGROUND

Equine influenza is a major health problem of equines worldwide. The polymerase genes of influenza virus have key roles in virus replication, transcription, transmission between hosts and pathogenesis. Hence, the comprehensive genetic and codon usage bias of polymerase genes of equine influenza virus (EIV) were analyzed to elucidate the genetic and evolutionary relationships in a novel perspective.

RESULTS

The group - specific consensus amino acid substitutions were identified in all polymerase genes of EIVs that led to divergence of EIVs into various clades. The consistent amino acid changes were also detected in the Florida clade 2 EIVs circulating in Europe and Asia since 2007. To study the codon usage patterns, a total of 281,324 codons of polymerase genes of EIV H3N8 isolates from 1963 to 2015 were systemically analyzed. The polymerase genes of EIVs exhibit a weak codon usage bias. The ENc-GC3s and Neutrality plots indicated that natural selection is the major influencing factor of codon usage bias, and that the impact of mutation pressure is comparatively minor. The methods for estimating host imposed translation pressure suggested that the polymerase acidic (PA) gene seems to be under less translational pressure compared to polymerase basic 1 (PB1) and polymerase basic 2 (PB2) genes. The multivariate statistical analysis of polymerase genes divided EIVs into four evolutionary diverged clusters - Pre-divergent, Eurasian, Florida sub-lineage 1 and 2.

CONCLUSIONS

Various lineage specific amino acid substitutions observed in all polymerase genes of EIVs and especially, clade 2 EIVs underwent major variations which led to the emergence of a phylogenetically distinct group of EIVs originating from Richmond/1/07. The codon usage bias was low in all the polymerase genes of EIVs that was influenced by the multiple factors such as the nucleotide compositions, mutation pressure, aromaticity and hydropathicity. However, natural selection was the major influencing factor in defining the codon usage patterns and evolution of polymerase genes of EIVs.

Codon usage analysis of photolyase encoding genes of cyanobacteria inhabiting diverse habitats

Nucleotide and amino acid compositions were studied to determine the genomic and structural relationship of photolyase gene in freshwater, marine and hot spring cyanobacteria. Among three habitats, photolyase encoding genes from hot spring cyanobacteria were found to have highest GC content. The genomic GC content was found to influence the codon usage and amino acid variability in photolyases. The third position of codon was found to have more effect on amino acid variability in photolyases than the first and second positions of codon. The variation of amino acids Ala, Asp, Glu, Gly, His, Leu, Pro, Gln, Arg and Val in photolyases of three different habitats was found to be controlled by first position of codon (G1C1). However, second position (G2C2) of codon regulates variation of Ala, Cys, Gly, Pro, Arg, Ser, Thr and Tyr contents in photolyases. Third position (G3C3) of codon controls incorporation of amino acids such as Ala, Phe, Gly, Leu, Gln, Pro, Arg, Ser, Thr and Tyr in photolyases from three habitats. Photolyase encoding genes of hot spring cyanobacteria have 85% codons with G or C at third position, whereas marine and freshwater cyanobacteria showed 82 and 60% codons, respectively, with G or C at third position. Principal component analysis (PCA) showed that GC content has a profound effect in separating the genes along the first major axis according to their RSCU (relative synonymous codon usage) values, and neutrality analysis indicated that mutational pressure has resulted in codon bias in photolyase genes of cyanobacteria.