The Evolutionary Dynamics of Human Influenza B Virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.     

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.     

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.     

Reconstruction of the Evolutionary Dynamics of A(H3N2) Influenza Viruses Circulating in Italy from 2004 to 2012

Background

Influenza A viruses are characterised by their rapid evolution, and the appearance of point mutations in the viral hemagglutinin (HA) domain causes seasonal epidemics. The A(H3N2) virus has higher mutation rate than the A(H1N1) virus. The aim of this study was to reconstruct the evolutionary dynamics of the A(H3N2) viruses circulating in Italy between 2004 and 2012 in the light of the forces driving viral evolution.

Methods

Phylodinamic analyses were made using a Bayesian method, and codon-specific positive selection acting on the HA coding sequence was evaluated.

Results

Global and local phylogenetic analyses showed that the Italian strains collected between 2004 and 2012 grouped into five significant Italian clades that included viral sequences circulating in different epidemic seasons. The time of the most recent common ancestor (tMRCA) of the tree root was between May and December 2003. The tMRCA estimates of the major clades suggest that the origin of a new viral strain precedes the effective circulation of the strain in the Italian population by 6–31 months, thus supporting a central role of global migration in seeding the epidemics in Italy. The study of selection pressure showed that four codons were under positive selection, three of which were located in antigenic sites. Analysis of population dynamics showed the alternation of periods of exponential growth followed by a decrease in the effective number of infections corresponding to epidemic and inter-epidemic seasons.

Conclusions

Our analyses suggest that a complex interaction between the immune status of the population, migrations, and a few selective sweeps drive the influenza A(H3N2) virus evolution. Our findings suggest the possibility of the year-round survival of local strains even in temperate zones, a hypothesis that warrants further investigation.     

Host-Specific and Segment-Specific Evolutionary Dynamics of Avian and Human Influenza A Viruses: A Systematic Review

Understanding the evolutionary dynamics of influenza viruses is essential to control both avian and human influenza. Here, we analyze host-specific and segment-specific Tajima’s D trends of influenza A virus through a systematic review using viral sequences registered in the National Center for Biotechnology Information. To avoid bias from viral population subdivision, viral sequences were stratified according to their sampling locations and sampling years. As a result, we obtained a total of 580 datasets each of which consists of nucleotide sequences of influenza A viruses isolated from a single population of hosts at a single sampling site within a single year. By analyzing nucleotide sequences in the datasets, we found that Tajima’s D values of viral sequences were different depending on hosts and gene segments. Tajima’s D values of viruses isolated from chicken and human samples showed negative, suggesting purifying selection or a rapid population growth of the viruses. The negative Tajima’s D values in rapidly growing viral population were also observed in computer simulations. Tajima’s D values of PB2, PB1, PA, NP, and M genes of the viruses circulating in wild mallards were close to zero, suggesting that these genes have undergone neutral selection in constant-sized population. On the other hand, Tajima’s D values of HA and NA genes of these viruses were positive, indicating HA and NA have undergone balancing selection in wild mallards. Taken together, these results indicated the existence of unknown factors that maintain viral subtypes in wild mallards.     

Evolutionary Dynamics of Multiple Sublineages of H5N1 Influenza Viruses in Nigeria from 2006 to 2008

