Genome-wide association analysis of the lipid and fatty acid metabolism regulatory network in the mesocarp of oil palm (Elaeis guineensis Jacq.) based on small noncoding RNA sequencing

Integrated analysis of miRNA, transcriptome, and degradome sequencing provides new insights into lipid metabolism in perilla seed

MicroRNAs (miRNA) are small noncoding RNAs that play a crucial as molecular regulators in lipid metabolism in various oil crops. Perilla (Perilla frutescens) is a specific oil crop known for its high alpha-linolenic acid (C18:3n3, ALA) content (>65 %) in their seed oils. In view of the regulatory mechanism of miRNAs in perilla remains unclear, we conducted miRNAs and transcriptome sequencing in two cultivars with distinct lipid compositions. A total of 525 unique miRNAs, including 142 differentially expressed miRNAs was identified in perilla seeds. The 318 miRNAs targeted 7,761 genes. Furthermore, we identified 112 regulated miRNAs and their 610 target genes involved in lipid metabolism. MiR159b and miR167a as the core nodes to regulate the expression of genes in oil biosynthesis (e.g., KAS, FATB, GPAT, FAD, DGK, LPAAT) and key regulatory TFs (e.g., MYB, ARF, DOF, SPL, NAC, TCP, and bHLH). The 1,219 miRNA-mRNA regulation modules were confirmed through degradome sequencing. Notably, pf-miR159b-MYBs and pf-miR167a-ARFs regulation modules were confirmed. They exhibited significantly different expression levels in two cultivars and believed to play important roles in oil biosynthesis in perilla seeds. This provides valuable insights into the functional analysis of miRNA-regulated lipid metabolism in perilla seeds.

Genome-Wide Identification, Evolution and Expression Profile Analysis of NAC Transcription Factor in Simmondsia chinensis

NAC transcription factors (TFs) are one of the largest plant-specific gene families and play important roles in plant growth, development, and the biotic and abiotic stress response. Although the sequencing of Jojoba (Simmondsia chinensis) has been completed, the genome-wide identification and analysis of its NAC TFs has not been reported. In this study, a total of 57 genes were identified in Jojoba, which were divided into eight groups based on phylogenetic analysis. The genes clustered in the same groups have a similar gene structure and motif distribution. Based on the analysis of cis-elements in NAC TFs, nine cis-acting elements were identified in the promoter region that involved in light response, hormonal response, and stress response. Synteny analysis showed a greater collinearity between Jojoba and V. vinifera than Arabidopsis thaliana. The 24 genes in the Jojoba NAC TFs are derived from fragment replication, which may be the main source of NAC amplification. Gene expression analysis identified seven genes that were highly expressed in seeds. The differential expression analysis of NAC TFs in cotyledon and embryonic axis tissues showed that the expression of 10 genes was up-regulated and 1 gene was down-regulated. This study provides more information on the classification, gene structure, conserved motif, and evolution of NAC TFs in Jojoba, facilitating further exploration of their specific functional analysis in Jojoba seed development.

The early response of oil palm (Elaeis guineensis Jacq.) plants to water deprivation: Expression analysis of miRNAs and their putative target genes, and similarities with the response to salinity stress

Oil palm (Elaeis guineensis Jacq.) is a oilseed crop of great economic importance drastically affected by abiotic stresses. MicroRNAs (miRNAs) play crucial roles in transcription and post-transcription regulation of gene expression, being essential molecules in the response of plants to abiotic stress. To better understand the molecular mechanisms behind the response of young oil palm plants to drought stress, this study reports on the prediction and characterization of miRNAs and their putative target genes in the apical leaf of plants subjected to 14 days of water deprivation. Then, the data from this study were compared to the data from a similar study that focused on salinity stress. Both, the drought-and salt-responsive miRNAs and their putative target genes underwent correlation analysis to identify similarities and dissimilarities among them. Among the 81 identified miRNAs, 29 are specific for oil palm, including two (egu-miR28ds and egu-miR29ds) new ones - described for the first time. As for the expression profile, 62 miRNAs were significantly differentially expressed under drought stress, being five up-regulated (miR396e, miR159b, miR529b, egu-miR19sds, and egu-miR29ds) and 57 down-regulated. Transcription factors, such as MYBs, HOXs, and NF-Ys, were predicted as putative miRNA-target genes in oil palm under water deprivation; making them the most predominant group of such genes. Finally, the correlation analysis study revealed a group of putative target genes with similar behavior under salt and drought stresses. Those genes that are upregulated by these two abiotic stresses encode lncRNAs and proteins linked to stress tolerance, stress memory, modulation of ROS signaling, and defense response regulation to abiotic and biotic stresses. In summary, this study provides molecular evidence for the possible involvement of miRNAs in the drought stress response in oil palm. Besides, it shows that, at the molecular level, there are many similarities in the response of young oil palm plants to these two abiotic stresses.

