ObjectiveWe aimed to characterise the liver tissue bacterial metataxonomic signature in two independent cohorts of patients with biopsy-proven non-alcoholic fatty liver disease (NAFLD) diagnosis, as differences in the host phenotypic features—from moderate to severe obesity—may be associated with significant changes in the microbial DNA profile.Design and methodsLiver tissue samples from 116 individuals, comprising of 47 NAFLD overweight or moderately obese patients, 50 NAFLD morbidly obese patients elected for bariatric surgery and 19 controls, were analysed using high-throughput 16S rRNA gene sequencing.ResultsLiver bacterial DNA profile significantly differs between morbidly obese and non-morbidly obese patients with NAFLD. Bacteroidetes (p=1.8e-18) and Firmicutes (p=0.0044) were over-represented in morbidly obese patients and Proteobacteria (p=5.2e-10)—specifically Gammaproteobacteria and Alphaproteobacteria, and Deinococcus-Thermus (p=0.00012)—were over-represented in the non-morbidly obese cohort. Cohort-specific analysis of liver microbial DNA signatures shows patterns linked to obesity. The imbalance in Proteobacteria (Alpha or Gamma) among non-morbidly obese patients, and Peptostreptococcaceae, Verrucomicrobia, Actinobacteria and Gamma Proteobacteria DNA among morbidly obese patients was associated with histological severity. Decreased amounts of bacterial DNA from the Lachnospiraceae family were associated with more severe histological features. Proteobacteria DNA was consistently associated with lobular and portal inflammation scores. Microbial DNA composition corresponded to predicted functional differences.ConclusionThis is the first comprehensive study showing that the liver tissue of NAFLD patients contains a diverse repertoire of bacterial DNA (up to 2.5×104 read counts). The liver metataxonomic signature may explain differences in the NAFLD pathogenic mechanisms as well as physiological functions of the host.

The human transcriptome comprises a myriad of non protein-coding RNA species, including long noncoding RNAs (lncRNAs), which have a remarkable role in transcriptional and epigenetic regulation. We hypothesized that variants in lncRNAs influence the susceptibility to nonalcoholic fatty liver disease (NAFLD). Using next generation sequencing, we performed a survey of genetic variation associated with randomly selected lncRNA-genomic regions located within both experimentally validated and computationally predicted regulatory elements. We used a two-stage (exploratory, n = 96 and replication, n = 390) case-control approach that included well-characterized patients with NAFLD diagnosed by liver biopsy. We sequenced > 263 megabase pairs at quality score > Q17, in a total of 2,027,565 reads, including 170 lncRNA-genomic regions. In the sequencing analysis and the validated dataset, we found that the rs2829145 A/G located in a lncRNA (lnc-JAM2-6) was associated with NAFLD and the disease severity. Prediction of regulatory elements in lnc-JAM2-6 showed potential sequence-specific binding motifs of oncogenes MAFK and JUND, and the transcription factor CEBPB that is involved in inflammatory response. The A-allele was significantly associated with NAFLD as disease trait (p = 0.0081) and the disease severity (NASH-nonalcoholic steatohepatitis vs. controls: OR 2.36 [95% CI: 1.54−3.62], p = 0.000078). The A-allele carriers also have significantly higher body mass index and glucose-related traits compared with homozygous GG. Hence, our results suggest that variation in lncRNAs contributes to NAFLD severity, while pointing toward the complexity of the genetic component of NAFLD, which involves still unexplored regulatory regions of the genome.

Background Human body microbiotas are influenced by several factors, including the interaction of the host with the environment and dietary preferences. The role of host genetics in modulating the liver microbiota in the context of NAFLD remains unknown. To address this gap, we examined the interplay between the liver metataxonomic profile and host genetics.Methods We obtained 16S rRNA gene sequences from liver biopsies and genotypes by Taqman-assays in 116 individuals. We compared taxon abundance at the genus level across host genotypes using dominant models of inheritance. We focused the analysis on variants influencing the risk/ protection against NAFLD-histological severity (PNPLA3-rs738409, TM6SF2-rs58542926, MBOAT7-rs641738, and HSD17B13-rs72613567) and a variant influencing macronutrient intake (FGF21-rs838133). We also explored the variants' combined effect via a polygenic risk score (PRS).Findings We identified at least 18 bacterial taxa associated with variants in the selected loci. Members of the Gammaproteobacteria class were significantly enriched in carriers of the rs738409 and rs58542926 risk-alleles, including Enterobacter (fold change [FC]=6.2) and Pseudoalteromonas (FC=2) genera, respectively. Lawsonella (1.6-FC), Prevo-tella_9 (FC=1.5), and Staphylococcus (FC=1.3) genera were enriched in rs838133-minor allele carriers, which is linked to sugar consumption and carbohydrate intake. Tyzzerella abundance (FC=2.64) exhibited the strongest association (p = 0.0019) with high PRS values (>4 risk alleles). The percentage of genus-level taxa variation explained by the PRS was »7.4%, independently of liver steatosis score and obesity.Interpretation We provided evidence that genetic variation may influence the liver microbial DNA composition. These observations may represent potentially actionable mechanisms of disease.