Comparing the effect of imputation reference panel composition in four distinct Latin American cohorts

Genome-wide association studies have been useful in identifying genetic risk factors for various phenotypes. These studies rely on imputation and many existing panels are largely composed of individuals of European ancestry, resulting in lower levels of imputation quality in underrepresented populations. We aim to analyze how the composition of imputation reference panels affects imputation quality in four target Latin American cohorts. We compared imputation quality for chromosomes 7 and X when altering the imputation reference panel by: 1) increasing the number of Latin American individuals; 2) excluding either Latin American, African, or European individuals, or 3) increasing the Indigenous American (IA) admixture proportions of included Latin Americans. We found that increasing the number of Latin Americans in the reference panel improved imputation quality in the four populations; however, there were differences between chromosomes 7 and X in some cohorts. Excluding Latin Americans from analysis resulted in worse imputation quality in every cohort, while differential effects were seen when excluding Europeans and Africans between and within cohorts and between chromosomes 7 and X. Finally, increasing IA-like admixture proportions in the reference panel increased imputation quality at different levels in different populations. The difference in results between populations and chromosomes suggests that existing and future reference panels containing Latin American individuals are likely to perform differently in different Latin American populations.

X-Chromosome Association Study in Latin American Cohorts Identifies New Loci in Parkinson's Disease

BACKGROUND

Sex differences in Parkinson's disease (PD) risk are well-known. However, the role of sex chromosomes in the development and progression of PD is still unclear.

OBJECTIVE

The objective of this study was to perform the first X-chromosome-wide association study for PD risk in a Latin American cohort.

METHODS

We used data from three admixed cohorts: (1) Latin American Research consortium on the Genetics of Parkinson's Disease (n = 1504) as discover cohort, and (2) Latino cohort from International Parkinson Disease Genomics Consortium (n = 155) and (3) Bambui Aging cohort (n = 1442) as replication cohorts. We also developed an X-chromosome framework specifically designed for admixed populations.

RESULTS

We identified eight linkage disequilibrium regions associated with PD. We replicated one of these regions (top variant rs525496; discovery odds ratio [95% confidence interval]: 0.60 [0.478-0.77], P = 3.13 × 10-5 replication odds ratio: 0.60 [0.37-0.98], P = 0.04). rs5525496 is associated with multiple expression quantitative trait loci in brain and non-brain tissues, including RAB9B, H2BFM, TSMB15B, and GLRA4, but colocalization analysis suggests that rs5525496 may not mediate risk by expression of these genes. We also replicated a previous X-chromosome-wide association study finding (rs28602900), showing that this variant is associated with PD in non-European populations.

CONCLUSIONS

Our results reinforce the importance of including X-chromosome and diverse populations in genetic studies. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Factors Associated With Survival Disparities Between Non-Hispanic Black and White Patients With Uterine Cancer

Importance

Disparities in survival exist between non-Hispanic Black (hereafter, Black) and non-Hispanic White (hereafter, White) patients with uterine cancer.

Objective

To investigate factors associated with racial disparities in survival between Black and White patients with uterine cancer.

Design, Setting, and Patients

This cohort study used data from the National Cancer Database on 274 838 Black and White patients who received a diagnosis of uterine cancer from January 1, 2004, to December 31, 2017, with follow-up through December 2020. Statistical analysis was performed in July 2022.

Main Outcomes and Measures

Overall survival by self-reported race and evaluation of explanatory study factors associated with hazard ratio (HR) reduction for Black vs White patients. A propensity scoring approach was applied sequentially to balance racial differences in demographic characteristics, comorbidity score, neighborhood income, insurance status, histologic subtype, disease stage, and treatment.

Results