Characterization and functional analysis of the MADS-box EgAGL9 transcription factor from the mesocarp of oil palm (Elaeis guineensis Jacq.)

Oil palm (Elaeis guineensis Jacq.) is one of the most important oil crops in the world, and compared to all oil crops, it has the highest productive efficiency. In the present study, a MADS-box transcription factor of the AGAMOUS class, named EgAGL9, was identified by expression profile analysis in the different developmental stages of oil palm mesocarp. Real-time quantitative PCR results confirmed that the expression of EgAGL9 increased rapidly during the last stages of oil palm mesocarp development. Then, three downstream genes, including EgSAD (Stearoyl-ACP desaturase), EgTSA (Tryptophan synthase) and EgSDH (Succinate dehydrogenase), were screened by ChIP-Seq and data analysis. EMSA analysis verified that EgAGL9 interacted with the promoter regions of EgSAD, EgTSA and EgSDH. Moreover, the expression levels of EgSAD, EgTSA and EgSDH were downregulated in EgAGL9-overexpressing protoplasts and calli of oil palm. Compared to WT, the total lipid content and ratio of unsaturated fatty acids in transgenic calli (including oleic acid, linoleic acid and linolenic acid) were significantly decreased. Together, these results revealed that these three EgAGL9-regulated genes are involved in regulatory pathways in the oil palm mesocarp. Compared with previous studies, the present study provides a new research strategy for understanding of the molecular regulatory pathways of lipid metabolism in mesocarp of oil palm. The obtained results will bring a new perspective for a comprehensive understanding of the regulation of the metabolic accumulation in the oil palm mesocarp.

EgMYB108 regulates very long-chain fatty acid (VLCFA) anabolism in the mesocarp of oil palm

EgMYB108 regulates VLCFA anabolism in oil palm. Very long-chain fatty acids (VLCFAs), which are fatty acids with more than 18 C, can not only be used as a form of triglyceride (TAG) but also provide precursors for the biosynthesis of cuticle wax, and they exist in plant epidermal cells in the form of wax in higher plants. However, which and how transcriptional factors (TFs) regulate this process is largely unknown in oil palm. In this study, a MYB transcription factor (EgMYB108) with high expression in the mesocarp of oil palm fruit was characterized. Overexpression of EgMYB108 promoted not only total lipid content but also VLCFA accumulation in oil palm embryoids. Subsequently, transient transformation in protoplasts and qRT-PCR analysis indicated that the EgKCS5 and EgLACS4 genes were significantly increased with the overexpression of EgMYB108. Furthermore, yeast one‑hybrid assays, dual-luciferase assays and EMSAs demonstrated that EgMYB108 binds to the promoters of EgKCS5 and EgLACS4 and regulates their transcription. Finally, EgMYB108 interacts with the promoters of EgLACS and EgKCS simultaneously and finally improves the VLCFA and total lipid contents; a pathway summarizing this interaction was depicted.. The results provide new insight into the mechanism by which EgMYB108 regulates lipid and VLCFA accumulation in oil palm.

The reference genome of camellia chekiangoleosa provides insights into camellia evolution and tea oil biosynthesis

Camellia oil extracted from Camellia seeds is rich in unsaturated fatty acids (UFAs) and secondary metabolites beneficial to human health. However, no oil-tea tree genome has yet been published, which is a major obstacle to investigating the heredity improvement of oil-tea trees. Here, using both Illumina and PicBio sequencing technologies, we present the first chromosome-level genome sequence of the oil-tea tree species Camellia chekiangoleosa Hu. (CCH). The assembled genome consists of 15 pseudochromosomes with a genome size of 2.73 Gb and a scaffold N50 of 185.30 Mb. At least 2.16 Gb of the genome assembly consists of repetitive sequences, and the rest involves a high-confidence set of 64 608 protein-coding gene models. Comparative genomic analysis revealed that the CCH genome underwent a whole-genome duplication (WGD) event shared across the Camellia genus at ~57.48 MYA and a γ-WGT event shared across all core eudicot plants at ~120 MYA. Gene family clustering revealed that the genes involved in terpenoid biosynthesis have undergone rapid expansion. Furthermore, we determined the expression patterns of oleic acid accumulation- and terpenoid biosynthesis-associated genes in six tissues. We found that these genes tend to be highly expressed in leaves, pericarp tissues, roots, and seeds. The first chromosome-level genome of oil-tea trees will provide valuable resources for determining Camellia evolution and utilizing the germplasm of this taxon.

