Influenza Sequences Research Articles

MOLECULAR CHARACTERIZATION ON THE BASIS OF HA, NA AND M GENE REVEALED CHANGES IN CRITICAL AMINO ACID POSITIONS OF INFLUENZA A (H3N2) VIRUS CIRCULATING IN INDIA DURING 2011-2013

Influenza a (H1N1) virus is responsible for acute respiratory infection in human, which occasionally causes epidemic and rarely a pandemic. Full length DNA sequencing of HA, NA and M gene of H3N2 virus from year 2011-2013 was performed for determination of the circulating strains in India along with monitoring status of various mutations associated with drug resistance and virulence. Genetic analysis of HA gene of study samples in comparison to reference strain showed maximum amino acid changes in epitopes of high neutralization efficiencies (antigenic sites A, B and D) along with gain of glycosylation sites at position 45 (except single sample) and loss at position 126 (7 samples). NA gene has four amino acid changes in antibody binding region at position S367N, K369T, S332F and S334I along with gain and loss of glycosylation site at position 367 and 402 respectively whereas no change was observed in its active site. While comparison of deduced amino acid sequences of study samples with Indian strains from earlier years showed six amino acid changes in three antigenic sites on HA protein and three amino acid changes in NA protein epitope sequence. All these changes indicate that the H3N2 virus is undergoing antigenic drift and which may result in the emergence of new antigenic variants. We want to report that our study provided for the first time full length matrix gene sequence of influenza a (H3N2) virus from India.

Sequencing artifacts in the type A influenza databases and attempts to correct them

BackgroundThere are over 276 000 influenza gene sequences in public databases, with the quality of the sequences determined by the contributor.ObjectiveAs part of a high school class project, influenza sequences with possible errors were identified in the public databases based on the size of the gene being longer than expected, with the hypothesis that these sequences would have an error. Students contacted sequence submitters alerting them of the possible sequence issue(s) and requested they the suspect sequence(s) be correct as appropriate.MethodsType A influenza viruses were screened, and gene segments longer than the accepted size were identified for further analysis. Attention was placed on sequences with additional nucleotides upstream or downstream of the highly conserved non-coding ends of the viral segments.Results and ConclusionsA total of 1081 sequences were identified that met this criterion. Three types of errors were commonly observed: non-influenza primer sequence wasn't removed from the sequence; PCR product was cloned and plasmid sequence was included in the sequence; and Taq polymerase added an adenine at the end of the PCR product. Internal insertions of nucleotide sequence were also commonly observed, but in many cases it was unclear if the sequence was correct or actually contained an error. A total of 215 sequences, or 22.8% of the suspect sequences, were corrected in the public databases in the first year of the student project. Unfortunately 138 additional sequences with possible errors were added to the databases in the second year. Additional awareness of the need for data integrity of sequences submitted to public databases is needed to fully reap the benefits of these large data sets.

Full-Genome Sequences of Influenza A(H1N1)pdm09 Viruses Isolated from Finnish Patients from 2009 to 2013

Here we report full-length sequencing of the first large set of influenza A(H1N1)pdm09 virus genomes isolated in Finland between the years 2009 and 2013 and discuss the advantages and needs of influenza virus sequencing efforts.

Deep sequencing approach for genetic stability evaluation of influenza A viruses

Assessment of genetic stability of viruses could be used to monitor manufacturing process of both live and inactivated viral vaccines. Until recently such studies were limited by the difficulty of detecting and quantifying mutations in heterogeneous viral populations. High-throughput sequencing technologies (deep sequencing) can generate massive amounts of genetic information and could be used to reveal and quantify mutations. Comparison of different approaches for deep sequencing of the complete influenza A genome was performed to determine the best way to detect and quantify mutants in attenuated influenza reassortant strain A/Brisbane/59/2007 (H1N1) and its passages in different cell substrates. Full-length amplicons of influenza A virus segments as well as multiple overlapping amplicons covering the entire viral genome were subjected to several ways of DNA library preparation followed by deep sequencing using Solexa (Illumina) and pyrosequencing (454 Life Science) technologies. Sequencing coverage (the number of times each nucleotide was determined) of mutational profiles generated after 454-pyrosequencing of individually synthesized overlapping amplicons were relatively low and insufficiently uniform. Amplification of the entire genome of influenza virus followed by its enzymatic fragmentation, library construction, and Illumina sequencing resulted in high and uniform sequencing coverage enabling sensitive quantitation of mutations. A new bioinformatic procedure was developed to improve the post-alignment quality control for deep-sequencing data analysis.

