Simultaneous inference of parental admixture proportions and admixture times from unphased local ancestry calls

Population genetic analyses of local ancestry tracts routinely assume that the ancestral admixture process is identical for both parents of an individual, an assumption that may be invalid when considering recent admixture. Here we present Parental Admixture Proportion Inference (PAPI), a Bayesian tool for inferring the admixture proportions and admixture times for each parent of a single admixed individual. PAPI analyzes unphased local ancestry tracts and has two components models: a binomial model that exploits the informativeness of homozygous ancestry regions to infer parental admixture proportions, and a hidden Markov model (HMM) that infers admixture times from tract lengths. Crucially, the HMM employs an approximation to the pedigree crossover dynamics that accounts for unobserved within-ancestry recombination, enabling inference of parental admixture times. We compared the accuracy of PAPI's admixture proportion estimates with those of ANCESTOR in simulated admixed individuals and found that PAPI outperforms ANCESTOR by an average of 46% in a representative set of simulation scenarios, with PAPI's estimates deviating from the ground truth by 0.047 on average. Moreover, PAPI's admixture time estimates were strongly correlated with the ground truth in these simulations (R = 0.76), but have an average downward bias of 1.01 generations that is partly attributable to inaccuracies in local ancestry inference. As an illustration of its utility, we ran PAPI on real African Americans from the PAGE study (N = 5,786) and found strong evidence of assortative mating by ancestry proportion: couples' ancestry proportions are closer to each other than expected by chance (P<10-6), and are highly correlated (R = 0.87). We anticipate that PAPI will be useful in studying the population dynamics of admixture and will also be of interest to individuals seeking to learn about their personal genealogies.

Dysregulation of DNA Methylation and Epigenetic Clocks in Prostate Cancer among Puerto Rican Men

In 2021, approximately 248,530 new prostate cancer (PCa) cases are estimated in the United States. Hispanic/Latinos (H/L) are the second largest racial/ethnic group in the US. The objective of this study was to assess DNA methylation patterns between aggressive and indolent PCa along with ancestry proportions in 49 H/L men from Puerto Rico (PR). Prostate tumors were classified as aggressive (n = 17) and indolent (n = 32) based on the Gleason score. Genomic DNA samples were extracted by macro-dissection. DNA methylation patterns were assessed using the Illumina EPIC DNA methylation platform. We used ADMIXTURE to estimate global ancestry proportions. We identified 892 differentially methylated genes in prostate tumor tissues as compared with normal tissues. Based on an epigenetic clock model, we observed that the total number of stem cell divisions (TNSC) and stem cell division rate (SCDR) were significantly higher in tumor than adjacent normal tissues. Regarding PCa aggressiveness, 141 differentially methylated genes were identified. Ancestry proportions of PR men were estimated as African, European, and Indigenous American; these were 24.1%, 64.2%, and 11.7%, respectively. The identification of DNA methylation profiles associated with risk and aggressiveness of PCa in PR H/L men will shed light on potential mechanisms contributing to PCa disparities in PR population.

Development and Evaluation of the Ancestry Informative Marker Panel of the VISAGE Basic Tool

We detail the development of the ancestry informative single nucleotide polymorphisms (SNPs) panel forming part of the VISAGE Basic Tool (BT), which combines 41 appearance predictive SNPs and 112 ancestry predictive SNPs (three SNPs shared between sets) in one massively parallel sequencing (MPS) multiplex, whereas blood-based age analysis using methylation markers is run in a parallel MPS analysis pipeline. The selection of SNPs for the BT ancestry panel focused on established forensic markers that already have a proven track record of good sequencing performance in MPS, and the overall SNP multiplex scale closely matched that of existing forensic MPS assays. SNPs were chosen to differentiate individuals from the five main continental population groups of Africa, Europe, East Asia, America, and Oceania, extended to include differentiation of individuals from South Asia. From analysis of 1000 Genomes and HGDP-CEPH samples from these six population groups, the BT ancestry panel was shown to have no classification error using the Bayes likelihood calculators of the Snipper online analysis portal. The differentiation power of the component ancestry SNPs of BT was balanced as far as possible to avoid bias in the estimation of co-ancestry proportions in individuals with admixed backgrounds. The balancing process led to very similar cumulative population-specific divergence values for Africa, Europe, America, and Oceania, with East Asia being slightly below average, and South Asia an outlier from the other groups. Comparisons were made of the African, European, and Native American estimated co-ancestry proportions in the six admixed 1000 Genomes populations, using the BT ancestry panel SNPs and 572,000 Affymetrix Human Origins array SNPs. Very similar co-ancestry proportions were observed down to a minimum value of 10%, below which, low-level co-ancestry was not always reliably detected by BT SNPs. The Snipper analysis portal provides a comprehensive population dataset for the BT ancestry panel SNPs, comprising a 520-sample standardised reference dataset; 3445 additional samples from 1000 Genomes, HGDP-CEPH, Simons Foundation and Estonian Biocentre genome diversity projects; and 167 samples of six populations from in-house genotyping of individuals from Middle East, North and East African regions complementing those of the sampling regimes of the other diversity projects.

