Population-based change-point detection for the identification of homozygosity islands.

MOTIVATION: This work is motivated by the problem of identifying homozygosity islands on the genome of individuals in a population. Our method directly tackles the issue of identification of the homozygosity islands at the population level, without the need of analysing single individuals and then combine the results, as is made nowadays in state-of-the-art approaches. RESULTS: We propose regularized offline change-point methods to detect changes in the parameters of a multidimensional distribution when we have several aligned, independent samples of fixed resolution. We present a penalized maximum likelihood approach that can be efficiently computed by a dynamic programming algorithm or approximated by a fast binary segmentation algorithm. Both estimators are shown to converge almost surely to the set of change-points without the need of specifying a priori the number of change-points. In simulation, we observed similar performances from the exact and greedy estimators. Moreover, we provide a new methodology for the selection of the regularization constant which has the advantage of being automatic, consistent, and less prone to subjective analysis. AVAILABILITY AND IMPLEMENTATION: The data used in the application are from the Human Genome Diversity Project (HGDP) and is publicly available. Algorithms were implemented using the R software R Core Team (R: A Language and Environment for Statistical Computing. Vienna (Austria): R Foundation for Statistical Computing, 2020.) in the R package blockcpd, found at https://github.com/Lucas-Prates/blockcpd.

Fourteen-year trends in overweight, general obesity, and abdominal obesity in Amazonian indigenous peoples.

BACKGROUND: Available data show that the epidemiological profile of most indigenous Brazilian populations is characterized by the coexistence of long-standing health problems (high prevalence of infectious and parasitic diseases, malnutrition, and deficiency diseases, such as anemia in children and women of reproductive age), associated with new health problems, especially those related to obesity (hypertension, type 2 diabetes mellitus and dyslipidemia). Based on this scenario, this study analyzed the nutritional profile of the adult population of seven indigenous peoples from the Brazilian Amazon in the years 2007 and 2021. METHODS: A total of 598 adults individuals were analyzed in 2007 (319 women and 279 men) and 924 in 2021 (483 women and 441 men), from seven indigenous peoples located in the state of Pará, who were assisted during health actions carried out in 2007 and in 2021. Body mass index classification used the World Health Organization criteria for adults: low weight, < 18.5 kg/m2; normal weight, ≥ 18.5 and < 25 kg/m2); overweight, ≥ 25 and < 30 kg/m2, and obesity, ≥ 30 kg/m2. A waist circumference (WC) < 90 cm in men and < 80 cm in women was considered normal. RESULTS: The data revealed heterogeneous anthropometric profiles, with a low prevalence of nutritional changes in the Araweté, Arara and Parakanã peoples, and high proportions of excess weight and abdominal obesity in the Kararaô, Xikrin do Bacajá, Asurini do Xingu and Gavião peoples, similar to or even higher than the national averages. CONCLUSION: Different stages of nutritional transition were identified in the indigenous peoples analyzed, despite apparently having been subjected to the same environmental pressures that shaped their nutritional profile in recent decades, which may indicate different genetic susceptibilities to nutritional changes. The evidence shown in this study strongly suggests the need to investigate in greater depth the genetic and environmental factors associated with the nutritional profile of Brazilian indigenous peoples, with assessment of diet, physical activity and sociodemographic and socioeconomic variables that enable the development of appropriate prevention and monitoring measures.

Population Histories and Genomic Diversity of South American Natives.

South America is home to one of the most culturally diverse present-day native populations. However, the dispersion pattern, genetic substructure, and demographic complexity within South America are still poorly understood. Based on genome-wide data of 58 native populations, we provide a comprehensive scenario of South American indigenous groups considering the genomic, environmental, and linguistic data. Clear patterns of genetic structure were inferred among the South American natives, presenting at least four primary genetic clusters in the Amazonian and savanna regions and three clusters in the Andes and Pacific coast. We detected a cline of genetic variation along a west-east axis, contradicting a hard Andes-Amazon divide. This longitudinal genetic variation seemed to have been shaped by both serial population bottlenecks and isolation by distance. Results indicated that present-day South American substructures recapitulate ancient macroregional ancestries and western Amazonia groups show genetic evidence of cultural exchanges that led to language replacement in precontact times. Finally, demographic inferences pointed to a higher resilience of the western South American groups regarding population collapses caused by the European invasion and indicated precontact population reductions and demic expansions in South America.