Comparison of Mutation Patterns in Full-Genome A/H3N2 Influenza Sequences Obtained Directly from Clinical Samples and the Same Samples after a Single MDCK Passage

Human influenza viruses can be isolated efficiently from clinical samples using Madin-Darby canine kidney (MDCK) cells. However, this process is known to induce mutations in the virus as it adapts to this non-human cell-line. We performed a systematic study to record the pattern of MDCK-induced mutations observed across the whole influenza A/H3N2 genome. Seventy-seven clinical samples collected from 2009-2011 were included in the study. Two full influenza genomes were obtained for each sample: one from virus obtained directly from the clinical sample and one from the matching isolate cultured in MDCK cells. Comparison of the full-genome sequences obtained from each of these sources showed that 42% of the 77 isolates had acquired at least one MDCK-induced mutation. The presence or absence of these mutations was independent of viral load or sample origin (in-patients versus out-patients). Notably, all the five hemagglutinin missense mutations were observed at the hemaggutinin 1 domain only, particularly within or proximal to the receptor binding sites and antigenic site of the virus. Furthermore, 23% of the 77 isolates had undergone a MDCK-induced missense mutation, D151G/N, in the neuraminidase segment. This mutation has been found to be associated with reduced drug sensitivity towards the neuraminidase inhibitors and increased viral receptor binding efficiency to host cells. In contrast, none of the neuraminidase sequences obtained directly from the clinical samples contained the D151G/N mutation, suggesting that this mutation may be an indicator of MDCK culture-induced changes. These D151 mutations can confound the interpretation of the hemagglutination inhibition assay and neuraminidase inhibitor resistance results when these are based on MDCK isolates. Such isolates are currently in routine use in the WHO influenza vaccine and drug-resistance surveillance programs. Potential data interpretation miscalls can therefore be avoided by careful exclusion of such D151 mutants after further sequence analysis.

Rapid sequencing of influenza A virus vRNA, cRNA and mRNA non-coding regions

Characterizing the genomic sequences of influenza A viruses is important for pathophysiological and evolutionary studies. Noncoding regions (NCR) of influenza A virus have been shown to play critical roles in replication and transcription but their sequences are infrequently determined. In this study, a method employing poly(A) addition and SMART (switching mechanism at 5′ end of RNA transcript) technology is described for directly determining and discriminating both NCR ends of viral RNA (vRNA), complementary RNA (cRNA), or NCR and cap sequences from viral mRNA. This modified method may also be used to characterize the NCRs of influenza A virus samples in which the RNA has been degraded.

Differing Epidemiological Dynamics of Influenza B Virus Lineages in Guangzhou, Southern China, 2009-2010

The epidemiological and evolutionary dynamics of the two cocirculating lineages of influenza B virus, Victoria and Yamagata, are poorly understood, especially in tropical or subtropical areas of Southeast Asia. We performed a phylogenetic analysis of the hemagglutinin (HA) and neuraminidase (NA) sequences of influenza B viruses isolated in Guangzhou, a southern Chinese city, during 2009 to 2010 and compared the demographic and clinical features of infected patients. We identified multiple viral introductions of Victoria strains from both Chinese and international sources, which formed two phylogenetically and antigenically distinct clades (Victoria 1 and 2), some of which persisted between seasons. We identified one dominant Yamagata introduction from outside China during 2009. Our phylogenetic analysis reveals the occurrence of reassortment events among the Victoria and Yamagata lineages and also within the Victoria lineage. We found no significant difference in clinical severity by influenza B lineage, with the exceptions that (i) the Yamagata lineage infected older people than either Victoria lineage and (ii) fewer upper respiratory tract infections were caused by the Victoria 2 than the Victoria 1 clade. Overall, our study reveals the complex epidemiological dynamics of different influenza B lineages within a single geographic locality and has implications for vaccination policy in southern China.