Abstract P075: Genetic European Ancestry And Incident Diabetes Mellitus In Black Individuals-Insights From The SPRINT Trial

Background: A metabolic paradox of lower triglyceride (TG) and higher high-density lipoprotein cholesterol (HDL-C) levels but a higher incidence of diabetes has been described in Blacks. We evaluated the association of genetic ancestry-related variants with the risk of incident diabetes in self-identified Black individuals in the Systolic Blood Pressure (BP) Intervention Trial (SPRINT). Methods: The genetic European ancestry proportions were estimated using 106 biallelic genotype markers in non-diabetic Black participants. Participants were stratified into tertiles (T1-T3) of European ancestry proportions. Multivariable-adjusted association of baseline metabolic syndrome indices (fasting plasma glucose [FPG], TG, HDL-C, body mass index [BMI], and BP) with the ancestry proportion was assessed. Multivariable-adjusted Cox regression was used to assess the risk of incident diabetes (FPG ≥126 mg/dL or self-reported diabetes treatment) across tertiles. Results: In 2,466 participants with median European ancestry 19% (13 - 27%), those in T1 were relatively younger, with lower Framingham risk score, and worse renal function (p<0.05 for all). At baseline, higher European ancestry proportion was associated with higher TG and lower HDL-C levels after controlling for clinical and demographic factors (p<0.05 for both). FPG, BMI, and BP were not significantly associated with ancestry proportion. Compared with T1, those in second (HR: 0.64 [95% CI: 0.45-0.90]) and third tertiles (HR: 0.61 [0.44-0.89]) had a lower risk of incident diabetes after controlling for covariates. Specifically, a 5% increment in European ancestry was associated with a 29% lower risk of incident diabetes (HR: 0.71 [0.55-0.93]). No interaction was observed between intensive BP control strategy and European ancestry proportion tertiles on incident diabetes. Conclusions: Genetic factors may account for racial differences in TG and HDL-C levels. The higher risk of incident diabetes in Blacks may similarly be genetically determined.

Evidence That Geographic Variation in Genetic Ancestry Associates With Uterine Fibroids

Uterine fibroids disproportionately impact African American (AA) women. Evidence suggests AA women have earlier onset and higher cumulative risk. This risk disparity may be due an imbalance of risk alleles in one parental geographic subpopulation relative to others. We investigated ancestry proportions for the 1000 Genomes phase 3 populations clustered into 6 geographic groups for association with fibroid traits in AA women (n=583 cases, 797 controls) and European American (EA) women (n=1,195 cases, 1,164 controls). Global ancestry proportions were estimated using ADMIXTURE. Dichotomous (fibroids status and multiple fibroid status) and continuous outcomes (volume and largest dimension) were modeled for association with ancestry proportions using logistic and linear regression adjusting for age. Effect estimates are reported per 10% increase in genetically inferred ancestry proportion. Among AAs, West African (WAFR) ancestry was associated with fibroid risk, East African ancestry was associated with risk of multiple fibroids, Northern European (NEUR) ancestry was protective for multiple fibroids, Southern European ancestry was protective for fibroids and multiple fibroids, and South Asian (SAS) ancestry was positively associated with volume and largest dimension. In EAs, NEUR ancestry was protective for fibroids, SAS ancestry was associated with fibroid risk, and WAFR ancestry was positively associated with volume and largest dimension. These results suggest that a proportion of fibroid risk and fibroid trait racial disparities are due to genetic differences between geographic groups. Further investigation at the local ancestry and single variant levels may yield novel insights about disease architecture and genetic mechanisms underlying ethnic disparities in fibroid risk.

“Immunogenetics of resistance to SARS-CoV-2 infection in discordant couples”

BackgroundDespite the high number of individuals infected by SARS-CoV-2 who develop COVID-19 symptoms worldwide, many exposed individuals remain asymptomatic and/or stay uninfected. This could be explained by a combination of environmental (exposure, previous infection), epigenetic, and genetic factors. Aiming to identify genetic variants involved in SARS-CoV-2 resistance, we analyzed 86 discordant Brazilian couples where one was infected and symptomatic while the partner remained asymptomatic and seronegative despite sharing the same bedroom during the infection. The discordant partners had comparable ages, and genetic ancestry proportions.MethodsWhole-exome sequencing followed by a state-of-the-art method to call genotypes and haplotypes across the highly polymorphic MHC and LRC.ResultsWe observed a minor impact in antigen-presentation genes and KIR genes associated with resistance. Interestingly, genes related to immune modulation, mainly involved in NK cell killing activation/inhibition harbor variants potentially contributing to infection resistance. We hypothesize that individuals prone to produce higher amounts of MICA (possibly soluble), LILRB1, LILRB2, and low amounts of MICB, would be more susceptible to infection.ConclusionAccording to this hypothesis, quantitative differences in these NK activity-related molecules could contribute to resistance to COVID-19 down regulating NK cell cytotoxic activity in infected individuals but not in resistant partners.

