Transcriptomic Analysis Reveals the High-Oleic Acid Feedback Regulating the Homologous Gene Expression of Stearoyl-ACP Desaturase 2 (SAD2) in Peanuts

Genomic and co-expression network analyses reveal candidate genes for oil accumulation based on an introgression population in Upland cotton (Gossypium hirsutum)

KEY MESSAGE

Integrated QTL mapping and WGCNA condense the potential gene regulatory network involved in oil accumulation. A glycosyl hydrolases gene (GhHSD1) for oil biosynthesis was confirmed in Arabidopsis, which will provide useful knowledge to understand the functional mechanism of oil biosynthesis in cotton. Cotton is an economical source of edible oil for the food industry. The genetic mechanism that regulates oil biosynthesis in cottonseeds is essential for the genetic enhancement of oil content (OC). To explore the functional genomics of OC, this study utilized an interspecific backcross inbred line population to dissect the quantitative trait locus (QTL) interlinked with OC. In total, nine OC QTLs were identified, four of which were novel, and each QTL explained 3.62-34.73% of the phenotypic variation of OC. The comprehensive transcript profiling of developing cottonseeds revealed 3,646 core genes differentially expressed in both inbred parents. Functional enrichment analysis determined 43 genes were annotated with oil biosynthesis processes. Implementation of weighted gene co-expression network analysis showed that 803 differential genes had a significant correlation with the OC phenotype. Further integrated analysis identified seven important genes located in OC QTLs. Of which, the GhHSD1 gene located in stable QTL qOC-Dt3-1 exhibited the highest functional linkages with the other network genes. Phylogenetic analysis showed significant evolutionary differences in the HSD1 sequences between oilseed- and starch- crops. Furthermore, the overexpression of GhHSD1 in Arabidopsis yielded almost 6.78% higher seed oil. This study not only uncovers important genetic loci for oil accumulation in cottonseed, but also provides a set of new candidate genes that potentially influence the oil biosynthesis pathway in cottonseed.

Genome-wide comprehensive characterization and transcriptomic analysis of AP2/ERF gene family revealed its role in seed oil and ALA formation in perilla (Perilla frutescens)

Perilla (Perilla frutescens) is a potential specific oilseed crop with an extremely high α-linolenic acid (ALA) content in its seeds. AP2/ERF transcription factors (TFs) play important roles in multiple biological processes. However, limited information is known about the regulatory mechanism of the AP2/ERF family in perilla's oil accumulation. In this research, we identified 212 AP2/ERF family members in the genome of perilla, and their domain characteristics, collinearity, and sub-genome differentiation were comprehensively analyzed. Transcriptome sequencing revealed that genes encoding key enzymes involved in oil biosynthesis (e.g., ACCs, KASII, GPAT, PDAT and LPAAT) were up-regulated in the high-oil variety. Moreover, the endoplasmic reticulum-localized FAD2 and FAD3 were significantly up-regulated in the high-ALA variety. To investigate the roles of AP2/ERFs in lipid biosynthesis, we conducted a correlation analysis between non-redundant AP2/ERFs and key lipid metabolism genes using WGCNA. A significant correlation was found between 36 AP2/ERFs and 90 lipid metabolism genes. Among them, 12 AP2/ERFs were identified as hub genes and showed significant correlation with lipid synthase genes (e.g., FADs, GPAT and ACSL) and key regulatory TFs (e.g., LEC2, IAA, MYB, UPL3). Furthermore, gene expression analysis identified three AP2/ERFs (WRI, ABI4, and RAVI) potentially playing an important role in the regulation of oil accumulation in perilla. Our study suggests that PfAP2/ERFs are important regulatory TFs in the lipid biosynthesis pathway, providing a foundation for the molecular understanding of oil accumulation in perilla and other oilseed crops.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

scRNA-seq Reveals the Mechanism of Fatty Acid Desaturase 2 Mutation to Repress Leaf Growth in Peanut (Arachis hypogaea L.)

