Comparative analysis of codon usage patterns in Rift Valley fever virus.

Rift Valley fever virus (RVFV) is a vector-borne pathogen and is the most widely known virus in the genus Phlebovirus. Since it was first reported, RVFV has spread to western Africa, Egypt and Madagascar from its traditional endemic region, and infections continue to occur in new areas. In this study, we analyzed genomic patterns according to the infection properties of RVFV. Among the four segments of RVFV, the nucleotide composition, overall GC content and the difference of GC composition in the third position of the codons (%GC3) between groups were the largest in the S (NP) segment, showing that more diverse codons were used than in other segments. Furthermore, the results of CAI analysis of the S (NP) segment showed that viruses isolated from regions where no previous infections had been reported had the highest values, indicating greater adaptability to human hosts compared with other viruses. This result suggests that mutations in the S (NP) segment co-evolve with the infected hosts and may lead to expansion of the geographic range. The distinctive codon usage patterns observed in specific genomic regions of a group with similar infection properties may be related to the increasing likelihood of RVFV infections in new areas.

Codon usage patterns of LT-Ag genes in polyomaviruses from different host species.

BACKGROUND: Polyomaviruses (PyVs) have a wide range of hosts, from humans to fish, and their effects on hosts vary. The differences in the infection characteristics of PyV with respect to the host are assumed to be influenced by the biochemical function of the LT-Ag protein, which is related to the cytopathic effect and tumorigenesis mechanism via interaction with the host protein. METHODS: We carried out a comparative analysis of codon usage patterns of large T-antigens (LT-Ags) of PyVs isolated from various host species and their functional domains and sequence motifs. Parity rule 2 (PR2) and neutrality analysis were applied to evaluate the effects of mutation and selection pressure on codon usage bias. To investigate evolutionary relationships among PyVs, we carried out a phylogenetic analysis, and a correspondence analysis of relative synonymous codon usage (RSCU) values was performed. RESULTS: Nucleotide composition analysis using LT-Ag gene sequences showed that the GC and GC3 values of avian PyVs were higher than those of mammalian PyVs. The effective number of codon (ENC) analysis showed host-specific ENC distribution characteristics in both the LT-Ag gene and the coding sequences of its domain regions. In the avian and fish PyVs, the codon diversity was significant, whereas the mammalian PyVs tended to exhibit conservative and host-specific evolution of codon usage bias. The results of our PR2 and neutrality analysis revealed mutation bias or highly variable GC contents by showing a narrow GC12 distribution and wide GC3 distribution in all sequences. Furthermore, the calculated RSCU values revealed differences in the codon usage preference of the LT-AG gene according to the host group. A similar tendency was observed in the two functional domains used in the analysis. CONCLUSIONS: Our study showed that specific domains or sequence motifs of various PyV LT-Ags have evolved so that each virus protein interacts with host cell targets. They have also adapted to thrive in specific host species and cell types. Functional domains of LT-Ag, which are known to interact with host proteins involved in cell proliferation and gene expression regulation, may provide important information, as they are significantly related to the host specificity of PyVs.

Analysis of the codon usage pattern of the RdRP gene of mycovirus infecting Aspergillus spp.

BACKGROUND: Mycoviruses that infect fungi generally do not have a significant effect on the host and, instead, reduce the toxicity of the fungi. However, recent studies have shown that polymycovirus-1, a mycovirus that infects Aspergillus species known to cause disease in humans, is related to increased virulence of the fungus. METHODS: Comparative analysis was performed of RdRP gene codon usage patterns of Aspergillus fumigatus polymycovirus-1 (AfuPmV-1) and other mycoviruses known to infect Aspergillus spp. to examine the genetic characteristics of AfuPmV-1. In addition, codon usage analysis was performed to determine whether the nucleotide composition and codon usage characteristics of AfuPmV-1 were also present in other polymycoviruses and hypervirulence-related mycoviruses. Phylogenetic analysis was also performed to investigate their evolutionary relationship. RESULTS: Analysis of nucleotide composition indicated that AfuPmV-1 had the highest GC content among analyzed mycoviruses and relative synonymous codon usage analysis indicated that all of the codons preferred by AfuPmV-1 ended with C or G, while codons ending with A or U were not observed. Moreover, the effective number of codons, the codon adaptation index, and correspondence analysis showed that AfuPmV-1 had greater codon preference compared with other mycoviruses and that AfuPmV-1 had relatively high adaptability to humans and fungi. These results were generally similar among polymycoviruses. CONCLUSIONS: The codon usage pattern of AfuPmV-1 differs from other mycoviruses that infect Aspergillus spp. This difference may be related to the hypervirulence effect of AfuPmV-1. Analysis of AfuPmV-1 codon usage patterns could contribute to the identification and prediction of virulence effects of mycoviruses with similar genetic characteristics.

