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.     

南京市2011年乙型流感血凝素基因分子特征分析

[目的]分析2011年南京市乙型流感病毒的血凝素(HA)分子学特征.[方法]选择7株2011年南京市不同时间段有代表性的乙型流感毒株进行HA基因序列测定,通过生物信息学方法对HA分子学特征进行分析.[结果]7株乙型流感毒株分为两个系,4株为Victoria,3株为Yamagata;与2011年度疫苗株相比,Victoria和Yamagata系毒株分别在抗原位点146、197和116、198发生了氨基酸替换;其中197和198位点分别是Victoria和Yamagata毒株的受体结合位点,由于上述位点的替换使得Victoria系/Yamagata系毒株分别在197/196位增加了一个潜在的糖基化位点.[结论]2011年南京市乙型流感Victoria 系和Yamagata系病毒同时存在,Victoria/Yamagata毒株197/198位点的氨基酸替换,值得做进一步的探讨.     

The Wilhelmine E. Key 1999 Invitational lecture. Predicting the evolution of human influenza A

We studied the evolution of the HA1 domain of the H3 hemagglutinin gene from human influenza virus type A. The phylogeny of these genes showed a single dominant lineage persisting over time. We tested the hypothesis that the progenitors of this single evolutionarily successful lineage were viruses carrying mutations at codons at which prior mutations had helped the virus to avoid human immune surveillance. We found evidence that eighteen hemagglutinin codons appeared to have been under positive selection to change the amino acid they encoded in the past. Retrospective tests show that viral lineages undergoing the greatest number of mutations in the positively selected codons were the progenitors of future H3 lineages in nine of eleven recent influenza seasons. Codons under positive selection were associated with antibody combining sites A or B or the sialic acid receptor binding site. However, not all codons in these sites had predictive value. Monitoring new H3 isolates for additional changes in positively selected codons might help identify the most fit extant viral strains that arise during antigenic drift.     

Estimating the Lineage Dynamics of Human Influenza B Viruses

The prediction of the lineage dynamics of influenza B viruses for the next season is one of the largest obstacles for constructing an appropriate influenza trivalent vaccine. Seasonal fluctuation of transmissibility and epidemiological interference between the two major influenza B lineages make the lineage dynamics complicated. Here we construct a parsimonious model describing the lineage dynamics while taking into account seasonal fluctuation of transmissibility and epidemiological interference. Using this model we estimated the epidemiological and evolutional parameters with the time-series data of the lineage specific isolates in Japan from the 2010–2011 season to the 2014–2015 season. The basic reproduction number is similar between Victoria and Yamagata, with a minimum value during one year as 0.82 (95% highest posterior density (HPD): 0.77–0.87) for the Yamagata and 0.83 (95% HPD: 0.74–0.92) for Victoria, the amplitude of seasonal variation of the basic reproduction number is 0.77 (95% HPD:0.66–0.87) for Yamagata and 1.05 (95% HPD: 0.89–1.02) for Victoria. The duration for which the acquired immunity is effective against infection by the Yamagata lineage is shorter than the acquired immunity for Victoria, 424.1days (95% HPD:317.4–561.5days). The reduction rate of susceptibility due to immune cross-reaction is 0.51 (95% HPD: 0.084–0.92) for the immunity obtained from the infection with Yamagata against the infection with Victoria and 0.62 (95% HPD: 0.42–0.80) for the immunity obtained from the infection with Victoria against the infection with Yamagata. Using estimated parameters, we predicted the dominant lineage in 2015–2016 season. The accuracy of this prediction is 68.8% if the emergence timings of the two lineages are known and 61.4% if the emergence timings are unknown. Estimated seasonal variation of the lineage specific reproduction number can narrow down the range of emergence timing, with an accuracy of 64.6% if the emergence times are assumed to be the time at which the estimated reproduction number exceeds one.     

Epidemiological and Evolutionary Dynamics of Influenza B Viruses in Malaysia, 2012-2014

