DOI: 10.5281/zenodo.8397632
title: Grandmaternal allergen sensitization reprograms epigenetic and airway responses to allergen in second generation offspring

description: <strong>SUPPLEMENTARY FIGURE LEGEND</strong> <strong>Supplemental Figure 1. Grandmaternal allergen sensitization during pregnancy regulates differentially methylated cytosines in airway epithelium of second-generation offspring at baseline and after allergen sensitization</strong>. Grandmaternal HDM exposure resulted in 1,390 DMCs in airway epithelium of F2 offspring at baseline (F0 HDM•F2 Base vs F0 Veh•F2 Base). F2 HDM sensitization alone resulted in 1,467 DMCs (F0Veh•F2Sens vs F0 Veh•F2 Base). In contrast, HDM sensitization in F2 mice from HDM-exposed founders resulted in 3,756 total DMCs (F0 HDM•F2 Sens vs F0 HDM•F2 Base). N=8 total for each group (n=4 males and n=4 females). <strong>Supplemental Figure 2. Gene ontology analysis of DMCs in airway epithelium of second-generation offspring</strong>. Gene ontology analysis of DMCs in airway epithelium of F2 mice identified enriched biological processes (individual small circles), which were then grouped based on functional annotation using Cytoscape Enrichment Map (larger, labeled, black circles). Each small circle represents an individual biological process and size is proportional to number of involved genes, while connecting edges represent shared genes between individual biological processes. (A) Gene ontology analysis showing enriched biological processes due to grandmaternal HDM exposure in F2 offspring at baseline (F0 HDM•F2 Base versus F0 VEH•F2 Base). (B) Gene ontology analysis showing effects of HDM sensitization alone on enriched biological processes in F2 offspring from vehicle-exposed F0 mice (F0 Veh•F2 Sens vs F0 VEH•F2 Base). (C) Gene ontology analysis showing potentiating effects of F0 grandmaternal HDM exposure on enriched biological processes after HDM sensitization in F2 mice (F0 HDM•F2 Sens vs F0 HDM•F2 Base). N=8 total for each group (n=4 males and n=4 females). <strong>Supplementary Figure 3. Differential methylation is linked to altered gene expression in airway epithelium of second-generation offspring. </strong>A subset of F2 airway epithelial genes identified during differentially-methylated pathway analyses were analyzed by real-time RT-PCR to assess for expression. Grandmaternal HDM exposure significantly increased Smad3 expression at baseline and potentiated its response after F2 HDM sensitization. Paxillin expression was also increased after HDM sensitization compared to non-sensitized F2 descendants from HDM exposed founders, although this difference did not meet statistical significance (<em>p</em> = 0.06). No differences were observed between Ets3, Cdh3, or Cacna1d expression between groups. *<em>p</em> &lt; 0.05. Delta CT values were analyzed with two-way ANOVA and Tukey’s post-test. N=8 total for each group (n=4 males “open circles” and n=4 females “closed circles”), except for F0 HDM•F2 Base, n=5. <strong>Supplementary Figure 4. Grandmaternal allergen exposure during pregnancy regulates vagal ganglia cytosine methylation in second-generation offspring at baseline and after allergen sensitization</strong>. Grandmaternal HDM sensitization during pregnancy resulted in 901 DMCs in F2 mice at baseline (F0 HDM•F2 Base vs F0 Veh•F2 Base). F2 HDM sensitization alone resulted in 1,467 DMCs (F0Veh•F2Sens vs F0 Veh•F2 Base). In contrast, HDM sensitization in F2 mice from HDM-sensitized F0 grandmaternal founders resulted in 1,787 total DMCs (F0 HDM•F2 Sens vs F0 HDM•F2 Base). N=8 total for each group (n=4 males and n=4 females). <strong>Supplementary Figure 5. Grandmaternal allergen sensitization during pregnancy modifies epigenetic signatures in second-generation offspring in a tissue-specific manner. </strong>Only rarely were differentially methylated genes shared between vagal ganglia neurons and airway epithelium in all groups. Specifically, neurons and epithelium in non-sensitized F2 mice from grandmaternal HDM-sensitized F0 founders shared only 15 genes in common (light blue circles). Neurons and epithelium in HDM-sensitized F2 mice from Vehicle-exposed grandmaternal founders shared only 9 genes in common (grey circles). Neurons and epithelium in HDM-sensitized F2 mice from HDM-sensitized grandmaternal F0 founders shared only 41 genes in common. N=8 total for each group (n=4 males and n=4 females). <strong>Supplementary Figure 6. Gene ontology analysis of biological process enrichment in vagal ganglia from second-generation offspring</strong>. Gene ontology analysis of DMCs in vagal ganglia of F2 mice identified enriched biological processes (individual small circles), which were then grouped based on functional annotation using Cytoscape Enrichment Map (larger, labeled, black circles). Each small circle represents an individual biological process and size is proportional to number of involved genes, while connecting edges represent shared genes between individual biological processes. (A) Gene ontology analysis showing effects of HDM sensitization on biological process enrichment in F2 offspring from vehicle-exposed F0 mice (F0 Veh•F2 Sens vs F0 VEH•F2 Base). (B) Gene ontology analysis showing potentiating effects of F0 grandmaternal HDM sensitization on biological process enrichment in F2 mice after HDM sensitization (F0 HDM•F2 Sens vs F0 HDM•F2 Base). N=8 total for each group (n=4 males and n=4 females). <strong>Supplemental Figure 7. Quantification of airway sensory nerve density using 3D confocal microscopy.</strong> Airway epithelial substance P-positive sensory nerves were analyzed in whole mount specimens using immunofluorescence and 3-dimensional confocal microscopy. Total nerve length and the number of branch points were quantified from nerve models generated based on substance P-positive voxels. Top) Representative images of substance P-expressing nerves (purple, top left) and nerve modeling (white filaments, top right). Bottom) Nerve length (left) and the number of branch points (right) were similar between untreated F2 mice from HDM exposed grandmothers and vehicle exposed grandmothers. N=6 both groups (n=3 males “open circles” and n=3 females “closed circles”). <strong>SUPPLEMENTAL TABLES</strong> <strong>Supplemental Table S1.</strong> Pathways identified by DMC analysis that were not present on DMR analysis in airway epithelium of second-generation F2 mice. All DMR-associated pathways were also identified by DMC analysis, as indicated in Table 1 in main text. <strong>Supplemental Table S2. </strong>Over-represented transcription factor binding sites in DMRs of airway epithelium and vagal ganglia from F2 mice before and after allergen sensitization. <strong>Supplemental Table S3. </strong>Pathway identified by DMC analysis in vagal ganglia of second-generation F2 mice. All DMR-associated pathways were also identified by DMC analysis, as indicated in Table 3 in main text. <strong>SUPPLEMENTARY DATA FILES</strong> <strong>Lebold&Cook_DMC_Gene List.xlsx. </strong> <strong>Lebold&Cook_DMR Gene List.xlsx.</strong> Differentially methylated cytosines (DMC) and regions (DMR) source data files.