Assisting the analysis of insertions and deletions using regional allele frequencies.

Accurate estimation of population allele frequency (AF) is crucial for gene discovery and genetic diagnostics. However, determining AF for frameshift-inducing small insertions and deletions (indels) faces challenges due to discrepancies in mapping and variant calling methods. Here, we propose an innovative approach to assess indel AF. We developed CRAFTS-indels (Calculating Regional Allele Frequency Targeting Small indels), an algorithm that combines AF of distinct indels within a given region and provides "regional AF" (rAF). We tested and validated CRAFTS-indels using three independent datasets: gnomAD v2 (n=125,748 samples), an internal dataset (IGM; n=39,367), and the UK BioBank (UKBB; n=469,835). By comparing rAF against standard AF, we identified rare indels with rAF exceeding standard AF (sAF≤10-4 and rAF>10-4) as "rAF-hi" indels. Notably, a high percentage of rare indels were "rAF-hi", with a higher proportion in gnomAD v2 (11-20%) and IGM (11-22%) compared to the UKBB (5-9% depending on the CRAFTS-indels' parameters). Analysis of the overlap of regions based on their rAF with low complexity regions and with ClinVar classification supported the pertinence of rAF. Using the internal dataset, we illustrated the utility of CRAFTS-indel in the analysis of de novo variants and the potential negative impact of rAF-hi indels in gene discovery. In summary, annotation of indels with cohort specific rAF can be used to handle some of the limitations of current annotation pipelines and facilitate detection of novel gene disease associations. CRAFTS-indels offers a user-friendly approach to providing rAF annotation. It can be integrated into public databases such as gnomAD, UKBB and used by ClinVar to revise indel classifications.

The expanded spectrum of human disease associated with GREB1L likely includes complex congenital heart disease.

OBJECTIVE: GREB1L has been linked prenatally to Potter's sequence, as well as less severe anomalies of the kidney, uterus, inner ear, and heart. The full phenotypic spectrum is unknown. The purpose of this study was to characterize known and novel pre- and postnatal phenotypes associated with GREB1L. METHODS: We solicited cases from the Fetal Sequencing Consortium, screened a population-based genomic database, and conducted a comprehensive literature search to identify disease cases associated with GREB1L. We present a detailed phenotypic spectrum and molecular changes. RESULTS: One hundred twenty-seven individuals with 51 unique pathogenic or likely pathogenic GREB1L variants were identified. 24 (47%) variants were associated with isolated kidney anomalies, 19 (37%) with anomalies of multiple systems, including one case of hypoplastic left heart syndrome, five (10%) with isolated sensorineural hearing loss, two (4%) with isolated uterine agenesis; and one (2%) with isolated tetralogy of Fallot. CONCLUSION: GREB1L may cause complex congenital heart disease (CHD) in humans. Clinicians should consider GREB1L testing in the setting of CHD, and cardiac screening in the setting of GREB1L variants.

Clinical and Genetic Characteristics of CKD Patients with High-Risk APOL1 Genotypes.

SIGNIFICANCE STATEMENT: APOL1 high-risk genotypes confer a significant risk of kidney disease, but variability in patient outcomes suggests the presence of modifiers of the APOL1 effect. We show that a diverse population of CKD patients with high-risk APOL1 genotypes have an increased lifetime risk of kidney failure and higher eGFR decline rates, with a graded risk among specific high-risk genotypes. CKD patients with high-risk APOL1 genotypes have a lower diagnostic yield for monogenic kidney disease. Exome sequencing revealed enrichment of rare missense variants within the inflammasome pathway modifying the effect of APOL1 risk genotypes, which may explain some clinical heterogeneity. BACKGROUND: APOL1 genotype has significant effects on kidney disease development and progression that vary among specific causes of kidney disease, suggesting the presence of effect modifiers. METHODS: We assessed the risk of kidney failure and the eGFR decline rate in patients with CKD carrying high-risk ( N =239) and genetically matched low-risk ( N =1187) APOL1 genotypes. Exome sequencing revealed monogenic kidney diseases. Exome-wide association studies and gene-based and gene set-based collapsing analyses evaluated genetic modifiers of the effect of APOL1 genotype on CKD. RESULTS: Compared with genetic ancestry-matched patients with CKD with low-risk APOL1 genotypes, those with high-risk APOL1 genotypes had a higher risk of kidney failure (Hazard Ratio [HR]=1.58), a higher decline in eGFR (6.55 versus 3.63 ml/min/1.73 m 2 /yr), and were younger at time of kidney failure (45.1 versus 53.6 years), with the G1/G1 genotype demonstrating the highest risk. The rate for monogenic kidney disorders was lower among patients with CKD with high-risk APOL1 genotypes (2.5%) compared with those with low-risk genotypes (6.7%). Gene set analysis identified an enrichment of rare missense variants in the inflammasome pathway in individuals with high-risk APOL1 genotypes and CKD (odds ratio=1.90). CONCLUSIONS: In this genetically matched cohort, high-risk APOL1 genotypes were associated with an increased risk of kidney failure and eGFR decline rate, with a graded risk between specific high-risk genotypes and a lower rate of monogenic kidney disease. Rare missense variants in the inflammasome pathway may act as genetic modifiers of APOL1 effect on kidney disease.

Genomic Disorders in CKD across the Lifespan.

SIGNIFICANCE STATEMENT: Pathogenic structural genetic variants, also known as genomic disorders, have been associated with pediatric CKD. This study extends those results across the lifespan, with genomic disorders enriched in both pediatric and adult patients compared with controls. In the Chronic Renal Insufficiency Cohort study, genomic disorders were also associated with lower serum Mg, lower educational performance, and a higher risk of death. A phenome-wide association study confirmed the link between kidney disease and genomic disorders in an unbiased way. Systematic detection of genomic disorders can provide a molecular diagnosis and refine prediction of risk and prognosis. BACKGROUND: Genomic disorders (GDs) are associated with many comorbid outcomes, including CKD. Identification of GDs has diagnostic utility. METHODS: We examined the prevalence of GDs among participants in the Chronic Kidney Disease in Children (CKiD) cohort II ( n =248), Chronic Renal Insufficiency Cohort (CRIC) study ( n =3375), Columbia University CKD Biobank (CU-CKD; n =1986), and the Family Investigation of Nephropathy and Diabetes (FIND; n =1318) compared with 30,746 controls. We also performed a phenome-wide association analysis (PheWAS) of GDs in the electronic MEdical Records and GEnomics (eMERGE; n =11,146) cohort. RESULTS: We found nine out of 248 (3.6%) CKiD II participants carried a GD, replicating prior findings in pediatric CKD. We also identified GDs in 72 out of 6679 (1.1%) adult patients with CKD in the CRIC, CU-CKD, and FIND cohorts, compared with 199 out of 30,746 (0.65%) GDs in controls (OR, 1.7; 95% CI, 1.3 to 2.2). Among adults with CKD, we found recurrent GDs at the 1q21.1, 16p11.2, 17q12, and 22q11.2 loci. The 17q12 GD (diagnostic of renal cyst and diabetes syndrome) was most frequent, present in 1:252 patients with CKD and diabetes. In the PheWAS, dialysis and neuropsychiatric phenotypes were the top associations with GDs. In CRIC participants, GDs were associated with lower serum magnesium, lower educational achievement, and higher mortality risk. CONCLUSION: Undiagnosed GDs are detected both in children and adults with CKD. Identification of GDs in these patients can enable a precise genetic diagnosis, inform prognosis, and help stratify risk in clinical studies. GDs could also provide a molecular explanation for nephropathy and comorbidities, such as poorer neurocognition for a subset of patients.

Education and electronic medical records and genomics network, challenges, and lessons learned from a large-scale clinical trial using polygenic risk scores.

Polygenic risk scores (PRS) have potential to improve health care by identifying individuals that have elevated risk for common complex conditions. Use of PRS in clinical practice, however, requires careful assessment of the needs and capabilities of patients, providers, and health care systems. The electronic Medical Records and Genomics (eMERGE) network is conducting a collaborative study which will return PRS to 25,000 pediatric and adult participants. All participants will receive a risk report, potentially classifying them as high risk (∼2-10% per condition) for 1 or more of 10 conditions based on PRS. The study population is enriched by participants from racial and ethnic minority populations, underserved populations, and populations who experience poorer medical outcomes. All 10 eMERGE clinical sites conducted focus groups, interviews, and/or surveys to understand educational needs among key stakeholders-participants, providers, and/or study staff. Together, these studies highlighted the need for tools that address the perceived benefit/value of PRS, types of education/support needed, accessibility, and PRS-related knowledge and understanding. Based on findings from these preliminary studies, the network harmonized training initiatives and formal/informal educational resources. This paper summarizes eMERGE's collective approach to assessing educational needs and developing educational approaches for primary stakeholders. It discusses challenges encountered and solutions provided.

Mouse and human studies support DSTYK loss of function as a low-penetrance and variable expressivity risk factor for congenital urinary tract anomalies.

PURPOSE: Previous work identified rare variants in DSTYK associated with human congenital anomalies of the kidney and urinary tract (CAKUT). Here, we present a series of mouse and human studies to clarify the association, penetrance, and expressivity of DSTYK variants. METHODS: We phenotypically characterized Dstyk knockout mice of 3 separate inbred backgrounds and re-analyzed the original family segregating the DSTYK c.654+1G>A splice-site variant (referred to as "SSV" below). DSTYK loss of function (LOF) and SSVs were annotated in individuals with CAKUT, epilepsy, or amyotrophic lateral sclerosis vs controls. A phenome-wide association study analysis was also performed using United Kingdom Biobank (UKBB) data. RESULTS: Results demonstrate ∼20% to 25% penetrance of obstructive uropathy, at least, in C57BL/6J and FVB/NJ Dstyk-/- mice. Phenotypic penetrance increased to ∼40% in C3H/HeJ mutants, with mild-to-moderate severity. Re-analysis of the original family segregating the rare SSV showed low penetrance (43.8%) and no alternative genetic causes for CAKUT. LOF DSTYK variants burden showed significant excess for CAKUT and epilepsy vs controls and an exploratory phenome-wide association study supported association with neurological disorders. CONCLUSION: These data support causality for DSTYK LOF variants and highlights the need for large-scale sequencing studies (here >200,000 cases) to accurately assess causality for genes and variants to lowly penetrant traits with common population prevalence.

Kidney disease progression and collider bias in GWAS.

New genome-wide meta-analysis for longitudinal kidney function decline identified several genetic loci related to kidney disease progression. The study illustrated the complexity of modeling longitudinal traits in genome-wide association studies and highlighted the issue of a collider bias that can be introduced when a kidney disease progression phenotype is adjusted for baseline kidney function. Herein, we briefly outline the key findings of this study, their limitations, and implications for future studies.

Medical Records-Based Genetic Studies of the Complement System.

BACKGROUND: Genetic variants in complement genes have been associated with a wide range of human disease states, but well-powered genetic association studies of complement activation have not been performed in large multiethnic cohorts. METHODS: We performed medical records-based genome-wide and phenome-wide association studies for plasma C3 and C4 levels among participants of the Electronic Medical Records and Genomics (eMERGE) network. RESULTS: In a GWAS for C3 levels in 3949 individuals, we detected two genome-wide significant loci: chr.1q31.3 (CFH locus; rs3753396-A; ß=0.20; 95% CI, 0.14 to 0.25; P=1.52x10-11) and chr.19p13.3 (C3 locus; rs11569470-G; ß=0.19; 95% CI, 0.13 to 0.24; P=1.29x10-8). These two loci explained approximately 2% of variance in C3 levels. GWAS for C4 levels involved 3998 individuals and revealed a genome-wide significant locus at chr.6p21.32 (C4 locus; rs3135353-C; ß=0.40; 95% CI, 0.34 to 0.45; P=4.58x10-35). This locus explained approximately 13% of variance in C4 levels. The multiallelic copy number variant analysis defined two structural genomic C4 variants with large effect on blood C4 levels: C4-BS (ß=-0.36; 95% CI, -0.42 to -0.30; P=2.98x10-22) and C4-AL-BS (ß=0.25; 95% CI, 0.21 to 0.29; P=8.11x10-23). Overall, C4 levels were strongly correlated with copy numbers of C4A and C4B genes. In comprehensive phenome-wide association studies involving 102,138 eMERGE participants, we cataloged a full spectrum of autoimmune, cardiometabolic, and kidney diseases genetically related to systemic complement activation. CONCLUSIONS: We discovered genetic determinants of plasma C3 and C4 levels using eMERGE genomic data linked to electronic medical records. Genetic variants regulating C3 and C4 levels have large effects and multiple clinical correlations across the spectrum of complement-related diseases in humans.

Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux.

BACKGROUND: Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood. METHODS: A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry. RESULTS: Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis. CONCLUSIONS: These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

Experimental evidence of pathogenic role of IgG autoantibodies in IgA nephropathy.

BACKGROUND: IgA nephropathy is thought to be an autoimmune disease wherein galactose-deficient IgA1 (Gd-IgA1) is recognized by IgG autoantibodies, resulting in formation and renal accumulation of nephritogenic immune complexes. Although this hypothesis is supported by recent findings that, in renal immunodeposits of IgA nephropathy patients, IgG is enriched for Gd-IgA1-specific autoantibodies, experimental proof is still lacking. METHODS: IgG isolated from sera of IgA nephropathy patients or produced as a recombinant IgG (rIgG) was mixed with human Gd-IgA1 to form immune complexes. IgG from healthy individuals served as a control. Nude and SCID mice were injected with human IgG and Gd-IgA1, in immune complexes or individually, and their presence in kidneys was ascertained by immunofluorescence. Pathologic changes in the glomeruli were evaluated by quantitative morphometry and exploratory transcriptomic profiling was performed by RNA-Seq. RESULTS: Immunodeficient mice injected with Gd-IgA1 mixed with IgG autoantibodies from patients with IgA nephropathy, but not Gd-IgA1 mixed with IgG from healthy individuals, displayed IgA, IgG, and mouse complement C3 glomerular deposits and mesangioproliferative glomerular injury with hematuria and proteinuria. Un-complexed Gd-IgA1 or IgG did not induce pathological changes. Moreover, Gd-IgA1-rIgG immune complexes injected into immunodeficient mice induced histopathological changes characteristic of human disease. Exploratory transcriptome profiling of mouse kidney tissues indicated that these immune complexes altered gene expression of multiple pathways, in concordance with the changes observed in kidney biopsies of patients with IgA nephropathy. CONCLUSIONS: This study provides the first in vivo evidence for a pathogenic role of IgG autoantibodies specific for Gd-IgA1 in the pathogenesis of IgA nephropathy.

iMEGES: integrated mental-disorder GEnome score by deep neural network for prioritizing the susceptibility genes for mental disorders in personal genomes.

BACKGROUND: A range of rare and common genetic variants have been discovered to be potentially associated with mental diseases, but many more have not been uncovered. Powerful integrative methods are needed to systematically prioritize both variants and genes that confer susceptibility to mental diseases in personal genomes of individual patients and to facilitate the development of personalized treatment or therapeutic approaches. METHODS: Leveraging deep neural network on the TensorFlow framework, we developed a computational tool, integrated Mental-disorder GEnome Score (iMEGES), for analyzing whole genome/exome sequencing data on personal genomes. iMEGES takes as input genetic mutations and phenotypic information from a patient with mental disorders, and outputs the rank of whole genome susceptibility variants and the prioritized disease-specific genes for mental disorders by integrating contributions from coding and non-coding variants, structural variants (SVs), known brain expression quantitative trait loci (eQTLs), and epigenetic information from PsychENCODE. RESULTS: iMEGES was evaluated on multiple datasets of mental disorders, and it achieved improved performance than competing approaches when large training dataset is available. CONCLUSION: iMEGES can be used in population studies to help the prioritization of novel genes or variants that might be associated with the susceptibility to mental disorders, and also on individual patients to help the identification of genes or variants related to mental diseases.

Participant-guided development of bilingual genomic educational infographics for Electronic Medical Records and Genomics Phase IV study.

OBJECTIVE: Developing targeted, culturally competent educational materials is critical for participant understanding of engagement in a large genomic study that uses computational pipelines to produce genome-informed risk assessments. MATERIALS AND METHODS: Guided by the Smerecnik framework that theorizes understanding of multifactorial genetic disease through 3 knowledge types, we developed English and Spanish infographics for individuals enrolled in the Electronic Medical Records and Genomics Network. Infographics were developed to explain concepts in lay language and visualizations. We conducted iterative sessions using a modified "think-aloud" process with 10 participants (6 English, 4 Spanish-speaking) to explore comprehension of and attitudes towards the infographics. RESULTS: We found that all but one participant had "awareness knowledge" of genetic disease risk factors upon viewing the infographics. Many participants had difficulty with "how-to" knowledge of applying genetic risk factors to specific monogenic and polygenic risks. Participant attitudes towards the iteratively-refined infographics indicated that design saturation was reached. DISCUSSION: There were several elements that contributed to the participants' comprehension (or misunderstanding) of the infographics. Visualization and iconography techniques best resonated with those who could draw on prior experiences or knowledge and were absent in those without. Limited graphicacy interfered with the understanding of absolute and relative risks when presented in graph format. Notably, narrative and storytelling theory that informed the creation of a vignette infographic was most accessible to all participants. CONCLUSION: Engagement with the intended audience who can identify strengths and points for improvement of the intervention is necessary to the development of effective infographics.

Cell subtype-specific effects of genetic variation in the Alzheimer's disease brain.

The relationship between genetic variation and gene expression in brain cell types and subtypes remains understudied. Here, we generated single-nucleus RNA sequencing data from the neocortex of 424 individuals of advanced age; we assessed the effect of genetic variants on RNA expression in cis (cis-expression quantitative trait loci) for seven cell types and 64 cell subtypes using 1.5 million transcriptomes. This effort identified 10,004 eGenes at the cell type level and 8,099 eGenes at the cell subtype level. Many eGenes are only detected within cell subtypes. A new variant influences APOE expression only in microglia and is associated with greater cerebral amyloid angiopathy but not Alzheimer's disease pathology, after adjusting for APOEε4, providing mechanistic insights into both pathologies. Furthermore, only a TMEM106B variant affects the proportion of cell subtypes. Integration of these results with genome-wide association studies highlighted the targeted cell type and probable causal gene within Alzheimer's disease, schizophrenia, educational attainment and Parkinson's disease loci.

Risk factors affecting polygenic score performance across diverse cohorts.

Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed effects of covariate stratification and interaction on body mass index (BMI) PGS (PGSBMI) across four cohorts of European (N=491,111) and African (N=21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R2 being nearly double between best and worst performing quintiles for certain covariates. 28 covariates had significant PGSBMI-covariate interaction effects, modifying PGSBMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R2 differences among strata and interaction effects - across all covariates, their main effects on BMI were correlated with their maximum R2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGSBMI individuals have highest R2 and increase in PGS effect. Using quantile regression, we show the effect of PGSBMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGSBMI performance and effects, we investigated ways to increase model performance taking into account non-linear effects. Machine learning models (neural networks) increased relative model R2 (mean 23%) across datasets. Finally, creating PGSBMI directly from GxAge GWAS effects increased relative R2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGSBMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.

Clinical associations with a polygenic predisposition to benign lower white blood cell counts.

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30-0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69-0.88], p = 4.0 × 10-5) or immunosuppressant (n = 354, HR = 0.61 [0.38-0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44-0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

BIGKnock: fine-mapping gene-based associations via knockoff analysis of biobank-scale data.

We propose BIGKnock (BIobank-scale Gene-based association test via Knockoffs), a computationally efficient gene-based testing approach for biobank-scale data, that leverages long-range chromatin interaction data, and performs conditional genome-wide testing via knockoffs. BIGKnock can prioritize causal genes over proxy associations at a locus. We apply BIGKnock to the UK Biobank data with 405,296 participants for multiple binary and quantitative traits, and show that relative to conventional gene-based tests, BIGKnock produces smaller sets of significant genes that contain the causal gene(s) with high probability. We further illustrate its ability to pinpoint potential causal genes at [Formula: see text] of the associated loci.

CARMA is a new Bayesian model for fine-mapping in genome-wide association meta-analyses.

Fine-mapping is commonly used to identify putative causal variants at genome-wide significant loci. Here we propose a Bayesian model for fine-mapping that has several advantages over existing methods, including flexible specification of the prior distribution of effect sizes, joint modeling of summary statistics and functional annotations and accounting for discrepancies between summary statistics and external linkage disequilibrium in meta-analyses. Using simulations, we compare performance with commonly used fine-mapping methods and show that the proposed model has higher power and lower false discovery rate (FDR) when including functional annotations, and higher power, lower FDR and higher coverage for credible sets in meta-analyses. We further illustrate our approach by applying it to a meta-analysis of Alzheimer's disease genome-wide association studies where we prioritize putatively causal variants and genes.