The epidemiological signature of influenza B virus and its B/Victoria and B/Yamagata lineages in the 21st century

Authors: Saverio Caini aff001;  Gabriela Kusznierz aff002;  Verònica Vera Garate aff002;  Sonam Wangchuk aff003;  Binay Thapa aff003;  Francisco José de Paula Júnior aff004;  Walquiria Aparecida Ferreira de Almeida aff004;  Richard Njouom aff005;  Rodrigo A. Fasce aff006;  Patricia Bustos aff006;  Luzhao Feng aff007;  Zhibin Peng aff007;  Jenny Lara Araya aff008;  Alfredo Bruno aff009;  Doménica de Mora aff009;  Mónica Jeannette Barahona de Gámez aff011;  Richard Pebody aff012;  Maria Zambon aff012;  Rocio Higueros aff013;  Rudevelinda Rivera aff014;  Herman Kosasih aff015;  Maria Rita Castrucci aff016;  Antonino Bella aff017;  Hervé A. Kadjo aff018;  Coulibaly Daouda aff019;  Ainash Makusheva aff020;  Olga Bessonova aff021;  Sandra S. Chaves aff022;  Gideon O. Emukule aff023;  Jean-Michel Heraud aff024;  Norosoa H. Razanajatovo aff024;  Amal Barakat aff025;  Fatima El Falaki aff025;  Adam Meijer aff026;  Gé A. Donker aff001;  Q. Sue Huang aff027;  Tim Wood aff027;  Angel Balmaseda aff028;  Rakhee Palekar aff029;  Brechla Moreno Arévalo aff030;  Ana Paula Rodrigues aff031;  Raquel Guiomar aff032;  Vernon Jian Ming Lee aff033;  Li Wei Ang aff033;  Cheryl Cohen aff034;  Florette Treurnicht aff034;  Alla Mironenko aff036;  Olha Holubka aff036;  Joseph Bresee aff037;  Lynnette Brammer aff022;  Mai T. Q. Le aff038;  Phuong V. M. Hoang aff038;  Clotilde El Guerche-Séblain aff039;  John Paget aff001
Authors place of work: Netherlands Institute for Health Services Research (Nivel), Utrecht, The Netherlands aff001;  National Institute of Respiratory Diseases "Emilio Coni", Santa Fe, Argentina aff002;  Royal Centre for Disease Control, Department of Public Health, Ministry of Health, Thimphu, Bhutan aff003;  Ministry of Health, Department of Surveillance of Transmissible Diseases, Brasília/DF, Brazil aff004;  Virology Department, Centre Pasteur of Cameroon, Yaoundé, Cameroon aff005;  Sub-Department of Viral Diseases, Instituto de Salud Pública de Chile, Santiago, Chile aff006;  Division of Infectious Diseases, Chinese Center for Disease Control and Prevention, Beijing, P.R. China aff007;  National Influenza Center, Ministry of Health, San José, Costa Rica aff008;  National Institute of Public Health Research (INSPI), National Reference Centre for Influenza and Other Respiratory Viruses, Guayaquil, Ecuador aff009;  Agricultural University of Ecuador, Guayaquil, Ecuador aff010;  National Influenza Center, Ministry of Health, San Salvador, El Salvador aff011;  Public Health England, London, England, United Kingdom aff012;  National Influenza Center, Ministry of Health, Guatemala City, Guatemala aff013;  National Influenza Center, Ministry of Health, Tegucigalpa, Honduras aff014;  US Naval Medical Research Unit No.2, Jakarta, Indonesia aff015;  National Influenza Center, Department of Infectious Diseases, National Institute of Health, Rome, Italy aff016;  Department of Infectious Diseases, National Institute of Health, Rome, Italy aff017;  Department of Epidemic Virus, Institut Pasteur, Abidjan, Côte d'Ivoire aff018;  Service of Epidemiological Diseases Surveillance, National Institute of Public Hygiene, Abidjan, Côte d'Ivoire aff019;  National Center of Expertise, Committee of Public Health Protection, Ministry of Health, Astana, Kazakhstan aff020;  National Center of Expertise, Committee of Public Health Protection, Ministry of Health, Uralsk City, Kazakhstan aff021;  Influenza Division, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America aff022;  Influenza Program, Centers for Disease Control and Prevention, Nairobi, Kenya aff023;  National Influenza Center, Virology Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar aff024;  National Influenza Center, Institut National d'Hygiène, Ministry of Health, Rabat, Morocco aff025;  National Institute for Public Health and the Environment, Centre for Infectious Diseases Research, Diagnostics and Laboratory Surveillance, Bilthoven, The Netherlands aff026;  Institute of Environmental Science and Research, Weillngton, New Zealand aff027;  National Influenza Center, Ministry of Health, Managua, Nicaragua aff028;  Pan American Health Organization, Washington, District of Columbia, United States of America aff029;  National Influenza Center, IC Gorgas, Panama City, Panama aff030;  Department of epidemiology, National Institute of Health Doutor Ricardo Jorge, Lisbon, Portugal aff031;  National Influenza Reference Laboratory, National Institute of Health Doutor Ricardo Jorge, Lisbon, Portugal aff032;  Public Health Group, Ministry of Health, Singapore, Singapore aff033;  Centre for Respiratory Diseases and Meningitis, National Institute for Communicable Diseases of the National Health Laboratory Service, Johannesburg, South Africa aff034;  School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa aff035;  L.V.Gromashevsky Institute of Epidemiology and Infectious Diseases, National Academy of Medical Science of Ukraine, Department of Respiratory and other Viral Infections, Kyiv, Ukraine aff036;  Influenza Division, National Center for Immunizations and Respiratory Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of America aff037;  National Institute of Hygiene and Epidemiology, Hanoi, Vietnam aff038;  Global Vaccine Epidemiology and Modeling Department (VEM), Franchise Epidemiologist, Sanofi Pasteur, Lyon, France aff039
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0222381


