SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2012

Transcriptome-Wide Changes in Chlamydomonas reinhardtii Gene Expression Regulated by Carbon Dioxide and the CO2-Concentrating Mechanism Regulator CIA5/CCM1

We used RNA sequencing to query the Chlamydomonas reinhardtii transcriptome for regulation by CO2 and by the transcription regulator CIA5 (CCM1). Both CO2 and CIA5 are known to play roles in acclimation to low CO2 and in induction of an essential CO2-concentrating mechanism (CCM), but less is known about their interaction and impact on the whole transcriptome. Our comparison of the transcriptome of a wild type versus a cia5 mutant strain under three different CO2 conditions, high CO2 (5%), low CO2 (0.03 to 0.05%), and very low CO2 (< 0.02%), provided an entry into global changes in the gene expression patterns occurring in response to the interaction between CO2 and CIA5. We observed a massive impact of CIA5 and CO2 on the transcriptome, affecting almost 25% of all Chlamydomonas genes, and we discovered an array of gene clusters with distinctive expression patterns that provide insight into the regulatory interaction between CIA5 and CO2. Several individual clusters respond primarily to either CIA5 or CO2, providing access to genes regulated by one factor but decoupled from the other. Three distinct clusters clearly associated with CCM-related genes may represent a rich source of candidates for new CCM components, including a small cluster of genes encoding putative inorganic carbon transporters.

Fang, W↗

DOE OSTI2024

Cell‐type‐specific transcriptomics uncovers spatial regulatory networks in bioenergy sorghum stems

SUMMARY Bioenergy sorghum is a low‐input, drought‐resilient, deep‐rooting annual crop that has high biomass yield potential enabling the sustainable production of biofuels, biopower, and bioproducts. Bioenergy sorghum's 4–5 m stems account for ~80% of the harvested biomass. Stems accumulate high levels of sucrose that could be used to synthesize bioethanol and useful biopolymers if information about cell‐type gene expression and regulation in stems was available to enable engineering. To obtain this information, laser capture microdissection was used to isolate and collect transcriptome profiles from five major cell types that are present in stems of the sweet sorghum Wray. Transcriptome analysis identified genes with cell‐type‐specific and cell‐preferred expression patterns that reflect the distinct metabolic, transport, and regulatory functions of each cell type. Analysis of cell‐type‐specific gene regulatory networks (GRNs) revealed that unique transcription factor families contribute to distinct regulatory landscapes, where regulation is organized through various modes and identifiable network motifs. Cell‐specific transcriptome data was combined with known secondary cell wall (SCW) networks to identify the GRNs that differentially activate SCW formation in vascular sclerenchyma and epidermal cells. The spatial transcriptomic dataset provides a valuable source of information about the function of different sorghum cell types and GRNs that will enable the engineering of bioenergy sorghum stems, and an interactive web application developed during this project will allow easy access and exploration of the data ( https://mc‐lab.shinyapps.io/lcm‐dataset/ ).

09 BIOMASS FUELS↗

DOE OSTI2015

Genomics and Transcriptomics Analyses of the Oil-Accumulating Basidiomycete Yeast Trichosporon oleaginosus : Insights into Substrate Utilization and Alternative Evolutionary Trajectories of Fungal Mating Systems

Microbial fermentation of agro-industrial waste holds great potential for reducing the environmental impact associated with the production of lipids for industrial purposes from plant biomass. However, the chemical complexity of many residues currently prevents efficient conversion into lipids, creating a high demand for strains with the ability to utilize all energy-rich components of agricultural residues. Here, we present results of genome and transcriptome analyses of Trichosporon oleaginosus. This oil-accumulating yeast is able to grow on a wide variety of substrates, including pentoses and N-acetylglucosamine, making it an interesting candidate for biotechnological applications. Transcriptomics shows specific changes in gene expression patterns under lipid-accumulating conditions. Furthermore, gene content and expression analyses indicate that T. oleaginosus is well-adapted for the utilization of chitin-rich biomass. We also focused on the T. oleaginosus mating type, because this species is a member of the Tremellomycetes, a group that has been intensively analyzed as a model for the evolution of sexual development, the best-studied member being Cryptococcus neoformans. The structure of the T. oleaginosus mating-type regions differs significantly from that of other Tremellomycetes and reveals a new evolutionary trajectory paradigm. Comparative analysis shows that recruitment of developmental genes to the ancestral tetrapolar mating-type loci occurred independently in the Trichosporon and Cryptococcus lineages, supporting the hypothesis of a trend toward larger mating-type regions in fungi. IMPORTANCE: Finite fossil fuel resources pose sustainability challenges to society and industry. Microbial oils are a sustainable feedstock for biofuel and chemical production that does not compete with food production. We describe genome and transcriptome analyses of the oleaginous yeast Trichosporon oleaginosus, which can accumulate up to 70% of its dry weight as lipids. In contrast to conventional yeasts, this organism not only shows an absence of diauxic effect while fermenting hexoses and pentoses but also effectively utilizes xylose and N-acetylglucosamine, which are building blocks of lignocellulose and chitin, respectively. Transcriptome analysis revealed metabolic networks that govern conversion of xylose or N-acetylglucosamine as well as lipid accumulation. These data form the basis for a targeted strain optimization strategy. Furthermore, analysis of the mating type of T. oleaginosus supports the hypothesis of a trend toward larger mating-type regions in fungi, similar to the evolution of sex chromosomes in animals and plants.

10 SYNTHETIC FUELS↗

DOE OSTI2013

Global transcriptome analysis of Clostridium thermocellum ATCC 27405 during growth on dilute acid pretreated Populus and switchgrass

Background The thermophilic anaerobe Clostridium thermocellum is a candidate consolidated bioprocessing (CBP) biocatalyst for cellulosic ethanol production. The aim of this study was to investigate C. thermocellum genes required to ferment biomass substrates and to conduct a robust comparison of DNA microarray and RNA sequencing (RNA-seq) analytical platforms. Results C. thermocellum ATCC 27405 fermentations were conducted with a 5 g/L solid substrate loading of either pretreated switchgrass or Populus. Quantitative saccharification and inductively coupled plasma emission spectroscopy (ICP-ES) for elemental analysis revealed composition differences between biomass substrates, which may have influenced growth and transcriptomic profiles. High quality RNA was prepared for C. thermocellum grown on solid substrates and transcriptome profiles were obtained for two time points during active growth (12 hours and 37 hours postinoculation). A comparison of two transcriptomic analytical techniques, microarray and RNA-seq, was performed and the data analyzed for statistical significance. Large expression differences for cellulosomal genes were not observed. We updated gene predictions for the strain and a small novel gene, Cthe_3383, with a putative AgrD peptide quorum sensing function was among the most highly expressed genes. RNAseq data also supported different small regulatory RNA predictions over others. The DNA microarray gave a greater number (2,351) of significant genes relative to RNA-seq (280 genes when normalized by the kernel density mean of M component (KDMM) method) in an analysis of variance (ANOVA) testing method with a 5 % false discovery rate (FDR). When a 2-fold difference in expression threshold was applied, 73 genes were significantly differentially expressed in common between the two techniques. Sulfate and phosphate uptake/utilization genes, along with genes for a putative efflux pump system were some of the most differentially regulated transcripts when profiles for C. thermocellum grown on either pretreated switchgrass or Populus were compared. Conclusions Our results suggest that a high degree of agreement in differential gene expression measurements between transcriptomic platforms is possible, but choosing an appropriate normalization regime is essential.

BESC↗

DOE OSTI2021

Optimal dimensionality selection for independent component analysis of transcriptomic data

