Gut microbiota profile in patients with nonalcoholic fatty liver disease and presumed nonalcoholic steatohepatitis
Zahra Mohammadi1, Hossein Poustchi1, Azita Hekmatdoost2, Arash Etemadi3, Sareh Eghtesad1, Maryam Sharafkhah1, Delisha Stewart4, Reza Ghanbari5, George Edward Chlipala6, Faraz Bishehsari7, Shahin Merat8, Reza Malekzadeh5
1 Liver and Pancreatobiliary Diseases Research Center, Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran
2 Department of Clinical Nutrition, School of Nutritional Sciences and Dietetics, Shahid Beheshti University of Medical Sciences, Tehran, Iran
3 Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, Maryland, USA
4 Nutrition Research Institute, University of North Carolina at Chapel Hill, Kannapolis, USA
5 Digestive Oncology Research Center, Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran
6 Research Informatics Core, Research Resources Center, University of Illinois at Chicago, Chicago, Illinois, USA
7 Department of Internal Medicine, Division of Gastroenterology, Rush University Medical Center, Chicago, Illinois, USA
8 Digestive Disease Research Center, Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran
|Date of Submission||02-Aug-2021|
|Date of Decision||06-Jan-2022|
|Date of Acceptance||31-Jan-2022|
|Date of Web Publication||29-Jul-2022|
Dr. Faraz Bishehsari
Department of Internal Medicine, Division of Gastroenterology, Rush University Medical Center, Chicago, Illinois
Digestive Diseases Research Institute, Shariati Hospital, North Kargar Street, Tehran 1411713135
Dr. Shahin Merat
Digestive Disease Research Center, Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran
Source of Support: None, Conflict of Interest: None
Background: The main composition of intestinal microbiota in nonalcoholic fatty liver disease (NAFLD) and nonalcoholic steatohepatitis (NASH) patients has not yet been elucidated. In this, case-control study, we identified differences of intestinal microbiota in male patients with NAFLD, presumed NASH, and healthy controls. Materials and Methods: We compared gut microbial composition of 25 patients with NAFLD, 13 patients with presumed NASH, and 12 healthy controls. Demographic information as well as clinical, nutritional, and physical activity data was gathered. Stool and blood samples were collected to perform the laboratory analysis. The taxonomic composition of gut microbiota was assessed using V4 regions of microbial small subunit ribosomal Ribonucleic acid genes sequencing of stool samples. Results: Firmicutes, Actinobacteria, and Bacteroidetes were the most frequently phyla in all groups. Our results revealed that Veillonella was the only genus with significantly different amounts in presumed NASH patients compared with patients with NAFLD (P = 2.76 × 10−6, q = 2.07 × 10−4, logFC = 5.52). Conclusion: This pilot study was the first study to compare gut microbial composition in patients with NAFLD and presumed NASH in the Middle East. Given the potential effects of gut microbiota on the management and prevention of NAFLD, larger, prospective studies are recommended to confirm this study's findings.
Keywords: Gut microbiota, nonalcoholic fatty liver disease, presumed nonalcoholic steatohepatitis
|How to cite this article:|
Mohammadi Z, Poustchi H, Hekmatdoost A, Etemadi A, Eghtesad S, Sharafkhah M, Stewart D, Ghanbari R, Chlipala GE, Bishehsari F, Merat S, Malekzadeh R. Gut microbiota profile in patients with nonalcoholic fatty liver disease and presumed nonalcoholic steatohepatitis. J Res Med Sci 2022;27:54
|How to cite this URL:|
Mohammadi Z, Poustchi H, Hekmatdoost A, Etemadi A, Eghtesad S, Sharafkhah M, Stewart D, Ghanbari R, Chlipala GE, Bishehsari F, Merat S, Malekzadeh R. Gut microbiota profile in patients with nonalcoholic fatty liver disease and presumed nonalcoholic steatohepatitis. J Res Med Sci [serial online] 2022 [cited 2023 Mar 25];27:54. Available from: https://www.jmsjournal.net/text.asp?2022/27/1/54/352944
| Introduction|| |
Nonalcoholic fatty liver disease (NAFLD) has been recently recognized as the most prevalent liver disease worldwide affecting over 25% of the population. This disease is characterized by fat deposition in the liver cells and can progress to nonalcoholic steatohepatitis (NASH), liver fibrosis, cirrhosis, and even hepatocellular carcinoma.
The prevalence of NAFLD is rising parallel to other metabolic morbidities such as obesity, insulin resistance, metabolic syndrome, and dyslipidemia. Various risk factors and pathophysiologies have been suggested for NAFLD, such as environmental, nutritional, genetic and immunological factors.,
The role of intestinal dysbiosis in NAFLD/NASH has also been debated by researchers., The human intestinal tract houses over 1014 bacteria involved in food digestion and interactions with the immune system, and many studies have highlighted their role in the pathogenesis of NAFLD., There are several mechanisms, which propose that microbiota dysregulation may affect NAFLD. Impaired intestinal permeability allows an increased flow of microbiota derivatives such as lipopolysaccharides (LPS) into the blood stream. Blood from the intestines is received by the liver through the portal vein, exposing it to the maximum concentration of these metabolites. While liver Kupffer cells destroy most of these intruders, the liver is still exposed to higher than normal levels of toxins, microbes and fatty compounds such as LPS, subsequently causing inflammation and damage to liver cells.,
Another mechanism by which intestinal microbiota can cause NAFLD is through ethanol-producing bacteria which produce ethanol by fermenting ingested sugar. Ethanol metabolism induces fatty acid synthesis which is deposited in the liver causing inflammation and damage.
Some studies have shown dysbiosis to cause NAFLD by over-activation of de novo lipogenesis (DNL), a regular metabolic process in which excess carbohydrates consumed are converted to fatty acids and stored in adipose tissue as triglycerides to be used for energy production through beta-oxidation when needed. Greater fatty acid production caused by the over-activation of DNL leads to fat accumulation in the liver, causing inflammation through oxidative stress and ultimately the development of NAFLD.,,
Finding the main variations in gut microbiota of patients with NAFLD, NASH, and healthy people may lead to novel strategies in the management of the disease., Thus, we aimed to compare the profile of gut microbiota in NAFLD and presumed NASH patients as well as healthy controls in this case–control study.
| Materials and Methods|| |
Fifty men 18–60 years of age referring to the Tehran Gastroenterology and Hepatology Clinic, were enrolled in the study. Given that previous studies have shown sex-related differences in gut microbiome, only men were included. In this study, we had 25 NAFLD, 13 presumed NASH, and 12 healthy participants as control. NAFLD was diagnosed by transient elastography (FibroScan, Echosence, France) and defined as having a controlled attenuation parameter (CAP) score above 260 dB/m. Presumed NASH was defined as having a CAP score above 260 dB/m along with a serum alanine transaminase (ALT) level >45 IU/L. Men with a CAP score below 260 dB/m who had normal ALT levels were chosen as controls.,
Exclusion criteria included having any of the following conditions: Hepatitis, autoimmune disorders, advanced liver disease, cancers, irritable bowel syndrome, inflammatory bowel disease, chronic diarrhea, liver enzymes 10 times above the normal values, any gastrointestinal surgeries, alcohol use >40 g per week, use of corticosteroids, probiotics, Vitamin E and fish oil supplements within 6 months, use of antibiotics within 6 weeks and dieting within 1 month.
The study protocol was approved by the Tehran University of Medical Sciences ethical committee and written consent was obtained from all participants.
An interviewer-administered questionnaire was completed for each participant obtaining information on demographics, past medical and surgical history, medication history, alcohol use, and smoking. To assess dietary intake, a 90-item, nonquantitative food frequency questionnaire was completed. Physical activity was assessed using the 7-question, International Physical Activity Questionnaire, through which the type, duration and difficulty level of different activities performed are questioned yielding a metabolic equivalent of task score for each individual. The Pittsburgh Sleep Quality Index was also measured using the standard validated questionnaire. Height, weight and waist and hip circumferences were measured using the National Institute of Health protocols.
Fasting blood was collected and serum levels of blood sugar, aspartate and ALT, cholesterol, triglycerides, very low-density lipoprotein (VLDL) were measured by BT 3000 Auto Analyzer. Antinuclear antibody, insulin, C-reactive protein, and hemoglobin A1C were performed by Immunoturbidimetric.
Stool collection kits were given to all participants and participants were instructed to return their specimen within the 24 h. The stool was collected directly into a sterile 20 ml polypropylene fecal container with a spoon attached to the snap-on lid and stored in-80°C.
DNA was extracted from stool samples using FavorPrep TM Stool DNA Isolation Mini Kits (FAVORGEN, Taiwan). DNA concentration was evaluated by Nanodrop (IMPLEN, Germany). Genomic DNA was polymerase chain reaction (PCR) amplified with primers (modified from the primer set employed by the Earth Microbiome Project (GTGYCAGCMGCCGCGGTAA and GGACTACNVGGGTWTCTAAT) targeting the V4 regions of microbial small subunit ribosomal Ribonucleic acid (RNA) genes. Amplicons were generated using a two-stage “targeted amplicon sequencing” protocol. The primers contained 5' common sequence tags (CS1 and CS2). First and second stages PCR amplifications were performed using MyTaq HS 2X mastermix (Bioline). PCR conditions for first PCR amplification was 95°C for 5 min, followed by 28 cycles of 95°C for 30”, 55°C for 45” and 72°C for 60”.
In second PCR amplification, each well received a separate primer pair with a unique 10-base barcode, obtained from the Access Array Barcode Library for Illumina (Fluidigm, South San Francisco, CA; Item# 100-4876). These primers contained the CS1 and CS2 linkers at the 3' ends. Cycling conditions were: 95°C for 5 min, followed by 8 cycles of 95°C for 30”, 60°C for 30” and 72°C for 30”.
Samples were then pooled in equal volume using an EpMotion5075 liquid handling robot (Eppendorf, Hamburg, Germany). The pooled library was purified using an AMPure XP cleanup protocol (0.6X, vol/vol; Agencourt, Beckmann-Coulter) to remove fragments smaller than 300 bp. The pooled libraries, with a 20% phiX spike-in, were loaded onto an Illumina MiniSeq mid-output flow cell. Based on the distribution of reads per barcode, the amplicons were re-pooled to generate a more balanced distribution of reads. The re-pooled library was purified using AMPure XP cleanup. The re-pooled libraries, with a 20% phiX spike-in, were loaded onto a Miniseq flow cell and sequenced. Fluidigm sequencing primers, targeting the CS1 and CS2 linker regions, were used to initiate sequencing. De-multiplexing of reads was performed on the instrument. Library preparation, pooling, and sequencing were performed at the University of Illinois at Chicago Genome Research Core within the Research Resources Center.
Forward and reverse reads were merged using PEAR. Merged reads were trimmed to remove ambiguous nucleotides, primer sequences, and trimmed based on the quality threshold of P = 0.01. Reads that lacked either primer sequence or any sequences <225 bp were discarded. Chimeric sequences were identified and removed using the USEARCH algorithm with a comparison to Silva v132 reference sequence database.,
The standard QIIME pipeline was modified to generate taxonomic summaries using sub-OTU resolution of the sequence dataset., Briefly, the resulting sequence files were then merged with sample information. All sequences were then dereplicated to produce a list of unique sequences. All sequences that had an abundance of at least 10 counts were designated seed sequences. USEARCH was then used to find the nearest seed sequence for any non-seed sequence with a minimum identity threshold of 97%. For any non-seed sequence that matched a seed sequence, its counts were merged with the seed sequence counts. For any non-seed sequence that did not match a seed sequence, it would remain an independent sequence.
Taxonomic annotations for seed and unmatched non-seed sequences were assigned using the USEARCH and Silva v132 reference with a minimum similarity threshold of 90%., In order to improve depth of annotation, the standard QIIME assignment algorithm was modified to only consider hits at each taxonomic level that had an assigned name. Furthermore, any hits in the reference database must have a minimum identity of 97% or 99% to be considered for genus or species level assignment, respectively. Taxonomic annotations and sequence abundance data were then merged into a single sequence table.
Analyses method for base-line characteristics
Quantile–quantile plot was used to assess normality of data. Mean ± standard deviation and median (range) were calculated for normal and skewed variables, respectively.
One-way analysis of variance (ANOVA) and Kruskal–Wallis tests were performed for normal and nonnormal variables, respectively. A P < 0.05 was considered to be statistically significant. Statistical analysis was performed using Stata version 12(StataCorp, Texas, USA).
Differential analysis of amplicon sequence data
Differential analyses of taxa as compared with experimental covariates were performed using the software package edge R on raw sequence counts. Prior to analysis, the data were filtered to remove any sequences that were annotated as chloroplast or mitochondria in origin as well as removing taxa that accounted for <0.1% of the total sequence counts. Data were normalized as counts per million. Normalized data were then fit using a negative binomial generalized linear model (GLM) using experimental covariates, and statistical tests were performed using a likelihood ratio test. Adjusted P values were calculated using the Benjamini–Hochberg false discovery rate (FDR) correction. Significant taxa were determined based on an FDR threshold of 5% (0.05).
Alpha and beta diversity analysis of amplicon sequence data
Shannon indices were calculated with default parameters in R using the vegan library. Prior to analysis, the data were rarefied to a depth of 6500 counts per sample. The resulting Shannon indices were then modeled with the sample covariates using a GLM assuming a Gaussian distribution. Significance of the model ANOVA was tested using the F test. Post hoc, pair wise tests were performed using Mann–Whitney test. Plots were generated in R using the ggplot2 library.
Bray-Curtis indices were calculated with default parameters in R using the vegan library. Prior to analysis the normalized data were square root transformed. The resulting dissimilarity indices were modeled and tested for significance with the sample covariates using the ADONIS test. Additional comparisons of the individual covariates were also performed using ANOSIM. Plots were generated in R using the ggplot2 library.
| Results|| |
Overall, 50 men (41.3 ± 8.9 years) entered this study, 12, 25, and 13 of whom were in the control, NAFLD, and presumed NASH groups, respectively [Table 1]. The number of individuals overweight or obese (body mass index [BMI] >25) was significantly greater in the presumed NASH group compared to NAFLD and in the presumed NASH and NAFLD groups, compared to controls. Serum triglycerides, VLDL and aspartate transaminase were also significantly greater in the presumed NASH and NAFLD groups (P < 0.05).
|Table 1: Demographic and paraclinic characteristics of study participants|
Click here to view
Alpha and beta diversity
Beta diversity was not significantly different at the genus level using ADONIS (P = 0.466, R2 = 0.041) and ANOSIM (P = 0.573, R = −0.0121) methods [Figure 1]. Alpha diversity assessed at the genus level via the Simpson (A), Evenness (B) and Richness (C) methods did not yield significant differences either, P = 0.482, P = 0.573 and P = 0.464, respectively [Figure 2].
|Figure 1: Nonmetric multidimensional scaling plot of Bray-Curtis dissimilarity indices computed using the taxonomic summary data at the genus level. Points are colored based on sample group|
Click here to view
|Figure 2: Dot and box plots of (a) Shannon, (b) Pielou's eveness and (c) richness diversity indices as compared with sample group|
Click here to view
The relative average phyla distribution of gut microbiomes in the NAFLD, presumed NASH, and control groups are reported in [Figure 3]. Firmicutes, Actinobacteria, and Bacteroidetes were dominant in all study groups; differences in these phyla were insignificant. At the phyla level, there were no significant differences. At the Genus level, only Veillonella were found to be significantly different among the study groups [Figure 4]. Veillonella was more abundant in the presumed NASH group than the NAFLD group (P = 2.76 × 10−6, q = 2.07 × 10−4, logFC = 5.52).
|Figure 3: Relative sequence abundance of main bacterial phyla as compared with sample group|
Click here to view
|Figure 4: Relative sequence abundance of genus Veillonella as compared with sample group. *Veillonella was significantly different between presumed nonalcoholic steatohepatitis and nonalcoholic fatty liver disease group|
Click here to view
| Discussion|| |
Intestinal microbial composition affects host metabolism. Accumulating evidence suggests a relationship between microbial composition and fatty liver pathogenesis., In this study, we compared gut microbiota in presumed NASH, NAFLD and healthy individuals at phyla, class, order, family and genus levels, and accounted for factors (diet, sleep quality and physical activity) that were shown to be associated with microbial composition as well as NAFLD in the prior studies.
Overall, our results indicate no significant difference in alpha and beta diversity, among the study groups. Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria have been recognized as the four dominant bacterial phyla. Some studies have reported Firmicutes and Bacteroidetes to be the dominant phyla,,, while our study, found Firmicutes, followed by Actinobacteria and Bacteroidetes to be dominant. These variations can be explained by ethnicity, genetics, diet, and lifestyles.
Recent systematic reviews in NAFLD/NASH patients have also shown heterogeneity in biodiversity at both phyla and genus levels., While some have reported significant differences at the phyla level between NASH patients and controls,,, others observed no significant differences., Among those reporting significant differences, the composition in the cases and controls did not follow a similar pattern, consistent with our findings; in some studies, Bacteroidetes increased in NASH patients and Firmicutes decreased,,, while in others it was the opposite. In this study, while Firmicutes increased in presumed NASH patients compared to the NAFLD and controls, no significant differences were observed at the phyla level.
Inconsistent results may be explained by technical heterogeneity such as stool sampling and storage, as well as DNA extraction methods. Ethnicity and genetics have also been proposed to influence gut microbial composition. We used 16s rRNA analysis as the standard method to study taxonomic and phylogenetic composition of microbiota., This study is the first evaluation of gut microbiota in NAFLD/presumed NASH patients in Middle East Region; hence, different results compared to studies in other countries were expected. Interestingly, some studies have shown that even within a single nation with similar lifestyles and dietary habits, microbial composition can vary significantly.
BMI >25 was found to be significantly different among NAFLD/presumed NASH patients and controls. Prior to adjusting, microbial composition varied significantly, so the association of BMI with NAFLD/NASH may not be independent from microbial composition.
Veillonella was the only genus exhibiting a significant difference between the NAFLD and presumed NASH groups, with greater abundance among those with presumed NASH. Given that some studies have shown Veillonella to be greater in cirrhotic patients,, it is possible that this genus increases as liver disease progresses.
Veillonella, highly recognized for its involvement in lactate fermentation, causing the release of acetate and propionate. Buildup of high amounts of these compounds has been previously shown to trigger gluconeogenesis and lipogenesis, increasing lipid storage in the liver and body tissues. Veillonella has been shown to play a role in small intestinal bacterial overgrowth (SIBO), which occurs more frequently in overweight/obese individuals. In addition, SIBO is correlated with increased TLR4 expression and Interleukin 8 secretions, both of which affect inflammatory pathways involved in NAFLD pathogenesis.
In this study, many factors, known to affect microbial composition in NAFLD/NASH patients were measured and controlled. For example, given that even small alterations in diet can affect the microbial composition, only individuals who had stable, routine diets in the month prior to study recruitment were eligible to enter the study and dietary information was obtained for all individuals to control for variations. As expected, individuals with NAFLD and NASH consumed higher overall calories as well as calories from fat compared to healthy controls.
| Conclusion|| |
The dominant phyla in this study population were different from those of many other populations; however, these results were not significantly different among the healthy individuals and those with fatty liver. Larger, longitudinal cohort studies are needed to better control for all factors affecting NAFLD/presumed NASH development and capture long-term changes in microbial composition of this patient population.
This study was performed at the Digestive Diseases Research Institute, Tehran University of Medical Sciences (TUMS), Tehran, Iran (96-02-37-35077). The study protocol was approved by the TUMS Ethical Committee (IR.TUMS.VCR.REC.1396.3322). Sequencing in the project was performed by the University of Illinois at Chicago Genome Research Core within the Research Resources Center. Basic processing of the raw data was performed by the University of Illinois at Chicago Research Informatics Core.
Financial support and sponsorship
Funding for this study was provided by the Digestive Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran (96-02-37-35077).
Conflicts of interest
There are no conflicts of interest.
| References|| |
Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liverdisease – Meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology 2016;64:73-84.
Sattar N, Forrest E, Preiss D. Non-alcoholic fatty liver disease. BMJ 2014;349:g4596.
Vanni E, Bugianesi E, Kotronen A, De Minicis S, Yki-Järvinen H, Svegliati-Baroni G. From the metabolic syndrome to NAFLD or vice versa? Dig Liver Dis 2010;42:320-30.
Rotman Y, Koh C, Zmuda JM, Kleiner DE, Liang TJ; NASH CRN. The association of genetic variability in patatin-like phospholipase domain-containing protein 3 (PNPLA3) with histological severity of nonalcoholic fatty liver disease. Hepatology 2010;52:894-903.
Marchesini G, Petta S, Dalle Grave R. Diet, weight loss, and liver health in nonalcoholic fatty liver disease: Pathophysiology, evidence, and practice. Hepatology 2016;63:2032-43.
Wieland A, Frank DN, Harnke B, Bambha K. Systematic review: Microbial dysbiosis and nonalcoholic fatty liver disease. Aliment Pharmacol Ther 2015;42:1051-63.
Mohammadi Z, Poustchi H, Motamed-Gorji N, Eghtesad S, Hekmatdoost A, Saniee P, et al.
Fecal microbiota in non-alcoholic fatty liver disease and non-alcoholic steatohepatitis: A systematic review. Arch Iran Med 2020;23:44-52.
Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host-bacterial mutualism in the human intestine. Science 2005;307:1915-20.
Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, et al.
Metagenomic analysis of the human distal gut microbiome. Science 2006;312:1355-9.
Mazzotti A, Caletti MT, Sasdelli AS, Brodosi L, Marchesini G. Pathophysiology of nonalcoholic fatty liver disease: Lifestyle-gut-gene interaction. Dig Dis 2016;34 Suppl 1:3-10.
Poeta M, Pierri L, Vajro P. Gut-liver axis derangement in non-alcoholic fatty liver disease. Children (Basel) 2017;4:E66.
Matsushita N, Osaka T, Haruta I, Ueshiba H, Yanagisawa N, Omori-Miyake M, et al.
Effect of lipopolysaccharide on the progression of non-alcoholic fatty liver disease in high caloric diet-fed mice. Scand J Immunol 2016;83:109-18.
Reid DT, McDonald B, Khalid T, Vo T, Schenck LP, Surette MG, et al.
Unique microbial-derived volatile organic compounds in portal venous circulation in murine non-alcoholic fatty liver disease. Biochim Biophys Acta 2016;1862:1337-44.
Elshaghabee FM, Bockelmann W, Meske D, de Vrese M, Walte HG, Schrezenmeir J, et al.
Ethanol production by selected intestinal microorganisms and lactic acid bacteria growing under different nutritional conditions. Front Microbiol 2016;7:47.
Paglialunga S, Dehn CA. Clinical assessment of hepatic de novo
lipogenesis in non-alcoholic fatty liver disease. Lipids Health Dis 2016;15:159.
Sanders FW, Griffin JL. De novo
lipogenesis in the liver in health and disease: More than just a shunting yard for glucose. Biol Rev Camb Philos Soc 2016;91:452-68.
Basaranoglu M, Basaranoglu G, Bugianesi E. Carbohydrate intake and nonalcoholic fatty liver disease: Fructose as a weapon of mass destruction. Hepatobiliary Surg Nutr 2015;4:109-16.
Kim YS, Unno T, Kim BY, Park MS. Sex differences in gut microbiota. World J Mens Health 2020;38:48-60.
Jamali R, Pourshams A, Amini S, Deyhim MR, Rezvan H, Malekzadeh R. The upper normal limit of serum alanine aminotransferase in Golestan Province, northeast Iran. Arch Iran Med 2008;11:602-7.
Kabir A, Pourshams A, Khoshnia M, Malekzadeh F. Normal limit for serum alanine aminotransferase level and distribution of metabolic factors in old population of Kalaleh, Iran. Hepat Mon 2013;13:e10640.
Vasheghani-Farahani A, Tahmasbi M, Asheri H, Ashraf H, Nedjat S, Kordi R. The Persian, last 7-day, long form of the International Physical Activity Questionnaire: Translation and validation study. Asian J Sports Med 2011;2:106-16.
Farrahi Moghaddam J, Nakhaee N, Sheibani V, Garrusi B, Amirkafi A. Reliability and validity of the Persian version of the Pittsburgh Sleep Quality Index (PSQI-P). Sleep Breath 2012;16:79-82.
Naqib A, Poggi S, Wang W, Hyde M, Kunstman K, Green SJ. Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol. In: Gene Expression Analysis. New York, NY: Humana Press; 2018. p. 149-69.
Moonsamy PV, Williams T, Bonella P, Holcomb CL, Höglund BN, Hillman G, et al.
High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ system for simplified amplicon library preparation. Tissue Antigens 2013;81:141-9.
Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: A fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics (Oxford, England) 2014;30:614-20.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010;26:2460-1.
Glöckner FO, Yilmaz P, Quast C, Gerken J, Beccati A, Ciuprina A, et al
. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J Biotechnol 2017;261:169-76.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al.
QIIME allows analysis of high-throughput community sequencing data. Nat Meth 2010;7:335-6.
Tikhonov M, Leach RW, Wingreen NS. Interpreting 16S metagenomic data without clustering to achieve sub-OTU resolution. ISME J 2015;9:68-80.
Khanna S, Tosh PK. A clinician's primer on the role of the microbiome in human health and disease. Mayo Clin Proc 2014;89:107-14.
Zhu L, Baker SS, Gill C, Liu W, Alkhouri R, Baker RD, et al.
Characterization of gut microbiomes in nonalcoholic steatohepatitis (NASH) patients: A connection between endogenous alcohol and NASH. Hepatology 2013;57:601-9.
Jiang W, Wu N, Wang X, Chi Y, Zhang Y, Qiu X, et al.
Dysbiosis gut microbiota associated with inflammation and impaired mucosal immune function in intestine of humans with non-alcoholic fatty liver disease. Sci Rep 2015;5:8096.
Sobhonslidsuk A, Chanprasertyothin S, Pongrujikorn T, Kaewduang P, Promson K, Petraksa S, et al.
The association of gut microbiota with nonalcoholic steatohepatitis in thais. Biomed Res Int 2018;2018:9340316.
Rinninella E, Raoul P, Cintoni M, Franceschi F, Miggiano GA, Gasbarrini A, et al.
What is the healthy gut microbiota composition? A changing ecosystemacross age, environment, diet, and diseases. Microorganisms 2019;7:E14.
Vernekar M, Singhal R, Joshi K, Amarapurkar D. Variation in the plasma levels of polyunsaturated fatty acids in control vis-avis nonalcoholic fatty liver disease subjects and its possible association with gut microbiome. Metab Syndr Relat Disord 2018;16:329-35.
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:120-7.
Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature 2012;486:207-14.
Claesson MJ, O'Toole PW. Evaluating the latest high-throughput molecular techniques for the exploration of microbial gut communities. Gut Microbes 2010;1:277-8.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al.
A human gut microbial gene catalogue established by metagenomic sequencing. Nature 2010;464:59-65.
Chen Y, Ji F, Guo J, Shi D, Fang D, Li L. Dysbiosis of small intestinal microbiota in liver cirrhosis and its association with etiology. Sci Rep 2016;6:34055.
Ponziani FR, Bhoori S, Castelli C, Putignani L, Rivoltini L, Del Chierico F, et al.
Hepatocellular carcinoma is associated with gut microbiota profile and inflammation in nonalcoholic fatty liver disease. Hepatology 2019;69:107-20.
Mashima I, Nakazawa F. The influence of oral Veillonella species on biofilms formed by Streptococcus species. Anaerobe 2014;28:54-61.
Sohail MU, Althani A, Anwar H, Rizzi R, Marei HE. Role of the gastrointestinal tract microbiome in the pathophysiology of diabetes mellitus. J Diabetes Res 2017;2017:9631435.
Wigg AJ, Roberts-Thomson IC, Dymock RB, McCarthy PJ, Grose RH, Cummins AG. The role of small intestinal bacterial overgrowth, intestinal permeability, endotoxaemia, and tumour necrosis factor alpha in the pathogenesis of non-alcoholic steatohepatitis. Gut 2001;48:206-11.
Sabaté JM, Jouët P, Harnois F, Mechler C, Msika S, Grossin M, et al.
High prevalence of small intestinal bacterial overgrowth in patients with morbid obesity: A contributor to severe hepatic steatosis. Obes Surg 2008;18:371-7.
Shanab AA, Scully P, Crosbie O, Buckley M, O'Mahony L, Shanahan F, et al
. Small intestinal bacterial overgrowth in nonalcoholic steatohepatitis: Association with toll-like receptor 4 expression and plasma levels of interleukin 8. Dig Dis Sci 2011;56:1524-34.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]