Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gene Expression Profiling”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Translational Line of Sight: a multi-omics longitudinal study of the murine retinal response to spaceflight hazard analogs

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗
NASA NTRSDate not supplied

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and/or Simulated Microgravity

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗
DOE OSTI2026

Carbon source-induced changes in the transcriptional landscape of Novacetimonas hansenii

Novacetimonas hansenii is a model organism for bacterial cellulose (BC) production, a biopolymer with broad industrial applications. To better understand its carbon source-dependent metabolism, we conducted the first time-resolved transcriptomic analysis of N. hansenii ATCC 23769 grown on five carbon sources (glucose, fructose, mannitol, galactose, and sucrose) across five time points (24–120h). Glucose, fructose, and mannitol supported growth and efficient BC production, whereas galactose and sucrose resulted in poor yields, likely due to limited transport and enzymatic capacity. Solid-state NMR revealed that cellulose crystallinity and polymorphism were carbon source-dependent, with reduced bcsD expression correlating with lower crystallinity. Transcriptomic profiling showed distinct metabolic responses. Cultures containing glucose and mannitol activated core metabolic pathways, fructose induced a delayed adaptive shift, and galactose and sucrose triggered stress-related gene expression. These results highlight the metabolic plasticity of N. hansenii and offer insights for enhancing BC production through carbon source optomization.

Bacterial cellulose↗
DOE OSTI2021

SMAD4 target genes are part of a transcriptional network that integrates the response to BMP and SHH signaling during early limb bud patterning

SMAD4 regulates gene expression in response to BMP and TGFβ signal transduction, and is required for diverse morphogenetic processes, but its target genes have remained largely elusive. Here, we identify the SMAD4 target genes in mouse limb buds using an epitope-tagged Smad4 allele for ChIP-seq analysis in combination with transcription profiling. This analysis shows that SMAD4 predominantly mediates BMP signal transduction during early limb bud development. Unexpectedly, the expression of cholesterol biosynthesis enzymes is precociously downregulated and intracellular cholesterol levels are reduced in Smad4-deficient limb bud mesenchymal progenitors. Most importantly, our analysis reveals a predominant function of SMAD4 in upregulating target genes in the anterior limb bud mesenchyme. Analysis of differentially expressed genes shared between Smad4- and Shh-deficient limb buds corroborates this function of SMAD4 and also reveals the repressive effect of SMAD4 on posterior genes that are upregulated in response to SHH signaling. This analysis uncovers opposing trans-regulatory inputs from SHH- and SMAD4-mediated BMP signal transduction on anterior and posterior gene expression during the digit patterning and outgrowth in early limb buds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Spatial distribution of marker gene activity in the mouse lung during alveolarization

Here, this data is a curated collection of visual images of gene expression patterns from the pre- and post-natal mouse lung, accompanied by associated mRNA probe sequences and RNA-Seq expression profiles. Mammalian lungs undergo significant growth and cellular differentiation before and after the transition to breathing air. Documenting normal lung development is an important step in understanding abnormal lung development, as well as the challenges faced during a preterm birth. Images in this dataset indicate the spatial distribution of mRNA transcripts for over 500 different genes that are active during lung development, as initially determined via RNA-Seq. Images were systematically acquired using high-throughput in situ hybridization with non-radioactive digoxigenin-labeled mRNA probes across mouse lungs from developmental time points E16.5, E18.5, P7, and P28.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Compartmentalized microbial composition, oxygen gradients and nitrogen fixation in the gut of Odontotaenius disjunctus

Coarse woody debris is an important biomass pool in forest ecosystems that numerous groups of insects have evolved to take advantage of. These insects are ecologically important and represent useful natural analogs for biomass to biofuel conversion. Using a range of molecular approaches combined with microelectrode measurements of oxygen, here we have characterized the gut microbiome and physiology of Odontotaenius disjunctus, a wood-feeding beetle native to the eastern United States. We hypothesized that morphological and physiological differences among gut regions would correspond to distinct microbial populations and activities. In fact, significantly different communities were found in the foregut (FG), midgut (MG)/posterior hindgut (PHG) and anterior hindgut (AHG), with Actinobacteria and Rhizobiales being more abundant toward the FG and PHG. Conversely, fermentative bacteria such as Bacteroidetes and Clostridia were more abundant in the AHG, and also the sole region where methanogenic Archaea were detected. Although each gut region possessed an anaerobic core, micron-scale profiling identified radial gradients in oxygen concentration in all regions. Nitrogen fixation was confirmed by 15 N 2 incorporation, and nitrogenase gene (nifH) expression was greatest in the AHG. Phylogenetic analysis of nifH identified the most abundant transcript as related to Ni–Fe nitrogenase of a Bacteroidetes species, Paludibacter propionicigenes. Overall, we demonstrate not only a compartmentalized microbiome in this beetle digestive tract but also sharp oxygen gradients that may permit aerobic and anaerobic metabolism to occur within the same regions in close proximity. We provide evidence for the microbial fixation of N 2 that is important for this beetle to subsist on woody biomass.

