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 397 records · Page 22

Metabolic flux, metabolite, and transcript analysis uncover reprogramming of metabolism toward higher seed oil

Overexpression of WRINKLED1 (WRI1), a master regulator of glycolysis and fatty acid biosynthesis, together with DIACYLGLYCEROL ACYLTRANSFERASE1 (DGAT1), which catalyzes the final step of triacylglycerol assembly, is a promising strategy for enhancing seed oil content. However, how these regulators coordinate system-wide metabolic reprogramming at the levels of gene expression, metabolite pools, and fluxes remains poorly understood. To address this, we performed 13 C-metabolic flux analysis, metabolomics, and transcriptomics on in vitro cultured pennycress (Thlaspi arvense L.) embryos overexpressing the native WRI1 and DGAT1 homologs. Here, in cultured embryos, WRI1/DGAT1 overexpression increased triacylglycerol accumulation by 28% while reducing protein content by 34%, relative to the wild type. Embryos showed ∼20-fold and 50-fold upregulation of WRI1 and DGAT1 along with induction of WRI1 target genes in glycolysis and fatty acid biosynthesis. Genes associated with photosynthesis and Calvin cycle functions were also upregulated, whereas genes encoding ribosomal proteins and seed storage proteins were strongly repressed, consistent with the observed lipid–protein tradeoff. Flux analysis revealed that enhanced triacylglycerol biosynthesis is supported by increased flux through the Rubisco shunt and cytosolic pyruvate kinase, while the oxidative pentose phosphate pathway and malic enzyme contributed little to NADPH or pyruvate supply. Metabolomic profiling revealed extensive perturbations in glycolytic intermediates, tricarboxylic acid cycle metabolites, and amino acids. In plant grown seeds, WRI1/DGAT1 lines also showed a modest but significant increase in total lipid content. Collectively, these findings reveal how WRI1 and DGAT1 reprogram central metabolism to enhance oil accumulation, with relevance to mature seeds.

59 BASIC BIOLOGICAL SCIENCES↗

Patterns of gene expression in atrophying skeletal muscles: response to food deprivation

During fasting and many systemic diseases, muscle undergoes rapid loss of protein and functional capacity. To define the transcriptional changes triggering muscle atrophy and energy conservation in fasting, we used cDNA microarrays to compare mRNAs from muscles of control and food-deprived mice. Expression of >94% of genes did not change, but interesting patterns emerged among genes that were differentially expressed: 1) mRNAs encoding polyubiquitin, ubiquitin extension proteins, and many (but not all) proteasome subunits increased, which presumably contributes to accelerated protein breakdown; 2) a dramatic increase in mRNA for the ubiquitin ligase, atrogin-1, but not most E3s; 3) a significant suppression of mRNA for myosin binding protein H (but not other myofibrillar proteins) and IGF binding protein 5, which may favor cell protein loss; 4) decreases in mRNAs for several glycolytic enzymes and phosphorylase kinase subunits, and dramatic increases in mRNAs for pyruvate dehydrogenase kinase 4 and glutamine synthase, which should promote glucose sparing and gluconeogenesis. During fasting, metallothionein mRNA increased dramatically, mRNAs for extracellular matrix components fell, and mRNAs that may favor cap-independent mRNA translation rose. Significant changes occurred in mRNAs for many growth-related proteins and transcriptional regulators. These transcriptional changes indicate a complex adaptive program that should favor protein degradation and suppress glucose oxidation in muscle. Similar analysis of muscles atrophying for other causes is allowing us to identify a set of atrophy-specific changes in gene expression.

Non-NASA Center↗

Effects of spaceflight and simulated weightlessness on longitudinal bone growth

Indirect measurements have suggested that spaceflight impairs bone elongation in rats. To test this possibility, our laboratory measured, by the fluorochrome labeling technique, bone elongation that occurred during a spaceflight experiment. The longitudinal growth rate (LGR) in the tibia of rats in spaceflight experiments (Physiological Space Experiments 1, 3, and 4 and Physiological-Anatomical Rodent Experiment 3) and in two models of skeletal unloading (hind-limb elevation and unilateral sciatic neurotomy) were calculated. The effects of an 11 day spaceflight on gene expression of cartilage matrix proteins in rat growth plates were also determined by northern analysis and are reported for the first time in this study. Measurements of longitudinal growth indicate that skeletal unloading generally did not affect LGR, regardless of age, strain, gender, duration of unloading, or method of unloading. There was, however, one exception with 34% suppression in LGR detected in slow-growing, ovariectomized rats skeletally unloaded for 8 days by hind-limb elevation. This detection of reduced LGR by hind-limb elevation is consistent with changes in steady-state mRNA levels for type II collagen (-33%) and for aggrecan (-53%) that were detected in rats unloaded by an 11 day spaceflight. The changes detected in gene expression raise concern that spaceflight may result in changes in the composition of extracellular matrix, which could have a negative impact on conversion of growth-plate cartilage into normal cancellous bone by endochondral ossification.

manned↗

Multi-Omics Reveals Temporal Scales of Carbon Metabolism in Synechococcus Elongatus PCC 7942 Under Light Disturbance

Central carbon metabolism in model cyanobacteria involves multiple pathways to adapt to energy-light limitations across diel cycles. However, the success in mechanistic modeling for phenotypic prediction of the protein regulators in the metabolic state depends on capturing the vast possibilities emerging from multiple regulatory pathways in complex biological processes. Here, we developed a physics-informed machine learning approach based on energy-landscape concepts to predict regulatory proteins responding to cyclic circadian and unforeseen light perturbations in cyanobacterial metabolic networks. Our approach provides interpretable de novo models for inferring gene expression dynamics from Synechococcus elongatus over diel cycles and using redox proteome analysis to distinguish immediate light-responsive elements from circadian-regulated processes in carbon metabolism pathways. We identified distinct temporal signatures with the analysis of the redox proteome: there was an immediate shift in cysteine redox states accompanied by a limited change in protein abundance under constant illumination and after 2 hours of darkness. This discovery indicates that the generation of reductants coordinates photoinduced electron transport with redox metabolic pathways in two discernable molecular mechanisms: fast redox-based protein modifications occur immediately after the light disturbance, followed by slow transcriptional regulations across networks. This temporal regulation reveals how metabolic networks integrate rapid light responses with programmed circadian rhythms to maintain cellular homeostasis under the light-energy limitations over the diel cycle.

Biomolecular & subcellular processes↗

Molecular and biochemical evidence for the involvement of calcium/calmodulin in auxin action

The use of (35)S-labeled calmodulin (CaM) to screen a corn root cDNA expression library has led to the isolation of a CaM-binding protein, encoded by a cDNA with sequence similarity to small auxin up RNAs (SAURs), a class of early auxin-responsive genes. The cDNA designated as ZmSAUR1 (Zea mays SAURs) was expressed in Escherichia coli, and the recombinant protein was purified by CaM affinity chromatography. The CaM binding assay revealed that the recombinant protein binds to CaM in a calcium-dependent manner. Deletion analysis revealed that the CaM binding site was located at the NH(2)-terminal domain. A synthetic peptide of amino acids 20-45, corresponding to the potential CaM binding region, was used for calcium-dependent mobility shift assays. The synthetic peptide formed a stable complex with CaM only in the presence of calcium. The CaM affinity assay indicated that ZmSAUR1 binds to CaM with high affinity (K(d) approximately 15 nM) in a calcium-dependent manner. Comparison of the NH(2)-terminal portions of all of the characterized SAURs revealed that they all contain a stretch of the basic alpha-amphiphilic helix similar to the CaM binding region of ZmSAUR1. CaM binds to the two synthetic peptides from the NH(2)-terminal regions of Arabidopsis SAUR-AC1 and soybean 10A5, suggesting that this is a general phenomenon for all SAURs. Northern analysis was carried out using the total RNA isolated from auxin-treated corn coleoptile segments. ZmSAUR1 gene expression began within 10 min, increased rapidly between 10 and 60 min, and peaked around 60 min after 10 microM alpha-naphthaleneacetic acid treatment. These results indicate that ZmSAUR1 is an early auxin-responsive gene. The CaM antagonist N-(6-aminohexyl)5-chloro-1-naphthalenesulfonamide hydrochloride inhibited the auxin-induced cell elongation but not the auxin-induced expression of ZmSAUR1. This suggests that calcium/CaM do not regulate ZmSAUR1 at the transcriptional level. CaM binding to ZmSAUR1 in a calcium-dependent manner suggests that calcium/CaM regulate ZmSAUR1 at the post-translational level. Our data provide the first direct evidence for the involvement of calcium/CaM-mediated signaling in auxin-mediated signal transduction.

