Kelch-motif containing acyl-CoA binding proteins AtACBP4 and AtACBP5 are differentially expressed and function in floral lipid metabolism

A promoter polymorphism defines distinct roles in anther development for Col-0 and Ler-0 alleles of Arabidopsis ACYL-COA BINDING PROTEIN3

Acyl-CoA-Binding Proteins (ACBPs) bind acyl-CoA esters and function in lipid metabolism. Although acbp3-1, the ACBP3 mutant in Arabidopsis thaliana ecotype Col-0, displays normal floral development, the acbp3-2 mutant from ecotype Ler-0 characterized herein exhibits defective adaxial anther lobes and improper sporocyte formation. To understand these differences and identify the role of ERECTA in ACBP3 function, the acbp3 mutants and acbp3-erecta (er) lines were analyzed by microscopy for anther morphology and high-performance liquid chromatography for lipid composition. Defects in Landsberg anther development were related to the ERECTA-mediated pathway because the progenies of acbp3-2 × La-0 and acbp3-1 × er-1 in Col-0 showed normal anthers, contrasting to that of acbp3-2 in Ler-0. Polymorphism in the regulatory region of ACBP3 enabled its function in anther development in Ler-0 but not Col-0 which harbored an AT-repeat insertion. ACBP3 expression and anther development in acbp3-2 were restored using ACBP3pro (Ler)::ACBP3 not ACBP3pro (Col)::ACBP3. SPOROCYTELESS (SPL), a sporocyte formation regulator activated ACBP3 transcription in Ler-0 but not Col-0. For anther development, the ERECTA-related role of ACBP3 is required in Ler-0, but not Col-0. The disrupted promoter regulatory region for SPL binding in Col-0 eliminates the role of ACBP3 in anther development.

Characterization and expression analysis of chalcone synthase gene family members suggested their roles in the male sterility of a wheat temperature-sensitive genic male sterile (TGMS) line

Chalcone synthase (CHS), also known as the plants-specific type III polyketide synthases (PKSs), catalyzes the first key step in the biosynthesis of plant flavonoids. Flavonoids are one of the most important secondary metabolites which participate in flower pigmentation and pollen fertility. Recent reports have demonstrated the role of the CHS family in plant pollen exine formation. This study focused on the potential roles of CHS in the pollen exine formation of wheat. In the present study, a genome-wide investigation of the CHS family was carried out, and 87 CHS genes were identified in wheat. TaCHS3, TaCHS10, and TaCHS13 are wheat orthologs of Arabidopsis LESS ADHESIVE POLLEN (LAP5); TaCHS58, TaCHS64, and TaCHS67 are wheat orthologs of AtLAP6. TaCHS3, TaCHS10, and TaCHS67 showed anther-specific patterns. The expression of TaCHS3, TaCHS10, and TaCHS67 was positively co-expressed with sporopollenin biosynthetic genes, including TaCYP703A2, TaCYP704B1, TaDRL1, TaTKPR2, and TaMS2. Coincidently, the expression of TaCHS3, TaCHS10, and TaCHS67, together with those sporopollenin biosynthetic genes, were repressed at the tetrads and uninucleate stages in the temperature-sensitive genic male-sterile (TGMS) line BS366 under sterile conditions. Wheat anther-specific CHS genes might participate in the exine formation of BS366 through co-expressing with sporopollenin biosynthetic genes, which will undoubtedly provide knowledge of the roles of CHS in wheat pollen development.

Roles of acyl-CoA-binding proteins in plant reproduction

Acyl-CoA-binding proteins (ACBPs) constitute a well-conserved family of proteins in eukaryotes that are important in stress responses and development. Past studies have shown that ACBPs are involved in maintaining, transporting and protecting acyl-CoA esters during lipid biosynthesis in plants, mammals, and yeast. ACBPs show differential expression and various binding affinities for acyl-CoA esters. Hence, ACBPs can play a crucial part in maintaining lipid homeostasis. This review summarizes the functions of ACBPs during the stages of reproduction in plants and other organisms. A comprehensive understanding on the roles of ACBPs during plant reproduction may lead to opportunities in crop improvement in agriculture.

Phylogeny and subcellular localization analyses reveal distinctions in monocot and eudicot class IV acyl-CoA-binding proteins

