Chromosome 10q24.32 Variants Associate With Brain Arterial Diameters in Diverse Populations: A Genome-Wide Association Study

BACKGROUND

Brain arterial diameters (BADs) are novel imaging biomarkers of cerebrovascular disease, cognitive decline, and dementia. Traditional vascular risk factors have been associated with BADs, but whether there may be genetic determinants of BADs is unknown.

METHODS AND RESULTS

The authors studied 4150 participants from 6 geographically diverse population-based cohorts (40% European, 14% African, 22% Hispanic, 24% Asian ancestries). Brain arterial diameters for 13 segments were measured and averaged to obtain a global measure of BADs as well as the posterior and anterior circulations. A genome-wide association study revealed 14 variants at one locus associated with global BAD at genome-wide significance (P<5×10-8) (top single-nucleotide polymorphism, rs7921574; β=0.06 [P=1.54×10-8]). This locus mapped to an intron of CNNM2. A trans-ancestry genome-wide association study meta-analysis identified 2 more loci at NT5C2 (rs10748839; P=2.54×10-8) and AS3MT (rs10786721; P=4.97×10-8), associated with global BAD. In addition, 2 single-nucleotide polymorphisms colocalized with expression of CNNM2 (rs7897654; β=0.12 [P=6.17×10-7]) and AL356608.1 (rs10786719; β=-0.17 [P=6.60×10-6]) in brain tissue. For the posterior BAD, 2 variants at one locus mapped to an intron of TCF25 were identified (top single-nucleotide polymorphism, rs35994878; β=0.11 [P=2.94×10-8]). For the anterior BAD, one locus at ADAP1 was identified in trans-ancestry genome-wide association analysis (rs34217249; P=3.11×10-8).

CONCLUSIONS

The current study reveals 3 novel risk loci (CNNM2, NT5C2, and AS3MT) associated with BADs. These findings may help elucidate the mechanism by which BADs may influence cerebrovascular health.

Admixture mapping identifies novel Alzheimer's disease risk regions in African Americans

BACKGROUND

This study used admixture mapping to prioritize the genetic regions associated with Alzheimer's disease (AD) in African American (AA) individuals, followed by ancestry-aware regression analysis to fine-map the prioritized regions.

METHODS

We analyzed 10,271 individuals from 17 different AA datasets. We performed admixture mapping and meta-analyzed the results. We then used regression analysis, adjusting for local ancestry main effects and interactions with genotype, to refine the regions identified from admixture mapping. Finally, we leveraged in silico annotation and differential gene expression data to prioritize AD-related variants and genes.

RESULTS

Admixture mapping identified two genome-wide significant loci on chromosomes 17p13.2 (p = 2.2 × 10-5 ) and 18q21.33 (p = 1.2 × 10-5 ). Our fine mapping of the chromosome 17p13.2 and 18q21.33 regions revealed several interesting genes such as the MINK1, KIF1C, and BCL2.

DISCUSSION

Our ancestry-aware regression approach showed that AA individuals have a lower risk of AD if they inherited African ancestry admixture block at the 17p13.2 locus.

HIGHLIGHTS

We identified two genome-wide significant admixture mapping signals: on chromosomes 17p13.2 and 18q21.33, which are novel in African American (AA) populations. Our ancestry-aware regression approach showed that AA individuals have a lower risk of Alzheimer's disease (AD) if they inherited African ancestry admixture block at the 17p13.2 locus. We found that the overall proportion of African ancestry does not differ between the cases and controls that suggest African genetic ancestry alone is not likely to explain the AD prevalence difference between AA and non-Hispanic White populations.

Using source-associated mobile genetic elements to identify zoonotic extraintestinal E. coli infections