We describe the epidemiological characteristics, pattern of circulation, and geographical distribution of influenza B viruses and its lineages using data from the Global Influenza B Study. We included over 1.8 million influenza cases occurred in thirty-one countries during 2000–2018. We calculated the proportion of cases caused by influenza B and its lineages; determined the timing of influenza A and B epidemics; compared the age distribution of B/Victoria and B/Yamagata cases; and evaluated the frequency of lineage-level mismatch for the trivalent vaccine. The median proportion of influenza cases caused by influenza B virus was 23.4%, with a tendency (borderline statistical significance, p = 0.060) to be higher in tropical vs. temperate countries. Influenza B was the dominant virus type in about one every seven seasons. In temperate countries, influenza B epidemics occurred on average three weeks later than influenza A epidemics; no consistent pattern emerged in the tropics. The two B lineages caused a comparable proportion of influenza B cases globally, however the B/Yamagata was more frequent in temperate countries, and the B/Victoria in the tropics (p = 0.048). B/Yamagata patients were significantly older than B/Victoria patients in almost all countries. A lineage-level vaccine mismatch was observed in over 40% of seasons in temperate countries and in 30% of seasons in the tropics. The type B virus caused a substantial proportion of influenza infections globally in the 21st century, and its two virus lineages differed in terms of age and geographical distribution of patients. These findings will help inform health policy decisions aiming to reduce disease burden associated with seasonal influenza.


Medicine and health sciences – Infectious diseases – Viral diseases – Influenza – Infectious disease control – Vaccines – Pathology and laboratory medicine – Pathogens – Influenza B virus – Biology and life sciences – Organisms – Viruses – RNA viruses – Orthomyxoviruses – Influenza viruses – Influenza A virus – Microbiology – Medical microbiology – Microbial pathogens – Viral pathogens – Population biology – Population metrics – Age distribution – Earth sciences – Geography – Geographic areas – Northern Hemisphere – Southern Hemisphere


Influenza causes a major burden of disease on populations globally [1]. The impact of seasonal and pandemic influenza on population health in high-income countries has been described extensively [2]. In Europe, influenza ranks first among infectious diseases in terms of burden, accounting for 30% of the total disability-adjusted life years (DALYs) lost due to infectious diseases annually [3]. More recently, evidence has been accumulating that the burden of disease associated with influenza is high in low- and middle-income countries as well [4]. Type A viruses cause most influenza cases and are responsible for pandemics, but influenza B is an important cause of morbidity and mortality during interpandemic periods, and its prevention represents an important public health priority globally [56]. Influenza B viruses split into two antigenically distinct phylogenetic lineages (B/Victoria/2/87 representative, abbreviated B/Victoria, and B/Yamagata/16/88 representative, abbreviated B/Yamagata) in the early 1980s [7]. Whilst the circulation of the Victoria lineage was geographically limited to eastern Asia for most of the 1990s, the two lineages have co-circulated globally in the 21st century [89].

The mainstay of influenza prevention is vaccination, aimed primarily at individuals who are at greater risk of developing complications when infected, like the elderly, pregnant women, and people with underlying medical conditions, and more recently at children to provide both direct protection to the children and indirect protection to the wider population [10]. Both the composition and period of administration of the vaccine are crucial to ensure adequate immunity against influenza; thus, investigating the timing of influenza epidemics and the patterns of circulation of the different influenza virus types, subtypes and lineages is of great importance from a public health perspective [11]. Until 2012, trivalent influenza vaccines (TIV) were available that contained only one type B influenza virus, belonging to either the B/Victoria or B/Yamagata lineage. Each year, the WHO issued recommendations for the B virus strain to be included in the TIV, but the frequency of lineage-level vaccine mismatch has been high (around 50%) [5]. This often resulted in unsatisfactory protection, considering the limited efficacy against the mismatched lineage [12], and substantial health impact [13]. Quadrivalent influenza vaccine (QIV), containing one B/Victoria and one B/Yamagata virus strain, were first approved by the US Food and Drug Administration (US-FDA) in 2012. This prompted the establishment of several studies that aimed to describe the epidemiology, clinical characteristics and burden of disease of influenza B in more depth, and to compare the cost-effectiveness of QIV vs. TIV in different settings and populations.

Despite significant advances in our understanding of influenza B [14], important knowledge gaps persist, in particular concerning the pattern of circulation and geographical distribution of B/Victoria and B/Yamagata viruses, and the age distribution of patients who are most often infected with either lineage. These aspects are of great importance, however, for instance in order to evaluate the cost-effectiveness of alternative influenza prevention strategies and thus enable more appropriate public health choices (e.g. around the choice of tri- and quadrivalent vaccines) [15]. In fact, while some reports have addressed this topic in single countries [1622], a global overview is still lacking. Here, we conducted a global analysis of the epidemiology of influenza B virus and its lineages using the database of the Global Influenza B Study (GIBS).


The Global Influenza B Study