Plant class IV ACBPs diverged with the split of monocots and eudicots. Difference in the subcellular localization supported the functional variation of plant class IV ACBP. Acyl-CoA-binding proteins (ACBPs) are divided into class I-IV in plants. Class IV ACBPs are kelch motif containing proteins that are specific to plants. The currently known subcellular localizations of plant class IV ACBPs are either in the cytosol (Arabidopsis) or in the peroxisomes (rice). However, it is not clear whether peroxisomal localization of class IV ACBP is a shared character that distinguishes eudicots and monocots. Here, the phylogeny of class IV ACBPs from 73 plant species and subcellular localization of class IV ACBPs from six monocots and eudicots were conducted. Phylogenetic analysis of 112 orthologues revealed that monocot class IV ACBPs were basal to the monophyletic clade formed by eudicots and basal angiosperm. Transient expression of GFP fusions in onion epidermal cells demonstrated that monocot maize (Zea mays), wheat (Triticum aestivum), and sorghum (Sorghum bicolor) and eudicot poplar (Populus trichocarpa) all contained at least one peroxisomal localized class IV ACBP, while orthologues from cucumber (Cucumis sativus L.) and soybean (Glycine max) were all cytosolic. Combining the location of Arabidopsis and rice class IV ACBPs, it indicates that maintaining at least one peroxisomal class IV ACBP could be a shared feature within the tested monocots, while cytosolic class IV ACBPs would be preferred in the tested eudicots. Furthermore, the interaction between OsACBP6 and peroxisomal ATP-binding cassette (ABC) transporter provided clues for the functional mechanism of OsACBP6.

A kelch-repeat superfamily gene, ZmNL4, controls leaf width in maize (Zea mays L.)

Leaf width (LW) is an important component of plant architecture that extensively affects both light capture during photosynthesis and grain yield, particularly under dense planting conditions. However, the genetic and molecular mechanisms regulating LW remain largely elusive in maize (Zea mays L.). In this study, qLW4a, a major quantitative trait locus controlling LW, was identified in a population constructed with maize inbred lines PH6WC, with wide leaves, and Lin387, with narrow leaves. Map-based cloning revealed that ZmNL4, a kelch-repeat superfamily gene, emerged to be the candidate for qLW4a, and a single-base deletion in the conserved SMC_prok_B domain of ZmNL4 in Lin387 caused a frame shift, leading to premature termination. Consistently, the knockout of ZmNL4 by CRISPR/Cas9 editing significantly reduced the LW, which was attributed to a reduction in the cell number instead of cell size, indicating a role of ZmNL4 in regulating cell division. Transcriptomic comparison of ZmNL4 knockout lines with the wild type B73-329 revealed that ZmNL4 might participate in cell wall biogenesis, asymmetric cell division, metabolic processes, transmembrane transport and response to external stimulus, etc. These results provide insights into the genetic and molecular mechanisms of ZmNL4 in controlling LW and could potentially contribute to optimizing plant architecture for maize breeding.

CRISPR/Cas9-Mediated Knockout of GmFATB1 Significantly Reduced the Amount of Saturated Fatty Acids in Soybean Seeds

Soybean (Glycine max) oil is one of the most widely used vegetable oils across the world. Breeding of soybean to reduce the saturated fatty acid (FA) content, which is linked to cardiovascular disease, would be of great significance for nutritional improvement. Acyl-acyl carrier protein thioesterases (FATs) can release free FAs and acyl-ACP, which ultimately affects the FA profile. In this study, we identified a pair of soybean FATB coding genes, GmFATB1a and GmFATB1b. Mutants that knock out either or both of the GmFATB1 genes were obtained via CRISPR/Cas9. Single mutants, fatb1a and fatb1b, showed a decrease in leaf palmitic and stearic acid contents, ranging from 11% to 21%. The double mutant, fatb1a:1b, had a 42% and 35% decrease in palmitic and stearic acid content, displayed growth defects, and were male sterility. Analysis of the seed oil profile revealed that fatb1a and fatb1b had significant lower palmitic and stearic acid contents, 39-53% and 17-37%, respectively, while that of the unsaturated FAs were the same. The relative content of the beneficial FA, linoleic acid, was increased by 1.3-3.6%. The oil profile changes in these mutants were confirmed for four generations. Overall, our data illustrate that GmFATB1 knockout mutants have great potential in improving the soybean oil quality for human health.

