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 comprehensive characterization and transcriptomic analysis of AP2/ERF gene family revealed its role in seed oil and ALA formation in perilla (Perilla frutescens)

Perilla (Perilla frutescens) is a potential specific oilseed crop with an extremely high α-linolenic acid (ALA) content in its seeds. AP2/ERF transcription factors (TFs) play important roles in multiple biological processes. However, limited information is known about the regulatory mechanism of the AP2/ERF family in perilla's oil accumulation. In this research, we identified 212 AP2/ERF family members in the genome of perilla, and their domain characteristics, collinearity, and sub-genome differentiation were comprehensively analyzed. Transcriptome sequencing revealed that genes encoding key enzymes involved in oil biosynthesis (e.g., ACCs, KASII, GPAT, PDAT and LPAAT) were up-regulated in the high-oil variety. Moreover, the endoplasmic reticulum-localized FAD2 and FAD3 were significantly up-regulated in the high-ALA variety. To investigate the roles of AP2/ERFs in lipid biosynthesis, we conducted a correlation analysis between non-redundant AP2/ERFs and key lipid metabolism genes using WGCNA. A significant correlation was found between 36 AP2/ERFs and 90 lipid metabolism genes. Among them, 12 AP2/ERFs were identified as hub genes and showed significant correlation with lipid synthase genes (e.g., FADs, GPAT and ACSL) and key regulatory TFs (e.g., LEC2, IAA, MYB, UPL3). Furthermore, gene expression analysis identified three AP2/ERFs (WRI, ABI4, and RAVI) potentially playing an important role in the regulation of oil accumulation in perilla. Our study suggests that PfAP2/ERFs are important regulatory TFs in the lipid biosynthesis pathway, providing a foundation for the molecular understanding of oil accumulation in perilla and other oilseed crops.

Assessing the Potential for Genome-Assisted Breeding in Red Perilla Using Quantitative Trait Locus Analysis and Genomic Prediction

Red perilla is an important medicinal plant used in Kampo medicine. The development of elite varieties of this species is urgently required. Medicinal compounds are generally considered target traits in medicinal plant breeding; however, selection based on compound phenotypes (i.e., conventional selection) is expensive and time consuming. Here, we propose genomic selection (GS) and marker-assisted selection (MAS), which use marker information for selection, as suitable selection methods for medicinal plants, and we evaluate the effectiveness of these methods in perilla breeding. Three breeding populations generated from crosses between one red and three green perilla genotypes were used to elucidate the genetic mechanisms underlying the production of major medicinal compounds using quantitative trait locus analysis and evaluating the accuracy of genomic prediction (GP). We found that GP had a sufficiently high accuracy for all traits, confirming that GS is an effective method for perilla breeding. Moreover, the three populations showed varying degrees of segregation, suggesting that using these populations in breeding may simultaneously enhance multiple target traits. This study contributes to research on the genetic mechanisms of the major medicinal compounds of red perilla, as well as the breeding efficiency of this medicinal plant.

Development of an SNP marker set for marker-assisted backcrossing using genotyping-by-sequencing in tetraploid perilla

High-quality molecular markers are essential for marker-assisted selection to accelerate breeding progress. Compared with diploid species, recently diverged polyploid crop species tend to have highly similar homeologous subgenomes, which is expected to limit the development of broadly applicable locus-specific single-nucleotide polymorphism (SNP) assays. Furthermore, it is particularly challenging to make genome-wide marker sets for species that lack a reference genome. Here, we report the development of a genome-wide set of kompetitive allele specific PCR (KASP) markers for marker-assisted recurrent selection (MARS) in the tetraploid minor crop perilla. To find locus-specific SNP markers across the perilla genome, we used genotyping-by-sequencing (GBS) to construct linkage maps of two F2 populations. The two resulting high-resolution linkage maps comprised 2326 and 2454 SNP markers that spanned a total genetic distance of 2133 cM across 16 linkage groups and 2169 cM across 21 linkage groups, respectively. We then obtained a final genetic map consisting of 22 linkage groups with 1123 common markers from the two genetic maps. We selected 96 genome-wide markers for MARS and confirmed the accuracy of markers in the two F2 populations using a high-throughput Fluidigm system. We confirmed that 91.8% of the SNP genotyping results from the Fluidigm assay were the same as the results obtained through GBS. These results provide a foundation for marker-assisted backcrossing and the development of new varieties of perilla.

