A computational framework for influenza antigenic cartography

Influenza viruses have been responsible for large losses of lives around the world and continue to present a great public health challenge. Antigenic characterization based on hemagglutination inhibition (HI) assay is one of the routine procedures for influenza vaccine strain selection. However, HI assay is only a crude experiment reflecting the antigenic correlations among testing antigens (viruses) and reference antisera (antibodies). Moreover, antigenic characterization is usually based on more than one HI dataset. The combination of multiple datasets results in an incomplete HI matrix with many unobserved entries. This paper proposes a new computational framework for constructing an influenza antigenic cartography from this incomplete matrix, which we refer to as Matrix Completion-Multidimensional Scaling (MC-MDS). In this approach, we first reconstruct the HI matrices with viruses and antibodies using low-rank matrix completion, and then generate the two-dimensional antigenic cartography using multidimensional scaling. Moreover, for influenza HI tables with herd immunity effect (such as those from Human influenza viruses), we propose a temporal model to reduce the inherent temporal bias of HI tables caused by herd immunity. By applying our method in HI datasets containing H3N2 influenza A viruses isolated from 1968 to 2003, we identified eleven clusters of antigenic variants, representing all major antigenic drift events in these 36 years. Our results showed that both the completed HI matrix and the antigenic cartography obtained via MC-MDS are useful in identifying influenza antigenic variants and thus can be used to facilitate influenza vaccine strain selection. The webserver is available at http://sysbio.cvm.msstate.edu/AntigenMap.     

Patterns of predicted T-cell epitopes associated with antigenic drift in influenza H3N2 hemagglutinin

Antigenic drift allowing escape from neutralizing antibodies is an important feature of transmission and survival of influenza viruses in host populations. Antigenic drift has been studied in particular detail for influenza A H3N2 and well defined antigenic clusters of this virus documented. We examine how host immunogenetics contributes to determination of the antibody spectrum, and hence the immune pressure bringing about antigenic drift. Using uTOPE™ bioinformatics analysis of predicted MHC binding, based on amino acid physical property principal components, we examined the binding affinity of all 9-mer and 15-mer peptides within the hemagglutinin 1 (HA1) of 447 H3N2 virus isolates to 35 MHC-I and 14 MHC-II alleles. We provide a comprehensive map of predicted MHC-I and MHC-II binding affinity for a broad array of HLA alleles for the H3N2 influenza HA1 protein. Each HLA allele exhibited a characteristic predicted binding pattern. Cluster analysis for each HLA allele shows that patterns based on predicted MHC binding mirror those described based on antibody binding. A single amino acid mutation or position displacement can result in a marked difference in MHC binding and hence potential T-helper function. We assessed the impact of individual amino acid changes in HA1 sequences between 10 virus isolates from 1968–2002, representative of antigenic clusters, to understand the changes in MHC binding over time. Gain and loss of predicted high affinity MHC-II binding sites with cluster transitions were documented. Predicted high affinity MHC-II binding sites were adjacent to antibody binding sites. We conclude that host MHC diversity may have a major determinant role in the antigenic drift of influenza A H3N2.     

Identification of Low- and High-Impact Hemagglutinin Amino Acid Substitutions That Drive Antigenic Drift of Influenza A(H1N1) Viruses

Determining phenotype from genetic data is a fundamental challenge. Identification of emerging antigenic variants among circulating influenza viruses is critical to the vaccine virus selection process, with vaccine effectiveness maximized when constituents are antigenically similar to circulating viruses. Hemagglutination inhibition (HI) assay data are commonly used to assess influenza antigenicity. Here, sequence and 3-D structural information of hemagglutinin (HA) glycoproteins were analyzed together with corresponding HI assay data for former seasonal influenza A(H1N1) virus isolates (1997–2009) and reference viruses. The models developed identify and quantify the impact of eighteen amino acid substitutions on the antigenicity of HA, two of which were responsible for major transitions in antigenic phenotype. We used reverse genetics to demonstrate the causal effect on antigenicity for a subset of these substitutions. Information on the impact of substitutions allowed us to predict antigenic phenotypes of emerging viruses directly from HA gene sequence data and accuracy was doubled by including all substitutions causing antigenic changes over a model incorporating only the substitutions with the largest impact. The ability to quantify the phenotypic impact of specific amino acid substitutions should help refine emerging techniques that predict the evolution of virus populations from one year to the next, leading to stronger theoretical foundations for selection of candidate vaccine viruses. These techniques have great potential to be extended to other antigenically variable pathogens.     