The evolutionary and genetic patterns of African swine fever virus.

African swine fever (ASF) is a serious animal disease, and has spread to Africa, Europe and Asia, causing massive economic losses. African swine fever virus (ASFV) is transmitted from a reservoir host (warthog) to domestic pigs via a sylvatic cycle (transmission between warthogs and soft ticks) and a domestic cycle (transmission between domestic pigs) and survives by expressing a variety of genes related to virus-host interactions. We evaluated differences in codon usage patterns among ASFV genotypes and clades and explored the common and specific evolutionary and genetic characteristics of ASFV sequences. We analysed the evolutionary relationships, nucleotide compositions, codon usage patterns, selection pressures (mutational pressure and natural selection) and viral adaptation to host codon usage based on the coding sequences (CDS) of key functional genes of ASFV. AT bias was detected in the six genes analysed, irrespective of clade. The AT bias of genes (A224L, A179L, EP153R) encoding proteins involved in interaction with host cells after infection was high; among them, the AT bias of EP153R was the greatest at 78.3%. A large number of overrepresented codons were identified in EP153R, whereas there were no overrepresented codons with a relative synonymous codon usage (RSCU) value of ≥3 in B646L. In most genes, the pattern of selection pressure was similar for each clade, but in EP153R, diverse patterns of selection pressure were captured within the same clade and genotype. As a result of evaluating host adaptation based on the codon adaptation index (CAI), for B646L, E183L, CP204L and A179L, the codon usage patterns in all sequences were more similar to tick than domestic pig or wild boar. However, EP153R showed the lowest average CAI value of 0.52 when selecting tick as a reference set. The genes analysed in this study showed different magnitudes of selection pressure at the clade and genotype levels, which is likely to be related to the function of the encoded proteins and may determine key evolutionary traits of viruses, such as the level of genetic variation and host range. The diversity of codon adaptations at the genetic level in ASFV may account for differences in translational selection in ASFV hosts and provides insight into viral host adaptation and co-evolution.

Utilization of systematic error-assessment software to improve phylogenetic accuracy.

Unlike classical systems based on the use of morphological data, modern phylogenetic analyses use genetic information to construct phylogenetic trees. Ongoing research in the field of phylogenetics is evaluating the accuracy of phylogenetic estimation results and the reliability of phylogenetic trees to explain evolutionary relationships. Recently, the probability of stochastic errors in large-scale phylogenetic datasets has decreased, while the probability of systematic errors has increased. Therefore, before constructing a phylogenetic tree, it is necessary to assess the causes of systematic bias to improve the accuracy of phylogenetic estimates. We performed analyses of three datasets (Terebelliformia, Daphniid, and Glires clades) using bioinformatics software to assess systematic error and improve phylogenetic tree accuracy. Then, we proposed a combination of systematic biases capable of discerning the most suitable gene markers within a series of taxa and generating conflicting phylogenetic topologies. Our findings will help improve the reliability of phylogenetic software to estimate phylogenies more accurately by exploiting systematic bias.

Anti-adhesive functions of CD43 expressed on colon carcinoma cells through the modulation of integrins.

CD43 has conflicting roles in both pro- and anti-adhesive function in cell-to-cell adhesion in hematopoietic cells. We examined the role of CD43 glycoprotein in a colorectal carcinoma cell line. We expressed human CD43 antigen on HT-29 cells, a colon adenocarcinoma cell line, and compared the adhesion to the extracellular matrix with that of mock-transduced cells in vitro. CD43 expression inhibited the adhesion to extracellular matrix, such as collagen type IV and laminin. As the expression of ß1 integrin was downregulated in CD43-expressing HT-29 cells, the anti-adhesive effect of CD43 might be implicated in its expression. Our findings suggest that the anti-adhesive function of CD43 in colon carcinoma cells plays a role in the tumorigenesis and metastasis of colorectal carcinoma cells.

Analysis of evolutionary and genetic patterns in structural genes of primate lentiviruses.

