Genomic analysis of codon usage shows influence of mutation pressure, natural selection, and host features on Marburg virus evolution.

BACKGROUND: The Marburg virus (MARV) has a negative-sense single-stranded RNA genome, belongs to the family Filoviridae, and is responsible for several outbreaks of highly fatal hemorrhagic fever. Codon usage patterns of viruses reflect a series of evolutionary changes that enable viruses to shape their survival rates and fitness toward the external environment and, most importantly, their hosts. To understand the evolution of MARV at the codon level, we report a comprehensive analysis of synonymous codon usage patterns in MARV genomes. Multiple codon analysis approaches and statistical methods were performed to determine overall codon usage patterns, biases in codon usage, and influence of various factors, including mutation pressure, natural selection, and its two hosts, Homo sapiens and Rousettus aegyptiacus. RESULTS: Nucleotide composition and relative synonymous codon usage (RSCU) analysis revealed that MARV shows mutation bias and prefers U- and A-ended codons to code amino acids. Effective number of codons analysis indicated that overall codon usage among MARV genomes is slightly biased. The Parity Rule 2 plot analysis showed that GC and AU nucleotides were not used proportionally which accounts for the presence of natural selection. Codon usage patterns of MARV were also found to be influenced by its hosts. This indicates that MARV have evolved codon usage patterns that are specific to both of its hosts. Moreover, selection pressure from R. aegyptiacus on the MARV RSCU patterns was found to be dominant compared with that from H. sapiens. Overall, mutation pressure was found to be the most important and dominant force that shapes codon usage patterns in MARV. CONCLUSIONS: To our knowledge, this is the first detailed codon usage analysis of MARV and extends our understanding of the mechanisms that contribute to codon usage and evolution of MARV.

A 2-year retrospective study of viral and host-associated risk factors in Pakistani hepatocellular carcinoma patients.

OBJECTIVE: Persistent chronic hepatitis C (CHC) infection is associated strongly with serious complications such as hepatitis C virus-associated liver cirrhosis (HCV-LC) and hepatitis C virus-associated hepatocellular carcinoma (HCV-HCC). The aim of this study was to assess the distribution of hepatitis C virus (HCV) genotypes among HCV-positive patients and examine the potential associations between viral and host-associated factors with the risk of developing HCV-HCC. PATIENTS AND METHODS: HCV-positive patients (n = 300) were enrolled and divided into three groups: CHC (n = 171), HCV-LC (n = 51), and HCV-HCC (n = 78). RESULTS: HCV genotype 3a showed the highest prevalence among HCV-positive individuals (66% of patients), followed by genotype 1a (15% of patients). The proportion of individuals infected with mixed HCV genotypes was higher among HCV-HCC patients. Interestingly, there were a significantly higher proportion of women (54/78; 69.2%) among HCV-HCC patients compared with CHC patients (89/171 or 52%; χ = 6.47; P=1 × 10). Women with HCV had two-fold higher odds of developing HCV-HCC (odds ratio = 2.07, 95% confidence interval: 1.18-3.71). In comparison with CHC patients, significantly more HCV-HCC patients were 50 years of age or older (59/78 or 75.6% of HCV-HCC patients and 61/171 or 35.7% of CHC patients; χ = 34.27; P < 0.0001), suggesting that HCV-positive patients aged 50 years or older had an ~five-fold higher risk of developing HCV-HCC (odds ratio = 5.6, 95% confidence interval: 3.02-10.01). CONCLUSION: In summary, HCV genotype 3a had the highest prevalence in the studied HCV-positive population, and women and older patients were at a higher risk of developing HCV-LC and HCV-HCC following CHC infections.

Evolution of codon usage in Zika virus genomes is host and vector specific.