Independent component analysis is an unsupervised machine learning algorithm that separates a set of mixed signals into a set of statistically independent source signals. Applied to high-quality gene expression datasets, independent component analysis effectively reveals both the source signals of the transcriptome as co-regulated gene sets, and the activity levels of the underlying regulators across diverse experimental conditions. Two major variables that affect the final gene sets are the diversity of the expression profiles contained in the underlying data, and the user-defined number of independent components, or dimensionality, to compute. Availability of high-quality transcriptomic datasets has grown exponentially as high-throughput technologies have advanced; however, optimal dimensionality selection remains an open question. We computed independent components across a range of dimensionalities for four gene expression datasets with varying dimensions (both in terms of number of genes and number of samples). We computed the correlation between independent components across different dimensionalities to understand how the overall structure evolves as the number of user-defined components increases. We then measured how well the resulting gene clusters reflected known regulatory mechanisms, and developed a set of metrics to assess the accuracy of the decomposition at a given dimension. We found that over-decomposition results in many independent components dominated by a single gene, whereas under-decomposition results in independent components that poorly capture the known regulatory structure. From these results, we developed a new method, called OptICA, for finding the optimal dimensionality that controls for both over- and under-decomposition. Specifically, OptICA selects the highest dimension that produces a low number of components that are dominated by a single gene. We show that OptICA outperforms two previously proposed methods for selecting the number of independent components across four transcriptomic databases of varying sizes. OptICA avoids both over-decomposition and under-decomposition of transcriptomic datasets resulting in the best representation of the organism’s underlying transcriptional regulatory network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

De novo transcriptome analysis of white teak (Gmelina arborea Roxb) wood reveals critical genes involved in xylem development and secondary metabolism

Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees. For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR. We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2014

Comparative transcriptomics between Synechococcus PCC 7942 and Synechocystis PCC 6803 provide insights into mechanisms of adaptation to stress.

Background: Synechococcus sp. PCC 7942 and Synechocystis sp. PCC 6803 are model cyanobacteria from which the metabolism and adaptive responses of other cyanobacteria are inferred. Here we report the gene expression response of these two strains to a variety of nutrient and environmental stresses of varying duration, using transcriptomics. Our data comprise both stranded and 5' enriched libraries in order to elucidate many aspects of the transcriptome. Results: Both organisms were exposed to stress conditions due to nutrient deficiency (inorganic carbon) or change of environmental conditions (salinity, temperature, pH, light) sampled at 1 and 24 hours after the application of stress. The transcriptome profile of each strain revealed similarities and differences in gene expression for photosynthetic and respiratory electron transport chains and carbon fixation. Transcriptome profiles also helped us improve the structural annotation of the genome and identify possible missed genes (including anti-sense) and determine transcriptional units (operons). Finally, we predicted association of proteins of unknown function biochemical pathways by associating them to well-characterized ones based on their transcript levels correlation. Conclusions: Overall, this study results an informative annotation of those species and the comparative analysis of the response of the two organisms revealed similarities but also significant changes in the way they respond to external stress and the duration of the response

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Plant genotype and rhizobia strain combinations strongly influence the transcriptome under heavy metal stress conditions in Medicago truncatula

Heavy metals such as cadmium (Cd) and mercury (Hg) pose significant threats to plant health and food safety as they are absorbed from the environment. Legumes are generally considered sensitive to heavy metals but possess standing genetic variation for accumulation and tolerance to toxic ions. We conducted a transcriptomic analysis on hydroponically and soil grown Medicago truncatula plants to investigate gene expression responses to Cd and Hg exposure in roots, leaves, and nodules. By using plant genotypes with varying metal tolerance or accumulation levels, we observed distinct clustering of gene ontologies, indicating tissue-specific, genotype-specific, and metal-specific gene expression patterns. Considering the symbiotic relationship between legumes and nitrogen-fixing bacteria, we further examined plant phenotypes and transcriptomes of plant genotypes with contrasting Hg accumulation levels and inoculated them with high or low Hg-tolerant Sinorhizobium medicae strains that have presence-absence variation for a mercury reductase (Mer) operon. Host plants inoculated with the Hg-tolerant rhizobia strain possessing a Mer operon exhibited less reduction in nodule number and plant biomass. A smaller reduction in iron (Fe) distribution in nodules after Hg stress was measured using X-ray Fluorescence (XRF) imaging. Dual transcriptome (host plant and bacteria) analysis of nodules revealed a remarkable decrease in the number of differentially expressed genes (DEGs) and clustering of gene ontologies in plants inoculated with the Hg-tolerant rhizobia strain, including symbiosis related genes. This finding suggests that the Hg-tolerant rhizobia strain has the potential to mitigate Hg stress in host plants. Furthermore, we observed genotype by-genotype interactions between the high Hg accumulating plant genotype and the Hg-tolerant rhizobia strain. These findings provide insights into enhancing plant resilience in contaminated environments through optimizing legume-rhizobia interactions for heavy metal tolerance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Parallel measurement of transcriptomes and proteomes from same single cells using nanodroplet splitting

