Ancestry, admixture and fitness in Colombian genomes

Prevalence of axial spondyloarthritis in Colombia: data from the National Health Registry 2017-2021

INTRODUCTION

Registries allow ascertaining the epidemiology of chronic diseases such as axial spondyloarthritis (axSpA). The Colombian Ministry of Health has implemented a National Health Registry (SISPRO) that collects data from each medical contact in the system, which provides close to universal coverage (around 98%).

OBJECTIVE

To establish the 5-year prevalence of axSpA in Colombia, and to describe its demographics, using data from January 1st, 2017, to December 31st, 2021.

METHODS

We performed an observational, cross-sectional study using the International Statistical Classification of Diseases and Related Health Problems as search terms related to ax-SpA, based on SISPRO data. We estimated the prevalence using three approaches: (1) ankylosing spondylitis (AS) diagnoses; (2) diagnoses compatible with axSpA; and (3) diagnoses compatible with axSpA, including sacroiliitis. We calculated prevalence per 100,000 inhabitants.

RESULTS

Based on our three approaches, patients with a primary diagnosis compatible with ax-SpA ranged between 12,684 and 117,648, with an estimated 5-year adjusted prevalence between 26.3 and 244 cases per 100,000 inhabitants (0.03-0.2%). The male-to-female ratio ranged between 1.2:1 and 0.4:1, which was markedly skewed towards a higher prevalence in women when we included the code for sacroiliitis. We found the highest frequency of cases in the 50-54 years group. A differential prevalence was observed between different regions in our country, particularly in regions known to have European ancestors.

CONCLUSION

This is the first study that describes demographic characteristics of ax-SpA in Colombia and offers valuable information for stakeholders. Key Points • Using the official country-level health database, the prevalence of axSpA in Colombia ranges between 26.3 and 244 cases per 100,000 inhabitants (0.03% - 0.2%) • The prevalence of axSpA peaked among the 50-54 years patient group, suggesting an increased survival • Nations with a substantial admixture, such as Colombia, may present a differential prevalence of axSpA among regions within the country • Including the ICD-10 code for sacroiliitis (M46.1) in epidemiological studies probably overestimates the frequency of axSpA.

Genomic and demographic processes differentially influence genetic variation across the human X chromosome

Many forces influence genetic variation across the genome including mutation, recombination, selection, and demography. Increased mutation and recombination both lead to increases in genetic diversity in a region-specific manner, while complex demographic patterns shape patterns of diversity on a more global scale. While these processes act across the entire genome, the X chromosome is particularly interesting because it contains several distinct regions that are subject to different combinations and strengths of these forces: the pseudoautosomal regions (PARs) and the X-transposed region (XTR). The X chromosome thus can serve as a unique model for studying how genetic and demographic forces act in different contexts to shape patterns of observed variation. We therefore sought to explore diversity, divergence, and linkage disequilibrium in each region of the X chromosome using genomic data from 26 human populations. Across populations, we find that both diversity and substitution rate are consistently elevated in PAR1 and the XTR compared to the rest of the X chromosome. In contrast, linkage disequilibrium is lowest in PAR1, consistent with the high recombination rate in this region, and highest in the region of the X chromosome that does not recombine in males. However, linkage disequilibrium in the XTR is intermediate between PAR1 and the autosomes, and much lower than the non-recombining X. Finally, in addition to these global patterns, we also observed variation in ratios of X versus autosomal diversity consistent with population-specific evolutionary history as well. While our results were generally consistent with previous work, two unexpected observations emerged. First, our results suggest that the XTR does not behave like the rest of the recombining X and may need to be evaluated separately in future studies. Second, the different regions of the X chromosome appear to exhibit unique patterns of linked selection across different human populations. Together, our results highlight profound regional differences across the X chromosome, simultaneously making it an ideal system for exploring the action of evolutionary forces as well as necessitating its careful consideration and treatment in genomic analyses.

Inferring multi-locus selection in admixed populations

Admixture, the exchange of genetic information between distinct source populations, is thought to be a major source of adaptive genetic variation. Unlike mutation events, which periodically generate single alleles, admixture can introduce many selected alleles simultaneously. As such, the effects of linkage between selected alleles may be especially pronounced in admixed populations. However, existing tools for identifying selected mutations within admixed populations only account for selection at a single site, overlooking phenomena such as linkage among proximal selected alleles. Here, we develop and extensively validate a method for identifying and quantifying the individual effects of multiple linked selected sites on a chromosome in admixed populations. Our approach numerically calculates the expected local ancestry landscape in an admixed population for a given multi-locus selection model, and then maximizes the likelihood of the model. After applying this method to admixed populations of Drosophila melanogaster and Passer italiae, we found that the impacts between linked sites may be an important contributor to natural selection in admixed populations. Furthermore, for the situations we considered, the selection coefficients and number of selected sites are overestimated in analyses that do not consider the effects of linkage among selected sites. Our results imply that linkage among selected sites may be an important evolutionary force in admixed populations. This tool provides a powerful generalized method to investigate these crucial phenomena in diverse populations.

Riding the merry-go-round of racial disparities in ADRD research

INTRODUCTION

With the rapid expansion of the aging population, the burden of Alzheimer's disease related dementias (ADRD) is anticipated to increase in racialized and minoritized groups who are at disproportionately higher risk. To date, research emphasis has been on further characterizing the existence of racial disparities in ADRD through comparisons to groups racialized as White that are assumed to be normative. Much of the literature on this comparison insinuates that racialized and minoritized groups experience poorer outcomes due to genetics, culture, and/or health behaviors.

METHODS

This perspective shines a light on a category of ADRD research that employs ahistorical methodological approaches to describe racial disparities in ADRD that puts us on a merry-go-round of research with no benefits to society.

METHODS

This commentary provides historical context for the use of race in ADRD research and justification for the study of structural racism. The commentary concludes with recommendations to guide future research.

Localizing Post-Admixture Adaptive Variants with Object Detection on Ancestry-Painted Chromosomes

Gene flow between previously differentiated populations during the founding of an admixed or hybrid population has the potential to introduce adaptive alleles into the new population. If the adaptive allele is common in one source population, but not the other, then as the adaptive allele rises in frequency in the admixed population, genetic ancestry from the source containing the adaptive allele will increase nearby as well. Patterns of genetic ancestry have therefore been used to identify post-admixture positive selection in humans and other animals, including examples in immunity, metabolism, and animal coloration. A common method identifies regions of the genome that have local ancestry "outliers" compared with the distribution across the rest of the genome, considering each locus independently. However, we lack theoretical models for expected distributions of ancestry under various demographic scenarios, resulting in potential false positives and false negatives. Further, ancestry patterns between distant sites are often not independent. As a result, current methods tend to infer wide genomic regions containing many genes as under selection, limiting biological interpretation. Instead, we develop a deep learning object detection method applied to images generated from local ancestry-painted genomes. This approach preserves information from the surrounding genomic context and avoids potential pitfalls of user-defined summary statistics. We find the method is robust to a variety of demographic misspecifications using simulated data. Applied to human genotype data from Cabo Verde, we localize a known adaptive locus to a single narrow region compared with multiple or long windows obtained using two other ancestry-based methods.

Determination of TMPRSS2-ERG, SPOP, FOXA1, and IDH1 prostate cancer molecular subtypes in Colombian patients and their possible implications for prognosis