Ceja-Navarro, Javier A.↗
DOE OSTI2022

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass)

Background: MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass ( Panicum virgatum ), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass. Results: Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession “Alamo”, upland accession “Dacotah”, and two upland/lowland F 1 hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors. Conclusion: Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Membrane lipid and expression responses of Saccharolobus islandicus REY15A to acid and cold stress

Archaea adjust the number of cyclopentane rings in their glycerol dibiphytanyl glycerol tetraether (GDGT) membrane lipids as a homeostatic response to environmental stressors such as temperature, pH, and energy availability shifts. However, archaeal expression patterns that correspond with changes in GDGT composition are less understood. Here we characterize the acid and cold stress responses of the thermoacidophilic crenarchaeon Saccharolobus islandicus REY15A using growth rates, core GDGT lipid profiles, transcriptomics and proteomics. We show that both stressors result in impaired growth, lower average GDGT cyclization, and differences in gene and protein expression. Transcription data revealed differential expression of the GDGT ring synthase grsB in response to both acid stress and cold stress. Although the GDGT ring synthase encoded by grsB forms highly cyclized GDGTs with ≥5 ring moieties, S. islandicus grsB upregulation under acidic pH conditions did not correspond with increased abundances of highly cyclized GDGTs. Our observations highlight the inability to predict GDGT changes from transcription data alone. Broader analysis of transcriptomic data revealed that S. islandicus differentially expresses many of the same transcripts in response to both acid and cold stress. These included upregulation of several biosynthetic pathways and downregulation of oxidative phosphorylation and motility. Transcript responses specific to either of the two stressors tested here included upregulation of genes related to proton pumping and molecular turnover in acid stress conditions and upregulation of transposases in cold stress conditions. Overall, our study provides a comprehensive understanding of the GDGT modifications and differential expression characteristic of the acid stress and cold stress responses in S. islandicus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Aristolochic acid-induced genotoxicity and toxicogenomic changes in rodents

Aristolochic acid (AA) is a group of structurally related nitrophenanthrene carboxylic acids found in many plants that are widely used by many cultures as traditional herbal medicines. AA is a causative agent for Chinese herbs nephropathy, a term replaced later by AA nephropathy. Evidence indicates that AA is nephrotoxic, genotoxic, and carcinogenic in humans; and it also induces tumors in the forestomach, kidney, renal pelvis, urinary bladder, and lung of rats and mice. Therefore, plants containing AA have been classified as carcinogenic to humans (Group 1) by the International Agency for Research on Cancer. In our laboratories, we have conducted a series of genotoxicity and toxicogenomic studies in the rats exposed to AA of 0.1–10 mg/kg for 12 weeks. Our results demonstrated that AA treatments induced DNA adducts and mutations in the kidney, liver, and spleen of rats, as well as significant alteration of gene expression in both its target and nontarget tissues. AA treatments altered mutagenesis- or carcinogenesis-related microRNA expression in rat kidney and resulted in significant changes in protein expression profiling. We also applied benchmark dose (BMD) modeling to the 3-month AA-induced genotoxicity data. The obtained BMDL10 (the lower 95% confidence interval of the BMD 10 that is a 10% increase over the background level) for AA-induced mutations in the kidney of rats was about 7 μg/kg body weight per day. This review constitutes an overview of our investigations on AA-induced genotoxicity and toxicogenomic changes including gene expression, microRNA expression, and proteomics; and presents updated information focused on AA-induced genotoxicity in rodents.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2000

Genetic Network Inference: From Co-Expression Clustering to Reverse Engineering

