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.

Introduction

Tuberculosis (TB) is an infectious disease caused by the Mycobacterium tuberculosis bacillus (MTB). In 2018, 10 million incident cases of TB were registered, resulting in 1.2 million TB deaths, a number higher than that caused by Human Immunodeficiency Virus [1,2]. While the incidence of TB and the total number of deaths, worldwide, have declined slightly in recent years, TB continues to be a major endemic concern for the amerindian people of the Brazilian Amazon [3,4].

In amerindian population tend to have much higher rates of TB in comparison with the general (non-indigenous) population [5-7]. This disparity has been attributed to multiple factors, including environmental and socioeconomic variables, co-infections with other pathogens, and the genetics of the host [8]. The profile of genomic ancestry is known to be a determining factor in the susceptibility of some groups to certain diseases, as well as their response to treatment [9].

Brazil is one of the most heterogeneous countries in the world with a population formed primarily by the admixture of three continental groups: European (EUR), African (AFR), and Amerindian, AMR [10]. Overall, EUR ancestry predominates in Brazil (62.4%), while AFR ancestry is the second most prevalent (22.6%), and AMR, the least common, at 14.7% [11]. In northern Brazil, however, the two minor ancestries are reversed, with AMR increasing to 28%, while EUR decreases to 60% and AFR to 12% [12-15].

Genetic studies have indicated that populations of different ancestry, in particular AMR, may be relatively more susceptible to certain diseases than other populations [15-18]. In Mexican and Peruvian populations, for example, AMR ancestry was associated with increased risk for systemic lupus erythematosus [16]. In northern Brazil, AMR ancestry has been associated with susceptibility to the development of leprosy, including the severe clinical form [17], as well as an increased prevalence of hepatitis B and C [18], and an increased risk of developing acute lymphoblastic leukemia [15].

While these studies have focused on different diseases and populations from northern Brazil and other regions of Latin America, no data are available on the incidence of TB or its association with genetic factors in the population of the Amazon region. Given this, the present study investigated the potential role of AMR genetic ancestry in the susceptibility to TB infection of a highly admixed Amazonian population.

Material and methods

Study populations

This is a case-control study was conducted between March 2008 and March 2010 at João de Barros Barreto Hospital, in Belém city, Brazil. Two hundred eighty TB patients were enrolled in case group and one hundred and thirty-eight in control group. The case group was composed of participants aged over 18 years with confirmed pulmonary TB according to the World Health Organization (WHO). Exclusion criteria was individuals with confirmed extra-pulmonary TB. The control group was composed of participants who were employees from the Hospital in contact with active TB patients for at least five years with neither symptoms nor diagnostic confirmed of TB.

Ethical aspects

The protocol adopted in the present study was approved by the Committee for Research Ethics of the Federal University of Pará (approval no. 350507). Informed consent was obtained from all the participants of the study.

Extraction and quantification of the DNA

The genomic DNA was isolated from peripheral blood leukocytes using the Biopur mini spin plus kit (Biometrix Inc., San Francisco, CA, USA) and quantified with a Thermo Scientific NanoDrop 1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE).

Ancestry analysis

A set of 61 Ancestry-Informative Markers (AIMs) was used to estimate the proportion of three different ancestries: European (EUR), African (AFR), and Amerindian (AMR). These indel-type markers were previously described [12,19]. The AIMs were genotyped by Multiplex PCR, followed by capillary electrophoresis with fragment analysis. The multiplex amplifications were run using the QIAGEN Multiplex PCR kit (QIAGEN, Germany), with the PCR being run in an ABI Veriti thermocycler (Life Technologies, USA), followed by the capillary electrophoresis protocol. The DNA fragments were separated using an ABI PRISM 3130 Genetic Analyzer and peak reads were obtained in GeneMapper ID v3.2 software (Life Technologies). The experiment was replicated twice.

Statistical analysis

The samples were analyzed in STRUCTURE v2.2 to determine the relative contribution of EUR, AFR and AMR ancestries to the genetic profile of each individual [20,21]. Additional statistical analyses were made using the software RStudio v.3.6.3 with the basic statistics package (stats) for the R language. The analysis of a categorical non-parametric clinical data (sex) was performed with the Fisher test, while the analysis of an ordinal non-parametric clinical data (age) was performed with Mann-Whitney test. Ancestry informative markers (AIM) was analyzed using the multivariate logistic regression model. Clinical variables with significant results were used as covariates. All tests were two-tailed and p-value ≤ 0.05 was considered statistically significant.

Results