Prostate cancer (PCa) is one of cancer with of the highest incidence and mortality worldwide. Current disease prognostic markers do not differentiate aggressive from indolent PCa with sufficient certainty, and characterization by molecular subtypes has been sought to allow a better classification. TMPRSS2-ERG, SPOP, FOXA1, and IDH1 molecular subtypes have been described, but the association of these subtypes with prognosis in PCa is unclear; their frequency in Colombian patients is also unknown. Formalin-fixed and paraffin-embedded samples of radical prostatectomy from 112 patients with PCa were used. The TMPRSS2-ERG subtype was assessed with fluorescent in situ hybridization. The mutations in SPOP, FOXA1, and IDH1 in hot-spot regions were evaluated using Sanger sequencing. Fusion was detected in 71 patients (63.4%). No statistically significant differences were found between the state of fusion and the variables analyzed. In the 41 fusion-negative cases (36.6%), two patients (4.9%) had missense mutations in SPOP (p.F102C and p.F133L), representing a 1.8% of the overall cohort. The low frequency of this subtype in Colombians could be explained by the reported variability in the frequency of these mutations according to the population (5%-20%). No mutations were found in FOXA1 in the cases analyzed. The synonym SNP rs11554137 IDH1^105GGT was found in tumor tissue but not in the normal tissue in one case. A larger cohort of Colombian PCa patients is needed for future studies to validate these findings and gain a better understanding of the molecular profile of this cancer in our population and if there are any differences by Colombian regions.

Sleep problems in low income, urban pediatric populations living at different altitudes in Colombia

OBJECTIVE

To determine the frequency of sleep problems in low-income, urban pediatric populations in cities at different altitudes in Colombia.

METHODS

A descriptive, cross-sectional population-based observational study was conducted in children aged between 2 and 12 years in the low income, urban areas of three cities in Colombia (Santa Marta, Bucaramanga, and Bogotá) located at 15, 959, and 2640 m above sea level, respectively. Sociodemographic data were collected, and the Spanish version of the Pediatric Sleep Questionnaire was used.

RESULTS

1989 children were surveyed, distributed as follows: Santa Marta (32.0%), Bucaramanga (33.4%), and Bogotá (34.6%). The overall prevalence of sleep problems was 39.0%. Children from Santa Marta had the highest frequency of parasomnias (58.0%); those from Bucaramanga had the highest frequency of attention deficit symptoms (4.0%) and apneic pauses witnessed by parents or caregivers (5.7%). Finally, Bogotá, the only high-altitude location, had the highest frequency of sleep disordered breathing (17.2%).

CONCLUSIONS

The study found a high frequency of sleep problems in the pediatric population, especially at higher altitudes when compared to lower altitude settings. Sleep disorders warrant early detection and timely therapeutic intervention.

SALAI-Net: species-agnostic local ancestry inference network

MOTIVATION

Local ancestry inference (LAI) is the high resolution prediction of ancestry labels along a DNA sequence. LAI is important in the study of human history and migrations, and it is beginning to play a role in precision medicine applications including ancestry-adjusted genome-wide association studies (GWASs) and polygenic risk scores (PRSs). Existing LAI models do not generalize well between species, chromosomes or even ancestry groups, requiring re-training for each different setting. Furthermore, such methods can lack interpretability, which is an important element in each of these applications.

RESULTS

We present SALAI-Net, a portable statistical LAI method that can be applied on any set of species and ancestries (species-agnostic), requiring only haplotype data and no other biological parameters. Inspired by identity by descent methods, SALAI-Net estimates population labels for each segment of DNA by performing a reference matching approach, which leads to an interpretable and fast technique. We benchmark our models on whole-genome data of humans and we test these models' ability to generalize to dog breeds when trained on human data. SALAI-Net outperforms previous methods in terms of balanced accuracy, while generalizing between different settings, species and datasets. Moreover, it is up to two orders of magnitude faster and uses considerably less RAM memory than competing methods.

AVAILABILITY AND IMPLEMENTATION

We provide an open source implementation and links to publicly available data at github.com/AI-sandbox/SALAI-Net. Data is publicly available as follows: https://www.internationalgenome.org (1000 Genomes), https://www.simonsfoundation.org/simons-genome-diversity-project (Simons Genome Diversity Project), https://www.sanger.ac.uk/resources/downloads/human/hapmap3.html (HapMap), ftp://ngs.sanger.ac.uk/production/hgdp/hgdp_wgs.20190516 (Human Genome Diversity Project) and https://www.ncbi.nlm.nih.gov/bioproject/PRJNA448733 (Canid genomes).

SUPPLEMENTARY INFORMATION

Supplementary data are available from Bioinformatics online.

Killer-Cell Immunoglobulin-like Receptor Diversity in an Admixed South American Population

Natural Killer (NK) cells are innate immune cells that mediate antiviral and antitumor responses. NK cell activation and induction of effector functions are tightly regulated by the integration of activating and inhibitory receptors such as killer immunoglobulin-like receptors (KIR). KIR genes are characterized by a high degree of diversity due to presence or absence, gene copy number and allelic polymorphism. The aim of this study was to establish the distribution of KIR genes and genotypes, to infer the most common haplotypes in an admixed Colombian population and to compare these KIR gene frequencies with some Central and South American populations and worldwide. A total of 161 individuals from Medellin, Colombia were included in the study. Genomic DNA was used for KIR and HLA genotyping. We analyzed only KIR gene-content (presence or absence) based on PCR-SSO. The KIR genotype, most common haplotypes and combinations of KIR and HLA ligands frequencies were estimated according to the presence or absence of KIR and HLA genes. Dendrograms, principal component (PC) analysis and Heatmap analysis based on genetic distance were constructed to compare KIR gene frequencies among Central and South American, worldwide and Amerindian populations. The 16 KIR genes analyzed were distributed in 37 different genotypes and the 7 most frequent KIR inferred haplotypes. Importantly, we found three new genotypes not previously reported in any other ethnic group. Our genetic distance, PC and Heatmap analysis revealed marked differences in the distribution of KIR gene frequencies in the Medellin population compared to worldwide populations. These differences occurred mainly in the activating KIR isoforms, which are more frequent in our population, particularly KIR3DS1. Finally, we observed unique structural patterns of genotypes, which evidences the potential diversity and variability of this gene family in our population, and the need for exhaustive genetic studies to expand our understanding of the KIR gene complex in Colombian populations.

The immunogenetic impact of European colonization in the Americas

The introduction of pathogens originating from Eurasia into the Americas during early European contact has been associated with high mortality rates among Indigenous peoples, likely contributing to their historical and precipitous population decline. However, the biological impacts of imported infectious diseases and resulting epidemics, especially in terms of pathogenic effects on the Indigenous immunity, remain poorly understood and highly contentious to this day. Here, we examine multidisciplinary evidence underpinning colonization-related immune genetic change, providing contextualization from anthropological studies, paleomicrobiological evidence of contrasting host-pathogen coevolutionary histories, and the timings of disease emergence. We further summarize current studies examining genetic signals reflecting post-contact Indigenous population bottlenecks, admixture with European and other populations, and the putative effects of natural selection, with a focus on ancient DNA studies and immunity-related findings. Considering current genetic evidence, together with a population genetics theoretical approach, we show that post-contact Indigenous immune adaptation, possibly influenced by selection exerted by introduced pathogens, is highly complex and likely to be affected by multifactorial causes. Disentangling putative adaptive signals from those of genetic drift thus remains a significant challenge, highlighting the need for the implementation of population genetic approaches that model the short time spans and complex demographic histories under consideration. This review adds to current understandings of post-contact immunity evolution in Indigenous peoples of America, with important implications for bettering our understanding of human adaptation in the face of emerging infectious diseases.

Association Between the LZTFL1 rs11385942 Polymorphism and COVID-19 Severity in Colombian Population

Genetic and non-genetic factors are responsible for the high interindividual variability in the response to SARS-CoV-2. Although numerous genetic polymorphisms have been identified as risk factors for severe COVID-19, these remain understudied in Latin-American populations. This study evaluated the association of non-genetic factors and three polymorphisms: ACE rs4646994, ACE2 rs2285666, and LZTFL1 rs11385942, with COVID severity and long-term symptoms by using a case-control design. The control group was composed of asymptomatic/mild cases (n = 61) recruited from a private laboratory, while the case group was composed of severe/critical patients (n = 63) hospitalized in the Hospital Universitario Mayor-Méderi, both institutions located in Bogotá, Colombia. Clinical follow up and exhaustive revision of medical records allowed us to assess non-genetic factors. Genotypification of the polymorphism of interest was performed by amplicon size analysis and Sanger sequencing. In agreement with previous reports, we found a statistically significant association between age, male sex, and comorbidities, such as hypertension and type 2 diabetes mellitus (T2DM), and worst outcomes. We identified the polymorphism LZTFL1 rs11385942 as an important risk factor for hospitalization (p < 0.01; OR = 5.73; 95% CI = 1.2–26.5, under the allelic test). Furthermore, long-term symptoms were common among the studied population and associated with disease severity. No association between the polymorphisms examined and long-term symptoms was found. Comparison of allelic frequencies with other populations revealed significant differences for the three polymorphisms investigated. Finally, we used the statistically significant genetic and non-genetic variables to develop a predictive logistic regression model, which was implemented in a Shiny web application. Model discrimination was assessed using the area under the receiver operating characteristic curve (AUC = 0.86; 95% confidence interval 0.79–0.93). These results suggest that LZTFL1 rs11385942 may be a potential biomarker for COVID-19 severity in addition to conventional non-genetic risk factors. A better understanding of the impact of these genetic risk factors may be useful to prioritize high-risk individuals and decrease the morbimortality caused by SARS-CoV2 and future pandemics.

The Bank Vole (Clethrionomys glareolus) as a Model System for Adaptive Phylogeography in the European Theater

The legacy of climatic changes during the Pleistocene glaciations allows inferences to be made about the patterns and processes associated with range expansion/colonization, including evolutionary adaptation. With the increasing availability of population genomic data, we have the opportunity to examine these questions in detail and in a variety of non-traditional model species. As an exemplar, here we review more than two decades of work by our group and others that illustrate the potential of a single “non-model model” mammal species - the bank vole (Clethrionomys glareolus), which is particularly well suited to illustrate the complexities that may be associated with range expansion and the power of genomics (and other datasets) to uncover them. We first summarize early phylogeographic work using mitochondrial DNA and then describe new phylogeographic insights gained from population genomic analysis of genome-wide SNP data to highlight the bank vole as one of the most compelling examples of a forest mammal, that survived in cryptic extra-Mediterranean (“northern”) glacial refugia in Europe, and as one of the species in which substantial replacement and mixing of lineages originating from different refugia occurred during end-glacial colonization. Our studies of bank vole hemoglobin structure and function, as well as our recent ecological niche modeling study examining differences among bank vole lineages, led us to develop the idea of “adaptive phylogeography.” This is what we call the study of the role of adaptive differences among populations in shaping phylogeographic patterns. Adaptive phylogeography provides a link between past population history and adaptation that can ultimately help predict the potential of future species responses to climate change. Because the bank vole is part of a community of organisms whose range has repeatedly contracted and then expanded in the past, what we learn from the bank vole will be useful for our understanding of a broad range of species.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

The genomic signatures of natural selection in admixed human populations

Admixture has been a pervasive phenomenon in human history, extensively shaping the patterns of population genetic diversity. There is increasing evidence to suggest that admixture can also facilitate genetic adaptation to local environments, i.e., admixed populations acquire beneficial mutations from source populations, a process that we refer to as "adaptive admixture." However, the role of adaptive admixture in human evolution and the power to detect it remain poorly characterized. Here, we use extensive computer simulations to evaluate the power of several neutrality statistics to detect natural selection in the admixed population, assuming multiple admixture scenarios. We show that statistics based on admixture proportions, F_adm and LAD, show high power to detect mutations that are beneficial in the admixed population, whereas other statistics, including iHS and F_ST, falsely detect neutral mutations that have been selected in the source populations only. By combining F_adm and LAD into a single, powerful statistic, we scanned the genomes of 15 worldwide, admixed populations for signatures of adaptive admixture. We confirm that lactase persistence and resistance to malaria have been under adaptive admixture in West Africans and in Malagasy, North Africans, and South Asians, respectively. Our approach also uncovers other cases of adaptive admixture, including APOL1 in Fulani nomads and PKN2 in East Indonesians, involved in resistance to infection and metabolism, respectively. Collectively, our study provides evidence that adaptive admixture has occurred in human populations whose genetic history is characterized by periods of isolation and spatial expansions resulting in increased gene flow.

A neurodegenerative disease landscape of rare mutations in Colombia due to founder effects

BACKGROUND

The Colombian population, as well as those in other Latin American regions, arose from a recent tri-continental admixture among Native Americans, Spanish invaders, and enslaved Africans, all of whom passed through a population bottleneck due to widespread infectious diseases that left small isolated local settlements. As a result, the current population reflects multiple founder effects derived from diverse ancestries.

METHODS

We characterized the role of admixture and founder effects on the origination of the mutational landscape that led to neurodegenerative disorders under these historical circumstances. Genomes from 900 Colombian individuals with Alzheimer's disease (AD) [n = 376], frontotemporal lobar degeneration-motor neuron disease continuum (FTLD-MND) [n = 197], early-onset dementia not otherwise specified (EOD) [n = 73], and healthy participants [n = 254] were analyzed. We examined their global and local ancestry proportions and screened this cohort for deleterious variants in disease-causing and risk-conferring genes.

RESULTS

We identified 21 pathogenic variants in AD-FTLD related genes, and PSEN1 harbored the majority (11 pathogenic variants). Variants were identified from all three continental ancestries. TREM2 heterozygous and homozygous variants were the most common among AD risk genes (102 carriers), a point of interest because the disease risk conferred by these variants differed according to ancestry. Several gene variants that have a known association with MND in European populations had FTLD phenotypes on a Native American haplotype. Consistent with founder effects, identity by descent among carriers of the same variant was frequent.

CONCLUSIONS

Colombian demography with multiple mini-bottlenecks probably enhanced the detection of founder events and left a proportionally higher frequency of rare variants derived from the ancestral populations. These findings demonstrate the role of genomically defined ancestry in phenotypic disease expression, a phenotypic range of different rare mutations in the same gene, and further emphasize the importance of inclusiveness in genetic studies.

The Impact of Ethnicity and Genetic Ancestry on Disease Prevalence and Risk in Colombia