BACKGROUND: Primate lentiviruses (HIV1, HIV2, and Simian immunodeficiency virus [SIV]) cause immune deficiency, encephalitis, and infectious anemia in mammals such as cattle, cat, goat, sheep, horse, and puma. OBJECTIVE: This study was designed and conducted with the main purpose of confirming the overall codon usage pattern of primate lentiviruses and exploring the evolutionary and genetic characteristics commonly or specifically expressed in HIV1, HIV2, and SIV. METHODS: The gag, pol, and env gene sequences of HIV1, HIV2, and SIV were analyzed to determine their evolutionary relationships, nucleotide compositions, codon usage patterns, neutrality, selection pressure (influence of mutational pressure and natural selection), and viral adaptation to human codon usage. RESULTS: A strong 'A' bias was confirmed in all three structural genes, consistent with previous findings regarding HIV. Notably, the ENC-GC3s plot and neutral evolution analysis showed that all primate lentiviruses were more affected by selection pressure than by mutation caused by the GC composition of the gene, consistent with prior reports regarding HIV1. The overall codon usage bias of pol was highest among the structural genes, while the codon usage bias of env was lowest. The virus groups showing high codon bias in all three genes were HIV1 and SIVcolobus. The codon adaptation index (CAI) and similarity D(A, B) values indicated that although there was a high degree of similarity to human codon usage in all three structural genes of HIV, this similarity was not caused by translation pressure. In addition, compared with HIV1, the codon usage of HIV2 is more similar to the human codon usage, but the overall codon usage bias is lower. CONCLUSION: The origin viruses of HIV (SIVcpz_gor and SIVsmm) exhibit greater similarity to human codon usage in the gag gene, confirming their robust adaptability to human codon usage. Therefore, HIV1 and HIV2 may have evolved to avoid human codon usage by selection pressure in the gag gene after interspecies transmission from SIV hosts to humans. By overcoming safety and stability issues, information from codon usage analysis will be useful for attenuated HIV1 vaccine development. A recoded HIV1 variant can be used as a vaccine vector or in immunotherapy to induce specific innate immune responses. Further research regarding HIV1 dinucleotide usage and codon pair usage will facilitate new approaches to the treatment of AIDS.

Analysis of protein determinants of host-specific infection properties of polyomaviruses using machine learning.

BACKGROUND: The large tumor antigen (LT-Ag) and major capsid protein VP1 are known to play important roles in determining the host-specific infection properties of polyomaviruses (PyVs). OBJECTIVE: The objective of this study was to investigate the physicochemical properties of amino acids of LT-Ag and VP1 that have important effects on host specificity, as well as classification techniques used to predict PyV hosts. METHODS: We collected and used reference sequences of 86 viral species for analysis. Based on the clustering pattern of the reconstructed phylogenetic tree, the dataset was divided into three groups: mammalian, avian, and fish. We then used random forest (RF), naïve Bayes (NB), and k-nearest neighbors (kNN) algorithms for host classification. RESULTS: Among the three algorithms, classification accuracy using kNN was highest in both LT-Ag (ACC = 98.83) and VP1 (ACC = 96.51). The amino acid physicochemical property most strongly correlated with host classification was charge, followed by solvent accessibility, polarity, and hydrophobicity in LT-Ag. However, in VP1, amino acid composition showed the highest correlation with host classification, followed by charge, normalized van der Waals volume, and solvent accessibility. CONCLUSIONS: The results of the present study suggest the possibility of determining or predicting the host range and infection properties of PyVs at the molecular level by identifying the host species of active and emerging PyVs that exhibit different infection properties among diverse host species. Structural and biochemical differences of LT-Ag and VP1 proteins in host species that reflect these amino acid properties can be considered primary factors that determine the host specificity of PyV.

Simulation Model of Bacterial Resistance to Antibiotics Using Individual-Based Modeling.

We designed and implemented simulation models of bacterial growth and antibiotic resistance to determine the appropriate antibiotics to use against antibiotic-resistant bacteria. Simulation models were designed using individual-based modeling, and a simulation tool, ARSim, was developed to conduct experiments using the models. Simulations of bacterial growth were conducted by virtually growing Klebsiella pneumoniae bacteria in a virtual environment with predefined parameters. Other experiments included predicting the effects of antibiotics when added to two different groups, one group of nonresistant bacteria and another group of both resistant and nonresistant bacteria. Carbapenem class antibiotics such as Imipenem were used for the simulation. The simulation results showed that the biological principles of bacteria and their antibiotic resistance mechanisms were correctly designed and implemented. Using the computational approaches developed in this study, we hope to provide researchers with a more effective method for finding new ways to fight antibiotic resistance.

