Mucosal Genomics Implicate Lymphocyte Activation and Lipid Metabolism in Refractory Environmental Enteric Dysfunction

Yael Haberman, Najeeha T Iqbal, Sudhir Ghandikota, Indika Mallawaarachchi, Tzipi Braun, Phillip J Dexheimer, Najeeb Rahman, Rotem Hadar, Kamran Sadiq, Zubair Ahmad, Romana Idress, Junaid Iqbal, Sheraz Ahmed, Aneeta Hotwani, Fayyaz Umrani, Lubaina Ehsan, Greg Medlock, Sana Syed, Chris Moskaluk, Jennie Z Ma, Anil G Jegga, Sean R Moore, Syed Asad Ali, Lee A Denson, Yael Haberman, Najeeha T Iqbal, Sudhir Ghandikota, Indika Mallawaarachchi, Tzipi Braun, Phillip J Dexheimer, Najeeb Rahman, Rotem Hadar, Kamran Sadiq, Zubair Ahmad, Romana Idress, Junaid Iqbal, Sheraz Ahmed, Aneeta Hotwani, Fayyaz Umrani, Lubaina Ehsan, Greg Medlock, Sana Syed, Chris Moskaluk, Jennie Z Ma, Anil G Jegga, Sean R Moore, Syed Asad Ali, Lee A Denson

Abstract

Background & aims: Environmental enteric dysfunction (EED) limits the Sustainable Development Goals of improved childhood growth and survival. We applied mucosal genomics to advance our understanding of EED.

Methods: The Study of Environmental Enteropathy and Malnutrition (SEEM) followed 416 children from birth to 24 months in a rural district in Pakistan. Biomarkers were measured at 9 months and tested for association with growth at 24 months. The duodenal methylome and transcriptome were determined in 52 undernourished SEEM participants and 42 North American controls and patients with celiac disease.

Results: After accounting for growth at study entry, circulating insulin-like growth factor-1 (IGF-1) and ferritin predicted linear growth, whereas leptin correlated with future weight gain. The EED transcriptome exhibited suppression of antioxidant, detoxification, and lipid metabolism genes, and induction of anti-microbial response, interferon, and lymphocyte activation genes. Relative to celiac disease, suppression of antioxidant and detoxification genes and induction of antimicrobial response genes were EED-specific. At the epigenetic level, EED showed hyper-methylation of epithelial metabolism and barrier function genes, and hypo-methylation of immune response and cell proliferation genes. Duodenal coexpression modules showed association between lymphocyte proliferation and epithelial metabolic genes and histologic severity, fecal energy loss, and wasting (weight-for-length/height Z < -2.0). Leptin was associated with expression of epithelial carbohydrate metabolism and stem cell renewal genes. Immune response genes were attenuated by giardia colonization.

