Search NASA⌕ Search

SEARCH · Search NASA

Results for “methylation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Low rotational barriers for the most dynamically active methyl groups in the proposed antiviral drugs for treatment of SARS-CoV-2, apilimod and tetrandrine

A recent screening study highlighted a molecular compound, apilimod, for its efficacy against the SARS-CoV-2 virus, while another compound, tetrandrine, demonstrated a remarkable synergy with the benchmark antiviral drug, remdesivir. In this work, we find that because of significantly reduced potential energy barriers, which also give rise to pronounced quantum effects, the rotational dynamics of the most dynamically active methyl groups in apilimod and tetrandrine are much faster than those in remdesivir. Because dynamics of methyl groups are essential for biochemical activity, screening studies based on the computed potential energy profiles may help identify promising candidates within a given class of drugs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

High-pressure ignition delay time measurements of a four-component gasoline surrogate and its high-level blends with ethanol and methyl acetate

Ignition delay times of a gasoline surrogate (iso-octane/toluene/n-heptane/1-hexene at 55/25/15/5% by liquid volume, developed in McCormick et al. 2017) with high-level bioblendstock-gasoline surrogate blends (50% and 85% biofuel by liquid volume of each ethanol and methyl acetate) were collected behind reflected shock waves. Post-reflected-shock temperatures ranged from 968 to 1361 K at pressures of about 4, 10, 25, and 50 atm. Data were collected for real fuel–air mixtures at fuel lean and stoichiometric conditions (φ = 0.5, 1.0) with focus on the more-practical, fuel-lean conditions. Ignition delay times were measured from OH* chemiluminescence around 307 nm using an endwall diagnostic. The RON and MON of the surrogate were 90.3 and 84.7, respectively. Ethanol was chosen due to its wide use in flex-fuel vehicles and availability, while also providing an increased octane rating. Methyl acetate was chosen for its especially high octane rating to investigate the effect of an extreme case. To validate the gasoline surrogate, the data are compared to real gasoline (RD387) and several gasoline surrogate experiments from the literature. Using the wide range of pressures studied, the gasoline surrogate’s pressure dependence was quantified to account for test-to-test variations using regression analysis. Similarly, a global correlation for gasoline and its surrogates was developed using all available data from the literature. Two modern chemical kinetics models targeting gasoline and its surrogates are compared to the ignition delay time measurements. These new high-pressure, high-bioblendstock concentration tests provide required chemical kinetic data for optimizing fuel and engine design.

42 ENGINEERING↗
DOE OSTI2025

Ecology of methyl-coenzyme M reductase encoding Thermoproteota

The recent demonstration that members of at least three classes of archaea affiliated with the Thermoproteota superphylum are involved in the production of the climate-active gas methane has sparked discussions about how well we understand the diversity of methanogens. Here, in this work, we show that members of all three of these lineages, as well as several other, yet uncultured and physiologically uncharacterized groups within the Thermoproteota that encode the key enzyme of anaerobic methane cycling, methyl-coenzyme M reductase (MCR), are widely distributed in anoxic ecosystems. We postulate that the taxonomic, metabolic, and ecological diversity of methanogenic and MCR-encoding Thermoproteota are poorly understood, and that the contribution of methylotrophic and thermoproteotal methanogenesis to methane production is largely unknown. We hypothesize that thermoproteotal methanogens could contribute, potentially substantially, to methane emissions in many anoxic environments that harbor methylated precursors, including wetlands, sediments, peat, rice paddies, wastewater sludge, and geothermal systems. We highlight the necessity to experimentally test the (eco)physiology of these widely distributed archaea using both culture-dependent ( in vitro ) and culture-independent ( in situ ) approaches to assess their potential contribution to methane emissions. Last, we stress the importance of remaining agnostic about the physiology of MCR-encoding Thermoproteota in the absence of experimental data because most of these archaea also carry the genetic potential to grow non-methanogenically.

MCR↗
DOE OSTI2021

Separation, speciation, and mechanism of astatine and bismuth extraction from nitric acid into 1-octanol and methyl anthranilate

We report a detailed study of At and Bi extraction from nitric acid media into conventional solvents, namely 1-octanol and methyl anthranilate, has been performed. The analysis includes a mathematical modeling which allows the fitting of experimental data and determination of extraction constants of the two above mentioned elements. Also, this approach helped to estimate a stability constant of a weak AtO(NO3) complex along with thermodynamic constants describing the redox process of At species in the acidic solution and formation of an adduct of Bi in the presence of methyl anthranilate. The results of the fitting have been used to calculate corresponding separation factors of Bi and At as well. Moreover, a computational study has been performed to evaluate At interaction with the above mentioned solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Conformational Changes Induced by Methyl Side-Chains in Protonated Tripeptides Containing Glycine and Alanine Residues