Computational Identification of Antigenicity-Associated Sites in the Hemagglutinin Protein of A/H1N1 Seasonal Influenza Virus

The antigenic variability of influenza viruses has always made influenza vaccine development challenging. The punctuated nature of antigenic drift of influenza virus suggests that a relatively small number of genetic changes or combinations of genetic changes may drive changes in antigenic phenotype. The present study aimed to identify antigenicity-associated sites in the hemagglutinin protein of A/H1N1 seasonal influenza virus using computational approaches. Random Forest Regression (RFR) and Support Vector Regression based on Recursive Feature Elimination (SVR-RFE) were applied to H1N1 seasonal influenza viruses and used to analyze the associations between amino acid changes in the HA1 polypeptide and antigenic variation based on hemagglutination-inhibition (HI) assay data. Twenty-three and twenty antigenicity-associated sites were identified by RFR and SVR-RFE, respectively, by considering the joint effects of amino acid residues on antigenic drift. Our proposed approaches were further validated with the H3N2 dataset. The prediction models developed in this study can quantitatively predict antigenic differences with high prediction accuracy based only on HA1 sequences. Application of the study results can increase understanding of H1N1 seasonal influenza virus antigenic evolution and accelerate the selection of vaccine strains.     

Antigenic drift in influenza virus H3 hemagglutinin from 1968 to 1980: multiple evolutionary pathways and sequential amino acid changes at key antigenic sites

Surveys of the antigenic properties of a wide range of variants of the H3N2 (Hong Kong) influenza virus subtype have revealed complex patterns of variants cocirculating during each of the main epidemic eras of the subtype. We determined hemagglutinin (HA) gene sequences for 14 isolates chosen to give the wildest possible spread of variant types. The addition of these data to existing HA gene sequence information for other variants provides a comprehensive picture of HA gene evolution during antigenic drift among H3N2 subtype viruses. The data reveal the existence of multiple evolutionary pathways during at least one period of development of the subtype and strikingly demonstrate that amino acid changes are limited to a small number of locations on the HA molecule during antigenic drift. The occurrence of sequential amino acid changes at key positions within these variable regions suggests that the HA structure has remained constant during subtype evolution so that only limited possibilities remain for further antigenic drift among H3N2 viruses.     

The antigenic property of the H5N1 avian influenza viruses isolated in central China

ABSTRACT: BACKGROUND: Three influenza pandemics outbroke in the last century accompanied the viral antigen shift and drift, resulting in the change of antigenic property and the low cross protective ability of the existed antibody to the newly emerged pandemic virus, and eventually the death of millions of people. The antigenic characterizations of the viruses isolated in central China in 2004 and 2006-2007 were investigated in the present study. RESULTS: Hemagglutinin inhibition assay and neutralization assay displayed differential antigenic characteristics of the viruses isolated in central China in two periods (2004 and 2006-2007). HA genes of the viruses mainly located in two branches in phylogeny analysis. 53 mutations of the deduced amino acids of the HA genes were divided into 4 patterns. Mutations in pattern 2 and 3 showed the main difference between viruses isolated in 2004 and 2006-2007. Meanwhile, most amino acids in pattern 2 and 3 located in the globular head of the HA protein, and some of the mutations evenly distributed at the epitope sites. CONCLUSIONS: The study demonstrated that a major antigenic drift had occurred in the viruses isolated in central China. And monitoring the antigenic property should be the priority in preventing the potential pandemic of H5N1 avian influenza virus.     

