Search NASA⌕ Search

SEARCH · Search NASA

Results for “sequence 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Five key aspects of metaproteomics as a tool to understand functional interactions in host-associated microbiomes

Host-associated microbial communities (microbiomes) play critical roles in human, animal, and plant health and development. However, interactions between the host, members of the microbiome, and invading pathogens are in most cases still poorly understood. Such interactions are multidimensional and can alter the taxonomic composition and/or the functional metabolic activities of the microbiome in response to disease or treatment conditions. For example, after 2 days of antibiotic treatment, the mouse gut microbiome is altered and more susceptible to invasion by the pathogen Clostridioides difficile. Studies of these multidimensional interactions have been fueled by the ability to use high-throughput sequencing of phylogenetic marker genes to profile microbial community composition and shotgun metagenomics to profile functional potential. However, many protein-coding genes predicted from metagenomes are not necessarily expressed under a given condition, and thus, it is difficult to assess the activities and functional interactions in microbial communities based on DNA sequencing data alone. The physiological and pathological processes expressed in these communities under specific conditions are better reflected by the abundances of transcripts or proteins. In this Pearl, we provide a brief introduction to metaproteomics, which is a tool for the large-scale analysis of proteins in microbiomes that allows researchers to address a diversity of questions related to functions and interactions in microbiomes. The term “metaproteomics” was first used in 2004 for “the large-scale characterization of the entire protein complement of environmental microbiota at a given point in time”, and since then, a large array of metaproteomics approaches have been developed. Our objective in this Pearl is to highlight what we feel are 5 essential elements to be considered for a metaproteomics research campaign and to introduce nonexpert readers to the topic without going into too much technical detail.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

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

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

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The parallelism motifs of genomic data analysis

Genomic datasets are growing dramatically as the cost of sequencing continues to decline and small sequencing devices become available. Enormous community databases store and share these data with the research community, but some of these genomic data analysis problems require large-scale computational platforms to meet both the memory and computational requirements. These applications differ from scientific simulations that dominate the workload on high-end parallel systems today and place different requirements on programming support, software libraries and parallel architectural design. For example, they involve irregular communication patterns such as asynchronous updates to shared data structures. We consider several problems in high-performance genomics analysis, including alignment, profiling, clustering and assembly for both single genomes and metagenomes. We identify some of the common computational patterns or ‘motifs’ that help inform parallelization strategies and compare our motifs to some of the established lists, arguing that at least two key patterns, sorting and hashing, are missing. This article is part of a discussion meeting issue ‘Numerical algorithms for high-performance computational science’.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
DOE OSTI2019

Limited genetic diversity across pathogen populations responsible for the global emergence of boxwood blight identified using SSR s

Boxwood blight is an emerging disease of ornamental as well as native boxwood. The disease became widely established in Europe at the beginning of the 21st century, prior to its more recent discovery in North America and Asia. Two sister‐species of fungi cause the disease, Calonectria pseudonaviculata (Cps) and C. henricotiae (Che). Prior efforts to quantify intraspecific genetic polymorphisms of Cps and Che have yielded little information, limiting the ability to understand the evolution and migration of these pathogens. This study describes the development and implementation of simple sequence repeat ( SSR ) markers to analyse genetic diversity from a global collection of Cps and Che isolates, representing major blight outbreaks since the disease was first identified in the UK in the late 1990s. Analysis of the Cps CB 002 genome sequence identified 180 single copy SSR loci using stringent search criteria, 11 of which were polymorphic and used to screen a global sample of 306 isolates. Fourteen multilocus genotypes of Cps and two multilocus genotypes of Che were identified. Twelve of the 14 Cps genotypes differed from each other by a single allele. The most common Cps genotype was found on all continents where boxwood blight is confirmed. Based on measurement of linkage disequilibrium, Cps showed no evidence of sexual recombination. Further in silico analysis identified 1594 SSR s using relaxed SSR definition criteria. Comparison of these SSR ‐containing loci with Cps and Che genome sequences representing three different genotypes demonstrated that single nucleotide polymorphisms might serve as informative genetic markers for future studies.