The study included 32 230 Black female patients (mean [SD] age at diagnosis, 63.8 [10.0] years) and 242 608 White female patients (mean [SD] age at diagnosis, 63.5 [10.5] years) and had a median follow-up of 74.0 months (range, 43.5-113.8 months). Black patients were more likely than White patients to have low income (44.1% vs 14.0%), be uninsured (5.7% vs 2.6%), present with nonendometrioid histologic characteristics (46.1% vs 21.6%), have an advanced disease stage (34.1% vs 19.8%), receive first-line chemotherapy (33.8% vs 18.2%), and have worse 5-year survival (58.6% vs 78.5%). Among patients who received a diagnosis at younger than 65 years of age, the HR for death for Black vs White patients was 2.43 (95% CI, 2.34-2.52) in a baseline demographic-adjusted model and 1.29 (95% CI, 1.23-1.35) after balancing other factors. Comorbidity score, neighborhood income, insurance status, histologic subtype, disease stage, treatment, and unexplained factors accounted for 0.8%, 7.2%, 11.5%, 53.1%, 5.8%, 1.2%, and 20.4%, respectively, of the excess relative risk (ERR) among the younger Black vs White patients. Among patients 65 years or older, the HR for death for Black vs White patients was 1.87 (95% CI, 1.81-1.93) in the baseline model and 1.14 (95% CI, 1.09-1.19) after balancing other factors. Comorbidity score, neighborhood income, insurance status, histologic subtype, disease stage, treatment, and unexplained factors accounted for 3.0%, 7.5%, 0.0%, 56.2%, 10.6%, 6.9%, and 15.8%, respectively, of the ERR among Black vs White patients aged 65 years or older.

Conclusions and Relevance

This study suggests that histologic subtype was the dominant factor associated with racial survival disparity among patients with uterine cancer, while insurance status represented the main modifiable factor for women younger than 65 years. Additional studies of interactions between biology and social determinants of health are merited.

X-Chromosome Association Study in Latin American Cohorts Identifies New Loci in Parkinson Disease

Sex differences in Parkinson Disease (PD) risk are well-known. However, it is still unclear the role of sex chromosomes in the development and progression of PD. We performed the first X-chromosome Wide Association Study (XWAS) for PD risk in Latin American individuals. We used data from three admixed cohorts: (i) Latin American Research consortium on the GEnetics of Parkinson's Disease (n=1,504) as discover cohort and (ii) Latino cohort from International Parkinson Disease Genomics Consortium (n = 155) and (iii) Bambui Aging cohort (n= 1,442) as replication cohorts. After developing a X-chromosome framework specifically designed for admixed populations, we identified eight linkage disequilibrium regions associated with PD. We fully replicated one of these regions (top variant rs525496; discovery OR [95%CI]: 0.60 [0.478 - 0.77], p = 3.13 × 10 ^-5 ; replication OR: 0.60 [0.37-0.98], p = 0.04). rs525496 is an expression quantitative trait loci for several genes expressed in brain tissues, including RAB9B, H2BFM, TSMB15B and GLRA4 . We also replicated a previous XWAS finding (rs28602900), showing that this variant is associated with PD in non-European populations. Our results reinforce the importance of including X-chromosome and diverse populations in genetic studies.

Exome Array Analysis of 9721 Ischemic Stroke Cases from the SiGN Consortium

Recent genome wide association studies have identified 89 common genetic variants robustly associated with ischemic stroke and primarily located in non-coding regions. To evaluate the contribution of coding variants, which are mostly rare, we performed an exome array analysis on 106,101 SNPs for 9721 ischemic stroke cases from the SiGN Consortium, and 12,345 subjects with no history of stroke from the Health Retirement Study and SiGN consortium. We identified 15 coding variants significantly associated with all ischemic stroke at array-wide threshold (i.e., p < 4.7 × 10-7), including two common SNPs in ABO that have previously been associated with stroke. Twelve of the remaining 13 variants were extremely rare in European Caucasians (MAF < 0.1%) and the associations were driven by African American samples. There was no evidence for replication of these associations in either TOPMed Stroke samples (n = 5613 cases) or UK Biobank (n = 5874 stroke cases), although power to replicate was very low given the low allele frequencies of the associated variants and a shortage of samples from diverse ancestries. Our study highlights the need for acquiring large, well-powered diverse cohorts to study rare variants, and the technical challenges using array-based genotyping technologies for rare variant genotyping.

Tracing the Distribution of European Lactase Persistence Genotypes Along the Americas

