SEARCH · Engineering Papers

Results for “DNA 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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Accelerated epigenetic aging as a risk factor for chronic obstructive pulmonary disease and decreased lung function in two prospective cohort studies

Chronic obstructive pulmonary disease (COPD) is a frequent diagnosis in older individuals and contributor to global morbidity and mortality. Given the link between lung disease and aging, we need to understand how molecular indicators of aging relate to lung function and disease. Using data from the population-based KORA (Cooperative Health Research in the Region of Augsburg) surveys, we associated baseline epigenetic (DNA methylation) age acceleration with incident COPD and lung function. Models were adjusted for age, sex, smoking, height, weight, and baseline lung disease as appropriate. Associations were replicated in the Normative Aging Study. Of 770 KORA participants, 131 developed incident COPD over 7 years. Baseline accelerated epigenetic aging was significantly associated with incident COPD. The change in age acceleration (follow-up – baseline) was more strongly associated with COPD than baseline aging alone. The association between the change in age acceleration between baseline and follow-up and incident COPD replicated in the Normative Aging Study. Associations with spirometric lung function parameters were weaker than those with COPD, but a meta-analysis of both cohorts provide suggestive evidence of associations. Accelerated epigenetic aging, both baseline measures and changes over time, may be a risk factor for COPD and reduced lung function.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Gene content evolution in the arthropods

Background: Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods. Results: Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception. Conclusions: These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Arthropods↗

DOE OSTI2020

Biochemical and structural basis for YTH domain of human YTHDC1 binding to methylated adenine in DNA

The recently characterized mammalian writer (methyltransferase) and eraser (demethylase) of the DNA N6-methyladenine (N6mA) methyl mark act on single-stranded (ss) and transiently-unpaired DNA. As YTH domain-containing proteins bind N6mA-containing RNA in mammalian cells, we investigated whether mammalian YTH domains are also methyl mark readers of N6mA DNA. Here, we show that the YTH domain of YTHDC1 (known to localize in the nucleus) binds ssDNA containing N6mA, with a 10 nM dissociation constant. This binding is stronger by a factor of 5 than in an RNA context, tested under the same conditions. However, the YTH domains of YTHDF2 and YTHDF1 (predominantly cytoplasmic) exhibited the opposite effect with ~1.5–2× stronger binding to ssRNA containing N6mA than to the corresponding DNA. We determined two structures of the YTH domain of YTHDC1 in complex with N6mA-containing ssDNA, which illustrated that YTHDC1 binds the methylated adenine in a single-stranded region flanked by duplexed DNA. We discuss the hypothesis that the writer-reader-eraser of N6mA-containining ssDNA is associated with maintaining genome stability. Structural comparison of YTH and SRA domains (the latter a DNA 5-methylcytosine reader) revealed them to be diverse members of a larger family of DNA/RNA modification readers, apparently having originated from bacterial modification-dependent restriction enzymes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Mechanistic Evaluation of Black Cohosh Extract-Induced Genotoxicity in Human Cells

Abstract Black cohosh extract (BCE) is marketed to women as an alternative to hormone replacement therapy for alleviating menopausal symptoms. Previous studies by the National Toxicology Program revealed that BCE induced micronuclei (MN) and a nonregenerative macrocytic anemia in rats and mice, likely caused by disruption of the folate metabolism pathway. Additional work using TK6 cells showed that BCE induced aneugenicity by destabilizing microtubules. In the present study, BCE-induced MN were confirmed in TK6 and HepG2 cells. We then evaluated BCE-induced DNA damage using the comet assay at multiple time points (0.5–24 h). Following a 0.5-h exposure, BCE induced significant, concentration-dependent increases in %tail DNA in TK6 cells only. Although DNA damage decreased in TK6 cells over time, likely due to repair, small but statistically significant levels of DNA damage were observed after 2 and 4 h exposures to 250 µg/ml BCE. A G1/S arrest in TK6 cells exposed to 125 µg/ml BCE (24 h) was accompanied by apoptosis and increased expression of γH2A.X, p-Chk1, p-Chk2, p53, and p21. Conditioning TK6 cells to physiological levels of folic acid (120 nM) did not increase the sensitivity of cells to BCE-induced DNA damage. BCE did not alter global DNA methylation in TK6 and HepG2 cells cultured in standard medium. Our results suggest that BCE induces acute DNA strand breaks which are quickly repaired in TK6 cells, whereas DNA damage seen at 4 and 24 h may reflect apoptosis. The present study supports that BCE is genotoxic mainly by inducing MN with an aneugenic mode of action.