Two hundred and eighty tuberculosis patients and one hundred thirty-eight control individuals were analyzed. Table 1 presents the demographic and clinical characteristics of the case and control group. The results showed significant differences when the cases were compared against controls in relation to sex and African, European and Amerindian ancestry. The female gender was predominant in the case and control group with 52% and 85% respectively (p = 1.33x10-11).

Table 1

Demographic and clinical characteristics of the case and control groups selected for the present study.

VARIABLE	TUBERCULOSIS (n = 280)	CONTROL (n = 138)	p-value
Age (years)	47 ± 21.1	52 ± 14.4	0.052 a
Sex (male/female)	134 (48%) / 146 (52%)	21 (15%) / 117 (85%)	1.33x10-11 b
AFR ancestry	0.254 ± 0.12	0.216 ± 0.11	0.004 c
EUR ancestry	0.409 ± 0.13	0.499 ± 0.17	7.56x10-7 c
AMR ancestry	0.337 ± 0.13	0.285 ± 0.14	0.001 c

aMann-Whitney test

bFisher's exact test

cMultivariate regression logistic with sex as a covariate. The data are presented as the median ± standard deviation.

The ethnic composition of the case group was 25.4%% African, 40.9% European, and 33.7% Amerindian, while in the control group composition was 21.6% African, 49.9% European, 28.5% Amerindian (Table 1). These results suggested an increased contribution from Amerindian ancestry (p = 0.001) and a loss of contribution from European ancestry (p <0.001) in the case group compared to the control group.

Distribution of the three principal ancestries, EUR, AFR, and AMR, in the case and control groups was presented in Fig 1. The ancestry of the case-control groups was investigated by comparing interdependent components. These results indicated that EUR ancestry was less amply distributed in the TB group, whereas the AMR and AFR ancestries of this group were higher than those of the control group (Fig 1).

Fig 1

Distribution of the three principal ancestries (EUR, AFR, and AMR) in the case (TB) and control (non-TB) groups analyzed in the present study.

Significance determined by multivariate regression logistic with sex as a covariate.

The individual interethnic profiles of the study participants were presented in Fig 2. In this plot, the TB patients were represented by orange dots and the control subjects by green dots, and their respective positions within the vertices of the triangle correspond to the relative contribution of EUR, AFR and AMR ancestry. This analysis provided a visual representation of the different ancestries of the subjects of the two groups. Most of the TB patients were closer to the lower part of the triangle, which reflects their greater average AMR ancestry. By contrast, most individuals in the control group were closer to the upper portion of the triangle, reflecting their greater EUR ancestry.

Fig 2

Estimates of genetic admixture in the case and control individuals analyzed in the present study.

The influence of the range of AMR ancestry on the susceptibility of individuals to TB was assessed using multivariate regression logistic model with sex as a covariate (Fig 3). This analysis indicated that AMR ancestry in the 20% to 60% range was a risk factor for susceptibility to TB, with the risk peaking in the upper portion of this range (40–50%). When the contribution of AMR ancestry was in the 20–30% range, the estimated risk effect was 1.9 (OR = 1.92, 95%CI = 1.03–3.63; p = 0.041), in the 30–40% was 1.92 (OR = 1.91, 95% CI = 1.01–3.66, p = 0.047). Risk peaked in the 40–50% range, where the Odds Ratio for TB infection was 3.25 (95%CI = 1.56–7.01, p = 0.002) and decreased in the 50–60% range where the OR was 2.72 (95% CI = 1.09–7.30, p = 0.037). The lack of any significant risk effect at the extremes of this range (below 20% AMR ancestry and above 60%) was probably due to the restrictively small sample sizes for these levels of ancestry, which limited the discriminative power of the statistical tests.

Fig 3

Variation in the odds ratios recorded for different percentages of Amerindian ancestry as a risk factor for tuberculosis.

Significance as determined by multivariate regression logistic with sex as a covariate.

Discussion

In the present study, we investigated the potential role of Amerindian genomic ancestry in the susceptibility to tuberculosis of individuals from the Brazilian Amazon region. This population is admixed from three groups, of European (EUR), African (AFR), and Amerindian (AMR) ancestry. Brazil’s interethnic mix is ​​one of the most heterogeneous anywhere in the world, resulting in a uniquely diverse population [22]. The population of Brazil is ethnically distinct from the populations of both neighboring countries and other continents [23].

