Molecular evolution of the antiretroviral <Emphasis Type="Italic">TRIM5</Emphasis> gene

In 2004, the first report of TRIM5α as a cellular antiretroviral factor triggered intense interest among virologists, particularly because some primate orthologs of TRIM5α have activity against HIV. Since that time, a complex and eventful evolutionary history of the TRIM5 locus has emerged. A review of the TRIM5 literature constitutes a veritable compendium of evolutionary phenomena, including elevated rates of nonsynonymous substitution, divergence in subdomains due to short insertions and deletions, expansions and contractions in gene copy number, pseudogenization, balanced polymorphism, trans-species polymorphism, convergent evolution, and the acquisition of new domains by exon capture. Unlike most genes, whose history is dominated by long periods of purifying selection interspersed with rare instances of genetic innovation, analysis of restriction factor loci is likely to be complicated by the unpredictable and more-or-less constant influence of positive selection. In this regard, the molecular evolution and population genetics of restriction factor loci most closely resemble patterns that have been documented for immunity genes, such as class I and II MHC genes, whose products interact directly with microbial targets. While the antiretroviral activity encoded by TRIM5 provides plausible mechanistic hypotheses for these unusual evolutionary observations, evolutionary analyses have reciprocated by providing significant insights into the structure and function of the TRIM5α protein. Many of the lessons learned from TRIM5 should be applicable to the study of other restriction factor loci, and molecular evolutionary analysis could facilitate the discovery of new antiviral factors, particularly among the many TRIM genes whose functions remain as yet unidentified.     

Strepsirhine dichotomy from a human perspective (Primates: Lorisoidea, Lemuroidea)

Using monospecific, polyclonal antisera against 69 human plasma proteins, 128 antigenic determinants from 40 cross-reacting homologues were characterized in representatives of the prosimian genera Lemur, Eulemur, Varecia (Lemuridae) and Otolemur (Galagidae). Seventeen determinants from 16 different proteins were absent in homologues of the gagalo but were shared by lemurs and several platyrrhines, cercopithecids, and hominoids. Smaller locus samples for potto(Perodicticus potto) and slow loris (Nycticebus coucang) confirmed the more distant immunological relationship of lorises than of lemurs to anthropoids. If evolutionary rates are constant and equal in lorisiform and lemuriform prosimians, this patern of character distribution indicates strepsirhine paraphyly, lorises diverging earlier (possibly some 6 times 10⁶ years) than lemurs from anthropoid ancestors. If this is so, lorises rather than true lemurs should be elected to root the polarity of character evolution in Primates. As an alternative, galagine proteins evolve more rapidly than lemurine homologues.     

Large-scale comparative mapping of the MHC class I region of predominant haplotypes in Japanese

 In order to further our understanding of major histocompatibility complex (MHC) class I gene organization, we began a comparative analysis of the large scale organization of the class I region in diverse haplotypes. For these studies, the MHC in healthy Japanese donors who have the predominant MHC haplotypes and/or HLA-A or -B alleles was examined by pulsed field gel electrophoresis and Southern analysis using probes spanning the class I region. Hybridization with probes from the HLA-A to HLA-G region revealed that individuals expressing HLA-A30, -A31, or -A33 have an approximately 70 kilobase (kb) insertion near the HLA-A gene as compared with haplotypes containing the HLA-A11 or -A26 allele. Conversely, HLA-A24-containing haplotypes appear to have an approximately 50 kb deletion from the same region. Further, it appears that chromosomes carrying closely related alleles are similar to each other in this region, consistent with their presumed evolutionary relationship. While little is known about the gene content between the HLA-A and HLA-G region, it will be interesting to examine the prospect that functional genes do in fact reside within the inserted or deleted portions, thereby raising the possibility that distinct functional differences are conferred by different haplotypes. Overall, the results reported here should contribute to furthering our understanding of the association between diseases and HLA as well as provide new insights into the evolution of the MHC. Received: 11 December 1996     