Fatty Acid Desaturase 2 (FAD2) controls the conversion of oleic acids into linoleic acids. Mutations in FAD2 not only increase the high-oleic content, but also repress the leaf growth. However, the mechanism by which FAD2 regulates the growth pathway has not been elucidated in peanut leaves with single-cell resolution. In this study, we isolated fad2 mutant leaf protoplast cells to perform single-cell RNA sequencing. Approximately 24,988 individual cells with 10,249 expressed genes were classified into five major cell types. A comparative analysis of 3495 differentially expressed genes (DEGs) in distinct cell types demonstrated that fad2 inhibited the expression of the cytokinin synthesis gene LOG in vascular cells, thereby repressing leaf growth. Further, pseudo-time trajectory analysis indicated that fad2 repressed leaf cell differentiation, and cell-cycle evidence displayed that fad2 perturbed the normal cell cycle to induce the majority of cells to drop into the S phase. Additionally, important transcription factors were filtered from the DEG profiles that connected the network involved in high-oleic acid accumulation (WRKY6), activated the hormone pathway (WRKY23, ERF109), and potentially regulated leaf growth (ERF6, MYB102, WRKY30). Collectively, our study describes different gene atlases in high-oleic and normal peanut seedling leaves, providing novel biological insights to elucidate the molecular mechanism of the high-oleic peanut-associated agronomic trait at the single-cell level.

Analysis of Delta(9) fatty acid desaturase gene family and their role in oleic acid accumulation in Carya cathayensis kernel

Carya cathayensis, commonly referred to as Chinese hickory, produces nuts that contain high-quality edible oils, particularly oleic acid (18:1). It is known that stearoyl-ACP desaturase (SAD) is the first key step converting stearic acid (C18:0, SA) to oleic acid (C18:1, OA) in the aminolevulinic acid (ALA) biosynthetic pathway and play an important role in OA accumulation. Thus far, there is little information about SAD gene family in C. cathayensis and the role of individual members in OA accumulation. This study searched the Chinese Hickory Genome Database and identified five members of SAD genes, designated as CcSADs, at the whole genome level through the comparison with the homologous genes from Arabidopsis. RNA-Seq analysis showed that CcSSI2-1, CcSSI2-2, and CcSAD6 were highly expressed in kernels. The expression pattern of CcSADs was significantly correlated with fatty acid accumulation during the kernel development. In addition, five full-length cDNAs encoding SADs were isolated from the developing kernel of C. cathayensis. CcSADs-green fluorescent protein (GFP) fusion construct was infiltrated into tobacco epidermal cells, and results indicated their chloroplast localization. The catalytic function of these CcSADs was further analyzed by heterologous expression in Saccharomyces cerevisiae, Nicotiana benthamiana, and walnut. Functional analysis demonstrated that all CcSADs had fatty acid desaturase activity to catalyze oleic acid biosynthesis. Some members of CcSADs also have strong substrate specificity for 16:0-ACP to synthesize palmitoleic acid (C16:1, PA). Our study documented SAD gene family in C. cathayensis and the role of CcSSI2-1, CcSSI2-2, and CcSAD6 in OA accumulation, which could be important for future improvement of OA content in this species via genetic manipulation.

Effects of Low Nighttime Temperature on Fatty Acid Content in Developing Seeds from Brassica napus L. Based on RNA-Seq and Metabolome

Brassica napus L. is a vital plant oil resource worldwide. The fatty acid biosynthesis and oil accumulation in its seeds are controlled by several genetic and environmental factors, including daytime and nighttime temperatures. We analyzed changes in oleic and erucic acid content in two double haploid (DH) lines, DH0729, a weakly temperature-sensitive line, and DH0815, a strongly temperature-sensitive line, derived from B. napus plants grown at different altitudes (1600, 1800, 2000, 2200, and 2400 m a.s.l., 28.85° N, 112.35° E) and nighttime temperatures (20/18, 20/16, 20/13 and 20/10 °C, daytime/nighttime temperature). Based on medium- and long-chain fatty acid metabolites, the total oleic acid content 35 and 43 days after flowering was significantly lower in low nighttime temperature (LNT, 20/13 °C) plants than in high nighttime temperature (HNT, 20/18 °C) plants (HNT: 58-62%; LNT: 49-54%; an average decrease of 9%), and the total erucic acid content was significantly lower in HNT than in LNT plants (HNT: 1-2%; LNT: 8-13%; an average increase of 10%). An RNA-seq analysis showed that the expression levels of SAD (LOC106366808), ECR (LOC106396280), KCS (LOC106419344), KAR (LOC106367337), HB1(LOC106430193), and DOF5 (LOC111211868) in STSL seeds increased under LNT conditions. In STSL seeds, a base mutation in the cis-acting element involved in low-temperature responsiveness (LTR), the HB1 and KCS promoter caused loss of sensitivity to low temperatures, whereas that of the KCS promoter caused increased sensitivity to low temperatures.

