Search NASA⌕ Search

SEARCH · Search NASA

Results for “gene expression analysis”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sorghum bicolor BTx623 Nitrogen Grown Conditions Set2 Gene Expression Profiling

Dhurrin, a cyanogenic glucoside, plays an important role in Sorghum bicolor physiology and defense. The concentration of dhurrin in sorghum is influenced by both nitrogen status and stage of plant organ development. While nitrogen resupply activates the expression of genes for dhurrin biosynthesis, the molecular mechanisms underlying this regulation remain unclear. In this study, we investigated the transcriptional response of sorghum to nitrogen resupply following growth under nitrogen-limiting conditions. Using a time-course design, we measured hydrogen cyanide potential (HCNp), growth, and nitrate content at 0-, 2-, 6-, 12-, 24-, 36-, 48-, and 60-h after resupply and collected tissue for RNAseq analysis in parallel for analysis of gene expression and construction of gene regulatory networks (GRNs). HCNp (mg g−1 DW) increased significantly in leaf and stem tissues following nitrogen resupply, with increases in the leaf partially driven by continued declines in controls under ongoing nitrogen stress. Expression of the dhurrin pathway genes was upregulated in leaves from 24 h after nitrogen resupply, with diel expression patterns observable over the remaining time points. No upregulation was observed in roots or stems, suggesting that developmental context overrides environmental cues. GRN analysis identified candidate transcription factors regulating dhurrin biosynthesis genes, including members of the MYB, bZIP, and GARP-type transcription factor families. Some of these candidate transcription factors may be involved in relieving senescence-associated suppression of dhurrin biosynthesis and link nitrogen signaling to pathway activation. These findings provide new insight into the nitrogen-responsive regulation of dhurrin in sorghum, highlighting candidate regulators for future functional characterization.

cyanogenic glucoside↗
DOE OSTI2020

Perimatrix of middle ear cholesteatoma: A granulation tissue with a specific transcriptomic signature

Objectives/Hypothesis: To establish comprehensive transcriptomic profiles of cholesteatoma perimatrix tissue and granulation tissue from chronic otitis media (COM) that did not develop cholesteatoma, which can indicate molecular pathways involved in the cholesteatoma perimatrix pathology and invasiveness. Study Design: Retrospective Case Series. Methods: Transcriptome data were obtained from cholesteatoma perimatrix tissue and COM granulation tissue by an Illumina iScan microarray. Differentially expressed genes (DEGs) were subsequently analyzed using both bioinformatical functional annotation and network analysis. Expression of candidate genes (MMP9 and LCN2) was validated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) on a larger group of samples. Results: Analysis of the transcriptome led to the identification of 169 differentially expressed genes between investigated tissues. Bioinformatic analysis suggested that most significant biological processes involving DEGs were previously described in cholesteatoma pathology. Network analysis identified ERBB2, TFAP2A, and TP63 as major hubs of the DEGs molecular network. Furthermore, it was observed that the cellular component most significantly enriched in DEGs was extracellular space containing 47 DEGs. Using qRT-PCR, it was confirmed that mRNA levels of the major extracellular hub (MMP9) are increased, whereas its interacting molecule (LCN2) mRNA levels were decreased in cholesteatoma perimatrix tissue compared to COM granulation tissue. Conclusions: The current study approach offers an overall look at molecular mechanisms that describe the cholesteatoma entity by focusing exclusively on the perimatrix processes in comparison to COM granulation tissue. The observed differences in gene expression between cholesteatoma perimatrix and COM granulation tissue could suggest novel markers potentially influenced by the perimatrix–matrix molecular interplay, which is not present in COM without cholesteatoma. Level of Evidence: NA. Laryngoscope, 130:E220–E227, 2020. © 2019 The American Laryngological, Rhinological and Otological Society, Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Transcriptomic and Proteomic Insights into Host Immune Responses in Pediatric Severe Malarial Anemia: Dysregulation in HSP60-70-TLR2/4 Signaling and Altered Glutamine Metabolism