Allelic and interlocus comparison of the PERB11 multigene family in the MHC

 The major histocompatibility complex (MHC) contains at least a hundred genes over 4 megabases of DNA. Within the MHC there are several new multigene families which have been recently described. PERB11 is a multigene family which occurs over the class I and central region of the MHC. Two members of the family have been shown to be functional and share domains with members of the supergene family including HLA class I, FcRn, and Zn-α2-glycoprotein molecules. The two functional members are contained within an area of the MHC which has been associated with increased susceptibility to autoimmune diseases such as insulin-dependent diabetes mellitus and also rapid progression to AIDS following HIV-1 infection. Intralocus and interlocus differences between PERB11.1 and PERB11.2 include: (1) several nucleotide substitutions leading to amino acid changes; (2) presence and absence of potential glycosylation sites; (3) insertions and deletions leading to a frame shift resulting in diversity at the amino acid level and an early termination signal. There are ten different alleles of PERB11.1 including one allele which contains a frame shift in the transmembrane region causing a putative truncated molecule lacking the cytoplasmic tail. The significance of this polymorphism in disease associations is under investigation. The most divergent domain is the transmembrane region when PERB11.1 and PERB11.2 are compared. The results suggest that these two molecules may have different functions. Received: 23 July 1996 / Revised: 17 September 1996     

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     

Signatures of historical demography and pathogen richness on <Emphasis Type="Italic">MHC</Emphasis> class I genes

The extreme polymorphism of MHC class I has been argued to be driven by balancing selection from pathogens, with the prediction that populations exposed to a wider variety of diseases should have higher diversity. We assembled a global database of allotype frequencies for MHC class I genes and investigated possible drivers of genetic diversity, measured in different ways. We first looked for a decline in diversity with distance from Africa (a consequence of drift during human expansions) and then investigated the link with pathogen richness once the effect of drift had been corrected for. Using heterozygosity, we recovered a clear decline in diversity from Africa and confirmed the positive relationship between genetic diversity and pathogen richness for all three classical MHC class I genes. However, when we considered a sequence-based measure of genetic diversity, the correlation with geographic distance from Africa vanished for HLA-C, and the correlations with pathogen richness for the three MHC class I genes were much weaker. HLA-C is known to consist of two functional classes of allotypes (classified with respect to the 80th residue), which interact with different KIR receptors. While this separation provided some improvement in the fit between genetic diversity and distance from Africa for one class, much clearer and consistent patterns were recovered when we used the 90th residue to separate HLA-C allotypes into two new classes. This suggests that this residue, which is also involved in the binding of KIR, might have had an important evolutionary role that has been overlooked.     

Locus-specific conservation of the HLA class I introns by intra-locus homogenization

 The protein-coding sequences of the major histocompatibility complex (MHC) genes are characterized by extraordinarily high polymorphism, apparently maintained by balancing selection, which favors diversity in the peptide-binding domains of the MHC glycoproteins. Here we report that the introns flanking the polymorphic exons of the human MHC class I loci HLA-A, -B, and -C genes have been relatively conserved and have become locus-specific apparently as a result of recombination and subsequent genetic drift, leading to homogenization within loci over evolutionary time. Thus, HLA class I genes have been shaped by contrasting evolutionary forces maintaining polymorphism in the exons and leading to conservation in the introns. This study provides the first extensive analysis of the introns of a highly polymorphic gene family. Received: 10 April 1997 / Revised: 15 July 1997     

Polymorphism of the HLA class II loci in Siberian populations

