Identification of Mx gene nucleotide dimorphism (G/A) as genetic marker for antiviral activity in Egyptian chickens

Egyptian chickens, representing 2 breeds and 7 strains, were genotyped using the PCR-RFLP and sequencing techniques for detection of a non-synonymous dimorphism (G/A) in exon 14 of chicken Myxovirus resistance (Mx) gene. This dimorphic position is responsible for altering Mx protein’s antiviral activity. Polymerase Chain reactions were performed using Egyptian chickens DNA and specific primer set to amplify Mx DNA fragments of 299 or 301?bp, containing the dimorphic position. Amplicons were cut with restriction enzyme Hpy81. Genotype and allele frequencies for the resistant allele A and sensitive allele G were calculated in all the tested chickens. Results of PCR-RFLP were confirmed by sequencing. The three genotypes AA, AG, GG at the target nucleotide position in Mx gene were represented in all the studied Egyptian chicken breeds and strains except Baladi strain which showed only one genotype AA. The average allele frequency of the resistant A allele in the tested birds (0.67) was higher than the sensitive G allele average frequency in the same birds (0.33). Appling PCR-RFLP technique in the breeding program can be used to select chickens carrying the A allele with high frequencies. This will help in improving poultry breeding in Egypt by producing infectious disease-resistant chickens.     

Different genetic patterns in avian Toll-like receptor (TLR)5 genes

Toll-like receptors (TLRs) mediate immune response via recognition of pathogen-associated molecular patterns (PAMPs), thus play important roles in host defense. Polymorphisms of TLR5 may affect their recognition of bacterial flagellin, leading to varied host resistance to pathogenic infections. Here, we cloned TLR5 genes from Common Pheasant, Guinea fowl and 9 Chicken breeds and analyzed their sequences. The open reading frames of TLR5 were sequenced. Amino acid analysis indicated that TLR5 from Chicken breeds shared 99.4–99.9% homology. The amino acid homology of TLR5 ranged from 92.1 to 92.5% between Chickens and Guinea fowl, 95.7–96.1% between Chickens and Turkey, 94.3–94.7% between Chickens and Common Pheasant, and 79.9–80.1% between Chickens and Zebra-finch. Different genetic patterns were determined among Chickens, Common Pheasant, Guinea fowl, Turkey and Zebra-finch. It was found that there were 92 amino acid polymorphic sites, among which 5 sites in chicken TLR5, 63 sites in Guinea fowl TLR5 and 44 sites in Common Pheasant TLR5. Our data indicate that the positive Darwinian selection occurred in avian TLR5 genes since frequency of non-synonymous (d _N ) > frequency of synonymous (d _S ). These results also demonstrate that avian TLR5 genes are polymorphic among avian breeds, suggesting a varied resistance among breeds of avian. This information might be of help to improve the health of avian by breeding and vaccination.     

Natural selection in the <Emphasis Type="Italic">TLR-</Emphasis>related genes in the course of primate evolution

The innate immune system constitutes the front line of host defense against pathogens. Toll-like receptors (TLRs) recognize molecules derived from pathogens and play crucial roles in the innate immune system. Here, we provide evidence that the TLR-related genes have come under natural selection pressure in the course of primate evolution. We compared the nucleotide sequences of 16 TLR-related genes, including TLRs (TLR1–10), MYD88, TILAP, TICAM1, TICAM2, MD2, and CD14, among seven primate species. Analysis of the non-synonymous/synonymous substitution ratio revealed the presence of both strictly conserved and rapidly evolving regions in the TLR-related genes. The genomic segments encoding the intracellular Toll/interleukin 1 receptor domains, which exhibited lower rates of non-synonymous substitution, have undergone purifying selection. In contrast, TLR4, which carried a high proportion of non-synonymous substitutions in the part of extracellular domain spanning 200 amino acids, was found to have been the suggestive target of positive Darwinian selection in primate evolution. However, sequence analyses from 25 primate species, including eight hominoids, six Old World monkeys, eight New World monkeys, and three prosimians, showed no evidence that the pressure of positive Darwinian selection has shaped the pattern of sequence variations in TLR4 among New World monkeys and prosimians. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Three types of polymorphisms in exon 14 in porcine Mx1 gene