We present a systematic study of the conformational and isomeric populations in gas-phase protonated tripeptides containing glycine and alanine residues using infrared predissociation spectroscopy of cryogenically cooled ions. Specifically, the protonated forms of Gly-Gly-Gly, Ala-Gly-Gly, Gly-Ala-Gly, Gly-Gly-Ala, Ala-Ala-Gly, Ala-Gly-Ala, Gly-Ala-Ala, and Ala-Ala-Ala allow us to sample all permutations of the methyl side-chain position, providing a comprehensive view of the effects of this simple side-chain on the 3-D structure of the peptide. The individual structural populations for all but one of these peptide species are determined via conformer-specific IR–IR double-resonance spectroscopy and comparison with electronic structure predictions. The observed structures can be classified into three main families defined by the protonation site and the number of internal hydrogen bonds. The relative contribution of each structural family is highly dependent on the exact amino acid sequence of the tripeptide. Here, these observed changes in structural population can be rationalized in terms of the electron-donating effect of the methyl side-chain modulating the local proton affinities of the amine and various carbonyl groups in the tripeptide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Bimolecular Reaction of Methyl-Ethyl-Substituted Criegee Intermediate with SO 2

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO 2 , a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10 –10 cm 3 s –1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO 3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO 2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO 2 (k syn = 4.42 × 10 11 T –7.80 exp(-1401/T) cm 3 s –1 and k anti = 1.26 × 10 11 T –7.55 exp(-1397/T) cm 3 s –1 ). In conclusion, our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO 2 is significantly faster than MVK-oxide + SO 2 , demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modeling–Experiment–Theory Analysis of Reactions Initiated from Cl + Methyl Formate