Seo, Ji-Eun↗

DOE OSTI2021

RACGAP1 modulates ECT2-Dependent mitochondrial quality control to drive breast cancer metastasis

Most cancer deaths are due to the colonization of tumor cells in distant organs. More evidence indicates that overexpression of RACGAP1 plays a critical role in cancer metastasis. However, the underlying mechanism still remains poorly understood. Here we found that RACGAP1 promoted breast cancer metastasis through regulating mitochondrial quality control. Overexpression of RACGAP1 in breast cancer cells led to the fragmentation of mitochondria, increased mitophagy intensity, mitochondrial turnover, and aerobic glycolysis ATP production. We showed that RACGAP1 promoted mitochondrial fission through recruiting ECT2 during anaphase and subsequently had activated ERK-DRP1 pathway. We further demonstrated the phosphorylation of RACGAP1 is essential for its ability of binding with ECT2 and its downstream effects. RACGAP1 overexpression also increased the expression of PGC-1a, a key mitochondrial biogenesis regulator, presumably by the increased mitophagy intensity induced by RACGAP1. PGC-1a increased the enrichment of DNMT1 in mitochondria, mitochondrial DNMT1 augmented mitochondrial DNA methylation and upregulated mitochondrial genome transcription. Our data indicated that RACGAP1 simultaneously facilitated mitophagy and mitochondrial biogenesis through regulating DRP1 phosphorylation and PGC-1a expression, eventually improved mitochondrial quality control in breast cancer cells. Our study provided a new angle in understanding the RACGAP1-overexpression related malignancy in breast cancer patients.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Differential localization patterns of pyruvate kinase isoforms in murine naïve, formative, and primed pluripotent states

Highlights: • PKM1/2 protein abundance is greater in formative mEpiLCs compared to naïve mESCs or primed mEpiSCs. The ratio of PKM1/2 is maintained. • PKM1/2 are localized in both nuclear and cytoplasmic regions relative to GAPDH and OCT4 across the pluripotent continuum. • PKM1 localization is strongly correlated to OCT4 localization, and moderately to GAPDH in formative mEpiLCs. • PKM1/2 localization is moderately correlated to OCT4 and GAPDH localization in naïve mESCs. • PKM1/2 localization are moderately correlated to GAPDH localization in primed mEpiSCs. Mouse embryonic stem cells (mESCs) and mouse epiblast stem cells (mEpiSCs) represent opposite ends of the pluripotency continuum, referred to as naïve and primed pluripotent states, respectively. These divergent pluripotent states differ in several ways, including growth factor requirements, transcription factor expression, DNA methylation patterns, and metabolic profiles. Naïve cells employ both glycolysis and oxidative phosphorylation (OXPHOS), whereas primed cells preferentially utilize aerobic glycolysis, a trait shared with cancer cells referred to as the Warburg Effect. Until recently, metabolism has been regarded as a by-product of cell fate, however, evidence now supports metabolism as being a driver of stem cell state and fate decisions. Pyruvate kinase muscle isoforms (PKM1 and PKM2) are important for generating and maintaining pluripotent stem cells (PSCs) and mediating the Warburg Effect. Both isoforms catalyze the final, rate limiting step of glycolysis, generating adenosine triphosphate and pyruvate, however, the precise role(s) of PKM1/2 in naïve and primed pluripotency is not well understood. The primary objective of this study was to characterize the cellular expression and localization patterns of PKM1 and PKM2 in mESCs, chemically transitioned epiblast-like cells (mEpiLCs) representing formative pluripotency, and mEpiSCs using immunoblotting and confocal microscopy. The results indicate that PKM1 and PKM2 are not only localized to the cytoplasm, but also accumulate in differential subnuclear regions of mESC, mEpiLCs, and mEpiSCs as determined by a quantitative confocal microscopy employing orthogonal projections and airyscan processing. Importantly, we discovered that the subnuclear localization of PKM1/2 changes during the transition from mESCs, mEpiLCs, and mEpiSCs. Finally, we have comprehensively validated the appropriateness and power of the Pearson's correlation coefficient and Manders's overlap coefficient for assessing nuclear and cytoplasmic protein colocalization in PSCs by immunofluorescence confocal microscopy. We propose that nuclear PKM1/2 may assist with distinct pluripotency state maintenance and lineage priming by non-canonical mechanisms. These results advance our understanding of the overall mechanisms controlling naïve, formative, and primed pluripotency.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

DNA transfer between two different species mediated by heterologous cell fusion in Clostridium coculture