Currently, the vast majority of genomic research cohorts are made up of participants with European ancestry. Genomic medicine will only reach its full potential when genomic studies become more broadly representative of global populations. We are working to support the establishment of genomic medicine in developing countries in Latin America via studies of ethnically and ancestrally diverse Colombian populations. The goal of this study was to analyze the effect of ethnicity and genetic ancestry on observed disease prevalence and predicted disease risk in Colombia. Population distributions of Colombia's three major ethnic groups - Mestizo, Afro-Colombian, and Indigenous - were compared to disease prevalence and socioeconomic indicators. Indigenous and Mestizo ethnicity show the highest correlations with disease prevalence, whereas the effect of Afro-Colombian ethnicity is substantially lower. Mestizo ethnicity is mostly negatively correlated with six high-impact health conditions and positively correlated with seven of eight common cancers; Indigenous ethnicity shows the opposite effect. Malaria prevalence in particular is strongly correlated with ethnicity. Disease prevalence co-varies across geographic regions, consistent with the regional distribution of ethnic groups. Ethnicity is also correlated with regional variation in human development, partially explaining the observed differences in disease prevalence. Patterns of genetic ancestry and admixture for a cohort of 624 individuals from Medellín were compared to disease risk inferred via polygenic risk scores (PRS). African genetic ancestry is most strongly correlated with predicted disease risk, whereas European and Native American ancestry show weaker effects. African ancestry is mostly positively correlated with disease risk, and European ancestry is mostly negatively correlated. The relationships between ethnicity and disease prevalence do not show an overall correspondence with the relationships between ancestry and disease risk. We discuss possible reasons for the divergent health effects of ethnicity and ancestry as well as the implication of our results for the development of precision medicine in Colombia.

Genetic ancestry and ethnic identity in Ecuador

We investigated the ancestral origins of four Ecuadorian ethnic groups-Afro-Ecuadorian, Mestizo, Montubio, and the Indigenous Tsáchila-in an effort to gain insight on the relationship between ancestry, culture, and the formation of ethnic identities in Latin America. The observed patterns of genetic ancestry are largely concordant with ethnic identities and historical records of conquest and colonization in Ecuador. Nevertheless, a number of exceptional findings highlight the complex relationship between genetic ancestry and ethnicity in Ecuador. Afro-Ecuadorians show far less African ancestry, and the highest levels of Native American ancestry, seen for any Afro-descendant population in the Americas. Mestizos in Ecuador show high levels of Native American ancestry, with substantially less European ancestry, despite the relatively low Indigenous population in the country. The recently recognized Montubio ethnic group is highly admixed, with substantial contributions from all three continental ancestries. The Tsáchila show two distinct ancestry subgroups, with most individuals showing almost exclusively Native American ancestry and a smaller group showing a Mestizo characteristic pattern. Considered together with historical data and sociological studies, our results indicate the extent to which ancestry and culture interact, often in unexpected ways, to shape ethnic identity in Ecuador.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians

Background

Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors).

Methods

In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors.

Results

The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors.

Conclusions

Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01900-y.

HLA-B Allele, Genotype, and Haplotype Frequencies in a Group of Healthy Individuals in Colombia

BACKGROUND

The sequencing of alleles of the HLA-B, a human leukocyte antigen (HLA) class I gene, was established as the most polymorphic of chromosome 6 and of the entire human genome. In this locus, the HLA-B*27 allele is highly polymorphic and has clinical relevance. Literature about the subtypes and singular frequency of these alleles in Colombia's healthy population is scarce.

OBJECTIVE

The aim of this study was to establish the HLA-B allele, genotype, and haplotype frequencies in a healthy Colombian population and analyze their association with the sex and geographical distribution of the individuals studied.

METHODS

This is a nonexperimental and descriptive study. The data from whole-blood samples whose HLA genes were genotyped by protocol with the Luminex 100/200 xMAP technology were evaluated. HLA-B*27 positivity was confirmed by the new-generation sequencing technology. The associations between the HLA-B alleles and demographic variables were evaluated by χ2 and Fisher exact tests.

RESULTS

Twenty-seven HLA-B genotypes were identified in 255 individuals, with the highest frequencies for HLA-B*35 (44.7%), B*40 (19.6%), and B*44 (16.8%). Additionally, 89 HLA-B alleles were found; the most common were HLA-B*35:01 (6.7%) and B*40:02 (6.5%). Nine individuals tested positive for the HLA-B*27 allele with genotype and allele frequencies of 3.5% and 1.8%, respectively; the HLA-B*27:05:02 subtype predominated.

CONCLUSIONS

Here, we report the most common HLA-B allele, genotype, and haplotype frequencies in a healthy Colombian population group and analyzed their association with the sex and geographical distribution of the individuals studied. Results for the HLA-B*27 allele confirm racial mixing in Colombia with a high degree of Caucasian influence, as well as the repopulation of Colombia's central region, attributed to the migration phenomena. Results agree with data published in Colombia that was obtained from cord blood samples.

Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP) and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the Major Histocompatibility Complex region (OR = 0.74, 95%CI = 0.66-0.83, lowest p-value = 4.53x10-8). The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95%CI = 0.90-0.97, p-value = 3.54x10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the Human Leukocyte Antigen (HLA) genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils, and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Continental-scale genomic analysis suggests shared post-admixture adaptation in the Americas

American populations are one of the most interesting examples of recently admixed groups, where ancestral components from three major continental human groups (Africans, Eurasians and Native Americans) have admixed within the last 15 generations. Recently, several genetic surveys focusing on thousands of individuals shed light on the geography, chronology and relevance of these events. However, even though gene flow could drive adaptive evolution, it is unclear whether and how natural selection acted on the resulting genetic variation in the Americas.

In this study, we analysed the patterns of local ancestry of genomic fragments in genome-wide data for ~ 6000 admixed individuals from 10 American countries. In doing so, we identified regions characterized by a divergent ancestry profile (DAP), in which a significant over or under ancestral representation is evident.

Our results highlighted a series of genomic regions with DAPs associated with immune system response and relevant medical traits, with the longest DAP region encompassing the human leukocyte antigen locus. Furthermore, we found that DAP regions are enriched in genes linked to cancer-related traits and autoimmune diseases. Then, analysing the biological impact of these regions, we showed that natural selection could have acted preferentially towards variants located in coding and non-coding transcripts and characterized by a high deleteriousness score.

Taken together, our analyses suggest that shared patterns of post admixture adaptation occurred at a continental scale in the Americas, affecting more often functional and impactful genomic variants.

The Phenotypic Consequences of Genetic Divergence between Admixed Latin American Populations: Antioquia and Chocó, Colombia

Genome-wide association studies have uncovered thousands of genetic variants that are associated with a wide variety of human traits. Knowledge of how trait-associated variants are distributed within and between populations can provide insight into the genetic basis of group-specific phenotypic differences, particularly for health-related traits. We analyzed the genetic divergence levels for 1) individual trait-associated variants and 2) collections of variants that function together to encode polygenic traits, between two neighboring populations in Colombia that have distinct demographic profiles: Antioquia (Mestizo) and Chocó (Afro-Colombian). Genetic ancestry analysis showed 62% European, 32% Native American, and 6% African ancestry for Antioquia compared with 76% African, 10% European, and 14% Native American ancestry for Chocó, consistent with demography and previous results. Ancestry differences can confound cross-population comparison of polygenic risk scores (PRS); however, we did not find any systematic bias in PRS distributions for the two populations studied here, and population-specific differences in PRS were, for the most part, small and symmetrically distributed around zero. Both genetic differentiation at individual trait-associated single nucleotide polymorphisms and population-specific PRS differences between Antioquia and Chocó largely reflected anthropometric phenotypic differences that can be readily observed between the populations along with reported disease prevalence differences. Cases where population-specific differences in genetic risk did not align with observed trait (disease) prevalence point to the importance of environmental contributions to phenotypic variance, for both infectious and complex, common disease. The results reported here are distributed via a web-based platform for searching trait-associated variants and PRS divergence levels at http://map.chocogen.com (last accessed August 12, 2020).

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