The codon usage patterns of viruses reflect the evolutionary changes that allow them to optimize their survival and adapt their fitness to the external environment and, most importantly, their hosts. Here we report the genotype-specific codon usage patterns of Zika virus (ZIKV) strains from the current and previous outbreaks. Several genotype-specific and common codon usage traits were noted in the ZIKV coding sequences, indicating their independent evolutionary origins from a common ancestor. The overall influence of natural selection was more profound than that of mutation pressure, acting on a specific set of viral genes in the Asian-genotype ZIKV strains from the recent outbreak. An interplay between codon adaptation and deoptimization may have allowed the virus to adapt to multiple host and vectors and is reported for the first time in ZIKV genomes. Combining our codon analysis with geographical data on Aedes populations in the Americas suggested that ZIKV has evolved host- and vector-specific codon usage patterns to maintain successful replication and transmission chains within multiple hosts and vectors.

Genome-wide analysis of codon usage and influencing factors in chikungunya viruses.

Chikungunya virus (CHIKV) is an arthropod-borne virus of the family Togaviridae that is transmitted to humans by Aedes spp. mosquitoes. Its genome comprises a 12 kb single-strand positive-sense RNA. In the present study, we report the patterns of synonymous codon usage in 141 CHIKV genomes by calculating several codon usage indices and applying multivariate statistical methods. Relative synonymous codon usage (RSCU) analysis showed that the preferred synonymous codons were G/C and A-ended. A comparative analysis of RSCU between CHIKV and its hosts showed that codon usage patterns of CHIKV are a mixture of coincidence and antagonism. Similarity index analysis showed that the overall codon usage patterns of CHIKV have been strongly influenced by Pan troglodytes and Aedes albopictus during evolution. The overall codon usage bias was low in CHIKV genomes, as inferred from the analysis of effective number of codons (ENC) and codon adaptation index (CAI). Our data suggested that although mutation pressure dominates codon usage in CHIKV, patterns of codon usage in CHIKV are also under the influence of natural selection from its hosts and geography. To the best of our knowledge, this is first report describing codon usage analysis in CHIKV genomes. The findings from this study are expected to increase our understanding of factors involved in viral evolution, and fitness towards hosts and the environment.

Comparative genomics analysis of Mycobacterium ulcerans for the identification of putative essential genes and therapeutic candidates.

Mycobacterium ulcerans, the causative agent of Buruli ulcer, is the third most common mycobacterial disease after tuberculosis and leprosy. The present treatment options are limited and emergence of treatment resistant isolates represents a serious concern and a need for better therapeutics. Conventional drug discovery methods are time consuming and labor-intensive. Unfortunately, the slow growing nature of M. ulcerans in experimental conditions is also a barrier for drug discovery and development. In contrast, recent advancements in complete genome sequencing, in combination with cheminformatics and computational biology, represent an attractive alternative approach for the identification of therapeutic candidates worthy of experimental research. A computational, comparative genomics workflow was defined for the identification of novel therapeutic candidates against M. ulcerans, with the aim that a selected target should be essential to the pathogen, and have no homology in the human host. Initially, a total of 424 genes were predicted as essential from the M. ulcerans genome, via homology searching of essential genome content from 20 different bacteria. Metabolic pathway analysis showed that the most essential genes are associated with carbohydrate and amino acid metabolism. Among these, 236 proteins were identified as non-host and essential, and could serve as potential drug and vaccine candidates. Several drug target prioritization parameters including druggability were also calculated. Enzymes from several pathways are discussed as potential drug targets, including those from cell wall synthesis, thiamine biosynthesis, protein biosynthesis, and histidine biosynthesis. It is expected that our data will facilitate selection of M. ulcerans proteins for successful entry into drug design pipelines.

Mycoplasma genitalium: a comparative genomics study of metabolic pathways for the identification of drug and vaccine targets.

