No observed effect of homologous recombination on influenza C virus evolution

Improved Algorithmic Complexity for the 3SEQ Recombination Detection Algorithm

Identifying recombinant sequences in an era of large genomic databases is challenging as it requires an efficient algorithm to identify candidate recombinants and parents, as well as appropriate statistical methods to correct for the large number of comparisons performed. In 2007, a computation was introduced for an exact nonparametric mosaicism statistic that gave high-precision P values for putative recombinants. This exact computation meant that multiple-comparisons corrected P values also had high precision, which is crucial when performing millions or billions of tests in large databases. Here, we introduce an improvement to the algorithmic complexity of this computation from O(mn3) to O(mn2), where m and n are the numbers of recombination-informative sites in the candidate recombinant. This new computation allows for recombination analysis to be performed in alignments with thousands of polymorphic sites. Benchmark runs are presented on viral genome sequence alignments, new features are introduced, and applications outside recombination analysis are discussed.

Bioinformatics studies of Influenza A hemagglutinin sequence data indicate recombination-like events leading to segment exchanges

BACKGROUND

The influenza genome is highly variable due primarily to two mechanisms: antigenic drift and antigenic shift. A third mechanism for genetic change, known as copy choice or template switching, can arise during replication when, if two viral strains infect a cell, a part of a gene from the second viral strain can be copied into the growing progeny of a gene of the first viral strain as replacement leading to a new variety of the virus. This template switching between the same genes of the two strains is known as homologous recombination. While genetic drift and shift are well-understood, the presence or absence of intra-segment homologous recombination in influenza genomes is controversial.

CONTEXT AND PURPOSE OF STUDY

We are interested to study the possibility of subunit-wise homologous recombination. The idea is that where well-defined subunits are separated by consensus sequences, it might be possible for template switching to take place at such junctions. The influenza hemagglutinin gene has basically two subunits, HA1 and HA2, with HA1 being mostly surface exposed and containing the active site for binding to cells, while HA2 secures the hemagglutinin to the viral coat. We undertook a thorough search of the major human infecting influenza hemagglutinin gene sequences, viz., the H1N1, H5N1, H3N2 and H7N9 subtypes, over the period 2010-2014 in Asia to determine if certain sequences could be identified that had HA1 from a previous strain and HA2 from another.

RESULTS

Our search yielded several instances where sequence identities between segments of various strains could be interpreted as indicating possibilities of segment exchange. In some cases, on closer examination they turn out to differ by a few mutations in each segment, due perhaps to the short time span of our database.

CONCLUSIONS AND POTENTIAL IMPLICATIONS

The study reported here, and in combination with our earlier observations on the neuraminidase, shows that subunit-wise recombination-like events in the influenza genes may be occurring more often than have been accounted for and merits further detailed studies.

Detection of influenza C virus but not influenza D virus in Scottish respiratory samples

BACKGROUND

A newly proposed genus of influenza virus (influenza D) is associated with respiratory disease in pigs and cattle. The novel virus is most closely related to human influenza C virus and can infect ferrets but infection has not been reported in humans.

OBJECTIVES

To ascertain if influenza D virus can be detected retrospectively in patient respiratory samples.

STUDY DESIGN

3300 human respiratory samples from Edinburgh, Scotland, covering the period 2006-2008, were screened in pools of 10 by RT-PCR using primers capable of detecting both influenza C and D viruses.

RESULTS

Influenza D was not detected in any sample. Influenza C was present in 6 samples (0.2%), compared with frequencies of 3.3% and 0.9% for influenza A and B viruses from RT-PCR testing of respiratory samples over the same period. Influenza C virus was detected in samples from individuals <2 years or >45 years old, with cases occurring throughout the year. Phylogenetic analysis of nearly complete sequences of all seven segments revealed the presence of multiple, reassortant lineages.

CONCLUSION

We were unable to detect viruses related to influenza D virus in human respiratory samples. Influenza C virus was less prevalent than influenza A and B viruses, was associated with mild disease in the young (<2 years) and old (>45 years) and comprised multiple, reassortant lineages. Inclusion of influenza C virus as part of a diagnostic testing panel for respiratory infections would be of limited additional value.

Published sequences do not support transfer of oseltamivir resistance mutations from avian to human influenza A virus strains

Background