Genomic insight into the origins and dispersal of the Brazilian coastal natives.

In the 15th century, ∼900,000 Native Americans, mostly Tupí speakers, lived on the Brazilian coast. By the end of the 18th century, the coastal native populations were declared extinct. The Tupí arrived on the east coast after leaving the Amazonian basin ∼2,000 y before present; however, there is no consensus on how this migration occurred: toward the northern Amazon and then directly to the Atlantic coast, or heading south into the continent and then migrating to the coast. Here we leveraged genomic data from one of the last remaining putative representatives of the Tupí coastal branch, a small, admixed, self-reported Tupiniquim community, as well as data of a Guaraní Mbyá native population from Southern Brazil and of three other native populations from the Amazonian region. We demonstrated that the Tupiniquim Native American ancestry is not related to any extant Brazilian Native American population already studied, and thus they could be considered the only living representatives of the extinct Tupí branch that used to settle the Atlantic Coast of Brazil. Furthermore, these data show evidence of a direct migration from Amazon to the Northeast Coast in pre-Columbian time, giving rise to the Tupí Coastal populations, and a single distinct migration southward that originated the Guaraní people from Brazil and Paraguay. This study elucidates the population dynamics and diversification of the Brazilian natives at a genomic level, which was made possible by recovering data from the Brazilian coastal population through the genomes of mestizo individuals.

Haplotypes of single cancer driver genes and their local ancestry in a highly admixed long-lived population of Northeast Brazil.

Admixed populations have not been examined in detail in cancer genetic studies. Here, we inferred the local ancestry of cancer-associated single nucleotide polymorphisms (SNPs) and haplotypes of a highly admixed Brazilian population. SNP array was used to genotype 73 unrelated individuals aged 80-102 years. Local ancestry inference was performed by merging genotyped regions with phase three data from the 1000 Genomes Project Consortium using RFmix. The average ancestry tract length was 9.12-81.71 megabases. Strong linkage disequilibrium was detected in 48 haplotypes containing 35 SNPs in 10 cancer driver genes. All together, 19 risk and eight protective alleles were identified in 23 out of 48 haplotypes. Homozygous individuals were mainly of European ancestry, whereas heterozygotes had at least one Native American and one African ancestry tract. Native-American ancestry for homozygous individuals with risk alleles for HNF1B, CDH1, and BRCA1 was inferred for the first time. Results indicated that analysis of SNP polymorphism in the present admixed population has a high potential to identify new ancestry-associated alleles and haplotypes that modify cancer susceptibility differentially in distinct human populations. Future case-control studies with populations with a complex history of admixture could help elucidate ancestry-associated biological differences in cancer incidence and therapeutic outcomes.

Common BMI and diabetes-related genetic variants: A pilot study among indigenous people in the Brazilian Amazon.

This study was carried out to investigate the frequency of genetic variants related to body mass index (BMI) and type 2 diabetes (T2D) and evaluating the potential impact of risk alleles on susceptibility to these disorders in six indigenous peoples from Brazilian Amazon region. The majority of Fst values for pairwise population comparisons among the indigenous groups are low or moderate. The indigenous people show high values of differentiation with Africans, Europeans and Southeast Asians and moderate values with East Asian and American populations, as expected. The allelic frequencies among indigenous indicate that the majority of associations observed with T2D in continental populations can be replicated in native Amazonians. The genetic risk scores calculated for T2D in indigenous are high and similar to those calculated for Americans and East Asians, while the estimates obtained for obesity are low, probably due to the low frequencies of the risk allele of the FTO gene found in our samples. ADRB3-rs4994 and ABCC8-rs1799854 genes showed a significant association with BMI and waist circumference, and the KCNJ11-rs5219 gene with hyperglycemia. These results emphasize the importance of knowing the genetic variability underlying complex genetic diseases in indigenous peoples and the search for particular or rare variants.