Increasing emergence of antibiotic-resistant pathogenic microorganisms is one of the biggest challenges for biomedical research and drug development. Traditional drug discovery methods are time-consuming, expensive and often yield few drug targets. In contrast, advances in complete genome sequencing, bioinformatics and cheminformatics represent an attractive alternative approach to identify drug targets worthy of experimental follow-up. Mycoplasma genitalium is a human parasitic pathogen that is associated with several sexually transmitted diseases. Recently, emergence of treatment-resistant isolates has been reported, which raises serious concern and a need for identification of additional drug targets. In the present study, a systematic workflow consisting of comparative genomics, metabolic pathways analysis and additional drug prioritizing parameters was defined for the identification of novel drug and vaccine targets that are essential for M. genitalium, but absent in its human host. In silico analyses and manual mining identified 79 proteins of M. genitalium, which showed no similarity to human proteins. Among these, 67 proteins were identified as non-homologous essential proteins that could serve as potential drug and vaccine targets. Subcellular localization, molecular weight, and three-dimensional structural characteristics that could facilitate filtering of attractive drug targets were also calculated for the non-homologous essential proteins. Enzymes from thiamine biosynthesis, protein biosynthesis, and folate biosynthesis were identified as attractive candidates for drug development. Furthermore, druggability of each of the identified drug targets was also evaluated by the DrugBank database. Results from this study could facilitate selection of M. genitalium proteins for entry into drug design and vaccine production pipelines.

Computational identification of interplay between phosphorylation and O-ß-glycosylation of human occludin as potential mechanism to impair hepatitis C virus entry.

Hepatitis C virus (HCV) is one of the leading causes of liver diseases. Several host factors that facilitate the attachment and entry of HCV have been discovered, of which human occludin seems to be the most promising. Studies have shown that activity of occludin is dependent upon its phosphorylation status, and that during HCV infection deregulation of phosphorylated occludin collectively leads to a reduction in tight junction (TJ) integrity of hepatocytes and favors HCV entry. However, detailed information of the posttranslational modifications (PTMs) of occludin still remains largely unknown. In addition to phosphorylation, serine/threonine residues of several proteins are also regulated by a unique type of modification known as O-ß-glycosylation and this crosstalk serves as a functional switch. To identify the O-ß-glycosylation potential and how interplay between phosphorylation and O-ß-glycosylation can be exploited for the inhibition of HCV entry, here we report a computational analysis of PTMs of human occludin. Several conserved phosphorylation residues and kinases that can alter the ability of occludin to regulate the integrity of TJs were identified. In addition to previously reported Tyr residues, two additional Tyr residues (Tyr29 and Tyr287) were identified as target sites of Src kinase. To our knowledge, this is the first study to report the O-ß-GlcNAc potential of occludin and target sites of ERK (Ser8, Ser310, and Thr345), GSK-3 (Ser8, Ser341) and Cdk5 (Thr376). Furthermore, based on findings from this study, a potential novel interplay between phosphorylation and O-ß-glycosylation at the two Yin Yang sites (Ser408 and Ser490) is also proposed.

Comparative sequence, antigenic and phylogenetic analysis of avian influenza (H9N2) surface proteins isolated in Pakistan between 1999 and 2008.

INTRODUCTION: Influenza A viruses possess a unique genomic structure which leads to genetic instability, especially in products of neuraminidase and hemagglutinin genes. These surface proteins play major roles in viral entry and release, and in the activation of the host immune system. METHODOLOGY: This study involved an in silico sequence, phylogenetic and antigenic analyses of hemagglutinin and neuraminidase proteins of avian influenza A (H9N2) strains that circulated in Pakistan's poultry flocks from 1999 to 2008 and determined variations among these sequences at different levels. RESULTS: Sequence and phylogenetic analysis revealed a large number of similar substitution mutations and close evolutionary relation among sequences of both proteins. Changes were observed in the N-glycosylation sites of both surface proteins, along with the appearance of a new glycosylation site in the neuraminidase sequence isolated in 2007. Epitopes for hemagglutinin remained conserved, whereas for neuraminidase, epitopes from older strains reappeared in present sequences. CONCLUSIONS: Because of the rapid mutating nature of avian influenza subtype H9N2, constant surveillance of annual sequence variations is important. Preventive measures and vaccine products can be evaluated by keeping track of changes that may lead to reassortment among different circulating strains in Pakistan's commercial poultry flocks or in humans.