SEARCH · Search PubMed

Results for “synonymous substitution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

PubMed2002

Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes.

The nonsynonymous to synonymous substitution rate ratio (omega = d(N)/d(S)) provides a sensitive measure of selective pressure at the protein level, with omega values <1, =1, and >1 indicating purifying selection, neutral evolution, and diversifying selection, respectively. Maximum likelihood models of codon substitution developed recently account for variable selective pressures among amino acid sites by employing a statistical distribution for the omega ratio among sites. Those models, called random-sites models, are suitable when we do not know a priori which sites are under what kind of selective pressure. Sometimes prior information (such as the tertiary structure of the protein) might be available to partition sites in the protein into different classes, which are expected to be under different selective pressures. It is then sensible to use such information in the model. In this paper, we implement maximum likelihood models for prepartitioned data sets, which account for the heterogeneity among site partitions by using different omega parameters for the partitions. The models, referred to as fixed-sites models, are also useful for combined analysis of multiple genes from the same set of species. We apply the models to data sets of the major histocompatibility complex (MHC) class I alleles from human populations and of the abalone sperm lysin genes. Structural information is used to partition sites in MHC into two classes: those in the antigen recognition site (ARS) and those outside. Positive selection is detected in the ARS by the fixed-sites models. Similarly, sites in lysin are classified into the buried and solvent-exposed classes according to the tertiary structure, and positive selection was detected at the solvent-exposed sites. The random-sites models identified a number of sites under positive selection in each data set, confirming and elaborating the results of the fixed-sites models. The analysis demonstrates the utility of the fixed-sites models, as well as the power of previous random-sites models, which do not use the prior information to partition sites.

Alleles↗

PubMed2003

Molecular epidemiological study of dengue virus type 1 in Taiwan.

Taiwan has experienced several major outbreaks of dengue (DEN) virus since 1981. The predominant virus type involved has been dengue virus type one (DEN-1), which first appeared in 1987. To understand the molecular epidemiology of this virus, 15 strains of DEN-1 isolated during 1987-1991 and 1994-1995, including 11 epidemic strains, two sporadic strains, and two imported strains have been studied. Fragments of 490 nucleotides (nt) from the E/NS1 junction were amplified by reverse transcription-polymerase chain reaction and the nt sequences were determined. Of the 490 nt of the E/NS1 junction, 240 nt (nt 2282-2521) were aligned and compared. Nucleotide substitutions were found at 54 positions among 15 isolates. Most nt changes were synonymous substitutions, and only three amino acid changes were found. A total of 61 strains isolated worldwide were analyzed by the Neighbor-joining method, and separated phylogenetically into three distinct genotypes, I-III. Genotype I comprised isolates from Japan and Hawaii collected in the 1940s. Genotype II included most strains isolated from Asia in 1977-1995. Genotype III consisted of isolates from three continents in 1964-1995: Asia, the Americas, and Africa. Genotype III was divided further into two subgenotypes, IIIA and IIIB. Most recent isolates from Taiwan, except for the sporadic strain isolated in 1995, were similar genetically and have been classified as Genotype II.

Amino Acid Sequence↗

PubMed1995

Comparison of the evolution of ribulose-1, 5-biphosphate carboxylase (rbcL) and atpB-rbcL noncoding spacer sequences in a recent plant group, the tribe Rubieae (Rubiaceae).

Plastid sequences of the atpB-rbcL spacer and rbcL gene itself were used to evaluate their respective potential in reconstructing the phylogeny of 15 taxa from the tribe Rubieae (Rubiaceae). From our previous analyses using the atpB-rbcL spacer, the 15 selected taxa represent most of the variability of the tribe. Since this group is considered to be relatively recent (Upper Tertiary), it should allow the study of early dynamics of nucleotide substitutions in recent divergences. The results show that the spacer and rbcL inferred phylogenies are not totally congruent; the spacer trees are more similar to interpretations of morphological data. A comparative analysis of the pattern of nucleotide substitution of these two sequences in the Rubieae shows that (1) the overall rate of substitution is similar in the spacer and in rbcL, and the rate of synonymous substitution in rbcL is much higher; (2) the level of homoplasy is higher in rbcL than in the spacer matrix which shows a higher phylogenetic structure; and (3) the pattern of transition and transversion substitutions is different in the two sequences, and is not linear in rbcL. As a result of these observations, we suggest that (1) the spacer is evolving relatively slowly because of unsuspected, and phylogenetically important, selective constraints on its sequence; and (2) in the rbcL sequence, many sites, free of constraint, are changing at high rate, and some of these sites seem to have undergone multiple substitutions even in this recent tribe. This could explain the high level of homoplasy found in Rubieae rbcL sequences.

