Heirs of the jaguar and the anaconda: HLA, conquest and disease in the indigenous populations of the Americas

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Variable NK cell receptors and their MHC class I ligands in immunity, reproduction and human evolution

Natural killer (NK) cells have roles in immunity and reproduction that are controlled by variable receptors that recognize MHC class I molecules. The variable NK cell receptors found in humans are specific to simian primates, in which they have progressively co-evolved with MHC class I molecules. The emergence of the MHC-C gene in hominids drove the evolution of a system of NK cell receptors for MHC-C molecules that is most elaborate in chimpanzees. By contrast, the human system of MHC-C receptors seems to have been subject to different selection pressures that have acted in competition on the immunological and reproductive functions of MHC class I molecules. We suggest that this compromise facilitated the development of the bigger brains that enabled archaic and modern humans to migrate out of Africa and populate other continents.

Analytical methods for immunogenetic population data

In this chapter, we describe analyses commonly applied to immunogenetic population data, along with software tools that are currently available to perform those analyses. Where possible, we focus on tools that have been developed specifically for the analysis of highly polymorphic immunogenetic data. These analytical methods serve both as a means to examine the appropriateness of a dataset for testing a specific hypothesis, as well as a means of testing hypotheses. Rather than treat this chapter as a protocol for analyzing any population dataset, each researcher and analyst should first consider their data, the possible analyses, and any available tools in light of the hypothesis being tested. The extent to which the data and analyses are appropriate to each other should be determined before any analyses are performed.

The Origin of Amerindians and the Peopling of the Americas According to HLA Genes: Admixture with Asian and Pacific People

The classical three-waves theory of American peopling through Beringia was based on a mixed anthropological and linguistic methodology. The use of mtDNA, Y chromosome and other DNA markers offers different results according to the different markers and methodologies chosen by different authors. At present, the peopling of Americas remains uncertain, regarding: time of population, number of peopling waves and place of peopling entrance among other related issues. In the present review, we have gathered most available HLA data already obtained about First Native American populations, which raise some doubts about the classical three waves of American peopling hypothesis. In summary, our conclusions are: 1) North West Canadian Athabaskans have had gene flow with: a) close neighboring populations, b) Amerindians, c) Pacific Islanders including East Australians and d) Siberians; 2) Beringia was probably not the only entrance of people to America: Pacific Ocean boat trips may have contributed to the HLA genetic American profile (or the opposite could also be true); 3) Amerindians entrance to America may have been different to that of Athabaskans and Eskimos and Amerindians may have been in their lands long before Athabaskans and Eskimos because they present and altogether different set of HLA-DRB1 allele frequencies; 4) Amerindians show very few “particular alleles”, almost all are shared with other Amerindians, Athabaskans and Pacific Islanders, including East Australians and Siberians; 5) Our results do not support the three waves model of American peopling, but another model where the people entrance is not only Beringia, but also Pacific Coast. Reverse migration (America to Asia) is not discarded and different movements of people in either direction in different times are supported by the Athabaskan population admixture with Asian-Pacific population and with Amerindians, 6) HLA variability is more common than allele veriability in Amerindians. Finally, it is shown that gene genealogy analises should be completed with allele frequency analyses in population relatednes and migrations studies.

Epidemiology of spondyloarthritis in North America

Many challenges have made it difficult to determine the prevalence of spondyloarthritis (SpA) in North America. They include the ethnic heterogeneity of the population, the lack of feasibility of applying current criteria (such as requirements for human leukocyte antigen-B27 testing and imaging studies such as pelvic radiographs and magnetic resonance imaging scanning) and the transient nature of some SpA symptoms (ie, peripheral arthritis and enthesitis). Current estimates of the prevalence of SpA in the United States range between 0.2% and 0.5% for ankylosing spondylitis, 0.1% for psoriatic arthritis, 0.065% for enteropathic peripheral arthritis, between 0.05% and 0.25% for enteropathic axial arthritis and an overall prevalence of SpA as high as >1%. With newer population-based instruments becoming available, the availability of the widely validated European Spondyloarthropathy Study Group criteria and the lower cost and greater feasibility of genetic testing, opportunities for true population-based studies of SpA are possible and will likely soon ensue.

