Genome-wide association study identifies five risk loci for pernicious anemia and implicates the role of HLA-DR15 haplotype
===========================================================================================================================

* Triin Laisk
* Maarja Lepamets
* Reedik Mägi

## Abstract

Pernicious anemia is a rare condition characterized by vitamin B12 deficiency anemia due to lack of intrinsic factor, often caused by autoimmune gastritis. Patients with pernicious anemia have a higher incidence of other autoimmune disorders, such as type 1 diabetes, vitiligo and autoimmune thyroid issues. Therefore, the disease has a clear autoimmune basis, although the genetic susceptibility factors have thus far remained poorly studied. We conducted a genome-wide association study meta-analysis in 2,166 cases and 659,516 European controls from population-based biobanks and identified genome-wide significant signals in or near the *PTPN22 (*rs6679677, p=1.91 ⨯ 10−24, OR=1.63*), PNPT1* (rs12616502, p=3.14 ⨯ 10−8, OR=1.70), *HLA-DQB1* (rs28414666, p=1.40 ⨯ 10−16, OR=1.38), *IL2RA* (rs2476491, p=1.90 ⨯ 10−8, OR=1.22) and *AIRE* (rs74203920, p=2.33 ⨯ 10−9, OR=1.83) genes, thus providing the first robust associations between pernicious anemia and genetic risk factors. We further mapped the susceptibility in the HLA region to the HLA-DR15 haplotype. Analysis of associated diagnoses and disease trajectories confirm the association between pernicious anemia and thyroid issues, vitiligo, gastritis, stomach cancer, osteoporosis and other diagnoses.

Keywords
*   genome-wide association study
*   pernicious anemia
*   B12 deficiency

## Main text

B12 deficiency anemia due to intrinsic factor deficiency, also known as pernicious anemia, is characterized by impaired B12 uptake caused by lack of intrinsic factor, a substance produced by epithelial cells of the stomach lining. Intrinsic factor normally binds B12 and facilitates absorption in the intestinal tract. Pernicious anemia is often caused by autoimmune damage to the stomach lining (autoimmune gastritis) in which case the gastric epithelial lining is damaged or destroyed1. The prevalence of pernicious anemia is around 0.1% in populations of European ancestry; however, it is more common in older people (∼2% in >60-year-olds), and believed to be less prevalent in Asian populations2. Symptoms of pernicious anemia range from fatigue to megaloblastic anemia and neurological abnormalities (peripheral numbness, paresthesia, ataxia) in more serious cases3.

Pernicious anemia is a complex disease with familial clustering, with a clear autoimmune basis and higher incidence of other autoimmune diseases, such as autoimmune thyroid conditions4, vitiligo5, type 1 diabetes6 in both patients with pernicious anemia and their relatives7. Although studies focusing on HLA serotypes have been conducted for pernicious anemia, the results have been conflicting7 and currently there is no clear consensus on the HLA alleles or other genetic risk factors predisposing to pernicious anemia.