A one-health perspective may provide new and actionable information about Escherichia coli transmission. E. coli colonizes a broad range of vertebrates, including humans and food-production animals, and is a leading cause of bladder, kidney, and bloodstream infections in humans. Substantial evidence supports foodborne transmission of pathogenic E. coli strains from food animals to humans. However, the relative contribution of foodborne zoonotic E. coli (FZEC) to the human extraintestinal disease burden and the distinguishing characteristics of such strains remain undefined. Using a comparative genomic analysis of a large collection of contemporaneous, geographically-matched clinical and meat-source E. coli isolates (n = 3111), we identified 17 source-associated mobile genetic elements - predominantly plasmids and bacteriophages - and integrated them into a novel Bayesian latent class model to predict the origins of clinical E. coli isolates. We estimated that approximately 8 % of human extraintestinal E. coli infections (mostly urinary tract infections) in our study population were caused by FZEC. FZEC strains were equally likely to cause symptomatic disease as non-FZEC strains. Two FZEC lineages, ST131-H22 and ST58, appeared to have particularly high virulence potential. Our findings imply that FZEC strains collectively cause more urinary tract infections than does any single non-E. coli uropathogenic species (e.g., Klebsiella pneumoniae). Our novel approach can be applied in other settings to identify the highest-risk FZEC strains, determine their sources, and inform new one-health strategies to decrease the heavy public health burden imposed by extraintestinal E. coli infections.

Chromosome 10q24.32 Variants Associate with Brain Arterial Diameters in Diverse Populations: A Genome-Wide Association Study

Background

Brain arterial diameters are novel imaging biomarkers of cerebrovascular disease, cognitive decline and dementia. Traditional vascular risk factors have been associated with brain arterial diameters but whether there may be genetic determinants of brain arterial diameters is unknown.

Results

We studied 4150 participants from six geographically diverse population-based cohorts (40% European, 14% African, 22% Hispanic, 24% Asian ancestries). We measured brain arterial diameters for 13 segments and averaged them to obtain a global measure of brain arterial diameters as well as the posterior and anterior circulations. A genome-wide association study (GWAS) revealed 14 variants at one locus associated with global brain arterial diameter at genome-wide significance (P<5×10-8) (top SNP, rs7921574; β =0.06, P=1.54×10-8). This locus mapped to an intron of CNNM2. A trans-ancestry GWAS meta-analysis identified two more loci at NT5C2 (rs10748839; P=2.54×10-8) and at AS3MT (rs10786721; P=4.97×10-8), associated with global brain arterial diameter. In addition, two SNPs co-localized with expression of CNNM2 (rs7897654, β=0.12, P=6.17×10-7) and AL356608.1 (rs10786719, β =-0.17, P=6.60×10-6) in brain tissue. For the posterior brain arterial diameter, two variants at one locus mapped to an intron of TCF25 were identified (top SNP, rs35994878; β =0.11, P=2.94×10-8). For the anterior brain arterial diameter, one locus at ADAP1 was identified in trans-ancestry genome-wide association analysis (rs34217249; P=3.11×10-8).

Conclusion

Our study reveals three novel risk loci (CNNM2, NT5C2 and AS3MT) associated with brain arterial diameters. Our finding may elucidate the mechanisms by which brain arterial diameters influence the risk of stroke and dementia.

Integration of GWAS and brain transcriptomic analyses in a multiethnic sample of 35,245 older adults identifies DCDC2 gene as predictor of episodic memory maintenance

Identifying genes underlying memory function will help characterize cognitively resilient and high-risk declining subpopulations contributing to precision medicine strategies. We estimated episodic memory trajectories in 35,245 ethnically diverse older adults representing eight independent cohorts. We conducted apolipoprotein E (APOE)-stratified genome-wide association study (GWAS) analyses and combined individual cohorts' results via meta-analysis. Three independent transcriptomics datasets were used to further interpret GWAS signals. We identified DCDC2 gene significantly associated with episodic memory (Pmeta = 3.3 x 10^-8 ) among non-carriers of APOE ε4 (N = 24,941). Brain transcriptomics revealed an association between episodic memory maintenance and (1) increased dorsolateral prefrontal cortex DCDC2 expression (P = 3.8 x 10^-4 ) and (2) lower burden of pathological Alzheimer's disease (AD) hallmarks (paired helical fragment tau P = .003, and amyloid beta load P = .008). Additional transcriptomics results comparing AD and cognitively healthy brain samples showed a downregulation of DCDC2 levels in superior temporal gyrus (P = .007) and inferior frontal gyrus (P = .013). Our work identified DCDC2 gene as a novel predictor of memory maintenance.

