HLA class II diversity in seven Amerindian populations. Clues about the origins of the Ache

Association between HLA Class II Alleles/Haplotypes and Genomic Ancestry in Brazilian Patients with Type 1 Diabetes: A Nationwide Exploratory Study

We aimed to identify HLA-DRB1, -DQA1, and -DQB1 alleles/haplotypes associated with European, African, or Native American genomic ancestry (GA) in admixed Brazilian patients with type 1 diabetes (T1D). This exploratory nationwide study enrolled 1599 participants. GA percentage was inferred using a panel of 46 ancestry informative marker-insertion/deletion. Receiver operating characteristic curve analysis (ROC) was applied to identify HLA class II alleles related to European, African, or Native American GA, and showed significant (p < 0.05) accuracy for identifying HLA risk alleles related to European GA: for DRB1*03:01, the area under the curve was (AUC) 0.533; for DRB1*04:01 AUC = 0.558, for DRB1*04:02 AUC = 0.545. A better accuracy for identifying African GA was observed for the risk allele DRB1*09:01AUC = 0.679 and for the protective alleles DRB1*03:02 AUC = 0.649, DRB1*11:02 AUC = 0.636, and DRB1*15:03 AUC = 0.690. Higher percentage of European GA was observed in patients with risk haplotypes (p < 0.05). African GA percentage was higher in patients with protective haplotypes (p < 0.05). Risk alleles and haplotypes were related to European GA and protective alleles/haplotypes to African GA. Future studies with other ancestry markers are warranted to fill the gap in knowledge regarding the genetic origin of T1D in highly admixed populations such as that found in Brazil.

Natural selection of immune and metabolic genes associated with health in two lowland Bolivian populations

A growing body of work has addressed human adaptations to diverse environments using genomic data, but few studies have connected putatively selected alleles to phenotypes, much less among underrepresented populations such as Amerindians. Studies of natural selection and genotype-phenotype relationships in underrepresented populations hold potential to uncover previously undescribed loci underlying evolutionarily and biomedically relevant traits. Here, we worked with the Tsimane and the Moseten, two Amerindian populations inhabiting the Bolivian lowlands. We focused most intensively on the Tsimane, because long-term anthropological work with this group has shown that they have a high burden of both macro and microparasites, as well as minimal cardiometabolic disease or dementia. We therefore generated genome-wide genotype data for Tsimane individuals to study natural selection, and paired this with blood mRNA-seq as well as cardiometabolic and immune biomarker data generated from a larger sample that included both populations. In the Tsimane, we identified 21 regions that are candidates for selective sweeps, as well as 5 immune traits that show evidence for polygenic selection (e.g., C-reactive protein levels and the response to coronaviruses). Genes overlapping candidate regions were strongly enriched for known involvement in immune-related traits, such as abundance of lymphocytes and eosinophils. Importantly, we were also able to draw on extensive phenotype information for the Tsimane and Moseten and link five regions (containing PSD4, MUC21 and MUC22, TOX2, ANXA6, and ABCA1) with biomarkers of immune and metabolic function. Together, our work highlights the utility of pairing evolutionary analyses with anthropological and biomedical data to gain insight into the genetic basis of health-related traits.

Distribution of a novel CYP2C haplotype in Native American populations

The CYP2C19 gene, located in the CYP2C cluster, encodes the major drug metabolism enzyme CYP2C19. This gene is highly polymorphic and no-function (CYP2C19*2 and CYP2C19*3), reduced function (CYP2C19*9) and increased function (CYP2C19*17) star alleles (haplotypes) are commonly used to predict CYP2C19 metabolic phenotypes. CYP2C19*17 and the genotype-predicted rapid (RM) and ultrarapid (UM) CYP2C19 metabolic phenotypes are absent or rare in several Native American populations. However, discordance between genotype-predicted and pharmacokinetically determined CYP2C19 phenotypes in Native American cohorts have been reported. Recently, a haplotype defined by rs2860840T and rs11188059G alleles in the CYP2C cluster has been shown to encode increased rate of metabolism of the CYP2C19 substrate escitalopram, to a similar extent as CYP2C19*17. We investigated the distribution of the CYP2C:TG haplotype and explored its potential impact on CYP2C19 metabolic activity in Native American populations. The study cohorts included individuals from the One Thousand Genomes Project AMR superpopulation (1 KG_AMR), the Human Genome Diversity Project (HGDP), and from indigenous populations living in Brazil (Kaingang and Guarani). The frequency range of the CYP2C:TG haplotype in the study cohorts, 0.469 to 0.598, is considerably higher than in all 1 KG superpopulations (range: 0.014-to 0.340). We suggest that the high frequency of the CYP2C:TG haplotype might contribute to the reported discordance between CYP2C19-predicted and pharmacokinetically verified CYP2C19 metabolic phenotypes in Native American cohorts. However, functional studies involving genotypic correlations with pharmacokinetic parameters are warranted to ascertain the importance of the CYP2C:TG haplotype.

Human leukocyte antigen Class II alleles associated with acral lentiginous melanoma in Mexican Mestizo patients: A case-control study

Background Melanoma is an aggressive cutaneous cancer. Acral lentiginous melanoma is a melanoma subtype arising on palms, soles, and nail-units. The incidence, prevalence and prognosis differ among populations. The link between expression of major histocompatibility complex Class II alleles and melanoma progression is known. However, available studies report variable results regarding the association of melanoma with specific HLA Class II loci. Aims The aim of the study was to determine HLA Class II allele frequencies in acral lentiginous melanoma patients and healthy Mexican Mestizo individuals. Methods Eighteen patients with acral lentiginous melanoma and 99 healthy controls were recruited. HLA Class II typing was performed based on the sequence-specific oligonucleotide method. Results Three alleles were associated with increased susceptibility to develop acral lentiginous melanoma, namely: HLA-DRB1*13:01; pC = 0.02, odds ratio = 6.1, IC95% = 1.4-25.5, HLA-DQA1*01:03; pC = 0.001, odds ratio = 9.3, IC95% = 2.7-31.3 and HLA-DQB1*02:02; pC = 0.01, odds ratio = 3.7, IC95% = 1.4-10.3. Limitations The small sample size was a major limitation, although it included all acral lentiginous melanoma patients seen at the dermatology department of Dr. Manuel Gea González General Hospital during the study period. Conclusion HLA-DRB1*13:01, HLA-DQB1*02:02 and HLA-DQA*01:03 alleles are associated with increased susceptibility to develop acral lentiginous melanoma in Mexican Mestizo patients.