The populations that colonized Siberia diverged from one another in the Paleolithic and evolved in isolation until today. These populations are therefore a rich source of information about the conditions under which the initial divergence of modern humans occurred. In the present study we used the HLA system, first, to investigate the evolution of the human major histocompatibility complex (MHC) itself, and second, to reveal the relationships among Siberian populations. We determined allelic frequencies at five HLA class II loci (DRB1, DQA1, DQB1, DPA1, and DPB1) in seven Siberian populations (Ket, Evenk, Koryak, Chukchi, Nivkh, Udege, and Siberian Eskimo) by the combination of single-stranded conformational polymorphism and DNA sequencing analysis. We then used the gene frequency data to deduce the HLA class II haplotypes and their frequencies. Despite high polymorphism at four of the five loci, no new alleles could be detected. This finding is consistent with a conserved evolution of human class II MHC genes. We found a high number of HLA class II haplotypes in Siberian populations. More haplotypes have been found in Siberia than in any other population. Some of the haplotypes are shared with non-Siberian populations, but most of them are new, and some represent “forbidden” combinations of DQA1 and DQB1 alleles. We suggest that a set of “public” haplotypes was brought to Siberia with the colonizers but that most of the new haplotypes were generated in Siberia by recombination and are part of a haplotype pool that is turning over rapidly. The allelic frequencies at the DRB1 locus divide the Siberian populations into eastern and central Siberian branches; only the former shows a clear genealogical relationship to Amerinds. Received: 18 August 1997 / Accepted: 6 October 1997     

Composite origin of major histocompatibility complex genes

Major histocompatibility complex WHO genes have now been cloned from representatives of all vertebrate classes except Agnatha. The recent accumulation of sequence data has given great insight into the course of evolution of these genes. Although the primary structure of the MHC genes varies greatly from class to class and also within the individual classes, the general features of the tertiary and quaternary structure have been conserved remarkably well during more than 400 million years.of evolution. The ancestral MHC genes may have been assembled from at least three structural elements derived from different gene families. Class II MHC genes appear to have been assembled first, and then to have given rise to class I genes.     

Current perspectives on the intensity of natural selection of MHC loci

Polymorphism of genes in the major histocompatibility complex (MHC) is believed to be maintained by balancing selection. However, direct evidence of selection has proven difficult to demonstrate. In 1994, Satta and colleagues estimated the selection intensity of the human MHC (human leukocyte antigen (HLA)) loci; however, at that time the number of HLA sequences was limited. By comparing five different methods, this study demonstrated the best way to calculate the selection coefficient, through a computer simulation study. Since the study, many HLA nucleotide sequences have been made available. Our new analysis takes advantage of these newly available sequences and compares new estimates with those of the previous study. Generally, our new results are consistent with those of the 1994 study. Our results show that, even after 20 years of exhaustive sequencing of human HLA, the number of dominant HLA alleles, on which our original estimate of selection intensity depended, appears to be conserved. Indeed, according to the frequency distribution for each HLA allele, most sequences in the database were minor or private alleles; therefore, we conclude that the selection intensities of HLA loci are at most 4.4 % even though the HLA is the prominent example on which the natural selection has been operating.     

Nucleotide sequencing analysis of the swine 433-kb genomic segment located between the non-classical and classical<Emphasis Type="Italic"> SLA</Emphasis> class I gene clusters

Genome analysis of the swine leukocyte antigen (SLA) region is needed to obtain information on the MHC genomic sequence similarities and differences between the swine and human, given the possible use of swine organs for xenotransplantation. Here, the genomic sequences of a 433-kb segment located between the non-classical and classical SLA class I gene clusters were determined and analyzed for gene organization and contents of repetitive sequences. The genomic organization and diversity of this swine non-class I gene region was compared with the orthologous region of the human leukocyte antigen (HLA) complex. The length of the fully sequenced SLA genomic segment was 433 kb compared with 595 kb in the corresponding HLA class I region. This 162-kb difference in size between the swine and human genomic segments can be explained by indel activity, and the greater variety and density of repetitive sequences within the human MHC. Twenty-one swine genes with strong sequence similarity to the corresponding human genes were identified, with the gene order from the centromere to telomere of HCR - SPR1 - SEEK1 - CDSN - STG - DPCR1 - KIAA1885 - TFIIH - DDR - IER3 - FLOT1 - TUBB - KIAA0170 - NRM - KIAA1949 - DDX16 - FLJ13158 - MRPS18B - FB19 - ABCFI - CAT56. The human SEEK1 and DPCR1 genes are pseudogenes in swine. We conclude that the swine non-class I gene region that we have sequenced is highly conserved and therefore homologous to the corresponding region located between the HLA-C and HLA-E genes in the human.The nucleotide sequence data reported in this paper have been submitted to DDBJ, EMBL and GenBank databases under accession numbers AB113354, AB113355, AB113356, AB113357     