Base Sequence↗

PubMed2000

Neutral evolution of the sex-determining gene transformer in Drosophila.

The amino acid sequence of the transformer (tra) gene exhibits an extremely rapid rate of evolution among Drosophila species, although the gene performs a critical step in sex determination. These changes in amino acid sequence are the result of either natural selection or neutral evolution. To differentiate between selective and neutral causes of this evolutionary change, analyses of both intraspecific and interspecific patterns of molecular evolution of tra gene sequences are presented. Sequences of 31 tra alleles were obtained from Drosophila americana. Many replacement and silent nucleotide variants are present among the alleles; however, the distribution of this sequence variation is consistent with neutral evolution. Sequence evolution was also examined among six species representative of the genus Drosophila. For most lineages and most regions of the gene, both silent and replacement substitutions have accumulated in a constant, clock-like manner. In exon 3 of D. virilis and D. americana we find evidence for an elevated rate of nonsynonymous substitution, but no statistical support for a greater rate of nonsynonymous relative to synonymous substitutions. Both levels of analysis of the tra sequence suggest that, although the gene is evolving at a rapid pace, these changes are neutral in function.

Animals↗

PubMed2003

Neutral substitutions occur at a faster rate in exons than in noncoding DNA in primate genomes.

Point mutation rates in exons (synonymous sites) and noncoding (introns and intergenic) regions are generally assumed to be the same. However, comparative sequence analyses of synonymous substitutions in exons (81 genes) and that of long intergenic fragments (141.3 kbp) of human and chimpanzee genomes reveal a 30%-60% higher mutation rate in exons than in noncoding DNA. We propose a differential CpG content hypothesis to explain this fundamental, and seemingly unintuitive, pattern. We find that the increased exonic rate is the result of the relative overabundance of synonymous sites involved in CpG dinucleotides, as the evolutionary divergence in non-CpG sites is similar in noncoding DNA and synonymous sites of exons. Expectations and predictions of our hypothesis are confirmed in comparisons involving more distantly related species, including human-orangutan, human-baboon, and human-macaque. Our results suggest an underlying mechanism for higher mutation rate in GC-rich genomic regions, predict nonlinear accumulation of mutations in pseudogenes over time, and provide a possible explanation for the observed higher diversity of single nucleotide polymorphisms (SNPs) in the synonymous sites of exons compared to the noncoding regions.

Animals↗

PubMed2005

Multiple major increases and decreases in mitochondrial substitution rates in the plant family Geraniaceae.

BACKGROUND: Rates of synonymous nucleotide substitutions are, in general, exceptionally low in plant mitochondrial genomes, several times lower than in chloroplast genomes, 10-20 times lower than in plant nuclear genomes, and 50-100 times lower than in many animal mitochondrial genomes. Several cases of moderate variation in mitochondrial substitution rates have been reported in plants, but these mostly involve correlated changes in chloroplast and/or nuclear substitution rates and are therefore thought to reflect whole-organism forces rather than ones impinging directly on the mitochondrial mutation rate. Only a single case of extensive, mitochondrial-specific rate changes has been described, in the angiosperm genus Plantago. RESULTS: We explored a second potential case of highly accelerated mitochondrial sequence evolution in plants. This case was first suggested by relatively poor hybridization of mitochondrial gene probes to DNA of Pelargonium hortorum (the common geranium). We found that all eight mitochondrial genes sequenced from P. hortorum are exceptionally divergent, whereas chloroplast and nuclear divergence is unexceptional in P. hortorum. Two mitochondrial genes were sequenced from a broad range of taxa of variable relatedness to P. hortorum, and absolute rates of mitochondrial synonymous substitutions were calculated on each branch of a phylogenetic tree of these taxa. We infer one major, approximately 10-fold increase in the mitochondrial synonymous substitution rate at the base of the Pelargonium family Geraniaceae, and a subsequent approximately 10-fold rate increase early in the evolution of Pelargonium. We also infer several moderate to major rate decreases following these initial rate increases, such that the mitochondrial substitution rate has returned to normally low levels in many members of the Geraniaceae. Finally, we find unusually little RNA editing of Geraniaceae mitochondrial genes, suggesting high levels of retroprocessing in their history. CONCLUSION: The existence of major, mitochondrial-specific changes in rates of synonymous substitutions in the Geraniaceae implies major and reversible underlying changes in the mitochondrial mutation rate in this family. Together with the recent report of a similar pattern of rate heterogeneity in Plantago, these findings indicate that the mitochondrial mutation rate is a more plastic character in plants than previously realized. Many molecular factors could be responsible for these dramatic changes in the mitochondrial mutation rate, including nuclear gene mutations affecting the fidelity and efficacy of mitochondrial DNA replication and/or repair and--consistent with the lack of RNA editing--exceptionally high levels of "mutagenic" retroprocessing. That the mitochondrial mutation rate has returned to normally low levels in many Geraniaceae raises the possibility that, akin to the ephemerality of mutator strains in bacteria, selection favors a low mutation rate in plant mitochondria.