Application of the multiensemble sampling to the equilibrium folding of proteins.

MOTIVATION: Conventional Monte Carlo and molecular dynamics simulations of proteins in the canonical ensemble are of little use, because they tend to get trapped in states of energy local minima at low temperatures. One way to surmount this difficulty is to use a non-Boltzmann sampling method in which conformations are sampled upon a general weighting function instead of the conventional Boltzmann weighting function. The multiensemble sampling (MES) method is a non-Boltzmann sampling method that was originally developed to estimate free energy differences between systems with different potential energies and/or at different thermodynamic states. The method has not yet been applied to studies of complex molecular systems such as proteins. RESULTS: MES Monte Carlo simulations of small proteins have been carried out using a united-residue force field. The proteins at several temperatures from the unfolded to the folded states were simulated in a single MC run at a time and their equilibrium thermodynamic properties were calculated correctly. The distributions of sampled conformations clearly indicate that, when going through states of energy local minima, the MES simulation did not get trapped in them but escaped from them so quickly that all the relevant parts of conformation space could be sampled properly. A two-step folding process consisting of a collapse transition followed by a folding transition is observed. This study demonstrates that the use of MES alleviates the multiple-minima problem greatly. AVAILABILITY: Available on request from the authors.

Protein backbone torsion angle-based structure comparison and secondary structure database web server.

Structural information has been a major concern for biological and pharmaceutical studies for its intimate relationship to the function of a protein. Three-dimensional representation of the positions of protein atoms is utilized among many structural information repositories that have been published. The reliability of the torsional system, which represents the native processes of structural change in the structural analysis, was partially proven with previous structural alignment studies. Here, a web server providing structural information and analysis based on the backbone torsional representation of a protein structure is newly introduced. The web server offers functions of secondary structure database search, secondary structure calculation, and pair-wise protein structure comparison, based on a backbone torsion angle representation system. Application of the implementation in pair-wise structural alignment showed highly accurate results. The information derived from this web server might be further utilized in the field of ab initio protein structure modeling or protein homology-related analyses.

Correlation analysis for the incubation period of prion disease.

Previous studies have shown that genetic quantitative trait loci (QTL), strain barriers, inoculation dose and inoculation method modulate the incubation period of prion diseases. We examined the relationship between a diverse set of physical, genetic and immunological characteristics and the incubation period of prion disease using correlation analyses. We found that incubation period was highly correlated with brain weight. In addition, mean corpuscular volume and cell size were strongly correlated with incubation period, indicating that the physical magnitude of prion-infected organs or individual cells may be important in determining the incubation period. Given the same prion inoculation dose, animals with a lower brain weight, mean corpuscular volume or cell size may experience more virulent disease, as the effective concentration of abnormal prion, which might regulate the attachment rate of prions to aggregates, is increased with smaller capacity of brains and cells. This is partly consistent with previous theoretical modeling. The strong correlations between incubation period and physical properties of the brain and cells in this study suggest that the mechanism underlying prion disease pathology may be physical, indicating that the incubation process is governed by simple chemical stoichiometry.

[The current trend of avian influenza viruses in bioinformatics research].

OBJECTIVES: Since the first human infection from avian influenza was reported in Hong Kong in 1997, many Asian countries have confirmed outbreaks of highly pathogenic H5N1 avian influenza viruses. In addition to Asian countries, the EU authorities also held an urgent meeting in February 2006 at which it was agreed that Europe could also become the next target for H5N1 avian influenza in the near future. In this paper, we provide the general and applicable information on the avian influenza in the bioinformatics field to assist future studies in preventive medicine. METHODS: We introduced some up-to-date analytical tools in bioinformatics research, and discussed the current trends of avian influenza outbreaks. Among the bioinformatics methods, we focused our interests on two topics: pattern analysis using the secondary database of avian influenza, and structural analysis using the molecular dynamics simulations in vaccine design. RESULTS: Use of the public genome databases available in the bioinformatics field enabled intensive analysis of the genetic patterns. Moreover, molecular dynamic simulations have also undergone remarkable development on the basis of the high performance supercomputing infrastructure these days. CONCLUSIONS: The bioinformatics techniques we introduced in this study may be useful in preventive medicine, especially in vaccine and drug discovery.

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.

