Allelic diversity at the Mhc-DQA locus in cotton rats (Sigmodon hispidus) and a comparison of DQA sequences within the family Muridae (Mammalia: Rodentia)

 The cotton rat (Sigmodon hispidus) is a common murid rodent of the southern United States, Mexico, and Central America. Using single-stranded conformation polymorphism analysis and DNA sequencing techniques, 11 DQA exon 2 alleles were detected among 180 S. hispidus from Caddo County, Oklahoma, USA. The alleles represent a single locus exhibiting a high level of polymorphism. Nucleotide and amino acid distance values among DQA alleles of S. hispidus were higher than those within Mus musculus and species of Rattus. Although the distribution of polymorphic amino acid residues among alleles of S. hispidus was similiar to that of Mus and Rattus, some residues of the α-helix region were more variable in S. hispidus. Comparisons of nonsynonymous and synonymous substitutions indicated a trend toward higher numbers of nonsynonymous substitutions; however, this difference was not significant statistically among S. hispidus alleles. To examine evolution of DQA alleleswithin Muridae, we performed a phylogenetic analysis that included DQA alleles from S. hispidus, Peromyscus leucopus, M. musculus, R. norvegicus, and six Australian species of Rattus. Results depicted monophyly for each genus, and this concordance between species and gene trees represents a lack of evidence for trans-species persistence of alleles among these genera. Received: 22 October 1998 / Revised: 17 March 1999     

Evolution of a new nonclassical MHC class I locus in two Old World primate species

 HLA-G is a nonclassical major histocompatibility complex (MHC) class I molecule that is expressed only in the human placenta, suggesting that it plays an important role at the fetal-maternal interface. In rhesus monkeys, which have similar placentation to humans, the HLA-G orthologue is a pseudogene. However, rhesus monkeys express a novel placental MHC class I molecule, Mamu-AG, which has HLA-G-like characteristics. Phylogenetic analysis of AG alleles in two Old World primate species, the baboon and the rhesus macaque, revealed limited diversity characteristic of a nonclassical MHC class I locus. Gene trees constructed using classical and nonclassical primate MHC class I alleles demonstrated that the AG locus was most closely related to the classical A locus. Interestingly, gene tree analyses suggested that the AG alleles were most closely related to a subset of A alleles which are the products of an ancestral interlocus recombination event between the A and B loci. Calculation of the rates of synonymous and nonsynonymous substitution at the AG locus revealed that positive selection was not acting on the codons encoding the peptide binding region. In exon 4, however, the rate of nonsynonymous substitution was significantly lower than the rate of synonymous substitution, suggesting that negative selection was acting on these codons. Received: 22 April 1998 / Revised: 15 July 1998     

A statistical examination of hypervariability in complementarity-determining regions of immunoglobulins

To determine the relative importance of gene conversion followed by natural selection and of natural selection for point mutation in generating variability in immunoglobulins, the numbers of synonymous and nonsynonymous substitutions in immunoglobulin sequences of various subgroups were estimated for complementarity-determining regions (CDRs) and for framework regions (FRs). Both the number of synonymous substitutions and the number of nonsynonymous substitutions in the CDR were found to exceed the corresponding numbers in the FR. Therefore, gene conversion is likely to be an important mechanism for providing variability in the CDR of immunoglobulins. The correlation coefficients between the number of synonymous substitutions and the number of nonsynonymous substitutions and between the substitution number in the CDR and that in the FR were found to be very low. Again, gene conversion is thought to be responsible for this finding.     