Base Sequence↗

PubMed1992

Evolution of Ac and Dsl elements in select grasses (Poaceae).

We present data on evolution of the Ac/Ds family of transposable elements in select grasses (Poaceae). An Ac-like element was cloned from a DNA library of the grass Pennisetum glaucum (pearl millet) and 2387 bp of it have been sequenced. When the pearl millet Ac-like sequence is aligned with the corresponding region of the maize Ac sequence, it is found that all sequences corresponding to intron II in maize Ac are absent in pearl millet Ac. Kimura's evolutionary distance between maize and pearl millet Ac sequences is estimated to be 0.429 +/- 0.020 nucleotide substitutions per site. This value is not significantly different from the average number of synonymous substitutions for coding regions of the Adh1 gene between maize and pearl millet, which is 0.395 +/- 0.051 nucleotide substitutions per site. If we can assume Ac and Adh1 divergence times are equivalent between maize and pearl millet, then the above calculations suggest Ac-like sequences have probably not been strongly constrained by natural selection. The level of DNA sequence divergence between maize and pearl millet Ac sequences, the estimated date when maize and pearl millet diverged (25-40 million years ago), coupled with their reproductive isolation/lack of current genetic exchange, all support the theory that Ac-like sequences have not been recently introduced into pearl millet from maize. Instead, Ac-like sequences were probably present in the progenitor of maize and pearl millet, and have thus existed in the grasses for at least 25 million years. Ac-like sequences may be widely distributed among the grasses. We also present the first 2 Ds1 controlling element sequences from teosinte species: Zea luxurians and Zea perennis. A total of 10 Ds1 elements had previously been sequenced from maize and a distant maize relative, Tripsacum. When a maximum likelihood network of genetic relationships is constructed for all 12 sequenced Ds1 elements, the 2 teosinte Ds1 elements are as distant from most maize Ds1 elements and from each other, as the maize Ds1 elements are from one another. Our new teosinte sequence data support the previous conclusion that Ds1 elements have been accumulating mutations independently since maize and Tripsacum diverged. We present a scenario for the origin of Ds1 elements.

Amino Acid Sequence↗

PubMed1994

Molecular evolutionary characterization of an Activator (Ac)-like transposable element sequence from pearl millet (Pennisetum glaucum) (Poaceae).