Methyl formate (MF; CH 3 OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH 3 OCO and CH 2 OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl 2 was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl 2 , HCl, CO 2 , CH 3 , CH 3 Cl, CH 2 O, and CH 2 ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH 3 OCO (R1a) and Cl + MF → HCl + CH 2 OCHO (R1b): k 1a,theory = 6.71 × 10 –15 ·T 1.14 ·exp(—606/T) cm 3 /molecule·s; k 1b,theory = 4.67 × 10 –18 ·T 2.21 ·exp(—245/T) cm 3 /molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH 3 OCO and CH 2 OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH 3 + CO 2 (R3) and CH 2 O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k 3,theory (P = ∞) = 2.94 × 10 9 ·T 1.21 ·exp(—6209/T) s –1 and k 5,theory (P = ∞) = 8.45 × 10 8 ·T 1.39 ·exp(—15132/T) s –1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH 2 OCHO + Cl, CH 2 OCHO + Cl 2 , and CH 3 + Cl 2 were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dual Feedstock Upcycling of α-Methylstyrene-Doped Poly(methyl methacrylate) and Biomass via the Telescope of Depolymerization and Diels–Alder Reaction

Nearly 90% of poly(methyl methacrylate) (PMMA) is not recycled and instead ends up in landfills. Conventional pyrolysis of PMMA recovers impure methyl methacrylate (MMA) with low economic value. Here, we present a telescoped dual upcycling strategy that integrates PMMA depolymerization, Diels–Alder cycloaddition, and aromatization to convert AMS-doped PMMA and biomass-derived 2,5-dimethylfuran (DMF) into 1,2,4-trimethylbenzene (pseudocumene), a valuable chemical feedstock. BBr 3 proved effective in promoting the challenging Diels–Alder reaction between MMA and DMF under high pressure of argon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Restricted Rotational Flexibility of the C5α-Methyl-Substituted Carbapenem NA-1-157 Leads to Potent Inhibition of the GES-5 Carbapenemase

Carbapenem antibiotics are used as a last-resort treatment for infections caused by multidrug-resistant bacteria. The wide spread of carbapenemases in Gram-negative bacteria has severely compromised the utility of these drugs and represents a serious public health threat. To combat carbapenemase-mediated resistance, new antimicrobials and inhibitors of these enzymes are urgently needed. Here, we describe the interaction of the atypically C5α-methyl-substituted carbapenem, NA-1-157, with the GES-5 carbapenemase. MICs of this compound against Escherichia coli, Klebsiella pneumoniae, and Acinetobacter baumannii producing the enzyme were reduced 4–16-fold when compared to MICs of the commercial carbapenems, reaching clinically sensitive breakpoints. When NA-1-157 was combined with meropenem, a strong synergistic effect was observed. Kinetic and ESI-LC/MS studies demonstrated that NA-1-157 is a potent inhibitor of GES-5, with a high inactivation efficiency of (2.9 ± 0.9) × 105 M –1 s –1 . Acylation of GES-5 by NA-1-157 was biphasic, with the fast phase completing within seconds, and the slow phase taking several hours and likely proceeding through a reversible tetrahedral intermediate. Deacylation was extremely slow (k 3 = (2.4 ± 0.3) × 10 –7 s –1 ), resulting in a residence time of 48 ± 6 days. MD simulation of the GES-5-meropenem and GES-5-NA-1-157 acyl-enzyme complexes revealed that the C5α-methyl group in NA-1- 157 sterically restricts rotation of the 6α-hydroxyethyl group preventing ingress of the deacylating water into the vicinity of the scissile bond of the acyl-enzyme intermediate. Furthermore, these data demonstrate that NA-1-157 is a potent irreversible inhibitor of the GES-5 carbapenemase.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Electrochemically Initiated Depolymerization of Poly(Methyl Methacrylate)

Efficient depolymerization of polymers with all-carbon backbones under mild conditions would be valuable in chemically recycling commodity plastics. Poly(methyl methacrylate) (PMMA) is a commodity thermoplastic that is currently depolymerized under temperatures in excess of 400 °C. Herein, we lower the temperatures needed to achieve depolymerization of PMMA by performing radical generation and depropagation with orthogonal stimuli. This first demonstration of electrochemically initiated PMMA depolymerization relies on reduction of phthalimide esters that, upon subsequent decarboxylation, generate polymer-centered radicals. These radicals then spontaneously unzip the polymer back to its monomeric constituents at temperatures as low as 105 °C. We studied the mechanism and efficiency of this transformation as a function of phthalimide ester placement, incorporation density, and polymer molecular weight. We found that chain-end activation is effective for modest molecular weights but suffers diminished efficiency at higher degrees of polymerization. In contrast, pendent-group activation is more effective for depolymerizing higher molecular weight species. Integrating higher molar amounts of phthalimide ester pendants leads to more effective depolymerization, with >95% depolymerization in copolymers with 5 mol% phthalimide ester incorporation. We leveraged this understanding to create a custom electro-distillation apparatus that allowed us to simultaneously electrochemically depolymerize PMMA and directly distill methyl methacrylate in >22% yield, which could ultimately be repolymerized. These findings establish electrochemistry as a versatile and orthogonal stimulus for vinyl polymer depolymerization and provide a foundation for closed-loop electrochemical recycling of widely used plastics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

rRNA methylation by Spb1 regulates the GTPase activity of Nog2 during 60S ribosomal subunit assembly

Biogenesis of the large ribosomal (60S) subunit involves the assembly of three rRNAs and 46 proteins, a process requiring approximately 70 ribosome biogenesis factors (RBFs) that bind and release the pre-60S at specific steps along the assembly pathway. The methyltransferase Spb1 and the K-loop GTPase Nog2 are essential RBFs that engage the rRNA A-loop during sequential steps in 60S maturation. Spb1 methylates the A-loop nucleotide G2922 and a catalytically deficient mutant strain ( spb 1 D52A ) has a severe 60S biogenesis defect. However, the assembly function of this modification is currently unknown. Here, we present cryo-EM reconstructions that reveal that unmethylated G2922 leads to the premature activation of Nog2 GTPase activity and capture a Nog2-GDP-AlF 4 - transition state structure that implicates the direct involvement of unmodified G2922 in Nog2 GTPase activation. Genetic suppressors and in vivo imaging indicate that premature GTP hydrolysis prevents the efficient binding of Nog2 to early nucleoplasmic 60S intermediates. We propose that G2922 methylation levels regulate Nog2 recruitment to the pre-60S near the nucleolar/nucleoplasmic phase boundary, forming a kinetic checkpoint to regulate 60S production. Our approach and findings provide a template to study the GTPase cycles and regulatory factor interactions of the other K-loop GTPases involved in ribosome assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Methylation enables the use of fluorine-free ether electrolytes in high-voltage lithium metal batteries

Lithium metal batteries represent a promising technology for next-generation energy storage, but they still suffer from poor cycle life due to lithium dendrite formation and cathode cracking. Fluorinated solvents can improve battery longevity by improving LiF content in the solid-electrolyte interphase; however, the high cost and environmental concerns of fluorinated solvents limit battery viability. Here, for this work, we designed a series of fluorine-free solvents through the methylation of 1,2-dimethoxyethane, which promotes inorganic LiF-rich interphase formation through anion reduction and achieves high oxidation stability. The anion-derived LiF interphases suppress lithium dendrite growth on the lithium anode and minimize cathode cracking under high-voltage operation. The Li + -solvent structure is investigated through in situ techniques and simulations to draw correlations between the interphase compositions and electrochemical performances. The methylation strategy provides an alternative pathway for electrolyte engineering towards high-voltage electrolytes while reducing dependence on expensive fluorinated solvents.

25 ENERGY STORAGE↗
DOE OSTI2020

Structure determination of the HgcAB complex using metagenome sequence data: insights into microbial mercury methylation

Bacteria and archaea possessing the hgcAB gene pair methylate inorganic mercury (Hg) to form highly toxic methylmercury. HgcA consists of a corrinoid binding domain and a transmembrane domain, and HgcB is a dicluster ferredoxin. However, their detailed structure and function have not been thoroughly characterized. We modeled the HgcAB complex by combining metagenome sequence data mining, coevolution analysis, and Rosetta structure calculations. In addition, we overexpressed HgcA and HgcB in Escherichia coli, confirmed spectroscopically that they bind cobalamin and [4Fe-4S] clusters, respectively, and incorporated these cofactors into the structural model. Surprisingly, the two domains of HgcA do not interact with each other, but HgcB forms extensive contacts with both domains. The model suggests that conserved cysteines in HgcB are involved in shuttling HgII, methylmercury, or both. These findings refine our understanding of the mechanism of Hg methylation and expand the known repertoire of corrinoid methyltransferases in nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Experimental and theoretical studies of the doubly substituted methyl-ethyl Criegee intermediate: Infrared action spectroscopy and unimolecular decay to OH radical products

In this work the infrared (IR) action spectrum of the doubly substituted methyl-ethyl Criegee intermediate (MECI) is observed in the CH stretch overtone region with detection of OH products. The MECI exhibits four conformers, all of which undergo unimolecular decay via a 1,4 H-atom transfer mechanism, followed by the rapid release of OH products. Conformers with different orientations of the carbonyl oxide group with respect to the methyl and ethyl substituents (i.e., anti and syn ) decay via distinct transition state barriers (16.1 kcal mol –1 and 15.4 kcal mol –1 , respectively). The observed IR action spectrum is in good agreement with the predicted anharmonic IR absorption spectrum, but exhibits significant congestion, which is attributed to couplings between spectroscopic bright states and nearby dark states. Energy-dependent OH appearance rates are measured upon IR excitation of the strongest features in the IR action spectrum and are found to be on the order of 10 6 –10 7 s –1 . The experimental rates are in good agreement with computed Rice–Ramsperger–Kassel–Marcus rates for the unimolecular decay of MECI at these energies, which incorporate quantum mechanical tunneling and sophisticated hindered rotor treatments, as well as high-level theoretical calculations of the TS barrier heights, rovibrational properties, and torsional barriers associated with the MECI conformers. Master equation modeling is used to predict thermal rates for the unimolecular decay of anti - and syn -MECI of 473 s –1 and 660 s –1 , respectively. Comparison with other previously studied Criegee intermediate systems provides insights into substituent effects on unimolecular decay under both energy-dependent and thermal conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Photodissociation dynamics of methyl vinyl ketone oxide: A four-carbon unsaturated Criegee intermediate from isoprene ozonolysis

The electronic spectrum of methyl vinyl ketone oxide (MVK-oxide), a four-carbon Criegee intermediate derived from isoprene ozonolysis, is examined on its second π* ← π transition, involving primarily the vinyl group, at UV wavelengths (λ) below 300 nm. A broad and unstructured spectrum is obtained by a UV-induced ground state depletion method with photoionization detection on the parent mass (m/z 86). Electronic excitation of MVK-oxide results in dissociation to O ( 1 D) products that are characterized using velocity map imaging. Electronic excitation of MVK-oxide on the first π* ← π transition associated primarily with the carbonyl oxide group at λ > 300 nm results in a prompt dissociation, and yields broad total kinetic energy release (TKER) and anisotropic angular distributions for the O ( 1 D) + methyl vinyl ketone products. By contrast, electronic excitation at λ ≤ 300 nm results in bimodal TKER and angular distributions, indicating two distinct dissociation pathways to O ( 1 D) products. Furthermore, one pathway is analogous to that at λ > 300 nm, while the second pathway results in very low TKER and isotropic angular distributions indicative of internal conversion to the ground electronic state and statistical unimolecular dissociation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Knock-limited combustion of ethanol-, isobutanol-, and 2-methyl-3-buten-2-ol-gasoline blends in a direct-injected spark-ignition engine

Knock-limited combustion of alcohol-gasoline fuel blends was studied in a direct-injection spark-ignition engine. Ethanol, isobutanol (2-methyl-1-propanol), and methylbutenol (2-methyl-3-buten-2-ol) were splash-blended with a blendstock for oxygenate blending (BOB) gasoline. Ethanol was blended in fractions of 10%, 20%, and 30% by volume (E10, E20, and E30). Isobutanol and methylbutenol were blended to match the oxygen weight percentage of the ethanol blends, resulting in blends of 16%, 32%, and 49% (I16, I32, and I49) isobutanol by volume and 18%, 37%, and 56% (M18, M37, and M56) methylbutenol by volume. Neat BOB and gasoline primary reference fuels (PRF) with octane numbers (ON) of 87 (PRF87) and 100 (PRF100) were included for reference. The engine was operated at fixed speed, equivalence ratio, and injection timing while knock-limited spark advance was located across a range of intake pressures. Low-level alcohol blends E10 and M18 had similar knocking behavior and appeared to be slightly more knock-resistant than I16. For mid-level blends, E20 showed further improvement over I32 while M37 showed noticeably better knock resistance, able to match the knock-resistance of PRF100 for loads above 1000 kPa gross indicated mean effective pressure (GIMEP). For high-level blends, the difference between E30 and I49 was similar to their mid-level counterparts. However, M56 showed significant improvement over E30 and I49 throughout the entire range of loads tested. M56 exceeded the knock-resistance of PRF100 for loads above 700 kPa GIMEP. The improved knock-resistance combined with the increased volumetric energy density of methylbutenol resulted in improved fuel consumption and indicated efficiency. The overall results of this study indicate that the order of knock-resistance of these three alcohols, for direct-injection operation with injection during the intake stroke, is isobutanol < ethanol < methylbutenol when compared on an equal oxygen weight basis.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Genome-wide profiling of histone (H3) lysine 4 (K4) tri-methylation (me3) under drought, heat, and combined stresses in switchgrass

Background: Switchgrass (Panicum virgatum L.) is a warm-season perennial (C4) grass identified as an important biofuel crop in the United States. It is well adapted to the marginal environment where heat and moisture stresses predominantly affect crop growth. However, the underlying molecular mechanisms associated with heat and drought stress tolerance still need to be fully understood in switchgrass. The methylation of H3K4 is often associated with transcriptional activation of genes, including stress-responsive. Therefore, this study aimed to analyze genome-wide histone H3K4-tri-methylation in switchgrass under heat, drought, and combined stress. Results: In total, ~ 1.3 million H3K4me3 peaks were identified in this study using SICER. Among them, 7,342; 6,510; and 8,536 peaks responded under drought (DT), drought and heat (DTHT), and heat (HT) stresses, respectively. Most DT and DTHT peaks spanned 0 to + 2000 bases from the transcription start site [TSS]. By comparing differentially marked peaks with RNA-Seq data, we identified peaks associated with genes: 155 DT-responsive peaks with 118 DT-responsive genes, 121 DTHT-responsive peaks with 110 DTHT-responsive genes, and 175 HT-responsive peaks with 136 HT-responsive genes. We have identified various transcription factors involved in DT, DTHT, and HT stresses. Gene Ontology analysis using the AgriGO revealed that most genes belonged to biological processes. Most annotated peaks belonged to metabolite interconversion, RNA metabolism, transporter, protein modifying, defense/immunity, membrane traffic protein, transmembrane signal receptor, and transcriptional regulator protein families. Further, we identified significant peaks associated with TFs, hormones, signaling, fatty acid and carbohydrate metabolism, and secondary metabolites. qRT-PCR analysis revealed the relative expressions of six abiotic stress-responsive genes (transketolase, chromatin remodeling factor-CDH3, fatty-acid desaturase A, transmembrane protein 14C, beta-amylase 1, and integrase-type DNA binding protein genes) that were significantly (P < 0.05) marked during drought, heat, and combined stresses by comparing stress-induced against un-stressed and input controls. Conclusion: Our study provides a comprehensive and reproducible epigenomic analysis of drought, heat, and combined stress responses in switchgrass. Significant enrichment of H3K4me3 peaks downstream of the TSS of protein-coding genes was observed. In addition, the cost-effective experimental design, modified ChIP-Seq approach, and analyses presented here can serve as a prototype for other non-model plant species for conducting stress studies.

59 BASIC BIOLOGICAL SCIENCES↗