Polymorphism and transcription of Mhc class I genes in a passerine bird, the great reed warbler

 The class I genes of the major histocompatibility complex (Mhc) are here investigated for the first time in a passerine bird. The great reed warbler is a rare species in Sweden with a few semi-isolated populations. Yet, we found extensive Mhc class I variation in the study population. The variable exon 3, corresponding to the α₂ domain, was amplified from genomic DNA with degenerated primers. Seven different genomic class I sequences were detected in a single individual. One of the sequences had a deletion leading to a shift in the reading frame, indicating that it was not a functional gene. A randomly selected clone was used as a probe for restriction fragment length polymorphism (RFLP) studies in combination with the restriction enzyme Pvu II. The RFLP pattern was complex with 21–25 RFLP fragments per individual and extensive variation. Forty-nine RFLP genotypes were detected in 55 tested individuals. To study the number of transcribed genes, we isolated 14 Mhc class I clones from a cDNA library from a single individual. We found eight different sequences of four different lengths (1.3–2.2 kilobases), suggesting there are at least four transcribed loci. The number of nonsynonymous substitutions (d _N ) in the peptide binding region of exon 3 were higher than the number of synonymous substitutions (d _S ), indicating balancing selection in this region. The number of transcribed genes and the numerous RFLP fragments found so far suggest that the great reed warbler does not have a "minimal essential Mhc" as has been suggested for the chicken. Received: 13 May 1998 / Revised: 18 August 1998     

Sequence variation at the major histocompatibility complex DRB loci in beluga (Delphinapterus leucas) and narwhal (Monodon monoceros)

 The variation at loci with similarity to DRB class II major histocompatibility complex loci was assessed in 313 beluga collected from 13 sampling locations across North America, and 11 narwhal collected in the Canadian high Arctic. Variation was assessed by amplification of exon 2, which codes for the peptide binding region, via the polymerase chain reaction, followed by either cloning and DNA sequencing or single-stranded conformation polymorphism analysis. Two DRB loci were identified in beluga: DRB1, a polymorphic locus, and, DRB2, a monomorphic locus. Eight alleles representing five distinct lineages (based on sequence similarity) were found at the beluga DRB1 locus. Although the relative number of alleles is low when compared with terrestrial mammals, the amino acid variation found among the lineages is moderate. At the DRB1 locus, the average number of nonsynonymous substitutions per site is greater than the average number of synonymous substitutions per site (0.0806 : 0.0207, respectively;P<0.01). Most of the 31 amino acid substitutions do not conserve the physiochemical properties of the residue, and 21 of these are located at positions implicated as forming pockets responsible for the selective binding of foreign peptide side chains. Only DRB1 variation was examined in 11 narwhal, revealing a low amount of variation. These data are consistent with an important role for the DRB1 locus in the cellular immune response of beluga. In addition, the ratio of nonsynonymous to synonymous substitutions is similar to that among primate alleles, arguing against a reduction in the balancing selection pressure in the marine environment. Two hypotheses may explain the modest amount of Mhc variation when compared with terrestrial mammals: small population sizes at speciation or a reduced neutral substitution rate in cetaceans. Received: 15 July 1997 / Revised: 24 March 1998     

Diversity and locus specificity of chicken MHC B class I sequences

The major histocompatibility complex B (MHC B) region in a standard haplotype of Leghorn chickens contains two closely linked class I loci, B-FI and B-FIV. Few sequences of B-FI alleles are available, and therefore alleles of the two loci have not been compared with regard to sequence diversity or locus specificity. Here, we report eight new B-F alpha 1/alpha 2-coding sequences from broiler chicken MHC B haplotypes, and a unique recombinant between the two B-F loci. The new sequences were combined with existing B-F sequences from Leghorn and broiler haplotypes for analysis. On the basis of phylogenetic analysis and conserved sequence motifs, B-F sequences separated into two groups (Groups A and B), corresponding to B-FIV and B-FI locus, respectively. Every broiler haplotype had one B-F sequence in Group A and the second B-F sequence, if it existed, clustered in Group B. Group B (presumptive B-FI locus) sequences identified in broiler haplotypes resembled the human MHC class I HLA-C locus in their distinctive pattern of allelic polymorphism. Compared with B-FIV, B-FI alleles were less polymorphic and possessed a conserved locus-specific motif in the alpha1 helix, but nevertheless demonstrated evidence of diversifying selection. One B-FI alpha 1/alpha 2-coding nucleotide sequence was completely conserved in four different broiler haplotypes, but each allele differed in the exon encoding the alpha 3 domain.     

Sequence variability analysis on major histocompatibility complex class II DRB alleles in three felines

