Next Article in Journal
Considering Predictive Factors in the Diagnosis of Clinically Significant Prostate Cancer in Patients with PI-RADS 3 Lesions
Previous Article in Journal
Differential Regulation of Myocardial E3 Ligases and Deubiquitinases in Ischemic Heart Failure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts

1
Key Laboratory for Forest Resources Conservation and Utilization in the Southwest Mountains of China, Ministry of Education, Southwest Landscape Architecture Engineering Research Center of National Forestry and Grassland Administration, Southwest Forestry University, Kunming 650224, China
2
Lincang Academy of Forestry, Lincang 677009, China
3
Yuxi Sannong Plateau Characteristic Modern Agriculture Co., Ltd., Chengjiang 652599, China
4
Guangxi South Subtropical Agricultural Science Research Institute, Longzhou 532415, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Life 2021, 11(12), 1431; https://doi.org/10.3390/life11121431
Submission received: 10 November 2021 / Revised: 7 December 2021 / Accepted: 15 December 2021 / Published: 18 December 2021
(This article belongs to the Section Plant Science)

Abstract

:
Macadamia nuts are considered a high-quality oil crop worldwide. To date, the lipid diversity and the genetic factors that mediate storage lipid biosynthesis in Macadamia ternifolia are poorly known. Here, we performed a comprehensive transcriptomic and lipidomic data analysis to understand the mechanism of lipid biosynthesis by using young, medium-aged, and mature fruit kernels. Our lipidomic analysis showed that the M. ternifolia kernel was a rich source of unsaturated fatty acids. Moreover, different species of triacylglycerols, diacylglycerol, ceramides, phosphatidylethanolamine, and phosphatidic acid had altered accumulations during the developmental stages. The transcriptome analysis revealed a large percentage of differently expressed genes during the different stages of macadamia growth. Most of the genes with significant differential expression performed functional activity of oxidoreductase and were enriched in the secondary metabolite pathway. The integration of lipidomic and transcriptomic data allowed for the identification of glycerol-3-phosphate acyltransferase, diacylglycerol kinase, phosphatidylinositols, nonspecific phospholipase C, pyruvate kinase 2, 3-ketoacyl-acyl carrier protein reductase, and linoleate 9S-lipoxygenase as putative candidate genes involved in lipid biosynthesis, storage, and oil quality. Our study found comprehensive datasets of lipidomic and transcriptomic changes in the developing kernel of M. ternifolia. In addition, the identification of candidate genes provides essential prerequisites to understand the molecular mechanism of lipid biosynthesis in the kernel of M. ternifolia.

1. Introduction

Macadamia is a nut tree with high-quality kernel oil. It belongs to the family Proteaceae, known as Australian or sometimes Hawaiian nuts, and has been grown worldwide in tropical and subtropical regions [1,2]. The genus Macadamia has four types of species that are commercially grown to produce nuts. These species are Macadamia integrifolia, M. tetraphylla, M. ternifolia, and M. jansenii [3]. However, M. integrifolia and M. tetraphylla have more significance due to their high-quality oils [4]. Approximately 44,000 metric tons of macadamia kernels are produced per year in the world and 14,100 metric tons per year are produced alone in Australia [5]. Because of its simple production technology, survival in low-fertility soil, high profitability, and low-temperature resistance compared to other traditional crops of tropical regions [6], the planting area of macadamia promptly increased over the past few years worldwide. China is a world-leading country in terms of planting area and is expected to reach more than 15,000 ha in the southern regions. About 90% of China’s total production comes from just Yunnan and Guangxi province [7]. The research knowledge related to breeding tools, nutrition management, disease management, processing methods, and product development is fundamental for the rapid devolvement of the macadamia oil industry around the world [8].
The macadamia is an evergreen nut tree and produces fruits 5–6 years after planting [9]; the fruits are commonly called nuts, and consist of a husk, shell, and kernel [10]. The kernel is a quality product, the edible part, and contains more than 60% lipids [5]. The kernel is consumed fresh, fried, roasted, or caramelized [11]. Interestingly, macadamia nuts are comprised of a high amount of unsaturated fats (70 to 80%) as compared with other edible nuts, specifically oleic and palmitoleic acids [12,13]. The high rate of palmitoleic acids can be utilized as biofuel to reduce issues of energy availability [14]. However, fuel properties, such as the cetane number, kinematic viscosity, oxidative stability, cold flow, and lubricity, in M. ternifolia still need exploration. The high percentage of unsaturated fats maintained lower blood cholesterol and ultimately prevent chronic diseases, e.g., cardiovascular diseases, cancer, and dyslipidemia in humans [15,16,17]. Because of these health benefits, macadamia nuts, together with walnuts, peanuts, and almonds, are named as heart-healthy foods by the Food and Drug Administration [18]. Moreover, macadamia nuts are a vital source of minerals (calcium, iron, phosphorus, magnesium, and potassium) [4], vitamins (thiamine, riboflavin, retinol, and niacin), carbohydrates, protein, and fiber [13]. However, the nut nutrient characteristics generally depend on the species, cultivar, growing conditions, location, and management practices [19]. The refined oil obtained from macadamia nut is pale yellow and is widely utilized in the cake and cookies industry [20]. Besides this, macadamia oil is used in the pharmaceutical industry but also has potential applications in the cosmetics industry to produce skin-cleaning and -lightening agents [21].
In terms of various interests, plant lipids are the primary source of energy and essential fatty acids and they have useful industry applications and a significant impact on the economy and human diet around the world [22]. Natural fatty acids are divided into saturated or unsaturated forms. Most saturated fatty acids have a hydrocarbon chain with a single bond, high melting point, and are solid at room temperature. In contrast, unsaturated fatty acids have a hydrocarbon chain with multiple bonds, low melting points, and are liquid at room temperature [23]. In recent years, high-throughput sequencing has become a routine task to identify putative candidate genes for a specific trait or character [24,25]. The lipids profile is a major determinant of oil quality during storage and processing. In this regard, high-quality RNA-seq research was utilized to understand the genetic mechanisms of lipid regulation in Arabidopsis [26], soybean [27], Brassica napus [28], jatropha [29], peanut [30], peach [31], castor bean [32], and pecan [33]. Research in oilseed crops identified that genes associated with 3-ketoacyl-acyl carrier protein synthase (KAS), beta-ketoacyl reductase (KAR), wrinkledI1 (WRI1), glycerol-3-phosphate acyltransferase (GPAT), lysophosphatidic acid acyltransferas (LPAT), diacylglycerol acyltransferase (DGAT), long-chain acyl-CoA synthetase (LACS), phosphatidic acid phosphatase (PAP), enoyl-acp reductase (ENR), and fatty acid desaturase (FADs) enzymes are the key regulators of plant lipids [34,35,36,37,38]. To date, the genomes of M. integrifolia [39] and M. tetraphylla [1] have been sequenced. However, the research progress has been very slow in M. ternifolia. The identification of genes involved in lipid biosynthesis and characterization of their expression patterns are two critical prerequisites to reveal the genetic mechanism of lipid biosynthesis in the M. ternifolia fruit kernel. Therefore, we designed this study regarding developing the fruit kernels of M. ternifolia with the following aims: (1) examine the diversity of lipids content, (2) identify key genes involved in the regulation of lipids, and (3) characterize their expression patterns. This study provides comprehensive insights into lipid composition and associated transcripts in developing fruit kernels. Our data resources are an important foundation to clarify the genetic mechanism of storage lipid accumulation in M. ternifolia.

2. Results

2.1. RNA Sequencing, De Novo Assembly, and Annotations

