Analysis of MRI-derived spleen iron in the UK Biobank identifies genetic variation linked to iron homeostasis and erythrocyte morphology
========================================================================================================================================

* Elena P. Sorokin
* Nicolas Basty
* Brandon Whitcher
* Yi Liu
* Jimmy D. Bell
* Robert L. Cohen
* Madeleine Cule
* E. Louise Thomas

## Abstract

The spleen plays a key role in iron homeostasis. It is the largest filter of the blood and performs iron reuptake from old or damaged erythrocytes. Despite this role, spleen iron concentration has not been measured in a large, population-based cohort. In this study, we quantify spleen iron in 41,764 participants of the UK Biobank using magnetic resonance imaging, and provide the first reference range for spleen iron in an unselected population. Through genome-wide association study, we identify associations between spleen iron and regulatory variation at two hereditary spherocytosis genes, *ANK1* and *SPTA1*. Spherocytosis-causing coding mutations in these genes are associated with lower reticulocyte volume and increased reticulocyte percentage, while these novel common alleles are associated with increased expression of *ANK1* and *SPTA1* in blood and with larger reticulocyte volume and reduced reticulocyte percentage. As genetic modifiers, these common alleles may explain mild spherocytosis phenotypes that have been observed clinically. Our genetic study also identifies a signal which co-localizes with a splicing quantitative trait locus for *MS4A7*, and we show this gene is abundantly expressed in the spleen and in macrophages. The combination of deep learning and efficient image processing enables non-invasive measurement of spleen iron and, in turn, characterization of genetic factors related to iron recycling and erythrocyte morphology.

Keywords
*   Spleen iron
*   Magnetic Resonance Imaging
*   Imaging-derived phenotype
*   quantitative traits
*   human iron economy
*   iron homeostasis
*   hereditary spherocytosis
*   hemolytic anemia
*   genome-wide association study
*   rare variant association study

## Introduction

In normal human physiology, iron is recycled much faster than new dietary iron is absorbed, and this iron economy is regulated by the spleen.1–3 The spleen plays a critical role in removing senescent erythrocytes from the blood, and does so via a population of splenic macrophages as well as by the action of the protein ferroportin, which transports iron back to the plasma.4 Given the large iron flux through the normal spleen due to erythrocyte recycling, measurement of spleen iron has the potential to reflect activity of iron salvage pathways.

Spleen iron also may reflect erythrocyte biology and dysfunction. Hereditary spherocytosis (HS) is a relatively common erythrocyte membrane disorder, occurring at a prevalence of 1:1000-2500 in European populations.5 Patients have defects in membrane and cytoskeletal genes that contribute to erythrocyte membrane integrity and deformability, including *SPTA1 (*encoding the filamentous protein alpha-spectrin) and *ANK1 (*encoding a protein that tethers spectrin filaments to erythrocyte membranes*)*.*6,7* The protein products of these genes interact in the formation of the mature erythrocyte cytoskeleton, a process accompanying cellular remodeling during reticulocytosis.8 In HS, the resulting erythrocytes are spherically shaped and lose deformability, and are ultimately trapped and ingested by macrophages within the red pulp of the spleen, resulting in an enlarged spleen and anemia despite brisk reticulocytosis. When measured directly with radiolabeled cells, erythrocyte turnover is dramatically accelerated in patients with severe HS.9 While laboratory measures such as sphered cell and reticulocyte volume can diagnose HS10, less is understood about the phenotypic and genetic heterogeneity of the disease.11–13

Although spleen iron has been investigated in specific disease groups14–20, few studies were performed in unascertained cohorts21 partly due to limitations in detecting and quantifying low levels of tissue iron in non-overloaded populations.22,23 The UK Biobank (UKBB) is a prospective study of half a million adults in the UK,24 with genetic and phenotypic data including magnetic resonance imaging (MRI).25 In this study, we applied computer vision techniques to quantify spleen iron non-invasively and at scale, by repurposing the dedicated liver MRI acquisition. Spleen iron is only moderately correlated with other measures of iron stores in the body, making it a novel trait. Through genome-wide association study (GWAS), we characterized associations between spleen iron and genes involved in the human iron economy including common regulatory variation in the ferroportin gene, *SLC40A1*. Our GWAS identified a signal which co-localizes with a splicing quantitative trait locus for *MS4A7*, which we found to be abundantly expressed in the spleen and in macrophages. Regulatory variation in the genes encoding alpha-spectrin (*SPTA1*) and ankyrin (*ANK1*) were also linked to spleen iron, increased mRNA expression of these genes, and furthermore exhibited effects on reticulocyte and erythrocyte parameters opposite to the effects observed in HS patients. Non-invasive imaging coupled with genetic analysis enabled us to develop a novel trait and illuminated new aspects of iron homeostasis and erythrocyte biology.

## Results

### Characterization of spleen iron in a large, population-based cohort

We quantified spleen iron concentration (spleen iron hereafter) in 41,764 UK Biobank (UKBB) participants with both the 3D neck-to-knee and the quantitative liver single-slice MRI sequences available. We opportunistically measured spleen iron concentration by segmenting the spleen from the neck-to-knee image26 (**Figure 1A**) and subsequently extracting a 2D mask where the spleen volume intersects with the quantitative 2D liver slice (**Figure 1B**).45 The average spleen iron was 0.92±0.32 mg/g, significantly lower than liver iron of 1.24±0.29 mg/g (**Table 1, Figure 1C**; paired t-test p < 2.2e-16). While there is no accepted normal range of spleen iron, 1.98 mg/g has been suggested as an upper cut-off, with 2.74 mg/g reported to be pathological.27 Using 1.98 mg/g as the threshold, 1.04% (n=435) of this cohort had elevated spleen iron, while 0.32% (n=137) had spleen iron above the 2.74 mg/g threshold. 95% of the population fell into the range of 0.53 to 1.67 mg/g, and we propose this as a possible reference range in an unselected population.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F1)

Figure 1: High-throughput spleen iron quantification.
**A)** Example neck-to-knee abdominal MRI acquisition, with the liver slice acquistion shown as a red plane, and example deep learning-derived spleen segmentation shown in green. **B)** Axial view of iron concentration (mg/g) obtained through deep learning and image processing, with liver shown at left and the spleen at right. **C)** Distribution of spleen iron (mg/g), as compared to liver iron, in the UK Biobank (n=41,764). **D)** Spleen iron increases with age for both males and females. Standardized residuals are shown, adjusted for study center, scan date and scan time.

