Optimization of codon usage of the envelope protein E2 gene from various genotypes of hepatitis C virus to predict the expression level in <Emphasis Type="Italic">Pichia pastoris</Emphasis>

Hepatitis C virus infection (HCV) alarmingly increases worldwide; it causes chronic hepatitis, liver cirrhosis and hepatocellular carcinoma, so there is urgent need of developing effective and sufficient quantity of vaccine. HCV envelope protein E2 is the main target for developing as a vaccine candidate. Presently recombinant proteins can successfully be used as a vaccine for many diseases. This concern, it is challenging to produce sufficient quantities of many recombinant proteins from their expression hosts. One of the main factors affecting the success of expression of foreign genes in heterologous hosts is the divergence of codon usage of the target gene from that used in the expression system. In this study, we optimized the various genotypes of HCV envelope protein E2 gene according to the codon usage of Pichia pastoris and predicted the expression level. Synonymous codon usage of E2 adapted to that used by P. pastoris was estimated using the relative synonymous codon usage value (RSCU), codon adaptation index (CAI) and effective number of codon (ENC). The CAI of optimized HCV E2 sequences was enhanced from 0.638 to 0.833 and %GC was decreased from 56.05 to 44.05; this was significantly (p < 0.01) different from the native sequences. Codon with RSCU value less than one was replaced with most preferred synonymous codons. The ENC values of optimized HCV E2 sequences varied from 47.00 to 47.50, with a mean value of 47.15 and an SD of 0.14. Our study suggested that, from the measured values of predicted expression level, the codon optimized HCV E2 protein could be produced in sufficient quantity in the expression host; knowledge of the codon usage patterns of E2 of various genotypes facilitate the production of a promising unique vaccine candidate for HCV.     

Revelation of Influencing Factors in Overall Codon Usage Bias of Equine Influenza Viruses

Equine influenza viruses (EIVs) of H3N8 subtype are culprits of severe acute respiratory infections in horses, and are still responsible for significant outbreaks worldwide. Adaptability of influenza viruses to a particular host is significantly influenced by their codon usage preference, due to an absolute dependence on the host cellular machinery for their replication. In the present study, we analyzed genome-wide codon usage patterns in 92 EIV strains, including both H3N8 and H7N7 subtypes by computing several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis disclosed bias of preferred synonymous codons towards A/U-ended codons. The overall codon usage bias in EIVs was slightly lower, and mainly affected by the nucleotide compositional constraints as inferred from the RSCU and effective number of codon (ENc) analysis. Our data suggested that codon usage pattern in EIVs is governed by the interplay of mutation pressure, natural selection from its hosts and undefined factors. The H7N7 subtype was found less fit to its host (horse) in comparison to H3N8, by possessing higher codon bias, lower mutation pressure and much less adaptation to tRNA pool of equine cells. To the best of our knowledge, this is the first report describing the codon usage analysis of the complete genomes of EIVs. The outcome of our study is likely to enhance our understanding of factors involved in viral adaptation, evolution, and fitness towards their hosts.     

Synonymous Codon Usage in TTSuV2: Analysis and Comparison with TTSuV1

Two species of the DNA virus Torque teno sus virus (TTSuV), TTSuV1 and TTSuV2, have become widely distributed in pig-farming countries in recent years. In this study, we performed a comprehensive analysis of synonymous codon usage bias in 41 available TTSuV2 coding sequences (CDS), and compared the codon usage patterns of TTSuV2 and TTSuV1. TTSuV codon usage patterns were found to be phylogenetically conserved. Values for the effective number of codons (ENC) indicated that the overall extent of codon usage bias in both TTSuV2 and TTSuV1 was not significant, the most frequently occurring codons had an A or C at the third codon position. Correspondence analysis (COA) was performed and TTSuV2 and TTSuV1 sequences were located in different quadrants of the first two major axes. A plot of the ENC revealed that compositional constraint was the major factor determining the codon usage bias for TTSuV2. In addition, hierarchical cluster analysis of 41 TTSuV2 isolates based on relative synonymous codon usage (RSCU) values suggested that there was no association between geographic distribution and codon bias of TTSuV2 sequences. Finally, the comparison of RSCU for TTSuV2, TTSuV1 and the corresponding host sequence indicated that the codon usage pattern of TTSuV2 was similar to that of TTSuV1. However the similarity was low for each virus and its host. These conclusions provide important insight into the synonymous codon usage pattern of TTSuV2, as well as better understangding of the molecular evolution of TTSuV2 genomes.     

