Admixed Ancestry Research Articles

Leveraging diverse genomic data to guide equitable carrier screening: Insights from gnomAD v.4.1.0.

Analysis of exome data from the latest release of the Genome Aggregation Database (gnomAD v.4.1.0) revealed a significant carrier burden of pathogenic/likely pathogenic (P/LP) variants in genes associated with autosomal-recessive conditions across diverse ancestral populations. Carrier screening panels are routinely offered to reproductive partners to inform their risk of having an affected child. Current guidelines from the American College of Medical Genetics and Genomics (ACMG) recommend screening for genes with a carrier frequency of at least 1/200 and associated with moderate/severe conditions. Here, we systematically analyzed >700,000 gnomAD v.4.1.0 exomes spanning eight ancestries to estimate the carrier frequency of P/LP variants in 2,987 genes associated with autosomal-recessive conditions. After expert curation for clinical severity, we identified 286 genes meeting the criteria for carrier screening. The number of genes exceeding the 1/200 threshold varied across populations, with 40 in the South Asian ancestry and up to 119 in the Ashkenazi Jewish population. Simulations showed that pan-ethnic screening panels offer advantages for individuals of diverse or admixed ancestry, while ancestry-specific panels may be preferable for genetically homogeneous populations. This study leveraged the most comprehensive genomic dataset to date to provide an updated candidate gene list for equitable carrier screening across diverse populations. Our findings highlight the need for continued expansion of genomic resources to better understand rare disease risk and inform screening efforts in underrepresented groups.

Biobank-scale characterization of Alzheimer's disease and related dementias identifies potential disease-causing variants, risk factors, and genetic modifiers across diverse ancestries.

Alzheimer's disease and related dementias (AD/ADRDs) pose a significant global public health challenge, underscored by the intricate interplay of genetic and environmental factors that differ across ancestries. To effectively implement equitable, personalized therapeutic interventions on a global scale, it is essential to identify disease-causing mutations and genetic risk and resilience factors across diverse ancestral backgrounds. Exploring genetic-phenotypic correlations across the globe enhances the generalizability of research findings, contributing to a more inclusive and universal understanding of disease. This study leveraged biobank-scale data to conduct the largest multi-ancestry whole-genome sequencing characterization of AD/ADRDs. We aimed to build a valuable catalog of potential disease-causing, genetic risk and resilience variants impacting the etiology of these conditions. We thoroughly characterized genetic variants from key genes associated with AD/ADRDs across 11 genetic ancestries, utilizing data from All of Us, UK Biobank, 100,000 Genomes Project, Alzheimer's Disease Sequencing Project, and the Accelerating Medicines Partnership in Parkinson's Disease, including a total of 25,001 cases and 93,542 controls. We prioritized 116 variants possibly linked to disease, including 18 known pathogenic and 98 novel variants. We detected previously described disease-causing variants among controls, leading us to question their pathogenicity. Notably, we showed a higher frequency of APOE ε4/ε4 carriers among individuals of African and African Admixed ancestry compared to other ancestries, confirming ancestry-driven modulation of APOE -associated AD/ADRDs. A thorough assessment of APOE revealed a disease-modifying effect conferred by the TOMM40 :rs11556505, APOE :rs449647, 19q13.31 :rs10423769, NOCT :rs13116075, CASS4 :rs6024870, and LRRC37A :rs2732703 variants among APOE ε4 carriers across different ancestries. In summary, we compiled the most extensive catalog of established and novel genetic variants in known genes increasing risk or conferring resistance to AD/ADRDs across diverse ancestries, providing clinical insights into their genetic-phenotypic correlations. The findings from this investigation hold significant implications for potential clinical trials and therapeutic interventions on a global scale. Finally, we present an accessible and user-friendly platform for the AD/ADRDs research community to help inform and support basic, translational, and clinical research on these debilitating conditions ( https://niacard.shinyapps.io/MAMBARD_browser/ ).

Prevalence of germline variants in Brazilian pancreatic carcinoma patients