Genetic substructure and complex demographic history of South African Bantu speakers

South Eastern Bantu-speaking (SEB) groups constitute more than 80% of the population in South Africa. Despite clear linguistic and geographic diversity, the genetic differences between these groups have not been systematically investigated. Based on genome-wide data of over 5000 individuals, representing eight major SEB groups, we provide strong evidence for fine-scale population structure that broadly aligns with geographic distribution and is also congruent with linguistic phylogeny (separation of Nguni, Sotho-Tswana and Tsonga speakers). Although differential Khoe-San admixture plays a key role, the structure persists after Khoe-San ancestry-masking. The timing of admixture, levels of sex-biased gene flow and population size dynamics also highlight differences in the demographic histories of individual groups. The comparisons with five Iron Age farmer genomes further support genetic continuity over ~400 years in certain regions of the country. Simulated trait genome-wide association studies further show that the observed population structure could have major implications for biomedical genomics research in South Africa.

Alpha thalassemia and alpha-MRE haplotypes in Uruguayan patients with microcytosis and hypochromia without anemia

Alpha thalassemia is the most common genetic disorder across the world, being the α-3.7 deletion the most frequent mutation. In order to analyze the spectrum and origin of alpha thalassemia mutations in Uruguay, we obtained a sample of 168 unrelated outpatients with normal hemoglobin levels with microcytosis and hypochromia from two cities: Montevideo and Salto. The presence of α-thalassemia mutations was investigated by gap-PCR, restriction endonucleases analysis and HBA2 and HBA1 genes sequencing, whereas the alpha-MRE haplotypes were investigated by sequencing. We found 55 individuals (32.7%) with α-thalassemia mutations, 51(30.4%) carrying the -α3.7 deletion, one with the -α4.2 deletion and three having the rare punctual mutation HBA2:c.-59C>T. Regarding alpha-MRE analysis, we observed a significant higher frequency of haplotype D, characteristic of African populations, in the sample with the -α3.7 deletion. These results show that α-thalassemia mutations are an important determinant of microcytosis and hypochromia in Uruguayan patients with microcytosis and hypochromia without anemia, mainly due to the -α3.7 deletion. The alpha-MRE haplotypes and the α-thalassemia mutations spectrum suggest a predominant, but not exclusive, African origin of these mutations in Uruguay.

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.

Rapid adaptation to malaria facilitated by admixture in the human population of Cabo Verde

Humans have undergone large migrations over the past hundreds to thousands of years, exposing ourselves to new environments and selective pressures. Yet, evidence of ongoing or recent selection in humans is difficult to detect. Many of these migrations also resulted in gene flow between previously separated populations. These recently admixed populations provide unique opportunities to study rapid evolution in humans. Developing methods based on distributions of local ancestry, we demonstrate that this sort of genetic exchange has facilitated detectable adaptation to a malaria parasite in the admixed population of Cabo Verde within the last ~20 generations. We estimate that the selection coefficient is approximately 0.08, one of the highest inferred in humans. Notably, we show that this strong selection at a single locus has likely affected patterns of ancestry genome-wide, potentially biasing demographic inference. Our study provides evidence of adaptation in a human population on historical timescales.

Association of variants in IL1B, TLR9, TREM1, IL10RA, and CD3G and Native American ancestry on malaria susceptibility in Colombian populations

Host genetics is an influencing factor in the manifestation of infectious diseases. In this study, the association of mild malaria with 28 variants in 16 genes previously reported in other populations and/or close to ancestry-informative markers (AIMs) selected was evaluated in an admixed 736 Colombian population sample. Additionally, the effect of genetic ancestry on phenotype expression was explored. For this purpose, the ancestral genetic composition of Turbo and El Bagre was determined. A higher Native American ancestry trend was found in the population with lower malaria susceptibility [odds ratio (OR) = 0.416, 95% confidence interval (95% CI) = 0.234-0.740, P = 0.003]. Three AIMs presented significant associations with the disease phenotype (MID1752, MID921, and MID1586). The first two were associated with greater malaria susceptibility (D/D, OR = 2.23, 95% CI = 1.06-4.69, P = 0.032 and I/D-I/I, OR = 2.14, 95% CI = 1.18-3.87, P = 0.011, respectively), and the latter has a protective effect on the appearance of malaria (I/I, OR = 0.18, 95% CI = 0.08-0.40, P < 0.0001). After adjustment by age, sex, municipality, and genetic ancestry, genotype association analysis showed evidence of association with malaria susceptibility for variants in or near IL1B, TLR9, TREM1, IL10RA, and CD3G genes: rs1143629-IL1B (G/A-A/A, OR = 0.41, 95% CI = 0.21-0.78, P = 0.0051), rs352139-TLR9 (T/T, OR = 0.28, 95% CI = 0.11-0.72, P = 0.0053), rs352140-TLR9 (C/C, OR = 0.41, 95% CI = 0.20-0.87, P = 0.019), rs2234237-TREM1 (T/A-A/A, OR = 0.43, 95% CI = 0.23-0.79, P = 0.0056), rs4252246-IL10RA (C/A-A/A, OR = 2.11, 95% CI = 1.18-3.75, P = 0.01), and rs1561966-CD3G (A/A, OR = 0.20, 95% CI = 0.06-0.69, P = 0.0058). The results showed the participation of genes involved in immunological processes and suggested an effect of ancestral genetic composition over the traits analyzed. Compared to the paisa population (Antioquia), Turbo and El Bagre showed a strong decrease in European ancestry and an increase in African and Native American ancestries. Also, a novel association of two single nucleotide polymorphisms with malaria susceptibility was identified in this study.

Canalization and Robustness in Human Genetics and Disease

Canalization refers to the evolution of populations such that the number of individuals who deviate from the optimum trait, or experience disease, is minimized. In the presence of rapid cultural, environmental, or genetic change, the reverse process of decanalization may contribute to observed increases in disease prevalence. This review starts by defining relevant concepts, drawing distinctions between the canalization of populations and robustness of individuals. It then considers evidence pertaining to three continuous traits and six domains of disease. In each case, existing genetic evidence for genotype-by-environment interactions is insufficient to support a strong inference of decanalization, but we argue that the advent of genome-wide polygenic risk assessment now makes an empirical evaluation of the role of canalization in preventing disease possible. Finally, the contributions of both rare and common variants to congenital abnormality and adult onset disease are considered in light of a new kerplunk model of genetic effects.

Minimal Erythema Dose: Correlation with Fitzpatrick Skin Type and Concordance Between Methods of Erythema Assessment in a Patient Sample in Colombia

BACKGROUND AND OBJECTIVE

The minimal erythema dose (MED), an essential measurement in studies of skin photosensitivity, requires establishing MED values for specific populations, given genetic variation. Different ways to assess erythema are also relevant. We aimed to determine MED values in a sample of Colombian patients and correlations between MED and Fitzpatrick skin type. We also studied concordance correlation between MEDs and two alternative ways to assess erythema.