The rationale and methodology of the Global Influenza B Study has been described in detail elsewhere [14,2324]. The GIBS database encompasses epidemiological and virological influenza surveillance data from thirty-one countries around the world (Fig 1). Participating countries were classified as located in the Northern hemisphere (n = 11), inter-tropical belt (n = 15), or Southern hemisphere (n = 5) based on the latitude of the country’s population centre [25]. Participating countries were required to make available influenza surveillance data for as many consecutive years as possible (from 2000 onwards). Data include the weekly number of laboratory-confirmed influenza cases broken down by virus type (A, B), type A subtype (H1N1, 2009 pandemic H1N1, H3N2, A not subtyped), and type B lineage (Victoria, Yamagata, B not characterized) and weekly influenza-like illness/acute respiratory infection rates (per 100,000 population or 100 consultations, depending on country). Information on age (exact age or age groups) was also available. In most participating countries, the influenza surveillance system covers the whole national territory, and a sample of influenza virus positive specimens or virus isolates is routinely sent to a WHO collaborating centre for further investigations. However, countries differ in other important aspects of their surveillance system, e.g. the clinical case definition in use, the distribution of outpatients and inpatients being sampled, and the availability of population denominators [14,24]. The GIBS database was initially assembled in 2013–14 by merging influenza surveillance data up to as late as December 2013 from thirty countries [14], but was updated for the purposes of this paper with data until 2015 to 2018 (depending on the country) for twenty-five countries. The Netherlands joined the project in 2018, bringing the total number of participating countries to thirty-one.

<h2>Countries participating in the Global Influenza B Study (GIBS), 2000–2018.</h2>
Fig. 1.

Countries participating in the Global Influenza B Study (GIBS), 2000–2018.

Data for China were provided separately for the Northern and Southern parts of the country, which were therefore entered separately in the analyses (for brevity, however, we will us the term “country” to refer to a whole country or a part of it henceforth). Brazil also provided data stratified by sub-national regions (i.e. five administrative regions: north, north-east, central-west, south, and south-east); however, geographically-stratified data were too sparse in some years, thus we opted to treat Brazil as a single country in the analysis.

Statistical methods

Similar to previous GIBS publications [14,2324], the unit of analysis was the “season”: this corresponded to the calendar year in tropical countries and countries located in the Southern hemisphere, and was defined as the period between the 27th week of a year and the 26th week of the following year for countries located in the Northern hemisphere. As previously explained [23], the purpose of this methodological approach is to give each “season” an equal weighting in the analysis, thus limiting the impact on results arising from any differences in reporting between countries (e.g. high- vs. low-resources countries) and over time within the same country (e.g. before vs. after the 2009 pandemic).

For each country, we included in the analyses the seasons with at least 20 weeks of reported data and at least 100 laboratory-confirmed influenza cases. In each season, we calculated the percentage of influenza cases that were due to influenza B. We then determined the percentage of cases that were caused by viruses belonging to the B/Victoria or B/Yamagata lineage: in this analysis, to increase the stability of the results, we only included (for each country) the subset of seasons with 50 or more characterized influenza B viruses. The non-parametric Kruskal-Wallis test was applied to compare the median proportion of influenza A and B cases (over all influenza cases) and B/Victoria and B/Yamagata cases (over all cases with characterized influenza B viruses) between countries based on their latitude (Northern hemisphere, inter-tropical belt, or Southern hemisphere) and population age structure (median age below 30 years, between 30 and 35 years, or above 35 years) [25].

We compared the timing of the primary peak of influenza A and B epidemics in each country using the EPIPOI software [26]. For this analysis, we excluded data from the 2009 season (2009–2010 in Northern hemisphere countries) because of the markedly atypical timing of the influenza A(H1N1)pdm09 pandemics. Time series were first standardized by dividing the weekly number of influenza cases by the maximum weekly number (per country by defined season). Next, a periodic annual function (PAF) was generated by summing up the annual, semi-annual and quarterly harmonics (obtained by Fourier decomposition). The timing of the primary peak (defined as the month where the PAF reaches its maximum value) was then compared for A and B epidemics in each country.

We aimed to evaluated the frequency of influenza B lineage-level vaccine mismatch, which was defined as a mismatch between the lineage that caused the majority (>50%) of influenza B cases in a given season and country, and the lineage included in the trivalent vaccine. Information on influenza vaccine formulation was retrieved from the WHO website [27]. The above definition is straightforward for temperate countries in the Northern and Southern hemisphere, but much less so for countries of the inter-tropical belt. While the WHO issues recommendations on the composition of influenza vaccines to be used in tropical countries [28], these countries often adopt either formulation based on local considerations or make no specific recommendations [29]. To cope with this and consider all possible scenarios (and also for consistency with previous GIBS publications [14]), the frequency of B lineage-level vaccine mismatch in tropical countries was calculated by assuming that: (i) all countries situated north of the equator use the northern hemisphere vaccine, and all countries south of the equator use the southern hemisphere vaccine formulation; (ii) all countries use the northern hemisphere vaccine formulation; or (iii) all countries use the southern hemisphere vaccine formulation.

The recent update of the GIBS database and the resulting increase in the number of cases with characterized influenza B viruses made it possible to compare the age distribution of influenza B/Victoria and B/Yamagata cases in a similar way to what had been done previously for influenza B vs. A (and its subtypes) [14,24]. This analysis was limited to countries for which information on the exact age was available for at least 50 B/Victoria and at least 50 B/Yamagata influenza cases over all seasons available. The virus-lineage-specific age distribution of influenza B cases was visualized, separately in each country, using histograms with 5-age-year-wide bars. Since the distributions were skewed to the right, we used the non-parametric Wilcoxon rank sum test to compare the median age of influenza B cases by virus lineage.

Statistical analyses were conducted using Stata software version 14 (StataCorp LP, College Station, TX, USA). All statistical tests were two-sided, and p-values were considered statistically significant when below 0.05.