Investigations of natural multispecies microbiomes and synthetic microbial cocultures are attracting renewed interest for their potential application in biotechnology, ecology, and medical fields. Previously, we have shown the syntrophic coculture of C. acetobutylicum and C. ljungdahlii undergoes heterologous cell-to-cell fusion, which facilitates the exchange of cytoplasmic protein and RNA between the two organisms. We now show that heterologous cell fusion between the two Clostridium organisms can facilitate the exchange of DNA. By applying selective pressures to this coculture system, we isolated clones of wild-type C. acetobutylicum which acquired the erythromycin resistance (erm) gene from the C. ljungdahlii strain carrying a plasmid with the erm gene. Single-molecule real-time sequencing revealed that the erm gene was integrated into the genome in a mosaic fashion. Our data also support the persistence of hybrid C. acetobutylicum / C. ljungdahlii cells displaying hybrid DNA-methylation patterns.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2014

Does Simulated Weightlessness Alter the Methylation Status of Gene Promoters During Bone Remodeling?

The skeleton interacts with its environment in a way that resembles a mechanostat - through a controlled process of bone remodeling, namely local formation and resorption, to maintain a healthy structure. During weightlessness, astronauts lose structure in weight-bearing bones due to decreased formation by osteoblasts and increased resorption by osteoclasts. In contrast, increased mechanical loading through exercise targets bone remodeling to remove and repair microdamage, improving structural integrity. In fact, recent advances in astronaut exercise regimens have prevented the deleterious changes in skeletal structure during spaceflight. However, knowledge of the molecular underpinnings of the skeletal response to spaceflight and to mechanical stimulation is limited. We propose that epigenetic modification, specifically DNA methylation, may influence osteoblast differentiation and activity during spaceflight and exercise. We hypothesize that simulated weightlessness hypermethylates pro-osteoblastogenic gene promoters and decreases expression of osteoblastogenic genes. Oppositely, we hypothesize that mechanical loading hypomethylates pro-osteoblastogenic gene promoters and increases expression of osteoblastogenic genes.

Stevick, Rebecca J.↗

DOE OSTI2020

The genomic and epigenomic evolutionary history of papillary renal cell carcinomas

Intratumor heterogeneity (ITH) and tumor evolution have been well described for clear cell renal cell carcinomas (ccRCC), but they are less studied for other kidney cancer subtypes. Here we investigate ITH and clonal evolution of papillary renal cell carcinoma (pRCC) and rarer kidney cancer subtypes, integrating whole-genome sequencing and DNA methylation data. In 29 tumors, up to 10 samples from the center to the periphery of each tumor, and metastatic samples in 2 cases, enable phylogenetic analysis of spatial features of clonal expansion, which shows congruent patterns of genomic and epigenomic evolution. In contrast to previous studies of ccRCC, in pRCC, driver gene mutations and most arm-level somatic copy number alterations (SCNAs) are clonal. These findings suggest that a single biopsy would be sufficient to identify the important genetic drivers and that targeting large-scale SCNAs may improve pRCC treatment, which is currently poor. While type 1 pRCC displays near absence of structural variants (SVs), the more aggressive type 2 pRCC and the rarer subtypes have numerous SVs, which should be pursued for prognostic significance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Expanded encyclopaedias of DNA elements in the human and mouse genomes

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Mycena genomes resolve the evolution of fungal bioluminescence

Mushroom-forming fungi in the order Agaricales represent an independent origin of bioluminescence in the tree of life; yet the diversity, evolutionary history, and timing of the origin of fungal luciferases remain elusive. We sequenced the genomes and transcriptomes of five bonnet mushroom species ( Mycena spp.), a diverse lineage comprising the majority of bioluminescent fungi. Two species with haploid genome assemblies ~150 Mb are among the largest in Agaricales, and we found that a variety of repeats between Mycena species were differentially mediated by DNA methylation. We show that bioluminescence evolved in the last common ancestor of mycenoid and the marasmioid clade of Agaricales and was maintained through at least 160 million years of evolution. Analyses of synteny across genomes of bioluminescent species resolved how the luciferase cluster was derived by duplication and translocation, frequently rearranged and lost in most Mycena species, but conserved in the Armillaria lineage. Luciferase cluster members were coexpressed across developmental stages, with the highest expression in fruiting body caps and stipes, suggesting fruiting-related adaptive functions. Our results contribute to understanding a de novo origin of bioluminescence and the corresponding gene cluster in a diverse group of enigmatic fungal species.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Alternative splicing and allosteric regulation modulate the chromatin binding of UHRF1