LeBlanc, N.↗
DOE OSTI2020

Predicting the chemical space of fungal polyketides by phylogeny-based bioinformatics analysis of polyketide synthase-nonribosomal peptide synthetase and its modification enzymes

Fungal polyketide synthase (PKS)–nonribosomal peptide synthetase (NRPS) hybrids are key enzymes for synthesizing structurally diverse hybrid natural products (NPs) with characteristic biological activities. Predicting their chemical space is of particular importance in the field of natural product chemistry. However, the unexplored programming rule of the PKS module has prevented prediction of its chemical structure based on amino acid sequences. Here, we conducted a phylogenetic analysis of 884 PKS–NRPS hybrids and a modification enzyme analysis of the corresponding biosynthetic gene cluster, revealing a hidden relationship between its genealogy and core structures. This unexpected result allowed us to predict 18 biosynthetic gene cluster (BGC) groups producing known carbon skeletons (number of BGCs; 489) and 11 uncharacterized BGC groups (171). The limited number of carbon skeletons suggests that fungi tend to select PK skeletons for survival during their evolution. The possible involvement of a horizontal gene transfer event leading to the diverse distribution of PKS–NRPS genes among fungal species is also proposed. This study provides insight into the chemical space of fungal PKs and the distribution of their biosynthetic gene clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

VRC34-Antibody Lineage Development Reveals How a Required Rare Mutation Shapes the Maturation of a Broad HIV-Neutralizing Lineage

Rare mutations have been proposed to restrict the development of broadly neutralizing antibodies against HIV-1, but this has not been explicitly demonstrated. We hypothesized that such rare mutations might be identified by comparing broadly neutralizing and non-broadly neutralizing branches of an antibody-developmental tree. Because sequences of antibodies isolated from the fusion peptide (FP)-targeting VRC34-antibody lineage suggested it might be suitable for such rare mutation analysis, we carried out next-generation sequencing (NGS) on B cell transcripts from donor N123, the source of the VRC34 lineage, and functionally and structurally characterized inferred intermediates along broadly neutralizing and poorly neutralizing developmental branches. The broadly neutralizing VRC34.01 branch required the rare heavy-chain mutation Y33P to bind FP, whereas the early bifurcated VRC34.05 branch did not require this rare mutation and evolved less breadth. Our results demonstrate how a required rare mutation can restrict development and shape the maturation of a broad HIV-1-neutralizing antibody lineage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Transcriptome sequencing of three Pseudo-nitzschia species reveals comparable gene sets and the presence of Nitric Oxide Synthase genes in diatoms

Diatoms are among the most diverse eukaryotic microorganisms on Earth, they are responsible for a large fraction of primary production in the oceans and can be found in different habitats. Pseudo-nitzschia are marine planktonic diatoms responsible for blooms in coastal and oceanic waters. We analyzed the transcriptome of three species, Pseudo-nitzschia arenysensis, Pseudo-nitzschia delicatissima and Pseudo-nitzschia multistriata, with different levels of genetic relatedness. These species have a worldwide distribution and the last one produces the neurotoxin domoic acid. We were able to annotate about 80% of the sequences in each transcriptome and the analysis of the relative functional annotations allowed comparison of the main metabolic pathways, pathways involved in the biosynthesis of isoprenoids (MAV and MEP pathways) and pathways putatively involved in domoic acid synthesis. The search for homologous transcripts among the target species and other congeneric species resulted in the discovery of a sequence annotated as Nitric Oxide Synthase (NOS), found uniquely in Pseudo-nitzschia multistriata. The predicted protein product contained all the domains of the canonical metazoan sequence. Putative NOS sequences were found in other available diatom datasets, supporting a role for nitric oxide as signaling molecule in this group of microalgae.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Genomic analysis of six new Geobacillus strains reveals highly conserved carbohydrate degradation architectures and strategies