View this table:
[Table 1:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/T1)

Table 1: Summary of the UKBB cohort at recruitment, image acquisition, and quantification of liver iron and spleen iron.
Spleen iron measurements (mg/g) are provided as mean (standard deviation) and range. For other quantitative values, the mean (standard deviation) is given. For binary values, mean [95% confidence interval] is given.

Spleen iron differed by age and sex. Men had higher spleen iron than women (men: 0.96±0.34 mg/g, women: 0.87±0.29 mg/g) (p=6.3e-219; **Supplementary Table 2**). Increasing spleen iron was associated with age (0.0044 mg/g/year or 0.012 SD/year) (**Figure 1D**). In women, menopause was associated with 0.12 mg/g higher spleen iron [95% CI 0.08-0.16].

In a phenome-wide association study with over 3,200 quantitative traits and disease outcomes, spleen iron was correlated with erythrocyte parameters: reticulocyte percentage (beta=0.091; p=3.7e-66), reticulocyte count (beta=0.087; p=6.6e-64), and high light scatter reticulocyte count (beta=0.089; p=5.8e-64). Spleen iron was also associated with lifestyle factors, including consumption of lamb (beta=0.150; p=2.1e-44) and beef (beta=0.143; p=1.3e-43), and negatively associated with alcohol consumption (beta=-0.087; p=1.7e-18). Spleen iron correlated with brain iron content, specifically with T2* (inversely proportional to iron) in the caudate (beta=-0.063; p=4.3e-16) and putamen (beta=-0.061, p=4.1e-15).28 Spleen iron was associated with myeloid leukemia (beta=0.386; p=9.6e-10), chronic dermatitis (beta=0.328; p=4.1e-07), hypokalemia (beta =0.286, p=5.8e-6), and glaucoma (beta=0.186, p=1.0e-6) (**Figure 2, Supplementary Tables 3-4**). Spleen iron was negatively correlated with iron-deficiency anemia, but this did not achieve Bonferroni significance (beta=-0.134, p=0.002), though this diagnosis may not have been fully captured by medical billing codes (n=580 cases, n=35,316 controls).

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F2)

Figure 2: Phenome-wide association between spleen iron and complex traits in the UK Biobank.
The horizontal line shows a significant association after Bonferroni correction. The top three associations in each category are annotated. ML: Myeloid leukaemia. HLS reticulocyte: High light scatter reticulocyte.

### Genome-wide association study of spleen iron identifies DNA polymorphisms linked to global iron homeostasis

In a common variant genome-wide association study (GWAS) of spleen iron, seven loci reached genome-wide significance (p < 5e-8; **Figure 3, Table 2**). Conditional analysis yielded no secondary signals. We estimated the narrow-sense heritability of spleen iron to be 16.7% (s.e. 1.64%). Spleen iron was moderately genetically correlated with ferritin (rg=0.56) and mean corpuscular hemoglobin concentration (rg=0.42), but the genetic correlations with other iron measures, including liver iron, were not significant (**Supplementary Tables 5-6)**.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F3)

Figure 3. Genome-wide association study of spleen iron.
Genome-wide significant signals are annotated with the variant at each locus. Dotted horizontal line marks genome-wide significance (p=5e-8).

View this table:
[Table 2:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/T2)

Table 2: Fine-mapped lead SNPs from genome-wide association study of spleen iron.
Genome-wide significant associations (p<5e-8) are shown by locus and lead SNP after fine-mapping. A seventh association at the major histocompatibility locus (MHC) could not be fine-mapped.

We observed a signal on chromosome 2 at *SLC40A1*, which encodes the iron transporter ferroportin (lead variant: rs13008848[G], beta=-0.057, p=2.9e-10, **Figure 4A**). To test whether the spleen signal was shared across other body iron traits, we re-analyzed previous genetic studies of serum ferritin, serum iron, and liver iron26,29, and observed evidence of co-localization at *SLC40A1*(posterior probability≥0.99) (**Supplementary Table 7**). Since the lead variant lies upstream of the *SLC40A1* open reading frame, we tested for a shared effect between spleen iron and *SLC40A1* expression. Across 56 human tissues of the GTEx Consortium,30 we observed evidence for regional co-localization with a quantitative trait locus for *SLC40A1* mRNA in many tissues, including whole blood (posterior probability≥0.99) (**Figure 4D**; **Supplementary Table 8**). As expected, spleen iron was associated with *SLC40A1*, and furthermore, we showed that this locus likely influences tissue iron levels through *SLC40A1* mRNA abundance.

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F4)

Figure 4: Common variants in *SLC40A1, SPTA1*, and *ANK1* loci are associated with spleen iron, and colocalize with cis-regulatory variation.
**(A)** Fine-mapped locus near *SLC40A1* is associated with spleen iron (lead SNP: rs13008848). **(B)** Fine-mapped genome-wide significant locus near *SPTA1* is associated with spleen iron (lead SNP rs2479868). **(C)** Fine-mapped locus near *ANK1* is associated with spleen iron (lead SNP: 4373010). **(D)** Co-localization of cis-regulatory variation in the *SLC40A1* locus and spleen iron. Co-localization at this locus was observed in multiple tissues (posterior probability≥0.99 for each). Blood cis-eQTL data is shown. **(E)** Cis-regulatory variation in the *SPTA1* locus co-localizes with spleen iron. A signal was observed in multiple tissues but did not meet significance threshold in whole blood (posterior probability=0.72). **(F)** Cis regulatory variation at the *ANK1* locus co-localizes with spleen iron. Co-localization was observed in multiple tissues (posterior probability≥0.99). Grey dashed line indicates genome-wide significance (p=5e-8). For all eQTL signals, a threshold of FDR<5% was used in GTEx38. Linkage disequilibrium was calculated using 1000 Genomes Phase 3. Gene models are shown at bottom in GRCh38 coordinates.

