Molecular characterization of MHC class I genes in Four species of Turdidae family to assess genetic diversity and selection.

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deal with intracellular pathogens (Virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene ﬂow. Here, we evaluated the MHC class I gene in four species of Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequence, including 47 putatively functional sequences and 30 pseudo sequences from 80 individuals. Method based on analysis of cloned amplicons (n=25) for each species, we found two and seven MHC I sequence per individual indicates more than one MHC I loci identiﬁed in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of diﬀerent selection patterns among the domains of PBR and non-PBR. Alleles found to be divergent with overall polymorphic sites per species ranging between 58 to 70 (out of 291 sites). Moreover, trans-species alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for existence of species. Our study signiﬁes a valuable tool for the characterization of evolutionary relevant diﬀerence across a population of birds with high conservational concerns.


INTRODUCTION
Major Histocompatibility Complex (MHC) are molecules encoded by certain genes that are most polymorphic to have been described in vertebrates' genomes [1].Two types of MHC gene families, class I and class II, which are useful cell surface glycoproteins that regulate the immune response.MHC class II molecules are heterodimers consists on α chain and a β chain, both contribute to presenting peptides from the processing of extracellular pathogens such as bacteria to the CD4+ T-helper cells [2].Heterodimers molecules of MHC class I are made up of an α chain and a non-MHC molecule, the β2 microglobulin.The α chain constitutes a cytoplasmic tail, a transmembrane domain, and three extracellular domains named α1, α2 and α3 [3] that are encoded by exons 2, 3 and 4. The MHC class I molecules are expressed in almost all somatic cells and trigger an adaptive immune response by presenting endogenously derived peptides of viral protein and an individual's own body cells to CD8+ cytotoxic T-cells [4].Polymorphism is largely confined within the region encoding the ABS (Antigen Binding Site) of the MHC class I [5].Maintenance of surprising diversity is supposed to be take place by two types of selection: heterozygote advantage and frequency-dependent selection.Heterozygotes could recognize a broader range of antigens from multiple pathogens and therefore have more fitness than either individual having homozygote [6].Other is frequency-dependent selection, in which rare alleles deliver a selective advantage where pathogens have found a means to escape against common immune defensive alleles in the population.Thus alteration in pathogens community with time and locality result in MHC variation in the host population.Generally, an individual possessing huge numbers and diverse MHC alleles; more pathogens it can recognize [1].
Structural diversity and immune response has been explored in numerous research, including genomics [7,8], ailment [9][10][11], mate choice [12][13][14].Sequence similarity at PBR based assignment to the locus is frequently hampered by various evolutionary indicators due to current recombination, duplication and/or concerted evolution as well as positive selection mediated by a variety of pathogens [15].Thus, numerous studies emphasized MHC genes as important markers to evaluate the adaptive potential and evolutionary status of threatened population [16].
The emerging scenario inspires researchers to collect statistics from a group of wild taxa to enlarge our understanding of the evolution of MHC gene [17].Despite significant efforts, protocols for locus-specific MHC genotyping in avian are still difficult to achieve and remarkably rare [18].MHC studies in population of wild birds remain neglected possibly due to complications in amplifying gene sequences from bird species not closely related to systematically studied chicken [19,20].
A significant decline in habitats and fragmentation of available habitats are predisposing factors for dramatic deterioration in population sizes [21].The avian genus Turdus is one of the broadly distributed passerine genera, with 65 documented extant species.The genus is listed wild territorial birds that are beneficial to china having economic and research value.Birds of this genus are strongly migratory thus experiencing a variety of environments.Up to the present, there are no studies on MHC class I genes in Turdidaespecies, which is the first step towards exploring the role of selection mediated by pathogens in the maintenance of MHC class I diversity.Precisely, this study aims to (1) characterize locus-specific variation in MHC I exon 3 genes across the Turdidae family and to examine the mode of evolution by which such variation comes about.To achieve this, we have measured the diversity and selection at MHC I genes to make available the variations that exist across Turdidae family.(2) We investigate the numbers of alleles possessed by each species and the general features of alleles in terms of functional genetic diversity.(3) Phylogenetic analyses to assess evolutionary relationships and processes driving avian MHC I diversity among four species of Turdidae family and other avian species.Results of this study provide us an insight into the evolutionary processes involved in avian MHC diversity, an important path for detailed studies on avian MHC evolution.

Study Population
The study population was non-sympatrically distributed 80 individuals of four species of genus Turdus of the Turdidae family.Samples include two to three contour feathers, tissue from breast and liver of birds accidentally injured or died during migratory season of 2017-19 in autumn and deposit in State key laboratory of wildlife detection center in northeast forestry university, stored at 4 °C.The sample material, voucher identity, and geographical location are presented (Table.S1; Figure 1).

Extraction of genomic DNA
Region of calamus to the rachis of contour feathers was excised, tissues from skeletal muscles were minced, placed into a 1.5 ml Eppendorf tube containing TNE buffer (10mM Tris-HCl (pH 8.0), 150mM NaCl, 2mM EDTA, 1% SDS).Total genomic DNA was extracted with AxyPrep Multisource Genomic DNA Miniprep Kit (AXYGEN, China) according to the manufacturer's instructions.The DNA concentration was measured with Nanopore Spectrophotometer at 260 nm absorbance.Samples above 100ng/μl concentration were used for further analysis.

PCR, cloning and sequencing
Polymerase chain reaction was conducted using motif specific primers designed for the amplification of MHC class I genes in great reed warbler.The forward primers HN36 5'-TCCCCACAGGTCTCCACACAGT-3' and HN46 reverse 5'-ATCCCAAATTCCCACCCACCTT-3' correspond to exon 3 region in the flanking introns, the region coding most of the peptide-binding site in MHC molecules (subunit α2) [22][23][24].The primers were used due to their successful amplification in many passerine species.A list of primers trialed is given in Table S2.Amplification was performed in the reaction mixture containing 20ng DNA template, 0.2μM of each primer, 25μl 2× EasyTaq® PCR SuperMix (+dye) (Trans, China), and water (deionized) to reach 50μl as final volume.Thermal cycling for MHC class I amplification began with one cycle at 94°C for 5 min, followed by 30 cycles of denaturation consisting of sequential steps of 94°C for 30s, 52°C for 30s, and 72°C for 30s, ending with a single extension step at 72°C for 5 min.Purification was carried out with AxyPrep TM DNA Gel Extraction Kit in accordance with the manufacturer's protocol.Purified PCR product was cloned using pEASY ®-T5 Zero Cloning Kit containing Trans1-T1 Phage resistant chemically competent cells (Transgen Biotech).PCRs were performed for positive clones using M13 forward and reverse primers.Several colonies (20)(21)(22)(23)(24)(25) per individual were selected and used as a template for sequencing directionally on an automatic sequencer (ABI PRISM 3730; Invitrogen Biotechnology Co. Ltd.).

Definition of Allele
Since few artifacts introduced during the recombination of PCR products in cloning [25,26].Amplification, cloning and sequencing were performed twice.Sequences were verified and referred to as an Allele; either minimum of three sequences have the same nucleotide composition or repeated in both events.The sequences which showed any deletion, insertion, or premature stop codons within exons were identified as presumed pseudogene sequence, and others were considered as putative functional allele (PFA) [27].All sequences appropriate to our criteria have been deposited into the GenBank (Accession No: MN849308-54).

Sequence Analysis
Chromatogram signals of all sequencing were examined with chromas 2.2.6.Sequences without ambiguous signals were selected.Vector sequence from the MHC class I gene was removed using seqMan in the DNA-Star7.1 package.Sequence editing and organization were done with BioEdit [28].Sequences were aligned individually and then altogether four sampled species using CLUSTAL X [29].The unique alleles were named according to the nomenclature for MHC in non-human species [30].NCBI BLAST [31] was used for sequences confirmation representing close identity to passerine species previously published MHC class I exon 3 sequences.Sequences having at least one stop codon (shift in the reading frame due to indels or nonsense sequences) were classified as pseudogenes.Based upon sequences found to be translatable, a minimum number of functional loci MHC class I was estimated using a conservational approach that all Loci from samples species' individual were in heterozygote state.

Tests for Selection
For selection, we conduct a priori classification of peptide binding region (PBR) and non-peptide region upon inferred passerine PBR sequences [40,41] homology sites with chicken MHC [42,43] and human HLA [44].The identification of sites subjected to selection in MHC class I Exon 3 was performed using various methods.The first standard selection test (Tajima's D, Fu & Li's F* and Fu & Li's D*) were calculated using DnaSP 5.0 [32].Second method was the calculation of parameter (ω) for functional alleles.It was carried out an overall estimation ofd N /d S of MHC class I Exon 3 and the other was codons comprising only PBR and non-PBR, which was calculated with MEGA 7.0 according to the Nei-Gojobori method [45] with the Jukes and Cantor correction.Standard error estimates were derived from 1000 bootstrap replicates.Z test of historical positive selection [46] was calculated in MEGA 7.0.Third, the Maximum likelihood implemented in codeml in PAML 4.9 was used for identification of sites involved in the positive selection, which are indicated where the ratio ω (d N /d S ) larger than 1 [47].Two different models corresponding ω were tested: M7 (beta), M8 (β and ω).To find whether the alternative model (M8) provided better fitter than the M7, we performed Likelihood ratio tests to compare twice the difference of the log-likelihood ratios (2ΔlnL) using a distributionχ 2 .PSSs in the M8 model was identified by PP more than 95% using the Bayes empirical Bayes procedure.Positively selected sites were verified at each codon site separately using many complementary approaches implemented in datamonkey (http://www.datamonkey.org/)[48] in addition to afore mention methods.Specifically, we used MEME [49], FEL, SALC [50] and FUBER [51]

Phylogenetic Analysis:
To assess the phylogenetic relationship, we construct two phylogenies (One for sampled species and other representing MHC class I sequences of related passerines plus sampled species) using Bayesian inference.We find the GTR+T nucleotide substitution model [52] that fits our data using MrModeltest [35] through the Akaike Information Criterion (AICc) [53].Bayesian Markov chain Monte Carlo (MCMC) was run for two million generations and sampling every 1,000 generations to ascertain when log Likelihood reached stationary phase.The phylogenetic tree was summarized in MrBayes v3.1.2[54] and the first 25% of the tree as burn-in was removed.Fig tree was used for visualization of the consensus tree.Exploration of relation between sampled species and related avian species, we conducted a maximum likelihood (ML) analysis with MEGA 7.0 [55].The data were analyzed with the T92+G model.We conducted 1000 bootstrap replicates to estimate the support.Values greater than 75% were indicated in the ML phylogenetic trees.The species covered are mainly fromPasseridae, Acrocephalidae, Paridae, Motacillidae, Muscicapidae, Hirundinidae, Phylloscopidae, Fringillidae, Cardinalidae, Sturnidae .Representative sequences from samples species and other selected for phylogenetic analysis were summarized in (See Table S3) with respective studied reference from GenBank.
To further identify allelic lineages among sampled species and related avian species, we conducted the Neighbor-Net algorithm in SplitsTree 4.14.8.Neighbor-Net networks were based on uncorrected P-distances and carried out 1000 bootstrap replicates to estimate nodal support.Nodal support values (> 75%) were displayed.

Characterization of Alleles
We successfully and selectively amplified MHC class I exon 3 genes across 80 individuals from four species of Turdidae family using HN36 and HN46 primers.An average of 22.7 clones per individual was sequenced.Sequences varied between 459 and 579 base pairs.The multiple sequence alignments of all sampled species were 411 base pair long.The final aligned MHC class I dataset included 285-291bp (Primers not include).Analysis of gDNA alignment revealed a total of 77 distinct Haplotypes/alleles including 47 PFA.Frequency, identity, distribution and the proposed name of haplotypes recovered from the individual of Turdidae family are given in (Supplementary Material Table S4).Each sequence was confirmed to exhibit similarity (81%-93%) with earlier reported passerine MHC class 1 sequences based upon BLAST search.The numbers of PFA sequences found in a single individual ranged from one to five, indicating that one to three loci exist in three of the four species of Turdidae family.However, the number of putative functional alleles found in a single individual ranged from two to seven in Turdus atrogularisexhibiting two to four loci.Number of the individual tested, number of PFA and pseudogene retrieved, the minimum number of functional loci estimated is given in Table.Of the 77 sequences, 30 were non-translatable due to indels or the presence of stop codons resulted changes in the reading frame.Sequences were thus presumed to be pseudogenes.The number of identified pseudogenes within the four species ranged between three to five in most individuals of study population and six of the thirteen pseudogene sequences were found to be identical in three individuals from the population of sampled species.We cannot ignore the likelihood that some of the identified pseudogene sequences may be due to PCR or sequencing artifacts, as such events would more often result in nonfunctional sequences.The nucleotide deletion result in loss of 3 amino acids was obvious in Tuna-MHCI*PS 07-9 and Tueu-MHCI*PS 01-04 andTueu-MHCI*PS 08.Both nucleotide deletion, frame shift mutation and premature stop codons were detected in Turu-MHCI*PS 01,03 andMHCI*PS 09 at amino acid 33 encoding Exon 3. Loss of 3 amino acids was at position 78 was detected in Tuat-MHCI*PS 05 andTuat-MHCI*PS06.

Analysis of genetic diversity
Overall we find an elevated genetic diversity (π) within exon 3 alleles repertoire among individuals of Turdus atrogularis was (0.151) than Turdus enumus (0.113).The average number of nucleotides difference (K ) varied between (43.95) in Turdus atrogularis to ( 32.32) in Turdus eunomus.TABLE 1 : Amplification success and genetic diversity within each of the four species of Turdidae family investigated.MHC class I Exon 3 gene size (L), the overall number of polymorphic sites per allele repertoire (S), and the average number of nucleotide differences (k ), Nucleotide diversity π at all sites, PBR and non-PBR.

Analysis of Recombination
The recombination detection program not only analyzes brake points but also identify parent sequences.We ran the test of recombination by pooling all putative functional alleles recovered from four species ofTurdidae family.We only find one potential recombination event in Tuna-MHCI*PFA06 inTurdus nau-manni at two recombinant breakpoints at position 148 and 253.Tuna-MHCI*PFA02 as major and Tuna-MHCI*PFA011 minor parent.Likewise, a single recombination was significant inTueu-MHCI*PFA07.We detected no recombination among other alleles.However, these recombinations were only significant in two out of seven tests and not consistent with recombination breakpoint identified by GARD, hence the results represent that overall recombination is not likely to have any prominent effects on tests for positively selected sites (Table .2).The recombination breakpoints identified by these two programs are often inconsistent, probably because they use different computational methods.+ Bold represents significant results.
Application of Likelihood models represents that the model M8 allows for positive selection provides a better than the neutral evolution models M7.Sites being positively selected were recognized, are given in (Table 4).In total, we find 12 codons under positive selection in sampled species, of which three sites (25%) match homologues codons found positively selected in other avian species and one (8.3%)matched human peptide binding region (Table 4).-847.53 -849.78P = 0.83, q = 0.211, P 0 = 0.585 P = 0.173, q = 0.279, P 1 = 0.17  Across all tests for positive selection, four codons (9,29,65 and 88) were frequently identified by all methods as having under positive selection.Of these, codons (42,59) were corresponding to PBR in human and codons 9, 29, 64 and 88 also match homology to PBR, known as positively selected among passerine in general [56] (Figure 4).The ten most frequent MHC class I alleles retrieved from sampled species displayed 87%-91% sequence similarity to 18 sequences from five other passerine families(Acrocephalidae, Passeridae, Muscicapidae, Paridae, Passerellidae).None of the 77 alleles studied had 100% sequence similarity to other published sequences to GenBank; thus, it establishes no allelic pair in the study population that was 100 % sequence likeness shared by another species.

Phylogenetic Analysis
In phylogenetic analysis, we observed that sampled species form a well-supported monophyletic clade with Erithacus rubeculs members of Turdidae family in maximum likelihood analysis (Figure S1).Bayesian analysis represents that most of the alleles shared amongTurdus atrogularis and Turdus reficollis .This pattern was almost consistent among Turdus naumanni and Turdus eunomuspresented in

DISCUSSION
In this study, we have for the first time characterize MHC Class I gene in four species of Turdidae family in the order Passeriformes from the wide geographical area of Northeast china.Analysis of MHC class I sequences revealed a total of 77 distinct Haplotypes/alleles including 47 putative functional alleles ever reported in passerine species, a group which is reported to have surprising MHC diversity [58,59].According to our findings based on MHC class I sequences, the functional loci in an individual ranged from one to three in three of the four species, which was consistent with findings from other passerine species studied till now [60].In addition, we detected a large number of presumed pseudogene sequences in the sampled population as it retains important information about the evolution of MHC.This is not surprising, as it is consistent with the expectation of evolution by birth-and-death [61].We made a significant effort to characterize the variation in regions of MHC class I exon 3 in our study population, we find that the primers would make some unlikely bias in allelic variations among individuals.Hence, MHC class I alleles variations per individual should, largly be due to copy number of genes variation among individuals, which has been confirmed in other birds [62].Few MHC class I alleles were shared between Turdus naumanni andTurdus eunomus as well as among individuals of Turdus ruficollis and Turdus atrogularis is indicating allelic sharing due to common ancestors or challenging common pathogens, as this event is frequent in numerous avian species such as owls, ardeid birds, penguins and passerines [63][64][65].
Generally, abundant variation in genetic material in a species is an indicator of the capacity to adapt to numerous environmental changes by that species.Rapidly evolving environmental pathogens would cause MHC genes to exhibit enlarged genetic diversity in species [66,67].Collectively, in our study, we find elevated genetic diversity among functional sequences and significant divergence, whereas pseudogene has low genetic variation and limited divergence.Similar results also have been described in other passerine species, including common yellow throat [68], great reed warbler and the great tit [69].The allelic variation described in our study could be due to increased immunological defense against the internal pathogen since these are highly unlikely to adapt to novel, infrequent variant [15].
Recombination has been considered an important mechanism that influences allelic diversity and driving evolution of the MHC gene [70,71] We only find one potential recombination event in Tuna-MHCI*PFA06 at two recombinant breakpoints at position 148 and 253 identified with recombination detection program.Similarly, single recombination was significant in Tueu-MHCI*PFA07.Recombination pattern was also restricted two out of seven tests; hence our finding indicate recombination is unlikely to have any significant influence on tests for PSs.Though we could not find any substantial recombination among other alleles, qualitatively our result suggests a role for recombination during the evolution of MHC class I in our species studied.Our finding is consistent with, that micro-recombination is frequently observed in MHC genes [57].
Further study of recombinant function in the future will contribute to a detailed understanding of its role in the evolution of the MHC gene.
Positive selection is the maintainer of alleles having the advantageous mutation that maintain fitness of an individual.In our study, the classical test of selection Tajima's D, Fu & Li's D* and Fu & Li's F* showed no deviation from neutral selection or balance selection.Considering the level of variation, conventional methods used to find selection are not influential [72].As sites positively selected are likely to accumulate more non-synonymous than synonymous substitutions, influencing amino acid variation to result in functional modifications in proteins [73] It should be noted that the pooling of all alleles across loci will mostly reduce selection detection tests, so the outcomes might be conservative, but will be less prone to false positives [77,78].Therefore, attention should be given while inferencing about the detected diversity in MHC and the possible effects of selection on individual loci.Our results suggested that α2 domain of MHC class I exon 3 of all species are under positive selection pressure.Pronounced positive selection at antigen-binding sites permits a species or population to present a larger repertoire of peptides (antigens), thus increase the defensive ability against parasitic and pathogenic infections.
Finally, phylogenetic clustering of MHC class I data set of sampled species when pooled with other passerine species produces a contrasting pattern.In general, the MHC class I sequence of Turdidae family clustered together with sequences from congeneric species.We found increased sequences similarities between same species rather than within species (trans specific likenesses), is usually described with trans species polymorphism (TSP), which occurs due to alleles passage from ancestral to the decedent via partial arrangement of lineages [79].Although trans specific similarities can be described with convergent evolution due to the results of similar environmental selective pressure.Studies indicate that TSP is a primary mechanism responsible for clustering of alleles at avian MHC class I [80].

CONCLUSION
Our study shows that species of Turdidae family has retained significant MHC class I diversity, which supports high conservational value and contribute to evolution of MHC class I genes.Importantly, we specifically amplify the exon 3 locus, provide an opportunity to avoid chimera formation during molecular characterization of hypervariable genes of immunity.At the same time, our study is the first to validate contrasting patterns of allelic diversity and positive selection upon inferred PBR and non-PBR codons supported the hypothesis that different mechanisms can shape evolutionary paths of MHC class I.

FIGURE 1 :
FIGURE 1: Geographical locations of samples included in our study.Square represents actual site of sample and size of square represents the approximate diameter of the sample's geographical range.

10 *
The log-likelihood values and parameters estimated were computed using codeml implemented in PAML 4.9.PSSs were inferred in model M8 by BEB with posterior probabilities > 95%.Usually consistent with the above finding, every substitute test (MEME, SALC, FEL, and FUBAR) for positive selection implemented in online adoptive evolutionary server Datamonkey (Weaver et al., 2018) identify numerous codons under positive selection (Figure2).

FIGURE 2 :
FIGURE 2: Positively selected sites using online adoptive evolutionary server Datamonkey with (a) MEME.(b) FEL.(c) FUBAR and (d) SALC.Substitution tests identify numerous codons showed the signature of positive selection.

Figure 3 :
Figure 3: Alignment of deduced amino acid sequences of all allele retrieved from four species of Turdidae family.Dots (.) indicate identity with the reference sequence.Light gray color represents codons presumably coding for peptide binding regions upon alignment with human and other avian species.Light brown region indicate flanking introns.(-) indicate missing nucleotides.( * ) Represents codons positively selected in almost all of test perform for selection analysis.(! ) indicates homologues region also positively selected in most of avian species.(# ) Inferred homology with human.

Figure 4 .
The Net network of putative functional and pseudogene MHC class I exon 3 sequences in Turdidae family with other passerines indicate that allelic distribution among them is almost congruent with limited divergence.For instance, Tueu-MHCI*PFA02 and Tuna-MHCI*PFA07 networks formed a monophyletic clade in the phylogenetic network of exon 3. Three alleles were shared amongTurdus naumanni and Turdus eunomus two among Turdus rufficollis and Turdus atrogularis.The clustering of the sequences among species could be due to trans-species polymorphism or orthology[57].

FIGURE 4 :
FIGURE 4: Bayesian phylogenetic reconstruction of MHC class I exon 3 of four species of Turdidae family.All the nodes are well supported (PP>0.90%)unless indicated otherwise.AB268885 Gallus gallus was used as an outgroup.

FIGURE 5 :
FIGURE 5: The phylogenetic networks of MHC class I of four species of turdidae family along with homologues sequences from passerines species.GenBank accession number are provided.Species name are mention lower right side with branch colors.Neighbor-Net networks based on uncorrected P -distances and carried out 1000 bootstrap replicates to estimate nodal support.Nodal support values (> 75%) were displayed.

TABLE 2 :
Recombinants detected in Turdidae family MHC class I alleles; parent sequences, breakpoints detected by the recombination detection program (RDP) and the genetic algorithm for recombination detection (GARD); and the RDP analyses.
Note: NS indicates not significant.Maj P and Min P represent as major and minor parents.BP denote Breakpoint.The numbers in parentheses are BP locations in the recombinant nucleotide sequences without gaps.The values after the semicolon are Max Chi values for that BPs.**indicateP < 0.01.Analysis of SelectionConsidering that the evolutionary history of each domain might have been different, we tested each domain separately for evidence of positive selection.Selection statistics by traditional methods did not disclose any statistical significant signal of selection that deviate from neutral expectations for turdus enumus (Tajima's D: -0.87309, P > 0.10; Fu & Li's D* test statistic: 0.36, P > 0.10; Fu & Li's F* test statistic: 0.03, P > 0.10) and turdus atrogularis (Tajima's D: -0.86107, P > 0.10 Fu & Li's D* test statistic: 0.19, P > 0.10; Fu & Li's F* test statistic: -0.077, P > 0.10).Still, overall d N value was significantly higher statistically than dS in Turdus atrogularis ( 1.687) and ratio d N/ ds was more pronounced at codons presumably coding PBR (1.994) than codons not involved in such activity (0.884) is presented in (Table3).

TABLE 3 :
The average rates of nonsynonymous (d N ), synonymous (ds) substitutions and the result of Z-test, the average nucleotide distances (d nt ) and amino acid distances (d aa ) for PBR, non-PBR and all sites in MHC class I of Turdidae family.

TABLE 4 :
Estimation of d N andd S substitution rates for sites positively selected and their ratio for codons chosen a priori to constitute (PBR) and not involved (non-PBR).
[76]r study revealed differential expression of selection pattern in functional sequences on regions related with PBR and non-PBR of the MHC class I gene.Codons involved in peptide binding region revealed more non-synonymous substitution than synonymous (d N /d S =1.99) in Turdus atrogularis as compared to non-peptide binding region (d N /d S =0.884), pattern was consistent among all species tested, which might be enlightened that stronger selection pressure from intracellular pathogens than extracellular pathogens[74].Evidence of positive selection at PBR of MHC has been reported in the house sparrow (PBR d N /d S = 1.55 vs non-PBR d N /d S = 0.51)[75]and golden pheasant (PBR d N /d S = 1.45 vs non-PBR d N /d S =0.91)[76].Of the 12 codons in total among species tested exhibit positive selection with Likelihood methods using PAML, 9,29, 64 and 88 match homologues codons found positively selected in other passerine species.