The Genome of the Korean Island-Originated Perilla citriodora 'Jeju17' Sheds Light on Its Environmental Adaptation and Fatty Acid and Lipid Production Pathways

Perilla is a key component of Korean food. It contains several plant-specialized metabolites that provide medical benefits. In response to an increased interest in healthy supplement food from the public, people are focusing on the properties of Perilla. Nevertheless, unlike rice and soybeans, there are few studies based on molecular genetics on Perilla, so it is difficult to systematically study the molecular breed. The wild Perilla, Perilla citriodora 'Jeju17', was identified a decade ago on the Korean island of Jeju. Using short-reads, long-reads, and Hi-C, a chromosome-scale genome spanning 676 Mbp, with high contiguity, was assembled. Aligning the 'Jeju17' genome to the 'PC002' Chinese species revealed significant collinearity with respect to the total length. A total of 31,769 coding sequences were predicted, among which 3331 were 'Jeju17'-specific. Gene enrichment of the species-specific gene repertoire highlighted environment adaptation, fatty acid metabolism, and plant-specialized metabolite biosynthesis. Using a homology-based approach, genes involved in fatty acid and lipid triacylglycerol biosynthesis were identified. A total of 22 fatty acid desaturases were found and comprehensively characterized. Expression of the FAD genes in 'Jeju17' was examined at the seed level, and hormone signaling factors were identified. The results showed that the expression of FAD genes in 'Jeju17' at the seed level was high 25 days after flowering, and their responses of hormones and stress were mainly associated with hormone signal transduction and abiotic stress via cis-elements patterns. This study presents a chromosome-level genome assembly of P. citriodora 'Jeju17', the first wild Perilla to be sequenced from the Korean island of Jeju. The analyses provided can be useful in designing ALA-enhanced Perilla genotypes in the future.

RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed

Perilla frutescen is used as traditional food and medicine in East Asia. Its seeds contain high levels of α-linolenic acid (ALA), which is important for health, but is scarce in our daily meals. Previous reports on RNA-seq of perilla seed had identified fatty acid (FA) and triacylglycerol (TAG) synthesis genes, but the underlying mechanism of ALA biosynthesis and its regulation still need to be further explored. So we conducted Illumina RNA-sequencing in seven temporal developmental stages of perilla seeds. Sequencing generated a total of 127 million clean reads, containing 15.88 Gb of valid data. The de novo assembly of sequence reads yielded 64,156 unigenes with an average length of 777 bp. A total of 39,760 unigenes were annotated and 11,693 unigenes were found to be differentially expressed in all samples. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 486 unigenes were annotated in the “lipid metabolism” pathway. Of these, 150 unigenes were found to be involved in fatty acid (FA) biosynthesis and triacylglycerol (TAG) assembly in perilla seeds. A coexpression analysis showed that a total of 104 genes were highly coexpressed (r > 0.95). The coexpression network could be divided into two main subnetworks showing over expression in the medium or earlier and late phases, respectively. In order to identify the putative regulatory genes, a transcription factor (TF) analysis was performed. This led to the identification of 45 gene families, mainly including the AP2-EREBP, bHLH, MYB, and NAC families, etc. After coexpression analysis of TFs with highly expression of FAD2 and FAD3 genes, 162 TFs were found to be significantly associated with two FAD genes (r > 0.95). Those TFs were predicted to be the key regulatory factors in ALA biosynthesis in perilla seed. The qRT-PCR analysis also verified the relevance of expression pattern between two FAD genes and partial candidate TFs. Although it has been reported that some TFs are involved in seed development, more direct evidence is still needed to verify their function. However, these findings can provide clues to reveal the possible molecular mechanisms of ALA biosynthesis and its regulation in perilla seed.