Remarkably low KIR and HLA diversity in Amerindians reveals signatures of strong purifying selection shaping the centromeric KIR region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Amerindian genetic ancestry as a risk factor for tuberculosis in an amazonian population

In recent years, the incidence of tuberculosis (TB) has declined worldwide, although this disease still occurs at relatively high rates in Amerindian populations. This suggests that the genetic ancestry of Amerindians may be an important factor in the development of infections, and may account for at least some of the variation in infection rates in the different populations. The present study investigated the potential influence of Amerindian genetic ancestry on susceptibility to tuberculosis in an Amazon population. The study included 280 patients diagnosed with tuberculosis and 138 asymptomatic hospital employees with no history of TB, but who were in contact with bacterially active TB patients. Ancestry analysis was run on a set of 61 Ancestry-Informative Markers to estimate European, African, and Amerindian genetic ancestry using STRUCTURE v2.2. The TB group had significantly higher Amerindian ancestry in comparison with the control group, and significantly lower European ancestry. Amerindian ancestry in the 20–60% range was found to be the principal risk factor for increased susceptibility to TB. The results of the study indicate that Amerindian ancestry is an important risk factor for susceptibility to TB in the admixed population of the Brazilian Amazon region.

The immunogenetic diversity of the HLA system in Mexico correlates with underlying population genetic structure

