Extent to which array genotyping and imputation with large reference panels approximate deep whole-genome sequencing

The predictive capacity of polygenic risk scores for disease risk is only moderately influenced by imputation panels tailored to the target population

MOTIVATION

Polygenic risk scores (PRSs) predict individuals' genetic risk of developing complex diseases. They summarize the effect of many variants discovered in genome-wide association studies (GWASs). However, to date, large GWASs exist primarily for the European population and the quality of PRS prediction declines when applied to other ethnicities. Genetic profiling of individuals in the discovery set (on which the GWAS was performed) and target set (on which the PRS is applied) is typically done by SNP arrays that genotype a fraction of common SNPs. Therefore, a key step in GWAS analysis and PRS calculation is imputing untyped SNPs using a panel of fully sequenced individuals. The imputation results depend on the ethnic composition of the imputation panel. Imputing genotypes with a panel of individuals of the same ethnicity as the genotyped individuals typically improves imputation accuracy. However, there has been no systematic investigation into the influence of the ethnic composition of imputation panels on the accuracy of PRS predictions when applied to ethnic groups that differ from the population used in the GWAS.

RESULTS

We estimated the effect of imputation of the target set on prediction accuracy of PRS when the discovery and the target sets come from different ethnic groups. We analyzed binary phenotypes on ethnically distinct sets from the UK Biobank and other resources. We generated ethnically homogenous panels, imputed the target sets, and generated PRSs. Then, we assessed the prediction accuracy obtained from each imputation panel. Our analysis indicates that using an imputation panel matched to the ethnicity of the target population yields only a marginal improvement and only under specific conditions.

AVAILABILITY AND IMPLEMENTATION

The source code used for executing the analyses is this paper is available at https://github.com/Shamir-Lab/PRS-imputation-panels.

The hazards of genotype imputation when mapping disease susceptibility variants

BACKGROUND

The cost-free increase in statistical power of using imputation to infer missing genotypes is undoubtedly appealing, but is it hazard-free? This case study of three type-2 diabetes (T2D) loci demonstrates that it is not; it sheds light on why this is so and raises concerns as to the shortcomings of imputation at disease loci, where haplotypes differ between cases and reference panel.

RESULTS

T2D-associated variants were previously identified using targeted sequencing. We removed these significantly associated SNPs and used neighbouring SNPs to infer them by imputation. We compared imputed with observed genotypes, examined the altered pattern of T2D-SNP association, and investigated the cause of imputation errors by studying haplotype structure. Most T2D variants were incorrectly imputed with a low density of scaffold SNPs, but the majority failed to impute even at high density, despite obtaining high certainty scores. Missing and discordant imputation errors, which were observed disproportionately for the risk alleles, produced monomorphic genotype calls or false-negative associations. We show that haplotypes carrying risk alleles are considerably more common in the T2D cases than the reference panel, for all loci.

CONCLUSIONS

Imputation is not a panacea for fine mapping, nor for meta-analysing multiple GWAS based on different arrays and different populations. A total of 80% of the SNPs we have tested are not included in array platforms, explaining why these and other such associated variants may previously have been missed. Regardless of the choice of software and reference haplotypes, imputation drives genotype inference towards the reference panel, introducing errors at disease loci.

Genotype imputation accuracy and the quality metrics of the minor ancestry in multi-ancestry reference panels

Large-scale imputation reference panels are currently available and have contributed to efficient genome-wide association studies through genotype imputation. However, whether large-size multi-ancestry or small-size population-specific reference panels are the optimal choices for under-represented populations continues to be debated. We imputed genotypes of East Asian (180k Japanese) subjects using the Trans-Omics for Precision Medicine reference panel and found that the standard imputation quality metric (Rsq) overestimated dosage r2 (squared correlation between imputed dosage and true genotype) particularly in marginal-quality bins. Variance component analysis of Rsq revealed that the increased imputed-genotype certainty (dosages closer to 0, 1 or 2) caused upward bias, indicating some systemic bias in the imputation. Through systematic simulations using different template switching rates (θ value) in the hidden Markov model, we revealed that the lower θ value increased the imputed-genotype certainty and Rsq; however, dosage r2 was insensitive to the θ value, thereby causing a deviation. In simulated reference panels with different sizes and ancestral diversities, the θ value estimates from Minimac decreased with the size of a single ancestry and increased with the ancestral diversity. Thus, Rsq could be deviated from dosage r2 for a subpopulation in the multi-ancestry panel, and the deviation represents different imputed-dosage distributions. Finally, despite the impact of the θ value, distant ancestries in the reference panel contributed only a few additional variants passing a predefined Rsq threshold. We conclude that the θ value substantially impacts the imputed dosage and the imputation quality metric value.

Populational pan-ethnic screening panel enabled by deep whole genome sequencing