Genome-Wide Identification and Expression of FAR1 Gene Family Provide Insight Into Pod Development in Peanut (Arachis hypogaea)

The far-red-impaired response 1 (FAR1) transcription family were initially identified as important factors for phytochrome A (phyA)-mediated far-red light signaling in Arabidopsis; they play crucial roles in controlling the growth and development of plants. The reported reference genome sequences of Arachis, including A. duranensis, A. ipaensis, A. monticola, and A. hypogaea, and its related species Glycine max provide an opportunity to systematically perform a genome-wide identification of FAR1 homologous genes and investigate expression patterns of these members in peanut species. Here, a total of 650 FAR1 genes were identified from four Aarchis and its closely related species G. max. Of the studied species, A. hypogaea contained the most (246) AhFAR1 genes, which can be classified into three subgroups based on phylogenic relationships. The synonymous (Ks) and non-synonymous (Ka) substitution rates, phylogenetic relationship and synteny analysis of the FAR1 family provided deep insight into polyploidization, evolution and domestication of peanut AhFAR1 genes. The transcriptome data showed that the AhFAR1 genes exhibited distinct tissue- and stage-specific expression patterns in peanut. Three candidate genes including Ahy_A10g049543, Ahy_A06g026579, and Ahy_A10g048401, specifically expressed in peg and pod, might participate in pod development in the peanut. The quantitative real-time PCR (qRT-PCR) analyses confirmed that the three selected genes were highly and specifically expressed in the peg and pod. This study systematically analyzed gene structure, evolutionary characteristics and expression patterns of FAR1 gene family, which will provide a foundation for the study of genetic and biological function in the future.

Transcriptome analysis and identification of genes associated with oil accumulation in upland cotton

Cotton is not only the most important fiber crop but also the fifth most important oilseed crop in the world because of its oil-rich seeds as a byproduct of fiber production. By comparative transcriptome analysis between two germplasms with diverse oil accumulation, we reveal pieces of the gene expression network involved in the process of oil synthesis in cottonseeds. Approximately, 197.16 Gb of raw data from 30 RNA sequencing samples with 3 biological replicates were generated. Comparison of the high-oil and low-oil transcriptomes enabled the identification of 7682 differentially expressed genes (DEGs). Based on gene expression profiles relevant to triacylglycerol (TAG) biosynthesis, we proposed that the Kennedy pathway (diacylglycerol acyltransferase-catalyzed diacylglycerol to TAG) is the main pathway for oil production, rather than the phospholipid diacylglycerol acyltransferase-mediated pathway. Using weighted gene co-expression network analysis, 5312 DEGs were obtained and classified into 14 co-expression modules, including the MEblack module containing 10 genes involved in lipid metabolism. Among the DEGs in the MEblack module, GhCYSD1 was identified as a potential key player in oil biosynthesis. The overexpression of GhCYSD1 in yeast resulted in increased oil content and altered fatty acid composition. This study may not only shed more light on the underlying molecular mechanism of oil accumulation in cottonseed oil, but also provide a set of new gene for potential enhancement of oil content in cottonseeds.

A Comparison of Lipid Contents in Different Types of Peanut Cultivars Using UPLC-Q-TOF-MS-Based Lipidomic Study

Peanuts are a rich dietary source of lipids, which are essential for human health. In this study, the lipid contents of 13 peanut cultivars were analyzed using UPLC-Q-TOF-MS and GC–MS. The OXITEST reactor was used to test their lipid oxidation stabilities. A total of 27 subclasses, 229 individual lipids were detected. The combined analysis of lipid and oxidation stability showed that lipid unsaturation was inversely correlated with oxidation stability. Moreover, lipid profiles differed significantly among the different peanut cultivars. A total of 11 lipid molecules (TG 18:2/18:2/18:2, TG 24:0/18:2/18:3, TG 20:5/14:1/18:2, TG 18:2/14:1/18:2, PE 17:0/18:2, BisMePA 18:2/18:2, PG 38:5, PMe 18:1/18:1, PC 18:1/18:1, MGDG 18:1/18:1, TG 10:0/10:1/18:1) might be employed as possible indicators to identify high oleic acid (OA) and non-high OA peanut cultivars, based on the PLS-DA result of lipid molecules with a VIP value greater than 2. This comprehensive analysis will help in the rational selection and application of peanut cultivars.