HLA-E polymorphism in Amerindians from Mexico (Mazatecans), Colombia (Wayu) and Chile (Mapuches): evolution of MHC-E gene

Human leukocyte antigen (HLA)-E is a nonclassical class I (Ib) gene with a restricted polymorphism. Only eight DNA alleles and three proteins of this gene have been described and their frequencies analyzed in Caucasian, Oriental, Asian Indian, and Negroid populations. In the present study, HLA-E polymorphism has been analyzed in six Amerindian and Mestizo populations from North and South America and compared with previously described populations. HLA-E*0101 is the most frequent allele found in all populations except in Afrocolombian and Wayu Amerindians, in which blood group analyses show a high admixture with Caucasian and African populations. Mazatecan and Mapuche (two Amerindian groups from North and South America, respectively) presented similar HLA-E frequencies, whereas Wayu Indians are more similar to the Afrocolombian population. The Mexican and Colombian Mestizo show similar allele frequencies to Amerindians with high frequencies of HLA-E*0101 and HLA-E*010302 alleles. Also, frequencies in Negroids and Asian Indians present a similar distribution of HLA-E alleles. These data are in agreement with worldwide restricted polymorphism of HLA-E because no new allele was detected in the six populations studied. The allelic frequencies show differences among Caucasian, Oriental, Mestizo and Indian populations. Ape major histocompatibility complex-E allelism is also very restricted: common chimpanzee (one allele); bonobo (two alleles); gorilla (two alleles); orangutan (one allele); rhesus monkey (eight alleles); cynomolgus monkey (two alleles); and green monkey (two alleles).

Effects of race on survival after stem cell transplantation

Effects of race or ethnicity on survival after high-dose chemoradiation followed by stem cell transplantation (SCT) have not been thoroughly evaluated. We analyzed survival according to racial/ethnic categories for 3587 consecutive patients who had SCT at a single US institution between July 1992 and December 2000. Among 1366 patients who received autologous SCT, race or ethnicity was not significantly associated with survival. In contrast, among 2221 patients who received allogeneic SCT from HLA-matched unrelated or sibling donors, blacks had a significantly greater mortality than whites (unadjusted hazard ratio, 1.65; 95% confidence interval, 1.21-2.25). Mortality among other racial or ethnic groups was not significantly different from that among whites. The greater mortality hazard among blacks persisted after controlling for donor type, pretransplantation risk category, patient age, donor/patient sex, and cytomegalovirus exposure (hazard ratio, 1.71; 95% confidence interval, 1.25-2.34). SCT from both HLA-matched unrelated and HLA-identical sibling donors was associated with more severe acute graft-versus-host disease and higher nonrelapse mortality among blacks compared with whites. Furthermore, blacks who received SCT for chronic myeloid leukemia had longer diagnosis-to-transplantation intervals than whites. A matched-cohort analysis showed that the higher mortality among blacks could not be explained by obvious socioeconomic differences. The higher incidence of severe graft-versus-host disease among blacks compared with whites, both with HLA-identical sibling donors, might be related to yet-unidentified "immune-enhancing" genetic polymorphisms. We cannot exclude the possibility that the increased mortality risk among blacks after discharge from the transplant center might in part be related to unidentified sociocultural differences that influence medical care.

HLA alleles and haplotypes among the lakota sioux: report of the ASHI minority workshops, part III