PATIENTS AND METHODS

Cross-sectional study of 113 individuals in Bogotá, Colombia. We used a solar simulator to measure UV-A radiation and combined UV-A and UV-B (UVA+UVB) radiation, o se podría suprimir este término porque UVA y UVB son términos conocidos for MED calculation. Narrowband UV-B (NBUVB) radiation was measured in a phototherapy cabin. Erythema was assessed visually and with a Mexameter MX 18 device.

RESULTS

The median MEDs of UVA+UVB radiation were 22mJ/cm² for Fitzpatrick skin typesI andII, and 33 and 43mJ/cm², respectively, for typesIII andIV. The MEDs of UV-A radiation were 22, 42, 86, and 100J/cm² for typesI, II, III, andIV, respectively. The MEDs of NBUVB light were 390, 550, 770, and 885mJ/cm² for the 4 skin types. The correlation between MEDs and skin types ranged from 0.5 to 0.69. Lin's concordance correlation coefficients between visual and Mexameter assessments of erythema were greater than 0.8 in all cases.

CONCLUSION

This study allowed us to understand MED values for UV-A, UVA+UVB, and NBUVB according to different skin types in the Colombian population. Concordance correlation coefficients between the different methods of erythema assessment were very good. Correlations between MEDs and skin types were moderate to good.

Human Alveolar and Splenic Macrophage Populations Display a Distinct Transcriptomic Response to Infection With Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) infects alveolar macrophages (AMs), causing pulmonary tuberculosis (PTB), the most common form of the disease. Less frequently, Mtb is disseminated to many other organs and tissues, resulting in different extrapulmonary forms of TB. Nevertheless, very few studies have addressed the global mRNA response of human AMs, particularly from humans with the active form of the disease. Strikingly, almost no studies have addressed the response of human extrapulmonary macrophages to Mtb infection. In this pilot study, using microarray technology, we examined the transcriptomic ex vivo response of AMs from PTB patients (AMTBs) and AMs from control subjects (AMCTs) infected with two clinical isolates of Mtb. Furthermore, we also studied the infection response of human splenic macrophages (SMs) to Mtb isolates, as a model for extrapulmonary infection, and compared the transcriptomic response between AMs and SMs. Our results showed a striking difference in global mRNA profiles in response to infection between AMs and SMs, implicating a tissue-specific macrophage response to Mtb.

Development of a small panel of SNPs to infer ancestry in Chileans that distinguishes Aymara and Mapuche components

Background

Current South American populations trace their origins mainly to three continental ancestries, i.e. European, Amerindian and African. Individual variation in relative proportions of each of these ancestries may be confounded with socio-economic factors due to population stratification. Therefore, ancestry is a potential confounder variable that should be considered in epidemiologic studies and in public health plans. However, there are few studies that have assessed the ancestry of the current admixed Chilean population. This is partly due to the high cost of genome-scale technologies commonly used to estimate ancestry. In this study we have designed a small panel of SNPs to accurately assess ancestry in the largest sampling to date of the Chilean mestizo population (n = 3349) from eight cities. Our panel is also able to distinguish between the two main Amerindian components of Chileans: Aymara from the north and Mapuche from the south.

Results

A panel of 150 ancestry-informative markers (AIMs) of SNP type was selected to maximize ancestry informativeness and genome coverage. Of these, 147 were successfully genotyped by KASPar assays in 2843 samples, with an average missing rate of 0.012, and a 0.95 concordance with microarray data. The ancestries estimated with the panel of AIMs had relative high correlations (0.88 for European, 0.91 for Amerindian, 0.70 for Aymara, and 0.68 for Mapuche components) with those obtained with AXIOM LAT1 array. The country’s average ancestry was 0.53 ± 0.14 European, 0.04 ± 0.04 African, and 0.42 ± 0.14 Amerindian, disaggregated into 0.18 ± 0.15 Aymara and 0.25 ± 0.13 Mapuche. However, Mapuche ancestry was highest in the south (40.03%) and Aymara in the north (35.61%) as expected from the historical location of these ethnic groups. We make our results available through an online app and demonstrate how it can be used to adjust for ancestry when testing association between incidence of a disease and nongenetic risk factors.

Conclusions

We have conducted the most extensive sampling, across many different cities, of current Chilean population. Ancestry varied significantly by latitude and human development. The panel of AIMs is available to the community for estimating ancestry at low cost in Chileans and other populations with similar ancestry.

Evolutionary and population (epi)genetics of immunity to infection

Immune response is one of the functions that have been more strongly targeted by natural selection during human evolution. The evolutionary genetic dissection of the immune system has greatly helped to distinguish genes and functions that are essential, redundant or advantageous for human survival. It is also becoming increasingly clear that admixture between early Eurasians with now-extinct hominins such as Neanderthals or Denisovans, or admixture between modern human populations, can be beneficial for human adaptation to pathogen pressures. In this review, we discuss how the integration of population genetics with functional genomics in diverse human populations can inform about the changes in immune functions related to major lifestyle transitions (e.g., from hunting and gathering to farming), the action of natural selection to the evolution of the immune system, and the history of past epidemics. We also highlight the need of expanding the characterization of the immune system to a larger array of human populations-particularly neglected human groups historically exposed to different pathogen pressures-to fully capture the relative contribution of genetic, epigenetic, and environmental factors to immune response variation in humans.

Evaluating and sharing global genetic ancestry in biomedical datasets

Genetic ancestry is a critical co-factor to study phenotype-genotype associations using cohorts of human subjects. Most publicly available molecular datasets are, however, missing this information or only share self-reported race and ethnicity, representing a limitation to identify and repurpose datasets to investigate the contribution of ancestry to diseases and traits. We propose an analytical framework to enrich the metadata from publicly available cohorts with genetic ancestry information and a resulting diversity score at continental resolution, calculated directly from the data. We illustrate this framework using The Cancer Genome Atlas datasets searched through the DataMed Data Discovery Index. Data repositories and contributors can use this framework to provide genetic diversity measurements for controlled access datasets, minimizing the work involved in requesting a dataset that may ultimately prove inadequate for a researcher's purpose. With the increasing global scale of human genetics research, studies on disease risk and susceptibility would benefit greatly from the adequate estimation and sharing of genetic diversity in publicly available datasets following a framework such as the one presented.

Admixture-enabled selection for rapid adaptive evolution in the Americas

Background

Admixture occurs when previously isolated populations come together and exchange genetic material. We hypothesize that admixture can enable rapid adaptive evolution in human populations by introducing novel genetic variants (haplotypes) at intermediate frequencies, and we test this hypothesis through the analysis of whole genome sequences sampled from admixed Latin American populations in Colombia, Mexico, Peru, and Puerto Rico.

Results

Our screen for admixture-enabled selection relies on the identification of loci that contain more or less ancestry from a given source population than would be expected given the genome-wide ancestry frequencies. We employ a combined evidence approach to evaluate levels of ancestry enrichment at single loci across multiple populations and multiple loci that function together to encode polygenic traits. We find cross-population signals of African ancestry enrichment at the major histocompatibility locus on chromosome 6, consistent with admixture-enabled selection for enhanced adaptive immune response. Several of the human leukocyte antigen genes at this locus, such as HLA-A, HLA-DRB51, and HLA-DRB5, show independent evidence of positive selection prior to admixture, based on extended haplotype homozygosity in African populations. A number of traits related to inflammation, blood metabolites, and both the innate and adaptive immune system show evidence of admixture-enabled polygenic selection in Latin American populations.