We evaluated the prevalence of pathogenic/likely pathogenic germline variants (PGV) in Brazilian pancreatic adenocarcinoma (PC) patients, that represent a multiethnic population, in a cross-sectional study. We included 192 PC patients unselected for family history of cancer. We evaluated a panel of 113 cancer genes, through genomic DNA sequencing and 46 ancestry-informative markers, through multiplex PCR. The median age was 61 years; 63.5% of the patients presented disease clinical stages III or IV; 8.3% reported personal history of cancer; 4.7% and 16.1% reported first-degree relatives with PC or breast and/or prostate cancer, respectively. Although the main ancestry was European, there was considerable genetic composition admixture. Twelve patients (6.25%) were PGV carriers in PC predisposition genes (ATM, BRCA1, BRCA2, CDKN2A, MSH2, PALB2) and another 25 (13.0%) were PGV carriers in genes with a limited association or not previously associated with PC (ACD, BLM, BRIP1, CHEK2, ERCC4, FANCA, FANCE, FANCM, GALNT12, MITF, MRE11, MUTYH, POLE, RAD51B, RAD51C, RECQL4, SDHA, TERF2IP). The most frequently affected genes were CHEK2, ATM and FANC. In tumor samples from PGV carriers in ACD, BRIP1, MRE11, POLE, SDHA, TERF2IP, which were examined through exome sequencing, the main single base substitutions (SBS) mutational signature was SBS1+5+18, probably associated with age, tobacco smoking and reactive oxygen species. SBS3 associated with homologous repair deficiency was also represented, but on a lower scale. There was no difference in the frequency of PGV carriers between: (a) patients with or without first-degree relatives with cancer; and (b) patients with admixed ancestry versus those with predominantly European ancestry. Furthermore, there was no difference in overall survival between PGV carriers and non-carriers. Therefore, genetic testing should be offered to all Brazilian pancreatic cancer patients, regardless of their ancestry. Genes with limited or previously unrecognized associations with pancreatic cancer should be further investigated to clarify their role in cancer risk.

Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species.

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Complex trait susceptibilities and population diversity in a sample of 4,145 Russians

The population of Russia consists of more than 150 local ethnicities. The ethnic diversity and geographic origins, which extend from eastern Europe to Asia, make the population uniquely positioned to investigate the shared properties of inherited disease risks between European and Asian ancestries. We present the analysis of genetic and phenotypic data from a cohort of 4,145 individuals collected in three metro areas in western Russia. We show the presence of multiple admixed genetic ancestry clusters spanning from primarily European to Asian and high identity-by-descent sharing with the Finnish population. As a result, there was notable enrichment of Finnish-specific variants in Russia. We illustrate the utility of Russian-descent cohorts for discovery of novel population-specific genetic associations, as well as replication of previously identified associations that were thought to be population-specific in other cohorts. Finally, we provide access to a database of allele frequencies and GWAS results for 464 phenotypes.

Assessing the Risk Stratification of Breast Cancer Polygenic Risk Scores in a Brazilian Cohort

Polygenic risk scores (PRSs) for breast cancer have a clear clinical utility in risk prediction. PRS transferability across populations and ancestry groups is hampered by population-specific factors, ultimately leading to differences in variant effects, such as linkage disequilibrium and differences in variant frequency (allele frequency differences). Thus, locally sourced population-based phenotypic and genomic data sets are essential to assess the validity of PRSs derived from signals detected across populations. This study assesses the transferability of a breast cancer PRS composed of 313 risk variants (313-PRS) in a Brazilian trihybrid admixed ancestries (European, African, and Native American) whole-genome sequenced cohort, the Rare Genomes Project. 313-PRS was computed in the Rare Genomes Project (n=853) using the UK Biobank (UKBB; n=264,307) as reference. The Brazilian cohorts have a high European ancestry (EA) component, with allele frequency differences and to a lesser extent linkage disequilibrium patterns similar to those found in EA populations. The 313-PRS distribution was found to be inflated when compared with that of the UKBB, leading to potential overestimation of PRS-based risk if EA is taken as a standard. However, case controls lead to equivalent predictive power when compared with UKBB-EA samples with area under the receiver operating characteristic curve values of 0.66 to 0.62 compared with 0.63 for UKBB.

Genome-wide study of half a million individuals with major depression identifies 697 independent associations, infers causal neuronal subtypes and biological targets for novel pharmacotherapies.