NASA Program Fundamental Space Biology↗

Microarray Data Analysis of Space Grown Arabidopsis Leaves for Genes Important in Vascular Patterning

Venation patterning in leaves is a major determinant of photosynthesis efficiency because of its dependency on vascular transport of photoassimilates, water, and minerals. Arabidopsis thaliana grown in microgravity show delayed growth and leaf maturation. Gene expression data from the roots, hypocotyl, and leaves of A. thaliana grown during spaceflight vs. ground control analyzed by Affymetrix microarray are available through NASAs GeneLab (GLDS-7). We analyzed the data for differential expression of genes in leaves resulting from the effects of spaceflight on vascular patterning. Two genes were found by preliminary analysis to be upregulated during spaceflight that may be related to vascular formation. The genes are responsible for coding an ARGOS like protein (potentially affecting cell elongation in the leaves), and an F-boxkelch-repeat protein (possibly contributing to protoxylem specification). Further analysis that will focus on raw data quality assessment and a moderated t-test may further confirm upregulation of the two genes and/or identify other gene candidates. Plants defective in these genes will then be assessed for phenotype by the mapping and quantification of leaf vascular patterning by NASAs VESsel GENeration (VESGEN) software to model specific vascular differences of plants grown in spaceflight.

Weitzeal, A. J.↗

Microarray Data Analysis of Space Grown Arabidopsis Leaves for Genes Important in Vascular Patterning

Venation patterning in leaves is a major determinant of photosynthesis efficiency because of its dependency on vascular transport of photoassimilates, water, and minerals. Arabidopsis thaliana grown in microgravity show delayed growth and leaf maturation. Gene expression data from the roots, hypocotyl, and leaves of A. thaliana grown during spaceflight vs. ground control analyzed by Affymetrix microarray are available through NASA's GeneLab (GLDS-7). We analyzed the data for differential expression of genes in leaves resulting from the effects of spaceflight on vascular patterning. Two genes were found by preliminary analysis to be upregulated during spaceflight that may be related to vascular formation. The genes are responsible for coding an ARGOS like protein (potentially affecting cell elongation in the leaves), and an F-box/kelch-repeat protein (possibly contributing to protoxylem specification). Further analysis that will focus on raw data quality assessment and a moderated t-test may further confirm upregulation of the two genes and/or identify other gene candidates. Plants defective in these genes will then be assessed for phenotype by the mapping and quantification of leaf vascular patterning by NASA's VESsel GENeration (VESGEN) software to model specific vascular differences of plants grown in spaceflight.

Weitzeal, A. J.↗

Comparative genomics analysis of drought response between obligate CAM and C 3 photosynthesis plants

Crassulacean acid metabolism (CAM) plants exhibit elevated drought and heat tolerance compared to C 3 and C 4 plants through an inverted pattern of day/night stomatal closure and opening for CO 2 assimilation. However, the molecular responses to water-deficit conditions remain unclear in obligate CAM species. In this study, we presented genome-wide transcription sequencing analysis using leaf samples of an obligate CAM species Kalanchoë fedtschenkoi under moderate and severe drought treatments at two-time points of dawn (2-h before the start of light period) and dusk (2-h before the dark period). Differentially expressed genes were identified in response to environmental drought stress and a whole genome wide co-expression network was created as well. We found that the expression of CAM-related genes was not regulated by drought stimuli in K. fedtschenkoi. Our comparative analysis revealed that CAM species (K. fedtschenkoi) and C 3 species (Arabidopsis thaliana, Populus deltoides ‘WV94’) share some common transcriptional changes in genes involved in multiple biological processes in response to drought stress, including ABA signaling and biosynthesis of secondary metabolites.

59 BASIC BIOLOGICAL SCIENCES↗

Providing biological context for GWAS results using eQTL regulatory and co‐expression networks in Populus

