Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gene expression”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiple highly expressed phospho enol pyruvate carboxylase genes have divergent enzyme kinetic properties in two C4 grasses

Abstract Background and Aims Phosphoenolpyruvate (PEP) carboxylase (PEPC) catalyses the irreversible carboxylation of PEP with bicarbonate to produce oxaloacetate. This reaction powers the carbon-concentrating mechanism (CCM) in plants that perform C4 photosynthesis. This CCM is generally driven by a single PEPC gene product that is highly expressed in the cytosol of mesophyll cells. We found two C4 grasses, Panicum miliaceum and Echinochloa colona, that each have two highly expressed PEPC genes. We characterized the kinetic properties of the two most abundant PEPCs in E. colona and P. miliaceum to better understand how the enzyme’s amino acid structure influences its function. Methods Coding sequences of the two most abundant PEPC proteins in E. colona and P. miliaceum were synthesized by GenScript and were inserted into bacteria expression plasmids. Point mutations resulting in substitutions at conserved amino acid residues (e.g. N-terminal serine and residue 890) were created via site-directed PCR mutagenesis. The kinetic properties of semi-purified plant PEPCs from Escherichia coli were analysed using membrane-inlet mass spectrometry and a spectrophotometric enzyme-coupled reaction. Key Results The two most abundant P. miliaceum PEPCs (PmPPC1 and PmPPC2) have similar sequence identities (>95 %), and as a result had similar kinetic properties. The two most abundant E. colona PEPCs (EcPPC1 and EcPPC2) had identities of ~78 % and had significantly different kinetic properties. The PmPPCs and EcPPCs had different responses to allosteric inhibitors and activators, and substitutions at the conserved N-terminal serine and residue 890 resulted in significantly altered responses to allosteric regulators. Conclusions The two, significantly expressed C4Ppc genes in P. miliaceum were probably the result of genomes combining from two closely related C4Panicum species. We found natural variation in PEPC’s sensitivity to allosteric inhibition that seems to bypass the conserved 890 residue, suggesting alternative evolutionary pathways for increased malate tolerance and other kinetic properties.

DiMario, Robert J. (ORCID:0000000250566868)↗
DOE OSTI2024

Enhancers in Plant Development, Adaptation and Evolution

Understanding plant responses to developmental and environmental cues is crucial for studying morphological divergence and local adaptation. Gene expression changes, governed by cis-regulatory modules (CRMs) including enhancers, are a major source of plant phenotypic variation. However, while genome-wide approaches have revealed thousands of putative enhancers in mammals, far fewer have been identified and functionally characterized in plants. This review provides an overview of how enhancers function to control gene regulation, methods to predict DNA sequences that may have enhancer activity, methods utilized to functionally validate enhancers and the current knowledge of enhancers in plants, including how they impact plant development, response to environment and evolutionary adaptation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Spray‐induced gene silencing to identify powdery mildew gene targets and processes for powdery mildew control

Abstract Spray‐induced gene silencing (SIGS) is an emerging tool for crop pest protection. It utilizes exogenously applied double‐stranded RNA to specifically reduce pest target gene expression using endogenous RNA interference machinery. In this study, SIGS methods were developed and optimized for powdery mildew fungi, which are widespread obligate biotrophic fungi that infect agricultural crops, using the known azole‐fungicide target cytochrome P450 51 (CYP51) in the Golovinomyces orontii–Arabidopsis thaliana pathosystem. Additional screening resulted in the identification of conserved gene targets and processes important to powdery mildew proliferation: apoptosis‐antagonizing transcription factor in essential cellular metabolism and stress response; lipid catabolism genes lipase a , lipase 1 , and acetyl‐CoA oxidase in energy production ; and genes involved in manipulation of the plant host via abscisic acid metabolism ( 9‐cis‐epoxycarotenoid dioxygenase , xanthoxin dehydrogenase , and a putative abscisic acid G‐protein coupled receptor ) and secretion of the effector protein, effector candidate 2 . Powdery mildew is the dominant disease impacting grapes and extensive powdery mildew resistance to applied fungicides has been reported. We therefore developed SIGS for the Erysiphe necator–Vitis vinifera system and tested six successful targets identified using the G. orontii–A. thaliana system. For all targets tested, a similar reduction in powdery mildew disease was observed between systems. This indicates screening of broadly conserved targets in the G. orontii–A. thaliana pathosystem identifies targets and processes for the successful control of other powdery mildew fungi. The efficacy of SIGS on powdery mildew fungi makes SIGS an exciting prospect for commercial powdery mildew control.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The DYRKP1 kinase regulates cell wall degradation in Chlamydomonas by inducing matrix metalloproteinase expression