The GWAS of spleen iron identified other loci relevant to iron homeostasis, including a locus on chromosome 9 (**Supplementary Figure 5B**). The lead SNP rs41276777[A] (beta=0.17, p=2.1e-9) occurred in the 5’ untranslated region of *PRPF4* and *CDC26* (**Supplementary Figure 5A**), and was associated with a regulatory locus affecting expression of both *PRPF4* and *CDC26* mRNA in whole blood and the spleen (**Supplementary Table 8**), suggesting bi-directional regulation of gene expression. This region colocalized with signals for serum ferritin, serum iron, and other blood iron traits (**Supplementary Table 7**), suggesting roles not only in spleen iron but body iron more broadly.

### Elevated spleen iron colocalizes with a splicing quantitative trait locus for *MS4A7* and monocyte traits

Our GWAS with spleen iron found a significant association at the *MS4A7* locus (lead variant: rs950802[A], beta=0.086, p=7.7e-26; **Table 2, Figure 5A**). rs950802[A] causes a synonymous mutation at Leu57 in the third exon of *MS4A7* and is also a variant in the first intron of *MS4A14*. This signal colocalized with a signal of serum ferritin29 (posterior probability≥0.99. Using splicing quantitative trait data30, we identified an alternative splicing event in *MS4A7*, which colocalized with spleen iron (**Figure 5B**). Based on an analysis of the open reading frame of *MS4A7*, this alternative splicing event was predicted to increase skipping of the second exon and therefore interrupt a conserved CD20-like transmembrane domain in *MS4A7* (**Figure 5C)**. We thus identified a novel association between *MS4A7* and spleen iron, and found a plausible molecular mechanism by which this variant disrupts *MS4A7* function.

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F5)

Figure 5: Characterization of the *MS4A7* locus.
**(A)** The *MS4A7* locus on chromosome 11 is associated with an increase in spleen iron. Purple triangle displays the lead SNP rs950802[G] **(B)** This locus co-localizes with splicing quantitative trait locus of *MS4A7*, whole blood (posterior probability≥0.99). **(C)** This splicing quantitative trait locus is associated with exon skipping events in the *MS4A7* locus. Corresponding betas and p-value are shown. The *MS4A7* gene model is displayed with the conserved CD20-like domain shown in green. **(D)** Genetic associations with a panel of blood cell traits show associations with monocyte count, monocyte percentage and mean platelet volume, and negative associations with platelet count and platelet crit. **(E)** Expression of *MS4A7* was enriched in monocytes, including non-classical, intermediate, and classical forms. **(F)** Expression of *MS4A7* in human tissues displayed enrichment for lymphoid tissue, notably spleen and lymph nodes.

To explore the functional consequences of regulatory variation in *MS4A7*, we examined its association with hematological parameters in the UKBB, and observed a significant association with monocyte count and percentage, as well as platelet count and crit (**Figure 5D**). We replicated these findings in the Blood Cell Consortium31, confirming a positive association with monocyte count (beta=0.011, p=1.5e-8) and a negative association with platelet count (beta=-0.012, p=2.7e-10) in a European cohort (**Supplementary Table 10**). Further, we found enrichment of *MS4A7* mRNA in monocytes compared to sixteen other hematopoietic cell types32 (**Figure 5E**), and also higher expression of *MS4A7* mRNA in the spleen relative to other tissues (**Figure 5F**). A novel common allele at *MS4A7* was thus associated with elevated spleen iron and monocyte traits. *MS4A7* mRNA was enriched in the spleen and in monocytes, suggesting a role for this gene in iron recycling in the spleen.

### GWAS of spleen iron identifies novel common alleles in *SPTA1* and *ANK1* linked to increased gene expression and erythrocyte function

In addition to signals linked to iron homeostasis, the GWAS of spleen iron revealed associations in *SPTA1* and *ANK1*, encoding structural components of erythrocytes (**Table 2, Supplementary Figure 6**). The lead variant rs2479868[T] (beta=-0.083, p=2.5e-22) on chromosome 1 was located in the 3’ untranslated region of *SPTA1*, and the lead SNP rs4737010[A] (beta=-0.077, p = 1.0e-17) on chromosome 8 was the first intron of *ANK1* (**Figure 4B-C**). We tested for a shared signal with mRNA expression levels of *SPTA1* and *ANK1*, and observed co-localization with cis-regulatory variation for *ANK1* in multiple tissues (**Figure 4E-F**).

Variation at each locus was associated with increased mRNA expression and mean reticulocyte volume; and decreased mean corpuscular hemoglobin concentration, reticulocyte percentage, and spleen volume (**Figure 6**). The magnitude and directionality for the changes in red cell parameters associated with the lead SNPs at *SPTA1* and *ANK1* were similar, but a regression model testing for interaction between these two loci found that the effects are independent (beta=0.016, p=0.28). These signals did not co-localize with iron traits in other tissues such as serum or liver, suggesting specific effects to the spleen (**Supplementary Table 7**).

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/19/2021.11.16.21266431/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2022/01/19/2021.11.16.21266431/F6)

Figure 6: Associations of molecular, cellular, organ, and tissue traits with lead SNPs at *SPTA1* and *ANK1* loci indicate protective effects on red blood cell (erythrocyte) function.
Both loci were associated with increased mRNA expression, beneficial effects on red blood cell parameters, decreases in spleen iron and spleen volume, and no effects on other measures of body iron or anemia diagnosis. Standardized effect sizes are shown. Error bars represent 95% confidence intervals.

Using data from the Blood Cell Consortium31, we replicated the associations and directions of effect between *SPTA1* and *ANK1* lead variants and mean corpuscular hemoglobin concentration and mean corpuscular volume. In a cohort of European ancestry, we replicated all findings (p<5e-10). Additionally, we observed associations for both loci and MCHC in East Asian and Hispanic/Latino ancestry cohorts (p<0.01) (**Supplementary Table 10**).

### Common alleles in *SPTA1* and *ANK1* show effects on erythrocytes opposite to effects observed with rare deleterious alleles and in hereditary spherocytosis