Highly pathogenic A/H5N1 avian influenza (HPAI H5N1) viruses have seriously affected the Nigerian poultry industry since early 2006. Previous studies have identified multiple introductions of the virus into Nigeria and several reassortment events between cocirculating lineages. To determine the spatial, evolutionary, and population dynamics of the multiple H5N1 lineages cocirculating in Nigeria, we conducted a phylogenetic analysis of whole-genome sequences from 106 HPAI H5N1 viruses isolated between 2006 and 2008 and representing all 25 Nigerian states and the Federal Capital Territory (FCT) reporting outbreaks. We identified a major new subclade in Nigeria that is phylogenetically distinguishable from all previously identified sublineages, as well as two novel reassortment events. A detailed analysis of viral phylogeography identified two major source populations for the HPAI H5N1 virus in Nigeria, one in a major commercial poultry area (southwest region) and one in northern Nigeria, where contact between wild birds and backyard poultry is frequent. These findings suggested that migratory birds from Eastern Europe or Russia may serve an important role in the introduction of HPAI H5N1 viruses into Nigeria, although virus spread through the movement of poultry and poultry products cannot be excluded. Our study provides new insight into the genesis and evolution of H5N1 influenza viruses in Nigeria and has important implications for targeting surveillance efforts to rapidly identify the spread of the virus into and within Nigeria.Since its emergence in 1996 in Guangdong, China, highly pathogenic avian influenza virus of the H5N1 subtype (HPAI H5N1 virus) has disseminated widely across Asia, Europe, and Africa, infecting a range of domestic and wild avian species and sporadically spilling over into humans and other mammals (4, 35). Over time, the HPAI H5N1 virus has diversified into multiple phylogenetically distinct lineages, classified as clades 0 to 9 according to the unified nomenclature system (33). The H5N1 lineage currently circulating in central Asia, the Middle East, Europe, and Africa is referred to as clade 2.2 (33) and has also been described as “EMA” or Qinghai-like in previous publications (4, 17, 27). This clade originated in April 2005 during a large outbreak of a phylogenetically distinct H5N1 virus among wild bird populations at Qinghai Lake in western China (4, 17) and rapidly spread west through central Asia and Europe, eventually reaching Africa in 2006 (27). Clade 2.2 has further diversified, forming the genetic third-order clade 2.2.1 (32) and three genetically distinct sublineages (I, II, and III) (2, 19, 28), all of which are found in Africa.Since 2006 HPAI H5N1 viruses belonging to clade 2.2 have disseminated across multiple countries in western, eastern, and northern Africa: Egypt, Niger, Cameroon, Sudan, Burkina Faso, Djibouti, Ivory Coast, Ghana, Togo, Benin, and Nigeria (2). With a large poultry industry, estimated at 140 million birds (11), Nigeria has experienced several major outbreaks of HPAI H5N1 virus, posing a serious threat to food security and public health in Africa. The first case of HPAI H5N1 virus in Nigeria (sublineage I) occurred in January 2006 in the state of Kaduna, and the virus subsequently was detected in Ghana, Burkina Faso, Ivory Coast, and Sudan (2). In February 2006 sublineage II was reported in Nigeria, and it disseminated widely across the country during 2006 and 2007, also appearing in Togo (2). Clade 2.2.1, which has been prevalent in Egypt, Israel, and the Gaza Strip from 2006 to 2008, was also detected in Nigeria in 2006 (10).By the end of 2007, outbreaks of HPAI H5N1 virus in Nigeria appeared to have been successfully controlled by measures such as “stamping out with compensation,” restrictions on movement of poultry, and enhanced surveillance (13). However, in July 2008 new cases of HPAI H5N1 from a sublineage never previously detected in Africa (sublineage III) were registered in the Nigerian states of Kano and Katsina and in live bird markets in Gombe and Kebbi states (13, 21). Hence, Nigeria is the only African country where viruses belonging to clade 2.2.1 and to three different sublineages (I, II, and III) of clade 2.2 have all been detected. At least three different reassortment events between sublineages have been documented in Nigeria. Salzberg et al. identified the first reassortant strain (which we refer to as “R1”), in which four genome segments (hemagglutinin [HA], NP, NS, and PB1) belong to sublineage I and the other four segments (NA, MP, PA, and PB2) are derived from sublineage II (27). Subsequently, phylogenetic analysis showed that a 2007 reassortant strain (which we refer to as “R3”) contained the HA and NS segments from sublineage I and the other six segments from sublineage II (19, 22). Another reassortant virus (which we refer to as “R5”) contained only the NS gene segment from sublineage I, while the other seven segments were derived from sublineage II (22).Although the genetic diversity of the Nigerian HPAI H5N1 virus population has been well characterized, including multiple introductions of the virus into Nigeria and several reassortment events, little is known about the evolutionary and population growth dynamics of the virus within Nigeria. Particularly understudied are the spatial movements of individual sublineages among Nigeria''s vast poultry population. To explore the spatial, evolutionary, and population dynamics of the multiple H5N1 lineages cocirculating in Nigeria, we conducted a phylogenetic analysis of whole-genome sequences from 106 HPAI H5N1 viruses isolated between 2006 and 2008 and representing all 25 Nigerian states and the Federal Capital Territory (FCT) reporting outbreaks. Using the exact date and location of collection for each viral isolate, we inferred from their phylogenetic relationships the directionality of viral gene flow among Nigerian states and identified critical regions that are likely to serve as key sources for the H5N1 virus in Nigeria.     

The Evolutionary and Epidemiological Dynamics of the Paramyxoviridae