The overexpression of OsACBP5 protects transgenic rice against necrotrophic, hemibiotrophic and biotrophic pathogens

The most devastating diseases in rice (Oryza sativa) are sheath blight caused by the fungal necrotroph Rhizoctonia solani, rice blast by hemibiotrophic fungus Magnaporthe oryzae, and leaf blight by bacterial biotroph Xanthomonas oryzae (Xoo). It has been reported that the Class III acyl-CoA-binding proteins (ACBPs) such as those from dicots (Arabidopsis and grapevine) play a role in defence against biotrophic pathogens. Of the six Arabidopsis (Arabidopsis thaliana) ACBPs, AtACBP3 conferred protection in transgenic Arabidopsis against Pseudomonas syringae, but not the necrotrophic fungus, Botrytis cinerea. Similar to Arabidopsis, rice possesses six ACBPs, designated OsACBPs. The aims of this study were to test whether OsACBP5, the homologue of AtACBP3, can confer resistance against representative necrotrophic, hemibiotrophic and biotrophic phytopathogens and to understand the mechanisms in protection. Herein, when OsACBP5 was overexpressed in rice, the OsACBP5-overexpressing (OsACBP5-OE) lines exhibited enhanced disease resistance against representative necrotrophic (R. solani & Cercospora oryzae), hemibiotrophic (M. oryzae & Fusarium graminearum) and biotrophic (Xoo) phytopathogens. Progeny from a cross between OsACBP5-OE9 and the jasmonate (JA)-signalling deficient mutant were more susceptible than the wild type to infection by the necrotroph R. solani. In contrast, progeny from a cross between OsACBP5-OE9 and the salicylic acid (SA)-signalling deficient mutant was more susceptible to infection by the hemibiotroph M. oryzae and biotroph Xoo. Hence, enhanced resistance of OsACBP5-OEs against representative necrotrophs appears to be JA-dependent whilst that to (hemi)biotrophs is SA-mediated.

Lipid Metabolism: Critical Roles in Male Fertility and Other Aspects of Reproductive Development in Plants

Fatty acids and their derivatives are essential building blocks for anther cuticle and pollen wall formation. Disruption of lipid metabolism during anther and pollen development often leads to genic male sterility (GMS). To date, many lipid metabolism-related GMS genes that are involved in the formation of anther cuticle, pollen wall, and subcellular organelle membranes in anther wall layers have been identified and characterized. In this review, we summarize recent progress on characterizing lipid metabolism-related genes and their roles in male fertility and other aspects of reproductive development in plants. On the basis of cloned GMS genes controlling biosynthesis and transport of anther cutin, wax, sporopollenin, and tryphine in Arabidopsis, rice, and maize as well as other plant species, updated lipid metabolic networks underlying anther cuticle development and pollen wall formation were proposed. Through bioinformatics analysis of anther RNA-sequencing datasets from three maize inbred lines (Oh43, W23, and B73), a total of 125 novel lipid metabolism-related genes putatively involved in male fertility in maize were deduced. More, we discuss the pathways regulating lipid metabolism-related GMS genes at the transcriptional and post-transcriptional levels. Finally, we highlight recent findings on lipid metabolism-related genes and their roles in other aspects of plant reproductive development. A comprehensive understanding of lipid metabolism, genes involved, and their roles in plant reproductive development will facilitate the application of lipid metabolism-related genes in gene editing, haploid and callus induction, molecular breeding and hybrid seed production in crops.

Overexpression of a Monocot Acyl-CoA-Binding Protein Confers Broad-Spectrum Pathogen Protection in a Dicot