Epidemiological and evolutionary dynamics of influenza B Victoria and Yamagata lineages remained poorly understood in the tropical Southeast Asia region, despite causing seasonal outbreaks worldwide. From 2012–2014, nasopharyngeal swab samples collected from outpatients experiencing acute upper respiratory tract infection symptoms in Kuala Lumpur, Malaysia, were screened for influenza viruses using a multiplex RT-PCR assay. Among 2,010/3,935 (51.1%) patients infected with at least one respiratory virus, 287 (14.3%) and 183 (9.1%) samples were tested positive for influenza A and B viruses, respectively. Influenza-positive cases correlate significantly with meteorological factors—total amount of rainfall, relative humidity, number of rain days, ground temperature and particulate matter (PM10). Phylogenetic reconstruction of haemagglutinin (HA) gene from 168 influenza B viruses grouped them into Yamagata Clade 3 (65, 38.7%), Yamagata Clade 2 (48, 28.6%) and Victoria Clade 1 (55, 32.7%). With neuraminidase (NA) phylogeny, 30 intra-clade (29 within Yamagata Clade 3, 1 within Victoria Clade 1) and 1 inter-clade (Yamagata Clade 2-HA/Yamagata Clade 3-NA) reassortants were identified. Study of virus temporal dynamics revealed a lineage shift from Victoria to Yamagata (2012–2013), and a clade shift from Yamagata Clade 2 to Clade 3 (2013–2014). Yamagata Clade 3 predominating in 2014 consisted of intra-clade reassortants that were closely related to a recent WHO vaccine candidate strain (B/Phuket/3073/2013), with the reassortment event occurred approximately 2 years ago based on Bayesian molecular clock estimation. Malaysian Victoria Clade 1 viruses carried H274Y substitution in the active site of neuraminidase, which confers resistance to oseltamivir. Statistical analyses on clinical and demographic data showed Yamagata-infected patients were older and more likely to experience headache while Victoria-infected patients were more likely to experience nasal congestion and sore throat. This study describes the evolution of influenza B viruses in Malaysia and highlights the importance of continuous surveillance for better vaccination policy in this region.     

Evolution of Influenza B Virus in Kuala Lumpur,Malaysia, between 1995 and 2008

Influenza B virus causes significant disease but remains understudied in tropical regions. We sequenced 72 influenza B viruses collected in Kuala Lumpur, Malaysia, from 1995 to 2008. The predominant circulating lineage (Victoria or Yamagata) changed every 1 to 3 years, and these shifts were associated with increased incidence of influenza B. We also found poor lineage matches with recommended influenza virus vaccine strains. While most influenza B virus lineages in Malaysia were short-lived, one circulated for 3 to 4 years.     

Molecular Epidemiology and Phylogenetic Analyses of Influenza B Virus in Thailand during 2010 to 2014

Influenza B virus remains a major contributor to the seasonal influenza outbreak and its prevalence has increased worldwide. We investigated the epidemiology and analyzed the full genome sequences of influenza B virus strains in Thailand between 2010 and 2014. Samples from the upper respiratory tract were collected from patients diagnosed with influenza like-illness. All samples were screened for influenza A/B viruses by one-step multiplex real-time RT-PCR. The whole genome of 53 influenza B isolates were amplified, sequenced, and analyzed. From 14,418 respiratory samples collected during 2010 to 2014, a total of 3,050 tested positive for influenza virus. Approximately 3.27% (471/14,418) were influenza B virus samples. Fifty three isolates of influenza B virus were randomly chosen for detailed whole genome analysis. Phylogenetic analysis of the HA gene showed clusters in Victoria clades 1A, 1B, 3, 5 and Yamagata clades 2 and 3. Both B/Victoria and B/Yamagata lineages were found to co-circulate during this time. The NA sequences of all isolates belonged to lineage II and consisted of viruses from both HA Victoria and Yamagata lineages, reflecting possible reassortment of the HA and NA genes. No significant changes were seen in the NA protein. The phylogenetic trees generated through the analysis of the PB1 and PB2 genes closely resembled that of the HA gene, while trees generated from the analysis of the PA, NP, and M genes showed similar topology. The NS gene exhibited the pattern of genetic reassortment distinct from those of the PA, NP or M genes. Thus, antigenic drift and genetic reassortment among the influenza B virus strains were observed in the isolates examined. Our findings indicate that the co-circulation of two distinct lineages of influenza B viruses and the limitation of cross-protection of the current vaccine formulation provide support for quadrivalent influenza vaccine in this region.     

