Influenza B viruses with site-specific mutations introduced into the HA gene.

We have succeeded in engineering changes into the genome of influenza B virus. First, model RNAs containing the chloramphenicol acetyltransferase gene flanked by the noncoding sequences of the HA or NS genes of influenza B virus were transfected into cells which were previously infected with an influenza B helper virus. Like those of the influenza A viruses, the termini of influenza B virus genes contain cis-acting signals which are sufficient to direct replication, expression, and packaging of the RNA. Next, a full-length copy of the HA gene from influenza B/Maryland/59 virus was cloned. Following transfection of this RNA, we rescued transfectant influenza B viruses which contain a point mutation introduced into the original cDNA. A series of mutants which bear deletions or changes in the 5' noncoding region of the influenza B/Maryland/59 virus HA gene were constructed. We were able to rescue viruses which contained deletions of 10 or 33 nucleotides at the 5' noncoding region of the HA gene. The viability of these viruses implies that this region of the genome is flexible in sequence and length.     

A/swine/H1N1流感病毒中国分离株血凝素进化分析

2009年3月以来,甲型H1N1病毒已在全球包括中国造成了巨大的危害,所以利用生物信息技术对其进行研究显得十分必要。从NCBI数据库下载中国大陆境内A/swine/H1N1病毒HA基因核酸序列及其编码蛋白序列,利用MEGA4.0软件对核苷酸编码序列构建系统进化树,利用BioEdit软件对蛋白序列进行比对,分析重要抗原位点变异情况,结果显示2010年在广东流行的病毒,从其他地方传播到广东的,而非早期广东流行的的病毒变异而来。2008年福建,山东,北京等地区的病毒传播比较迅速.这些分析结果阐明了中国大陆境内A/swine/H1N1病毒血凝素(HA)基因的进化关系和变异趋向,对于研究A/swine/H1N1病毒具有重要的参考价值。     

Insights from modeling protein evolution with context-dependent mutation and asymmetric amino acid selection

We develop an approximate maximum likelihood method to estimate flanking nucleotide context-dependent mutation rates and amino acid exchange-dependent selection in orthologous protein-coding sequences and use it to analyze genome-wide coding sequence alignments from mammals and yeast. Allowing context-dependent mutation provides a better fit to coding sequence data than simpler (context-independent or CpG "hotspot") models and significantly affects selection parameter estimates. Allowing asymmetric (nonreciprocal) selection on amino acid exchanges gives a better fit than simple dN/dS or symmetric selection models. Relative selection strength estimates from our models show good agreement with independent estimates derived from human disease-causing and engineered mutations. Selection strengths depend on local protein structure, showing expected biophysical trends in helical versus nonhelical regions and increased asymmetry on polar-hydrophobic exchanges with increased burial. The more stringent selection that has previously been observed for highly expressed proteins is primarily concentrated in buried regions, supporting the notion that such proteins are under stronger than average selection for stability. Our analyses indicate that a highly parameterized model of mutation and selection is computationally tractable and is a useful tool for exploring a variety of biological questions concerning protein and coding sequence evolution.     

Detecting overlapping coding sequences with pairwise alignments

MOTIVATION: Overlapping gene coding sequences (CDSs) are particularly common in viruses but also occur in more complex genomes. Detecting such genes with conventional gene-finding algorithms can be difficult for several reasons. If an overlapping CDS is on the same read-strand as a known CDS, then there may not be a distinct promoter or mRNA. Furthermore, the constraints imposed by double-coding can result in atypical codon biases. However, these same constraints lead to particular mutation patterns that may be detectable in sequence alignments. RESULTS: In this paper, we investigate several statistics for detecting double-coding sequences with pairwise alignments--including a new maximum-likelihood method. We also develop a model for double-coding sequence evolution. Using simulated sequences generated with the model, we characterize the distribution of each statistic as a function of sequence composition, length, divergence time and double-coding frame. Using these results, we develop several algorithms for detecting overlapping CDSs. The algorithms were tested on known overlapping CDSs and other overlapping open reading frames (ORFs) in the hepatitis B virus (HBV), Escherichia coli and Salmonella typhimurium genomes. The algorithms should prove useful for detecting novel overlapping genes--especially short coding ORFs in viruses. AVAILABILITY: Programs may be obtained from the authors. SUPPLEMENTARY INFORMATION: http://biochem.otago.ac.nz/double.html.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.     