The variation of the exon 2 of the major histo-compatibility complex (MHC) class II gene DRB locus in three feline species were examined on clouded leopard (Neofelis nebulosa), leopard (Panthera pardus) and Amur tiger (Panthera tigris altaica). A pair of degenerated primers was used to amplify DRB locus covering almost the whole exon 2. Exon 2 encodes the β1 domain which is the most variable fragments of the MHC class II molecule. Single-strand conformational polymorphism (SSCP) analysis was applied to detect different MHC class II DRB haplotypes. Fifteen recombinant plasmids for each individual were screened out, isolated, purified and sequenced finally. Totally eight distinct haplotypes of exon 2 were obtained in four individuals. Within 237 bp nucleotide sequences from four samples, 30 variable positions were found, and 21 putative peptide-binding positions were disclosed in 79 amino acid residues. The ratio of nonsynonymous substitutions (d _N ) was much higher than that of synonymous substitutions (d _S ), which indicated that balancing selection probably maintain the variation of exon 2. MEGA neighbor joining (NJ) and PAUP maximum parsimony (MP) methods were used to reconstruct phylogenetic trees among species, respectively. Results displayed a more close relationship between leopard and tiger; however, clouded leopard has a comparatively distant relationship form the other two. __________ Translated from Zoological Research, 2006, 27(2): 181-C188 [译自:动物学研究]     

Allelic repertoire of the humanMHC class IMICA gene

 The hallmark of the classical major histocompatibility complex (MHC) class I molecules is their astonishing level of polymorphism, a characteristic not shared by the nonclassical MHC class I genes. A distinct family of MHC class I genes has been recently identified within the human MHC class I region. The MICA (MHC class I chain-related A) gene in this family is a highly divergent member of the MHC class I family and has a unique pattern of tissue expression. We have sequenced exons encoding the extracellular α1, α2, and α3 domains of the MICA gene from twenty HLA homozygous typing cell lines and four unrelated individuals. We report the identification of eleven new alleles defined by a total of twenty-two amino acid substitutions. Thus, the total number of MICA alleles is sixteen. Interestingly, a tentative superimposition of MICA variable residues on the HLA-A2 structure reveals a unique pattern of distribution, concentrated primarily on the outer edge of the MICA putative antigen binding cleft, apparently bordering an invariant ligand binding site. Received: 13 May 1996 / Revised: 29 May 1996     

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.     

Evidence of the domestication history of flax (Linum usitatissimum L.) from genetic diversity of the sad2 locus

A phylogenetic analysis was conducted on 34 alleles of 2.5 kb sized stearoyl-ACP desaturase II (sad2), obtained from 30 accessions of cultivated and pale flax (Linum spp.), to elucidate the history of flax domestication. The analysis supports a single domestication origin for extant cultivated flax. The phylogenetic evidence indicates that flax was first domesticated for oil, rather than fibre. The genetic diversity of the sad2 locus in cultivated flax is low when compared to that of the pale flax assayed. An absolute archaeological date could be applied to the synonymous substitution rate of sad2 in cultivated flax, yielding a high estimate of 1.60–1.71×10⁻⁷ substitutions/site/year. The occurrence of nonsynonymous substitutions at conserved positions of the third exon in alleles from cultivated flax suggests that the locus may have been subjected to an artificial selection pressure. The elevated synonymous substitution rate is also compatible with a population expansion of flax since domestication, followed by a population decline in historic times. These findings provide new insight into flax domestication and are significant for the continuous exploration of the flax germplasm for utilization.     