Human leukocyte antigen (HLA) class I and II alleles were defined for 302 Lakota Sioux American Indians as part of the American Society for Histocompatibility and Immunogenetics coordinated studies on minority populations. The study group was comprised of adult volunteers from the Cheyenne River and Ogala Sioux tribes residing, respectively, on the Cheyenne River and Pine Ridge Reservations in South Dakota. Of the participants, 263 (87%) claimed full American Indian ancestry through both maternal and paternal grandparents. The study group included 25 nuclear families that were informative for genotyping. HLA phenotypes from 202 adults with no other known first-degree relative included in the study were used for calculation of allele and haplotype frequencies by maximum likelihood estimation. HLA-A, -B, and -Cw alleles were found to be in Hardy Weinberg equilibrium. Deviation from equilibrium was observed for DRB1 alleles (p=0.01), but could be attributed to the sample size and the occurrence of some genotypes with low expected frequencies. Polymorphism among the Sioux was limited with four to seven alleles comprising >80% of those observed at each locus. Several alleles were found at high frequency (0.05-0.30) among the Sioux that are also prevalent in other Native Americans and Alaska Natives, including: A*2402, *3101, and *0206; B*3501,*3901, *5101, and *2705; Cw*0702, *0404, and *03041; DRB1*0407, *0404, *1402, and *16021; and DQB1*0301, *0302, and *0402. DRB1*0811, which has been only previously described in Navajo and Tlingit Indians, was found to occur at a frequency of 0.119 among the Sioux. Two new alleles were defined among the Sioux: Cw*0204 and DRB1*040703, which were found in two and four individuals, respectively. In the haplotype analyses, significant linkage disequilibrium (p<0.00001) was seen in all pairwise comparisons of loci and numerous two and three locus haplotypes were found to have strong, positive linkage disequilibrium values. The two most common extended haplotypes among the Sioux, determined by maximum likelihood estimation and genotyping were: A*31012, B*3501, Cw*0404, DRB1*0407; and A*24021, B*3501, Cw*0404, DRB1*0404.

Sequence-based typing of HLA class I alleles in Alaskan Yupik Eskimo

In comparison to South America, native North Americans tend to be less diverse in their repertoire of HLA class I alleles. Based upon this observation, we hypothesized that the Yupik Eskimo would exhibit a limited number of previously identified class I HLA alleles. To test this hypothesis, sequence-based typing was performed at the HLA-A, -B and -C loci for 99 Central Yupik individuals from southwestern Alaska. Two new class I alleles, A*2423 and Cw*0806, were identified. While A*2423 was observed in only one sample, Cw*0806 was present in 26 of the 99 individuals and all of the Cw*0806 samples contained B*4801. Allele Cw*0806 differs from Cw*0803 by a single nucleotide substitution such that Cw*0803 may be the progenitor of Cw*0806. Allele Cw*0803 was originally characterized as unique to South America, but detection of Cw*0803 in the Yupik indicates that Cw*0803 was a founding allele of the Americas. The presence of new alleles and previously unrecognized founding alleles in the Yupik population show that natives of North America are more diverse than previously envisioned.

Novel HLA allele, HLA-B*4013, identified from a potential bone marrow recipient. Putative gene conversion events with donor sequence

A new B40 allele was identified in a leukemic Caucasian patient. This allele, designated B*4013, differs in alpha 1 domain from B*4002 at six amino acidic positions: 67, 77, 80, 81, 82 and 83. Most of this substitutions could alter the antigen binding site of the HLA-B molecule. B*4013 may have originated by gene conversion or reciprocal recombination involving B*4002 as the recipient allele of sequence donated by B*4406. The new allele was serologically typed as a "blank" associated with the Bw4 epitope.

Diversity is demonstrated in class I HLA-A and HLA-B alleles in Cameroon, Africa: description of HLA-A*03012, *2612, *3006 and HLA-B*1403, *4016, *4703

To examine the genetic diversity in west Africa, class I HLA-A and HLA-B alleles of 92 unrelated individuals from two areas in the Cameroon, the capital Yaoundé and the village of Etoa, were identified by direct automated DNA sequencing of exons 2 and 3 of the HLA-B locus alleles and sequence-specific oligonucleotide probe (SSOP) and/or sequencing of the HLA-A locus alleles. HLA-A*2301 (18.7%), A*2902 (10.4%), B*5301 (10.9%), and B*5802 (10.9%) were the most frequently detected alleles, present in at least 10% of the population. A total of 30 HLA-A locus and 33 HLA-B locus alleles, including six novel alleles, were detected. The novel alleles were HLA-A*03012, A*2612, A*3006 and HLA-B*1403, B*4016, and B*4703. HLA-B*4703 contains a novel amino acid sequence that is a combination of the first 5 amino acids of the Bw6 epitope and the last 2 residues of the Bw4 epitope. The addition of 6 alleles to the ever-expanding number of known class I HLA alleles supports our hypothesis that extensive genetic diversity, including previously undescribed alleles, would be observed in this African population. In the Yaoundé population, the allele frequency distribution at the HLA-A locus is consistent with distributions indicative of balancing selection. Extensive HLA-A-B haplotypes were observed in this population suggesting that only a fraction of the Cameroon HLA-A-B haplotype diversity has been observed.

