If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Corresponding authors. Addresses: Institute of Health Informatics, Faculty of Population Health Sciences, University College London, NW12DA. Tel.: +(44)7899786998
Perspectum Diagnostics Ltd., Oxford, UKOxford Centre for Clinical Magnetic Resonance Research, Division of Cardiovascular Medicine, Oxford NIHR Biomedical Research Centre, University of Oxford, Oxford, UK
Research Centre for Optimal Health, School of Life Sciences, University of Westminster, London, UKGenetics of Complex Traits, College of Medicine and Health, University of Exeter, Exeter, UKDivision of Medical Sciences, Department of Health Sciences, Luleå University of Technology, Luleå, Sweden
Variants in metal ion transporter and NAFLD genes are associated with liver MRI-derived cT1, a steatohepatitis and fibrosis proxy.
•
cT1 is highly heritable, and is correlated with BMI, NAFLD and VLDL, and inversely correlated with HDL.
•
Insulin resistance, NAFLD and higher BMI are genetically linked to higher liver cT1, whilst favourable adiposity is linked to lower cT1.
Background & Aims
MRI-based corrected T1 (cT1) is a non-invasive method to grade the severity of steatohepatitis and liver fibrosis. We aimed to identify genetic variants influencing liver cT1 and use genetics to understand mechanisms underlying liver fibroinflammatory disease and its link with other metabolic traits and diseases.
Methods
First, we performed a genome-wide association study (GWAS) in 14,440 Europeans, with liver cT1 measures, from the UK Biobank. Second, we explored the effects of the cT1 variants on liver blood tests, and a range of metabolic traits and diseases. Third, we used Mendelian randomisation to test the causal effects of 24 predominantly metabolic traits on liver cT1 measures.
Results
We identified 6 independent genetic variants associated with liver cT1 that reached the GWAS significance threshold (p <5×10-8). Four of the variants (rs759359281 in SLC30A10, rs13107325 in SLC39A8, rs58542926 in TM6SF2, rs738409 in PNPLA3) were also associated with elevated aminotransferases and had variable effects on liver fat and other metabolic traits. Insulin resistance, type 2 diabetes, non-alcoholic fatty liver and body mass index were causally associated with elevated cT1, whilst favourable adiposity (instrumented by variants associated with higher adiposity but lower risk of cardiometabolic disease and lower liver fat) was found to be protective.
Conclusion
The association between 2 metal ion transporters and cT1 indicates an important new mechanism in steatohepatitis. Future studies are needed to determine whether interventions targeting the identified transporters might prevent liver disease in at-risk individuals.
Lay summary
We estimated levels of liver inflammation and scarring based on magnetic resonance imaging of 14,440 UK Biobank participants. We performed a genetic study and identified variations in 6 genes associated with levels of liver inflammation and scarring. Participants with variations in 4 of these genes also had higher levels of markers of liver cell injury in blood samples, further validating their role in liver health. Two identified genes are involved in the transport of metal ions in our body. Further investigation of these variations may lead to better detection, assessment, and/or treatment of liver inflammation and scarring.
It has come to our attention that there are errors in the Table 2 of the original manuscript “Genome-wide and Mendelian randomisation studies of liver MRI yield insights into the pathogenesis of steatohepatitis”.
Prevalence of clinically significant liver disease within the general population, as defined by non-invasive markers of liver fibrosis: a systematic review.
However, an important paradox in the history of liver fat accumulation exists; despite the large proportion of adults affected by simple steatosis (fatty liver), only a relatively small proportion (2.4–12.8%) will experience significant liver disease or liver-related death.
It is important to identify which individuals are at risk of developing the more inflammatory phenotype, steatohepatitis (a condition characterised by lipotoxicity and histological necroinflammation), which is considered to be the main pathophysiological driver of liver fibrosis and subsequent disease progression.
Prevalence of nonalcoholic fatty liver disease and nonalcoholic steatohepatitis among a largely middle-aged population utilizing ultrasound and liver biopsy: a prospective study.
T1 relaxation time reflects extracellular fluid, which is characteristic of fibrosis and inflammation. The presence of iron, which can be determined from T2∗ maps, has an opposing effect. Combining T2∗ and T1 values can correct for this opposing effect, from which cT1 (in milliseconds) is derived. Higher cT1 values are associated with both histological liver inflammation and fibrosis, although their relative contributions to the score are still unknown.
Understanding the underlying genetic susceptibility to steatohepatitis and fibrosis may provide new insights into the main pathophysiological mechanisms that contribute to chronic liver disease, helping in the identification of new drug targets. Genetic studies have so far been limited due to the phenotyping challenge. Liver biopsy is an invasive procedure with associated risks, significant sampling error and marked interobserver variance,
while routinely available liver blood tests such as aminotransferases, despite being useful in the identification of important liver disease susceptibility loci, are overall poor predictors of liver disease severity.
The diagnosis and management of non-alcoholic fatty liver disease: practice guideline by the American Gastroenterological Association, American Association for the Study of Liver Diseases, and American College of Gastroenterology.
Another challenging question is which metabolic traits cause steatohepatitis since treating causal factors can help prevent liver disease. Observational associations between steatohepatitis and other features of the metabolic syndrome might occur because they share common risk factors, rather than one causing the other. Mendelian randomisation is an established epidemiological approach that uses genetic studies to provide insight on causality.
Mendelian randomisation uses genetic variants associated with an exposure (e.g. body mass index [BMI], LDL cholesterol, insulin resistance) to assess their causal effect on an outcome of interest (e.g. cT1, steatohepatitis). Genetic markers of a risk factor are largely independent of confounders that may otherwise cause bias since genetic variants are randomly allocated before birth. Furthermore, the non-modifiable nature of genetic variants provides an analogy to randomised trials, in which exposure is allocated randomly and is non-modifiable by subsequent disease.
In this study, we aimed to (i) identify genetic variants influencing liver cT1 (ii) identify the effect of liver cT1 variants on other metabolic traits, (iii) investigate which metabolic traits are genetically correlated with cT1 measures and (iv) use Mendelian randomisation to investigate whether 24 metabolic traits and conditions are causally associated with cT1. We performed the first genome-wide association study (GWAS) on MRI liver cT1 in 14,440 European individuals from the UK Biobank. Finally, to investigate whether there are shared variants between liver cT1 and liver fat, we carried out a GWAS on MRI determined liver proton density fat fraction (PDFF) in the same cohort.
Materials and methods
UK Biobank participants
UK Biobank is a prospective cohort study that consists of over 500,000 individuals aged 37–73 years (99.5% were between 40 and 69 years of age) who were recruited between 2006 and 2010 from across the UK.
This research has been conducted using the data obtained via UK Biobank Access Application number 9914. The UK Biobank has approval from the North West Multi-Centre Research Ethics Committee (ref: 11/NW/0382) and obtained written informed consent from all participants prior to the study.
Imaging protocol and analysis
Invitation to the UK Biobank imaging study is based only on proximity to one of the main imaging sites. Participants were invited and scanned at the UK Biobank Imaging Centre in Cheadle (UK) using a Siemens 1.5T Magnetom Aera as previously described.
Reference range of liver corrected T1 values in a population at low risk for fatty liver disease-a UK Biobank sub-study, with an appendix of interesting cases.
Medical conditions were not taken into account except from those which would exclude the participant from being able to have an MRI (e.g. if they had an implanted defibrillator or metal implant).
Characterisation of cT1 in the UK Biobank cohort, alongside normal values and inter- and intra-reader variability have previously been published.
Reference range of liver corrected T1 values in a population at low risk for fatty liver disease-a UK Biobank sub-study, with an appendix of interesting cases.
Briefly, 2 sequences were used to acquire data: a shortened modified look locker inversion (ShMOLLI) to quantify liver T1, and a multiecho-spoiled gradient-echo, to quantify liver iron and fat (PDFF). In both cases, data was acquired as a single transverse slice captured through the centre of the liver superior to the porta hepatis. Acquisition was performed in end-expiration breath-hold and without the aid of any contrast agent injection. The slice-based methodology has previously been shown to correlate well with histology and predict liver-related outcomes.
Reference range of liver corrected T1 values in a population at low risk for fatty liver disease-a UK Biobank sub-study, with an appendix of interesting cases.
All values reported in this work are the median, for each metric, of all usable voxels in the liver within the image slice. T1 relaxation time reflects extracellular fluid and is characteristic of fibrosis and inflammation. The presence of iron, which can be determined from T2∗ maps, has an opposing effect on the T1, and algorithms have been formed to correct for the resulting bias.
All processed data are available through application to the UK Biobank. Fig. 1A illustrates the 3 MRI scans with different levels of cT1 in 3 participants.
From an initial collection of 20,386 imaging sessions (each of a unique individual), 691 did not have all the necessary imaging data, 1,354 were run with an early flawed protocol, 1,717 did not correctly trigger the sequence, 126 had more than half of their liver excluded due to poor model fitting and motion artefacts, leaving 16,498 for human quality control.
From these, a further 959 were removed through a combination of fat/water swaps, erroneous overcorrection of iron, misplacement of the image slice, segmentation failure, field artefacts, and cysts within the image slice that prevented reasonable quantification of parenchyma, leaving 15,539 participants.
Genetic data
Protocols for the participant genotyping, data collection, and quality control have previously been described in detail.
Briefly, participants were genotyped using 1 of 2 purpose-designed arrays (UK BiLEVE Axiom Array [n = 50,520] and UK Biobank Axiom Array [n = 438,692]) with 95% marker overlap. We excluded individuals who were identified by the UK Biobank as outliers based on either genotyping missingness rate or heterogeneity, or whose sex inferred from the genotypes did not match their self-reported sex. We removed individuals with a missingness >5% across variants which passed our quality control procedure. We used the latest release which included imputed data using 2 reference panels: a combined UK10K and 1000 Genomes panel and the Haplotype Reference Consortium panel. We limited our analysis to genetic variants with a minimum minor allele frequency (MAF) >1% and imputation quality score >0.3.
To define “white European” ancestry, we first used data from 1000 genomes samples to generate ancestry informative principal components (PCs). We then used these PCs in UK Biobank participants and employed K-means clustering to identify samples clustered with the 3 main 1000 genomes populations (European, African, and South Asian). Those clustered with the 1000 genomes’ “European” cluster were classified as having European ancestry.
In total, after image analysis and quality control steps, liver cT1 and PDFF measures were available for 14,440 white European individuals who also had genetic data available and were classified as white European.
Genome-wide association analysis
We used BOLT-LMM v2.3.4 to conduct a linear mixed model GWAS which accounts for population structure and relatedness. We increased our power by including all related individuals of European descent (n = 14,440). The relatedness matrix was computed using common (MAF >5%) genotyped variants that passed quality control in all 106 batches and were present on both genotyping arrays. Prior to association testing, liver cT1 and PDFF were inverse-normal transformed. We used age, sex, centre and genotyping arrays as covariates in the model.
Sensitivity analyses
We performed 6 sensitivity analyses (Table S1). We carried out GWASs and adjusted for (i) BMI and (ii) alcohol units consumed. We derived an alcohol units per day variable from the UK Biobank as previously suggested.
In summary, a 125 ml glass of wine (red, white, or sparkling) was considered to be 1.5 units, a pint of beer or cider was considered to be 2.8 units, other alcoholic drinks (e.g. alcopops) were considered to be 1.5 units, and a measure of spirit was considered to be 1 unit. We further adjusted for (iii) MRI-determined liver fat and (iv) liver iron to rule out the confounding effects of these 2 traits in our image processing pipeline. Finally, we carried out GWASs in (v) males and (vi) females separately to detect sex-specific associations.
Association of cT1 variants with liver biomarkers and metabolic traits and diseases
To further understand the role of each cT1 variant in the pathophysiology of liver disease, and also as a positive control, we tested the association between each variant and liver biomarkers in white European participants from the UK Biobank. We measured the following liver biomarkers: liver enzymes (alanine aminotransferase [ALT], aspartate aminotransferase [AST], gamma glutamyltransferase, alkaline phosphatase in up to 378,821 individuals), MRI-derived liver PDFF (n = 14,440), and MRI-derived liver iron (to understand if the correction of T1 measures for liver iron content has caused any bias; n = 14,440). The protocols for the derivation of MRI PDFF and liver iron have previously been published.
To validate the associations with aminotransferases in a non-UK Biobank dataset, we looked up the effects of cT1 variants in an existing GWAS of ALT and AST levels in up to 61,089 individuals.
To understand the effect of cT1 variants on cardiometabolic traits and diseases, we tested their associations with 15 predominantly metabolic traits including BMI, HDL-cholesterol, LDL-cholesterol, triglycerides, systolic blood pressure, diastolic blood pressure, type 2 diabetes, and coronary artery disease in up to n = 379,308 white European UK Biobank participants.
LD score regression and cross-trait genetic correlation analysis
We used LD Hub to conduct linkage disequilibrium (LD) score regression and heritability analysis. LD Hub is a centralised database of summary level GWAS for >500 diseases and traits from publicly available resources/consortia and uses a web interface that automates LD score regression, heritability and cross-trait genetic correlation analysis.
LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis.
We ran heritability analysis as well as genetic correlation analysis across 120 potentially relevant traits. Single-nucleotide polymorphism (SNP)-based heritability (h2SNP) is the proportion of total variation in liver cT1 measures due to the additive genetic variation between individuals in our study population.
Liver cirrhosis variants
To investigate the effect of liver cirrhosis variants on cT1 measures, and also as a positive control, we used variants associated with all-cause cirrhosis including rs2642438 (in or near MARC1), rs72613567 (HSD17B13), rs58542926 (TM6SF2), rs738409 (PNPLA3), rs1800562 (HFE), and rs28929474 (SERPINA).
We used the inverse variance weighted approach (IVW) as our main analysis, and Mendelian randomisation-Egger and penalised weighted median as sensitivity analyses in order to detect unidentified pleiotropy of our genetic instruments. Genetic instruments were constructed by using the independent genome-wide significant genetic variants (R2 <0.1) of the exposure of interest from previous GWASs. For more information on Mendelian randomisation and genetic instrument selection please see the Supplementary Material.
Results
The characteristics of liver cT1 cohort
In our discovery cohort, median age was 57 years (interquartile range (IQR) 50–62) for males and 55 years (IQR 48–60) for females. The median liver cT1 was 694 ms (IQR 662–730) in males and 676 ms (IQR 647–710) in females (Fig. S1); 5.3% of males (299/5,595) and 2.6% of females (169/6,455) had values above 800 ms, a threshold that has been set in current clinical trials as a cut-off for steatohepatitis,
and is under evaluation by the FDA and EMA as a diagnostic enrichment biomarker for NASH. Baseline characteristics were comparable to the rest of the UK Biobank cohort who did not participate in the imaging study except BMI, waist circumference and diabetes prevalence which were lower in both males and females in the liver cT1 cohort compared to the rest of the UK Biobank (Table 1). Although invitation was not based on any medical information, MRI exclusion criteria (e.g. metal or electrical implants, surgery 6 weeks prior to appointment, severe hearing or breathing problems) and the imaging site location (Cheadle, UK) may have contributed to a slightly healthier cohort.
BMI (Mann-Whitney U test, p = 1×10−80), waist circumference (Mann-Whitney U test, p = 1×10−100), diabetes prevalence (Pearson's chi squared test, p = 1×10−27) were lower in the imaging subset compared to the rest of UK Biobank. Levels of significance for all tests: (p <0.05).
BMI (Mann-Whitney U test, p = 1×10−80), waist circumference (Mann-Whitney U test, p = 1×10−100), diabetes prevalence (Pearson's chi squared test, p = 1×10−27) were lower in the imaging subset compared to the rest of UK Biobank. Levels of significance for all tests: (p <0.05).
BMI (Mann-Whitney U test, p = 1×10−80), waist circumference (Mann-Whitney U test, p = 1×10−100), diabetes prevalence (Pearson's chi squared test, p = 1×10−27) were lower in the imaging subset compared to the rest of UK Biobank. Levels of significance for all tests: (p <0.05).
26.6 (24.5–28.8)
25 (22.9–28)
27.3 (25–30.1)
26.1 (23.5–30)
∗ BMI (Mann-Whitney U test, p = 1×10−80), waist circumference (Mann-Whitney U test, p = 1×10−100), diabetes prevalence (Pearson's chi squared test, p = 1×10−27) were lower in the imaging subset compared to the rest of UK Biobank. Levels of significance for all tests: (p <0.05).
Genetic variants in 6 loci show association with liver cT1
In our GWAS of liver cT1 in individuals of European ancestry, variants in 6 independent loci (Table 2) reached genome-wide significance. Genomic inflation was low (λGC = 1.006, Fig. S2). We observed the strongest association with a missense variant, rs13107325, located in an exon of SLC39A8 (Fig. 1B). The minor allele (T; allele frequency 7%) of rs13107325 was associated with 0.54 SD increase in cT1 (p = 1.2×10−133). The mean cT1 was 692 ms in individuals with no risk allele, 727 ms in heterozygotes, and 772 ms in risk allele homozygotes (Fig. S3).
Table 2The association between 6 independent genetic variants and liver cT1. A linear mixed model was used for genetic associations (levels of significance: p <5×10−8).
Other independent variants included an intronic variant (rs759359281-CA >C) in SLC30A10 (p = 2.8×10−8), a missense variant (rs111723834-G >A) in PCK2 (p = 3.0×10−11), a missense variant (rs4820268-A >G) in TMPRSS6 (p = 1.6×10−9), and 2 known cirrhosis variants (rs58542926-A >G) in TM6SF2 (p = 1.4×10−8) and (rs738409-C >G) in PNPLA3 (p = 9.6×10−13). The 6 variants together explained 5.38% of variation in cT1 measures in white European UK Biobank participants with the SLC39A8 variant explaining most of this variation (3.95%) (Table 2). We estimated the SNP-based heritability (h2SNP) of liver cT1 to be 20%. This is higher than the heritability estimated for conditions and traits such as coronary artery disease (7%),
We did not detect any sex-specific associations and the effects were similar between men and women (Table S1). Sensitivity analyses that further controlled for alcohol unit intake and BMI did not identify any additional signals and did not significantly change the effect size (Table S1). Sensitivity analyses that controlled for liver PDFF removed the effects of rs58542926 in TM6SF2 and rs738409 in PNPLA3, suggesting that the effects of these variants on cT1 measures are mediated through liver fat accumulation (Table S1). The cT1 increasing allele (G) at TMPRSS6-rs4820268 is associated with lower plasma iron levels and lower liver iron.
The effect of this variant on cT1 may be due to its effect on liver iron concentrations since iron has an opposing effect to T1 relaxation time. However, sensitivity analyses that controlled for liver iron only slightly attenuated its effect on cT1 (from beta = 0.066, p = 2×10−9 to beta = 0.054, p = 7×10−7) suggesting that other mechanisms are involved and that this is a true signal.
Genetic variants in 4 loci show association with liver MRI-determined PDFF
In our GWAS of liver PDFF in 14,440 individuals of European ancestry missense variants in 4 independent loci reached genome-wide significance (rs1260326-C >T in GCKR, p = 3.9×10−8, rs58542926-C >T in TM6SF2, p = 6.3×10−37, rs429358-C >T in APOE, p = 5.6×10−11, rs738409-C >G in PNPLA3, p = 5.4×10−66 (Table S2, Fig. S4). Genomic inflation was low (λGC = 1.04). Two of the 4 variants (rs738409 in PNPLA3, rs58542926 in TM6SF2) were shared between PDFF and cT1 in our GWASs.
Four of the cT1 variants are associated with higher levels of aminotransferases and demonstrate variable effects on metabolic traits and diseases
To validate these variants and further understand their role in other metabolic traits and diseases, we investigated their association with liver blood tests, MRI-determined liver iron and liver PDFF, lipids, blood pressure, BMI and cardiometabolic disease outcomes (Fig. 2, Table S3). cT1-increasing alleles at 4 variants (in SLC30A10, SLC39A8, TM6SF2, and PNPLA3) were associated with higher ALT and AST (all with p values <2×10−5) and higher risk of type 2 diabetes (all with p <0.002, except the SLC30A10 variant). None of cT1 variants were associated with cardiovascular disease risk, whilst their effects on other metabolic traits including lipids and blood pressure were variable (Fig. 2). Among the novel identified and replicated variants (rs759359281 in SLC30A10, and rs13107325 in SLC39A8), only the latter was available in a non-UK Biobank cohort with available liver blood tests. The cT1-increasing allele in rs13107325 showed a similar direction of effect on ALT (n = 46,316, beta = 0.01, p = 0.27) and AST (n = 39,015, beta = 0.014, p = 0.0005) levels in an independent cohort (Table S4).
Effects are in SD for continuous traits and log(OR) for disease outcomes per copy of the risk allele. A linear mixed model was used for genetic associations. Levels of significance: p <0.05. ALP, alkaline phosphatase; ALT, alanine aminotransferase; AST, aspartate aminotransferase; CAD, coronary artery disease; cT1, corrected T1; GGT, gamma-glutamyltransferase; HDL-C, HDL-cholesterol; LDL-C, LDL cholesterol; OR, odds ratio; T2DM, type 2 diabetes.
Liver cT1 measures correlate genetically with components of metabolic syndrome
We calculated genetic correlations using the GWAS summary statistics (120 predominantly metabolic traits/diseases) in LD score regression analysis (Fig. 3, Table S5). Measures of insulin resistance, triglycerides, VLDL, type 2 diabetes, coronary artery disease, body fat percentage, BMI and waist-to-hip ratio were positively genetically correlated with liver cT1 measures after correcting p values for multiple testing (false discovery rate <0.05). The most genetically correlated traits were homeostatic model for insulin resistance (HOMA-IR, rG = 0.53, p = 0.0004) and mean diameter of VLDL particles (rG = 0.52, p = 0.0004), whereas the strongest inverse correlation was seen with total cholesterol in very large HDL (rG = -0.62, p = 0.04).
Fig. demonstrating the significant genetic correlations (rg) between cT1 and metabolic traits following correction for multiple testing (levels of significance: p false discovery rate <0.05) among more than 120 traits. The colours correspond to significance of correlation (t test); red: p <1×10−8; orange: 1×10−6 <p <1×10−5; blue: 1×10−5 <p <1×10−4; green: 1×10−4 <p <1×10−3; yellow: 0.001 <p <0.01. Higher cT1 is positively genetically correlated with VLDL, type 2 diabetes, coronary artery disease, and inversely correlated with HDL. cT1, corrected T1; HOMA-IR, homeostatic model assessment of insulin resistance, HOMA-B, homeostatic model assessment of β cell function; LD, linkage disequilibrium. (This figure appears in color on the web.)
Association of liver cirrhosis variants with liver cT1
We investigated the effects of all-cause cirrhosis risk variants on cT1 values. Among 6 variants associated with all-cause cirrhosis in a recent GWAS of 5,770 cases and 572,850 controls,
4 variants (those in or near MARC1, HSD17B13, TM6SF2 and PNPLA3) demonstrated associations with cT1 (Table 3), where alleles associated with higher risk of liver cirrhosis were also associated with higher cT1. The HFE haemochromatosis risk allele (in rs1800562) was inversely associated with cT1, however this is to be expected since cT1 measures are corrected for liver iron content. Consistently, this association became remarkably attenuated (from beta = −0.11, p = 8×10−7 to beta = −0.055, p = 0.02) in our sensitivity analysis correcting for liver iron content. In the GWAS of all-cause cirrhosis, the effect of a1-antitrypsin risk variant (rs28929474 in SERPINA1) was very weak (p = 0.01) and present only when a recessive model was used (Table 3).
∗Indicates recessive models were run for the previously published all-cause cirrhosis GWAS; all other association analyses used additive models. Logistic regression was used for the genetic associations with cirrhosis; a linear mixed model was used for the genetic associations with cT1 (levels of significance: p <5x10−8, suggestive p <0.05).
† Beta cirrhosis is the effect on all-cause cirrhosis in log(OR) and Beta cT1 is the effect on cT1 in SD.
Mendelian randomisation analysis provides genetic evidence that non-alcoholic fatty liver, insulin resistance and obesity causally elevate liver cT1
Demonstrating causality using observational studies can be challenging due to the presence of confounders such as other features of metabolic syndrome and behaviours including smoking and alcohol intake.
In UK Biobank, we detected a strong correlation between cT1 and BMI (r2 = 0.36, p = 5×10−324) and also between cT1 and MRI-determined liver fat PDFF (r2 = 0.62, p = 5×10−324), and a weak but significant inverse correlation with liver iron (r2 = -0.069, p = 6.6×10−18), which is to be expected since cT1 measures were corrected for liver iron (Fig. S5). We used genetic methods (Mendelian randomisation, Fig. 4) that are generally free of biases such as confounding and reverse causation to examine the potential causal effect of metabolic traits on liver cT1. We found evidence of a causal association between insulin resistance (IVW p = 0.0001), non-alcoholic fatty liver (IVW p = 0.01), type 2 diabetes (IVW p = 0.004), BMI (IVW p = 0.002) and higher cT1. We also found evidence for a protective role of favourable adiposity variants (variants associated with higher adiposity but lower risk of cardiometabolic diseases and lower ectopic fat)
Genome-wide and abdominal MRI data provide evidence that a genetically determined favorable adiposity phenotype is characterized by lower ectopic liver fat and lower risk of type 2 diabetes, heart disease, and hypertension.
We used 2-sample Mendelian randomisation analysis to investigate the causal effects of metabolic traits on liver cT1. For full results, including sensitivity analyses, please see Table S4. The inverse variance weighted test was used as the main analysis. Levels of significance: p <0.05. 2hGlu, 2 hour glucose tolerance test; cT1, corrected T1; NAFLD, non-alcoholic fatty liver disease; WHR_BMI, waist hip ratio adjusted for BMI.
We identified associations between 6 independent genetic variants and MRI-based liver cT1, a non-invasive marker of liver inflammation and fibrosis, in 14,440 participants from the UK Biobank. These include 5 missense variants (in SLC39A8, PCK2, TM6SF2, PNPLA3, and TMPRSS6) and 1 intronic variant (in SLC30A10). The cT1-increasing alleles in 4 genes (SLC39A8, SLC30A10, PNPLA3, and TM6SF2) were also associated with higher AST (n = 360,731) and higher ALT (n = 361,940) in the UK Biobank and also in an independent GWAS of liver enzymes (except for SLC30A10 where data was not available).
In some earlier publications, cT1 was reported using the LIF (liver inflammation and fibrosis) score (Supplementary Material). The LIF score is a tri-linear mapping of cT1 onto a continuous scale from 0 to 4 based on the association of cT1 with histological fibrosis.
The LIF cut-off of 1.4 had a sensitivity of 91% and a specificity of 52% for the diagnosis of NASH vs. steatosis (AUROC = 0.80), and corresponds to a cT1 value of 780 ms; a slightly higher cut-off of 800 ms is used in clinical trials
Reply to: “Multiparametric magnetic resonance imaging to predict clinical outcomes in patients with chronic liver disease: a cautionary note on a promising technique”.
The LIF score is no longer used since the medical and MRI physics community is more familiar with T1 for the assessment of inflammation and fibrosis across all specialties including cardiology and neurology.
Reference range of liver corrected T1 values in a population at low risk for fatty liver disease-a UK Biobank sub-study, with an appendix of interesting cases.
Utility and variability of three non-invasive liver fibrosis imaging modalities to evaluate efficacy of GR-MD-02 in subjects with NASH and bridging fibrosis during a phase-2 randomized clinical trial.
Evaluation of age-related interstitial myocardial fibrosis with cardiac magnetic resonance contrast-enhanced T1 mapping: MESA (Multi-Ethnic Study of Atherosclerosis).
In this GWAS study, the cT1 values reported are standardised across the MRI scanner model and field strength, showing very high repeatability and reproducibility.
The missense variant (rs13107325-C >T) in SLC39A8 is predicted to be deleterious in both Polyphen-2 and SIFT, and is associated with lower expression of SLC39A8 in human liver.
ZIP8 is a divalent cation importer capable of transporting zinc, manganese, iron, cadmium and selinate; the substitution of C for T allele impairs the cellular uptake of metals by this protein.
It is not known which metal is involved in liver pathogenicity but there is evidence that hepatic ZIP8 regulates manganese metabolism in the liver, a metal ion that is hepatotoxic at high levels.
Zinc and selenium also have important roles in liver cellular injury, oxidative stress and dysregulated inflammation; dietary supplementation of both has shown benefit in animal models of liver disease.
The pathogenic role of SLC39A8 in liver inflammation and fibrosis is supported by studies in mice which provide mechanistic evidence for the critical role of ZIP8 in liver disease. Liu et al.
used 2 mouse models to study the function of SLC39A8 in the liver. In the first model, they studied the chronic effect of SLC39A8 knockdown. The SLC39A8(neo/neo) homozygous mice died before or immediately after birth. The SLC39A8(+/neo) heterozygous mice had moderate ZIP8 deficiency which led to disruption of normal hepatocellular architecture, necrosis, inflammation, fibrosis and development of liver tumours with histopathological features consistent with hepatocellular neoplasms.
In the second model, they studied liver-specific SLC39A8 knockdown by adenovirus-delivered short-hairpin RNA and demonstrated that liver damage in the chronic model is not due to some extrahepatic process. Liver-specific ZIP8 downregulation for 7 days resulted in substantial hepatomegaly, inflammation, proliferation, oxidative stress, liver injury and cell death.
The intronic variant in SLC30A10, a gene which codes a predominantly manganese metal ion transporter, was also associated with elevated cT1 measures in our study, as well as elevated aminotransferases in the UK Biobank. Manganese is an essential metal required for the adequate functioning of numerous enzymes, however it is toxic and induces cell death at elevated cellular levels.
Loss-of-function mutations in SLC30A10 have previously been associated with cirrhosis, higher manganese levels in liver biopsy samples and neurotoxicity including parkinsonian-like movement disorders.
The association between cT1 increasing alleles at the 2 novel loci (SLC39A8 and SLC30A10) and higher ALT and AST adds supportive evidence for their pathogenic role in the liver. The missense variant in SLC39A8 has previously been shown to be associated with multiple traits including alcohol intake, BMI, schizophrenia, Crohn's disease, lower brain grey matter volume and microbiome diversity;
we show for the first time a further novel association with higher diabetes and triglyceride levels, whilst highlighting variable effects on cholesterol levels. The associations of both variants with cT1 were independent of BMI, alcohol intake, liver fat percentage and liver iron content in our sensitivity models.
We identified a further 2 missense variants that were associated with cT1 but not with elevated aminotransferases; therefore, further research is required to validate these findings and explore their potential role in liver inflammation and fibrosis. The cT1-increasing allele in rs111723834 (missense variant in PCK2, also an intronic variant in NRL) was associated with lower aminotransferases, lower risk of type 2 diabetes, and lower triglycerides. PCK2 encodes a mitochondrial enzyme that catalyses the conversion of oxaloacetate to phosphoenolpyruvate and has a key role in hepatic gluconeogenesis. Mitochondrial phosphoenolpyruvate carboxykinase deficiency (M-PEPCKD) is a rare autosomal recessive disorder resulting from impaired gluconeogenesis, and clinical characteristics include hypotonia, hepatomegaly, failure to thrive, lactic acidosis and hypoglycaemia.
The missense variant in PCK2 is also an intronic variant in NRL, and it is unclear which gene is associated with elevated cT1 measures. NRL however encodes for neural retinal leucine zipper transcription factor that is specifically expressed in neuronal retina cells, making it an unlikely causal gene candidate for liver cT1. The cT1-increasing allele (rs4820268-A >G) in TMPRSS6 has previously been reported to be associated with lower plasma iron levels and lower liver iron content.
It is also associated with a dysmetabolic profile including higher LDL cholesterol, higher cardiovascular disease risk and hypertension (Fig. 2). Its effect on cT1 however remained significant even after correcting for liver iron content in sensitivity analyses, making it unlikely that the association was secondary to bias resulting from iron correction when calculating cT1. Previous Mendelian randomisation studies have shown that higher circulating iron may be cardioprotective,
The same mechanisms may explain why the allele associated with lower circulating iron levels is associated with higher cT1.
Known NAFLD and cirrhosis risk alleles in PNPLA3 and TM6SF2 were also associated with both elevated cT1 and MRI-derived PDFF in our cohort. These associations provide strong positive controls for our study and validate for the first time the association with MRI-determined liver PDFF. The risk alleles in these 2 genes were further associated with higher risk of type 2 diabetes, but with lower serum triglycerides, LDL cholesterol, and lower risk for cardiovascular disease, as previously described.
In our GWAS on PDFF, alongside PNPLA3 and TM6SF2, we further identified variants in GCKR (another known NAFLD variant which we have replicated) and APOE (apolipoprotein E, a gene which encodes a major cholesterol carrier).
The APOE risk allele (T) for PDFF is associated with a higher risk of diabetes, and lower risk of cardiovascular disease and elevated LDL cholesterol in independent GWASs.
were also associated with liver cT1 with the first 2 reaching genome-wide significance. The paradoxical inverse association between the liver iron-increasing allele in HFE may be due to overcorrection since cT1 measures are corrected for liver iron content and were inversely correlated in our cohort. Adjustment for liver iron content in our sensitivity analysis remarkably attenuated the association with cT1. The SERPINA1 variant was only associated with all-cause cirrhosis in a recessive model.
We did not have any homozygotes in our liver cT1 cohort to detect a recessive model of association with cT1.
Identifying causal mechanisms to steatohepatitis is crucial since interventions targeting these modifiable exposures may prevent liver disease progression. Our Mendelian randomisation study investigated 24 possible metabolic traits that may cause steatohepatitis. We provide genetic evidence that insulin resistance, non-alcoholic fatty liver and higher BMI causally increase cT1. Recent genetic studies have further identified variants associated with higher BMI but lower risk of type 2 diabetes, hypertension and heart disease.
These “favourable adiposity” variants are also associated with higher subcutaneous-to-visceral adipose tissue ratio and may protect from disease through higher adipose storage capacity, by promoting lipid deposition in subcutaneous tissue rather than within the circulation and ectopic places. The inverse link between favourable adiposity and steatohepatitis provides supportive evidence for the protective effects of this phenotype on a variety of cardiometabolic diseases, underlying mechanisms which can be further explored and pointing to future preventive and therapeutic strategies.
Our study had several limitations. We did not have any independent cohort to replicate our findings. To overcome this limitation, we investigated associations between cT1 variants and ALT and AST levels both in the UK Biobank and an independent GWAS of liver enzymes.
While MRI-derived cT1 is clinically available and is used to assess the severity of steatohepatitis, this measure is still novel, and further research is needed to determine the relative contributions of inflammation and fibrosis to cT1.
Whilst it would be useful to have histological reference data for cT1, pathologist-interpreted liver biopsies do not lend themselves to large studies of this nature because of the risk to patients and inter-rater variance in assessment of histology. This may be improved with advances in digitally processed histology and centralised collection of pathological data for large consortia like the European LITMUS study. While cT1 has demonstrated excellent repeatability
Only 2 of the 6 cT1 variants were associated with liver steatosis, which highlights the complementarity of cT1 and liver fat PDFF as biomarkers of liver status, and their potential to recognise different mechanisms predisposing individuals to liver disease.
Conclusion
cT1 and PDFF phenotypes share some but not all aetiopathogenic mechanisms. We identified novel associations between an MRI-derived measure of fibroinflammatory liver disease and variants in SLC30A10 and SLC39A8 that replicated with blood biomarkers of hepatocyte injury. These genes have a critical role in transporting heavy metal cofactors for a multitude of biological processes. Future studies may determine whether targeting SLC30A10 and SLC39A8 are possible therapeutic options to prevent liver disease in at-risk individuals. Our Mendelian randomisation study provides genetic evidence that addressing weight gain and insulin resistance are useful strategies in the prevention of steatohepatitis.
Abbreviations
ALP, alkaline phosphatase; ALT, alanine aminotransferase; AST, aspartate aminotransferase; BMI, body mass index; CAD, coronary artery disease; cT1, corrected T1; EA, effect allele; EAF, effect allele frequency; GGT, gamma-glutamyltransferase; GWAS, genome-wide association study; HDL-C, HDL-cholesterol; IVW, inverse variance weighted; LD, linkage disequilibrium; LDL-C, LDL cholesterol; LIF, liver inflammation and fibrosis; NAFLD, non-alcoholic fatty liver disease; NASH, non-alcoholic steatohepatitis; PC, principal components; PDFF, proton density fat fraction; OA, other allele; OR, odds ratio; ShMOLLI, shortened modified look locker inversion; SNP, single-nucleotide polymorphism; T2DM, type 2 diabetes; WHR_BMI, waist hip ratio adjusted for BMI.
Financial support
H.Y. is funded by a Diabetes UK RD Lawrence fellowship (17/0005594). C.A.P is funded by a Wellcome Trust Clinical PhD Programme (206274/Z/17/Z). H.W is funded by an Innovate UK Knowledge Transfer Partnership (KTP10271). H.H is a National Institute for Health Research Senior Investigator.
Authors' contributions
C.A.P, H.W, H.Y, R.B. were involved in the study conception and design, analysis and writing of the manuscript. J.M and R.C.N were involved in analysis. H.H, A.D.H, R.S.P, C.A.P. E.L.T, J.B. edited the manuscript. S.N, R.B, M.K. edited the manuscript and provided the infrastructure underlying the MRI liver cT1 measurements.
Ethical approval
This research has been conducted using data from the UK Biobank resource and carried out under UK Biobank project application number 9914. UK Biobank protocols were approved by the National Research Ethics Service Committee.
Patient consent
No participants were directly involved in our study, as we used data derived from the UK Biobank study, under project application number 9914. For the UK Biobank overall study, participants signed written informed consent, specifically applicable to health-related research. All ethical regulations were followed. No patients or participants were specifically or directly involved in setting the research question or the outcome measures or in developing plans for recruitment, design, or implementation of this study. No patients were asked to advise on interpretation or drafting of results. There are no specific plans to disseminate the research results to study participants, but the UK Biobank disseminates key findings from projects on its website.
Data availability
Full data including individual cT1 and PDFF measures will be returned to UK Biobank and made publicly available via application (amsportal.ukbiobank.ac.uk).
Conflict of interest
M.K, J.M, R.C.N and R.B. are employees and shareholders of Perspectum Diagnostics. H.W. and S.N. are shareholders in Perspectum Diagnostics.
Prevalence of clinically significant liver disease within the general population, as defined by non-invasive markers of liver fibrosis: a systematic review.
Prevalence of nonalcoholic fatty liver disease and nonalcoholic steatohepatitis among a largely middle-aged population utilizing ultrasound and liver biopsy: a prospective study.
The diagnosis and management of non-alcoholic fatty liver disease: practice guideline by the American Gastroenterological Association, American Association for the Study of Liver Diseases, and American College of Gastroenterology.
Reference range of liver corrected T1 values in a population at low risk for fatty liver disease-a UK Biobank sub-study, with an appendix of interesting cases.
LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis.
Genome-wide and abdominal MRI data provide evidence that a genetically determined favorable adiposity phenotype is characterized by lower ectopic liver fat and lower risk of type 2 diabetes, heart disease, and hypertension.
Reply to: “Multiparametric magnetic resonance imaging to predict clinical outcomes in patients with chronic liver disease: a cautionary note on a promising technique”.
Utility and variability of three non-invasive liver fibrosis imaging modalities to evaluate efficacy of GR-MD-02 in subjects with NASH and bridging fibrosis during a phase-2 randomized clinical trial.
Evaluation of age-related interstitial myocardial fibrosis with cardiac magnetic resonance contrast-enhanced T1 mapping: MESA (Multi-Ethnic Study of Atherosclerosis).
30194-X&id=fx1.jpg){kind=link}
30194-X&id=gr1.jpg){kind=link}
30194-X&id=gr2.jpg){kind=link}
30194-X&id=gr3.jpg){kind=link}
30194-X&id=gr4.jpg){kind=link}