The database of the Global Influenza B Study encompassed a total of 1.820.301 influenza cases between 2000 and 2018, of which 419,167 (23.0%) had type B influenza. Cases were unevenly distributed between countries, with USA and Australia contributing 54.1% and 25.3% of all influenza cases to the database, respectively. Information on age was not available for the USA; for the remaining countries, information on age was available for most (95.0%) influenza cases. After applying our exclusion criteria, 299 seasons with 100 or more reported influenza cases overall were included in the analysis, of which 110 were from Northern hemisphere countries, 131 from countries located in the inter-tropical belt, and 58 from Southern hemisphere countries. The median number of seasons per country was nine, ranging from three (for Turkey) to eighteen (for New Zealand). The number of overall and influenza B cases in each country and season included in the analysis, the proportion of influenza B cases of which the virus was characterized, and their breakdown into B/Victoria and B/Yamagata lineages, are provided in S1 Table.

Frequency of influenza B

Influenza B virus caused a median 23.4% (interquartile range [IQR] 9.3–38.9%) of reported influenza cases in a season. More specifically, the proportion of influenza cases caused by the B virus type was <20% in 118 seasons (42.3%), between 20% and 50% in 125 seasons (41.8%), and above 50% in 45 seasons (15.1%). The proportion of influenza B cases over all influenza cases reported in the season varied geographically (Fig 2): its median value tended to be higher (with borderline statistical significance, p = 0.060) in countries of the inter-tropical belt (27.4%, IQR 12.2–41.7%) compared to temperate countries of the Northern (21.0%, IQR 7.3–37.4%) and Southern (22.2%, IQR 9.1–34.5%) hemispheres. The proportion of seasons in which influenza B cases constituted between 20 and 50% or more than 50% of all reported influenza cases did not significantly differ (p = 0.275) between tropical countries (47.3% and 16.0%, respectively) and temperate countries of the two hemispheres (39.9% and 14.3%, respectively).

<h2>Distribution of influenza seasons by proportion of influenza B cases and geographical area (Inter-tropical belt, Northern hemisphere, Southern hemisphere).</h2>
Fig. 2.

Distribution of influenza seasons by proportion of influenza B cases and geographical area (Inter-tropical belt, Northern hemisphere, Southern hemisphere).

Red bars indicate 25th, 50th (median) and 75th percentiles. The Global Influenza B Study, 2000–2018.

Timing of influenza A and B epidemics

In countries of the Southern hemisphere, influenza A epidemics typically peaked in July-September, on average 1.1 month earlier than influenza B (August-September) (Table 1). Influenza epidemics also peaked in the winter months of Northern hemisphere countries, earlier for influenza A (January-February, except Ukraine, in March) than for influenza B (February-March, except England, in January), with an average difference of 0.6 months. The only exception was Bhutan (i.e. the southernmost countries among those situated in the Northern hemisphere), where both influenza A and B epidemics typically peaked in August. There was more variability in countries of the inter-tropical belt: the timing of the primary peak varied widely between countries and could take place in practically any month of the year, with a small difference in timing (0.3 months earlier on average) for the primary peak of influenza A and B epidemics.

Tab. 1.

Typical timing of the peak of influenza A and B epidemics according to countries’ latitude.

<h2>Typical timing of the peak of influenza A and B epidemics according to countries’ latitude.</h2>

The Global Influenza B Study, 2000–2018.

Frequency of B/Victoria and B/Yamagata lineages viruses

There were 84 seasons where the virus lineage was determined for 50 or more influenza B cases (out of a total 299 seasons, 28.1%) and this data covered 18 countries (out of 31). In these seasons, the B/Victoria and B/Yamagata lineage viruses caused a similar range of proportions of influenza B cases. More precisely, the median proportion of influenza B cases reported in a season that were caused by B/Victoria lineage viruses was 46.0% (IQR 5.3–89.6%). The range of the proportion of B/Victoria over all influenza B cases was 0%-100%, meaning that there were seasons in which all influenza B cases were caused by either virus lineage. Co-circulation of the two lineages was frequent: in 27 (32.1%) of 84 seasons, both lineages accounted for at least 20% of influenza B cases, while the B/Victoria and B/Yamagata lineages caused over 80% of all influenza B cases in 27 (32.1%) and 30 (35.8%) seasons, respectively.

The proportion of either lineage over all influenza B cases varied according to the age structure of countries, and consequently, according to the countries’ latitude as well (Fig 3). The proportion of influenza B cases caused by the B/Victoria lineage in a given season varied widely (ranging from 0% to nearly 100%) also in the analyses stratified by the country’s age structure and latitude. However, the B/Victoria lineage was more frequent (p-value 0.049) in countries with a median age below 30 years (median proportion over all influenza B cases = 67.0%) compared to countries with a median age was 30–35 years (42.0%) or above 35 years (36.2%). Likewise, the median proportion of B/Victoria over all characterized influenza B cases in a season was higher in countries of the inter-tropical belt (57.5%) compared to temperate countries of the Northern (27.9%) and the Southern hemisphere (43.8%) (p-value = 0.048).

Fig. 3.

Proportion of influenza B cases in a season that were caused by the B/Victoria lineage viruses, according to countries median age (top) and geographical area (bottom). The Global Influenza B Study, 2000–2018.

Age distribution of B/Victoria and B/Yamagata influenza cases

In Table 2, we compared the median age of influenza B cases according to virus lineage. In Australia, B/Victoria cases were non-significantly older than B/Yamagata cases: the median age was 11 vs. 8 years, respectively. In all other countries, B/Yamagata cases were older than B/Victoria cases, and the difference between the median age was always significant except in Madagascar (p-value 0.248) and in Ukraine (p-value 0.082).