Given the known associations of both *SPTA1* and *ANK1* gene defects with hereditary spherocytosis, we analyzed rare variation (minor allele frequency <0.001) predicted to cause loss-of-function. Starting with 167,246 exomes of European ancestry, we conducted rare variant association studies for reticulocyte percentage and volume and identified one significant gene, *SPTA1*, in both studies (pSKAT-O< 1e-24; **Supplementary Figure 7**). Performing a scan of the same loss-of-function rare variation in *SPTA1* across a hematology panel, we recapitulated clinical hallmarks of spherocytosis: *SPTA1* loss-of-function was significantly associated with increased reticulocyte percentage and increased MCHC, increased bilirubin, decreased mean spherical cell volume and decreased MRV (**Supplementary Figure 7)**. We did not observe a significant association with spleen iron in the imaging subcohort for either putative loss-of-function or deleterious missense variation, perhaps due to reduced statistical power of the imaging subcohort (n=18,420) compared to the exome cohort (n=167,246) (**Supplementary Figure 8**).

Finally, we genetically identified HS in the UKBB exome cohort using clinical assertions of pathogenicity and predicted high confidence loss-of-function in one of six HS genes, and estimated prevalence to be 1:389 [95% CI 1:427 - 1:354]. We asked whether the common alleles found via GWAS could modify the effects of rare deleterious alleles for hematology parameters relevant for HS, including mean reticulocyte volume and percentage. We estimated that carrying either of the *SPTA1* or *ANK1* lead GWAS SNPs skewed erythrocyte parameters toward beneficial effects in deleterious allele carriers and non-carriers, suggesting that the common GWAS alleles could modify the effects of the rare deleterious alleles (**Supplementary Figure 10**).

## Discussion

Most studies of spleen iron have not included healthy volunteers, or methods that reliably assess organ iron15,33–35 or only provided qualitative histological grading35,36 (**Supplementary Table 11**).14,17–19 In this study, we combined deep learning algorithms and efficient image processing to quantify spleen iron in 41,764 participants of the UKBB. This enabled us to estimate the first reference range for spleen iron in a large, unselected population.

Spleen iron was higher in men, and our analysis of pre-and post-menopausal women is consistent with prior reports of iron stores in other tissues.37 Spleen iron increased with age, is associated with red meat intake and inversely associated with alcohol consumption, extending prior observations in the liver.38-39,40Further, spleen iron was only modestly associated with measures of iron in other tissues. Unlike other iron-rich organs such as the liver, we discovered that spleen iron was associated with indicators of reticulocytosis and reticulocyte turnover.

Our genome-wide association study of spleen iron identifies regulatory loci in *SPTA1* and *ANK1*, which when combined with other evidence, suggests a model of low splenic turnover due to relatively large, long-lived erythrocytes. These alleles were associated with decreased spleen iron and increased mRNA expression of their respective *cis* gene, as well as larger reticulocyte volume, and reduced measures of reticulocytosis (**Figure 6**). Since the spleen is the major route of erythrocyte clearance and iron salvage, lower levels of reticulocytosis would be expected to result in lower spleen iron at steady state. Neither variant was associated with anemia, suggesting that they are not pathogenic. This reduced-turnover model also predicts that iron levels in other organs not involved in erythrocyte clearance (such as the pancreas) would be unaffected by these two loci, consistent with the data.

As genetic modifiers, these alleles in *SPTA1* and *ANK1*, which segregate frequently across global populations, may explain the variable penetrance and expressivity observed in hereditary spherocytosis (HS). In a recent study of patients with identified HS mutations, 64% involved *SPTA1* or *ANK1*, and the investigators observed multiple HS families with broad phenotypic variability, including a compelling example of dizygotic twins sharing the same pathogenic *ANK1* mutation presenting with mild disease in one case and severe disease requiring splenectomy in the other11. The variation could not be explained, leading to speculation that yet-unknown genetic factors may be contributing. We were able to identify the hallmarks of HS through rare loss-of-function variation in *SPTA1*, even within a cohort not ascertained for hemolytic anemias. Our analysis suggested that common genetic variants may modify the effects of rare deleterious alleles. To further support this, we note that both common alleles deviated from Hardy Weinberg equilibrium (p= 9.0e-8 for *SPTA1* lead SNP, p=7.28e-17 for *ANK1* lead SNP), suggesting that they could be under positive selection. The expression-increasing variants identified here may help to explain the heterogeneity that has been observed in HS patients.

In addition to genes specific to red cell turnover in the spleen, our genetic study of spleen iron also pointed to known and novel regulators of the human body’s iron economy. We characterized a spleen iron signal at *MS4A7*, which belongs to the CD20 family of membrane proteins, which are expressed within the hematopoietic lineage and largely uncharacterized.42,43 Here, we linked a splicing quantitative trait locus in *MS4A7* to spleen iron, macrophage abundance, and even found that this gene was enriched in the spleen and in macrophages. It is possible that excess circulating monocytes can be recruited to provide an expanded reservoir in the spleen, contributing to the body’s iron economy in addition to splenic macrophages44. Our findings regarding *MS4A7* potentially illuminate additional details of splenic erythrocyte clearance mechanisms. For example, the specific functional roles played by macrophages and monocytes in the splenic red pulp remain incompletely understood.44,9

This study has limitations. First, as the only study to quantify spleen iron in a large cohort, no replication cohort is available, though we are able to replicate all our findings on relevant blood cell traits in independent cohorts. Second, while we found limited evidence in the UK Biobank that spleen iron varies by ethnicity (**Supplementary Figures 3-4**), additional imaging studies are warranted to quantify spleen iron across populations. Third, experimental studies in the hematopoietic lineage of a model system will be needed to test the functional consequences of the regulatory variation observed here. Fourth, while we genetically identified HS in the UKBB and showed that the common alleles in *SPTA1* and *ANK1* act as genetic modifiers of erythrocyte parameters, additional clinical validation is needed to substantiate these findings in an HS patient cohort.

In summary, we have quantified spleen iron by repurposing liver MRI acquisitions, maximising use of the data at no extra cost to scanning or participant time. Our findings suggest that steady-state levels of spleen iron are sensitive to alterations in erythrocyte structure affecting cell turnover, as well as alterations in iron transport by macrophages. We identified common regulatory variation in HS genes showing effects on erythrocytes that are opposite of the effects observed in HS patients. Imaging-derived spleen iron is thus a novel quantitative biomarker that is tractable for genetic analysis, and has potential to contribute to future clinical characterization of hemolytic anemias.

## Methods

### Image Analysis

Image processing and segmentation pipelines are described in Liu et al.26 We extracted 2D masks from the 3D spleen segmentations at their intersection with the liver acquisition (**Figure 1A**).45 We applied one-pixel erosion before computing the median value within that mask. We excluded 2D masks that had < 1% of 3D volume or < 20 voxels. Additional methods are described in the Supplement.