Advances in molecular biological, analytical, and computational technologies are enabling us to systematically investigate the complex molecular processes underlying biological systems. In particular, using high-throughput gene expression assays, we are able to measure the output of the gene regulatory network. We aim here to review datamining and modeling approaches for conceptualizing and unraveling the functional relationships implicit in these datasets. Clustering of co-expression profiles allows us to infer shared regulatory inputs and functional pathways. We discuss various aspects of clustering, ranging from distance measures to clustering algorithms and multiple-duster memberships. More advanced analysis aims to infer causal connections between genes directly, i.e., who is regulating whom and how. We discuss several approaches to the problem of reverse engineering of genetic networks, from discrete Boolean networks, to continuous linear and non-linear models. We conclude that the combination of predictive modeling with systematic experimental verification will be required to gain a deeper insight into living organisms, therapeutic targeting, and bioengineering.

Dhaeseleer, Patrik↗
DOE OSTI2013

Transcription factors, sucrose, and sucrose metabolic genes interact to regulate potato phenylpropanoid metabolism

Much remains unknown about how transcription factors and sugars regulate phenylpropanoid metabolism in tuber crops like potato (Solanum tuberosum). Based on phylogeny and protein similarity to known regulators of phenylpropanoid metabolism, 15 transcription factors were selected and their expression was compared in white, yellow, red, and purple genotypes with contrasting phenolic and anthocyanin profiles. Red and purple genotypes had increased phenylalanine ammonia lyase (PAL) enzyme activity, markedly higher levels of phenylpropanoids, and elevated expression of most phenylpropanoid structural genes, including a novel anthocyanin O-methyltransferase. The transcription factors Anthocyanin1 (StAN1), basic Helix Loop Helix1 (StbHLH1), and StWD40 were more strongly expressed in red and purple potatoes. Expression of 12 other transcription factors was not associated with phenylpropanoid content, except for StMYB12B, which showed a negative relationship. Increased expression of AN1, bHLH1, and WD40 was also associated with environmentally mediated increases in tuber phenylpropanoids. Treatment of potato plantlets with sucrose induced hydroxycinnamic acids, flavonols, anthocyanins, structural genes, AN1, bHLH1, WD40, and genes encoding the sucrose-hydrolysing enzymes SUSY1, SUSY4, and INV2. Transient expression of StAN1 in tobacco leaves induced bHLH1, structural genes, SUSY1, SUSY4, and INV1, and increased phenylpropanoid amounts. StAN1 infiltration into tobacco leaves decreased sucrose and glucose concentrations. In silico promoter analysis revealed the presence of MYB and bHLH regulatory elements on sucrolytic gene promoters and sucrose-responsive elements on the AN1 promoter. These findings reveal an interesting dynamic between AN1, sucrose, and sucrose metabolic genes in modulating potato phenylpropanoids.

sucrose↗
DOE OSTI2013

An attenuated Lassa vaccine in SIV-infected rhesus macaques does not persist or cause arenavirus disease but does elicit Lassa virus-specific immunity

Background: Lassa hemorrhagic fever (LHF) is a rodent-borne viral disease that can be fatal for human beings. In this study, an attenuated Lassa vaccine candidate, ML29, was tested in SIV-infected rhesus macaques for its ability to elicit immune responses without instigating signs pathognomonic for arenavirus disease. ML29 is a reassortant between Lassa and Mopeia viruses that causes a transient infection in non-human primates and confers sterilizing protection from lethal Lassa viral challenge. However, since the LHF endemic area of West Africa also has high HIV seroprevalence, it is important to determine whether vaccination could be safe in the context of HIV infection. Results: SIV-infected and uninfected rhesus macaques were vaccinated with the ML29 virus and monitored for specific humoral and cellular immune responses, as well as for classical and non-classical signs of arenavirus disease. Classical disease signs included viremia, rash, respiratory distress, malaise, high liver enzyme levels, and virus invasion of the central nervous system. Non-classical signs, derived from profiling the blood transcriptome of virulent and non-virulent arenavirus infections, included increased expression of interferon-stimulated genes (ISG) and decreased expression of COX2, IL-1β, coagulation intermediates and nuclear receptors needed for stress signaling. All vaccinated monkeys showed ML29-specific antibody responses and ML29-specific cellmediated immunity. Conclusion: SIV-infected and uninfected rhesus macaques responded similarly to ML29 vaccination, and none developed chronic arenavirus infection. Importantly, none of the macaques developed signs, classical or nonclassical, of arenavirus disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Inter-strain heterogeneity in rat hepatic transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)