Evolution of HLA-class I compared to HLA-class II polymorphism in Terena, a South-American Indian tribe

We have studied the HLA alleles of 60 unrelated healthy Terena and 10 Terena families. They are members of an isolated Brazilian tribe located in Mato Grosso do Sul (South Central Brazil). Six novel alleles were found in this population: HLA-A*0219 (gf = 0.02), A*0222 (gf = 0.15), HLA-B* 3520 (gf = 0.01), B*3521 (gf = 0.03), B*3912 (gf = 0.03) and B*4803 (gf = 0.16). Five of the six novel alleles differ from their putative progenitors by amino acid replacements in residues that contribute to the pockets of the peptide-binding site. Many of the variants defined by molecular methods were not identified correctly by serological typing. We calculated heterozygosity values (H) for HLA-A, -B, -C, DRB1, DQB1 and DPB . The highest values were observed at the HLA-B locus, followed by HLA-A, -DRB1 and DQB1. Residue positions 9, 24, 45, 62, 67, 95, 114, 116, 156, and 163 of HLA class I showed heterozygosity values greater than 0.50. Nine of them contribute to the peptide-binding specificity pockets and one to the T cell receptor binding site. If HLA antigens are useful for defense against pathogenic agents, heterozygosity would offer an advantage by allowing binding of a larger repertoire of peptides to the class I molecules. Individuals that are heterozygous at these positions would probably have a wider repertoire of peptide presentation to T cells. The observed results including the presence of novel alleles in the class I HLA loci suggest a functionally significant, more rapid evolution of class I compared to class II loci in this South American isolated population.

Novel HLA-A and HLA-B alleles in South American Indians

The human leukocyte antigen (HLA) complex includes the most polymorphic genes in humans. More than 600 allelic variants have been described in different populations. The HLA-B locus has contributed the largest number of alleles. Although Native American populations display a restricted number of HLA-alleles, many novel HLA class I alleles have been identified in indigenous communities of Central and South America. We have studied 248 unrelated individuals from three tribes of North-East Argentina and one from South-West Brazil, as well as 80 related individuals from the Brazilian tribe. In the course of this work, we found 8 new B-locus alleles and 2 novel A-locus alleles in these populations. Here we report the nucleotide sequences of A*0219, A*0222, B*3519, B*3520, B*3521, B*3912, B*4009 and B*4803 and we show their relationship with similar alleles. The new alleles B*35092 and B*3518 have been described by us in a previous paper. The possible mechanisms that may have produced these alleles over evolutionary time are discussed.

The South Amerindian allotype HLA-B*3909 has the largest known similarity in peptide specificity and common natural ligands with HLA-B27

HLA-B*3909 has only been found among South Amerindians, and presumably arose locally in these populations. It differs from B*3901 by a single Tyr to Ser change at position 99. To analyze the influence of this polymorphism on peptide specificity, pool sequence analysis and sequencing of multiple individual ligands from B*3901 and B*3909 were carried out. Both allotypes bind peptides with Arg2 or His2 and nonpolar C-terminal residues. However, whereas His2 is the predominant B*3901 motif, a majority of the B*3909-bound peptides have Arg2. In addition, B*3909 binds peptides with Pro2, and also shows an increased preference for Pro3. In spite of their differences, both subtypes bind overlapping peptide repertoires, as indicated by the identification of several identical ligands from their respective peptide pools. B*3909 is significantly more similar in its peptide specificity to HLA-B27 than B*3901. This is due to the increased preference of B*3909 for Arg2 and to low suitability of HLA-B27 for His2. The similarity between HLA-B27 and B*3909 was confirmed by identification of a natural ligand common to both allotypes. In addition, multiple HLA-B27 ligands bound efficiently B*3909 in vitro. The results indicate that, among the HLA class I allotypes of known peptide specificity, B*3909 is the most similar in its peptide binding properties to HLA-B27, which is absent in South Amerindians. This may have implications for the susceptibility of individuals in these populations to spondyloarthropathies.