In a genome-wide association study (GWAS) of 685,808 individuals with major depression (MD) and 4,364,225 controls from 29 countries and across diverse and admixed ancestries, we identify 697 independent associations at 636 genetic loci, 293 of which are novel. Using fine-mapping and functional genomic datasets, we find 308 high-confidence gene associations and enrichment of postsynaptic density and receptor clustering. Leveraging new single-cell gene expression data, we conducted a causal neural cell type enrichment analysis that implicated excitatory and inhibitory midbrain and forebrain neurons, peptidergic neurons, and medium spiny neurons in MD. Critically, our findings are enriched for the targets of antidepressants and provide potential antidepressant repurposing opportunities (e.g., pregabalin and modafinil). Polygenic scores (PGS) from European ancestries explained up to 5.7% of the variance in liability to MD in European samples and PGS trained using either European or multi-ancestry data significantly predicted case control status across all four diverse ancestries. These findings represent a major advance in our understanding of MD across global populations. We provide evidence that MD GWAS reveals known and novel biological targets that may be used to target and develop pharmacotherapies addressing the considerable unmet need for effective treatment.

Does nuclear DNA support the recognition of three species within the Splendid Fairywren Malurus splendens?

ABSTRACT The Splendid Fairywren Malurus splendens has four distinct phenotypic forms currently recognised as subspecies, although three species have been recognised within the group in the past. Here, we use multilocus nuclear DNA to test the recent hypothesis that the Splendid Fairywren represents three allopatric taxa with little to no modern-day contact, despite phenotypic evidence of periods of past contact and gene flow. Our nuclear dataset is concordant with earlier mtDNA-based data in supporting three distinct genetic groups within M. splendens. These genetic groups align with the three phenotypically distinct forms recently hypothesised to be allopatric and previously described at species-level – western M. splendens, central M. callainus and eastern M. melanotus. Nuclear gene flow does not appear to be rampant between the three forms; however, several individuals show signals of admixed nuclear ancestry, which we cannot distinguish as having resulted from incomplete lineage sorting, poor phylogenetic signal, or gene flow. Given changing modern taxonomic interpretations of the significance of gene flow in speciation, we argue that recognising three distinct species within the M. splendens complex might better represent their evolutionary distinctiveness, rather than continuing to treat these genetically and morphologically distinct forms as subspecies of a widespread polytypic species. Our study highlights the interesting complexity of determining species status when taxa are in the ‘grey zone’ near the subspecies/species continuum, and the value gained from adding perspectives from nuclear DNA.

Data Harmonization Guidelines to Combine Multi-platform Genomic Data from Admixed Populations and Boost Power in Genome-Wide Association Studies.

Data harmonization involves combining data from multiple independent sources and processing the data to produce one uniform dataset. Merging separate genotypes or whole-genome sequencing datasets has been proposed as a strategy to increase the statistical power of association tests by increasing the effective sample size. However, data harmonization is not a widely adopted strategy due to the difficulties with merging data (including confounding produced by batch effects and population stratification). Detailed data harmonization protocols are scarce and are often conflicting. Moreover, data harmonization protocols that accommodate samples of admixed ancestry are practically non-existent. Existing data harmonization procedures must be modified to ensure the heterogeneous ancestry of admixed individuals is incorporated into additional downstream analyses without confounding results. Here, we propose a set of guidelines for merging multi-platform genetic data from admixed samples that can be adopted by any investigator with elementary bioinformatics experience. We have applied these guidelines to aggregate 1544 tuberculosis (TB) case-control samples from six separate in-house datasets and conducted a genome-wide association study (GWAS) of TB susceptibility. The GWAS performed on the merged dataset had improved power over analyzing the datasets individually and produced summary statistics free from bias introduced by batch effects and population stratification. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Processing separate datasets comprising array genotype data Alternate Protocol 1: Processing separate datasets comprising array genotype and whole-genome sequencing data Alternate Protocol 2: Performing imputation using a local reference panel Basic Protocol 2: Merging separate datasets Basic Protocol 3: Ancestry inference using ADMIXTURE and RFMix Basic Protocol 4: Batch effect correction using pseudo-case-control comparisons.

Genetic relationships and introgression events between wild and cultivated grapevines (Vitis vinifera L.): focus on Italian Lambruscos