Tuberculosis still occurs at high rates and is one of the leading causes of disease-related mortality, worldwide [24]. The northern (Amazon) region of Brazil has relatively high rates of TB in comparison with the other regions of the country [25]. That is probably because of factors such as human polygenic variability which is an important determinant of the outcome of infection with MTB [26]. Epidemiological studies have evaluated high incidence of TB disease in Amerindian [27,28] and North of Brazil has high contribution of AMR ancestry [13] so these factors may be associated and consequently have important clinical implications.

The present study evaluated the role of genetic ancestry in the susceptibility to TB in a highly miscegenated population. The results of the genetic analysis of the three ancestral groups, i.e., European (EUR), African (AFR), and Amerindian (AMR), found significant differences among these groups. In particular, the AMR ancestry of the TB group was significantly higher than that of the control group and increasing AMR ancestry also contributed to an increased predisposition for infection by TB.

The highly heterogenous population of Brazil may be relatively more susceptible to some diseases than other populations, and the results of the present study further reinforce the need for the identification of the EUR, AFR and AMR ancestries in a given population in order to best predict the health risks of each individual [29]. Norris et al. [30] reported significant variation in the pathways of immunity, metabolism, and disease in distinct populations. These factors may contribute to the health-related disparities found in substructured Latin American populations, in particular, that different ancestral profiles may have distinct predispositions to specific diseases.

Chimusa et al. [31] reported an association of ancestry with TB susceptibility in an African population and concluded that this ancestry conferred a high risk of contracting the disease. Daya et al. [32] also concluded that increased AFR ancestry and decreased EUR ancestry were associated with an increased risk of developing TB.

A number of studies in Brazil have also highlighted the importance of the role of ancestry in the susceptibility of a population to disease, including the study of Pinto et al. [17], which identified a set of potential genetic markers related to susceptibility to leprosy, influenced in particular by AFR and EUR genomic ancestry. Carvalho et al. [33] also found an association between AMR ancestry and an increased risk of developing acute lymphoblastic leukemia in the population of Brazilian Amazonia. Lopes et al. [34] also concluded that reduced EUR genetic ancestry and increased AFR ancestry had a significant effect on susceptibility to chronic periodontitis.

Amerindian populations are immunologically distinct from non-Amerindian ones, with reduced variability in some immune system genes such as HLA, KM, GM, and KIR [35-39]. This implies that Amerindian populations may have a distinct immune response from the general population, in particular, that the Th2-type immune response is predominant in this group [27]. The Th2-type response induces a humoral response, which is inefficient in the defense of the organism against intracellular pathogens, such as MTB [40,41].

The results of the present study show that Amerindian genomic ancestry is a risk factor for tuberculosis, although we believe that this influence does not arise from the neutral polymorphisms used to infer genomic ancestry, but rather, reflects the low genetic diversity of the immune response in Amerindian populations and, in turn, of the admixed population with a high AMR ancestry. This means that individuals who have a higher Amerindian genomic ancestry have a higher risk of developing tuberculosis.

The incidence of TB in Amerindian populations has been amply studied from both epidemiological and socioeconomic viewpoints [5,28,42,43], albeit with a deficit in genetic research. Given this, the present study was a pioneering investigation of the association of Amerindian genetic ancestry with TB in the population of Brazilian Amazonia, which highlighted the importance of disease-oriented genomic ancestry studies in admixed populations. Our findings indicate that Amerindian ancestry is an important risk factor for TB in the admixed populations of the Brazilian Amazon region. Further research into the relationship between genetic ancestry and the prevalence of infectious diseases may provide important practical insights for the diagnosis, prevention, and treatment of diseases in different populations.

There were some limitations in the study. First, the sample sizes were relatively small and the case and control groups were individuals with high admixture of ancestry and genetic background is different from other populations, so findings might not generalizable to the population of Brazil. Second, the study was unable to obtain sufficient data related to socio-economic factor, cigarette smoking, alcohol consumption, in order to perform a stratified analysis. Despite these drawbacks, it could detect association between genetic ancestry and tuberculosis. There is lack in the literature about the influence of Amerindian ancestry in the risk or protection on infectious disease as tuberculosis. In future work, there is an interest in whole exome-sequencing in admixture population with high proportion of Amerindian ancestry to investigate immune gene variants associated to susceptibility to pulmonary tuberculosis infection.

Distribution table of patients according to the Amerindian ancestry range, with p-value, odds ratio and 95% confidence interval.

All statistical data were performed using the multivariate logistic regression model with the sex control variable.

(DOCX)

Click here for additional data file.

Tabela de Dados Brutos.

(PDF)

Click here for additional data file.