Genetic polymorphism of the swine major histocompatibility complex (SLA) class I genes, SLA-1, -2 and -3

In order to identify and characterize genetic polymorphism of the swine major histocompatibility complex (Mhc: SLA) class I genes, RT-PCR products of the second and third exons of the three SLA classical class I genes, SLA-1, SLA-2 and SLA-3 were subjected to nucleotide determination. These analyses allowed the identification of four, eight and seven alleles at the SLA-1, SLA-2 and SLA-3 loci, respectively, from three different breeds of miniature swine and one mixed breed. Among them, 12 alleles were novel. Construction of a phylogenetic tree using the nucleotide sequences of those 19 alleles indicated that the SLA-1 and -2 genes are more closely related to each other than to SLA-3. Selective forces operating at single amino acid sites of the SLA class I molecules were analyzed by the Adaptsite Package program. Ten positive selection sites were found at the putative antigen recognition sites (ARSs). Among the 14 positively selected sites observed in the human MHC (HLA) classical class I molecules, eight corresponding positions in the SLA class I molecules were inferred as positively selected. On the other hand, four amino acids at the putative ARSs were identified as negatively selected in the SLA class I molecules. These results suggest that selective forces operating in the SLA class I molecules are almost similar to those of the HLA class I molecules, although several functional sites for antigen and cytotoxic T-lymphocyte recognition by the SLA class I molecules may be different from those of the HLA class I molecules.The DNA sequence data reported in this paper have been submitted to the DDBJ, EMBL and GenBank nucleotide databases and have been assigned the accession numbers, AB105379, AB105380, AB105381, AB105382, AB105383, AB105384, AB105385, AB105386, AB105388, AB105389, AB105390 and AB105391     

Polymorphic SVA retrotransposons at four loci and their association with classical HLA class I alleles in Japanese, Caucasians and African Americans

Polymorphic insertion frequencies of the retrotransposons known as the “SVA” elements were investigated at four loci in the MHC class I genomic region to determine their allele and haplotype frequencies and associations with the HLA-A, -B or -C genes for 100 Japanese, 100 African Americans, 174 Australian Caucasians and 66 reference cell lines obtained from different ethnic groups. The SVA insertions representing different subfamily members varied in frequency between none for SVA-HF in Japanese and 65% for SVA-HB in Caucasians or African Americans with significant differences in frequencies between the three populations at least at three loci. The SVA loci were in Hardy–Weinberg equilibrium except for the SVA-HA locus which deviated significantly in African Americans and Caucasians possibly because of a genomic deletion of this locus in individuals with the HLA-A*24 allele. Strong linkage disequilibria and high percentage associations between the human leucocyte antigen (HLA) class I gene alleles and some of the SVA insertions were detected in all three populations in spite of significant frequency differences for the SVA and HLA class I alleles between the three populations. The highest percentage associations (>86%) were between SVA-HB and HLA-B*08, -B*27, -B*37 to -B*41, -B*52 and -B*53; SVA-HC and HLA-B*07; SVA-HA and HLA-A*03, -A*11 and -A*30; and SVA-HF and HLA-A*03 and HLA-B*47. From pairwise associations in the three populations and the homozygous cell line results, it was possible to deduce the SVA and HLA class I allelic combinations (haplotypes), population differences and the identity by descent of several common HLA-A allelic lineages.     

Inter- and Intralocus Recombination Drive MHC Class IIB Gene Diversification in a Teleost, the Three-Spined Stickleback Gasterosteus aculeatus