Annotation of selection strengths in viral genomes

MOTIVATION: Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses. RESULTS: We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.     

Polymerase errors accumulating during natural evolution of the glycoprotein gene of vesicular stomatitis virus Indiana serotype isolates

We report the entire glycoprotein (G) gene nucleotide sequences of 26 vesicular stomatitis virus Indiana serotype (VSV IND) type 1 isolates from North and Central America. These sequences are also compared with partial G gene sequences of VSV IND type 2 (Cocal) and type 3 (Alagoas) viruses and the complete G gene sequences of the more distantly related VSV New Jersey (NJ) and Chandipura viruses. Phylogenetic analysis of the G gene sequences by maximum parsimony revealed four major lineages or subtypes within the classical VSV IND (type 1) viruses, each with a distinct geographic distribution. A high degree of VSV genetic diversity was found in Central America, with several virus subtypes of both VSV IND and NJ serotypes existing in this mainly enzootic disease region. Nineteen percent sequence variation but no deletions or insertions were evident within the 5' noncoding and the coding regions of the VSV IND type 1 G genes. In addition to numerous base substitutions, the 3' noncoding regions of these viruses also contained numerous base insertions and deletions. This resulted in striking variation in G gene sizes, with gene lengths ranging from 1,652 to 1,868 nucleotides. As the VSV IND type 1 subtypes have diverged from the common ancestor with the NJ subtypes, their G mRNAs have accumulated more 3' noncoding sequence inserts, ranging up to 303 nucleotides in length. These primarily consist of an imprecise reiteration of the sequence UUUUUAA, apparently generated by a unique polymerase stuttering error. Analysis of the deduced amino acid sequence differences among VSV IND type 1 viruses revealed numerous substitutions within defined antigenic epitopes, suggesting that immune selection may play a role in the evolution of these viruses.     

基于Web服务的流感病毒基因组自动化翻译注释系统

随着流感病毒基因组测序数据的急剧增加,深入挖掘流感病毒基因组大数据蕴含的生物学信息成为研究热点。基于中国流感病毒流行特征数据,建设一个集自动化、一体化和信息化的序列库系统,对于实现流感病毒基因组批量快速翻译、注释、存储、查询、分析具有重要的应用价值。本课题组通过集成一系列软件和工具包,并结合自主研发的其他功能,在底层维护的2个关键的参考数据集基础上另外追加了翻译注释信息最佳匹配的精细化筛选规则,构建具有流感病毒基因组信息存储、自动化翻译、蛋白序列精准注释、同源序列比对和进化树分析等功能的自动化系统。结果显示,通过Web端输入fasta格式的流感病毒基因序列,本系统可针对参考序列片段数据集(blastdb.fasta)进行Blast同源性检索,可以鉴定流感病毒的型别(A、B或C)、亚型和基因片段(1~8片段);在此基础上,通过查询数据库底层用于翻译、注释的基因片段参考数据集,可以获得一组肽段数据集,然后通过循环调用ProSplign软件对其进行预测。结合精细化的筛选准入规则,选出与输入序列匹配最好的翻译后产物,作为该输入序列的预测蛋白,输出为gbk,asn和fasta等通用格式的文件,给出序列长度、是否全长、病毒型别、亚型、片段等信息。基于以上工作,另外自主研发了系统其他的附加功能如进化树分析展示、基因组数据存储等功能,构建成基于Web服务的流感病毒基因组自动化翻译注释系统。本研究提示,系统高度集成系列软件以及自有的注释翻译数据库文件,实现从序列存储、翻译、注释到序列分析和展示的功能,可全面满足我国高通量基因检测数据共享化、本土化、一体化、自动化的需求。     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsynonymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in viruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected positions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

Genome sequencing and phylogenetic analysis of three avian influenza H9N2 subtypes in Guangxi

Three isolates of H9N2 Avian Influenza viruses (AIV) were isolated from chickens in Guangxi province. Eight pairs of specific primers were designed and synthesized according to the sequences of H9N2 at GenBank. phylogenetic analysis showed a high degree of homology between the Guangxi isolates and isolates from Guangdong and Jiangsu provinces, suggesting that the Guangxi isolates originated from the same source. However, the eight genes of the three isolates from Guangxi were not in the same sublineages in their respective phylogenetic trees, which suggests that they were products of natural reassortment between H9N2 avian influenza viruses from different sublineages. The 9 nucleotides ACAGAGATA which encode amino acids T, G, I were absent between nucleotide 205 and 214 in the open reading frame of the NA gene in the Guangxi isolates. AIV strains that infect human have, in their HA proteins, leucine at position 226. The analysis of deduced amino acid sequence of HA proteins showed that position 226 of these isolates contained glycine instead of leucine, suggesting that these three isolates differ from H9N2 AIV strains isolated from human infections.     