Paramyxoviruses are responsible for considerable disease burden in human and wildlife populations: measles and mumps continue to affect the health of children worldwide, while canine distemper virus causes serious morbidity and mortality in a wide range of mammalian species. Although these viruses have been studied extensively at both the epidemiological and the phylogenetic scales, little has been done to integrate these two types of data. Using a Bayesian coalescent approach, we infer the evolutionary and epidemiological dynamics of measles, mumps and canine distemper viruses. Our analysis yielded data on viral substitution rates, the time to common ancestry, and elements of their demographic history. Estimates of rates of evolutionary change were similar to those observed in other RNA viruses, ranging from 6.585 to 11.350 × 10⁻⁴ nucleotide substitutions per site, per year. Strikingly, the mean Time to the Most Recent Common Ancestor (TMRCA) was both similar and very recent among the viruses studied, ranging from only 58 to 91 years (1908 to 1943). Worldwide, the paramyxoviruses studied here have maintained a relatively constant level of genetic diversity. However, detailed heterchronous samples illustrate more complex dynamics in some epidemic populations, and the relatively low levels of genetic diversity (population size) in all three viruses is likely to reflect the population bottlenecks that follow recurrent outbreaks.     

Molecular Evolutionary and Epidemiological Dynamics of Genotypes 1G and 2B of Rubella Virus

Rubella Virus (RV), which causes measles-like rashes in children, puts millions of infants at risk of congenital defects across the globe. Employing phylogenetic approaches to the whole genome sequence data and E1 glycoprotein sequence data, the present study reports the substitution rates and dates of emergence of all thirteen previously described rubella genotypes, and gains important insights into the epidemiological dynamics of two geographically widely distributed genotypes 1G and 2B. The overall nucleotide substitution rate of this non-vector-borne RV is in the order of 10⁻³ substitutions/site/year, which is considerably higher than the substitution rates previously reported for the vector-borne alphaviruses within the same family. Currently circulating strains of RV share a common ancestor that existed within the last 150 years, with 95% Highest Posterior Density values ranging from 1868 to 1926 AD. Viral strains within the respective genotypes began diverging between the year 1930 s and 1980 s. Both genotype 1G and 2B have shown a decline in effective number of infections since 1990 s, a period during which mass immunization programs against RV were adapted across the globe. Although both genotypes showed some extent of spatial genetic structuring, the analyses also depicted an inter-continental viral dispersal. Such a viral dispersal pattern could be related to the migration of infected individuals across the regions coupled with a low coverage of MMR vaccination.     

Isolation of Influenza A and B Viruses in HeLa Cells

The HeLa cell line which is one of the most popular cell lines was shown to be suitable for isolation of types A (H₃N₂) and B influenza viruses from throat washings of patients. Sixty-nine and 67 out of 147 throat washings taken from patients during the period from January to April, 1994, were positive for influenza A virus in HeLa cells and MDCK cells, respectively. Seven out of 10 throat washings taken between January and March, 1993, were positive for influenza B virus in MDCK. Of these 7, 4 were also positive for HeLa cells.     

A、B、C三型流感病毒病毒学、流行病学、临床特征和流感疫苗

２０世纪人类遭受了４次流感大流行,数千万人失去了生命,流感病毒分Ａ、Ｂ、Ｃ三型,对其病毒学、流行病学和临床特征,以及流感病毒传统疫苗－－灭活疫苗和新型疫苗－－核酸疫苗的研究进展作了论述。     

Epidemiological and Virological Characteristics of Influenza Viruses Circulating in Cambodia from 2009 to 2011

Background

The Cambodian National Influenza Center (NIC) monitored and characterized circulating influenza strains from 2009 to 2011.

Methodology/Principal Findings

Sentinel and study sites collected nasopharyngeal specimens for diagnostic detection, virus isolation, antigenic characterization, sequencing and antiviral susceptibility analysis from patients who fulfilled case definitions for influenza-like illness, acute lower respiratory infections and event-based surveillance. Each year in Cambodia, influenza viruses were detected mainly from June to November, during the rainy season. Antigenic analysis show that A/H1N1pdm09 isolates belonged to the A/California/7/2009-like group. Circulating A/H3N2 strains were A/Brisbane/10/2007-like in 2009 before drifting to A/Perth/16/2009-like in 2010 and 2011. The Cambodian influenza B isolates from 2009 to 2011 all belonged to the B/Victoria lineage represented by the vaccine strains B/Brisbane/60/2008 and B/Malaysia/2506/2004. Sequences of the M2 gene obtained from representative 2009–2011 A/H3N2 and A/H1N1pdm09 strains all contained the S31N mutation associated with adamantanes resistance except for one A/H1N1pdm09 strain isolated in 2011 that lacked this mutation. No reduction in the susceptibility to neuraminidase inhibitors was observed among the influenza viruses circulating from 2009 to 2011. Phylogenetic analysis revealed that A/H3N2 strains clustered each year to a distinct group while most A/H1N1pdm09 isolates belonged to the S203T clade.