In this work we report the whole genome sequences of six new Geobacillus xylanolytic strains along with the genomic analysis of their capability to degrade carbohydrates. The six sequenced Geobacillus strains described here have a range of GC contents from 43.9% to 52.5% and clade with named Geobacillus species throughout the entire genus. We have identified a ~200 kb unique super-cluster in all six strains, containing five to eight distinct carbohydrate degradation clusters in a single genomic region, a feature not seen in other genera. The Geobacillus strains rely on a small number of secreted enzymes located within distinct clusters for carbohydrate utilization, in contrast to most biomass-degrading organisms which contain numerous secreted enzymes located randomly throughout the genomes. All six strains are able to utilize fructose, arabinose, xylose, mannitol, gluconate, xylan, and α-1,6-glucosides. The gene clusters for utilization of these seven substrates have identical organization and the individual proteins have a high percent identity to their homologs. The strains show significant differences in their ability to utilize inositol, sucrose, lactose, α-mannosides, α-1,4-glucosides and arabinan.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Development of peptoid-based ligands for the removal of cadmium from biological media

Cadmium poisoning poses a serious health concern due to cadmium's increasing industrial use, yet there is currently no recommended treatment. The selective coordination of cadmium in a biological environment—i.e. in the presence of serum ions, small molecules, and proteins—is a difficult task. To address this challenge, a combinatorial library of peptoid-based ligands has been evaluated to identify structures that selectively bind to cadmium in human serum with minimal chelation of essential metal ions. Eighteen unique ligands were identified in this screening procedure, and the binding affinity of each was measured using metal titrations monitored by UV-vis spectroscopy. To evaluate the significance of each chelating moiety, sequence rearrangements and substitutions were examined. Analysis of a metal–ligand complex by NMR spectroscopy highlighted the importance of particular residues. Depletion experiments were performed in serum mimetics and human serum with exogenously added cadmium. These depletion experiments were used to compare and demonstrate the ability of these peptoids to remove cadmium from blood-like mixtures. In one of these depletion experiments, the peptoid sequence was able to deplete the cadmium to a level comparable to the reported acute toxicity limit. Evaluation of the metal selectivity in buffered solution and in human serum was performed to verify minimal off-target binding. These studies highlight a screening platform for the identification of metal–ligands that are capable of binding in a complex environment. They additionally demonstrate the potential utility of biologically-compatible ligands for the treatment of heavy metal poisoning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Cacao pod transcriptome profiling of seven genotypes identifies features associated with post-penetration resistance to Phytophthora palmivora

Abstract The oomycete Phytophthora palmivora infects the fruit of cacao trees ( Theobroma cacao ) causing black pod rot and reducing yields. Cacao genotypes vary in their resistance levels to P. palmivora , yet our understanding of how cacao fruit respond to the pathogen at the molecular level during disease establishment is limited. To address this issue, disease development and RNA-Seq studies were conducted on pods of seven cacao genotypes (ICS1, WFT, Gu133, Spa9, CCN51, Sca6 and Pound7) to better understand their reactions to the post-penetration stage of P. palmivora infection. The pod tissue- P. palmivora pathogen assay resulted in the genotypes being classified as susceptible (ICS1, WFT, Gu133 and Spa9) or resistant (CCN51, Sca6 and Pound7). The number of differentially expressed genes (DEGs) ranged from 1625 to 6957 depending on genotype. A custom gene correlation approach identified 34 correlation groups. De novo motif analysis was conducted on upstream promoter sequences of differentially expressed genes, identifying 76 novel motifs, 31 of which were over-represented in the upstream sequences of correlation groups and associated with gene ontology terms related to oxidative stress response, defense against fungal pathogens, general metabolism and cell function. Genes in one correlation group (Group 6) were strongly induced in all genotypes and enriched in genes annotated with defense-responsive terms. Expression pattern profiling revealed that genes in Group 6 were induced to higher levels in the resistant genotypes. An additional analysis allowed the identification of 17 candidate cis -regulatory modules likely to be involved in cacao defense against P. palmivora . This study is a comprehensive exploration of the cacao pod transcriptional response to P. palmivora spread after infection. We identified cacao genes, promoter motifs, and promoter motif combinations associated with post-penetration resistance to P. palmivora in cacao pods and provide this information as a resource to support future and ongoing efforts to breed P. palmivora -resistant cacao.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Establishing Data Analysis Pipeline for Bulk ATAC-Seq Datasets