Polymorphism and signature of selection in the MHC class I genes of the three-spined stickleback Gasterosteus aculeatus

The role and intensity of positive selection maintaining the polymorphism of major histocompatibility complex (MHC) class I genes in the three-spined stickleback Gasterosteus aculeatus was investigated. The highly polymorphic set of MHC class I genes found was organized in a single linkage group. Between 5 and 14 sequence variants per individual were identified by single-stranded conformation polymorphism (SSCP) analysis. Segregation analysis studied in 10 three-spined stickleback families followed the expected pattern of Mendelian inheritance. The gamete fusion in three-spined stickleback thus seems to be random with respect to the MHC class I genes. The DNA sequence analyses showed that the expressed MHC class I loci are under strong selection pressure, possibly mediated by parasites. Codons that were revealed to be under positive selection are potentially important in antigen binding. MHC class I sequences did not form significant supported clusters within a phylogenetic tree. Analogous to MHC class II genes, it was not possible to assign the class I sequences to a specific locus, suggesting that the class I genes may have been generated by recent gene duplication.     

The variable codons of H5N1 avian influenza A virus haemagglutinin genes

We investigated the selection pressures on the haemagglutinin genes of H5N1 avian influenza viruses using fixed effects likelihood models. We found evidence of positive selection in the sequences from isolates from 1997 to 2007, except viruses from 2000. The haemagglutinin sequences of viruses from southeast Asia, Hong Kong and mainland China were the most polymorphic and had similar nonsyn-onymous profiles. Some sites were positively selected in viruses from most regions and a few of these sites displayed different amino acid patterns. Selection appeared to produce different outcomes in vi-ruses from Europe, Africa and Russia and from different host types. One position was found to be positively selected for human isolates only. Although the functions of some positively selected posi-tions are unknown, our analysis provided evidence of different temporal, spatial and host adaptations for H5N1 avian influenza viruses.     

Transfection-mediated recombination of influenza A virus.

Several mechanisms, including a high mutation rate and reassortment of genes, have been found to be responsible for the variability of influenza A viruses. RNA recombination would be another mechanism leading to genetic variation; however, recombination has only rarely been reported to occur in influenza viruses. During ribonucleoprotein transfection experiments designed to generate viable influenza viruses from in vitro-synthesized RNA, we discovered several viruses which must have originated from recombination events. The ribonucleoprotein transfection system may enhance the formation of viruses which result from jumping of the viral polymerase between RNAs or from ligation of different viral RNAs. Five different recombinant viruses are described. Two of these, REC1 and REC2, contain a neuraminidase (NA) gene whose defective polyadenylation signal has been repaired via intergenic recombination; 124 and 95 nucleotides have been added, respectively. Another virus, REC5, must have originated by multiple recombination events since it contains a mosaic gene with sequences derived from the NA gene of influenza A/WSN/33 virus and the matrix, polymerase protein PB1, and NA genes of influenza A/PR/8/34 virus.     

Complete nucleotide sequence of the influenza A/PR/8/34 virus NS gene and comparison with the NS genes of the A/Udorn/72 and A/FPV/Rostock/34 strains.

The nucleotide sequence of the NS gene of the human influenza virus A/PR/8/34 was determined and found to be the same length (890 nucleotides) as the NS gene of another human influenza virus A/Udorn/72 and of the avian isolate A/FPV/Rostock/34. Comparison of the sequences of the NS genes of the two human influenza viruses shows an 8.9% difference whereas the NS gene of the avian isolate differs by only 8% from that of the human strain A/PR/8/34. The extensive sequence similarity among these three genes does not support the notion of species specific homology groups among NS genes of avian and human influenza virus strains. The primary sequence of the A/PR/8/34 NS gene is consistent with the findings that the influenza virus NS gene may code for two overlapping polypeptides. In addition, an open reading frame potentially coding for a polypeptide 167 amino acids in length was found in the negative strand RNA of the A/PR/8/34 virus NS gene.     

Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes

Phylogenetic profiles of the genes coding for the hemagglutinin (HA) protein, nucleoprotein (NP), matrix (M) protein, and nonstructural (NS) proteins of influenza B viruses isolated from 1940 to 1998 were analyzed in a parallel manner in order to understand the evolutionary mechanisms of these viruses. Unlike human influenza A (H3N2) viruses, the evolutionary pathways of all four genes of recent influenza B viruses revealed similar patterns of genetic divergence into two major lineages. Although evolutionary rates of the HA, NP, M, and NS genes of influenza B viruses were estimated to be generally lower than those of human influenza A viruses, genes of influenza B viruses demonstrated complex phylogenetic patterns, indicating alternative mechanisms for generation of virus variability. Topologies of the evolutionary trees of each gene were determined to be quite distinct from one another, showing that these genes were evolving in an independent manner. Furthermore, variable topologies were apparently the result of frequent genetic exchange among cocirculating epidemic viruses. Evolutionary analysis done in the present study provided further evidence for cocirculation of multiple lineages as well as sequestering and reemergence of phylogenetic lineages of the internal genes. In addition, comparison of deduced amino acid sequences revealed a novel amino acid deletion in the HA1 domain of the HA protein of recent isolates from 1998 belonging to the B/Yamagata/16/88-like lineage. It thus became apparent that, despite lower evolutionary rates, influenza B viruses were able to generate genetic diversity among circulating viruses through a combination of evolutionary mechanisms involving cocirculating lineages and genetic reassortment by which new variants with distinct gene constellations emerged.     

Inferring parameters of mutation, selection and demography from patterns of synonymous site evolution in Drosophila

Selection acting on codon usage can cause patterns of synonymous evolution to deviate considerably from those expected under neutrality. To investigate the quantitative relationship between parameters of mutation, selection, and demography, and patterns of synonymous site divergence, we have developed a novel combination of population genetic models and likelihood methods of phylogenetic sequence analysis. Comparing 50 orthologous gene pairs from Drosophila melanogaster and D. virilis and 27 from D. melanogaster and D. simulans, we show considerable variation between amino acids and genes in the strength of selection acting on codon usage and find evidence for both long-term and short-term changes in the strength of selection between species. Remarkably, D. melanogaster shows no evidence of current selection on codon usage, while its sister species D. simulans experiences only half the selection pressure for codon usage of their common ancestor. We also find evidence for considerable base asymmetries in the rate of mutation, such that the average synonymous mutation rate is 20-30% higher than in noncoding regions. A Bayesian approach is adopted to investigate how accounting for selection on codon usage influences estimates of the parameters of mutation.     

Integrating historical, clinical and molecular genetic data in order to explain the origin and virulence of the 1918 Spanish influenza virus.

The Spanish influenza pandemic of 1918-1919 caused acute illness in 25-30% of the world's population and resulted in the death of 40 million people. The complete genomic sequence of the 1918 influenza virus will be deduced using fixed and frozen tissues of 1918 influenza victims. Sequence and phylogenetic analyses of the complete 1918 haemagglutinin (HA) and neuraminidase (NA) genes show them to be the most avian-like of mammalian sequences and support the hypothesis that the pandemic virus contained surface protein-encoding genes derived from an avian influenza strain and that the 1918 virus is very similar to the common ancestor of human and classical swine H1N1 influenza strains. Neither the 1918 HA genes nor the NA genes possessed mutations that are known to increase tissue tropicity, which accounts for the virulence of other influenza strains such as A/WSN/33 or fowl plague viruses. The complete sequence of the nonstructural (NS) gene segment of the 1918 virus was deduced and tested for the hypothesis that the enhanced virulence in 1918 could have been due to type I interferon inhibition by the NS1 protein. The results from these experiments were inconclusive. Sequence analysis of the 1918 pandemic influenza virus is allowing us to test hypotheses as to the origin and virulence of this strain. This information should help to elucidate how pandemic influenza strains emerge and what genetic features contribute to their virulence.     

Phylogenetic Analysis of the Entire Genome of Influenza A (H3N2) Viruses from Japan: Evidence for Genetic Reassortment of the Six Internal Genes