Hepatoblastomas exhibit marked NNMT downregulation driven by promoter DNA hypermethylation.

Hepatoblastomas exhibit the lowest mutational burden among pediatric tumors. We previously showed that epigenetic disruption is crucial for hepatoblastoma carcinogenesis. Our data revealed hypermethylation of nicotinamide N-methyltransferase, a highly expressed gene in adipocytes and hepatocytes. The expression pattern and the role of nicotinamide N-methyltransferase in pediatric liver tumors have not yet been explored, and this study aimed to evaluate the effect of nicotinamide N-methyltransferase hypermethylation in hepatoblastomas. We evaluated 45 hepatoblastomas and 26 non-tumoral liver samples. We examined in hepatoblastomas if the observed nicotinamide N-methyltransferase promoter hypermethylation could lead to dysregulation of expression by measuring mRNA and protein levels by real-time quantitative polymerase chain reaction, immunohistochemistry, and Western blot assays. The potential impact of nicotinamide N-methyltransferase changes was evaluated on the metabolic profile by high-resolution magic angle spinning nuclear magnetic resonance spectroscopy. Significant nicotinamide N-methyltransferase downregulation was revealed in hepatoblastomas, with two orders of magnitude lower nicotinamide N-methyltransferase expression in tumor samples and hepatoblastoma cell lines than in hepatocellular carcinoma cell lines. A specific TSS1500 CpG site (cg02094283) of nicotinamide N-methyltransferase was hypermethylated in tumors, with an inverse correlation between its methylation level and nicotinamide N-methyltransferase expression. A marked global reduction of the nicotinamide N-methyltransferase protein was validated in tumors, with strong correlation between gene and protein expression. Of note, higher nicotinamide N-methyltransferase expression was statistically associated with late hepatoblastoma diagnosis, a known clinical variable of worse prognosis. In addition, untargeted metabolomics analysis detected aberrant lipid metabolism in hepatoblastomas. Data presented here showed the first evidence that nicotinamide N-methyltransferase reduction occurs in hepatoblastomas, providing further support that the nicotinamide N-methyltransferase downregulation is a wide phenomenon in liver cancer. Furthermore, this study unraveled the role of DNA methylation in the regulation of nicotinamide N-methyltransferase expression in hepatoblastomas, in addition to evaluate the potential effect of nicotinamide N-methyltransferase reduction in the metabolism of these tumors. These preliminary findings also suggested that nicotinamide N-methyltransferase level may be a potential prognostic biomarker for hepatoblastoma.

HW_TEST, a program for comprehensive HARDY-WEINBERG equilibrium testing.

This article deals with a Windows (© Microsoft Inc.) executable, user-friendly program that tests the hypothesis of Hardy-Weinberg (HW) proportions from autosomal multiallelic data using different methods that include parametric, nonparametric and exact bootstrap tests, the latter obtained through computer simulations. The program can be obtained free of charge directly from the internet repository https://github.com/Lemes-RenanB/HardyWeinbergTesting.

Estimation of inbreeding and substructure levels in African-derived Brazilian quilombo populations.