We conducted a genome-wide association study (GWAS) of vitamin B12 deficiency due to lack of intrinsic factor to evaluate the contribution of genetic variation to the etiology of this disease in a combined dataset of 2,166 cases and 659,516 controls from three large population based biobanks - Estonian Biobank (EstBB)8, UK Biobank (UKBB)9, and FinnGen study. In the EstBB, individuals with pernicious anemia were identified using the ICD-10 code D51.0, resulting in 378 cases and 138,207 controls for analysis. Association testing was carried out with SAIGE 0.38 software10, adjusting for sex, year of birth and 10 PCs. Analysis in the EstBB was carried out under ethical approval 1.1-12/624 from the Estonian Committee on Bioethics and Human Research and data release N05 from the EstBB. GWAS summary statistics for the UKBB analysis including White British participants were downloaded from the UKBB PheWeb ([http://pheweb.sph.umich.edu/SAIGE-UKB/about](http://pheweb.sph.umich.edu/SAIGE-UKB/about)). Similarly, cases had been identified using the ICD-10 code D51.0 (754 cases, 390026 controls) and SAIGE had been used for association testing, adjusting for sex, birth year and the first 4 PCs. Summary statistics for FinnGen study were obtained from the publicly available R3 release PheWeb ([https://www.finngen.fi/en/access_results](https://www.finngen.fi/en/access_results)). In the FinnGen study, we used the endpoint ‘Vitamin B12 deficiency anemia’, which included all the subcodes in the ICD10 D51 diagnosis group (1,034 cases, 131,283 controls). Similarly to other cohorts, SAIGE had been used for association testing, adjusting for sex, age, 10 PCs and genotyping batch. For meta-analysis, we used fixed-effects meta-analysis implemented in GWAMA11. Additional details available in **Methods** section.

In our GWAS meta-analysis we identify five genome-wide significant (p<5 × 10−8) associations (Table 1, Figure 1, Supplementary Figure 1) on 1p13.2 (lead signal rs6679677, p=1.91 × 10−24), 2p16.1 (rs12616502, p=3.14 × 10−8), 6p21.32 (rs28414666, p=1.40 × 10−16), 10p15.1 (rs2476491, p=1.90 × 10−8), and 21q22.3 (rs74203920, p=2.33 × 10−9), with similar effect estimates (effect heterogeneity was measured using Cochran’s test, p-values ranging from 0.05 to 0.83) in all analyzed cohorts (Figure 1), except for the lead variant on chromosome 2, which is absent in the Finnish data.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/10/14/2020.10.13.20211912/T1)

Table 1. Summary of GWAS meta-analysis results for pernicious anemia

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.10.13.20211912/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/10/14/2020.10.13.20211912/F1)

Figure 1. Results of the pernicious anemia GWAS meta-analysis.
A) Manhattan and B) QQ plot; C) Forest plot of effect estimates for lead variants associated with pernicious anemia. The odds ratios (dots) and 95% confidence intervals (error bars) are shown for all included cohorts and meta-analysis. The size of the dot is proportional to the effective sample size, (calculated as 4/((1/N_cases)+(1/N_controls)). Pernicious anemia is defined as ICD10 code D51.0 in EstBB and UKBB and as D51 (vitamin B12 deficiency anemia) in FinnGen.

The lead SNP on chr1, rs6679677, is in high LD (r2=0.96) with a non-synonymous variant in exon 12 of the *PTPN22* gene (rs2476601, p=2.82 × 10−24). *PTPN22* is a known immune regulator gene and this particular variant has been associated with several autoimmune diseases, including rheumatoid arthritis, Crohn’s disease, systemic lupus erythematosus, vitiligo, autoimmune thyroid conditions, type 1 diabetes and others (Supplementary Table 1 and 3, Supplementary Figure 2).

On chromosome 2, the lead signal rs12616502 is in high LD with exonic non-synonymous variants in *CCDC104* (r2=0.94, rs1045920) and *PNPT1* (r2=0.94, rs7594497), and this region has been previously associated with both vitiligo, hypothyroidism and myelodysplastic syndrome (Supplementary Tables 1 and 3, Supplementary Figure 2). Colocalization analysis demonstrated that pernicious anemia has a common causal variant with the expression of *PTPN1* in thyroid tissue in GTEx v8 dataset (Figure 2, Supplementary Table 2). *PTPN1* is further supported as a candidate causal gene in this locus by data from mouse knockouts, as *Pnpt1tm1a(KOMP)Wtsi/Pnpt1+* mice exhibit increased mean corpuscular volume (MCV) together with increased mean corpuscular hemoglobin (MCH)12. Increased red blood cell MCV is a common feature in macrocytic anemias, both megaloblastic (caused by B12 deficiency and pernicious anemia) and nonmegaloblastic (caused by diseases such as myelodysplastic syndrome and hypothyroidism)13.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.10.13.20211912/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/10/14/2020.10.13.20211912/F2)

Figure 2. GWAS signature for pernicious anemia co-localizes with the eQTL signal for *PNPT1* in thyroid tissue in GTEx v8 dataset (posterior probability for shared causal variant PP4=0.87).
Left – comparison on -log10 p-values from pernicious anemia GWAS-meta-analysis (x-axis) and eQTL analysis(y-axis). Right – locus plots for the *PNPT1* locus from GWAS meta-analysis (top) and eQTL analysis (bottom). For all plots, LD is colored with respect to the GWAS lead signal rs12616502 (labeled).