To obtain the genome-wide transcriptional changes in the fruit kernel of M. ternifolia, kernels from three different fruit developmental stages (young (S1-1,1-2,1-3), medium-aged (S2-1,2-2,2-3), and mature (S3-1,3-2,3-3) fruit) were selected to perform RNA sequencing. The mean raw reads in S1 were 54,363,432, whereas an average of 48,083,672 and 52,648,740 raw reads were obtained in the S2 and S3 stages, respectively (Table 1). The filtering of low-quality reads produced an average of 52,621,607, 45,606,890, and 50,793,072 clean reads in the libraries of S1, S2, and S3, respectively. A total of 7.9 Gb of clean data was obtained from the S1 stage, while 6.8 and 7.6 Gb were acquired from stages S2 and S3, respectively. The Q30% was 93% and the GC content percentage was almost 45 in our sequencing. This result suggested that we had produced high-quality data that was suitable for performing transcriptome downstream analysis. Trinity software was used to assemble the 307,336,019 total bases. The total number of unigenes assembled was 279,972. The average length was 1098 nt, with N50 and N90 length values of 1633 and 432 nt, respectively (Supplementary Materials Table S1).
The length distribution range of the unigenes was 200 to 2000 nt (Figure 1A). Due to the lack of a whole reference genome of M. ternifolia, all unigenes were blasted by using Nr, Swiss-Prot, KEGG, KOG, GO, and Trembl public databases to obtain annotations (Figure 1B). The majority of unigenes showed functional annotation in the Nr (163,325) and Trembl (162,849) databases. There were 119,590, 109,208, 96,830, 133,334, and 110,274 unigenes that were annotated in KEGG, SwissProt, KOG, GO, and Pfam, respectively. The abundance of each unigene was measured with fragments per kilobase per million reads (FPKM) values. The experiment and sampling suitability was determined with principal component analysis (PCA) and correlation matrixes. The PCA showed that the S1, S2, and S3 fruit kernel stages had significant variations with three distinct clusters (Figure 2A). However, a weak variation among replicates of each stage was observed. The difference was equally covered by PC1 and PC2, which explained 22.9 and 22.6% of the total variation, respectively. Moreover, a strong correlation (R2 > 0.90) was observed between the replicates of each stage (Figure 2B). In contrast, a weaker correlation was determined between different fruit kernel stages. The correlation value between S1 and S2 was around 0.50, whereas S1 and S3 showed a much weaker correlation (R2 = 0.15). The lower degree of correlation between the three stages indicated strong transcriptional changes.

2.2. Transcriptome Changes in Different Fruit Kernel Stages

To explore the gene expression changes in developing fruit kernels, we compared the young (S1), medium-aged (S2), and mature (S3) fruit kernel stages. Differentially expressed genes (DEGs) were obtained for two compared samples using p ≤ 0.05 and fold changes (ratio ≥ 1 or ratio ≤ −1) as thresholds (Figure 3A). The comparison between S1 and S2 resulted in 12,128 downregulated and 12,672 upregulated DEGs. A total of 36,135 DEGs were determined in the comparison of S1 and S3. From these, 18,546 genes were downregulated and 17,589 genes were upregulated. When we compared S2 and S3, 14,177 genes were downregulated and 11,421 genes were upregulated, with a total of 25,598 DEGs. The distribution of DEGs revealed that a major portion of the DEGs was overlapped but less stage-specific DEGs were identified among the different fruit kernel stages (Figure 3B). For example, the number of overlapped DEGs in S1 vs. S2 and S1 vs. S3 was 13,185. In addition, 3316 DEGs were overlapped among the S1, S2, and S3 fruit kernel stages. Globally, these results indicate major gene expression dynamics for fruit kernel development in M. ternifolia.

2.3. Functional Enrichment Analysis of DEGs

To obtain the complete functional information about transcripts that regulate fruit kernel development, GO and KEGG enrichment analyses were performed by using p ≤ 0.05 as the criterion for significant enrichment. In the GO analysis, the majority of DEGs in the S1 and S2 stages were involved in the molecular function of beta-glucosidase and oxidoreductase activity (Figure 4A). Meanwhile, in pathway enrichment, metabolic and biosynthesis of secondary metabolites had more significance and total genes. A total of 1595 DEGs were enriched in each of these pathways (Figure 4B). High numbers of oxidoreductase- and peroxidase-activity-associated genes were observed in the DEGs from the comparison of S1 vs. S3 (Figure 5A). DEGs with annotations regarding metabolism and biosynthesis of secondary metabolites exhibited predominant pathway enrichment (Figure 5B). Functional enrichment analysis of the S2 vs. S3 stages revealed that the majority of DEGs performed the biological function of cellulose biosynthesis, the molecular function of beta-glucosidase, and the cellular function of the kinesin complex (Figure 6A). The majority of DEGs in the S2 vs. S3 stages had the same pathway enrichment as the S1 vs. S2 and S1 vs. S3 stages (Figure 6B). These findings suggested that genes involved in the functional activity of beta-glucosidase and oxidoreductase made significant contributions to the fruit kernel development in M. ternifolia.

2.4. Overview of the Lipid Composition at Different Fruit Kernel Developmental Stages

The kernel of M. ternifolia has more than 60% oil content. To understand the lipid diversity and content in the fruit kernel, we analyzed the lipid profiles at different fruit kernel growth stages. In total, 408 different types of lipids were detected in our study (Supplementary Materials Table S2). Most of these lipids were categorized as glycerolipids, sterol lipids, sphingolipids, glycerophospholipids, and fatty acyls. These categories comprised major classes, including 118 triacylglycerols (TAGs), 31 diacylglycerols (DAGs), 50 ceramides (CERs), 4 sulfoquinovosyl diacylglycerols (SQDGs), 4 sphingomyelins (SMs), 4 stigmasterol hexosides (SGs), 14 acylated sterol glycosides (ASGs), 11 phosphatidylmethanols (PMeOHs), 19 phosphatidylcholine (PCs), 35 phosphatidylethanolamines (Pess), 11 phosphatidylglycerols (PGs), 10 phosphatidylinositols (PIs), 3 phosphatidic acids (PAs), 6 monogalactosyldiacylglycerols (MGDGs), 20 hexosyl ceramides (HexCers), 22 fatty acyls (FAs), 5 diacylglycerol glucuronides (DGGAs), and 12 digalactosyl diacylglycerols (DGDGs). It was observed that the contents of 155 subclasses of lipid were altered in the S2 stage compared with S1. Among all of these subclasses, 58 had high accumulation at S2 and 97 lipids exhibited a lower concentration during the S2 stage compared with S1 (Supplementary Materials Table S3). Moreover, three CERs, namely, 40:0, 42:2, and 44:1 subclasses, showed the highest fold change values (Figure 7A). All these lipids had a high concentration in the S2 stage of the fruit kernel. In addition, DAGs 32:1, 32:2, 40:3, 40:4, 42:3; HexCer 42:2; DGGA 34:1; and TG 50:5 lipid species were upregulated in the S2 stage. The top ten lipid subclasses with the lowest fold change values in the S2 stage were three DGDGs (36:3, 36:6, 38:6), four PEs (36:4, 38:1, 38:4, 38:5), MGDG 36:6, PC 36:1, and TAG 50:5 (Figure 7A).
The comparative analysis between the S1 and S3 stages revealed 261 lipid subclasses with differential concentrations (Supplementary Materials Table S4). From these, 187 had a high accumulation in S3 compared to S1 and only 74 subclasses of lipid showed a lower content in S3. Quantification showed that TAGs, including 49:3, 52:2, 52:3, 54:2, 56:2, 56:3, 56:4, 58:1, 60:2, and 62:2, had their highest levels in the S3 stage as compared with S1 (Figure 7B). Meanwhile, the top ten lipids with the lowest accumulation content in the S3 stage belonged to different subclasses (Figure 7B). For example, PEs 35:3, 36:5, 38:5, 40:3, 42:3; ASG 28:1; DGDG 38:6; DGGA 34:1; PA 36:3; and PC 36:5 exhibited higher contents during the S1 stage than S3. In total, 404 lipid subclasses had differential content accumulation in S3 compared to S2 (Supplementary Materials Table S5). A total of 203 lipids were obtained with a higher level, whereas 201 had a lower level in the S3 stage. Further results suggested that TAGs, such as 52:2, 52:3, 54:2, 56:1, 56:2;1O, 56:2, 56:3, 58:2, 60:2, and 61:2 had a significantly increased level in S3 than S2 (Figure 7C). In contrast, PEs 32:3, 34:4, 35:3, 36:5, 40:3, 42:4; ASGs 28:1, 29:1; DG 42:3; and PC 35:4 were highly reduced in the S3 stage (Figure 7C). It was noticed that there were few stage-specific lipids but there were more overlapped or common lipid profiles among the kernel developmental stages (Supplementary Materials Figure S2). Overall, the results of the lipid profiling revealed that during fruit kernel growth, triacylglycerols and diacylglycerol related lipids accumulation were increased. In contrast, phosphatidylethanolamine-associated lipids showed decreased accumulation.

2.5. Lipid-Related Gene Expression Changes during Fruit Kernel Development

