Search NASA⌕ Search

SEARCH · Search NASA

Results for “comparative transcriptomics”

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 73 records · Page 4

Diterpenoid Alkaloid Biosynthesis in the Delphinium and Aconitum Genera

Diterpenoid alkaloids are known for their range of potential medicinal applications and are found primarily in the Delphinium and Aconitum genera. We sequenced RNA isolated from different tissue types of Delphinium grandiflorum, Aconitum plicatum, and Aconitum lycoctonum to identify six entry steps in their biosynthesis.

biosynthesis↗

Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae)

Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2 , a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C 3 and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.

54 ENVIRONMENTAL SCIENCES↗

divergence in basal transcriptomes of "Wild Type" Chlamydomonas reinhardtii strains

This is a study to compare the basal transcriptomes of several widely used laboratory strains of the Chlorophyte alga, Chlamydomonas reinhardtii. Given that there is a high degree of genetic diversity among the closely-related laboratory strains, we wished to examine how much variation there is at the transcriptome level. A panel of WT strains (CC-124, CC-125, CC-1009, CC-1690, CC-1691), all believed to be descended from a single zygospore isolated in 1945, were chosen based on their representing the oldest lineages among the standard laboratory strains. Additionally, CC-4532, which was the source for the current (v6) reference assembly, and CC-4533, which is the initial parental strain of the CLiP library collection of mutant strains, were also included in this study based on their significance to the Chlamydomonas community. All strains were grown in liquid cultures under identical, mixotrophic conditions (light + acetate) to mid-log phase before collecting mRNA for RNA-Seq analysis. Overall design: Pre-cultures of each strain (CC-124, CC-125, CC-1009, CC-1690, CC-1691, CC-4532, and CC-4533) were used to innoculate flasks of Tris-Acetate-Phosphate (TAP) media supplemented with Kropat's trace metals at a concentration of 1.5 x 10^4 cells/mL. Flasks were placed on a shaking platform at 180 RPM under 50-70 µmol x m^-2 x s^-1 of light. Cultures were grown until they reached ~2-3 x 10^6 cells/mL, about three days, before total RNA was collected and subjected RNA-Seq analysis.

Source record↗

A tissue‐resolved, network‐based transcriptomic framework for abiotic stress responses in sorghum

Developing climate‐resilient crops requires a detailed understanding of stress‐induced gene expression dynamics, as maladaptive responses can compromise their productivity and survival. Sorghum, a globally important cereal with exceptional tolerance to multiple abiotic stresses, provides a powerful system for investigating these dynamics. However, how stress type, tissue specificity, and temporal progression jointly shape transcriptomic responses in crops remains poorly understood. Here, we present a comparative, time‐resolved transcriptomic atlas of sorghum responses to drought, heat, and salinity stress across shoot and root tissues. Integrative analyses revealed that tissue specificity is the dominant determinant of abiotic stress‐induced gene reprogramming across all three stresses. Building on these global comparisons, we focused on heat stress, as it elicited the most coherent and pronounced transcriptional and regulatory responses, enabling deeper network‐level interrogation. Co‐expression network analysis identified tissue‐specific modules enriched for phytohormone‐responsive genes, while gene regulatory network (GRN) mapping and cistrome analyses uncovered transcription factors (TFs) controlling key hub genes within these modules. Together, this study provides a foundational transcriptomic and network‐based resource for dissecting the regulatory architecture of abiotic stress responses in sorghum and offers prioritized candidates for future functional validation and engineering of climate‐resilient crops.

abiotic stress↗

Polar algae flaunt their zinc assets