Summix: A method for detecting and adjusting for population structure in genetic summary data

Publicly available genetic summary data have high utility in research and the clinic including prioritizing putative causal variants, polygenic scoring, and leveraging common controls. However, summarizing individual-level data can mask population structure resulting in confounding, reduced power, and incorrect prioritization of putative causal variants. This limits the utility of publicly available data, especially for understudied or admixed populations where additional research and resources are most needed. While several methods exist to estimate ancestry in individual-level data, methods to estimate ancestry proportions in summary data are lacking. Here, we present Summix, a method to efficiently deconvolute ancestry and provide ancestry-adjusted allele frequencies from summary data. Using continental reference ancestry, African (AFR), Non-Finnish European (EUR), East Asian (EAS), Indigenous American (IAM), South Asian (SAS), we obtain accurate and precise estimates (within 0.1%) for all simulation scenarios. We apply Summix to gnomAD v2.1 exome and genome groups and subgroups finding heterogeneous continental ancestry for several groups including African/African American (∼84% AFR, ∼14% EUR) and American/Latinx (∼4% AFR, ∼5% EAS, ∼43% EUR, ∼46% IAM). Compared to the unadjusted gnomAD AFs, Summix’s ancestry-adjusted AFs more closely match respective African and Latinx reference samples. Even on modern, dense panels of summary statistics, Summix yields results in seconds allowing for estimation of confidence intervals via block bootstrap. Given an accompanying R package, Summix increases the utility and equity of public genetic resources, empowering novel research opportunities.

Identification of ancestry proportions in admixed groups across the Americas using clinical pharmacogenomic SNP panels

We evaluated the performance of three PGx panels to estimate biogeographical ancestry: the DMET panel, and the VIP and Preemptive PGx panels described in the literature. Our analysis indicate that the three panels capture quite well the individual variation in admixture proportions observed in recently admixed populations throughout the Americas, with the Preemptive PGx and DMET panels performing better than the VIP panel. We show that these panels provide reliable information about biogeographic ancestry and can be used to guide the implementation of PGx clinical decision-support (CDS) tools. We also report that using these panels it is possible to control for the effects of population stratification in association studies in recently admixed populations, as exemplified with a warfarin dosing GWA study in a sample from Brazil.

Differential DNA Methylation in Prostate Tumors from Puerto Rican Men

In 2020, approximately 191,930 new prostate cancer (PCa) cases are estimated in the United States (US). Hispanic/Latinos (H/L) are the second largest racial/ethnic group in the US. This study aims to assess methylation patterns between aggressive and indolent PCa including DNA repair genes along with ancestry proportions. Prostate tumors classified as aggressive (n = 11) and indolent (n = 13) on the basis of the Gleason score were collected. Tumor and adjacent normal tissue were annotated on H&E (Haemotoxylin and Eosin) slides and extracted by macro-dissection. Methylation patterns were assessed using the Illumina 850K DNA methylation platform. Raw data were processed using the Bioconductor package. Global ancestry proportions were estimated using ADMIXTURE (k = 3). One hundred eight genes including AOX1 were differentially methylated in tumor samples. Regarding the PCa aggressiveness, six hypermethylated genes (RREB1, FAM71F2, JMJD1C, COL5A3, RAE1, and GABRQ) and 11 hypomethylated genes (COL9A2, FAM179A, SLC17A2, PDE10A, PLEKHS1, TNNI2, OR51A4, RNF169, SPNS2, ADAMTSL5, and CYP4F12) were identified. Two significant differentially methylated DNA repair genes, JMJD1C and RNF169, were found. Ancestry proportion results for African, European, and Indigenous American were 24.1%, 64.2%, and 11.7%, respectively. The identification of DNA methylation patterns related to PCa in H/L men along with specific patterns related to aggressiveness and DNA repair constitutes a pivotal effort for the understanding of PCa in this population.

Inference of recent admixture using genotype data

The inference of biogeographic ancestry (BGA) has become a focus of forensic genetics. Mis-inference of BGA can have profound unwanted consequences for investigations and society. We show that recent admixture can lead to misclassification and erroneous inference of ancestry proportions, using state of the art analysis tools with (i) simulations, (ii) 1000 genomes project data, and (iii) two individuals analyzed using the ForenSeq DNA Signature Prep Kit. Subsequently, we extend existing tools for estimation of individual ancestry (IA) by allowing for different IA in both parents, leading to estimates of parental individual ancestry (PIA), and a statistical test for recent admixture. Estimation of PIA outperforms IA in most scenarios of recent admixture. Furthermore, additional information about parental ancestry can be acquired with PIA that may guide casework.