Severe malarial anemia (SMA, Hb < 6.0 g/dL) is a leading cause of childhood morbidity and mortality in holoendemic Plasmodium falciparum transmission zones. This study explored the entire expressed human transcriptome in whole blood from 66 Kenyan children with non-SMA (Hb ≥ 6.0 g/dL, n = 41) and SMA (n = 25), focusing on host immune response networks. RNA-seq analysis revealed 6862 differentially expressed genes, with equally distributed up-and down-regulated genes, indicating a complex host immune response. Deconvolution analyses uncovered leukocytic immune profiles indicative of a diminished antigenic response, reduced immune priming, and polarization toward cellular repair in SMA. Weighted gene co-expression network analysis revealed that immune-regulated processes are central molecular distinctions between non-SMA and SMA. A top dysregulated immune response signaling network in SMA was the HSP60-HSP70-TLR2/4 signaling pathway, indicating altered pathogen recognition, innate immune activation, stress responses, and antigen recognition. Validation with high-throughput gene expression from a separate cohort of Kenyan children (n = 50) with varying severities of malarial anemia (n = 38 non-SMA and n = 12 SMA) confirmed the RNA-seq findings. Proteomic analyses in 35 children with matched transcript and protein abundance (n = 19 non-SMA and n = 16 SMA) confirmed dysregulation in the HSP60-HSP70-TLR2/4 signaling pathway. Additionally, glutamine transporter and glutamine synthetase genes were differentially expressed, indicating altered glutamine metabolism in SMA. This comprehensive analysis underscores complex immune dysregulation and novel pathogenic features in SMA.

Microbiology↗
DOE OSTI2023

Comparative transcriptome analysis of rainbow trout gonadal cells (RTG-2) infected with U and J genogroup infectious hematopoietic necrosis virus

Infectious hematopoietic necrosis virus (IHNV) is the causative pathogen of infectious hematopoietic necrosis, outbreaks of which are responsible for significant losses in rainbow trout aquaculture. Strains of IHNV isolated worldwide have been classified into five major genogroups, J, E, L, M, and U. To date, comparative transcriptomic analysis has only been conducted individually for the J and M genogroups. In this study, we compared the transcriptome profiles in U genogroup and J genogroup IHNV-infected RTG-2 cells with mock-infected RTG-2 cells. The RNA-seq results revealed 17,064 new genes, of which 7,390 genes were functionally annotated. Differentially expressed gene (DEG) analysis between U and J IHNV-infected cells revealed 2,238 DEGs, including 1,011 downregulated genes and 1,227 upregulated genes. Among the 2,238 DEGs, 345 new genes were discovered. The DEGs related to immune responses, cellular signal transduction, and viral diseases were further analyzed. RT-qPCR validation confirmed that the changes in expression of the immune response-related genes trpm2 , sting , itgb7 , ripk2 , and irf1 , cellular signal transduction-related genes irl , cacnb2 , bmp2l , gadd45α , and plk2 , and viral disease-related genes mlf1, mtor, armc5 , pik3r1 , and c-myc were consistent with the results of transcriptome analysis. Taken together, our findings provide a comprehensive transcriptional analysis of the differential virulence of the U and J genogroups of IHNV, and shed new light on the pathogenic mechanisms of IHNV strains.

Zhao, Jing-Zhuang↗
DOE OSTI2020

Insights of Molecular Mechanism of Xylem Development in Five Black Poplar Cultivars

Black poplar (Populus deltoides, P. nigra, and their hybrids) is the main poplar cultivars in China. It offers interesting options of large-scale biomass production for bioenergy due to its rapid growth and high yield. Poplar wood properties were associated with chemical components and physical structures during wood formation. In this study, five poplar cultivars, P. euramericana ‘Zhonglin46’ (Pe1), P. euramericana ‘Guariento’ (Pe2), P. nigra ‘N179’ (Pn1), P. deltoides ‘Danhong’ (Pd1), and P. deltoides ‘Nanyang’ (Pd2), were used to explore the molecular mechanism of xylem development. We analyzed the structural differences of developing xylem in the five cultivars and profiled the transcriptome-wide gene expression patterns through RNA sequencing. The cross sections of the developing xylem showed that the cell wall thickness of developed fiber in Pd1 was thickest and the number of xylem vessels of Pn1 was the least. A total of 10,331 differentially expressed genes were identified among 10 pairwise comparisons of the five cultivars, most of them were related to programmed cell death and secondary cell wall thickening. K-means cluster analysis and Gene Ontology enrichment analysis showed that the genes highly expressed in Pd1 were related to nucleotide decomposition, metabolic process, transferase, and microtubule cytoskeleton; whereas the genes highly expressed in Pn1 were involved in cell wall macromolecule decomposition and polysaccharide binding processes. Based on a weighted gene co-expression network analysis, a large number of candidate regulators for xylem development were identified. And their potential regulatory roles to cell wall biosynthesis genes were validated by a transient overexpression system. This study provides a set of promising candidate regulators for genetic engineering to improve feedstock and enhance biofuel conversion in the bioenergy crop Populus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Transcript and metabolite network perturbations in lignin biosynthetic mutants of Arabidopsis