Research efforts on genomic structure and ecology of wild populations of Vitis vinifera L. offer insights on grape domestication processes and on the assortment evolution of the cultivated forms. Attention is also paid to the origin of traditional, long-cultivated varieties, often producing renowned and valuable wines. The genetic relationships between 283 Vitis vinifera cultivated varieties (subsp. sativa) and 65 individuals from 9 populations of the sylvestris subspecies mainly from northern Italy were explored by means of molecular markers (27 nuclear and 4 chloroplastic microsatellites). Several episodes of contamination of the wild germplasm by the pollen of specific grape cultivars were detected, implying concern for maintaining the purity of the wild form. At the same time, events of introgression from the wild subspecies resulted playing a crucial role in the emergence of several cultivated varieties with a clear admixed genome ancestry sativa-sylvestris. These included Lambruscos originated from the flat areas crossed by the Po and Adige rivers in northern Italy, while other cultivars still called Lambrusco but typical of hilly areas did not show the same admixed genome. Historical and ecological evidences suggesting an adaptative recent post-domestication process in the origin of several Italian Lambruscos are discussed.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic β-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.

Identifying Rare Genetic Determinants for Improved Polygenic Risk Prediction of Bone Mineral Density and Fracture Risk.

Osteoporosis and fractures severely impact the elderly population. Polygenic risk scores for bone mineral density have demonstrated potential clinical utility. However, the value of rare genetic determinants in risk prediction has not been assessed. With whole-exome sequencing data from 436,824 UK Biobank participants, we assigned White British ancestry individuals into a training data set (n = 317,434) and a test data set (n = 74,825). In the training data set, we developed a common variant-based polygenic risk score for heel ultrasound speed of sound (SOS). Next, we performed burden testing to identify genes harboring rare determinants of bone mineral density, targeting influential rare variants with predicted high deleteriousness. We constructed a genetic risk score, called ggSOS, to incorporate influential rare variants in significant gene burden masks into the common variant-based polygenic risk score. We assessed the predictive performance of ggSOS in the White British test data set, as well as in populations of non-White British European (n = 18,885), African (n = 7165), East Asian (n = 2236), South Asian (n = 9829), and other admixed (n = 1481) ancestries. Twelve genes in pivotal regulatory pathways of bone homeostasis harbored influential rare variants associated with SOS (p < 5.5 × 10-7 ), including AHNAK, BMP5, CYP19A1, FAM20A, FBXW5, KDM5B, KREMEN1, LGR4, LRP5, SMAD6, SOST, and WNT1. Among 4013 (5.4%) individuals in the test data set carrying these variants, a one standard deviation decrease in ggSOS was associated with 1.35-fold (95% confidence interval [CI] 1.16-1.57) increased hazard of major osteoporotic fracture. However, compared with a common variant-based polygenic risk score (C-index = 0.641), ggSOS had only marginally improved prediction accuracy in identifying at-risk individuals (C-index = 0.644), with overlapping confidence intervals. Similarly, ggSOS did not demonstrate substantially improved predictive performance in non-European ancestry populations. In summary, modeling the effects of rare genetic determinants may assist polygenic prediction of fracture risk among carriers of influential rare variants. Nonetheless, improved clinical utility is not guaranteed for population-level risk screening. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Genetic divergences and hybridization within the Sebastes inermis complex.

The Sebastes inermis complex includes three sympatric species (Sebastes cheni, viz Sebastes inermis, and Sebastes ventricosus) with clear ecomorphological differences, albeit incomplete reproductive isolation. The presence of putative morphological hybrids (PMH) with plausibly higher fitness than the parent species indicates the need to confirm whether hybridization occurs within the complex. In this sense, we assessed the dynamics of genetic divergence and hybridization within the species complex using a panel of 10 microsatellite loci, and sequences of the mitochondrial control region (D-loop) and the intron-free rhodopsin (RH1) gene. The analyses revealed the presence of three distinct genetic clusters, large genetic distances using D-loop sequences, and distinctive mutations within the RH1 gene. These results are consistent with the descriptions of the three species. Two microsatellite loci had signatures of divergent selection, indicating that they are linked to genomic regions that are crucial for speciation. Furthermore, nonsynonymous mutations within the RH1 gene detected in S. cheni and "Kumano" (a PMH) suggest dissimilar adaptations related to visual perception in dim-light environments. The presence of individuals with admixed ancestry between two species confirmed hybridization. The presence of nonsynonymous mutations within the RH1 gene and the admixed ancestry of the "Kumano" morphotype highlight the potential role of hybridization in generating novelties within the species complex. We discuss possible outcomes of hybridization within the species complex, considering hybrid fitness and assortative mating. Overall, our findings indicate that the genetic divergence of each species is maintained in the presence of hybridization, as expected in a scenario of speciation-with-gene-flow.

Environment predicts the maintenance of reproductive isolation in a mosaic hybrid zone of rubber rabbitbrush.

Widely distributed plants of western North America experience divergent selection across environmental gradients, have complex histories shaped by biogeographic barriers and distributional shifts and often illustrate continuums of reproductive isolation. Rubber rabbitbrush (Ericameria nauseosa) is a foundational shrub species that occurs across diverse environments of western North America. Its remarkable phenotypic diversity is currently ascribed to two subspecies-Ericameria nauseosa nauseosa and Ericameria nauseosa consimilis-and 22 named varieties. To understand how genetic variation is partitioned across subspecies, varieties, and environments, we used high throughput sequencing of reduced representation libraries. We found clear evidence for divergence between the two subspecies, despite largely sympatric distributions. Numerous locations exhibiting admixed ancestry were not geographically localized but were widely distributed across a mosaic hybrid zone. The occurrence of hybrid and subspecific ancestries was strongly predicted by environmental variables as well as the proximity to major ecotones between ecoregions. Although this repeatability illustrates the importance of environmental factors in shaping reproductive isolation, variability in the prevalence of hybridization also indicates these factors likely differ across ecological contexts. There was mixed evidence for the evolutionary cohesiveness of varieties, but several genetically distinct and narrow endemic varieties exhibited admixed subspecific ancestries, hinting at the possibility for transgressive hybridization to contribute to phenotypic novelty and the colonization of new environments in E. nauseosa.

Identification of genetic risk loci and causal insights associated with Parkinson's disease in African and African admixed populations: a genome-wide association study

An understanding of the genetic mechanisms underlying diseases in ancestrally diverse populations is an important step towards development of targeted treatments. Research in African and African admixed populations can enable mapping of complex traits, because of their genetic diversity, extensive population substructure, and distinct linkage disequilibrium patterns. We aimed to do a comprehensive genome-wide assessment in African and African admixed individuals to better understand the genetic architecture of Parkinson's disease in these underserved populations. We performed a genome-wide association study (GWAS) in people of African and African admixed ancestry with and without Parkinson's disease. Individuals were included from several cohorts that were available as a part of the Global Parkinson's Genetics Program, the International Parkinson's Disease Genomics Consortium Africa, and 23andMe. A diagnosis of Parkinson's disease was confirmed clinically by a movement disorder specialist for every individual in each cohort, except for 23andMe, in which it was self-reported based on clinical diagnosis. We characterised ancestry-specific risk, differential haplotype structure and admixture, coding and structural genetic variation, and enzymatic activity. We included 197 918 individuals (1488 cases and 196 430 controls) in our genome-wide analysis. We identified a novel common risk factor for Parkinson's disease (overall meta-analysis odds ratio for risk of Parkinson's disease 1·58 [95% CI 1·37-1·80], p=2·397 × 10-14) and age at onset at the GBA1 locus, rs3115534-G (age at onset β=-2·00 [SE=0·57], p=0·0005, for African ancestry; and β=-4·15 [0·58], p=0·015, for African admixed ancestry), which was rare in non-African or non-African admixed populations. Downstream short-read and long-read whole-genome sequencing analyses did not reveal any coding or structural variant underlying the GWAS signal. The identified signal seems to be associated with decreased glucocerebrosidase activity. Our study identified a novel genetic risk factor in GBA1 in people of African ancestry, which has not been seen in European populations, and it could be a major mechanistic basis of Parkinson's disease in African populations. This population-specific variant exerts substantial risk on Parkinson's disease as compared with common variation identified through GWAS and it was found to be present in 39% of the cases assessed in this study. This finding highlights the importance of understanding ancestry-specific genetic risk in complex diseases, a particularly crucial point as the Parkinson's disease field moves towards targeted treatments in clinical trials. The distinctive genetics of African populations highlights the need for equitable inclusion of ancestrally diverse groups in future trials, which will be a valuable step towards gaining insights into novel genetic determinants underlying the causes of Parkinson's disease. This finding opens new avenues towards RNA-based and other therapeutic strategies aimed at reducing lifetime risk of Parkinson's disease. The Global Parkinson's Genetics Program, which is funded by the Aligning Science Across Parkinson's initiative, and The Michael J Fox Foundation for Parkinson's Research.

Causal effects on complex traits are similar for common variants across segments of different continental ancestries within admixed individuals.

Individuals of admixed ancestries (for example, African Americans) inherit a mosaic of ancestry segments (local ancestry) originating from multiple continental ancestral populations. This offers the unique opportunity of investigating the similarity of genetic effects on traits across ancestries within the same population. Here we introduce an approach to estimate correlation of causal genetic effects (radmix) across local ancestries and analyze 38 complex traits in African-European admixed individuals (N = 53,001) to observe very high correlations (meta-analysis radmix = 0.95, 95% credible interval 0.93-0.97), much higher than correlation of causal effects across continental ancestries. We replicate our results using regression-based methods from marginal genome-wide association study summary statistics. We also report realistic scenarios where regression-based methods yield inflated heterogeneity-by-ancestry due to ancestry-specific tagging of causal effects, and/or polygenicity. Our results motivate genetic analyses that assume minimal heterogeneity in causal effects by ancestry, with implications for the inclusion of ancestry-diverse individuals in studies.

Hybridization in the absence of an ecotone favors hybrid success in woodrats (Neotoma spp.).

Hybridization is a common process that has broadly impacted the evolution of multicellular eukaryotes; however, how ecological factors influence this process remains poorly understood. Here, we report the findings of a 3-year recapture study of the Bryant's woodrat (Neotoma bryanti) and desert woodrat (Neotoma lepida), two species that hybridize within a creosote bush (Larrea tridentata) shrubland in Whitewater, CA, USA. We used a genotype-by-sequencing approach to characterize the ancestry distribution of individuals across this hybrid zone coupled with Cormack-Jolly-Seber modeling to describe demography. We identified a high frequency of hybridization at this site with ~40% of individuals possessing admixed ancestry, which is the result of multigenerational backcrossing and advanced hybrid-hybrid crossing. F1, F2, and advanced generation hybrids had apparent survival rates similar to parental N. bryanti, while parental and backcross N. lepida had lower apparent survival rates and were far less abundant. Compared to bimodal hybrid zones where hybrids are often rare and selected against, we find that hybrids at Whitewater are common and have comparable survival to the dominant parental species, N. bryanti. The frequency of hybridization at Whitewater is therefore likely limited by the abundance of the less common parental species, N. lepida, rather than selection against hybrids.

Abstract A006: Increased germline mutational burden in individuals of African ancestry: Implications for interpretation of tumor mutation burden

Abstract Introduction: Tumor mutation burden (TMB), defined as the number of somatic gene mutations per megabase in a tumor genome, is used clinically to identify cancer patients that may respond to immune checkpoint inhibitors. Recent studies suggest that patient ancestry can influence TMB, where individuals of African ancestry were found to have elevated TMB values based on tumor-only sequencing analysis. However, the impact of patient ancestry on germline mutational burden and implications for interpretation of tumor-only mutational burden data is largely unexplored. Methods: We examined the influence of patient ancestry on germline and tumor mutation burden using tumor-normal whole exome sequencing (WES) data from a pan-cancer cohort of 1228 individuals from a single institution. Genetic ancestry was estimated from constitutional WES data using single nucleotide polymorphism weights from external reference panels. Variant calling was performed using constitutional, tumor-only, and paired tumor-normal workflows. Total and loss-of-function burden scores were calculated for each participant from each workflow based on the total number of mutations detected in each sample and the total number of predicted loss-of-function mutations, based on snpEff annotations, respectively. Results: Genetic ancestry analysis found that one-third of this pan-cancer cohort was of non-European ancestry. 9.4% of individuals were of Eastern Asian ancestry, 5.3% of African ancestry, 1.0% of South Asian ancestry, 0.2% of Native American ancestry, and 17.1% of admixed ancestry, predominantly European and Native American admixed ancestry (15.6%). Total and loss-of-function germline burden scores varied across ancestral groups, with individuals of African ancestry showing significantly increased germline burden scores compared to other ancestral groups (p&amp;lt;0.0001). Tumor-only analysis also showed increased mutation burden for individuals of African ancestry. However, when private and rare germline variation was taken into account using paired tumor-normal analysis, tumor mutation burden did not differ based on patient ancestry, suggesting the differences seen in tumor-only analysis are due to germline variation rather than differences in true somatic mutation burden. Conclusion: Our study reveals that the paucity of knowledge of ancestry-specific reference genomes leads to unacceptable health disparities in cancer, specifically in patients with African ancestry. We show that germline mutational burden varies by ancestral background, with individuals of African ancestry displaying increased genetic variation compared to other ancestral populations. These results suggest that patient ancestry should be considered when interpreting tumor mutational burden values, particularly from tumor-only analysis. Approaches utilizing paired tumor-normal analysis or ancestry-specific reference genomes can aid in more accurate assessment of TMB, thereby avoiding futile treatments targeted at presumed high mutational burden tumors specifically in African American populations. Citation Format: Janith Don, Sara A. Byron, Guangfa Zhang, Tyler Izatt, Jiaming Zhang, Bethany Davis, Bryce Turner, Jonathan J. Keats, Jeffrey M. Trent, Lorna Rodriguez-Rodriguez, Nicholas J. Schork. Increased germline mutational burden in individuals of African ancestry: Implications for interpretation of tumor mutation burden [abstract]. In: Proceedings of the 15th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2022 Sep 16-19; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr A006.

Local Ancestry-Informed Candidate Pathway Analysis of Warfarin Stable Dose in Latino Populations.

Accuracy of warfarin dose prediction algorithms may be improved by including data from diverse populations in genetic studies of dose variability. Here, we surveyed single nucleotide polymorphisms in vitamin K-related genetic pathways for association with warfarin dose requirements in two admixed Latino populations in standard-principal component adjusted and contemporary-local ancestry adjusted regression models. A total of five variants from vitamin K-related genes/pathways were associated with warfarin dose in both cohorts (P < 0.0125) in standard models. Local ancestry-adjusted analysis unveiled 35 associated variants with absolute effects ranging from β=9.04 ( ±2.23) to 39.18 ( ±10.89) per ancestral allele in the discovery cohort and β=6.47 (± 2.02) to 17.82 (± 6.83) in the replication cohort. Importantly, we demonstrate the technical validity of the Tractor model in cohorts with admixed ancestry from three founder populations and bring attention to the technical hurdles obstructing the inclusion of diverse, especially admixed, populations in pharmacogenomic research.

Admixed ancestral composition with Amerindian predominance at the Peruvian Alzheimer Disease Initiative (PeADI)

AbstractBackgroundCurrent genetic studies for AD and other dementias are making efforts to incorporate underrepresented populations including admixed Latinos. Peruvian population is characterized by admixed ancestry with a significant Amerindian component, which varies according to specific regions across Peru. The Peruvian Alzheimer Disease initiative (PeADI) was developed to ascertain a cohort for AD and other related dementias for genetic studies in Peru. We aim to determine the patterns of continental ancestry by regions across Peru.MethodsOver the last 3 years, The PeADI study has recruited 212 unrelated cognitive participants through collaborative health centers and community outreach ascertainment strategies. Cases were assessed by neurologists following NINDS‐ADRDA criteria. Controls were screened using MMSE, Clock drawing test and Pfeffer functional activities questionnaire. Genome‐wide genotyping was performed by Illumina screening array. PC‐AiR and model‐based. The cohort was divided into five regions (Northern, Southern, Lima, Central Highlands, and Amazonian) based on place of birth of the participant or the ancient known ancestor.ResultsThe global Admixture analysis showed that Peruvians have a substantial Amerindian component (63.6%), followed by European (35.9%), African (2.5%) and East Asian (2.1%) components. When analyzed by regions, we found that the Central region concentrates the highest Amerindian ancestry (72.7%), followed by the Northern region (58.2%), Lima and Callao (57%), the Amazonian region (55.6%), and the Southern region (55.1%). The highest European ancestry is located in the Amazonian region (42.1%), while the Central region has the lowest (26.6%). African and East Asian ancestry has little influence in Peru, being the Northern with the higher African component (4.1%), and Lima and Callao the region with the higher East Asian component (3.5%). There is no significant for Amerindian component across five regions (p=0.054); however, Central highlands region has higher Amerindian component compared to Lima (p=0.005), the Northern (p=0.005) and Southern(p=0.017) regions.ConclusionOur results confirmed the ancestry admixture of the Peruvian population with predominance of the Amerindian component. The Central region concentrates the highest Amerindian ancestry compared with other regions across Peru.