Sequence variation of HIV and bioinformatics

The envelope glycoprotein of human immunodeficiency virus type 1 (HIV-1) interacts with receptors on the target cell and mediates virus entry by fusing the viral and cell membranes. To maintain the viral infectivity, amino acids that interact with receptors are expected to be more conserved than the other sites on the protein surface. In contrast to the functional constraint of amino acids for the receptor binding, some amino acid changes in this protein may produce antigenic variations that enable the virus to escape from recognition of the host immune system. Therefore, both positive selection (higher fitness) and negative selection (lower fitness) against amino acid changes are taking place during evolution of surface proteins of parasites To elucidate the evolutionary mechanisms of the whole HIV-1 gp120 envelope glycoprotein at the single site level, we collected and analyzed all available sequence data for the protein. By analyzing 186 sequences of the HIV-1 gp120 (subtype B), we reevaluated amino acid variability at the single site level, and estimated the numbers of synonymous and nonsynonymous substitutions at each codon position to detect positive and negative selection. We identified 33 amino acid positions which may be under positive selection. Some of these positions may form discontinuous epitopes. We also analyzed amino acid sequences to find amino acid positions responsible for usage of the second receptor. We found that, in addition to the V3 loop, amino acid variation at residue 440 in C4 region is clearly linked with the usage of CXCR 4.     

Analysis of synonymous codon usage bias in <Emphasis Type="Italic">helicase</Emphasis> gene from <Emphasis Type="Italic">Autographa californica</Emphasis> multiple nucleopolyhedrovirus

The helicase gene of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is not only involved in viral DNA replication, but also plays a role in viral host range. To identify the codon usage bias of helicase of AcMNPV, the codon usage bias of helicase was especially studies in AcMNPV and 41 reference strains of baculoviruses by calculating the codon adaptation index (CAI), effective number of codon (ENc), relative synonymous codon usage (RSCU), and other indices. The helicase of baculovirus is less biased (mean ENc?=?50.539?>?40; mean CAI?=?0.246). AcMNPV helicase has a strong bias toward the synonymous codons with G and C at the third codon position (GC3s?=?53.6%). The plot of GC3s against ENc values revealed that GC compositional constraints are the main factor that determines the codon usage bias of major of helicase. Several indicators supported that the codon usage pattern of helicase is mainly subject to mutation pressure. Analysis of variation in codon usage and amino acid composition indicated AcMNPV helicase shows the significant preference for one or more postulated codons for each amino acid. A cluster analysis based on RSCU values suggested that AcMNPV is evolutionarily closer to members of group I alphabaculovirus. Comparison of the codon usage pattern among E. coli, yeast, mouse, human and AcMNPV showed that yeast is a suitable expression system for AcMNPV helicase. AcMNPV helicase shows weak codon usage bias. This study may help in elucidating the functional mechanism of AcMNPV helicase and the evolution of baculovirus helicases.     

Synonymous codon usage pattern among the S,M, and L segments in Crimean-congo hemorrhagic fever virus

Crimean-Congo hemorrhagic fever (CCHF) virus is one among the major zoonosis viral diseases that use the Hyalomma ticks as their transmission vector to cause viral infection to the human and mammalian community. The fatality of infectious is high across the world especially in Africa, Asia, Middle East, and Europe. This study regarding codon usage bias of S, M, and L segments of the CCHF virus pertaining to the host Homo sapiens, reveals in-depth information about the evolutionary characteristics of CCHFV. Relative Synonymous Codon Usage (RSCU), Effective number of codons (ENC) were calculated, to determine the codon usage pattern in each segment. Correlation analysis between Codon adaptation index (CAI), GRAVY (Hydrophobicity), AROMO (Aromaticity), and nucleotide composition revealed bias in the codon usage pattern. There was no strong codon bias found among any segments of the CCHF virus, indicating both the factors i.e., natural selection and mutational pressure shapes the codon usage bias.     