Nucleotide sequences of all eight RNA segments of 10 human H3N2 influenza viruses isolated during a 5-year period from 1993 to 1997 were determined and analyzed phylogenetically in order to define the evolutionary pathways of all genes in a parallel fashion. It was evident that the hemagglutinin and neuraminidase genes of these viruses evolved essentially in a single lineage and that amino acid changes accumulated sequentially with respect to time. In contrast, amino acid differences in the internal proteins were erratic and did not accumulate over time. Parallel analysis of the phylogenetic patterns of all genes revealed that the evolutionary pathways of the six internal genes were not linked to the surface glycoproteins. Genes coding for the basic polymerase-1, nucleoprotein, and matrix proteins of 1997 isolates were closest phylogenetically to those of earlier isolates of 1993 and 1994. Furthermore, all six internal genes of four viruses isolated in the 1995 epidemic season consistently divided into two distinct branch clusters, and two 1995 isolates contained PB2 genes apparently originating from those of viruses before 1993. It was apparent that the lack of correlation between the topologies of the phylogenetic trees of the genes coding for the surface glycoproteins and internal proteins was a reflection of genetic reassortment among human H3N2 viruses. This is the first evidence demonstrating the occurrence of genetic reassortment involving the internal genes of human H3N2 viruses. Furthermore, internal protein variability coincided with marked increases in the activity of H3N2 viruses in 1995 and 1997.     

Evolution of the H3 influenza virus hemagglutinin from human and nonhuman hosts.

The nucleotide and amino acid sequences of 40 influenza virus hemagglutinin genes of the H3 serotype from mammalian and avian species and 9 genes of the H4 serotype were compared, and their evolutionary relationships were evaluated. From these relationships, the differences in the mutational characteristics of the viral hemagglutinin in different hosts were examined and the RNA sequence changes that occurred during the generation of the progenitor of the 1968 human pandemic strain were examined. Three major lineages were defined: one containing only equine virus isolates; one containing only avian virus isolates; and one containing avian, swine, and human virus isolates. The human pandemic strain of 1968 was derived from an avian virus most similar to those isolated from ducks in Asia, and the transfer of this virus to humans probably occurred in 1965. Since then, the human viruses have diverged from this progenitor, with the accumulation of approximately 7.9 nucleotide and 3.4 amino acid substitutions per year. Reconstruction of the sequence of the hypothetical ancestral strain at the avian-human transition indicated that only 6 amino acids in the mature hemagglutinin molecule were changed during the transition between an avian virus strain and a human pandemic strain. All of these changes are located in regions of the molecule known to affect receptor binding and antigenicity. Unlike the human H3 influenza virus strains, the equine virus isolates have no close relatives in other species and appear to have diverged from the avian viruses much earlier than did the human virus strains. Mutations were estimated to have accumulated in the equine virus lineage at approximately 3.1 nucleotides and 0.8 amino acids per year. Four swine virus isolates in the analysis each appeared to have been introduced into pigs independently, with two derived from human viruses and two from avian viruses. A comparison of the coding and noncoding mutations in the mammalian and avian lineages showed a significantly lower ratio of coding to total nucleotide changes in the avian viruses. Additionally, the avian virus lineages of both the H3 and H4 serotypes, but not the mammalian virus lineages, showed significantly greater conservation of amino acid sequence in the internal branches of the phylogenetic tree than in the terminal branches. The small number of amino acid differences between the avian viruses and the progenitor of the 1968 pandemic strain and the great phenotypic stability of the avian viruses suggest that strains similar to the progenitor strain will continue to circulate in birds and will be available for reintroduction into humans.     

Antigenic drift in the influenza A virus (H3N2) nucleoprotein and escape from recognition by cytotoxic T lymphocytes

Viruses exploit different strategies to escape immune surveillance, including the introduction of mutations in cytotoxic T-lymphocyte (CTL) epitopes. The sequence of these epitopes is critical for their binding to major histocompatibility complex (MHC) class I molecules and recognition by specific CTLs, both of which interactions may be lost by mutation. Sequence analysis of the nucleoprotein gene of influenza A viruses (H3N2) isolated in The Netherlands from 1989 to 1999 revealed two independent amino acid mutations at the anchor residue of the HLA-B27-specific CTL epitope SRYWAIRTR (383 to 391). A R384K mutation was found in influenza A viruses isolated during the influenza season 1989-1990 but not in subsequent seasons. In the influenza season 1993-1994, a novel mutation in the same CTL epitope at the same position was introduced. This R384G mutation proved to be conserved in all influenza A viruses isolated from 1993 onwards. Both mutations R384K and R384G abrogated MHC class I presentation and allowed escape from recognition by specific CTLs.