How the ratio of nonsynonymous to synonymous pseudogene substitutions can be less than one

Cheetah paradigm revisited: MHC diversity in the world's largest free-ranging population

For more than two decades, the cheetah (Acinonyx jubatus) has been considered a paradigm of disease vulnerability associated with low genetic diversity, particularly at the immune genes of the major histocompatibility complex (MHC). Cheetahs have been used as a classic example in numerous conservation genetics textbooks as well as in many related scientific publications. However, earlier studies used methods with low resolution to quantify MHC diversity and/or small sample sizes. Furthermore, high disease susceptibility was reported only for captive cheetahs, whereas free-ranging cheetahs show no signs of infectious diseases and a good general health status. We examined whether the diversity at MHC class I and class II-DRB loci in 149 Namibian cheetahs was higher than previously reported using single-strand conformation polymorphism analysis, cloning, and sequencing. MHC genes were examined at the genomic and transcriptomic levels. We detected ten MHC class I and four class II-DRB alleles, of which nine MHC class I and all class II-DRB alleles were expressed. Phylogenetic analyses and individual genotypes suggested that the alleles belong to four MHC class I and three class II-DRB putative loci. Evidence of positive selection was detected in both MHC loci. Our study indicated that the low number of MHC class I alleles previously observed in cheetahs was due to a smaller sample size examined. On the other hand, the low number of MHC class II-DRB alleles previously observed in cheetahs was further confirmed. Compared with other mammalian species including felids, cheetahs showed low levels of MHC diversity, but this does not seem to influence the immunocompetence of free-ranging cheetahs in Namibia and contradicts the previous conclusion that the cheetah is a paradigm species of disease vulnerability.

More than 80R2R3-MYB regulatory genes in the genome of Arabidopsis thaliana

Transcription factors belonging to the R2R3-MYB family contain the related helix-turn-helix repeats R2 and R3. The authors isolated partial cDNA and/or genomic clones of 78 R2R3-MYB genes from Arabidopsis thaliana and found accessions corresponding to 31 Arabidopsis genes of this class in databanks, seven of which were not represented in the authors' collection. Therefore, there are at least 85, and probably more than 100, R2R3-MYB genes present in the Arabidopsis thaliana genome, representing the largest regulatory gene family currently known in plants. In contrast, no more than three R2R3-MYB genes have been reported in any organism from other phyla. DNA-binding studies showed that there are differences but also frequent overlaps in binding specificity among plant R2R3-MYB proteins, in line with the distinct but often related functions that are beginning to be recognized for these proteins. This large-sized gene family may contribute to the regulatory flexibility underlying the developmental and metabolic plasticity displayed by plants.

Sequence analysis of a polymorphic Mhc class II gene in Pacific salmon

Polymorphism of the nucleotide sequences encoding 149 amino acids of linked major histocompatibility complex (Mhc) class II B1 and B2 peptides, and of the intervening intron (548-773 base pairs), was examined within and among seven Pacific salmon (Oncorhynchus) species. Levels of nucleotide diversity were higher for the B1 sequence than for B2 or the intron in comparisons both within and between species. For the codons of the peptide binding region of the B1 sequence, the level of nonsynonymous nucleotide substitution (dN) exceeded the level of synonymous substitution (dS) by a factor of ten for within-species comparisons, and by a factor of four for between-species comparisons. The excess of dN indicates that balancing selection maintains diversity at this salmonid Mhc class II locus, as is common for Mhc loci in other vertebrates. Levels of nucleotide diversity for both the exon and intron sequences were greater among than within species, and there were numerous species-specific nucleotides present in both the coding and noncoding regions. Thus, neighbor-joining analysis of both the intron and exon regions provided phylogenies in which the sequences clustered strongly by species. There was little evidence of shared ancestral (trans-species) polymorphism in the exon phylogeny, and the intron phylogeny depicted standard relationships among the Pacific salmon species. The lack of shared allelic B1 lineages in these closely related species may result from severe bottlenecks that occurred during speciation or during the ice ages that glaciated the rim of the north Pacific Ocean approximately every 100000 years in the Pleistocene.

Evolutionary relationship of HLA-DRB genes inferred from intron sequences

