Genome-wide identification and analysis of long non-coding RNAs involved in fatty acid biosynthesis in young soybean pods

Long non-coding RNAs (lncRNAs) are non-coding RNAs of more than 200 nucleotides. To date, the roles of lncRNAs in soybean fatty acid synthesis have not been fully studied. Here, the low-linolenic acid mutant ‘MT72′ and the wild-type control ‘JN18′ were used as materials. The lncRNAs in young pods at 30 and 40 days (d) after flowering were systematically identified and analyzed using transcriptome sequencing technology combined with bioinformatics tools. A total of 39,324 lncRNAs and 561 differentially expressed lncRNAs were identified. A lncRNAs-miRNAs-protein-coding genes (mRNAs) network was constructed, and 46 lncRNAs, 46 miRNAs and 137 mRNAs were found to be correlated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of 12 targeted mRNAs in the competing endogenous RNA network showed that these lncRNAs may be involved in the biological processes of fatty acid transport, lipid synthesis and cell division. Finally, the expression levels of differentially expressed lncRNAs, miRNAs and mRNAs were verified using qRT-PCR. The expression patterns of most genes were consistent with the sequencing results. In conclusion, new information was provided for the study of fatty acid synthesis by lncRNAs in young soybean pods.