Next Article in Journal
A Comparative Transcriptome Analysis Reveals Physiological Maturation Properties of Mycelia in Pleurotus tuoliensis
Previous Article in Journal
Overcoming Fish Defences: The Virulence Factors of Yersinia ruckeri
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 9
10.3390/genes10090701
genes-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

The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection

by
Quan Wang
1,2,
Guang Li
1,2,
Kaijie Zheng
1,
Xiaobin Zhu
3,
Jingjing Ma
1,2,
Dongmei Wang
1,2,
Kuanqiang Tang
1,2,
Xingxing Feng
1,2,
Jiantian Leng
1,
Hui Yu
1,
Suxin Yang
1 and
Xianzhong Feng
1,*
1
Key Laboratory of Soybean Molecular Design Breeding, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China
2
School of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
3
School of Life Science, Jilin Agricultural University, Changchun 130118, China
*
Author to whom correspondence should be addressed.
Genes 2019, 10(9), 701; https://doi.org/10.3390/genes10090701
Submission received: 8 July 2019 / Revised: 18 August 2019 / Accepted: 9 September 2019 / Published: 11 September 2019
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Laccase is a widely used industrial oxidase for food processing, dye synthesis, paper making, and pollution remediation. At present, laccases used by industries come mainly from fungi. Plants contain numerous genes encoding laccase enzymes that show properties which are distinct from that of the fungal laccases. These plant-specific laccases may have better potential for industrial purposes. The aim of this work was to conduct a genome-wide search for the soybean laccase genes and analyze their characteristics and specific functions. A total of 93 putative laccase genes (GmLac) were identified from the soybean genome. All 93 GmLac enzymes contain three typical Cu-oxidase domains, and they were classified into five groups based on phylogenetic analysis. Although adjacent members on the tree showed highly similar exon/intron organization and motif composition, there were differences among the members within a class for both conserved and differentiated functions. Based on the expression patterns, some members of laccase were expressed in specific tissues/organs, while some exhibited a constitutive expression pattern. Analysis of the transcriptome revealed that some laccase genes might be involved in providing resistance to oomycetes. Analysis of the selective pressures acting on the laccase gene family in the process of soybean domestication revealed that 10 genes could have been under artificial selection during the domestication process. Four of these genes may have contributed to the transition of the soft and thin stem of wild soybean species into strong, thick, and erect stems of the cultivated soybean species. Our study provides a foundation for future functional studies of the soybean laccase gene family.

1. Introduction

Laccases, p-diphenol:oxygen oxidoreductases (EC 1.10.3.2), represent the largest subgroup of multicopper-containing oxidases (MCOs), also known as “blue” oxidases [1,2,3]. They are widespread not only in fungi and plants, but also in bacteria and insects [4,5,6,7]. Laccase was initially discovered in Japanese lacquer tree by Yoshida in 1883 [8]. Subsequently, it was defined as a metal-containing oxidase [9]. Commonly, the laccase protein sequence has three conserved domains (Cu-oxidase, Cu-oxidase_2, and Cu-oxidase_3) that are used to identify canonical laccases. The active center of laccases contains four copper atoms: a type-1 copper, a type-2 copper, and two type-3 copper atoms, which are responsible for electron transfer in the catalytic process [10,11]. The catalytic capacity of laccases is actually nonspecific. In most cases, laccase substrates are aromatic compounds; however, it also has the ability to catalytically oxidize non-aromatic compounds [12,13,14,15]. Precisely, because of their nonspecific catalytic ability, laccases were reported to be involved in various physiological activities in living organisms and are utilized in diverse biotechnological applications in industrial production.
Laccases used to be classified as fungal, plant, or bacterial laccases [2,16]. Basic characteristics and functions of laccases are diverse both within and across kingdoms. Delignification is probably the most studied function of fungal laccases [17], as it can be identified in nearly all wood-decomposing fungi [18], although the detailed mechanism of this process is still unclear [19,20]. In addition to the delignification function, fungal laccases are also connected with many other physiological processes including the regulation of morphology, such as fruiting body formation [21,22,23], pathogenicity when the fungi are plant pathogens [24,25], promotion of special compound synthesis and metabolism such as pigment formation [26], and flavonoid metabolism [27].
Unlike fungal laccases, plant laccases have a higher glycosylation degree, which is closely related to secretion, thermal stability, and enzyme activity [28]. The isoelectric point (pI) of fungal laccase is around 4.0, whereas the isoelectric point of plant laccases is usually higher than this value (usually above 5.0) [28,29]. These observations suggest that plant and fungal laccases may have functional differences. In plants, laccases are reported to play a role in forming protective barriers. For example, the secretion of laccases in lacquer tree helps wound healing to prevent infection [30,31]. Plant laccases were also reported to regulate plant organ development [32]. In Arabidopsis, SKU5 encoding a laccase was localized to both the plasma membrane and the cell wall. Loss of SKU5 function exhibited a left-handed rotation bias in roots and etiolated hypocotyls [33]. Plant laccases were suggested to contribute toward pathogen defense [34,35]. The most well-known function of plant laccases is probably the involvement of the enzyme in the lignification process [36,37]. For example, an increased lignin level was observed in transgenic poplar plants following over-expression of the cotton laccase gene GaLAC1 [38].
Although the presence of bacterial laccase is widely recognized [4,39], reports of characterized and purified bacterial laccase are still very limited. There are relatively few studies on the function of bacterial laccases. Unlike plant and fungal laccases, which are often extracellular, bacterial laccases are intracellular [5]. A well-studied bacterial laccase is CotA, identified in Bacillus subtilis and suggested to participate in the biosynthesis of the brown spore pigment to prevent ultraviolet (UV) damage [40]. Laccases were also identified in some marine bacteria. For example, Lac15 identified from the South China Sea provides chloride tolerance and dye discoloration ability to the bacteria and may have specific industrial applications [41]. In general, laccase is a widely distributed enzyme with some functional differentiation across the kingdoms, providing diversity to laccases for industrial usages. In addition, laccase was always a topic of study ever since it was discovered because of its wide range of potential functions for the fields of environmental protection, biological detection, paper making, and food industries. Thus, identification of novel laccases is always meaningful to facilitate efficient industrial production.
In plants, although laccases were identified and classified in many species such as Arabidopsis thaliana [3,20], sorghum [42], rice [43], Populus [28], and cotton [44], an investigation of the laccase gene family in soybean is yet to be conducted. In this study, we identified all possible laccase-coding genes from the soybean reference genome (Wm82.a2.v1). We then predicted the gene and protein structures, studied the expression patterns of the gene family with or without Phytophthora sojae infection, and investigated the selection pressure exerted on the gene family during the domestication process. This study also generates a strong base for future soybean laccase research.

2. Materials and Methods

2.1. Identification of Laccase Gene Family Members in Soybean