### Epidemiological modelling of spleen iron risk factors

Associations of spleen iron with age, genetic sex, and self-reported ethnicity were performed in R v3.6.3. using the linear and logistic models for spleen iron as a quantitative or binary trait, respectively, after adjusting for covariates including imaging center, date and time.

### Phenome-wide association study

We generated a list of variables derived from raw data used PHESANT46, and removed procedural metrics (e.g. measurement date), duplicates, and raw measures, resulting in 1,824 traits (**Supplementary Table 2**). We used PheWAS47 to combine ICD10 codes (Field 41270) into distinct phenotype codes or phecodes (**Supplementary Table 3**). We performed linear (quantitative traits) or logistic regression (binary traits) on spleen iron, adjusting for imaging center, date, scan time, age, sex, BMI, height, and ethnicity.

### Genome-wide association study and replication analysis

We conducted a GWAS in n=35324 participants using methods as previously described.48 Additional methods are described in the Supplement. Summary statistics from the Blood Cell Consortium31 were downloaded from [http://www.mhi-humangenetics.org/en/resources/](http://www.mhi-humangenetics.org/en/resources/) and harmonized using dbSNP build 151 (GRCh37). Using relevant blood cell traits, we tested for replication using all available ancestries: Europeans, East Asians, African-Americans, and Hispanic/Latinos.

### Conditional analysis and fine-mapping

We performed conditional analysis using GCTA49, considering variants within 500kb of an index variant. We constructed a reference LD panel of 5,000 randomly selected, unrelated European UKBB participants24. We excluded the major histocompatibility complex region due to the complexity of LD structure at this locus (GRCh37::6:28,477,797-33,448,354; see [https://www.ncbi.nlm.nih.gov/grc/human/regions/MHC](https://www.ncbi.nlm.nih.gov/grc/human/regions/MHC)). For each locus, we considered variants with locus-wide evidence of association (pjoint<10−6) to be conditionally independent. We followed an iterative procedure to determine credible sets of causal variants with 95% coverage.48

### Colocalization studies

For gene expression studies, we used summary statistics from GTEx v850. For disease and quantitative trait studies, we used UKBB summary statistics of phecodes51, normalized quantitative traits ([http://www.nealelab.is/blog/2017/7/19/rapid-gwas-of-thousands-of-phenotypes-for-337000-samples-in-the-uk-biobank](http://www.nealelab.is/blog/2017/7/19/rapid-gwas-of-thousands-of-phenotypes-for-337000-samples-in-the-uk-biobank)). We selected phenotypes with p<5×10−8 within 500kb of the index variant. We performed colocalization analysis using coloc52 with default priors and considered variants within 500kb of the index variant. We considerd two genetic signals to have strong evidence of colocalization if PP3+PP4≥0.99 and PP4/PP3≥5.53

## Supporting information

Supplementary Information [[supplements/266431_file02.pdf]](pending:yes)

Supplementary Tables [[supplements/266431_file03.xlsx]](pending:yes)

## Data Availability

Model weights for spleen segmentation are available via Github ([www.github.com/calico-ukbb-mri-sseg](http://www.github.com/calico-ukbb-mri-sseg)) and code for resampling at [www.github.com/recoh/pipeline](http://www.github.com/recoh/pipeline). All summary statistics will be made available from the GWAS Catalog ([www.ebi.ac.uk/gwas](http://www.ebi.ac.uk/gwas)) upon publication. All derived data will be made available from the UK Biobank ([www.ukbiobank.ac.uk](http://www.ukbiobank.ac.uk)).

## Author Contributions

NB, MC, EPS, JDB, ELT designed the study. NB, BW, and YL implemented the image processing methods. MC performed the data processing. EPS and MC performed genetic analysis. ELT, EPS, RLC, MC and NB drafted the manuscript. All authors edited, read, and approved the manuscript.

## Declaration of interests

EPS, RLC and MC are employees of Calico Life Sciences LLC. YL is a former employee of Calico Life Sciences LLC. NB, BW, JDB, and ELT have no competing interests.

## Data Availability

Model weights for spleen segmentation are available via Github ([www.github.com/calico-ukbb-mri-sseg](http://www.github.com/calico-ukbb-mri-sseg)) and code for resampling at [www.github.com/recoh/pipeline](http://www.github.com/recoh/pipeline). All summary statistics will be made available from the GWAS Catalog ([www.ebi.ac.uk/gwas](http://www.ebi.ac.uk/gwas)) upon publication. All derived data will be made available from the UK Biobank ([www.ukbiobank.ac.uk](http://www.ukbiobank.ac.uk)).

## Acknowledgments

The authors would like to thank David Botstein, Marcia Paddock, Amoolya Singh, Kevin Wright, and Bert van der Zwaag for helpful feedback and discussion. This work was made possible by the UK Biobank, including staff, funders, and study volunteers. This research has been conducted using the UKBB Resource under Application Number 44584 and was funded by Calico Life Sciences LLC.

*   Received November 16, 2021.
*   Revision received January 18, 2022.
*   Accepted January 19, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

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

## References

1.  1.Cesta, M.F. (2006). Normal Structure, Function, and Histology of the Spleen. Toxicologic Pathology 34, 455– 465.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/01926230600867743&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17067939&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241520700004&link_type=ISI) 

2.  2.Dzierzak, E., and Philipsen, S. (2013). Erythropoiesis: Development and Differentiation. Cold Spring Harbor Perspectives in Medicine 3, a011601–a011601.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNzaHBlcnNwZWN0bWVkIjtzOjU6InJlc2lkIjtzOjExOiIzLzQvYTAxMTYwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzE5LzIwMjEuMTEuMTYuMjEyNjY0MzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

3.  3.Mebius, R.E., and Kraal, G. (2005). Structure and function of the spleen. Nat. Rev. Immunol. 5, 606–616.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri1669&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16056254&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230880900012&link_type=ISI) 

4.  4.Ganz, T. (2013). Systemic iron homeostasis. Physiol. Rev. 93, 1721–1741.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.00008.2013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24137020&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

5.  5.Gallagher, P.G. (2005). Red cell membrane disorders. Hematology Am. Soc. Hematol. Educ. Program 13–18.
    
    

6.  6.Delaunay, J. (2002). Molecular basis of red cell membrane disorders. Acta Haematol. 108, 210–218.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000065657&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12432217&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179701300004&link_type=ISI) 

7.  7.Bennett, V., and Baines, A.J. (2001). Spectrin and ankyrin-based pathways: metazoan inventions for integrating cells into tissues. Physiol. Rev. 81, 1353–1392.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/physrev.2001.81.3.1353&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11427698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169570800010&link_type=ISI) 

8.  8.Ovchynnikova, E., Aglialoro, F., von Lindern, M., and van den Akker, E. (2018). The Shape Shifting Story of Reticulocyte Maturation. Frontiers in Physiology 9,.
    
    

9.  9.Klei, T., Van Bruggen, R., Dalimot, J., Veldthuis, M., Mul, E., Rademakers, T., Hoogenboezem, M., Nagelkerke, S.Q., Van Ijcken, W., Moestrup, S.K., et al. (2019). Hemolysis in the Spleen Drives Erythrocyte Turnover. Blood 134, 946–946.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMDoiMTM0LzEyLzk0NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzE5LzIwMjEuMTEuMTYuMjEyNjY0MzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

10. 10.Arora, R.D., Dass, J., Maydeo, S., Arya, V., Kotwal, J., and Bhargava, M. (2018). Utility of mean sphered cell volume and mean reticulocyte volume for the diagnosis of hereditary spherocytosis. Hematology 23, 413–416.
    
    

11. 11.van Vuren, A., van der Zwaag, B., Huisjes, R., Lak, N., Bierings, M., Gerritsen, E., van Beers, E., Bartels, M., and van Wijk, R. (2019). The Complexity of Genotype-Phenotype Correlations in Hereditary Spherocytosis: A Cohort of 95 Patients. HemaSphere 3, e276.
    
    

12. 12.Chonat, S., Risinger, M., Sakthivel, H., Niss, O., Rothman, J.A., Hsieh, L., Chou, S.T., Kwiatkowski, J.L., Khandros, E., Gorman, M.F., et al. (2019). The Spectrum of SPTA1-Associated Hereditary Spherocytosis. Frontiers in Physiology 10,.
    
    

13. 13.Tole, S., Dhir, P., Pugi, J., Drury, L.J., Butchart, S., Fantauzzi, M., Langer, J.C., Baker, J.M., Blanchette, V.S., Kirby-Allen, M., et al. (2020). Genotype–phenotype correlation in children with hereditary spherocytosis. British Journal of Haematology 191, 486–496.
    
    

14. 14.Kolnagou, A., Michaelides, Y., Kontoghiorghe, C.N., and Kontoghiorghes, G.J. (2013). The importance of spleen, spleen iron, and splenectomy for determining total body iron load, ferrikinetics, and iron toxicity in thalassemia major patients. Toxicol. Mech. Methods 23, 34–41.
    
    

15. 15.Brewer, C.J., Coates, T.D., and Wood, J.C. (2009). Spleen R2 and R2* in iron-overloaded patients with sickle cell disease and thalassemia major. J. Magn. Reson. Imaging 29, 357–364.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmri.21666&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19161188&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

16. 16.Papakonstantinou, O., Alexopoulou, E., Economopoulos, N., Benekos, O., Kattamis, A., Kostaridou, S., Ladis, V., Efstathopoulos, E., Gouliamos, A., and Kelekis, N.L. (2009). Assessment of iron distribution between liver, spleen, pancreas, bone marrow, and myocardium by means of R2 relaxometry with MRI in patients with beta-thalassemia major. J. Magn. Reson. Imaging 29, 853–859.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmri.21707&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19306409&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

17. 17.Aslan, E., Luo, J.W., Lesage, A., Paquin, P., Cerny, M., Chin, A.S.-L., Olivié, D., Gilbert, G., Soulières, D., and Tang, A. (2021). MRI-based R2* mapping in patients with suspected or known iron overload. Abdom Radiol (NY).
    
    

18. 18.Çetinçakmak, M.G., Hattapoglu, S., Söker, M., Ekici, F., Yilmaz, K., Göya, C., and Hamidi, C. (2020). Evaluation of the relationship between splenic iron overload and liver, heart and muscle features evident on T2*-weighted magnetic resonance imaging. Adv. Clin. Exp. Med. 29, 475–480.
    
    

19. 19.França, M., Martí-Bonmatí, L., Porto, G., Silva, S., Guimarães, S., Alberich-Bayarri, Á., Vizcaíno, J.R., and Pessegueiro Miranda, H. (2018). Tissue iron quantification in chronic liver diseases using MRI shows a relationship between iron accumulation in liver, spleen, and bone marrow. Clin. Radiol. 73, 215.e1–e215.e9.
    
    

20. 20.Adler, D.D., Glazer, G.M., and Aisen, A.M. (1986). MRI of the spleen: normal appearance and findings in sickle-cell anemia. AJR Am. J. Roentgenol. 147, 843–845.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3489394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986E249700044&link_type=ISI) 

21. 21.Schwenzer, N.F., Machann, J., Haap, M.M., Martirosian, P., Schraml, C., Liebig, G., Stefan, N., Häring, H.-U., Claussen, C.D., Fritsche, A., et al. (2008). T2* relaxometry in liver, pancreas, and spleen in a healthy cohort of one hundred twenty-nine subjects-correlation with age, gender, and serum ferritin. Invest. Radiol. 43, 854–860.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/RLI.0b013e3181862413&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19002057&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

22. 22.Angelucci, E., Giovagnoni, A., Valeri, G., Paci, E., Ripalti, M., Muretto, P., McLaren, C., Brittenham, G.M., and Lucarelli, G. (1997). Limitations of magnetic resonance imaging in measurement of hepatic iron. Blood 90, 4736– 4742.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMDoiOTAvMTIvNDczNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzE5LzIwMjEuMTEuMTYuMjEyNjY0MzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

23. 23.Saito, H. (2014). METABOLISM OF IRON STORES. Nagoya J. Med. Sci. 76, 235–254.
    
    

24. 24.Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L.T., Sharp, K., Motyer, A., Vukcevic, D., Delaneau, O., O’Connell, J., et al. (2018). The UK Biobank resource with deep phenotyping and genomic data. Nature 562, 203–209.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

25. 25.Littlejohns, T.J., Holliday, J., Gibson, L.M., Garratt, S., Oesingmann, N., Alfaro-Almagro, F., Bell, J.D., Boultwood, C., Collins, R., Conroy, M.C., et al. (2020). The UK Biobank imaging enhancement of 100,000 participants: rationale, data collection, management and future directions. Nat. Commun. 11, 2624.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-15948-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32457287&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

26. 26.Liu, Y., Basty, N., Whitcher, B., Bell, J.D., Sorokin, E.P., van Bruggen, N., Louise Thomas, E., and Cule, M. (2021). Genetic architecture of 11 organ traits derived from abdominal MRI using deep learning. eLife 10,.
    
    

27. 27.Henninger, B., Alustiza, J., Garbowski, M., and Gandon, Y. (2020). Practical guide to quantification of hepatic iron with MRI. Eur. Radiol. 30, 383–393.
    
    

28. 28.Elliott, L.T., Sharp, K., Alfaro-Almagro, F., Shi, S., Miller, K.L., Douaud, G., Marchini, J., and Smith, S.M. (2018). Genome-wide association studies of brain imaging phenotypes in UK Biobank. Nature 562, 210–216.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0571-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

29. 29.Bell, S., Rigas, A.S., Magnusson, M.K., Ferkingstad, E., Allara, E., Bjornsdottir, G., Ramond, A., Sørensen, E., Halldorsson, G.H., Paul, D.S., et al. (2021). A genome-wide meta-analysis yields 46 new loci associating with biomarkers of iron homeostasis. Commun Biol 4, 156.
    
    

30. 30.Consortium, T.G., and The GTEx Consortium (2020). The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 369, 1318–1330.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNjkvNjUwOS8xMzE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMTkvMjAyMS4xMS4xNi4yMTI2NjQzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

31. 31.Chen, M.-H., Raffield, L.M., Mousas, A., Sakaue, S., Huffman, J.E., Moscati, A., Trivedi, B., Jiang, T., Akbari, P., Vuckovic, D., et al. (2020). Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations. Cell 182, 1198–1213.e14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.06.045&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32888493&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

32. 32.Uhlen, M., Karlsson, M.J., Zhong, W., Tebani, A., Pou, C., Mikes, J., Lakshmikanth, T., Forsström, B., Edfors, F., Odeberg, J., et al. (2019). A genome-wide transcriptomic analysis of protein-coding genes in human blood cells. Science 366,.
    
    

33. 33.Sirlin, C.B., and Reeder, S.B. (2010). Magnetic resonance imaging quantification of liver iron. Magn. Reson. Imaging Clin. N. Am. 18, 359–381, ix.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mric.2010.08.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21094445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

34. 34.Fernandes, J.L. (2018). MRI for Iron Overload in Thalassemia. Hematol. Oncol. Clin. North Am. 32, 277–295.
    
    

35. 35.Arrivé, L., Thurnher, S., Hricak, H., and Price, D.C. (1990). Magnetic resonance imaging of splenic iron overload. Eur. J. Radiol. 10, 98–104.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2338103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

36. 36.Ooi, G.C., Chen, F.E., Chan, K.N., Tsang, K.W., Wong, Y.H., Liang, R., Chan, V., and Ngan, H. (1999). Qualitative and quantitative magnetic resonance imaging in haemoglobin H disease: screening for iron overload. Clin. Radiol. 54, 98–102.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0009-9260(99)91068-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10050737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078471600007&link_type=ISI) 

37. 37.Rushton, D.H., and Barth, J.H. (2010). What is the evidence for gender differences in ferritin and haemoglobin? Crit. Rev. Oncol. Hematol. 73, 1–9.
    
    

38. 38.Timmers, P.R.H.J., Wilson, J.F., Joshi, P.K., and Deelen, J. (2020). Multivariate genomic scan implicates novel loci and haem metabolism in human ageing. Nat. Commun. 11, 3570.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-17312-3&link_type=DOI) 

39. 39.McKay, A., Wilman, H.R., Dennis, A., Kelly, M., Gyngell, M.L., Neubauer, S., Bell, J.D., Banerjee, R., and Thomas, E.L. (2018). Measurement of liver iron by magnetic resonance imaging in the UK Biobank population. PLoS One 13, e0209340.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0209340&link_type=DOI) 

40. 40.Harrison-Findik, D.D. (2007). Role of alcohol in the regulation of iron metabolism. World Journal of Gastroenterology 13, 4925.
    
    

41. 41.Pietrangelo, A. (2017). Ferroportin disease: pathogenesis, diagnosis and treatment. Haematologica 102, 1972–1984.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiaGFlbWF0b2wiO3M6NToicmVzaWQiO3M6MTE6IjEwMi8xMi8xOTcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMTkvMjAyMS4xMS4xNi4yMTI2NjQzMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

42. 42.Silva-Gomes, R., Mapelli, S.N., Boutet, M.-A., Mattiola, I., Sironi, M., Grizzi, F., Colombo, F., Supino, D., Carnevale, S., Pasqualini, F., et al. (2021). Differential expression and regulation of MS4A family members in myeloid cells in physiological and pathological conditions. J. Leukoc. Biol.
    
    

43. 43.Kuek, L.E., Leffler, M., Mackay, G.A., and Hulett, M.D. (2016). The MS4A family: counting past 1, 2 and 3. Immunology & Cell Biology 94, 11–23.
    
    

44. 44.Swirski, F.K., Nahrendorf, M., Etzrodt, M., Wildgruber, M., Cortez-Retamozo, V., Panizzi, P., Figueiredo, J.-L., Kohler, R.H., Chudnovskiy, A., Waterman, P., et al. (2009). Identification of splenic reservoir monocytes and their deployment to inflammatory sites. Science 325, 612–616.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMjUvNTk0MC82MTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMS8xOS8yMDIxLjExLjE2LjIxMjY2NDMxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

45. 45.Basty, N., Liu, Y., Cule, M., Louise Thomas, E., Bell, J.D., and Whitcher, B. (2020). Automated Measurement of Pancreatic Fat and Iron Concentration Using Multi-Echo and T1-Weighted MRI Data. 2020 IEEE 17th International Symposium on Biomedical Imaging (ISBI).
    
    

46. 46.Millard, L.A.C., Davies, N.M., Gaunt, T.R., Davey Smith, G., and Tilling, K. (2018). Software Application Profile: PHESANT: a tool for performing automated phenome scans in UK Biobank. Int. J. Epidemiol. 47, 29–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyx204&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29040602&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

47. 47.Carroll, R.J., Bastarache, L., and Denny, J.C. (2014). R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment. Bioinformatics 30, 2375–2376.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btu197&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24733291&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000342746000018&link_type=ISI) 

