A complete map of potential pathogenicity markers of avian influenza virus subtype H5 predicted from 11 expressed proteins

Predictive analysis for pathogenicity classification of H5Nx avian influenza strains using machine learning techniques

Over the past decades, avian influenza (AI) outbreaks have been reported across different parts of the globe, resulting in large-scale economic and livestock loss and, in some cases raising concerns about their zoonotic potential. The virulence and pathogenicity of H5Nx (e.g., H5N1, H5N2) AI strains for poultry could be inferred through various approaches, and it has been frequently performed by detecting certain pathogenicity markers in their haemagglutinin (HA) gene. The utilization of predictive modeling methods represents a possible approach to exploring this genotypic-phenotypic relationship for assisting experts in determining the pathogenicity of circulating AI viruses. Therefore, the main objective of this study was to evaluate the predictive performance of different machine learning (ML) techniques for in-silico prediction of pathogenicity of H5Nx viruses in poultry, using complete genetic sequences of the HA gene. We annotated 2137 H5Nx HA gene sequences based on the presence of the polybasic HA cleavage site (HACS) with 46.33% and 53.67% of sequences previously identified as highly pathogenic (HP) and low pathogenic (LP), respectively. We compared the performance of different ML classifiers (e.g., logistic regression (LR) with the lasso and ridge regularization, random forest (RF), K-nearest neighbor (KNN), Naïve Bayes (NB), support vector machine (SVM), and convolutional neural network (CNN)) for pathogenicity classification of raw H5Nx nucleotide and protein sequences using a 10-fold cross-validation technique. We found that different ML techniques can be successfully used for the pathogenicity classification of H5 sequences with ∼99% classification accuracy. Our results indicate that for pathogenicity classification of (1) aligned deoxyribonucleic acid (DNA) and protein sequences, with NB classifier had the lowest accuracies of 98.41% (+/-0.89) and 98.31% (+/-1.06), respectively; (2) aligned DNA and protein sequences, with LR (L1/L2), KNN, SVM (radial basis function (RBF)) and CNN classifiers had the highest accuracies of 99.20% (+/-0.54) and 99.20% (+/-0.38), respectively; (3) unaligned DNA and protein sequences, with CNN's achieved accuracies of 98.54% (+/-0.68) and 99.20% (+/-0.50), respectively. ML methods show potential for regular classification of H5Nx virus pathogenicity for poultry species, particularly when sequences containing regular markers were frequently present in the training dataset.

Emergence of highly pathogenic H5N2 and H7N1 influenza A viruses from low pathogenic precursors by serial passage in ovo

Highly pathogenic (HPAI) strains emerge from their low pathogenic (LPAI) precursors and cause severe disease in poultry with enormous economic losses, and zoonotic potential. Understanding the mechanisms involved in HPAI emergence is thus an important goal for risk assessments. In this study ostrich-origin H5N2 and H7N1 LPAI progenitor viruses were serially passaged seventeen times in 14-day old embryonated chicken eggs and Ion Torrent ultra-deep sequencing was used to monitor the incremental changes in the consensus genome sequences. Both virus strains increased in virulence with successive passages, but the H7N1 virus attained a virulent phenotype sooner. Mutations V63M, E228V and D272G in the HA protein, Q357K in the nucleoprotein (NP) and H155P in the neuraminidase protein correlated with the increased pathogenicity of the H5N2 virus; whereas R584H and L589I substitutions in the polymerase B2 protein, A146T and Q220E in HA plus D231N in the matrix 1 protein correlated with increased pathogenicity of the H7N1 virus in embryos. Enzymatic cleavage of HA protein is the critical virulence determinant, and HA cleavage site motifs containing multibasic amino acids were detected at the sub-consensus level. The motifs PQERRR/GLF and PQRERR/GLF were first detected in passages 11 and 15 respectively of the H5N2 virus, and in the H7N1 virus the motifs PELPKGKK/GLF and PELPKRR/GLF were detected as early as passage 7. Most significantly, a 13 nucleotide insert of unknown origin was identified at passage 6 of the H5N2 virus, and at passage 17 a 42 nucleotide insert derived from the influenza NP gene was identified. This is the first report of non-homologous recombination at the HA cleavage site in an H5 subtype virus. This study provides insights into how HPAI viruses emerge from low pathogenic precursors and demonstrated the pathogenic potential of H5N2 and H7N1 strains that have not yet been implicated in HPAI outbreaks.