Single-cell RNA-seq describes the transcriptome landscape and identifies critical transcription factors in the leaf blade of the allotetraploid peanut (Arachis hypogaea L.)

Single-cell RNA-seq (scRNA-seq) has been highlighted as a powerful tool for the description of human cell transcriptome, but the technology has not been broadly applied in plant cells. Herein, we describe the successful development of a robust protoplast cell isolation system in the peanut leaf. A total of 6,815 single cells were divided into eight cell clusters based on reported marker genes by applying scRNA-seq. Further, a pseudo-time analysis was used to describe the developmental trajectory and interaction network of transcription factors (TFs) of distinct cell types during leaf growth. The trajectory enabled re-investigation of the primordium-driven development processes of the mesophyll and epidermis. These results suggest that palisade cells likely differentiate into spongy cells, while the epidermal cells originated earlier than the primordium. Subsequently, the developed method integrated multiple technologies to efficiently validate the scRNA-seq result in a homogenous cell population. The expression levels of several TFs were strongly correlated with epidermal ontogeny in accordance with obtained scRNA-seq values. Additionally, peanut AHL23 (AT-HOOK MOTIF NUCLEAR LOCALIZED PROTEIN 23), which is localized in nucleus, promoted leaf growth when ectopically expressed in Arabidopsis by modulating the phytohormone pathway. Together, our study displays that application of scRNA-seq can provide new hypotheses regarding cell differentiation in the leaf blade of Arachis hypogaea. We believe that this approach will enable significant advances in the functional study of leaf blade cells in the allotetraploid peanut and other plant species.

Transcriptomic Characterization of Nitrate-Enhanced Stevioside Glycoside Synthesis in Stevia (Stevia rebaudiana) Bertoni

Nitrogen forms (nitrate (NO₃⁻) or ammonium (NH₄⁺)) are vital to plant growth and metabolism. In stevia (Stevia rebaudiana), it is important to assess whether nitrogen forms can influence the synthesis of the high-value terpene metabolites-steviol glycosides (SGs), together with the underlying mechanisms. Field and pot experiments were performed where stevia plants were fertilized with either NO₃⁻ or NH₄⁺ nutrition to the same level of nitrogen. Physiological measurements suggested that nitrogen forms had no significant impact on biomass and the total nitrogen content of stevia leaves, but NO₃⁻-enhanced leaf SGs contents. Transcriptomic analysis identified 397 genes that were differentially expressed (DEGs) between NO₃⁻ and NH₄⁺ treatments. Assessment of the DEGs highlighted the responses in secondary metabolism, particularly in terpenoid metabolism, to nitrogen forms. Further examinations of the expression patterns of SGs synthesis-related genes and potential transcription factors suggested that GGPPS and CPS genes, as well as the WRKY and MYB transcription factors, could be driving N form-regulated SG synthesis. We concluded that NO₃⁻, rather than NH₄⁺, can promote leaf SG synthesis via the NO₃⁻-MYB/WRKY-GGPPS/CPS module. Our study suggests that insights into the molecular mechanism of how SG synthesis can be affected by nitrogen forms.

Screening and Interaction Analysis Identify Genes Related to Anther Dehiscence in Solanum melongena L

Anther indehiscence is an important form of functional male sterility that can facilitate the production of hybrid seeds. However, the molecular mechanisms of anther indehiscence-based male sterility in eggplant (Solanum melongena L.) have not been thoroughly explored. We performed transcriptome sequencing and real-time quantitative reverse transcription-PCR (qRT-PCR) assays to compare the fertile line (F142) and male sterile line (S12) eggplant. We identified 2,670 differentially expressed genes (DEGs) between lines. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses identified 31 DEGs related to hormone biosynthesis. We, therefore, measured phytohormone contents, such as jasmonic acid (JA), auxin (IAA), gibberellin (GA), and abscisic acid (ABA) in S12 and F142. There were differences in IAA, GA₃, and ABA levels between S12 and F142, while JA levels were significantly lower in S12 than in F142. Five key genes in the JA signaling pathway were differentially expressed in S12 vs. F142. Of these, SmJAZ1 and SmJAR1 were significantly upregulated and SmDAD1, SmLOX, and SmCOI1 were downregulated in S12 vs. F142. Protein-protein interaction studies identified a direct interaction between SmDAD1 and SmLOX, while SmDAD1 failed to interact with SmJAR1, SmCOI1, and SmJAZ1. The data represent a valuable resource for further exploration of regulatory mechanisms underlying anther dehiscence in eggplant.

A global survey of the gene network and key genes for oil accumulation in cultivated tetraploid cottons

To enrich our knowledge about gene network of fatty acid biosynthesis in cottonseed, we conducted comparative transcriptome to reveal the differences in gene expression between Gossypium hirsutum and Gossypium barbadense during cottonseed development. The prolonged expression period and increased expression abundance of oil-related genes are the main reasons for producing high seed oil content (SOC) in G. barbadense, which manifested as the bias of homeologous gene expression in Dt-subgenome after 25 day postanthesis (DPA). The dynamic expression profile showed that SAD6 and FATA are more important for oil biosynthesis in G. barbadense than that in G. hirsutum. Three key transcription factors, WRI1, NF-YB6 and DPBF2, showed their elite roles in regulating seed oil in cotton. We observed that sequence variations in the promoter region of BCCP2 genes might contribute to its divergence in expression level between the two species. Based on the quantitative trait loci (QTL) information of the seed oil content and utilizing additional G. barbadense introgression lines (ILs), we propose 21 candidate genes on the basis of their differential expression level, of which the GbSWEET and the GbACBP6 showed the potential functional to improve the oil content. Taken together, studying the different expression of oil-related genes and their genetic regulation mechanisms between G. hirsutum and G. barbadense provide new insights to understanding the mechanism of fatty acid biosynthesis network and fatty acid genetic improving breeding in cotton.

Characterization of the stearoyl-ACP desaturase gene (PoSAD) from woody oil crop Paeonia ostii var. lishizhenii in oleic acid biosynthesis

Paeonia ostii var. lishizhenii has been approved as a woody oil crop with the outstanding characteristic of abundant α-linolenic acid (C18:3, ALA) in its seed oil. The stearoyl-ACP desaturase gene (SAD) regulates the first key step from stearic acid (C18:0, SA) to oleic acid (C18:1, OA) in the ALA biosynthetic pathway, but its functional characterization in P. ostii var. lishizhenii is absent to date. In this study, a key PoSAD gene (1719 bp in length) was acquired from endosperm of P. ostii var. lishizhenii by transcriptome sequencing analysis and the RACE (rapid-amplification of cDNA ends) method. Bioinformatic analysis of the PoSAD protein showed high homology (ranging from 90.4% to 94.4%) and similar physical and chemical properties to SAD from other higher plants, indicating that it encodes a putative stearoyl-ACP desaturase. Analysis of cis-acting elements found several endosperm tissue-specific motifs; i.e., one Prolamin box, thirteen DOFCOREs and one RY repeat in its promoter. The results of the qRT-PCR experiments verified that PoSAD was most highly expressed in developing endosperm at 59 days after pollination (53.7 times that in shoots) compared with that in roots (1.4 times), stems (2.5 times), leaves (3.1 times), petals (13.1 times) and stamens (46.0 times). Meanwhile, the fatty acid contents in P. ostii var. lishizhenii endosperm at seven growth stages were compared with variation in PoSAD expression. Heterologous expression of PoSAD significantly decreased SA and increased OA content, which effectively reduced the ratios of SA to OA in Saccharomyces cerevisiae and Arabidopsis thaliana. However, contents and ratios of palmitic acid (C16:0) and palmitoleic acid (C16:1) were stable in transgenic yeast, and palmitoleic acid remained absent in transgenic A. thaliana seeds. These results illustrate that PoSAD plays an essential role in endosperm development of P. ostii var. lishizhenii, strictly in catalysis of SA desaturation and OA biosynthesis but without functioning in PA desaturation. The results contribute to our understanding of the characterization of PoSAD in OA biosynthesis in P. ostii var. lishizhenii.