The third signal on chromosome 6 is located in the HLA region, a common hub for autoimmune condition associations, downstream *HLA-DQB1*. To further clarify the association signal in the HLA region, we used data on a total of 249 imputed HLA alleles at four-digit level available for 183 cases and 67,505 controls in the EstBB. Three tested alleles passed the Bonferroni corrected threshold of association (0.05/249=2 × 10−4), including HLA-DRB1*15:01 (p=1.3 × 10−6), HLA-DQB1*06:02 (p=1.3 × 10−6) and HLA-DQA1*01:02 (p=2.8 × 10−5). The DRB1\*15:01-DQB1\*06:02-DQA1*01:02 combination forms the HLA-DR15 haplotype, which is a reported risk factor for multiple sclerosis14. HLA-DR15 belongs under the HLA-DR2 group, which has been associated with pernicious anemia in a 1981 study15.

On chromosome 10, the sentinel variant rs2476491 is intronic to *IL2RA. IL2RA* encodes the interleukin-2 receptor alpha chain, thus being involved in regulating regulatory T-cells and immune tolerance, as regulatory T cells suppress autoreactive T-cells. Accordingly, this locus has previously been associated with multiple sclerosis, juvenile idiopathic arthritis, vitiligo and hypothyroidism (Supplementary Table 1).

Finally, the association on chromosome 21, rs74203920, is a missense variant in the *AIRE* gene, a known autoimmune regulator. Interestingly, mutations in *AIRE* are a known cause of autoimmune polyendocrinopathy syndrome type 1 (APS-1), which is a rare autosomal recessive syndrome, that sometimes includes pernicious anemia among other components16.

We used the individual level data in EstBB to evaluate the association between pernicious anemia and other diseases (defined by ICD-10 codes). According to our analysis, individuals with pernicious anemia have more diagnoses of other anemias and vitamin deficiencies (Figure 3, Supplementary Table 4), thyroid problems (thyroiditis and hypothyroidism), gastrointestinal tract diagnoses (gastritis, malignant neoplasm of stomach, intestinal malabsorption, irritable bowel syndrome), but also of vitiligo, dermatitis, osteoporosis and spontaneous abortion. Evaluation of disease timeline shows that while there is considerable individual variation, on average, thyroid issues (ICD10 E03 other hypothyroidism and E06 thyroiditis, Figure 3), intestinal problems (K12 stomatitis, K29 gastritis and K58 irritable bowel syndrome), spontaneous abortion (O03), depression (F32) and vitiligo (L80) are diagnosed before pernicious anemia, while nutritional deficiencies are diagnosed after (D52 folate deficiency, E53 deficiency of other B vitamins, E55 vitamin D deficiency, E61 deficiency of other nutrient elements). Majority of these diagnoses reflect the etiology of pernicious anemia (gastritis) or symptoms (skin problems, syncope and collapse, depression, stomatitis), known comorbidities (vitiligo, thyroid issues4,5), or diseases where pernicious anemia is a known risk factor (such as osteoporosis17 and stomach cancer18). The association with spontaneous abortion is interesting, as although there is some evidence B12 deficiency and pernicious anemia could cause recurrent miscarriage19–21, the data is scarce and the link with spontaneous miscarriage has not been explored in depth.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/14/2020.10.13.20211912/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/10/14/2020.10.13.20211912/F3)

Figure 3. Association plot for ICD10 codes associated with pernicious anemia diagnosis in Estonian Biobank.
A) Each triangle in the plot represents one ICD10 main code and the direction of the triangle represents direction of effect – upward-pointing triangles show increased prevalence of diagnosis code in pernicious anemia cases. Red line – Bonferroni-corrected significance level (2.5 × 10−5). B) Time at first diagnosis of associated diagnosis code compared to time of pernicious anemia (D51.0) diagnosis. Plot shows mean difference in years, error bars correspond to standard deviation. Red dashed vertical line represents D51.0 diagnosis.