We present data on the evolution of the Ac/Ds family of transposable elements in select grasses (Poaceae). A defective Ac-like element was cloned from a DNA library of the grass Pennisetum glaucum (pearl millet) and its entire 4531 bp sequence has been determined. When the pearl millet Ac-like sequence is aligned with the maize Ac sequence, it is found that there is approximately 70% DNA similarity in the central region spanning most of maize Ac exon II and all of exon III. In addition, there are two smaller regions of similarity at the Ac terminii. Besides these three major structural similarities, Pennisetum Ac has two large regions, one 5' and one 3', that show little similarity to Zea Ac. Furthermore, most of the sequences corresponding to intron II in maize Ac are absent in pearl millet Ac. Kimura's evolutionary distance between the central region of maize and pearl millet Ac sequences is estimated to be 0.429 +/- 0.020 nucleotide substitutions per site. This value is not significantly different from the average number of synonymous substitutions for coding regions of the Adh1 gene between maize and pearl millet, which is 0.395 +/- 0.051 nucleotide substitutions per site. If we assume Ac and Adh1 divergence times are equivalent between maize and pearl millet, then the above calculations suggest Ac-like sequences have probably not been strongly constrained by natural selection. Conserved DNA and amino acid sequence motifs are also examined. The level of DNA sequence divergence between maize and pearl millet Ac sequences, the estimated date when maize and pearl millet diverged (25-40 million years ago), coupled with their reproductive isolation/lack of current genetic exchange, all support the theory that Ac-like sequences have not been recently introduced into pearl millet from maize. Instead, Ac-like sequences were probably present in the progenitor of maize and pearl millet and have thus existed in the grasses for at least 25 million years.

Amino Acid Sequence↗

PubMed2005

Substitution rate and structural divergence of 5'UTR evolution: comparative analysis between human and cynomolgus monkey cDNAs.

The substitution rate and structural divergence in the 5'-untranslated region (UTR) were investigated by using human and cynomolgus monkey cDNA sequences. Due to the weaker functional constraint in the UTR than in the coding sequence, the divergence between humans and macaques would provide a good estimate of the nucleotide substitution rate and structural divergence in the 5'UTR. We found that the substitution rate in the 5'UTR (K5UTR) averaged approximately 10%-20% lower than the synonymous substitution rate (Ks). However, both the K5UTR and nonsynonymous substitution rate (Ka) were significantly higher in the testicular cDNAs than in the brain cDNAs, whereas the Ks did not differ. Further, an in silico analysis revealed that 27% (169/622) of macaque testicular cDNAs had an altered exon-intron structure in the 5'UTR compared with the human cDNAs. The fraction of cDNAs with an exon alteration was significantly higher in the testicular cDNAs than in the brain cDNAs. We confirmed by using reverse transcriptase-polymerase chain reaction that about one-third (6/16) of in silico "macaque-specific" exons in the 5'UTR were actually macaque specific in the testis. The results imply that positive selection increased K5UTR and structural alteration rate of a certain fraction of genes as well as Ka. We found that both positive and negative selection can act on the 5'UTR sequences.

Animals↗

PubMed2001

MtDNA mutations in maternally inherited diabetes: presence of the 3397 ND1 mutation previously associated with Alzheimer's and Parkinson's disease.

Mutations in the mitochondrial tRNA(leu) (UUR) gene have been associated with diabetes mellitus and deafness. We screened for the presence of mtDNA mutations in the tRNA(leu) (UUR) gene and adjacent ND1 sequences in 12 diabetes mellitus pedigrees with a possible maternal inheritance of the disease. One patient carried a G to A substitution at nt 3243 (tRNA(leu) (UUR) gene) in heteroplasmic state. In a second pedigree a patient had an A to G substitution at nt 3397 in the ND1 gene. All maternal relatives of the proband had the 3397 substitution in homoplasmic state. This substitution was not present in 246 nonsymptomatic Caucasian controls. The 3397 substitution changes a highly conserved methionine to a valine at aa 31 and has previously been found in Alzheimer's (AD) and Parkinson's (PD) disease patients. Substitutions in the mitochondrial ND1 gene at aa 30 and 31 have associated with a number of different diseases (e.g. AD/PD, MELAS, cardiomyopathy and diabetes mellitus, LHON, Wolfram-syndrome and maternal inherited diabetes) suggesting that changes at these two codons may be associated with very diverse pathogenic processes. In a further attempt to search for mtDNA mutations outside the tRNAleu gene associated with diabetes, the whole mtDNA genome sequence was determined for two patients with maternally inherited diabetes and deafness. Except for substitutions previously reported as polymorphisms, none of the two patients showed any non-synonymous substitutions either in homoplasmic or heteroplasmic state. These results imply that the maternal inherited diabetes and deafness in these patients must result from alterations of nuclear genes and/or environmental factors.

Alzheimer Disease↗

PubMed2000

Sequence variations of the hypervariable region of hepatitis C virus and their clinical significance.