To find the key genes associated with fruit kernel development, further mining of metabolic genes was achieved through the selection of DEGs having a fold change value ≥ 5 or ratio ≤ −5 as thresholds. The Venn diagram revealed 282 stage-specific highly expressed DEGs between S1 and S2 (Supplementary Materials Figure S2). The S1 and S3 stages comparison had 246 stage-specific DEGs. It was observed that S2 vs. S3 had 412 stage-specific DEGs with a high fold change ratio. Moreover, 400 overlapped DEGs among all stages were considered as a regulator of the metabolic process during fruit kernel growth. The expression trends of annotated genes associated with lipid biosynthesis are shown in Table 2. Significant low expression levels of four genes related to 3-ketoacyl-acyl carrier protein reductase (unigene 73227, unigene 211387, unigene 73228, unigene 26385) were observed during the S3 fruit kernel developmental stage compared to the S1 stage. The expression levels of five beta-glucosidase genes, namely, unigene 170150, unigene 81638, unigene 81642, unigene 81652, and unigene 91684, were highest during the S2 stage. Similar results were obtained for unigene 150581, unigene 150580, and unigene 140441, which encode the enoyl-acyl-carrier-protein reductase (NADH) enzyme. In addition, the diacylglycerol-kinase-2-linked unigene 147364 had an unpredictable expression during the fruit kernel development stage. A significantly altered expression was also determined for pyruvate-kinase-2-annotated unigene 63669, unigene 63675, and unigene 86831 during fruit kernel growth.
Triacylglycerols are the key resources of energy in the seeds. Two unigenes encoding glycerol-3-phosphate acyltransferase and denoted as unigene 152617 and unigene 128368 were found with low expressions during the late stage of fruit kernel growth in our study. Moreover, a significant change in the expression trends of some genes related to GDSL esterase lipase was obtained during the fruit kernel developmental stages. For example, unigene 240247, unigene 82774, unigene 183098, unigene 38911, unigene 113245, and unigene 238268 showed increased expression levels during the S1 stage compared with S2 and S3. In contrast, unigene 213504, unigene 62756, unigene 6967, and unigene 103298 had higher levels of expression in the S2 stage than S1 and S3. Globally, the expression of linoleate 13S-lipoxygenase 2-1 gene family members, such as unigene 9629, unigene 66716, unigene 54520, unigene 59928, unigene 213589, unigene 37584, unigene 34267, unigene 13471, unigene 5404, and unigene 139752, were significantly downregulated during the late stages of fruit kernel growth. Various genes involved in phospholipids regulation were determined with altered expression levels during the fruit kernel growth. Among these, unigene 190439, unigene 227869, and unigene 171746 had higher expression values in the S2 stage compared with S1 and S3, whereas unigene 158347 and unigene 122230 showed upregulated expression during the S1 (early fruit kernel growth) stage. Overall, this suggests that the stage-specific regulation of genes related to lipid metabolism and catabolism caused the diversity of lipid compositions and concentrations in M. ternifolia nuts.

2.6. Validation of Key Genes Involved in Lipid Regulation Using qRT-PCR Analysis

Based on the transcript downstream data analysis, ten genes putatively involved in the lipid regulation were selected to perform a real-time qRT-PCR analysis in order to confirm the accuracy of the RNA-seq. The selected genes showed differential expressions at the different developmental stages of the fruit kernel in M. ternifolia (Supplementary Materials Figure S3). It was observed that the expression trends of all genes were consistent with RNA-seq, thus confirming the reliability of RNA-seq results in our study.

3. Discussion

Out of four cultivated species, M. integrifolia and M. tetraphylla are widely grown in commercial orchards [40]. There are many constraints that limit the productivity of M. ternifolia. The major challenge is a high level of cyanogenic glycosides present in the kernel, which are toxic and unsuitable for human consumption without steeping and cooking [4,41]. A detailed understanding of the genetic mechanisms underlying lipid composition in the kernel of M. ternifolia will be useful to promote its commercial production level. Through comprehensive lipidomics and transcriptomic data analyses, this study improved our understanding of the molecular factors that influence oil accumulation in M. ternifolia kernels.

3.1. Divergence of Lipid Composition in the Kernels of M. ternifolia

A lipid is a form of stored carbon in plant seeds, yields more energy than carbohydrates, and is an integral part of plant organs [29]. The lipid composition was significantly changed in the developing fruit kernels of M. ternifolia. In our study, a total of 407 different types of lipid species were determined at different fruit kernel growth stages. The majority of these lipids belonged to the category of glycerolipids including triacylglycerols (TAGs). TAGs are the primary components of Macadamia oil, accounting for more than 95% of total lipids. These influence the nutritional characteristics of nut oil through physical and biochemical changes [42]. TAGs are entirely acylated derivatives of glycerol and the most abundant neutral lipid bodies present in seeds [43], senescence leaves, and pollen grains [44,45]. In the mature fruit kernels of M. ternifolia, the most accumulated species of TAGs were long-chain unsaturated fatty acids, including TAG 52:2;1O, TAG 54:2;1O, TAG 56:3;1O, and TAG 56:2;1O. An increase in unsaturated TAGs was previously observed in the later stages of seed development in Paeonia ostii [46], Glycine max [47], and Brassica napus [48]. The stored level of TAGs has a functional linkage to seed germination [49] and stress management [50]. Moreover, the nuts of Macadamia species are usually consumed in roasted form, which can modify the lipid profiles due to oxidation [51]. Excessive oxidation in unsaturated TAGs produces free fatty acids due to esterases and lipases through the activity of hydrolytic enzymes [52]. Free fatty acids are more vulnerable to oxidation during nut roasting, which causes adverse chemical changes through rancidity and the deterioration process, which ultimately leads to rejection by consumers due to a lower quality of nuts [53]. In a recent study, it was observed that roasting had a non-significant impact on the TAGs composition and it ultimately improved the flavor, aroma, color, texture, and appearance of macadamia nuts [42]. Among other detected glycerolipids in the developing fruit kernels of M. ternifolia, different species of diacylglycerols (DAGs) have a high abundance. In particular, DAGs 40:4 and 42:3 species had a high level of accumulation. DAGs are precursors of glycolipids, are swiftly phosphorylated to produce phosphatidic acid (PA) by diacylglycerol kinase enzyme, and act as a potential secondary messenger in plants [54]. The high contents of TAGs and DAGs in the fruit kernel of M. ternifolia may be useful to improve the oil quality through the biosynthesis of other essential fatty acids and lipids.
Phosphatidylethanolamine (PE) is a polar lipid belonging to the phospholipid category, is dominant in the plasma membrane, and plays an important role in cold stress tolerance [55]. In our study, the quantitative level of many PE species was dramatically reduced at the mature fruit kernel stage. Specifically, PE 35:3, PE 36:5, PE 40:3, PE 42:4, and PA 36:3 had higher levels of accumulation in the young fruit kernel of M. ternifolia. Previous research showed that Pes, together with phosphatidylcholine species, alleviates chilling injury to protect cell membrane damage in cold conditions [56]. PA is a type of glycerophospholipid that mediates the biosynthesis of phospholipids and glycolipids in plants. Moreover, it is involved in secondary signal transduction during abiotic stress, including salinity [57], drought [58], and cold [59]. Our analysis further reported a differential level of ceramides (CERs) in the medium fruit kernel stage, in particular, high levels of Cer 44:1 and Cer 42:2 contents were detected. CERs are sphingolipids that are essential components of the membrane in plants. Functional studies showed that CERs are involved in responses to biotic and abiotic stimuli. However, little is known about their signaling mechanism during the programmed cell death of certain bacteria and fungi pathogens [60]. PEs, PAs, and CERs have a critical role in the defense against various stressful conditions. A similar function in developing the kernel of M. ternifolia might be linked with a higher level of lipid stores in the nut. The data of lipidomic analysis identified the presence of some species of acylated sterol glycosides, phosphatidylmethanol, diacylglyceryl glucuronide, and digalactosyldiacylglycerol in the kernel of M. ternifolia. However, their non-differential trend suggested a minor role for the oil quality index. In contrast, a higher abundance of oil species, such as TAGS, DAGs, PEs, and PCs, in the mature nut is of substantial importance to the improve oil quality in M. ternifolia nut.

3.2. Lipid Biosynthesis in Fruit Kernels of M. ternifolia