Tab. 2.

Age distribution of influenza B/Victoria and B/Yamagata patients by country (sorted by median age of the general population).

<h2>Age distribution of influenza B/Victoria and B/Yamagata patients by country (sorted by median age of the general population).</h2>

The Global Influenza B Study, 2000–2018.

The visual inspection of the graphs revealed a fairly consistent pattern in the age distribution of B/Victoria and B/Yamagata cases (S1 File). In most countries, influenza B/Victoria cases tended to be distributed according to their age along a unimodal curve, with a peak below 10 years of age, while the age distribution of B/Yamagata cases frequently followed a bimodal curve, with an earlier, larger peak below 10 years of age, and a smaller, yet discernible peak at older age, mostly between 25 and 50 years of age. The adult age peak of B/Yamagata cases was barely noticeable in countries where the majority of inhabitants are younger than 25 years (e.g. Madagascar and Ivory Coast), and became increasingly evident as the median age of the country population increased. A partial exception was Portugal, where the peak in the age distribution of B/Victoria cases was between 10 and 20 years of age, and B/Yamagata distributed along a unimodal curve peaking at above 50 years of age.

Frequency of lineage-level influenza B trivalent vaccine mismatches

The proportion of influenza B lineage-level vaccine mismatch was 54.2% in countries in the Southern hemisphere, and 42.9% in countries of the Northern hemisphere. For countries in the inter-tropical belt, the proportion of seasons with a lineage-level vaccine mismatch was 29.6% in the scenario in which these countries use the vaccine composition recommended for the hemisphere they are situated in; 40.7% assuming all these countries used the recommended vaccine composition for the Northern hemisphere; and 22.2% if all countries used the recommended vaccine composition for the Southern hemisphere.


Our analysis found that influenza B virus was responsible for nearly one fourth of all influenza cases in an average season between 2000 and 2018. More specifically, influenza B virus accounted for over 20% of all influenza cases in more than half of the seasons, and was the dominant virus type in about one out of every seven seasons. This proportion varied geographically, as the proportion of influenza B over all influenza cases tended to be higher on average in tropical countries compared to countries of the Northern and Southern hemispheres. These findings were broadly in line with reports from countries around the world [5,9,2021,3031] and with GIBS data published in 2015 [14], with moderate discrepancies (less than 5%) being most likely linked to the variability of the time span covered by different studies. Influenza B epidemics tended to peak on average three weeks later than influenza A epidemics during the winter period in temperate countries of both hemispheres, which is consistent with previous country reports [21,30,3234]. The timing of influenza A and B epidemics was different in tropical countries where they were highly heterogeneous and showed no consistent pattern in the timing of the different epidemics. Influenza viruses exhibit differential ability to induce temporary immunity against re-infection with other viruses in experimental settings [35]. This “viral hierarchy” may help explain the differences in timing of epidemics caused by different virus types and subtypes where the weather conditions that are most favourable to the occurrence of influenza epidemics (e.g. cold-dry) last for only a limited number of weeks each year.

The study of patterns of circulation of B/Victoria and B/Yamagata lineages revealed two main findings. Firstly, the proportion of B cases caused by either lineage in a given season ranged between 0% and nearly 100% in each geographical area (i.e. anywhere around the world seasons occurred where all influenza B cases were caused by either the B/Victoria or the B/Yamagata lineage), which emphasizes the challenges faced when trying to predict which lineage will cause most influenza B cases next season. This helps explain the very high frequency of lineage-level vaccine mismatch for the TIV which was around 50% (i.e. close to the scenario in which the choice of the vaccine lineage is made at random) in temperate countries of the two hemispheres, which has been reported in the past [30,3637].

The second important finding that emerges from the B/Victoria-B/Yamagata lineage analysis is that whilst the two lineages caused a comparable proportion of influenza B cases globally during the study period, their distribution varied geographically, as B/Victoria was relatively more frequent in tropical countries, while B/Yamagata was more frequent in temperate climate countries of the Southern and Northern hemispheres. The reason for this pattern of circulation of the two lineages is not clear. One possibility is that the survival and transmissibility of B/Victoria and B/Yamagata viruses is differentially affected by “cold-dry” and “humid-rainy” environmental conditions that drive the occurrence of influenza epidemics in temperate and tropical countries, respectively [38]. This hypothesis does not seem, however, to be supported by the available evidence, which suggests instead that influenza viruses of different types and subtypes (e.g. A(H3N2), A(H1N1)pdm09, and B) are similarly affected by weather conditions [39]. Based on our findings on the differential age distribution of B/Victoria and B/Yamagata infected cases, a more likely explanation for the unequal geographical distribution of the two B lineages lies is in the diverse demographic structure of world countries, with those located around the tropics having a lower median age, on average, than those in temperate climates.