Single-cell multiomics provides comprehensive insights into gene regulatory networks, cellular diversity, and temporal dynamics. Here, we introduce nanoSPLITS (nanodroplet SPlitting for Linked-multimodal Investigations of Trace Samples), an integrated platform that enables global profiling of the transcriptome and proteome from same single cells via RNA sequencing and mass spectrometry-based proteomics, respectively. Benchmarking of nanoSPLITS demonstrates high measurement precision with deep proteomic and transcriptomic profiling of single-cells. We apply nanoSPLITS to cyclin-dependent kinase 1 inhibited cells and found phospho-signaling events could be quantified alongside global protein and mRNA measurements, providing insights into cell cycle regulation. We extend nanoSPLITS to primary cells isolated from human pancreatic islets, introducing an efficient approach for facile identification of unknown cell types and their protein markers by mapping transcriptomic data to existing large-scale single-cell RNA sequencing reference databases. Accordingly, we establish nanoSPLITS as a multiomic technology incorporating global proteomics and anticipate the approach will be critical to furthering our understanding of biological systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Spatial co-transcriptomics reveals discrete stages of the arbuscular mycorrhizal symbiosis

Abstract The symbiotic interaction of plants with arbuscular mycorrhizal (AM) fungi is ancient and widespread. Plants provide AM fungi with carbon in exchange for nutrients and water, making this interaction a prime target for crop improvement. However, plant–fungal interactions are restricted to a small subset of root cells, precluding the application of most conventional functional genomic techniques to study the molecular bases of these interactions. Here we used single-nucleus and spatial RNA sequencing to explore both Medicago truncatula and Rhizophagus irregularis transcriptomes in AM symbiosis at cellular and spatial resolution. Integrated, spatially registered single-cell maps revealed infected and uninfected plant root cell types. We observed that cortex cells exhibit distinct transcriptome profiles during different stages of colonization by AM fungi, indicating dynamic interplay between both organisms during establishment of the cellular interface enabling successful symbiosis. Our study provides insight into a symbiotic relationship of major agricultural and environmental importance and demonstrates a paradigm combining single-cell and spatial transcriptomics for the analysis of complex organismal interactions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Transcriptomic and proteomic signatures of stemness and differentiation in the colon crypt

The intestinal epithelium is replaced weekly by non-quiescent stem cells with kinetics that rely on a rapid loss of stemness and choice for secretory or absorptive lineage differentiation. To determine how the cellular transcriptome and proteome changes during these transitions, we developed a new cell sorting method to purify stem cells, secretory and absorptive progenitor cells, and mature, differentiated cells. Transcriptome analyses revealed that as stem cells transition to the progenitor stage, alternative mRNA splicing and polyadenylation dominate changes in the transcriptome. In contrast, as progenitors differentiate into mature cell types, alterations in gene expression and mRNA levels drive the changes. RNA processing targets mRNAs encoding regulators of cell cycle, RNA regulators, cell adhesion, SUMOylation, and Wnt and Notch signaling. Additionally, global proteome profiling of sorted cell populations detected >2,800 proteins and revealed RNA:protein patterns of abundance and correlation. Paired together, these data highlight new potentials for autocrine and feedback regulation and provide new insights into cell state transitions in the crypt.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Compilation and utilization of a sorghum transcriptome compendium for gene regulatory network analysis and crop trait engineering