Conclusions/Significance

In Cambodia, from 2009 to 2011, influenza activity occurred throughout the year with peak seasonality during the rainy season from June to November. Seasonal influenza epidemics were due to multiple genetically distinct viruses, even though all of the isolates were antigenically similar to the reference vaccine strains. The drug susceptibility profile of Cambodian influenza strains revealed that neuraminidase inhibitors would be the drug of choice for influenza treatment and chemoprophylaxis in Cambodia, as adamantanes are no longer expected to be effective.     

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.     

Independence of Evolutionary and Mutational Rates after Transmission of Avian Influenza Viruses to Swine

In 1979, an H1N1 avian influenza virus crossed the species barrier, establishing a new lineage in European swine. Because there is no direct or serologic evidence of previous H1N1 strains in these pigs, these isolates provide a model for studying early evolution of influenza viruses. The evolutionary rates of both the coding and noncoding changes of the H1N1 swine strains are higher than those of human and classic swine influenza A viruses. In addition, early H1N1 swine isolates show a marked plaque heterogeneity that consistently appears after a few passages. The presence of a mutator mutation was postulated (C. Scholtissek, S. Ludwig, and W. M. Fitch, Arch. Virol. 131:237–250, 1993) to account for these observations and the successful establishment of an avian H1N1 strain in swine. To address this question, we calculated the mutation rates of A/Mallard/New York/6750/78 (H2N2) and A/Swine/Germany/2/81 (H1N1) by using the frequency of amantadine-resistant mutants. To account for the inherent variability of estimated mutation rates, we used a probabilistic model for the statistical analysis. The resulting estimated mutation rates of the two strains were not significantly different. Therefore, an increased mutation rate due to the presence of a mutator mutation is unlikely to have led to the successful introduction of avian H1N1 viruses in European swine.Influenza viruses undergo rapid variation in nature, thereby limiting prevention of epidemics and pandemics (20). Therefore, questions regarding their evolution remain important, and the answers may yield information useful for predicting further antigenic changes and for explaining the occurrence of new pandemic strains. In 1979, an H1N1 influenza virus of avian origin was transmitted to pigs in northern Europe, thereby introducing a new stable lineage (22, 26). Because there is no direct or serologic evidence of previous H1N1 strains in these pigs (35), we do not need to extrapolate the point of introduction from evolutionary changes over time. Therefore, the European swine influenza H1N1 viruses are a model for studying the early evolution of influenza virus strains in a new host after their introduction in the absence of reassortment.Several characteristics of these H1N1 swine viruses make them particularly interesting for studying the evolution of influenza viruses. Analysis of phylogenic data on their HA, NP, M, and NS genes revealed that the evolutionary rates of both the coding and noncoding changes of these strains are up to 54% higher than those of human and classic swine viruses (18, 33). Such higher evolutionary rates might be due to positive Darwinian selection in avian-like swine influenza viruses, reflecting adaptation to a new host; to differences in sampling frequencies (10) of the viral lineages compared; or to a higher mutation rate. The elevated rate of noncoding changes might be due to fixation of linked coding changes. In addition, early H1N1 swine viruses show a marked plaque heterogeneity and an unusually high escape rate in the presence of various monoclonal antibodies (18). These observations suggest that at least one other factor, independent from positive selection of advantageous variants, must have influenced the evolution of the swine H1N1 isolates since 1979. One candidate is the presence of a mutator mutation (24). An initial variant could have a mutation in its polymerase complex that caused it to be more error prone and to generate a broader spectrum of variants. Such a polymerase might be detrimental for an established strain but advantageous under stress conditions such as adaptation to a new environment. After that, variants with a less error-prone polymerase might again become predominant in the population (18).We wanted to determine whether increased mutation rates due to the presence of a mutator mutation might have contributed to the establishment of the European swine H1N1 viruses. To this end, we calculated the mutation rates of A/Swine/Germany/2/81, a well-characterized early H1N1 swine isolate, and A/Mallard/New York/6750/78 (H2N2), a well-established avian isolate; the evolution of avian influenza viruses is considerably slower than that of classic swine and human isolates (16). An avian H1N1 precursor virus was not available for study. To accommodate the inherent high variability of estimated mutation rates, we developed a new probabilistic model (see Appendix) to calculate the mutation rate and its standard deviation (SD) and screened several parallel clones. To facilitate rapid screening of parallel cultures, we evaluated the frequency with which amantadine-resistant mutants developed. Our results suggest that the mutation rates of A/Swine/Germany/2/81 and A/Mallard/New York/6750/78 are not significantly different and fail to reflect their very different evolutionary backgrounds.     

