Disclosing the genetic structure of Brazil through analysis of male lineages with highly discriminating haplotypes

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

Association of apoptosis-related variants to malaria infection and parasite density in individuals from the Brazilian Amazon

BACKGROUND

In malaria infection, apoptosis acts as an important immunomodulatory mechanism that leads to the elimination of parasitized cells, thus reducing the parasite density and controlling immune cell populations. Here, it was investigated the association of INDEL variants in apoptotic genes-rs10562972 (FAS), rs4197 (FADD), rs3834129 and rs59308963 (CASP8), rs61079693 (CASP9), rs4647655 (CASP3), rs11269260 (BCL-2), and rs17880560 (TP53)-and the influence of genetic ancestry with susceptibility to malaria and parasite density in an admixed population from the Brazilian Amazon.

METHODS

Total DNA was extracted from 126 malaria patients and 101 uninfected individuals for investigation of genetic ancestries and genotypic distribution of apoptosis-related variants by Multiplex PCR. Association analyses consisted of multivariate logistic regressions, considering the following comparisons: (i) DEL/DEL genotype vs. INS/DEL + INS/INS; and (ii) INS/INS vs. INS/DEL + DEL/DEL.

RESULTS

Individuals infected by Plasmodium falciparum had significantly higher African ancestry proportions in comparison to uninfected controls, Plasmodium vivax, and mixed infections. The INS/INS genotype of rs3834129 (CASP8) seemed to increase the risk for P. falciparum infection (P = 0.038; OR = 1.867; 95% CI 0.736-3.725), while the DEL/DEL genotype presented a significant protective effect against infection by P. falciparum (P = 0.049; OR = 0.446; 95% CI 0.185-0.944) and mixed infection (P = 0.026; OR = 0.545; 95% CI 0.281-0.996), and was associated with lower parasite density in P. falciparum malaria (P = 0.009; OR = 0.383; 95% CI 0.113-1.295). Additionally, the INS/INS genotype of rs10562972 (FAS) was more frequent among individuals infected with P. vivax compared to P. falciparum (P = 0.036; OR = 2.493; 95% CI 1.104-4.551), and the DEL/DEL genotype of rs17880560 (TP53) was significantly more present in patients with mono-infection by P. vivax than in individuals with mixed infection (P = 0.029; OR = 0.667; 95% CI 0.211-1.669).

CONCLUSIONS

In conclusion, variants in apoptosis genes are associated with malaria susceptibility and parasite density, indicating the role of apoptosis-related genetic profiles in immune responses against malaria infection.

Allelic frequencies distribution and forensic parameters of 23 autosomal short tandem repeats in the population of the State of Pernambuco, Brazil

Located in the Northeast Region, the Pernambuco State is one of the 27 federative units of Brazil. Here, we determined populational data for 23 short tandem repeat (STR) markers - CSF1PO, FGA, TH01, TPOX, vWA, D1S1656, D2S1338, D2S441, D3S1358, D5S818, D7S820, D8S1179, D10S1248, D12S391, D13S317, D16S539, D18S51, D19S433, D21S11, D22S1045, PENTA D, PENTA E and SE33 - of the Pernambuco population. The sample consisted of 767 healthy, adult, unrelated individuals (437 males, 330 females) born and resident in the State of Pernambuco. STRs were amplified using three multiplex kits, according to the availability: PowerPlex® Fusion 6C System (Promega Corporation), PowerPlex® Fusion System (Promega Corporation) and GlobalFiler™ Express (Thermo Fisher Scientific). Allelic frequencies, forensic parameters and Hardy-Weinberg equilibrium determinations were estimated for all the 23 loci. No deviations from the Hardy-Weinberg equilibrium were observed for any of the markers, after Bonferroni correction. We observed that the most and less informative markers were SE33 and TPOX, respectively. The combined power of discrimination (CPD) was 0.99999999999999999999999999999, and the combined power of exclusion (CPE) was 0.99999999997. The cumulative typical paternity index was 37,919,301,869.3021. Interpopulation analyses (Nei's genetic distance) based on the expanded CODIS core loci was performed between the Pernambuco sample and other global populations. Pernambuco was the closest Brazilian population to African group and stayed distant from the Native American group. This work contributed to show that a panel of 23 autosomal STR loci is very informative, being able for forensic applications related in this population.

The Search for Cancer Biomarkers: Assessing the Distribution of INDEL Markers in Different Genetic Ancestries

Cancer is a multifactorial group of diseases, being highly incident and one of the leading causes of death worldwide. In Brazil, there is a great variation in cancer incidence and impact among the different geographic regions, partly due to the genetic heterogeneity of the population in this country, composed mainly by European (EUR), Native American (NAM), African (AFR), and Asian (ASN) ancestries. Among different populations, genetic markers commonly present diverse allelic frequencies, but in admixed populations, such as the Brazilian population, data is still limited, which is an issue that might influence cancer incidence. Therefore, we analyzed the allelic and genotypic distribution of 12 INDEL polymorphisms of interest in populations from the five Brazilian geographic regions and in populations representing EUR, NAM, AFR, and ASN, as well as tissue expression in silico. Genotypes were obtained by multiplex PCR and the statistical analyses were done using R, while data of tissue expression for each marker was extracted from GTEx portal. We highlight that all analyzed markers presented statistical differences in at least one of the population comparisons, and that we found 39 tissues to be differentially expressed depending on the genotype. Here, we point out the differences in genotype distribution and gene expression of potential biomarkers for risk of cancer development and we reinforce the importance of this type of study in populations with different genetic backgrounds.

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.

Mutational data and population profiling of 23 Y-STRs in three Brazilian populations

Y-chromosomal STRs are important markers in forensic genetics, due to some peculiar characteristics. The absence of recombination makes them a useful tool to infer kinship in complex cases involving distant paternal relatives, or to infer paternal bio-geographic ancestry. The presence of a single copy, being transmitted from father to son, allow tracing mutational events in Y-STRs without ambiguity. For the statistical interpretation of forensic evidences based on Y-STR profiles, it is necessary to have estimates on both mutation rates and haplotype frequencies. In this work, 407 father-son duos from São Paulo and Rio de Janeiro states and 204 unrelated individuals from Manaus were analyzed. Haplotype frequencies and mutation rates for the Y-STRs from the PowerPlex Y23 commercial kit were estimated. Thirty-six mutations were observed in 15 of the 22 Y-STRs analyzed, for an average mutation rate of 3.84 × 10^-3 (95 % CI 2.69 × 10^-3 to 5.32 × 10^-3). All mutations in GAAA repeats occurred in alleles with 13 or more uninterrupted units. Mutations in GATA repeats were observed in alleles with 9-17 uninterrupted units. An analysis carried out in different father's age groups showed an increase of 2.48 times the mutation rate in the age group of 40-50 years, when compared to the 20-30 age group, in agreement with the described for autosomal STRs. A high haplotype diversity was found in the three Brazilian populations. Pairwise genetic distance analysis (F_ST) showed no significant differences between the three populations in this study, which were also close to populations with strong European influence. The highest distances among the Brazilian populations were with São Gabriel da Cachoeira, which has a high Native American ancestry.

New insights on intercontinental origins of paternal lineages in Northeast Brazil

Background

The current Brazilian population is the product of centuries of admixture between intercontinental founding groups. Although previous results have revealed a heterogeneous distribution of mitochondrial lineages in the Northeast region, the most targeted by foreign settlers during the sixteenth century, little is known about the paternal ancestry of this particular population. Considering historical records have documented a series of territorial invasions in the Northeast by various European populations, we aimed to characterize the male lineages found in Brazilian individuals in order to discover to what extent these migrations have influenced the present-day gene pool. Our approach consisted of employing four hierarchical multiplex assays for the investigation of 45 unique event polymorphisms in the non-recombining portion of the Y-chromosome of 280 unrelated men from several Northeast Brazilian states.

Results

Primary multiplex results allowed the identification of six major haplogroups, four of which were screened for downstream SNPs and enabled the observation of 19 additional lineages. Results reveal a majority of Western European haplogroups, among which R1b-S116* was the most common (63.9%), corroborating historical records of colonizations by Iberian populations. Nonetheless, F_ST genetic distances show similarities between Northeast Brazil and several other European populations, indicating multiple origins of settlers. Regarding Native American ancestry, our findings confirm a strong sexual bias against such haplogroups, which represented only 2.5% of individuals, highly contrasting previous results for maternal lineages. Furthermore, we document the presence of several Middle Eastern and African haplogroups, supporting a complex historical formation of this population and highlighting its uniqueness among other Brazilian regions.

Conclusions

We performed a comprehensive analysis of the major Y-chromosome lineages that form the most dynamic migratory region from the Brazilian colonial period. This evidence suggests that the ongoing entry of European, Middle Eastern, and African males in the Brazilian Northeast, since at least 500 years, was significantly responsible for the present-day genetic architecture of this population.

Demand for larger Y-STR reference databases in ethnic melting-pot countries: Argentina as a test case

The Y chromosome behaves as a single locus. Its genetic information is useful in forensic casework, deficiency kinship testing, and population genetics studies. Continuous increases of loci number within commercial kits forced modification of worldwide reference databases. In Pan American countries, like Argentina, diverse parental ethnic groups contributed to the extant admixed urban populations. We report 509 additional haplotypes of 23 Y-STRs from donors inhabiting urban areas of six Argentinean provinces: Buenos Aires, Santiago del Estero, Santa Cruz, Rio Negro, Santa Fe, and Formosa. To better understand the demographic landscape of the admixed urban paternal lineages, structural analysis was performed using published data from other Argentinean provinces. AMOVA by Rst distance and inferred haplogroups by two predictive online software methods based on haplotypes yielded complementary results with respect to detected population structure, probably due to the different proportions of the Native American Q3-M3 haplogroup in the studied samples. This situation, which is common to most North, Meso, and South American countries, underscores the need for the additional step of typing specific SNPs for haplogroup diagnosis. We propose organizing Y-STR haplotype reference databases according to the most frequent haplogroups detected in a given admixed population.

mtDNA structure: the women who formed the Brazilian Northeast

BACKGROUND

The distribution of mitochondrial DNA (mtDNA) lineages in Brazil is heterogeneous due to different regional colonization dynamics. Northeastern Brazil, although being an important region in terms of human imigration and ethnic admixture, has little information regarding its population mtDNA composition. Here, we determine which mitochondrial lineages contributed to the formation of the Northeastern Brazilian population. Our sample consisted of 767 individuals distributed as follows i) 550 individuals from eight Northeastern states (Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas, Sergipe, and Bahia) which were sequenced for mtDNA hypervariable segments I, II, and III; ii) 217 individuals from Alagoas and Pernambuco (previously published data). Data analysis was performed through sequence alignment and Haplogrep 2.0 haplogroup assignment tools. Furthermore, maternal ancestry distribution was contextualized and, when possible, related to historical events to better understand the biological interactions and population dynamics that occurred in this region since the beginning of colonization.