Novel bioinformatics strategies for prediction of directional sequence changes in influenza virus genomes and for surveillance of potentially hazardous strains

BackgroundWith the remarkable increase of microbial and viral sequence data obtained from high-throughput DNA sequencers, novel tools are needed for comprehensive analysis of the big sequence data. We have developed “Batch-Learning Self-Organizing Map (BLSOM)” which can characterize very many, even millions of, genomic sequences on one plane. Influenza virus is one of zoonotic viruses and shows clear host tropism. Important issues for bioinformatics studies of influenza viruses are prediction of genomic sequence changes in the near future and surveillance of potentially hazardous strains.MethodsTo characterize sequence changes in influenza virus genomes after invasion into humans from other animal hosts, we applied BLSOMs to analyses of mono-, di-, tri-, and tetranucleotide compositions in all genome sequences of influenza A and B viruses and found clear host-dependent clustering (self-organization) of the sequences.ResultsViruses isolated from humans and birds differed in mononucleotide composition from each other. In addition, host-dependent oligonucleotide compositions that could not be explained with the host-dependent mononucleotide composition were revealed by oligonucleotide BLSOMs. Retrospective time-dependent directional changes of mono- and oligonucleotide compositions, which were visualized for human strains on BLSOMs, could provide predictive information about sequence changes in newly invaded viruses from other animal hosts (e.g. the swine-derived pandemic H1N1/09).ConclusionsBasing on the host-dependent oligonucleotide composition, we proposed a strategy for prediction of directional changes of virus sequences and for surveillance of potentially hazardous strains when introduced into human populations from non-human sources. Millions of genomic sequences from infectious microbes and viruses have become available because of their medical and social importance, and BLSOM can characterize the big data and support efficient knowledge discovery.

Highly Pathogenic Avian Influenza A(H5N1) Mutants Transmissible by Air Are Susceptible to Human and Animal Neutralizing Antibodies

A laboratory-generated reassortant H5 hemagglutinin (HA)/influenza A(H1N1) strain containing 4 mutations in influenza A(H5N1) HA has become transmissible by air among mammals. Here, we constructed 15 influenza A(H5N1) pseudoviruses containing a single mutation or a combination of mutations and showed that the pseudoviruses were susceptible to neutralizing antibodies from patients with influenza A(H5N1) infection and from mice immunized with a vaccine containing the conserved HA1 sequence of influenza A(H5N1). These results indicate that antibodies in patients currently infected by influenza A(H5N1) and antibodies induced by vaccines containing conserved sequences in HA1 of wild-type influenza A(H5N1) are highly effective in cross-neutralizing future influenza A(H5N1) mutants with airborne transmissibility, suggesting that human influenza pandemics caused by these influenza A(H5N1) variants can be prevented.

Simultaneous and complete genome sequencing of influenza A and B with high coverage by Illumina MiSeq Platform

Active global surveillance and characterization of influenza viruses are essential for better preparation against possible pandemic events. Obtaining comprehensive information about the influenza genome can improve our understanding of the evolution of influenza viruses and emergence of new strains, and improve the accuracy when designing preventive vaccines. This study investigated the use of deep sequencing by the next-generation sequencing (NGS) Illumina MiSeq Platform to obtain complete genome sequence information from influenza virus isolates. The influenza virus isolates were cultured from 6 respiratory acute clinical specimens collected in Thailand and Nepal. DNA libraries obtained from each viral isolate were mixed and all were sequenced simultaneously. Total information of 2.6Gbases was obtained from a 455±14K/mm2 density with 95.76% (8,571,655/8,950,724 clusters) of the clusters passing quality control (QC) filters. Approximately 93.7% of all sequences from Read1 and 83.5% from Read2 contained high quality sequences that were ≥Q30, a base calling QC score standard. Alignments analysis identified three seasonal influenza A H3N2 strains, one 2009 pandemic influenza A H1N1 strain and two influenza B strains. The nearly entire genomes of all six virus isolates yielded equal or greater than 600-fold sequence coverage depth. MiSeq Platform identified seasonal influenza A H3N2, 2009 pandemic influenza A H1N1and influenza B in the DNA library mixtures efficiently.