Impact of antigenic and genetic drift on the serologic surveillance of H5N2 avian influenza viruses

Background  

Serologic surveillance of Avian Influenza (AI) viruses is carried out by the hemagglutination inhibition (HI) test using reference reagents. This method is recommended by animal health organizations as a standard test to detect antigenic differences (subtypes) between circulating influenza virus, vaccine- and/or reference- strains. However, significant discrepancies between reference antisera and field isolates have been observed during serosurveillance of influenza A viruses in pig and poultry farms. The objective of this study was to examine the effects of influenza virus genetic and antigenic drift on serologic testing using standard HI assays and reference reagents. Low pathogenic AI H5N2 viruses isolated in Mexico between 1994 and 2008 were used for phylogenetic analysis of AI hemagglutinin genes and for serologic testing using antisera produced with year-specific AI virus isolates.     

Minor changes in the hemagglutinin of influenza A(H1N1)2009 virus alter its antigenic properties

Background

The influenza A(H1N1)2009 virus has been the dominant type of influenza A virus in Finland during the 2009–2010 and 2010–2011 epidemic seasons. We analyzed the antigenic characteristics of several influenza A(H1N1)2009 viruses isolated during the two influenza seasons by analyzing the amino acid sequences of the hemagglutinin (HA), modeling the amino acid changes in the HA structure and measuring antibody responses induced by natural infection or influenza vaccination.

Methods/Results

Based on the HA sequences of influenza A(H1N1)2009 viruses we selected 13 different strains for antigenic characterization. The analysis included the vaccine virus, A/California/07/2009 and multiple California-like isolates from 2009–2010 and 2010–2011 epidemic seasons. These viruses had two to five amino acid changes in their HA1 molecule. The mutation(s) were located in antigenic sites Sa, Ca1, Ca2 and Cb region. Analysis of the antibody levels by hemagglutination inhibition test (HI) indicated that vaccinated individuals and people who had experienced a natural influenza A(H1N1)2009 virus infection showed good immune responses against the vaccine virus and most of the wild-type viruses. However, one to two amino acid changes in the antigenic site Sa dramatically affected the ability of antibodies to recognize these viruses. In contrast, the tested viruses were indistinguishable in regard to antibody recognition by the sera from elderly individuals who had been exposed to the Spanish influenza or its descendant viruses during the early 20^th century.

Conclusions

According to our results, one to two amino acid changes (N125D and/or N156K) in the major antigenic sites of the hemagglutinin of influenza A(H1N1)2009 virus may lead to significant reduction in the ability of patient and vaccine sera to recognize A(H1N1)2009 viruses.     

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.     

Effect of vaccine use in the evolution of Mexican lineage H5N2 avian influenza virus

An outbreak of avian influenza (AI) caused by a low-pathogenic H5N2 type A influenza virus began in Mexico in 1993 and several highly pathogenic strains of the virus emerged in 1994-1995. The highly pathogenic virus has not been reported since 1996, but the low-pathogenic virus remains endemic in Mexico and has spread to two adjacent countries, Guatemala and El Salvador. Measures implemented to control the outbreak and eradicate the virus in Mexico have included a widespread vaccination program in effect since 1995. Because this is the first case of long-term use of AI vaccines in poultry, the Mexican lineage virus presented us with a unique opportunity to examine the evolution of type A influenza virus circulating in poultry populations where there was elevated herd immunity due to maternal and active immunity. We analyzed the coding sequence of the HA1 subunit and the NS gene of 52 Mexican lineage viruses that were isolated between 1993 and 2002. Phylogenetic analysis indicated the presence of multiple sublineages of Mexican lineage isolates at the time vaccine was introduced. Further, most of the viruses isolated after the introduction of vaccine belonged to sublineages separate from the vaccine's sublineage. Serologic analysis using hemagglutination inhibition and virus neutralization tests showed major antigenic differences among isolates belonging to the different sublineages. Vaccine protection studies further confirmed the in vitro serologic results indicating that commercial vaccine was not able to prevent virus shedding when chickens were challenged with antigenically different isolates. These findings indicate that multilineage antigenic drift, which has not been observed in AI virus, is occurring in the Mexican lineage AI viruses and the persistence of the virus in the field is likely aided by its large antigenic difference from the vaccine strain.     