Conclusions

The results reported here, considered together with the ubiquity of admixture in human evolution, suggest that admixture serves as a fundamental mechanism that drives rapid adaptive evolution in human populations.

Type 1 diabetes loci display a variety of native American and African ancestries in diseased individuals from Northwest Colombia

BACKGROUND

Type 1 diabetes (T1D) is a complex disease with a higher incidence in Europeans than other populations. The Colombians Living in Medellin (CLM) is admixed with ancestry contributions from Europeans, Native Americans (NAT) and Africans (AFR).

AIM

Our aim was to analyze the genetic admixture component at candidate T1D loci in Colombian individuals with the disease.

METHODS

Seventy-four ancestry informative markers (AIMs), which tagged 41 T1D candidate loci/genes, were tested by studying a cohort of 200 Northwest Colombia diseased individuals. T1D status was classified by testing for glutamic acid decarboxylase (GAD-65 kDa) and protein tyrosine-like antigen-2 auto-antibodies in serum samples. Candidate loci/genes included HLA, INS, PTPN22, CTLA4, IL2RA, SUMO4, CLEC16A, IFIH1, EFR3B, IL7R, NRP1 and RNASEH1, amongst others. The 1,000 genome database was used to analyze data from 94 individuals corresponding to the reference CLM. As the data did not comply with a normal distribution, medians were compared between groups using the Mann-Whitney U-test.

RESULTS

Both T1D patients and individuals from CLM displayed mainly European ancestry (61.58 vs 62.06) followed by Native American (27.34 vs 27.46) and to a lesser extent the AFR ancestry (10.28 vs 10.65) components. However, compared to CLM, ancestry of T1D patients displayed a decrease of NAT ancestry at gene EFR3B (24.30 vs 37.10) and an increase at genes IFIH1 (32.07 vs 14.99) and IL7R (52.18 vs 39.18). Also, for gene NRP1 (36.67 vs 0.003), we observed a non-AFR contribution (attributed to NAT). Autoimmune patients (positive for any of two auto-antibodies) displayed lower NAT ancestry than idiopathic patients at the MHC region (20.36 vs 31.88). Also, late onset patients presented with greater AFR ancestry than early onset patients at gene IL7R (19.96 vs 6.17). An association analysis showed that, even after adjusting for admixture, an association exists for at least seven such AIMs, with the strongest findings on chromosomes 5 and 10 (gene IL7R, P = 5.56 × 10^-6 and gene NRP1, P = 8.70 × 10^-19, respectively).

CONCLUSION

Although Colombian T1D patients have globally presented with higher European admixture, specific T1D loci have displayed varying levels of Native American and AFR ancestries in diseased individuals.

Native American admixture recapitulates population-specific migration and settlement of the continental United States

European and African descendants settled the continental US during the 17th-19th centuries, coming into contact with established Native American populations. The resulting admixture among these groups yielded a significant reservoir of Native American ancestry in the modern US population. We analyzed the patterns of Native American admixture seen for the three largest genetic ancestry groups in the US population: African descendants, Western European descendants, and Spanish descendants. The three groups show distinct Native American ancestry profiles, which are indicative of their historical patterns of migration and settlement across the country. Native American ancestry in the modern African descendant population does not coincide with local geography, instead forming a single group with origins in the southeastern US, consistent with the Great Migration of the early 20th century. Western European descendants show Native American ancestry that tracks their geographic origins across the US, indicative of ongoing contact during westward expansion, and Native American ancestry can resolve Spanish descendant individuals into distinct local groups formed by more recent migration from Mexico and Puerto Rico. We found an anomalous pattern of Native American ancestry from the US southwest, which most likely corresponds to the Nuevomexicano descendants of early Spanish settlers to the region. We addressed a number of controversies surrounding this population, including the extent of Sephardic Jewish ancestry. Nuevomexicanos are less admixed than nearby Mexican-American individuals, with more European and less Native American and African ancestry, and while they do show demonstrable Sephardic Jewish ancestry, the fraction is no greater than seen for other New World Spanish descendant populations.

Genetic variability analysis in a population from Bogotá: Towards a haplotype map

Introducción.

Los proyectos del mapa de haplotipos (HapMap) y de los 1.000 genomas han sido fundamentales para la compresión del componente genético de las enfermedades comunes y los fenotipos normales. Sin embargo, la variabilidad genética colombiana incluida en estos proyectos no es representativa del país.

Objetivo.

Contribuir al conocimiento de la variabilidad genética de la población colombiana a partir del estudio genómico de una muestra de individuos de Bogotá.

Materiales y métodos.

Se genotipificaron 2’372.784 marcadores genéticos de 32 individuos nacidos en Bogotá y de padres originarios de la misma ciudad utilizando la plataforma Illumina™. Los niveles de variabilidad genética se determinaron y se compararon con los datos disponibles de otras poblaciones del proyecto de los 1.000 genomas.

Resultados.

Los individuos analizados presentaron una variabilidad genética semejante a la de poblaciones con las que comparten ancestros. No obstante, a pesar de la poca diferenciación genética detectada en la población de Bogotá y en la de Medellín, el análisis de los componentes principales sugiere una composición genética diferente en las dos poblaciones.

Conclusiones.

El análisis genómico de la muestra de Bogotá permitió detectar similitudes y diferencias con otras poblaciones americanas. El aumento de tamaño de la muestra bogotana y la inclusión de muestras de otras regiones del país permitirán una mejor compresión de la variabilidad genética en Colombia, lo cual es fundamental para los estudios de salud humana, y la prevención y el tratamiento de enfermedades comunes en el país.

Reference values for spirometric parameters in healthy children living in a Colombian city located at 2640 m altitude

OBJECTIVES

To identify the spirometric equations that are most appropriate for use in children and adolescents living in Bogota, Colombia after evaluating a set of relevant previously-developed equations, including the Global Lung Function Initiative (GLI) 2012 spirometry reference equations.

METHODS

Healthy children aged between 6 and 17 years that were attending two randomly-selected schools in Bogota were invited to participate in the study, from January 2017 to January 2018. All participants underwent spirometry, following the procedures recommended by the American Thoracic Society/European Respiratory Society (ATS/ERS) Task Force. To identify the model or group of models that best predict each spirometric parameter in our population, we performed extensive residuals analyses and constructed Bland-Altman plots.

RESULTS

Three hundred twenty-six spirometric tests (149 boys, 177 girls) formed the reference data set. Knudson and GLI-2012 spirometry reference equations proved to be the most accurate in predicting the majority of spirometry parameters in both sexes and both age groups, each providing the lowest median prediction error in the residual analyses or the narrowest limits of agreement in the Bland-Altman plots in approximately one-third of the spirometry parameters analyzed.

CONCLUSION

For the majority of spirometry parameters, we recommend the use of Knudson and GLI-2012 spirometry reference equations for evaluating the respiratory function of children living in Bogota, Colombia, a city located at an altitude of 2640 m. Future investigations should target additional spirometric equations from Latin American populations living at moderate to high altitude to improve the GLI-2012 equations.

Assortative Mating on Ancestry-Variant Traits in Admixed Latin American Populations