Differential trends in the codon usage patterns in HIV-1 genes

Host-pathogen interactions underlie one of the most complex evolutionary phenomena resulting in continual adaptive genetic changes, where pathogens exploit the host's molecular resources for growth and survival, while hosts try to eliminate the pathogen. Deciphering the molecular basis of host-pathogen interactions is useful in understanding the factors governing pathogen evolution and disease propagation. In host-pathogen context, a balance between mutation, selection, and genetic drift is known to maintain codon bias in both organisms. Studies revealing determinants of the bias and its dynamics are central to the understanding of host-pathogen evolution. We considered the Human Immunodeficiency Virus (HIV) type 1 and its human host to search for evolutionary signatures in the viral genome. Positive selection is known to dominate intra-host evolution of HIV-1, whereas high genetic variability underlies the belief that neutral processes drive inter-host differences. In this study, we analyze the codon usage patterns of HIV-1 genomes across all subtypes and clades sequenced over a period of 23 years. We show presence of unique temporal correlations in the codon bias of three HIV-1 genes illustrating differential adaptation of the HIV-1 genes towards the host preferred codons. Our results point towards gene-specific translational selection to be an important force driving the evolution of HIV-1 at the population level.     

Role of the human immunodeficiency virus type 1 envelope gene in viral fitness

A human host offers a variety of microenvironments to the infecting human immunodeficiency virus type 1 (HIV-1), resulting in various selective pressures, most of them directed against the envelope (env) gene. Therefore, it seems evident that the replicative capacity of the virus is largely related to viral entry. In this study we have used growth competition experiments and TaqMan real-time PCR detection to measure the fitness of subtype B HIV-1 primary isolates and autologous env-recombinant viruses in order to analyze the contribution of wild-type env sequences to overall HIV-1 fitness. A significant correlation was observed between fitness values obtained for wild-type HIV-1 isolates and those for the corresponding env-recombinant viruses (r = 0.93; P = 0.002). Our results suggest that the env gene, which is linked to a myriad of viral characteristics (e.g., entry into the host cell, transmission, coreceptor usage, and tropism), plays a major role in fitness of wild-type HIV-1. In addition, this new recombinant assay may be useful for measuring the contribution of HIV-1 env to fitness in viruses resistant to novel antiretroviral entry inhibitors.     

Comparative analysis of human coronaviruses focusing on nucleotide variability and synonymous codon usage patterns

The prevailing COVID-19 pandemic has drawn the attention of the scientific community to study the evolutionary origin of Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2). This study is a comprehensive quantitative analysis of the protein-coding sequences of seven human coronaviruses (HCoVs) to decipher the nucleotide sequence variability and codon usage patterns. It is essential to understand the survival ability of the viruses, their adaptation to hosts, and their evolution.The current analysis revealed a high abundance of the relative dinucleotide (odds ratio), GC and CT pairs in the first and last two codon positions, respectively, as well as a low abundance of the CG pair in the last two positions of the codon, which might be related to the evolution of the viruses. A remarkable level of variability of GC content in the third position of the codon among the seven coronaviruses was observed. Codons with high RSCU values are primarily from the aliphatic and hydroxyl amino acid groups, and codons with low RSCU values belong to the aliphatic, cyclic, positively charged, and sulfur-containing amino acid groups. In order to elucidate the evolutionary processes of the seven coronaviruses, a phylogenetic tree (dendrogram) was constructed based on the RSCU scores of the codons. The severe and mild categories CoVs were positioned in different clades. A comparative phylogenetic study with other coronaviruses depicted that SARS-CoV-2 is close to the CoV isolated from pangolins (Manis javanica, Pangolin-CoV) and cats (Felis catus, SARS(r)-CoV). Further analysis of the effective number of codon (ENC) usage bias showed a relatively higher bias for SARS-CoV and MERS-CoV compared to SARS-CoV-2. The ENC plot against GC3 suggested that the mutational bias might have a role in determining the codon usage variation among candidate viruses.A codon adaptability study on a few human host parasites (from different kingdoms), including CoVs, showed a diverse adaptability pattern. SARS-CoV-2 and SARS-CoV exhibit relatively lower but similar codon adaptability compared to MERS-CoV.     