Soybean genome was downloaded from Phytozome (https://phytozome.jgi.doe.gov). The genome of Arabidopsis was obtained from The Arabidopsis Information Resource (TAIR) database (http://arabidopsis.org). In total, 17 laccase sequences were obtained according to published research [3]. Arabidopsis laccase members contain Cu-oxidase, Cu-oxidase_2, and Cu-oxidase_3 (PF00394, PF07731, and PF07732) domains. The three domains were searched in the soybean genome using HAMMER 3.0. We took the intersection of the three BLAST results and manually eliminated genes with no open reading frame and multi-transcripts from the same gene. Additionally, identified genes were verified using the SMART domain database (http://smart.embl-heidelberg.de).

2.2. Physical Map of the Soybean Laccase Genes and Properties of the Laccase Proteins

Gene identifiers (IDs) and chromosomal locations were obtained from the Phytozome database (https://phytozome.jgi.doe.gov/). Mapchart 2.32 was used to visualize laccase gene (GmLac) distribution on soybean chromosomes [45]. Peptide size (amino acids, aa), isoelectric point (pI), and protein molecular weight (Da) were calculated using the online ExPASy tool (https://web.expasy.org/protparam/) [46]. Subcellular locations of the GmLac members were determined using the online software CELLO (http://cello.life.nctu.edu.tw/) [47]. Transmembrane topology and signal peptides were predicted using the online Phobius software (http://phobius.sbc.su.se/index.html) [48].

2.3. Sequence Alignment and Phylogenetic Analysis

An alignment of 93 GmLac and 17 AtLac peptides was performed using ClustalW in MEGA 7.0 [49]. All conserved sites were used to construct evolutionary relationships with the neighbor-joining (NJ) method (1000 replicates). The final phylogenetic tree was visualized and edited in iTOL (http://itol.embl.de/) [50].

2.4. Gene Duplication and Expansion Analysis

Intra-species synteny blocks were detected by the MCScanX algorithm [51], with following settings (five genes required to call a collinear block based on the previous all-to-all BLASTP result E-value ≤ 1 × 10−5). The ratio of nonsynonymous substitutions per nonsynonymous site (Ka) to synonymous substitutions per synonymous site (Ks) was computed using KaKs_Calculator 2.0 [52]. Then, the estimated date (mya, million years ago) of each duplication event was calculated using the mean Ks values (T = Ks/2λ), assuming clocklike rates (λ) of 6.1 × 10−9 [53]. Tandem repeat genes were calculated through the physical location on the chromosome, and segmental duplication genes were defined as those located in the same synteny blocks.

2.5. Gene Structure and Motif Compositions Analysis

A phylogenetic tree based only 93 GmLac peptides was constructed according to same method in Section 2.3. The exon/intron organizations of soybean laccase genes were identified with the Gene Structure Display Server, GSDS 3.0 (http://gsds.cbi.pku.edu.cn/) [54]. Conserved motifs of laccase proteins were identified statistically using MEME (http://meme-suite.org/tools/meme). Discovery mode was set as discriminative mode and we used 17 Arabidopsis laccase protein sequences as a control. The maximum number of motifs to find was set at 10 [55]. Visualization of motif compositions was executed with TBtools V0.52 (https://github.com/CJ-Chen/TBtools).

2.6. Expression Pattern of Soybean Laccase Family Genes

Nine soybean-tissue/organ RNA sequencing (RNA-Seq) data, FPKM (Fragments per kilobase of transcript per million fragments mapped) values, were obtained from the Phytozome database [56]. Average expression valued were calculated, and a three-fold increase compared to the average value was defined as a high-expression gene. Color gradations (light to dark green representing different expression abundances) were appended to the expression data.

2.7. Expression Patterns of GmLacs Following Phytophthora sojae Infection

RNA-Seq data of Williams 82 roots 0, 0.5, 3, 6, and 12 h post P. sojae inoculation were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA318321). The variation amount of expression values (ΔFPKM) ware counted through [n_hpi (FPKM + 0.001)/0_hpi (FPKM + 0.001] and then normalized to log2 (ΔFPKM). Visualization of the expression data was conducted by the pheatmap package in R.

2.8. qRT-PCR Analysis of GmLacs Expression Following Phytophthora sojae Infection

P. sojae strain (P7076) was posted on the roots of 72-h-old seedings (grown in darkness), collecting infected samples at 0, 0.5, 3, 6, and 12 h. Total RNA of these samples was extracted using the Easy-Pure Plant RNA kit (Transgen, ER301). The RNA quality was detected through an ultra-microspectrophotometer (IMPLEN, P330). RNA samples were reverse-transcribed with EasyScript One-Step gDNA Removal and complementary DNA (cDNA) Synthesis SuperMix (Transgen, AE311) strictly following the manufacturer’s instructions. The gene expression pattern was analyzed through qRT-PCR using SYBR® Green Master (Roche, 4913914001) in the Agilent Technologies Stratagene Mx3005P Detection Device. Three replicates were investigated to confirm the result. The relative expression levels were calculated via the 2−ΔΔCt method after taking GmPDF (serine/threonine-protein phosphatase 2A regulatory subunit A gene, Glyma.20G114000) as an internal control. The qRT-PCR analysis primers are listed in the Table S8.

2.9. Selective Pressure Analysis

SNP (Single nucleotide polymorphism) data gathered from resequencing of 302 soybean genomes were used for the genetic diversity analysis of GmLac genes in soybean [57]. The SNPs with missing data >10% or MAF (Minor Allele Frequency) <5% were filtered. The soybean accessions were divided into two populations: Glycine soja and G. max (landrace and improved cultivar). Nucleotide diversity (π) was calculated using a 20-k-10-k sliding window and VCFtools [58]. The ratio of diversity πG.max/πG.soja for each window was calculated. The πG.max/πG.soja ratio of a target region was compared to the mean πG.max/πG.soja ratio of the corresponding entire chromosome to determine whether the target region was under selection pressure or not [59]. The 5% windows with the lowest πG.max/πG.soja ratios values were compared with the target region πG.max/πG.soja ratios to determine the genes that were highly selected [60].

3. Results

3.1. Identification and Classification of Laccase Genes in Soybean

To identify the laccase genes in soybean, we downloaded the reference genome of Glycine max from Phytozome and searched the genome with HAMMER 3.0 software for Cu-oxidase, Cu-oxidase_2, and Cu-oxidase_3 domains (PF00394, PF07731, and PF07732). After elimination of the redundant genes with only one or two typical domains on the peptide sequence or with no integral open reading frames, 93 putative soybean laccase-coding genes GmLac1 through GmLac93 were identified (Table S1). Furthermore, we obtained the physical locations of these laccase genes and mapped them onto 19 of the 20 soybean chromosomes (Figure 1). Chromosome 18 had the greatest number of laccase genes, which was 12 (GmLac70–GmLac81). On the contrary, chromosomes 3, 9, 15, and 16 contained only one gene each. GmLac91 and GmLac92 mapped onto scaffold_27, while GmLac93 mapped onto scaffold_614. Basic information of all soybean laccases including gene name, gene locus, amino-acid length, signal peptide length, molecular weight, pI value, and subcellular localization was also determined and is presented in the Table S1. The amino-acid length ranged from 315 to 637 aa, with molecular weight ranging from 34,958.38 to 71,749.27 Daltons, and isoelectric points (pI) ranging from 4.95 to 9.96. Most soybean laccase pI values were above 7.00, whereas only 19 laccase pI values were below 7.00. The predicted subcellular locations of 76 soybean laccases included the extracellular space. Thirty-four of the laccases were predicted to localize only to the extracellular space. Thus, most laccases are probably secretory proteins. Laccases may also commonly be located in the plasma membrane and less commonly in lysosomes. Finally, very few laccases localized to other organelles including the mitochondria, peroxisomes, and vacuole. Specifically, GmLac21, GmLac39, and GmLac86 may be located in the mitochondria and peroxisomes, whereas GmLac77 may be located only in the mitochondria. GmLac90 was the only one predicted to be located in the vacuole. Transmembrane topology and signal peptide analysis indicated that nine laccases contain transmembrane domains. GmLac29 contains two transmembrane regions. The hydrophobic region (H-Region) of signal peptides of most laccase members ranged from 7–20 aa (Table S1).

3.2. Phylogenetic Analysis of the Soybean Laccase Gene Family

To further investigate the phylogenetic relationships of laccases between soybean and Arabidopsis, we constructed a neighbor-joining phylogenetic tree based on multiple sequence alignments of 93 putative soybean laccase and 17 Arabidopsis laccase peptide sequences (Figure 2). According to classification standard of Arabidopsis laccases and location on the phylogenetic tree [61], 93 soybean laccases were divided into five groups. The distribution of the laccase genes in each group was rather uneven; 15 members were classified as class I, 13 members as class Ⅱ, eight members as class Ⅲ, eight members as class Ⅳ, and 49 members as class Ⅴ. Unlike in Arabidopsis, soybean laccase members in class Ⅳ classified into two subgroups: IVa subgroup (Gm18 and Gm19) and IVb subgroup (the six other genes). None of the soybean laccases clustered into class Ⅵ. The specific gene list of each category is presented in the Table S2a–c.

3.3. Duplication and Expansion Events in Soybean Laccase Gene Family

Gene duplication is considered one of the major driving forces of genome evolution [62]. Duplicated genes are the source for creating novel genetic variation. It is generally accepted that the soybean genome went through three rounds of whole-genome duplication (WGD) events and retained more than 50% of the repeated segments [63,64]. They were the γ WGD event about ~130 to 240 million years ago (Mya), the legume WGD event about ~58 Mya, and the Glycine genus WGD event about ~13 Mya [65,66,67]. The general Ks rate range of soybean genome γ WGD events is bigger than 1.5, while that of the legume WGD event is 0.3 to 1.5, and that of the Glycine genus WGD event is 0 to 0.3 [40,44,47]. For an in-depth study the WGD events of laccase genes, we conducted collinear analysis using the base MCScanX software. In total, 85 orthologous gene pairs were identified within the pairwise syntenic blocks (Figure 3). We calculated the non-synonymous nucleotide substitution rate (Ks) and synonymous nucleotide substitution rate (Ka) between collinear gene pairs (Table S3). We calculated that about 25.9% of gene pairs (22 of 85) were separated during the Glycine genus WGD event, nearly 43.5% of gene pairs (37 of 85) were separated during the legume WGD event, and approximately 30.6% of gene pairs (26 of 85) were separated during the γ WGD event. Additionally, we also calculated the Ka/Ks ratios to study the gene divergence after the WGD events. We observed that all 85 synteny pairs underwent purifying selection with Ka/Ks < 1 [48].
The most common patterns of gene duplication include interspersed repeats, tandem repeats, and segmental duplication [68], while segmental and tandem duplications are usually considered as the main causes of large gene family expansion in plants [69]. To gain insight into the soybean laccase gene duplication pattern, we analyzed the tandem and segmental duplication events of laccase genes. To identify tandem duplication, we searched adjacent gene clusters on the same chromosome with no more than one interferential gene being inserted. We found that 43.0% (40 of 93) of genes in the soybean laccase gene family were tandem repeats (Figure 1 and Table S4).
Segmental duplications are blocks of genomic sequence ranging from 1–200 kb that map to different loci in a genome and share high sequence identity (often >90%) [70,71,72]. In this study, we defined segmental duplication as genes located on the same homologous blocks [73]. Thus, collinearity gene pairs can be regarded as segmental duplication members in the laccase gen family. In total 45.1% (42 of 93) of laccase genes involved segmental duplication (Table S4).

3.4. Gene Structure and Motif Compositions of GmLacs

To reveal the structural diversity of soybean laccase genes, we constructed the exon/intron organizations and searched for conservative motifs based on the phylogenetic tree of all soybean laccase alignments (Figure 4). Overall, nearly all the closest genes on the phylogenetic tree showed remarkably similar gene structures. However, there was still a small proportion of the gene clades that exhibited different intron/exon organizations. For instance, GmLac13 contained six exons, whereas its nearby paralogous genes, GmLac68, GmLac67, GmLac40, and GmLac5, had eight exons even though their evolutionary relationships reached 98–100% bootstrap values, respectively.
To further reveal the conserved motifs of the GmLac proteins, we analyzed 93 GmLac proteins using the MEME program and identified 10 motifs (Figure 4; Figure S1). As expected, the motif compositions of peer groups had similar structures and organizations. Some classes of GmLac proteins with completely the same motif composition may indicate the possibility of functional redundancies among those members. Some motifs were found in all 93 putative members. These conserved motifs may play critical roles in mediating GmLac function. Despite an overall similarity among the 93 GmLac proteins, there were some differences in the number of motifs. For example, the number of motifs varied from 2–10 among those GmLac proteins. Varying numbers of motifs across the GmLac proteins may suggest functional divergence among some of the GmLac proteins.

3.5. Expression Patterns of GmLac Genes

Gene expression patterns often reveal how, when, and where genes are active. To verify the expression profiles of GmLac genes, we analyzed publicly available RNA-Seq data in Phytozome. Various organs and tissues including root, root hair, nodules, pod, seed, stem, shoot apical meristem (SAM), leaves, and flower transcriptomes were downloaded (Table S5a,b). According to the fragments per kilobase million (FPKM) values of 90 GmLacs excluding the three members located in scaffolds, a differential transcriptional regulation pattern for the laccase genes was observed. GmLac25 has the highest expression abundance and is specifically expressed in flowers (Table S5a). The overall mean FPKM value of all 90 GmLac genes among the nine transcriptomes was about 10. We defined the genes with mean values three-fold higher than the overall mean FPKM value across all transcriptomes as highly expressed GmLac genes. With this criterion, we identified 11 genes that were highly expressed in multiple tissues (Table 1). Most of the highly expressed genes showed transcript abundance in the stem, indicating that they may be involved in stem cell-wall maintenance to provide mechanical strength for the aboveground portion of the plant. Conversely, the expression of some genes was restricted to a specific organ or tissue and very poorly expressed in others. For example, GmLac1 and GmLac24 specifically and highly expressed in the roots, GmLac77 did so in pods, GmLac2, GmLac8, and GmLac9 did so in the stem, GmLac87 did so in the SAM, and GmLac25 did so in the flower (Table 2). The promoters of these genes were active in specific organs or tissues. In addition, 11 genes were inactive in all tissues and organs studied (Table S5b). For the remaining members, there was either a mild expression pattern with moderate abundance or expression in several samples.
We also compared the expression pattern of collinear gene pairs (Table S6). In total, 64.7% (55 of 85) of collinear gene pair expression patterns were distinct. Only a fraction of gene pairs (23 of 85) maintained a similar expression pattern. These results additionally suggested a differentiation of the promoter elements of laccase genes and possibly their functions.

3.6. Expression Changes of GmLacs after Inoculation with Phytophthora sojae

It is well established that laccases are involved in lignin metabolism and plant defenses [34,35,36,74,75]. Phytophthora sojae is an oomycete pathogen that causes root and stem rot disease in soybean, which severely affects yield [76]. We were interested if any laccase genes were rapidly induced following P. sojae infection. An earlier transcriptomic dataset of soybean–P. sojae interaction is available on the NCBI database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA318321). We studied the responses of 93 soybean laccase genes following infection (Table S7a). Considering 0 h post infection (hpi) as a control, we calculated the expression change value of the GmLac genes at 0.5, 3, 6, and 12 hpi (Table S7b), and the results are depicted in a heatmap (Figure 5). According to the clustering, 93 putative laccase genes can be roughly divided into three categories based their expression status. Most laccase genes (72 of 93) did not show obvious expression changes following P. sojae infection. Nine of the remaining were induced, including Gmlac28, Gmlac91, Gmlac92, Gmlac79, Gmlac23, Gmlac57, Gmlac47, Gmlac37, and Gmlac36. Although they were induced, their expression levels peaked at different stages. Expression abundances of the remaining 12 genes were suppressed. GmLac32 and GmLac75 were strongly repressed at 12 hpi, while the other nine (GmLac78, GmLac71, GmLac25, GmLac89, GmLac15, GmLac90, GmLac45, GmLac48, GmLac76, and GmLac18) were potentially suppressed at multiple stages. By examining the expression profile of these 21 genes (Table S5a), we noticed that these genes were barely expressed in normal conditions. On the other hand, after infection, the expression of these genes changed significantly, which strongly implies a close tie between these genes and oomycete pathogen resistance.
To further test the expression pattern of these genes, we performed qRT-PCR on a Williams 82 root sample with a soybean Phytophthora sojae strain (P7076) [77]. The results (Figure 5) suggested that most (8/9) induced laccase genes were similarly induced with the exception of GmLac36. In particular, GmLac28, GmLac91, GmLac92, GmLac79, GmLac23, and GmLac57 expression abundances were significantly increased in both tests. Conversely, the quantitative results of repressed genes also displayed coincident RNA-seq data. GmLac78, GmLac75, GmLac32, GmLac32, GmLac45, GmLac89, and GmLac25 showed obvious repression in expression. GmLac48 showed no expression signal. Automatically, we suggest a putative defense function for the P. sojae-induced genes, while suppressed genes could equally be important, as P. sojae effectors may manage to suppress the expression of some defense-related GmLac genes to promote colonization.

3.7. Selection Pressure Analysis of Laccase Gene Family during Soybean Domestication

Cultivated soybean was domesticated from its relative Glycine soja in ancient China [78], which led to a series of “domestication syndromes” and subsequently resulted in markedly decreased nucleotide diversity (π) of the respective domestication genes [79]. Wild soybeans uniformly exhibit a softer and thinner stem, while cultivated soybeans uniformly show thicker, stronger, and erect stems, suggesting that these phenotypes are most likely domestication-related traits [80].
Lignin is abundant in the stem cortex and xylem cell wall, playing an important role in maintaining mechanical strength of stem [81]. Laccase is well known for its ability to facilitate lignin synthesis and promote plant tissue-specific lignification [82]. Therefore, we wanted to investigate whether the domestication of soybean was accompanied by changes in laccase gene nucleotide diversity. Thus, we performed selection pressure analysis of GmLac genes based on π ratio (πG.max/πG.soja) between 62 wild species and 240 (130 landraces and 110 improved cultivars) cultivated soybean accessions [57]. A comparison of the πG.max/πG.soja values of a chromosomal region containing GmLac genes with the corresponding mean value of the respective chromosome (Table S9a) revealed that 57 of 90 (63.3%) GmLac genes exhibited above average πG.max/πG.soja values. This observation indicates that the nucleic acid diversity of GmLac genes significantly decreased following the domestication of soybean, and also suggests that most laccase genes were under selective pressure during the domestication of cultivated soybean (Table S9b). When judging whether a gene underwent selection pressure by comparing with the top 5% of πG.max/πG.soja values according to previous research [57,60], we identified 10 genes located in the selective sweep regions (Figure 6; Table S9c). Of these 10 genes, GmLac8, GmLac12, GmLac66, and GmLac85 are highly expressed in the stem, and, more importantly, the expression of GmLac8 is highly specific to the stem (Table S9d). These genes are most likely involved in the transition from a soft and thin soybean stem to a thick and erect stem.
To further support our conclusion that 10 genes have underwent strong selection pressure during the domestication of soybean, we conducted phylogenetic analysis using SNPs of the ten genes among 302 resequenced soybean accessions (Supplementary Materials). The 302 accessions can be divided into five groups (Figure 7). In Groups Ⅰ through Ⅳ, all accessions except SRR1533282 are wild soybean species. In Group Ⅴ, all but SRR1533169 and SRR1533198 are soybean cultivars. These results also imply that the nucleotide diversity of the 10 genes is much higher in the wild accessions than in the cultivated accessions.

4. Discussion

Laccase enzymes are multifunctional proteins. They were shown to contribute toward cell morphology [32,33], pigment formation [83], secondary cell-wall biosynthesis [84,85], and resistance to biotic and abiotic stresses in plant [74,75]. Functional studies of laccases were already carried out in many model and crop plants, including Arabidopsis [86], Brachypodium distachyon [87], sorghum [42], Oryza sativa [43], Populus trichocarpa [88], Gossypium raimondii, and Gossypium arboretum [74]. However, as one of the most important crops around the world, the functions of soybean laccase are rarely reported. In this study, gene structures, protein properties, phylogenetic relationships, replication characteristics, expression profiles, and selection pressures were investigated for members of the soybean laccase gene family.
The soybean Williams 82 genome encompasses 105 putative copper oxidase domain-containing genes that can be regarded as putative candidate laccase genes. Typical laccases usually have three copper oxidase domains; therefore, we investigated 93 GmLac genes that contained three copper oxidase domains, representing 0.2% of the total number of soybean genes, which is higher than the laccase gene proportions in other plant species including Arabidopsis thaliana (0.06%) [89], Brachypodium distachyon (0.13%) [90], Sorghum bicolor (0.08%) [91], Oryza sativa (0.08%) [92], Populus trichocarpa (0.07%) [93], Gossypium raimondii (0.11%) [94], and Gossypium arboretum (0.11%) [95]. Predicted isoelectric points of most soybean laccases (74 of 93) are higher than 7, suggesting that they tend to function in alkaline environments. Predicted subcellular locations of the GmLac proteins suggest that the vast majority of laccases (76 of 93) are extracellular. This result implies that most soybean laccases are secretory proteins.
Through multiple sequence alignment of 17 Arabidopsis and 93 soybean laccases and construction of an N-J phylogenetic tree, we divided the soybean laccases into five classes according to the classification criteria applied in Arabidopsis thaliana. Class Ⅳ was divided into two different branches, Ⅳa and Ⅳb. Further genome-wide duplication (GWD), homologous, and historical duplication event analyses of the GmLac genes revealed that laccase genes are an ancient gene family that expanded as early as ~321 Mya ago during the γ WGD event that was accompanied by the evolution of rosids. Laccase genes continued to expand in subsequent legume WGD and Glycine WGD events. According to the Ks value, purifying selection was the major force driving the evolution of the laccase gene family. Synteny analysis and the physical location of GmLac genes show that tandem repeats (41.93%) and segmental duplications (43.01%) were the main mechanisms of laccase gene duplication.
Gene structure and protein motif analyses showed that the most closely related members in the phylogenetic tree had similar exon/intron organization and common motif compositions. This may suggest the existence of a possible functional similarities among the closely related GmLac proteins. However, within the five basic classes, exon/intron organization and motif compositions were highly diverse, suggesting the occurrence of functional differentiation in each subgroup of the soybean laccase gene family.
Investigation of the expression profiles of GmLac genes among nine organs/tissues revealed the possible role of transcriptional regulation in mediating the functions of GmLac proteins. Some members have high expression abundance, while some genes are barely expressed. Both organ/tissue-specific and constitutive types of expressions were observed among members of the gene family.
In response to P. sojae infection, only a few GmLac genes were either induced or repressed. Although it is possible to assign a putative defense function for the P. sojae-induced genes, downregulation of some genes could equally be important in Phytophthora resistance. Exchange of the promoter of such a gene was shown to enhance resistance in transgenic soybean plants against this pathogen [96]. P. sojae effectors may manage to suppress the expression of some defense-related GmLac genes to avoid lignified cell barriers. Laccase genes are able to enhance plant resistance through accelerating the polymerization of monolignol to lignin. On one hand, increased lignin content helps create stronger physical barriers to prevent invasion; on the other hand, lignin synthesis reinforces cell biosynthesis, which is a branch of the phenylpropyl pathway, and phenylpropanoid metabolism is one of the most important secondary metabolic pathways involved in plant defense against biotic and abiotic stresses [75].
Our study on the natural variations in the GmLac gene family across a collection of both wild and cultivated accessions revealed that the laccase gene family may have been subjected to selection pressure following domestication (Figure 6). Together with the differential expression patterns and selection pressure, four genes were found to be highly expressed in the stem. Based on the differences in stem characteristics between wild soybean and cultivated soybean and laccase’s role in lignification, we speculate that these four genes may have made an important contribution to the transition from a soft and thin stem in wild species to a hard, thick, and erect stem in present day soybean cultivars.

5. Conclusions

A total of 93 putative GmLac genes were identified from the William 82 genome. Based on the phylogenetic relationships between these genes and Arabidopsis thaliana laccase genes, they were divided into five classes (Figure 1 and Figure 2; Tables S1 and S2). Gene duplication and expansion analysis showed that, as expected, duplication of the members of the soybean laccase gene family occurred in all three WGD events in soybean (Figure 3; Table S3). Purifying selection was the major driving force in laccase gene family evolution. The main gene duplication mechanisms for GmLac gene family were tandem and segmental duplications. The analyses of conserved domains suggested that, in general, adjacent members in the phylogenetic tree had common motif compositions. However, exceptions were observed within a clade for motif compositions. Differences were recorded for functional conservation, as well as in the motif compositions within clades (Figure 4; Table S4). Transcriptomic analysis revealed both tissue- and organ-specific specific expression patterns, as well as constitutive expression patterns, of the soybean laccase gene family (Table 1 and Table 2; Table S5). This study also revealed that distinct GmLac gene sets are induced or repressed in response to invasions by oomycete pathogen P. sojae (Figure 5; Table S7). Selective pressure analysis revealed a possible role of a few members of the gene family in transitioning from the soft and thin stem of wild soybean species to the strong, thick, and erect stem of the cultivated soybean accessions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/10/9/701/s1, Table S1: Basic information of soybean laccases; Table S2a: Gene numbers in each phylogenetic class in soybean and Arabidopsis thaliana, Table S2b: Specific gene list of each phylogenetic class of Arabidopsis thaliana, Table S2c: Specific gene list of each phylogenetic class of soybean; Table S3: Collinear gene pairs and their Ka, Ks value and estimated divergence time; Table S4: List of tandem repeat and segmental repeat laccase genes; Table S5a: Expression pattern of soybean laccase gene family; Table S5b: Genes that are barely expressed in 9 detecteted materials; Table S6: Comparison of expression patterns between collinearity pairs; Table S7a: Soybean laccase gene expression pattern following P. sojae infection; Table S7b: Change value of laccase gene inoculation with P. sojae (△FPKM); Table S8: Primers used for qRT-PCR; Table S9a: Sliding window analysis of nucleotide diversity of 93 laccase gene sequence; Table S9b: List of selection GmLacs above average value; Table S9c: List of GmLacs experience sellection pressure above top 5% genome-wide level; Table S9d: Expression pattern of the 10 highly selected GmLacs; Figure S1: Schematic diagram of the motifs of soybean laccase proteins; Data S1: SNPs (filtered) of all laccase gene of 302 soybean accessions.

Author Contributions

X.F. and S.Y. designed the research; Q.W. and G.L. performed bioinformatics analyses; Q.W. performed experiments; K.Z., X.Z., J.M., D.W., K.T., X.F., J.L. and H.Y. helped with data analysis; Q.W. and X.F. wrote the manuscript.

Funding

This research was funded by the Key Research Program of the Chinese Academy of Sciences (Grant Nos. ZDRW-ZS-2019-2-02) and by programs (Grant nos. 31571692 and 31471743) from the National Natural Science Foundation of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yaropolov, A.I.; Skorobogat, O.V.; Vartanov, S.S.; Varfolomeyev, S.D. Laccase. Appl. Biochem. Biotechnol. 1994, 49, 257–280. [Google Scholar] [CrossRef]
  2. Alfred, M.M.; Richard, C.S. Laccase: New functions for an old enzyme. Phytochemistry 2002, 60, 551–565. [Google Scholar]
  3. 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] [PubMed] [Green Version]
  4. Alexandre, G.; Zhulin, I.B. Laccases are widespread in bacteria. Trends Biotechnol. 2000, 18, 41–42. [Google Scholar] [CrossRef]
  5. Diamantidis, G.; Effosse, A.; Potier, P.; René, B. Purification and characterization of the first bacterial laccase in the rhizospheric bacterium Azospirillum lipoferum. Soil Biol. Biochem. 2000, 32, 919–927. [Google Scholar] [CrossRef]
  6. Claus, H. Laccases and their occurrence in prokaryotes. Arch. Microbiol. 2003, 179, 145–150. [Google Scholar] [CrossRef] [PubMed]
  7. Ausec, L.; Zakrzewski, M.; Goesmann, A.; Schlueter, A.; Mandic-Mulec, I. Bioinformatic analysis reveals high diversity of bacterial genes for laccase-like enzymes. PLoS ONE 2011, 6, e25724. [Google Scholar] [CrossRef]
  8. Yoshida, H. Chemistry of lacquer (urushi). J. Chem. Soc. 1883, 43, 472–486. [Google Scholar] [CrossRef]
  9. Bertrand, G. Sur la laccase et sur le pouvoir oxydant de cette diastase. CR Acad. Sci. 1895, 120, 266–269. [Google Scholar]
  10. Kudanga, T.; Nyanhongo, G.S.; Guebitz, G.; Burton, S. Potential applications of laccase-mediated coupling and grafting reactions: A review. Enzyme Microb. Technol. 2011, 48, 195–208. [Google Scholar] [CrossRef]
  11. Bento, I.; Carrondo, M.A.; Lindley, P.F. Reduction of dioxygen by enzymes containing copper. J. Biol. Inorg. Chem. 2006, 11, 539–547. [Google Scholar] [CrossRef] [PubMed]
  12. Shiba, T.; Xiao, L.; Miyakoshi, T.; Chen, C.L. Oxidation of isoeugenol and coniferyl alcohol catalyzed by laccases isolated from rhus vernicifera stokes and pycnoporus coccineus. J. Mol. Catal. B Enzym. 2000, 10, 605–615. [Google Scholar] [CrossRef]
  13. Rodriguez, C.S.; Toca-Herrera, J.L. Industrial and biotechnological applications of laccases: A review. Biotechnol. Adv. 2006, 24, 500–513. [Google Scholar] [CrossRef] [PubMed]
  14. Chirivì, C.; Fontana, G.; Monti, D.; Ottolina, G.; Riva, S.; Danieli, B. The quest for new mild and selective modifications of natural structures: Laccase-catalysed oxidation of ergot alkaloids leads to unexpected stereoselective C-4 hydroxylation. Chem. Eur. J. 2012, 18, 10355–10361. [Google Scholar] [CrossRef]
  15. Claus, H. Laccases: Structure, reaction, distribution. Micron 2004, 35, 93–96. [Google Scholar] [CrossRef] [PubMed]
  16. Morozova, O.; Shumakovich, G.P.; Gorbacheva, M.; Yaropolov, A.I. “Blue” laccases. Biochemistry 2007, 72, 1136–1150. [Google Scholar] [CrossRef] [PubMed]
  17. Linares, N.; Fernández, F.; Loske, A.; Gómez-Lim, M.A. Enhanced delignification of lignocellulosic biomass by recombinant fungus phanerochaete chrysosporium overexpressing laccases and peroxidases. J. Mol. Microb. Biotechnol. 2018, 28, 1–13. [Google Scholar]
  18. Leonowicz, A.; Cho, N.S.; Luterek, J.; Jarosz-Wilkołazka, A.; Wojtas-Wasilewska, M.; Matuszewska, A.; Hofrichter, M.; Dobritzsch, D.; Rogalski, J. Fungal laccase: Properties and activity on lignin. J. Basic Microb. 2001, 41, 185–227. [Google Scholar] [CrossRef]
  19. Baldrian, P. Fungal laccases—Occurrence and properties. FEMS Microbiol. Rev. 2006, 30, 215–242. [Google Scholar] [CrossRef]
  20. Hoegger, P.J.; Kilaru, S.; James, T.Y.; Thacker, J.R.; Kües, U. Phylogenetic comparison and classification of laccase and related multicopper oxidase protein sequences. FEBS J. 2006, 273, 2308–2326. [Google Scholar] [CrossRef]
  21. Wösten, H.A.B.; Wessels, J.G.H. The Emergence of Fruiting Bodies in Basidiomycetes. In Growth, Differentiation and Sexuality; Ursula, K., Reinhard, F., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 393–414. [Google Scholar]
  22. Lettera, V.; Piscitelli, A.; Leo, G.; Birolo, L.; Pezzella, C.; Sannia, G. Identification of a new member of pleurotus ostreatus laccase family from mature fruiting body. Fungal Biol. UK 2010, 114, 724–730. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, J.; Chen, H.; Chen, M.; Ren, A.; Huang, J.; Wang, H.; Zhao, M.; Feng, Z. Cloning and functional analysis of a laccase gene during fruiting body formation in Hypsizygus marmoreus. Microbiol. Res. 2015, 179, 54–63. [Google Scholar] [CrossRef] [PubMed]
  24. Liu, N.; Cao, Z.; Cao, K.; Ma, S.; Gong, X.; Jia, H.; Dai, D.; Dong, J. Identification of laccase-like multicopper oxidases from the pathogenic fungus Setosphaeria turcica and their expression pattern during growth and infection. Eur. J. Plant Pathol. 2019, 153, 1149–1163. [Google Scholar] [CrossRef]
  25. Cloete, M.; Mostert, L.; Fischer, M.; Halleen, F. Pathogenicity of south african hymenochaetales taxa isolated from esca-infected grapevines. Phytopathol. Mediterr. 2015, 54, 368–379. [Google Scholar]
  26. Temp, U.; Eggert, C. Novel interaction between laccase and cellobiose dehydrogenase during pigment synthesis in the white rot fungus Pycnoporus cinnabarinus. Appl. Environ. Microb. 1999, 65, 389–395. [Google Scholar]
  27. Díaz, L.; Del Río, J.A.; Pérez-Gilabert, M.; Ortuño, A. Involvement of an extracellular fungus laccase in the flavonoid metabolism in Citrus fruits inoculated with Alternaria alternata. Plant Physiol. Biochem. 2015, 89, 11–17. [Google Scholar] [CrossRef]
  28. Berthet, S.; Thévenin, J.; Baratiny, D.; Demont-Caulet, N.; Debeaujon, I.; Bidzinski, P.; Leplé, J.-C.; Huis, R.; Hawkins, S.; Gomez, L.-D.; et al. Role of plant laccases in lignin polymerization. In Lignins—Biosynthesis, Biodegradation and Bioengineering; Lise, J., Catherine, L., Eds.; Elsevier Academic Press: Amsterdam, The Netherlands, 2012; Chapter 5; pp. 145–172. [Google Scholar]
  29. Rodgers, C.J.; Blanford, C.F.; Giddens, S.R.; Skamnioti, P.; Armstrong, F.A.; Gurr, S.J. Designer laccases: A vogue for high-potential fungal enzymes? Trends Biotechnol. 2010, 28, 63–72. [Google Scholar] [CrossRef]
  30. Meehan, T. The lacquer tree. Proc. Acad. Natl. Sci. Phila. 2014, 29, 279. [Google Scholar]
  31. Nitta, K.; Kataoka, K.; Sakurai, T. Primary structure of a Japanese lacquer tree laccase as a prototype enzyme of multicopper oxidases. J. Inorg. Biochem. 2002, 91, 125–131. [Google Scholar] [CrossRef]
  32. Berthet, S.; Demontcaulet, N.; Pollet, B.; Bidzinski, P.; Cézard, L.; Bris, P.L.; Borrega, N.; Hervé, J.; Blondet, E.; Balzergue, S. Disruption of Laccase4 and 17 results in tissue-specific alterations to lignification of arabidopsis thaliana stems. Plant Cell 2011, 23, 1124–1137. [Google Scholar] [CrossRef]
  33. Sedbrook, J.C.; Carroll, K.L.; Kai, F.H.; Masson, P.H.; Somerville, C.R. The arabidopsis SKU5 gene encodes an extracellular glycosyl phosphatidylinositol-anchored glycoprotein involved in directional root growth. Plant Cell 2002, 14, 1635–1648. [Google Scholar] [CrossRef] [PubMed]
  34. Walker, J.R.L.; Ferrar, P.H. Diphenol oxidases, enzyme-catalysed browning and plant disease resistance. Biotechnol. Genet. Eng. 1998, 15, 457–498. [Google Scholar] [CrossRef]
  35. Iqbal, M.J.; Ahsan, R.; Afzal, A.J.; Jamai, A.; Meksem, K.; Elshemy, H.A.; Lightfoot, D.A.; Elshemy, H.A. Multigeneic QTL: The laccase encoded within the soybean Rfs2/rhg1 locus inferred to underlie part of the dual resistance to cyst nematode and sudden death syndrome. Curr. Issues Mol. Biol. 2009, 11 (Suppl. S1), 11. [Google Scholar]
  36. Bao, W.; O’Malley, D.M.; Whetten, R.; Sederoff, R.R. A laccase associated with lignification in loblolly pine xylem. Science 1993, 260, 672–674. [Google Scholar] [CrossRef] [PubMed]
  37. Sterjiades, R.; Dean, J.F.; Eriksson, K.E. Laccase from sycamore maple (acer pseudoplatanus) polymerizes monolignols. Plant Physiol. 1992, 99, 1162–1168. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, J.; Wang, C.; Zhu, M.; Yu, Y.; Zhang, Y.; Wei, Z. Generation and characterization of transgenic poplar plant overexpressing a cotton laccase gene. Plant Cell Tiss. Org. 2008, 93, 303–310. [Google Scholar] [CrossRef]
  39. Sharma, P.; Goel, R.; Capalash, N. Bacterial Laccases. World J. Microb. Biot. 2007, 23, 823–832. [Google Scholar] [CrossRef]
  40. Hullo, M.F.; Moszer, I.; Danchin, A.; Martin, V.I. CotA of bacillus subtilis Is a copper-dependent laccase. J. Bacteriol. 2001, 183, 5426–5430. [Google Scholar] [CrossRef]
  41. Fang, Z.; Li, T.; Wang, Q.; Zhang, X.; Peng, H.; Fang, W.; Hong, Y.; Ge, H.; Xiao, Y. A bacterial laccase from marine microbial metagenome exhibiting chloride tolerance and dye decolorization ability. Appl. Microbiol. Biot. 2011, 89, 1103–1110. [Google Scholar] [CrossRef]
  42. Wang, J.; Feng, J.; Jia, W.; Fan, P.; Bao, H.; Li, S.; Li, Y. Genome-wide identification of sorghum bicolor laccases reveals potential targets for lignin modification. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef]
  43. Liu, Q.; Luo, L.; Wang, X.; Shen, Z.; Zheng, L. Comprehensive analysis of rice laccase gene (OsLAC) family and ectopic expression of OsLAC10 enhances tolerance to copper stress in arabidopsis. Int. J. Mol. Sci. 2017, 18, 209. [Google Scholar] [CrossRef] [PubMed]
  44. Balasubramanian, V.K.; Rai, K.M.; Thu, S.W.; Hii, M.M.; Mendu, V. Genome-wide identification of multifunctional laccase gene family in cotton (Gossypium spp.); Expression and biochemical analysis during fiber development. Sci. Rep. UK 2016, 6. [Google Scholar] [CrossRef] [PubMed]
  45. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed]
  46. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.E.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein Identification and Analysis Tools on the ExPASy Server; Humana Press: Totowa, NJ, USA, 2005. [Google Scholar]
  47. Yu, C.S.; Chen, Y.C.; Lu, C.H.; Hwang, J.K. Prediction of protein subcellular localization. Proteins 2006, 64, 643–651. [Google Scholar] [CrossRef] [PubMed]
  48. Käll, L.; Krogh, A.; Sonnhammer, E.L. Advantages of combined transmembrane topology and signal peptide prediction—The Phobius web server. Nucleic Acids Res. 2007, 35, W429–W432. [Google Scholar] [CrossRef]
  49. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, msw054. [Google Scholar] [CrossRef] [PubMed]
  50. Ivica, L.; Peer, B. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar]
  51. Wang, Y.; Li, J.; Paterson, A.H. MCScanX-transposed: Detecting transposed gene duplications based on multiple colinearity scans. Bioinformatics 2013, 29, 1458–1460. [Google Scholar] [CrossRef]
  52. Wang, D.; Zhang, Y.; Zhang, Z.; Zhu, J.; Yu, J. KaKs_Calculator 2.0: A toolkit incorporating γ-series methods and sliding window strategies. Genom. Proteom. Bioinf. 2010, 8, 77–80. [Google Scholar] [CrossRef]
  53. Lynch, M.; Conery, J.S. The evolutionary fate and consequences of duplicate genes. Science 2000, 290, 1151–1155. [Google Scholar] [CrossRef]
  54. Guo, A.; Zhu, Q.H.; Chen, X.; Luo, J.C. GSDS: A gene structure display server. Yi chuan = Hereditas/Zhongguo yi chuan xue hui bian ji 2007, 29, 1023–1026. [Google Scholar] [CrossRef]
  55. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef] [PubMed]
  56. Libault, M.; Farmer, A.; Joshi, T.; Takahashi, K.; Langley, R.J.; Franklin, L.D.; He, J.; Xu, D.; May, G.; Stacey, G. An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants. Plant J. 2010, 63, 86–99. [Google Scholar] [CrossRef] [PubMed]
  57. Zhou, Z.; Jiang, Y.; Wang, Z.; Gou, Z.; Lyu, J.; Li, W.; Yu, Y.; Shu, L.; Zhao, Y.; Ma, Y.; et al. Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improvement in soybean. Nat. Biotechnol. 2015, 33, 408. [Google Scholar] [CrossRef] [PubMed]
  58. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C.A.; Banks, E.; DePristo, M.A.; Handsaker, R.E.; Lunter, G.; Marth, G.T.; Sherry, S.T.; et al. The variant call format and VCFtools. Bioinformatics 2011, 27, 2156–2158. [Google Scholar] [CrossRef] [PubMed]
  59. Zhang, D.; Sun, L.; Li, S.; Wang, W.; Ding, Y.; Swarm, S.A.; Li, L.; Wang, X.; Tang, X.; Zhang, Z.; et al. Elevation of soybean seed oil content through selection for seed coat shininess. Nat. Plants 2018, 4, 30–35. [Google Scholar] [CrossRef] [PubMed]
  60. Wang, M.; Li, W.; Fang, C.; Xu, F.; Liu, Y.; Wang, Z.; Yang, R.; Zhang, M.; Liu, S.; Lu, S.; et al. Parallel selection on a dormancy gene during domestication of crops from multiple families. Nat. Genet. 2018, 50, 1435–1441. [Google Scholar] [CrossRef] [PubMed]
  61. McCaig, B.C.; Meagher, R.B.; Dean, J.F. Gene structure and molecular analysis of the laccase-like multicopper oxidase (LMCO) gene family in Arabidopsis thaliana. Planta 2005, 221, 619–636. [Google Scholar] [CrossRef] [PubMed]
  62. Moore, R.C.; Purugganan, M.D. The early stages of gene duplication. Proc. Natl. Acad. Sci. USA 2004, 100, 15682–15687. [Google Scholar] [CrossRef] [PubMed]
  63. Lavin, M.; Herendeen, P.; Wojciechowski, M. Evolutionary rates analysis of leguminosae implicates a rapid diversification of lineages during the tertiary. Syst. Biol. 2005, 54, 575–594. [Google Scholar] [CrossRef] [PubMed]
  64. Pfeil, B.E.; Schlueter, J.A.; Shoemaker, R.C.; Doyle, J.J. Placing paleopolyploidy in relation to taxon divergence: A phylogenetic analysis in legumes using 39 gene families. Syst. Biol. 2005, 54, 441–454. [Google Scholar] [CrossRef] [PubMed]
  65. Fan, C.; Wang, X.; Hu, R.; Wang, Y.; Xiao, C.; Jiang, Y.; Zhang, X.; Zheng, C.; Fu, Y.F. The pattern of phosphate transporter 1 genes evolutionary divergence in Glycine max L. BMC Plant Biol. 2013, 13, 48. [Google Scholar] [CrossRef] [PubMed]
  66. Schlueter, J.A.; Lin, J.Y.; Schlueter, S.D.; Vasylenko-Sanders, I.F.; Deshpande, S.; Yi, J.; O’Bleness, M.; Roe, B.A.; Nelson, R.L.; Scheffler, B.E.; et al. Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing. BMC Genom. 2007, 8, 330. [Google Scholar] [CrossRef] [PubMed]
  67. Dai, A.H.; Yang, S.X.; Zhou, H.K.; Tang, K.Q.; Li, G.; Leng, J.T.; Yu, H.; Zhang, Y.H.; Gao, J.S.; Yang, X.; et al. Evolution and expression divergence of the CYP78A subfamily genes in soybean. Genes 2018, 9, 611. [Google Scholar] [CrossRef] [PubMed]
  68. Kong, H.; Landherr, L.; Frohlich, M.; Leebens-Mack, J.; Ma, H.; dePamphilis, C.W. Patterns of gene duplication in the plant SKP1 gene family in angiosperms: Evidence for multiple mechanisms of rapid gene birth. Plant J. 2007, 50, 873–885. [Google Scholar] [CrossRef] [PubMed]
  69. 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. [Google Scholar] [CrossRef] [PubMed]
  70. Van de Peer, Y.; Maere, S.; Meyer, A. The evolutionary significance of ancient genome duplications. Nat. Rev. Genet. 2009, 10, 725–732. [Google Scholar] [CrossRef] [Green Version]
  71. Giannuzzi, G.; D’Addabbo, P.; Gasparro, M.; Martinelli, M.; Carelli, F.; Antonacci, D.; Ventura, M. Analysis of high-identity segmental duplications in the grapevine genome. BMC Genom. 2011, 12. [Google Scholar] [CrossRef]
  72. Lupski, J.R.; Stankiewicz, M.D. Genomic Disorders-The Genomic Basis of Disease; Humana Press: New York, NY, USA, 2006; pp. 73–88. [Google Scholar]
  73. Schmutz, J.; Cannon, S.B.; Schlueter, J.A.; Ma, J.X.; Mitros, T.; Nelson, W.; Hyten, D.L.; Song, Q.; Thelen, J.J.; Cheng, J.; et al. Genome sequence of the palaeopolyploid soybean. Nature 2010, 465, 120. [Google Scholar] [CrossRef]
  74. Hu, Q.; Min, L.; Yang, X.; Jin, S.; Zhang, L.; Li, Y.; Ma, Y.; Qi, X.; Li, D.; Liu, H.Y.; 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]
  75. Zhang, S.; Ju, C.L.; Wang, X. Arabidopsis laccase Gene AtLAC4 regulates plant growth and responses to abiotic stress. CBB 2013, 47, 357–365. [Google Scholar]
  76. Gijzen, M.; Qutob, D. Phytophthora Sojae and Soybean. In Oomycete Genetics and Genomics: Diversity, Interactions, and Research Tools; Kurt, L.S., Ed.; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 303–329. [Google Scholar]
  77. Forster, H. Phytophthora sojae Races have arisen by clonal evolution and by rare outcrosses. Mol. Plant Microbe Interact. 1994, 7, 780. [Google Scholar] [CrossRef]
  78. Dashiell, K. SOYBEANS: Improvement, Production, and Uses, 3rd ed.; American Society of Agronomy, Crop Science Society of America, Soil Science Society of America: Madison, WI, USA, 2005; pp. 110–111. [Google Scholar]
  79. Zohary, D.; Hopf, M. Domestication of plants in the old world. J. Appl. Ecol. 1993, 130, 365–372. [Google Scholar]
  80. Liu, B.; Fujita, T.; Yan, Z.; Sakamoto, S.; Xu, D.J. QTL mapping of domestication-related traits in soybean (Glycine max). Ann. Bot. Lond. 2007, 100, 1027–1038. [Google Scholar] [CrossRef] [PubMed]
  81. Luo, Z.; Li, P.; Cai, D.; Chen, Q.; Qin, P.; Tan, T.; Cao, H. Comparison of performances of corn fiber plastic composites made from different parts of corn stalk. Ind. Crop Prod. 2016, 95, 521–527. [Google Scholar] [CrossRef]
  82. Dean, J.F.; Eriksson, K.E. Laccase and the deposition of lignin in vascular plants. Holzforschung 1994, 48, 21–33. [Google Scholar] [CrossRef]
  83. Pourcel, L.; Routaboul, J.M.; Kerhoas, L.; Caboche, M.; Lepiniec, L.; Debeaujon, I. Transparent testa10 encodes a laccase-like enzyme involved in oxidative polymerization of flavonoids in arabidopsis seed coat. Plant Cell 2005, 17, 2966–2980. [Google Scholar] [CrossRef] [PubMed]
  84. Ranocha, P.; Chabannes, M.; Chamayou, S.; Danoun, S.; Jauneau, A.; Boudet, A.M.; Goffner, D. Laccase down-regulation causes alterations in phenolic metabolism and cell wall structure in poplar. Plant Physiol. 2002, 129, 145–155. [Google Scholar] [CrossRef] [PubMed]
  85. Zhao, Q.; Nakashima, J.; Chen, F.; Yin, Y.; Fu, C.; Yun, J.; Shao, H.; Wang, X.; Wang, Z.Y.; Dixon, R. Laccase is necessary and nonredundant with peioxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell 2013, 25, 3976–3987. [Google Scholar] [CrossRef]
  86. Wang, C.Y.; Zhang, S.; Yu, Y.; Luo, Y.C.; Liu, Q.; Ju, C.; Zhang, Y.C.; Qu, L.H.; Lucas, W.J.; Wang, X.; et al. MiR397b regulates both lignin content and seed number in Arabidopsis via modulating a laccase involved in lignin biosynthesis. Plant Biotechnol. J. 2014, 12, 1132–1142. [Google Scholar] [CrossRef] [Green Version]
  87. Wang, Y.; Bouchabké-Coussa, O.; Lenris, P.; Antelme, S.; Soulhat, C.; Gineau, E.; Dalmais, M.; Bendahmane, A.; Morin, H.; Mouille, G.; et al. Lassase 5 is required for lignification of the brachypodium distachyon culm. Plant Physiol. 2015. [Google Scholar] [CrossRef] [PubMed]
  88. Lu, S.; Li, Q.; Wei, H.; Chang, M.J.; Tunlaya-Anukit, S.; Kim, H.; Liu, J.; Song, J.; Sun, Y.H.; Yuan, L.; et al. Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Proc. Natl. Acad. Sci. USA 2013, 110, 10848–10853. [Google Scholar] [CrossRef] [PubMed]
  89. The Arabidopsis Genome, I. Analysis of the genome sequence of the flowering plant arabidopsis thaliana. Nature 2000, 408, 796–815. [Google Scholar] [CrossRef] [PubMed]
  90. The International Brachypodium Initiative; Vogel, J.P.; Garvin, D.F.; Mockler, T.C.; Schmutz, J.; Rokhsar, D.; Bevan, M.W.; Barry, K.; Lucas, S.; Harmon-Smith, M.; et al. Genome sequencing and analysis of the model grass brachypodium distachyon. Nature 2010, 463, 763–768. [Google Scholar] [CrossRef] [PubMed]
  91. Paterson, A.H.; Bowers, J.E.; Bruggmann, R.; Dubchak, I.; Grimwood, J.; Gundlach, H.; Haberer, G.; Hellsten, U.; Mitros, T.; Poliakov, A.; et al. The sorghum bicolor genome and the diversification of grasses. Nature 2009, 457, 551–556. [Google Scholar] [CrossRef] [PubMed]
  92. Kawahara, Y.; de la Bastide, M.; Hamilton, J.P.; Kanamori, H.; McCombie, W.R.; Ouyang, S.; Schwartz, D.C.; Tanaka, T.; Wu, J.; Zhou, S.; et al. Improvement of the oryza sativa nipponbare reference genome using next generation sequence and optical map data. Rice 2013, 6, 4. [Google Scholar] [CrossRef] [PubMed]
  93. Tuskan, G.A.; DiFazio, S.; Jansson, S.; Bohlmann, J.; Grigoriev, I.; Hellsten, U.; Putnam, N.; Ralph, S.; Rombauts, S.; Salamov, A.; et al. The genome of black cottonwood, populus trichocarpa (Torr. & Gray). Science 2006, 313, 1596–1604. [Google Scholar]
  94. Wang, K.; Wang, Z.; Li, F.; Ye, W.; Wang, J.; Song, G.; Yue, Z.; Cong, L.; Shang, H.; Zhu, S.; et al. The draft genome of a diploid cotton Gossypium raimondii. Nat. Genet. 2012, 44, 1098–1104. [Google Scholar] [CrossRef] [PubMed]
  95. Li, F.; Fan, G.; Wang, K.; Sun, F.; Yuan, Y.; Song, G.; Li, Q.; Ma, Z.; Lu, C.H.; Zou, C.; et al. Genome sequence of the cultivated cotton Gossypium arboreum. Nat. Genet. 2014, 46, 567–574. [Google Scholar] [CrossRef]
  96. Ngaki, M.N.; Wang, B.; Sahu, B.; Srivastava, S.; Farooqi, M.; Sekhar, K.; Swaminathan, S.; Bhattacharyya, M.K. Tanscriptomic study of the soybean-fusarium virguliforme interaction revealed a novel ankyrin-repeat containing defense gene, expression of whose during infection led to enhanced resistance to the fungal pathogen in transgenic soybean plants. PLoS ONE 2016, 11, e0163106. [Google Scholar] [CrossRef]
