Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transcriptome 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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

pmartR 2.0 : A Quality Control, Visualization, and Statistics Pipeline for Multiple Omics Datatypes

The pmartR (https://github.com/pmartR/pmartR) package was designed for the quality control (QC) and analysis of mass spectrometry (MS) data, tailored to specific characteristics of proteomic (isobaric or labelled), metabolomic, and lipidomic datasets. Since its initial release, the tool has been expanded to address the needs of its growing userbase and now includes QC and statistics for nuclear magnetic resonance (NMR) metabolomic data, and leverages the DESeq2, edgeR, and limma-voom R packages for some transcriptomic data analyses. These improvements have made progress towards a unified omics processing pipeline for ease of reporting and streamlined statistical purposes. The package’s statistics and visualization capabilities have also been expanded by adding support for paired data and by integrating pmartR with the trelliscopejs R package for the quick creation of trellis displays (https://github.com/hafen/trelliscopejs). Here, we present relevant examples of each of these enhancements to pmartR and highlight how each new feature benefits the omics community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Widespread Polycistronic Transcripts in Fungi Revealed by Single-Molecule mRNA Sequencing

Genes in prokaryotic genomes are often arranged into clusters and co-transcribed into polycistronic RNAs. Isolated examples of polycistronic RNAs were also reported in some higher eukaryotes but their presence was generally considered rare. Here we developed a longread sequencing strategy to identify polycistronic transcripts in several mushroom forming fungal species including Plicaturopsis crispa, Phanerochaete chrysosporium, Trametes versicolor, and Gloeophyllum trabeum. We found genome-wide prevalence of polycistronic transcription in these Agaricomycetes, involving up to 8% of the transcribed genes. Unlike polycistronic mRNAs in prokaryotes, these co-transcribed genes are also independently transcribed. We show that polycistronic transcription may interfere with expression of the downstream tandem gene. Further comparative genomic analysis indicates that polycistronic transcription is conserved among a wide range of mushroom forming fungi. In summary, our study revealed, for the first time, the genome prevalence of polycistronic transcription in a phylogenetic range of higher fungi. Furthermore, we systematically show that our long-read sequencing approach and combined bioinformatics pipeline is a generic powerful tool for precise characterization of complex transcriptomes that enables identification of mRNA isoforms not recovered via short-read assembly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Two related families of metal transferases, ZNG1 and ZNG2, are involved in acclimation to poor Zn nutrition in Arabidopsis

Metal homeostasis has evolved to tightly modulate the availability of metals within the cell, avoiding cytotoxic interactions due to excess and protein inactivity due to deficiency. Even in the presence of homeostatic processes, however, low bioavailability of these essential metal nutrients in soils can negatively impact crop health and yield. While research has largely focused on how plants assimilate metals, acclimation to metal-limited environments requires a suite of strategies that are not necessarily involved in metal transport across membranes. The identification of these mechanisms provides a new opportunity to improve metal-use efficiency and develop plant foodstuffs with increased concentrations of bioavailable metal nutrients. Here, we investigate the function of two distinct subfamilies of the nucleotide-dependent metallochaperones (NMCs), named ZNG1 and ZNG2, that are found in plants, using Arabidopsis thaliana as a reference organism. AtZNG1 (AT1G26520) is an ortholog of human and fungal ZNG1, and like its previously characterized eukaryotic relatives, localizes to the cytosol and physically interacts with methionine aminopeptidase type I (AtMAP1A). Analysis of At ZNG1 , At MAP1A , At MAP2A , and At MAP2B transgenic mutants are consistent with the role of Arabidopsis ZNG1 as a Zn transferase for AtMAP1A, as previously described in yeast and zebrafish. Structural modeling reveals a flexible cysteine-rich loop that we hypothesize enables direct transfer of Zn from AtZNG1 to AtMAP1A during GTP hydrolysis. Based on proteomics and transcriptomics, loss of this ancient and conserved mechanism has pleiotropic consequences impacting the expression of hundreds of genes, including those involved in photosynthesis and vesicle transport. Members of the plant-specific family of NMCs, ZNG2A1 (AT1G80480) and ZNG2A2 (AT1G15730), are also required during Zn deficiency, but their target protein(s) remain to be discovered. In conclusion, RNA-seq analyses reveal wide-ranging impacts across the cell when the genes encoding these plastid-localized NMCs are disrupted.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

'Omics analysis of low dose acetaminophen intake demonstrates novel response pathways in humans

Acetaminophen is the primary cause of acute liver toxicity in Europe/USA, which led the FDA to reconsider recommendations concerning safe acetaminophen dosage/use. Unfortunately, the current tests for liver toxicity are no ideal predictive markers for liver injury, i.e. they only measure acetaminophen exposure after profound liver toxicity has already occurred. Furthermore, these tests do not provide mechanistic information. Here, 'omics techniques (global analysis of metabolomic/gene-expression responses) may provide additional insight. To better understand acetaminophen-induced responses at low doses, we evaluated the effects of (sub-)therapeutic acetaminophen doses on metabolite formation and global gene-expression changes (including, for the first time, full-genome human miRNA expression changes) in blood/urine samples from healthy human volunteers. Many known and several new acetaminophen-metabolites were detected, in particular in relation to hepatotoxicity-linked, oxidative metabolism of acetaminophen. Transcriptomic changes indicated immune-modulating effects (2 g dose) and oxidative stress responses (4 g dose). For the first time, effects of acetaminophen on full-genome human miRNA expression have been considered and confirmed the findings on mRNA level. 'Omics techniques outperformed clinical chemistry tests and revealed novel response pathways to acetaminophen in humans. Although no definitive conclusion about potential immunotoxic effects of acetaminophen can be drawn from this study, there are clear indications that the immune system is triggered even after intake of low doses of acetaminophen. Also, oxidative stress-related gene responses, similar to those seen after high dose acetaminophen exposure, suggest the occurrence of possible pre-toxic effects of therapeutic acetaminophen doses. Possibly, these effects are related to dose-dependent increases in levels of hepatotoxicity-related metabolites. -- Highlights: ► 'Omics techniques outperformed classic clinical chemistry tests. ► Metabolomic analyses led to the detection of five new acetaminophen metabolites. ► Low dose APAP changed immune and oxidative stress related gene expression in blood. ► APAP-induced full-genome human blood miRNA profiles were assessed for the first time.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Lipidomic profiling of plasma samples from patients with mitochondrial disease

Highlights: • All 14 lysoPCs were dramatically decreased in MD patients. • Acylcarnitines, TGs, PIs, PEs, unsaturated lysoPEs and PCs were significantly increased in MD patients. • PI 38:6, lysoPCs (20:0, 19:0, 18:0, 17:0) may be potentially sensitive and specific biomarkers for MD. Mitochondrial disease (MD) is a rare mitochondrial respiratory chain disorder with a high mortality and extremely challenging to treat. Although genomic, transcriptomic, and proteomic analyses have been performed to investigate the pathogenesis of MD, the role of metabolomics in MD, particularly of lipidomics remains unclear. This study was undertaken to identify potential lipid biomarkers of MD. An untargeted lipidomic approach was used to compare the plasma lipid metabolites in 20 MD patients and 20 controls through Liquid Chromatography coupled to Mass Spectrometry. Volcano plot analysis was performed to identify the different metabolites. Receiver operating characteristic (ROC) curves were constructed and the area under the ROC curves (AUC) was calculated to determine the potentially sensitive and specific biomarkers. A total of 41 lipids were significantly different in MD patients and controls. ROC curve analysis showed the top 5 AUC values of lipids (phosphatidylinositols 38:6, lysoPC 20:0, 19:0, 18:0, 17:0) are more than 0.99. Multivariate ROC curve based exploratory analysis showed the AUC of combination of top 5 lipids is 1, indicating they may be potentially sensitive and specific biomarkers for MD. We propose combination of these lipid species may be more valuable in predicting the development and progression of MD, and this will have important implications for the diagnosis and treatment of MD.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Differential gene expression profiling of human adipose stem cells differentiating into smooth muscle-like cells by TGFβ1/BMP4

Regenerative repair of the vascular system is challenging from the perspectives of translational medicine and tissue engineering. There are fundamental hurdles in front of creating bioartificial arteries, which involve recaputilation of the three-layered structure under laboratory settings. Obtaining and maintaining smooth muscle characteristics is an important limitation, as the transdifferentiated cells fail to display mature phenotype. This study aims to shed light on the smooth muscle differentiation of human adipose stem cells (hASCs). To this end, we first acquired hASCs from lipoaspirate samples. Upon characterization, the cells were induced to differentiate into smooth muscle (SM)-like cells using a variety of inducer combinations. Among all, TGFβ1/BMP4 combination had the highest differentiation efficiency, based on immunohistochemical analyses. hSM-like cell samples were compared to hASCs and to the positive control, human coronary artery-smooth muscle cells (hCA-SMCs) through gene transcription profiling. Microarray findings revealed the activation of gene groups that function in smooth muscle differentiation, signaling pathways, extracellular modeling and cell proliferation. Our results underline the effectiveness of the growth factors and suggest some potential variables for detecting the SM-like cell characteristics. Evidence in transcriptome level was used to evaluate the TGFβ1/BMP4 combination as a previously unexplored effector for the smooth muscle differentiation of adipose stem cells. - Highlights: • Human adipose stem cells (hASCs) were isolated, characterized and cultured. • Growth factor combinations were evaluated for their effectiveness in differentiation using IHC. • hASCs were differentiated into smooth muscle (SM)-like cells using TGF-β1 and BMP4 combination. • Microarray analysis was performed for hASCs, SM-like cells and coronary artery-SMCs. • Microarray data was used to perform hierarchical clustering and interpretation of activated pathways.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Analysis of clock-regulated genes in Neurospora reveals widespread posttranscriptional control of metabolic potential

Significance Circadian clocks regulate gene expression levels to allow an organism to anticipate environmental conditions. These clocks reside in all the major branches of life and confer a competitive advantage to the organisms that maintain them. The clock in the fungus Neurospora crassa is an excellent model for basic understanding of core circadian architecture as well as for filamentous fungi. Here, we identify genes whose expression is clock regulated; indeed, as much as 40% of the transcriptome may be clock regulated, broadly directing daytime catabolism and nighttime growth. Both transcriptional control and posttranscriptional regulation play major roles in control of cycling transcripts such that DNA binding of transcription factors alone appears insufficient to set the phase of circadian transcription.

Hurley, Jennifer M.↗
DOE OSTI2013

Utilizing toxicogenomic data to understand chemical mechanism of action in risk assessment

The predominant role of toxicogenomic data in risk assessment, thus far, has been one of augmentation of more traditional in vitro and in vivo toxicology data. This article focuses on the current available examples of instances where toxicogenomic data has been evaluated in human health risk assessment (e.g., acetochlor and arsenicals) which have been limited to the application of toxicogenomic data to inform mechanism of action. This article reviews the regulatory policy backdrop and highlights important efforts to ultimately achieve regulatory acceptance. A number of research efforts on specific chemicals that were designed for risk assessment purposes have employed mechanism or mode of action hypothesis testing and generating strategies. The strides made by large scale efforts to utilize toxicogenomic data in screening, testing, and risk assessment are also discussed. These efforts include both the refinement of methodologies for performing toxicogenomics studies and analysis of the resultant data sets. The current issues limiting the application of toxicogenomics to define mode or mechanism of action in risk assessment are discussed together with interrelated research needs. In summary, as chemical risk assessment moves away from a single mechanism of action approach toward a toxicity pathway-based paradigm, we envision that toxicogenomic data from multiple technologies (e.g., proteomics, metabolomics, transcriptomics, supportive RT-PCR studies) can be used in conjunction with one another to understand the complexities of multiple, and possibly interacting, pathways affected by chemicals which will impact human health risk assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

The Arabidopsis CSLD 5 functions in cell plate formation in a cell cycle-dependent manner

In plants, the presence of a load-bearing cell wall presents unique challenges during cell division. Unlike other eukaryotes, which undergo contractile cytokinesis upon completion of mitosis, plants instead synthesize and assemble a new dividing cell wall to separate newly formed daughter cells. In this study, we mine transcriptome data from individual cell types in the Arabidopsis thaliana stomatal lineage and identify CSLD5, a member of the Cellulose Synthase Like-D family, as a cell wall biosynthesis enzyme uniquely enriched in rapidly dividing cell populations. We further show that CSLD5 is a direct target of SPEECHLESS, the master transcriptional regulator of these divisions during stomatal development. Using a combination of genetic analysis and in vivo localization of fluorescently tagged fusion proteins, we show that CSLD5 preferentially accumulates in dividing plant cells where it participates in the construction of newly forming cell plates. We show that CSLD5 is an unstable protein that is rapidly degraded upon completion of cell division and that the protein turnover characteristics of CSLD5 are altered in ccs52a2 mutants, indicating that CSLD5 turnover may be regulated by a cell cycle-associated E3-ubiquitin ligase, the anaphase-promoting complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Engineering Anaerobic Gut Fungi for Lignocellulose Breakdown

Biofuels derived from plant biomass (lignocellulose) are attractive alternatives to petroleum-based products, as they avoid many of the economic and environmental problems that plague traditional energy sources. Breakdown of plant waste into simple sugars would translate into the production of renewable fuels via microbial fermentation. However, existing technologies are insufficient to allow for industrial-scale production of these products due to difficulties associated with the recalcitrance of crude lignin-rich biomass, and the high cost/poor performance of known cellulolytic enzymes (cellulases). Therefore, there is a critical need to develop new technologies to break down lignocellulosic biomass into fermentable sugars for downstream fuel development. Towards this goal, much can be learned by studying how anaerobic microbes depolymerize lignocellulose in biomass-rich environments, such as the digestive tract of large herbivores. Anaerobic gut fungi are native to the gut and rumen of these animals, where they have evolved unique abilities to break down lignocellulosic biomass through invasive growth and the secretion of powerful enzymes. The goal of this project was to engineer anaerobic gut fungi as novel platform organisms for biofuel production from plant material. To accomplish this goal, a panel of anaerobic fungi was isolated form different herbivores and screened for their ability to degrade lignocellulose. Four novel strains were isolated and genetically sequenced, including Piromyces finnis, Neocallimastix californiae, Anaeromyces robustus, and Caecomyces churrovis. Transcriptomic and genomic characterization of these anaerobic fungi revealed that they harbor the largest collection of biomass-degrading enzymes of any sequenced fungus, and delineated the basic metabolic networks that govern lignocellulose hydrolysis within anaerobic fungi, RNA sequencing enabled the dynamics and mechanism associated with carbon catabolite repression to be delineated, which revealed how important enzyme groups are coordinated during breakdown. Through cluster analysis, a set of putative novel enzymes associated with lignocellulose breakdown was discovered, which also revealed a set of proteins critical to the formation of fungal cellulosomes (multi-enzyme complexes). Using this information, genetic transformation strategies to manipulate gut fungi were developed, which would endow them with enhanced functionality against a range of industrially relevant substrates. Collectively, this information will establish the molecular framework for anaerobic fungal hydrolysis, and will guide in the development of lignocellulosic biofuels.

09 BIOMASS FUELS↗
DOE OSTI2021

BETO 2021 Peer Review - Improving Formate Upgrading by Cupriavidus Necator 2.3.2.111

Formate can be generated by electrocatalytic reduction of CO2 and has been proposed as a soluble intermediate for the storage of carbon and energy. Biological systems capable of assimilating formate could enable conversion of formate generated from low-cost renewable energy and waste CO2 to myriad fuels and chemicals. To that end, the goal of this project is to develop the natural formatotroph Cupriavidus necator as a robust microbial chassis for efficient conversion of formate to value-added products. Using a combination of laboratory evolution and rational engineering, we have identified mutations and genetic deletions that improve its growth rate and biomass yield on formate by 1.2X. Using RNA-seq transcriptomics, we are elucidating the mechanisms that underly these improvements and anticipate this will contribute to our understanding of formate conversion. We have also engineered C. necator to convert formate to the potential polymer precursor 2-hydroxymuconate semialdehyde (2HMS) by introducing an exogenous pathway and deleting competing pathways and demonstrated production and robust growth using a pH-stat bioreactor. Final techno-economic analysis and life-cycle assessment on biological production of 2HMS from formate will enable the DOE and other stake holders to evaluate the economic and environmental impact of this project and related technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

Bioprocessing analysis of Pyrococcus furiosus strains engineered for CO 2 -based 3-hydroxypropionate production

In this paper, metabolically engineered strains of the hyperthermophile Pyrococcus furiosus (T opt 95–100°C), designed to produce 3-hydroxypropionate (3HP) from maltose and CO 2 using enzymes from the Metallosphaera sedula (T opt 73°C) carbon fixation cycle, were examined with respect to the impact of heterologous gene expression on metabolic activity, fitness at optimal and sub-optimal temperatures, gas-liquid mass transfer in gas-intensive bioreactors, and potential bottlenecks arising from product formation. Transcriptomic comparisons of wild-type P. furiosus , a genetically-tractable, naturally-competent mutant (COM1), and COM1-based strains engineered for 3HP production revealed numerous differences after being shifted from 95°C to 72°C, where product formation catalyzed by the heterologously-produced M. sedula enzymes occurred. At 72°C, significantly higher levels of metabolic activity and a stress response were evident in 3HP-forming strains compared to the non-producing parent strain (COM1). Gas–liquid mass transfer limitations were apparent, given that 3HP titers and volumetric productivity in stirred bioreactors could be increased over 10-fold by increased agitation and higher CO 2 sparging rates, from 18 mg/L to 276 mg/L and from 0.7 mg/L/h to 11 mg/L/h, respectively. 3HP formation triggered transcription of genes for protein stabilization and turnover, RNA degradation, and reactive oxygen species detoxification. Lastly, the results here support the prospects of using thermally diverse sources of pathways and enzymes in metabolically engineered strains designed for product formation at sub-optimal growth temperatures.

3-hydroxypropionate↗
DOE OSTI2017

Transcriptomic Analyses Elucidate Adaptive Differences of Closely Related Strains of Pseudomonas aeruginosa in Fuel

Pseudomonas aeruginosa can utilize hydrocarbons, but different strains have various degrees of adaptation despite their highly conserved genome. P. aeruginosa ATCC 33988 is highly adapted to hydrocarbons, while P. aeruginosa strain PAO1, a human pathogen, is less adapted and degrades jet fuel at a lower rate than does ATCC 33988. We investigated fuel-specific transcriptomic differences between these strains in order to ascertain the underlying mechanisms utilized by the adapted strain to proliferate in fuel. During growth in fuel, the genes related to alkane degradation, heat shock response, membrane proteins, efflux pumps, and several novel genes were upregulated in ATCC 33988. Overexpression of alk genes in PAO1 provided some improvement in growth, but it was not as robust as that of ATCC 33988, suggesting the role of other genes in adaptation. Expression of the function unknown gene PA5359 from ATCC 33988 in PAO1 increased the growth in fuel. Bioinformatic analysis revealed that PA5359 is a predicted lipoprotein with a conserved Yx(FWY)xxD motif, which is shared among bacterial adhesins. Overexpression of the putative resistance-nodulation-division (RND) efflux pump PA3521 to PA3523 increased the growth of the ATCC 33988 strain, suggesting a possible role in fuel tolerance. Interestingly, the PAO1 strain cannot utilize n -C 8 and n -C 10 . The expression of green fluorescent protein (GFP) under the control of alkB promoters confirmed that alk gene promoter polymorphism affects the expression of alk genes. Promoter fusion assays further confirmed that the regulation of alk genes was different in the two strains. Protein sequence analysis showed low amino acid differences for many of the upregulated genes, further supporting transcriptional control as the main mechanism for enhanced adaptation. Lastly, these results support that specific signal transduction, gene regulation, and coordination of multiple biological responses are required to improve the survival, growth, and metabolism of fuel in adapted strains. This study provides new insight into the mechanistic differences between strains and helpful information that may be applied in the improvement of bacterial strains for resistance to biotic and abiotic factors encountered during bioremediation and industrial biotechnological processes.

54 ENVIRONMENTAL 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 OSTI2024

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. Here, the standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

UnigeneFinder: An Automated Pipeline for Gene Calling From Transcriptome Assemblies Without a Reference Genome

ABSTRACT For most species, transcriptome data are much more readily available than genome data. Without a reference genome, gene calling is cumbersome and inaccurate because of the high degree of redundancy in de novo transcriptome assemblies. To simplify and increase the accuracy of de novo transcriptome assembly in the absence of a reference genome, we developed UnigeneFinder. Combining several clustering methods, UnigeneFinder substantially reduces the redundancy typical of raw transcriptome assemblies. This pipeline offers an effective solution to the problem of inflated transcript numbers, achieving a closer representation of the actual underlying genome. UnigeneFinder performs comparably or better, compared with existing tools, on plant species with varying genome complexities. UnigeneFinder is the only available transcriptome redundancy solution that fully automates the generation of primary transcript, coding region, and protein sequences, analogous to those available for high‐quality reference genomes. These features, coupled with the pipeline’s cross‐platform implementation, focus on automation, and an accessible, user‐friendly interface, make UnigeneFinder a useful tool for many downstream sequence‐based analyses in nonmodel organisms lacking a reference genome, including differential gene expression analysis, accurate ortholog identification, functional enrichments, and evolutionary analyses. UnigeneFinder also runs efficiently both on high‐performance computing (HPC) systems and personal computers, further reducing barriers to use.

Xue, Bo [Plant Resilience Institute Michigan State↗
DOE OSTI2023

Trichoderma harzianum transcriptome in response to the nematode Pratylenchus brachyurus

The root-lesion nematode Pratylenchus brachyurus causes extensive damage in several crops of economic importance. Fungi of Trichoderma genus have been highlighted as biopesticide agents in the control of several plant diseases. Although it is already widely used in agriculture, there are few studies, especially at the molecular level, that evaluate T. harzianum in the control of P. brachyurus . The aim of the present study was to investigate how interaction with the nematode P. brachyurus influences gene expression of T. harzianum , by using RNA-Seq analysis. Of the 13,932 predicted genes in the T. harzianum genome, 2,922 (21%) were differentially expressed in the presence of P. brachyurus , in relation to the absence of the nematode. Among the differentially expressed genes, we found genes encoding Carbohydrate Active EnZymes (CAZy), MEROPS peptidases, and proteins related to secondary metabolite synthesis. 118 pathways were identified as related to the biosynthesis of secondary metabolites. Additionally, the analysis identified 136 metabolic pathways related to these differentially expressed genes, among which we highlight: aminobenzoate degradation, xenobiotic metabolism by cytochrome P450, and sesquiterpenoid and triterpenoid biosynthesis. Our results contribute to a better understanding of the response of T. harzianum to the nematode P. bracyurus , the potential for biocontrol by this fungus.

59 BASIC BIOLOGICAL SCIENCES↗