Assortative mating is a universal feature of human societies, and individuals from ethnically diverse populations are known to mate assortatively based on similarities in genetic ancestry. However, little is currently known regarding the exact phenotypic cues, or their underlying genetic architecture, which inform ancestry-based assortative mating. We developed a novel approach, using genome-wide analysis of ancestry-specific haplotypes, to evaluate ancestry-based assortative mating on traits whose expression varies among the three continental population groups – African, European, and Native American – that admixed to form modern Latin American populations. Application of this method to genome sequences sampled from Colombia, Mexico, Peru, and Puerto Rico revealed widespread ancestry-based assortative mating. We discovered a number of anthropometric traits (body mass, height, and facial development) and neurological attributes (educational attainment and schizophrenia) that serve as phenotypic cues for ancestry-based assortative mating. Major histocompatibility complex (MHC) loci show population-specific patterns of both assortative and disassortative mating in Latin America. Ancestry-based assortative mating in the populations analyzed here appears to be driven primarily by African ancestry. This study serves as an example of how population genomic analyses can yield novel insights into human behavior.

Population Pharmacogenomics for Precision Public Health in Colombia

While genomic approaches to precision medicine hold great promise, they remain prohibitively expensive for developing countries. The precision public health paradigm, whereby healthcare decisions are made at the level of populations as opposed to individuals, provides one way for the genomics revolution to directly impact health outcomes in the developing world. Genomic approaches to precision public health require a deep understanding of local population genomics, which is still missing for many developing countries. We are investigating the population genomics of genetic variants that mediate drug response in an effort to inform healthcare decisions in Colombia. Our work focuses on two neighboring populations with distinct ancestry profiles: Antioquia and Chocó. Antioquia has primarily European genetic ancestry followed by Native American and African components, whereas Chocó shows mainly African ancestry with lower levels of Native American and European admixture. We performed a survey of the global distribution of pharmacogenomic variants followed by a more focused study of pharmacogenomic allele frequency differences between the two Colombian populations. Worldwide, we found pharmacogenomic variants to have both unusually high minor allele frequencies and high levels of population differentiation. A number of these pharmacogenomic variants also show anomalous effect allele frequencies within and between the two Colombian populations, and these differences were found to be associated with their distinct genetic ancestry profiles. For example, the C allele of the single nucleotide polymorphism (SNP) rs4149056 [Solute Carrier Organic Anion Transporter Family Member 1B1 (SLCO1B1)^∗5], which is associated with an increased risk of toxicity to a commonly prescribed statin, is found at relatively high frequency in Antioquia and is associated with European ancestry. In addition to pharmacogenomic alleles related to increased toxicity risk, we also have evidence that alleles related to dosage and metabolism have large frequency differences between the two populations, which are associated with their specific ancestries. Using these findings, we have developed and validated an inexpensive allele-specific PCR assay to test for the presence of such population-enriched pharmacogenomic SNPs in Colombia. These results serve as an example of how population-centered approaches to pharmacogenomics can help to realize the promise of precision medicine in resource-limited settings.

Genetic origin of a large family with a novel PSEN1 mutation (Ile416Thr)

INTRODUCTION

A small percentage of Alzheimer's disease (AD) cases are caused by genetic mutations with autosomal dominant inheritance. We report a family with a novel variant in PSEN1.

METHODS

We performed clinical and genetic evaluation of 93 related individuals from a Colombian admixed population. 31 individuals had whole-genome sequencing.

RESULTS

Genetic analysis revealed a missense variant in PSEN1 (NM_000021.3: c.1247T>C p.Ile416Thr), which originated on an African haplotype and segregated with AD logarithm of the odds score of 6. Their clinical phenotype is similar to sporadic AD except for earlier age at onset: the mean age at onset for mild cognitive impairment was 47.6 years (standard deviation 5.83) and for dementia 51.6 years (standard deviation 5.03).

DISCUSSION

Ile416Thr is a novel pathogenic variant that causes AD in the sixth decade of life. The history of the region that included slave importation and admixtures within a confined geographic locale represents a "mini-population bottleneck" and subsequent emergence of a rare dominant mutation.

Adaptation to milking agropastoralism in Chilean goat herders and nutritional benefit of lactase persistence

The genetic trait of lactase persistence (LP) evolved as an adaptation to milking pastoralism in the Old World and is a well-known example of positive natural selection in humans. However, the specific mechanisms conferring this selective advantage are unknown. To understand the relationship between milk drinking, LP, growth, reproduction, and survival, communities of the Coquimbo Region in Chile, with recent adoption of milking agropastoralism, were used as a model population. DNA samples and data on stature, reproduction, and diet were collected from 451 participants. Lactose tolerance tests were done on 41 of them. The European -13,910*T (rs4988235) was the only LP causative variant found, showing strong association (99.6%) with LP phenotype. Models of associations of inferred LP status and milk consumption, with fertility, mortality, height, and weight were adjusted with measures of ancestry and relatedness to control for population structure. Although we found no statistically significant effect of LP on fertility, a significant effect (P = 0.002) was observed of LP on body mass index (BMI) in males and of BMI on fertility (P = 0.003). These results fail to support a causal relationship between LP and fertility yet suggest the idea of a nutritional advantage of LP. Furthermore, the proportion of European ancestry around the genetic region of -13,910*T is significantly higher (P = 0.008) than the proportion of European ancestry genome-wide, providing evidence of recent positive selection since European-Amerindian admixture. This signature was absent in nonpastoralist Latin American populations, supporting the hypothesis of specific adaptation to milking agropastoralism in the Coquimbo communities.

Genetic ancestry, admixture and health determinants in Latin America

BACKGROUND

Modern Latin American populations were formed via genetic admixture among ancestral source populations from Africa, the Americas and Europe. We are interested in studying how combinations of genetic ancestry in admixed Latin American populations may impact genomic determinants of health and disease. For this study, we characterized the impact of ancestry and admixture on genetic variants that underlie health- and disease-related phenotypes in population genomic samples from Colombia, Mexico, Peru, and Puerto Rico.

RESULTS

We analyzed a total of 347 admixed Latin American genomes along with 1102 putative ancestral source genomes from Africans, Europeans, and Native Americans. We characterized the genetic ancestry, relatedness, and admixture patterns for each of the admixed Latin American genomes, finding a spectrum of ancestry proportions within and between populations. We then identified single nucleotide polymorphisms (SNPs) with anomalous ancestry-enrichment patterns, i.e. SNPs that exist in any given Latin American population at a higher frequency than expected based on the population's genetic ancestry profile. For this set of ancestry-enriched SNPs, we inspected their phenotypic impact on disease, metabolism, and the immune system. All four of the Latin American populations show ancestry-enrichment for a number of shared pathways, yielding evidence of similar selection pressures on these populations during their evolution. For example, all four populations show ancestry-enriched SNPs in multiple genes from immune system pathways, such as the cytokine receptor interaction, T cell receptor signaling, and antigen presentation pathways. We also found SNPs with excess African or European ancestry that are associated with ancestry-specific gene expression patterns and play crucial roles in the immune system and infectious disease responses. Genes from both the innate and adaptive immune system were found to be regulated by ancestry-enriched SNPs with population-specific regulatory effects.

CONCLUSIONS

Ancestry-enriched SNPs in Latin American populations have a substantial effect on health- and disease-related phenotypes. The concordant impact observed for same phenotypes across populations points to a process of adaptive introgression, whereby ancestry-enriched SNPs with specific functional utility appear to have been retained in modern populations by virtue of their effects on health and fitness.