In adulthood, the ability to digest lactose, the main sugar present in milk of mammals, is a phenotype (lactase persistence) observed in historically herder populations, mainly Northern Europeans, Eastern Africans, and Middle Eastern nomads. As the -13910∗T allele in the MCM6 gene is the most well-characterized allele responsible for the lactase persistence phenotype, the -13910C > T (rs4988235) polymorphism is commonly evaluated in lactase persistence studies. Lactase non-persistent adults may develop symptoms of lactose intolerance when consuming dairy products. In the Americas, there is no evidence of the consumption of these products until the arrival of Europeans. However, several American countries' dietary guidelines recommend consuming dairy for adequate human nutrition and health promotion. Considering the extensive use of dairy and the complex ancestry of Pan-American admixed populations, we studied the distribution of -13910C > T lactase persistence genotypes and its flanking haplotypes of European origin in 7,428 individuals from several Pan-American admixed populations. We found that the -13910∗T allele frequency in Pan-American admixed populations is directly correlated with allele frequency of the European sources. Moreover, we did not observe any overrepresentation of European haplotypes in the -13910C > T flanking region, suggesting no selective pressure after admixture in the Americas. Finally, considering the dominant effect of the -13910∗T allele, our results indicate that Pan-American admixed populations are likely to have higher frequency of lactose intolerance, suggesting that general dietary guidelines deserve further evaluation across the continent.

Genome-Wide Analysis of Copy Number Variation in Latin American Parkinson's Disease Patients

BACKGROUND

Parkinson's disease is the second most common neurodegenerative disorder and affects people from all ethnic backgrounds, yet little is known about the genetics of Parkinson's disease in non-European populations. In addition, the overall identification of copy number variants at a genome-wide level has been understudied in Parkinson's patients. The objective of this study was to understand the genome-wide burden of copy number variants in Latinos and its association with Parkinson's disease.

METHODS

We used genome-wide genotyping data from 747 Parkinson's disease patients and 632 controls from the Latin American Research Consortium on the Genetics of Parkinson's disease.

RESULTS

Genome-wide copy number burden analysis showed that patients were significantly enriched for copy number variants overlapping known Parkinson's disease genes compared with controls (odds ratio, 3.97; 95%CI, 1.69-10.5; P = 0.018). PRKN showed the strongest copy number burden, with 20 copy number variant carriers. These patients presented an earlier age of disease onset compared with patients with other copy number variants (median age at onset, 31 vs 57 years, respectively; P = 7.46 × 10-7 ).

CONCLUSIONS

We found that although overall genome-wide copy number variant burden was not significantly different, Parkinson's disease patients were significantly enriched with copy number variants affecting known Parkinson's disease genes. We also identified that of 250 patients with early-onset disease, 5.6% carried a copy number variant on PRKN in our cohort. Our study is the first to analyze genome-wide copy number variant association in Latino Parkinson's disease patients and provides insights about this complex disease in this understudied population. © 2020 International Parkinson and Movement Disorder Society.

Detecting geospatial patterns of Plasmodium falciparum parasite migration in Cambodia using optimized estimated effective migration surfaces

BACKGROUND

Understanding the genetic structure of natural populations provides insight into the demographic and adaptive processes that have affected those populations. Such information, particularly when integrated with geospatial data, can have translational applications for a variety of fields, including public health. Estimated effective migration surfaces (EEMS) is an approach that allows visualization of the spatial patterns in genomic data to understand population structure and migration. In this study, we developed a workflow to optimize the resolution of spatial grids used to generate EEMS migration maps and applied this optimized workflow to estimate migration of Plasmodium falciparum in Cambodia and bordering regions of Thailand and Vietnam.

METHODS

The optimal density of EEMS grids was determined based on a new workflow created using density clustering to define genomic clusters and the spatial distance between genomic clusters. Topological skeletons were used to capture the spatial distribution for each genomic cluster and to determine the EEMS grid density; i.e., both genomic and spatial clustering were used to guide the optimization of EEMS grids. Model accuracy for migration estimates using the optimized workflow was tested and compared to grid resolutions selected without the optimized workflow. As a test case, the optimized workflow was applied to genomic data generated from P. falciparum sampled in Cambodia and bordering regions, and migration maps were compared to estimates of malaria endemicity, as well as geographic properties of the study area, as a means of validating observed migration patterns.

RESULTS

Optimized grids displayed both high model accuracy and reduced computing time compared to grid densities selected in an unguided manner. In addition, EEMS migration maps generated for P. falciparum using the optimized grid corresponded to estimates of malaria endemicity and geographic properties of the study region that might be expected to impact malaria parasite migration, supporting the validity of the observed migration patterns.

