Next Article in Journal
Effect of Smut Infection on the Photosynthetic Physiological Characteristics and Related Defense Enzymes of Sugarcane
Next Article in Special Issue
Utilization of Genotyping-by-Sequencing (GBS) for Rice Pre-Breeding and Improvement: A Review
Previous Article in Journal
Atlas of Nervous System Vascular Malformations: A Systematic Review
Previous Article in Special Issue
Insight into the Roles of Proline-Rich Extensin-like Receptor Protein Kinases of Bread Wheat (Triticum aestivum L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 12
Issue 8
10.3390/life12081200
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Peroxidase and Laccase Gene Families and In Silico Identification of Potential Genes Involved in Upstream Steps of Lignan Formation in Sesame

by
Yedomon Ange Bovys Zoclanclounon
1,2,
Michael Rostás
3,
Nam-Jin Chung
2,
Youngjun Mo
2,
Petr Karlovsky
3 and
Komivi Dossa
3,*
1
Genomics Division, National Institute of Agricultural Sciences, Rural Development Administration, Jeonju 54874, Korea
2
Department of Crop Science and Biotechnology, Jeonbuk National University, Jeonju 54896, Korea
3
Molecular Phytopathology and Mycotoxin Research, Grisebachstrasse 6, Georg-August-University Goettingen, D-37077 Goettingen, Germany
*
Author to whom correspondence should be addressed.
Life 2022, 12(8), 1200; https://0-doi-org.brum.beds.ac.uk/10.3390/life12081200
Submission received: 26 July 2022 / Revised: 5 August 2022 / Accepted: 6 August 2022 / Published: 8 August 2022
(This article belongs to the Special Issue Research Advances in Plant Genomics: 2nd Edition)
Browse Figures
Versions Notes

Abstract

:
Peroxidases and laccases are oxidative enzymes involved in physiological processes in plants, covering responses to biotic and abiotic stress as well as biosynthesis of health-promoting specialized metabolites. Although they are thought to be involved in the biosynthesis of (+)-pinoresinol, a comprehensive investigation of this class of enzymes has not yet been conducted in the emerging oil crop sesame and no information is available regarding the potential (+)-pinoresinol synthase genes in this crop. In the present study, we conducted a pan-genome-wide identification of peroxidase and laccase genes coupled with transcriptome profiling of diverse sesame varieties. A total of 83 and 48 genes have been identified as coding for sesame peroxidase and laccase genes, respectively. Based on their protein domain and Arabidopsis thaliana genes used as baits, the genes were classified into nine and seven groups of peroxidase and laccase genes, respectively. The expression of the genes was evaluated using dynamic transcriptome sequencing data from six sesame varieties, including one elite cultivar, white vs black seed varieties, and high vs low oil content varieties. Two peroxidase genes (SiPOD52 and SiPOD63) and two laccase genes (SiLAC1 and SiLAC39), well conserved within the sesame pan-genome and exhibiting consistent expression patterns within sesame varieties matching the kinetic of (+)-pinoresinol accumulation in seeds, were identified as potential (+)-pinoresinol synthase genes. Cis-acting elements of the candidate genes revealed their potential involvement in development, hormonal signaling, and response to light and other abiotic triggers. Transcription factor enrichment analysis of promoter regions showed the predominance of MYB binding sequences. The findings from this study pave the way for lignans-oriented engineering of sesame with wide potential applications in food, health and medicinal domains.

1. Introduction

Sesame (Sesamum indicum L.), a member of Pedaliaceae family, is an oil crop whose seeds contain lignans including (+)- sesamin, (+)- sesamolin, and (+)- sesaminol [1,2]. Therapeutical properties of sesame lignans against neurodegenerative diseases [3,4], prostate, and breast cancers [5] have been reported. Besides, lignans represent an emerging perspective for health care and disease prevention as functional foods and nutraceuticals [6,7,8,9,10,11,12]. The lignan market is exploding and may reach over USD 610 Million by 2028 [13] with multiple applications covering food, pharmaceutics, and cosmetics industries. Meanwhile, several patents have been deposited concerning the extraction, purification, and transformation of lignans from sesame [14,15,16,17,18,19,20,21,22], showing the growing interest in this class of plant-specialized metabolites.
To date, tremendous works have been done to elucidate the lignan biosynthesis pathway in sesame by using wild (S. alatum and S. radiatum) and cultivated (S. indicum) material [1,23,24,25,26]. A total of six enzymes involved in lignan biosynthesis in sesame have been characterized, including two cytochrome P450 coding genes (CYP81Q1 and CYP92B14), and four glycosyltransferases (UGT71A9, UGT94D1, UGT94AG1, and UGT94AA2) [1,25]. The connections between the identified enzymes and their respective targets are depicted in Figure 1. At the initial step of the lignan biosynthesis pathway, an oxidative coupling reaction involving two molecules of conyferol alcohol takes place. With the help of a dirigent protein (DIR), the primary precursor of lignan, (+)-pinoresinol, is generated. The latter is then sequentially converted by CYP81Q1 to produce (+)-piperitol and (+)-sesamin. Further, (+)-sesamolin and (+)-sesaminol synthesis is guided by CYP92B14.
It is worth mentioning that the mechanism for the synthesis of the central precursor ie (+)-pinoresinol, very likely involves a combinatory action of oxidase enzymes (laccases and/or peroxidases) with a dirigent protein [27]. Oxidation by peroxidases and/or laccases followed by a stereo-selective radical coupling guided by a dirigent protein is assumed to take place at an early stage of the synthesis of all plant lignans [27,28,29,30]. In brief, one-electron oxidation on conyferyl alcohol is catalyzed by peroxides and/or laccasses. As a result, an intermediate molecule, a conyferyl alcohol (1)-derived free radical, is formed. Then an intramolecular cyclization guided by a dirigent protein induces the formation of (+)-pinoresinol molecule [27].
To the best of our knowledge, evidence of oxidase activity involved in the critical precursor step of sesame lignan biosynthesis has not been yet established. While a dirigent protein (XP_011080883) has been detected [2] in sesame by sequence homology, the proof of its function as well as the identity of the peroxidase or laccase responsible for the conversion of coniferyl alcohol to (+)-pinoresinol remain to be established.
Peroxidases and laccases are multifunctional enzymes playing a wide range of roles in plants covering biotic [31,32,33] and abiotic responses [34,35,36,37,38], and other biological processes such as fiber initiation [39], cell elongation [40], lignification [41,42], seed setting and panicle branching [43], pigmentation [44,45], and flavonoid oxidation [46]. Peroxidases and laccases are also used in therapeutic and industrial applications. For instance, peroxidases are used in industrial chemical synthesis, diagnostic tests, and enzyme immunoassays [47]. Meanwhile, laccase has great importance in the paper industry due to its capacity for delignification [48,49]. Moreover, laccases are also useful in ethanol production, wine clarification, industrial effluents treatment, herbicide degradation, dyes decoloration, and drug analysis [50].
Owing to the importance of the peroxidases and laccases in plants, the present study was carried out to comprehensively characterize these enzymes in sesame and outline candidate genes likely involved in the early steps of lignan biosynthesis.

2. Materials and Methods

2.1. Genome-Wide Identification of Peroxidase and Laccase Genes and Core Genes Inference

Putative peroxidase and laccase genes were searched for in the genome data of S. indicum var. Zhongzhi13 (NCBI RefSeq accession: GCF_000512975.1), S. indicum var. Goenbaek (https://zenodo.org/record/6350881, accessed on 15 March 2022), S. indicum cv Mishouzhima, S. indicum cv. Baizhima, S. indicum var. Yuzhi11, and S. indicum var. Swetha (http://www.sesame-bioinfo.org/pan-genome, accessed on 15 November 2018).
NCBI HMM accession TIGR03390.1 (EC 1.10.3.2) served for the identification of candidate laccase genes, whereas PFAM HMM accession PF00141.26 was employed to detect peroxidase genes. After a hit search using hmmsearch (-E 1 × 10−5—domE 1 × 10−5), a domain verification was executed with PfamScan v1.6 [51] to check the presence of the POD (peroxidase) and LAC (laccase) domains. Spurious genes were filtered out. An additional check of the presence of the POD and LAC domains was performed using the InterProScan v5 [52].
In order to infer the core conserved POD and LAC genes within the sesame pangenome dataset, Orthofinder v2.3.12 [53] was run with the default settings.
The identified POD and LAC genes from the reference genome of S. indicum var Zhongzhi13 were retained for downstream analysis.

2.2. Chromosome Location and Syntheny Analyses

Genome mapping of SiPOD (peroxidase genes) and SiLAC (laccase genes) was rendered using MG2C V2.1 [54] based on the annotation information. Sesame-to-sesame and sesame-to-Arabidopsis synteny blocks were investigated with the MCScanX toolkit [55]. Regarding the evolutionary origin of the duplicated genes, the duplicate genes classifier Perl script from MCScanX helped to distinguish between genes evolved by tandem or segmental duplication.

2.3. Phylogenetic Analysis

Peroxidase and laccase genes from A. thaliana were added to SiPOD and SiLAC genes for the construction of phylogenetic trees. Prior to the tree inference, the genes were aligned using MAFFT v7.464-0 [56]. The resulting alignment was trimmed with trimAl v1.4.1 [57]. Subsequently, the trees were constructed with IQ-TREE v1.6.12 [58]. Peroxidase and laccase maximum likelihood trees were inferred following models LG+R5 and LG+I+G4 with 1000 iterations, respectively. The tree models were selected using the ModelFinder package [59].

2.4. RNA-Seq Data Retrieval

A set of six sesame varieties (Zhongzhi13, Zhongfengzhi No.1, Zhongzhi No.33, ZZM4728, ZZM2161, and ZZM3495) were selected for the investigation of peroxidase and laccase gene expression using RNA sequencing data (Table S1). The criteria for the selection were (a) oil content (high oil-producing variety versus low oil-producing variety) and (b) seed color (white versus black). Thus, in addition to the multi-organs transcriptome data of the reference genome Zhongzhi13 [60], the seed RNA-Seq of two pure lines of sesame Zhongfengzhi No.1 and Zhongzhi No.33 [61], exhibiting white and black seed color respectively, were used. The seed transcriptome sequences of one high (59.1%) oil content (ZZM4728) and two low oil content ZZM2161 (48.4%) and ZZM3495 (51.0%) varieties [62], were also downloaded from NCBI. The sesame varieties were grown under identical growth conditions at the Oil Crops Research Institute (OCRI) experimental station in Hubei Province, Wuhan, China, at N 30.57°, E 114.30°, altitude 27 m [60,61,62]. The RNA samplings were performed at 10, 20, 25, and 30 days after anthesis for ZZM2161, ZZM3495, and ZZM4728, while Zhongfengzhi No.1 and Zhongzhi No.33 were sampled at 5, 8, 11, 14, 17, 20, 23, 26, and 30 days after anthesis. Detailed information regarding the materials and RNA-Seq raw data SRAs are provided in Table S1.

2.5. Expression Profile Analysis and Candidate Genes Selection

The RNA-Seq raw data were quality-checked using FastQC v0.11.2 [63]. Sequencing adapters and low-quality (Q < 30) reads were filtered out with Trimmomatic v0.36 [64]. The clean data were then mapped to the reference genome using HISAT v2.2.1 [65]. The gene expression profile in each tissue was assessed with the RSEM package v1.3.3 [66] as fragments per kilobase of transcript per million fragments mapped (FPKM). Heatmaps showing the expression in different tissues were plotted with TBTools v1.098746 [67] with log2(FPKM) values.
To select candidate genes, we applied three major filtering criteria. Firstly, from the reference genome dataset, the genes of interest should be preferentially expressed in the seed tissue compared to the root, leaf, stem, and capsule. Secondly, the candidate genes from the first filtering step have been checked for their expression in high versus low oil content varieties with a particular emphasis on those that were expressed early during seed development. The early-stage criterion was included since pinoresinol is the precursor of all lignans in sesame (Figure 1). Additionally, the kinetic of the biosynthesis of the pinoresinol described by Ono et al. [23] implied an increasing expression at the early stage of the seed development that enables later accumulation of downstream lignans in the seed at the maturity stage. Thirdly, after the second filtering step, the same approach was applied to the white versus black seed dataset. It is worth mentioning that the seed color was considered in this study because it has been reported that white sesame seeds contain higher levels of lignans (sesamolin and sesamin) compared to black sesame seeds [68,69].
Overall, the retained candidate genes belonged to the core gene repertoire from the sesame pan-genome; were preferentially expressed in seed tissues; and were strongly expressed at the early stage of the seed development.

2.6. Conserved Motifs, Gene Structure, GO Annotation, and Orthologs Detection of Candidate SiPOD and SiLAC Genes

The candidate genes from the three filtering steps were screened for conserved motifs using MEME Suite v5.0.4 [70] with a maximum number of motifs set to 20. A 2 Kbp promoter sequence of candidate genes were submitted to the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 6 February 2022) [71] to find out cis-acting regulatory elements. Besides, the InterProScan v5 [52]) was employed to find out molecular functions associated with the candidate genes. Meanwhile, their orthologs search was executed with SHOOT (https://www.shoot.bio/, accessed on 6 February 2022) [72] using the plant database option.

2.7. Transcription Factor Enrichment Analysis

Taking advantage of the Plant Transcriptional Regulatory Map (PlantRegMap) platform [73], we performed a transcription factor (TF) enrichment analysis (http://plantregmap.gao-lab.org/tf_enrichment.php, accessed on 8 February 2022) in order to estimate the most contributive TF families potentially involved in the regulation of peroxidase and laccase genes.

2.8. Gene Expression Assays Using qRT-PCR

The expression profile of the candidate genes was verified by performing a qRT-PCR in a LightCycler® 480II real-time PCR instrument (Roche Diagnostics, Rotkreuz, Switzerland). Prior to the PCR, five pairs of primers (Table S2) were designed using the primer 3 blast online tool (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/tools/primer-blast/, accessed on 11 April 2022). The qRT-PCR experiment was carried out according to Dossa et al. [74]. A total of three independent replicates were applied for each gene. The sesame actin gene (NCBI gene ID: LOC105159390) served as a positive control. The 2−ΔΔCT method [75] was employed to calculate the relative expression of the target genes.

3. Results and Discussion

3.1. Variability of Peroxidase and Laccase Genes in Sesame Pangenome and Phylogenetic Analyses

From the sesame pan-genome gene sets, 83, 82, 60, 59, 54, and 54 peroxidase genes were counted in Zongzhi13, Swetha, Mishuozima, Baizhima, Goenbaek, and Yuzhi11 genomes, respectively (Figure 2). Similarly to the peroxidase count variability observed at the intra-species level in sesame, high inter-species variability was observed with 138, 119, 102, 90, 73, and 47 peroxidase genes counted in Oryza sativa [76], Zea mays [77], Solanum tuberosum [78], Betula pendula [79], A. thaliana [80], and Vitis vinifera [81], respectively. Most of the gene clusters (25) were shared by all varieties while only Swetha exhibited species-specific gene clusters (7).
An average of 44 ± 7 laccase genes were identified in the sesame pangenome. The most abundant laccases were observed in Swetha (n = 56), followed by Zhongzhi13 (n = 48), Goenbaek (n = 45), Baizhima (n = 42), Yuzhi11 (n = 40), and Mishuozhima (n = 35). Similar laccases counts were observed in land plants Prunus persica (n = 48) [82], Panicum virgatum (n = 49) [83], and Solanum melongena (n = 42) [84]; a lower number of laccase genes was reported in A. thaliana (n = 17) [85].
The orthology analysis revealed that the core laccase genes were grouped into 20 clusters, while Swetha harbored 6 specific gene clusters. Except for Yuzhi11 and Goenbaek, the peroxidase and laccase genes were globally abundant in modern varieties (Zhongzhi13, Swetha) compared to landraces (Mishuozima and Baizhima); suggesting the influence of the human oil-oriented selection. In fact, at a whole-genome scale, landraces (Mishuozima and Baizhima) exhibited specific genes coding for environmental adaption, while modern varieties showed preferential genes with oil-related functional attributes [86].
Interestingly, only Swetha showed a specific gene cluster suggesting a unique peroxidase and laccase gene repertoire in this variety. However, this should be interpreted with caution since the Swetha chromosome-scale genome was constructed based on the reference genome of Zhongzhi13 [86] with short-reads assembly as an initial contigs-level assembly.
For downstream analyses, peroxidases and laccases genes from the reference genome Zhongzhi13 were used (Tables S3 and S4).
After locating the identified genes on the chromosomes (Figure 3), we assessed the evolutionary determinants of the distribution of SiPOD and SiLAC genes within the sesame genome. Two types of gene relationships were found (Figures S1 and S2): paralogous genes that are adjacent on the same chromosome (tandem-duplicated genes) and genes that are far away from each other, often located on different chromosomes (segmental duplicated genes). For SiPOD, the results highlighted 18 pairs paralogous genes that underwent segmental duplication, while one pair originated from tandem duplication (Figure S1). Tandem and segmental duplications are considered as an evolutionary driving force resulting in gene families’ expansion [87,88]. Peroxidase genes duplication via tandem and segmental duplications were extensively reported in plants, including soybean [36], cassava [89], carrot [90], potato [78], cotton [91], watermelon [92], maize [77], tomato [93], Chinese pear [94] and others. In our study, peroxidase gene expansion mainly resulted from segmental duplication, which is consistent with the findings of Cao et al. [94] in pear. However, in trees such as Betula pendula [79] and Populus trichocarpa [95], tandem duplication was the main driver of peroxidase gene expansion.
Similarly, a total of 20 SiLAC genes belonged exclusively to the segmental duplication group (Figure S2). In contrast, only a tandem duplication pattern was detected in Solanum melongena with 16 laccase genes [84]. In Panicum virgatum, both tandem and segmental duplications of laccase genes were found [83].
The syntheny analysis revealed that 31 (37%) SiPOD and 26 (54%) SiLAC genes showed syntheny with A. thalina respective gene sets (Figures S3 and S4), suggesting that they are conserved within these species.
To classify the identified peroxidase and laccase genes, phylogenetic trees were constructed using A. thaliana genes as baits. The results revealed nine and seven groups for peroxidase and laccase gene families, respectively (Figure 4 and Figure 5). The tree topology was in line with sesame-alone peroxidase and laccase genes trees (Figures S5 and S6).

3.2. Expression Profiles of SiPOD and SiLAC Genes in Different Tissues and Candidate Genes Selection

To identify potential (+)-pinoresinol synthase genes, a multi-varieties comparative transcriptome approach was utilized. RNA-Seq data from six sesame varieties (Table S1) were inspected following the three filtering steps including (i) preferential expression in seed, (ii) expression at least at the early stage of the seed development, and (ii) expression of the gene across all tested varieties.
For peroxidase genes (Figure 6), a set of eight SiPOD (SiPOD41, SiPOD42, SiPOD47, SiPOD50, SiPOD52, SiPOD8, SiPOD63, and SiPOD65) genes were preferentially expressed in the seeds of Zhongzhi 13 (Figure 6A).
The selected genes were mined for their expression in the seeds of low (LOS) versus high (HOS) oil content varieties (Figure 6B). The results pinpointed SiPOD42, SiPOD52, SiPOD58, SiPOD63, SiPOD65, and SiPOD50. SiPOD42 was most expressed in HOS at 10 days as compared to the two LOS varieties. Its expression was maintained within all the seed development stages preferentially in the HOS. Similar observations were noted for SiPOD52, SiPOD58, and SiPOD50 with an ascending expression from 10 to 20 days after anthesis followed by a decline starting from 25 days after anthesis. Interestingly, SiPOD63 and SiPOD65 were expressed at the early stage (10 days after anthesis) in all varieties before the expression fell in the following development stages.
The candidate genes from the second filtering round were screened for preferential high expression in seed or at least expression at an earlier stage of seed development regardless of seed coat color (Figure 6C). Thus, the genes SiPOD52, SiPOD63, SiPOD50, and SiPOD65 came up to be the potential candidates. SiPOD52 and SiPOD63 showed a higher expression in white seed at early stage (five and eight days after anthesis). Intriguingly, a higher expression of SiPOD50 was noted in black seed at five and eight days after anthesis. However, the expression of the gene was quite stable within all development stages, regardless of the seed color.
As for the laccase genes (Figure 7), SiLAC1, SiLAC12, and SiLAC39 came out on the top in the first filtering round (Figure 7A). From this set, SiLAC39 showed the peak of expression in both low and high oil content varieties at 20 days after anthesis. Similarly, SiLAC1 expression was higher at 25 days after anthesis regardless of the type of variety (Figure 7B). Furthermore, the expression of the two later genes was checked in the black versus white varieties (Figure 7C). Interestingly, SiLAC1 exhibited a higher expression in the white seed variety (Zhongzhi No.1) at the early stage of the seed development (8 days after anthesis) compared to the black one (Zhongzhi No.33). Similarly, SiLAC39 expression was differentially higher in white seed variety (Zhongzhi No.1) at both 11 and 17 days after anthesis.
It is worth noting that all candidate genes fall into the principle according to which they should belong to the core gene repertoire of the sesame pan-genome and also, follow the pinoresinol synthesis kinetic, as described by Ono et al. [23] and depicted in Figure 8.
Briefly, the kinetic of the lignans biosynthesis suggests the expression of oxidative enzymes (potentially here peroxidase and/or laccase) at an upstream stage (from +(−)coniferyl alcohol to +(−)pinoresinol synthesis) of the lignans biosynthesis. Therefore, pinoresinol synthase is supposed to be expressed at an early stage of the seed development to enable downstream biosynthesis, resulting in +(−) sesamin, +(−) sesamolin, and +(−) sesaminol at the maturity stage (~30 days after anthesis). Thus, by combining the pinoresinol content extracted from sesame seed by Ono et al. [23] and the transcript data sets (Figure 8), we selected two peroxidase (SiPOD52, SiPOD63) and two laccase (SiLAC1, SiLAC39) genes. Globally, the FPKM values of the four candidate genes were declining at 23 days after anthesis (approximatively one week before the maturation stage) (Figure 8). The SiPOD50 gene was filtered out since it was constantly expressed at all development stages of the seed (Figure 7C), which does not match the lignan biosynthesis kinetic. Therefore, in downstream bioinformatic analyses (cis-acting elements, orthologs identification and gene ontology analyses), only SiPOD52, SiPOD63, SiLAC1, and SiLAC39 genes were selected.

3.3. Cis-Acting Elements, Related Transcription Factors, and Functional Attributes

Although the candidate genes shared similar gene structures (Figure 9A–C), a wide diversity of regulation and functional characteristics were highlighted through the cis-element analysis (Figure 9D), including hormonal response, light response, abiotic stress response, and physiological development. This suggests that SiPOD and SiLAC genes might be involved in a broad spectrum of biological processes in sesame.
Knowing the important role of transcription factors in gene regulation, we performed a transcription factor-oriented enrichment analysis by using all SiPOD and SiLAC genes to identify candidate transcription factors (TF) potentially involved in peroxidase and laccase gene regulation.
A panel of TF families was predicted, among which v-myb avian myeloblastosis viral oncogene homolog (MYB), NAM (no apical meristem), ATAF1-2 (Arabidopsis thaliana activating factor), CUC2 (cup-shaped cotyledon) (NAC), Basic leucine zipper (bZIP), Heat shock factors (HSF), Homeodomain-leucine zipper (HD-ZIP), MIKC-type MADS-box (MIKC_MADS) were the most abundant. As depicted in Figure 9E,F, MYB is the most predominant TFs indicating their putative regulatory role in the expression of both peroxidases and laccases. Using transgenic A. thaliana lines, Shen et al. [96] showed that the sweet cherry (Prunus avium cv. Hong Deng) R2R3 MYB was able to alleviate salt stress and provide anti-bacterial resistance through the activation of peroxidases and accumulation of anthocyanin. Furthermore, co-expression of A. thaliana laccases (lac4 or lac17) with MYB63 genes is known to rescue dwarfism in A. thaliana mutant lines [97].
Looking at the orthologs of the candidate genes in other taxa by a phylogenetic approach with SHOOT [72], we found homologous sequences in oil- and non-oil crops including Solanum lycopersicum, Arabidopsis thaliana, Brassica oleracea, Gossypium raimondii, Glycine max, Triticum aestivum, Oryza sativa, and Zea mays (Table S5).
The GO annotation supported the hypothesis that the all-candidate peroxidase genes are involved in the hydrogen peroxidase catabolic process (GO:0042744) with a heme-binding (GO:0020037) and peroxidase activity (GO:0004601) as main molecular functions. Regarding laccase, the GO annotation indicated that SiLAC1 and SiLAC39 may be related to the degradation of lignin (GO:0046274), with copper ion binding (GO:0005507) and oxidase activity (GO:0016491) as major molecular functions. From the GO annotation results, both SiPOD and SiLAC genes were predicted to have an oxidative role, which is a key requirement for the transformation of (+)- coniferyl alcohol into (+)- pinoresinol.
To assess the validity of the RNA-seq data, a qRT-PCR experiment was performed using the selected genes (Figure S7A). The trends of changes in the relative expression were highly consistent with those of the transcriptome sequencing data (R2 = 0.8505) (Figure S7B); indicating that the RNA-seq data were reliable.
As matter of fact, peroxidase and laccase genes were reported to be able to oxidize coniferyl and p-coumaryl alcohols, acting as catalysts during cell wall lignification in Zinnia elegans [98,99,100,101,102]. Therefore, the suggested genes are a valuable candidate for functional validation and ultimately, for usage in the pharmaceutical and food industries through bioengineering. Since sesame is recalcitrant to genetic transformation using classical methods [2], the hairy roots method might be a valuable alternative path for in vitro production of sesame lignans as demonstrated by Ogasawara et al. [103].

4. Conclusions

Candidate genes for peroxidase (83 genes) and laccase (48 genes) were identified in the genome of sesame. Gene count for both genes varied considerably between varieties. Taking advantage of a large panel of transcriptome data and stringent filtering, four genes (SiPOD52, SiPOD63, SiLAC1, and SiLAC39) were proposed as candidates for encoding (+)-pinoresinol synthase. The genes were predicted to interact with a wide range of transcription factors, indicating that they are involved in diverse physiological processes. The findings of this study open a way for functional investigation of the candidate genes and eventually for the genetic improvement of sesame regarding the synthesis of lignans.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/life12081200/s1, Figure S1: Circos plot showing paralogous peroxidase genes exhibiting segmental duplication in sesame; Figure S2: Circos plot showing paralogous laccase genes belonging to segmental duplication pattern; Figure S3: Circos plot showing synthenic peroxidase genes between Sesamum indicum and Arabidopsis thaliana; Figure S4: Circos plot showing synthenic laccase genes between Sesamum indicum and Arabidopsis thaliana; Figure S5: Motif structure of sesame peroxidase genes; Figure S6: Motif structure of sesame laccase genes; Figure S7: qRT-PCR results of the target genes and correlation with RNA-seq data; Table S1: Information relative to the SRA accessions of RNA sequencing data used in the present study; Table S2: List of primers used in the quantitative real time-PCR analysis; Table S3: List of identified peroxidase genes in Sesamum indicum var. Zhongzhi13; Table S4. List of identified laccase genes in Sesamum indicum var. Zhongzhi13, Table S5. List of candidate genes, functional attributes from TAIR database and their orthologous sequences.

Author Contributions

Conceptualization, K.D., P.K. and M.R.; methodology, Y.A.B.Z. and K.D.; software, Y.A.B.Z.; validation, Y.A.B.Z. and K.D.; formal analysis, Y.A.B.Z.; data curation, Y.A.B.Z.; writing—original draft preparation, Y.A.B.Z.; writing—review and editing, Y.A.B.Z., N.-J.C., Y.M., K.D., M.R. and P.K.; visualization, Y.A.B.Z.; supervision, K.D.; project administration, K.D., M.R. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Alexander von Humboldt, grant number: Ref 3.4-1202788-SEN-GF-P.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The peroxidase and laccase genes presented in this study are available in Tables S2 and S3. The NCBI SRA accessions are listed in Table S1. The sesame reference genome assembly and annotation of Zhongzhi13 are available under the NCBI project number PRJNA186669.

Acknowledgments

The authors would like to thank Julien Alban Nguinkal for his technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Murata, J.; Ono, E.; Yoroizuka, S.; Toyonaga, H.; Shiraishi, A.; Mori, S.; Tera, M.; Azuma, T.; Nagano, A.J.; Nakayasu, M.; et al. Oxidative rearrangement of (+)-sesamin by CYP92B14 co-generates twin dietary lignans in sesame. Nat. Commun. 2017, 8, 2115. [Google Scholar] [CrossRef] [PubMed]
  2. Andargie, M.; Vinas, M.; Rathgeb, A.; Möller, E.; Karlovsky, P. Lignans of Sesame (Sesamum indicum L.): A Comprehensive Review. Molecules 2021, 26, 883. [Google Scholar] [CrossRef] [PubMed]
  3. Katayama, S.; Sugiyama, H.; Kushimoto, S.; Uchiyama, Y.; Hirano, M.; Nakamura, S. Effects of Sesaminol Feeding on Brain Aβ Accumulation in a Senescence-Accelerated Mouse-Prone 8. J. Agric. Food Chem. 2016, 64, 4908–4913. [Google Scholar] [CrossRef] [PubMed]
  4. Shimoyoshi, S.; Takemoto, D.; Ono, Y.; Kitagawa, Y.; Shibata, H.; Tomono, S.; Unno, K.; Wakabayashi, K. Sesame lignans suppress age-related cognitive decline in senescence-accelerated mice. Nutrients 2019, 11, 1582. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, Z.; Saarinen, N.M.; Thompson, L.U. Sesamin Is One of the Major Precursors of Mammalian Lignans in Sesame Seed (Sesamum indicum) as Observed In Vitro and in Rats. J. Nutr. 2006, 136, 906–912. [Google Scholar] [CrossRef] [Green Version]
  6. Wu, M.-S.; Aquino, L.B.B.; Barbaza, M.Y.U.; Hsieh, C.-L.; De Castro-Cruz, K.A.; Yang, L.-L.; Tsai, P.-W. Anti-Inflammatory and Anticancer Properties of Bioactive Compounds from Sesamum indicum L.—A Review. Molecules 2019, 24, 4426. [Google Scholar] [CrossRef] [Green Version]
  7. Peñalvo, J.L.; López-Romero, P. Urinary enterolignan concentrations are positively associated with serum HDL cholesterol and negatively associated with serum triglycerides in U.S. adults. J. Nutr. 2012, 142, 751–756. [Google Scholar] [CrossRef] [Green Version]
  8. Sun, Q.; Wedick, N.M.; Pan, A.; Townsend, M.K.; Cassidy, A.; Franke, A.A.; Rimm, E.B.; Hu, F.B.; Van Dam, R.M. Gut microbiota metabolites of dietary lignans and risk of type 2 diabetes: A prospective investigation in two cohorts of U.S. women. Diabetes Care 2014, 37, 1287–1295. [Google Scholar] [CrossRef] [Green Version]
  9. Rodríguez-García, C.; Sánchez-Quesada, C.; Toledo, E.; Delgado-Rodríguez, M.; Gaforio, J. Naturally Lignan-Rich Foods: A Dietary Tool for Health Promotion? Molecules 2019, 24, 917. [Google Scholar] [CrossRef] [Green Version]
  10. Durazzo, A.; Carcea, M.; Adlercreutz, H.; Azzini, E.; Polito, A.; Olivieri, L.; Zaccaria, M.; Meneghini, C.; Maiani, F.; Bausano, G.; et al. Effects of consumption of whole grain foods rich in lignans in healthy postmenopausal women with moderate serum cholesterol: A pilot study. Int. J. Food Sci. Nutr. 2014, 65, 637–645. [Google Scholar] [CrossRef]
  11. Zamora-Ros, R.; Agudo, A.; Luján-Barroso, L.; Romieu, I.; Ferrari, P.; Knaze, V.; Bueno-de-Mesquita, H.B.; Leenders, M.; Travis, R.C.; Navarro, C.; et al. Dietary flavonoid and lignan intake and gastric adenocarcinoma risk in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Am. J. Clin. Nutr. 2012, 96, 1398–1408. [Google Scholar] [CrossRef] [PubMed]
  12. Peterson, J.; Dwyer, J.; Adlercreutz, H.; Scalbert, A.; Jacques, P.; McCullough, M.L. Dietary lignans: Physiology and potential for cardiovascular disease risk reduction. Nutr. Rev. 2010, 68, 571–603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Anonymous. Lignans Market Share, Size, Trends, Industry Analysis Report; Polaris Market Research: San Franscisco, CA, USA, 2021; pp. 1–145. [Google Scholar]
  14. Kojima, A.; Yuasa, I.; Kiyomoto, K.; Omura, A. Composition for Promoting Collagen Production, for Promoting Elastin Production and/or for Promoting Keratinocyte Migration and Usage Therefor. US Patent No. 9,629,823 B2, 25 April 2017. [Google Scholar]
  15. Sok, D.-E.; Cui, H.; Kim, M. Isolation and Boactivities of Furfuran Type Lignan Compounds from Edible Plants. Recent Pat. Food Nutr. Agric. 2012, 1, 87–95. [Google Scholar] [CrossRef] [Green Version]
  16. Yamada, D.; Kato, M.; Ono, Y.; Tomimori, N.; Nishiumi, T.; Nakahara, K. Oil-in-Water Emulsions Containing Lignan-Class Compounds and Compositions Containing the Same. U.S. Patent No. 8,685,455 B2, 1 April 2014. [Google Scholar]
  17. You, J.W.; Rho, H.S.; Kim, D.H.; Chang, I.S.; Lee, O.S. Sesamol Derivatives or Their Salts, the Process for Preparing the Same, and the Skin External Composition Containing the Same. U.S. Patent No. 7,943,599 B2, 17 May 2011. [Google Scholar]
  18. Forse, A.R.; Chavali, S.R. Sesamol Inhibition of Δ-5-Desaturase Activity and Uses Therefor. U.S. Patent No. 6,172,106 B1, 9 January 2001. [Google Scholar]
  19. Markus, R.L. Process for Producing Sesamol. U.S. Patent No. 3,058,995, 16 October 1962. [Google Scholar]
  20. Hardwicke, J.E.; King, J.; Terrell, R.C. Synthesis of Sesamol Acetate and Sesamol. U.S. Patent No. 2,885,407, 5 May 1959. [Google Scholar]
  21. Chami, A.; Pillai, C.; Balachandran, M.; Rhesma, V.; Sundaresan, A.; Thomas, S.; Sukumar, D.; Saritha, S.K.S. Process for the Extraction of Bioactive Lignans with High Yield and Purity from Sesame Oil. U.S. Patent No. 8,350,066 B2, 8 January 2013. [Google Scholar]
  22. Namiki, M.; Kobayashi, T.; Hara, H. Process of Producing Sesame Lignans an/or Sasame Flavors. U.S. Patent No. 6,278,005 B1, 21 August 2001. [Google Scholar]
  23. Ono, E.; Nakai, M.; Fukui, Y.; Tomimori, N.; Fukuchi-Mizutani, M.; Saito, M.; Satake, H.; Tanaka, T.; Katsuta, M.; Umezawa, T.; et al. Formation of two methylenedioxy bridges by a Sesamum CYP81Q protein yielding a furofuran lignan, (+)-sesamin. Proc. Natl. Acad. Sci. USA 2006, 103, 10116–10121. [Google Scholar] [CrossRef] [Green Version]
  24. Noguchi, A.; Fukui, Y.; Iuchi-Okada, A.; Kakutani, S.; Satake, H.; Iwashita, T.; Nakao, M.; Umezawa, T.; Ono, E. Sequential glucosylation of a furofuran lignan, (+)-sesaminol, by Sesamum indicum UGT71A9 and UGT94D1 glucosyltransferases. Plant J. 2008, 54, 415–427. [Google Scholar] [CrossRef]
  25. Ono, E.; Waki, T.; Oikawa, D.; Murata, J.; Shiraishi, A.; Toyonaga, H.; Kato, M.; Ogata, N.; Takahashi, S.; Yamaguchi, M.; et al. Glycoside-specific glycosyltransferases catalyze regio-selective sequential glucosylations for a sesame lignan, sesaminol triglucoside. Plant J. 2019, 101, 1221–1233. [Google Scholar] [CrossRef]
  26. Harada, E.; Murata, J.; Ono, E.; Toyonaga, H.; Shiraishi, A.; Hideshima, K.; Yamamoto, M.P.; Horikawa, M. (+)-Sesamin-oxidising CYP92B14 shapes specialised lignan metabolism in sesame. Plant J. 2020, 104, 1117–1128. [Google Scholar] [CrossRef]
  27. Kim, K.W.; Smith, C.A.; Daily, M.D.; Cort, J.R.; Davin, L.B.; Lewis, N.G. Trimeric structure of (+)-pinoresinol-forming dirigent protein at 1.95 Å resolution with three isolated active sites. J. Biol. Chem. 2015, 290, 1308–1318. [Google Scholar] [CrossRef] [Green Version]
  28. Kim, M.K.; Jeon, J.H.; Fujita, M.; Davin, L.B.; Lewis, N.G. The western red cedar (Thuja plicata) 8-8′ DIRIGENT family displays diverse expression patterns and conserved monolignol coupling specificity. Plant Mol. Biol. 2002, 49, 199–214. [Google Scholar] [CrossRef]
  29. Davin, L.B.; Bin Wang, H.; Crowell, A.L.; Bedgar, D.L.; Martin, D.M.; Sarkanen, S.; Lewis, N.G. Stereoselective bimolecular phenoxy radical coupling by an auxiliary (dirigent) protein without an active center. Science 1997, 275, 362–366. [Google Scholar] [CrossRef]
  30. Pickel, B.; Constantin, M.-A.; Pfannstiel, J.; Conrad, J.; Beifuss, U.; Schaller, A. An Enantiocomplementary Dirigent Protein for the Enantioselective Laccase-Catalyzed Oxidative Coupling of Phenols. Angew. Chem. Int. Ed. 2010, 49, 202–204. [Google Scholar] [CrossRef] [PubMed]
  31. Bindschedler, L.V.; Dewdney, J.; Blee, K.A.; Stone, J.M.; Asai, T.; Plotnikov, J.; Denoux, C.; Hayes, T.; Gerrish, C.; Davies, D.R.; et al. Peroxidase-dependent apoplastic oxidative burst in Arabidopsis required for pathogen resistance. Plant J. 2006, 47, 851–863. [Google Scholar] [CrossRef] [Green Version]
  32. Hu, Q.; Min, L.; Yang, X.; Jin, S.; Zhang, L.; Li, Y.; Ma, Y.; Qi, X.; Li, D.; Liu, H.; et al. Laccase GhLac1 modulates broad-spectrum biotic stress tolerance via manipulating phenylpropanoid pathway and jasmonic acid synthesis. Plant Physiol. 2018, 176, 1808–1823. [Google Scholar] [CrossRef] [Green Version]
  33. Daudi, A.; Cheng, Z.; O’Brien, J.A.; Mammarella, N.; Khan, S.; Ausubel, F.M.; Paul Bolwell, G. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell 2012, 24, 275–287. [Google Scholar] [CrossRef] [Green Version]
  34. Kidwai, M.; Ahmad, I.Z.; Chakrabarty, D. Class III peroxidase: An indispensable enzyme for biotic/abiotic stress tolerance and a potent candidate for crop improvement. Plant Cell Rep. 2020, 39, 1381–1393. [Google Scholar] [CrossRef]
  35. Roy, J.; Blervacq, A.S.; Créach, A.; Huss, B.; Hawkins, S.; Neutelings, G. Spatial regulation of monolignol biosynthesis and laccase genes control developmental and stress-related lignin in flax. BMC Plant Biol. 2017, 17, 124. [Google Scholar] [CrossRef] [Green Version]
  36. Aleem, M.; Riaz, A.; Raza, Q.; Aleem, M.; Aslam, M.; Kong, K.; Atif, R.M.; Kashif, M.; Bhat, J.A.; Zhao, T. Genome-wide characterization and functional analysis of class III peroxidase gene family in soybean reveal regulatory roles of GsPOD40 in drought tolerance. Genomics 2022, 114, 45–60. [Google Scholar] [CrossRef]
  37. Cai, X.; Davis, E.J.; Ballif, J.; Liang, M.; Bushman, E.; Haroldsen, V.; Torabinejad, J.; Wu, Y. Mutant identification and characterization of the laccase gene family in Arabidopsis. J. Exp. Bot. 2006, 57, 2563–2569. [Google Scholar] [CrossRef] [Green Version]
  38. Li, W.; Huai, X.; Li, P.; Raza, A.; Mubarik, M.S.; Habib, M.; Faiz, S.; Zhang, B.; Pan, J.; Khan, R.S.A. Genome-wide characterization of glutathione peroxidase (GPX) gene family in rapeseed (brassica napus L.) revealed their role in multiple abiotic stress response and hormone signaling. Antioxidants 2021, 10, 1481. [Google Scholar] [CrossRef]
  39. Hu, Q.; Xiao, S.; Guan, Q.; Tu, L.; Sheng, F.; Du, X.; Zhang, X. The laccase gene GhLac1 modulates fiber initiation and elongation by coordinating jasmonic acid and flavonoid metabolism. Crop J. 2020, 8, 522–533. [Google Scholar] [CrossRef]
  40. Jackson, P.; Ricardo, C.P.P. The changing peroxidase polymorphism in Lupinus albus during vegetative development. Funct. Plant Biol. 1998, 25, 261. [Google Scholar] [CrossRef]
  41. Liang, M.; Davis, E.; Gardner, D.; Cai, X.; Wu, Y. Involvement of AtLAC15 in lignin synthesis in seeds and in root elongation of Arabidopsis. Planta 2006, 224, 1185–1196. [Google Scholar] [CrossRef] [PubMed]
  42. Gabaldón, C.; López-Serrano, M.; Pedreño, M.A.; Barceló, A.R. Cloning and molecular characterization of the basic peroxidase isoenzyme from Zinnia elegans, an enzyme involved in lignin biosynthesis. Plant Physiol. 2005, 139, 1138–1154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Zhang, Y.C.; Yu, Y.; Wang, C.Y.; Li, Z.Y.; Liu, Q.; Xu, J.; Liao, J.Y.; Wang, X.J.; Qu, L.H.; Chen, F.; et al. Overexpression of microRNA OsmiR397 improves rice yield by increasing grain size and promoting panicle branching. Nat. Biotechnol. 2013, 31, 848–852. [Google Scholar] [CrossRef]
  44. Fang, F.; Zhang, X.L.; Luo, H.H.; Zhou, J.J.; Gong, Y.H.; Li, W.J.; Shi, Z.W.; He, Q.; Wu, Q.; Li, L.; et al. An intracellular Laccase is responsible for epicatechin-mediated Anthocyanin degradation in litchi fruit Pericarp. Plant Physiol. 2015, 169, 2391–2408. [Google Scholar] [CrossRef] [Green Version]
  45. Ring, L.; Yeh, S.Y.; Hücherig, S.; Hoffmann, T.; Blanco-Portales, R.; Fouche, M.; Villatoro, C.; Denoyes, B.; Monfort, A.; Caballero, J.L.; et al. Metabolic interaction between anthocyanin and lignin biosynthesis is associated with peroxidase FaPRX27 in strawberry fruit. Plant Physiol. 2013, 163, 43–60. [Google Scholar] [CrossRef] [Green Version]
  46. Pourcel, L.; Routaboul, J.M.; Cheynier, V.; Lepiniec, L.; Debeaujon, I. Flavonoid oxidation in plants: From biochemical properties to physiological functions. Trends Plant Sci. 2007, 12, 29–36. [Google Scholar] [CrossRef]
  47. Yoshida, K.; Kaothien, P.; Matsui, T.; Kawaoka, A.; Shinmyo, A. Molecular biology and application of plant peroxidase genes. Appl. Microbiol. Biotechnol. 2003, 60, 665–670. [Google Scholar] [CrossRef]
  48. Poppius-Levlin, K.; Tamminen, T.; Kalliola, A.; Ohra-aho, T. Characterization of Residual Lignins in Pulps Delignified by Laccase/N -Hydroxyacetanilide; American Chemical Society: Washington, DC, USA, 2001; pp. 358–372. [Google Scholar]
  49. Hussain, A.; Bilal, M.; Rafeeq, H.; Jabeen, Z.; Afsheen, N.; Sher, F.; Kumar, V.; Bharagava, R.N.; Ferreira, L.F.R.; Iqbal, H.M.N. Role of laccase in the pulp and paper industry. In Nanotechnology in Paper and Wood Engineering; Elsevier: Amsterdam, The Netherlands, 2022; pp. 35–60. ISBN 9780323858359. [Google Scholar]
  50. Mayer, A.M.; Staples, R.C. Laccase: New functions for an old enzyme. Phytochemistry 2002, 60, 551–565. [Google Scholar] [CrossRef]
  51. Mistry, J.; Bateman, A.; Finn, R.D. Predicting active site residue annotations in the Pfam database. BMC Bioinform. 2007, 8, 298. [Google Scholar] [CrossRef] [Green Version]
  52. Jones, P.; Binns, D.; Chang, H.Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5: Genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Emms, D.M.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 2019, 20, 238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chao, J.; Li, Z.; Sun, Y.; Aluko, O.O.; Wu, X.; Wang, Q.; Liu, G. MG2C: A user-friendly online tool for drawing genetic maps. Mol. Hortic. 2021, 1, 16. [Google Scholar] [CrossRef]
  55. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  56. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef]
  58. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  59. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, L.; Yu, S.; Tong, C.; Zhao, Y.; Liu, Y.; Song, C.; Zhang, Y.; Zhang, X.; Wang, Y.; Hua, W.; et al. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. Genome Biol. 2014, 15, R39. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, L.; Dossou, S.S.K.; Wei, X.; Zhang, Y.; Li, D.; Yu, J.; Zhang, X. Transcriptome dynamics during black and white sesame (Sesamum indicum L.) seed development and identification of candidate genes associated with black pigmentation. Genes 2020, 11, 1399. [Google Scholar] [CrossRef]
  62. Wang, L.; Zhang, Y.; Li, D.; Dossa, K.; Wang, M.L.; Zhou, R.; Yu, J.; Zhang, X. Gene expression profiles that shape high and low oil content sesames. BMC Genet. 2019, 20, 45. [Google Scholar] [CrossRef] [PubMed]
  63. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data [Online]. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/2010 (accessed on 11 December 2021).
  64. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  68. Wang, L.; Zhang, Y.; Li, P.; Wang, X.; Zhang, W.; Wei, W.; Zhang, X. HPLC Analysis of Seed Sesamin and Sesamolin Variation in a Sesame Germplasm Collection in China. J. Am. Oil Chem. Soc. 2012, 89, 1011–1020. [Google Scholar] [CrossRef]
  69. Ghotbzadeh Kermani, S.; Saeidi, G.; Sabzalian, M.R.; Gianinetti, A. Drought stress influenced sesamin and sesamolin content and polyphenolic components in sesame (Sesamum indicum L.) populations with contrasting seed coat colors. Food Chem. 2019, 289, 360–368. [Google Scholar] [CrossRef]
  70. Bailey, T.L.; Johnson, J.; Grant, C.E.; Noble, W.S. The MEME Suite. Nucleic Acids Res. 2015, 43, W39–W49. [Google Scholar] [CrossRef] [Green Version]
  71. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van De Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  72. Emms, D.M.; Kelly, S. SHOOT: Phylogenetic gene search and ortholog inference. Genome Biol. 2022, 23, 85. [Google Scholar] [CrossRef]
  73. Tian, F.; Yang, D.C.; Meng, Y.Q.; Jin, J.; Gao, G. PlantRegMap: Charting functional regulatory maps in plants. Nucleic Acids Res. 2020, 48, D1104–D1113. [Google Scholar] [CrossRef]
  74. Dossa, K.; Wei, X.; Li, D.; Fonceka, D.; Zhang, Y.; Wang, L.; Yu, J.; Boshou, L.; Diouf, D.; Cissé, N.; et al. Insight into the AP2/ERF transcription factor superfamily in sesame and expression profiling of DREB subfamily under drought stress. BMC Plant Biol. 2016, 16, 171. [Google Scholar] [CrossRef] [Green Version]
  75. Audic, S.; Claverie, J.-M. The Significance of Digital Gene Expression Profiles. Genome Res. 1997, 7, 986–995. [Google Scholar] [CrossRef]
  76. Passardi, F.; Longet, D.; Penel, C.; Dunand, C. The class III peroxidase multigenic family in rice and its evolution in land plants. Phytochemistry 2004, 65, 1879–1893. [Google Scholar] [CrossRef]
  77. Wang, Y.; Wang, Q.; Zhao, Y.; Han, G.; Zhu, S. Systematic analysis of maize class III peroxidase gene family reveals a conserved subfamily involved in abiotic stress response. Gene 2015, 566, 95–108. [Google Scholar] [CrossRef]
  78. Yang, X.; Yuan, J.; Luo, W.; Qin, M.; Yang, J.; Wu, W.; Xie, X. Genome-Wide Identification and Expression Analysis of the Class III Peroxidase Gene Family in Potato (Solanum tuberosum L.). Front. Genet. 2020, 11, 593577. [Google Scholar] [CrossRef]
  79. Cai, K.; Liu, H.; Chen, S.; Liu, Y.; Zhao, X.; Chen, S. Genome-wide identification and analysis of class III peroxidases in Betula pendula. BMC Genom. 2021, 22, 314. [Google Scholar] [CrossRef]
  80. Tognolli, M.; Penel, C.; Greppin, H.; Simon, P. Analysis and expression of the class III peroxidase large gene family in Arabidopsis thaliana. Gene 2002, 288, 129–138. [Google Scholar] [CrossRef]
  81. Xiao, H.; Wang, C.; Khan, N.; Chen, M.; Fu, W.; Guan, L.; Leng, X. Genome-wide identification of the class III POD gene family and their expression profiling in grapevine (Vitis vinifera L). BMC Genom. 2020, 21, 444. [Google Scholar] [CrossRef]
  82. Qui, K.; Zhou, H.; Pan, H.; Sheng, Y.; Yu, H.; Xie, Q.; Chen, H.; Cai, Y.; Zhang, J.; He, J. Genome-wide identification and functional analysis of the peach (P. persica) laccase gene family reveal members potentially involved in endocarp lignification. Trees-Struct. Funct. 2022. [Google Scholar] [CrossRef]
  83. Li, R.; Zhao, Y.; Sun, Z.; Wu, Z.; Wang, H.; Fu, C. Genome-Wide Identification of Switchgrass Laccases Involved in Lignin Biosynthesis and Heavy-Metal Responses. Int. J. Mol. Sci. 2022, 23, 6530. [Google Scholar] [CrossRef]
  84. Wan, F.; Zhang, L.; Tan, M.; Wang, X.; Wang, G.L.; Qi, M.; Liu, B.; Gao, J.; Pan, Y.; Wang, Y. Genome-wide identification and characterization of laccase family members in eggplant (Solanum melongena L.). PeerJ 2022, 10, e12922. [Google Scholar] [CrossRef]
  85. Turlapati, P.V.; Kim, K.W.; Davin, L.B.; Lewis, N.G. The laccase multigene family in Arabidopsis thaliana: Towards addressing the mystery of their gene function(s). Planta 2011, 233, 439–470. [Google Scholar] [CrossRef]
  86. Yu, J.; Golicz, A.A.; Lu, K.; Dossa, K.; Zhang, Y.; Chen, J.; Wang, L.; You, J.; Fan, D.; Edwards, D.; et al. Insight into the evolution and functional characteristics of the pan-genome assembly from sesame landraces and modern cultivars. Plant Biotechnol. J. 2019, 17, 881–892. [Google Scholar] [CrossRef] [Green Version]
  87. Moore, R.C.; Purugganan, M.D. The early stages of duplicate gene evolution. Proc. Natl. Acad. Sci. USA 2003, 100, 15682–15687. [Google Scholar] [CrossRef] [Green Version]
  88. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [Green Version]
  89. Wu, C.; Ding, X.; Ding, Z.; Tie, W.; Yan, Y.; Wang, Y.; Yang, H.; Hu, W. The class III peroxidase (POD) gene family in cassava: Identification, phylogeny, duplication, and expression. Int. J. Mol. Sci. 2019, 20, 2730. [Google Scholar] [CrossRef] [Green Version]
  90. Meng, G.; Fan, W.; Rasmussen, S.K. Characterisation of the class III peroxidase gene family in carrot taproots and its role in anthocyanin and lignin accumulation. Plant Physiol. Biochem. 2021, 167, 245–256. [Google Scholar] [CrossRef]
  91. Duan, P.; Wang, G.; Chao, M.; Zhang, Z.; Zhang, B. Genome-wide identification and analysis of class III peroxidases in allotetraploid cotton (Gossypium hirsutum L.) and their responses to PK deficiency. Genes 2019, 10, 473. [Google Scholar] [CrossRef] [Green Version]
  92. Yang, T.; Zhang, P.; Pan, J.; Amanullah, S.; Luan, F.; Han, W.; Liu, H.; Wang, X. Genome-Wide Analysis of the Peroxidase Gene Family and Verification of Lignin Synthesis-Related Genes in Watermelon. Int. J. Mol. Sci. 2022, 23, 642. [Google Scholar] [CrossRef]
  93. Huang, Y.-L.; Zhang, L.-K.; Zhang, K.; Chen, S.-M.; Hu, J.-B.; Cheng, F. The impact of tandem duplication on gene evolution in Solanaceae species. J. Integr. Agric. 2022, 21, 1004–1014. [Google Scholar] [CrossRef]
  94. Cao, Y.; Han, Y.; Meng, D.; Li, D.; Jin, Q.; Lin, Y.; Cai, Y. Structural, evolutionary, and functional analysis of the class III peroxidase gene family in Chinese pear (pyrus bretschneideri). Front. Plant Sci. 2016, 7, 1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Ren, L.L.; Liu, Y.J.; Liu, H.J.; Qian, T.T.; Qi, L.W.; Wang, X.R.; Zeng, Q.Y. Subcellular relocalization and positive selection play key roles in the retention of duplicate genes of Populus class III peroxidase family. Plant Cell 2014, 26, 2404–2419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Shen, X.; Guo, X.; Guo, X.; Zhao, D.; Zhao, W.; Chen, J.; Li, T. PacMYBA, a sweet cherry R2R3-MYB transcription factor, is a positive regulator of salt stress tolerance and pathogen resistance. Plant Physiol. Biochem. 2017, 112, 302–311. [Google Scholar] [CrossRef] [PubMed]
  97. Perkins, M.L.; Schuetz, M.; Unda, F.; Smith, R.A.; Sibout, R.; Hoffmann, N.J.; Wong, D.C.J.; Castellarin, S.D.; Mansfield, S.D.; Samuels, L. Dwarfism of high-monolignol Arabidopsis plants is rescued by ectopic LACCASE overexpression. Plant Direct 2020, 4, e00265. [Google Scholar] [CrossRef] [PubMed]
  98. Barceló, A.R.; Ros, L.V.G.; Carrasco, A.E. Looking for syringyl peroxidases. Trends Plant Sci. 2007, 12, 486–491. [Google Scholar] [CrossRef]
  99. Ros Barceló, A.; Gómez Ros, L.V.; Gabaldón, C.; López-Serrano, M.; Pomar, F.; Carrión, J.S.; Pedreño, M.A. Basic peroxidases: The gateway for lignin evolution? Phytochem. Rev. 2004, 3, 61–78. [Google Scholar] [CrossRef]
  100. Novo-Uzal, E.; Fernández-Pérez, F.; Herrero, J.; Gutiérrez, J.; Gómez-Ros, L.V.; Bernal, M.Á.; Díaz, J.; Cuello, J.; Pomar, F.; Pedreño, M.Á. From Zinnia to Arabidopsis: Approaching the involvement of peroxidases in lignification. J. Exp. Bot. 2013, 64, 3499–3518. [Google Scholar] [CrossRef] [Green Version]
  101. Tugbaeva, A.; Ermoshin, A.; Plotnikov, D.; Wuriyanghan, H.; Kiseleva, I. Role of Class III Peroxidases in Stem Lignification of Zinnia elegans Jacq. Life Sci. Forum 2021, 4, 22. [Google Scholar] [CrossRef]
  102. Gavnholt, B.; Larsen, K. Molecular biology of plant laccases in relation to lignin formation. Physiol. Plant. 2002, 116, 273–280. [Google Scholar] [CrossRef]
  103. Ogasawara, T.; Chiba, K.; Tada, M. Sesamum indicum L. (Sesame): In Vitro Culture, and the Production of Naphthoquinone and Other Secondary Metabolites. In Medicinal and Aromatic Plants X; Springer: Berlin/Heidelberg, Germany, 1998; Volume 41, pp. 366–393. [Google Scholar]
Life 12 01200 g001 550
Figure 1. A simplified lignan biosynthesis pathway in sesame showing the route of synthesis of (+)- sesamin, (+)- sesamolin, and (+)- sesaminol. CYP81Q1 and CYP92B14 triggered the biosynthesis of (+)- sesamin, (+)- sesamolin and (+)- sesaminol, while UTG71A9, UGT94D1, UGT94AA2, and UGT94AG1 are suggested to catalyze the synthesis of (+)- sesaminol monoglucoside, (+)-sesaminol 2-O-β-D-glucosyl (1 2)-O-β-D-glucoside, (+)-sesaminol 2-O-β-D-glucosyl (1 6)-O-β-D-glucoside, and (+)-sesaminol triglucoside. The target step in the present study is marked by a yellow rectangle. The pathway is adapted from Ono et al. [23].
Figure 1. A simplified lignan biosynthesis pathway in sesame showing the route of synthesis of (+)- sesamin, (+)- sesamolin, and (+)- sesaminol. CYP81Q1 and CYP92B14 triggered the biosynthesis of (+)- sesamin, (+)- sesamolin and (+)- sesaminol, while UTG71A9, UGT94D1, UGT94AA2, and UGT94AG1 are suggested to catalyze the synthesis of (+)- sesaminol monoglucoside, (+)-sesaminol 2-O-β-D-glucosyl (1 2)-O-β-D-glucoside, (+)-sesaminol 2-O-β-D-glucosyl (1 6)-O-β-D-glucoside, and (+)-sesaminol triglucoside. The target step in the present study is marked by a yellow rectangle. The pathway is adapted from Ono et al. [23].
Life 12 01200 g001
Life 12 01200 g002 550
Figure 2. Gene count and conservation of peroxidase (A) and laccase (B) genes in sesame pangenome. Horizontal bar charts summarize the gene count. Upset plots show the core conserved count of peroxidases/laccases within sesame pan-genome.
Figure 2. Gene count and conservation of peroxidase (A) and laccase (B) genes in sesame pangenome. Horizontal bar charts summarize the gene count. Upset plots show the core conserved count of peroxidases/laccases within sesame pan-genome.
Life 12 01200 g002
Life 12 01200 g003 550
Figure 3. Chromosome location of peroxidase and laccase genes in sesame.
Figure 3. Chromosome location of peroxidase and laccase genes in sesame.
Life 12 01200 g003
Life 12 01200 g004 550
Figure 4. Unrooted maximum likelihood phylogenetic tree inferred from 73 A. thaliana and 83 S. indicum peroxidase genes. The tree was constructed following the LG+R5 model with a total of 1000 iterations. Square symbol represents genes from S. indicum while triangle symbol stands for A. thaliana genes. Yellow (41–80), and green (81–100) dot represents the clades support values. Arabidopsis-specific phylogenetic groups were designated A1 to A3. Groups containing SiPOD are designated G1 to G9.
Figure 4. Unrooted maximum likelihood phylogenetic tree inferred from 73 A. thaliana and 83 S. indicum peroxidase genes. The tree was constructed following the LG+R5 model with a total of 1000 iterations. Square symbol represents genes from S. indicum while triangle symbol stands for A. thaliana genes. Yellow (41–80), and green (81–100) dot represents the clades support values. Arabidopsis-specific phylogenetic groups were designated A1 to A3. Groups containing SiPOD are designated G1 to G9.
Life 12 01200 g004
Life 12 01200 g005 550
Figure 5. Unrooted maximum likelihood phylogenetic tree inferred from 17 A. thaliana and 48 S. indicum laccase genes. The tree was constructed following the LG+I+G4 model with a total of 1000 iterations. Square shape represents the genes from S. indicum while triangle shape stands for A. thaliana genes. Yellow (41–80), and green (81–100) dot represents the clades support values. An Arabidopsis-specific phylogenetic group was designated A1. Groups containing SiLAC are designated G1 to G7.
Figure 5. Unrooted maximum likelihood phylogenetic tree inferred from 17 A. thaliana and 48 S. indicum laccase genes. The tree was constructed following the LG+I+G4 model with a total of 1000 iterations. Square shape represents the genes from S. indicum while triangle shape stands for A. thaliana genes. Yellow (41–80), and green (81–100) dot represents the clades support values. An Arabidopsis-specific phylogenetic group was designated A1. Groups containing SiLAC are designated G1 to G7.
Life 12 01200 g005
Life 12 01200 g006 550
Figure 6. Expression profile of sesame peroxidase genes in different tissues: (A) seed, root, leaf, stem and capsule from Zhongzhi13; (B) Seed from ZZM4728 (HOS), ZZM3495 (LOSA), and ZZM2161 (LOSB); (C) Seed from black seed Zhongzhi No. 33 (B) and white seed Zhongzhi No.1 (W) varieties. Candidate genes for each set of tissues are highlighted with a red star symbol.
Figure 6. Expression profile of sesame peroxidase genes in different tissues: (A) seed, root, leaf, stem and capsule from Zhongzhi13; (B) Seed from ZZM4728 (HOS), ZZM3495 (LOSA), and ZZM2161 (LOSB); (C) Seed from black seed Zhongzhi No. 33 (B) and white seed Zhongzhi No.1 (W) varieties. Candidate genes for each set of tissues are highlighted with a red star symbol.
Life 12 01200 g006
Life 12 01200 g007 550
Figure 7. Expression profile of sesame laccase genes in different tissues: (A) seed, root, leaf, stem and capsule from Zhongzhi13; (B) Seed from ZZM4728 (HOS), ZZM3495 (LOSA), and ZZM2161 (LOSB); (C) Seed from black seed Zhongzhi No. 33 (B) and white seed Zhongzhi No.1 (W) varieties. Candidate genes for each set of tissues are highlighted with a red star symbol.
Figure 7. Expression profile of sesame laccase genes in different tissues: (A) seed, root, leaf, stem and capsule from Zhongzhi13; (B) Seed from ZZM4728 (HOS), ZZM3495 (LOSA), and ZZM2161 (LOSB); (C) Seed from black seed Zhongzhi No. 33 (B) and white seed Zhongzhi No.1 (W) varieties. Candidate genes for each set of tissues are highlighted with a red star symbol.
Life 12 01200 g007
Life 12 01200 g008 550
Figure 8. FPKM variation of candidate peroxidase and laccase genes following sesame seed development and pinoresinol content in white seed genotype Zhongfengzhi No1. Pinoresinol content values were obtained from Ono et al. [23] study. SiPOD52 (A), SiLAC39 (B), SiPOD63 (C), and SiLAC1 (D) variations following pinoresinol content were depicted during the seed development stages starting from five days after anthesis (W5D) to 30 days after anthesis (W30D). Pinoresinol content was colored in red. Peroxidase and laccase gene expression counts were colored in purple and green, respectively.
Figure 8. FPKM variation of candidate peroxidase and laccase genes following sesame seed development and pinoresinol content in white seed genotype Zhongfengzhi No1. Pinoresinol content values were obtained from Ono et al. [23] study. SiPOD52 (A), SiLAC39 (B), SiPOD63 (C), and SiLAC1 (D) variations following pinoresinol content were depicted during the seed development stages starting from five days after anthesis (W5D) to 30 days after anthesis (W30D). Pinoresinol content was colored in red. Peroxidase and laccase gene expression counts were colored in purple and green, respectively.
Life 12 01200 g008
Life 12 01200 g009 550
Figure 9. Bioinformatic analysis of candidate peroxidase and laccase genes. Protein motifs (A), domain (B) and gene structure (C), cis-acting elements functional attributes (D), and transcription factor enrichment word clouds for laccase (E) and peroxidase (F) candidate genes.
Figure 9. Bioinformatic analysis of candidate peroxidase and laccase genes. Protein motifs (A), domain (B) and gene structure (C), cis-acting elements functional attributes (D), and transcription factor enrichment word clouds for laccase (E) and peroxidase (F) candidate genes.
Life 12 01200 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zoclanclounon, Y.A.B.; Rostás, M.; Chung, N.-J.; Mo, Y.; Karlovsky, P.; Dossa, K. Characterization of Peroxidase and Laccase Gene Families and In Silico Identification of Potential Genes Involved in Upstream Steps of Lignan Formation in Sesame. Life 2022, 12, 1200. https://0-doi-org.brum.beds.ac.uk/10.3390/life12081200

AMA Style

Zoclanclounon YAB, Rostás M, Chung N-J, Mo Y, Karlovsky P, Dossa K. Characterization of Peroxidase and Laccase Gene Families and In Silico Identification of Potential Genes Involved in Upstream Steps of Lignan Formation in Sesame. Life. 2022; 12(8):1200. https://0-doi-org.brum.beds.ac.uk/10.3390/life12081200

Chicago/Turabian Style

Zoclanclounon, Yedomon Ange Bovys, Michael Rostás, Nam-Jin Chung, Youngjun Mo, Petr Karlovsky, and Komivi Dossa. 2022. "Characterization of Peroxidase and Laccase Gene Families and In Silico Identification of Potential Genes Involved in Upstream Steps of Lignan Formation in Sesame" Life 12, no. 8: 1200. https://0-doi-org.brum.beds.ac.uk/10.3390/life12081200

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop