Using the UK Biobank as a global reference of worldwide populations: application to measuring ancestry diversity from GWAS summary statistics

A novel method for cell deconvolution using DNA methylation in PCA space

BACKGROUND

In this study, we present a novel method for reference-based cell deconvolution using data from DNA methylation arrays. Different from existing methods like IDOL-Ext, which operate on probe-level data, our approach represents features in the principal component analysis (PCA) space for cell type deconvolution.

RESULTS

Our method's accuracy in estimating cell compositions is validated across various public datasets, including blood samples from glioma patients. It demonstrates precision comparable to IDOL-Ext, with R2 values ranging from 0.73 to 0.99 for most cell types, while offering improved discrimination between similar cell types, particularly T cell subtypes in glioma patient samples (R2 0.42-0.75 vs. 0.36-0.66 for IDOL-Ext). However, both methods showed lower accuracy for certain cell types, such as memory CD8 T cells in glioma patients (R2 0.42 vs. 0.36 for IDOL-Ext), highlighting the challenges in distinguishing closely related cell populations. We have made this method available as an R package "BloodCellDecon" on GitHub.

CONCLUSIONS

Our study confirms the efficacy of cell type deconvolution in PCA space. The results indicate wide-ranging applicability and potential for adaptation to other forms of genomic data.

scAI-SNP: a method for inferring ancestry from single-cell data

Collaborative efforts, such as the Human Cell Atlas, are rapidly accumulating large amounts of single-cell data. To ensure that single-cell atlases are representative of human genetic diversity, we need to determine the ancestry of the donors from whom single-cell data are generated. Self-reporting of race and ethnicity, although important, can be biased and is not always available for the datasets already collected. Here, we introduce scAI-SNP, a tool to infer ancestry directly from single-cell genomics data. To train scAI-SNP, we identified 4.5 million ancestry-informative single-nucleotide polymorphisms (SNPs) in the 1000 Genomes Project dataset across 3201 individuals from 26 population groups. For a query single-cell data set, scAI-SNP uses these ancestry-informative SNPs to compute the contribution of each of the 26 population groups to the ancestry of the donor from whom the cells were obtained. Using diverse single-cell data sets with matched whole-genome sequencing data, we show that scAI-SNP is robust to the sparsity of single-cell data, can accurately and consistently infer ancestry from samples derived from diverse types of tissues and cancer cells, and can be applied to different modalities of single-cell profiling assays, such as single-cell RNA-seq and single-cell ATAC-seq. Finally, we argue that ensuring that single-cell atlases represent diverse ancestry, ideally alongside race and ethnicity, is ultimately important for improved and equitable health outcomes by accounting for human diversity.

Characterizing substructure via mixture modeling in large-scale genetic summary statistics

Genetic summary data are broadly accessible and highly useful including for risk prediction, causal inference, fine mapping, and incorporation of external controls. However, collapsing individual-level data into groups masks intra- and inter-sample heterogeneity, leading to confounding, reduced power, and bias. Ultimately, unaccounted substructure limits summary data usability, especially for understudied or admixed populations. Here, we present Summix2, a comprehensive set of methods and software based on a computationally efficient mixture model to estimate and adjust for substructure in genetic summary data. In extensive simulations and application to public data, Summix2 characterizes finer-scale population structure, identifies ascertainment bias, and identifies potential regions of selection due to local substructure deviation. Summix2 increases the robust use of diverse publicly available summary data resulting in improved and more equitable research.

On the effect heterogeneity of established disease susceptibility loci for Alzheimer's disease across different genetic ancestries

INTRODUCTION

Genome-wide association studies have identified numerous disease susceptibility loci (DSLs) for Alzheimer's disease (AD). However, only a limited number of studies have investigated the dependence of the genetic effect size of established DSLs on genetic ancestry.

METHODS

We utilized the whole genome sequencing data from the Alzheimer's Disease Sequencing Project (ADSP) including 35,569 participants. A total of 25,459 subjects in four distinct populations (African ancestry, non-Hispanic White, admixed Hispanic, and Asian) were analyzed.