Exploring the Molecular Epidemiology and Evolutionary Dynamics of Influenza A Virus in Taiwan

The evolution and population dynamics of human influenza in Taiwan is a microcosm of the viruses circulating worldwide, which has not yet been studied in detail. We collected 343 representative full genome sequences of human influenza A viruses isolated in Taiwan between 1979 and 2009. Phylogenetic and antigenic data analysis revealed that H1N1 and H3N2 viruses consistently co-circulated in Taiwan, although they were characterized by different temporal dynamics and degrees of genetic diversity. Moreover, influenza A viruses of both subtypes underwent internal gene reassortment involving all eight segments of the viral genome, some of which also occurred during non-epidemic periods. The patterns of gene reassortment were different in the two subtypes. The internal genes of H1N1 viruses moved as a unit, separately from the co-evolving HA and NA genes. On the other hand, the HA and NA genes of H3N2 viruses tended to segregate consistently with different sets of internal gene segments. In particular, as reassortment occurred, H3HA always segregated as a group with the PB1, PA and M genes, while N2NA consistently segregated with PB2 and NP. Finally, the analysis showed that new phylogenetic lineages and antigenic variants emerging in summer were likely to be the progenitors of the epidemic strains in the following season. The synchronized seasonal patterns and high genetic diversity of influenza A viruses observed in Taiwan make possible to capture the evolutionary dynamic and epidemiological rules governing antigenic drift and reassortment and may serve as a “warning” system that recapitulates the global epidemic.

Molecular and phylogenetic analysis of influenza A H1N1 pandemic viruses in Cuba, May 2009 to August 2010

The influenza A(H1N1)pdm09 virus was detected in Cuba in May 2009. The introduction of a new virus with increased transmissibility into a population makes surveillance of the pandemic strain to the molecular level necessary. The aim of the present study was the molecular and phylogenetic analysis of pandemic influenza A(H1N1)pdm09 strains that circulated in Cuba between May 2009 and August 2010. Seventy clinical samples were included in the study. Nucleotide sequences from the hemagglutinin HA1 region segment were obtained directly from clinical samples. Genetic distances were calculated using MEGA v.5.05. A phylogenetic tree was constructed using MrBayes v.3.1.2 software. Potential N-glycosylation sites were predicted using NetNGlyc server 1.0. The 48 Cuban sequences of influenza A(H1N1)pdm09 obtained were similar to the A/California/07/2009 (H1N1) vaccine strain. Most of the Cuban strains belonged to clade 7. Cuban viruses showed amino acid changes, some of them located at three antigenic sites: Ca, Sa, and Sb. Two dominant mutations were detected: P83S (100%) and S203T (85.7%). Glycosylation site analysis revealed the gain of one site at position 162 in 13 sequences. The findings in this study contribute to our understanding of the progress of the influenza A(H1N1)pdm09 virus, since this virus is at the starting point of its evolution in humans.

Impact of influenza virus escape‐mutations on influenza detection by the rapid influenza diagnostic test

During the influenza pandemic of 2009-2010, rapid influenza diagnostic tests (RIDTs) were used to detect influenza viral infections because they are quick and simple to use. However, retrospective studies showed that RIDTs performed poorly when used to diagnose pandemic viral infections. Determining how amino acid sequence changes in pandemic or epidemic influenza viral antigens impact clinical value of RIDTs has not been possible, because the viral epitopes recognized by RIDTs have been not mapped. In this study, the effect of escape-variations or mutations in influenza viral antigens upon the sensitivity and specificity of an RIDT was investigated by characterizing the immunological properties of the antibodies used in the RIDT. Escape-mutants were generated by cultivating A/Korea/01/2009 in the presence of an excess of the same antibodies used in the RIDT. Escape-mutants not recognized by the RIDT were selected. Epitopes recognized by the RIDT were mapped by comparing the sequence and immunological analysis of the escape-variants and wild-type isolates. The RIDT antibodies recognized epitopes on the Sa antigenic site and in the F subdomain in hemagglutinin. Variants bearing mutations in these epitopes were not detected by the RIDT. The frequency of escape-variants emerging since the 2009-2010 pandemic was calculated as 1.27% using in silico surveillance of influenza sequence databases. These results suggest that mapping the relevant epitopes of RIDTs and making such information available to clinics would be helpful for determining whether RIDTs match newly emergent strains and subtypes prior to retrospective re-evaluation of the RIDTs using clinical specimens.