Conclusions: Children with reduced circulating IGF-1 are more likely to experience stunting. Leptin and a gene signature for lymphocyte activation and dysregulated lipid metabolism are implicated in wasting, suggesting new approaches for EED refractory to nutritional intervention. ClinicalTrials.gov, Number: NCT03588013. (https://ichgcp.net/clinical-trials-registry/NCT03588013).

Keywords: Anthropometrics; DNA Methylation; Intestine; RNA Sequencing.

Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.

Figures

Graphical abstract
Graphical abstract
Figure 1
Figure 1
The EED intestinal transcriptome and enriched biologic pathways. The core EED transcriptome was comprised of 1,262 genes (481 down- and 781 up-regulated) differentially expressed between 31 AKU-EED malnourished cases and 21 Cincinnati well-nourished controls in the training set (FDR A) demonstrating the averaged normalized expression in AKU-EED malnourished cases and Cincinnati well-nourished controls in the training and validation groups for the top differentially expressed genes (more detailed heatmap in Supplementary Figure 4). (B) Principal component analysis (PCA) using the 1,262 EED genes transcriptome (determined only using the training subset) showing separation of the AKU-EED malnourished cases and the well-nourished controls from Cincinnati in both the training and validation groups on the PC1 axis that explains 38% of the total variance in gene expression. Functional enrichment analyses of the 781 up- (C) and 481 down-regulated (D) genes between AKU-EED malnourished cases and Cincinnati well-nourished controls was performed using ToppGene/ToppCluster and was visualized using Cytoscape.
Figure 2
Figure 2
Shared and disease-specific immune and metabolic intestinal gene expression features of EED and celiac disease. (A) Representative hematoxylin and eosin stained duodenal biopsy specimens from a Cincinnati well-nourished control, a Cincinnati celiac disease patient (Marsh celiac disease score 3a; EED histology score of 12), a malnourished AKU-EED-1 case with EED histology score of 9, and a malnourished AKU-EED-2 case with EED histology score of 4 are shown. ∗Paneth cells in a Cincinnati well-nourished control. Arrow indicates villous blunting and arrowhead indicates intraepithelial lymphocytes in a patient from Cincinnati with celiac disease and a malnourished AKU-EED case. Bar equals 247 μm. (B) The Venn diagram shows the overlap between the 718 genes comprising the celiac disease transcriptome (differentially expressed genes between 17 patients from Cincinnati with celiac disease and 25 well-nourished controls from Cincinnati, FDR < 0.05 and fold change [FC] ≥ 1.5 using bulk RNASeq of duodenal RNA) and 1,262 genes comprising the EED transcriptome. This demonstrates 212 shared down- and 85 shared up-regulated genes. (C) Unsupervised hierarchical clustering heatmap with the top differentially expressed genes in the EED transcriptome demonstrating the averaged normalized expression across malnourished AKU-EED cases, patients from Cincinnati with celiac disease, and Cincinnati well-nourished controls. Functional enrichment analysis of the up- (D) and down-regulated (E) shared and unique genes in the EED and celiac disease transcriptomes was performed using ToppGene/ToppCluster and was visualized using Cytoscape. (F) Immunohistochemistry was performed using antibodies against DUOX2 (yellow chromogen) and LCN2 (teal chromogen) in a dual stain. Original magnification x200 for i & ii. (G) Data for the relative tissue area exhibiting staining for the analytes, normalized against the total area of tissue in each sample, are shown for controls (n = 10), celiac disease (n = 10), and EED (n = 57); Kruskal-Wallis test with Dunn multiple comparisons test; ∗∗P < .01; ∗P < .05.
Figure 3
Figure 3
Variation in DNAm associated with expression of immune and metabolic genes in EED. Genome-wide intestinal DNAm was profiled in DNA prepared from duodenal biopsy specimens using the Illumina Infinium MethylationEPIC BeadChip platform. (A) A Manhattan plot is shown displaying the overlapping DMRs associated with EED in 2 methylation profile batches including 31 malnourished AKU-EED cases compared with 21 well-nourished Cincinnati controls in batch 1, and 33 malnourished AKU-EED cases compared with 9 well-nourished Cincinnati controls in batch 2, of which 12 AKU-EED cases and 5 Cincinnati controls were tested in both batches. The corrected P values (−log10 Stouffer) of each DMR are plotted against their respective positions on each chromosome. (B) The Venn diagram shows the overlap between 481 down- and 781 up-regulated genes in the EED transcriptome and DMRs highlighting 453 rDMR including genes that show evidence for both differential methylation (DM) and differential expression (DE). Beta-value methylation levels of differentially methylated points within rDMR showing a significant relationship (P < 1E-6) between methylation levels and expression (TPM) of specific down- (C) and up-regulated (D) genes as indicated. We highlight genes that are expressed in intestinal epithelial cells based upon a previous isolated ileal epithelial cell dataset and single-cell datasets. The gray lines illustrate a linear model fit, whereas rho values indicate the Spearman correlation coefficients.
Figure 4
Figure 4
Gene coexpression modules are associated with EED diagnosis and measures of clinical and histologic severity. WGCNA was implemented to identify modules of coexpressed genes. For each module, the first principal component, referred to as the eigengene, was considered to be the module representative tested for association with phenotypic traits. (A) Heatmap representation of the WGCNA demonstrates gene coexpression modules (represented by module eigengenes), which were correlated with EED diagnosis (first column, EED), in an analysis that included 52 malnourished AKU-EED cases and 25 well-nourished Cincinnati controls and other clinical traits as shown in an analysis limited to the 52 malnourished AKU-EED cases. Seven gene coexpression modules identified based on the correlation strength with the EED diagnosis are shown, together with the results for correlations with clinical traits within the AKU-EED cases. Data are shown as the correlation coefficient and P value for each comparison. (B) A representation plot of hierarchically clustered selected top functionally enriched (FDR P < .05) biological processes, pathways, and cell types in each of the 7 gene coexpression modules is shown. The size of the circles and the intensity of the color is proportional to the enrichment strength. (C) Hub genes (ovals) are shown that were also differentially methylated together with functionally enriched pathways from 3 gene coexpression modules that were strongly correlated with WHZ (salmon, pink, and brown modules). (D) A heatmap of specific hub genes from (C) and their correlation with EED diagnosis and other clinical traits including WHZ around birth as indicated is shown.

References

    1. Victora C.G., Adair L., Fall C. Maternal and child undernutrition: consequences for adult health and human capital. Lancet. 2008;371:340–357.
    1. Lunn P.G., Northrop-Clewes C.A., Downes R.M. Intestinal permeability, mucosal injury, and growth faltering in Gambian infants. Lancet. 1991;338:907–910.
    1. Fasano A., Catassi C. Clinical practice. Celiac disease. N Engl J Med. 2012;367:2419–2426.
    1. Campbell D.I., Murch S.H., Elia M. Chronic T cell-mediated enteropathy in rural west African children: relationship with nutritional status and small bowel function. Pediatr Res. 2003;54:306–311.
    1. Chen R.Y., Kung V.L., Das S. Duodenal microbiota in stunted undernourished children with enteropathy. N Engl J Med. 2020;383:321–333.
    1. Uchiyama R., Kupkova K., Shetty S.J. Histone H3 lysine 4 methylation signature associated with human undernutrition. Proc Natl Acad Sci U S A. 2018;115:E11264–E11273.
    1. Schulze K.V., Swaminathan S., Howell S. Edematous severe acute malnutrition is characterized by hypomethylation of DNA. Nat Commun. 2019;10:5791.
    1. Iqbal N.T., Syed S., Sadiq K. Study of Environmental Enteropathy and Malnutrition (SEEM) in Pakistan: protocols for biopsy based biomarker discovery and validation. BMC Pediatr. 2019;19:247.
    1. Aguayo V.M., Badgaiyan N., Qadir S.S. Community Management of Acute Malnutrition (CMAM) programme in Pakistan effectively treats children with uncomplicated severe wasting. Matern Child Nutr. 2018;14(Suppl 4)
    1. Liu T.C., VanBuskirk K., Ali S.A. A novel histological index for evaluation of environmental enteric dysfunction identifies geographic-specific features of enteropathy among children with suboptimal growth. PLoS Negl Trop Dis. 2020;14
    1. Scanlon S.A., Murray J.A. Update on celiac disease - etiology, differential diagnosis, drug targets, and management advances. Clin Exp Gastroenterol. 2011;4:297–311.
    1. Haberman Y., Karns R., Dexheimer P.J. Ulcerative colitis mucosal transcriptomes reveal mitochondriopathy and personalized mechanisms underlying disease severity and treatment response. Nat Commun. 2019;10:38.
    1. Howell K.J., Kraiczy J., Nayak K.M. DNA methylation and transcription patterns in intestinal epithelial cells from pediatric patients with inflammatory bowel diseases differentiate disease subtypes and associate with outcome. Gastroenterology. 2018;154:585–598.
    1. Zhang B., Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4
    1. Moore S.R., Quinn L.A., Maier E.A. Oral alanyl-glutamine for inflammation, nutrition, & enteropathy: a randomized dose-response controlled trial in children. J Pediatr Gastroenterol Nutr. 2020;71:393–400.
    1. Iqbal N.T., Sadiq K., Syed S. Promising biomarkers of environmental enteric dysfunction: a prospective cohort study in Pakistani children. Sci Rep. 2018;8:2966.
    1. Haberman Y., Tickle T.L., Dexheimer P.J. Pediatric Crohn disease patients exhibit specific ileal transcriptome and microbiome signature. J Clin Invest. 2014;124:3617–3633.
    1. Bragde H., Jansson U., Jarlsfelt I. Gene expression profiling of duodenal biopsies discriminates celiac disease mucosa from normal mucosa. Pediatr Res. 2011;69:530–537.
    1. Cavalli G., Heard E. Advances in epigenetics link genetics to the environment and disease. Nature. 2019;571:489–499.
    1. Dolinoy D.C., Huang D., Jirtle R.L. Maternal nutrient supplementation counteracts bisphenol A-induced DNA hypomethylation in early development. Proc Natl Acad Sci U S A. 2007;104:13056–13061.
    1. Dominguez-Salas P., Moore S.E., Baker M.S. Maternal nutrition at conception modulates DNA methylation of human metastable epialleles. Nat Commun. 2014;5:3746.
    1. Smillie C.S., Biton M., Ordovas-Montanes J. Intra- and inter-cellular rewiring of the human colon during ulcerative colitis. Cell. 2019;178 714e22–730e22.
    1. Haberman Y., Minar P., Karns R. Mucosal inflammatory and wound healing gene programs reveal targets for stricturing behavior in pediatric Crohn’s disease. J Crohns Colitis. 2020;15:273–286.
    1. Chama M., Amadi B.C., Chandwe K. Transcriptomic analysis of enteropathy in Zambian children with severe acute malnutrition. EBioMedicine. 2019;45:456–463.
    1. Gehrig J.L., Venkatesh S., Chang H.W. Effects of microbiota-directed foods in gnotobiotic animals and undernourished children. Science. 2019;365
    1. Jones K.D., Hunten-Kirsch B., Laving A.M. Mesalazine in the initial management of severely acutely malnourished children with environmental enteric dysfunction: a pilot randomized controlled trial. BMC Med. 2014;12:133.
    1. Gazi M.A., Das S., Siddique M.A. Plasma kynurenine to tryptophan ratio is negatively associated with linear growth of children living in a slum of Bangladesh: results from a community-based intervention study. Am J Trop Med Hyg. 2021;104:766–773.
    1. Naylor C., Lu M., Haque R. Environmental enteropathy, oral vaccine failure and growth faltering in infants in Bangladesh. EBioMedicine. 2015;2:1759–1766.
    1. Hartmann W., Brunn M.L., Stetter N. Helminth infections suppress the efficacy of vaccination against seasonal influenza. Cell Rep. 2019;29:2243e4–2256e4.
    1. Zhang Y., Zhou J., Niu F. Characterizing early child growth patterns of height-for-age in an urban slum cohort of Bangladesh with functional principal component analysis. BMC Pediatr. 2017;17:84.
    1. Manger M.S., Taneja S., Strand T.A. Poor folate status predicts persistent diarrhea in 6- to 30-month-old north Indian children. J Nutr. 2011;141:2226–2232.
    1. Yu D.H., Gadkari M., Zhou Q. Postnatal epigenetic regulation of intestinal stem cells requires DNA methylation and is guided by the microbiome. Genome Biol. 2015;16
    1. da Silva A.V.A., Oliveira S.B.D., Di Rienzi S.C. Murine methyl donor deficiency impairs early growth in association with dysmorphic small intestinal crypts and reduced gut microbial community diversity. Curr Dev Nutr. 2018;3:nzy070.
    1. Kaimal V., Bardes E.E., Tabar S.C. ToppCluster: a multiple gene list feature analyzer for comparative enrichment clustering and network-based dissection of biological systems. Nucleic Acids Res. 2010;38:W96–W102.
    1. Saito R., Smoot M.E., Ono K. A travel guide to Cytoscape plugins. Nat Methods. 2012;9:1069–1076.

Source: PubMed

Patrocinadores e Colaboradores

Condições médicas

Intervenções de drogas

3
Se inscrever