In this regard, our study showed consistently at a global level (i.e. across countries that differ greatly under many aspects, including the type of influenza surveillance system) that B/Yamagata cases are on average older than B/Victoria cases, and in particular, that the former virus lineage infects adult individuals (≥25 years) more frequently than the latter. The only exception (i.e. the only country in which B/Victoria cases were older than B/Yamagata cases, albeit not significantly) was Australia. The reason for this divergent pattern is not clear; however, since the majority of B/Yamagata cases in the database from Australia occurred in 2008 (see S1 Table), this finding might be driven by a single season with an unusual age distribution of B/Yamagata patients. Importantly, our findings support earlier work showing that each influenza A subtype (H3N2, pre-2009-pandemic H1N1, and H1N1pdm09) tends to affect a different age group [24]. A recent study has shown that the B/Victoria viruses antigenic drift parallels that of A(H3N2), with limited cross-reactivity between phylogenetic clusters compared to B/Yamagata (which, in contrast, resembled A(H1N1) under this regard, with multiple variants circulating and greater levels of antigenic cross-reactivity) [40]. According to these observations, one would expect re-infection among adult individuals to occur more frequently for B/Victoria rather than B/Yamagata, but the observed pattern runs counter to this expectation. Vijaykrishna and colleagues suggested that the observed differential age distribution may be explained by either the higher effective reproductive number of B/Victoria; a broader response against B/Victoria viruses in older people; and/or differences in the prevalence of receptor binding molecules that help B/Victoria and B/Yamagata viruses infect respiratory tract cells in young children and adults [40]. Another hypothesis is that the bimodal age distribution of B/Yamagata cases may be the result of a shorter duration of acquired immunity developed after infection with B/Yamagata vs. B/Victoria viruses [41]. Finally, antigenic original sin (i.e. preferential response to influenza B viruses encountered in the past [42]) may also play a role. If this were true, the mean age of patients infected with either B lineage would be expected to drift over time: this has occurred in Australia between 2008–2011 [40], but we were unable to replicate these results in the GIBS database (results not shown).

These findings have important implications for assessing the impact of influenza B viruses in a population each season. Knowing early in the season (through surveillance activities) which influenza B lineage virus is circulating will allow one to predict which age groups will be more affected, with B/Yamagata seasons being generally characterized by more adults and elderly people being infected and, therefore, by a greater burden of disease compared to B/Victoria seasons. This appears to have been the case during the 2017/18 season in Europe, when B/Yamagata was dominant (and there was a lineage-level vaccine mismatch) and a high burden of influenza B was observed in many European countries [43].

We analyzed influenza surveillance data from around the world to determine the epidemiology of influenza B virus and its lineages in the 21st century. We have enriched the analysis of an earlier GIBS paper [14] by including more countries (31 vs. 26) and additional years (data to up to 2018): the overall number of reports analysed is now twice as large (>1.8 million vs. 935,000 cases). This allowed us to conduct a comprehensive comparative analysis of the epidemiological characteristics of B/Victoria and B/Yamagata virus lineages, thus achieving a better understanding of the epidemiology of influenza B globally. The most important limitation of our study is that it relies on data that were collected within national surveillance systems rather than for research purposes. Countries participating in GIBS inevitably differ in terms of a number of important characteristics, including percentage of out- and in-patients included, number of years with data and average number of influenza cases reported per year, proportion of influenza B viruses that were characterized availability of information on age of these cases and sampling strategies across age groups. In each individual country, changes in the aforementioned characteristics may have occurred over time, for instance before and after the 2009 pandemic. Another limitation of our study is the uneven distribution of cases among countries and seasons, which is partly due to differences among countries in terms of population, the intensity of surveillance activities, and the coverage of each country’s influenza surveillance system. Throughout this paper, we have used statistical methods that were intended to minimize the impact of this variability on the results. These include the use of the ‘season’ as a unit of analysis for the comparison of the frequency of B influenza across countries, the standardization offered by the EPIPOI software in the analysis of temporal patterns of influenza epidemics, and the use of medians to compare the age distribution of influenza cases between countries. However, we acknowledge that these sources of variability may have affected our findings to an extent that is difficult to quantify precisely, and a further improvement in the harmonization of influenza surveillance data remains desirable as it would allow researchers to conduct more detailed and in-depth analyses and obtain more reliable results.

In conclusion, our study provides important new information that describes the epidemiology of influenza B viruses and its lineages during the start of the 21st century. We showed that influenza B virus caused nearly one fourth (on average) of all influenza cases occurring annually in the world during 2000–2018, and that both lineages accounted each for a comparable proportion of all influenza B cases globally. We also observed that cases infected with the B/Victoria or B/Yamagata lineages differed consistently (i.e. in a similar fashion across countries) in terms of their age distribution, and suggested that this may help explain the uneven geographical distribution of the two lineages globally. Finally, we confirmed that the TIV suffers from a high frequency of lineage-level mismatch, especially in temperate countries of the two hemispheres. While expanding our knowledge of the epidemiology and pattern of circulation of influenza B virus, these findings will help inform health policy decisions aiming to reduce disease burden associated with seasonal influenza locally, regionally and globally.

Supporting information

S1 Table [a]
Laboratory-confirmed influenza cases reported in each country and season, and proportion of those caused by the influenza B virus type (overall and by lineage).

S1 File [pdf]
Age distribution of laboratory-confirmed influenza patients infected with type B viruses belonging to the Victoria or Yamagata lineage.

S2 File [doc]
Contact information of data owners.


1. Iuliano AD, Roguski KM, Chang HH, Muscatello DJ, Palekar R, Tempia S, et al. Estimates of global seasonal influenza-associated respiratory mortality: a modelling study. Lancet. 2018;391(10127):1285–1300. doi: 10.1016/S0140-6736(17)33293-2 29248255

2. Rolfes MA, Foppa IM, Garg S, Flannery B, Brammer L, Singleton JAet al. Annual estimates of the burden of seasonal influenza in the United States: A tool for strengthening influenza surveillance and preparedness. Influenza Other Respir Viruses. 2018;12(1):132–137. doi: 10.1111/irv.12486 29446233

3. Cassini A, Colzani E, Pini A, Mangen MJ, Plass D, McDonald SA, et al. Impact of infectious diseases on population health using incidence-based disability-adjusted life years (DALYs): results from the Burden of Communicable Diseases in Europe study, European Union and European Economic Area countries, 2009 to 2013. Euro Surveill. 2018;23(16).