Metal ions have been at the centre of pivotal points in the evolution of extant life. Oxygen-evolving photosynthesis, which irrevocably reshaped the geochemistry and biology of our planet, has an absolute requirement for metal ions to harvest light, split water and transfer electrons. Access to abundant oxygen then led to the propagation of organisms capable of oxidative metabolism, a process that is also dependent on metal ions for electron transfer and reduction of molecular oxygen. Because of the essential links between metal ions and the proteins that require them to function, as well as changes in metal bioavailability through time, metals have shaped the trajectories that evolution can take. Iron typically steals this show, but with access to whole-genome sequences and transcriptomes, the imprint that zinc has made on biology is coming into focus. Writing in Nature Ecology & Evolution, Ye and colleagues present new insights into the complex connections between zinc bioavailability, adaptation of algae to the polar oceans and the evolution of regulatory networks. By combining field and laboratory-based analyses, the authors suggest that expanded families of zinc-binding proteins have enabled the green alga Microglena sp. YARC to flourish in the harsh conditions of the polar Southern Ocean (Fig. 1). As the waters where this alga occurs are known for their enrichment of zinc, the authors further propose that availability of this metal ion was directly responsible for successful microalgal colonization of polar oceans. To test this hypothesis, the authors compare meta-transcriptomes collected from pole-to-pole and find positive correlations between higher copy numbers for transcripts encoding putative zinc-binding proteins, higher latitudes, lower surface temperatures and dissolved zinc.

59 BASIC BIOLOGICAL SCIENCES↗

Cell Population–resolved Multiomics Atlas of the Developing Lung

The lung is a vital organ that undergoes extensive morphological and functional changes during postnatal development. To disambiguate how different cell populations contribute to organ development, we performed proteomic and transcriptomic analyses of four sorted cell populations from the lung of human subjects aged 0 to 8 years-old with a focus on early life. The cell populations analyzed included epithelial, endothelial, mesenchymal, and immune cells. Our results revealed distinct molecular signatures for each of the sorted cell populations that enable the description of molecular shifts occurring in these populations during post-natal development. Here, we confirmed that the proteome of the different cell populations was distinct regardless of age and identified functions specific to each population. We identified a series of cell population protein markers, including those located at the cell surface, that show differential expression and distribution on RNA in situ hybridization and immunofluorescence imaging. We validated the spatial distribution of AT1 and endothelial cell surface markers. Temporal analyses of the proteome of each of the four populations revealed processes modulated during postnatal development and disambiguating results obtained on whole tissue proteome. Finally, the proteome was compared to a transcriptomics survey performed on the same lung samples to evaluate processes under post-transcriptional control.

59 BASIC BIOLOGICAL SCIENCES↗

Microbial inoculum effects on the rumen epithelial transcriptome and rumen epimural metatranscriptome in calves

Manipulation of the rumen microbial ecosystem in early life may affect ruminal fermentation and enhance the productive performance of dairy cows. The objective of this experiment was to evaluate the effects of dosing three different types of microbial inoculum on the rumen epithelium tissue (RE) transcriptome and the rumen epimural metatranscriptome (REM) in dairy calves. For this objective, 15 Holstein bull calves were enrolled in the study at birth and assigned to three different intraruminal inoculum treatments dosed orally once weekly from three to six weeks of age. The inoculum treatments were prepared from rumen contents collected from rumen fistulated lactating cows and were either autoclaved (control; ARF), processed by differential centrifugation to create the bacterial-enriched inoculum (BE), or through gravimetric separation to create the protozoal-enriched inoculum (PE). Calves were fed 2.5 L/d pasteurized waste milk 3x/d from 0 to 7 weeks of age and texturized starter until euthanasia at 9 weeks of age, when the RE tissues were collected for transcriptome and microbial metatranscriptome analyses, from four randomly selected calves from each treatment. The different types of inoculum altered the RE transcriptome and REM. Compared to ARF, 9 genes were upregulated in the RE of BE and 92 in PE, whereas between BE and PE there were 13 genes upregulated in BE and 114 in PE. Gene ontology analysis identified enriched GO terms in biological process category between PE and ARF, with no enrichment between BE and ARF. The RE functional signature showed different KEGG pathways related to BE and ARF, and no specific KEGG pathway for PE. We observed a lower alpha diversity index for RE microbiome in ARF (observed genera and Chao1 (p < 0.05)). Five microbial genera showed a significant correlation with the changes in host gene expression: Roseburia (25 genes), Entamoeba (two genes); Anaerosinus, Lachnospira, and Succiniclasticum were each related to one gene. sPLS-DA analysis showed that RE microbial communities differ among the treatments, although the taxonomic and functional microbial profiles show different distributions. Co-expression Differential Network Analysis indicated that both BE and PE had an impact on the abundance of KEGG modules related to acyl-CoA synthesis, type VI secretion, and methanogenesis, while PE had a significant impact on KEGGs related to ectoine biosynthesis and D-xylose transport. Our study indicated that artificial dosing with different microbial inocula in early life alters not only the RE transcriptome, but also affects the REM and its functions.