This article deals with the estimation of inbreeding and substructure levels in a set of 10 (later regrouped as eight) African-derived quilombo communities from the Ribeira River Valley in the southern portion of the state of São Paulo, Brazil. Inbreeding levels were assessed through F-values estimated from the direct analysis of genealogical data and from the statistical analysis of a large set of 30 molecular markers. The levels of population substructure found were modest, as was the degree of inbreeding: in the set of all communities considered together, F-values were 0.00136 and 0.00248 when using raw and corrected data from their complete genealogical structures, respectively, and 0.022 and 0.036 when using the information taken from the statistical analysis of all 30 loci and of 14 single-nucleotide polymorphic loci, respectively. The overall frequency of consanguineous marriages in the set of all communities considered together was ∼ 2%. Although modest, the values of the estimated parameters are much larger than those obtained for the overall Brazilian population and in general much smaller than the ones recorded for other Brazilian isolates. To circumvent problems related to heterogeneous sampling and virtual absence of reliable records of biological relationships, we had to develop or adapt several methods for making valid estimates of the prescribed parameters.

The seroconversion history to SARS-CoV-2 in Indigenous people from Brazil - the interplay between exposure, vaccination, and tuberculosis.

The COVID-19 pandemic caused a significant loss of human lives and a worldwide decline in quality of life. Although our understanding of the pandemic has improved significantly since the beginning, the natural history of COVID-19 and its impacts on under-represented populations, such as Indigenous people from America, remain largely unknown. We performed a retrospective serological survey with two Brazilian Indigenous populations (n=624), Tupiniquim and Guarani-Mbyá. Samples were collected between September 2020 and July 2021: a period comprising the dissemination of SARS-CoV-2 variants and the beginning of COVID-19 vaccination in Brazil. Seroconversions against S and N antigens were assessed using three different commercially available ELISA kits. Samples were also used to assess the prevalence of tuberculosis (TB) in the same population (n=529). Seroconversion against SARS-CoV-2 antigens was considered positive if at least one of the three ELISA kits detected levels of specific antibodies above the threshold specified by the manufacturer. In this sense, we report 56.0% (n=349/623) of seroconverted individuals. Relative seroconversion peaked after introduction of the Coronavac vaccine in February 2021. Vaccination increased the production of anti-S IgG from 3.9% to 48.6%. Our results also indicated that 11.0% (n=46/417) of all individuals were positive for TB. Seroconversion to SARS-CoV-2 was similar between individuals with positive tuberculosis test results to those with negative test results. Most vaccinated individuals seroconverted to SARS-CoV-2, indicating that Coronavac may be as protective in individuals from these indigenous groups as observed in the general Brazilian population. COVID-19 severity was minimal regardless of incomplete vaccine coverage, suggesting that vaccination may not be the only factor protecting individuals from severe COVID-19. Tuberculosis is highly prevalent and not associated with increased seroconversion to SARS-CoV-2.

Host genetics and the profile of COVID-19 in indigenous people from the Brazilian Amazon: A pilot study with variants of the ACE1, ACE2 and TMPRSS2 genes.

This pilot study aimed to investigate genetic factors that may have contributed to the milder clinical outcomes of COVID-19 in Brazilian indigenous populations. 263 Indigenous from the Araweté, Kararaô, Parakanã, Xikrin do Bacajá, Kayapó and Munduruku peoples were analyzed, 55.2% women, ages ranging from 10 to 95 years (average 49.5 ± 20.7). Variants in genes involved in the entry of SARS-CoV-2 into the host cell (ACE1 rs1799752 I/D, ACE2 rs2285666 C/T, ACE2 rs73635825 A/G and TMPRSS2 rs123297605 C/T), were genotyped in indigenous peoples from the Brazilian Amazon, treated during the SARS-CoV-2 pandemic between 2020 and 2021. The distribution of genotypes did not show any association with the presence or absence of IgG antibodies. Additionally, the influence of genetic variations on the severity of the disease was not examined extensively because a significant number of indigenous individuals experienced the disease with either mild symptoms or no symptoms. It is worth noting that the frequencies of risk alleles were found to be lower in Indigenous populations compared to both continental populations and Brazilians. Indigenous Brazilian Amazon people exhibited an ethnic-specific genetic profile that may be associated with a milder disease, which could explain the unexpected response they demonstrated to COVID-19, being less impacted than Brazilians.