Virus-like Particles Containing Multiple M2 Extracellular Domains Confer Improved Cross-protection Against Various Subtypes of Influenza Virus

The extracellular domain of M2 (M2e), a small ion channel membrane protein, is well conserved among different human influenza A virus strains. To improve the protective efficacy of M2e vaccines, we genetically engineered a tandem repeat of M2e epitope sequences (M2e5x) of human, swine, and avian origin influenza A viruses, which was expressed in a membrane-anchored form and incorporated in virus-like particles (VLPs). The M2e5x protein with the transmembrane domain of hemagglutinin (HA) was effectively incorporated into VLPs at a several 100-fold higher level than that on influenza virions. Intramuscular immunization with M2e5x VLP vaccines was highly effective in inducing M2e-specific antibodies reactive to different influenza viruses, mucosal and systemic immune responses, and cross-protection regardless of influenza virus subtypes in the absence of adjuvant. Importantly, immune sera were found to be sufficient for conferring protection in naive mice, which was long-lived and cross-protective. Thus, molecular designing and presenting M2e immunogens on VLPs provide a promising platform for developing universal influenza vaccines without using adjuvants.

Molecular surveillance of low pathogenic avian influenza viruses in wild birds across the United States: inferences from the hemagglutinin gene.

A United States interagency avian influenza surveillance plan was initiated in 2006 for early detection of highly pathogenic avian influenza viruses (HPAIV) in wild birds. The plan included a variety of wild bird sampling strategies including the testing of fecal samples from aquatic areas throughout the United States from April 2006 through December 2007. Although HPAIV was not detected through this surveillance effort we were able to obtain 759 fecal samples that were positive for low pathogenic avian influenza virus (LPAIV). We used 136 DNA sequences obtained from these samples along with samples from a public influenza sequence database for a phylogenetic assessment of hemagglutinin (HA) diversity in the United States. We analyzed sequences from all HA subtypes except H5, H7, H14 and H15 to examine genetic variation, exchange between Eurasia and North America, and geographic distribution of LPAIV in wild birds in the United States. This study confirms intercontinental exchange of some HA subtypes (including a newly documented H9 exchange event), as well as identifies subtypes that do not regularly experience intercontinental gene flow but have been circulating and evolving in North America for at least the past 20 years. These HA subtypes have high levels of genetic diversity with many lineages co-circulating within the wild birds of North America. The surveillance effort that provided these samples demonstrates that such efforts, albeit labor-intensive, provide important information about the ecology of LPAIV circulating in North America.

Incidence of respiratory viruses in a pediatric population: molecular and epidemiological aspects

Introduction: Respiratory infections are not well defined and the etiology is often unknown. Material and method: four hundred fortynine subjectrs were enrolled in the study; in all patientes there was a suspect of inflammatory diseases of the respiratory tract. At admission, a nasopharyngeal swab was made. A multiplex PCR was performed after extraction and reverse transcription of viral RNA. The amplified fragments were revealed by using an electrophoresis separation. Results: Two hundred and four patients (45.4%) were hospitalized for infection of the upper respiratory tract, 141 (31.4%) for lower respiratory infection and the remaining (23%) for other symptoms. One hundred fiftyseven (35%) patients were positive for human influenza A (H1N1 subtype) and 184 for other respiratory viruses,of which 59 (32%) gave a positive for respiratory syncytial virus, 42 (23%) for rhinovirus, 31 (17%) for parainfluenza virus, 12 (6.5%) for coronavirus, 28 (15%) for adenovirus and 6 (3%) for influenza B (3%) and 6 (3%) for metapneumovirus. The M1 gene sequence of influenza A H1N1 strains from 12 patients had a high identity with that of the reference virus. Conclusion: Furthermore H1N1 and RSV were the main causative agents of acute respiratory infection. A molecular approach provides an accurate and rapid aetiological diagnosis of viral respiratory infections. The molecolar features in the M1 gene suggested that the H1N1 influenza strains circulating in Apulia region had a conserved genetic make up.