Genes 10 00701 g001 550
Figure 1. Physical map of the soybean laccase genes. Chromosome size is indicated by its relative length. The scale bar is shown on the left, and the numbers represent the physical position of the genes on the chromosomes. Tandemly duplicated genes are represented by boxes with red outlines.
Figure 1. Physical map of the soybean laccase genes. Chromosome size is indicated by its relative length. The scale bar is shown on the left, and the numbers represent the physical position of the genes on the chromosomes. Tandemly duplicated genes are represented by boxes with red outlines.
Genes 10 00701 g001
Genes 10 00701 g002 550
Figure 2. The phylogenetic tree of the Arabidopsis and soybean laccase proteins. In total, 93 soybean and 17 Arabidopsis laccase proteins were aligned using the ClustalW program. The neighbor-joining tree was constructed using MEGA7.0. The bootstrap value was set as 1000 replicates. The length of each branch is shown next to the branches.
Figure 2. The phylogenetic tree of the Arabidopsis and soybean laccase proteins. In total, 93 soybean and 17 Arabidopsis laccase proteins were aligned using the ClustalW program. The neighbor-joining tree was constructed using MEGA7.0. The bootstrap value was set as 1000 replicates. The length of each branch is shown next to the branches.
Genes 10 00701 g002
Genes 10 00701 g003 550
Figure 3. Collinearity of the laccase gene (GmLac) gene pairs. The outermost scale represents the chromosomes. Genes are listed in the inner circle according to their chromosomal location. Synteny relationships between gene pairs are marked with a black line.
Figure 3. Collinearity of the laccase gene (GmLac) gene pairs. The outermost scale represents the chromosomes. Genes are listed in the inner circle according to their chromosomal location. Synteny relationships between gene pairs are marked with a black line.
Genes 10 00701 g003
Genes 10 00701 g004 550
Figure 4. Phylogenetic relationships, gene structure, and motif compositions of the soybean laccase gene family. (A) A neighbor-joining tree of 93 GmLac protein sequences constructed using MEGA v7.0. (B) The structure of the 93 GmLac genes. Yellow squares correspond to exons and linking black lines indicate introns, while the blue squares refer to the 5′ untranslated region (UTR) and 3′ UTR sequences. (C) Schematic motif composition of 93 GmLac proteins. The colored boxes represent the different motifs, indicated in the top right-hand corner. The scales at the bottom of the image indicate the estimated exon/intron length in kb and motif length in numbers of amino acids (aa).
Figure 4. Phylogenetic relationships, gene structure, and motif compositions of the soybean laccase gene family. (A) A neighbor-joining tree of 93 GmLac protein sequences constructed using MEGA v7.0. (B) The structure of the 93 GmLac genes. Yellow squares correspond to exons and linking black lines indicate introns, while the blue squares refer to the 5′ untranslated region (UTR) and 3′ UTR sequences. (C) Schematic motif composition of 93 GmLac proteins. The colored boxes represent the different motifs, indicated in the top right-hand corner. The scales at the bottom of the image indicate the estimated exon/intron length in kb and motif length in numbers of amino acids (aa).
Genes 10 00701 g004
Genes 10 00701 g005 550
Figure 5. Expression profiles of the 93 soybean laccase genes following Phytophthora sojae infection. (A) Transcriptome heatmap. The hierarchical cluster color code: the highest induction values vs. controls are displayed in red (hot), while the values displayed in blue (cool) indicate the most reduced expression levels vs. controls; intermediate expression levels are represented by lighter intermediate colors. The color scale represents △FPKM normalized to log2-transformed vales. The genes marked by the black brackets on the right in (A) are those obviously mis-regulated genes. (B) Expression levels of 21 GmLacs after P7076 infection. Bars represent average values of three replicates ± standard deviation (SD). All expression levels of the GmLac genes were normalized to the expression levels of GmPDF. ***, **, and * indicate significant differences compared to the control (0 h) at p < 0.001, p < 0.01, and p < 0.05, respectively.
Figure 5. Expression profiles of the 93 soybean laccase genes following Phytophthora sojae infection. (A) Transcriptome heatmap. The hierarchical cluster color code: the highest induction values vs. controls are displayed in red (hot), while the values displayed in blue (cool) indicate the most reduced expression levels vs. controls; intermediate expression levels are represented by lighter intermediate colors. The color scale represents △FPKM normalized to log2-transformed vales. The genes marked by the black brackets on the right in (A) are those obviously mis-regulated genes. (B) Expression levels of 21 GmLacs after P7076 infection. Bars represent average values of three replicates ± standard deviation (SD). All expression levels of the GmLac genes were normalized to the expression levels of GmPDF. ***, **, and * indicate significant differences compared to the control (0 h) at p < 0.001, p < 0.01, and p < 0.05, respectively.
Genes 10 00701 g005
Genes 10 00701 g006 550
Figure 6. Selective pressure analysis of laccase gene family. The outermost scale represents chromosomes, and the green bar chart in the middle represents πG.max/πG.soja values per chromosome. The top 5% selected laccase genes are listed according to their physical location.
Figure 6. Selective pressure analysis of laccase gene family. The outermost scale represents chromosomes, and the green bar chart in the middle represents πG.max/πG.soja values per chromosome. The top 5% selected laccase genes are listed according to their physical location.
Genes 10 00701 g006
Genes 10 00701 g007 550
Figure 7. Phylogenetic tree of 302 accessions. SNPs of the 10 selected genes among resequenced genomes of the 302 soybean accessions were used to construct the neighbor-joining (N-J) phylogenetic tree. The bootstrap value was set as 1000 replicates. Corresponding specie idnetifiers (IDs) of the accessions are available from Reference [57].
Figure 7. Phylogenetic tree of 302 accessions. SNPs of the 10 selected genes among resequenced genomes of the 302 soybean accessions were used to construct the neighbor-joining (N-J) phylogenetic tree. The bootstrap value was set as 1000 replicates. Corresponding specie idnetifiers (IDs) of the accessions are available from Reference [57].
Genes 10 00701 g007
Table
Table 1. Constitutively expressed laccase genes (GmLac).
Table 1. Constitutively expressed laccase genes (GmLac).
Gene NameGene IDRootRHNodulesPodSeedStemSAMLeavesFlowerScale Bar
GmLac5Glyma.01G18310029.8766.0433.9937.2314.5546.4124.142.006.780.00
GmLac10Glyma.04G01950011.5944.8911.6173.40111.7553.3823.712.2611.455.00
GmLac12Glyma.05G05610054.7365.7760.5223.3622.1647.6355.086.075.5210.00
GmLac16Glyma.06G30710036.5049.4710.29129.804.8199.05114.4014.2818.4315.00
GmLac40Glyma.11G05920033.5359.2840.5133.9210.2687.0843.946.109.0120.00
GmLac55Glyma.13G076900162.5050.7530.2312.383.0936.758.7712.194.6325.00
GmLac62Glyma.14G22300027.3253.2632.2221.1215.6833.048.421.002.1030.00
GmLac65Glyma.17G01230069.51127.5854.7326.1133.7326.3162.911.394.16
GmLac66Glyma.17G13830066.8356.0548.3656.5235.8884.1343.278.409.79
GmLac69Glyma.17G26150047.0965.8829.722.730.3628.5316.043.854.14
GmLac85Glyma.20G05190093.6567.1530.2832.4744.2148.5113.4912.042.55
Note: ID, identifier; RH, root hairs; SAM, stem apical meristem. Data in the table are FPKM (fragments per kilobase million) values obtained from Phytozome. Table colors are labeled according to scale bar.
Table
Table 2. Organ-/tissue-specific expression of GmLac genes.
Table 2. Organ-/tissue-specific expression of GmLac genes.
Gene NameGene IDRootRHNodulesPodSeedStemSAMLeavesFlowerScale Bar
GmLac1Glyma.01G10820041.644.946.820.060.000.040.010.030.140.00
GmLac2Glyma.01G1126000.946.3012.5810.450.4570.533.315.931.925.00
GmLac8Glyma.02G2616002.954.363.947.970.6932.580.560.981.9110.00
GmLac9Glyma.03G0779001.007.2815.9310.780.6175.273.255.622.1815.00
GmLac24Glyma.07G14260038.5311.6621.130.080.000.030.092.860.2120.00
GmLac25Glyma.07G2253000.200.000.000.010.001.060.200.03505.0325.00
GmLac77Glyma.18G1837000.161.652.1335.660.0516.900.301.160.3230.00
GmLac87Glyma.20G17260013.600.030.000.000.012.87117.910.002.54
Note: ID, identifier; RH, root hairs; SAM, stem apical meristem. Data in the table are FPKM (fragments per kilobase million) value obtained from Phytozome. Table colors are labeled according to scale bar.

Share and Cite

MDPI and ACS Style

Wang, Q.; Li, G.; Zheng, K.; Zhu, X.; Ma, J.; Wang, D.; Tang, K.; Feng, X.; Leng, J.; Yu, H.; et al. The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection. Genes 2019, 10, 701. https://doi.org/10.3390/genes10090701

AMA Style

Wang Q, Li G, Zheng K, Zhu X, Ma J, Wang D, Tang K, Feng X, Leng J, Yu H, et al. The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection. Genes. 2019; 10(9):701. https://doi.org/10.3390/genes10090701

Chicago/Turabian Style

Wang, Quan, Guang Li, Kaijie Zheng, Xiaobin Zhu, Jingjing Ma, Dongmei Wang, Kuanqiang Tang, Xingxing Feng, Jiantian Leng, Hui Yu, and et al. 2019. "The Soybean Laccase Gene Family: Evolution and Possible Roles in Plant Defense and Stem Strength Selection" Genes 10, no. 9: 701. https://doi.org/10.3390/genes10090701

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