A model for the dynamics of expanded CAG repeat alleles: ATXN2 and ATXN3 as prototypes.

Background: Spinocerebellar ataxia types 2 (SCA2) and 3 (SCA3/MJD) are diseases due to dominant unstable expansions of CAG repeats (CAGexp). Age of onset of symptoms (AO) correlates with the CAGexp length. Repeat instability leads to increases in the expanded repeats, to important AO anticipations and to the eventual extinction of lineages. Because of that, compensatory forces are expected to act on the maintenance of expanded alleles, but they are poorly understood. Objectives: we described the CAGexp dynamics, adapting a classical equation and aiming to estimate for how many generations will the descendants of a de novo expansion last. Methods: A mathematical model was adapted to encompass anticipation, fitness, and allelic segregation; and empirical data fed the model. The arbitrated ancestral mutations included in the model had the lowest CAGexp and the highest AO described in the literature. One thousand generations were simulated until the alleles were eliminated, fixed, or 650 generations had passed. Results: All SCA2 lineages were eliminated in a median of 10 generations. In SCA3/MJD lineages, 593 were eliminated in a median of 29 generations. The other ones were eliminated due to anticipation after the 650th generation or remained indefinitely with CAG repeats transitioning between expanded and unexpanded ranges. Discussion: the model predicted outcomes compatible with empirical data - the very old ancestral SCA3/MJD haplotype, and the de novo SCA2 expansions -, which previously seemed to be contradictory. This model accommodates these data into understandable dynamics and might be useful for other CAGexp disorders.

Extracellular matrix and vascular dynamics in the kidney of a murine model for Marfan syndrome.

Fibrillin-1 is a pivotal structural component of the kidney's glomerulus and peritubular tissue. Mutations in the fibrillin-1 gene result in Marfan syndrome (MFS), an autosomal dominant disease of the connective tissue. Although the kidney is not considered a classically affected organ in MFS, several case reports describe glomerular disease in patients. Therefore, this study aimed to characterize the kidney in the mgΔlpn-mouse model of MFS. Affected animals presented a significant reduction of glomerulus, glomerulus-capillary, and urinary space, and a significant reduction of fibrillin-1 and fibronectin in the glomerulus. Transmission electron microscopy and 3D-ultrastructure analysis revealed decreased amounts of microfibrils which also appeared fragmented in the MFS mice. Increased collagen fibers types I and III, MMP-9, and α-actin were also observed in affected animals, suggesting a tissue-remodeling process in the kidney. Video microscopy analysis showed an increase of microvessel distribution coupled with reduction of blood-flow velocity, while ultrasound flow analysis revealed significantly lower blood flow in the kidney artery and vein of the MFS mice. The structural and hemodynamic changes of the kidney indicate the presence of kidney remodeling and vascular resistance in this MFS model. Both processes are associated with hypertension which is expected to worsen the cardiovascular phenotype in MFS.

Genomic history of coastal societies from eastern South America.

Sambaqui (shellmound) societies are among the most intriguing archaeological phenomena in pre-colonial South America, extending from approximately 8,000 to 1,000 years before present (yr BP) across 3,000 km on the Atlantic coast. However, little is known about their connection to early Holocene hunter-gatherers, how this may have contributed to different historical pathways and the processes through which late Holocene ceramists came to rule the coast shortly before European contact. To contribute to our understanding of the population history of indigenous societies on the eastern coast of South America, we produced genome-wide data from 34 ancient individuals as early as 10,000 yr BP from four different regions in Brazil. Early Holocene hunter-gatherers were found to lack shared genetic drift among themselves and with later populations from eastern South America, suggesting that they derived from a common radiation and did not contribute substantially to later coastal groups. Our analyses show genetic heterogeneity among contemporaneous Sambaqui groups from the southeastern and southern Brazilian coast, contrary to the similarity expressed in the archaeological record. The complex history of intercultural contact between inland horticulturists and coastal populations becomes genetically evident during the final horizon of Sambaqui societies, from around 2,200 yr BP, corroborating evidence of cultural change.