Complete Genome Sequence of an Avian-Origin H3N2 Canine Influenza A Virus Isolated in Farmed Dogs in Southern China

We report here the complete genomic sequence of an avian-origin H3N2 canine influenza A virus containing multiple mutations in farmed dogs in southern China. Phylogenetic analyses of the sequences of all eight viral RNA segments demonstrated that these are wholly avian influenza viruses of the Asia lineage. To our knowledge, this is the first report of interspecies transmission of an avian H3N2 influenza virus to domestic farm dogs under natural conditions in Southern China. The amino acid information provided herein suggests that continued study is required to determine if this virus could be established in the farm dog population and pose potential threats to public health.

A novel empirical mutual information approach to identify co-evolving amino acid positions of influenza A viruses

Mutual information (MI) is an approach commonly used to estimate the evolutionary correlation of 2 amino acid sites. Although several MI methods exist, prior to our contribution no systematic method had been developed to assess their performance, or to establish numerical thresholds to detect co-evolving amino acid sites. The current study performed a Markov chain Monte Carlo (MCMC) algorithm on influenza viral sequences to capture their evolutionary characteristics. A consensus maximum clade credibility (MCC) tree was estimated from the samples, together with their amino acid substitution statistics, from which we generated synthetic sequences of known dependent and independent paired amino acid sites. A pair-to-pair and influenza-specific amino acid substitution matrix (P2PFLU) incorporated into Bayesian Evolutionary Analysis Sampling Trees (BEAST) enumerated these synthetic sequences. The sequences inherited evolutionary features and co-varying characteristics from the real viral sequences, rendering these synthetic data ideal for exploring their co-evolving features. For the MI measure, we proposed a novel metric called the empirical MI (MIEm), which outperformed other MI measures in analysis of receiver operating characteristics (ROC). We implemented our approach on 1086 all-time PB2 sequences of influenza A H5N1 viruses, in which we found 97 sites exhibiting co-evolutionary substitution of one or more amino acid sites. In particular, PB2 451, along with eight other PB2 sites of various MIEm scores, was found to co-evolve with PB2 627, a known species-associated amino acid residue which plays a critical role in influenza virus replication.

Deep sequencing of H1N1pdm Influenza A clinical samples from Luxembourg 2010/2011

Deep sequencing of H1N1pdm Influenza A clinical samples from Luxembourg 2010/2011

Receptor Binding Surveillance of Influenza Clinical Isolates

Aims: The outbreak of highly pathogenic avian influenza H5N1 in poultry and the influenza A pandemic have wreaked havoc on public health. Viruses alter their carbohydrate binding preferences, thereby causing pandemics all over the world. In this study, we tried to investigate the carbohydrate binding specificity of influenza clinical isolates. Materials & methods: Biotin-conjugated polyacrylamide-based glycan epitopes were immobilized on UltraBind™ membranes and used to survey the glycan-binding preference of influenza clinical isolates, including seasonal influenza A, A(H1N1)pdm09 and influenza B viruses. In addition, the DNA sequences of influenza B virus hemagglutinin were analyzed. Results: Human influenza A, especially the A(H1N1)pdm09 viruses, accepted α2,6 and α2,3 sialylated glycans, sulfated glycans and α2,8 sialosides. Although all influenza B clinical isolates bound strongly to NeuAc, 6´-sialyl lactose and sialyl biantennary N-glycan, some viruses also recognized sulfated and α2,3 sialylated glycans. According to the nucleotide sequences of viral hemagglutinin, influenza B viruses that exhibited weak interaction with sulfated and α2,3 sialylated glycans showed fewer charged amino acids. Conclusion: The substrate specificities of influenza clinical isolates were surveyed. Influenza A exhibited more complicated glycan-binding patterns than influenza B viruses. Our findings provided a systematic investigation of receptor-binding specificities for influenza clinical isolates, as well as useful information for exploring viral tropism.