48. 48.Sethi, A., Taylor, L., Graham Ruby, J., Venkataraman, J., Cule, M., and Melamud, E. Calcification of abdominal aorta is an underappreciated cardiovascular disease risk factor.
    
    

49. 49.Yang, J., Ferreira, T., Morris, A.P., Medland, S.E., Genetic Investigation of ANthropometric Traits (GIANT) Consortium, DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium, Madden, P.A.F., Heath, A.C., Martin, N.G., Montgomery, G.W., et al. (2012). Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat. Genet. 44, 369–375, S1–S3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2213&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22426310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

50. 50.GTEx Consortium, Laboratory, Data Analysis &Coordinating Center (LDACC)—Analysis Working Group, Statistical Methods groups—Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, NIH/NCI, NIH/NHGRI, NIH/NIMH, NIH/NIDA, Biospecimen Collection Source Site—NDRI, et al. (2017). Genetic effects on gene expression across human tissues. Nature 550, 204–213.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29022597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000412829500039&link_type=ISI) 

51. 51.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%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

52. 52.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%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

53. 53.Guo, H., Fortune, M.D., Burren, O.S., Schofield, E., Todd, J.A., and Wallace, C. (2015). Integration of disease association and eQTL data using a Bayesian colocalisation approach highlights six candidate causal genes in immune-mediated diseases. Hum. Mol. Genet. 24, 3305–3313.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddv077&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25743184&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