FMNL2 regulates gliovascular interactions and is associated with vascular risk factors and cerebrovascular pathology in Alzheimer's disease

Alzheimer's disease (AD) has been associated with cardiovascular and cerebrovascular risk factors (CVRFs) during middle age and later and is frequently accompanied by cerebrovascular pathology at death. An interaction between CVRFs and genetic variants might explain the pathogenesis. Genome-wide, gene by CVRF interaction analyses for AD, in 6568 patients and 8101 controls identified FMNL2 (p = 6.6 × 10-7). A significant increase in FMNL2 expression was observed in the brains of patients with brain infarcts and AD pathology and was associated with amyloid and phosphorylated tau deposition. FMNL2 was also prominent in astroglia in AD among those with cerebrovascular pathology. Amyloid toxicity in zebrafish increased fmnl2a expression in astroglia with detachment of astroglial end feet from blood vessels. Knockdown of fmnl2a prevented gliovascular remodeling, reduced microglial activity and enhanced amyloidosis. APP/PS1dE9 AD mice also displayed increased Fmnl2 expression and reduced the gliovascular contacts independent of the gliotic response. Based on this work, we propose that FMNL2 regulates pathology-dependent plasticity of the blood-brain-barrier by controlling gliovascular interactions and stimulating the clearance of extracellular aggregates. Therefore, in AD cerebrovascular risk factors promote cerebrovascular pathology which in turn, interacts with FMNL2 altering the normal astroglial-vascular mechanisms underlying the clearance of amyloid and tau increasing their deposition in brain.

Admixture Mapping of Alzheimer's disease in Caribbean Hispanics identifies a new locus on 22q13.1

Late-onset Alzheimer's disease (LOAD) is significantly more frequent in Hispanics than in non-Hispanic Whites. Ancestry may explain these differences across ethnic groups. To this end, we studied a large cohort of Caribbean Hispanics (CH, N = 8813) and tested the association between Local Ancestry (LA) and LOAD ("admixture mapping") to identify LOAD-associated ancestral blocks, separately for ancestral components (European [EUR], African [AFR], Native American[NA]) and jointly (AFR + NA). Ancestral blocks significant after permutation were fine-mapped employing multi-ethnic whole-exome sequencing (WES) to identify rare variants associated with LOAD (SKAT-O) and replicated in the UK Biobank WES dataset. Candidate genes were validated studying (A) protein expression in human LOAD and control brains; (B) two animal AD models, Drosophila and Zebrafish. In the joint AFR + NA model, we identified four significant ancestral blocks located on chromosomes 1 (p value = 8.94E-05), 6 (p value = 8.63E-05), 21 (p value = 4.64E-05) and 22 (p value = 1.77E-05). Fine-mapping prioritized the GCAT gene on chromosome 22 (SKAT-O p value = 3.45E-05) and replicated in the UK Biobank (SKAT-O p value = 0.05). In LOAD brains, a decrease of 28% in GCAT protein expression was observed (p value = 0.038), and GCAT knockdown in Amyloid-β42 Drosophila exacerbated rough eye phenotype (68% increase, p value = 4.84E-09). In zebrafish, gcat expression increased after acute amyloidosis (34%, p value = 0.0049), and decreased upon anti-inflammatory Interleukin-4 (39%, p value = 2.3E-05). Admixture mapping uncovered genomic regions harboring new LOAD-associated loci that might explain the observed different frequency of LOAD across ethnic groups. Our results suggest that the inflammation-related activity of GCAT is a response to amyloid toxicity, and reduced GCAT expression exacerbates AD pathology.