Reassortment and Insertion-Deletion Are Strategies for the Evolution of Influenza B Viruses in Nature

The evolution of influenza B viruses is poorly understood. Reassortment of influenza B viruses in nature as a means of genetic variation has not been considered to be a major contributor to their evolution. However, the current practice of assigning evolutionary relationships by antigenic analysis of the hemagglutinin of influenza B viruses would fail to detect reassortants. In this study, influenza B viruses isolated within the past 10 years from sites in the United States and China were studied by nucleotide sequencing of the hemagglutinin and neuraminidase genes and construction of phylogenetic trees to assess evolutionary relationships. A group of viruses represented by B/Houston/1/92 possess a hemagglutinin derived from a B/Yamagata/16/88-like strain and a neuraminidase derived from a B/Victoria/2/87-like strain. A second reassortment event between the hemagglutinin of a B/Yamagata/16/88-like virus closely related to the B/Beijing/184/93 strain and the neuraminidase of a B/Victoria/2/87-like strain is represented by a single virus, B/Memphis/3/93. The neuraminidase of the reassortant viruses is most closely related to that of B/Victoria/2/87-like viruses currently circulating in Nanchang, China. A pattern of insertions and deletions in the hemagglutinin and the neuraminidase of different strains of influenza B viruses is observed. Reassortment plays a role in the evolution of influenza B viruses and may necessitate a change in the methods used to assess and identify new influenza viruses.     

Synonymous Codon Usage,GC<Subscript>3</Subscript>, and Evolutionary Patterns Across Plastomes of Three Pooid Model Species: Emerging Grass Genome Models for Monocots

We have analyzed factors affecting the codon usage pattern of the chloroplasts genomes of representative species of pooid grass family. Correspondence analysis of relative synonymous codon usages (RSCU) showed that genes on secondary axis were correlated with their GC_3S values (all r > 0.3, p < 0.05), indicating mutational bias as an important selective force that shaped the variation in the codon usage among chloroplast genes. The Nc-plot showed that although a majority of the points with low-Nc values were lying below the expected curve, a few genes lied on the expected curve. Nc plot clearly showed that mutational bias plays a major role in codon biology across the monocot plastomes. The hydrophobicity and aromaticity of encoded proteins of each species were found to be other factors of codon usage variation. In the view of above light, besides natural selection, several other factors also likely to be involved in determining the selective constraints on codon bias in plastomes of pooid grass genomes. In addition, five codons (B. distachyon), seven codons (H. vulgare), and four codons (T. aestivum) were identified as optimal codons of the three grass chloroplasts. To identify genes evolving under positive selection, rates of nonsynonymous substitutions (Ka) and synonymous substitutions (Ks) were computed for all groups of orthologous gene pairs.     

Maximum Likelihood Estimation on Large Phylogenies and Analysis of Adaptive Evolution in Human Influenza Virus A

Algorithmic details to obtain maximum likelihood estimates of parameters on a large phylogeny are discussed. On a large tree, an efficient approach is to optimize branch lengths one at a time while updating parameters in the substitution model simultaneously. Codon substitution models that allow for variable nonsynonymous/synonymous rate ratios (ω=d _N/d _S) among sites are used to analyze a data set of human influenza virus type A hemagglutinin (HA) genes. The data set has 349 sequences. Methods for obtaining approximate estimates of branch lengths for codon models are explored, and the estimates are used to test for positive selection and to identify sites under selection. Compared with results obtained from the exact method estimating all parameters by maximum likelihood, the approximate methods produced reliable results. The analysis identified a number of sites in the viral gene under diversifying Darwinian selection and demonstrated the importance of including many sequences in the data in detecting positive selection at individual sites. Received: 25 April 2000 / Accepted: 24 July 2000     

Cross-lineage influenza B and heterologous influenza A antibody responses in vaccinated mice: immunologic interactions and B/Yamagata dominance