Plants are continuously infected by various pathogens throughout their lifecycle. Previous studies have reported that the expression of Class III acyl-CoA-binding proteins (ACBPs) such as the Arabidopsis ACBP3 and rice ACBP5 were induced by pathogen infection. Transgenic Arabidopsis AtACBP3-overexpressors (AtACBP3-OEs) displayed enhanced protection against the bacterial biotroph, Pseudomonas syringae, although they became susceptible to the fungal necrotroph Botrytis cinerea. A Class III ACBP from a monocot, rice (Oryza sativa) OsACBP5 was overexpressed in the dicot Arabidopsis. The resultant transgenic Arabidopsis lines conferred resistance not only to the bacterial biotroph P. syringae but to fungal necrotrophs (Rhizoctonia solani, B. cinerea, Alternaria brassicicola) and a hemibiotroph (Colletotrichum siamense). Changes in protein expression in R. solani-infected Arabidopsis OsACBP5-overexpressors (OsACBP5-OEs) were demonstrated using proteomic analysis. Biotic stress-related proteins including cell wall-related proteins such as FASCILIN-LIKE ARABINOGALACTAN-PROTEIN10, LEUCINE-RICH REPEAT EXTENSIN-LIKE PROTEINS, XYLOGLUCAN ENDOTRANSGLUCOSYLASE/HYDROLASE PROTEIN4, and PECTINESTERASE INHIBITOR18; proteins associated with glucosinolate degradation including GDSL-LIKE LIPASE23, EPITHIOSPECIFIER MODIFIER1, MYROSINASE1, MYROSINASE2, and NITRILASE1; as well as a protein involved in jasmonate biosynthesis, ALLENE OXIDE CYCLASE2, were induced in OsACBP5-OEs upon R. solani infection. These results indicated that upregulation of these proteins in OsACBP5-OEs conferred protection against various plant pathogens.

Characterization of Oil Palm Acyl-CoA-Binding Proteins and Correlation of Their Gene Expression with Oil Synthesis

Acyl-CoA-binding proteins (ACBPs) are involved in binding and trafficking acyl-CoA esters in eukaryotic cells. ACBPs contain a well-conserved acyl-CoA-binding domain. Their various functions have been characterized in the model plant Arabidopsis and, to a lesser extent, in rice. In this study, genome-wide detection and expression analysis of ACBPs were performed on Elaeis guineensis (oil palm), the most important oil crop in the world. Seven E. guineensis ACBPs were identified and classified into four groups according to their deduced amino acid domain organization. Phylogenetic analysis showed conservation of this family with other higher plants. All seven EgACBPs were expressed in most tissues while their differential expression suggests various functions in specific tissues. For example, EgACBP3 had high expression in inflorescences and stalks while EgACBP1 showed strong expression in leaves. Because of the importance of E. guineensis as an oil crop, expression of EgACBPs was specifically examined during fruit development. EgACBP3 showed high expression throughout mesocarp development, while EgACBP1 had enhanced expression during rapid oil synthesis. In endosperm, both EgACBP1 and EgACBP3 exhibited increased expression during seed development. These results provide important information for further investigations on the biological functions of EgACBPs in various tissues and, in particular, their roles in oil synthesis.

Arabidopsis cytosolic acyl-CoA-binding proteins function in determining seed oil composition

As plant seed oils provide animals with essential fatty acids (FAs), genes that regulate plant lipid metabolism have been used in genetic manipulation to improve dietary seed oil composition and benefit human health. Herein, the Arabidopsis thaliana cytosolic acyl-CoA-binding proteins (AtACBPs), AtACBP4, AtACBP5, and AtACBP6 were shown to play a role in determining seed oil content by analysis of atacbp (atacbp4, atacbp5, atacbp6, atacbp4atacbp5, atacbp4atacbp6, atacbp5atacbp6, and atacbp4atacbp5atacbp6) seed oil content in comparison with the Col-0 wild type (WT). Triacylglycerol (TAG) composition in electrospray ionization-mass spectrometer (ESI-MS) analysis on atacbp6 seed oil showed a reduction (-50%) of C58-TAGs in comparison with the WT. Investigations on fatty acid composition of atacbp mutants indicated that 18:2-FA accumulated in atacbp6 and 18:3-FA in atacbp4, both at the expense of 20:1-FA. As TAG composition can be modified by acyl editing through phosphatidylcholines (PC) and lysophosphatidylcholines (LPC), total PC and LPC content in atacbp6 mature seeds was determined and ESI-MS analysis revealed that LPC had increased (+300%) at the expense of PC. Among all the 14 tested PC species, all (34:1-, 34:2-, 34:3-, 34:4-, 34:5-, 34:6-, 36:2-, 36:3-, 36:5-, 36:6-, 38:2-, 38:3-, and 38:4-PCs) but 36:4-PC were lower in atacbp6 than the WT. In contrast, all LPC species (16:0-, 18:1-, 18:2-, 18:3-, and 20:1-LPC) examined were elevated in atacbp6. LPC abundance also increased in atacbp4atacbp5, but not atacbp4 and atacbp5. Interestingly, when LPC composition in atacbp4atacbp5 was compared with atacbp4 and atacbp5, significant differences were observed between atacbp4atacbp5 and each single mutant, implying that AtACBP4 and AtACBP5 play combinatory roles by affecting LPC (but not PC) biosynthesis. Furthermore, PC-related genes such as those encoding acyl-CoA:lysophphosphatidylcholine acyltransferase (LPCAT1) and phospholipase A2 alpha (PLA2α) were upregulated in atacbp6 developing seeds. A model on the role of AtACBP6 in modulating TAG through regulating LPCAT1 and PLA2α expression is proposed. Taken together, cytosolic AtACBPs appear to affect unsaturated TAG content and are good candidates for engineering oil crops to enhance seed oil composition.