Genetic determinants of intracranial large artery stenosis in the northern Manhattan study

BACKGROUND

Intracranial stenosis is one of the most common causes of stroke worldwide. Several single nucleotide polymorphisms have been associated with intracranial atherosclerosis, which is inferred to be the most common underlying cause of intracranial large artery stenosis (ILAS). We previously reviewed known genetic variants related to ILAS in predominantly Asian cohorts, but their prevalence and role in ILAS among western multiethnic populations are uncertain.

METHODS

We leveraged existing imaging and genetic data from the Northern Manhattan Study, a multiethnic prospective cohort study. Based on literature review, we selected adiponectin Q (ADIPOQ) rs2241767 and rs182052, ring finger protein 213 (RNF213) rs112735431, apolipoprotein E (APOE) rs429358, phosphodiesterase 4D (PDE4D) rs2910829, lipoprotein lipase (LPL) rs320, and aldosterone synthase (CYP11B2) rs1799998 variants as candidates to explore. We defined ILAS as luminal stenosis >50% in any intracranial large artery using time-of-flight magnetic resonance angiography (MRA).

RESULTS

We included 1109 participants (mean age 70 ± 9 years, 70% Hispanic, 60% women) in this study. ILAS was identified in 81 (7%) NOMAS participants. Logistic regression analyses adjusted for age, sex, principal components, and vascular risk factors showed ILAS prevalence associated with CYP11B2 rs1799998 under the dominant model (OR = 0.56, 95%CI: 0.35-0.89) and LPL rs320 heterozygote genotype (OR = 1.68, 95%CI: 1.05-2.71). The genotype distributions of ADIPOQ rs2241767 and rs182052, APOE rs429358 and CYP11B2 rs1799998 variants were significantly different among non-Hispanic white and Black, and Hispanic groups. When participants were further stratified by race/ethnicity, the estimates were consistent for CYP11B2 rs1799998 across race/ethnic groups but not for LPL rs320.

CONCLUSION

The CYP11B2 rs1799998 variant may be a protective genetic factor for ILAS across race/ethnic groups, but the risk of ILAS associated with LPL rs320 varies by race/ethnic group. Further functional studies may help elucidate the role that these variants play in the pathophysiology of ILAS.

Abstract WP205: Trans-ethnic Genome-wide Association Studies Meta-analysis Of Large Brain Arterial Diameters Identifies Two Novel Genetic Loci

Introduction: Larger brain arterial diameters, which may be a consequence of connective tissue disorders or hypertension, are associated with mortality, cognitive decline and risk of dementia. Whether there may be genetic determinants of large brain arterial diameters is unknown, but identifying such traits may shed light into possible mechanistic links between large brain arterial diameters and the observed outcomes.

Methods: We studied 3654 participants from four ethnically diverse population-based cohorts: The Atherosclerosis Risk in Communities study (n=1565: 74.3% white, 25.7% black), The Northern Manhattan Study (n=1092: 13.8% white, 16.6% black, 69.6% Hispanic), Epidemiology of Dementia in Singapore (n=647: 31.7% Chinese, 35.7% Malay, 32.6% Indian), and Memory Clinic Cohort in Singapore (n=350: 84.9% Chinese, 7.8% Malay, 7.3% Indian). We measured brain arterial diameters for 13 segments and normalized the averaged diameter per participant to obtain global brain arterial diameter as dependent variable. Genome-wide association studies (GWAS) was performed in each study using linear regression, adjusted for sex, age, head size, number of missing arteries and population substructure. Meta-analysis of association statistics was performed using the METAL to generate combined trait-specific effect estimates. We also performed trans-ethnic meta-regression by MR-MEGA to account for the heterogeneity in allelic effect that is correlated with ancestry.