Annexin A2 (ANXA2) interacts with nonstructural protein 1 and promotes the replication of highly pathogenic H5N1 avian influenza virus

BACKGROUND

Non-structural protein 1 (NS1) is a multifunctional protein and a crucial regulatory factor in the replication and pathogenesis of avian influenza virus (AIV). Studies have shown that NS1 can interact with a variety of host proteins to modulate the viral life cycle. We previously generated a monoclonal antibody against NS1 protein; In the current research study, using this antibody, we immunoprecipitated host proteins that interact with NS1 to better understand the roles played by NS1 in communications between virus and host.

RESULTS

Co-immunoprecipitation experiments identified annexin A2 (ANXA2) as a target molecule interacting with NS1. Results from confocal laser scanning microscopy indicated that NS1 co-localized with ANXA2 in the cell cytoplasm. Overexpression of ANXA2 significantly increased the titer of H5N1 subtype HPAIV, whereas siRNA-mediated knockdown of ANXA2 markedly inhibited the expression of viral proteins and reduced the progeny virus titer.

CONCLUSIONS

Our results indicate that ANXA2 interacts with NS1 and ANXA2 expression increases HPAIV replication.

Biochemical and structural characterization of the interface mediating interaction between the influenza A virus non-structural protein-1 and a monoclonal antibody

We have previously shown that a non-structural protein 1 (NS1)-binding monoclonal antibody, termed as 2H6, can significantly reduce influenza A virus (IAV) replication when expressed intracellularly. In this study, we further showed that 2H6 binds stronger to the NS1 of H5N1 than A/Puerto Rico/8/1934(H1N1) because of an amino acid difference at residue 48. A crystal structure of 2H6 fragment antigen-binding (Fab) has also been solved and docked onto the NS1 structure to reveal the contacts between specific residues at the interface of antibody-antigen complex. In one of the models, the predicted molecular contacts between residues in NS1 and 2H6-Fab correlate well with biochemical results. Taken together, residues N48 and T49 in H5N1 NS1 act cooperatively to maintain a strong interaction with mAb 2H6 by forming hydrogen bonds with residues found in the heavy chain of the antibody. Interestingly, the pandemic H1N1-2009 and the majority of seasonal H3N2 circulating in humans since 1968 has N48 in NS1, suggesting that mAb 2H6 could bind to most of the currently circulating seasonal influenza A virus strains. Consistent with the involvement of residue T49, which is well-conserved, in RNA binding, mAb 2H6 was also found to inhibit the interaction between NS1 and double-stranded RNA.

Application of Species Distribution Modeling for Avian Influenza surveillance in the United States considering the North America Migratory Flyways

Highly Pathogenic Avian Influenza (HPAI) has recently (2014-2015) re-emerged in the United States (US) causing the largest outbreak in US history with 232 outbreaks and an estimated economic impact of $950 million. This study proposes to use suitability maps for Low Pathogenic Avian Influenza (LPAI) to identify areas at high risk for HPAI outbreaks. LPAI suitability maps were based on wild bird demographics, LPAI surveillance, and poultry density in combination with environmental, climatic, and socio-economic risk factors. Species distribution modeling was used to produce high-resolution (cell size: 500m x 500m) maps for Avian Influenza (AI) suitability in each of the four North American migratory flyways (NAMF). Results reveal that AI suitability is heterogeneously distributed throughout the US with higher suitability in specific zones of the Midwest and coastal areas. The resultant suitability maps adequately predicted most of the HPAI outbreak areas during the 2014-2015 epidemic in the US (i.e. 89% of HPAI outbreaks were located in areas identified as highly suitable for LPAI). Results are potentially useful for poultry producers and stakeholders in designing risk-based surveillance, outreach and intervention strategies to better prevent and control future HPAI outbreaks in the US.