4. Lee VJ, Ho ZJM, Goh EH, Campbell H, Cohen C, Cozza V, et al. Advances in measuring influenza burden of disease. Influenza Other Respir Viruses. 2018;12(1):3–9. doi: 10.1111/irv.12533 29460425

5. Paul Glezen W, Schmier JK, Kuehn CM, Ryan KJ, Oxford J. The burden of influenza B: a structured literature review. Am J Public Health. 2013;103(3):e43–51. doi: 10.2105/AJPH.2012.301137 23327249

6. Tafalla M, Buijssen M, Geets R, Vonk Noordegraaf-Schouten M. A comprehensive review of the epidemiology and disease burden of Influenza B in 9 European countries. Hum Vaccin Immunother. 2016;12(4):993–1002. doi: 10.1080/21645515.2015.1111494 26890005

7. Rota PA, Wallis TR, Harmon MW, Rota JS, Kendal AP, Nerome K. Cocirculation of two distinct evolutionary lineages of influenza type B virus since 1983. Virology. 1990;175(1):59–68. doi: 10.1016/0042-6822(90)90186-u 2309452

8. Paiva TM, Benega MA, Silva DB, Santos KC, Cruz AS, Hortenci MF, et al. Evolutionary pattern of reemerging influenza B/Victoria lineage viruses in São Paulo, Brazil, 1996–2012: Implications for vaccine composition strategy. J Med Virol. 2013;85(11):1983–9. doi: 10.1002/jmv.23684 23926069

9. Jennings L, Huang QS, Barr I, Lee PI, Kim WJ, Buchy P, et al. Literature review of the epidemiology of influenza B disease in 15 countries in the Asia-Pacific region. Influenza Other Respir Viruses. 2018;12(3):383–411. doi: 10.1111/irv.12522 29127742

10. Mertz D, Kim TH, Johnstone J, Lam PP, Science M, Kuster SP, et al. Populations at risk for severe or complicated influenza illness: systematic review and meta-analysis. BMJ. 2013;347:f5061. doi: 10.1136/bmj.f5061 23974637

11. Caini S, Schellevis F, El-Guerche Séblain C, Paget J. Important changes in the timing of influenza epidemics in the WHO European Region over the past 20 years: virological surveillance 1996 to 2016. Euro Surveill. 2018;23(1).

12. Tricco AC, Chit A, Soobiah C, Hallett D, Meier G, Chen MH, et al. Comparing influenza vaccine efficacy against mismatched and matched strains: a systematic review and meta-analysis. BMC Med. 2013;11:153. doi: 10.1186/1741-7015-11-153 23800265

13. Heikkinen T, Ikonen N, Ziegler T. Impact of influenza B lineage-level mismatch between trivalent seasonal influenza vaccines and circulating viruses, 1999–2012. Clin Infect Dis. 2014;59(11):1519–24. doi: 10.1093/cid/ciu664 25139969

14. Caini S, Huang QS, Ciblak MA, Kusznierz G, Owen R, Wangchuk S, et al. Epidemiological and virological characteristics of influenza B: results of the Global Influenza B Study. Influenza Other Respir Viruses. 2015;9 Suppl 1:3–12.

15. de Boer PT, van Maanen BM, Damm O, Ultsch B, Dolk FCK, Crépey P, et al. A systematic review of the health economic consequences of quadrivalent influenza vaccination. Expert Rev Pharmacoecon Outcomes Res. 2017;17(3):249–265. doi: 10.1080/14737167.2017.1343145 28613092

16. Sočan M, Prosenc K, Učakar V, Berginc N. A comparison of the demographic and clinical characteristics of laboratory-confirmed influenza B Yamagata and Victoria lineage infection. J Clin Virol. 2014;61(1):156–60. doi: 10.1016/j.jcv.2014.06.018 25034374

17. Le TT, Pham TH, Pham TH, Nguyen le KH, Nguyen CT, Hoang VM, et al. Circulation of influenza B lineages in northern Viet Nam, 2007–2014. Western Pac Surveill Response J. 2015;6(4):17–23. doi: 10.5365/WPSAR.2015.6.1.022 26798557

18. Xu C, Chan KH, Tsang TK, Fang VJ, Fung RO, Ip DK, et al. Comparative epidemiology of influenza B Yamagata- and Victoria-lineage viruses in households. Am J Epidemiol. 2015;182(8):705–13. doi: 10.1093/aje/kwv110 26400854

19. Barr IG, Vijaykrishna D, Sullivan SG. Differential age susceptibility to influenza B/Victoria lineage viruses in the 2015 Australian influenza season. Euro Surveill. 2016;21(4).

20. Seleka M, Treurnicht FK, Tempia S, Hellferscee O, Mtshali S, Cohen AL, et al. Epidemiology of influenza B/Yamagata and B/Victoria lineages in South Africa, 2005–2014. PLoS One. 2017;12(5):e0177655. doi: 10.1371/journal.pone.0177655 28542324

21. Moa AM, Muscatello DJ, Turner RM, MacIntyre CR. Epidemiology of influenza B in Australia: 2001–2014 influenza seasons. Influenza Other Respir Viruses. 2017;11(2):102–109. doi: 10.1111/irv.12432 27650482

22. Lapinscki B, Pereira LA, Nogueira MB, Vidal LR, Riediger I, Debur MC, et al. Molecular epidemiology of influenza B virus and implications in immunization strategy, Southern Brazil. Vaccine. 2018;36(1):107–113. doi: 10.1016/j.vaccine.2017.11.033 29174679