Antigenic and genetic evolution of swine influenza A (H3N2) viruses in Europe

In the early 1970s, a human influenza A/Port Chalmers/1/73 (H3N2)-like virus colonized the European swine population. Analyses of swine influenza A (H3N2) viruses isolated in The Netherlands and Belgium revealed that in the early 1990s, antigenic drift had occurred, away from A/Port Chalmers/1/73, the strain commonly used in influenza vaccines for pigs. Here we show that Italian swine influenza A (H3N2) viruses displayed antigenic and genetic changes similar to those observed in Northern European viruses in the same period. We used antigenic cartography methods for quantitative analyses of the antigenic evolution of European swine H3N2 viruses and observed a clustered virus evolution as seen for human viruses. Although the antigenic drift of swine and human H3N2 viruses has followed distinct evolutionary paths, potential cluster-differentiating amino acid substitutions in the influenza virus surface protein hemagglutinin (HA) were in part the same. The antigenic evolution of swine viruses occurred at a rate approximately six times slower than the rate in human viruses, even though the rates of genetic evolution of the HA at the nucleotide and amino acid level were similar for human and swine H3N2 viruses. Continuous monitoring of antigenic changes is recommended to give a first indication as to whether vaccine strains may need updating. Our data suggest that humoral immunity in the population plays a smaller role in the evolutionary selection processes of swine H3N2 viruses than in human H3N2 viruses.     

一种含不同流感病毒血凝素抗原表位的核酸疫苗

流感病毒表面抗原血凝素( hemagglutinin,HA)是流感核酸疫苗重要的靶抗原,针对HA的保护性中和抗体主要由HA上的五个抗原表位诱导产生.在本文中,我们构建了一种以新甲型H1N1流感病毒HA1为骨架的含2个A/PR/8( H1N1)流感病毒HA抗原表位和3个新甲型H1N1流感病毒HA抗原表位的核酸疫苗,并在B...     

Antigenic drift in H5N1 avian influenza virus in poultry is driven by mutations in major antigenic sites of the hemagglutinin molecule analogous to those for human influenza virus

H5N1 highly pathogenic avian influenza virus has been endemic in poultry in Egypt since 2008, notwithstanding the implementation of mass vaccination and culling of infected birds. Extensive circulation of the virus has resulted in a progressive genetic evolution and an antigenic drift. In poultry, the occurrence of antigenic drift in avian influenza viruses is less well documented and the mechanisms remain to be clarified. To test the hypothesis that H5N1 antigenic drift is driven by mechanisms similar to type A influenza viruses in humans, we generated reassortant viruses, by reverse genetics, that harbored molecular changes identified in genetically divergent viruses circulating in the vaccinated population. Parental and reassortant phenotype viruses were antigenically analyzed by hemagglutination inhibition (HI) test and microneutralization (MN) assay. The results of the study indicate that the antigenic drift of H5N1 in poultry is driven by multiple mutations primarily occurring in major antigenic sites at the receptor binding subdomain, similarly to what has been described for human influenza H1 and H3 subtype viruses.     

长沙A(H1N1)亚型季节性流感暴发疫情的实验室诊断及其HA基因特性分析