Global Transcriptome and Correlation Analysis Reveal Cultivar-Specific Molecular Signatures Associated with Fruit Development and Fatty Acid Determination in Camellia oleifera Abel

Background

Oil-tea Camellia is a very important edible oil plant widely distributed in southern China. Tea oil extracted from the oil-tea Camellia seeds is beneficial to health and is considered as a health edible oil. We attempt to identify genes related to fatty acid biosynthesis in an oil-tea Camellia seed kernel, generated a comprehensive transcriptome analysis of the seed kernel at different developmental stages, and explore optimal picking time of fruit. Material and Methods. A gas chromatography-mass spectrometer was used to detect the content of various fatty acids in samples. Transcriptome analysis was performed to detect gene dynamics and corresponding functions.

Results

Multiple phenotypic data were counted in detail, including the oil content, oleic acid content, linoleic acid content, linolenic acid content, fruit weight, fruit height, fruit diameter, single seed weight, seed length, and seed width in different developmental stages, which indicate that a majority of indicators increased with the development of oil-tea Camellia. The transcriptomics was conducted to perform a comprehensive and system-level view on dynamic gene expression networks for different developmental stages. Short Time-series Expression Miner (STEM) analysis of XL106 (the 6 time points) and XL210 (8 time points) was performed to screen related fatty acid (FA) gene set, from which 1041 candidate genes related to FA were selected in XL106 and 202 related genes were screened in XL210 based on GO and KEGG enrichment. Then, candidate genes and trait dataset were combined to conduct correlation analysis, and 10 genes were found to be strongly connected with several key traits.

Conclusions

The multiple phenotypic data revealed the dynamic law of changes during the picking stage. Transcriptomic analysis identified a large number of potential key regulatory factors that can control the oil content of dried kernels, oleic acid, linoleic acid, linolenic acid, fresh seed rate, and kernel-to-seed ratio, thereby providing a new insight into the molecular networks underlying the picking stage of oil-tea Camellia, which provides a theoretical basis for the optimal fruit picking point.

Global transcriptome analysis of subterranean pod and seed in peanut (Arachis hypogaea L.) unravels the complexity of fruit development under dark condition

Peanut pods develop underground, which is the most salient characteristic in peanut. However, its developmental transcriptome remains largely unknown. In the present study, we sequenced over one billion transcripts to explore the developmental transcriptome of peanut pod using Illumina sequencing. Moreover, we identified and quantified the abundances of 165,689 transcripts in seed and shell tissues along with a pod developmental gradient. The dynamic changes of differentially expressed transcripts (DETs) were described in seed and shell. Additionally, we found that photosynthetic genes were not only pronouncedly enriched in aerial pod, but also played roles in developing pod under dark condition. Genes functioning in photomorphogenesis showed distinct expression profiles along subterranean pod development. Clustering analysis unraveled a dynamic transcriptome, in which transcripts for DNA synthesis and cell division during pod expansion were transitioning to transcripts for cell expansion and storage activity during seed filling. Collectively, our study formed a transcriptional baseline for peanut fruit development under dark condition.

Integrated Analysis of Comparative Lipidomics and Proteomics Reveals the Dynamic Changes of Lipid Molecular Species in High-Oleic Acid Peanut Seed

Modern peanut contains fatty acid desaturase 2 (FAD2) mutation, which is capable of producing high oleic acid for human health. However, the dynamic changes of the lipidome regarding fad2 remain elusive in peanut seed. In the present study, 547 lipid features were identified in high- and normal-oleic peanut seeds by utilizing the mass spectrometric approach. The fad2-induced differently expressed lipids (DELs) were polarly distributed at early and maturation stages during high-oleic acid (OA) seed development. Subsequently, integration of previously published proteomic data and lipidomic data revealed that 21 proteins and 149 DELs were annotated into the triacylglycerol assembly map, of which nine enzymes and 31 lipid species shared similar variation tendencies. Additionally, the variation tendencies of 17 acyl fatty acids were described in a hypothetical biosynthetic pathway. Collectively, the understanding of the lipid composition correlated with fad2 established a foundation for future high-OA peanut breeding based on lipidomic data.