Summary Our study utilized genome‐wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa , a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co‐expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation of a Novel Salt Stress-Responsive Pathway Mediated by Arabidopsis DEAD-Box RNA Helicase Gene AtRH17 Using RNA-Seq Analysis

Previously, we reported that overexpression of AtRH17, an Arabidopsis DEAD-box RNA helicase gene, confers salt stress-tolerance via a pathway other than the well-known salt stress-responsive pathways. To decipher the salt stress-responsive pathway in AtRH17-overexpressing transgenic plants (OXs), we performed RNA-Sequencing and identified 397 differentially expressed genes between wild type (WT) and AtRH17 OXs. Among them, 286 genes were upregulated and 111 genes were downregulated in AtRH17 OXs relative to WT. Gene ontology annotation enrichment and KEGG pathway analysis showed that the 397 upregulated and downregulated genes are involved in various biological functions including secretion, signaling, detoxification, metabolic pathways, catabolic pathways, and biosynthesis of secondary metabolites as well as in stress responses. Genevestigator analysis of the upregulated genes showed that nine genes, namely, LEA4-5, GSTF6, DIN2/BGLU30, TSPO, GSTF7, LEA18, HAI1, ABR, and LTI30, were upregulated in Arabidopsis under salt, osmotic, and drought stress conditions. In particular, the expression levels of LEA4-5, TSPO, and ABR were higher in AtRH17 OXs than in WT under salt stress condition. Taken together, our results suggest that a high AtRH17 expression confers salt stress-tolerance through a novel salt stress-responsive pathway involving nine genes, other than the well-known ABA-dependent and ABA-independent pathways.

59 BASIC BIOLOGICAL SCIENCES↗

Differential Responses to Mechanostimulation in Embryonic Stem Cells Versus the Embryoid Body Model of Development Assessed at Single Cell RNA-Seq Resolution

Mechanicalforces generated by gravity have shaped life on Earth and impact gene expression and morphogenesis during early development. In contrast disuse canreduce normal mechanical loading, resulting in altered cell and tissue function. Although loading in adult mammals is known to promote increased cell proliferation and differentiation, little is known about how cells respondto this stimulusduring early development. In this study we sought to understand, with single cell RNA-sequencing resolution, how a 60-minute pulse of 50xg hypergravity-generated 5kPa hydrostatic pressure, influences transcriptomic regulation of developmental processes in the Embryoid Body (EB) model. Our study included both day-9 EBs and progenitor mouse embryonic stem cells (ESCs) with or without the hydrostatic pressurepulse. Single cell tSNE mapping shows limited transcriptome shifts in response to thispulse in either ESCs or EBs; this pulse,however, induces greater positional shifts in EB mapping compared to ESCs, indicating the influence of mechanotransduction is more pronounced in later states of cell commitment within the developmental program.We assessed ESCs and EBs for differentially expressed (DE) genes with hydrostatic pressurepulse and found approximately 1/3 DE genes were shared. However, gene ontology (GO) pathway analysis show that EBs have choreographed responses associated with upregulation ofpathways formulticellular development, mechanical signal transduction, and DNA damage repair. Cluster transcriptome analysis of the EBs showsmechanostimulationpromotes maintenance of transitory cell phenotypes in early development,including EB cluster co-expression of markers for progenitor, post-implant epiblast and primitive endoderm phenotypes versus expression exclusivity in the non-pulsed clusters. Pseudotime analysisidentified three branching cell types susceptible tohydrostatic pressureinduction of cell fate decisions. In summary, this study provides novel evidence that ESC maintenance and EB development can be regulated by mechanostimulation,and that stem cells committed to a differentiation program are more sensitive to force-induced changes to their transcriptome.

Cassandra Juran↗

Microarray Data Analysis of Space Grown Arabidopsis Leaves for Genes Important in Vascular Patterning

Venation patterning in leaves is a major determinant of photosynthesis efficiency because of its dependency on vascular transport of photo-assimilates, water, and minerals. Arabidopsis thaliana grown in microgravity show delayed growth and leaf maturation. Gene expression data from the roots, hypocotyl, and leaves of A. thaliana grown during spaceflight vs. ground control analyzed by Affymetrix microarray are available through NASA's GeneLab (GLDS-7). We analyzed the data for differential expression of genes in leaves resulting from the effects of spaceflight on vascular patterning. Two genes were found by preliminary analysis to be up-regulated during spaceflight that may be related to vascular formation. The genes are responsible for coding an ARGOS (Auxin-Regulated Gene Involved in Organ Size)-like protein (potentially affecting cell elongation in the leaves), and an F-box/kelch-repeat protein (possibly contributing to protoxylem specification). Further analysis that will focus on raw data quality assessment and a moderated t-test may further confirm up-regulation of the two genes and/or identify other gene candidates. Plants defective in these genes will then be assessed for phenotype by the mapping and quantification of leaf vascular patterning by NASA's VESsel GENeration (VESGEN) software to model specific vascular differences of plants grown in spaceflight.

venation↗

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

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

Prachi Kothiyal↗

Functional analysis of Pogostemon cablin farnesyl pyrophosphate synthase gene and its binding transcription factor PcWRKY44 in regulating biosynthesis of patchouli alcohol

Farnesyl pyrophosphate synthase (FPPS) plays an important role in the synthesis of plant secondary metabolites, but its function and molecular regulation mechanism remain unclear in Pogostemon cablin . In this study, the full-length cDNA of the FPP synthase gene from P. cablin ( PcFPPS ) was cloned and characterized. The expressions of PcFPPS are different among different tissues (highly in P. cablin flowers). Subcellular localization analysis in protoplasts indicated that PcFPPS was located in the cytoplasm. PcFPPS functionally complemented the lethal FPPS deletion mutation in yeast CC25. Transient overexpression of PcFPPS in P. cablin leaves accelerated terpene biosynthesis, with an ~47% increase in patchouli alcohol. Heterologous overexpression of PcFPPS in tobacco plants was achieved, and it was found that the FPP enzyme activity was significantly up-regulated in transgenic tobacco by ELISA analysis. In addition, more terpenoid metabolites, including stigmasterol, phytol, and neophytadiene were detected compared with control by GC-MS analysis. Furthermore, with dual-LUC assay and yeast one-hybrid screening, we found 220 bp promoter of PcFPPS can be bound by the nuclear-localized transcription factor PcWRKY44. Overexpression of PcWRKY44 in P. cablin upregulated the expression levels of PcFPPS and patchoulol synthase gene ( PcPTS ), and then promote the biosynthesis of patchouli alcohol. Taken together, these results strongly suggest the PcFPPS and its binding transcription factor PcWRKY44 play an essential role in regulating the biosynthesis of patchouli alcohol.

Wang, Xiaobing↗

Dynamics of Global Gene Expression and Regulatory Elements in Growing Brachypodium Root System

Root systems are dynamic and adaptable organs that play critical roles in plant development. However, how roots grow and accumulate biomass during plant life cycle and in relation to shoot growth phenology remains understudied. A comprehensive time-dependent root morphological analysis integrated with molecular signatures is then required to advance our understanding of root growth and development. Here we studied Brachypodium distachyon rooting process by monitoring root morphology, biomass production, and C/N ratios during developmental stages. To provide insight into gene regulation that accompanies root growth, we generated comprehensive transcript profiles of Brachypodium whole-root system at four developmental stages. Our data analysis revealed that multiple biological processes including trehalose metabolism and various families of transcription factors (TFs) were differentially expressed in root system during plant development. In particular, the AUX/IAA, ERFs, WRKY, NAC, and MADS TF family members were upregulated as plant entered the booting/heading stage, while ARFs and GRFs were downregulated suggesting these TF families as important factors involved in specific phases of rooting, and possibly in regulation of transition to plant reproductive stages. We identified several Brachypodium candidate root biomass-promoting genes and cis -regulatory elements for further functional validations and root growth improvements in grasses.

59 BASIC BIOLOGICAL SCIENCES↗

Deletion of Carboxypeptidase E in β-Cells Disrupts Proinsulin Processing but Does Not Lead to Spontaneous Development of Diabetes in Mice