Two types of alcohol dehydrogenase from Perilla can form citral and perillaldehyde

Studies on the biosynthesis of oil compounds in Perilla will help in understanding regulatory systems of secondary metabolites and in elucidating reaction mechanisms for natural product synthesis. In this study, two types of alcohol dehydrogenases, an aldo-keto reductase (AKR) and a geraniol dehydrogenase (GeDH), which are thought to participate in the biosynthesis of perilla essential oil components, such as citral and perillaldehyde, were isolated from three pure lines of perilla. These enzymes shared high amino acid sequence identity within the genus Perilla, and were expressed regardless of oil type. The overall reaction from geranyl diphosphate to citral was performed in vitro using geraniol synthase and GeDH to form a large proportion of citral and relatively little geraniol as reaction products. The biosynthetic pathway from geranyl diphosphate to citral, the main compound of citral-type perilla essential oil, was established in this study.

Effects of SYJN, a Chinese herbal formula, on chronic unpredictable stress-induced changes in behavior and brain BDNF in rats

AIM OF THE STUDY

Suyu-Jiaonang (SYJN) is a Chinese herbal formula that contains four herbs: Bupleurum chinense DC, Curcuma aromatica Salisb., Perilla frutescens (Linn.) Britt., and Acorus tatarinowii Schott. Previous studies conducted in our laboratory have revealed an antidepressant-like effect of the formula in various mouse models of behavioral despair. The present study aimed to investigate whether SYJN could produce antidepressant-like effects in chronic unpredictable stress (CUS)-induced depression model in rats and its possible mechanism(s).

MATERIALS AND METHODS

Rats were subjected to an experimental setting of CUS. The effect of SYJN treatment on CUS-induced depression was examined using behavioral tests including the sucrose consumption and open field tests. The mechanism underlying the antidepressant-like action of SYJN was examined by measuring brain-derived neurotrophic factor (BDNF) protein and mRNA expression in brain tissues of CUS-exposed rats.

RESULTS

Exposure to CUS for 4 weeks caused depression-like behavior in rats, as indicated by significant decreases in sucrose consumption and locomotor activity (assessed in the open field test). In addition, it was found that BDNF protein and mRNA levels in the hippocampus and frontal cortex were lower in CUS-treated rats, as compared to controls. Daily intragastric administration of SYJN (1300 or 2600 mg/kg) during the 4-week period of CUS significantly suppressed behavioral changes and attenuated the CUS-induced decrease in BDNF protein and mRNA levels in the hippocampus and frontal cortex.

CONCLUSION

The results suggest that SYJN alleviates depression induced by CUS. The antidepressant-like activity of SYJN is likely mediated by the increase in BDNF expression in brain tissues.

Authentic seed-specific activity of the Perilla oleosin 19 gene promoter in transgenic Arabidopsis

The Perilla (Perilla frutescens L. cv. Okdong) oleosin gene, PfOle19, produces a 19-kDa protein that is highly expressed only in seeds. The activity of the -2,015 bp 5'-upstream promoter region of this gene was investigated in transgenic Arabidopsis plants using the fusion reporter constructs of enhanced green fluorescent protein (EGFP) and beta-glucuronidase (GUS). The PfOle19 promoter directs Egfp expression in developing siliques, but not in leaves, stems or roots. In the transgenic Arabidopsis, EGFP fluorescence and histochemical GUS staining were restricted to early seedlings, indehiscent siliques and mature seeds. Progressive 5'-deletions up to the -963 bp position of the PfOle19 promoter increases the spatial control of the gene expression in seeds, but reduces its quantitative levels of expression. Moreover, the activity of the PfOle19 promoter in mature seeds is 4- and 5-fold greater than that of the cauliflower mosaic virus 35S promoter in terms of both EGFP intensity and fluorometric GUS activity, respectively.