Sorghum bicolor (Sorghum) is a drought and heat tolerant C4 grass crop used to produce grain, forage, biofuels, and other bioproducts. Genetic improvement of sorghum hybrid crops is aided by a large and diverse germplasm, sorghum's diploid inbreeding genetics, and a relatively small genome that has facilitated genomic research. Over the past 20 years, the sorghum research community characterized the cytogenetic and recombinant landscapes of sorghum's 10 chromosomes, sequenced and annotated the sorghum genome, and used that information to identify genes/alleles that modulate flowering time, plant height, seed shattering, and other important traits. More recently, >1000 RNA-seq transcriptome profiles were collected from 15 sorghum genotypes to help understand the genetic basis of variation in growth and development of sorghum stems, tillers, roots, and leaves, and the regulation of biosynthetic pathways that produce epicuticular wax, dhurrin, and RFOs, compounds that contribute to sorghum's resilience. Transcriptome studies were designed to identify differentially expressed genes that are co-expressed during development or in response to a treatment to enable construction of gene regulatory networks. Co-expression and network analysis identified transcription factors and their cognate binding sites in target gene promoters and signaling pathways that modulate gene regulatory networks providing gene editing targets for further trait optimization. RNA-seq data from >20 experiments targeting sorghum organs, tissues, cell types, developmental stages, and responses to environmental conditions (i.e., diel, day-length, shading, water-deficit, temperature) has been compiled in a sorghum transcriptome compendium. The goal of this resource paper is to describe compendium content, accessibility, and a compendium data analysis pipeline and to illustrate the types of information that can be derived from the compendium with a focus on the elucidation of gene regulatory networks useful for guiding the improvement of sorghum traits through gene editing.

RNA-seq↗

DOE OSTI2021

Plant Cell Identity in the Era of Single-Cell Transcriptomics

High-throughput single-cell transcriptomic approaches have revolutionized our view of gene expression at the level of individual cells, providing new insights into their heterogeneity, identities, and functions. Recently, technical challenges to the application of single-cell transcriptomics to plants have been overcome, and many plant organs and tissues have now been subjected to analyses at single-cell resolution. In this review, we describe these studies and their impact on our understanding of the diversity, differentiation, and activities of plant cells. We particularly highlight their impact on plant cell identity, including unprecedented views of cell transitions and definitions of rare and novel cell types. Lastly, we also point out current challenges and future opportunities for the application and analyses of single-cell transcriptomics in plants.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Global transcriptome profile of the developmental principles of in vitro iPSC-to-motor neuron differentiation

Abstract Background Human induced pluripotent stem cells (iPSC) have opened new avenues for regenerative medicine. Consequently, iPSC-derived motor neurons have emerged as potentially viable therapies for spinal cord injuries and neurodegenerative disorders including Amyotrophic Lateral Sclerosis. However, direct clinical application of iPSC bears in itself the risk of tumorigenesis and other unforeseeable genetic or epigenetic abnormalities. Results Employing RNA-seq technology, we identified and characterized gene regulatory networks triggered by in vitro chemical reprogramming of iPSC into cells with the molecular features of motor neurons (MNs) whose function in vivo is to innervate effector organs. We present meta-transcriptome signatures of 5 cell types: iPSCs, neural stem cells, motor neuron progenitors, early motor neurons, and mature motor neurons. In strict response to the chemical stimuli, along the MN differentiation axis we observed temporal downregulation of tumor growth factor-β signaling pathway and consistent activation of sonic hedgehog, Wnt/β-catenin, and Notch signaling. Together with gene networks defining neuronal differentiation (neurogenin 2, microtubule-associated protein 2, Pax6, and neuropilin-1), we observed steady accumulation of motor neuron-specific regulatory genes, including Islet-1 and homeobox protein HB9. Interestingly, transcriptome profiling of the differentiation process showed that Ca 2+ signaling through cAMP and LPC was downregulated during the conversion of the iPSC to neural stem cells and key regulatory gene activity of the pathway remained inhibited until later stages of motor neuron formation. Pathways shaping the neuronal development and function were well-represented in the early motor neuron cells including, neuroactive ligand-receptor interactions, axon guidance, and the cholinergic synapse formation. A notable hallmark of our in vitro motor neuron maturation in monoculture was the activation of genes encoding G-coupled muscarinic acetylcholine receptors and downregulation of the ionotropic nicotinic acetylcholine receptors expression. We observed the formation of functional neuronal networks as spontaneous oscillations in the extracellular action potentials recorded on multi-electrode array chip after 20 days of differentiation. Conclusions Detailed transcriptome profile of each developmental step from iPSC to motor neuron driven by chemical induction provides the guidelines to novel therapeutic approaches in the re-construction efforts of muscle innervation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