PKB/Akt inhibits ceramide-induced apoptosis in neuroblastoma cells by blocking apoptosis-inducing factor (AIF) translocation.

Ceramide is a sphingolipid that is abundant in the plasma membrane of neuronal cells and is thought to have regulatory roles in cell differentiation and cell death. Ceramide is known to induce apoptosis in a variety of different cell types, whereas the physiological significance of gangliosides, another class of sphingolipids, in these processes is still unclear. We examined the mechanisms of ceramide-induced cell death using a human neuroblastoma cell line. Treatment of the human neuroblastoma cell line SH-SY5Y with ceramide induced dephosphorylation of the PKB/Akt kinase and subsequent mitochondrial dysfunction. In addition, ceramide-induced neuronal cell death was not completely blocked by inhibition of caspase activity. This incomplete inhibition appeared to be attributable to the translocation of apoptosis-inducing factor to the nucleus. Furthermore, overexpression of active PKB/Akt or Bcl-2 successfully blocked ceramide-induced neuronal cell death through inhibition of the translocation of apoptosis-inducing factor.

Genomic analysis of influenza A viruses, including avian flu (H5N1) strains.

This study was designed to conduct genomic analysis in two steps, such as the overall relative synonymous codon usage (RSCU) analysis of the five virus species in the orthomyxoviridae family, and more intensive pattern analysis of the four subtypes of influenza A virus (H1N1, H2N2, H3N2, and H5N1) which were isolated from human population. All the subtypes were categorized by their isolated regions, including Asia, Europe, and Africa, and most of the synonymous codon usage patterns were analyzed by correspondence analysis (CA). As a result, influenza A virus showed the lowest synonymous codon usage bias among the virus species of the orthomyxoviridae family, and influenza B and influenza C virus were followed, while suggesting that influenza A virus might have an advantage in transmitting across the species barrier due to their low codon usage bias. The ENC values of the host-specific HA and NA genes represented their different HA and NA types very well, and this reveals that each influenza A virus subtype uses different codon usage patterns as well as the amino acid compositions. In NP, PA and PB2 genes, most of the virus subtypes showed similar RSCU patterns except for H5N1 and H3N2 (A/HK/1774/1999) subtypes which were suspected to be transmitted across the species barrier, from avian and porcine species to human beings, respectively. This distinguishable synonymous codon usage patterns in non-human origin viruses might be useful in determining the origin of influenza A viruses in genomic levels as well as the serological tests. In this study, all the process, including extracting sequences from GenBank flat file and calculating codon usage values, was conducted by Java codes, and these bioinformatics-related methods may be useful in predicting the evolutionary patterns of pandemic viruses.

Structural characterization of the molten globule state of apomyoglobin by limited proteolysis and HPLC-mass spectrometry.

A method to characterize the structural conformation of an acidic molten globule apomyoglobin (apoMb) at pH 4.2 was developed using limited proteolysis and HPLC-mass spectrometry (HPLC-MS). Endoproteinase Glu-C, which has a double maximum activity at pH 4.0 and pH 7.8 toward glutamic acid (Glu), was used as a proteolytic enzyme. Using this method enabled us to compare the proteolytic cleavages of native apoMb (at pH 8.0) and molten globule (at pH 4.2) directly. Only the first cleavage event in each molecule was considered as reflecting original structural information since the original structure of the protein can be altered after the fist cleavage. Structural changes of apoMb in various pH conditions were studied here to elucidate the local helicity of molten globule apoMb. Among 13 Glu sites, only Glu83 and Glu85 in the F-helix were cleaved at pH 8.0, which confirms that only helix F is frayed upon removal of heme group. At acidic molten globule state, rapid cleavages at Glu38, Glu52, Glu54, Glu85, and Glu148 were detected, while the remaining eight sites were protected. Glu6 and Glu18 in the A-helix, and Glu105 in the G-helix were protected due to the helicity of the secondary structures. The cleavage at Glu38 and the protection at Glu41 in the C-helix indicate that the first half of the C-helix is frayed and the second half of the C-helix is structured. Cleavage at both Glu52 and Glu54 in the D-helix proves that the D-helix is disordered. The N-terminal end of the E-helix at Glu59 was protected, and the beginning of the F-helix was protected by aid of the pH-induced C-cap of the E-helix. The cleavage at Glu148 in H suggests that the C-terminal end of the H-helix is disordered. The A-helix and the first half of the B-helix were highly stable.