Microbiota, type 2 diabetes and non-alcoholic fatty liver disease: protocol of an observational study
Journal of Translational Medicine volume 17, Article number: 408 (2019)
Non-alcoholic fatty liver disease (NAFLD) is characterized by triglyceride accumulation in the hepatocytes in the absence of alcohol overconsumption, commonly associated with insulin resistance and obesity. Both NAFLD and type 2 diabetes (T2D) are characterized by an altered microbiota composition, however the role of the microbiota in NAFLD and T2D is not well understood. To assess the relationship between alteration in the microbiota and NAFLD while dissecting the role of T2D, we established a nested study on T2D and non-T2D individuals within the Cooperative Health Research In South Tyrol (CHRIS) study, called the CHRIS-NAFLD study. Here, we present the study protocol along with baseline and follow-up characteristics of study participants.
Among the first 4979 CHRIS study participants, 227 individuals with T2D were identified and recalled, along with 227 age- and sex-matched non-T2D individuals. Participants underwent ultrasound and transient elastography examination to evaluate the presence of hepatic steatosis and liver stiffness. Additionally, sampling of saliva and faeces, biochemical measurements and clinical interviews were carried out.
We recruited 173 T2D and 183 non-T2D participants (78% overall response rate). Hepatic steatosis was more common in T2D (63.7%) than non-T2D (36.3%) participants. T2D participants also had higher levels of liver stiffness (median 4.8 kPa, interquartile range (IQR) 3.7, 5.9) than non-T2D participants (median 3.9 kPa, IQR 3.3, 5.1). The non-invasive scoring systems like the NAFLD fibrosis score (NFS) suggests an increased liver fibrosis in T2D (mean − 0.55, standard deviation, SD, 1.30) than non-T2D participants (mean − 1.30, SD, 1.17).
Given the comprehensive biochemical and clinical characterization of study participants, once the bioinformatics classification of the microbiota will be completed, the CHRIS-NAFLD study will become a useful resource to further our understanding of the relationship between microbiota, T2D and NAFLD.
Non-alcoholic fatty liver disease (NAFLD) encompasses a spectrum of disorders characterized by hepatic triglyceride accumulation (hepatic steatosis) in the absence of alcohol overconsumption . Twenty-to-thirty % of the NAFLD patients progress to non-alcoholic steatohepatitis (NASH), implying liver inflammation and association with liver-related diseases such as fibrosis, cirrhosis, and hepatocellular carcinoma . NAFLD has an estimated prevalence of about 25% in the general population  and is expected to become the leading cause of liver transplantation over the next 20 years, with expanding costs for healthcare systems . In the presence of type 2 diabetes (T2D) and metabolic syndrome (MetS), NAFLD prevalence can rise up to 70% [5,6,7]. T2D may increase the risk of developing NAFLD and NASH, but also NAFLD itself may be a risk factor for T2D onset [8, 9]. The interplay between NAFLD, T2D, and MetS is complex, with NAFLD and MetS sharing clinical manifestations such as obesity, insulin resistance, T2D, dyslipidemia, and hypertension .
Triggering factors, such as the translocation of bacterial components and their products from the gut into the systemic circulation following alterations of the intestinal integrity, have been identified as an important mechanism of NAFLD onset [11, 12]. Animal studies suggest that bacterial components may also have a crucial role in NAFLD and NASH onset in humans .
The gut microbiota represents the community of microorganisms inhabiting the digestive tract, with vital functions in relation to vitamin biosynthesis, bile acid degradation, maintenance of the intestinal mucosal barrier integrity, and complex carbohydrate digestion. Microbiota composition can be altered by medications, environmental factors, and diet [14, 15]. Individuals with T2D [16,17,18] or MetS  may suffer a disrupted intestinal microbial composition, which then promotes an imbalance between protective and harmful effects of the microbiota on the host. Increasing evidence suggests an important role for the intestinal microbiota in the pathogenesis of T2D, regulating metabolic pathways and glucose hemostasis . For example, obesity is associated with a larger number of bacterial strains that ferment food components and increase the potential of the host to harvest energy . Obese and MetS-affected individuals may have an altered ratio of Firmicutes-to-Bacteroidetes specific phyla [22, 23]. Members of both phyla are involved in the production of short chain fatty acids (SCFAs) from dietary compounds not completely digested in the small intestine . Evidence suggests a role of SCFAs as key mediators of the cross talk between brain and gut in the pathogenesis of obesity . It has been shown that the total amount of SCFA produced is higher in obese subjects, suggesting that SCFA metabolism might play a considerable role in obesity . Recently, using bidirectional Mendelian randomization, the causal relationship between the gut microbiome and metabolic traits has been explored, providing evidence of a causal effect of the gut microbiome on metabolic traits . Several studies have shown alterations of the microbiota composition in NAFLD [28,29,30] and NASH  patients. Increased relative abundance of Bacteroides and Ruminococcus in the intestine has been associated with more severe histology in NAFLD patients . In contrast, NAFLD patients seem to have lower relative abundance of the Prevotella strain [32, 33]. Recently, 37 bacterial strains from the gut were identified that allowed discrimination between mild and severe hepatic fibrosis in biopsy-proven NAFLD patients .
Another important microbial habitat is the oral cavity . In addition to contributing to oral diseases [36, 37], the oral microbiota may represent a risk factor for systemic diseases such as T2D [38, 39]. NAFLD was associated with periodontitis, and some characteristics of periodontitis such as systemic inflammation and invasion of the commensal bacteria are involved in progression of liver fibrosis in NAFLD-affected individuals . Experimental models showed an association between Porphyromonas gingivalis and risk of NAFLD and NASH [41, 42]. The importance of microbial invasion from the oral cavity into the lower intestine compartments in patients with cirrhotic liver disease was recently demonstrated by a study showing that > 50% of bacterial species found to be enriched in the intestine of cirrhotic patients were of buccal origin . While the oral microbiota has gained much attention only recently, there remain many unanswered questions on the role of bacterial strains on specific pathologies of the liver.
To further improve the general understanding on the relationship between the microbiota, from gut and mouth, and NALFD, and to illuminate such relationships in the context of T2D, we carried out an observational study nested within the Cooperative Health Research in South Tyrol (CHRIS) study , called CHRIS-NAFLD. Here, we describe the study protocol, recruitment and measurement procedures, and offer a description of the epidemiological characteristics of study participants at the baseline and follow-up times.
The CHRIS-NAFLD study was setup in the context of the CHRIS study, a population-based study carried out in a rural Alpine context [44, 45]. CHRIS study participants were on overnight fasting, underwent blood and urine collection, anthropometric and blood pressure (BP) measurements, electrocardiographic analysis, and tremor assessment. Participants were also administered a series of interviewer- and self-administered questionnaires about their health status. The CHRIS-NAFLD study was designed after the recruitment of the first 4979 CHRIS study participants, which was carried out between 2011 and 2014 (baseline data) . From the 4979 participants, we selected for the CHRIS-NAFLD study all 227 individuals affected by T2D and an equal number of non-T2D individuals, matched on age (± 2 years tolerance) and sex. T2D was defined according to standard guidelines , as a positive response to the question “Has a doctor ever diagnosed you with diabetes?” or either having fasting plasma glucose levels of ≥ 126 mg/dl or glycated hemoglobin (HbA1c) levels of ≥ 6.5%. Participants with other types of diabetes were excluded from the selection. Among participants selected as non-T2D diabetes (negative response to the doctor’s diagnosed diabetes question), participants were further excluded in case of HbA1c levels beyond 5.6%, to prevent inclusion of potentially pre-diabetic participants.
Recruitment of study participants
Selected participants were informed about the objective and content of the CHRIS-NAFLD study by invitation letter, mailed between 2 and 3 weeks prior to enrollment. They were subsequently contacted by phone to arrange an appointment at the study center. The participants were requested not to eat anything from 8:00 pm of the night before the day of participation, and to abstain from drinking or smoking for at least 2 h prior to the visit. At the study center, participants underwent blood drawing, urine collection, saliva sampling, anthropometric and BP measurements, tremor assessment, and clinical examination by a medical doctor to assess the hepatic status. Body mass index (BMI), fat percentage, and visceral and subcutaneous fat were assessed using a body composition monitor (OMRON BF508). Waist and hip circumference were measured according to the WHO protocol .
Questionnaires and interview
Questionnaires concerning the change of participants’ health status since the baseline participation and their life-style were administered by an interviewer, and the medication intake in the last 7 days was documented using an electronic optical scan of their medication box bar codes according to the Anatomical Therapeutic Chemical (ATC) classification system. A Food Frequency Questionnaire (FFQ) based on the Global Allergy and Asthma European Network of Excellence study  was mailed to their homes prior to participation, in order to limit the time spent at the study center. The FFQ also asked about the average frequency of the consumption of alcoholic drinks during the last 12 months (rarely or never, 1–3/month, 1/week, 2–4/week, 5–6/week, 1/day, 2+/day), specifically of beer (200 ml), red wine (125 ml), white wine (125 ml), rosé wine (125 ml), liqueurs (50 ml), and spirits (50 ml). At the study center, participants answered a computer-assisted interviewer-administered questionnaire on smoking habits, based on the European Community Respiratory Health Survey II  from which we derived pack-years as a measure of cumulative smoking. Based on their smoking habits, participants were classified as never smokers (never smoked or smoked for < 1 year in their lifetime), past-smokers (smoked for ≥ 1 year in their lifetime but stop smoking ≥ 1 year before the interview), and current-smokers (currently smoking at the time of the interview or stopped smoking < 1 month before the interview).
Evaluation of the hepatic steatosis and fibrosis
To evaluate the presence of hepatic steatosis and fibrosis, participants underwent abdominal ultrasound examination (5-1 MHz Phased Array Transducer, iViz, SonoSite, USA) and transient elastography (Fibroscan®, Echosens, France) performed by a trained medical doctor (Christoph Grander).
To assess transient elastography (TE), participants were placed in a supine position with their right arm fully adducted and asked to hold their breath. At least ten independent resistance measurements were taken, starting always with an M + probe but using an XL + probe as a backup option when prompted by the automatic probe selection tool [50, 51]. TE values were defined as unreliable when the IQR to median ratio was > 30%. Fibrosis was then scaled into four stages, F0 to F4, based on the resistance levels of the liver . TE values of > 6.5 kPa were considered as diagnosis of fibrosis (≥ F1) .
Additionally, abdominal ultrasound scanning was performed in every participant after overnight fasting. Presence and severity of steatosis was evaluated as documented by Ballestri et al. . Steatosis was classified into three grades: normal or very slight increase in the echo pattern with normal visualization of vessels and diaphragm (grade 1); moderate increase in echogenicity with reduced visibility of portal veins and diaphragm (grade 2); or distinct increase in echo pattern with poor visibility of intrahepatic vessels and diaphragm (grade 3). Presence of gallstones, gallbladder size, and wall thickness were also assessed, as well as the visceral and subcutaneous fat thickness to gain insights into the participant’s metabolic risk profile .
The presence of NAFLD was defined as steatosis grade of ≥ 2, after the exclusion of other causes such as overt hepatitis due to virus infection, or hereditary liver disorders or other liver diseases . Three participants were excluded from further analyses because of potential drug induced steatosis by methotrexate. Participants were not excluded based on alcohol consumption levels.
Advanced fibrosis was defined as LSM values ≥ 6.5 kPa
For a deeper characterization of steatosis and fibrosis, additional surrogate markers have been calculated as reported in Box 1.
Biospecimen collection and biobanking
For the CHRIS-NAFLD study, blood (49 ml) and urine (30 ml) samples for laboratory analysis and biobanking were collected in the early morning, after an overnight fasting, following the same procedures previously described for the CHRIS study in terms of sample pre-analytical processing, transportation, and biobanking [44, 45]. In addition, two 1 ml aliquots of serum were stored at − 80 °C and sent in dry ice in a unique batch at the end of the recruitment to Synlab Italia Srl for insulin measurement. The CHRIS biobank was assigned a “Bioresource Research Impact Factor” code BRIF6107 [44, 66].
Stool and saliva collection for the metagenomic analysis
Stool collection tubes were shipped to the participant’s home some days before the enrollment with instructions for the sample collection. Participants were asked to collect the samples in the same morning of their visit or, failing that, within 24 h of the visit. Participants brought their stool samples to the study center at room temperature. Once at the study center, samples were stored at − 20 °C. Samples were then transported frozen to the biobank, where they were finally stored at − 80 °C. Information about the exact time of defecation was collected and the Bristol stool scale (BSS) was assessed . The BSS is used to classify stool consistence: it can be used as a surrogate marker for stool transit time  and is applied in both clinical and experimental fields .
Non-stimulated saliva samples were collected using the Omnigene oral collection device (OM-501, DNA Genotek, USA) at the study center. After collection, stabilized saliva samples were transported to the biobank, where they were stored at − 80 °C after splitting into 500 μl aliquots. Status of teeth and gums was assessed using items 3, 6, and 12 of the WHO’s Oral Health Questionnaire (OHQ) for Adults .
Microbiome extraction protocol
DNA extraction from the faeces was performed using a Chemagic Magnetic Separation Module I Dispenser (Perkin Elmer, USA) according to the Chemagen protocol (Chemagic DNA Feces 1 k drying prefilling H12 VD160617.che) using a blood kit (CMG-763-1, Perkin Elmer) supplemented with a lysis buffer specifically for stools (CMG-852, Perkin Elmer). Briefly, under a sterile hood, up to 1 g of each sample was taken and immediately immersed in 8 ml lysis buffer. The weight was annotated and the tube mixed thoroughly on a vortex. After adding 50 μl protease mix (provided in the kit), the sample was incubated for 20 min at 70 °C followed by 5 min inactivation at 95 °C. The sample was centrifuged and the supernatant transferred to a new tube that was further processed on the robot.
For the DNA extraction protocol from saliva, based on chemical lysis and purification for downstream applications, we followed the manufacturer’s protocol (CMG-1037, Chemagic DNA Saliva Kit special, Perkin Elmer). DNA was quantified with QuantiFluor (E2670, Promega) on an Envision plate reader (Perkin Elmer) and quality was assayed on a NanoDrop spectrophotometer (ThermoFisher, USA) and by running on a 0.5% agarose gel.
Metagenomic sequencing will be based on the amplification of the V3–V4 hypervariable region of the 16S RNA gene with specific primers which can then be used to incorporate unique indexes into the fragments which will further allow unique indexing of up to 384 samples that can then be pooled together. The resulting multiplexed pool will be run on the MiSeq System using the V3 chemistry 600 cycle kit (16S Metagenomic Sequencing Library Preparation System, Illumina, USA).
16S data processing and analysis
We will follow the data processing and the analytical pipeline developed by the MiBioGen consortium , comprising the following steps: 16S data processing, genotype data processing, and genome-wide association study (GWAS) . For the 16S data processing, we will use the Ribosomal Database Project (RDP) Classifier instead of OTU picking, since it leads to more consistent results and for the genotype imputation we will use the freely available Michigan Imputation Server . For the GWAS analysis we will follow the uniform analytical pipeline developed by the consortium.
Statistical analyses and power calculation
The association between microbiota composition and NAFLD will be assessed using the Fisher’s exact test for proportions. We built power scenarios using the “power two-proportions” command with the “test(fisher)” option implemented in Stata version 14. The impact of predictor variables on the presence of NAFLD or NAFLD-related symptoms and markers will be determined using univariable and multivariable logistic regression models. The role of T2D in the microbiome-NAFLD relationship will be investigated by means of interaction analyses. To assess the relationship between NAFLD severity and microbiota we will fit linear regression models. All models will be further controlled for participants’ relatedness to avoid biased estimates due to population structure. Finally, we will analyze the difference between microbiota measured from saliva and from stools in participants with NAFLD and in those with T2D.
The CHRIS-NAFLD study protocol was approved by the Ethical Committee of the Healthcare System of the Autonomous Province of Bolzano (Südtiroler Sanitätsbetrieb/Azienda Sanitaria dell’Alto Adige), protocol no. 85-2016 (19 Oct 2016). As it is nested within the CHRIS study, the CHRIS-NAFLD study follows the CHRIS’s ethics protocols with regard to the collection, use and access of data and biosamples, which are stored for a long-term (30 years) after the end of recruitment. Participants were recruited based on a previous consent that allowed re-call. The CHRIS study uses online dynamic consent procedures for empowering the autonomy and compliance of study participants. The dynamic consent options for re-contact, allowed us to re-invite the participants for the CHRIS-NAFLD study, collecting additional information and re-consent. Prior to participation, participants were informed about the objectives and extra procedures of this additional study, for which they provided written informed consent.
Recruitment took place between October 2016 and February 2017. Out of 454 invited individuals, 356 were recruited (78.4% participation rate). Participation rate was similar in the T2D (173 out of 227 invited participants) and non-T2D (183 out of 227 invited participants) groups. At baseline, the 356 individuals who accepted to participate in the CHRIS-NAFLD study, were on average 67.6 years old (standard deviation, SD = 10.6) and 177 (49.7%) were females. One participant selected in the non-T2D group reported an incident T2D diagnosis at the time of the CHRIS-NAFLD examination and was thus included in the T2D group. We collected saliva and stool samples of 354 and 350 participants, respectively. The DNA extraction and quality control were completed for all samples and they all amplified in PCR. DNA was normalized to ~ 100 ng/µl and stored at − 80 °C until further use.
Characteristics of the 356 participants at the time of participation to the CHRIS-NAFLD study are described in Table 1, following the stratification by T2D status that was used for recruitment. Briefly, 173 (48.6%) and 183 (51.4%) were T2D and non-T2D, respectively. Mean follow-up time was 3.56 years (SD = 0.82) and 3.91 years (0.88) for T2D and non-T2D participants, respectively. As expected, T2D and non-T2D participants had similar age and sex distributions. Of the T2D participants, 8.1% reported to have never consumed alcohol, while 22.5% drink daily (5.5% and 29.5% in the non-T2D group). Most of the participants were never (60.5% and 62.3% for T2D and non-T2D, respectively) or past smokers (32.0% and 32.2%).
Among T2D individuals, 111 (64.2%) reported a diagnosis of diabetes, 88 (51.8%) were on diabetic treatment at the time of the visit, 42.4% had HbA1c values in the pre-diabetic range (HbA1c 6.0 to 6.4%), while 44.2% had HbA1c values in the diabetic range (HbA1c ≥ 6.5%). Fifty percent of the T2D individuals and 18.3% of the non-T2D individuals were obese (BMI ≥ 30). Hypertension was observed at the time of the visit in 55.9% of the T2D individuals and 44.2% of the non-T2D individuals. Use of proton-pump inhibitors (PPIs), which were shown to potentially influence gut microbial composition [73, 74], had similar distribution in T2D and non-T2D participants, while statins were more common in T2D individuals (Table 1).
Hepatic steatosis was diagnosed in 43 of the 173 T2D individuals (24.9%) and 8 of the 183 non-T2D ones (4.4%). Furthermore, T2D participants showed higher values of TE than non-T2D participants (median = 4.8 (IQR = 3.7, 5.9) kPa vs. 3.9 (3.3, 5.1) kPa) (Fig. 1). These findings agreed with non-invasive scoring systems, like NAFLD fibrosis score, suggesting increased liver fibrosis in T2D individuals (Table 1).
According to the criteria for NAFLD calculation, we identified 172 as affected with NAFLD and 181 healthy individuals among the two groups of T2D and not-T2D participants. Power calculations are shown in Fig. 2. Assuming we wish to compare the relative bacterial abundance in participants affected with NAFLD versus healthy ones and analyzing data using a Fisher exact test, we anticipate a > 80% power to detect a difference in the proportion of relative abundance of single bacterial strains in the two groups higher than 15% (Fig. 2, scenario 1). Despite previous studies that did not identify more than 6 phyla, we predict finding 10 phyla (Fig. 2, scenario 2); in this scenario, the study will have 80% power to detect differences in the proportions equal or higher than 20%. In Fig. 2, scenario 3, we consider 250 bacterial subgroups, calculated assuming to find 10 phyla, 5 families per phyla and 5 genera per family. Under this scenario the study would have 80% power to detect differences in proportions higher than 22.5% in most cases. However, in terms of multiple comparisons, this should be considered as a pessimistic scenario as none of previous studies was able to identify phyla, families and genera at such a high level of detail.
CHRIS-NAFLD is a population-based study aimed to assess the relationship between gut and oral microbiota and NAFLD in individuals with and without T2D. Stratification by diabetes will enable assessment of whether the NAFLD-microbiota relationship is modified or moderated by, or is independent of, the T2D status. In addition to the cross-sectional nature of the investigation, the study will help to evaluate whether NAFLD assessed at the time of the CHRIS-NAFLD participation was associated with clinical or biochemical information collected at the time of the CHRIS baseline participation. Finally, the study will provide important information on the relationship between microbiota measured at different sites of the human body, specifically from samples of saliva and stool.
Currently, extensive efforts are being focused on the analyses of microbiota in relation to several diseases. In recent years, the importance of gut microbiota in NAFLD was demonstrated by several independent groups. NAFLD patients exhibit distinct changes in their intestinal flora, which impact the host metabolism . Patients show a higher abundance of bacterial strains which supply the host with nutrient source out of indigestible products, such as complex carbohydrates [21, 76]. However, studies conducted in the general population on the association between microbiota and NAFLD are few and with a smaller number of participants compared to ours [32, 77,78,79].
It is likely that the interaction of genetic and environmental factors with metabolic alterations accelerate NAFLD progression in T2D patients . NAFLD and T2D commonly co-exist and several studies have demonstrated that NAFLD might be found in up to 70% of patients with T2D [6, 81]. From this perspective, our study will allow the question as to whether the NAFLD-microbiota relationship is affected by the diabetic status to be answered.
In a landmark work, Qin and colleagues defined a distinct microbial composition in late stage liver disease. In liver cirrhosis patients, the specific beneficial bacterial strains, like Faecalibacterium prausnitzii, were diminished compared to healthy individuals. Additional analyses of buccal flora showed that, in the patients with liver cirrhosis, harmful bacterial strains are transferred from the oral cavity to the intestine, possibly contributing to the development of cirrhosis . In terms of “oralisation” of the intestinal microbiota, the widespread use of PPIs, which reduce the barrier function of gastric acid , received great attention in the recent years. Several studies analyzed the impact of PPIs on liver diseases such as hepatic encephalopathy  and alcoholic liver disease . Therefore, there is an urgent need to understand the pathophysiological mechanisms leading to NAFLD and to gain more insight into the role of intestinal and oral microbiota in NAFLD.
In the CHRIS-NAFLD study, we collected both saliva and faeces for microbiota analyses. The oral and gastrointestinal microbiome represents the bulk of the overall human microbial load. The correlation of oral microbiota and gut microbiota in NAFLD patients has not been evaluated yet. This will give novel insights into the composition of the microbiota in individuals with and without T2D, possibly identifying microbial transfer in NAFLD patients. Moreover, since saliva is easier to collect compared to stool, if we will observe a similar composition of salivary and stool microbiota, this would increase the compliance of individual participation in such microbiota studies.
In a preliminary descriptive analysis we observed a higher level of liver stiffness in T2D participants, a finding also reported in other cohorts . We also observed a higher prevalence of NAFLD-affected individuals in the T2D groups, as reported by other studies [6, 7].
Strengths of this study include comprehensive assessment of NAFLD, T2D and microbiota in individuals from the general population submitted to a comprehensive evaluation of their hepatic health via ultrasound and elastography examination. In addition to the described data, a wealth of additional genetic, molecular, clinical, environmental data and biological biobanked samples collected in the framework of the CHRIS study are available [44, 45, 85]. Our study has also potential limitations. Even though liver biopsy represents the gold standard for the diagnosis of fibrosis, we used TE for this assessment since the use of an invasive procedure, such as biopsy, would not be ethically justifiable in a population-based observational study. Furthermore, TE is considered a valid noninvasive alternative for this assessment, as reported previously . The ultrasound based method used for detection of steatosis has 85% sensitivity and 94% specificity for identifying a degree of ≥ 20–30% steatosis . Finally, we classified participants that had fasting HbA1c levels of ≥ 6.5% as having T2D since the recent International Expert Committee’s statements recommended using these HbA1c levels as diagnostic criteria for diabetes , but there are also some potential factors that can lead to altered HbA1c levels such as chronic salicylate intake in some individuals . Furthermore, the daily medication of participants was also evaluated in our study. PPIs were equally distributed in T2D and non-T2D individuals. In a large meta-analysis, PPIs were shown to be associated with increased enteric infection with Clostridium difficile . Statins, which were also proven to influence gut microbiota in mice [90, 91], were more commonly reported by T2D individuals. This difference could be explained by higher numbers of dyslipidemia in T2D participants.
In summary, by combining comprehensive bio-sampling with clinical characterization including detailed information on drug history of a large group of individuals with or without T2D and related NAFLD, the CHRIS-NAFLD study will help elucidate important questions on the relationship between microbiota and presence of NAFLD in patients with and without T2D.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
anatomical therapeutic chemical
body mass index
Bristol stool scale
Cooperative Health Research In South Tyrol
Food Frequency Questionnaire
fatty liver index
fasting plasma glucose
genome-wide association study
hepatic steatosis score
homeostatic model assessment-insulin resistance
impaired fasting glucose
NAFLD liver fat score
non-alcoholic fatty liver disease
NAFLD fibrosis score
Oral Health Questionnaire
ribosomal database project
short chain fatty acids
type 2 diabetes
visceral adiposity index
Vernon G, Baranova A, Younossi ZM. Systematic review: the epidemiology and natural history of non-alcoholic fatty liver disease and non-alcoholic steatohepatitis in adults. Aliment Pharmacol Ther. 2011;34(3):274–85.
Singh S, Allen AM, Wang Z, Prokop LJ, Murad MH, Loomba R. Fibrosis progression in nonalcoholic fatty liver vs nonalcoholic steatohepatitis: a systematic review and meta-analysis of paired-biopsy studies. Clin Gastroenterol Hepatol. 2015;13(4):643–54.e1–9 (quiz e39–40).
Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liver disease-meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. 2016;64(1):73–84.
Charlton MR, Burns JM, Pedersen RA, Watt KD, Heimbach JK, Dierkhising RA. Frequency and outcomes of liver transplantation for nonalcoholic steatohepatitis in the United States. Gastroenterology. 2011;141(4):1249–53.
Bellentani S, Saccoccio G, Masutti F, Croce LS, Brandi G, Sasso F, et al. Prevalence of and risk factors for hepatic steatosis in Northern Italy. Ann Intern Med. 2000;132(2):112–7.
Targher G, Bertolini L, Padovani R, Rodella S, Tessari R, Zenari L, et al. Prevalence of nonalcoholic fatty liver disease and its association with cardiovascular disease among type 2 diabetic patients. Diabetes Care. 2007;30(5):1212–8.
Leite NC, Salles GF, Araujo AL, Villela-Nogueira CA, Cardoso CR. Prevalence and associated factors of non-alcoholic fatty liver disease in patients with type-2 diabetes mellitus. Liver Int. 2009;29(1):113–9.
Adams LA, Waters OR, Knuiman MW, Elliott RR, Olynyk JK. NAFLD as a risk factor for the development of diabetes and the metabolic syndrome: an eleven-year follow-up study. Am J Gastroenterol. 2009;104(4):861–7.
Zelber-Sagi S, Lotan R, Shibolet O, Webb M, Buch A, Nitzan-Kaluski D, et al. Non-alcoholic fatty liver disease independently predicts prediabetes during a 7-year prospective follow-up. Liver Int. 2013;33(9):1406–12.
Neuschwander-Tetri BA, Caldwell SH. Nonalcoholic steatohepatitis: summary of an AASLD Single Topic Conference. Hepatology. 2003;37(5):1202–19.
Alisi A, Ceccarelli S, Panera N, Nobili V. Causative role of gut microbiota in non-alcoholic fatty liver disease pathogenesis. Front Cell Infect Microbiol. 2012;2:132.
De Gottardi A, McCoy KD. Evaluation of the gut barrier to intestinal bacteria in non-alcoholic fatty liver disease. J Hepatol. 2011;55(6):1181–3.
Csak T, Velayudham A, Hritz I, Petrasek J, Levin I, Lippai D, et al. Deficiency in myeloid differentiation factor-2 and toll-like receptor 4 expression attenuates nonalcoholic steatohepatitis and fibrosis in mice. Am J Physiol Gastrointest Liver Physiol. 2011;300(3):G433–41.
Wang BH, Yao MF, Lv LX, Ling ZX, Li LJ. The human microbiota in health and disease. Eng Prc. 2017;3(1):71–82.
Langdon A, Crook N, Dantas G. The effects of antibiotics on the microbiome throughout development and alternative approaches for therapeutic modulation. Genome Med. 2016;8(1):39.
Karlsson F, Tremaroli V, Nielsen J, Backhed F. Assessing the human gut microbiota in metabolic diseases. Diabetes. 2013;62(10):3341–9.
Khan MT, Nieuwdorp M, Backhed F. Microbial modulation of insulin sensitivity. Cell Metab. 2014;20(5):753–60.
Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60.
Le Chatelier E, Nielsen T, Qin J, Prifti E, Hildebrand F, Falony G, et al. Richness of human gut microbiome correlates with metabolic markers. Nature. 2013;500(7464):541–6.
Muscogiuri G, Balercia G, Barrea L, Cignarelli A, Giorgino F, Holst JJ, et al. Gut: a key player in the pathogenesis of type 2 diabetes? Crit Rev Food Sci Nutr. 2018;58(8):1294–309.
Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–31.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.
Million M, Maraninchi M, Henry M, Armougom F, Richet H, Carrieri P, et al. Obesity-associated gut microbiota is enriched in Lactobacillus reuteri and depleted in Bifidobacterium animalis and Methanobrevibacter smithii. Int J Obes (Lond). 2012;36(6):817–25.
Louis P, Scott KP, Duncan SH, Flint HJ. Understanding the effects of diet on bacterial metabolism in the large intestine. J Appl Microbiol. 2007;102(5):1197–208.
Barrea L, Muscogiuri G, Annunziata G, Laudisio D, Pugliese G, Salzano C, et al. From gut microbiota dysfunction to obesity: could short-chain fatty acids stop this dangerous course? Hormones (Athens). 2019;18:245–50.
Schwiertz A, Taras D, Schafer K, Beijer S, Bos NA, Donus C, et al. Microbiota and SCFA in lean and overweight healthy subjects. Obesity (Silver Spring). 2010;18(1):190–5.
Sanna S, van Zuydam NR, Mahajan A, Kurilshikov A, Vich Vila A, Vosa U, et al. Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet. 2019;51(4):600–5.
Raman M, Ahmed I, Gillevet PM, Probert CS, Ratcliffe NM, Smith S, et al. Fecal microbiome and volatile organic compound metabolome in obese humans with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. 2013;11(7):868–75.e1–3.
Mouzaki M, Comelli EM, Arendt BM, Bonengel J, Fung SK, Fischer SE, et al. Intestinal microbiota in patients with nonalcoholic fatty liver disease. Hepatology. 2013;58(1):120–7.
Zhang X, Shen D, Fang Z, Jie Z, Qiu X, Zhang C, et al. Human gut microbiota changes reveal the progression of glucose intolerance. PLoS ONE. 2013;8(8):e71108.
Wong VW, Tse CH, Lam TT, Wong GL, Chim AM, Chu WC, et al. Molecular characterization of the fecal microbiota in patients with nonalcoholic steatohepatitis—a longitudinal study. PLoS ONE. 2013;8(4):e62885.
Boursier J, Mueller O, Barret M, Machado M, Fizanne L, Araujo-Perez F, et al. The severity of nonalcoholic fatty liver disease is associated with gut dysbiosis and shift in the metabolic function of the gut microbiota. Hepatology. 2016;63(3):764–75.
Del Chierico F, Nobili V, Vernocchi P, Russo A, Stefanis C, Gnani D, et al. Gut microbiota profiling of pediatric nonalcoholic fatty liver disease and obese patients unveiled by an integrated meta-omics-based approach. Hepatology. 2017;65(2):451–64.
Loomba R, Seguritan V, Li W, Long T, Klitgord N, Bhatt A, et al. Gut microbiome-based metagenomic signature for non-invasive detection of advanced fibrosis in human nonalcoholic fatty liver disease. Cell Metab. 2017;25(5):1054–62.e5.
Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner AC, Yu WH, et al. The human oral microbiome. J Bacteriol. 2010;192(19):5002–17.
Belda-Ferre P, Alcaraz LD, Cabrera-Rubio R, Romero H, Simon-Soro A, Pignatelli M, et al. The oral metagenome in health and disease. ISME J. 2012;6(1):46–56.
Huang S, Yang F, Zeng X, Chen J, Li R, Wen T, et al. Preliminary characterization of the oral microbiota of Chinese adults with and without gingivitis. BMC Oral Health. 2011;11:33.
Casarin RC, Barbagallo A, Meulman T, Santos VR, Sallum EA, Nociti FH, et al. Subgingival biodiversity in subjects with uncontrolled type-2 diabetes and chronic periodontitis. J Periodontal Res. 2013;48(1):30–6.
Hintao J, Teanpaisan R, Chongsuvivatwong V, Ratarasan C, Dahlen G. The microbiological profiles of saliva, supragingival and subgingival plaque and dental caries in adults with and without type 2 diabetes mellitus. Oral Microbiol Immunol. 2007;22(3):175–81.
Alazawi W, Bernabe E, Tai D, Janicki T, Kemos P, Samsuddin S, et al. Periodontitis is associated with significant hepatic fibrosis in patients with non-alcoholic fatty liver disease. PLoS ONE. 2017;12(12):e0185902.
Yoneda M, Naka S, Nakano K, Wada K, Endo H, Mawatari H, et al. Involvement of a periodontal pathogen, Porphyromonas gingivalis on the pathogenesis of non-alcoholic fatty liver disease. BMC Gastroenterol. 2012;12:16.
Acharya C, Sahingur SE, Bajaj JS. Microbiota, cirrhosis, and the emerging oral-gut-liver axis. JCI Insight. 2017;2(19):e94416.
Qin N, Yang F, Li A, Prifti E, Chen Y, Shao L, et al. Alterations of the human gut microbiome in liver cirrhosis. Nature. 2014;513(7516):59–64.
Pattaro C, Gogele M, Mascalzoni D, Melotti R, Schwienbacher C, De Grandi A, et al. The Cooperative Health Research in South Tyrol (CHRIS) study: rationale, objectives, and preliminary results. J Transl Med. 2015;13:348.
Noce D, Gogele M, Schwienbacher C, Caprioli G, De Grandi A, Foco L, et al. Sequential recruitment of study participants may inflate genetic heritability estimates. Hum Genet. 2017;136(6):743–57.
Standards of Medical Care in Diabetes. American Diabetes Association. Diabetes care. 2014;37(Suppl 1):S14–80. https://doi.org/10.2337/dc14-S014.
Waist circumference and waist–hip ratio. Report of a WHO expert consultation, Geneva, 8–11 December 2008.
Garcia-Larsen V, Luczynska M, Kowalski ML, Voutilainen H, Ahlstrom M, Haahtela T, et al. Use of a common food frequency questionnaire (FFQ) to assess dietary patterns and their relation to allergy and asthma in Europe: pilot study of the GA2LEN FFQ. Eur J Clin Nutr. 2011;65(6):750–6.
Chinn S, Jarvis D, Melotti R, Luczynska C, Ackermann-Liebrich U, Anto JM, et al. Smoking cessation, lung function, and weight gain: a follow-up study. Lancet. 2005;365(9471):1629–35 (discussion 00-1).
Castera L. Assessing liver fibrosis. Expert Rev Gastroenterol Hepatol. 2008;2(4):541–52.
Vuppalanchi R, Siddiqui MS, Van Natta ML, Hallinan E, Brandman D, Kowdley K, et al. Performance characteristics of vibration-controlled transient elastography for evaluation of nonalcoholic fatty liver disease. Hepatology. 2018;67(1):134–44.
Wong VW, Vergniol J, Wong GL, Foucher J, Chan HL, Le Bail B, et al. Diagnosis of fibrosis and cirrhosis using liver stiffness measurement in nonalcoholic fatty liver disease. Hepatology. 2010;51(2):454–62.
Siddiqui MS, Vuppalanchi R, Van Natta ML, Hallinan E, Kowdley KV, Abdelmalek M, et al. Vibration-controlled transient elastography to assess fibrosis and steatosis in patients with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. 2019;17(1):156–63 e2.
Ballestri S, Lonardo A, Romagnoli D, Carulli L, Losi L, Day CP, et al. Ultrasonographic fatty liver indicator, a novel score which rules out NASH and is correlated with metabolic parameters in NAFLD. Liver Int. 2012;32(8):1242–52.
Sturm W, Sandhofer A, Engl J, Laimer M, Molnar C, Kaser S, et al. Influence of visceral obesity and liver fat on vascular structure and function in obese subjects. Obesity (Silver Spring). 2009;17(9):1783–8.
Abd El-Kader SM, El-Den Ashmawy EM. Non-alcoholic fatty liver disease: the diagnosis and management. World J Hepatol. 2015;7(6):846–58.
IDF Consensus Worldwide Definition of the Metabolic Syndrome.
Zimmet P, Magliano D, Matsuzawa Y, Alberti G, Shaw J. The metabolic syndrome: a global public health problem and a new definition. J Atheroscler Thromb. 2005;12(6):295–300.
Amato MC, Giordano C, Galia M, Criscimanna A, Vitabile S, Midiri M, et al. Visceral Adiposity Index: a reliable indicator of visceral fat function associated with cardiometabolic risk. Diabetes Care. 2010;33(4):920–2.
Kotronen A, Peltonen M, Hakkarainen A, Sevastianova K, Bergholm R, Johansson LM, et al. Prediction of non-alcoholic fatty liver disease and liver fat using metabolic and genetic factors. Gastroenterology. 2009;137(3):865–72.
Lee JH, Kim D, Kim HJ, Lee CH, Yang JI, Kim W, et al. Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. 2010;42(7):503–8. https://doi.org/10.1016/j.dld.2009.08.002.
Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, et al. The Fatty Liver Index: a simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006;6:33.
Angulo P, Hui JM, Marchesini G, Bugianesi E, George J, Farrell GC, et al. The NAFLD fibrosis score: a noninvasive system that identifies liver fibrosis in patients with NAFLD. Hepatology. 2007;45:846–54. https://doi.org/10.1002/hep.21496
Sterling RK, Lissen E, Clumeck N, Sola R, Correa MC, Montaner J, et al. Development of a simple noninvasive index to predict significant fibrosis in patients with HIV/HCV coinfection. Hepatology. 2006;43(6):1317–25.
Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28(7):412–9.
Mabile L, Dalgleish R, Thorisson GA, Deschenes M, Hewitt R, Carpenter J, et al. Quantifying the use of bioresources for promoting their sharing in scientific research. Gigascience. 2013;2(1):7.
O’Donnell LJ, Virjee J, Heaton KW. Detection of pseudodiarrhoea by simple clinical assessment of intestinal transit rate. BMJ. 1990;300(6722):439–40.
Chumpitazi BP, Self MM, Czyzewski DI, Cejka S, Swank PR, Shulman RJ. Bristol Stool Form Scale reliability and agreement decreases when determining Rome III stool form designations. Neurogastroenterol Motil. 2016;28(3):443–8.
Riegler G, Esposito I. Bristol scale stool form. A still valid help in medical practice and clinical research. Tech Coloproctol. 2001;5(3):163–4.
World Health Organization (WHO). Oral health surveys: basic methods, 5th edition. WHO
Wang J, Kurilshikov A, Radjabzadeh D, Turpin W, Croitoru K, Bonder MJ, et al. Meta-analysis of human genome-microbiome association studies: the MiBioGen consortium initiative. Microbiome. 2018;6(1):101.
Michigan Imputation Server. https://imputationserver.sph.umich.edu/index.html#!
Imhann F, Bonder MJ, Vich Vila A, Fu J, Mujagic Z, Vork L, et al. Proton pump inhibitors affect the gut microbiome. Gut. 2016;65(5):740–8.
Jackson MA, Goodrich JK, Maxan ME, Freedberg DE, Abrams JA, Poole AC, et al. Proton pump inhibitors alter the composition of the gut microbiota. Gut. 2016;65(5):749–56.
Abdul-Hai A, Abdallah A, Malnick SD. Influence of gut bacteria on development and progression of non-alcoholic fatty liver disease. World J Hepatol. 2015;7(12):1679–84.
Samuel BS, Gordon JI. A humanized gnotobiotic mouse model of host-archaeal-bacterial mutualism. Proc Natl Acad Sci USA. 2006;103(26):10011–6.
Lelouvier B, Servant F, Paisse S, Brunet AC, Benyahya S, Serino M, et al. Changes in blood microbiota profiles associated with liver fibrosis in obese patients: a pilot analysis. Hepatology. 2016;64(6):2015–27.
Spencer MD, Hamp TJ, Reid RW, Fischer LM, Zeisel SH, Fodor AA. Association between composition of the human gastrointestinal microbiome and development of fatty liver with choline deficiency. Gastroenterology. 2011;140(3):976–86.
Yun Y, Kim HN, Lee EJ, Ryu S, Chang Y, Shin H, et al. Fecal and blood microbiota profiles and presence of nonalcoholic fatty liver disease in obese versus lean subjects. PLoS ONE. 2019;14(3):e0213692.
Leite NC, Villela-Nogueira CA, Cardoso CR, Salles GF. Non-alcoholic fatty liver disease and diabetes: from physiopathological interplay to diagnosis and treatment. World J Gastroenterol. 2014;20(26):8377–92.
Williamson RM, Price JF, Glancy S, Perry E, Nee LD, Hayes PC, et al. Prevalence of and risk factors for hepatic steatosis and nonalcoholic Fatty liver disease in people with type 2 diabetes: the Edinburgh Type 2 Diabetes Study. Diabetes Care. 2011;34(5):1139–44.
Bian J, Wang A, Lin J, Wu L, Huang H, Wang S, et al. Association between proton pump inhibitors and hepatic encephalopathy: a meta-analysis. Medicine (Baltimore). 2017;96(17):e6723.
Llorente C, Jepsen P, Inamine T, et al. Gastric acid suppression promotes alcoholic liver disease by inducing overgrowth of intestinal Enterococcus. Nat Commun. 2017;8(1):837. https://doi.org/10.1038/s41467-017-00796-x
Koehler EM, Plompen EP, Schouten JN, Hansen BE, Darwish Murad S, Taimr P, et al. Presence of diabetes mellitus and steatosis is associated with liver stiffness in a general population: the Rotterdam study. Hepatology. 2016;63(1):138–47.
Melotti R, Ruscheweyh R, Pramstaller PP, Hicks AA, Pattaro C. Structural Consistency of the Pain Sensitivity Questionnaire in the Cooperative Health Research In South Tyrol (CHRIS) Population-Based Study. J Pain. 2018;19(12):1424–34.
Dhyani M, Anvari A, Samir AE. Ultrasound elastography: liver. Abdom Imaging. 2015;40(4):698–708.
Hernaez R, Lazo M, Bonekamp S, Kamel I, Brancati FL, Guallar E, et al. Diagnostic accuracy and reliability of ultrasonography for the detection of fatty liver: a meta-analysis. Hepatology. 2011;54(3):1082–90.
Radin MS. Pitfalls in hemoglobin A1c measurement: when results may be misleading. J Gen Intern Med. 2014;29(2):388–94.
Janarthanan S, Ditah I, Adler DG, Ehrinpreis MN. Clostridium difficile-associated diarrhea and proton pump inhibitor therapy: a meta-analysis. Am J Gastroenterol. 2012;107(7):1001–10.
Caparros-Martin JA, Lareu RR, Ramsay JP, Peplies J, Reen FJ, Headlam HA, et al. Statin therapy causes gut dysbiosis in mice through a PXR-dependent mechanism. Microbiome. 2017;5(1):95.
Khan TJ, Ahmed YM, Zamzami MA, Mohamed SA, Khan I, Baothman OAS, et al. Effect of atorvastatin on the gut microbiota of high fat diet-induced hypercholesterolemic rats. Sci Rep. 2018;8(1):662.
We thank all study participants, the study personnel, and Dr. Helmuth Weiss (Hospital of Silandro) for his collaboration and supervision; the Biobank’s laboratory technicians and Dr. Claudia Volpato of the Institute for Biomedicine for her collaboration in libraries’ preparation; the IT support and data manager of the Institute for Biomedicine; Dr. Andrew Hicks of the Institute for Biomedicine for the linguistic revision. We thank the Labormed Srl – GmbH, Bolzano (represented by Dr. Monica Cembran) for funding the insulin measurements, which were performed by Synlab Italia Srl, Castenedolo. The authors thank the Department of Innovation, research and University of autonomous Province of Bozen/Bolzano for covering the Open Access publication costs.
This study was supported by the Department of Innovation and Research and University of the Autonomous Province of Bolzano.
Ethics approval and consent to participate
The CHRIS-NAFLD study protocol was approved by the Ethical Committee of the Healthcare System of the Autonomous Province of Bolzano (Südtiroler Sanitätsbetrieb/Azienda Sanitaria dell’Alto Adige), protocol no. 85-2016 (19 Oct 2016). Participants were informed on the objectives of this study and provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Motta, B.M., Grander, C., Gögele, M. et al. Microbiota, type 2 diabetes and non-alcoholic fatty liver disease: protocol of an observational study. J Transl Med 17, 408 (2019). https://doi.org/10.1186/s12967-019-02130-z