RESULTS

Unexpectedly, Amerindian mitochondrial ancestry was the highest in the Northeastern region overall, followed by African, European and non-Amerindian Asian, unlike previous results for this region. Alagoas and Pernambuco states, however, showed a larger African mtDNA frequency. The Northeastern region showed an intraregional heterogeneous distribution regarding ancestral groups, in which states/mesoregions located to the north had a prevalent Amerindian ancestral frequency and those to the south had predominance of African ancestry. Moreover, results showed great diversity of European haplogroups and the presence of non-Amerindian Asian haplogroups.

CONCLUSIONS

Our findings are in disagreement with previous investigations that suggest African mitochondrial ancestry is the most prevalent in the Brazilian Northeast. The predominance of Amerindian lineages exemplifies the importance of indigenous women in the formation of the population, despite intense African slave entry and conflicts with European settlers. The variable distribution of ancestral groups observed in the Northeast is in accordance with historical records showing the similarities with colonization dynamics occurred in the Amazon region and the Brazilian Southeast. Moreover, the variety of European haplogroups suggests multiple origins of founding groups, specially those found in Western European populations.

HABP2 p.G534E variant in patients with family history of thyroid and breast cancer

Familial Papillary Thyroid Carcinoma (PTC) has been described as a hereditary predisposition cancer syndrome associated with mutations in candidate genes including HABP2. Two of 20 probands from families with history of PTC and breast carcinoma (BC) were evaluated by whole exome sequencing (WES) revealing HABP2 p.G534E. Sanger sequencing was used to confirm the involvement of this variant in three families (F1: 7 relatives; F2: 3 and F3: 3). The proband and his sister (with no malignant tumor so far) from F1 were homozygous for the variant whereas one relative with PTC from F2 was negative for the variant. Although the proband of the F3 with PTC was HABP2 wild type, three relatives presented the variant. Five of 170 healthy Brazilian individuals with no family history of BC or PTC and three of 50 sporadic PTC presented the p.G534E. These findings suggested no association of this variant with our familial PTC cases. Genes potentially associated with deregulation of the extracellular matrix organization pathway (CTSB, TNXB, COL4A3, COL16A1, COL24A1, COL5A2, NID1, LOXL2, MMP11, TRIM24 and MUSK) and DNA repair function (NBN and MSH2) were detected by WES, suggesting that other cancer-associated genes have pathogenic effects in the risk of familial PTC development.

Analysis of the genetic ancestry of patients with oral clefts from South American admixed populations

Increased susceptibility to cleft lip, with or without cleft palate (CL±P) has been observed in South America, as related to Amerindian ancestry, using epidemiological data, uniparental markers, and blood groups. In this study, it was evaluated whether this increased risk remains when Amerindian ancestry is estimated using autosomal markers and considered in the predictive model. Ancestry was estimated through genotyping 62 insertion and deletion (INDEL) markers in sample sets of patients with CL±P, patients with cleft palate (CP), and controls, from Patagonia in southern Argentina and Belém in northern Brazil. The Amerindian ancestry in patients from Patagonia with CL±P was greater than in controls although it did not reach statistical significance. The European ancestry in patients with CL±P from Belém and in patients with CP from Belém and Patagonia was higher than in controls and statistically significant for patients with CP who were from Belém. This high contribution of European genetic ancestry among patients with CP who were from Belém has not been previously observed in American populations. Our results do not corroborate the currently accepted risks for CL±P and CP estimated by epidemiological studies in the North American populations and probably reflect the higher admixture found in South American ethnic groups when compared with the same ethnic groups from the North American populations.

Male Lineages in Brazil: Intercontinental Admixture and Stratification of the European Background

The non-recombining nature of the Y chromosome and the well-established phylogeny of Y-specific Single Nucleotide Polymorphisms (Y-SNPs) make them useful for defining haplogroups with high geographical specificity; therefore, they are more apt than the Y-STRs to detect population stratification in admixed populations from diverse continental origins. Different Y-SNP typing strategies have been described to address issues of population history and movements within geographic territories of interest. In this study, we investigated a set of 41 Y-SNPs in 1217 unrelated males from the five Brazilian geopolitical regions, aiming to disclose the genetic structure of male lineages in the country. A population comparison based on pairwise FST genetic distances did not reveal statistically significant differences in haplogroup frequency distributions among populations from the different regions. The genetic differences observed among regions were, however, consistent with the colonization history of the country. The sample from the Northern region presented the highest Native American ancestry (8.4%), whereas the more pronounced African contribution could be observed in the Northeastern population (15.1%). The Central-Western and Southern samples showed the higher European contributions (95.7% and 93.6%, respectively). The Southeastern region presented significant European (86.1%) and African (12.0%) contributions. The subtyping of the most frequent European lineage in Brazil (R1b1a-M269) allowed differences in the genetic European background of the five Brazilian regions to be investigated for the first time.

Distribution of allelic and genotypic frequencies of IL1A, IL4, NFKB1 and PAR1 variants in Native American, African, European and Brazilian populations

BACKGROUND

The inflammatory response plays a key role at different stages of cancer development. Allelic variants of the interleukin 1A (IL1A), interleukin 4 (IL4), nuclear factor kappa B1 (NFKB1) and protease-activated receptor 1 (PAR1) genes may influence not only the inflammatory response but also susceptibility to cancer development. Among major ethnic or continental groups, these polymorphic variants present different allelic frequencies. In admixed populations, such as the Brazilian population, data on distribution of these polymorphisms are limited. Here, we collected samples of cancer-free individuals from the north, northeast, midwest, south and southeast regions of Brazil and from the three main groups that gave rise to the Brazilian population: Native Americans from the Brazilian Amazon, Africans and Europeans. We describe the allelic distributions of four IL1A (rs3783553), IL4 (rs79071878), NFKB1 (rs28362491) and PAR1 (rs11267092) gene polymorphisms, which the literature describes as polymorphisms with a risk of cancer or worse prognosis for cancer.

RESULTS

The genotypic distribution of the four polymorphisms was statistically distinct between Native Americans, Africans and Europeans. For the allelic frequency of these polymorphisms, the Native American population was the most distinct among the three parental populations, and it included the greatest number of alleles with a risk of cancer or worse prognosis for cancer. The PAR1 gene polymorphism allelic distribution was similar among all Brazilian regions. For the other three markers, the northern region population was statistically distinct from other Brazilian region populations.

CONCLUSION

The IL1A, IL4, NFKB1 and PAR1 gene polymorphism allelic distributions are homogeneous among the regional Brazilian populations, except for the northern region, which significantly differs from the other four Brazilian regions. Among the parental populations, the Native American population exhibited a higher incidence of alleles with risk of cancer or worse prognosis for cancer, which can indicate greater susceptibility to this disease. These genetic data may be useful for future studies on the association between these polymorphisms and cancer in the investigated populations.

Influence of Genetic Ancestry on INDEL Markers of NFKβ1, CASP8, PAR1, IL4 and CYP19A1 Genes in Leprosy Patients

BACKGROUND

Leprosy is an insidious infectious disease caused by the obligate intracellular bacteria Mycobacterium leprae, and host genetic factors can modulate the immune response and generate distinct categories of leprosy susceptibility that are also influenced by genetic ancestry.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the possible effects of CYP19A1 [rs11575899], NFKβ1 [rs28362491], IL1α [rs3783553], CASP8 [rs3834129], UGT1A1 [rs8175347], PAR1 [rs11267092], CYP2E1 [INDEL 96pb] and IL4 [rs79071878] genes in a group of 141 leprosy patients and 180 healthy individuals. The INDELs were typed by PCR Multiplex in ABI PRISM 3130 and analyzed with GeneMapper ID v3.2. The NFKβ1, CASP8, PAR1 and IL4 INDELs were associated with leprosy susceptibility, while NFKβ1, CASP8, PAR1 and CYP19A1 were associated with the MB (Multibacilary) clinical form of leprosy.

CONCLUSIONS/SIGNIFICANCE

NFKβ1 [rs28362491], CASP8 [rs3834129], PAR1 [rs11267092] and IL4 [rs79071878] genes are potential markers for susceptibility to leprosy development, while the INDELs in NFKβ1, CASP8, PAR1 and CYP19A1 (rs11575899) are potential markers for the severe clinical form MB. Moreover, all of these markers are influenced by genetic ancestry, and European contribution increases the risk to leprosy development, in other hand an increase in African contribution generates protection against leprosy.

Male lineage strata of Brazilian population disclosed by the simultaneous analysis of STRs and SNPs

Brazil has a large territory divided in five geographical regions harboring highly diverse populations that resulted from different degrees and modes of admixture between Native Americans, Europeans and Africans. In this study, a sample of 605 unrelated males was genotyped for 17 Y-STRs and 46 Y-SNPs aiming a deep characterization of the male gene pool of Rio de Janeiro and its comparison with other Brazilian populations. High values of Y-STR haplotype diversity (0.9999±0.0001) and Y-SNP haplogroup diversity (0.7589±0.0171) were observed. Population comparisons at both haplotype and haplogroup levels showed significant differences between Brazilian South Eastern and Northern populations that can be explained by differences in the proportion of African and Native American Y chromosomes. Statistical significant differences between admixed urban samples from the five regions of Brazil were not previously detected at haplotype level based on smaller size samples from South East, which emphasizes the importance of sample size to detected population stratification for an accurate interpretation of profile matches in kinship and forensic casework. Although not having an intra-population discrimination power as high as the Y-STRs, the Y-SNPs are more powerful to disclose differences in admixed populations. In this study, the combined analysis of these two types of markers proved to be a good strategy to predict population sub-structure, which should be taken into account when delineating forensic database strategies for Y chromosome haplotypes.

Revisiting the genetic ancestry of Brazilians using autosomal AIM-Indels

There are many different studies that contribute to the global picture of the ethnic heterogeneity in Brazilian populations. These studies use different types of genetic markers and are focused on the comparison of populations at different levels. In some of them, each geographical region is treated as a single homogeneous population, whereas other studies create different subdivisions: political (e.g., pooling populations by State), demographic (e.g., urban and rural), or ethnic (e.g., culture, self-declaration, or skin colour). In this study, we performed an enhanced reassessment of the genetic ancestry of ~ 1,300 Brazilians characterised for 46 autosomal Ancestry Informative Markers (AIMs). In addition, 798 individuals from twelve Brazilian populations representing the five geographical macro-regions of Brazil were newly genotyped, including a Native American community and a rural Amazonian community. Following an increasing North to South gradient, European ancestry was the most prevalent in all urban populations (with values up to 74%). The populations in the North consisted of a significant proportion of Native American ancestry that was about two times higher than the African contribution. Conversely, in the Northeast, Center-West and Southeast, African ancestry was the second most prevalent. At an intrapopulation level, all urban populations were highly admixed, and most of the variation in ancestry proportions was observed between individuals within each population rather than among population. Nevertheless, individuals with a high proportion of Native American ancestry are only found in the samples from Terena and Santa Isabel. Our results allowed us to further refine the genetic landscape of Brazilians while establishing the basis for the effective application of an autosomal AIM panel in forensic casework and clinical association studies within the highly admixed Brazilian populations.