Identification of combinatorial host-specific signatures with a potential to affect host adaptation in influenza A H1N1 and H3N2 subtypes

BACKGROUND

The underlying strategies used by influenza A viruses (IAVs) to adapt to new hosts while crossing the species barrier are complex and yet to be understood completely. Several studies have been published identifying singular genomic signatures that indicate such a host switch. The complexity of the problem suggested that in addition to the singular signatures, there might be a combinatorial use of such genomic features, in nature, defining adaptation to hosts.

RESULTS

We used computational rule-based modeling to identify combinatorial sets of interacting amino acid (aa) residues in 12 proteins of IAVs of H1N1 and H3N2 subtypes. We built highly accurate rule-based models for each protein that could differentiate between viral aa sequences coming from avian and human hosts. We found 68 host-specific combinations of aa residues, potentially associated to host adaptation on HA, M1, M2, NP, NS1, NEP, PA, PA-X, PB1 and PB2 proteins of the H1N1 subtype and 24 on M1, M2, NEP, PB1 and PB2 proteins of the H3N2 subtypes. In addition to these combinations, we found 132 novel singular aa signatures distributed among all proteins, including the newly discovered PA-X protein, of both subtypes. We showed that HA, NA, NP, NS1, NEP, PA-X and PA proteins of the H1N1 subtype carry H1N1-specific and HA, NA, PA-X, PA, PB1-F2 and PB1 of the H3N2 subtype carry H3N2-specific signatures. M1, M2, PB1-F2, PB1 and PB2 of H1N1 subtype, in addition to H1N1 signatures, also carry H3N2 signatures. Similarly M1, M2, NP, NS1, NEP and PB2 of H3N2 subtype were shown to carry both H3N2 and H1N1 host-specific signatures (HSSs).

CONCLUSIONS

To sum it up, we computationally constructed simple IF-THEN rule-based models that could distinguish between aa sequences of avian and human IAVs. From the rules we identified HSSs having a potential to affect the adaptation to specific hosts. The identification of combinatorial HSSs suggests that the process of adaptation of IAVs to a new host is more complex than previously suggested. The present study provides a basis for further detailed studies with the aim to elucidate the molecular mechanisms providing the foundation for the adaptation process.

Prevalence of the C-terminal truncations of NS1 in avian influenza A viruses and effect on virulence and replication of a highly pathogenic H7N1 virus in chickens

Highly pathogenic (HP) avian influenza viruses (AIV) evolve from low pathogenic (LP) precursors after circulation in poultry by reassortment and/or single mutations in different gene segments including that encoding NS1. The carboxyl terminal end (CTE) of NS1 exhibits deletions between amino acid 202 and 230 with still unknown impact on virulence of AIV in chickens. In this study, NS1 protein sequences of all AIV subtypes in birds from 1902 to 2015 were analyzed to study the prevalence and distribution of CTE truncation (ΔCTE). Thirteen different ΔCTE forms were observed in NS1 proteins from 11 HA and 8 NA subtypes with high prevalences in H9, H7, H6 and H10 and N9, N2, N6 and N1 subtypes particularly in chickens and minor poultry species. With 88% NS217 lacking amino acids 218-230 was the most common ΔCTE form followed by NS224 (3.6%). NS217 was found in 10 and 8 different HA and NA subtypes, respectively, whereas NS224 was detected exclusively in the Italian HPAIV H7N1 suggesting relevance for virulence. To test this assumption, 3 recombinant HPAIV H7N1 were constructed carrying wild-type HP NS1 (Hp-NS224), NS1 with extended CTE (Hp-NS230) or NS1 from LPAIV H7N1 (Hp-NSLp), and tested in-vitro and in-vivo. Extension of CTE in Hp NS1 significantly decreased virus replication in chicken embryo kidney cells. Truncation in the NS1 decreased the tropism of Hp-NS224 to the endothelium, central nervous system and respiratory tract epithelium without significant difference in virulence in chickens. This study described the variable forms of ΔCTE in NS1 and indicated that CTE is not an essential virulence determinant particularly for the Italian HPAIV H7N1 but may be a host-adaptation marker required for efficient virus replication.