23. Caini S, Andrade W, Badur S, Balmaseda A, Barakat A, Bella A, et al. Temporal patterns of influenza A and B in tropical and temperate countries: what are the lessons for influenza vaccination? PLoS One. 2016;11(3):e0152310. doi: 10.1371/journal.pone.0152310 27031105

24. Caini S, Spreeuwenberg P, Kusznierz GF, Rudi JM, Owen R, Pennington K, et al. Distribution of influenza virus types by age using case-based global surveillance data from twenty-nine countries, 1999–2014. BMC Infect Dis. 2018;18(1):269. doi: 10.1186/s12879-018-3181-y 29884140

25. Central Intelligence Agency. The World Factbook. Available from Cited 20 June 2018

26. Alonso WJ, McCormick BJ. EPIPOI: a user-friendly analytical tool for the extraction and visualization of temporal parameters from epidemiological time series. BMC Public Health. 2012;12(1):982.

27. World Health Organization. Recommendations on the composition of influenza virus vaccines. Cited 18 June 2018

28. World Health Organization. Which vaccine formulation to use–Northern or Southern hemisphere? Cited 18 June 2018

29. Hirve S, Lambach P, Paget J, Vandemaele K, Fitzner J, Zhang W. Seasonal influenza vaccine policy, use and effectiveness in the tropics and subtropics—a systematic literature review. Influenza Other Respir Viruses. 2016;10(4):254–67. doi: 10.1111/irv.12374 26842617

30. Mosnier A, Caini S, Daviaud I, Bensoussan JL, Stoll-Keller F, Bui TT, et al. Ten influenza seasons in France: distribution and timing of influenza A and B circulation, 2003–2013. BMC Infect Dis. 2015;15:357. doi: 10.1186/s12879-015-1056-z 26289794

31. Barros EN, Cintra O, Rossetto E, Freitas L, Colindres R. Patterns of influenza B circulation in Brazil and its relevance to seasonal vaccine composition. Braz J Infect Dis. 2016;20(1):81–90. doi: 10.1016/j.bjid.2015.09.009 26626166

32. Pan Y, Zhang Y, Yang P, Qian H, Shi W, Wu S, et al. Epidemiological and phylogenetic characteristics of influenza B infection in severe acute respiratory infection cases in Beijing, 2014 to 2015. Medicine (Baltimore). 2015;94(52):e2399.

33. van Asten L, Bijkerk P, Fanoy E, van Ginkel A, Suijkerbuijk A, van der Hoek W, et al. Early occurrence of influenza A epidemics coincided with changes in occurrence of other respiratory virus infections. Influenza Other Respir Viruses. 2016;10(1):14–26. doi: 10.1111/irv.12348 26369646

34. Caini S, Alonso WJ, Séblain CE, Schellevis F, Paget J. The spatiotemporal characteristics of influenza A and B in the WHO European Region: can one define influenza transmission zones in Europe? Euro Surveill. 2017;22(35).

35. Laurie KL, Guarnaccia TA, Carolan LA, Yan AW, Aban M, Petrie S, et al. Interval between infections and viral hierarchy are determinants of viral interference following influenza virus infection in a ferret model. J Infect Dis. 2015;212(11):1701–10. doi: 10.1093/infdis/jiv260 25943206

36. Belshe RB, Coelingh K, Ambrose CS, Woo JC, Wu X. Efficacy of live attenuated influenza vaccine in children against influenza B viruses by lineage and antigenic similarity. Vaccine. 2010;28(9):2149–56. doi: 10.1016/j.vaccine.2009.11.068 20003926

37. de Boer PT, Kelso JK, Halder N, Nguyen TP, Moyes J, Cohen C, et al. The cost-effectiveness of trivalent and quadrivalent influenza vaccination in communities in South Africa, Vietnam and Australia. Vaccine. 2018;36(7):997–1007. doi: 10.1016/j.vaccine.2017.12.073 29373192

38. Tamerius JD, Shaman J, Alonso WJ, Bloom-Feshbach K, Uejio CK, Comrie A, et al. Environmental predictors of seasonal influenza epidemics across temperate and tropical climates. PLoS Pathog. 2013;9(3):e1003194. doi: 10.1371/journal.ppat.1003194 23505366

39. Lowen AC, Steel J. Roles of humidity and temperature in shaping influenza seasonality. J Virol. 2014;88(14):7692–5. doi: 10.1128/JVI.03544-13 24789791

40. Vijaykrishna D, Holmes EC, Joseph U, Fourment M, Su YC, Halpin R, et al. The contrasting phylodynamics of human influenza B viruses. Elife. 2015;4:e05055. doi: 10.7554/eLife.05055 25594904

41. Nyirenda M, Omori R, Tessmer HL, Arimura H, Ito K. Estimating the lineage dynamics of human influenza B viruses. PLoS One. 2016;11(11):e0166107. doi: 10.1371/journal.pone.0166107 27829058

42. Skowronski DM, Hottes TS, De Serres G, Ward BJ, Janjua NZ, Sabaiduc S, et al. Influenza B/Victoria antigen induces strong recall of B/Yamagata but lower B/Victoria response in children primed with two doses of B/Yamagata. Pediatr Infect Dis J. 2011;30(10):833–9. doi: 10.1097/INF.0b013e31822db4dc 21857263

43. Adlhoch C, Snacken R, Melidou A, Ionescu S, Penttinen; The European Influenza Surveillance Network. Dominant influenza A(H3N2) and B/Yamagata virus circulation in EU/EEA, 2016/17 and 2017/18 seasons, respectively. Euro Surveill. 2018;23(13).

Článok vyšiel v časopise


2019 Číslo 9