The biochemical and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) have been the subject of intense study for decades. It is now clear that essentially all TCDD-induced toxicities are mediated by DNA–protein interactions involving the Aryl Hydrocarbon Receptor (AHR). Nevertheless, it remains unknown which AHR target genes cause TCDD toxicities. Several groups, including our own, have developed rodent model systems to probe these questions. mRNA expression profiling of these model systems has revealed significant inter-species heterogeneity in rodent hepatic responses to TCDD. It has remained unclear if this variability also exists within a species, amongst rodent strains. To resolve this question, we profiled the hepatic transcriptomic response to TCDD of diverse rat strains (L-E, H/W, F344 and Wistar rats) and two lines derived from L-E × H/W crosses, at consistent age, sex, and dosing (100 μg/kg TCDD for 19 h). Using this uniquely consistent dataset, we show that the majority of TCDD-induced alterations in mRNA abundance are strain/line-specific: only 11 genes were affected by TCDD across all strains, including well-known dioxin-responsive genes such as Cyp1a1 and Nqo1. Our analysis identified two novel universally dioxin-responsive genes as well as 4 genes induced by TCDD in dioxin-sensitive rats only. These 6 genes are strong candidates to explain TCDD-related toxicities, so we validated them using 152 animals in time-course (0 to 384 h) and dose–response (0 to 3000 μg/kg) experiments. This study reveals that different rat strains exhibit dramatic transcriptional heterogeneity in their hepatic responses to TCDD and that inter-strain comparisons can help identify candidate toxicity-related genes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Proteome-wide Light/Dark Modulation of Thiol Oxidation in Cyanobacteria Revealed by Quantitative Site-Specific Redox Proteomics

Reversible protein thiol oxidation is an essential regulatory mechanism of photosynthesis, metabolism, and gene expression in photosynthetic organisms. Herein, we present the first proteome-wide quantitative site-specific profiling of in vivo dynamics of thiol oxidation modulated by light/dark in cyanobacterium Synechocystis sp. PCC 6803, an oxygenic photosynthetic prokaryote. Our proteomic approach integrates resin-assisted enrichment with isobaric labeling to enable site-specific and quantitative measurements of reversibly oxidized thiols. The redox dynamics of ~2,200 cysteine sites from 1,060 proteins under different conditions (light, dark, and in the presence of a photosystem II inhibitor DCMU) were quantified. The results revealed broad proteome-wide changes in thiol oxidation in many key biological processes including photosynthesis, carbon fixation, and glycolysis. Moreover, the redox sensitivity data provides a good prediction of potential functional cysteine sites for proteins of interest. Together, our findings provided significant novel insights into the broad redox regulation of photosynthetic organisms.

redox proteomics, photosynthesis, Synechocystis 68↗
NASA NTRS2022

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and Simulated Microgravity

The space environment includes unique hazards like radiation and microgravity which can adversely affect biological systems. We assessed a multi-omics NASA GeneLab dataset where mice were hindlimb unloaded and/or gamma irradiated for 21 days followed by retinal analysis at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles resulting in a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation difference was seen across different conditions at the same time point. Nucleotide metabolism biological processes were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. A total of 23 genes showed significant changes in methylation and expression compared to unexposed controls, including genes involved in retinal function and inflammatory response. This multi-omics analysis interrogates the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination and highlights important molecular mechanisms at different post-exposure stages.

Prachi Kothiyal↗
DOE OSTI2014

Investigating the roles of MicroRNAs in biotic stress response induced by Rhizoctonia solani in rice

Sheath blight disease, caused by Rhizoctonia solani 1802/KB was screened on two rice varieties, Oryza sativaindica cultivar MR219 and Oryza sativa indica cultivar UKMRC9. The disease symptom was severe in MR219 compared to UKMRC9. Total RNA from R. solani 1802/KB, infected rice leaves of MR219 and infected rice leaves of UKMRC9 were extracted using TRIzol reagent, purified and sent for small RNA sequencing. Three miRNA libraries were generated and analyzed. The libraries generated 65 805, 78 512 and 81 325 known miRNAs respectively. The structure of miRNA of these samples was predicted. The up-regulated and down-regulated of miRNAs target gene prediction and its target functions were discovered and were mainly related to the growth and development of metabolism, protein transport, transcriptional regulation, stress response, and hormone signaling and electron transfer. Sheath blight-induced differential expression of known miRNAs tends to targetMYB transcription factor, F-box proteins, NBS-LRR, leucine-rich repeat receptor protein kinases and zinc finger proteins. Detecting new miRNAs and measuring the expression profiles of known miRNAs is an important tasks required for a better understanding of various biological conditions. Therefore, further analysis using Gene Ontology Slim will be conducted to deduce some biological information from the datasets obtained.

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↗