Genomic loci mispositioning in Tmem120a knockout mice yields latent lipodystrophy

Little is known about how the observed fat-specific pattern of 3D-spatial genome organisation is established. Here we report that adipocyte-specific knockout of the gene encoding nuclear envelope transmembrane protein Tmem120a disrupts fat genome organisation, thus causing a lipodystrophy syndrome. Tmem120a deficiency broadly suppresses lipid metabolism pathway gene expression and induces myogenic gene expression by repositioning genes, enhancers and miRNA-encoding loci between the nuclear periphery and interior. Tmem120a-/- mice, particularly females, exhibit a lipodystrophy syndrome similar to human familial partial lipodystrophy FPLD2, with profound insulin resistance and metabolic defects that manifest upon exposure to an obesogenic diet. Interestingly, similar genome organisation defects occurred in cells from FPLD2 patients that harbour nuclear envelope protein encoding LMNA mutations. Our data indicate TMEM120A genome organisation functions affect many adipose functions and its loss may yield adiposity spectrum disorders, including a miRNA-based mechanism that could explain muscle hypertrophy in human lipodystrophy.

Crosstalk Between Polygonatum kingianum, the miRNA, and Gut Microbiota in the Regulation of Lipid Metabolism

Objectives: Polygonatum kingianum is a medicinal herb used in various traditional Chinese medicine formulations. The polysaccharide fraction of P. kingianum can reduce insulin resistance and restore the gut microbiota in a rat model of aberrant lipid metabolism by down regulating miR-122. The aim of this study was to further elucidate the effect of P. kingianum on lipid metabolism, and the roles of specific miRNAs and the gut microbiota. Key findings: P. kingianum administration significantly altered the abundance of 29 gut microbes and 27 differentially expressed miRNAs (DEMs). Several aberrantly expressed miRNAs closely related to lipid metabolism were identified, of which some were associated with specific gut microbiota. MiR-484 in particular was identified as the core factor involved in the therapeutic effects of P. kingianum. We hypothesize that the miR-484-Bacteroides/Roseburia axis acts as an important bridge hub that connects the entire miRNA-gut microbiota network. In addition, we observed that Parabacteroides and Bacillus correlated significantly with several miRNAs, including miR-484, miR-122-5p, miR-184 and miR-378b. Summary: P. kingianum alleviates lipid metabolism disorder by targeting the network of key miRNAs and the gut microbiota.

Expression analysis of miRNAs and their putative target genes confirm a preponderant role of transcription factors in the early response of oil palm plants to salinity stress

BACKGROUND

Several mechanisms regulating gene expression contribute to restore and reestablish cellular homeostasis so that plants can adapt and survive in adverse situations. MicroRNAs (miRNAs) play roles important in the transcriptional and post-transcriptional regulation of gene expression, emerging as a regulatory molecule key in the responses to plant stress, such as cold, heat, drought, and salt. This work is a comprehensive and large-scale miRNA analysis performed to characterize the miRNA population present in oil palm (Elaeis guineensis Jacq.) exposed to a high level of salt stress, to identify miRNA-putative target genes in the oil palm genome, and to perform an in silico comparison of the expression profile of the miRNAs and their putative target genes.

RESULTS

A group of 79 miRNAs was found in oil palm, been 52 known miRNAs and 27 new ones. The known miRNAs found belonged to 28 families. Those miRNAs led to 229 distinct miRNA-putative target genes identified in the genome of oil palm. miRNAs and putative target genes differentially expressed under salinity stress were then selected for functional annotation analysis. The regulation of transcription, DNA-templated, and the oxidation-reduction process were the biological processes with the highest number of hits to the putative target genes, while protein binding and DNA binding were the molecular functions with the highest number of hits. Finally, the nucleus was the cellular component with the highest number of hits. The functional annotation of the putative target genes differentially expressed under salinity stress showed several ones coding for transcription factors which have already proven able to result in tolerance to salinity stress by overexpression or knockout in other plant species.