Much is known about the antiviral activity of Mx proteins in species such as mouse and human. In the mouse, loss of resistibility to influenza virus has been shown to be due to specific polymorphisms in the Mx gene. This gene is therefore an interesting candidate gene for disease resistance in farm animals. The porcine Mx1 gene has already been identified and characterized based on its homology with mouse Mx1; however, until now no evidence of polymorphisms in the porcine gene has been reported. In this study, we have found two new polymorphisms in exon 14 of porcine Mx1 by DNA sequencing and confirmed their presence in different breeds, using polymerase chain reaction (PCR)–restriction fragment length polymorphisms (RFLP) with NarI and NaeI restriction enzymes. On the basis of the deduced amino acid sequence, one allele contains a deletion that may result in a frameshift to yield several amino acid substitutions and extension of the carboxyl terminal region of Mx1 protein. The deletion allele, Mx1 ^c, was found to be segregating in Landrace, Berkshire, Duroc, Hampshire, and Yucatan miniature pig. A second point mutation, Mx1 ^b, was detected in Meishan and two Vietnamese native pig breeds. All other breeds tested were fixed for the Mx1 ^a allele that is identical to the sequence reported previously. It will be interesting to determine if the Mx1 ^c deletion is associated with variation in resistance to the myxovirus family in the pig.     

Flesh color association with polymorphism of the tyrosinase gene in different Chinese chicken breeds

In order to study genetic variation of tyrosinase gene in four different flesh color chicken breeds selected from special districts including Guyuan, Wenchang, Tibetan and Hisex chicken, five loci of the TYR gene exon-1 and one locus of 5′ flanking region were analyzed in PCR-SSCP and DNA sequencing. The results indicated that there were polymorphisms only at TYR1 and TYR3 locus. At TYR1 locus located in exon-1, there were three genotypes (TT, CC, TC), respectively, in three Chinese chicken breeds, and Genotype CC had not been detected in Hisex chicken. At TYR3 locus located in 5′ flanking region, there were three genotypes (GG, AA and GA) in Chinese local chicken breeds and genotype AA had not been detected in Hisex chicken breed. It was concluded that there were many variations of TYR gene in Chinese local chicken breeds. DNA sequencing of PCR products for different genotypes showed that there were two mutation sites, respectively, C to T at TYR1 locus and G to A at TYR3 locus. Mutation at TYR1 locus did not cause any amino acid variation. The chi-square analysis revealed that there were significant statistical differences generally between flesh color and the two loci among four chicken populations (P < 0.01). Our results suggested that the flesh color was related to genotype of TYR gene in Chinese chicken breeds. This study provided original information for elucidating the possible roles of exon-1 of TYR gene and 5′ flanking region in chickens with different flesh color chicken.     

Molecular evolution of bovine Toll-like receptor 2 suggests substitutions of functional relevance

Background  

There is accumulating evidence that polymorphism in Toll-like receptor (TLR) genes might be associated with disease resistance or susceptibility traits in livestock. Polymorphic sites affecting TLR function should exhibit signatures of positive selection, identified as a high ratio of non-synonymous to synonymous nucleotide substitutions (ω). Phylogeny based models of codon substitution based on estimates of ω for each amino acid position can therefore offer a valuable tool to predict sites of functional relevance. We have used this approach to identify such polymorphic sites within the bovine TLR2 genes from ten Bos indicus and Bos taurus cattle breeds. By analysing TLR2 gene phylogeny in a set of mammalian species and a subset of ruminant species we have estimated the selective pressure on individual sites and domains and identified polymorphisms at sites of putative functional importance.     

Purifying selection and positive selection on the myxovirus resistance gene in mammals and chickens