54. 54.Reeder, S.B., Pineda, A.R., Wen, Z., Shimakawa, A., Yu, H., Brittain, J.H., Gold, G.E., Beaulieu, C.H., and Pelc, N.J. (2005). Iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL): Application with fast spin-echo imaging. Magnetic Resonance in Medicine 54, 636–644.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/mrm.20624&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16092103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

55. 55.Hong, C.W., Hamilton, G., Hooker, C., Park, C.C., Tran, C.A., Henderson, W.C., Hooker, J.C., Fazeli Dehkordy, S., Schwimmer, J.B., Reeder, S.B., et al. (2019). Measurement of spleen fat on MRI-proton density fat fraction arises from reconstruction of noise. Abdom Radiol (NY) 44, 3295–3303.
    
    

56. 56.Wood, J.C., Enriquez, C., Ghugre, N., Tyzka, J.M., Carson, S., Nelson, M.D., and Coates, T.D. (2005). MRI R2 and R2* mapping accurately estimates hepatic iron concentration in transfusion-dependent thalassemia and sickle cell disease patients. Blood 106, 1460–1465.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMDoiMTA2LzQvMTQ2MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzE5LzIwMjEuMTEuMTYuMjEyNjY0MzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

57. 57.Labranche, R., Gilbert, G., Cerny, M., Vu, K.-N., Soulières, D., Olivié, D., Billiard, J.-S., Yokoo, T., and Tang, A. (2018). Liver Iron Quantification with MR Imaging: A Primer for Radiologists. Radiographics 38, 392–412.
    
    