CONCLUSIONS

Our findings provide new insights into the early response of young oil palm plants to salinity stress and confirm an expected preponderant role of transcription factors - such as NF-YA3, HOX32, and GRF1 - in this response. Besides, it points out potential salt-responsive miRNAs and miRNA-putative target genes that one can utilize to develop oil palm plants tolerant to salinity stress.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Plant miR397 and its functions

MicroRNAs (miRNAs) are noncoding, small RNAs of 20-24 nucleotides (nt) and function critically at the post-transcriptional level to regulate gene expression through cleaving mRNA targets or interfering with translation of the target mRNAs. They are broadly involved in many biological processes in plants. The miR397 family in plants contains several conserved members either in 21-nt or in 22-nt that mainly target the laccase (LAC) genes functioning in lignin synthesis and are involved in the development of plants under various conditions. Recent findings showed that miR397b in Arabidopsis could also target to Casein Kinase II Subunit Beta 3 (CKB3) and mediate circadian regulation and plant flowering. This review aims to summarise recent updates on miR397 and provides the available basis for understanding the functional mechanisms of miR397 in plant growth and development regulation and in response to external adverse stimulation.

Integration of miRNAs, Degradome, and Transcriptome Omics Uncovers a Complex Regulatory Network and Provides Insights Into Lipid and Fatty Acid Synthesis During Sesame Seed Development

Sesame (Sesamum indicum L.) has always been known as a health-promoting oilseed crop because of its nutrient-rich oil. In recent years, studies have focused on lipid and fatty acid (FA) biosynthesis in various plants by high-throughput sequencing. Here, we integrated transcriptomics, small RNAs, and the degradome to establish a comprehensive reserve intensive on key regulatory micro RNA (miRNA)-targeting circuits to better understand the transcriptional and translational regulation of the oil biosynthesis mechanism in sesame seed development. Deep sequencing was performed to differentially express 220 miRNAs, including 65 novel miRNAs, in different developmental periods of seeds. GO and integrated KEGG analysis revealed 32 pairs of miRNA targets with negatively correlated expression profiles, of which 12 miRNA-target pairs were further confirmed by RT-PCR. In addition, a regulatory co-expression network was constructed based on the differentially expressed gene (DEG) profiles. The FAD2, LOC10515945, LOC105161564, and LOC105162196 genes were clustered into groups that regulate the accumulation of unsaturated fatty acid (UFA) biosynthesis. The results provide a unique advanced molecular platform for the study of lipid and FA biosynthesis, and this study may serve as a new theoretical reference to obtain increased levels of UFA from higher-quality sesame seed cultivars and other plants.

Functional Role of Long-Chain Acyl-CoA Synthetases in Plant Development and Stress Responses

Fatty acids (FAs) play vital roles in plants as components of lipid membranes that demarcate cells and organelles, as sources of stored energy in the form of neutral lipids, and as signaling molecules that elicit plant responses to adverse conditions. The activation of FAs through the formation of acyl-CoA intermediates by acyl-CoA synthetase (ACS) family enzymes is required for their synthesis and degradation. Long-chain ACSs (LACSs) represent a small subgroup of ACS enzymes that specifically convert long-chain or very-long-chain FAs into corresponding thioesters for multiple lipid-associated processes. Alteration of LACS activity often results in pleiotropic phenotypes such as male sterility, organ fusion, aberrant cuticular structure, delayed seed germination, altered seed oil content, and plant capacity to respond to various environmental stresses. This review provides a comprehensive analysis of LACS family enzymes including substrate specificity, tissue-specific expression patterns, and distinct subcellular localization highlighting their specific roles in lipid synthesis and degradation, the effects of altered LACS activity on plant development, the relationship between LACS activity and stress resistance, and the regulation of LACS activity. Finally, we pose several major questions to be addressed, which would advance our current understanding of LACS function in plants.

Integrated Analysis of Seed microRNA and mRNA Transcriptome Reveals Important Functional Genes and microRNA-Targets in the Process of Walnut (Juglans regia) Seed Oil Accumulation