We studied HLA class I (HLA-A, -B) and class II (HLA-DRB1, -DQB1) allele groups and alleles by PCR-SSP based typing in a total of 15,318 mixed ancestry Mexicans from all the states of the country divided into 78 sample sets, providing information regarding allelic and haplotypic frequencies and their linkage disequilibrium, as well as admixture estimates and genetic substructure. We identified the presence of 4268 unique HLA extended haplotypes across Mexico and find that the ten most frequent (HF > 1%) HLA haplotypes with significant linkage disequilibrium (Δ'≥0.1) in Mexico (accounting for 20% of the haplotypic diversity of the country) are of primarily Native American ancestry (A*02~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*08~DQB1*04, A*68~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*14~DQB1*03:01, A*24~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*04~DQB1*03:02, A*02~B*40:02~DRB1*04~DQB1*03:02, A*68~B*35~DRB1*04~DQB1*03:02, A*02~B*15:01~DRB1*04~DQB1*03:02). Admixture estimates obtained by a maximum likelihood method using HLA-A/-B/-DRB1 as genetic estimators revealed that the main genetic components in Mexico as a whole are Native American (ranging from 37.8% in the northern part of the country to 81.5% in the southeastern region) and European (ranging from 11.5% in the southeast to 62.6% in northern Mexico). African admixture ranged from 0.0 to 12.7% not following any specific pattern. We were able to detect three major immunogenetic clusters correlating with genetic diversity and differential admixture within Mexico: North, Central and Southeast, which is in accordance with previous reports using genome-wide data. Our findings provide insights into the population immunogenetic substructure of the whole country and add to the knowledge of mixed ancestry Latin American population genetics, important for disease association studies, detection of demographic signatures on population variation and improved allocation of public health resources.

Beyond the HLA polymorphism: A complex pattern of genetic susceptibility to pemphigus

Pemphigus is a group of autoimmune bullous skin diseases that result in significant morbidity. As for other multifactorial autoimmune disorders, environmental factors may trigger the disease in genetically susceptible individuals. The goals of this review are to summarize the state of knowledge about the genetic variation that may affect the susceptibility and pathogenesis of pemphigus vulgaris and pemphigus foliaceus – both the endemic and the sporadic forms –, to compare and discuss the possible meaning of the associations reported, and to propose recommendations for new research initiatives. Understanding how genetic variants translate into pathogenic mechanisms and phenotypes remains a mystery for most of the polymorphisms that contribute to disease susceptibility. However, genetic studies provide a strong foundation for further developments in this field by generating testable hypotheses. Currently, results still have limited influence on disease prevention and prognosis, drug development, and clinical practice, although the perspectives for future applications for the benefit of patients are encouraging. Recommendations for the continued advancement of our understanding as to the impact of genetic variation on pemphigus include these partially overlapping goals: (1) Querying the functional effect of genetic variants on the regulation of gene expression through their impact on the nucleotide sequence of cis regulatory DNA elements such as promoters and enhancers, the splicing of RNA, the structure of regulatory RNAs and proteins, binding of these regulatory molecules to regulatory DNA elements, and alteration of epigenetic marks; (2) identifying key cell types and cell states that are implicated in pemphigus pathogenesis and explore their functional genomes; (3) integrating structural and functional genomics data; (4) performing disease-progression longitudinal studies to disclose the causal relationships between genetic and epigenetic variation and intermediate disease phenotypes; (5) understanding the influence of genetic and epigenetic variation in the response to treatment and the severity of the disease; (6) exploring gene-gene and genotype-environment interactions; (7) developing improved pemphigus-prone and non-prone animal models that are appropriate for research about the mechanisms that link genotypes to pemphigus. Achieving these goals will demand larger samples of patients and controls and multisite collaborations.

Association between NDO-LID and PGL-1 for leprosy and class I and II human leukocyte antigen alleles in an indigenous community in Southwest Amazon

The frequencies of the Human leukocyte antigen (HLA) alleles in the Puyanawa indigenous reserve population and their association with the NDO-LID and ELISA PGL-1 rapid serological test was assessed. This was a cross-sectional study with an epidemiological clinical design conducted in two indigenous communities in the state of Acre, Brazil. Blood was collected in a tube with EDTA to identify HLA alleles and perform serological tests. DNA was obtained using the salting out procedure. The LabType™ technique (One-Lambda-USA) was used for HLA class I (loci A*, B* and C*) and II (loci DRB1*, DQA1* and DQB1*) typing. Allele frequency was obtained by direct count, and the chi-square test was used to assess the association with the NDO-LID and PGL-1 tests. The most frequent alleles in the two communities were: HLA-A*02:01, HLA-B*40:02, HLA-DRB1*16:02, HLA-DQA1*05:05 and HLA-DQB1*03:01. The allele HLA-C*04:01 was the most common in the Barão community, and the allele HLA-C*07:01 in Ipiranga. Among individuals who presented seropositivity to the NDO-LID test, the association with alleles HLA-A*02 (43.18% vs 24.8%, p = 0.03, OR = 2.35) and HLA-B*53 (6.83% vs 0.0%, p = 0.03, OR = 8.95) was observed in the Barão community. HLA-B*15 was associated with non-seroconversion to the NDO-LID test in Ipiranga. In both communities, HLA-B*40 and HLA-C*03 were associated with positive serological response to ELISA PGL-1. The HLA class I and II alleles most frequently found in this study have already been described among Terena indigenous groups, and HLA class I contributes to seroconversion to NDO-LID and PGL-1 tests in inhabitants of the Barão and Ipiranga communities.

High-resolution characterization of 12 classical and non-classical HLA loci in Southern Brazilians

The human leukocyte antigen (HLA) are the most polymorphic genes in the human genome. Because of their importance for antigen recognition, HLA molecules play a central role in host defense and graft rejection upon transplantation. The aim of this study was to characterize allelic diversity of the classical HLA genes HLA-A, -B, -C, -DRA, -DRB1, -DQA1, -DQB1, -DPA1, -DPB1, and the non-classical class I genes HLA-E, -F and -G at high-resolution for a population of predominantly European ancestry from Curitiba, Brazil. Genotyping of 108 individuals was performed by next-generation sequencing on the MiSeq platform and also by Sanger sequencing. The genotype distributions of all loci were in accordance with Hardy-Weinberg equilibrium (P > 0.05) and a total of 202 HLA variants at second field resolution were observed for the 12 loci. The strongest linkage disequilibrium (r² = 1.0, P < 10^-5 ) was observed for the following pairs of alleles: HLA-B*42:01:01 ~ HLA-DRB1*03:02:01; HLA-B*14:02:01 ~ HLA-C*08:02:01; B*42:01:01 ~ HLA-C*17:01:01; HLA-DRB1*03:01:01 ~ HLA-DQB1*02:01:01 ~ DRB1*03:01:01 ~ HLA-DQB1*02:01:01; DRB1*13:01:01~ HLA-DQB1*06:03:01 and HLA-DRB1*09:01:02 ~ HLA-DQA1*03:02. This is the first study to characterize all 12 HLA genes at high resolution in a single population. On the basis of the allelic frequencies of worldwide populations and principal component analysis, we confirmed the similarity of the study population to European and other Euro-descendant populations.

HLA diversity in Brazil

Brazil is the fifth largest country in the world in area and the fifth most populous. The Brazilian voluntary Bone Marrow Donor Registry is the third largest in terms of number of donors in the world, being a valuable source of HLA genetics to characterize the donor population of Brazil as well. The genetic background of the Brazilian population is quite heterogeneous, resulting from 5 centuries of admixture among Native Americans, Europeans and Africans, making the Brazilian population unique in terms of genetic ancestry. The unique characteristics of populations in different Brazilian regions make them an exciting focus for genetic diversity studies. Studies on HLA genetic diversity of Brazilian populations have been conducted since the late 1980s and, in this review, we highlight the main findings from studies carried out in Brazil based on classical HLA. In addition, we calculated the genetic distance from the molecular data of the studies included in this review in order to have a broader view of the HLA diversity in Brazilian populations. We emphasize that characterization of HLA diversity is not only important for transplantation programs, but can shed a light on ancestry, history and other demographic patterns with or without association with autoimmune disease.

Unveiling the Diversity of Immunoglobulin Heavy Constant Gamma (IGHG) Gene Segments in Brazilian Populations Reveals 28 Novel Alleles and Evidence of Gene Conversion and Natural Selection

Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene (IGH) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2, and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 (p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations (p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.

Impact of HLA-DRB1 allele polymorphisms on control of HIV infection in a Peruvian MSM cohort

Associations between HLA class II polymorphisms and HIV control were assessed in a Peruvian MSM cohort. Among 233 treatment naïve HIV+ individuals, DRB1*13:02 was linked to elevated viral loads (P = .044) while DRB1*12:01 showed significantly lower viral set points (P = .015) and restricted a dominant T cell response to HIV Gag p24 (P = .038). The present work contributes to a better knowledge of the Peruvian immunogenetics and supports the important role of HLA class II restricted T cells in HIV control.

HLA-A, B, DRB1, DQA1, DQB1 alleles and haplotype frequencies in Dene and Cree cohorts in Manitoba, Canada

BACKGROUND

First Nations in the Canadian province of Manitoba have disproportionately high rates of epidemic and endemic TB. Gene polymorphisms that modulate HLA Class I and II antigens are among the risk markers for TB, along with other biologic, and social determinants of health. HLA-A, B, DRB1, DQA1, DQB1 were typed in two Manitoba First Nation indigenous groups to identify and compare the frequency of gene polymorphisms that may influence susceptibility or resistance to TB.

METHODS

Participants who self-identified as either Dene or Cree enrolled into the study from two First Nation communities in Manitoba, Canada. Genomic DNA was extracted from blood samples collected with informed consent from Dene (N=63) and Cree (N=42) First Nation study participants. Participants self-reported having treated active TB, treated latent TB or no TB. HLA Class I and II molecules were typed using sequence-specific oligonucleotide (SSO) probes from commercially available kits.

RESULTS

The rates of treated active and latent TB were marginally higher among the Dene than the Cree participants (p=0.112). Class I and II HLA loci were in Hardy-Weinberg equilibrium in both the Dene and Cree groups. In this exploratory analysis of TB and HLA allele frequencies in Dene and Cree cohorts HLA-A*03 and HLA-DQB1*05:03 were significantly associated with TB.

CONCLUSIONS

The high incidence of TB in both Dene and Cree populations in Canada requires both biomedical and socioeconomic prevention and control measures. Among the former, an understanding of HLA diversity among First Nations groups may aid the development of new effective vaccine and therapeutic modalities that depend on the interaction between small molecules and specific HLA epitopes.

Association between HLA genes and dust mite sensitivity in a Brazilian population

BACKGROUND

Type I hypersensitivity, also known as IgE-mediated allergy, is a complex, multifactorial condition whose onset and severity are influenced by both genetic and environmental factors. Mite allergens stimulate the production of humoral response (IgE), especially in children, which is closely involved in atopic asthma and rhinitis.

OBJECTIVE

This study aimed to investigate the association between HLA class I (-A, -B, and -C), and HLA class II (-DRB1) genes in individuals sensitive to dust mites (Dermatophagoides farinae, Dermatophagoides pteronyssinus, or Blomia tropicalis) and mite-insensitive controls.

METHODS

396 participants were grouped as mite-sensitive and mite-insensitive according to immediate hypersensitivity as determined by skin-prick tests, and to HLA genotyping by polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO).

RESULTS

After chi-square heterogeneity testing no significant differences were observed in HLA-A, B, and C genes, except for the HLA-DRB1 locus, which, showed a negative association for DRB1∗04, between mite-sensitive and mite-insensitive individuals. In high resolution, DRB1∗04:11 allele was significantly different from all other results (P=0.0042, OR=0.26, and 95%CI=0.09-0.70). The analysis stratified by etiologic agent confirmed these associations.

CONCLUSION

Our results suggest a possible association between HLA-DRB1 genes and hypersensitivity to dust mites.

The role of variants from the innate immune system genes in tuberculosis and skin test response in a Native American population

Native American populations show higher tuberculosis (TB) mortality and infectivity rates than non-Native populations. Variants in the innate immune system seem to have an important role on TB susceptibility. The role of some innate immune system variants in TB susceptibility and/or skin test response (PPD) were investigated in the Aché, a Native American population. Complement receptor 1 and toll like receptor 9 variants were associated with anergy to PPD and protection to TB, respectively. These findings demonstrate an important role of the innate immune system variants in TB susceptibility.

Genetics, genomics, and evolutionary biology of NKG2D ligands

Human and mouse NKG2D ligands (NKG2DLs) are absent or only poorly expressed by most normal cells but are upregulated by cell stress, hence, alerting the immune system in case of malignancy or infection. Although these ligands are numerous and highly variable (at genetic, genomic, structural, and biochemical levels), they all belong to the major histocompatibility complex class I gene superfamily and bind to a single, invariant, receptor: NKG2D. NKG2D (CD314) is an activating receptor expressed on NK cells and subsets of T cells that have a key role in the recognition and lysis of infected and tumor cells. Here, we review the molecular diversity of NKG2DLs, discuss the increasing appreciation of their roles in a variety of medical conditions, and propose several explanations for the evolutionary force(s) that seem to drive the multiplicity and diversity of NKG2DLs while maintaining their interaction with a single invariant receptor.

Alleles of HLA-DRB1*04 Associated with Pulmonary Tuberculosis in Amazon Brazilian Population

Immunogenetic host factors are associated with susceptibility or protection to tuberculosis (TB). Strong associations of HLA class II genes with TB are reported. We analyzed the HLA-DRB1*04 alleles to identify subtypes associated with pulmonary TB and their interaction with risk factors such as alcohol, smoking, and gender in 316 pulmonary TB patients and 306 healthy individuals from the Brazilian Amazon. The HLA-DRB1*04 was prevalent in patients with pulmonary TB (p<0.0001; OR = 2.94; 95% CI = 2.12 to 4.08). Direct nucleotide sequencing of DRB1 exon 2 identified nine subtypes of HLA-DRB1*04. The subtype HLA-DRB1*04:11:01 (p = 0.0019; OR = 2.23; 95% CI = 1.34 to 3.70) was associated with susceptibility to pulmonary TB while DRB1*04:07:01 (p<0.0001; OR = 0.02; 95% CI = 0.001 to 0.33) to protection. Notably, the interaction between alcohol and HLA-DRB1*04:11:01 increased the risk for developing pulmonary TB (p = 0.0001; OR = 51.3; 95% CI = 6.81 to 386). Multibacillary pulmonary TB, the clinical presentation of disease transmission, was strongly associated with interaction to alcohol (p = 0.0026; OR = 11.1; 95% CI = 3.99 to 30.9), HLA-DRB1*04:11:01 (p = 0.0442; OR = 2.01; 95% CI = 1.03 to 3.93) and DRB1*04:92 (p = 0.0112; OR = 8.62; 95% CI = 1.63 to 45.5). These results show that HLA-DRB1*04 are associated with pulmonary TB. Interestingly, three subtypes, DRB1*04:07:01, DRB1*04:11:01 and DRB1*04:92 of the HLA-DRB1*04 could be potential immunogenetic markers that may help to explain mechanisms involved in disease development.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology

OBJECTIVES

The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations.

MATERIALS AND METHODS

We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database.

RESULTS

Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure.

DISCUSSION

Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Genetic Variations in the TP53 Pathway in Native Americans Strongly Suggest Adaptation to the High Altitudes of the Andes

The diversity of the five single nucleotide polymorphisms located in genes of the TP53 pathway (TP53, rs1042522; MDM2, rs2279744; MDM4, rs1563828; USP7, rs1529916; and LIF, rs929271) were studied in a total of 282 individuals belonging to Quechua, Aymara, Chivay, Cabanaconde, Yanke, Taquile, Amantani, Anapia, Uros, Guarani Ñandeva, and Guarani Kaiowá populations, characterized as Native American or as having a high level (> 90%) of Native American ancestry. In addition, published data pertaining to 100 persons from five other Native American populations (Surui, Karitiana, Maya, Pima, and Piapoco) were analyzed. The populations were classified as living in high altitude (≥ 2,500 m) or in lowlands (< 2,500 m). Our analyses revealed that alleles USP7-G, LIF-T, and MDM2-T showed significant evidence that they were selected for in relation to harsh environmental variables related to high altitudes. Our results show for the first time that alleles of classical TP53 network genes have been evolutionary co-opted for the successful human colonization of the Andes.

HLA class II alleles in the Otomi population of the Mezquital Valley: a genetic approach to the history of interethnic migrations in the Mexican Central Plateau

From a historical and genetic point of view, the Otomi of the Mezquital Valley are a frontier people that have played an important role in the population dynamics of the Mexican Central Plateau. Due to the antiquity of their presence in the area, the Otomi may be bearers of ancient genetic variability, shared mainly today with other groups belonging to the Otomanguean linguistic family and with the Nahua. In this study we analyzed the HLA class II allele frequencies reported in Mexican indigenous populations, in order to provide an intraregional-level historical perspective of the genetic relationships between the Otomi of the Mezquital Valley and indigenous populations from other regions of Mexico. We examined genetic variation in HLA-DRB1 and -DQB1 loci in 66 nonrelated individuals belonging to seven indigenous communities from the Ixmiquilpan municipality in the Mezquital Valley, in the State of Hidalgo, Mexico. The variability of the HLA-DRB1 gene among the Otomi of the Mezquital Valley was mainly concentrated in five alleles: -DRB1*08:02 (31.06%), -DRB1*04:07 (25.77%), -DRB1*14:06 (7.55%), -DRB1*14:02 (6.06%), and -DRB1*16:02 (4.55%); these alleles have been previously described in other indigenous populations. The most frequent alleles at the HLA-DQB1 locus were -DQB1*03:02 (34.09%), -DQB1*04:02 (31.03%), and -DQB1*03:01 (19.7%). Furthermore, the HLA-DQB1*02:02 allele was found in the Otomi group with a frequency of 2.27%; this allele has not been reported in Mexican indigenous populations. In conclusion, the genetic constitution of the Otomi population is intermediate to the northern groups and the genetic variability shared by the peoples of the central regions of Mexico. Furthermore, HLA-DRB1 and -DQB1 allelic variability among the Otomi provides insight into the historical processes implied in the biological admixture with European, Asian, and African populations as well as in the admixture with the population of Mexico City associated with long-standing migratory processes.

Polymorphisms of HLA-DRB1, -DQA1 and -DQB1 in inhabitants of Astana, the capital city of Kazakhstan

BACKGROUND

Kazakhstan has been inhabited by different populations, such as the Kazakh, Kyrgyz, Uzbek and others. Here we investigate allelic and haplotypic polymorphisms of human leukocyte antigen (HLA) genes at DRB1, DQA1 and DQB1 loci in the Kazakh ethnic group, and their genetic relationship between world populations.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 157 unrelated Kazakh ethnic individuals from Astana were genotyped using sequence based typing (SBT-Method) for HLA-DRB1, -DQA1 and -DQB1 loci. Allele frequencies, neighbor-joining method, and multidimensional scaling analysis have been obtained for comparison with other world populations. Statistical analyses were performed using Arlequin v3.11. Applying the software PAST v. 2.17 the resulting genetic distance matrix was used for a multidimensional scaling analysis (MDS). Respectively 37, 17 and 19 alleles were observed at HLA-DRB1, -DQA1 and -DQB1 loci. The most frequent alleles were HLA-DRB1*07:01 (13.1%), HLA-DQA1*03:01 (13.1%) and HLA-DQB1*03:01 (17.6%). In the observed group of Kazakhs DRB1*07:01-DQA1*02:01-DQB1*02:01 (8.0%) was the most common three loci haplotype. DRB1*10:01-DQB1*05:01 showed the strongest linkage disequilibrium. The Kazakh population shows genetic kinship with the Kazakhs from China, Uyghurs, Mongolians, Todzhinians, Tuvinians and as well as with other Siberians and Asians.

CONCLUSIONS/SIGNIFICANCE

The HLA-DRB1, -DQA1 and -DQB1 loci are highly polymorphic in the Kazakh population, and this population has the closest relationship with other Asian and Siberian populations.

Frequencies of allele groups HLA-A, HLA-B and HLA-DRB1 in a population from the northwestern region of São Paulo State, Brazil

The aim of this study was to estimate the HLA-A, HLA-B and HLA-DRB1 allele groups frequencies in a population of 1559 volunteer bone marrow donors from the northwestern region of São Paulo State grouped according to ethnicity. An additional objective was to compare the allele frequencies of the current study with data published for other Brazilian populations. The allele groups were characterized by the PCR-rSSO method using Luminex(®) technology. Twenty HLA-A, 32 HLA-B and 13 HLA-DRB1 allele groups were identified. The most common allele groups in European descent and mixed African and European descent samples were HLA-A*02, HLA-B*35 and HLA-DRB1*13, while HLA-A*02, HLA-B*35 and HLA-DRB1*11 were more common in African descent samples. The HLA-A*23, HLA-A*36, HLA-B*58 and HLA-B*81 allele groups were more common in sample from African descent than European descent, and the HLA-DRB1*08 was more common in mixed African and European descent than in European descent. Allele group frequencies were compared with samples from other Brazilian regions. The HLA-A*30 and HLA-A*23 were more common in this study than in the populations of Rio Grande do Sul and Paraná; and the HLA-A*01, HLA-B*18, HLA-B*57 and HLA-DRB1*11 were more common in this study than in the population of Piauí. The least frequent allele groups were HLA-A*31, HLA-B*15, HLA-B*40 and HLA-DRB1*08 for the population of Piauí, HLA-A*01 and HLA-A*11 for Parana, HLA-A*02 and -A*03 for Rio Grande do Sul and HLA-DRB1*04 for Paraná, Rio Grande do Sul and Piauí. These data provide an overview on the knowledge on HLA diversity in the population of the northwestern region of São Paulo State and show that the genes of this system are useful to distinguish different ethnic groups.

Genetic diversity of the HLA-G coding region in Amerindian populations from the Brazilian Amazon: a possible role of natural selection

HLA-G has an important role in the modulation of the maternal immune system during pregnancy, and evidence that balancing selection acts in the promoter and 3'UTR regions has been previously reported. To determine whether selection acts on the HLA-G coding region in the Amazon Rainforest, exons 2, 3 and 4 were analyzed in a sample of 142 Amerindians from nine villages of five isolated tribes that inhabit the Central Amazon. Six previously described single-nucleotide polymorphisms (SNPs) were identified and the Expectation-Maximization (EM) and PHASE algorithms were used to computationally reconstruct SNP haplotypes (HLA-G alleles). A new HLA-G allele, which originated in Amerindian populations by a crossing-over event between two widespread HLA-G alleles, was identified in 18 individuals. Neutrality tests evidenced that natural selection has a complex part in the HLA-G coding region. Although balancing selection is the type of selection that shapes variability at a local level (Native American populations), we have also shown that purifying selection may occur on a worldwide scale. Moreover, the balancing selection does not seem to act on the coding region as strongly as it acts on the flanking regulatory regions, and such coding signature may actually reflect a hitchhiking effect.

HLA class I and class II conserved extended haplotypes and their fragments or blocks in Mexicans: implications for the study of genetic diversity in admixed populations

Major histocompatibility complex (MHC) genes are highly polymorphic and informative in disease association, transplantation, and population genetics studies with particular importance in the understanding of human population diversity and evolution. The aim of this study was to describe the HLA diversity in Mexican admixed individuals. We studied the polymorphism of MHC class I (HLA-A, -B, -C), and class II (HLA-DRB1, -DQB1) genes using high-resolution sequence based typing (SBT) method and we structured the blocks and conserved extended haplotypes (CEHs) in 234 non-related admixed Mexican individuals (468 haplotypes) by a maximum likelihood method. We found that HLA blocks and CEHs are primarily from Amerindian and Caucasian origin, with smaller participation of African and recent Asian ancestry, demonstrating a great diversity of HLA blocks and CEHs in Mexicans from the central area of Mexico. We also analyzed the degree of admixture in this group using short tandem repeats (STRs) and HLA-B that correlated with the frequency of most probable ancestral HLA-C/−B and -DRB1/−DQB1 blocks and CEHs. Our results contribute to the analysis of the diversity and ancestral contribution of HLA class I and HLA class II alleles and haplotypes of Mexican admixed individuals from Mexico City. This work will help as a reference to improve future studies in Mexicans regarding allotransplantation, immune responses and disease associations.

MHC Class II haplotypes of Colombian Amerindian tribes

We analyzed 1041 individuals belonging to 17 Amerindian tribes of Colombia, Chimila, Bari and Tunebo (Chibcha linguistic family), Embera, Waunana (Choco linguistic family), Puinave and Nukak (Maku-Puinave linguistic families), Cubeo, Guanano, Tucano, Desano and Piratapuyo (Tukano linguistic family), Guahibo and Guayabero (Guayabero Linguistic Family), Curripaco and Piapoco (Arawak linguistic family) and Yucpa (Karib linguistic family). for MHC class II haplotypes (HLA-DRB1, DQA1, DQB1). Approximately 90% of the MHC class II haplotypes found among these tribes are haplotypes frequently encountered in other Amerindian tribes. Nonetheless, striking differences were observed among Chibcha and non-Chibcha speaking tribes. The DRB1*04:04, DRB1*04:11, DRB1*09:01 carrying haplotypes were frequently found among non-Chibcha speaking tribes, while the DRB1*04:07 haplotype showed significant frequencies among Chibcha speaking tribes, and only marginal frequencies among non-Chibcha speaking tribes. Our results suggest that the differences in MHC class II haplotype frequency found among Chibcha and non-Chibcha speaking tribes could be due to genetic differentiation in Mesoamerica of the ancestral Amerindian population into Chibcha and non-Chibcha speaking populations before they entered into South America.

KIR gene content in amerindians indicates influence of demographic factors

Although the KIR gene content polymorphism has been studied worldwide, only a few isolated or Amerindian populations have been analyzed. This extremely diverse gene family codifies receptors that are expressed mainly in NK cells and bind HLA class I molecules. KIR-HLA combinations have been associated to several diseases and population studies are important to comprehend their evolution and their role in immunity. Here we analyzed, by PCR-SSP (specific sequencing priming), 327 individuals from four isolated groups of two of the most important Brazilian Amerindian populations: Kaingang and Guarani. The pattern of KIR diversity among these and other ten Amerindian populations disclosed a wide range of variation for both KIR haplotypes and gene frequencies, indicating that demographic factors, such as bottleneck and founder effects, were the most important evolutionary factors in shaping the KIR polymorphism in these populations.

HLA-A, B and DRB1 allele and haplotype frequencies in volunteer bone marrow donors from the north of Parana State

Background

Knowledge of allele and haplotype frequencies of the human leukocyte antigen (HLA) system is important in the search for unrelated bone marrow donors. The Brazilian population is very heterogeneous and the HLA system is highly informative of populations because of the high level of polymorphisms.

Aim

The aim of this study was to characterize the immunogenetic profile of ethnic groups (Caucasians, Afro-Brazilians and Asians) in the north of Parana State.

Methods

A study was carried out of 3978 voluntary bone marrow donors registered in the Brazilian National Bone Marrow Donor Registry and typed for the HLA-A, B and DRB1 (low resolution) loci. The alleles were characterized by the polymerase chain reaction sequence-specific oligonucleotides method using the LabType SSO kit (One Lambda, CA, USA). The ARLEQUIN v.3.11 computer program was used to calculate allele and haplotype frequencies

Results

The most common alleles found in Caucasians were HLA-A*02, 24, 01; HLA-B*35, 44, 51; DRB1*11, 13, 07; for Afro-Brazilians they were HLA-A*02, 03, 30; HLA-B*35, 15, 44; DRB1*13, 11, 03; and for Asians they were: HLA-A*24, 02, 26; HLA-B*40, 51, 52; DRB1*04, 15, 09. The most common haplotype combinations were: HLA-A*01, B*08, DRB1*03 and HLA-A*29, B*44, DRB1*07 for Caucasians; HLA-A*29, B*44, DRB1*07 and HLA-A*01, B*08 and DRB1*03 for Afro-Brazilians; and HLA-A*24, B*52, DRB1*15 and HLA-A*24, B*40 and DRB1*09 for Asians.

Conclusion

There is a need to target and expand bone marrow donor campaigns in the north of Parana State. The data of this study may be used as a reference by the Instituto Nacional de Cancer/Brazilian National Bone Marrow Donor Registry to evaluate the immunogenetic profile of populations in specific regions and in the selection of bone marrow donors

Celiac disease in native Indians from Brazil: A clinical and epidemiological survey

Background:

Celiac disease has been described in populations from around the world, with recent data emphasizing the occurrence of the disease in ethnic minorities. There are only a few studies evaluating celiac disease in native Indians.

Aims:

This study aimed to screen the anti-endomysial antibody (IgA-EmA) in Kaingang and Guarani Indians from southern Brazil, in order to establish a clinical serological evaluation of celiac disease in these individuals.

Material and Methods:

Serum samples from 321 individuals (125 male and 196 female; 4-86 years old) from Mangueirinha Indigenous Reserve, State of Parana, Brazil, and 180 non-Indigenous healthy individuals (62 male and 118 female; 2-81 years old) were analysed to the presence of anti-endomysial antibody class IgA by indirect immunofluorescence assay. Amongst the Indians, 158 were Kaingang, 98 Guarani and 65 of mixed race. Indians presenting complaints of diarrhea (N=12) were also evaluated to the IgG class of anti-endomisyal antibody.

Results:

None of the individuals showed positive results either to IgA or IgG anti-endomysial antibodies.

Conclusions:

Although the results indicate an absence of celiac disease in Kaingang and Guarani Indians, the authors call attention to the importance of following up indigenous children or adults presenting gastrointestinal complaints or other symptoms related to the disease. Consideration should be given to the genetic background of these individuals, allied to the inter ethnic marriages and the changing habits or occupational activities, that have gradually introduced diseases previously not described in indigenous populations.

Impact of population diversity on the prediction of 7-SNP NAT2 phenotypes using the tagSNP rs1495741 or paired SNPs

A novel NAT2 tagSNP (rs1495741) and a 2-SNP genotype (rs1041983 and rs1801280) have been recently shown to accurately predict the NAT2 acetylator phenotypes in populations of exclusive or predominant European/White ancestry. We confirmed the accuracy of the tagSNP approach in White Brazilians, but not in Brown or Black Brazilians, sub-Saharan Mozambicans, and Guarani Amerindians. The combined rs1041983 and rs1801280 genotypes provided considerably better prediction of the NAT2 phenotype in Guarani, but no consistent improvement in Brown or Black Brazilians and Mozambicans. Best predictions of the NAT2 phenotype in Mozambicans using NAT2 SNP pairs were obtained with rs1801280 and rs1799930, but the accuracy of the estimates remained inadequate for clinical use or for investigations in this sub-Saharan group or in Brazilians with considerable African ancestry. In conclusion, the rs1495741 tagSNP cannot be applied to predict the NAT2 acetylation phenotype in Guarani and African-derived populations, whereas 2-SNP genotypes may accurately predict NAT2 phenotypes in Guarani, but not in Africans.

HLA polymorphism in a Guarani-Indian population from Paraguay and its usefulness for the Hispano-Indian admixture study in Paraguay

In this study we investigated the human leucocyte antigen-A (HLA-A), -B and DRB1 polymorphism of Native American population of Paraguay, the Guarani Indians. We found that the HLA variability consisted of 5 HLA-A, 7 HLA-B and 6 HLA-DRB1 groups of alleles and of several specific alleles (B*1504, B*3505, B*3912, B*4004, B*5104, DRB1*0411, DRB1*1413) common in other Native American populations. The comparison of the HLA polymorphism of the Guaranis from Paraguay with the «Mestizos» of Paraguay and the Spaniards showed that the «Mestizos» of Paraguay are genetically very distant from the Guarani Indians of Paraguay but much more close to the Spaniards. This can be explained, at least in part, by the history of the country. Our results are of importance in transplantation, in particular in the search for an unrelated donor for a Paraguayan patient requiring hematopoietic stem cell transplantation.

Variability of the BCHE gene in Amerindians from Paraná, Brazil

OBJECTIVE

This study aims to better understand the variability of the BCHE gene, and regions adjacent to it, in different populations. More specifically, it also aims to analyze the diversity of this DNA segment among two Amerindian populations, whose distinct characteristics provide us an extremely interesting material for genetic, evolutionary and anthropological studies, since they diverge genetically and culturally in spite of a very close living relationship in space and time.

METHODS

The two Amerindian populations analysed were from Guarani-Mbyá and Kaingang ethnic groups, both from the Rio das Cobras indian reservation area (Paraná, Brazil). The other populations data were obtained from the HapMap Project data, NCBI data and previously studies. The variability of seven SNPs was analyzed: two downstream, three within and two upstream of the BCHE gene.

RESULTS

The number of different BCHE variants in Amerindians is lower than in other populations. The present data are in accordance with the genetic distance between Guarani and Kaingang, and from them to other populations. Some SNPs (rs3495, rs7624915 and rs4387996) presented an interesting frequency pattern, with a higher frequency of the ancestral allele in the region of the modern human origin and its frequency decreasing gradually, from Africa toward the Americas.

CONCLUSIONS

Despite the geographical coexistence throughout their dispersion, each of these two Amerindian groups retains a high degree of genetic identity, probably maintained by cultural and social isolation.

Genetic link among Hani, Bulang and other Southeast Asian populations: evidence from HLA -A, -B, -C, -DRB1 genes and haplotypes distribution

Human leucocyte antigen (HLA) genetic characteristic of different ethnic minorities would be useful for tracing the origin of modern human HLA matching in transplantation and disease associations. In this study, we reported HLA-A, -B, -C and -DRB1 alleles and haplotypes in the Bulang and Hani populations of southwestern China using a high-resolution polymerase chain reaction-Luminex (PCR-Luminex) typing method. A total of nine HLA-A, 22 HLA-B, 13 HLA-C and 18 HLA-DRB1 alleles were identified in the Bulang population, and 11 HLA-A, 25 HLA-B, 14 HLA-C and 19 HLA-DRB1 alleles were in the Hani population. Compared with other populations, the predominant A*1101-B*1502-DRB1*1202 haplotype in the Bulang and Hani populations was also common in Jinuo, Wa, Dai, Maonan and Vietnamese populations. The distribution of HLA genes indicate that Bulang and Hani populations belong to the Southeast Asia group and they have maintained their original genetic characteristics and kept a long genetic distance from other populations owing to founder effects and subsequent geographic isolation. In addition, the close relationship among ethnic groups in Yunnan province and the Thai and Vietnamese populations were confirmed, whereas the Mon-Kmer-speaking populations shared other common HLA alleles and haplotypes compared with other linguistic groups.

Dopamine receptor D4 allele distribution in Amerindians: a reflection of past behavior differences?

The DRD4 variable number of tandem repeats (VNTR) allele distribution of 172 Guarani (Kaiowá and Ñandeva subgroups) and Kaingang Brazilian Amerindians is reported. These results are integrated with those previously obtained for this ethnic group. Allele frequencies for the three populations are within the interval observed for 15 other Native American populations and show intermediate values between those observed in Amazonia and Patagonia. Significant differences in allele distribution between recent past hunter-gatherer and agriculturalist populations are observed, with an increase of the 7R allele among hunter-gatherers (P < 0.001). Analysis of molecular variance (AMOVA) and pairwise F(ST) data suggest three distinct sectors for the genetic landscape of Native South America: Andes, Center/Southeast region, and Amazonia. Common traits among hunter-gatherers such as novelty-seeking temperament, hyperactivity, and impulsivity could have been important and advantageous in new environments during America's prehistoric colonization.

Variations in human herpesvirus type 8 seroprevalence in Native Americans, South America

To determine the epidemiology of human herpesvirus type 8 (HHV-8) among non-Amazonian native populations, we conducted a cross-sectional study in Brazil, Bolivia, and Paraquay. Our data show striking ethnic and geographic variations in the distribution of HHV-8 seroprevalences in Amazonian (77%) and non-Amazonian native populations (range 0%–83%).

Distribution of human leukocyte antigen alleles and haplotypes in Oroqen and Ewenki nationality minority in Inner Mongolia Autonomous Region of China

The frequencies of the human leukocyte antigen alleles HLA-A,-B, DRB1 and the A-B, A-DRB1, B-DRB1, A-B-DRB1 haplotypes were investigated through means of PCR-based reverse line-strip sequence specific oligonucleotide hybridization on 108 Oroqen and 104 Ewenki nationality unrelated healthy individuals from the Inner Mongolia Autonomous Region of China. A total of thirteen different HLA-A alleles, 21 different HLA-B alleles and 13 different HLA-DRB1 alleles were detected in the Oroqen ethnic group and the most frequent HLA alleles found were A*24(35.65%), B*15(17.92%), and DRB1*09(17.59%), respectively. The common HLA-A-B-DRB1 haplotypes were A*24-B*40-DRB1*09(5.09%), A*24-B*48-DRB1*12(2.78%) and A*24-B*51-DRB1*04(2.78%); and the HLA-A*33-B*58, A*30-B*13, A*01-B*37, A*33-DRB1*03, A*01-DRB1*10, A*30-DRB1*07, B*37-DRB1*10, B*58-DRB1*03, B*38-DRB1*08, B*13-DRB1*07 were significant positive linkage disequilibrium in the Oroqen nationality group. In total, 14 different HLA-A alleles, 27 B alleles and 12 DRB1 alleles were found in Ewenki nationality group, and the most frequent HLA alleles found were A*24(24.49%), B*40(17.35%), and DRB1*04(14.80%), respectively. The common HLA-A-B-DRB1 haplotypes were A*33-B*58-DRB1*03(6.25%), A*01-B*51-DRB1*11(2.88%) and A*24-B*40-DRB1*09(2.88%); the HLA-A*33-B*58, A*29-B*44, A*03-B*52, A*33-DRB1*03, A*29-DRB1*07, A*24-DRB1*09, B*58-DRB1*03, B*08-DRB1*03, B*46-DRB1*09 were significant positive linkage disequilibrium in Ewenki nationality group. The distribution of HLA A,-B, DRB1, alleles haplotypes frequencies and phylogenetic tree indicated that the Oroqen and Ewenki population groups belongs to northern group of China, together as a group cluster.