Tamiflu (oseltamivir phosphate ester, OE) is a widely used antiviral active against influenza A virus. Its active metabolite, oseltamivir carboxylate (OC), is chemically stable and secreted into wastewater treatment plants. OC contamination of natural habitats of waterfowl might induce OC resistance in influenza viruses persistently infecting waterfowl, and lead to transfer of OC-resistance from avian to human influenza. The aim of this study was to evaluate whether such has occurred.

Methods

A genomics approach including phylogenetic analysis and probability calculations for homologous recombination was applied on altogether 19,755 neuraminidase (N1 and N2) genes from virus sampled in humans and birds, with and without resistance mutations.

Results

No evidence for transfer of OE resistance mutations from avian to human N genes was obtained, and events suggesting recombination between human and avian influenza virus variants could not be traced in the sequence material studied.

Conclusions

The results indicate that resistance in influenza viruses infecting humans is due to the selection pressure posed by the global OE administration in humans rather than transfer from avian influenza A virus strains carrying mutations induced by environmental exposure to OC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-015-0860-9) contains supplementary material, which is available to authorized users.

Reassortment between influenza B lineages and the emergence of a coadapted PB1-PB2-HA gene complex

Influenza B viruses make a considerable contribution to morbidity attributed to seasonal influenza. Currently circulating influenza B isolates are known to belong to two antigenically distinct lineages referred to as B/Victoria and B/Yamagata. Frequent exchange of genomic segments of these two lineages has been noted in the past, but the observed patterns of reassortment have not been formalized in detail. We investigate interlineage reassortments by comparing phylogenetic trees across genomic segments. Our analyses indicate that of the eight segments of influenza B viruses only segments coding for polymerase basic 1 and 2 (PB1 and PB2) and hemagglutinin (HA) proteins have maintained separate Victoria and Yamagata lineages and that currently circulating strains possess PB1, PB2, and HA segments derived entirely from one or the other lineage; other segments have repeatedly reassorted between lineages thereby reducing genetic diversity. We argue that this difference between segments is due to selection against reassortant viruses with mixed-lineage PB1, PB2, and HA segments. Given sufficient time and continued recruitment to the reassortment-isolated PB1-PB2-HA gene complex, we expect influenza B viruses to eventually undergo sympatric speciation.

Evidence of selection pressures of neuraminidase gene (NA) of influenza A virus subtype H5N1 on different hosts in Guangxi Province of China

In the present study, the possible evidence of positive selection was analyzed for the neuraminidase (NA) sequences of Guangxi H5N1 strains of China. Based on an overall site-specific positive selection analysis, it was found that NA gene of H5N1 Guangxi strains underwent purifying selection and no significant positively selected sites were identified. For the branch-specific positive selection analysis, there was no positive selection evidence for the branches leading to different poultry hosts (chicken, duck and goose). Conclusively, positive selection seems not possible (if not rare) for the NA gene in influenza H5N1 subtype, at least for the samples found in Guangxi Province of China.

Tempo and mode in the molecular evolution of influenza C

Abstract:

Influenza C contributes to economic damage caused by working days lost through absence or inefficiency and may occasionally cause an acute respiratory illness in a paediatric setting. All Influenza C sequences from the NCBI Influenza Virus Resource were examined to determine the date of the most recent common ancestor (t-MRCA), the average nucleotide substitution rate, and the location of residues under positive selection, for each of the seven genome segments of this virus. The segment with the deepest phylogeny was found to be segment 4, encoding the haemagglutinin-esterase protein (HE) with mean t-MRCA at 1890 of the common era (AD), at a 95% highest posterior density (HPD) of 1857-1924 AD. Other genome segments have slightly more recent common ancestors, ranging from mean t-MRCAs of 1916 AD (HPD 1891-1937) for segment 7, encoding the two non-structural proteins (NS) to 1944 AD (HPD 1940-1948) for segment 2 encoding the type 1 basic polymerase (PB1). On the basis of the Bayesian analysis a reclassification of lineages within genome segments is proposed. Some evidence for positive selection was found in the receptor-binding domain of the haemagglutinin-esterase protein. However, average ω (omega) values ranged from 0.05 for polymerase basic protein 2 (PB2) to 0.38 for non-structural protein 2 (NS2), suggesting that strong to moderate purifying selection is the main trend. Characteristic combinations of segment lineages were identified (genome constellations) and shown to have a relatively short life-span before being broken up by reassortment.