对2009 年长沙麓山国际学校流感暴发疫情进行实验室诊断, 并探索新分离的A(H1N1)亚型流感病毒血凝素(HA)的基因特性。对流感暴发疫情的25 份鼻/咽拭子标本进行RT-PCR检测和流感病毒分离, 然后利用CEQ?8000 Genetic Analysis System对病毒分离株(A/Yuelu/314/2009)进行测序, 测序结果提交至GenBank(登录号: FJ912843)并用ClustalX和Mega4.1软件进行序列分析。结果显示, 分离出A(H1N1)亚型流感毒株18株, 检出21份A(H1N1)亚型流感病毒核酸阳性; A/Yuelu/314/2009(H1N1) HA基因序列与2008~2009 年疫苗株(A/Brisbane/59/2007)比较显示: 核苷酸和氨基酸同源性均为99%, 有6个位点的氨基酸发生了变异(V148A、S158N、G202A、I203D、A206T、W435R), 其中一个S158N氨基酸变异位于B抗原表位, HA基因序列上共有潜在糖基化位点9 个(27、28、40、71、151、176、303、497、536), 与A/Brisbane/59/2007相同且氨基酸序列保守。本实验诊断出此次流感暴发疫情的病原体为A(H1N1)型季节性流感病毒, 研究还发现A/Yuelu/314/2009(H1N1)长沙分离株与A/Brisbane/59/2007 疫苗株基因序列比较显示并未形成一个新的变种, 推测是由于分离株与疫苗株之间基因特性的改变和人群对A(H1N1)亚型流感病毒免疫力降低导致了此次长沙麓山国际学校A(H1N1)亚型流感的暴发。     

Role of Positive Selection Pressure on the Evolution of H5N1 Hemagglutinin

The surface glycoprotein hemagglutinin (HA) helps the influenza A virus to evade the host immune system by antigenic variation and is a major driving force for viral evolution. In this study, the selection pressure on HA of H5N1 influenza A virus was analyzed using bioinformatics algorithms. Most of the identified positive selection (PS) sites were found to be within or adjacent to epitope sites. Some of the identified PS sites are consistent with previous experimental studies, providing further support to the biological significance of our findings. The highest frequency of PS sites was observed in recent strains isolated during 2005-2007. Phylogenetic analysis was also conducted on HA sequences from various hosts. Viral drift is almost similar in both avian and human species with a progressive trend over the years. Our study reports new mutations in functional regions of HA that might provide markers for vaccine design or can be used to predict isolates of pandemic potential.     

Molecular evolution of influenza A (H3N2) viruses circulated in Fujian Province, China during the 1996–2004 period

We studied the genetic and epidemic characteristics of influenza A (H3N2) viruses circulated in human in Fujian Province, south of China from 1996 to 2004. Phylogenetic analysis was carried out for genes encoding hemagglutinin1 (HA1) of influenza A virus (14 new and 11 previously reported reference se-quences). Our studies revealed that in the 8 flu seasons, the mutations of HA1 genes occurred from time to time, which were responsible for about four times of antigenic drift of influenza H3N2 viruses in Fujian, China. The data demonstrated that amino acid changes were limited to some key codons at or near antibody binding sites A through E on the HA1 molecule. The changes at the antibody binding site B or A or sialic acid receptor binding site 226 were critical for antigenic drift. But the antigenic sites might change and the key codons for antigenic drift might change as influenza viruses evolve. It seems important to monitor new H3 isolates for mutations in the positively selected codons of HA1 gene in south of Asia.     

Molecular characterization and phylogenetic analysis of H3N2 human influenza A viruses in Cheongju, South Korea

To investigate the genetic characteristics of human influenza viruses circulating in Chungbuk province, we tested 510 clinical samples of nasopharyngeal suction from pediatric patients diagnosed with respiratory illness between June 2007 and June 2008. Genetic characterization of the HA genes of H3N2 isolates indicated the relative higher similarity to A/Virginia/04/07 (99.6%) rather than that of A/Wisconsin/67/2005 (98.4%), a Northern Hemisphere 2007∼2008 vaccine strain, based on amino acid sequences. We found several altered amino acids at the H3 HA1 antigenic sites compared with the vaccine strain; K140I at site A, K158R at site B, and K173N (H471) or K173Q, and S262N at site E, but there was no antigenic shift among the H3N2 viruses. Interestingly, A/Cheongju/H383/08 and A/Cheongju/H407/08 isolates had single amino acid substitution at D151G on the catalytic site of the N2 NA while A/Cheongju/H412/08 and A/Cheongju/ H398/07 isolates had one amino acid deletion at residue 146. Furthermore, we found that 25% (3 out of 12 isolates) of the H3N2 subtype viruses had the amino acid substitution at position 31 on the M2 protein (Aspartic acid to Asparagine) and confirmed their drug-resistance by biological assays. Taken together, the results of this study demonstrated continuous evolutions of human H3N2 viruses by antigenic drift and also highlighted the need to closely monitor antigenic drug resistance in influenza A viruses to aid in the early detection of potentially pandemic strains, as well as underscore the need for new therapeutics.     