We developed an analysis pipeline for transposase-accessible chromatin sequencing (ATAC-Seq) data derived from bulk samples, which brings together publicly available R packages in addition to command-line tools designed for analysis of bulk ATAC-Seq data and can be run on any computer running a Linux-like operating system such as Ubuntu or Apple OSX.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Emerging Trends and Technologies Used for the Identification, Detection, and Characterisation of Plant-Parasitic Nematode Infestation in Crops

Accurate identification and estimation of the population densities of microscopic, soil-dwelling plant-parasitic nematodes (PPNs) are essential, as PPNs cause significant economic losses in agricultural production systems worldwide. This study presents a comprehensive review of emerging techniques used for the identification of PPNs, including morphological identification, molecular diagnostics such as polymerase chain reaction (PCR), high-throughput sequencing, meta barcoding, remote sensing, hyperspectral analysis, and image processing. Classical morphological methods require a microscope and nematode taxonomist to identify species, which is laborious and time-consuming. Alternatively, quantitative polymerase chain reaction (qPCR) has emerged as a reliable and efficient approach for PPN identification and quantification; however, the cost associated with the reagents, instrumentation, and careful optimisation of reaction conditions can be prohibitive. High-throughput sequencing and meta-barcoding are used to study the biodiversity of all tropical groups of nematodes, not just PPNs, and are useful for describing changes in soil ecology. Convolutional neural network (CNN) methods are necessary to automate the detection and counting of PPNs from microscopic images, including complex cases like tangled nematodes. Remote sensing and hyperspectral methods offer non-invasive approaches to estimate nematode infestations and facilitate early diagnosis of plant stress caused by nematodes and rapid management of PPNs. This review provides a valuable resource for researchers, practitioners, and policymakers involved in nematology and plant protection. It highlights the importance of fast, efficient, and robust identification protocols and decision-support tools in mitigating the impact of PPNs on global agriculture and food security.

Plant Sciences↗
DOE OSTI2022

Microbial community diversity changes during voltage reversal repair in a 12-unit microbial fuel cell

Microbial fuel cell stacks (MFC-Stack) are often confronted with voltage reversals, likely due to an interplay between microbial community dynamics and insufficient electric circuit balancing. Herein, we provide new insight into voltage reversals by examining the microbiomes of twelve MFC units of a 12-liter Pilot-MFC-Stack during repair. Different biofilm repair methods (self-healing, electrostimulation, and re-acclimatization upon cross-inoculation) were used to evaluate the microbial community response. In addition, MFC-Stack simulation was performed based on Kirchhoff’s Second Law to predict values for source potentials and post-evaluate internal resistances. Analysis of the 16S rRNA amplicon sequencing data suggests that the biofilm repair methods could slowly heal damaged biofilms. Notably, severely voltage reversed MFC units had low electrogen relative abundances (18%) and positive anode potentials, while strong bioanodes and contained more than 50% electrogens and had negative anode potentials. Between-community analyses (beta diversity ordination and multinomial regression) of the voltage reversed MFC units revealed differences among biofilms in contrast to healthy/strong MFC units. Permutational multivariate analysis of variance (PERMANOVA) confirmed that reversed biofilms were, indeed, significantly (p < 0.05) different from stronger ones. Overall, these analyses demonstrated the utility of combining electrotechnical and microbial community analyses, especially beta diversity ordination and multinomial regression, to understand problematic MFC units and the potential success of a biofilm repair method. Finally, thicker biofilms were usually healthier and stronger, although thickness was no guarantee for proper structure and power function as all factors were interdependent. There was an evolutionary trend that strong anodes became stronger/healthier and others weaker. This spontaneous trend has to be considered to avoid irreversible voltage reversals and to repair electrogenic biofilms in an MFC-Stack.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Large‐scale insertional mutagenesis of Chlamydomonas supports phylogenomic functional prediction of photosynthetic genes and analysis of classical acetate‐requiring mutants