MHC class I genes of the channel catfish: sequence analysis and expression

 Four cDNAs encoding the major histocompatibility complex (MHC) class I α chain were isolated from a channel catfish clonal B-cell cDNA library. Sequence analysis suggests these cDNAs represent three different MHC class I loci. All cDNAs encoded conserved residues characteristic of the MHC class I α chain: namely, those involved in peptide binding, salt bridges, disulfide bond formation, and glycosylation. Southern blot analyses of individual outbred and second-generation gynogenetic fish indicated the existence of both polygenic and polymorphic loci. Northern blot studies demonstrated that catfish B, T, and macrophage cell lines transcribed markedly higher levels of class I α and β₂-microglobulin (β₂m) mRNA than fibroblast cell lines. In addition, immunoprecipitation data showed that a 41 000 M _r glycoprotein (presumably class I α) was associated with β₂m on the surface of catfish B cells. This latter finding is the first direct evidence for the cell surface association of β₂m with the MHC class I α chain on teleost cells and supports the notion that functional MHC class I proteins exist in teleosts. Received: 25 March 1998 / Revised: 28 July 1998     

中国地方鸡种MHC B-LBⅡ新等位基因的遗传多态性研究

为研究鸡MHC B-LBⅡ基因的遗传多态性,首先在8个中国地方鸡种(藏鸡、仙居鸡、北京油鸡、固始鸡、斗鸡、丝羽乌骨鸡、白耳鸡和狼山鸡)B-LBⅡ基因第二外显子扩增了一长度为 175 bp 的 DNA 片段并进行 SSCP 基因型分析;在8 个地方鸡种共 467 个个体中检测到 37 个 PCR-SSCP 基因型;从被检样品中筛选出不同基因型的个体,并在其 B-LBⅡ基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,通过克隆和测序获得了该片段的核苷酸序列。经序列分析,在前述地方鸡种被筛选出的 30 个无血缘关系的个体中发现了 31 个 B-LBⅡ新等位基因,并参照哺乳动物 MHC II 类 B 等位基因命名规则进行了命名。对这 31 个 B-LBⅡ新等位基因长度为 374 bp 的 DNA 片段进行比对表明,在其第二外显子序列上共有 68 个多态性变异位点,其中简约性信息位点 51 个,单变异位点 17 个,具有丰富的遗传多态性。在这些多态性变异位点中,出现在遗传密码子第一和第二位上的碱基替换率分别为 36.76% 和 35.29%。等位基因序列间的相似性估测为 90.6%-99.5%;B-LBⅡ基因第二外显子的错义替换率和同义替换率分别为 14.64±2.67%和 2.92±0.94%。结果表明,B-LBⅡ基因的丰富遗传多态性主要是由基因重组和平衡选择效应所引起的。对 B-LBⅡ等位基因第二外显子所编码的 B-LBⅡ分子β1 结构域氨基酸序列比对发现,31 个 B-LBⅡ新等位基因属于 26 个等位基因主型;在β1结构域氨基酸序列的 33个变异位点上,存在 6 个同义替换和 27 个错义替换。分析认为,那些发生在多肽结合位点上的氨基酸错义替换与鸡 MHC B-LBⅡ分子的免疫特异性有关。该结果可为鸡的抗病育种研究提供分子生物学依据。     

Characterization of Mhc-DRB allelic diversity in white-tailed deer (Odocoileus virginianus) provides insight into Mhc-DRB allelic evolution within Cervidae

 Although white-tailed deer (Odocoileus virginianus) are one of North America's best studied mammals, no information is available concerning allelic diversity at any locus of the major histocompatibility complex in this taxon. Using the polymerase chain reaction, single-stranded conformation polymorphism analysis, and DNA sequencing techniques, 15 DRB exon 2 alleles were identified among 150 white-tailed deer from a single population in southeastern Oklahoma. These alleles represent a single locus and exhibit a high degree of nucleotide and amino acid polymorphism, with most amino acid variation occurring at positions forming the peptide binding sites. Furthermore, twenty-seven amino acid residues unique to white-tailed deer DRB alleles were detected, with 19 of these occurring at residues forming contact points of the peptide binding region. Significantly higher rates of nonsynonymous than synonymous substitutions were detected among these DRB alleles. In contrast to other studies of Artiodactyla DRB sequences, interallelic recombination does not appear to be playing a significant role in the generation of allelic diversity at this locus in white-tailed deer. To examine evolution of white-tailed deer (Odvi-DRB) alleles within Cervidae, we performed a phylogenetic analysis of all published red deer (Ceel-DRB), roe deer (Caca-DRB), and moose (Alal-DRB) DRB alleles. The phylogenetic tree clearly shows a trans-species persistence of DRB lineages among these taxa. Moreover, this phylogenetic tree provides insight into evolution of DRB allelic lineages within Cervidae and may aid in assignment of red deer DRB alleles to specific loci. Received: 25 June 1998 / Revised: 2 September 1998     