The annually reformulated trivalent inactivated influenza vaccine (TIV) includes both influenza A/subtypes (H3N2 and H1N1) but only one of two influenza B/lineages (Yamagata or Victoria). In a recent series of clinical trials to evaluate prime-boost response across influenza B/lineages, influenza-na?ve infants and toddlers originally primed with two doses of 2008-09 B/Yamagata-containing TIV were assessed after two doses of B/Victoria-containing TIV administered in the subsequent 2009-10 and 2010-11 seasons. In these children, the Victoria-containing vaccines strongly recalled antibody to the initiating B/Yamagata antigen but induced only low B/Victoria antibody responses. To further evaluate this unexpected pattern of cross-lineage vaccine responses, we conducted additional immunogenicity assessment in mice. In the current study, mice were primed with two doses of 2008-09 Yamagata-containing TIV and subsequently boosted with two doses of 2010-11 Victoria-containing TIV (Group-Yam/Vic). With the same vaccines, we also assessed the reverse order of two-dose Victoria followed by two-dose Yamagata immunization (Group-Vic/Yam). The Group-Yam/Vic mice showed strong homologous responses to Yamagata antigen. However, as previously reported in children, subsequent doses of Victoria antigen substantially boosted Yamagata but induced only low antibody response to the immunizing Victoria component. The reverse order of Group-Vic/Yam mice also showed low homologous responses to Victoria but subsequent heterologous immunization with even a single dose of Yamagata antigen induced substantial boost response to both lineages. For influenza A/H3N2, homologous responses were comparably robust for the differing TIV variants and even a single follow-up dose of the heterologous strain, regardless of vaccine sequence, substantially boosted antibody to both strains. For H1N1, two doses of 2008-09 seasonal antigen significantly blunted response to two doses of the 2010-11 pandemic H1N1 antigen. Immunologic interactions between influenza viruses considered antigenically distant and in particular the cross-lineage influenza B and dominant Yamagata boost responses we have observed in both human and animal studies warrant further evaluation.     

Uncovering the mutation-fixation correlation in short lineages

Background  

We recently reported a highly unexpected positive correlation between the fixation probability of nonsynonymous mutations (estimated by ω) and neutral mutation rate (estimated by K _s) in mammalian lineages. However, this positive correlation was observed for lineages with relatively long divergence time such as the human-mouse lineage, and was not found for very short lineages such as the human-chimpanzee lineage. It was previously unclear how to interpret this discrepancy. It may indicate that the positive correlation between ω and K _s in long lineages is a false finding. Alternatively, it may reflect a biologically meaningful difference between various lineages. Finally, the lack of positive correlation in short lineages may be the result of methodological artifacts.     

Minor shift in background substitutional patterns in the Drosophila saltans and willistoni lineages is insufficient to explain GC content of coding sequences

Background  

Several lines of evidence suggest that codon usage in the Drosophila saltans and D. willistoni lineages has shifted towards a less frequent use of GC-ending codons. Introns in these lineages show a parallel shift toward a lower GC content. These patterns have been alternatively ascribed to either a shift in mutational patterns or changes in the definition of preferred and unpreferred codons in these lineages.     

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (= d_N/d_S), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.     

Spatial,Temporal and Genetic Dynamics Characteristics of Influenza B Viruses in China, 1973–2018

Annual influenza B virus epidemics and outbreaks cause severe influenza diseases in humans and pose a threat to public health. China is an important epidemic area of influenza B viruses. However, the spatial, temporal transmission pathways and the demography history of influenza B viruses in China remain unknown. We collected the haemagglutinin gene sequences sampled of influenza B virus in China between 1973 and 2018. A Bayesian Markov chain Monte Carlo phylogeographic discrete approach was used to infer the spatial and temporal phylodynamics of influenza B virus. The Bayesian phylogeographic analysis of influenza B viruses showed that the North subtropical and South subtropical zones are the origins of the Victoria and Yamagata lineage viruses, respectively. Furthermore, the South temperate and North subtropical zones acted as transition nodes in the Victoria lineage virus dispersion network and that the North subtropical and Mid subtropical zones acted as transition nodes in the Yamagata lineage virus dispersion network. Our findings contribute to the knowledge regarding the spatial and temporal patterns of influenza B virus outbreaks in China.     