UHRF1 is an important epigenetic regulator associated with apoptosis and tumour development. It is a multidomain protein that integrates readout of different histone modification states and DNA methylation with enzymatic histone ubiquitylation activity. Emerging evidence indicates that the chromatin-binding and enzymatic modules of UHRF1 do not act in isolation but interplay in a coordinated and regulated manner. Here, we compared two splicing variants (V1, V2) of murine UHRF1 (mUHRF1) with human UHRF1 (hUHRF1). We show that insertion of nine amino acids in a linker region connecting the different TTD and PHD histone modification-binding domains causes distinct H3K9me3-binding behaviour of mUHRF1 V1. Structural analysis suggests that in mUHRF1 V1, in contrast to V2 and hUHRF1, the linker is anchored in a surface groove of the TTD domain, resulting in creation of a coupled TTD-PHD module. This establishes multivalent, synergistic H3-tail binding causing distinct cellular localization and enhanced H3K9me3-nucleosome ubiquitylation activity. In contrast to hUHRF1, H3K9me3-binding of the murine proteins is not allosterically regulated by phosphatidylinositol 5-phosphate that interacts with a separate less-conserved polybasic linker region of the protein. Our results highlight the importance of flexible linkers in regulating multidomain chromatin binding proteins and point to divergent evolution of their regulation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Transcriptome and DNA methylome divergence of inflorescence development between 2 ecotypes in Panicum hallii

The morphological diversity of the inflorescence determines flower and seed production, which is critical for plant adaptation. Hall's panicgrass (Panicum hallii, P. hallii) is a wild perennial grass that has been developed as a model to study perennial grass biology and adaptive evolution. Highly divergent inflorescences have evolved between the 2 major ecotypes in P. hallii, the upland ecotype (P. hallii var hallii, HAL2 genotype) with compact inflorescence and large seed and the lowland ecotype (P. hallii var filipes, FIL2 genotype) with an open inflorescence and small seed. Here we conducted a comparative analysis of the transcriptome and DNA methylome, an epigenetic mark that influences gene expression regulation, across different stages of inflorescence development using genomic references for each ecotype. Global transcriptome analysis of differentially expressed genes (DEGs) and co-expression modules underlying the inflorescence divergence revealed the potential role of cytokinin signaling in heterochronic changes. Comparing DNA methylome profiles revealed a remarkable level of differential DNA methylation associated with the evolution of P. hallii inflorescence. We found that a large proportion of differentially methylated regions (DMRs) were located in the flanking regulatory regions of genes. Intriguingly, we observed a substantial bias of CHH hypermethylation in the promoters of FIL2 genes. The integration of DEGs, DMRs, and $K_a$/$K_s$ ratio results characterized the evolutionary features of DMR-associated DEGs that contribute to the divergence of the P. hallii inflorescence. This study provides insights into the transcriptome and epigenetic landscape of inflorescence divergence in P. hallii and a genomic resource for perennial grass biology.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

35S rDNA evolution↗

DOE OSTI2022

Unified epigenomic, transcriptomic, proteomic, and metabolomic taxonomy of Alzheimer’s disease progression and heterogeneity

Alzheimer’s disease (AD) is a heterogeneous disorder with abnormalities in multiple biological domains. In an advanced machine learning analysis of postmortem brain and in vivo blood multi-omics molecular data ( N = 1863), we integrated epigenomic, transcriptomic, proteomic, and metabolomic profiles into a multilevel biological AD taxonomy. We obtained a personalized multilevel molecular index of AD dementia progression that predicts severity of neuropathologies, and identified three robust molecular-based subtypes that explain much of the pathologic and clinical heterogeneity of AD. These subtypes present distinct patterns of alteration in DNA methylation, RNA, proteins, and metabolites, identifiable in the brain and subsequently in blood. In addition, the genetic variations that predispose to the various AD subtypes in brain predict distinct spatial patterns of alteration in cell types, suggesting a unique influence of each putative AD variant on neuropathological mechanisms. These observations support that an individually tailored multi-omics molecular taxonomy of AD may represent distinct targets for preventive or treatment interventions.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Pathogenic SIV infection is associated with acceleration of epigenetic age in rhesus macaques