This GWAS meta-analysis focused on individuals with pernicious anemia identified from population-based biobanks, thus more detailed analysis on the effect these genetic risk factors have on disease severity or subphenotypes is unfortunately not possible. Although the used summary statistics for FinnGen cohort are for a broader vitamin B12 deficiency phenotype definition compared toother two cohorts, we see no significant differences in effect estimates for the reported lead variants. The analyses included participants of European ancestry, however, analyses in other populations are warranted, since the prevalence of pernicious anemia differs depending on racial background, being more common in people with African or European ancestry, especially from Scandinavia and UK2. The analysis of associated phenotypes could potentially provide clinically useful insights into pernicious anemia disease trajectories and offer information for patient management; however, at the moment it only includes participants of the EstBB, with limited follow-up time, therefore further studies are needed.

In summary, our analysis of 2,166 cases and 659,516 controls identifies robust risk loci for pernicious anemia in or near candidate genes with a known role in autoimmune conditions (*PTPN22, HLA, IL2RA, AIRE*). We further narrow the signal in the HLA region to the HLA-DR15 haplotype and propose *PNPT1* as the most likely causal gene in the 2p16.1 locus. The associations between the identified loci and other autoimmune conditions, such as type 1 diabetes, vitiligo, and autoimmune thyroid conditions help to clarify the link between pernicious anemia and its common comorbidities. Analysis of associated diagnoses and disease trajectories confirm the association between pernicious anemia and thyroid issues, vitiligo, gastritis, stomach cancer, osteoporosis and other diagnoses, but also between pernicious anemia and spontaneous abortion.

## Methods

### Cohorts

#### Estonian Biobank

The Estonian Biobank is a population-based biobank with over 200,000 participants. The 150K data freeze was used for the analyses described in this paper. All biobank participants have signed a broad informed consent form. Individuals with pernicious anemia were identified using the ICD-10 code D51.0, and all biobank participants who did not have this diagnosis were considered as controls. Information on ICD codes is obtained via regular linking with the national Health Insurance Fund and other relevant databases8.