Lipid biosynthesis in nuts is predominantly altered by the activity of genes that are involved in the de novo biosynthesis of fatty acids, synthesis of triacylglycerol, and formation of oil bodies [35]. The genetic mechanism of fatty acid and lipid synthesis leading to TAG is well known in Arabidopsis [61]. The substrate for all synthesized fatty acid pools is acetyl-CoA. The enzyme pyruvate kinase (PK) in the plastids generates pyruvate, which is catalyzed by the plastidial pyruvate dehydrogenase complex to form acetyl-CoA for the de novo production of fatty acids [62]. The altered expression of PK2 genes in young fruit kernels of M. ternifolia suggests their vital role in the generation of diverse fatty acids. Functional genomics experiments on PK2 genes in rice stated their importance in fatty acid and starch metabolism [63]. Among the fatty acid synthase system, genes associated with 3-ketoacyl-acyl carrier protein reductase (KAR) or synthase and enoyl-acyl-carrier-protein reductase (NADH) (ENR) regulate the synthesis of long-chain unsaturated fatty acids in plants seeds [64]. In the formation of long-chain fatty acids, the KAR enzyme reduces 1,3-ketobutyryl-ACP to 3-hydroxybutyryl-ACP, which extends a hydrocarbon chain of the acyl group up to a 16- or 18-carbon chain through a series of cycle reactions [65]. Similarly, ENR is predicted to be involved in the last step of fatty acid elongation in plants, bacteria, and mammals [66]. The altered expression of genes involved in the functional activity of KAR and ENR enzymes probably contributes to the development of unsaturated fatty acids in M. ternifolia. In this way, it might be interlinked with the unique biochemical and physical properties of nut oil. However, further research evidence is needed to support this statement. In our lipidomic data analysis, TAGs, DAGs, PEs, and PAs are the most synthesized forms of lipids present in the fruit kernel of M. ternifolia. The de novo lipid biosynthesis of M. ternifolia most probably starts with the catalytic activity of glycerol-3-phosphate acyltransferase (GPAT) enzyme. The GPAT initiates the preliminary reaction of TAG and phospholipids biosynthesis in the endoplasmic reticulum. It further converts glyerol-3-phosphate into lysophosphatidic acid, which produces phosphatidic acid. The dephosphorylation of phosphatidic acid forms DAG and finally TAG [67]. Furthermore, the oleoyl CoA transferred from plastids is used as a substrate to produce membrane lipids (PCs, PEs, and PIs) or the storage lipid TAG in M. ternifolia. However, further metabolic studies will be essential to prove these findings.
Several studies have reported the role of phosphoglycerolipids in lipid signaling [68,69]. The key genes involved in lipid signal transduction include diacylglycerol kinase (DGK), phosphatidylinositol (PI), nonspecific phospholipase C (NPC), phospholipase C (PLC), phospholipase D (PLD), and PA [70]. Interestingly, we observed substantially altered transcription levels of GPAT, DGK, PI, and NPC in developing fruit kernel M. ternifolia. This indicates their key roles not only in higher TAGs, DAGs, PE, and PA accumulation but also in lipid signaling during the growth of M. ternifolia kernel. It was reported that enzyme GDSL-type esterase/lipase can hydrolyze thioesters, aryl esters, phospholipids, and amino acids. In this way, they altered the mechanism of secondary metabolism [71]. A significant change in expression trends of genes related to the GDSL esterase lipase enzyme probably performed a similar function in developing the M. ternifolia fruit kernel. Furthermore, our research also identified several genes associated with the linoleate 9S-lipoxygenase (LOX) enzyme, particularly LOX6, which had a significantly modified expression in the M. ternifolia fruit kernel. This enzyme utilizes linoleic acid/linolenic acid as a substrate to produce unsaturated fatty acids, which improves the oil content and quality for human health [72]. In addition, LOX genes play a critical role in the enzymatic oxidation of polyunsaturated fatty acids to prevent germination during seed aging [73]. The β-glucosidases (BGlu) enzymes have numerous roles in plants, such as defense, cell wall lignification, and phytohormone activation. These enzymes are also responsible for the release of scented volatile compounds in plants [74,75]. The increased expression of BGlu genes in the developing fruit kernel of M. ternifolia suggests an association with the diversity of functions that directly or indirectly affect oil quality.
Concisely, our comparative lipidomic and transcriptomic data analysis identified key regulatory genes associated with oil content and quality in M. ternifolia kernels. Due to the wide public acceptance of Macadamia nuts as a portion of quality food, these results are vital for promoting the nutritional and commercial value of M. ternifolia kernels. Our results provide a baseline for further transgenic and genetic engineering research with the aim to identify oil candidate genes and their regulators. Hence, further works are needed to boost the application of M. ternifolia nuts in the food, pharmaceutical, and cosmetic industries.

4. Materials and Methods

4.1. Plant Materials, RNA Extraction, and Sequencing

The plant material used in this study was collected from M. ternifolia plants cultivated in the Yunnan province of China. The developing fruit kernels comprised of young (S1), medium-aged (S2), and mature (S3) fruiting stages were picked in three biological repeats. The harvested samples were immediately placed in liquid nitrogen and stored at −80 °C before further processing. Total RNA was isolated from all nine samples with a Trizol reagent (Invitrogen, San Diego, CA, USA) following the manufacturer’s standard protocol. RNase-free DNase I (TaKaRa, Kyoto, Japan) was used to purify the total RNA. The RNA quality and quantity were determined by 1% agarose gel and Nanodrop spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA), respectively. Sequencing libraries were prepared with the Illumina TruSeq Standard mRNA library preparation kit following the manufacturer’s standard protocol. After the quality tests, the Illumina HiSeq Ten X platform (Illumina Inc., San Diego, CA, USA) was used according to the recommended protocol to perform de novo pair-end RNA sequencing.

4.2. Transcript Assembly, Annotation, and Analysis of DEGs

By removing the adaptors, ambiguous bases, and low-quality reads, the original sequencing data in the form of raw reads were filtered to obtain high-quality clean reads. Because of the non-availability of the reference genome of M. ternifolia, high-quality clean reads were assembled into unigenes for subsequent analysis with Trinity 2.6.6 software [76]. The BLAST software was used to compare the unigenes sequence with the Kyoto Encyclopedia of Genes and Genomes database (KEGG), NCBI non-redundant protein sequences (NR), Gene Ontology (GO), euKaryotic Ortholog Groups (KOG), Protein family (Pfam), Trembl, and Swiss-Prot databases to retrieve the unigenes’ functional annotations [77]. The expression levels of the mapped transcript were determined in the form of fragments per kilobase of transcript per million fragments mapped (FPKM) values. The R package DESeq2 version 1.22.2 [78] was used to obtain the total number of differential genes and the number of upregulated and downregulated genes among two different samples. The significant criteria for DEGs included log2 fold change ≥1 and ≤−1 and a false discovery rate < 0.05. The KEGG for linking genomes to life and the environment platform was used to perform the pathway enrichment analysis [79].

4.3. Real-Time qRT-PCR Analysis

To perform the qRT-PCR analysis for ten selected genes, the cDNAs were produced with iScript cDNA Synthesis Kit (BIO-RAD). Primer Premier 5.0 was used to design the gene-specific reverse and forward primers (Table S6). The reaction mixture in three technical repeats was prepared with iTaq Universal SYBR Green Supermix 50 mL (BIO-RAD) kit. The quantitative PCR analysis was conducted in triplicate biological replicates on a BIO-RAD iCycler iQ5 real-time PCR system. The Actin2 gene was used as an internal control. The running and data analysis protocol of the qRT-PCR was followed as detailed in a previous report [80].

4.4. Lipid Extraction and Analysis Using UPLC-MS/MS

Different fruit kernels consisting of young (S1), medium-aged (S2), and mature (S3) fruiting stages were utilized to extract lipids. In brief, a 15 mg sample was taken to make a fine powder. First, 2 mL of precooled (−20 °C) methanol was added, followed by 4 mL dichloromethane, and the mixture was vortexed for 1 h. After 10 min of incubation at 4 °C and sonication for 10 min in an ultra-sonication bath, 1.6 mL double-distilled water was added. Samples were then centrifuged for 5 min. Two phases were produced after the centrifugation. The upper clean and residual-free lipophilic phase was carefully removed. The lipophilic phase was then dried under a vacuum. Later, 1 mL isopropyl alcohol was added in the lipophilic phase and used for further lipidomic analysis. The lipid profiling was executed on a 6600 plus AccurateMass Q-TOF (AB Sciex TripleTOF® 6600, Framingham, MA, USA) mass spectrometer system. The process of liquid chromatography was completed using UPLC with an autosampler LC-30A (Shimadzu Corporation) column (100 × 2.1 mm, 2.6 m). The lipid measurement protocol was followed as in the description of Wan et al. [81]. The MSDIAL ver. 4.0, PeakView2.1, and MultiQuant 3.0 were used to perform the quantitative and qualitative data analyses of the lipids [82].

5. Conclusions