RESULTS

We found that nine DSLs showed significant heterogeneity across populations. Single nucleotide polymorphism (SNP) rs2075650 in translocase of outer mitochondrial membrane 40 (TOMM40) showed the largest heterogeneity (Cochran's Q = 0.00, I2 = 90.08), followed by other SNPs in apolipoprotein C1 (APOC1) and apolipoprotein E (APOE). Two additional loci, signal-induced proliferation-associated 1 like 2 (SIPA1L2) and solute carrier 24 member 4 (SLC24A4), showed significant heterogeneity across populations.

DISCUSSION

We observed substantial heterogeneity for the APOE-harboring 19q13.32 region with TOMM40/APOE/APOC1 genes. The largest risk effect was seen among African Americans, while Asians showed a surprisingly small risk effect.

Adjusting for principal components can induce spurious associations in genome-wide association studies in admixed populations

Principal component analysis (PCA) is widely used to control for population structure in genome-wide association studies (GWAS). Top principal components (PCs) typically reflect population structure, but challenges arise in deciding how many PCs are needed and ensuring that PCs do not capture other artifacts such as regions with atypical linkage disequilibrium (LD). In response to the latter, many groups suggest performing LD pruning or excluding known high LD regions prior to PCA. However, these suggestions are not universally implemented and the implications for GWAS are not fully understood, especially in the context of admixed populations. In this paper, we investigate the impact of pre-processing and the number of PCs included in GWAS models in African American samples from the Women's Women's Health Initiative SNP Health Association Resource and two Trans-Omics for Precision Medicine Whole Genome Sequencing Project contributing studies (Jackson Heart Study and Genetic Epidemiology of Chronic Obstructive Pulmonary Disease Study). In all three samples, we find the first PC is highly correlated with genome-wide ancestry whereas later PCs often capture local genomic features. The pattern of which, and how many, genetic variants are highly correlated with individual PCs differs from what has been observed in prior studies focused on European populations and leads to distinct downstream consequences: adjusting for such PCs yields biased effect size estimates and elevated rates of spurious associations due to the phenomenon of collider bias. Excluding high LD regions identified in previous studies does not resolve these issues. LD pruning proves more effective, but the optimal choice of thresholds varies across datasets. Altogether, our work highlights unique issues that arise when using PCA to control for ancestral heterogeneity in admixed populations and demonstrates the importance of careful pre-processing and diagnostics to ensure that PCs capturing multiple local genomic features are not included in GWAS models.

Hunter-gatherer genetics research: Importance and avenues

Major developments in the field of genetics in the past few decades have revolutionised notions of what it means to be human. Although currently only a few populations around the world practise a hunting and gathering lifestyle, this mode of subsistence has characterised members of our species since its very origins and allowed us to migrate across the planet. Therefore, the geographical distribution of hunter-gatherer populations, dependence on local ecosystems and connections to past populations and neighbouring groups have provided unique insights into our evolutionary origins. However, given the vulnerable status of hunter-gatherers worldwide, the development of the field of anthropological genetics requires that we reevaluate how we conduct research with these communities. Here, we review how the inclusion of hunter-gatherer populations in genetics studies has advanced our understanding of human origins, ancient population migrations and interactions as well as phenotypic adaptations and adaptability to different environments, and the important scientific and medical applications of these advancements. At the same time, we highlight the necessity to address yet unresolved questions and identify areas in which the field may benefit from improvements.

Inferring disease architecture and predictive ability with LDpred2-auto

LDpred2 is a widely used Bayesian method for building polygenic scores (PGSs). LDpred2-auto can infer the two parameters from the LDpred model, the SNP heritability h2 and polygenicity p, so that it does not require an additional validation dataset to choose best-performing parameters. The main aim of this paper is to properly validate the use of LDpred2-auto for inferring multiple genetic parameters. Here, we present a new version of LDpred2-auto that adds an optional third parameter α to its model, for modeling negative selection. We then validate the inference of these three parameters (or two, when using the previous model). We also show that LDpred2-auto provides per-variant probabilities of being causal that are well calibrated and can therefore be used for fine-mapping purposes. We also introduce a formula to infer the out-of-sample predictive performance r2 of the resulting PGS directly from the Gibbs sampler of LDpred2-auto. Finally, we extend the set of HapMap3 variants recommended to use with LDpred2 with 37% more variants to improve the coverage of this set, and we show that this new set of variants captures 12% more heritability and provides 6% more predictive performance, on average, in UK Biobank analyses.

Long-term adverse influence of smoking during pregnancy on height and body size of offspring at ten years old in the UK Biobank cohort

Background

To explore the long-term relationship between maternal smoking during pregnancy and early childhood growth in the UK Biobank cohort.

Methods

To estimate the effect of maternal smoking during pregnancy on offspring height and body size at ten years old, we performed binary logistic analyses and reported odds ratios (OR) as well as 95% confidence intervals (95%CIs). We also implemented the cross-contextual comparison study to examine whether such influence could be repeatedly observed among three different ethnicities in the UK Biobank cohort (n = 22,140 for White, n = 7094 for South Asian, and n = 5000 for Black). In particular, we conducted the sibling cohort study in White sibling cohort (n = 9953 for height and n = 7239 for body size) to control for unmeasured familial confounders.

Results

We discovered that children whose mothers smoked during pregnancy had greater risk of being shorter or plumper at age ten in the full UK Biobank White cohort, with 15.3% (95% CIs: 13.0%∼17.7%) higher risk for height and 32.4% (95%CIs: 29.5%∼35.4%) larger risk for body size. Similar associations were identified in the South Asian and Black ethnicities. These associations were robust and remained significant in the White sibling cohort (12.6% [95%CIs: 5.0%∼20.3%] for height and 36.1% [95%CIs: 26.3%∼45.9%] for body size) after controlling for family factors.

Conclusion

This study robustly confirms that maternal smoking during pregnancy can promote height deficit and obesity for offspring at ten years old. Our findings strongly encourage mothers to quit smoking during pregnancy for improving growth and development of offspring.

AI pitfalls and what not to do: mitigating bias in AI

Various forms of artificial intelligence (AI) applications are being deployed and used in many healthcare systems. As the use of these applications increases, we are learning the failures of these models and how they can perpetuate bias. With these new lessons, we need to prioritize bias evaluation and mitigation for radiology applications; all the while not ignoring the impact of changes in the larger enterprise AI deployment which may have downstream impact on performance of AI models. In this paper, we provide an updated review of known pitfalls causing AI bias and discuss strategies for mitigating these biases within the context of AI deployment in the larger healthcare enterprise. We describe these pitfalls by framing them in the larger AI lifecycle from problem definition, data set selection and curation, model training and deployment emphasizing that bias exists across a spectrum and is a sequela of a combination of both human and machine factors.

Neural ADMIXTURE for rapid genomic clustering

Characterizing the genetic structure of large cohorts has become increasingly important as genetic studies extend to massive, increasingly diverse biobanks. Popular methods decompose individual genomes into fractional cluster assignments with each cluster representing a vector of DNA variant frequencies. However, with rapidly increasing biobank sizes, these methods have become computationally intractable. Here we present Neural ADMIXTURE, a neural network autoencoder that follows the same modeling assumptions as the current standard algorithm, ADMIXTURE, while reducing the compute time by orders of magnitude surpassing even the fastest alternatives. One month of continuous compute using ADMIXTURE can be reduced to just hours with Neural ADMIXTURE. A multi-head approach allows Neural ADMIXTURE to offer even further acceleration by calculating multiple cluster numbers in a single run. Furthermore, the models can be stored, allowing cluster assignment to be performed on new data in linear time without needing to share the training samples.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

ADAMS project: a genetic Association study in individuals from Diverse Ancestral backgrounds with Multiple Sclerosis based in the UK

PURPOSE

Genetic studies of multiple sclerosis (MS) susceptibility and severity have focused on populations of European ancestry. Studying MS genetics in other ancestral groups is necessary to determine the generalisability of these findings. The genetic Association study in individuals from Diverse Ancestral backgrounds with Multiple Sclerosis (ADAMS) project aims to gather genetic and phenotypic data on a large cohort of ancestrally-diverse individuals with MS living in the UK.

PARTICIPANTS

Adults with self-reported MS from diverse ancestral backgrounds. Recruitment is via clinical sites, online (https://app.mantal.co.uk/adams) or the UK MS Register. We are collecting demographic and phenotypic data using a baseline questionnaire and subsequent healthcare record linkage. We are collecting DNA from participants using saliva kits (Oragene-600) and genotyping using the Illumina Global Screening Array V.3.

FINDINGS TO DATE

As of 3 January 2023, we have recruited 682 participants (n=446 online, n=55 via sites, n=181 via the UK MS Register). Of this initial cohort, 71.2% of participants are female, with a median age of 44.9 years at recruitment. Over 60% of the cohort are non-white British, with 23.5% identifying as Asian or Asian British, 16.2% as Black, African, Caribbean or Black British and 20.9% identifying as having mixed or other backgrounds. The median age at first symptom is 28 years, and median age at diagnosis is 32 years. 76.8% have relapsing-remitting MS, and 13.5% have secondary progressive MS.

FUTURE PLANS

Recruitment will continue over the next 10 years. Genotyping and genetic data quality control are ongoing. Within the next 3 years, we aim to perform initial genetic analyses of susceptibility and severity with a view to replicating the findings from European-ancestry studies. In the long term, genetic data will be combined with other datasets to further cross-ancestry genetic discoveries.

Genome-wide association and multi-trait analyses characterize the common genetic architecture of heart failure

Heart failure is a leading cause of cardiovascular morbidity and mortality. However, the contribution of common genetic variation to heart failure risk has not been fully elucidated, particularly in comparison to other common cardiometabolic traits. We report a multi-ancestry genome-wide association study meta-analysis of all-cause heart failure including up to 115,150 cases and 1,550,331 controls of diverse genetic ancestry, identifying 47 risk loci. We also perform multivariate genome-wide association studies that integrate heart failure with related cardiac magnetic resonance imaging endophenotypes, identifying 61 risk loci. Gene-prioritization analyses including colocalization and transcriptome-wide association studies identify known and previously unreported candidate cardiomyopathy genes and cellular processes, which we validate in gene-expression profiling of failing and healthy human hearts. Colocalization, gene expression profiling, and Mendelian randomization provide convergent evidence for the roles of BCKDHA and circulating branch-chain amino acids in heart failure and cardiac structure. Finally, proteome-wide Mendelian randomization identifies 9 circulating proteins associated with heart failure or quantitative imaging traits. These analyses highlight similarities and differences among heart failure and associated cardiovascular imaging endophenotypes, implicate common genetic variation in the pathogenesis of heart failure, and identify circulating proteins that may represent cardiomyopathy treatment targets.

Identifying and correcting for misspecifications in GWAS summary statistics and polygenic scores

Publicly available genome-wide association studies (GWAS) summary statistics exhibit uneven quality, which can impact the validity of follow-up analyses. First, we present an overview of possible misspecifications that come with GWAS summary statistics. Then, in both simulations and real-data analyses, we show that additional information such as imputation INFO scores, allele frequencies, and per-variant sample sizes in GWAS summary statistics can be used to detect possible issues and correct for misspecifications in the GWAS summary statistics. One important motivation for us is to improve the predictive performance of polygenic scores built from these summary statistics. Unfortunately, owing to the lack of reporting standards for GWAS summary statistics, this additional information is not systematically reported. We also show that using well-matched linkage disequilibrium (LD) references can improve model fit and translate into more accurate prediction. Finally, we discuss how to make polygenic score methods such as lassosum and LDpred2 more robust to these misspecifications to improve their predictive power.