Abstract Lignin, one of the most abundant polymers in plants, is derived from the phenylpropanoid pathway, which also gives rise to an array of metabolites that are essential for plant fitness. Genetic engineering of lignification can cause drastic changes in transcription and metabolite accumulation with or without an accompanying development phenotype. To understand the impact of lignin perturbation, we analyzed transcriptome and metabolite data from the rapidly lignifying stem tissue in 13 selected phenylpropanoid mutants and wild-type Arabidopsis (Arabidopsis thaliana). Our dataset contains 20,974 expressed genes, of which over 26% had altered transcript levels in at least one mutant, and 18 targeted metabolites, all of which displayed altered accumulation in at least one mutant. We found that lignin biosynthesis and phenylalanine supply via the shikimate pathway are tightly co-regulated at the transcriptional level. The hierarchical clustering analysis of differentially expressed genes (DEGs) grouped the 13 mutants into 5 subgroups with similar profiles of mis-regulated genes. Functional analysis of the DEGs in these mutants and correlation between gene expression and metabolite accumulation revealed system-wide effects on transcripts involved in multiple biological processes.

Plant Sciences↗
DOE OSTI2023

Single-cell quantification of ribosome occupancy in early mouse development

Translation regulation is critical for early mammalian embryonic development. However, previous studies had been restricted to bulk measurements, precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N 6 -methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Decomprolute is a benchmarking platform designed for multiomics-based tumor deconvolution

Tumor deconvolution is a reliable way to disentangle the diverse cell types that comprise solid tumors. To date, however, both the algorithms developed to deconvolve tumor samples, and the gold standard datasets used to assess the algorithms are geared toward the analysis of gene expression (e.g., RNA-seq) rather than protein levels in tumor cells. While gene expression is less expensive to measure, protein levels provide a more accurate view of immune markers. To facilitate the development as well as improve the reproducibility and reusability of multi-omic deconvolution algorithms, we introduce Decomprolute, a Common Workflow Language framework that leverages containerization to compare tumor deconvolution algorithms across multiomic data sets. Decomprolute incorporates the large-scale multiomic data sets produced by the Clinical Proteomic Tumor Analysis Consortium (CPTAC), which include matched mRNA expression and proteomic data from thousands of tumors across multiple cancer types to build a fully open-source, containerized proteogenomic tumor deconvolution benchmarking platform. The platform consists of modular architecture and it comes with well-defined input and output formats at each module. As a result, it is robust and extendable easily with additional algorithms or analyses. The platform is available for access and use at http://pnnl-compbio.github.io/decomprolute.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Phylogeny and Expression Atlas of the NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY in Agave

Agave species are widely cultivated crassulacean acid metabolism (CAM) plants for alcoholic beverages, food and fiber production. Among these, the Agave hybrid H11648 ((A. amaniensis × A. angustifolia) × A. amaniensis) is the main cultivar for sisal fiber in the tropical areas of Brazil, China, and African countries. The plants of Agave hybrid H11648 have a long life cycle and large leaves, which require a huge amount of nitrogen nutrient. However, the molecular basis of nitrogen transport and allocation has not been well understood in agave. In this study, we identified 19 NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY(NPF) genes (called AhNPFs) with full-length coding sequences in Agave hybrid H11648. Our analysis of gene expression in various types of tissues revealed the tissue-specific expression pattern of AhNPFs. We further examined their expression patterns at different leaf developmental stages, under abiotic/biotic stresses and nutrient deficiency. The results reveal several candidate regulators in the agave NPF family, including AhNPF4.3/5.2/7.1. We first characterized the NPF genes in agave based on published leaf transcriptome datasets and emphasized their potential functions. The study will benefit future studies related to nitrogen nutrient in agave.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

De novo transcriptome in roots of switchgrass ( Panicum virgatum L. ) reveals gene expression dynamic and act network under alkaline salt stress