OBJECTIVE: To understand the clinical significance of sequence variations in the hypervariable region (HVR) of hepatitis C virus during infection. METHODS: Eight patients with acute hepatitis C and 20 patients with chronic hepatitis C were followed up for two years. Blood samples were taken at intervals of six months for analysis of HCV-HVR sequences by reverse transcription-polymerase chain reaction (RT-PCR) and direct sequencing methods. RESULTS: HCV-HVR sequences of the 28 patients changed in various degrees. 92% of these nucleotide substitutions led to changes of corresponding aminoacid sequences. Only 8% of changed nucleotide were synonymous substitutions. Of 27 amino acids variation of amino acid ranged from 1 to 20 (mean 8, 30%). The most common nucleotide substitution (62%) occurred in the first position of codon, 31% in the second and the rest in the third. HVR variation rate was 0.89 x 10(-1) per genome site per year in acute hepatitis C, compared with 2.31 x 10(-1) per genome site per year in chronic hepatitis C (P < 0.05), but had no relation to HCV subtype. Variation of HVR in the flare up type (ALT > 150 mu/L) was much more than that in the quiescent type (ALT < 100 mu/L). CONCLUSION: Our results suggested that sequence variation of HVR during HCV chronic infection seems to be an adaptive response to HCV to evade the host immune pressure and might play a major role in the establishment of persistent infection as well as in the flare-up of hepatitis.

Acute Disease↗

PubMed1994

Evolution of pathogenic viruses with special reference to the rates of synonymous and nonsynonymous substitutions.

For pathogenic viruses such as human immunodeficiency virus (HIV), hepatitis B virus (HBV), hepatitis C virus (HCV), human influenza A virus, and human T-cell leukemia virus type I (HTLV-I), the evolutionary features were briefly reviewed with special reference to the rates of synonymous and nonsynonymous substitutions. In particular, these rates were discussed in connection with the neutral theory of molecular evolution. It was common to all the five pathogenic viruses that the rate of synonymous substitution was higher than that of nonsynonymous substitution particularly when the entire gene regions were compared between different isolates. This suggests that the viral proteins are quite conservative to functional and structural changes even though most of these viral genomes are evolving at a speed extraordinarily higher than their host genomes. Thus, this feature is consistent with the neutral theory. However, it is also pointed out that positive selection may be operating on some specific sites such as antigenic sites in order for the pathogenic viruses to escape from the host immune system.

Base Sequence↗

PubMed1989

Human glyceraldehyde-3-phosphate dehydrogenase pseudogenes: molecular evolution and a possible mechanism for amplification.