The myxovirus resistance gene (Mx) expresses antiviral activity in many species, e.g. mouse, human and chicken. It is not clear if the antiviral activity of Mx has evolved in these species to inhibit a set of species-specific pathogens, nor what factors drive Mx evolution in different animal lineages. Therefore, it is important to determine the evolutionary pattern of Mx and positively selected sites which affect the antiviral activity of the Mx gene in mammals and birds. We used sequence comparisons among species to detect positively selected sites by conducting phylogenetic analysis. The two-ratio model was significantly better than the one-ratio model in four species (mouse, rat, chicken and duck, p<0.05). Although selection pressure varied among different lineages, Mx had strong purifying selection in mammals and positive selection in chicken and duck lineages. Relative rate test revealed that Mx evolved faster in chickens than in ducks (Tajima's relative rate test, chi(2)=7.17, p<0.01). In the further analysis using a branch-site model A test, 8 sites were positively selected in the chicken lineage while no positive selection signals were observed for any site in the other lineages. The branch-site model A test had a omega value of 4.374 for the chicken lineage (2Deltal=14.20, d.f.=1, p<0.001). Comparisons of all currently available Mx mRNA sequences showed that these predicted positively selected sites had been fixed in the chicken lineage, suggesting that the chicken Mx gene evolved within the species to resist newly challenging environments. There is an increased selection constraint leading to mammals, while positive selection has acted on the chicken Mx.     

Bovine and water buffalo <Emphasis Type="Italic">Mx2</Emphasis> genes: polymorphism and antiviral activity

Millennia-long selective pressure of single-strand RNA viruses on the bovine Mx locus has increased the advantages of using the bovine Mx protein to evaluate the ultimate significance of the antiviral role of Mx proteins. The conclusions of research based only on the bovine Mx1 protein showed the need for comprehensive studies that demonstrate the role of all isoforms, individually or together, especially in the presence of a second isoform, the bovine Mx2 gene. This study provides information about bovine and water buffalo Mx2 genes, as well as their allelic polymorphism and basic antiviral potential. Observation of an Mx2 cDNA sequence (2,381 bp) obtained from 15 animals from 11 breeds using primers based on a previous sequence (NCBI accession no. AF335147) revealed several nucleotide substitutions, with eight different alleles and two amino acid exchanges: Gly to Ser at position 302 and Ile to Val at position 354, though the latter was found only in the NCBI database. A water buffalo Mx2 cDNA sequence was identified for the first time, revealing 46 nucleotide substitutions with 12 amino acid variations, in addition to a 9-bp insertion in the 5′ untranslated region UTR, compared with the bovine Mx2 cDNA. Transfected 3T3 cells expressing bovine Mx2 mRNAs coding Gly or Ser at position 302, water buffalo Mx2 mRNA, positive control bovine Mx1 mRNA-expressing cells, and negative control parental 3T3 were subjected to infection with recombinant vesicular stomatitis virus (VSVΔG*-G), as were empty pCI-neo vector-transfected cells. The positive control and all cells expressing Mx2 mRNAs displayed significantly higher levels of antiviral activity against VSVΔG*-G (P < 0.01) than did the negative controls.     

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

MHC class II variation in the endangered European mink <Emphasis Type="Italic">Mustela lutreola</Emphasis> (L. 1761)—consequences for species conservation

The polymorphic major histocompatibility complex (MHC) has gained a specific relevance in pathogen resistance and mate choice. Particularly the antigen-binding site (ABS), encoded by exon 2 of the DRB class II gene, exhibits numerous alleles and extensive sequence variations between alleles. A lack of MHC variability has attributed to instances such as bottleneck effects or relaxed selection pressure and has a certain impact on the long-term viability of the species concerned. As a result of seriously decreased population density during the last century, the current population of the endangered European mink (Mustela lutreola, L. 1761) has suffered from geographic isolation. In this study, we amplified a partial sequence of the MHC class II DRB exon 2 (229 bp), assessed the degree of genetic variation and compared the variability with those of other Mustelidae. As a result, nine alleles were detected in 20 investigated individuals, which differ from each other by four to 25 nucleotide substitutions (two to 11 amino acid substitutions). Whilst an equal ratio for synonymous and non-synonymous substitutions was found inside the ABS, synonymous substitutions were significantly higher than non-synonymous substitutions in the non-ABS region. Results might indicate that no positive selection exists within the ex situ population of M. lutreola, at least in the analysed fragment. In addition, phylogenetic analyses support the trans-species model of evolution. Becker and Nieberg have contributed equally to this work.     

藏鸡主要组织相容性复合体B-LBⅡ基因第二外显子序列遗传变异分析

根据鸡主要组织相容性复合体B-LBⅡ基因序列设计特异性引物,在藏鸡基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,并通过克隆和PCR直接测序获得了该片段的核苷酸序列。发现了15个B-LBⅡ新等位基因。对18个B-LBⅡ等位基因核苷酸序列和其所编码的MHCB-LBⅡ分子β1结构域的氨基酸序列分析显示,第二外显子核苷酸序列遗传多态性异常丰富,存在着62个多态变异位点(共包括80个突变),其中41个为简约性多态位点;衡量该序列遗传多样性的π值为0.0718;反映其群体内遗传变异度的平均遗传距离为0.056±0.008,低于在5个外来品种所估算的平均遗传距离。该编码区核苷酸相对异义替换率(15.61±2.69%)显著高于其同义替换率(3.25±0.94%),进一步分析表明,基因重组和平衡选择机制可能是引起B-LBⅡ基因序列变异的主要因素。在β1结构域氨基酸序列中,存在11个同义替换和27个异义替换;在24个肽结合位点中有12个变异位点;与其他6个中国地方鸡品种和一个外来品种比较发现,有11个异义氨基酸替换仅出现在藏鸡群体中,并被认为与藏鸡的免疫特异性有关,可为鸡的抗病力研究提供分子依据。     

Genomic structure, promoter analysis, and expression of the porcine (Sus scrofa) Mx1 gene

Allelic polymorphisms at the mouse Mx1 locus affect the probability of survival after experimental influenzal disease, raising the possibility that marker-assisted selection using the homologous locus could improve the innate resistance of pigs to natural influenza infections. Several issues need to be resolved before efficient large scale screening of the allelic polymorphism at the porcine (Sus scrofa) Mx1 locus can be implemented. First, the Mx1 genomic structure has to be established and sufficient flanking intronic sequences have to be gathered to enable simple PCR amplification of the coding portions of the gene. Then, a basic knowledge of the promoter region needs to be obtained as an allelic variation there can significantly alter absolute levels and/or tissue-specificity of MX protein expression. The results gathered here show that the porcine Mx1 gene and promoter share the major structural and functional characteristics displayed by their homologs described in cattle, mouse, chicken, and man. The crucial function of the proximal interferon-sensitive response elements motif for gene expression is also demonstrated. The sequence data compiled here will allow an extensive analysis of the polymorphisms present among the widest spectrum possible of porcine breeds with the aim to identify an Mx1 allele providing antiviral resistance.     

Reproductive tissue-specific expression profiling and genetic variation across a 19 gene bovine β-defensin cluster

β-defensins are small cationic peptides, with potent immunoregulatory and antimicrobial activity which are produced constitutively and inducibly by eukaryotic cells. This study profiles the expression of a cluster of 19 novel defensin genes which spans 320 kb on chromosome 13 in Bos taurus. It also assesses the genetic variation in these genes between two divergently selected cattle breeds. Using quantitative real-time PCR (qRT-PCR), all 19 genes in this cluster were shown to be expressed in the male genital tract and 9 in the female genital tract, in a region-specific manner. These genes were sequenced in Norwegian Red (NR) and Holstein-Friesian (HF) cattle for population genetic analysis. Of the 17 novel single nucleotide polymorphisms (SNPs) identified, 7 were non-synonymous, 6 synonymous and 4 outside the protein coding region. Significant frequency differences in SNPs in bovine β-defensins (BBD) 115, 117, 121, and 122 were detected between the two breeds, which was also reflected at the haplotype level (P < 0.05). There was clear segregation of the haplotypes into two blocks on chromosome 13 in both breeds, presumably due to historical recombination. This study documents genetic variation in this β-defensin gene cluster between Norwegian Red and Holstein-Friesian cattle which may result from divergent selection for production and fertility traits in these two breeds. Regional expression in the epididymis and fallopian tube suggests a potential reproductive-immunobiology role for these genes in cattle.     

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

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

Evolution of caprine and ovine β-defensin genes

Defensins comprise an important family of anti-microbial peptides. Among vertebrates, numerous defensin genes have been detected, but their evolutionary background is still discussed. We investigated the molecular evolution and variability of β-defensins of Caprini via sequence analyses of defensin introns. Screening of several domestic and wild species of Caprini revealed a total of 13 discrete β-defensin coding sequences, with three of them described before this study. Phylogenetic analyses revealed that the array of newly described defensin genes is of monophyletic origin and has arisen in numerous independent duplication events after separation of the ancestral defensins. As a result of that scenario, recent defensin genes are distributed in a species-specific manner. Values of synonymous and non-synonymous substitutions demonstrated that both modes of evolutionary pressure, positive as well as negative selection, have acted. In addition, conservation of some β-defensin exons is demonstrated. Discrimination of certain β-defensin genes was possible only due to intron-specific differences. Therefore, sequence analyses restricted to the exons would result in underestimation of the number of β-defensin genes. Our study shows that for reconstruction of the phylogenetic history data of defensin introns are more appropriated. Comparisons among the amino acid sequences show moderate substitutions without changing the net charge of the mature peptides. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

Genetic Variation within the Mx Gene of Commercially Selected Chicken Lines Reveals Multiple Haplotypes,Recombination and a Protein under Selection Pressure

The Mx protein is one of the best-characterized interferon-stimulated antiviral mediators. Mx homologs have been identified in most vertebrates examined; however, their location within the cell, their level of activity, and the viruses they inhibit vary widely. Recent studies have demonstrated multiple Mx alleles in chickens and some reports have suggested a specific variant (S631N) within exon 14 confers antiviral activity. In the current study, the complete genome of nine elite egg-layer type lines were sequenced and multiple variants of the Mx gene identified. Within the coding region and upstream putative promoter region 36 SNP variants were identified, producing a total of 12 unique haplotypes. Each elite line contained from one to four haplotypes, with many of these haplotypes being found in only one line. Observation of changes in haplotype frequency over generations, as well as recombination, suggested some unknown selection pressure on the Mx gene. Trait association analysis with either individual SNP or haplotypes showed a significant effect of Mx haplotype on several egg production related traits, and on mortality following Marek''s disease virus challenge in some lines. Examination of the location of the various SNP within the protein suggests synonymous SNP tend to be found within structural or enzymatic regions of the protein, while non-synonymous SNP are located in less well defined regions. The putative resistance variant N631 was found in five of the 12 haplotypes with an overall frequency of 47% across the nine lines. Two Mx recombinants were identified within the elite populations, indicating that novel variation can arise and be maintained within intensively selected lines. Collectively, these results suggest the conflicting reports in the literature describing the impact of the different SNP on chicken Mx function may be due to the varying context of haplotypes present in the populations studied.     

Discovery of a new nucleotide substitution in the MC1R gene and haplotype distribution in native and non-Japanese chicken breeds

Melanocortin 1-receptor (MC1R) is one of the major genes that controls chicken plumage colour. In this study, we investigated the sequence and haplotype distribution of the MC1R gene in native Japanese chickens, along with non-Japanese chicken breeds. In total, 732 and 155 chickens from 30 Japanese and eight non-Japanese breeds respectively were used. Three synonymous and 11 non-synonymous nucleotide substitutions were detected, resulting in 15 haplotypes (H0–H14). Of these, three were newly found haplotypes (H9, H13 and H14), of which one (H9) was composed of known substitutions C69T, T212C, G274A and G636A. The second one (H13) possessed newly found non-synonymous substitution C919G, apart from the known substitutions C69T, G178A, G274A, G636A and T637C. The third one (H14) comprised a newly discovered substitution C919G in addition to the known C69T, G274A and G409A substitutions. The homozygote for this new haplotype exhibited wt like plumage despite the presence of G274A. In addition to discovering a new nucleotide substitution (C919G) and three new haplotypes, we defined the plumage colour of the bird that was homozygous for the A644C substitution (H5 haplotype) as wheaten-like for the first time; although the substitution has been already reported, its effect was not revealed. Besides detecting the new plumage colour, we also confirmed that the A427G and G274A substitutions contribute in expressing brownish and black plumage colour respectively, as reported by the previous studies. Moreover, we confirmed that the buttercup allele does not express black plumage despite possessing a G274A substitution, under the suppression effect of A644C. In contrast, the birds homozygous for the birchen allele presented solid black plumage, which was contradictory to the previous reports. In conclusion, we revealed a large diversity in the MC1R gene of native Japanese chicken breeds, along with the discovery of a new non-synonymous nucleotide substitution (C919G) and three novel haplotypes (H9, H13 and H14).     

Patterns of variation of the major histocompatibility complex class IIB loci in Chinese goose (<Emphasis Type="Italic">Anser Cygnoides</Emphasis>)

In order to understand the variations of genomic organization of the major histocompatibility complex (MHC) and provide data for the studies on disease resistance of avian species, the MHC class II polymorphism in Chinese Z-goose was investigated for the first time in the present study. Eight alleles, which probably came from different loci, were found in six different geese with only one obvious band in the restriction fragment length polymorphism data. The numbers of nonsynonymous substitutions (dN) in peptide binding region of exon 2 were higher than that of synonymous substitutions (dS), and no stop codons or frameshift mutations were found in this region, indicating that balance selection was in operation, and the sequences are not likely to be pseudogenes. In addition, we successfully obtained five different long MHC class II fragments (about 1,162 bp) in six geese and found that the length of intron 1 was longer than that in chicken and some other birds, but intron 2 seemed to be intermediate in length. The phylogenetic tree appeared to branch in an order consistent with accepted evolutionary pathway. X. Zhou and C. Li contributed equally to this work and shared the co-first authorship.     

Heterogeneous Selective Pressure Acting on Influenza B Victoria- and Yamagata-Like Hemagglutinins

As a consequence of immune pressure, influenza virus hemagglutinin presents some of its amino acids under positive selection. Several authors have reported the existence of influenza A hemagglutinin codons under positive selective pressure (PSP). In this framework, the present work objectives were to demonstrate the presence of PSP and evaluate its effects on Victoria- and Yamagata-like influenza B viruses. Methodology adopted consisted in estimating the acceptance rate of nonsynonymous substitutions (ω = d_N/d_S) that describe the strength of selective pressure and identifying codons that may be positively selected, applying a set of continuous-time Markov chain codon-substitution models. Two groups of HA1 sequences (140 from Yamagata and 60 from Victoria lineage) were used. All the model maximum-likelihood estimates were obtained using codeml software application (PAML 3.15). The hypothesis of no existence of sites under PSP was rejected for both lineages (p < 0.001), using likelihood ratio tests. These results demonstrate the presence of positive selection acting on hemagglutinin of both Yamagata- and Victoria-like influenza B viruses. Several different sites were identified to be under PSP on Yamagata and Victoria hemagglutinins. Sites found with a posterior probability > 0.95 were codons 197 and 199 in both lineages, codon 75 in the Yamagata lineage, and codon 129 in the Victoria lineage. The detected amino acids are located at or near antigenic sites in influenza A virus H3 hemagglutinin. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.