58. 58.Loh, P.-R., Tucker, G., Bulik-Sullivan, B.K., Vilhjálmsson, B.J., Finucane, H.K., Salem, R.M., Chasman, D.I., Ridker, P.M., Neale, B.M., Berger, B., et al. (2015). Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat. Genet. 47, 284–290.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3190&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

59. 59.Yang, J., Benyamin, B., McEvoy, B.P., Gordon, S., Henders, A.K., Nyholt, D.R., Madden, P.A., Heath, A.C., Martin, N.G., Montgomery, G.W., et al. (2010). Common SNPs explain a large proportion of heritability for human height. Nat. Genet. 42, 565.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.608&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20562875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279242400007&link_type=ISI) 

60. 60.Bulik-Sullivan, B.K., Finucane, H.K., Anttila, V., Gusev, A., Day, F.R., Loh, P.-R., ReproGen Consortium, Psychiatric Genomics Consortium, Genetic Consortium for Anorexia Nervosa of the Wellcome Trust Case Control Consortium 3, Duncan, L., et al. (2015). An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47, 1236–1241.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

61. 61.Regier, A.A., Farjoun, Y., Larson, D.E., Krasheninina, O., Kang, H.M., Howrigan, D.P., Chen, B.-J., Kher, M., Banks, E., Ames, D.C., et al. (2018). Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects. Nat. Commun. 9, 4038.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-06159-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30279509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

62. 62.Karczewski, K.J., Francioli, L.C., Tiao, G., Cummings, B.B., Alföldi, J., Wang, Q., Collins, R.L., Laricchia, K.M., Ganna, A., Birnbaum, D.P., et al. (2020). The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

63. 63.Kircher, M., Witten, D.M., Jain, P., O’Roak, B.J., Cooper, G.M., and Shendure, J. (2014). A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46, 310–315.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2892&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24487276&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F19%2F2021.11.16.21266431.atom) 

64. 64.Zhou, W., Zhao, Z., Nielsen, J.B., Fritsche, L.G., LeFaive, J., Gagliano Taliun, S.A., Bi, W., Gabrielsen, M.E., Daly, M.J., Neale, B.M., et al. (2020). Scalable generalized linear mixed model for region-based association tests in large biobanks and cohorts. Nat. Genet. 52, 634–639.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-0621-6&link_type=DOI) 

65. 65.St Pierre, T.G., Clark, P.R., Chua-anusorn, W., Fleming, A.J., Jeffrey, G.P., Olynyk, J.K., Pootrakul, P., Robins, E., and Lindeman, R. (2005). Noninvasive measurement and imaging of liver iron concentrations using proton magnetic resonance. Blood 105, 855–861.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czo5OiIxMDUvMi84NTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMS8xOS8yMDIxLjExLjE2LjIxMjY2NDMxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)