Abstract The cell wall of plants and algae is an important cell structure that protects cells from changes in the external physical and chemical environment. This extracellular matrix, composed of polysaccharides and glycoproteins, must be constantly remodeled throughout the life cycle. However, compared to matrix polysaccharides, little is known about the mechanisms regulating the formation and degradation of matrix glycoproteins. We report here that a plant kinase belonging to the dual-specificity tyrosine phosphorylation-regulated kinase (DYRKP1) family present in all eukaryotes regulates cell wall degradation after mitosis of Chlamydomonas reinhardtii by inducing the expression of matrix metalloproteinases. Without DYRKP1, daughter cells cannot disassemble parental cell walls and remain trapped inside for more than 10 days. On the other hand, the dual-specificity tyrosine phosphorylation-regulated kinase complementation lines show normal degradation of the parental cell wall. Transcriptomic and proteomic analyses indicate a marked downregulation of MMP gene expression and accumulation, respectively, in the dyrkp1 mutants. The mutants deficient in matrix metalloproteinases retain palmelloid structures for a longer time than the background strain, like dyrkp1 mutants. Our findings show that dual-specificity tyrosine phosphorylation-regulated kinase, by ensuring timely MMP expression, enables the successful execution of the cell cycle. Altogether, this study provides insight into the life cycle regulation in plants and algae.

Kim, Minjae (ORCID:0000000223561295)↗
DOE OSTI2024

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders

Abstract Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN)cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genesMAFB, PHOX2A, CHN1, andEBF3– as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

Science & Technology - Other Topics↗
DOE OSTI2024

Greater aperture counteracts effects of reduced stomatal density on WUE: a case study on sugarcane and meta-analysis

Stomata regulate CO 2 and water vapor exchange between leaves and the atmosphere. Stomata are a target for engineering to improve crop intrinsic water use efficiency (iWUE). One example is by expressing genes that lower stomatal density (SD) and reduce stomatal conductance (g sw ). However, the quantitative relationship between reduced SD, g sw , and the mechanisms underlying it is poorly understood. We addressed this knowledge gap using low-SD sugarcane (Saccharum spp. hybrid) as a case study alongside a meta-analysis of data from 10 species. Transgenic expression of EPIDERMAL PATTERNING FACTOR 2 from Sorghum bicolor (SbEFP2) in sugarcane reduced SD by 26-38% but did not affect gsw compared to wildtype. Further, no changes occurred in stomatal complex size or proxies for photosynthetic capacity. Measurements of gas exchange at low CO 2 concentrations that promote complete stomatal opening to normalize aperture size between genotypes were combined with modeling of maximum gsw from anatomical data. These data suggest that increased stomatal aperture is the only possible explanation for maintaining gsw when SD is reduced. Meta-analysis across C 3 dicots, C 3 monocots, and C 4 monocots revealed engineered reductions in SD are strongly correlated with lower gsw (r 2 =0.60-0.98), but this response is damped relative to the change in anatomy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Bioinformatic and Mechanistic Analysis of the Palmerolide PKS-NRPS Biosynthetic Pathway From the Microbiome of an Antarctic Ascidian

Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the ~75 kbp trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicihabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression.

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

Magnetic resonance control of spin-correlated radical pair dynamics in vivo