The major histocompatibility complex (Mhc) consists of class I and class II genes. In the human Mhc (HLA) class II genes, nine DRB loci have been identified. To elucidate the origin of these duplicated loci and allelic divergences at the most polymorphic DRB1 locus, introns 4 and 5 as well as the 3' untranslated region (altogether approximately 1,000 base pairs) of seven HLA-DRB loci, three HLA-DRB1 alleles, and nine nonhuman primate DRB genes were examined. It is shown that there were two major diversification events in HLA-DRB genes, each involving gene duplications and allelic divergences. Approximately 50 million years (my) ago, DRB1*04 and an ancestor of the DRB1*03 cluster (DRB1*03, DRB1*15, and DRB3) diverged from each other and DRB5, DRB7, DRB8, and an ancestor of the DRB2 cluster (DRB2, DRB4, and DRB6) arose by gene duplication. Later, about 25 my ago, DRB1*15 diverged from DRB1*03, and DRB3 was duplicated from DRB1*03. Then, some 20 my ago, the lineage leading to the DRB2 cluster produced two new loci, DRB4 and DRB6. The DRB1*03 and DRB1*04 allelic lineages are extraordinarily old and have persisted longer than some duplicated genes. The orthologous relationships of DRB genes between human and Old World monkeys are apparent, but those between Catarrhini and New World monkeys are equivocal because of a rather rapid expansion and contraction of primate DRB genes by duplication and deletion.

Evolution of Mhc class II polymorphism: the rise and fall of class II gene function in primates

The substitution rate at the codons implicated at ARS of Mhc class II genes has previously been shown to be heavily biased towards nonsynonymous substitutions, indicative of positive selection for polymorphism. Based on our analysis of the number of synonymous changes at codons outside putative ARS in primates, the average age of the polymorphism at class II loci was found to increase in the following order: DPB1, DRB3, DRB5, DRB1, DRB4, DQB1, DQA1. For DRB loci, nonsynonymous changes were found to exceed synonymous changes at HLA-DRB1, DRB3 and DRB5, while no evidence of deviations from equal rates of synonymous and nonsynonymous substitutions were found for DRB6. The pattern of substitutions at the DRB loci of most Catarrhini species indicates constant positive selection at ARS codons over the evolutionary period examined. An exception to the relatively stable selection pattern between species exhibited by most loci is the appearance of polymorphism under positive selection at DRB4 only in the regular chimpanzee. The ds/dn ratios for DQA1 and DQB1 alleles are lower than for the most polymorphic DRB genes. Since the dn/ds ratio of ARS codons may be positively correlated to the ds for non-ARS codons, at least for DQB1, caution must be exercised in interpreting the low ratio for the DQ genes as an indication of weaker selection. The DQA1 allelic lineages show different dn/ds ratios, consistent with the hypothesis that the lineages are constrained from evolving in relation to the diversity of the interacting DQB1 alleles. In contrast to all other class II loci, DPB1 appears to have been subjected to strong positive selection only in the human lineage, and may represent the most conspicuous example of an Mhc locus acquiring an altered function in antigen presentation.

Intensity of natural selection at the major histocompatibility complex loci

Long persistence of allelic lineages, prevalence of nonsynonymous over synonymous substitutions in the peptide-binding region (PBR), and deviation from neutrality of the expected gene identity parameter F all indicate indirectly that balancing selection is operating at functional major histocompatibility complex (MHC) loci. Direct demonstrations of the existence of balancing selection at MHC loci are, however, either lacking or not fully convincing. To define the conditions under which balancing selection could be demonstrated, we estimated its intensity from the mean number of nonsynonymous substitutions, KB, at the PBR and the mutation rate mu. We compared the five available methods for estimating KB by computer simulation and chose the most reliable ones for estimation of selection intensity. For the human MHC, the selection coefficients of the HLA-A, -B, -C, -DRB1, -DQB1, -DQA1, and -DPB1 loci are 0.015, 0.042, 0.0026, 0.019, 0.0085, 0.0028, and 0.0007, respectively. This low selection intensity places severe restrictions on the possibility of measuring selection directly in vertebrate populations.

Structure, function, and evolution of mouse TL genes, nonclassical class I genes of the major histocompatibility complex

In contrast to well-studied "classical" class I genes of the major histocompatibility complex (MHC), the biology of nonclassical class I genes remains largely unexamined. The mouse TL genes constitute one of the best defined systems among nonclassical class I genes in the T region of the MHC. To elucidate the function and the evolution of TL genes and their relationship to classical class I genes, seven TL DNA sequences, including one from a Japanese wild mouse, were examined and compared with those of several mouse and human classical class I genes. The TL genes differ from either classical class I genes or pseudogenes in the extent and pattern of nucleotide substitutions. Natural selection appears to have operated so as to preserve the function of TL, which might have been acquired in an early stage of its evolution. In a putative peptide-binding region encoded by TL genes, the rate of nonsynonymous (amino acid replacing) substitution is considerably lower than that of synonymous substitution. This conservation is completely opposite that in classical class I genes, in which the peptide-binding region has evolved to diversify amino acid sequences so as to recognize a variety of antigens. Thus, it is suggested that the function of TL antigens is distinct from that of classical class I antigens and is related to the recognition of a relatively restricted repertoire of antigens and their presentation to T-cell receptors.