Arabidopsis acyl-CoA-binding proteins regulate the synthesis of lipid signals

Plant lipid signals are crucial developmental modulators and stress response mediators. A family of acyl-CoA-binding proteins (ACBPs) participates in the lipid trafficking of these signals. Isoform-specific functions can arise from differences in their subcellular distribution, tissue-specificity, stress-responsiveness, and ligand selectivity. In lipid-mediated cell signaling, plant ACBPs are not merely transporters but are also important regulators via their interaction with lipid-metabolic enzymes and precursor lipids. In this Insight, the regulatory roles of plant ACBPs in the synthesis of various signaling lipids, including phosphatidic acid, sterols, oxylipins, and sphingolipids, are reviewed. We focus on the functional significance of these lipid signals in plant development and stress responses with an overview of recent work using reverse genetics and transgenic Arabidopsis.

Functional and Structural Diversity of Acyl-coA Binding Proteins in Oil Crops

Diversities in structure and function of ACBP were discussed in this review. ACBP are important proteins that could transport newly synthesized fatty acid, activated into -coA, from plastid to endoplasmic reticulum, where oil in the form of triacylglycerol occurs. ACBP were detected in various animal and plants species, which indicated their importance in biological function. In fact, involvement of ACBP in important process such as lipid metabolism, regulation of enzyme and gene expression, and in response to plant stresses has been proven in several studies. In this review, findings on ACBP of 11 well-known oil crops were reviewed to comprehend diversity, comparative analyses on ACBP structure were made, and link between structure and function, tissue expression and subcellular location of ACBP were also observed. Incomplete reports in some species were mentioned, which might be encouraging to start or to perform deeper studies. Similar characteristics were found in paralogs ACBP, and orthologs ACBP had different functions, despite the high identity in amino acid sequence. At the end, it is confirmed that ortholog proteins could not necessarily display the same function, even from closely related species.

Arabidopsis ACYL-COA-BINDING PROTEIN1 interacts with STEROL C4-METHYL OXIDASE1-2 to modulate gene expression of homeodomain-leucine zipper IV transcription factors

Fatty acids (FAs) and sterols constitute building blocks of eukaryotic membranes and lipid signals. Co-regulation of FA and sterol synthesis is mediated by sterol regulatory element-binding proteins in animals but remains elusive in plants. We reported recently that Arabidopsis ACYL-COA-BINDING PROTEIN1 (ACBP1) modulates sterol synthesis via protein-protein interaction with STEROL C4-METHYL OXIDASE1-1 (SMO1-1). Herein, ACBP1 was demonstrated to co-express and interact with SMO1-2 by yeast two-hybrid, co-localization, pull-down, co-immunoprecipitation and β-glucuronidase assays. SMO1-2 silenced in acbp1 was used in phenotyping, GC-MS and expression profiling. ACBP1 co-expressed with SMO1-2 in embryo sacs, pollen and trichomes, corroborating with cooperative tissue-specific functions unseen with SMO1-1. SMO1-2 silencing in acbp1 impaired seed development, male and female gamete transmission, and pollen function. Genes encoding homeodomain-leucine zipper IV transcription factors (HDG5, HDG10, HDG11 and GLABRA2), which potentially bind phospholipids/sterols, were transcribed aberrantly. GLABRA2 targets (MYB23, MUM4 and PLDα1) were misregulated, causing glabra2-resembling trichome, seed coat mucilage and oil-accumulating phenotypes. Together with altered sterol and FA compositions upon ACBP1 mutation and/or SMO1-2 silencing, ACBP1-SMO1 interaction appears to mediate homeostatic co-regulation of FAs and sterols, which serve as lipid modulators for gene expression of homeodomain-leucine zipper IV transcription factors.