Limited Compatibility of Polymerase Subunit Interactions in Influenza A and B Viruses

Despite their close phylogenetic relationship, natural intertypic reassortants between influenza A (FluA) and B (FluB) viruses have not been described. Inefficient polymerase assembly of the three polymerase subunits may contribute to this incompatibility, especially because the known protein-protein interaction domains, including the PA-binding domain of PB1, are highly conserved for each virus type. Here we show that substitution of the FluA PA-binding domain (PB1-A_1–25) with that of FluB (PB1-B_1–25) is accompanied by reduced polymerase activity and viral growth of FluA. Consistent with these findings, surface plasmon resonance spectroscopy measurements revealed that PA of FluA exhibits impaired affinity to biotinylated PB1-B_1–25 peptides. PA of FluB showed no detectable affinity to biotinylated PB1-A_1–25 peptides. Consequently, FluB PB1 harboring the PA-binding domain of FluA (PB1-AB) failed to assemble with PA and PB2 into an active polymerase complex. To regain functionality, we used a single amino acid substitution (T6Y) known to confer binding to PA of both virus types, which restored polymerase complex formation but surprisingly not polymerase activity for FluB. Taken together, our results demonstrate that the conserved virus type-specific PA-binding domains differ in their affinity to PA and thus might contribute to intertypic exclusion of reassortants between FluA and FluB viruses.     

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.     

Microarray Multiplex Assay for the Simultaneous Detection and Discrimination of Influenza A and Influenza B Viruses

In this study, we present a microarray approach for the typing of influenza A and B viruses, and the subtyping of H1 and H3 subtypes. We designed four pairs of specific multiplex RT-PCR primers and eight specific oligonucleotide probes and prepared microarrays to identify the specific subtype of influenza virus. Through amplification and fluorescent marking of the multiplex RT-PCR products on the M gene of influenza A and B viruses and the HA gene of subtypes H1 and H3, the PCR products were hybridized with the microarray, and the results were analyzed using a microarray scanner. The results demonstrate that the chip developed by our research institute can detect influenza A and B viruses specifically and identify the subtypes H1 and H3 at a minimum concentration of 1 × 10² copies/μL of viral RNA. We tested 35 clinical samples and our results were identical to other fluorescent methods. The microarray approach developed in this study provides a reliable method for the monitoring and testing of seasonal influenza.     

Epidemiological Characteristics of Influenza A and B in Macau, 2010–2018

Influenza is one of the major respiratory diseases in humans. Macau is a tourist city with high density of population and special population mobility. The study on the epidemiological characteristics of influenza in Macau should bring great value for preventing influenza in tourist cities like Macau in the world. In this study, we collected a total of 104,874 samples with influenza-like illness (ILI) in Macau from 2010 to 2018. Chi-square test and binary multivariable logistic regression were used to investigate the epidemiological characteristics of influenza A and B in Macau. Among these ILI samples, the overall positive rate is 17.17% for influenza A and 6.97% for influenza B. The epidemics of influenza in three years (i.e., 2012, 2017 and 2018) differ from the remaining years (i.e., normal years). In a normal year, influenza A occurs year-round whereas influenza B is seasonal. Our research shows significant differences in influenza infections between different age groups in normal years. Interestingly, our analysis shows no significant difference between locals and tourists in influenza A and B infection in a normal year, whereas the odds of influenza A in tourists were significantly higher than those in locals in July 2017 and the odds of influenza B in tourists were significantly higher than those in locals in January–February 2012 and January–February 2018. This is possibly attributed by the policy of free vaccination to everyone in Macau. These findings should be valuable for preventing influenza in not only Macau but also the world.Supplementary InformationThe online version contains supplementary material available at 10.1007/s12250-021-00388-6.