A transcriptome software comparison for the analyses of treatments expected to give subtle gene expression responses

Background: In this comparative study we evaluate the performance of four software tools: DNAstar-D (DESeq2), DNAstar-E (edgeR), CLC Genomics and Partek Flow for identification of differentially expressed genes (DEGs) using a transcriptome of E. coli. The RNA-seq data are from the effect of below-background radiation 5.5 nGy total dose (0.2nGy/hr) on E. coli grown shielded from natural radiation 655 m below ground in a pre-World War II steel vault. The gene expression response to three supplemented sources of radiation designed to mimic natural background, 1952 – 5720 nGy in total dose (71–208 nGy/hr), are compared to this “radiation-deprived” treatment. In addition, RNA-seq data of Caenorhabditis elegans nematode from similar radiation treatments was analyzed by three of the software packages. Results: In E. coli, the four software programs identified one of the supplementary sources of radiation (KCl) to evoke about 5 times more transcribed genes than the minus-radiation treatment (69–114 differentially expressed genes, DEGs), and so the rest of the analyses used this KCl vs “Minus” comparison. After imposing a 30-read minimum cutoff, one of the DNAStar options shared two of the three steps (mapping, normalization, and statistic) with Partek Flow (they both used median of ratios to normalize and the DESeq2 statistical package), and these two programs identified the highest number of DEGs in common with each other (53). In contrast, when the programs used different approaches in each of the three steps, between 31 and 40 DEGs were found in common. Regarding the extent of expression differences, three of the four programs gave high fold-change results (15–178 fold), but one (DNAstar’s DESeq2) resulted in more conservative fold-changes (1.5–3.5). In a parallel study comparing three qPCR commercial validation software programs, these programs also gave variable results as to which genes were significantly regulated. Similarly, the C. elegans analysis showed exaggerated fold-changes in CLC and DNAstar’s edgeR while DNAstar-D was more conservative. Conclusions: Regarding the extent of expression (fold-change), and considering the subtlety of the very low level radiation treatments, in E. coli three of the four programs gave what we consider exaggerated fold-change results (15 – 178 fold), but one (DNAstar’s DESeq2) gave more realistic fold-changes (1.5–3.5). When RT-qPCR validation comparisons to transcriptome results were carried out, they supported the more conservative DNAstar-D’s expression results. When another model organism’s (nematode) response to these radiation differences was similarly analyzed, DNAstar-D also resulted in the most conservative expression patterns. Therefore, we would propose DESeq2 (“DNAstar-D”) as an appropriate software tool for differential gene expression studies for treatments expected to give subtle transcriptome responses.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

Proteomics and transcriptomics analyses of Arabidopsis floral buds uncover important functions of ARABIDOPSIS SKP1-LIKE1

The ARABIDOPSIS SKP1-LIKE1 (ASK1) protein functions as a subunit of SKP1-CUL1-F-box (SCF) E3 ubiquitin ligases. Previous genetic studies showed that ASK1 plays important roles in Arabidopsis flower development and male meiosis. However, the molecular impact of ASK1-containing SCF E3 ubiquitin ligases (ASK1-E3s) on the floral proteome and transcriptome is unknown. Here we identified proteins that are potentially regulated by ASK1-E3s by comparing floral bud proteomes of wild-type and the ask1 mutant plants. More than 200 proteins were detected in the ask1 mutant but not in wild-type and >300 were detected at higher levels in the ask1 mutant than in wild-type, but their RNA levels were not significantly different between wild-type and ask1 floral buds as shown by transcriptomics analysis, suggesting that they are likely regulated at the protein level by ASK1-E3s. Integrated analyses of floral proteomics and transcriptomics of ask1 and wild-type uncovered several potential aspects of ASK1-E3 functions, including regulation of transcription regulators, kinases, peptidases, and ribosomal proteins, with implications on possible mechanisms of ASK1-E3 functions in floral development. In conclusion, our results suggested that ASK1-E3s play important roles in Arabidopsis protein degradation during flower development. This study opens up new possibilities for further functional studies of these candidate E3 substrates.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Genome sequencing and transcriptome analysis of Trichoderma reesei QM9978 strain reveals a distal chromosome translocation to be responsible for loss of vib1 expression and loss of cellulase induction

