Integrated Analysis of lncRNA–mRNA Regulatory Networks Related to Lipid Metabolism in High-Oleic-Acid Rapeseed

A high oleic acid content is considered an essential characteristic in the breeding of high-quality rapeseed in China. Long-chain non-coding RNA (lncRNA) molecules play an important role in the plant’s growth and its response to stress. To better understand the role of lncRNAs in regulating plant reproductive development, we analyzed whole-transcriptome and physiological data to characterize the dynamic changes in lncRNA expression during the four representative times of seed development of high- and low-oleic-acid rapeseed in three regions. We identified 21 and 14 lncRNA and mRNA modules, respectively. These modules were divided into three types related to region, development stages, and material. Next, we analyzed the key modules related to the oil content and the oleic acid, linoleic acid, and linolenic acid contents with physiological data and constructed the key functional network analysis on this basis. Genes related to lipid metabolism, such as 3-ketoacyl-CoA synthase 16 (KCS16) and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1), were present in the co-expression network, suggesting that the effect of these genes on lipid metabolism might be embodied by the expression of these lncRNAs. Our results provide a fresh insight into region-, development-stage-, and material-biased changes in lncRNA expression in the seeds of Brassica napus. Some of these lncRNAs may participate in the regulatory network of lipid accumulation and metabolism, together with regulated genes. These results may help elucidate the regulatory system of lncRNAs in the lipid metabolism of high-oleic-acid rapeseed seeds.

[1]  Chuanyi Yao,et al.  Computationally Guided Enzymatic Studies on Schizochytrium-Sourced Malonyl-CoA:ACP Transacylase. , 2022, Journal of agricultural and food chemistry.

[2]  M. Guan,et al.  Comprehensive evaluation of high-oleic rapeseed (Brassica napus) based on quality, resistance, and yield traits: A new method for rapid identification of high-oleic acid rapeseed germplasm , 2022, PLoS ONE.

[3]  Boas Pucker,et al.  Mapping-by-Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus , 2022, bioRxiv.

[4]  Hui Huang,et al.  Differences in Fatty Acids and Lipids of Massive and Branching Reef-Building Corals and Response to Environmental Changes , 2022, Frontiers in Marine Science.

[5]  R. Oelmüller,et al.  Genome-Wide Identification and Expression Profiling of KCS Gene Family in Passion Fruit (Passiflora edulis) Under Fusarium kyushuense and Drought Stress Conditions , 2022, Frontiers in Plant Science.

[6]  Nan Xu,et al.  Identification and Structure Analysis of KCS Family Genes Suggest Their Reponding to Regulate Fiber Development in Long-Staple Cotton Under Salt-Alkaline Stress , 2022, Frontiers in Genetics.

[7]  B. J. Harris,et al.  Simulated Docking Predicts Putative Channels for the Transport of Long-Chain Fatty Acids in Vibrio cholerae , 2022, bioRxiv.

[8]  R. Sandberg,et al.  Transcriptional kinetics and molecular functions of long non-coding RNAs , 2020, bioRxiv.

[9]  M. Lodha,et al.  Functions of long non-coding RNA in Arabidopsis thaliana , 2021, Plant signaling & behavior.

[10]  Hao Lin,et al.  Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage , 2021, Horticulture Research.

[11]  P. LaFayette,et al.  Effects of type I Diacylglycerol O-acyltransferase (DGAT1) genes on soybean (Glycine max L.) seed composition , 2021, Scientific Reports.

[12]  Xin Feng,et al.  Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening , 2021, Scientific Reports.

[13]  N. Kav,et al.  Identification of lncRNAs in response to infection by Plasmodiophora brassicae in Brassica napus and development of lncRNA-based SSR markers. , 2020, Genome.

[14]  K. Kaufmann,et al.  Long non-coding RNAs in plants: emerging modulators of gene activity in development and stress responses , 2020, Planta.

[15]  K. Mysore,et al.  The 3-ketoacyl-CoA synthase WFL is involved in lateral organ development and cuticular wax synthesis in Medicago truncatula , 2020, Plant molecular biology.

[16]  Gao Caixia,et al.  Roles of lncRNAs in Rice: Advances and Challenges , 2020 .

[17]  Xueping Zhou,et al.  Transcriptome Analysis of Rice Reveals the lncRNA–mRNA Regulatory Network in Response to Rice Black-Streaked Dwarf Virus Infection , 2020, Viruses.

[18]  Yingfen Jiang,et al.  Association of Clubroot Resistance Locus PbBa8.1 With a Linkage Drag of High Erucic Acid Content in the Seed of the European Turnip , 2020, Frontiers in Plant Science.

[19]  S. Spasibionek,et al.  Marker assisted selection of new high oleic and low linolenic winter oilseed rape (Brassica napus L.) inbred lines revealing good agricultural value , 2020, PloS one.

[20]  F. Ariel,et al.  The Arabidopsis lncRNA ASCO modulates the transcriptome through interaction with splicing factors , 2020, EMBO reports.

[21]  S. Singer,et al.  The Role of Triacylglycerol in Plant Stress Response , 2020, Plants.

[22]  Liping Wang,et al.  On the Role of DGAT1 in Seed Glycerolipid Metabolic Network and Critical Stages of Plant Development in Arabidopsis. , 2020, Lipids.

[23]  Liyong Hu,et al.  Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering , 2020, BMC Plant Biology.

[24]  A. Żbikowska,et al.  Effect of high-oleic rapeseed oil oleogels on the quality of short-dough biscuits and fat migration , 2019, Journal of Food Science and Technology.

[25]  P. Seo,et al.  Interaction of DGAT1 and PDAT1 to enhance TAG assembly in Arabidopsis , 2018, Plant signaling & behavior.