Comparative study of the hemagglutinin and neuraminidase genes of influenza A virus H3N2, H9N2, and H5N1 subtypes using bioinformatics techniques

To investigate the genomic patterns of influenza A virus subtypes, such as H3N2, H9N2, and H5N1, we collected 1842 sequences of the hemagglutinin and neuraminidase genes from the NCBI database and parsed them into 7 categories: accession number, host species, sampling year, country, subtype, gene name, and sequence. The sequences that were isolated from the human, avian, and swine populations were extracted and stored in a MySQL database for intensive analysis. The GC content and relative synonymous codon usage (RSCU) values were calculated using JAVA codes. As a result, correspondence analysis of the RSCU values yielded the unique codon usage pattern (CUP) of each subtype and revealed no extreme differences among the human, avian, and swine isolates. H5N1 subtype viruses exhibited little variation in CUPs compared with other subtypes, suggesting that the H5N1 CUP has not yet undergone significant changes within each host species. Moreover, some observations may be relevant to CUP variation that has occurred over time among the H3N2 subtype viruses isolated from humans. All the sequences were divided into 3 groups over time, and each group seemed to have preferred synonymous codon patterns for each amino acid, especially for arginine, glycine, leucine, and valine. The bioinformatics technique we introduce in this study may be useful in predicting the evolutionary patterns of pandemic viruses.     

Codon usage bias in human cytomegalovirus and its biological implication

Human cytomegalovirus (HCMV) infection, a worldwide contagion, causes a serious disorder in infected individuals. Analysis of codon usage can reveal much molecular information about this virus. The effective number of codon (ENC) values, relative synonymous codon usage (RSCU) values, codon adaptation index (CAI), and nucleotide contents was investigated in approximately 160 coding sequences (CDS) among 17 human cytomegalovirus genomes using the software CodonW. Linear regression analysis and logistic regression were performed to explore the preliminary data. The results showed that, overall, HCMV genomes had low codon usage bias (mean ENC = 47.619). However, the ENC of individual CDS varied widely and was distributed unevenly between host-related genes and viral-self-function genes (P = 0.002, odds ratio (OR) = 3.194), as did the GC content (P = 0.016, OR = 2.178). The ENC values correlated with CAI, GC content, and the nucleotide composing at the 3rd codon position (GC3s) (P < 0.001). There was a significant variation in the codon preference that depended on the RSCU data. The predicted ENC curve suggested that mutational pressure, rather than natural selection, was one of the main factors that determined the codon usage bias in HCMV. Among 123 genes with known function, the genes related to viral self-replication and viral–host interaction showed different ENC and CAI values, and GC and GC3s contents. In conclusion, the detailed codon usage bias theoretically revealed information concerning HCMV evolution and could be a valuable additional parameter for HCMV gene function research.     

Synthesis and processing of human immunodeficiency virus type 1 envelope proteins encoded by a recombinant human adenovirus.

A recombinant adenovirus was constructed by inserting the human immunodeficiency virus type 1 (HIV-1) envelope gene downstream from the early region 3 (E3) promoter of adenovirus type 5 (Ad5), replacing the coding sequences of E3. The recombinant virus replicated as efficiently as the parent virus in all cell lines tested. Human cells infected with the recombinant virus synthesized the HIV-1 envelope precursor gp160, which was efficiently processed to the envelope glycoproteins gp120 and gp41. A human T-lymphoblast line (Molt-4) infected with the recombinant virus expressed HIV-1 envelope glycoproteins on the cell surface, leading to syncytium formation. The envelope gene was expressed from the E3 promoter at early times after infection and at late times from the major late promoter. When cotton rats were infected with the recombinant virus, antibodies against the HIV-1 envelope glycoproteins could be expressed in an immunoreactive form by the recombinant adenovirus, further illustrating the usefulness of adenoviruses as expression vectors.     