Depletion of Arabidopsis ACYL-COA-BINDING PROTEIN3 Affects Fatty Acid Composition in the Phloem

Oxylipins are crucial components in plant wound responses that are mobilised via the plant vasculature. Previous studies have shown that the overexpression of an Arabidopsis acyl-CoA-binding protein, AtACBP3, led to an accumulation of oxylipin-containing galactolipids, and AtACBP3pro::BETA-GLUCURONIDASE (GUS) was expressed in the phloem of transgenic Arabidopsis. To investigate the role of AtACBP3 in the phloem, reverse transcription-polymerase chain reaction and western blot analysis of phloem exudates from the acbp3 mutant and wild type revealed that the AtACBP3 protein, but not its mRNA, was detected in the phloem sap. Furthermore, micrografting demonstrated that AtACBP3 expressed from the 35S promoter was translocated from shoot to root. Subsequently, AtACBP3 was localised to the companion cells, sieve elements and the apoplastic space of phloem tissue by immunogold electron microscopy using anti-AtACBP3 antibodies. AtACBP3pro::GUS was induced locally in Arabidopsis leaves upon wounding, and the expression of wound-responsive jasmonic acid marker genes (JASMONATE ZIM-DOMAIN10, VEGETATIVE STORAGE PROTEIN2, and LIPOXYGENASE2) increased more significantly in both locally wounded and systemic leaves of the wild type in comparison to acbp3 and AtACBP3-RNAi. Oxylipin-related fatty acid (FA) (C18:2-FA, C18:3-FA and methyl jasmonate) content was observed to be lower in acbp3 and AtACBP3-RNAi than wild-type phloem exudates using gas chromatography-mass spectrometry. Experiments using recombinant AtACBP3 in isothermal titration calorimetry analysis showed that medium- and long-chain acyl-CoA esters bind (His)₆-AtACBP3 with K_D values in the micromolar range. Taken together, these results suggest that AtACBP3 is likely to be a phloem-mobile protein that affects the FA pool and jasmonate content in the phloem, possibly by its binding to acyl-CoA esters.

Arabidopsis thaliana Acyl-CoA-binding protein ACBP6 interacts with plasmodesmata-located protein PDLP8

In Arabidopsis thaliana, six acyl-CoA-binding proteins (ACBPs), designated as AtACBP1 to AtACBP6, have been identified to function in various events related to plant stress and development. The 10-kDa AtACBP6 is the smallest in this protein family, and recombinant AtACBP6 interacts with lipids in vitro by binding to acyl-CoA esters and phosphatidylcholine. Using anti-AtACBP6 antibodies in immunoelectron microscopy, we have localized AtACBP6 in the Arabidopsis phloem. The detection of immunogold grains in the plasmodesmata suggested that AtACBP6 could move from the companion cells to the sieve elements via the plasmodesmata. As AtACBP6 has been identified in a membrane-based interactome analysis to be a potential protein partner of Plasmodesmata-Localized Protein, PDLP8, AtACBP6-PDLP8 interaction was investigated herein utilizing isothermal titration calorimetry, as well as pull-down and bimolecular fluorescence complementation assays (BiFC). Notably, BiFC data revealed that AtACBP6-PDLP8 interaction occurred at the plasma membrane, which was unexpected as AtACBP6 has been previously identified in the cytosol. AtACBP6 expression was generally higher than PDLP8 in β-glucuronidase (GUS) assays on transgenic Arabidopsis transformed with AtACBP6 or PDLP8 promoter-driven GUS, consistent with qRT-PCR and microarray results. Furthermore, western blot analysis using anti-AtACBP6 antibodies showed a reduction in AtACBP6 expression in the pdlp8 T-DNA insertional mutant, suggesting that PDLP8 may possibly influence AtACBP6 accumulation in the sieve elements, probably in the plasmodesmata, where PDLP8 is confined and to where AtACBP6 has been immunodetected.

Regulation of fatty acid and flavonoid biosynthesis by miRNAs in Lonicera japonica

This is the first report on miRNAs in different varieties of honeysuckle, the results of which can further facilitate the discovery of functional regulatory miRNAs involved in the secondary metabolism in L. japonica.