Nonalcoholic Fatty Liver Disease (NAFLD) is the most common liver disease worldwide and closely associated with insulin resistance and metabolic syndrome1. The prevalence of NAFLD in the global population is about 25% and the incidence is increasing rapidly in parallel with the westernized diet, lack of exercise, lifestyle changes, and obesity2. It encompasses a spectrum of liver disease that ranges from simple accumulation of fat in liver cells to necroinflammation, liver fibrosis and cirrhosis, with the risk of hepatocellular carcinoma3. Nonalcoholic steatohepatitis (NASH) is considered to be a more severe part in the NAFLD spectrum, which can promote the progression of liver fibrosis faster (7 years per fibrosis stage) than without it (14 years)4; and liver fibrosis5, 6 is considered to be independently associated with long-term outcomes with a high risk of death from liver-related events and all-cause mortality (e.g. cardiovascular disease), and increased with the increased stage of fibrosis7.
In the United States, NAFLD has become the second major indication for liver transplantation and the third major cause of hepatocellular carcinoma, and is still growing 8, 9 . As the clinical consequences of NAFLD increase, there can be severe economic burdens. The annual burden is estimated to be $103 billion ($1,613 per patient) in the United States and €35 billion (€354 to €1,163 per patient) in European countries 10 . Therefore, it is extremely important to fully understand the risk factors and pathogenesis of NAFLD and disease progression. The "multiple strike" hypothesis which include insulin resistance, hormones secreted from the adipose tissue, nutritional factors, gut microbiota and genetic and epigenetic factors largely explains the pathogenesis and progression of NAFLD, but knowledge on the mechanisms of NAFLD still remains incomplete 11 .
Given thyroid hormones play an important role in regulating body metabolism, there has been a heated discussion about the relationship between thyroid dysfunction and NAFLD. Several studies had demonstrated that hyperthyroidism (both subclinical and clinical) is a risk factor contributing to the development of steatosis, liver biopsy-proven NASH and advanced fibrosis12–14. A reasonable explanation is that the thyroid hormones is associated with body fat distribution, metabolic syndrome, and insulin resistance15, 16. Recently, some studies have shown that variation in thyroid hormones in the reference range may also have negative effects on health, same to subclinical and clinical hypothyroidism17. Meanwhile, the relationship between normal thyroid function and NAFLD has also gained attention. However, there is no consensus. In the previous findings, based on different sample sizes and diagnosis methods of NAFLD, they showed that low-normal thyroid function (higher plasma TSH level [2.5 to 4.5 mIU/L] with a normal thyroid hormone level)18, low serum free thyroxine level(fT4)19, high serum free triiodothyronine level(fT3)20, 21, and high fT3/fT4 ratio22 were risk factors for NAFLD respectively.
For this reason, we performed a study in euthyroid subjects with NAFLD aimed at exploring the independent relationship between thyroid function parameters (i.e., triiodothyronine (T3), thyroxine (T4), fT3, fT4, fT3/fT4 ratio and thyroid-stimulating hormone (TSH) levels) and biopsy-proven NASH, significant fibrosis respectively after adjusting other well-identified risk factors (e.g. metabolic risk factors and insulin resistance).
Patients and study design
No informed consent was required because all the data were anonymized. The protocol was in accordance with the Helsinki Declaration and was approved by the ethics committee of the First Affiliated Hospital of Wenzhou Medical University. From our liver biopsy database at The First Affiliated Hospital of Wenzhou Medical University, the clinical records of 1371 subjects during January 2016 and July 2019 were reviewed. Strict exclusion criteria are designed and implemented: (1) ≤ 18-year-old, (2) excessive alcohol consumption (> 140 g/week for men and 70 g/week for women) evaluated by a questionnaire, (3) the steatosis (≤ 5% of liver cells) at histology and history of viral hepatitis, autoimmune hepatitis, or other forms of chronic liver disease, (4) history of malignancy, (5) history of thyroid disease, including clinical hyperthyroidism and hypothyroidism, thyroidectomy, radiofrequency ablation of thyroid gland, (6) with thyroid dysfunction defined as serum TSH > 4.5 mIU/l or < 0.5 mIU/l and/ or fT4 > 14.41 pmol/l or < 7.86 pmol/l, (7) insufficient clinical date. In the end, 307 strictly screened subjects were enrolled in this study. Data was collected from the time of liver biopsy. All variables in this study were objective results stored in the hospital computer system. The presence of diabetes mellitus (fasting blood glucose > = 7.0 mmol/L or treatment with antidiabetic drugs) was recorded. Hypertension was defined as systolic blood pressure ≥ 140 mmHg/diastolic blood pressure ≥ 90 mmHg and/or the current use of anti-hypertensive medication. Smokers were defined as those who had smoked at least one cigarette per day during the previous year.
Anthropometric and laboratory measurements
The body height and weight of the subjects were measured, while they were barefoot and wearing light clothing. The height was recorded to the nearest 1 cm, and body weight was measured to the nearest 0.1 kg. Body mass index (BMI) was calculated by dividing the weight in kilograms with the square of height in meters. Venous blood sampling was collected after overnight fasting for at least 8–12 hours and measured at the hospital Clinical Sample Test Room. Laboratory assays included albumin (g/L), total cholesterol (TC)(mmol/L), triglyceride (TG)(mmol/L), low-density lipoprotein cholesterol (LDL-c)(mmol/L), high-density lipoprotein cholesterol (HDL-c)(mmol/L), aspartate aminotransferase (ALT)(U/L), alanine aminotransferase (AST)(U/L), gamma-glutamyl transferase (γGT)(U/L), alkaline phosphatase (ALP)(U/L), type IV collagen(ng/ml), type III procollagen(ng/ml), hyaluronic acid(ng/ml), total bilirubin (TB)(µmol/L), thyroid stimulating hormone (TSH)(mIU/L), thyroxine (T4)(nmol/L), triiodothyronine (T3)(nmol/L), free thyroxine (fT4)(nmol/L), free triiodothyronine (fT3)(pmol/L), glucose(mmol/L), insulin(pmol/L), glycosylated hemoglobin (GHb)(%), platelet(*10^9), uric acid(µmol/L), creatinine(µmol/L). The fT3/fT4 ratio is the value of fT3 divided by fT4. Insulin resistance was evaluated using the homeostasis model assessment of insulin resistance (HOMA-IR): fasting blood glucose (mmol/l) * insulin (mU/l) / 22.523.
Liver biopsy
Liver biopsy was performed by senior operators using 16-gauge Hepafix needles under ultrasonography positioning. The liver specimens were fixed in 10% formalin, and was scored by experienced hepatologists who were blinded to the clinical data, treatment allocation, and imaging findings. A scoring system published by Kleiner et al24 was used. The histological NAFLD Active Score (NAS) is defined as the unweighted sum of the scores for steatosis (0–3), lobular inflammation (0–3), and hepatocellular ballooning (0–2); scores therefore ranged from 0 to 8. Subjects with scores of 5 or greater were diagnosed as NASH. Fibrosis was staged as follows: stage 0 = no fibrosis; stage 1 = perisinusoidal or periportal fibrosis with 3 different patterns: 1a = mild, zone 3, perisinusoidal; 1b = moderate, zone 3, perisinusoidal fibrosis, and 1c = portal/periportal fibrosis; stage 2 = perisinusoidal and portal/periportal fibrosis; stage 3 = bridging fibrosis; stage 4 = cirrhosis. In this study, we pooled the subtype 1a, 1b, 1c of fibrosis into a single F1 score. Significant fibrosis (SF) was defined as stage 2 or greater (≥ 2).
Statistical Analysis
Continuous data were presented as median (1st quartile, 3rd quartile), and categorical variables were expressed in frequency or as a percentage. First, the univariate analysis (student t-test, Mann-Whitney U test, chi-square test) were used to infer the difference between the 2 groups. It was noted that the fT3/fT4 ratio was too small; therefore, we expanded it 10 times and labelled per 0.1 change (henceforth fT3/fT4 ratio (per 0.1 change)). Second, significant variables from the univariate analysis (P < 0.05) were then subjected to stepwise logistic regression analysis (Probability for Stepwise enter: 0.05, remove: 0.1) to evaluate the risk factors for NASH, significant fibrosis respectively. Third, according to the recommendation of the STROBE statement25, we simultaneously showed the results from unadjusted, minimally adjusted analyses and those from fully adjusted analyses. Fourth, Spearman’s correlation analysis was performed to assess the relationship between thyroid function parameter and histological features. Fifth, subgroup analyses were performed using stratified linear regression models. The modifications and interactions of subgroups were inspected by likelihood ration tests. All of the analyses were performed with the statistical software packages R (http://www.R-project.org, The R Foundation) and EmpowerStats (http://www.empowerstats.com, X&Y Solutions, Inc., Boston, MA). P values less than 0.05 (two-sided) were considered statistically significant.