Background: Soil salinization is a major limiting factor for crop cultivation. Switchgrass is a perennial rhizomatous bunchgrass that is considered an ideal plant for marginal lands, including sites with saline soil. Here we investigated the physiological responses and transcriptome changes in the roots of Alamo (alkaline-tolerant genotype) and AM314/MS-155 (alkaline-sensitive genotype) under alkaline salt stress. Results: Alkaline salt stress significantly affected the membrane, osmotic adjustment and antioxidant systems in switchgrass roots, and the ASTTI values between Alamo and AM-314/MS-155 were divergent at different time points. A total of 108,319 unigenes were obtained after reassembly, including 73,636 unigenes in AM-314/MS-155 and 65,492 unigenes in Alamo. A total of 10,219 DEGs were identified, and the number of upregulated genes in Alamo was much greater than that in AM-314/MS-155 in both the early and late stages of alkaline salt stress. The DEGs in AM-314/MS-155 were mainly concentrated in the early stage, while Alamo showed greater advantages in the late stage. These DEGs were mainly enriched in plant-pathogen interactions, ubiquitin-mediated proteolysis and glycolysis/gluconeogenesis pathways. We characterized 1480 TF genes into 64 TF families, and the most abundant TF family was the C2H2 family, followed by the bZIP and bHLH families. A total of 1718 PKs were predicted, including CaMK, CDPK, MAPK and RLK. WGCNA revealed that the DEGs in the blue, brown, dark magenta and light steel blue 1 modules were associated with the physiological changes in roots of switchgrass under alkaline salt stress. The consistency between the qRT-PCR and RNA-Seq results confirmed the reliability of the RNA-seq sequencing data. A molecular regulatory network of the switchgrass response to alkaline salt stress was preliminarily constructed on the basis of transcriptional regulation and functional genes. Conclusions: Alkaline salt tolerance of switchgrass may be achieved by the regulation of ion homeostasis, transport proteins, detoxification, heat shock proteins, dehydration and sugar metabolism. These findings provide a comprehensive analysis of gene expression dynamic and act network induced by alkaline salt stress in two switchgrass genotypes and contribute to the understanding of the alkaline salt tolerance mechanism of switchgrass and the improvement of switchgrass germplasm.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Mercury methylation by metabolically versatile and cosmopolitan marine bacteria

Microbes transform aqueous mercury (Hg) into methylmercury (MeHg), a potent neurotoxin that accumulates in terrestrial and marine food webs, with potential impacts on human health. This process requires the gene pair hgcAB, which encodes for proteins that actuate Hg methylation, and has been well described for anoxic environments. However, recent studies report potential MeHg formation in suboxic seawater, although the microorganisms involved remain poorly understood. In this study, we conducted large-scale multi-omic analyses to search for putative microbial Hg methylators along defined redox gradients in Saanich Inlet, British Columbia, a model natural ecosystem with previously measured Hg and MeHg concentration profiles. Analysis of gene expression profiles along the redoxcline identified several putative Hg methylating microbial groups, including Calditrichaeota, SAR324 and Marinimicrobia, with the last the most active based on hgc transcription levels. Marinimicrobia hgc genes were identified from multiple publicly available marine metagenomes, consistent with a potential key role in marine Hg methylation. Computational homology modelling predicts that Marinimicrobia HgcAB proteins contain the highly conserved amino acid sites and folding structures required for functional Hg methylation. Furthermore, a number of terminal oxidases from aerobic respiratory chains were associated with several putative novel Hg methylators. Our findings thus reveal potential novel marine Hg-methylating microorganisms with a greater oxygen tolerance and broader habitat range than previously recognized.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Discovery of widespread transcription initiation at microsatellites predictable by sequence-based deep neural network

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A clinically and genomically annotated nerve sheath tumor biospecimen repository

Nerve sheath tumors occur as a heterogeneous group of neoplasms in patients with neurofibromatosis type 1 (NF1). The malignant form represents the most common cause of death in people with NF1, and even when benign, these tumors can result in significant disfigurement, neurologic dysfunction, and a range of profound symptoms. Lack of human tissue across the peripheral nerve tumors common in NF1 has been a major limitation in the development of new therapies. To address this unmet need, we have created an annotated collection of patient tumor samples, patient-derived cell lines, and patient-derived xenografts, and carried out high-throughput genomic and transcriptomic characterization to serve as a resource for further biologic and preclinical therapeutic studies. In this work, we release genomic and transcriptomic datasets comprised of 55 tumor samples derived from 23 individuals, complete with clinical annotation. All data are publicly available through the NF Data Portal and at http://synapse.org/jhubiobank.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