Results: Genome-wide significant (P<5x10 -8 ) associations were observed for global brain arterial diameter with three SNPs: rs7921574(P=1.2x10 -8 ) mapping to an intron of CNNM2; rs2075858(P=2.47x10 -8 ) and rs2075857 (P=2.37x10 -8 ) mapping to introns of MUC5B. MR-MEGA analysis showed suggestive evidence of association (P=5.05x10 -8 ) for rs7921574. These loci comprise genes involved in response to magnesium transport, neuronal transmission, hypertension, mucociliary clearance and host defense.

Conclusion: Our study reveals two novel risk loci for large brain arterial diameters. We have begun to replicate these associations and to identify functional causal variants in independent studies. Identifying these loci may help us to elucidate underlying disease mechanisms.

KYNA/Ahr Signaling Suppresses Neural Stem Cell Plasticity and Neurogenesis in Adult Zebrafish Model of Alzheimer's Disease

Neurogenesis decreases in Alzheimer's disease (AD) patients, suggesting that restoring the normal neurogenic response could be a disease modifying intervention. To study the mechanisms of pathology-induced neuro-regeneration in vertebrate brains, zebrafish is an excellent model due to its extensive neural regeneration capacity. Here, we report that Kynurenic acid (KYNA), a metabolite of the amino acid tryptophan, negatively regulates neural stem cell (NSC) plasticity in adult zebrafish brain through its receptor, aryl hydrocarbon receptor 2 (Ahr2). The production of KYNA is suppressed after amyloid-toxicity through reduction of the levels of Kynurenine amino transferase 2 (KAT2), the key enzyme producing KYNA. NSC proliferation is enhanced by an antagonist for Ahr2 and is reduced with Ahr2 agonists or KYNA. A subset of Ahr2-expressing zebrafish NSCs do not express other regulatory receptors such as il4r or ngfra, indicating that ahr2-positive NSCs constitute a new subset of neural progenitors that are responsive to amyloid-toxicity. By performing transcriptome-wide association studies (TWAS) in three late onset Alzheimer disease (LOAD) brain autopsy cohorts, we also found that several genes that are components of KYNA metabolism or AHR signaling are differentially expressed in LOAD, suggesting a strong link between KYNA/Ahr2 signaling axis to neurogenesis in LOAD.

Polygenic Risk Score for Alzheimer's Disease in Caribbean Hispanics

OBJECTIVE

Polygenic risk scores (PRSs) assess the individual genetic propensity to a condition by combining sparse information scattered across genetic loci, often displaying small effect sizes. Most PRSs are constructed in European-ancestry populations, limiting their use in other ethnicities. Here we constructed and validated a PRS for late-onset Alzheimer's Disease (LOAD) in Caribbean Hispanics (CH).

METHODS

We used a CH discovery (n = 4,312) and independent validation sample (n = 1,850) to construct an ancestry-specific PRS ("CH-PRS") and evaluated its performance alone and with other predictors using the area under curve (AUC) and logistic regression (strength of association with LOAD and statistical significance). We tested if CH-PRS predicted conversion to LOAD in a subsample with longitudinal data (n = 1,239). We also tested the CH-PRS in an independent replication CH cohort (n = 200) and brain autopsy cohort (n = 33). Finally, we tested the effect of ancestry on PRS by using European and African American discovery cohorts to construct alternative PRSs ("EUR-PRS", "AA-PRS").

RESULTS

The full model (LOAD ~ CH-PRS + sex + age + APOE-ɛ4), achieved an AUC = 74% (ORCH-PRS  = 1.51 95%CI = 1.36-1.68), raising to >75% in APOE-ɛ4 non-carriers. CH-PRS alone achieved an AUC = 72% in the autopsy cohort, raising to AUC = 83% in full model. Higher CH-PRS was significantly associated with clinical LOAD in the replication CH cohort (OR = 1.61, 95%CI = 1.19-2.17) and significantly predicted conversion to LOAD (HR = 1.93, CI = 1.70-2.20) in the longitudinal subsample. EUR-PRS and AA-PRS reached lower prediction accuracy (AUC = 58% and 53%, respectively).

