Cross talk between the liver microbiome and epigenome in patients with metabolic dysfunction-associated steatotic liver disease.

BACKGROUND: The pathogenesis of MASLD (metabolic dysfunction-associated steatotic liver disease), including its severe clinical forms, involves complex processes at all levels of biological organization. This study examined the potential link between the liver microbiome profile and epigenetic factors. METHODS: Liver microbial DNA composition was analysed using high throughput 16S rRNA gene sequencing in 116 individuals, with 55% being female, across the spectrum of liver disease severity. Total activity of histone deacetylases (HDACs) and acetyltransferases (HATs) was assayed in nuclear extracts from fresh liver samples. In addition, we measured the global 5-hydroxymethylcytosine (5-hmC) levels of liver DNA. FINDINGS: Patients with MASLD showed a 2.07-fold increase (p = 0.013) in liver total HAT activity. Moreover, a correlation was observed between liver total HAT activity and the score for histological steatosis (Spearman's R = 0.60, p = 1.0E-3) and disease severity (R = 0.40, p = 2.0E-2). Liver HAT and HDAC activities also showed associations with the abundance of several liver bacterial DNAs. Additionally, liver global levels of 5-hmC showed negative correlation with the read number of Bacteroidetes (R = -0.62, p = 9.3E-4) and Gammaproteobacteria (R = -0.43, p = 3.2E-2), while it was positively correlated with the abundance of Acidobacteria (R = 0.42, p = 4.1E-2) and Actinobacteria (R = 0.47, p = 1.8E-2). INTERPRETATION: The host liver epigenome, including the activity of enzymes involved in maintaining the balance between protein acetylation and deacetylation and the global DNA hydroxy-methylation status, may be the target of microbial signals. FUNDING: Agencia Nacional de Promoción Científica y Tecnológica, FonCyT.

Computational Pipeline for Next-Generation Sequencing (NGS) Studies in Genetics of NASH.

High-throughput sequencing (HTS) technologies have contributed to expand current knowledge of the biology of complex diseases, including nonalcoholic fatty liver disease (NAFLD). Genome-wide association studies, whole exome sequencing, and sequencing of entire genes are used to identify variants and/or mutations that predispose to the disease pathogenesis. Here, we present a tutorial that may guide readers to manage high volume of genetics data in the context of Next-Generation Sequencing (NGS) studies.

The influence of host genetics on liver microbiome composition in patients with NAFLD.

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), Prevotella_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. FUNDING: This study was partially supported by grants PICT 2018-889, PICT 2019-0528, PICT2016-0135 and PICT 2018-0620 (Agencia Nacional de Promoción Científica y Tecnológica, FONCyT), CONICET Proyectos Unidades Ejecutoras 2017, PUE 0055.

Pleiotropy within gene variants associated with nonalcoholic fatty liver disease and traits of the hematopoietic system.

Genome-wide association studies of complex diseases, including nonalcoholic fatty liver disease (NAFLD), have demonstrated that a large number of variants are implicated in the susceptibility of multiple traits - a phenomenon known as pleiotropy that is increasingly being explored through phenome-wide association studies. We focused on the analysis of pleiotropy within variants associated with hematologic traits and NAFLD. We used information retrieved from large public National Health and Nutrition Examination Surveys, Genome-wide association studies, and phenome-wide association studies based on the general population and explored whether variants associated with NAFLD also present associations with blood cell-related traits. Next, we applied systems biology approaches to assess the potential biological connection/s between genes that predispose affected individuals to NAFLD and nonalcoholic steatohepatitis, and genes that modulate hematological-related traits-specifically platelet count. We reasoned that this analysis would allow the identification of potential molecular mediators that link NAFLD with platelets. Genes associated with platelet count are most highly expressed in the liver, followed by the pancreas, heart, and muscle. Conversely, genes associated with NAFLD presented high expression levels in the brain, lung, spleen, and colon. Functional mapping, gene prioritization, and functional analysis of the most significant loci (P < 1 × 10-8) revealed that loci involved in the genetic modulation of platelet count presented significant enrichment in metabolic and energy balance pathways. In conclusion, variants in genes influencing NAFLD exhibit pleiotropic associations with hematologic traits, particularly platelet count. Likewise, significant enrichment of related genes with variants influencing platelet traits was noted in metabolic-related pathways. Hence, this approach yields novel mechanistic insights into NAFLD pathogenesis.

Intrahepatic bacterial metataxonomic signature in non-alcoholic fatty liver disease.

OBJECTIVE: We 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 METHODS: Liver 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. RESULTS: Liver 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. CONCLUSION: This 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.

Genetic variation in long noncoding RNAs and the risk of nonalcoholic fatty liver disease.

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.

The response to postnatal stress: amino acids transporters and PKC activity.

It is well known that animals exposed to stressful stimuli during their early life develop different neurological disorders when they become adults. In this study, we evaluated the effect of acute cold stress on gamma-aminobutyric acid (GABA) and L-Serine (L-Ser) transporters in vitro, using the uptake of [(3)H]-GABA and [(3)H]L-Ser by synaptosomes-enriched fractions isolated from rat cerebral cortex during postnatal development. GABA and L-Ser uptake studies in vitro will be used in this investigation as a colateral evidence of changes in the expression of transporters of GABA and L-Ser. We observed that the maximum velocity (V (max)) in L-Ser and GABA uptake after stress session increased in all stages studied. In contrast, K (m) values of L-Ser uptake enhancent in almost age calculated, excluding at PD21 after cold stress during development, at the same time as K (m) (uptake affinity) values of GABA increased in just about age considered but not at PD5 compared with the control group. Finally we investigated the mechanism by which cells regulate the substrate affinity of L-Ser and GABA transporters. We demonstrated a significantly increase in total PKC activity to PD5 from PD21. Pretreatment with PKC inhibitor: staurosporine (SP) led to a restoration of control uptake in several postnatal-days suggesting a relationship between amino acids system and PKC activation. These findings suggest that a single exposure to postnatal cold stress at different periods after birth modifies both GABA and L-Ser transporters and the related increase in total PKC activity could be intracellular events that participate in neuronal plasticity by early life stress, which could be relevant to function of transporters in the adult rat brain.