Global mtDNA genetic structure and hypothesized invasion history of a major pest of citrus,Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid Diaphorina citri Kuwayama is a key pest of citrus as the vector of the bacterium causing the “huanglongbing” disease (HLB). To assess the global mtDNA population genetic structure, and possible dispersal history of the pest, we investigated genetic variation at the COI gene collating newly collected samples with all previously published data. Our dataset consists of 356 colonies from 106 geographic sites worldwide. High haplotype diversity (H‐mean = 0.702 ± 0.017), low nucleotide diversity (π‐mean = 0.003), and significant positive selection (Ka/Ks = 32.92) were observed. Forty‐four haplotypes (Hap) were identified, clustered into two matrilines: Both occur in southeastern and southern Asia, North and South America, and Africa; lineages A and B also occur in eastern and western Asia, respectively. The most abundant haplotypes were Hap4 in lineage A (35.67%), and Hap9 in lineage B (41.29%). The haplotype network identified them as the ancestral haplotypes within their respective lineages. Analysis of molecular variance showed significant genetic structure (F_ST = 0.62, p < .0001) between the lineages, and population genetic analysis suggests geographic structuring. We hypothesize a southern and/or southeastern Asia origin, three dispersal routes, and parallel expansions of two lineages. The hypothesized first route involved the expansion of lineage B from southern Asia into North America via West Asia. The second, the expansion of some lineage A individuals from Southeast Asia into East Asia, and the third involved both lineages from Southeast Asia spreading westward into Africa and subsequently into South America. To test these hypotheses and gain a deeper understanding of the global history of D. citri, more data‐rich approaches will be necessary from the ample toolkit of next‐generation sequencing (NGS). However, this study may serve to guide such sampling and in the development of biological control programs against the global pest D. citri.     

Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance

Diversifying selection on metabolic pathways can reduce intraspecific gene flow and promote population divergence. An opportunity to explore this arises from mitonuclear discordance observed in an Australian bird Eopsaltria australis. Across >1500 km, nuclear differentiation is low and latitudinally structured by isolation by distance, whereas two highly divergent, parapatric mitochondrial lineages (>6.6% in ND2) show a discordant longitudinal geographic pattern and experience different climates. Vicariance, incomplete lineage sorting and sex‐biased dispersal were shown earlier to be unlikely drivers of the mitonuclear discordance; instead, natural selection on a female‐linked trait was the preferred hypothesis. Accordingly, here we tested for signals of positive, divergent selection on mitochondrial genes in E. australis. We used codon models and physicochemical profiles of amino acid replacements to analyse complete mitochondrial genomes of the two mitochondrial lineages in E. australis, its sister species Eopsaltria griseogularis, and outgroups. We found evidence of positive selection on at least five amino acids, encoded by genes of two oxidative phosphorylation pathway complexes NADH dehydrogenase (ND4 and ND4L) and cytochrome bc₁ (cyt‐b) against a background of widespread purifying selection on all mitochondrial genes. Three of these amino acid replacements were fixed in ND4 of the geographically most widespread E. australis lineage. The other two replacements were fixed in ND4L and cyt‐b of the geographically more restricted E. australis lineage. We discuss whether this selection may reflect local environmental adaptation, a by‐product of other selective processes, or genetic incompatibilities, and propose how these hypotheses can be tested in future.     

Positive Selection in the N-Terminal Extramembrane Domain of Lung Surfactant Protein C (SP-C) in Marine Mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (ω) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Richard Kliman     

Molecular data provide strong evidence of natural hybridization between native and introduced lineages of <Emphasis Type="Italic">Phragmites australis</Emphasis> in North America

Hybridization, both within and between taxa, can be an important evolutionary stimulus for bioinvasions. Novel intra-taxon hybridizations may arise either between formerly allopatric introduced lineages, or between native and introduced lineages. The latter can occur following a cryptic invasion of a non-native lineage, such as the nineteenth century introduction to North America of a European lineage of the common reed Phragmites australis. Previous studies found no evidence of natural hybridization between native and introduced lineages of P. australis, but produced some F₁ hybrids under experimental conditions when the seed parent was native and the pollen parent was introduced. In this study we used microsatellite data to compare genotypes of P. australis along a transect of approximately 2,000 km in eastern North America. Although hybridization appears uncommon, simulations and principle component analysis of genetic data provided strong evidence for natural hybridization at two sites adjacent to Lake Erie in which native and introduced lineages were sympatric. The seed parent was the native lineage in some cases, and the introduced lineage in other cases. There is now the potential for P. australis hybrids to become increasingly invasive, and managers should consider as a priority the removal of introduced stands from sites where they co-exist with native stands.