Cytological location and further characterization of the Third chromosome resistance gene in Drosophila melanogaster

Mutations in the Third chromosome resistance (Tcr; 3-39.6) gene confer dominant resistance to α-methyl dopa and suggest the gene is involved in catecholamine metabolism. Evidence for involvement in catecholamine metabolism comes from the three phenotypes associated with the mutant Tcr chromosomes dominant resistance, dominant rescue of partially complementing l(2)amd alleles, and recessive lethal phenotypes. Only dominant resistance to αs-methyl dopa, however, was mapped to the Tcr locus. Both recessive lethality and dominant rescue of l(2)amd alleles have now been mapped to the Tcr gene and, through the isolation of a new deletion in the region, we demonstrate these phenotypes are due to a loss of Tcr function. This deletion places the Tcr gene in the 69B4-5 to 69C8-11 region. Additionally, we have tested and verified three predictions of the biochemical model proposed by Bishop, Sherald, and Wright (1989) for the function of the Tcr protein. This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

The immune genes of the major histocompatibility complex (MHC) are classical examples for high levels of genetic diversity and non-neutral evolution. This is particularly true for the regions containing the antigen-binding sites as, for instance, in the exon 2 of the MHC class II gene DRB. We surveyed, for the first time in the order Chiroptera, the genetic diversity within this exon in the sac-winged bat Saccopteryx bilineata. We detected 11 alleles among 85 bats, of which 79 were sampled in one population. Pairwise comparisons revealed that interallelic sequence differences ranged between 3 and 22%, although nucleotide substitutions were not evenly distributed along the exon sequence. This was most probably the result of intragenic recombination. High levels of sequence divergence and significantly more nonsynonymous than synonymous substitutions (d(N)/d(S)>1) suggest long-term balancing selection. Thus, the data are consistent with the hypothesis that recombination gives rise to new alleles at the DRB locus of the sac-winged bat, and these are maintained in the population through balancing selection. In this respect, the sac-winged bat closely resembles other mammalian species.     

Population study of allelic diversity in the human MHC class I-related MIC-A gene

 The polymorphism of major histocompatibility complex (MHC) class I HLA-A, -B, and -C molecules may have evolved through pathogen-driven selection of alleles with diverse peptide-binding specificities. Two MHC-encoded molecules that are distantly related to class I, MIC-A and MIC-B, do not function in the presentation of pathogen-derived peptides to T cells with αβ T-cell receptors (TCRs), but are broadly recognized by intraepithelial T cells with γδ TCRs. However, both MIC-A and MIC-B are polymorphic, displaying an unusual distribution of a number of variant amino acids in their extracellular α1, α2, and α3 domains. In order to further define the polymorphism of MIC-A, we examined its alleles among 275 individuals with common and rare HLA genotypes. Of 16 previously defined alleles, 12 were confirmed and 5 new alleles were identified. A two-by-two analysis of MIC-A and HLA-B alleles uncovered a number of statistically significant associations. These results confirm and extend previous knowledge on the polymorphism of MIC-A. The strong positive linkage of certain MIC-A and HLA-B alleles may have implications for studies related to MHC-associated diseases and transplantation. Received: 5 August 1998 / Revised: 26 October 1998     