Magnetic fields can influence reactions involving spin-correlated radical pairs (SCRPs). This provides a mechanism by which both static and time-varying magnetic fields can affect living systems at the biomolecular level. However, an engineered SCRP system conferring magnetic sensitivity to a non-native biochemical process in a multicellular organism has not yet been demonstrated. Here, in this study, we demonstrate control of SCRP dynamics using magnetic resonance in a live transgenic animal. We show that the emission of various red fluorescent proteins (RFPs), in the presence of a flavin cofactor, can be modified by a combination of static and radiofrequency magnetic fields applied near the electron spin resonance frequency. This effect was measured at room temperature both in vitro and in the nematode Caenorhabditis elegans, genetically modified to express the RFP mScarlet. These observations suggest that the magnetic field effects measured in RFP-flavin systems are due to quantum-correlated radical pairs with a coherence time larger than 4 ns. Our experiments demonstrate that radiofrequency magnetic fields can influence dynamics of reactions involving SCRPs in vivo, potentially enabling new methods for remotely controlling biomolecular processes, such as gene expression, and suggest broader potential for quantum tools in biology.

Burd, Shaun C. [Stanford Univ., CA (United States)↗
DOE OSTI2021

Suppression of hepatitis B virus through therapeutic activation of RIG-I and IRF3 signaling in hepatocytes

Hepatitis B virus (HBV) mediates persistent infection, chronic hepatitis, and liver disease. HBV covalently closed circular (ccc)DNA is central to viral persistence such that its elimination is considered the cornerstone for HBV cure. Inefficient detection by pathogen recognition receptors (PRRs) in the infected hepatocyte facilitates HBV persistence via avoidance of innate immune activation and interferon regulatory factor (IRF)3 induction of antiviral gene expression. We evaluated a small molecule compound, F7, and 5'-triphosphate-poly-U/UC pathogen-associated-molecular-pattern (PAMP) RNA agonists of RIG-I, a PRR that signals innate immunity, for ability to suppress cccDNA. F7 and poly-U/UC PAMP treatment of HBV-infected cells induced RIG-I signaling of IRF3 activation to induce antiviral genes for suppression of cccDNA formation and accelerated decay of established cccDNA, and were additive to the actions of entecavir. Our study shows that activation of the RIG-I pathway and IRF3 to induce innate immune actions offers therapeutic benefit toward elimination of cccDNA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

HAND transcription factors cooperatively specify the aorta and pulmonary trunk

Congenital heart defects (CHDs) affecting the cardiac outflow tract (OFT) constitute a significant cause of morbidity and mortality. The OFT develops from migratory cell populations which include the cardiac neural crest cells (cNCCs) and secondary heart field (SHF) derived myocardium and endocardium. The related transcription factors HAND1 and HAND2 have been implicated in human CHDs involving the OFT. Although Hand1 is expressed within the OFT, Hand1 NCC-specific conditional knockout mice (H1CKOs) are viable. Here we show that these H1CKOs present a low penetrance of OFT phenotypes, whereas SHF-specific Hand1 ablation does not reveal any cardiac phenotypes. Further, HAND1 and HAND2 appear functionally redundant within the cNCCs, as a reduction/ablation of Hand2 on an NCC-specific H1CKO background causes pronounced OFT defects. Double conditional Hand1 and Hand2 NCC knockouts exhibit persistent truncus arteriosus (PTA) with 100% penetrance. NCC lineage-tracing and Sema3c in situ mRNA expression reveal that Sema3c-expressing cells are mis-localized, resulting in a malformed septal bridge within the OFTs of H1CKO;H2CKO embryos. Interestingly, Hand1 and Hand2 also genetically interact within the SHF, as SHF H1CKOs on a heterozygous Hand2 background exhibit Ventricular Septal Defects (VSDs) with incomplete penetrance. Previously, we identified a BMP, HAND2, and GATA-dependent Hand1 OFT enhancer sufficient to drive reporter gene expression within the nascent OFT and aorta. Using these transcription inputs as a probe, we identify a novel Hand2 OFT enhancer, suggesting that a conserved BMP-GATA dependent mechanism transcriptionally regulates both HAND factors. These findings support the hypothesis that HAND factors interpret BMP signaling within the cNCCs to cooperatively coordinate OFT morphogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Systematic identification of transcriptional activation domains from non-transcription factor proteins in plants and yeast

Transcription factors can promote gene expression through activation domains. Whole-genome screens have systematically mapped activation domains in transcription factors but not in non-transcription factor proteins (e.g., chromatin regulators and coactivators). To fill this knowledge gap, we employed the activation domain predictor PADDLE to analyze the proteomes of Arabidopsis thaliana and Saccharomyces cerevisiae. We screened 18,000 predicted activation domains from >800 non-transcription factor genes in both species, confirming that 89% of candidate proteins contain active fragments. Our work enables the annotation of hundreds of nuclear proteins as putative coactivators, many of which have never been ascribed any function in plants. Analysis of peptide sequence compositions reveals how the distribution of key amino acids dictates activity. Finally, we validated short, "universal" activation domains with comparable performance to state-of-the-art activation domains used for genome engineering. Our approach enables the genome-wide discovery and annotation of activation domains that can function across diverse eukaryotes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A riboswitch separated from its ribosome-binding site still regulates translation

Abstract Riboswitches regulate downstream gene expression by binding cellular metabolites. Regulation of translation initiation by riboswitches is posited to occur by metabolite-mediated sequestration of the Shine-Dalgarno sequence (SDS), causing bypass by the ribosome. Recently, we solved a co-crystal structure of a prequeuosine1-sensing riboswitch from Carnobacterium antarcticum that binds two metabolites in a single pocket. The structure revealed that the second nucleotide within the gene-regulatory SDS, G34, engages in a crystal contact, obscuring the molecular basis of gene regulation. Here, we report a co-crystal structure wherein C10 pairs with G34. However, molecular dynamics simulations reveal quick dissolution of the pair, which fails to reform. Functional and chemical probing assays inside live bacterial cells corroborate the dispensability of the C10–G34 pair in gene regulation, leading to the hypothesis that the compact pseudoknot fold is sufficient for translation attenuation. Remarkably, the C. antarcticum aptamer retained significant gene-regulatory activity when uncoupled from the SDS using unstructured spacers up to 10 nucleotides away from the riboswitch—akin to steric-blocking employed by sRNAs. Accordingly, our work reveals that the RNA fold regulates translation without SDS sequestration, expanding known riboswitch-mediated gene-regulatory mechanisms. The results infer that riboswitches exist wherein the SDS is not embedded inside a stable fold.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

AtDGCR14L contributes to salt-stress tolerance via regulating pre-mRNA splicing in Arabidopsis

In plants, the pre-mRNA alternative splicing has been demonstrated to be a crucial tier that regulates gene expression in response to salt stress. However, the underlying mechanisms remain elusive. Here, in this study, we studied the roles of DIGEORGE-SYNDROME CRITICAL REGION 14-like (AtDGCR14L) in regulating pre-mRNA splicing and salt stress tolerance. We discovered that Arabidopsis AtDGCR14L is required for maintaining plant salt stress tolerance and the constitutively spliced and active isoforms of important stress- and/or abscisic acid (ABA)-responsive genes. We also identified the interaction between AtDGCR14L and splicing factor U1-70k, which needs a highly conserved three amino acid (TWG) motif in DGCR14. Different from wild-type AtDGCR14L, the overexpression of TWG-substituted AtDGCR14L mutant did not change salt stress tolerance or pre-mRNA splicing of stress/ABA-responsive genes. Additionally, SWITCH3A (SWI3A) is a core subunit of the SWI/SUCROSE NONFERMENTING (SWI/SNF) chromatin-remodeling complexes. We found that SWI3A, whose splicing depends on AtDGCR14L, actively enhances salt stress tolerance. These results revealed that AtDGCR14L may play an essential role in crosstalk between plant salt-stress response and pre-mRNA splicing mechanisms. We also unveiled the potential role of SWI3A in controlling salt stress tolerance. The TWG motif in the intrinsically disordered region of AtDGCR14L is highly conserved and crucial for DGCR14 functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The human aortic endothelium undergoes dose-dependent DNA methylation in response to transient hyperglycemia

Glycemic control is a strong predictor of long-term cardiovascular risk in patients with diabetes mellitus, and poor glycemic control influences long-term risk of cardiovascular disease even decades after optimal medical management. This phenomenon, termed glycemic memory, has been proposed to occur due to stable programs of cardiac and endothelial cell gene expression. This transcriptional remodeling has been shown to occur in the vascular endothelium through a yet undefined mechanism of cellular reprogramming.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Polysaccharide niche partitioning of distinct Polaribacter clades during North Sea spring algal blooms

Massive releases of organic substrates during marine algal blooms trigger growth of many clades of heterotrophic bacteria. Algal polysaccharides represent the most diverse and structurally complex class of these substrates, yet their role in shaping the microbial community composition is poorly understood. We investigated, whether polysaccharide utilization capabilities contribute to niche differentiation of Polaribacter spp. (class Flavobacteriia; known to include relevant polysaccharidedegraders) that were abundant during 2009–2012 spring algal blooms in the southern North Sea. We identified six distinct Polaribacter clades using phylogenetic and phylogenomic analyses, quantified their abundances via fluorescence in situ hybridization, compared metagenome-assembled genomes, and assessed in situ gene expression using metaproteomics. Four clades with distinct polysaccharide niches were dominating. Polaribacter 2-a comprised typical first responders featuring small genomes with limited polysaccharide utilization capacities. Polaribacter 3-a were abundant only in 2010 and possessed a distinct sulfated α-glucoronomannan degradation potential. Polaribacter 3-b responded late in blooms and had the capacity to utilize sulfated xylan. Polaribacter 1-a featured high numbers of glycan degradation genes and were particularly abundant following Chattonella algae blooms. These results support the hypothesis that sympatric Polaribacter clades occupy distinct glycan niches during North Sea spring algal blooms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

X-ray Crystallographic Study of Preferred Spacing by the NF-κB p50 Homodimer on κB DNA

Though originally characterized as an inactive or transcriptionally repressive factor, the NF-κB p50 homodimer has become appreciated as a physiologically relevant driver of specific target gene expression. By virtue of its low affinity for cytoplasmic IκB protein inhibitors, p50 accumulates in the nucleus of resting cells, where it is a binding target for the transcriptional co-activator IκBζ. In this study, we employed X-ray crystallography to analyze the structure of the p50 homodimer on κB DNA from the promoters of human interleukin-6 (IL-6) and neutrophil-gelatinase-associated lipocalin (NGAL) genes, both of which respond to IκBζ. The NF-κB p50 homodimer binds 11-bp on IL-6 κB DNA, while, on NGAL κB DNA, the spacing is 12-bp. This begs the question: what DNA binding mode is preferred by NF-κB p50 homodimer? To address this, we engineered a “Test” κB-like DNA containing the core sequence 5'-GGGGAATTCCCC-3' and determined its X-ray crystal structure in complex with p50. This revealed that, when presented with multiple options, NF-κB p50 homodimer prefers to bind 11-bp, which necessarily imposes asymmetry on the complex despite the symmetry inherent in both the protein and its target DNA, and that the p50 dimerization domain can contact DNA via distinct modes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A t FUT4 and A t FUT6 Are Arabinofuranose-Specific Fucosyltransferases

The bulk of plant biomass is comprised of plant cell walls, which are complex polymeric networks, composed of diverse polysaccharides, proteins, polyphenolics, and hydroxyproline-rich glycoproteins (HRGPs). Glycosyltransferases (GTs) work together to synthesize the saccharide components of the plant cell wall. The Arabidopsis thaliana fucosyltransferases (FUTs), At FUT4, and At FUT6, are members of the plant-specific GT family 37 (GT37). At FUT4 and At FUT6 transfer fucose (Fuc) onto arabinose (Ara) residues of arabinogalactan (AG) proteins (AGPs) and have been postulated to be non-redundant AGP-specific FUTs. At FUT4 and At FUT6 were recombinantly expressed in mammalian HEK293 cells and purified for biochemical analysis. We report an updated understanding on the specificities of At FUT4 and At FUT6 that are involved in the synthesis of wall localized AGPs. Our findings suggest that they are selective enzymes that can utilize various arabinogalactan (AG)-like and non-AG-like oligosaccharide acceptors, and only require a free, terminal arabinofuranose. We also report with GUS promoter-reporter gene studies that AtFUT4 and AtFUT6 gene expression is sub-localized in different parts of developing A. thaliana roots.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Abstract Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Azodi, Christina B.↗