Our results revealed lipid and transcript profiles at various developing fruit kernel stages of M. ternifolia nuts. Moreover, the key lipids and their putative candidate genes were identified. In particular, GPAT, DGK, PI, NPC, PK, KAR, and LOX enzymes encoding genes were most probably rate-limiting for the biosynthesis of glycerolipids and phospholipid in M. ternifolia nuts. These results help to understand lipid diversity and provide first insights into the genetic mechanism of lipid accumulation. Gene functional characterization and further transcript data from different genotypes are mandatory to explore the regulatory network of lipid biosynthesis in M. ternifolia nuts.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11121431/s1, Table S1: Statistics of assembly for the different fruit kernel stages of M. ternifolia. Table S2: Total diversity of lipids and their content level in the different fruit kernel stages of M. ternifolia. Table S3: Total contents of lipids differentially accumulation in the comparison of the S1 vs. S2 stages. Table S4: Total contents of lipids differentially accumulation in the comparison of the S1 vs. S3 stages. Table S5: Total contents of lipids differentially accumulation in the comparison of the S2 vs. S3 stages. Table S6: Primer sequences used in the qRT-PCR experiment. Figure S1: Distribution of the overlapped and unique number of lipid species among the different fruit kernel stages of M. ternifolia. Figure S2: Distribution of the overlapped and unique highly expressed metabolic DEGs among the different fruit kernel stages of M. ternifolia. Figure S3: qRT-PCR validation of selected differentially expressed genes. The Actin2 gene was used as the internal control and the transcript levels in S2 and S3 were expressed as compared with S1. Young (S1), medium-aged (S2), and mature (S3) fruiting stages.

Author Contributions