59 BASIC BIOLOGICAL SCIENCES↗

Multi‐Omics Analyses Reveal Divergent Molecular Mechanisms Underlying Plant Biomass Conversion by Five Fungi

Fungal plant biomass conversion (FPBC) is of great importance to the global carbon cycle and has been increasingly applied for the production of biofuel and biochemicals from lignocellulose. However, the comprehensive understanding of relevant molecular mechanisms in different fungi remains challenging. Here, we comparatively analyzed the transcriptome, proteome and metabolome profile of four ascomycetes and one basidiomycete fungi during their growth on two common agricultural feedstocks (soybean hulls and corn stover). We revealed strong time‐, substrate‐ and species‐specific responses at multi‐omics levels for the tested fungi, highlighting species‐specific carbon utilization approaches and evolutionary adaptation to environmental niches. Notably, a remarkable expressional diversity of lignocellulose degrading enzymes, sugar transporter and metabolic genes, as well as industrially relevant metabolites were identified across different fungi and cultivation conditions. The findings improves our understanding of complex molecular networks underlying FPBC and fungal ecological roles, offering novel insights that can guide future genetic engineering of fungi for valorization of agriculture waste into value‐added bioproducts.

CAZy↗

Convergent evolution of aerobic fermentation through divergent mechanisms acting on key shared glycolytic genes

As the tree of life becomes increasingly accessible to molecular investigations, describing mechanisms underlying evolutionary convergence and constraint will be crucial to understanding diversification. The lineage including the model yeast Saccharomyces cerevisiae evolved aerobic fermentation in part through an ancient whole genome duplication and retention of glycolytic genes. To evaluate glycolytic rates across diverse yeasts, we developed and deployed an extracellular acidification rates (ECAR) assay on 299 species that span more than 400 million years of evolution and identified a clade in the genus Saturnispora that convergently evolved aerobic fermentation. Through comparative genomics and transcriptomics, we found that several glycolytic genes had higher expression and novel cis-regulatory elements in aerobically fermenting Saturnispora species. When the transcription factor required for their activation was deleted in Saturnispora dispora, the mutants had reduced glycolytic rates and increased respiration. Intriguingly, many of the upregulated genes are orthologous to duplicated glycolytic genes in S. cerevisiae. These divergent genetic mechanisms affecting the same set of genes suggest that there are strong evolutionary constraints on how aerobic fermentation can arise.