We screened two human genomic libraries and isolated 14 different clones, designated lambda G1 and EG1-EG13, homologous to human glyceraldehyde-3-phosphate dehydrogenase (GAPD) cDNA. Subcloning and sequencing these recombinant phages led us to classify them as five different pseudogenes (psi G1-psi G5). All these sequences show such features typical of processed pseudogenes as numerous mutations, insertions, and deletions. The identity of numerous mutated sites among these pseudogenes and the presence of two Alu sequences flanking both ends of psi G1 suggest that GAPD pseudogenes originated from a unique reverse transcribed mRNA followed by gene duplication. The rate of nucleotide substitutions per site per year for known GAPD functional genes is low both for the synonymous substitutions (1.87 x 10(-9] and for the nonsynonymous substitutions (0.12 x 10(-9] and indicates that the GAPD cDNA sequence is well conserved not only at the amino acid level, but also at the nucleotide level. The rate of nucleotide substitutions per site per year for GAPD pseudogenes shows a higher value (5.9 x 10(-9] and suggests that these pseudogenes do not have any functional role.

Base Sequence↗

PubMed2000

Identification of 142 single nucleotide polymorphisms in 41 candidate genes for rheumatoid arthritis in the Japanese population.

Single-nucleotide polymorphisms (SNPs) can make an important contribution to our understanding of genetic backgrounds that may influence medical conditions and ethnic diversity. We undertook a systematic survey of genomic DNA for SNPs located not only in coding sequences but also in non-coding regions (e.g., introns and 5' flanking regions) of selected genes. Using DNA samples from 48 Japanese patients with rheumatoid arthritis (RA) as templates, we surveyed 41 genes that represent candidates for RA, screening a total of 104 kb of DNA (30 kb of coding sequences and 74 kb of non-coding DNA). Within this 104 kb of genomic sequences we identified 163 polymorphisms (1 per 638 bases on average), of which 142 were single-nucleotide substitutions and the remainder, insertions or deletions. Of the coding SNPs, 52% were non-synonymous substitutions, and non-conservative amino acid changes were observed in a quarter of those. Sixty-nine polymorphisms showed high frequencies for minor alleles (more than 15%) and 20 revealed low frequencies (<5%). Our results indicated a greater average distance between SNPs than others have reported, but this disparity may reflect the type of genes surveyed and/or the relative ethnic homogeneity of our test population.

Alleles↗

PubMed1995

Variation of nitrate reductase genes in selected grass species.

In order to study the variation of nitrate reductase (NR) genes among grass species, gene number, intron size and number, and the heme-hinge fragment sequence of 25 grass species were compared. Genomic DNA cut with six restriction enzymes and hybridized with the barley NAD(P)H and NADH NR gene probes revealed a single NAD(P)H NR gene copy and two or more NADH NR gene copies per haploid genome in most of the species examined. Major exceptions were Hordeum vulgare, H. vulgare ssp. spontaneum, and Avena strigosa, which appeared to have a single NADH NR gene copy. The NADH NR gene intron number and lengths were examined by polymerase chain reaction amplification. Introns I and III appeared to be absent in at least one of the NADH NR genes in the grass species, while intron II varied from 0.8 to 2.4 kilobases in length. The NADH NR gene heme-hinge regions were amplified and sequenced. The estimated average overall nucleotide substitution rate in the sequenced region was 7.8 x 10(-10) substitutions/site per year. The synonymous substitution rate was 2.11 x 10(-9) substitutions/synonymous site per year and the nonsynonymous substitution rate was 4.10 x 10(-10) substitutions/nonsynonymous site per year. Phylogenetic analyses showed that all of the wild Hordeum species examined clustered in a group separate from H. vulgare and H. vulgare ssp. spontaneum.

Base Sequence↗

PubMed1978

The proportion of all point mutations which are unacceptable: an estimate based on hemoglobin amino acid and nucleotide sequences.

Statistical analysis of the distribution of 156 kinds of human hemoglobin beta (Hbbeta) chain variants suggests that mutations are essentially random in their location. Thus differential fitness, not differential mutability, is the principal source of nonrandom distribution of interspecies differences in Hbbeta amino acid sequence. Similar analyses of both the location and the kind of interspecies differences detected among primates support this viewpoint and lead us to estimate that at least 95% of all amino acid subsitutions,i.e., nonsynonymous mutations, in Hbbeta are functionally unacceptable in homozygous state. Through the combined use of this estimate and the number of nonsynonymous and synonymous substitutions per nucleotide site inferred from comparisons of entire human and rabbit HbbetamRNA nucleotide sequences, we calculate (a) approximately 70% of synonymous Hbbeta mutations are adaptively undersirable and (b) the mutation rate underlying all changes is lesser than or equal to 10(-8) nucleotide substitutions per nucleotide site per year. Apart from such calculations, analyses of nucleotide patterns in HbbetamRNA as well as in rat preproinsulin mRNA reinforce the notion that a large portion of synonymous mutations are functionally unacceptable and rendered so by selective constraint, at a pretranslational level, of the abundance of particular nucleotide doublets such as CpG.

Amino Acid Sequence↗

PubMed1998

Synonymous and nonsynonymous substitutions in mammalian genes: intragenic correlations.

Previous investigations indicated that synonymous and nonsynonymous substitution rates are correlated in mammalian genes. In the present work, this correlation has been studied at the intragenic level using a dataset of 48 orthologous genes from species belonging to at least four different mammalian orders. The results obtained show that the intragenic variability in synonymous rates is correlated with that of nonsynonymous rates. Moreover, the variation in GC level (and especially of C level) of silent positions along each gene is correlated with the variation in synonymous rate. These results reinforce the previous conclusions that synonymous and nonsynonymous rates as well as GC levels of silent positions are to some extent under common selective constraints.

Animals↗

PubMed1999

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 alpha2 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 (dN) in the peptide binding region of exon 3 were higher than the number of synonymous substitutions (dS), 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.

3' Untranslated Regions↗