Background: The hydrolysis of biomass to simple sugars used for the production of biofuels in biorefneries requires the action of cellulolytic enzyme mixtures. During the last 50 years, the ascomycete Trichoderma reesei, the main source of industrial cellulase and hemicellulase cocktails, has been subjected to several rounds of classical mutagenesis with the aim to obtain higher production levels. During these random genetic events, strains unable to produce cellulases were generated. Here, whole genome sequencing and transcriptomic analyses of the cellulase-negative strain QM9978 were used for the identifcation of mutations underlying this cellulase-negative phenotype. Results: Sequence comparison of the cellulase-negative strain QM9978 to the reference strain QM6a identifed a total of 43 mutations, of which 33 were located either close to or in coding regions. From those, we identifed 23 single-nucleotide variants, nine InDels, and one translocation. The translocation occurred between chromosomes V and VII, is located upstream of the putative transcription factor vib1, and abolishes its expression in QM9978 as detected during the transcriptomic analyses. Ectopic expression of vib1 under the control of its native promoter as well as overexpression of vib1 under the control of a strong constitutive promoter restored cellulase expression in QM9978, thus confrming that the translocation event is the reason for the cellulase-negative phenotype. Gene deletion of vib1 in the moderate producer strain QM9414 and in the high producer strain Rut-C30 reduced cellulase expression in both cases. Overexpression of vib1 in QM9414 and Rut-C30 had no efect on cellulase production, most likely because vib1 is already expressed at an optimal level under normal conditions. Conclusion: We were able to establish a link between a chromosomal translocation in QM9978 and the cellulasenegative phenotype of the strain. We identifed the transcription factor vib1 as a key regulator of cellulases in T. reesei whose expression is absent in QM9978. We propose that in T. reesei, as in Neurospora crassa, vib1 is involved in cellulase induction, although the exact mechanism remains to be elucidated. The data presented here show an example of a combined genome sequencing and transcriptomic approach to explain a specifc trait, in this case the QM9978 cellulase-negative phenotype, and how it helps to better understand the mechanisms during cellulase gene regulation. When focusing on mutations on the single base-pair level, changes on the chromosome level can be easily overlooked and through this work we provide an example that stresses the importance of the big picture of the genomic landscape during analysis of sequencing data.

09 BIOMASS FUELS↗

DOE OSTI2018

Revealing the transcriptomic complexity of switchgrass by PacBio long-read sequencing

Switchgrass ( Panicum virgatum L.) is an important bioenergy crop widely used for lignocellulosic research. While extensive transcriptomic analyses have been conducted on this species using short read-based sequencing techniques, very little has been reliably derived regarding alternatively spliced (AS) transcripts. Here, we present an analysis of transcriptomes of six switchgrass tissue types pooled together, sequenced using Pacific Biosciences (PacBio) single-molecular long-read technology. Our analysis identified 105,419 unique transcripts covering 43,570 known genes and 8795 previously unknown genes. 45,168 are novel transcripts of known genes. A total of 60,096 AS transcripts are identified, 45,628 being novel. We have also predicted 1549 transcripts of genes involved in cell wall construction and remodeling, 639 being novel transcripts of known cell wall genes. Most of the predicted transcripts are validated against Illumina-based short reads. Specifically, 96% of the splice junction sites in all the unique transcripts are validated by at least five Illumina reads. Comparisons between genes derived from our identified transcripts and the current genome annotation revealed that among the gene set predicted by both analyses, 16,640 have different exon-intron structures. Overall, substantial amount of new information is derived from the PacBio RNA data regarding both the transcriptome and the genome of switchgrass.

09 BIOMASS FUELS↗