The mutational mechanism underlying the striking diversity in MHC (major histocompatibility complex) genes in vertebrates is still controversial. In order to evaluate the role of inter- and intragenic recombination in MHC gene diversification, we examined patterns of nucleotide polymorphism across an exon/intron boundary in a sample of 31 MHC class IIB sequences of three-spined stickleback (Gasterosteus aculeatus). MHC class IIB genes of G. aculeatus were previously shown to be under diversifying (positive) selection in mate choice and pathogen selection experiments. Based on recoding of alignment gaps, complete intron 2 sequences were grouped into three clusters using maximum-parsimony analysis. Two of these groups had >90% bootstrap support and were tentatively assigned single locus status. Intron nucleotide diversity within and among loci was low (p-distance within and among groups = 0.016 and 0.019, respectively) and fourfold lower than the rate of silent mutations in exon 2, suggesting that noncoding regions are homogenized by frequent interlocus recombination. A substitution analysis using GENECONV revealed as many intergenic conversion events as intragenic ones. Recombination between loci may explain the occurrence of sequence variants that are particularly divergent, as is the case in three-spined stickleback, with nucleotide diversity attaining d_N = 0.39 (peptide-binding residues only). For both MHC class II loci we also estimated the amount of intragenic recombination as population rate (4N_er) under the coalescent and found it to be approximately three times higher compared to point mutations (Watterson estimate per gene, 4N_eμ). Nonindependence of molecular evolution across loci and frequent recombination suggest that MHC class II genes of bony fish may follow different evolutionary dynamics than those of mammals. Our finding of widespread recombination suggests that phylogenies of MHC genes should not be based on coding segments but rather on noncoding introns. [Reviewing Editor: Dr. Richard Kliman]     

A class I jumping clone places the HLA-G gene approximately 100 kilobases from HLA-H within the HLA-A subregion of the human MHC

By the combination of cosmid cloning, chromosomal jumping, and pulsed-field gel electrophoresis (PFGE), we have fine-mapped the HLA-A subregion of the human major histocompatibility complex (MHC). Through the isolation of a class I jumping clone, the Qa-like HLA-G class I gene has been placed within 100 kb of HLA-H. The tight physical linkage of these class I genes has been further supported by hybridizing PFGE blots with locus-specific probes. It has been found that both of the above class I genes are linked to HLA-A, with HLA-H residing no more than 200 kb from the HLA-A gene. These data support the possible existence of a Qa-like subregion composed of nonclassical HLA class I genes within the human MHC linked telomerically to the HLA-A locus.     

Characterization and Phylogenetic Relationship of Prosimian MHC Class I Genes

MHC class I cDNA sequences from the most divergent primate group of extant primates compared to human, the suborder Strepsirrhini (prosimians), are described. The sequences are derived from the gray mouse lemur (Microcebus murinus) and the ring-tailed lemur (Lemur catta), which are members of the malagasy Lemuriformes, as well as from the pygmy slow loris (Nycticebus pygmaeus), a prosimian from East Asia. The M. murinus sequences have been analyzed in detail. Analysis of the expression level, G/C content, and synonymous vs. nonsynonymous substitution rates in the peptide-binding region codons suggests that these cDNA clones represent classical class I (class Ia) genes. According to Southern blot analysis, the genome of the gray mouse lemur might contain about 10 class I genes. In gene tree analysis, the strepsirrhine class Ia genes described here cluster significantly separately from the known class I genes of Catarrhini (humans, apes, Old World monkeys) and Platyrrhini (New World monkeys) species, suggesting that the class I loci of Simiiformes arose by gene duplications which occurred after the divergence of prosimians.     

Comparative genomic analysis of the major histocompatibility complex class I region in the teleost genus <Emphasis Type="Italic">Oryzias</Emphasis>