A mechanistic model and therapeutic interventions for COVID-19 involving a RAS-mediated bradykinin storm

Neither the disease mechanism nor treatments for COVID-19 are currently known. Here, we present a novel molecular mechanism for COVID-19 that provides therapeutic intervention points that can be addressed with existing FDA-approved pharmaceuticals. The entry point for the virus is ACE2, which is a component of the counteracting hypotensive axis of RAS. Bradykinin is a potent part of the vasopressor system that induces hypotension and vasodilation and is degraded by ACE and enhanced by the angiotensin 1-9 produced by ACE2. Here, we perform a new analysis on gene expression data from cells in bronchoalveolar lavage fluid (BALF) from COVID-19 patients that were used to sequence the virus. Comparison with BALF from controls identifies a critical imbalance in RAS represented by decreased expression of ACE in combination with increases in ACE2, renin, angiotensin, key RAS receptors, kinogen and many kallikrein enzymes that activate it, and both bradykinin receptors. This very atypical pattern of the RAS is predicted to elevate bradykinin levels in multiple tissues and systems that will likely cause increases in vascular dilation, vascular permeability and hypotension. These bradykinin-driven outcomes explain many of the symptoms being observed in COVID-19.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Deletion analysis of the itaconic acid biosynthesis gene cluster components in Aspergillus pseudoterreus ATCC32359

The filamentous fungus Aspergillus terreus has been successfully used for industrial production of itaconic acid (IA) for many years. The IA biosynthesis pathway has recently been characterized at a molecular genetic level as an IA gene cluster (IA cluster) by a clone-based transcriptomic approach. The cluster consists of four genes, including genes for cis-aconitic acid decarboxylase (cadA), a predicted transcription factor (tf), mitochondrial organic acid transporter (mttA), and MFS (Major Facilitator Superfamily) type transporter (mfsA). In this research, we performed EST (Expressed Sequence Tag) analysis and systematic gene deletions to investigate the role of those genes in IA biosynthesis in A. pseudoterreus ATCC32359. EST analysis showed a similar expression pattern among those four genes that was distinct from neighbouring genes and further confirmed that they belong to the same biosynthesis cluster. Systematic gene deletion analysis demonstrated that tf, cadA, mttA, and mfsA genes in the cluster are essential for IA production. Deletion of any of them will either completely abolish the IA production or dramatically decrease amount of IA produced. More importantly, a significant amount of aconitic acid was detected in the cadA deletion strain but not in the other deletion strains. Therefore, the cadA deletion strain is a novel microorganism host for production of aconitic acid and other value-added chemicals from lignocellulosic biomass.

Deng, Shuang↗
DOE OSTI2024

A Sorghum F-Box Protein Induces an Oxidative Burst in the Defense Against Colletotrichum sublineola

The hemibiotrophic fungal pathogen Colletotrichum sublineola is the causal agent of anthracnose in sorghum ( Sorghum bicolor), resulting in leaf blight, stalk rot, and head blight in susceptible genotypes, with yield losses of up to 50%. The development of anthracnose-resistant cultivars can reduce reliance on fungicides and provide a more sustainable and economical means for disease management. A previous genome-wide association study of the sorghum association panel identified the candidate resistance gene Sobic.005G172300 encoding an F-box protein. To better understand the role of this gene in the defense against C. sublineola, gene expression following infection with C. sublineola was monitored by RNA sequencing in seedlings of sorghum accession SC110, which harbored the resistance allele, and three accessions that harbored a susceptible allele. Only in SC110 did the expression of Sobic.005G172300 increase during the biotrophic phase of infection. Subsequent transcriptome analysis, gene co-expression networks, and gene regulatory networks of inoculated and mock-inoculated seedlings of resistant and susceptible accessions suggest that the increase in expression of Sobic.005G172300 induces an oxidative burst by lowering the concentration of ascorbic acid during the biotrophic phase of infection. Based on gene regulatory network analysis, the protein encoded by Sobic.005G172300 is proposed to target proteins involved in the biosynthesis of ascorbic acid for polyubiquitination through the SCF E3 ubiquitin ligase, causing their degradation via the proteasome.

Plant 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 OSTI2022

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts—twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

transcription start and end sites↗