Carboxypeptidase E (CPE) facilitates the conversion of prohormones into mature hormones and is highly expressed in multiple neuroendocrine tissues. Carriers of CPE mutations have elevated plasma proinsulin and develop severe obesity and hyperglycemia. We aimed to determine whether loss of Cpe in pancreatic β-cells disrupts proinsulin processing and accelerates development of diabetes and obesity in mice. Pancreatic β-cell–specific Cpe knockout mice (βCpeKO; Cpe fl/fl x Ins1 Cre/+ ) lack mature insulin granules and have elevated proinsulin in plasma; however, glucose-and KCl-stimulated insulin secretion in βCpeKO islets remained intact. High-fat diet–fed βCpeKO mice showed weight gain and glucose tolerance comparable with those of Wt littermates. Notably, β-cell area was increased in chow-fed βCpeKO mice and β-cell replication was elevated in βCpeKO islets. Transcriptomic analysis of βCpeKO β-cells revealed elevated glycolysis and Hif1α-target gene expression. On high glucose challenge, β-cells from βCpeKO mice showed reduced mitochondrial membrane potential, increased reactive oxygen species, reduced MafA, and elevated Aldh1a3 transcript levels. Following multiple low-dose streptozotocin injections, βCpeKO mice had accelerated development of hyperglycemia with reduced β-cell insulin and Glut2 expression. In conclusion, these findings suggest that Cpe and proper proinsulin processing are critical in maintaining β-cell function during the development of hyperglycemia.

59 BASIC BIOLOGICAL SCIENCES↗

Exploring drought-responsive crucial genes in Sorghum

Drought severely affects global food production. Sorghum is a typical drought-resistant model crop. Based on RNA-seq data for Sorghum with multiple time points and the gray correlation coefficient, this paper firstly selects candidate genes via mean variance test and constructs weighted gene differential co-expression networks (WGDCNs); then, based on guilt-by-rewiring principle, the WGDCNs and the hidden Markov random field model, drought-responsive crucial genes are identified for five developmental stages respectively. Enrichment and sequence alignment analysis reveal that the screened genes may play critical functional roles in drought responsiveness. A multilayer differential co-expression network for the screened genes reveals that Sorghum is very sensitive to pre-flowering drought. Furthermore, a crucial gene regulatory module is established, which regulates drought responsiveness via plant hormone signal transduction, MAPK cascades, and transcriptional regulations. The proposed method can well excavate crucial genes through RNA-seq data, which have implications in breeding of new varieties with improved drought tolerance.

60 APPLIED LIFE SCIENCES↗

Discovery of FoTO1 and Taxol genes enables biosynthesis of baccatin III

Abstract Plants make complex and potent therapeutic molecules 1,2 , but sourcing these molecules from natural producers or through chemical synthesis is difficult, which limits their use in the clinic. A prominent example is the anti-cancer therapeutic paclitaxel (sold under the brand name Taxol), which is derived from yew trees (Taxusspecies) 3 . Identifying the full paclitaxel biosynthetic pathway would enable heterologous production of the drug, but this has yet to be achieved despite half a century of research 4 . WithinTaxus’ large, enzyme-rich genome 5 , we suspected that the paclitaxel pathway would be difficult to resolve using conventional RNA-sequencing and co-expression analyses. Here, to improve the resolution of transcriptional analysis for pathway identification, we developed a strategy we term multiplexed perturbation × single nuclei (mpXsn) to transcriptionally profile cell states spanning tissues, cell types, developmental stages and elicitation conditions. Our data show that paclitaxel biosynthetic genes segregate into distinct expression modules that suggest consecutive subpathways. These modules resolved seven new genes, allowing a de novo 17-gene biosynthesis and isolation of baccatin III, the industrial precursor to Taxol, inNicotiana benthamianaleaves, at levels comparable with the natural abundance inTaxusneedles. Notably, we found that a nuclear transport factor 2 (NTF2)-like protein, FoTO1, is crucial for promoting the formation of the desired product during the first oxidation, resolving a long-standing bottleneck in paclitaxel pathway reconstitution. Together with a new β-phenylalanine-CoA ligase, the eight genes discovered here enable the de novo biosynthesis of 3’-N-debenzoyl-2’-deoxypaclitaxel. More broadly, we establish a generalizable approach to efficiently scale the power of co-expression analysis to match the complexity of large, uncharacterized genomes, facilitating the discovery of high-value gene sets.

Science & Technology - Other Topics↗