The major histocompatibility complex (MHC) class I region of teleosts harbors a tight cluster of the class IA genes and several other genes directly involved in class I antigen presentation. Moreover, the dichotomous haplotypic lineages (termed d- and N- lineages) of the proteasome subunit beta genes, PSMB8 and PSMB10, are present in this region of the medaka, Oryzias latipes. To understand the evolution of the Oryzias MHC class I region at the nucleotide sequence level, we analyzed bacterial artificial chromosome clones covering the MHC class I region containing the d- lineage of Oryzias luzonensis and the d- and N- lineages of Oryzias dancena. Comparison among these three elucidated sequences and the published sequences of the d- and N- lineages of O. latipes indicated that the order and orientation of the encoded genes were completely conserved among these five genomic regions, except for the class IA genes, which showed species-specific variation in copy number. The PSMB8 and PSMB10 genes showed trans-species dimorphism. The remaining regions flanking the PSMB10, PSMB8, and class IA genes showed high degrees of sequence conservation at interspecies as well as intraspecies levels. Thus, the three independent evolutionary patterns under apparently distinctive selective pressures are recognized in the Oryzias MHC class I region. Electronic Supplementary Material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A possible role of plantations for primate conservation in Madagascar

The utilization of eucalyptus plantations by seven sympatric species of prosimians was studied in the eastern rainforest of Madagascar. The species were Avahi laniger, Cheirogaleus major, Hapalemur griseus, Indri indri, Lemur fulvus, Lepilemur mustelinus, and Microcebus rufus. None of the lemurs was ever found in young eucalyptus plantations with little undergrowth. This was mainly due to the lack of travel opportunities within the shrub layer and between the shrubs and the canopy. Food (mainly berries) is seasonally available in the shrub layer but cannot be exploited because frugivorous lemurs cannot reach it. Old eucalyptus plantations with dense undergrowth are used by all prosimian species. They provide food as well as travel and resting facilities. Mixed tree plantations in the western part of Madagascar were used by groups of Lemur fulvus, Lepilemur mustelinus, and Propithecus verreauxi. According to these results, old eucalyptus plantations and mixed tree species plantations could be used to provide firewood and construction wood for the human population. They also might extend the habitat for lemurs and serve as buffers against human disturbance.     

A continuous restriction map from HLA-E to HLA-F. Structural comparison between different HLA-A haplotypes

The class I region of the human major histocompatibility complex contains genes encoding the classical transplantation antigens (HLA-A, B, and C), at least three new class I genes (HLA-E, F, and G) and many class I pseudogenes (including HLA-H). By pulse field gel electrophoresis and using five rare cutter enzymes, we have constructed a precise and continuous map of 1200 kilobases (kb) around HLA-A. The blots were hybridized with HLA-A, E, and F-specific probes and with new probes derived from yeast artificial chromosomes and cosmids of the class I region. We have compared the genomic organization of the same 1200 kb in three homozygous lymphoblastoid cell lines corresponding to three different HLA haplotypes (A3, A24, and A31). The differences in size observed may have been caused by insertions and deletions and may prove valuable in understanding the evolution of the HLA chromosomal region.     

Ancestral haplotypes reveal the role of the central MHC in the immunogenetics of IDDM

The major histocompatibility complex (MHC) contains multiple and diverse genes which may be relevant to the induction adn regulation of autoimmune responses in insulin dependent diabetes mellitus (IDDM). In addition to HLA class I and II, the possible candidates include TNF, C4, and several other poorly defined polymorphic genes in the central MHC region. This study describes two approaches which take advantage of the fact that the relevant genes are carried by highly conserved ancestral haplotypes such as 8.1 (HLA-B8, TNFS, C4AQO, C4B1, DR3, DQ2). First, three diabetogenic haplotypes (two Caucasoid and one Mongoloid) have been compared and it has been shown that all three share a rare allele of BAT3 as well as sharing DR3, DQ2. In 43 sequential patients with IDDM the cross product ration for BAT3S was 4.8 (p<0.01) and 6.9 for HLA-B8 plus BAT3S (p<0.001). Second, partial or recombinant ancestral haplotypes with either HLA class I (HLA-B8) or II (HLA-DR3, DQ2) alleles were identified. Third, using haplotypic polymorphisms such as the one in BAT3, we have shown that all the patients carrying recombinants of the 8.1 ancestral haplotype share the central region adjacent to HLA-B. These findings suggest that both HLA and non-HLA genes are involved in conferring susceptibility to IDDM, and that the region between HLA-B and BAT3 contains some of the relevant genes. By contrast, similar approaches suggest that protective genes map to the HLA class II region.