Walnut (Juglans regia) is known as a promising woody oil crop with abundant polyunsaturated fatty acids in its kernel. However, the regulation mechanism of walnut oil accumulation and fatty acid metabolism is still poorly understood, which restricted the breeding and genetic improvement of high-quality oil-bearing walnuts. To reveal the molecular mechanism of walnut oil accumulation, considering the potential regulation of microRNA (miRNA) in seed development, in this study, the oil content of walnut kernel on the 80th, 100th and 120th day after flowering (DAF) was tested and the corresponding proportions are 11.51%, 40.40% and 53.20%. Between DAF of 80th~120th, the content of stearic acid and oleic acid tended to increase, but the proportion of other fatty acids tended to decrease. Meanwhile, comparative transcriptome and sRNA-seq analysis on three stages (80th, 100th and 120th DAF), found 204 conserved miRNAs and 554 novel miRNAs in walnut kernels, among which 104 key genes related to walnut oil accumulation were screened. The phospholipid:diacylglycerol acyltransferase metabolic pathway may contribute more to oil accumulation in walnut. 16 miRNA-mRNA regulatory modules related to walnut oil accumulation and fatty acid synthesis were constructed. 8 known miRNAs and 9 novel miRNAs regulate 28 genes involved in fatty acid (FA) metabolism and lipid synthesis. Among them, jre-miRn105, jre-miRn434, jre-miR477d and jre-miR156a.2 are key miRNAs that regulate walnut FA synthesis. Jre-miRn411 and jre-miR399a.1 are closely related to oil accumulation. These data provide new insights and lay the foundation for subsequent studies on walnut FA synthesis and oil accumulation.

A MADS-box gene, EgMADS21, negatively regulates EgDGAT2 expression and decreases polyunsaturated fatty acid accumulation in oil palm (Elaeis guineensis Jacq.)

EgMADS21 regulates PUFA accumulation in oil palm. Oil palm (Elaeis guineensis Jacq.) is the most productive world oil crop, accounting for 36% of world plant oil production. However, the molecular mechanism of the transcriptional regulation of fatty acid accumulation and lipid synthesis in the mesocarp of oil palm by up- or downregulating the expression of genes involved in related pathways remains largely unknown. Here, an oil palm MADS-box gene, EgMADS21, was screened in a yeast one-hybrid assay using the EgDGAT2 promoter sequence as bait. EgMADS21 is preferentially expressed in early mesocarp developmental stages in oil palm fruit and presents a negative correlation with EgDGAT2 expression. The direct binding of EgMADS21 to the EgDGAT2 promoter was confirmed by electrophoretic mobility shift assay. Subsequently, transient expression of EgMADS21 in oil palm protoplasts revealed that EgMADS21 not only binds to the EgDGAT2 promoter but also negatively regulates the expression of EgDGAT2. Furthermore, EgMADS21 was stably overexpressed in transgenic oil palm embryoids by Agrobacterium-mediated transformation. In three independent transgenic lines, EgDGAT2 expression was significantly suppressed by the expression of EgMADS21. The content of linoleic acid (C18:2) in the three transgenic embryoids was significantly decreased, while that of oleic acid (C18:1) was significantly increased. Combined with the substrate preference of EgDGAT2 identified in previous research, the results demonstrate the molecular mechanism by which EgMADS21 regulates EgDGAT2 expression and ultimately affects fatty acid accumulation in the mesocarp of oil palm.

Identification and Analysis of microRNAs in Chlorella sorokiniana Using High-Throughput Sequencing

Chlorella is a popular microalga with robust physiological and biochemical characteristics, which can be cultured under various conditions. The exploration of the small RNA content of Chlorella could improve strategies for the enhancement of metabolite production from this microalga. In this study, stress was introduced to the Chlorella sorokiniana culture to produce high-value metabolites such as carotenoids and phenolic content. The small RNA transcriptome of C. sorokiniana was sequenced, focusing on microRNA (miRNA) content. From the analysis, 98 miRNAs were identified in cultures subjected to normal and stress conditions. The functional analysis result showed that the miRNA targets found were most often involved in the biosynthesis of secondary metabolites, followed by protein metabolism, cell cycle, and porphyrin and chlorophyll metabolism. Furthermore, the biosynthesis of secondary metabolites such as carotenoids, terpenoids, and lipids was found mostly in stress conditions. These results may help to improve our understanding of regulatory mechanisms of miRNA in the biological and metabolic process of Chlorella species. It is important and timely to determine the true potential of this microalga species and to support the potential for genetic engineering of microalgae as they receive increasing focus for their development as an alternative source of biofuel, food, and health supplements.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.