Immunodominance of Antigenic Site B over Site A of Hemagglutinin of Recent H3N2 Influenza Viruses

H3N2 influenza viruses have now circulated in the human population for 43 years since the pandemic of 1968, accumulating sequence changes in the hemagglutinin (HA) and neuraminidase (NA) that are believed to be predominantly due to selection for escape from antibodies. Examination of mutations that persist and accumulate led to identification of antigenically significant mutations that are contained in five antigenic sites (A-E) mapped on to the H3 HA. In early H3N2 isolates, antigenic site A appeared to be dominant while in the 1990s site B seemed more important. To obtain experimental evidence for dominance of antigenic sites on modern H3 HAs, we have measured antibodies in plasma of human subjects who received the 2006-07 trivalent subunit influenza vaccine (H3 component A/Wisconsin/67/05) or the 2008-09 formulation (H3 component A/Uruguay/716/07). Plasmas were tested against expressed HA of Wisconsin-like influenza A/Oklahoma/309/06 and site-directed mutants in antigenic site A (NNES121-124ITEG, N126T, N133D, TSSS135-138GSNA, K140I, RSNNS142-146PGSG), and antigenic site B (HL156-157KS, KFK158-160GST, NDQI189-192QEQT, A196V). "Native ELISA" analysis and escape mutant selection with two human monoclonal antibodies demonstrated that antibody E05 binds to antigenic site A and 1_C02 binds to site B. We find that most individuals, after vaccination in seasons 2006-07 and/or 2008-09, showed dominance of antigenic site B recognition over antigenic site A. A minority showed dominance of site A in 2006 but these were reduced in 2008 when the vaccine virus had a site A mutation. A better understanding of immunodominance may allow prediction of future antigenic drift and assist in vaccine strain selection.     

Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes

Phylogenetic profiles of the genes coding for the hemagglutinin (HA) protein, nucleoprotein (NP), matrix (M) protein, and nonstructural (NS) proteins of influenza B viruses isolated from 1940 to 1998 were analyzed in a parallel manner in order to understand the evolutionary mechanisms of these viruses. Unlike human influenza A (H3N2) viruses, the evolutionary pathways of all four genes of recent influenza B viruses revealed similar patterns of genetic divergence into two major lineages. Although evolutionary rates of the HA, NP, M, and NS genes of influenza B viruses were estimated to be generally lower than those of human influenza A viruses, genes of influenza B viruses demonstrated complex phylogenetic patterns, indicating alternative mechanisms for generation of virus variability. Topologies of the evolutionary trees of each gene were determined to be quite distinct from one another, showing that these genes were evolving in an independent manner. Furthermore, variable topologies were apparently the result of frequent genetic exchange among cocirculating epidemic viruses. Evolutionary analysis done in the present study provided further evidence for cocirculation of multiple lineages as well as sequestering and reemergence of phylogenetic lineages of the internal genes. In addition, comparison of deduced amino acid sequences revealed a novel amino acid deletion in the HA1 domain of the HA protein of recent isolates from 1998 belonging to the B/Yamagata/16/88-like lineage. It thus became apparent that, despite lower evolutionary rates, influenza B viruses were able to generate genetic diversity among circulating viruses through a combination of evolutionary mechanisms involving cocirculating lineages and genetic reassortment by which new variants with distinct gene constellations emerged.