[26]  R. D. Lawrence,et al.  Metabolic regulation , 2019, Prokaryotic Metabolism and Physiology.

[27]  Qing Liu,et al.  Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus , 2018, BMC Genomics.

[28]  Guangsheng Zhou,et al.  Higher density planting benefits mechanical harvesting of rapeseed in the Yangtze River Basin of China , 2018 .

[29]  Jia Gu,et al.  fastp: an ultra-fast all-in-one FASTQ preprocessor , 2018, bioRxiv.

[30]  Chaofu Lu,et al.  Mutagenesis of the FAE1 genes significantly changes fatty acid composition in seeds of Camelina sativa. , 2018, Plant physiology and biochemistry : PPB.

[31]  P. Khavari,et al.  The functions and unique features of long intergenic non-coding RNA , 2017, Nature Reviews Molecular Cell Biology.

[32]  R. Jetter,et al.  Arabidopsis ketoacyl-CoA synthase 16 (KCS16) forms C36 /C38 acyl precursors for leaf trichome and pavement surface wax. , 2017, Plant, cell & environment.

[33]  Dawei Li,et al.  Whole transcriptome sequencing of Pseudomonas syringae pv. actinidiae-infected kiwifruit plants reveals species-specific interaction between long non-coding RNA and coding genes , 2017, Scientific Reports.

[34]  Li Wang,et al.  Identification and characterization of novel lncRNAs in Arabidopsis thaliana. , 2017, Biochemical and biophysical research communications.

[35]  T. Sang,et al.  Systematic comparison of lncRNAs with protein coding mRNAs in population expression and their response to environmental change , 2017, BMC Plant Biology.

[36]  Z. Yang,et al.  Characterization of long non-coding RNAs involved in cadmium toxic response in Brassica napus , 2016 .

[37]  Swati Megha,et al.  Genome Wide Identification and Functional Prediction of Long Non-Coding RNAs Responsive to Sclerotinia sclerotiorum Infection in Brassica napus , 2016, PloS one.

[38]  Fenghong Huang,et al.  A Study on Triacylglycerol Composition and the Structure of High-Oleic Rapeseed Oil , 2016 .

[39]  M. Özcan,et al.  Some rape/canola seed oils: fatty acid composition and tocopherols , 2016, Zeitschrift fur Naturforschung. C, Journal of biosciences.

[40]  Mukhlesur Rahman,et al.  Chapter 15 – Designer Oil Crops , 2016 .

[41]  Steven L Salzberg,et al.  HISAT: a fast spliced aligner with low memory requirements , 2015, Nature Methods.

[42]  S. Salzberg,et al.  StringTie enables improved reconstruction of a transcriptome from RNA-seq reads , 2015, Nature Biotechnology.

[43]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[44]  Tingting Guo,et al.  Selection and evaluation of novel reference genes for quantitative reverse transcription PCR (qRT-PCR) based on genome and transcriptome data in Brassica napus L. , 2014, Gene.

[45]  Y. Hang,et al.  Evolutionary Pattern of the FAE1 Gene in Brassicaceae and Its Correlation with the Erucic Acid Trait , 2013, PloS one.

[46]  Yi Zhao,et al.  Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts , 2013, Nucleic acids research.

[47]  J. Ohlrogge,et al.  Biochemical pathways in seed oil synthesis. , 2013, Current opinion in plant biology.

[48]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[49]  J. Browse,et al.  The pathway of triacylglycerol synthesis through phosphatidylcholine in Arabidopsis produces a bottleneck for the accumulation of unusual fatty acids in transgenic seeds. , 2011, The Plant journal : for cell and molecular biology.

[50]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[51]  J. Ohlrogge,et al.  DGAT1 and PDAT1 Acyltransferases Have Overlapping Functions in Arabidopsis Triacylglycerol Biosynthesis and Are Essential for Normal Pollen and Seed Development[W][OA] , 2009, The Plant Cell Online.

[52]  Steve Horvath,et al.  WGCNA: an R package for weighted correlation network analysis , 2008, BMC Bioinformatics.

[53]  P. Moreau,et al.  The VLCFA elongase gene family in Arabidopsis thaliana: phylogenetic analysis, 3D modelling and expression profiling , 2008, Plant Molecular Biology.

[54]  Yong Zhang,et al.  CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine , 2007, Nucleic Acids Res..

[55]  J. Jaworski,et al.  Substrate specificity of Arabidopsis 3-ketoacyl-CoA synthases. , 2006, Biochemical and biophysical research communications.

[56]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[57]  M. Pfaffl,et al.  A new mathematical model for relative quantification in real-time RT-PCR. , 2001, Nucleic acids research.

[58]  U. Grossniklaus,et al.  FIDDLEHEAD, a gene required to suppress epidermal cell interactions in Arabidopsis, encodes a putative lipid biosynthetic enzyme. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[59]  D. Taylor,et al.  CUT1, an Arabidopsis Gene Required for Cuticular Wax Biosynthesis and Pollen Fertility, Encodes a Very-Long-Chain Fatty Acid Condensing Enzyme , 1999, Plant Cell.

[60]  K. Krolikowski,et al.  Developmental regulation of cell interactions in the Arabidopsis fiddlehead-1 mutant: a role for the epidermal cell wall and cuticle. , 1997, Developmental biology.

[61]  D. W. James,et al.  Directed tagging of the Arabidopsis FATTY ACID ELONGATION1 (FAE1) gene with the maize transposon activator. , 1995, The Plant cell.

[62]  L Hartman,et al.  Rapid preparation of fatty acid methyl esters from lipids. , 1973, Laboratory practice.