Unraveling the Genetic Architecture of Hepatoblastoma Risk: Birth Defects and Increased Burden of Germline Damaging Variants in Gastrointestinal/Renal Cancer Predisposition and DNA Repair Genes.

The ultrarare hepatoblastoma (HB) is the most common pediatric liver cancer. HB risk is related to a few rare syndromes, and the molecular bases remain elusive for most cases. We investigated the burden of rare damaging germline variants in 30 Brazilian patients with HB and the presence of additional clinical signs. A high frequency of prematurity (20%) and birth defects (37%), especially craniofacial (17%, including craniosynostosis) and kidney (7%) anomalies, was observed. Putative pathogenic or likely pathogenic monoallelic germline variants mapped to 10 cancer predisposition genes (CPGs: APC, CHEK2, DROSHA, ERCC5, FAH, MSH2, MUTYH, RPS19, TGFBR2 and VHL) were detected in 33% of the patients, only 40% of them with a family history of cancer. These findings showed a predominance of CPGs with a known link to gastrointestinal/colorectal and renal cancer risk. A remarkable feature was an enrichment of rare damaging variants affecting different classes of DNA repair genes, particularly those known as Fanconi anemia genes. Moreover, several potentially deleterious variants mapped to genes impacting liver functions were disclosed. To our knowledge, this is the largest assessment of rare germline variants in HB patients to date, contributing to elucidate the genetic architecture of HB risk.

DNA methylation as a key epigenetic player for hepatoblastoma characterization.

BACKGROUND: Hepatoblastoma (HB) is a rare embryonal liver tumor of children. Although intrinsic biological differences between tumors can affect prognosis, few groups have studied these differences. Given the recent increased attention to epigenetic mechanisms in the genesis and progression of these tumors, we aimed to classify HB samples according to the stages of liver development and DNA methylation machinery. BASIC PROCEDURES: We evaluated the expression of 24 genes associated with DNA methylation and stages of hepatocyte differentiation and global DNA methylation. Using bioinformatics tools and expression data, we propose a stratification model for HB. MAIN FINDINGS: Tumors clustered into three groups that presented specific gene expression profiles of the panel of DNA methylation enzymes and hepatocyte differentiation markers. In addition to reinforcing these embryonal tumors' molecular heterogeneity, we propose that a panel of 13 genes can stratify HBs (TET1, TET2, TET3, DNMT1, DNMT3A, UHRF1, ALB, CYP3A4, TDO2, UGT1A1, AFP, HNF4A, and FOXA2). DNA methylation machinery participates in the characterization of HBs, directly reflected in diverse DNA methylation content. The data suggested that a subset of HBs were similar to differentiated livers, with upregulation of mature hepatocyte markers, decreased expression of DNA methylation enzymes, and higher global methylation levels; these findings might predict worse outcomes. CONCLUSIONS: HBs are heterogeneous tumors. Despite using a small cohort of 21 HB samples, our findings reinforce that DNA methylation is a robust biomarker for this tumor type.

Hyperkyphosis is not dependent on bone mass and quality in the mouse model of Marfan syndrome.

Marfan syndrome (MFS) is an autosomal dominant disease affecting cardiovascular, ocular and skeletal systems. It is caused by mutations in the fibrillin-1 (FBN1) gene, leading to structural defects of connective tissue and increased activation of TGF-ß. Angiotensin II (ang-II) is involved in TGF-ß activity and in bone mass regulation. Inhibition of TGF-ß signaling by blockage of the ang-II receptor 1 (AT1R) via losartan administration leads to improvement of cardiovascular and pulmonary phenotypes, but has no effect on skeletal phenotype in the haploinsufficient mouse model of MFS mgR, suggesting a distinct mechanism of pathogenesis in the skeletal system. Here we characterized the skeletal phenotypes of the dominant-negative model for MFS mgΔlpn and tested the effect of inhibition of ang-II signaling in improving those phenotypes. As previously shown, heterozygous mice present hyperkyphosis, however we now show that only males also present osteopenia. Inhibition of ang-II production by ramipril minimized the kyphotic deformity, but had no effect on bone microstructure in male mutant animals. Histological analysis revealed increased thickness of the anterior longitudinal ligament (ALL) of the spine in mutant animals (25.8 ± 6.3 vs. 29.7 ± 7.7 µm), coupled with a reduction in type I (164.1 ± 8.7 vs. 139.0 ± 4.4) and increase in type III (86.5 ± 10.2 vs. 140.4 ± 5.6) collagen in the extracellular matrix of this ligament. In addition, we identified in the MFS mice alterations in the erector spinae muscles which presented thinner muscle fibers (1035.0 ± 420.6 vs. 655.6 ± 239.5 µm2) surrounded by increased area of connective tissue (58.17 ± 6.52 vs. 105.0 ± 44.54 µm2). Interestingly, these phenotypes were ameliorated by ramipril treatment. Our results reveal a sex-dependency of bone phenotype in MFS, where females do not present alterations in bone microstructure. More importantly, they indicate that hyperkyphosis is not a result of osteopenia in the MFS mouse model, and suggest that incompetent spine ligaments and muscles are responsible for the development of that phenotype.

The structure of first-cousin marriages in Brazil.

This paper deals with the frequency and structure of first-cousin marriages, by far the most important and frequent type of consanguineous mating in human populations. Based on the analysis of large amounts of data from the world literature and from large Brazilian samples recently collected, we suggest some explanations for the asymmetry of sexes among the parental sibs of first-cousin marriages. We suggest also a simple manner to correct the method that uses population surnames to assess the different Wright fixation indexes FIS, FST and FIT taking into account not only alternative methods of surname transmission, but also the asymmetries that are almost always observed in the distribution of sexes among the parental sibs of first-cousins.

Inbreeding estimates in human populations: Applying new approaches to an admixed Brazilian isolate.

The analysis of genomic data (~400,000 autosomal SNPs) enabled the reliable estimation of inbreeding levels in a sample of 541 individuals sampled from a highly admixed Brazilian population isolate (an African-derived quilombo in the State of São Paulo). To achieve this, different methods were applied to the joint information of two sets of markers (one complete and another excluding loci in patent linkage disequilibrium). This strategy allowed the detection and exclusion of markers that biased the estimation of the average population inbreeding coefficient (Wright's fixation index FIS), which value was eventually estimated as around 1% using any of the methods we applied. Quilombo demographic inferences were made by analyzing the structure of runs of homozygosity (ROH), which were adapted to cope with a highly admixed population with a complex foundation history. Our results suggest that the amount of ROH <2Mb of admixed populations should be somehow proportional to the genetic contribution from each parental population.

Selection scan reveals three new loci related to high altitude adaptation in Native Andeans.

The Andean Altiplano has been occupied continuously since the late Pleistocene, ~12,000 years ago, which places the Andean natives as one of the most ancient populations living at high altitudes. In the present study, we analyzed genomic data from Native Americans living a long-time at Andean high altitude and at Amazonia and Mesoamerica lowland areas. We have identified three new candidate genes - SP100, DUOX2 and CLC - with evidence of positive selection for altitude adaptation in Andeans. These genes are involved in the TP53 pathway and are related to physiological routes important for high-altitude hypoxia response, such as those linked to increased angiogenesis, skeletal muscle adaptations, and immune functions at the fetus-maternal interface. Our results, combined with other studies, showed that Andeans have adapted to the Altiplano in different ways and using distinct molecular strategies as compared to those of other natives living at high altitudes.