Conceptualization, R.S., Z.X., B.L. and H.B.; methodology, R.S. and W.W.; software, R.S., H.B. and B.L.; validation, Z.X., Z.Y., W.W., B.L. and H.B.; formal analysis, H.B., Z.X., Z.Y., B.L. and W.W.; investigation, B.L., Z.X., Z.Y., W.W. and H.B.; resources, R.S. and C.L.; data curation, W.W. writing—original draft preparation, R.S., B.L. and H.B.; writing—review and editing, R.S., B.L.; visualization, Z.X.; supervision, Z.X. and W.W.; project administration, C.L. and B.L.; funding acquisition, Z.X. and W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Guangxi Natural Science Foundation (2019GXNSFBA185011), macadamia nut research team (2021YT156), Guangxi Academy of Agricultural Sciences Science and Technology Development Fund (2020YM141, 2020YM46, 2019ZX2015, 2019ZX17), Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (no. 1630062017034), Major Science and Technology Project of Yunnan (202102AE090042, 202002AA10007, 2019ZG00901), Major Science and Technology Project of Kunming (2021JH002), High-End Foreign Experts Program of Yunnan (202105AQ130011, 2019013), and Yunnan Provincial Financial Forestry Science and Technology Promotion Demonstration Special Project (2020) no. TS09.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw transcriptome data were submitted to NCBI SRA under the project number PRJNA761018.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Niu, Y.-F.; Li, G.-H.; Ni, S.-B.; He, X.-Y.; Zheng, C.; Liu, Z.-Y.; Gong, L.-D.; Liu, J. Genome assembly and annotation of Macadamia tetraphylla. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  2. Toft, B.D.; Alam, M.; Topp, B. Estimating genetic parameters of architectural and reproductive traits in young macadamia cultivars. Tree Genet. Genomes 2018, 14, 1–10. [Google Scholar] [CrossRef]
  3. Mast, A.R.; Willis, C.L.; Jones, E.H.; Downs, K.M.; Weston, P.H. A smaller Macadamia from a more vagile tribe: Inference of phylogenetic relationships, divergence times, and diaspore evolution in Macadamia and relatives (tribe Macadamieae; Proteaceae). Am. J. Bot. 2008, 95, 843–870. [Google Scholar] [CrossRef]
  4. Hu, W.; Fitzgerald, M.; Topp, B.; Alam, M.; O’Hare, T.J. A review of biological functions, health benefits, and possible de novo biosynthetic pathway of palmitoleic acid in macadamia nuts. J. Funct. Foods 2019, 62, 103520. [Google Scholar] [CrossRef]
  5. Navarro, S.L.; Rodrigues, C.E. Macadamia oil extraction methods and uses for the defatted meal byproduct. Trends Food Sci. Technol. 2016, 54, 148–154. [Google Scholar] [CrossRef]
  6. Insanu, M.; Hartati, R.; Bajri, F.; Fidrianny, I. Macadamia Genus: An Updated Review of Phytochemical Compounds and Pharmacological Activities. Biointerface Res. Appl. Chem. 2021, 11, 14898–14905. [Google Scholar]
  7. Hong, M.Y.; Groven, S.; Marx, A.; Rasmussen, C.; Beidler, J.J.M. Anti-inflammatory, antioxidant, and hypolipidemic effects of mixed nuts in atherogenic diet-fed rats. Molecules 2018, 23, 3126. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Hardner, C.M.; Wall, M.; Cho, A.J.H. Global macadamia science: Overview of the special section. HortScience 2019, 54, 592–595. [Google Scholar] [CrossRef] [Green Version]
  9. Hardner, C.; e Silva, J.C.; Williams, E.; Meyers, N.; McConchie, C.J.H. Breeding new cultivars for the Australian macadamia industry. HortScience 2019, 54, 621–628. [Google Scholar] [CrossRef] [Green Version]
  10. Penoni, E.d.S.; Pio, R.; Rodrigues, F.A.; Maro, L.A.C.; Costa, F.C.J.C.R. Analysis of fruits and nuts of macadamia walnut cultivars. Ciência Rural 2011, 41, 2080–2083. [Google Scholar] [CrossRef] [Green Version]
  11. Piza, P.L.B.T.; Moriya, L.M. Cultivo da macadâmia no Brasil. Rev. Bras. Frutic. 2014, 36, 39–45. [Google Scholar] [CrossRef] [Green Version]
  12. Venkatachalam, M.; Sathe, S.K. Chemical composition of selected edible nut seeds. J. Agric. Food Chem. 2006, 54, 4705–4714. [Google Scholar] [CrossRef]
  13. Rengel, A.; Pérez, E.; Piombo, G.; Ricci, J.; Servent, A.; Tapia, M.S.; Gibert, O.; Montet, D. Lipid profile and antioxidant activity of macadamia nuts (Macadamia integrifolia) cultivated in Venezuela. Nat. Sci. 2015, 12, 535–547. [Google Scholar]
  14. Knothe, G. Biodiesel derived from a model oil enriched in palmitoleic acid, macadamia nut oil. Energy Fuels 2010, 24, 2098–2103. [Google Scholar] [CrossRef]
  15. Wang, D.D.; Hu, F.B. Dietary fat and risk of cardiovascular disease: Recent controversies and advances. Annu. Rev. Nutr. 2017, 37, 423–446. [Google Scholar] [CrossRef] [PubMed]
  16. Garg, M.L.; Blake, R.J.; Wills, R.B. Macadamia nut consumption lowers plasma total and LDL cholesterol levels in hypercholesterolemic men. J. Nutr. 2003, 133, 1060–1063. [Google Scholar] [CrossRef] [Green Version]
  17. Kris-Etherton, P.M.; Yu-Poth, S.; Sabaté, J.; Ratcliffe, H.E.; Zhao, G.; Etherton, T.D. Nuts and their bioactive constituents: Effects on serum lipids and other factors that affect disease risk. Am. J. Clin. Nutr. 1999, 70, 504s–511s. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Vadivel, V.; Kunyanga, C.N.; Biesalski, H.K. Health benefits of nut consumption with special reference to body weight control. Nutrition 2012, 28, 1089–1097. [Google Scholar] [CrossRef] [PubMed]
  19. Kaijser, A.; Dutta, P.; Savage, G. Oxidative stability and lipid composition of macadamia nuts grown in New Zealand. Food Chem. 2000, 71, 67–70. [Google Scholar] [CrossRef]
  20. Kochhar, S.P.; Henry, C.J.K. Oxidative stability and shelf-life evaluation of selected culinary oils. Int. J. Food Sci. Nutr. 2009, 60, 289–296. [Google Scholar] [CrossRef]
  21. Yoon, W.-J.; Kim, M.-J.; Moon, J.-Y.; Kang, H.-J.; Kim, G.-O.; Lee, N.H.; Hyun, C.-G. Effect of palmitoleic acid on melanogenic protein expression in murine b16 melanoma. J. Oleo Sci. 2010, 59, 315–319. [Google Scholar] [CrossRef] [Green Version]
  22. Kuo, T.M.; Gardner, H. Perspectives on the Production of Industrial Oils in Genetically Engineered Oilseeds. In Lipid Biotechnology; CRC Press: Boca Raton, FL, USA, 2002; pp. 107–117. [Google Scholar]
  23. Drackley, J. Lipid Metabolism; CABI: Wallingford, UK, 2000; Volume 1, pp. 97–119. [Google Scholar]
  24. Mutz, K.-O.; Heilkenbrinker, A.; Lönne, M.; Walter, J.-G.; Stahl, F. Transcriptome analysis using next-generation sequencing. Curr. Opin. Biotechnol. 2013, 24, 22–30. [Google Scholar] [CrossRef] [PubMed]
  25. Shahzad, K.; Zhang, X.; Guo, L.; Qi, T.; Tang, H.; Zhang, M.; Zhang, B.; Wang, H.; Qiao, X.; Feng, J. Comparative transcriptome analysis of inbred lines and contrasting hybrids reveals overdominance mediate early biomass vigor in hybrid cotton. BMC Genom. 2020, 21, 1–16. [Google Scholar]
  26. Beisson, F.; Koo, A.J.; Ruuska, S.; Schwender, J.; Pollard, M.; Thelen, J.J.; Paddock, T.; Salas, J.J.; Savage, L.; Milcamps, A. Arabidopsis genes involved in acyl lipid metabolism. A 2003 census of the candidates, a study of the distribution of expressed sequence tags in organs, and a web-based database. Plant Physiol. 2003, 132, 681–697. [Google Scholar] [PubMed] [Green Version]
  27. Severin, A.J.; Woody, J.L.; Bolon, Y.-T.; Joseph, B.; Diers, B.W.; Farmer, A.D.; Muehlbauer, G.J.; Nelson, R.T.; Grant, D.; Specht, J.E. RNA-Seq Atlas of Glycine max: A guide to the soybean transcriptome. BMC Plant Biol. 2010, 10, 1–16. [Google Scholar] [CrossRef] [Green Version]
  28. O’Hara, P.; Slabas, A.R.; Fawcett, T. Fatty acid and lipid biosynthetic genes are expressed at constant molar ratios but different absolute levels during embryogenesis. Plant Physiol. 2002, 129, 310–320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Gu, K.; Yi, C.; Tian, D.; Sangha, J.S.; Hong, Y.; Yin, Z. Expression of fatty acid and lipid biosynthetic genes in developing endosperm of Jatropha curcas. Biotechnol. Biofuels 2012, 5, 1–15. [Google Scholar] [CrossRef] [Green Version]
  30. Gupta, K.; Kayam, G.; Faigenboim-Doron, A.; Clevenger, J.; Ozias-Akins, P.; Hovav, R. Gene expression profiling during seed-filling process in peanut with emphasis on oil biosynthesis networks. Plant Sci. 2016, 248, 116–127. [Google Scholar] [CrossRef]
  31. Zhu, Y.; Wang, K.; Wu, C.; Zhao, Y.; Yin, X.; Zhang, B.; Grierson, D.; Chen, K.; Xu, C. Effect of ethylene on cell wall and lipid metabolism during alleviation of postharvest chilling injury in peach. Cells 2019, 8, 1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Brown, A.P.; Kroon, J.T.; Swarbreck, D.; Febrer, M.; Larson, T.R.; Graham, I.A.; Caccamo, M.; Slabas, A.R. Tissue-specific whole transcriptome sequencing in castor, directed at understanding triacylglycerol lipid biosynthetic pathways. PLoS ONE 2012, 7, e30100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Xu, Z.; Ni, J.; Shah, F.A.; Wang, Q.; Wang, Z.; Wu, L.; Fu, S. Transcriptome analysis of pecan seeds at different developing stages and identification of key genes involved in lipid metabolism. PLoS ONE 2018, 13, e0195913. [Google Scholar] [CrossRef]
  34. Bates, P.D.; Stymne, S.; Ohlrogge, J. Biochemical pathways in seed oil synthesis. Curr. Opin. Plant Biol. 2013, 16, 358–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dussert, S.; Guerin, C.; Andersson, M.; Joët, T.; Tranbarger, T.J.; Pizot, M.; Sarah, G.; Omore, A.; Durand-Gasselin, T.; Morcillo, F. Comparative transcriptome analysis of three oil palm fruit and seed tissues that differ in oil content and fatty acid composition. Plant Physiol. 2013, 162, 1337–1358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Yang, W.; Simpson, J.P.; Li-Beisson, Y.; Beisson, F.; Pollard, M.; Ohlrogge, J.B. A land-plant-specific glycerol-3-phosphate acyltransferase family in Arabidopsis: Substrate specificity, sn-2 preference, and evolution. Plant Physiol. 2012, 160, 638–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Nandi, A.; Welti, R.; Shah, J. The Arabidopsis thaliana dihydroxyacetone phosphate reductase gene Suppressor of fatty acid desaturase deficiency1 is required for glycerolipid metabolism and for the activation of systemic acquired resistance. Plant Cell 2004, 16, 465–477. [Google Scholar] [CrossRef] [Green Version]
  38. Li, R.-J.; Gao, X.; Li, L.-M.; Liu, X.-L.; Wang, Z.-Y.; Lü, S.-Y. De novo assembly and characterization of the fruit transcriptome of Idesia polycarpa reveals candidate genes for lipid biosynthesis. Front. Plant Sci. 2016, 7, 801. [Google Scholar] [CrossRef] [Green Version]
  39. Nock, C.J.; Baten, A.; Barkla, B.J.; Furtado, A.; Henry, R.J.; King, G.J. Genome and transcriptome sequencing characterises the gene space of Macadamia integrifolia (Proteaceae). BMC Genom. 2016, 17, 1–12. [Google Scholar] [CrossRef] [Green Version]
  40. Howlett, B.G.; Nelson, W.R.; Pattemore, D.E.; Gee, M. Pollination of macadamia: Review and opportunities for improving yields. Sci. Hortic. 2015, 197, 411–419. [Google Scholar] [CrossRef]
  41. Monaghan, E.K. Chemical Composition and Protein Antigenicity–Almond (Prunus dulcis) and Macadamia Nut (Macadamia integrifolia) Seeds; The Florida State University: Tallahassee, FL, USA, 2008. [Google Scholar]
  42. Tu, X.-H.; Wu, B.-f.; Xie, Y.; Xu, S.-L.; Wu, Z.-Y.; Lv, X.; Wei, F.; Du, L.-Q.; Chen, H. A comprehensive study of raw and roasted macadamia nuts: Lipid profile, physicochemical, nutritional, and sensory properties. Food Sci. Nutr. 2021, 9, 1688–1697. [Google Scholar] [CrossRef] [PubMed]
  43. Xiu, Y.; Wu, G.; Tang, W.; Peng, Z.; Bu, X.; Chao, L.; Yin, X.; Xiong, J.; Zhang, H.; Zhao, X. Oil biosynthesis and transcriptome profiles in developing endosperm and oil characteristic analyses in Paeonia ostii var. lishizhenii. J. Plant Physiol. 2018, 228, 121–133. [Google Scholar] [CrossRef]
  44. Kaup, M.T.; Froese, C.D.; Thompson, J.E. A role for diacylglycerol acyltransferase during leaf senescence. Plant Physiol. 2002, 129, 1616–1626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kim, H.U.; Hsieh, K.; Ratnayake, C.; Huang, A.H. A novel group of Oleosins is present inside the pollen ofArabidopsis. J. Biol. Chem. 2002, 277, 22677–22684. [Google Scholar] [CrossRef] [Green Version]
  46. Yu, S.-Y.; Zhang, Y.; Lyu, Y.-P.; Yao, Z.-J.; Hu, Y.-H. Lipidomic profiling of the developing kernel clarifies the lipid metabolism of Paeonia ostii. Sci. Rep. 2021, 11, 1–12. [Google Scholar] [CrossRef] [PubMed]
  47. Li, M.; Butka, E.; Wang, X. Comprehensive quantification of triacylglycerols in soybean seeds by electrospray ionization mass spectrometry with multiple neutral loss scans. Sci. Rep. 2014, 4, 1–11. [Google Scholar] [CrossRef] [Green Version]
  48. Woodfield, H.K.; Cazenave-Gassiot, A.; Haslam, R.P.; Guschina, I.A.; Wenk, M.R.; Harwood, J.L. Using lipidomics to reveal details of lipid accumulation in developing seeds from oilseed rape (Brassica napus L.). Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2018, 1863, 339–348. [Google Scholar] [CrossRef] [PubMed]
  49. Poxleitner, M.; Rogers, S.W.; Lacey Samuels, A.; Browse, J.; Rogers, J.C. A role for caleosin in degradation of oil-body storage lipid during seed germination. Plant J. 2006, 47, 917–933. [Google Scholar] [CrossRef]
  50. Volk, G.M.; Crane, J.; Caspersen, A.M.; Hill, L.M.; Gardner, C.; Walters, C. Massive cellular disruption occurs during early imbibition of Cuphea seeds containing crystallized triacylglycerols. Planta 2006, 224, 1415–1426. [Google Scholar] [CrossRef]
  51. Lin, J.-T.; Liu, S.-C.; Hu, C.-C.; Shyu, Y.-S.; Hsu, C.-Y.; Yang, D.-J. Effects of roasting temperature and duration on fatty acid composition, phenolic composition, Maillard reaction degree and antioxidant attribute of almond (Prunus dulcis) kernel. Food Chem. 2016, 190, 520–528. [Google Scholar] [CrossRef] [PubMed]
  52. Wall, M.M. Functional lipid characteristics, oxidative stability, and antioxidant activity of macadamia nut (Macadamia integrifolia) cultivars. Food Chem. 2010, 121, 1103–1108. [Google Scholar] [CrossRef]
  53. Colzato, M.; Scramin, J.; Forato, L.; Colnago, L.; Assis, O. 1h Nmr investigation of oil oxidation in macadamia nuts coated with zein-based films. J. Food Process. Preserv. 2011, 35, 790–796. [Google Scholar] [CrossRef]
  54. Ha, K.S.; Thompson, G.A., Jr. Diacylglycerol metabolism in the green alga Dunaliella salina under osmotic stress: Possible role of diacylglycerols in phospholipase C-mediated signal transduction. Plant Physiol. 1991, 97, 921–927. [Google Scholar] [CrossRef] [Green Version]
  55. Degenkolbe, T.; Giavalisco, P.; Zuther, E.; Seiwert, B.; Hincha, D.K.; Willmitzer, L. Differential remodeling of the lipidome during cold acclimation in natural accessions of Arabidopsis thaliana. Plant J. 2012, 72, 972–982. [Google Scholar] [CrossRef]
  56. Chen, M.; Thelen, J.J. Acyl-lipid desaturase2 is required for chilling and freezing tolerance in Arabidopsis. Plant Cell 2013, 25, 1430–1444. [Google Scholar] [CrossRef] [Green Version]
  57. Zhang, Q.; Lin, F.; Mao, T.; Nie, J.; Yan, M.; Yuan, M.; Zhang, W. Phosphatidic acid regulates microtubule organization by interacting with MAP65-1 in response to salt stress in Arabidopsis. Plant Cell 2012, 24, 4555–4576. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Li, Y.; Tan, Y.; Shao, Y.; Li, M.; Ma, F. Comprehensive genomic analysis and expression profiling of diacylglycerol kinase gene family in Malus prunifolia (Willd.) Borkh. Gene 2015, 561, 225–234. [Google Scholar] [CrossRef]
  59. Ruelland, E.; Cantrel, C.; Gawer, M.; Kader, J.-C.; Zachowski, A. Activation of phospholipases C and D is an early response to a cold exposure in Arabidopsis suspension cells. Plant Physiol. 2002, 130, 999–1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Liang, H.; Yao, N.; Song, J.T.; Luo, S.; Lu, H.; Greenberg, J.T. Ceramides modulate programmed cell death in plants. Genes Dev. Cell 2003, 17, 2636–2641. [Google Scholar] [CrossRef] [Green Version]
  61. Li-Beisson, Y.; Shorrosh, B.; Beisson, F.; Andersson, M.X.; Arondel, V.; Bates, P.D.; Baud, S.; Bird, D.; DeBono, A.; Durrett, T.P. Acyl-lipid metabolism. In The Arabidopsis Book; American Society of Plant Biologists: Rockville, MD, USA, 2013; Volume 11. [Google Scholar]
  62. Andre, C.; Froehlich, J.E.; Moll, M.R.; Benning, C. A heteromeric plastidic pyruvate kinase complex involved in seed oil biosynthesis in Arabidopsis. Plant Cell 2007, 19, 2006–2022. [Google Scholar] [CrossRef] [Green Version]
  63. Cai, Y.; Li, S.; Jiao, G.; Sheng, Z.; Wu, Y.; Shao, G.; Xie, L.; Peng, C.; Xu, J.; Tang, S. Os PK 2 encodes a plastidic pyruvate kinase involved in rice endosperm starch synthesis, compound granule formation and grain filling. Plant Biotechnol. J. 2018, 16, 1878–1891. [Google Scholar] [CrossRef] [Green Version]
  64. Moreno-Pérez, A.J.; Sánchez-García, A.; Salas, J.J.; Garcés, R.; Martínez-Force, E. Acyl-ACP thioesterases from macadamia (Macadamia tetraphylla) nuts: Cloning, characterization and their impact on oil composition. Plant Physiol. Biochem. 2011, 49, 82–87. [Google Scholar] [CrossRef]
  65. Bhatla, S.C.; Lal, M.A. Plant Physiology, Development and Metabolism; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  66. Massengo-Tiassé, R.P.; Cronan, J.E. Diversity in enoyl-acyl carrier protein reductases. Cell. Mol. Life Sci. 2009, 66, 1507–1517. [Google Scholar] [CrossRef] [Green Version]
  67. Cagliari, A.; Margis-Pinheiro, M.; Loss, G.; Mastroberti, A.A.; de Araujo Mariath, J.E.; Margis, R. Identification and expression analysis of castor bean (Ricinus communis) genes encoding enzymes from the triacylglycerol biosynthesis pathway. Plant Sci. 2010, 179, 499–509. [Google Scholar] [CrossRef]
  68. Eichmann, T.O.; Lass, A. DAG tales: The multiple faces of diacylglycerol-stereochemistry, metabolism, and signaling. Cell. Mol. Life Sci. 2015, 72, 3931–3952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Liu, Y.; Su, Y.; Wang, X. Phosphatidic acid-mediated signaling. Lipid Mediat. Protein Signal. 2013, 991, 159–176. [Google Scholar]
  70. Wang, X.; Chapman, K.D. Lipid signaling in plants. Front. Plant Sci. 2013, 4, 216. [Google Scholar] [CrossRef] [Green Version]
  71. Akoh, C.C.; Lee, G.-C.; Liaw, Y.-C.; Huang, T.-H.; Shaw, J.-F. GDSL family of serine esterases/lipases. Prog. Lipid Res. 2004, 43, 534–552. [Google Scholar] [CrossRef] [PubMed]
  72. Liu, Y.; Zhou, F.; Huang, X.; Wang, W.; Zhang, S.; Feng, F. Identification and integrated analysis of mRNAs, lncRNAs, and microRNAs of developing seeds in high oleic acid sunflower (Helianthus annuus L.). Acta Physiol. Plant. 2021, 43, 1–13. [Google Scholar] [CrossRef]
  73. Oenel, A.; Fekete, A.; Krischke, M.; Faul, S.C.; Gresser, G.; Havaux, M.; Mueller, M.J.; Berger, S. Enzymatic and non-enzymatic mechanisms contribute to lipid oxidation during seed aging. Plant Cell Physiol. 2017, 58, 925–933. [Google Scholar]
  74. Opassiri, R.; Pomthong, B.; Onkoksoong, T.; Akiyama, T.; Esen, A.; Cairns, J.R.K. Analysis of rice glycosyl hydrolase family 1 and expression of Os4bglu12 β-glucosidase. BMC Plant Biol. 2006, 6, 1–19. [Google Scholar] [CrossRef] [Green Version]
  75. Ketudat Cairns, J.R.; Esen, A. β-Glucosidases. Cell. Mol. Life Sci. 2010, 67, 3389–3405. [Google Scholar]
  76. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. Trinity: Reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 2011, 29, 644. [Google Scholar] [CrossRef] [Green Version]
  77. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  78. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [Green Version]
  79. Kanehisa, M.; Araki, M.; Goto, S.; Hattori, M.; Hirakawa, M.; Itoh, M.; Katayama, T.; Kawashima, S.; Okuda, S.; Tokimatsu, T. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2007, 36, D480–D484. [Google Scholar] [CrossRef]
  80. Dossa, K.; Zhou, R.; Li, D.; Liu, A.; Qin, L.; Mmadi, M.A.; Su, R.; Zhang, Y.; Wang, J.; Gao, Y. A novel motif in the 5′-UTR of an orphan gene ‘Big Root Biomass’ modulates root biomass in sesame. Plant Biotechnol. J. 2021, 19, 1065–1079. [Google Scholar] [CrossRef] [PubMed]
  81. Wan, H.; Liu, H.; Zhang, J.; Lyu, Y.; Li, Z.; He, Y.; Zhang, X.; Deng, X.; Brotman, Y.; Fernie, A.R. Lipidomic and transcriptomic analysis reveals reallocation of carbon flux from cuticular wax into plastid membrane lipids in a glossy “Newhall” navel orange mutant. Hortic. Res. 2020, 7, 1–18. [Google Scholar] [CrossRef] [Green Version]
  82. Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Kanazawa, M.; VanderGheynst, J.; Fiehn, O.; Arita, M. MS-DIAL: Data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat. Methods 2015, 12, 523–526. [Google Scholar] [CrossRef]