CONCLUSIONS

Optimized grids reduce spatial uncertainty in the EEMS contours that can result from user-defined parameters, such as the resolution of the spatial grid used in the model. This workflow will be useful to a broad range of EEMS users as it can be applied to analyses involving other organisms of interest and geographic areas.

Ancestral characterization of 1018 cancer cell lines highlights disparities and reveals gene expression and mutational differences

Background

Although cell lines are an essential resource for studying cancer biology, many are of unknown ancestral origin, and their use may not be optimal for evaluating the biology of all patient populations.

Methods

An admixture analysis was performed using genome‐wide chip data from the Catalogue of Somatic Mutations in Cancer (COSMIC) Cell Lines Project to calculate genetic ancestry estimates for 1018 cancer cell lines. After stratifying the analyses by tissue and histology types, linear models were used to evaluate the influence of ancestry on gene expression and somatic mutation frequency.

Results

For the 701 cell lines with unreported ancestry, 215 were of East Asian origin, 30 were of African or African American origin, and 453 were of European origin. Notable imbalances were observed in ancestral representation across tissue type, with the majority of analyzed tissue types having few cell lines of African American ancestral origin, and with Hispanic and South Asian ancestry being almost entirely absent across all cell lines. In evaluating gene expression across these cell lines, expression levels of the genes neurobeachin line 1 (NBEAL1), solute carrier family 6 member 19 (SLC6A19), HEAT repeat containing 6 (HEATR6), and epithelial cell transforming 2 like (ECT2L) were associated with ancestry. Significant differences were also observed in the proportions of somatic mutation types across cell lines with varying ancestral proportions.

Conclusions

By estimating genetic ancestry for 1018 cancer cell lines, the authors have produced a resource that cancer researchers can use to ensure that their cell lines are ancestrally representative of the populations they intend to affect. Furthermore, the novel ancestry‐specific signal identified underscores the importance of ancestral awareness when studying cancer.

Preclinical cancer cell line research is often conducted without an awareness of ancestral background, which results in incongruities between the genetic backgrounds of used cell lines and the patient populations they are intended to represent. By calculating genetic ancestry for 1018 common cancer cell lines and identifying ancestry‐specific expression and somatic mutation patterns, the importance of ancestral awareness is emphasized, and a resource is provided that can be used by cancer researchers to ensure that their cell lines are ancestrally representative of the populations they aim to impact.

Identifying tagging SNPs for African specific genetic variation from the African Diaspora Genome

A primary goal of The Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to develop an 'African Diaspora Power Chip' (ADPC), a genotyping array consisting of tagging SNPs, useful in comprehensively identifying African specific genetic variation. This array is designed based on the novel variation identified in 642 CAAPA samples of African ancestry with high coverage whole genome sequence data (~30× depth). This novel variation extends the pattern of variation catalogued in the 1000 Genomes and Exome Sequencing Projects to a spectrum of populations representing the wide range of West African genomic diversity. These individuals from CAAPA also comprise a large swath of the African Diaspora population and incorporate historical genetic diversity covering nearly the entire Atlantic coast of the Americas. Here we show the results of designing and producing such a microchip array. This novel array covers African specific variation far better than other commercially available arrays, and will enable better GWAS analyses for researchers with individuals of African descent in their study populations. A recent study cataloging variation in continental African populations suggests this type of African-specific genotyping array is both necessary and valuable for facilitating large-scale GWAS in populations of African ancestry.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

Fine-scale patterns of population stratification confound rare variant association tests

Advances in next-generation sequencing technology have enabled systematic exploration of the contribution of rare variation to Mendelian and complex diseases. Although it is well known that population stratification can generate spurious associations with common alleles, its impact on rare variant association methods remains poorly understood. Here, we performed exhaustive coalescent simulations with demographic parameters calibrated from exome sequence data to evaluate the performance of nine rare variant association methods in the presence of fine-scale population structure. We find that all methods have an inflated spurious association rate for parameter values that are consistent with levels of differentiation typical of European populations. For example, at a nominal significance level of 5%, some test statistics have a spurious association rate as high as 40%. Finally, we empirically assess the impact of population stratification in a large data set of 4,298 European American exomes. Our results have important implications for the design, analysis, and interpretation of rare variant genome-wide association studies.