Birth defect is a global threat to the public health systems. Mitigating neonatal anomalies is hampered by elusive molecular mechanisms of pathogenic mutations and poor subsequent translation into preventative measures. Applying appropriate strategies in China to promote reproductive health is particularly challenging, as the Chinese population compromises complex genomic diversity due to the inclusion of many ethnic groups with distinct genetic backgrounds. To investigate and evaluate the feasibility of implementing a pan-ethnic screening strategy, and guide future reproductive counselling, high-quality variants associated with autosome recessive (AR) diseases derived from the largest publicly available cohort of the Chinese population were re-analysed using a bottom-up approach. The analyses of gene carrier rates (GCRs) across distinct ethnic groups revealed that substantial heterogeneity existed potentially due to diverse evolutionary selection. The sampling population, sequencing coverage and underlying population structure contributed to the differential variants observed between ChinaMAP and the East Asian group in gnomAD. Beyond characteristics of GCR, potential druggable targets were additionally explored according to genomic features and functional roles of investigated genes, demonstrating that phase separation could be a therapeutic target for autosomal recessive diseases. A further examination of estimated GCR across ethnic groups indicated that most genes shared by at least two populations could be utilised to direct the design of a pan-ethnic screening application once sequencing and interpreting costs become negligible. To this end, a list of autosomal recessive disease genes is proposed based on the prioritised rank of GCR to formulate a tiered screening strategy.

KiT-GENIE, the French genetic biobank of kidney transplantation

KiT-GENIE is a monocentric DNA biobank set up to consolidate the very rich and homogeneous DIVAT French cohort of kidney donors and recipients (D/R) in order to explore the molecular factors involved in kidney transplantation outcomes. We collected DNA samples for kidney transplantations performed in Nantes, and we leveraged GWAS genotyping data for securing high-quality genetic data with deep SNP and HLA annotations through imputations and for inferring D/R genetic ancestry. Overall, the biobank included 4217 individuals (n = 1945 D + 2,272 R, including 1969 D/R pairs), 7.4 M SNPs and over 200 clinical variables. KiT-GENIE represents an accurate snapshot of kidney transplantation clinical practice in Nantes between 2002 and 2018, with an enrichment in living kidney donors (17%) and recipients with focal segmental glomerulosclerosis (4%). Recipients were predominantly male (63%), of European ancestry (93%), with a mean age of 51yo and 86% experienced their first graft over the study period. D/R pairs were 93% from European ancestry, and 95% pairs exhibited at least one HLA allelic mismatch. The mean follow-up time was 6.7 years with a hindsight up to 25 years. Recipients experienced biopsy-proven rejection and graft loss for 16.6% and 21.3%, respectively. KiT-GENIE constitutes one of the largest kidney transplantation genetic cohorts worldwide to date. It includes homogeneous high-quality clinical and genetic data for donors and recipients, hence offering a unique opportunity to investigate immunogenetic and genetic factors, as well as donor-recipient interactions and mismatches involved in rejection, graft survival, primary disease recurrence and other comorbidities.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

A Large-Scale Exome-Wide Association Study Identifies Novel Germline Mutations in Lung Cancer

Rationale: Genome-wide association studies have identified common variants of lung cancer. However, the contribution of rare exome-wide variants, especially protein-coding variants, to cancers remains largely unexplored. Objectives: To evaluate the role of human exomes in genetic predisposition to lung cancer. Methods: We performed exome-wide association studies to detect the association of exomes with lung cancer in 30,312 patients and 652,902 control subjects. A scalable and accurate implementation of a generalized mixed model was used to detect the association signals for loss-of-function, missense, and synonymous variants and gene-level sets. Furthermore, we performed association and Bayesian colocalization analyses to evaluate their relationships with intermediate exposures. Measurements and Main Results: We systematically analyzed 216,739 single-nucleotide variants in the human exome. The loss-of-function variants exhibited the most notable effects on lung cancer risk. We identified four novel variants, including two missense variants (rs202197044TET3 [Pmeta (P values of meta-analysis) = 3.60 × 10-8] and rs202187871POT1 [Pmeta = 2.21 × 10-8]) and two synonymous variants (rs7447927TMEM173 [Pmeta = 1.32 × 10-9] and rs140624366ATRN [Pmeta = 2.97 × 10-9]). rs202197044TET3 was significantly associated with emphysema (odds ratio, 3.55; Pfdr = 0.015), whereas rs7447927POT1 was strongly associated with telomere length (β = 1.08; Pfdr (FDR corrected P value) = 3.76 × 10-53). Functional evidence of expression of quantitative trait loci, splicing quantitative trait loci, and isoform expression was found for the four novel genes. Gene-level association tests identified several novel genes, including POT1 (protection of telomeres 1), RTEL1, BSG, and ZNF232. Conclusions: Our findings provide insights into the genetic architecture of human exomes and their role in lung cancer predisposition.

Microglial efferocytosis: Diving into the Alzheimer's disease gene pool

Genome-wide association studies and functional genomics studies have linked specific cell types, genes, and pathways to Alzheimer's disease (AD) risk. In particular, AD risk alleles primarily affect the abundance or structure, and thus the activity, of genes expressed in macrophages, strongly implicating microglia (the brain-resident macrophages) in the etiology of AD. These genes converge on pathways (endocytosis/phagocytosis, cholesterol metabolism, and immune response) with critical roles in core macrophage functions such as efferocytosis. Here, we review these pathways, highlighting relevant genes identified in the latest AD genetics and genomics studies, and describe how they may contribute to AD pathogenesis. Investigating the functional impact of AD-associated variants and genes in microglia is essential for elucidating disease risk mechanisms and developing effective therapeutic approaches.