[SNP marker and allele-specific diagnostic PCR for authenticating herbs of Perilla]

AIM

To authenticate all the varieties of Perilla (single-species genus), to analyze sequences of rDNA ITS regions and single nucleotide polymorphism (SNP) within them and based on these, to design allele-specific diagnostic PCR primers.

METHODS

The rDNA ITS regions of the perilla varieties were sequenced and analyzed by Clustal X 1.8, MEGA 3.0. Allele-specific diagnostic PCR primers that can authenticate all the perilla varieties were designed based on SNPs loci.

RESULTS

The length of rDNA ITS sequences of perilla varieties ranged from 612 to 615 bp in size, including ITS1 (230 -232 bp), 5.8S (179 bp) and ITS2 (203 -204 bp). The GC content is about 61.5% - 61.9%. There is not only SNPs in non-coding region ITS1 and ITS2 (ncSNP), but also three coding SNPs (cSNP) loci in the conservative region of 5.8S. All the SNPs have only two allele loci polymorphism. The cSNP in 5.8S is related to the morphology variation among the varieties. Allele-specific diagnostic PCR primers have been designed according to SNPs loci to authenticate accurately all the seeds and leaves of Perilla varieties.

CONCLUSION

SNPs in rDNA ITS region can be used as an effective molecular markers to authenticate all the varieties of Perilla.

Limonene production in tobacco with Perilla limonene synthase cDNA

Limonene synthase (LS) catalyses the stereo-specific cyclization of geranyl diphosphate (GPP) to form a monocyclic monoterpene, limonene. In an attempt to engineer monoterpene biosynthesis, three expression constructs of LS cDNA of Perilla frutescens, which were designed to be localized in either the plastid, the cytosol or the endoplasmic reticulum (ER), were introduced into tobacco in order to examine differences in enzyme activity and the productivity of limonene. High and moderate enzyme activity, respectively, was observed for plastid- and cytosol-localized LS, whereas no enzyme activity was seen for ER-localized LS, suggesting that the plastid is the preferred compartment for LS, while LS may also have an active form in the cytosol. The formation of limonene in vivo was confirmed by gas chromatography-mass spectrometry (GC-MS) in leaf extracts of both plastid- and cytosol-localized LS transgenic plants. The amount of limonene in plastid-localized LS transgenic plants was 143 ng g-1 fresh wt, whereas that in the cytosol-type was 40 ng g-1 fresh wt, and these limonene contents increased by 2.7-fold and 3.0-fold, respectively, with the addition of methyl jasmonate. The headspace analyses showed that the plastid- and the cytosol-localized LS transgenic plants (12 cm high) emitted 390 ng and 515 ng limonene per month, respectively. The possibility of genetically engineering monoterpene production is discussed.

A WD-repeat-containing putative regulatory protein in anthocyanin biosynthesis in Perilla frutescens

The WD-repeat proteins are found in eukaryotes and play an important role in the regulation of a wide variety of cellular functions such as signal transduction, transcription, and proliferation. In this study, we have isolated a cDNA encoding a novel WD-repeat protein, PFWD, from the anthocyanin-pigmented leaves of Perilla frutescens using AN11 cDNA from Petunia hybrida as the probe. The C-terminal region of PFWD contains a WD repeat that is highly conserved in homologous proteins from a variety of organisms that do not produce anthocyanin such as yeast, nematodes and mammals. Transgenic Arabidopsis plants overexpressing PFWD exhibited phenotypic changes including enhancement of anthocyanin production and reduced viability. A study of the interaction between PFWD and anthocyanin regulatory proteins using a yeast two-hybrid system showed strong interaction between PFWD and MYC-RP, a MYC-like protein from P. frutescens. PFWD fusion proteins transiently expressed in onion epidermal cells were localized in the cytosol under both dark and light conditions. However, co-expression of PFWD and MYC-RP fusion proteins resulted in nuclear localization of PFWD. We propose a model of genetic regulation in which the PFWD protein acts in signal transduction process in a variety of pathways through protein interaction with MYC proteins.