All EstBB participants have been genotyped at the Core Genotyping Lab of the Institute of Genomics, University of Tartu, using Illumina GSAv1.0, GSAv2.0, and GSAv2.0\_EST arrays. Samples were genotyped and PLINK format files were created using Illumina GenomeStudio v2.0.4. Individuals were excluded from the analysis if their call-rate was < 95% or if sex defined based on heterozygosity of X chromosome did not match sex in phenotype data. Before imputation, variants were filtered by call-rate < 95%, HWE p-value < 1e-4 (autosomal variants only), and minor allele frequency < 1%. Variant positions were updated to b37 and all variants were changed to be from TOP strand using GSAMD-24v1-0_20011747_A1-b37.strand.RefAlt.zip files from [https://www.well.ox.ac.uk/~wrayner/strand/](https://www.well.ox.ac.uk/~wrayner/strand/) webpage. Prephasing was done using Eagle v2.3 software22 (number of conditioning haplotypes Eagle2 uses when phasing each sample was set to: --Kpbwt=20000) and imputation was done using Beagle v.28Sep18.79323 with effective population size ne=20,000. Population specific imputation reference of 2297 WGS samples was used.

Association analysis was carried out using SAIGE (v0.38)10 software implementing mixed logistic regression model with LOCO option, using sex, year of birth and 10 PCs as covariates in step I.

#### UK Biobank

The UK Biobank (UKBB) is a prospective cohort of 502,637 individuals aged 37-73 recruited in 2006-2010 from across the UK, who completed detailed questionnaires regarding socio-demographic and lifestyle characteristics and their medical history and had a clinical assessment. Additional information about medical conditions (both existing at baseline and occurring during follow-up) has been obtained through linking with hospital admission and mortality data. Full details of the study have been reported in Sudlow et al9. Publicly available GWAS summary statistics downloaded from the UKBB PheWeb [[http://pheweb.sph.umich.edu/SAIGE-UKB/about](http://pheweb.sph.umich.edu/SAIGE-UKB/about)] were used for the analysis. Briefly, the PheWeb includes GWAS summary statistics for ICD code-based traits extracted from electronic health records. Phenotypes have been classified into 1,403 broad PheWAS codes, including pernicious anemia (PheCode 281.11), defined using the ICD-10 code D51.0 and excluding other anemias under the PheCodes 280-285.99. Genetic analyses have been carried out using SAIGE10.

#### FinnGen

FinnGen is a public-private partnership project combining data from Finnish biobanks and electronic health records from different registries. After a one-year embargo, the FinnGenv summary stats are available for download. In this study, we used the results from the FinnGen release R3, which includes data from 135,638 individuals and more than 1,800 disease endpoints. FinnGen individuals have been genotyped with Illumina and Affymetrix arrays and imputed to the population-specific SISu v3 importation reference panel. Genetic association testing has been carried out with SAIGE10. The FinnGen disease endpoint “Vitamin B12 deficiency anemia” included all individuals with the ICD10 D51 diagnosis as cases. For more information on genotype data, disease endpoints and GWAS analyses, please see [https://finngen.gitbook.io/documentation/](https://finngen.gitbook.io/documentation/).

### GWAS meta-analysis

We extracted all genetic variants with a rs-number from the summary statistics of the three participating cohorts and conducted an inverse of variance weighted fixed-effects meta-analysis without genomic control using GWAMA11. A total of 30,907,385 variants were included in the meta-analysis. Genome-wide significance was set to p < 5 × 10−8.

### HLA allele imputation in the EstBB

Imputation of HLA alleles from SNP data was carried out at the Broad Institute using the SNP2HLA tool24.

### Colocalization

We conducted colocalization analyses to detect shared causal variants between pernicious anemia and gene expression using COLOC (v.3.2.1) R package25 and GWAS meta-analysis summary statistics. We set the prior probabilities to p1=1 × 10−4, p2=1 × 10−4, p12=5 × 10−6 as suggested by Wallace (2020)26 and used the COLOC version which takes regression coefficients and their variance into account.

In the analysis we compared our significant GWAS loci to all eQTL Catalogue ([https://www.ebi.ac.uk/eqtl/](https://www.ebi.ac.uk/eqtl/)) RNA-seq datasets (excluding Lepik et al. 201727 due to sample overlap) containing QTLs for gene expression, exon expression, transcript usage and txrevise event usage; eQTL Catalogue microarray datasets containing QTLs for gene expression (excluding Kasela et al 201728); and GTEx v8 datasets containing QTLs for gene expression (see Methods: [https://www.ebi.ac.uk/eqtl/Methods/](https://www.ebi.ac.uk/eqtl/Methods/)).

We lifted the GWAS summary statistics over to hg38 build to match the eQTL Catalogue. For each genome-wide significant (p<5 × 10−8) GWAS locus we extracted the 1Mbp radius of its top hit from QTL datasets and ran the colocalization analysis for those eQTL Catalogue traits that had at least one cis-QTL within this region with p< 1 × 10−6. We considered two signals to colocalize if the posterior probability for a shared causal variant was 0.8 or higher. All results with a PP4 > 0.8 can be found in Supplementary Table 2.

Results were visualized with the *LocusCompareR* library30.

### Look-up of phenome-wide associations in GWAS catalog and with PhenoScanner v2

FUMA v1.3.6a29 was used to compare the genome-wide significant lead signals and markers in high LD with these markers against the results in the GWAS catalogue. The results of this look-up are presented in Supplementary Table 1.

PhenoScanner v231,32 was used for look-up of phenotype associations for the GWAS lead variants in previous GWAS studies. PhenoScanner query was done using the rsid-s of GWAS lead variants and the *phenoscanner* R package ([https://github.com/phenoscanner/phenoscanner](https://github.com/phenoscanner/phenoscanner)). Query results were filtered to keep one association per variant per trait, keeping studies from newer or larger studies. Descriptions of Experimental Factor Ontology (EFO) terms and classification of EFO broad categories were obtained from the GWAS Catalog. Missing categories were added by manually searching the EMBL-EBI EFO webpage ([www.ebi.ac.uk/efo/](https://www.ebi.ac.uk/efo/)). For visualization of PhenoScanner results, parent categories with fewer results were grouped into larger categories and a heatmap was created using the *pheatmap* library in R 3.6.1. and a modified script from ([https://github.com/LappalainenLab/spiromics-covid19-eqtl/blob/master/eqtl/summary\_phenoscanner_lookup.Rmd](https://github.com/LappalainenLab/spiromics-covid19-eqtl/blob/master/eqtl/summary_phenoscanner_lookup.Rmd)). The results of this look-up are presented in Supplementary Table 3.

### Mouse phenotypes

We used the Mouse Genome Database12 ([http://www.informatics.jax.org](http://www.informatics.jax.org)) to evaluate the *PNPT1* effect on phenotype in mouse models.

### Analysis of associated phenotypes in EstBB

Using the individual level data in the EstBB, we conducted an analysis to find ICD10 diagnosis codes associated with the D51.0 diagnosis. We tested the association between pernicious anemia status (defined as ICD10 D51.0) and other ICD10 codes using logistic regression and adjusting for sex, age and 10 PCs. Bonferroni correction was applied to select statistically significant associations (Number of tested ICD main codes – 1,944, corrected p-value threshold – 2.5 × 10− 5). Results were visualised using the *PheWas* library ([https://github.com/PheWAS/PheWAS](https://github.com/PheWAS/PheWAS)). All analyses were carried out in R 3.6.1. Diagnosis timeline of pernicious anemia and significantly associated diagnoses was evaluated by calculating the time interval (in years) between the first date of pernicious anemia diagnosis (D51.0) and first date of associated disease for each individual. Results were visualized by plotting the means and standard deviations of calculated intervals using R 3.6.1. The results of this analysis are presented in Supplementary Table 4.

## Supporting information

Supplementary tables [[supplements/211912_file02.xlsx]](pending:yes)

Supplementary Figure 1 [[supplements/211912_file03.pdf]](pending:yes)

Supplementary Figure 2 [[supplements/211912_file04.pdf]](pending:yes)

## Data Availability

Used UKBB and FinnGen summary statistics can be browsed and downloaded from UKBB PheWeb and FinnGen PheWeb, respectively. Full meta-analysis summary statistics will be made available upon publication. The analyses in this manuscript also included data from the Mouse Genome Database and the GTEx Portal.

[http://pheweb.sph.umich.edu/SAIGE-UKB/about](http://pheweb.sph.umich.edu/SAIGE-UKB/about) 

[http://r3.finngen.fi](http://r3.finngen.fi) 

[http://www.informatics.jax.org](http://www.informatics.jax.org) 

[https://www.ebi.ac.uk/eqtl/](https://www.ebi.ac.uk/eqtl/) 

## Data availability

Used UKBB and FinnGen summary statistics can be browsed and downloaded from UKBB PheWeb and FinnGen PheWeb, respectively. Full meta-analysis summary statistics will be made available upon publication.

## Web resources

UKBB PheWeb: [http://pheweb.sph.umich.edu/SAIGE-UKB/about](http://pheweb.sph.umich.edu/SAIGE-UKB/about); FinnGen Freeze 3 PheWeb: [http://r3.finngen.fi](http://r3.finngen.fi); Mouse Genome Database: [http://www.informatics.jax.org](http://www.informatics.jax.org); GTEx Portal: [https://gtexportal.org/home/](https://gtexportal.org/home/); eQTL Catalogue: [https://www.ebi.ac.uk/eqtl/](https://www.ebi.ac.uk/eqtl/)

## Declaration of Interests

The authors declare no competing interests.

## Acknowledgements

TL is supported by the Estonian Research Council grant MOBTP155. ML and RM are supported by the Estonian Research Council grant PRG687. Computations were performed in the High Performance Computing Center, University of Tartu. We want to acknowledge the participants and investigators of the FinnGen and UKBB studies. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from the GTEx Portal on 10/01/20.

*   Received October 13, 2020.
*   Revision received October 13, 2020.
*   Accepted October 14, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Lenti, M.V., Rugge, M., Lahner, E., Miceli, E., Toh, B.H., Genta, R.M., De Block, C., Hershko, C., and Di Sabatino, A. (2020). Autoimmune gastritis. Nat. Rev. Dis. Prim. 6, 1–19.
    
    

2.  2.Stabler, S.P., and Allen, R.H. (2004). Vitamin B12 deficiency as a worldwide problem. Annu. Rev. Nutr. 24, 299–326.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.nutr.24.012003.132440&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15189123&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222766200014&link_type=ISI) 

3.  3.Bizzaro, N., and Antico, A. (2014). Diagnosis and classification of pernicious anemia. Autoimmun. Rev. 13, 565–568.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.autrev.2014.01.042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24424200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

4.  4.Ness-Abramof, R., Nabriski, D.A., Braverman, L.E., Shilo, L., Weiss, E., Reshef, T., Shapiro, M.S., and Shenkman, L. (2006). Prevalence and evaluation of B12 deficiency in patients with autoimmune thyroid disease. Am. J. Med. Sci. 332, 119–122.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16969140&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

5.  5.Dahir, A.M., and Thomsen, S.F. (2018). Comorbidities in vitiligo: comprehensive review. Int. J. Dermatol. 57, 1157–1164.
    
    

6.  6.Kahaly, G.J., and Hansen, M.P. (2016). Type 1 diabetes associated autoimmunity. Autoimmun. Rev. 15, 644–648.
    
    

7.  7.Banka, S., Ryan, K., Thomson, W., and Newman, W.G. (2011). Pernicious anemia - Genetic insights. Autoimmun. Rev. 10, 455–459.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.autrev.2011.01.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21296191&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

8.  8.Leitsalu, L., Haller, T., Esko, T., Tammesoo, M.-L., Alavere, H., Snieder, H., Perola, M., Ng, P.C., Mägi, R., Milani, L., et al. (2015). Cohort Profile: Estonian Biobank of the Estonian Genome Center, University of Tartu. Int. J. Epidemiol. 44, 1137–1147.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyt268&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24518929&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

9.  9.Sudlow, C., Gallacher, J., Allen, N., Beral, V., Burton, P., Danesh, J., Downey, P., Elliott, P., Green, J., Landray, M., et al. (2015). UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age. PLoS Med. 12, e1001779.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001779&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25826379&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

10. 10.Zhou, W., Nielsen, J.B., Fritsche, L.G., Dey, R., Gabrielsen, M.E., Wolford, B.N., LeFaive, J., VandeHaar, P., Gagliano, S.A., Gifford, A., et al. (2018). Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies. Nat. Genet. 50, 1335–1341.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0184-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30104761&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

11. 11.Mägi, R., and Morris, A.P. (2010). GWAMA: software for genome-wide association meta-analysis. BMC Bioinformatics 11, 288.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-11-288&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20509871&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

12. 12.Bult, C.J., Blake, J.A., Smith, C.L., Kadin, J.A., Richardson, J.E., Anagnostopoulos, A., Asabor, R., Baldarelli, R.M., Beal, J.S., Bello, S.M., et al. (2019). Mouse Genome Database (MGD) 2019. Nucleic Acids Res. 47, D801–D806.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky1056&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30407599&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

13. 13.Nagao, T., and Hirokawa, M. (2017). Diagnosis and treatment of macrocytic anemias in adults. J. Gen. Fam. Med. 18, 200–204.
    
    

14. 14.Schmidt, H., Williamson, D., and Ashley-Koch, A. (2007). HLA-DR15 haplotype and multiple sclerosis: A HuGE review. Am. J. Epidemiol. 165, 1097–1109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwk118&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17329717&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246120000001&link_type=ISI) 

15. 15.Ungar, B., Mathews, J.D., Tait, B.D., and Cowling, D.C. (1981). HLA-DR patterns in pernicious anaemia. Br. Med. J. (Clin. Res. Ed). 282, 768–770.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIyODIvNjI2Ni83NjgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMC8xNC8yMDIwLjEwLjEzLjIwMjExOTEyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

16. 16.Bruserud, Ø., Oftedal, B.E., Wolff, A.B., and Husebye, E.S. (2016). AIRE-mutations and autoimmune disease. Curr. Opin. Immunol. 43, 8–15.
    
    

17. 17.Goerss, J.B., Kim, C.H., Atkinson, E.J., Eastell, R., O’Fallon, W.M., and Melton, L.J. (1992). Risk of fractures in patients with pernicious anemia. J. Bone Miner. Res. 7, 573–579.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1615763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HV08200013&link_type=ISI) 

18. 18.Ye, W., and Nyrén, O. (2003). Risk of cancers of the oesophagus and stomach by histology or subsite patients hospitalised for pernicious anaemia. Gut 52, 938–941.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjUyLzcvOTM4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMTQvMjAyMC4xMC4xMy4yMDIxMTkxMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

19. 19.Bennett, M. (2001). Vitamin B 12 deficiency, infertility and recurrent fetal loss. J. Reprod. Med. Obstet. Gynecol. 46, 209–212.
    
    

20. 20.Reznikoff-Etiévant, M.F., Zittoun, J., Vaylet, C., Pernet, P., and Milliez, J. (2002). Low Vitamin B 12 level as a risk factor for very early recurrent abortion. Eur. J. Obstet. Gynecol. Reprod. Biol. 104, 156–159.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0301-2115(02)00100-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12206930&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000178413300013&link_type=ISI) 

21. 21.Hübner, U., Alwan, A., Jouma, M., Tabbaa, M., Schorr, H., and Herrmann, W. (2008). Low serum vitamin B12 is associated with recurrent pregnancy loss in Syrian women. Clin. Chem. Lab. Med. 46, 1265–1269.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18636794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

22. 22.Loh, P.R., Danecek, P., Palamara, P.F., Fuchsberger, C., Reshef, Y.A., Finucane, H.K., Schoenherr, S., Forer, L., McCarthy, S., Abecasis, G.R., et al. (2016). Reference-based phasing using the Haplotype Reference Consortium panel. Nat. Genet. 48, 1443–1448.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3679&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27694958&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

23. 23.Browning, S.R., and Browning, B.L. (2007). Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 81, 1084–1097.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/521987&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17924348&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250480900018&link_type=ISI) 

24. 24.Jia, X., Han, B., Onengut-Gumuscu, S., Chen, W.M., Concannon, P.J., Rich, S.S., Raychaudhuri, S., and de Bakker, P.I.W. (2013). Imputing Amino Acid Polymorphisms in Human Leukocyte Antigens. PLoS One 8,.
    
    

25. 25.Giambartolomei, C., Vukcevic, D., Schadt, E.E., Franke, L., Hingorani, A.D., Wallace, C., and Plagnol, V. (2014). Bayesian Test for Colocalisation between Pairs of Genetic Association Studies Using Summary Statistics. PLoS Genet. 10, e1004383.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

26. 26.Wallace, C. (2020). Eliciting priors and relaxing the single causal variant assumption in colocalisation analyses. PLoS Genet. 16,.
    
    

27. 27.Lepik, K., Annilo, T., Kukuškina, V., Kisand, K., Kutalik, Z., Peterson, P., and Peterson, H. (2017). C-reactive protein upregulates the whole blood expression of CD59 - an integrative analysis. PLoS Comput. Biol. 13,.
    
    

28. 28.Kasela, S., Kisand, K., Tserel, L., Kaleviste, E., Remm, A., Fischer, K., Esko, T., Westra, H.J., Fairfax, B.P., Makino, S., et al. (2017). Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4 + versus CD8 + T cells. PLoS Genet. 13,.
    
    

29. 29.Watanabe, K., Taskesen, E., van Bochoven, A., and Posthuma, D. (2017). Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01261-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

30. 30.Liu, B., Gloudemans, M.J., Rao, A.S., Ingelsson, E., and Montgomery, S.B. (2019). Abundant associations with gene expression complicate GWAS follow-up. Nat. Genet. 51, 768–769.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0404-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31043754&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

31. 31.Staley, J.R., Blackshaw, J., Kamat, M.A., Ellis, S., Surendran, P., Sun, B.B., Paul, D.S., Freitag, D., Burgess, S., Danesh, J., et al. (2016). PhenoScanner: A database of human genotype-phenotype associations. Bioinformatics 32, 3207–3209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw373&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27318201&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom) 

32. 32.Kamat, M.A., Blackshaw, J.A., Young, R., Surendran, P., Burgess, S., Danesh, J., Butterworth, A.S., Staley, J.R., and Kelso, J. (2019). PhenoScanner V2: An expanded tool for searching human genotype-phenotype associations. Bioinformatics 35, 4851–4853.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F14%2F2020.10.13.20211912.atom)