INTERPRETATION

Enriching diversity in genetic studies is critical to provide an effective PRS in profiling LOAD risk across populations. ANN NEUROL 2021;90:366-376.

1627. Outbreaks of Klebsiella pneumoniae in Special Care Nurseries (SCN) in Jamaica: Role of Whole-Genome Sequencing

Background

Klebsiella pneumoniae is a frequent cause of neonatal sepsis and carries a high mortality rate in lower and middle-income countries (LMICs). From March-November 2015, two Jamaican hospitals experienced K. pneumoniae outbreaks in their Special Care Nurseries (SCNs). New admissions to both SCNs were temporarily halted while additional infection control strategies were implemented. 31 babies were infected, of which 15 died. International collaboration was requested to help investigate if the sepsis cases were nosocomial transmission, repeated introductions from the community, or both using whole-genome sequencing

Methods

We sequenced DNA from 19 outbreak isolates (n = 13 from Hospital A, n = 6 from Hospital B) on an Illumina HiSeq2500 instrument and assembled short-reads using SPAdes. We used ResFinder v3.1.0 to screen resistance genes and assigned MLSTs using in-house scripts. To compare the outbreak isolates, we selected a reference genome from among the assembled isolates, aligned raw reads using the Burrows–Wheeler Aligner (BWA), identified SNPs using GATK UnifiedGenotyper, and removed the recombined regions using Gubbins v2.3.4. We further contextualized the 19 outbreak isolates against a global collection of more than 300 K. pneumoniae genomes.

Results

All 13 isolates from Hospital A appeared to be from a single source. All were ST45 and encoded blaCTX-M-15, which confers extended-spectrum β-lactam (ESBL) resistance. Five of 6 isolates from Hospital B appeared to be from a separate, single source. These 5 isolates were ST268 and susceptible to most antibiotics. 1 isolate from Hospital B was ST628, encoded blaCTX-M-15, and grouped separately from other Hospital B outbreak isolates. Hospital A and B outbreak isolates formed independent, unique clades within a global K. pneumoniae collection.

Conclusion

Our findings indicate nosocomial transmission was responsible for both neonatal K. pneumoniae outbreaks, rather than repeat introductions from the community. The main sequence types we detected (ST45 and ST268) are not known pandemic clones and may circulate regionally. Multifaceted infection control measures were implemented for effectively halting outbreaks.

Disclosures

All authors: No reported disclosures.

Association of Variants in PINX1 and TREM2 With Late-Onset Alzheimer Disease

Importance

Genetic causes of late-onset Alzheimer disease (LOAD) are not completely explained by known genetic loci. Whole-exome and whole-genome sequencing can improve the understanding of the causes of LOAD and provide initial steps required to identify potential therapeutic targets.

Objective

To identify the genetic loci for LOAD across different ethnic groups.

Design, Setting, and Participants

This multicenter cohort study was designed to analyze whole-exome sequencing data from a multiethnic cohort using a transethnic gene-kernel association test meta-analysis, adjusted for sex, age, and principal components, to identify genetic variants associated with LOAD. A meta-analysis was conducted on the results of 2 independent studies of whole-exome and whole-genome sequence data from individuals of European ancestry. This group of European American, African American, and Caribbean Hispanic individuals participating in an urban population-based study were the discovery cohort; the additional cohorts included affected individuals and control participants from 2 publicly available data sets. Replication was achieved using independent data sets from Caribbean Hispanic families with multiple family members affected by LOAD and the International Genetics of Alzheimer Project.

Main Outcomes and Measures

Late-onset Alzheimer disease.

Results