Allelic variants of the human MHC class I chain-related B gene (MICB)

 The human major histocompatibility complex (MHC) is located within a 4 megabase segment on chromosome 6p21.3. Recently, a highly divergent MHC class I chain-related gene family, MIC was identified within the class I region. The MICA and MICB genes in this family have unique patterns of tissue expression. The MICA gene is highly polymorphic, with more than 20 alleles identified to date. To elucidate the extent of MICB allelic variations, we sequenced exons 2 (α1), 3 (α2), 4 (α3), and 5 (transmembrane) as well as introns 2 and 4 of this gene in 46 HLA homozygous B-cell lines. We report the identification of eleven alleles based on seven non-synonymous, two synonymous, and four intronic nucleotide variations. Interestingly, one allele has a nonsense mutation resulting in a premature termination codon in the α2 domain. Thus, MICB appears to have fewer alleles than MICA, not unlike the allelic ratio between the HLA-C and -B loci. A preliminary linkage analysis of the MICB alleles with those of the closely located MICA and HLA-B genes revealed no conspicuous linkage disequilibrium between them, implying the presence of a potential recombination hotspot between the MICB and MICA genes. Received: 16 April 1997 / Revised: 19 May 1997     

MHC class II DRB variability and parasite load in the striped mouse (Rhabdomys pumilio) in the Southern Kalahari

Major histocompatibility complex (MHC) variability is believed to be maintained by pathogen-driven selection, mediated either through heterozygous advantage or frequency-dependent selection. However, empirical support for these hypotheses under natural conditions is rare. In this study, we investigated the genetic constitution of the functionally important MHC class II gene (DRB exon 2) and the parasite load in a population of the striped mouse (Rhabdomys pumilio) in the Southern Kalahari. Fifty-eight individuals were genetically examined and the endoparasite load was quantified by counting fecal helminth eggs by using a modified McMaster technique. Thirty-four animals (58.6%) were infected. We identified 20 different MHC alleles with high levels of sequence divergence between alleles. Particularly, the antigen-binding sites revealed a significant higher rate of nonsynonymous substitutions (d(N)) than synonymous substitutions (d(S)), giving strong evidence of balancing selection. Heterozygosity did influence the infection status (being infected or not) and the individual fecal egg count (FEC) value with significantly higher values observed in homozygous individuals. Furthermore, a positive relationship was found between specific alleles and parasite load. The allele Rhpu-DRB*1 significantly occurred more frequently in infected individuals and in individuals with high FEC values (high parasite load). Individuals with the allele Rhpu-DRB*1 had a 1.5-fold higher chance of being infected than individuals without this allele (odds ratio test, P < 0.05). Contrarily, the allele Rhpu-DRB*8 significantly occurred more frequent in individuals with low FEC values. Our results support the hypotheses that MHC polymorphism in R. pumilio is maintained through pathogen-driven selection acting by both heterozygosity advantage and frequency-dependent selection.     

Gene conversion as a source of nucleotide diversity in Plasmodium falciparum

Examination of polymorphisms in the Plasmodium falciparum gene for falcipain 2 revealed that this gene is one of two paralogs separated by 10.8 kb in chromosome 11. We designate the annotated gene denoted chr11.gen_424 as encoding falcipain 2A and the annotated gene denoted chr11.gen_427 as encoding falcipain 2B. The paralogs are 96% identical at the nucleotide level and 93% identical at the amino acid level. The consensus sequences differ in 31/309 synonymous sites and 45/1140 nonsynonymous sites, including three amino acid replacements (V393I, A400P, and Q414E) that are near the catalytic site and that may affect substrate affinity or specificity. In six reference isolates, among 36 synonymous sites and 46 nonsynonymous sites that are polymorphic in the gene for falcipain 2A, falcipain 2B, or both, significant spatial clustering is observed. All but one of the polymorphisms appear to result from gene conversion between the paralogs. The estimated rate of gene conversion between the paralogs may be as many as 1,400 to 1,700 times greater than the rate of mutation. Owing to gene conversion, one of the falcipain 2A alleles is more similar to the falcipain 2B alleles than it is to other falcipain 2A alleles. Divergence among the synonymous sites suggests that the paralogous genes last shared a common ancestor 15.2 MYA, with a range of 8.8 to 20.6 MYA. During this period, the paralogs have acquired 0.10 synonymous substitutions per synonymous site in the coding region. The 5' and 3' flanking regions differ in 47.7% and 39.8% of the nucleotide sites, respectively. Hence synonymous sites and flanking regions are not conserved in sequence in spite of their high AT content and T skew.