Summary Chlamydomonas reinhardtii is a unicellular green alga that is a key model organism in the study of photosynthesis and oxidative stress. Here we describe the large‐scale generation of a population of insertional mutants that have been screened for phenotypes related to photosynthesis and the isolation of 459 flanking sequence tags from 439 mutants. Recent phylogenomic analysis has identified a core set of genes, named GreenCut2, that are conserved in green algae and plants. Many of these genes are likely to be central to the process of photosynthesis, and they are over‐represented by sixfold among the screened insertional mutants, with insertion events isolated in or adjacent to 68 of 597 GreenCut2 genes. This enrichment thus provides experimental support for functional assignments based on previous bioinformatic analysis. To illustrate one of the uses of the population, a candidate gene approach based on genome position of the flanking sequence of the insertional mutant CAL 027_01_20 was used to identify the molecular basis of the classical C. reinhardtii mutation ac17 . These mutations were shown to affect the gene PDH 2 , which encodes a subunit of the plastid pyruvate dehydrogenase complex. The mutants and associated flanking sequence data described here are publicly available to the research community, and they represent one of the largest phenotyped collections of algal insertional mutants to date.

Dent, Rachel M.↗
DOE OSTI2017

Optimization of SAMtools sorting using OpenMP tasks

SAMtools is a widely-used genomics application for post-processing high-throughput sequence alignment data. Such sequence alignment data are commonly sorted to make downstream analysis more efficient. However, this sorting process itself can be computationally- and I/O-intensive: high-throughput sequence alignment files in the de facto standard binary alignment/map (BAM) format can be many gigabytes in size, and may need to be decompressed before sorting and compressed afterwards. As a result, BAM-file sorting can be a bottleneck in genomics workflows. Here, this paper describes a case study on the performance analysis and optimization of SAMtools for sorting large BAM files. OpenMP task parallelism and memory optimization techniques resulted in a speedup of 5.9X versus the upstream SAMtools 1.3.1 for an internal (in-memory) sort of 24.6 GiB of compressed BAM data (102.6 GiB uncompressed) with 32 processor cores, while a 1.98X speedup was achieved for an external (out-of-core) sort of a 271.4 GiB BAM file.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Tissue-specific differences in metabolites and transcripts contribute to the heterogeneity of ricinoleic acid accumulation in Ricinus communis L. (castor) seeds

Castor ( Ricinus communis L.) seeds are valued for their production of oils which can comprise up to 90% hydroxy-fatty acids (ricinoleic acid). Castor oil contains mono-, di- and tri- ricinoleic acid containing triacylglycerols (TAGs). Although the enzymatic synthesis of ricinoleic acid is well described, the differential compartmentalization of these TAG molecular species has remained undefined. To examine the distribution of hydroxy fatty acid accumulation within the endosperm and embryo tissues of castor seeds. Matrix assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) was used to map the distribution of triacylglycerols in tissue sections of castor seeds. In addition, the endosperm and embryo (cotyledons and embryonic axis) tissues were dissected and extracted for quantitative lipidomics analysis and Illumina-based RNA deep sequencing. Results: This study revealed an unexpected heterogeneous tissue distribution of mono-, di- and tri- hydroxy-triacylglycerols in the embryo and endosperm tissues of castor seeds. Pathway analysis based on transcript abundance suggested that distinct embryo- and endosperm-specific mechanisms may exist for the shuttling of ricinoleic acid away from phosphatidylcholine (PC) and into hydroxy TAG production. The embryo-biased mechanism appears to favor removal of ricinoleic acid from PC through phophatidylcholine: diacylglycerol acylytransferase (PDAT) while the endosperm pathway appears to remove ricinoleic acid from the PC pool by preferences of phospholipase A (PLA 2 α) and/or phosphatidylcholine: diacylglycerol cholinepshphotransferase PDCT. In conclusion, collectively, a combination of lipidomics and transcriptomics analyses revealed previously undefined spatial aspects of hydroxy fatty acid metabolism in castor seeds. These studies underscore a need for tissue-specific studies as a means to better understand the regulation of triacylglycerol accumulation in oilseeds.

59 BASIC BIOLOGICAL SCIENCES↗