The discovery cohort included 3595 affected individuals, while the additional cohorts included 5931 individuals with LOAD and 5504 control participants. Of 3916 individuals in the discovery cohort, we included 3595 individuals (1397 with LOAD and 2198 cognitively healthy controls; 2451 [68.2%] women; mean [SD] age, 80.3 [6.83] years). Another 321 individuals (8.2%) were excluded because of non-LOAD diagnosis, age younger than 60 years, missing covariates, duplicate data, or genetic outlier status. Gene-based tests that compared affected individuals (n = 7328) and control participants (n = 7702) and included only rare and uncommon variants annotated as having moderate-high functional effect supported PINX1 (8p23.1) as a locus with gene-wide significance (P = 2.81 × 10-6) after meta-analysis across the 3 studies. The PINX1 finding was replicated using data from the family-based study and the International Genetics of Alzheimer Project. Full meta-analysis of discovery and replication cohorts reached a P value of 6.16 × 10-7 for PINX1 (in 7620 affected individuals vs 7768 control participants). We also identified TREM2 in an annotation model that prioritized highly deleterious variants with a combined annotation dependent depletion greater than 20 (P= 1.0 × 10-7).

Conclusions and Relevance

This gene-based, transethnic approach identified PINX1, a gene involved in telomere integrity, and TREM2, a gene with a product of an immune receptor found in microglia, as associated with LOAD. Both genes have well-established roles in aging and neurodegeneration.

Rare Variants Imputation in Admixed Populations: Comparison Across Reference Panels and Bioinformatics Tools

Background

Imputation has become a standard approach in genome-wide association studies (GWAS) to infer in silico untyped markers. Although feasibility for common variants imputation is well established, we aimed to assess rare and ultra-rare variants’ imputation in an admixed Caribbean Hispanic population (CH).

Methods

We evaluated imputation accuracy in CH (N = 1,000), focusing on rare (0.1% ≤ minor allele frequency (MAF) ≤ 1%) and ultra-rare (MAF < 0.1%) variants. We used two reference panels, the Haplotype Reference Consortium (HRC; N = 27,165) and 1000 Genome Project (1000G phase 3; N = 2,504) and multiple phasing (SHAPEIT, Eagle2) and imputation algorithms (IMPUTE2, MACH-Admix). To assess imputation quality, we reported: (a) high-quality variant counts according to imputation tools’ internal indexes (e.g., IMPUTE2 “Info” ≥ 80%). (b) Wilcoxon Signed-Rank Test comparing imputation quality for genotyped variants that were masked and imputed; (c) Cohen’s kappa coefficient to test agreement between imputed and whole-exome sequencing (WES) variants; (d) imputation of G206A mutation in the PSEN1 (ultra-rare in the general population an more frequent in CH) followed by confirmation genotyping. We also tested ancestry proportion (European, African and Native American) against WES-imputation mismatches in a Poisson regression fashion.

Results

SHAPEIT2 retrieved higher percentage of imputed high-quality variants than Eagle2 (rare: 51.02% vs. 48.60%; ultra-rare 0.66% vs. 0.65%, Wilcoxon p-value < 0.001). SHAPEIT-IMPUTE2 employing HRC outperformed 1000G (64.50% vs. 59.17%; 1.69% vs. 0.75% for high-quality rare and ultra-rare variants, respectively, Wilcoxon p-value < 0.001). SHAPEIT-IMPUTE2 outperformed MaCH-Admix. Compared to 1000G, HRC-imputation retrieved a higher number of high-quality rare and ultra-rare variants, despite showing lower agreement between imputed and WES variants (e.g., rare: 98.86% for HRC vs. 99.02% for 1000G). High Kappa (K = 0.99) was observed for both reference panels. Twelve G206A mutation carriers were imputed and all validated by confirmation genotyping. African ancestry was associated with higher imputation errors for uncommon and rare variants (p-value < 1e-05).

Conclusion

Reference panels with larger numbers of haplotypes can improve imputation quality for rare and ultra-rare variants in admixed populations such as CH. Ethnic composition is an important predictor of imputation accuracy, with higher African ancestry associated with poorer imputation accuracy.