Genetic Imprint of Vaccination on Simian/Human Immunodeficiency Virus Type 1 Transmitted Viral Genomes in Rhesus Macaques

Understanding the genetic, antigenic and structural changes that occur during HIV-1 infection in response to pre-existing immunity will facilitate current efforts to develop an HIV-1 vaccine. Much is known about HIV-1 variation at the population level but little with regard to specific changes occurring in the envelope glycoprotein within a host in response to immune pressure elicited by antibodies. The aim of this study was to track and map specific early genetic changes occurring in the viral envelope gene following vaccination using a highly controlled viral challenge setting in the SHIV macaque model. We generated 449 full-length env sequences from vaccinees, and 63 from the virus inoculum. Analysis revealed a different pattern in the distribution and frequency of mutations in the regions of the envelope gene targeted by the vaccine as well as different patterns of diversification between animals in the naïve control group and vaccinees. Given the high stringency of the model it is remarkable that we were able to identify genetic changes associated with the vaccination. This work provides insight into the characterization of breakthrough viral populations in less than fully efficacious vaccines and illustrates the value of HIV-1 Env SHIV challenge model in macaques to unravel the mechanisms driving HIV-1 envelope genetic diversity in the presence of vaccine induced-responses.     

Comparison of Correspondence Analysis Methods for Synonymous Codon Usage in Bacteria

Synonymous codon usage varies both between organisms and among genes within a genome, and arises due to differences in G + C content, replication strand skew, or gene expression levels. Correspondence analysis (CA) is widely used to identify major sources of variation in synonymous codon usage among genes and provides a way to identify horizontally transferred or highly expressed genes. Four methods of CA have been developed based on three kinds of input data: absolute codon frequency, relative codon frequency, and relative synonymous codon usage (RSCU) as well as within-group CA (WCA). Although different CA methods have been used in the past, no comprehensive comparative study has been performed to evaluate their effectiveness. Here, the four CA methods were evaluated by applying them to 241 bacterial genome sequences. The results indicate that WCA is more effective than the other three methods in generating axes that reflect variations in synonymous codon usage. Furthermore, WCA reveals sources that were previously unnoticed in some genomes; e.g. synonymous codon usage related to replication strand skew was detected in Rickettsia prowazekii. Though CA based on RSCU is widely used, our evaluation indicates that this method does not perform as well as WCA.Key words: correspondence analysis, synonymous codon usage, horizontal gene transfer, strand-specific mutational bias, translational selection     

N-linked glycosylation of the HIV type-1 gp120 envelope glycoprotein as a major determinant of CCR5 and CXCR4 coreceptor utilization

The variable V1V2 and V3 regions of the human immunodeficiency virus type-1 (HIV-1) envelope glycoprotein (gp120) can influence viral coreceptor usage. To substantiate this we generated isogenic HIV-1 molecularly cloned viruses that were composed of the HxB2 envelope backbone containing the V1V2 and V3 regions from viruses isolated from a patient progressing to disease. We show that the V3 amino acid charge per se had little influence on altering the virus coreceptor phenotype. The V1V2 region and its N-linked glycosylation degree were shown to confer CXCR4 usage and provide the virus with rapid replication kinetics. Loss of an N-linked glycosylation site within the V3 region had a major influence on the virus switching from the R5 to X4 phenotype in a V3 charge-dependent manner. The loss of this V3 N-linked glycosylation site was also linked with the broadening of the coreceptor repertoire to incorporate CCR3. By comparing the amino acid sequences of primary HIV-1 isolates, we identified a strong association between high V3 charge and the loss of this V3 N-linked glycosylation site. These results demonstrate that the N-linked glycosylation pattern of the HIV-1 envelope can strongly influence viral coreceptor utilization and the R5 to X4 switch.     

Convergent host-parasite codon usage between honeybee and bee associated viral genomes