HIV infection accelerates biological aging, but the contribution of the host’s age to this process is unknown. We investigated the influence of SIV infection in macaques (SIVmac) on the risk of comorbidities and aging in young and old rhesus macaques (RMs) by assessing pathogenesis markers, DNA methylation–based epigenetic age (EA), and EA acceleration (EAA) in blood and tissues. Initially, upon SIV infection, the young RMs showed greater resilience to CD4 + T cell depletion, better control of T cell activation, hypercoagulation, and excessive inflammation, yet this resilience was progressively lost in the advanced stages of infection. During the late stages of infection, the young RMs, but not the aged ones, showed an increase in EA in PBMCs; also, EAA in the cerebellum and heart of young RMs was higher compared with old RMs. SIV infection was more pathogenic in aged animals in early stages, leading to a more rapid disease progression; however, accelerated aging mostly affected young animals, so that the levels of multiple key pathogenesis markers in the young RMs converged toward those specific to aged ones in the late stages of infection. We conclude that SIV infection–driven age acceleration is tissue specific, and that host age influences the susceptibility of different tissues to enhanced aging.

Research & Experimental Medicine↗

DOE OSTI2022

Comparative genomic and transcriptomic analyses of trans-kingdom pathogen Fusarium solani species complex reveal degrees of compartmentalization

Background: The Fusarium solani species complex (FSSC) comprises fungal pathogens responsible for mortality in a diverse range of animals and plants, but their genome diversity and transcriptome responses in animal pathogenicity remain to be elucidated. We sequenced, assembled and annotated six chromosome-level FSSC clade 3 genomes of aquatic animal and plant host origins. We established a pathosystem and investigated the expression data of F. falciforme and F. keratoplasticum in Chinese softshell turtle (Pelodiscus sinensis) host. Results: Comparative analyses between the FSSC genomes revealed a spectrum of conservation patterns in chromosomes categorised into three compartments: core, fast-core (FC), and lineage-specific (LS). LS chromosomes contribute to variations in genomes size, with up to 42.2% of variations between F. vanettenii strains. Each chromosome compartment varied in structural architectures, with FC and LS chromosomes contain higher proportions of repetitive elements with genes enriched in functions related to pathogenicity and niche expansion. We identified differences in both selection in the coding sequences and DNA methylation levels between genome features and chromosome compartments which suggest a multi-speed evolution that can be traced back to the last common ancestor of Fusarium. We further demonstrated that F. falciforme and F. keratoplasticum are opportunistic pathogens by inoculating P. sinensis eggs and identified differentially expressed genes also associated with plant pathogenicity. These included the most upregulated genes encoding the CFEM (Common in Fungal Extracellular Membrane) domain. Conclusions: The high-quality genome assemblies provided new insights into the evolution of FSSC chromosomes, which also serve as a resource for studies of fungal genome evolution and pathogenesis. This study also establishes an animal model for fungal pathogens of trans-kingdom hosts.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Temporal change in chromatin accessibility predicts regulators of nodulation in Medicago truncatula

Symbiotic associations between bacteria and leguminous plants lead to the formation of root nodules that fix nitrogen needed for sustainable agricultural systems. Symbiosis triggers extensive genome and transcriptome remodeling in the plant, yet an integrated understanding of the extent of chromatin changes and transcriptional networks that functionally regulate gene expression associated with symbiosis remains poorly understood. In particular, analyses of early temporal events driving this symbiosis have only captured correlative relationships between regulators and targets at mRNA level. Here, we characterize changes in transcriptome and chromatin accessibility in the model legume Medicago truncatula, in response to rhizobial signals that trigger the formation of root nodules. We profiled the temporal chromatin accessibility (ATAC-seq) and transcriptome (RNA-seq) dynamics of M. truncatula roots treated with bacterial small molecules called lipo-chitooligosaccharides that trigger host symbiotic pathways of nodule development. Using a novel approach, dynamic regulatory module networks, we integrated ATAC-seq and RNA-seq time courses to predict cis-regulatory elements and transcription factors that most significantly contribute to transcriptomic changes associated with symbiosis. Regulators involved in auxin (IAA4-5, SHY2), ethylene (EIN3, ERF1), and abscisic acid (ABI5) hormone response, as well as histone and DNA methylation (IBM1), emerged among those most predictive of transcriptome dynamics. RNAi-based knockdown of EIN3 and ERF1 reduced nodule number in M. truncatula validating the role of these predicted regulators in symbiosis between legumes and rhizobia. Our transcriptomic and chromatin accessibility datasets provide a valuable resource to understand the gene regulatory programs controlling the early stages of the dynamic process of symbiosis. The regulators identified provide potential targets for future experimental validation, and the engineering of nodulation in species is unable to establish that symbiosis naturally.

59 BASIC BIOLOGICAL SCIENCES↗