Plasticity of gene expression in the nervous system by exposure to environmental odorants that inhibit HDACs
Abstract
Eukaryotes respond to secreted metabolites from the microbiome. However, little is known about the effects of exposure to volatiles emitted by microbes or in the environment that we are exposed to over longer durations. Using Drosophila melanogaster, we evaluated a yeast-emitted volatile, diacetyl, found at high levels around fermenting fruits where they spend long periods of time. Exposure to the diacetyl molecules in headspace alters gene expression in the antenna. In vitro experiments demonstrated that diacetyl and structurally related volatiles inhibited conserved histone deacetylases (HDACs), increased histone-H3K9 acetylation in human cells, and caused changes in gene expression in both Drosophila and mice. Diacetyl crosses the blood–brain barrier and exposure caused modulation of gene expression in the mouse brain, therefore showing potential as a neuro-therapeutic. Using two separate disease models previously known to be responsive to HDAC inhibitors, we evaluated the physiological effects of volatile exposure. Diacetyl exposure halted proliferation of a neuroblastoma cell line in culture. Exposure to diacetyl vapors slowed progression of neurodegeneration in a Drosophila model for Huntington’s disease. These changes strongly suggest that certain volatiles in the surroundings can have profound effects on histone acetylation, gene expression, and physiology in animals.
eLife assessment
This interesting and important work shows that diacetyl, a volatile organic compound released by yeast in fermenting fruit, can act as a histone deacetylase (HDAC) inhibitor and trigger wide changes in gene expression, together with suppression neurotoxicity in a Drosophila model of Huntington's disease. While the effects on gene expression changes and degenerative phenotypes are convincingly shown, further studies are required to determine whether and how olfactory sensory neurons and odorant receptors mediate the effects of diacetyl described by the authors.
https://doi.org/10.7554/eLife.86823.3.sa0Introduction
Multicellular organisms evolve in the presence of microbial volatiles and metabolites (McFall-Ngai et al., 2013; Dorrestein et al., 2014), often exposed for prolonged periods of time. Human-associated microbes have incredible metabolic diversity, encoding 100× as many genes as the human genome itself (McFall-Ngai et al., 2013). It is estimated that about half of the thousands of metabolites found in mammalian blood are either produced or altered by microbes (Wikoff et al., 2009). Microbial metabolites are an important mechanism by which the microbiome influences host immunity amongst other things (Dorrestein et al., 2014; Rooks and Garrett, 2016). Apart from microbial volatiles, many volatiles are also generated from the metabolism of a variety of food components, like triglycerides, sugars, and amino acids (Shibamoto, 2014). Some chemicals are also produced by microorganisms such as yeasts and lactic bacteria during fermentation in many foods and beverages (Shibamoto, 2014).
In the model system Drosophila melanogaster, prolonged exposure to CO2 showed an increase in volume of the CO2-sensing glomerulus in the antennal lobe, which lead to adaptations in fly behavior (Devaud et al., 2001; Sachse et al., 2007). Carbon dioxide is a common by-product of yeast fermentation in overripe fruits, the preferred food source. In order to study the physiological effects of long-term exposure to other volatiles emitted during fermentation, we exposed flies for 5 d to diacetyl, a buttery odorant made at high levels by yeast in fermenting fruits (1350 ppb in ripe banana), that is also found in many commonly consumed foods such as butter, yoghurt, wine, fruits, beer, and popcorn (Shibamoto, 2014; Martineau et al., 1995; Nykanen and Nykanen, 1991; Hughes and Baxter, 2001). Diacetyl was selected for two reasons: first, D. melanogaster use fermenting fruits as their favored food source, oviposition site, and spend long durations of time exposed to it. Second, diacetyl is known to be an inhibitory odorant of the CO2 receptor and therefore an opposite effect is expected on the V-glomerulus from the previous study (Turner and Ray, 2009). Follow-up experiments revealed that diacetyl acts as a histone deacetylase (HDAC) inhibitor in vitro. We found that volatiles structurally like diacetyl have similar HDAC inhibitory properties as well. Subsequent analyses with two different model systems – D. melanogaster and Mus musculus – showed gene expression changes in tissues caused by odor exposure from a distant source. We found expression of hundreds of genes to be modulated upon exposure to the volatile, as would be expected from an inhibitor of HDACs. HDACs are histone-modifying enzymes involved in the removal of acetyl groups from lysine residues and the remodeling of chromatin structure to regulate gene expression (Shahbazian and Grunstein, 2007; Gräff and Tsai, 2013). With large-scale roles in gene regulation, HDACs are promising targets in drug development for many diseases such as cancers and neurodegenerative disorders (Minucci and Pelicci, 2006; Kazantsev and Thompson, 2008; Chuang et al., 2009). Several classes of orally administered HDAC inhibitors have been found to attenuate the progression of certain cancers and neurodegenerative diseases, including Alzheimer’s disease and Huntington’s disease (Chuang et al., 2009). We determined that exposure to diacetyl volatiles substantially altered gene expression in the mouse brain, presumably by crossing the blood–brain barrier. Diacetyl exposure also slowed degeneration of photoreceptor cells in a Huntington’s disease model in Drosophila. Last, exposure to diacetyl caused dramatic reduction in proliferation of a neuroblastoma cell line in culture. Our discovery of a family of volatile odorants that also act as HDAC inhibitors highlights the need to understand how small molecules present in our environment interact with and alter eukaryotic gene expression.
Results
Activity-dependent expression modulation of chemosensory receptors
Previous studies demonstrated that prolonged exposure of ~5 d to elevated CO2 levels causes an increase in volume of the CO2-sensitive V-glomerulus in the antennal lobe of Drosophila, which lead to adaptations in fly behavior (Devaud et al., 2001; Sachse et al., 2007). To assess the effect of the CO2 response inhibitor diacetyl on the antennal lobe, flies were exposed to the odorant during the first 5 days immediately following eclosion. The CO2-detecting ab1C neuron was labeled using the promoters of the CO2 receptors, Gr63a-Gal4; UAS-mcd8GFP and Gr21a-Gal4; UAS-mcd8GFP, and the flies were exposed to headspace above a 1% (V/V) diacetyl solution in an air-tight container. The flies were then imaged for expression of GFP in the ab1C neurons of the antenna, as well as the V glomerulus, which receives axonal input from ab1C neurons in the antennal lobe (Jones et al., 2007). Surprisingly, the Gr63a promoter-driven GFP signal in the V glomerulus was completely abolished after 5 days of exposure to the odor for all brains imaged (Figure 1A). However, there was also a concomitant decrease in GFP signal in ab1C neurons in the antenna (Figure 1A and B). Similar results were seen for Gr21a-driven expression, where expression of GFP was lost in both the V-glomerulus and ab1C neurons (Figure 1—figure supplement 1A and B). Interestingly, after 5 d of recovery in clean air, the GFP expression was regained in the V-glomerulus as well as in ab1C neurons in the antenna, indicating that both Gr21a and Gr63a expression had been restored (Figure 1—figure supplement 2). Because olfactory neurons do not regenerate in the fruit fly (Berdnik et al., 2006; Sigg et al., 1997), these results suggest that changes in GFP levels are likely due to changes in the promoter expression and not due to neuronal cell death.
In order to test whether expression of Gr63a and Gr21a is being downregulated by diacetyl, we used quantitative-PCR for testing along with a few members of the Odorant receptor (Or) family genes. Apart from Gr63a and Gr21a, we also evaluated receptors where diacetyl has no activity (Or47a, Or88a) and the broadly expressed co-receptor (Or83b). Wild-type flies that were exposed to diacetyl for 6 d showed dramatic downregulation of all genes tested. However, after a 5 d subsequent recovery period, expression of all but one (Or88a) of the genes was recovered, including Gr63a (Figure 1—figure supplement 1C). This form of plasticity in gene expression is unexpected and therefore mechanistically interesting to investigate.
Diacetyl and structurally related odorants act as HDAC inhibitors in vitro
Diacetyl is structurally related to a soluble metabolite secreted by the gut microbiome and liver, β-hydroxybutyrate, which is known to inhibit zinc-dependent HDACs (Shimazu et al., 2013). It is also structurally related to another known HDAC inhibitor sodium butyrate. Diacetyl is present widely in nature as a pH-neutral fermentation product of microorganisms such as yeasts in over-ripe fruit, the favorite food for Drosophila, and also from lactic acid bacteria, and is found in many foods and beverages (Shibamoto, 2014). In humans, diacetyl is also produced by microbes on the skin, in the oral cavity, and is found in breath (Whiteson et al., 2014). Interestingly, this volatile is known to traverse the cell membrane (Krogerus and Gibson, 2013).
Soluble metabolites from the gut microbiome, such as short-chain fatty acids, are known to affect histone-modifying enzymes, gene expression, and physiology in various parts of the host (McFall-Ngai et al., 2013; Krautkramer et al., 2017). It is presumed that these dissolved compounds absorbed in the blood circulate to various parts of the body. Microbes in our skin microbiome, mouth microbiome, and environment also emit a diverse repertoire of volatile metabolites that can traverse much greater distances away from the source (James et al., 2013). To examine whether diacetyl can directly modulate HDACs like β-hydroxybutyrate, we performed in vitro acetylation assays with purified recombinant HDACs. We found that, indeed, diacetyl inhibited all four purified human class I HDACs (HDAC1, 2, 3, and 8) and the one class II HDAC (HDAC6) tested. The inhibition occurred in a dose-dependent manner. The IC50 values for HDAC1, 2, 3, 8, and 6 were 7.3, 23.1, 7.5, 14.3, and 24.5 mM, respectively (Figure 1C). The levels of inhibition as observed from IC50 values for HDAC1 and HDAC3 were comparable to those of β-hydroxybutyrate whose IC50 values were previously reported as 5.3 and 2.4 mM, respectively (Shimazu et al., 2013).
To test other odorants for their ability to inhibit HDACs, we identified a list of volatile odorants structurally similar to diacetyl using the similarity search feature on PubChem (Figure 1D) and tested them on four purified human HDACs, from class I (HDAC1 and HDAC3) and class II (HDAC4 and HDAC6). Each compound was tested at 15 mM, the concentration at which diacetyl inhibits HDAC1, 3, and 6 with more than 50% efficacy (Figure 1C and D). All compounds tested reduced the activity of at least one of the HDACs tested (Figure 1D). Ethyl pyruvate, methyl pyruvate, 2,3-hexanedione, 2,3-heptandione, 2,3-pentandione inhibited HDAC1, 3, and 6 with more than 70% efficacy. Allyl butyrate also inhibited those HDACs but with lower efficacy. 1-Acetoxyacetone and 2,3-butanediol inhibited HDAC6 more than 70% but inhibitory effects against HDAC1 and 4 were lower. Moreover, propyl formate strongly inhibited HDAC 6 but did not inhibit HDAC1, 3, and 4. Taken together, these results indicate that microbial volatiles and structurally related odorants can inhibit human HDACs, and that each compound has a specific HDAC inhibitory spectrum.
Microbial volatile increases histone acetylation in human cell nuclei
In order to directly examine histone acetylation at the cellular level, we used human HEK293T cells, which offer a tractable system to prepare nuclear extracts, and picked the volatile diacetyl to test in more detail. We exposed the cells to different doses of diacetyl for 2 or 6 hr and monitored histone acetylation levels within nuclear extracts in western blot assays. Compared to the mock treatment, 10 mM diacetyl significantly increased H3K9 acetylation levels within 2 hr of treatment, whereas the acetylation levels of H3K14 and H4K5 were not affected (Figure 2A). This specificity for an increase of the H3K9 mark is consistent with previous observations for β-hydroxybutyrate. Preferential acetylation of H3K9 was also observed previously in HEK293T cells with treatment of the structurally related β-hydroxybutyrate (Shimazu et al., 2013). After 6 hr of treatment, the H3K9 acetylation induced by 10 mM diacetyl was further increased (Figure 2B). The increase in H3K9 acetylation with diacetyl treatment was dependent both on the duration of exposure and concentration of the odorant.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/86823/elife-86823-fig2-data1-v1.zip
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/86823/elife-86823-fig2-data2-v1.zip
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/86823/elife-86823-fig2-data3-v1.zip
-
Figure 2—source data 4
- https://cdn.elifesciences.org/articles/86823/elife-86823-fig2-data4-v1.zip
Organisms are constantly exposed to volatiles commonly found in their food and environment for prolonged periods of time. In order to test the effect of a lower concentration of the volatiles, we selected a 5 d exposure time at a 100-fold lower concentration. When we treated HEK293T cells with this lower dose of diacetyl (100 μM), H3K9 acetylation level increased after 96 hr of exposure and reached significantly higher levels than control after 120 hr (Figure 2C and D). These results demonstrate that even prolonged exposure to low levels of diacetyl can greatly impact the epigenetic landscape inside the nucleus. More importantly, a 5 d exposure was sufficient to alter the epigenetic state of nuclei at concentrations that are present in some food sources. Taken together, these results demonstrate that diacetyl can act as an HDAC inhibitor, with the potential of causing broad modulations of gene expression and histone acetylation in cells.
Transcriptional response to odor exposure is conserved in invertebrates and vertebrates
We next performed in vivo experiments to determine whether eukaryotes alter gene expression in response to diacetyl exposure, as would be expected with HDAC inhibition. Epigenetic changes occur slowly over days, and to test its effects on gene expression we used two model systems for eukaryotes: invertebrates (D. melanogaster) and vertebrates (M. musculus). We placed Drosophila adult males in vials and housed these vials within a closed container exposed to headspace from a 1% diacetyl solution (V/V in paraffin oil) for 5 d, similar to a previous long-term odor-exposure study (Sachse et al., 2007; Figure 3A). The transcriptome of the primary olfactory organ, the antenna, was compared with that of the control group of age-matched flies that were exposed to the solvent paraffin oil (PO). The antennal transcriptional profile of diacetyl-exposed flies showed substantial changes in gene expression when compared to the solvent control (Figure 3B). We identified 1234 differentially expressed genes (DEGs) (false discovery rate [FDR] < 0.05) in the antennal transcriptome of diacetyl-exposed flies compared to control animals. Of these, 645 genes were significantly upregulated (log2 fold-change > 1; red dots in Figure 3B) and 589 genes were significantly downregulated (log2 fold-change < –1; blue dots in Figure 3B). A broad range of genes were significantly altered, with several biological process GO terms significantly enriched in the upregulated gene list including ‘response to biotic stimulus’ (p<3 × 10–7), ‘response to bacterium’ (p<3 × 10–7), and ‘defense response’ (p<3 × 10–7). In the downregulated gene set, the GO term most enriched was ‘sensory perception of chemical stimulus’ (p<6 × 10–72). The GO molecular functions that were most common in the upregulated sets included enzymes like hydrolases and oxidoreductases, as well as nucleic acid binding genes (Figure 3C). Downregulated genes include receptors, hydrolases, and transporters (Figure 3C).
The first tissue that would have contact with airborne volatiles in terrestrial vertebrates are the airways and lungs. We therefore performed transcriptome analyses on lung tissue of mice exposed to diacetyl headspace at different doses for a period of 5 d, as was done in Drosophila (Figure 3D). The dose of the volatile in the chamber was measured in parallel experiments by tracking volatile loss and airflow and determined to range between a peak 1920 ppm at the end of hour 2, tapering off to 28 ppm by hour 120 (Figure 3—figure supplement 1B). Indeed, expression of a substantial number of genes was modulated in the diacetyl-exposed lungs compared to the control. The changes were dose-dependent and more pronounced in mice exposed to headspace from a 1% diacetyl source compared to those exposed similarly to 0.1% (v/v) diacetyl source (Figure 3E, Figure 3—figure supplement 1C). Among these diverse sets of regulated genes in the lung tissue, a significant overlap was found between 1% (v/v) and 0.1% (v/v) headspace for both upregulated genes (p=3 × 10–3) and downregulated genes (p=6 × 10–3, Figure 3F), further supporting a dose-dependent effect of diacetyl on gene expression in the mouse lung.
The GO enrichment analysis of the lung transcriptome from mice exposed to diacetyl revealed several interesting sets of genes that were significantly altered (Figure 3G and H). Among these genes is the SARS-CoV-2 entry receptor ortholog, ACE2. The level of this gene in the lungs of diacetyl-exposed mice showed a reduction in expression levels (log2 fold-change = - 0.62, FDR = 2.02 × 10–5). The level of ACE2 in the nasal tissue, lungs, and other organs plays a critical role in infection by the SARS-CoV-2 virus, suggesting that epigenetic mechanisms may contribute to expression. While we do not know whether expression of ACE2 in humans is regulated similarly, these results highlight the importance of understanding volatile-induced modulation of gene expression in the lungs and other tissues.
HDAC family members are highly conserved across eukaryotes, spanning both animal and insects. Volatile microbial metabolites have been present throughout the evolution of animals, insects, and plants, and have potential as signals for multidomain communication because they can travel wide ranges; plants detect root infections and time their immune response using volatiles (Effantin et al., 2011), and microbial volatiles elicit olfactory behaviors in insects and nematodes. Our results indicate that a volatile compound, capable of directly inhibiting HDACs and altering acetylation levels, can cause significant changes in gene expression in a variety of organisms following exposure from a distance, through the air.
DEGs overlap those for known HDAC inhibitors and upregulation is partially reversible
One of the hallmarks of epigenetic regulators which makes them attractive candidates for therapeutics is plasticity. Upon removal of the inhibitors, we expect that some of the changes in gene expression are reversible. To test for reversibility, we performed a recovery experiment following diacetyl exposure using the Drosophila model. We maintained 5 d diacetyl-exposed flies in clean air for five additional days (Figure 4A). In parallel, we performed age-matched mock experiments with paraffin oil solvent exposure alone. A large number of genes were downregulated following the recovery in comparison to the 10-day-old flies in the mock condition (Figure 4B). Interestingly, there was a significant overlap of these downregulated genes with the set that was upregulated in the diacetyl treatment (Figure 4C). These results suggest that the effects of HDAC inhibitory odorant exposure are not permanent but dynamic, and removal of the odorant leads to subsequent changes in gene expression of the upregulated set.
Drugs that inhibit HDACs are in various stages of approval as treatments for several diseases. The approved drugs include sodium butyrate and valproic acid. To compare the gene-regulatory effects of these drugs to diacetyl, we performed RNA-seq after raising the flies on food containing sodium butyrate (SB) or valproic acid (VA) (Steffan et al., 2001) and compared gene expression in the antennae of flies raised on untreated food for 5 d. We next compared the upregulated gene profiles following each treatment to the one induced by exposure to diacetyl. As expected, feeding SB and VA induced significant changes in expression levels of several genes (Figure 4D and E). Interestingly, we found that 113 of diacetyl upregulated genes were also upregulated in either SB, VA, or both treatment conditions (Figure 4F). Pairwise statistical analysis of each gene set revealed a significant overlap of diacetyl-induced genes with SB-induced genes (p=6 × 10–11) and with VA-induced genes (p=2 × 10–65) (Figure 4F). There was, as expected, also a significant overlap between SB- and VA- induced genes (p=1 × 10–52) (Figure 4F). This highly significant overlap among upregulated genes lends further support to our model that diacetyl vapors act as an HDAC inhibitor in vivo. As expected, each of the three treatments also modulated a substantial number of unique genes (Figure 4G and H), suggesting that differences in oral vs. vapor delivery, molecular structure, and inhibition profile across the repertoire of HDACs may contribute to differences in gene regulation. These results suggest that some volatile odorants, or the microbes that emit them, may impart therapeutic effects like other HDAC inhibitors from a distance, through the air.
Volatile diacetyl alters gene expression in the mouse brain
One of the categories of challenging diseases that HDAC inhibitors are being tested for drugs are neurodegenerative diseases (Fischer, 2010), where volatile odorants may provide a new class of inhaled therapeutics. A major design challenge for nervous system drugs is the ability to cross the blood–brain barrier. Due to the volatility and small size of volatile odorants, there is a possibility that they could diffuse through the intranasal route to the brain directly (Chauhan and Chauhan, 2015). We cannot tell a priori whether odorants like diacetyl can travel across the nasal membrane to the brain. In order to test whether cells in the brain of mice respond when presented with diacetyl vapors from a distance, we used gene expression as the final readout to measure it. We performed RNA-seq experiments on mice exposed only to aroma of diacetyl in the air for 5 d as before. Littermate controls were exposed in a similar manner to the solvent (PO) in the air. Interestingly, several genes were differentially expressed upon exposure to headspace from 0.1% (v/v) diacetyl passaged through the air in the holding cage where it is further diluted (49 upregulated, 32 downregulated, |log2 fold-change|>1, FDR < 0.05) or to headspace from 1% (v/v) diacetyl diluted in the holding cage (748 upregulated, 1031 downregulated, |log2 fold-change|> 1, FDR < 0.05) (Figure 5B and C). GO analysis of the regulated genes in the exposed mouse brain transcriptome revealed several interesting sets of genes were significantly altered in each set (Figure 5E). These results indicate that a volatile odorant like diacetyl can reach the brain in mammals and presumably alter gene expression in neuronal cells. Although the overall DEG sets were different across the lungs and brain, there was a statistically significant overlap, particularly in the 1% exposure, as would be expected from using the same mechanism of action. There was a highly significant overlap (p=2 × 10–101) between genes upregulated in the brain compared with the lungs (248 of 748), as well as for downregulated genes (p=3 × 10–163; 477 of 1031 genes) (Figure 5D).
Diacetyl prevents proliferation of neuroblastoma cells
Among the significantly downregulated genes was MYCN, which is known to show increased expression and plays a central role in many types of cancers, suggesting that epigenetic regulation may affect tumor cells. Broad-spectrum HDAC inhibitors have been approved by the FDA for treatment of a variety of cancers. To test the exciting possibility that a volatile HDAC inhibitor, which affects gene expression in the brain, would have an effect against cancers as well, we performed a cell-based assay. We picked neuroblastoma as a candidate since it has a strong dependence on MYCN and develop from fetal adrenal neuroblasts (Jansky et al., 2021). Several genetic and ‘omics’ studies have identified genes that are upregulated in neuroblastomas, some of which are considered to be key causative factors (Pugh et al., 2013). We evaluated the DEGs from the diacetyl-exposed mouse brain transcriptome for several of these key genes. Three of the six key genes identified across multiple studies were significantly altered, two of which were significantly downregulated, MYCN (p=2.06E-13, FDR = 1.76E-12) and Alk (p=3.3E-5, FDR = 9.48E-5) (Figure 5F). An unbiased computational study evaluated omics datasets and identified four genes that were the best predictors of clinical outcome (Zhong et al., 2018), all of which are significantly different in the brains of the diacetyl-exposed mice, strikingly three of which were significantly downregulated, NCAN (p=1.33E-16, FDR = 1.59E-15), STK33 (p=4.6E-7, FDR = 1.79E-6), and ERCC6L2 (p=4.57E-5, FDR = 0.00012) (Figure 5F). These findings are extremely provocative, given that the downregulated genes, particularly MYCN, play a role in neuroblastomas and in several other cancers.
We tested diacetyl for the ability to inhibit proliferation of an SH-SY5Y neuroblastoma cell line from a metastatic bone tumor. A significant reduction in cell count was seen in diacetyl treatments in the SH-SY5Y cancer cell lines (Figure 5G). To validate the initial result, we retested diacetyl at two concentrations across three different cancer cell lines. Diacetyl showed a significant reduction in cell numbers for the neuroblastoma cell line specifically, and not any other cell lines, relative to the control (Figure 5G). These results demonstrate that the effect is specific to the SH-SY5Y cancer cell line and unlikely to be caused by a nonspecific general toxic effect. It also suggests that the potential epigenetic regulation of cancer-promoting genes may be caused by exposure to volatiles like diacetyl in the environment.
Testing therapeutic effect in Huntington’s neurodegeneration model
Since diacetyl can affect changes in a mammalian brain by crossing the blood–brain barrier, we wanted to test whether it can impart therapeutic effects in disorders of the brain like neurodegeneration for which HDAC inhibitors have been proposed as therapeutics. In order to test diacetyl in the context of neurodegenerative disease, we used the well-established Drosophila model of Huntington’s disease, which has previously been used to demonstrate the efficacy of sodium butyrate in slowing neurodegeneration (Steffan et al., 2001; Marsh and Thompson, 2006). The human Huntingtin protein with expanded poly-Q repeats is expressed in the neurons of the compound eye in Drosophila, causing progressive degeneration of the seven visible photoreceptor rhabdomere cells in each ommatidium (Jackson et al., 1998). We selected this model in particular since polyglutamine disorders are well suited for targeting by inhibitors of HDAC1 and HDAC3 (Thomas, 2014) such as diacetyl (Figure 1B). Moreover, previous studies have shown that orally administered HDAC inhibitors such as sodium butyrate and SAHA can significantly reduce photoreceptor degeneration in this model (Steffan et al., 2001). When the transgenic flies expressing two copies of the human Huntingtin with poly-Q repeats (HTTQ120) under control of the eye-specific GMR promoter were raised at 18°C, the number of rhabdomeres in each ommatidium was similar to that normally observed in wild-type flies immediately post-eclosion (day 1, Figure 6A–C). When these HTTQ120 flies were moved to 25°C following eclosion (Figure 6A), they showed dramatic degeneration of rhabdomeres over a period of 10 d (Figure 6B and D–G). The mean number of rhabdomeres was reduced from 7 to ~1 by day 10. The dose of the volatile in the chamber was measured in parallel experiments by tracking volatile loss and airflow (0.5 L/min) and determined to range from 260 to 95 ppm, an average of 125 ppm at the end of hour 1, tapering off to 0 ppm by hour 24. The odorant was replaced each day, leading to a dosage profile that would seem like a once-a-day treatment regimen (Figure 3—figure supplement 1A).
Remarkably, when the Huntingtin (HTTQ120)-expressing flies were exposed immediately after eclosion to volatile headspace of 1% diacetyl (Figure 6A), they showed a substantial (~50%) reduction of rhabdomere loss (Figure 6B and D–G). The majority of ommatidia retained 6–7 rhabdomeres at day 5 (Figure 6D and F). Even after 10 d, the majority of the ommatidia still had 2–3 rhabdomeres left in the odor-exposed flies, while the solvent controls had only ~1–2 rhabdomeres (Figure 6E and G). These results demonstrate that exposure to volatile diacetyl in air slows down the photoreceptor degeneration caused in these Huntington’s disease model flies. Thus, odorants like diacetyl may have potential as a prophylactic against this neurodegenerative disorder. Therefore, the volatile odorant can both significantly alter gene expression from a distance as well as cause a phenotypic change in the fly with a human disease gene.
Discussion
In this study, we discovered that eukaryotic cells could alter gene expression in response to an odorant in a manner independent of traditional neuronal activity-induced pathways. Members of a conserved family of HDACs detected the concentration of diacetyl differentially. Diacetyl exposure increased levels of H3K9 acetylation in the nucleus in a dose-dependent manner and led to changes in gene expression. The odorant diacetyl occurs naturally and is generated from the metabolism of a variety of food components, including triglycerides, sugars, and amino acids (Shibamoto, 2014). Moreover, the chemical is produced by the activity of microorganisms such as yeasts and lactic bacteria during fermentation in many foods and beverages (Shibamoto, 2014; Figure 3—source data 1). Although beyond the scope of this study, the discovery that cells can alter gene expression upon prolonged exposure to a common naturally occurring odorant raises important questions about the potential physiological consequences of the expression changes. Given our repeated exposure to particular flavors and fragrances, the findings outlined here highlight a new consideration for evaluating the safety of certain volatile chemicals that can cross the cell membrane.
In mammals, upregulation of circulating β-hydroxybutyrate during fasting or calorie restriction induces changes in the expression of a set of genes (Shimazu et al., 2013). The acetylation mark specificity and IC50 of β-hydroxybutyrate is similar to that of diacetyl. Any beneficial physiological effects of odor detection at lower odor concentrations would also be countered by studies of potential risks at higher concentrations. In fact, the deleterious effect of exposure to high levels to diacetyl causing bronchiolitis obliterans, or ‘popcorn lung,’ and its toxicity in cultured cells is already known (More et al., 2012). Even so it is present in several foods we eat and is on the Generally Regarded as Safe list of Flavor and Extract Manufacturers Association for use as a flavoring ingredient at low concentrations. While the specific mechanism of ‘popcorn lung’ is unknown, it has previously been proposed that carbonyl groups of diacetyl can react to modify amino acid side chains on proteins or generate reactive dicarbonyl and reactive oxygen species, leading to excessive cytokine production and inflammation (Starek-Swiechowicz and Starek, 2014). Our results raise the possibility that the molecular mechanisms underlying this disorder could also be partially attributed to the unusually high HDAC-mediated genetic response to diacetyl by cells in the lungs. Several studies using rodent animal models have shown toxicity of high levels of diacetyl in the respiratory tracts including the nose and lungs (Morgan et al., 2008; Palmer et al., 2011). The toxic effect of higher dosages has been observed for many types of HDAC inhibitors. For example, in one Drosophila study, feeding of a high dosage of 4-phenylbutyrate reduces survival rate, while a lower dosage extends longevity (Kang et al., 2002). Another concern is that one of the highest known levels of diacetyl exposure to humans is through smoking cigarettes, where exposure from even a single cigarette (250–361 ppm) is high (Rigler and Longo, 2010; Pierce et al., 2014): Salem (128 μg/cigarette), Camel (128.1 μg/cigarette), and Carlton (45.6 μg/cigarette) (Pierce et al., 2014). Different types of diets and lifestyle choices have been known to have effects on health and well-being due to differences in balance of different health-promoting ingredients, yet little is known about the role of such flavors, fragrances, and metabolites acting on epigenetic marks and altering gene expression. Our findings raise the possibility that diacetyl and structurally related odors in the environment may affect gene expression directly through absorption/transport into these cells. In mammals, the epithelium-lined lungs and nasal/oral pathway is directly exposed to odorants, presenting an opportunity to affect gene expression and physiology.
Eukaryotes primarily detect volatile odorants in the environment using olfactory neurons that express a variety of transmembrane receptors, a family of genes that have independently evolved multiple times in different phyla. Examples of these unrelated genes include the ionotropic 7-transmembrane (TM) odorant receptor (Or) family, which is insect-specific; 3-TM ionotropic receptors, which are present across most arthropods; the nematode-divergent 7-TM GPCRs belonging to the str, sra, srg, srw, srz, srbc, srsx, and srr families; the mammalian 7-TM GPCR olfactory receptor (OR) family; and the trace amine-associated receptor family (Buck and Axel, 1991; Clyne et al., 1999; Troemel et al., 1995; Benton et al., 2009; Vosshall et al., 1999). The activation of these receptors induces neuronal action potentials, and this information is conveyed to higher brain centers where olfactory perception is generated (Knaden and Hansson, 2014). This second response pathway that we have discovered is much slower and has an enzymatic mechanism that draws parallels with the multiple pathways in place to detect light. In the specialized cells of the eye, detection of light occurs via rhodopsins (7-transmembrane GPCRs) and leads to neuronal activity and behavioral responses. However, other cells, including in plants, are also able to respond to changes in light intensity of certain wavelengths using the ancient, conserved cryptochrome proteins that are related to photolyases (Fogle et al., 2015). This suggests that vital sensory modalities can be detected via multiple pathways. Our observations raise the question whether HDACs represent a conserved detection mechanism for odorants like diacetyl that link odor detection to a specific gene expression response. HDAC proteins, however, are an ancient family of genes that predate even histones themselves (Gregoretti et al., 2004; Postberg et al., 2010). It is conceivable that odor detection mechanisms that emerged in ancient forms may not have involved specialized transmembrane receptors, or neurons, or resulted in rapid behavioral movements.
Our data suggests that an HDAC inhibitory odorant also affects gene expression in the lungs and brain of mice, including significantly decreasing the expression level of the ACE2 receptor, whose ortholog in humans is the entry receptor for the SARS-CoV-2 virus that causes COVID-19. The applied impact for the future could be far-reaching as HDACs are known to affect a wide variety of pathways in eukaryotes, and inhibitors are candidates for novel inhaled therapeutics as treatment of many conditions such as COVID-19, cancers, neurodegeneration, and a variety of infectious diseases (Witt and Lindemann, 2009; Wang et al., 2013; Wang et al., 2012; Yeung et al., 2012; Park et al., 2012; Cai et al., 2010; Shankar and Srivastava, 2008; Platta et al., 2007; Komatsu et al., 2006).
Our studies also show that one of the volatiles is protective against the neurodegenerative effects of Huntington’s disease model and neuroblastoma. Exposure to diacetyl volatiles in the fly model of Huntington’s disease reduces cell degeneration, as has been previously observed with orally administered HDAC inhibitors like sodium butyrate and SAHA in this genetic model (Steffan et al., 2001). Previous studies indicate that the inhibition of HDACs counters the acetyltransferase inhibitory activity of the polyglutamine-domain of the Htt protein, which binds to p300, P/CAF, and CBP (Steffan et al., 2001).
This newly discovered activity in this class of volatile chemicals will also allow us to understand how their presence in our environment, diet, and microbiome shapes the epigenetic landscape in the lungs and affects inflammation and immune responses. Our approach to investigate consequences in model systems allows for a foundation for future research and allow identification of key molecular biomarkers that can be associated with observation of epigenetic effects, toxic or not, in individuals. The public health impact of these findings can be far-reaching, both from the perspective of toxicity and therapeutics. Investigating the health safety concerns of exposure to such odorants at high doses or prolonged periods is essential to understand potential for damage. Conversely, at nontoxic concentrations, HDAC inhibitors can be used for the treatment of many debilitating conditions in the lungs such as COPD, inflammation, and COVID-19 infection (Shakespear et al., 2011).
Taken together, our discovery that cells can alter gene expression in response to an odor that acts as an HDAC inhibitor and reprogram gene expression promises the pursuit of new types of odor-based therapeutics already approved for human consumption as prophylactics for a myriad of diseases. Simultaneously, they raise the concern of a new type of environmental agent that can reprogram gene expression and have widespread effects that are yet unknown.
Materials and methods
Drosophila stocks and manipulations
Request a detailed protocolFly stocks were maintained on conventional fly food under a 12 hr light:12 hr dark cycle at 18°C or 25°C. The fly strain of w1118 backcrossed five times to Canton-S (wCS) was used in all the Drosophila transcriptome experiments. P{GMR-HTT.Q120}2.4 (Bloomington # 8533) were used for neurodegeneration experiments. For imaging Gr63a-Gal4; UAS-mcd8:GFP, Gr21a-Gal4; UAS-mcd8:GFP stocks were obtained as a kind gift from the Carlson lab, Yale.
Odor exposure experiments for imaging
Request a detailed protocolNewly eclosed male wCS flies were sorted in groups of 30 into vials with standard cornmeal medium. Vials were closed with two overlapping 3″ × 3″ polypropylene mesh squares and tied with cotton twine. Vials were used within 6 hr. Four vials were placed in a 1000 ml Nalgene straight-sided jar with a 10 ml glass beaker with 10 ml of diacetyl 10–2 in paraffin oil. Jars were closed for exposure periods ranging from 2 to 6 d. In recovery experiments, vials were removed from the jars and were placed in a 25°C incubator for the remainder of the recovery period. Flies were anesthetized on ice, sacrificed in ethanol, and then immediately put into 1× phosphate-buffered saline with 0.3% Triton-X (PTX). Brains or antenna were put in 4% paraformaldehyde in PTX and incubated for 30 min while rotating at 25°C. Samples were washed five times in PTX and blocked in 5% natural goat serum in PTX for 1 hr while rotating at 25°C. Samples were then incubated in primary antibody with 5% goat serum in PTX, mouse-nc82 (1:20) (Development Studies Hybridoma Bank, University of Iowa) and rabbit-antiGFP (1:150) (Invitrogen) in PTX for 48 hr in 4°C while rotating. After washing five times in PTX, samples were incubated in secondary antibody with 5% goat serum in PTX, rabbit-anti-Alexa488 (1:400) (Invitrogen), and mouse-anti-Alexa568 (1:400) (Invitrogen) for 48 hr in 4°C. Samples are washed five times in PTX and stored in 70% glycerol in PTX. Images were taken using a Ziess 510 laser scanning confocal microscope, and image analysis was done using ImageJ software.
QRT-PCR
Request a detailed protocolDrosophila heads (with antenna) were collected on liquid nitrogen and stored at –80°C until processing. A total of 300 heads were collected for each sample and ground on liquid nitrogen. mRNA was extracted using a PolyATract kit (Promega) and was reverse transcribed using Superscript III (Invitrogen) according to the manufacturer’s protocol. The quality of the cDNA was assessed on a 1.5% agarose gel. cDNA was added to individual reactions of SYBR green master mix (Bio-Rad) and run on a Bio-Rad-MyQ thermocycler. The program began with a single cycle for 3 min at 95°C, followed by 40 cycles of 15 s at 95°C, 30 s at 58°C and 30 s at 72°C. Afterward, the PCR products were heated to 95°C for 1 min and cooled to 55°C for 1 min to measure the dissociation curves. The efficiency of each primer set was first validated by constructing a standard curve and three 10× serial dilutions of first-strand cDNA. For each serial dilution, the CT value was plotted against the log (template dilution) and the slope and r2 value of each regression line calculated. Expression of each Or or Gr gene was assessed in triplicate. Dissociation curves were used to assess the purity of the PCR reactions. Or and Gr transcript levels were normalized by using the transcript levels of ribosomal protein 49. Expression levels were calculated using the Pfaffl equation.
In general, primers were designed by using coding sequences close to the 3′ end of the gene, and where possible, primers spanned an intron.
Odor exposure protocol for transcriptome analysis
Request a detailed protocolFlies were exposed to diacetyl (B85307, Sigma-Aldrich, St. Louis, MO) by placing them in vials in a cylindrical closed container (112 mm diameter × 151 mm height) along with an odor-containing glass vial. The odorant was dissolved in 5 ml paraffin oil at 1% dilution. For a given exposure protocol, two groups of flies were prepared: those exposed to 1% diacetyl headspace and those exposed to paraffin oil headspace alone (control flies). Adult male flies aged 1 d were transferred to fly vials containing fresh medium and put into the container with the odor vial. At the end of the fifth day of exposure, flies were collected, and their antennae were dissected for RNA extraction. All treatments and experiments were performed at room temperature. For the recovery experiment, flies were transferred to a container with a glass vial of paraffin oil after 5 d of diacetyl exposure. At the end of the fifth day of recovery, flies were collected, and their antennae were dissected for total RNA extraction. The second and third antennal segments from 40 to 60 male flies after treatment were carefully hand-dissected from the head and collected in 1.5 ml microfuge tubes kept cold in liquid nitrogen. Antennae were mechanically crushed with disposable RNAse-free plastic pestles, and total RNA was isolated using a Trizol-based protocol. cDNA libraries were prepared from total RNA using the Illumina TruSeq RNA Sample Preparation Kit (v2) and 50 bps single- and paired-end sequencing was done using the HighSeq2000. Two biological replicates were sequenced for each condition, with an average of 27 million reads/ replicate, and with an average of 84% mapped.
Two-month-old C57BL/6 male mice (2–3 for each condition in a single cage) were continually exposed to air flowing over headspace of paraffin oil (solvent control), 0.1%, or 1% diacetyl over a period of 5 d, then euthanized for collecting the lung and brain tissues and processing for mRNA isolation. All protocols for animal use and euthanasia were approved by the Institutional Animal Care and Use Committee (#20150028) and were in accordance with the provisions established by the Animal Welfare Act and the Public Health Services (PHS) Policy. In the transcriptome analysis, two replicates were performed for each condition, with an average of 123,687,411 reads/replicate, with an average of 88% mapped. Multiplexed libraries were made from total RNA input using the Illumina TruSeq RNA sample preparation kit (v2) and 50 bps single-end sequencing was done using the NextSeq500.
Calculation ppm in air based on weight loss of odorant over time
Request a detailed protocolThe odor solution in PO in each setup (Drosophila and mouse) was weighed using a sensitive balance. PO is nonvolatile and therefore weight loss is attributable to diacetyl volatilization. Volatile Headspace exposure (ppm) is calculated through odor cartridge weight change from pure odor compound volatizing out (Inamdar et al., 2012).
where ppm is the concentration of tested VOC in parts per million (v/v), w is the weight of the tested chemical in grams, MW is the molecular weight of the tested VOC (g/mole), v is the total volume of the respective exposure chambers, and Vm is the molecular volume under the tested condition.
where 24.45 is the gram molecular volume, under standard lab conditions of 760 mmHg, 25°C; P is the ambient pressure, in mmHg; and t is the ambient temperature, in°C. Odor cartridges, consisting of varying odor volume/volume dilutions, were weighed periodically under standard ambient temperature and pressure. Under these conditions, the pure volatile odor compounds volatize out of the odor cartridge and the evaporation rates were used for our exposure calculations. Based on ideal gas law and experimental conditions, ppm is calculated based on volatized odor compound and experimental chamber volume based on airflow.
HDAC inhibitor treatment protocol for transcriptome analysis
Request a detailed protocolSodium butyrate (B5887, Sigma-Aldrich) or valproic acid (P4543, Sigma-Aldrich) was dissolved in normal fly food medium at the final concentration of 10 mM. Three groups of flies were prepared: those treated with one of the HDAC inhibitors and those without HDAC inhibitor treatment (control flies). Adult flies aged 1 d were transferred to fly vials containing medium with or without a HDAC inhibitor. At the end of the fifth day of treatment, flies were collected, and their antennae were dissected for RNA extraction. All treatments and experiments were performed at room temperature. Two biological replicates with 60 flies/replicate were performed for each condition, with an average of 23 million reads/replicate, and with an average of 92% mapped. Multiplexed libraries were made from total RNA input using the Illumina TruSeq RNA sample preparation kit (v2) and 50 bps paired-end sequencing was done using the HighSeq2000.
Bioinformatic analysis of RNA-seq experiments
Request a detailed protocolReads aligned to the latest release of each of the genomes used (dm6 for the Drosophila genome, GRCm38 for the Mus musculus genome) and quantified with kallisto (version 0.43.1) (Bray et al., 2016) Only libraries for which we obtained >75% alignment were used for downstream analysis. Transcript counts were summarized to gene-level using tximport package (version 1.4.0) (Soneson et al., 2015). For any instances of detected batch effects, we removed unwanted variation using RuvR in the RuvSeq package (version 1.10.0) (Risso et al., 2014). DEG analysis was performed with the edgeR package (version 3.18.1) (Robinson et al., 2010) using low count filtering (cpm > 0.5) and TMM normalization. Protein classification analysis was performed with PANTHER (version 13.1) (Mi et al., 2017). All significance analyses of gene overlap were done using the GeneOverlap package in R package (version 1.14.0). GO enrichment analysis was performed using clusterProfiler (version 3.6.0) (Yu et al., 2012).
HDAC activity assays
Request a detailed protocolHDAC activity of class I HDACs (HDAC1, 2, 3, and 8) was measured with the fluorometric HDAC Activity Assay kit: HDAC1 (10011563, Cayman Chemical, Ann Arbor, MI), HDAC2, HDAC 3, and HDAC 8 (50062, 50073 and 50068, BPS Bioscience, San Diego, CA), according to the manufacturer’s instructions. HDAC activity of class II HDACs (HDAC4 and HDAC6) were measured with the fluorometric HDAC Activity Assay kit: HDAC4, HDAC 4 (50064 and 50076, BPS Bioscience), according to the manufacturer’s instructions. HDAC activity in these 96-well format kits was directly correlated with fluorescence produced by the enzyme and substrate interaction. Each test compound was tested in duplicates, minimum. Due to interference, concentrations of >30 mM were not used in the IC50 calculation. To account for any potential absorbance or fluorescence contributed by the test chemicals interacting with the assay substrate, all baseline background measurements of test compound + substrate control wells (everything, but no HDAC enzyme) were performed and subtracted against the respective test compound wells. It is fair to note that no significant changes in background fluorescence were seen with tested compounds in comparison to control background fluorescence. HDAC activity was calculated based on averaged maximum fluorescence from negative control wells and the averaged fluorescence in test compound wells.
HDAC % inhibition was directly calculated from obtained HDAC activity %.
Cell culture and treatment
Request a detailed protocolHuman embryonic kidney 293 (HEK293) cells were grown in 100 mm cell culture dishes with Dulbecco’s modified Eagle’s medium (DMEM) (10-013, Corning, Manassas, VA), supplemented with 10% fetal bovine serum (FBS) (26140-079, Gibco, Carlsbad, CA) at 37°C with 5% CO2. Cells that were ~80% confluent were treated with freshly prepared medium supplemented with diacetyl at the concentrations indicated. The cells for mock controls were handled in the same manner without adding diacetyl to the medium. To prevent diffusion of diacetyl odor from the treatment dishes to the ones of mock control, the cell culture dishes in different conditions were cultured in separate CO2 incubators.
Cell proliferation assays were performed to assess the effects of compounds against cancer cell lines. Cells were maintained in a 37°C incubator with 5% CO2 throughout the whole experiment. Compound stocks are dissolved in 10% DMSO and final concentrations made fresh in complete media as needed. Then, 10% DMSO solvent is used as control in complete media. Cells were seeded onto 12-well plates and allowed 12–24 hr until treatment to allow cells to adhere to wells. This is followed by minimum 5 d treatment of compounds, changing the media with treatment every 48 hr. Cell lines A549 and SK-MEL-5 were seeded at 8000 cells/well, which was determined to provide 80–100% confluency after 6 d growth. SH-SY5Y was seeded at 8000 cells/well and 16,000 cells/well due to their known slower division rate and given at least 10 d for sufficient growth. After allotted growth times, cell counts were performed using the Countess automated cell counter. All assays had a minimum of three replicates for all conditions and analyzed through GraphPad Prism one-way ANOVA for significance. Cell lines were from ATCC: HEK-293T, A549, SH-SY5Y, and SK-MEL-5. All four cell lines were evaluated to be negative for Mycoplasma using the following tests: A549: Hoechst DNA stain (indirect) method – negative, Agar culture (direct) method – negative, and PCR-based assay – negative. The HEK-293T, SH-SY5Y, and SK-MEL-5 were determined to be negative based on PCR-based assay.
Odor exposure protocol for Huntington’s disease model flies
Request a detailed protocolFlies were exposed to diacetyl in a cylindrical container (112 mm diameter × 151 mm height). Each container was tightly closed but had two holes, one of which connected to an air suction port, and the other to a vial containing either of 5 ml paraffin oil or 5 ml 1% diacetyl in paraffin oil. A gentle suction was applied to pull the headspace from the odor or paraffin oil vials into the cylindrical structure. pGMR-HTTQ120 flies were maintained at 18°C. Adult flies aged 1 d were transferred to fly vials containing fresh medium and put into the odor-filled container at room temperature. Paraffin oil and 1% diacetyl solution were prepared and replaced every day. At the end of the fifth day of exposure, half of the flies were collected and subjected to pseudopupil analysis as before (Steffan et al., 2001). The remaining flies were transferred to fresh medium and exposed to the odors for an additional 5 d. All treatments and experiments were performed at room temperature. The number of rhabdomeres in each ommatidium was counted using the pseudopupil analysis, and statistical significance was determined by unpaired t-test (two-tailed) against PO at each time point.
Preparation of nuclear extracts from HEK293 cells
Request a detailed protocolNuclear extracts of HEK293 cells were prepared according to a protocol described previously (Andrews and Faller, 1991), with minor modifications. In brief, HEK293 cells were washed twice with cold phosphate-buffered saline (PBS) and lysed with hypotonic buffer (10 mM HEPES-KOH [pH 7.9], 1.5 mM MgCl2, 10 mM KCl, protease inhibitor cocktail [04693159001, Roche, Indianapolis, IN], 1 mM DTT, 1 mM TSA). Following a brief centrifugation, the pellet was resuspended in hypertonic buffer (20 mM HEPES-KOH [pH 7.9], 25% glycerol, 420 mM NaCl 1.5 mM MgCl2, 0.2 mM EDTA, protease inhibitor cocktail, 1 mM DTT, 1 mM TSA). The supernatant was recovered as nuclear extract.
Western blot analysis
Request a detailed protocolProteins in the nuclear extracts (60 μg protein) were separated by SDS-PAGE gels (456-1043, Bio-Rad, Hercules, CA), transferred onto PVDF membranes (162-0174, Bio-Rad), and incubated with anti-histone antibodies: acetylated H3K9 (1/2000: ab4441, abcam, Cambridge, MA), acetylated H3K14 (1/5000: 06-911, EMD Millipore, Billerica, MA), and acetylated H4K5 (1/2000: 07-327, EMD Millipore). Bound antibody was detected by horseradish peroxidase-conjugated anti-rabbit secondary antibody (1/20,000: 1705046, Bio-Rad) and developed using Clarity Western ECL Substrate (1705060, Bio-Rad). Signals were detected and captured using ImageQuant LAS 4000 mini (GE Healthcare, Pittsburgh, PA), and band intensities were quantified with ImageJ software. H3K9 acetylation intensity in individual lanes was reported relative to the normalized Mock treatment and calculated using this formula: relative H3K9ace intensity for each timepoint = (H3K9ace/PCNA)/averaged (Mock H3K9ace/Mock PCNA). Statistical significance was determined by unpaired t-test (two-tailed) against mock at each time point.
Materials availability
Request a detailed protocolAll newly created reagents in this publication such as primers are available upon request.
Data availability
Sequencing data have been deposited in GEO under accession code GSE116502.
-
NCBI Gene Expression OmnibusID GSE116502. A conserved odor detection pathway via modulation of chromatin and cellular gene expression.
References
-
Wiring stability of the adult Drosophila olfactory circuit after lesionThe Journal of Neuroscience 26:3367–3376.https://doi.org/10.1523/JNEUROSCI.4941-05.2006
-
Near-optimal probabilistic RNA-seq quantificationNature Biotechnology 34:525–527.https://doi.org/10.1038/nbt.3519
-
Brain uptake of neurotherapeutics after intranasal versus intraperitoneal delivery in miceJournal of Neurology and Neurosurgery 2:009.
-
Multiple roles of HDAC inhibition in neurodegenerative conditionsTrends in Neurosciences 32:591–601.https://doi.org/10.1016/j.tins.2009.06.002
-
HDAC inhibitors as therapy for neural disorders. Discovery of a new therapyPharmazie in Unserer Zeit 39:204–209.https://doi.org/10.1002/pauz.201000367
-
Histone acetylation: molecular mnemonics on the chromatinNature Reviews. Neuroscience 14:97–111.https://doi.org/10.1038/nrn3427
-
Molecular evolution of the histone deacetylase family: functional implications of phylogenetic analysisJournal of Molecular Biology 338:17–31.https://doi.org/10.1016/j.jmb.2004.02.006
-
BookFlavour determinants of beer qualityIn: Hughes PS, Baxter ED, editors. Beer: Quality, Safety and Nutritional Aspects. The Royal Society od Chemistry. pp. 40–73.https://doi.org/10.1039/9781847550224
-
Microbiological and biochemical origins of human axillary odourFEMS Microbiology Ecology 83:527–540.https://doi.org/10.1111/1574-6941.12054
-
Therapeutic application of histone deacetylase inhibitors for central nervous system disordersNature Reviews. Drug Discovery 7:854–868.https://doi.org/10.1038/nrd2681
-
Mapping odor valence in the brain of flies and miceCurrent Opinion in Neurobiology 24:34–38.https://doi.org/10.1016/j.conb.2013.08.010
-
SAHA, a HDAC inhibitor, has profound anti-growth activity against non-small cell lung cancer cellsOncology Reports 15:187–191.
-
Chemical signaling between gut microbiota and host chromatin: What is your gut really saying?The Journal of Biological Chemistry 292:8582–8593.https://doi.org/10.1074/jbc.R116.761577
-
Influence of valine and other amino acids on total diacetyl and 2,3-pentanedione levels during fermentation of brewer’s wortApplied Microbiology and Biotechnology 97:6919–6930.https://doi.org/10.1007/s00253-013-4955-1
-
Reassessment of the Influence of Malolactic Fermentation on the Concentration of Diacetyl in WinesAmerican Journal of Enology and Viticulture 46:385–388.https://doi.org/10.5344/ajev.1995.46.3.385
-
The butter flavorant, diacetyl, forms a covalent adduct with 2-deoxyguanosine, uncoils DNA, and leads to cell deathJournal of Agricultural and Food Chemistry 60:3311–3317.https://doi.org/10.1021/jf300180e
-
Respiratory toxicity of diacetyl in C57BL/6 miceToxicological Sciences 103:169–180.https://doi.org/10.1093/toxsci/kfn016
-
BookDistilled BeveragesIn: Nykanen L, Nykanen I, editors. Volatile Compounds in Foods and Beverages. Taylor & Francis Group. pp. 548–580.https://doi.org/10.1201/9780203734285
-
The HDAC inhibitor trichostatin A inhibits growth of small cell lung cancer cellsThe Journal of Surgical Research 142:219–226.https://doi.org/10.1016/j.jss.2006.12.555
-
The genetic landscape of high-risk neuroblastomaNature Genetics 45:279–284.https://doi.org/10.1038/ng.2529
-
Emission of diacetyl (2,3 butanedione) from natural butter, microwave popcorn butter flavor powder, paste, and liquid productsInternational Journal of Occupational and Environmental Health 16:291–302.https://doi.org/10.1179/107735210799160237
-
Normalization of RNA-seq data using factor analysis of control genes or samplesNature Biotechnology 32:896–902.https://doi.org/10.1038/nbt.2931
-
Gut microbiota, metabolites and host immunityNature Reviews. Immunology 16:341–352.https://doi.org/10.1038/nri.2016.42
-
Functions of site-specific histone acetylation and deacetylationAnnual Review of Biochemistry 76:75–100.https://doi.org/10.1146/annurev.biochem.76.052705.162114
-
Histone deacetylases as regulators of inflammation and immunityTrends in Immunology 32:335–343.https://doi.org/10.1016/j.it.2011.04.001
-
Histone deacetylase inhibitors: mechanisms and clinical significance in cancer: HDAC inhibitor-induced apoptosisAdvances in Experimental Medicine and Biology 615:261–298.https://doi.org/10.1007/978-1-4020-6554-5_13
-
Diacetyl: occurrence, analysis, and toxicityJournal of Agricultural and Food Chemistry 62:4048–4053.https://doi.org/10.1021/jf500615u
-
Activity-dependent changes to the brain and behavior of the honey bee, Apis mellifera (L.)The Journal of Neuroscience 17:7148–7156.https://doi.org/10.1523/JNEUROSCI.17-18-07148.1997
-
Diacetyl exposure as a pneumotoxic factor: a reviewRoczniki Panstwowego Zakladu Higieny 65:87–92.
-
Sirtinol, a class III HDAC inhibitor, induces apoptotic and autophagic cell death in MCF-7 human breast cancer cellsInternational Journal of Oncology 41:1101–1109.https://doi.org/10.3892/ijo.2012.1534
-
HDAC inhibitor M344 suppresses MCF-7 breast cancer cell proliferationBiomedicine & Pharmacotherapy = Biomedecine & Pharmacotherapie 66:232–236.https://doi.org/10.1016/j.biopha.2011.06.007
-
Identification of potential prognostic genes for neuroblastomaFrontiers in Genetics 9:589.https://doi.org/10.3389/fgene.2018.00589
Article and author information
Author details
Funding
National Institute of Allergy and Infectious Diseases (R01AI087785)
- Sachiko Haga-Yamanaka
- Rogelio Nunez-Flores
- Christi A Scott
- Sarah Perry
- Stephanie Turner Chen
- Crystal Pontrello
- Anandasankar Ray
National Institute of Allergy and Infectious Diseases (R01AI153195)
- Meera G Nair
University of California, Riverside
- Anandasankar Ray
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Drs. Maurizio Pellecchia, Parima Udompholkul, Francis Sladek, and Xuan Liu for cell lines. We would like to thank Dr. Larry Marsh for advice about the fly eye assay and HDAC inhibitors, Dr. Joel Kowalewski for help with critical comments, Dr. Thomas Girke for guidance in bioinformatics analysis, Drs. Anupama Dahanukar, Naoki Yamanaka, Chih-cheng Tsai, Karine LeRoch. and Adler R Dillman for sharing equipment and reagents. This work was partially supported by funds provided by the National Institutes of Health, NIAID (R01AI153195 to MGN and R01AI087785 to SHY, RNF, CAS, SP, STC, CP and AR), and HATCH funds from the Agricultural Experimentation Station UC Riverside to AR.
Ethics
All protocols for animal use and euthanasia were approved by the Institutional Animal Care and Use Committee (#20150028) and were in accordance with the provisions established by the Animal Welfare Act and the Public Health Services (PHS) Policy.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.86823. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Haga-Yamanaka, Nunez-Flores, Scott et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,459
- views
-
- 215
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling due to steric similarities to endogenous hormones. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) and differentially expressed genes (DEGs) that predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many EDCs persists over multiple generations. Many studies have described direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDC-induced epimutations are introduced or subsequently propagated, whether there are cell type-specific susceptibilities to the same EDC, and whether this correlates with differential expression of relevant hormone receptors. We exposed cultured pluripotent (iPS), somatic (Sertoli and granulosa), and primordial germ cell-like (PGCLC) cells to BPS and found that differential incidences of BPS-induced epimutations and DEGs correlated with differential expression of relevant hormone receptors inducing epimutations near relevant hormone response elements in somatic and pluripotent, but not germ cell types. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and DEGs was largely retained, however, >90% of the specific epimutations and DEGs were replaced by novel epimutations and DEGs. These results suggest a unique mechanism by which an EDC-induced epimutated state may be propagated transgenerationally.
-
- Genetics and Genomics
Microbial secondary metabolites are a rich source for pharmaceutical discoveries and play crucial ecological functions. While tools exist to identify secondary metabolite clusters in genomes, precise sequence-to-function mapping remains challenging because neither function nor substrate specificity of biosynthesis enzymes can accurately be predicted. Here, we developed a knowledge-guided bioinformatic pipeline to solve these issues. We analyzed 1928 genomes of Pseudomonas bacteria and focused on iron-scavenging pyoverdines as model metabolites. Our pipeline predicted 188 chemically different pyoverdines with nearly 100% structural accuracy and the presence of 94 distinct receptor groups required for the uptake of iron-loaded pyoverdines. Our pipeline unveils an enormous yet overlooked diversity of siderophores (151 new structures) and receptors (91 new groups). Our approach, combining feature sequence with phylogenetic approaches, is extendable to other metabolites and microbial genera, and thus emerges as powerful tool to reconstruct bacterial secondary metabolism pathways based on sequence data.