By correlating the codon usage in four insects (the honeybee, red flour beetle, mosquito and fruit fly) with six honeybee host specific viruses, we found that the codon usage patterns of the bee viruses were strongly related to that of the honeybee and only weakly related to the red flour beetle. The insects shared varying degrees of codon usage similarity which roughly follow the known phylogenetic relatedness. All of the codon usage similarity can be described by relatedness-by-descent except for the high codon usage similarity between the honeybee and honeybee associated viruses. This evidence for the convergent evolution of the honeybee viruses toward the codon usage of the honeybee suggests that small host specific viral genomes have the freedom to quickly optimize codon usage to successfully parasitize their preferred host. The codon usage co-evolution of the six host specific honeybee viruses towards the codon usage of the honeybee described in this paper is the first evidence for codon usage correlation between an insect host and a single stranded RNA virus.     

Human immunodeficiency virus type 1 (HIV-1) diversity at time of infection is not restricted to certain risk groups or specific HIV-1 subtypes

African women frequently acquire several genetically distinct human immunodeficiency virus type 1 (HIV-1) variants from a heterosexual partner, whereas the acquisition of multiple variants appears to be rare in men. To determine whether newly infected individuals in other risk groups acquire genetically diverse viruses, we examined the viral envelope sequences in plasma samples from 13 women and 4 men from the United States infected with subtype B viruses and 10 men from Kenya infected with non-subtype B viruses. HIV-1 envelope sequences differed by more than 2% in three U.S. women, one U.S. man, and one Kenyan man near the time of seroconversion. These findings suggest that early HIV-1 genetic diversity is not exclusive to women from Africa or to infection with any particular HIV-1 subtype.     

子囊菌延长因子1 alpha基因密码子偏好性分析(英文)

文中对子囊菌代表类群的延伸因子1 alpha基因密码子的使用模式进行了研究。结果表明:该基因的密码子使用偏好性不仅与核酸碱基组成密切相关,也受到其他选择性压力的影响。统计分析揭示了子囊菌各类群该基因的密码子组成和编码特点,在同义密码子的选择模式上,酵母纲(Saccharomycetes)的成员具有较独特的偏好性。基于密码子用法分歧度的聚类分析方法较合理地反映了大部分类群的分类学地位,但在各个纲的内部,密码子偏好性的变化程度存在差异。     

Dendritic cells preferentially transfer CXCR4-using human immunodeficiency virus type 1 variants to CD4+ T lymphocytes in trans

Human immunodeficiency virus type 1 (HIV-1) preferentially utilizes the CCR5 coreceptor for target cell entry in the acute phase of infection, while later in disease progression the virus switches to the CXCR4 coreceptor in approximately 50% of patients. In response to HIV-1 the adaptive immune response is triggered, and antibody (Ab) production is elicited to block HIV-1 entry. We recently determined that dendritic cells (DCs) can efficiently capture Ab-neutralized HIV-1, restore infectivity, and transmit infectious virus to target cells. Here, we tested the effect of Abs on trans transmission of CCR5 or CXCR4 HIV-1 variants. We observed that transmission of HIV-1 by immature as well as mature DCs was significantly higher for CXCR4- than CCR5-tropic viral strains. Additionally, neutralizing Abs directed against either the gp41 or gp120 region of the envelope such as 2F5, 4E10, and V3-directed Abs inhibited transmission of CCR5-tropic HIV-1, whereas Ab-treated CXCR4-tropic virus demonstrated unaltered or increased transmission. To further study the effects of coreceptor usage we tested molecularly cloned HIV-1 variants with modifications in the envelope that were based on longitudinal gp120 V1 and V3 variable loop sequences from a patient progressing to AIDS. We observed that DCs preferentially facilitated infection of CD4⁺ T lymphocytes of viral strains with an envelope phenotype found late in disease. Taken together, our results illustrate that DCs transmit CXCR4-tropic HIV-1 much more efficiently than CCR5 strains; we hypothesize that this discrimination could contribute to the in vivo coreceptor switch after seroconversion and could be responsible for the increase in viral load.