Horianopoulos, Linda C. [Great Lakes Bioenergy Res↗

Gene and genome duplications have contrasting impacts on biosynthetic and flower developmental pathways in California poppy

Benzylisoquinoline alkaloids (BIAs) represent a vast group of specialized plant metabolites with diverse pharmaceutical applications, synthesized by a variety of gene families. Among the multiple plant lineages that produce BIAs, the most notable is the poppy family (Papaveraceae), with California poppy (Eschscholzia californica) emerging as a model organism. Here, we report a haplotype-resolved genome assembly, in combination with a high-density expression atlas, for California poppy. Genome analyses reveal recent diversification of BIA biosynthesis genes in poppy through localized duplications. Furthermore, we demonstrate that the degree of phylogenetic relatedness among paralogs within BIA biosynthesis-associated gene families correlates with similarities in gene expression. In contrast, gene families involved in carotenoid biosynthesis, which contributes to the intense orange petal pigmentation, are not phylogenetically clustered, and floral developmental regulators exhibit a high degree of retention of gene duplicates associated with ancient polyploidy events. These findings illustrate alternative roles for gene and genome duplications as drivers of trait evolution. Given the position of California poppy in the angiosperm phylogeny, the high-quality genomic resources generated for this work constitute a valuable resource for comparative genomic and transcriptomic analyses for poppies and flowering plants more generally.

Rössner, Le-Han [Justus-Liebig University, Giessen↗

The biosynthesis of the anti–microbial diterpenoid leubethanol in Leucophyllum frutescens proceeds via an all– cis prenyl intermediate

Serrulatane diterpenoids are natural products found in plants from a subset of genera within the figwort family (Scrophulariaceae). Many of these compounds have been characterized as having anti–microbial properties and share a common diterpene backbone. One example, leubethanol from Texas sage (Leucophyllum frutescens) has demonstrated activity against multi–drug–resistant tuberculosis. Leubethanol is the only serrulatane diterpenoid identified from this genus; however, a range of such compounds have been found throughout the closely related Eremophila genus. Despite their potential therapeutic relevance, the biosynthesis of serrulatane diterpenoids has not been previously reported. Here we leverage the simple product profile and high accumulation of leubethanol in the roots of L. frutescens and compare tissue–specific transcriptomes with existing data from Eremophila serrulata to decipher the biosynthesis of leubethanol. A short–chain cis–prenyl transferase (LfCPT1) first produces the rare diterpene precursor nerylneryl diphosphate, which is cyclized by an unusual plastidial terpene synthase (LfTPS1) into the characteristic serrulatane diterpene backbone. Final conversion to leubethanol is catalyzed by a cytochrome P450 (CYP71D616) of the CYP71 clan. This pathway documents the presence of a short–chain cis–prenyl diphosphate synthase, previously only found in Solanaceae, which is likely involved in the biosynthesis of other known diterpene backbones in Eremophila. LfTPS1 represents neofunctionalization of a compartment–switching terpene synthase accepting a novel substrate in the plastid. Biosynthetic access to leubethanol will enable pathway discovery to more complex serrulatane diterpenoids which share this common starting structure and provide a platform for the production and diversification of this class of promising anti–microbial therapeutics in heterologous systems.

59 BASIC BIOLOGICAL SCIENCES↗

Combined analysis of microbial metagenomic and metatranscriptomic sequencing data to assess in situ physiological conditions in the premature infant gut

Microbes alter their transcriptomic profiles in response to the environment. The physiological conditions experienced by a microbial community can thus be inferred using meta-transcriptomic sequencing by comparing transcription levels of specifically chosen genes. However, this analysis requires accurate reference genomes to identify the specific genes from which RNA reads originate. In addition, such an analysis should avoid biases in transcript counts related to differences in organism abundance. In this study we describe an approach to address these difficulties. Sample-specific meta-genomic assembled genomes (MAGs) were used as reference genomes to accurately identify the origin of RNA reads, and transcript ratios of genes with opposite transcription responses were compared to eliminate biases related to differences in organismal abundance, an approach hereafter named the “diametric ratio” method. We used this approach to probe the environmental conditions experienced by Escherichia spp. in the gut of 4 premature infants, 2 of whom developed necrotizing enterocolitis (NEC), a severe inflammatory intestinal disease. We analyzed twenty fecal samples taken from four premature infants (4–6 time points from each infant), and found significantly higher diametric ratios of genes associated with low oxygen levels in samples of infants later diagnosed with NEC than in samples without NEC. We also show this method can be used for examining other physiological conditions, such as exposure to nitric oxide and osmotic pressure. These study results should be treated with caution, due to the presence of confounding factors that might also distinguish between NEC and control infants. Nevertheless, together with benchmarking analyses, we show here that the diametric ratio approach can be applied for evaluating the physiological conditions experienced by microbes in situ. Results from similar studies can be further applied for designing diagnostic methods to detect NEC in its early developmental stages.

59 BASIC BIOLOGICAL SCIENCES↗

Oligosaccharide production and signaling correlate with delayed flowering in an Arabidopsis genotype grown and selected in high [CO 2 ]

Since industrialization began, atmospheric CO 2 ([CO 2 ]) has increased from 270 to 415 ppm and is projected to reach 800–1000 ppm this century. Some Arabidopsis thaliana (Arabidopsis) genotypes delayed flowering in elevated [CO 2 ] relative to current [CO 2 ], while others showed no change or accelerations. To predict genotype-specific flowering behaviors, we must understand the mechanisms driving flowering response to rising [CO 2 ]. [CO 2 ] changes alter photosynthesis and carbohydrates in plants. Plants sense carbohydrate levels, and exogenous carbohydrate application influences flowering time and flowering transcript levels. We asked how organismal changes in carbohydrates and transcription correlate with changes in flowering time under elevated [CO 2 ]. We used a genotype (SG) of Arabidopsis that was selected for high fitness at elevated [CO 2 ] (700 ppm). SG delays flowering under elevated [CO 2 ] (700 ppm) relative to current [CO 2 ] (400 ppm). We compared SG to a closely related control genotype (CG) that shows no [CO 2 ]-induced flowering change. We compared metabolomic and transcriptomic profiles in these genotypes at current and elevated [CO 2 ] to assess correlations with flowering in these conditions. While both genotypes altered carbohydrates in response to elevated [CO 2 ], SG had higher levels of sucrose than CG and showed a stronger increase in glucose and fructose in elevated [CO 2 ]. Both genotypes demonstrated transcriptional changes, with CG increasing genes related to fructose 1,6-bisphosphate breakdown, amino acid synthesis, and secondary metabolites; and SG decreasing genes related to starch and sugar metabolism, but increasing genes involved in oligosaccharide production and sugar modifications. Genes associated with flowering regulation within the photoperiod, vernalization, and meristem identity pathways were altered in these genotypes. Elevated [CO 2 ] may alter carbohydrates to influence transcription in both genotypes and delayed flowering in SG. Changes in the oligosaccharide pool may contribute to delayed flowering in SG. This work extends the literature exploring genotypic-specific flowering responses to elevated [CO 2 ].

59 BASIC BIOLOGICAL SCIENCES↗

EvoNet: A phylogenomic and systems biology approach to identify genes underlying plant survival in marginal, low‐N soils

The DOE‐BER “EvoNet” project investigates the genetic and molecular basis of plant resilience in extreme environments. We do this by identifying key genes that enable “extreme survivor” species to thrive in the nitrogen-poor soils of Chile’s hyper-arid Atacama Desert. Our collections focus on 32 Atacama extremophile species, including seven grass species with potential biofuel applications. To identify genes-of-importance to survival we compared genomic and transcriptomic profiles of extremophile species that thrive in the Atacama to those of closely related “sister” species from nitrogen-rich arid and mesic regions of California. Deep RNA sequencing and de novo transcriptome assembly across these triplet species sets supported a phylogenomic framework for identifying positively selected genes associated with adaptive divergence. Our integrative analysis combined ecological and environmental data, metagenomics, evolutionary and systems biology, and metabolomics. This enabled us to create an unprecedented framework for systematically understanding how non-model plants have adapted to survive in extreme conditions. Our resulting database of positively selected ortholog groups in the extremophile plants offers promising targets for engineering crop and biofuel species with enhanced resilience to drought and extreme weather. Additionally, our newest dataset explores and exploits a complementary metabolomic approach. This new aspect provides innovative strategies to manipulate plant cell metabolism, further supporting efforts to improve agricultural productivity in the face of extreme climates. Importantly, our combined evolutionary- and metabolomic-based strategies focused on convergent patterns of adaptation, providing a genetic and metabolomic toolkit for improving crop and biofuel resilience across diverse plant species. Finally, our novel exploration of ecological and evolutionary dynamics delivered to the community a phylogenomic computational pipeline called “PhyloGeneious.” Our continued adaptations of this pipeline are publicly available to expedite evolutionary genomic research for future scientific discoveries. In total, our DOE-BER has provided genomic, metabolomic, and computational strategies to understand how extremophile plants provide evolutionary and physiological targets for improving agricultural and biofuel production.

59 BASIC BIOLOGICAL SCIENCES↗

A consensus-based ensemble approach to improve transcriptome assembly

Systems-level analyses, such as differential gene expression analysis, co-expression analysis, and metabolic pathway reconstruction, depend on the accuracy of the transcriptome. Multiple tools exist to perform transcriptome assembly from RNAseq data. However, assembling high quality transcriptomes is still not a trivial problem. This is especially the case for non-model organisms where adequate reference genomes are often not available. Different methods produce different transcriptome models and there is no easy way to determine which are more accurate. Furthermore, having alternative-splicing events exacerbates such difficult assembly problems. While benchmarking transcriptome assemblies is critical, this is also not trivial due to the general lack of true reference transcriptomes. In this study, we first provide a pipeline to generate a set of the simulated benchmark transcriptome and corresponding RNAseq data. Using the simulated benchmarking datasets, we compared the performance of various transcriptome assembly approaches including both de novo and genome-guided methods. The results showed that the assembly performance deteriorates significantly when alternative transcripts (isoforms) exist or for genome-guided methods when the reference is not available from the same genome. To improve the transcriptome assembly performance, leveraging the overlapping predictions between different assemblies, we present a new consensus-based ensemble transcriptome assembly approach, ConSemble. Without using a reference genome, ConSemble using four de novo assemblers achieved an accuracy up to twice as high as any de novo assemblers we compared. When a reference genome is available, ConSemble using four genome-guided assemblies removed many incorrectly assembled contigs with minimal impact on correctly assembled contigs, achieving higher precision and accuracy than individual genome-guided methods. Furthermore, ConSemble using de novo assemblers matched or exceeded the best performing genome-guided assemblers even when the transcriptomes included isoforms. We thus demonstrated that the ConSemble consensus strategy both for de novo and genome-guided assemblers can improve transcriptome assembly. The RNAseq simulation pipeline, the benchmark transcriptome datasets, and the script to perform the ConSemble assembly are all freely available from: http://bioinfolab.unl.edu/emlab/consemble/.

59 BASIC BIOLOGICAL SCIENCES↗

MOSAIC: a joint modeling methodology for combined circadian and non-circadian analysis of multi-omics data

Abstract Motivation Circadian rhythms are approximately 24-h endogenous cycles that control many biological functions. To identify these rhythms, biological samples are taken over circadian time and analyzed using a single omics type, such as transcriptomics or proteomics. By comparing data from these single omics approaches, it has been shown that transcriptional rhythms are not necessarily conserved at the protein level, implying extensive circadian post-transcriptional regulation. However, as proteomics methods are known to be noisier than transcriptomic methods, this suggests that previously identified arrhythmic proteins with rhythmic transcripts could have been missed due to noise and may not be due to post-transcriptional regulation. Results To determine if one can use information from less-noisy transcriptomic data to inform rhythms in more-noisy proteomic data, and thus more accurately identify rhythms in the proteome, we have created the Multi-Omics Selection with Amplitude Independent Criteria (MOSAIC) application. MOSAIC combines model selection and joint modeling of multiple omics types to recover significant circadian and non-circadian trends. Using both synthetic data and proteomic data from Neurospora crassa, we showed that MOSAIC accurately recovers circadian rhythms at higher rates in not only the proteome but the transcriptome as well, outperforming existing methods for rhythm identification. In addition, by quantifying non-circadian trends in addition to circadian trends in data, our methodology allowed for the recognition of the diversity of circadian regulation as compared to non-circadian regulation. Availability and implementation MOSAIC’s full interface is available at https://github.com/delosh653/MOSAIC. An R package for this functionality, mosaic.find, can be downloaded at https://CRAN.R-project.org/package=mosaic.find. Supplementary information Supplementary data are available at Bioinformatics online.

De los Santos, Hannah↗

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↗