Figure 1. Length distribution and functional annotations of unigenes in M. ternifolia fruit kernels: (A) Overall length distribution of sequences and (B) functional annotations for unigenes in public databases.
Figure 1. Length distribution and functional annotations of unigenes in M. ternifolia fruit kernels: (A) Overall length distribution of sequences and (B) functional annotations for unigenes in public databases.
Life 11 01431 g001
Figure 2. Principal component analysis and correlation of M. ternifolia fruit kernels: (A) principal component analysis and (B) correlation analysis between the samples.
Figure 2. Principal component analysis and correlation of M. ternifolia fruit kernels: (A) principal component analysis and (B) correlation analysis between the samples.
Life 11 01431 g002
Figure 3. Distribution of total DEGs in different developing stages of M. ternifolia fruit kernels: (A) division of up- and downregulated DEGs and (B) division of overlapped and unique DEGs.
Figure 3. Distribution of total DEGs in different developing stages of M. ternifolia fruit kernels: (A) division of up- and downregulated DEGs and (B) division of overlapped and unique DEGs.
Life 11 01431 g003
Figure 4. Functional and pathway enrichment analysis for DEGs in the young and medium-aged (S1 vs. S2) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Figure 4. Functional and pathway enrichment analysis for DEGs in the young and medium-aged (S1 vs. S2) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Life 11 01431 g004
Figure 5. Functional and pathway enrichment analysis for DEGs in the young and mature (S1 vs. S3) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Figure 5. Functional and pathway enrichment analysis for DEGs in the young and mature (S1 vs. S3) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Life 11 01431 g005
Figure 6. Functional and pathway enrichment analysis for DEGs in the medium-aged and mature (S2 vs. S3) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Figure 6. Functional and pathway enrichment analysis for DEGs in the medium-aged and mature (S2 vs. S3) stages of M. ternifolia fruit kernels: (A) functional enrichment and (B) pathway enrichment.
Life 11 01431 g006
Figure 7. Heatmaps of the lipid content that had the highest differential accumulation among the different stages of the M. ternifolia fruit kernels: (A) in the comparison of the young (S1) and medium-aged (S2) stages, (B) in the comparison of the young (S1) and mature (S3) stages, and (C) in the comparison of the medium-aged (S1) and mature (S3) stages.
Figure 7. Heatmaps of the lipid content that had the highest differential accumulation among the different stages of the M. ternifolia fruit kernels: (A) in the comparison of the young (S1) and medium-aged (S2) stages, (B) in the comparison of the young (S1) and mature (S3) stages, and (C) in the comparison of the medium-aged (S1) and mature (S3) stages.
Life 11 01431 g007
Table 1. Overview of the transcriptome sequencing dataset and quality check.
Table 1. Overview of the transcriptome sequencing dataset and quality check.
SampleRaw ReadsClean ReadsClean Base (Gb)Mapped Reads (%)Q20 (%)Q30 (%)GC (%)
S1-155,604,26853,433,1468.08597.993.844.8
S1-256,740,04854,978,8588.38397.592.845.4
S1-350,745,98049,452,8187.48397.993.745.2
S2-149,242,86245,820,0886.88297.893.544.1
S2-246,406,95844,626,7006.78297.893.544.3
S2-348,601,19646,373,8827.08197.993.744.1
S3-155,107,46653,593,0488.07997.893.645.7
S3-246,080,42844,529,8226.77897.893.645.2
S3-356,758,32654,256,3488.17698.094.244.0
Table 2. Expression profile of genes that encoded enzymes involved in lipid biosynthesis in the M. ternifolia kernels.
Table 2. Expression profile of genes that encoded enzymes involved in lipid biosynthesis in the M. ternifolia kernels.
IDS1S2S3Annotations
Unigene 7322761.91.80.03-ketoacyl-acyl carrier protein reductase
Unigene 21138711.50.20.03-ketoacyl-acyl carrier protein reductase
Unigene 7322810.60.00.03-ketoacyl-acyl carrier protein reductase
Unigene 1632000.01.20.03-ketoacyl-acyl carrier protein reductase
Unigene 2638541.60.00.03-ketoacyl-acyl carrier protein reductase
Unigene 170150108.4640.63.6Beta-glucosidase 11
Unigene 4108527.40.20.3Beta-glucosidase 11
Unigene 816380.318.40.0Beta-glucosidase 12
Unigene 816521.322.10.0Beta-glucosidase 12
Unigene 916840.00.80.0Beta-glucosidase 44
Unigene 350022.60.10.0Beta-glucosidase 41
Unigene 1018470.00.01.7Choline/ethanolaminephosphotransferase 1
Unigene 8164256.3479.20.1Cyanogenic beta-glucosidase
Unigene 1473640.05.82.7Diacylglycerol kinase 2
Unigene 1505810.00.70.0Enoyl-acyl-carrier-protein reductase (NADH)
Unigene 1505800.01.40.0Enoyl-acyl-carrier-protein reductase (NADH)
Unigene 1404410.03.00.0Enoyl-acyl-carrier-protein reductase (NADH)
Unigene 2402474.60.00.1GDSL esterase lipase
Unigene 2135043.710.10.0GDSL esterase lipase
Unigene 827747.55.70.0GDSL esterase lipase
Unigene 1830986.60.90.0GDSL esterase lipase
Unigene 627560.03.20.1GDSL esterase lipase
Unigene 389111.40.00.0GDSL esterase lipase
Unigene 1132451.40.00.0GDSL esterase lipase
Unigene 23826812.80.50.0GDSL esterase lipase
Unigene 6967118.540.20.5GDSL esterase lipase APG
Unigene 1032985.311.00.0GDSL esterase lipase APG
Unigene 1526171.50.00.0Glycerol-3-phosphate acyltransferase
Unigene 1283681.10.00.0Glycerol-3-phosphate acyltransferase
Unigene 1632401.81.00.0Hydroxymethylglutaryl-CoA synthase
Unigene 12223071.42.060.9Inositol-3-phosphate synthase
Unigene 96290.70.00.0Linoleate 13S-lipoxygenase 2-1
Unigene 667169.10.10.0Linoleate 13S-lipoxygenase 2-1
Unigene 545201.30.00.0Linoleate 13S-lipoxygenase 2-1
Unigene 5992820.70.10.0Linoleate 13S-lipoxygenase 2-1
Unigene 2135895.50.00.0Linoleate 13S-lipoxygenase 2-1
Unigene 3758419.10.30.0Linoleate 13S-lipoxygenase 2-1
Unigene 3426715.00.10.0Linoleate 13S-lipoxygenase 2-1
Unigene 134712.40.10.0Linoleate 9S-lipoxygenase A
Unigene 54040.90.00.0Linoleate 9S-lipoxygenase A
Unigene 1397520.50.00.5Lipoxygenase 4
Unigene 2278694.611.80.1Phosphatidylinositol 4-kinase alpha 1
Unigene 17174672.595.50.6Phospholipase A2-alpha
Unigene 562850.00.00.8Phospholipid hydroperoxide
Unigene 6366924.820.00.5Pyruvate kinase 2
Unigene 636752.42.40.0Pyruvate kinase 2
Unigene 868312.81.00.0Pyruvate kinase 2
Unigene 1583474.210.000.26Non-specific phospholipase C6; EC
Unigene 1904395.876.860.08Non-specific phospholipase C2; EC
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shi, R.; Bai, H.; Li, B.; Liu, C.; Ying, Z.; Xiong, Z.; Wang, W. Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts. Life 2021, 11, 1431. https://doi.org/10.3390/life11121431

AMA Style

Shi R, Bai H, Li B, Liu C, Ying Z, Xiong Z, Wang W. Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts. Life. 2021; 11(12):1431. https://doi.org/10.3390/life11121431

Chicago/Turabian Style

Shi, Rui, Haidong Bai, Biao Li, Can Liu, Zhiping Ying, Zhi Xiong, and Wenlin Wang. 2021. "Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts" Life 11, no. 12: 1431. https://doi.org/10.3390/life11121431

APA Style

Shi, R., Bai, H., Li, B., Liu, C., Ying, Z., Xiong, Z., & Wang, W. (2021). Combined Transcriptome and Lipidomic Analyses of Lipid Biosynthesis in Macadamia ternifolia Nuts. Life, 11(12), 1431. https://doi.org/10.3390/life11121431

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop