- 1Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, China
- 2Function Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
- 3CNRS, UMR 8227, Integrative Biology of Marine Models, Station Biologique de Roscoff, Roscoff, France
- 4Sorbonne Universités, UPMC Univ Paris 06, UMR 8227, Integrative Biology of Marine Models, Station Biologique de Roscoff, Roscoff, France
Lipoxygenase (LOX) plays important roles in fatty acid oxidation and lipid mediator biosynthesis. In this study, we give first insights into brown algal LOX evolution. Whole genome searches revealed four, three, and eleven LOXs in Ectocarpus siliculosus, Cladosiphon okamuranus, and Saccharina japonica, respectively. In phylogenetic analyses, LOXs from brown algae form a robust clade with those from prokaryotes, suggesting an ancestral origin and slow evolution. Brown algal LOXs were divided into two clades, C1 and C2 in a phylogenetic tree. Compared to the two species of Ectocarpales, LOX gene expansion occurred in the kelp S. japonica through tandem duplication and segmental duplication. Selection pressure analysis showed that LOX genes in brown algae have undergone strong purifying selection, while the selective constraint in the C2 clade was more relaxed than that in the C1 clade. Furthermore, within each clade, LOXs of S. japonica evolved under more relaxed selection constraints than E. siliculosus and C. okamuranus. Structural modeling showed that unlike LOXs of plants and animals, which contain a β barrel in the N-terminal part of the protein, LOXs in brown algae fold into a single domain. Analysis of previously published transcriptomic data showed that LOXs in E. siliculosus are responsive to hyposaline, hypersaline, oxidative, and copper stresses. Moreover, clear divergence of expression patterns was observed among different life stages, as well as between duplicate gene pairs. In E. siliculosus, all four LOXs are male-biased in immature gametophytes, and mature gametophytes showed significantly higher LOX mRNA levels than immature gametophytes and sporophytes. In S. japonica, however, our RNA-Seq data showed that most LOXs are highly expressed in sporophytes. Even the most recently duplicated gene pairs showed divergent expression patterns, suggesting that functional divergence has likely occurred since LOX genes duplicated, which potentially contributes to the production of various oxylipins in brown algae.
Introduction
Lipid oxidation is an essential biological process in all living organisms and is subject to both developmental and environmental regulation (Liavonchanka and Feussner, 2006). Lipoxygenases (LOX; EC 1.13.11.12) constitute an important group of enzymes responsible for this kind of lipid metabolism. They are non-haem iron-containing dioxygenases that catalyze the addition of oxygen to polyunsaturated fatty acids (PUFAs) to produce hydroperoxides, which are further converted into a series of biologically active compounds collectively named oxylipins (Fatima et al., 2017). The oxidation of PUFAs via LOXs and the subsequent reactions are collectively called LOX-pathway (Blée, 2002). LOX enzymes are encoded by a multigenic family. They are found in a wide variety of plants, fungi, animals, as well as algae, protozoa, and bacteria (Kalms et al., 2017). The abundant linoleic acid and linolenic acid are major substrates of LOXs in plants, while the principal substrate in animals is arachidonic acid (Porta and Rocha-Sosa, 2001). Derived oxylipins, such as jasmonic acid in plants or prostaglandins and leukotrienes in animals, are important bio-regulators and participate in a wide range of physiological functions, including the regulation of growth, development, and senescence, the mediation of stress responses (Bueno et al., 2001; Hou et al., 2015), or the maintenance of cell homeostasis (Molina and Donaire, 2002; Ackermann et al., 2017).
Lipoxygenases contain a highly conserved region rich in His residues. Five His residues are arranged in the form of His-(X)4-His-(X)4-His-(X)17-His-(X)8-His, which is important for iron atom binding and catalytic activity (Minor et al., 1996). The molecular mass of the LOX proteins identified to date ranges from 75 to 80 kDa in animals, 94–105 kDa in plants and fungi, and 49–75 kDa in bacteria (Feng et al., 2010; Garreta et al., 2013). In both animals and plants, the enzymes possess a small N-terminal domain PLAT/LH2 (homologous to polycystin-1, lipoxygenase, α-toxin, or lipoxygenase) and a major C-terminal domain. With the rapid progress in functional genomics, LOXs have been identified and analyzed in many land plants. At least six LOX genes in Arabidopsis thaliana have been identified and expression of AtLOX1 was stimulated by the hormones abscisic acid and methyl jasmonate (MeJA) (Melan et al., 1993). Chen et al. identified 20 LOXs in the poplar genome and showed that many of them have various expression patterns across different tissues and stress treatments (Chen et al., 2015). Chen et al. also investigated the molecular evolution of LOX in modern rosid plants and found that strong purifying selection plays a critical role in LOX evolution (Chen et al., 2016). Lipid peroxidation caused by LOX is of particular importance in fruit ripening and forming volatile flavor aroma, which is a critical aspect of fruit quality. Twelve LOX genes in peach were differentially expressed during fruit ripening under different post-harvest treatments (Guo et al., 2017). In melon and grape berry, members of LOXs also exhibit complex patterns of gene expression in different berry tissues, or during different stages of ripening (Podolyan et al., 2010; Zhang et al., 2014). In pear and apple, an expansion of the LOX family was observed, accompanying a whole genome duplication event that occurred in the common ancestor of the two species. Expression patterns of 18 out of 23 LOXs in pear fruit samples at six different stages after flowering were found to correspond to the changes in the volatile components (Li et al., 2014). These differential expression patterns indicate rapid functional diversification after gene duplication. LOXs are also found in prokaryotes, though in this kingdom of life they are present at much lower frequency (Kalms et al., 2017). The first crystal structure of a prokaryote LOX was obtained in the opportunistic pathogen Pseudomonas aeruginosa (Garreta et al., 2013). This LOX was suggested to play a role in the interaction with the host cell membrane. Another study later reported this LOX exhibits a phospholipid oxygenase activity and destroys red blood cells (Banthiya et al., 2015).
LOX-pathway derived oxylipins are also known to be produced in many marine macroalgae. Some of them belong to the prostaglandins and share great similarities with the products of LOXs in animals (Imbs et al., 2001; Potin et al., 2002). The contents of hydroxyl-oxylipins have been quantified in 40 macroalgae, including green, brown, and red algae (Kumari et al., 2014). In the red alga Rhodymenia pertusa, four new oxylipins were identified, supporting the existence of a 5R-LOX in this species (Jiang et al., 2000). Furthermore, 8-, 9- and 12-LOX pathways were found in other red algae (Gerwick et al., 1999). Brown algae, together with diatoms and oomycetes, are part of an independent lineage, the stramenopiles that evolved from a secondary endosymbiosis event between a protist and an ancestral unicellular red alga. Oxylipins have been reported in several bloom-forming diatoms (Wichard et al., 2005). Polyunsaturated aldehydes have, for instance, been detected in the marine diatom Thalassiosira rotula, where LOX-mediated fatty acid transformation might function as chemical defense by reducing the hatching success from eggs of copepods (Pohnert, 2002), whereas the synthesis of oxylipins in T. rotula relies on an unusual heterolytic mechanism involving a multifunctional LOX to perform both hydroperoxidation and cleavage of hydroperoxide (Barofsky and Pohnert, 2007). Another example is the pennate diatom Pseudo-nitzschia delicatissima, which produces three different species of oxylipins via the 15-LOX pathway (d’Ippolito et al., 2009). Kousaka et al. (2003) also found a multitude of oxylipin metabolites in the brown alga Eisenia bicyclis.
Oxylipin studies in macroalgae have been mainly focused on their isolation, structural characterization, and biological properties. The past decade has seen a few studies on deciphering the roles of macroalgal oxylipins in defense and immunity. Brown algae contain high amounts of both C18 and C20 fatty acids (Hofmann and Eichenberger, 1997; Küpper et al., 2009). In the brown alga Laminaria digitata, arachidonic acid, linolenic acid, and MeJA were strong triggers of an oxidative burst, and brought about changes in fatty acid and oxylipin profiles (Küpper et al., 2009), which suggests that the lipid oxidation pathway plays a major role as a defense mechanism in brown algae. Oxylipins derived from oxidation of polyunsaturated fatty acids (PUFAs) in brown algae were responsive to external stimuli or mezo-grazers (Pohnert and Boland, 2002). Prostaglandin A2 has been reported to trigger an oxidative burst in Laminaria sp. and may constitute a novel defense inducer in brown algae (Zambounis et al., 2012). Copper stress in the brown alga L. digitata induced the production of as many as 23 oxylipins, and treatment with the inhibitors of LOX enzymes efficiently reduced oxylipin production (Ritter et al., 2008). Ritter et al. (2014) also reported the detection of diverse oxylipins produced via the LOX pathways in the copper-stressed Ectocarpus siliculosus.
Although LOX products have been described from Phaeophyceae, to the best of our knowledge, there has been no comprehensive study of LOX genes in brown algae available up to now, and their phylogenetic relationships with LOXs of other phyla remain poorly understood. Brown macroalgae are the dominant vegetation in many benthic marine habitats with high ecological and economic significance (Xu et al., 2017). Particularly, Saccharina japonica is one of the most economically important seaweed in marine aquaculture. The completion of genome sequencing projects for brown algae Ectocarpus (Cock et al., 2010), S. japonica (Ye et al., 2015) and Cladosiphon okamuranus (Nishitsuji et al., 2016) greatly facilitate exhaustive inventories of gene families and cross phyla comparisons. Furthermore, comparative genomics between closely related species provides a way to understand evolution and function of genes in multigene families and their roles in adaptation to ecological niches (Good et al., 2014). In this study, we characterized the brown algal LOX family and attempted to determine their expansion mechanism in S. japonica. Such knowledge will increase our understanding of the roles of LOXs and oxylipins in brown algal development as well as facilitate the utilization of these algal enzymes as tools in synthetic chemistry for the production of oxylipins.
Materials and Methods
Database Search and LOX Sequence Retrieval
Protein and CDS sequences of E. siliculosus V2016 version were downloaded from the website http://bioinformatics.psb.ugent.be/orcae/overview/Ectsi. Those of S. japonica were downloaded from NCBI under the accession code JXRI00000000. Sequences of C. okamuranus were downloaded from http://marinegenomics.oist.jp/algae/. Other genomes were downloaded from the following sites: green plants, Phytozome V12.0 https://phytozome.jgi.doe.gov/pz/portal.html. Emiliania huxleyi, oomycetes, http://genome.jgi.doe.gov/pages/tree-of-life.jsf. Other genes were downloaded at NCBI https://www.ncbi.nlm.nih.gov/. The LOX domain PF00305 was downloaded from the Pfam website http://pfam.xfam.org/ (Finn et al., 2015). HMMER was used to search for this domain in the predicted proteins of each species. All of the obtained LOXs were further manually examined using the NCBI online BLAST tool https://blast.ncbi.nlm.nih.gov/Blast.cgi to confirm the genes’ annotation. The online InterProScan program http://www.ebi.ac.uk/interpro/search/sequence-search was used to confirm the presence of LOX domain in each sequence.
Phylogenetic Analysis
Sequence alignments were performed on both the full-length proteins and the conserved PF00305 domain using MUSCLE embedded in MEGA 7.0 (Kumar et al., 2016). The maximum likelihood (ML) and Neighbor-joining (NJ) phylogenetic trees were constructed with best substitution models and rate parameter WAG + G + I found by MEGA 7.0. Bootstrap with 1000 bootstrap replicates was performed to obtain the confidence support level. Bayesian inference of phylogeny was carried out using BEAST2 (Bouckaert et al., 2014). Analyses were run for 100 million generations, with sampling every 10000th generation using the same substitution model as above and Yule Model as a prior. The first 10000 samples were discarded as burn-in.
Sequence Structure Analysis
Intron and exon information of the three brown algal LOX genes was obtained from their GFF files and was used for the graphic display using the Gene Structure Display Server of Peking University (Hu et al., 2015). The protein transmembrane helices were predicted by CCTOP1 (Dobson et al., 2015). Structural motifs were searched for using the MEME program http://meme-suite.org/ (Bailey et al., 2009). The motif number was set to 15. The most conserved motif of LOXs was displayed on the WebLogo website http://weblogo.berkeley.edu/logo.cgi (Crooks et al., 2004). Signal peptides were predicted using the website http://www.cbs.dtu.dk/services/SignalP/ (Petersen et al., 2011). Subcellular protein localization was predicted using ASAFind (Gruber et al., 2015). Mass weight (MW) and isoelectric point (pI) was calculated using the website http://www.expasy.org/.
Molecular Evolution Analysis
In total, 18 protein sequences of LOXs in three brown algae Ectocarpus, S. japonica, C. okamuranus were aligned using MUSCLE. A ML tree was constructed with MEGA 7.0 as described above. Brown algae LOXs were divided into two clades C1 and C2. In each clade, the branches were divided again into Ectocarpales clades (C-E) and Laminariales clades (C-S). To evaluate variation in selective pressures between different clades (C1-E vs. C1-S; C2-E vs. C2-S), a branch model in PAML v4.8 (Yang, 2007) was used to estimate the ω (dN/dS) ratio under two assumptions: a one ratio model (M = 0) that assumes the same ω in C1 and C2 clade, respectively; a two-ratio model (M = 2) that assumes that Ectocarpales clade (C-E) and Saccharina clade (C-S) have different ω ratios. To test which of the models best fit the data, likelihood ratio tests (LRTs) were performed by comparing twice the difference in log-likelihood values using a χ2 distribution with one degree of freedom. We also ran site models (M = 0, NSsites = M0, M3, M7, M8) in Codeml to test if there were sites of positive selection in brown algal LOXs. LRTs were performed (M0 vs. M3, M7 vs. M8) to verify which model fit the data using a χ2 distribution with two degrees of freedom. To further explore the evolutionary dynamics of expanded LOXs in S. japonica, pairwise dN, dS and dN/dS values between duplicated pairs of S. japonica LOXs were calculated using the online PAL2NAL program http://www.bork.embl.de/pal2nal/index.cgi (Suyama et al., 2006).
Expression of LOX Genes in E. siliculosus and S. japonica
Previous microarray data of the E. siliculosus transcriptome (Dittami et al., 2009; Ritter et al., 2014) was used to explore the expression levels of LOXs in response to abiotic stressors, notably hyposaline stress, hypersaline stress, oxidative stress, and Cu stress. The expression level was determined by averaging the expression values (previously quantile normalized by Roche NimbleGen, Madison, WI, United States) of four replicates for each experimental condition. Furthermore, expression levels of LOXs were examined in different life cycle stages. In E. siliculosus we compared previously published RNA-Seq data of parthenosporphytes, male gametophytes, and female gametophytes (Lipinska et al., 2015). Expression levels of S. japonica LOXs were examined in sporophytes, male gametophytes, and female gametophytes. Their transcriptomes were sequenced by our group using an Illumina HiSeq 2500, with three biological replicates in each life stage. Gene expression levels were estimated by FPKM (Fragments Per Kilobase of transcript per Million fragments mapped). Differential expression analysis among life stages was performed using the DESeq (Love et al., 2014) R package using a model based on the negative binominal distribution. The resulting p-values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate (Benjamini and Hochberg, 1995). Genes with an adjusted p-value < 0.05 found by DESeq were assigned as differentially expressed. The raw sequence data was deposited at NCBI database. The SRA accession numbers are SRR5860560, SRR5860561, SRR5860562, SRR5860563, SRR5860564, SRR5860565, SRR5860566, SRR5860567, and SRR5860568.
Structure Modeling
To predict the 3D structures of S. japonica LOX proteins, the X-ray structure of a P. aeruginosa LOX protein (PDB accession numbers 4g33.1, 5lc8.1, 4rpe.1) was used as template, as it has the highest coverage and similarity with S. japonica LOXs. Each LOX was aligned with the template and the structure was automatically built using SWISS-MODEL https://www.swissmodel.expasy.org/ (Biasini et al., 2014). Secondary structure was displayed using the Espript web server http://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi (Gouet et al., 2003).
Results
Phylogenetic Analysis and Structures of the LOX Gene Family in Brown Algae
A comprehensive genome-wide search identified 12, 4, and 3 LOX genes in genomes of the S. japonica, E. siliculosus, and C. okamuranus, respectively. The number of LOX genes in S. japonica is approximately three-fold that in E. siliculosus and C. okamuranus. The accuracy of these gene models was manually checked by comparison with publicly available RNA-Seq data. We found a large number of reads spanning the SJ21862 and SJ21863, suggesting they correspond to a single transcript encoding a 518 residue protein. Furthermore, SJ16633 and SJ16634 were truncated at the N-terminal end and manually extended increasing their length by 25 and 91 residues, respectively. Additionally, SJ08291 and SJ05859 have the shortest protein length, and no expression data available, suggesting they are probably pseudogenes. Detailed information about the LOX family in brown algae is listed in Table 1. The length of these proteins ranged from 346 to 747 amino acid residues. The molecular masses ranged from 39 to 82 kDa, with SJ08291 and Ec-03_000010 being the smallest and largest, respectively. The pI values of these proteins are similar, ranging from 4.65 to 5.79. In S. japonica, E. siliculosus, and C. okamuranus, three, one, and two LOXs were predicted to be located in the chloroplast, respectively. In order to perform phylogenetic analysis, we also searched for the LOX genes in other organisms, including oomycetes, Phaeodactylum tricornutum, Thalassiosira pseudonana, Fragilariopsis cylindrus, Nannochloropsis gaditana, and Emiliania huxleyi, all of which, together with brown algae belong to the Stramenopiles and Haptophytes. Besides, LOX in prokaryotes, model plants and animal were included. Sequences corresponding to the LOX domain (PF00305) were aligned using MUSCLE imbedded in MEGA 7.0. Then ML and NJ trees were constructed using MEGA 7.0, and a Bayesian tree was constructed using BEAST 2. The three trees exhibited nearly the same topology. As shown in Figure 1, brown algal LOXs formed one robust clade, with 100% bootstrap support. It is interesting that LOXs of brown algae are most closely related to those of prokaryotes with high bootstrap support but not with other Stramenopiles, plants, or animals. From this tree topology, we propose that LOXs in both brown algae and prokaryotes, as well as LOXs in plant/animal or haptophytes/oomycetes originate from several different ancestral LOXs, and that each lineage retained different copies.
FIGURE 1. Unrooted ML tree of 76 members of the LOX family from prokaryotes and 14 eukaryote species. Support values correspond to bootstrap values (1000 replicates) for ML and NJ trees and Bayesian posterior probabilities in this order; they are shown only for nodes with ≥50% support with all three methods. Nodes with 100% in all methods are shown by a solid dot. The tree was generated using MEGA 7.0 and WAG + G + I model for ML and NJ analyses and BEAST2 for Bayesian methods.
To gain further insight into the diversity of brown algal LOXs, we constructed a separate phylogenetic tree using only the brown algal LOX protein sequences. The tree showed that the 18 LOXs were divided into two distinct clades (C1 and C2, Figure 2). C1 and C2 contained 9 genes each. The C1 clade was further divided into two subclades (C1-S and C1-E). Similarly, the C2 clade also contained two subclades (C2-S and C2-E), in which the ‘S’ and ‘E’ denote Saccharina and Ectocarpales, respectively. Subclades C1-E and C2-E contained LOXs from E. siliculosus and C. okamuranus, whereas subclasses C1-S and C2-S contained only LOXs from S. japonica. Striking gene structure conservation was found within 18 LOXs of brown algae in terms of either intron numbers or exon length. A total of 15 out of 18 genes contain from 15 to 20 introns, of which ten genes contain 17 introns with exons of similar length. Fewest exons were seen in SJ05859 and SJ08291, with only ten exons, which further supported that they are possibly pseudogenes. Intron length in S. japonica varies from 39 bp in SJ16634 to 9,722 bp in SJ12567 whereas that in E. siliculosus varies from 159 bp in Ec-03_000010 to 6,520 bp in Ec-20_003610. Pairwise comparison of the 18 LOX full-length protein sequences revealed some features (Figure 3). C1 clade proteins showed 25.4–85.5% pairwise sequence identity. C2 clade showed 14.1–78.1% pairwise sequence identity. However, the protein sequences between the two LOX clades were significantly more different (Mann–Whitney test, p < 0.01) and the pairwise sequence identity was 8.2–58.2% (Supplementary Table S1).
FIGURE 2. Phylogenetic relationship and gene structure of brown algal LOX genes. The phylogenetic tree was constructed by the ML method with 1000 bootstrap replicates. The two LOX clades were designated as C1 and C2, in which the subclades were designated as C1-S, C1-E, C2-S, C2-E, respectively. Exons and introns are represented by green boxes and black lines, respectively.
FIGURE 3. Pairwise sequence identity of full-length LOX proteins. C1 and C2 represent pairwise sequence identities of clade 1 and clade 2 proteins, respectively. C1–C2 represent pairwise sequence identities between the C1 and C2 clades. The boxplot shows the median (black line), interquartile range (box), and maximum and minimum values (whiskers) for each data set. Circles are outliers. Characters above the figure indicate significant differences between groups at p = 0.01, Mann–Whitney test.
To explore the detailed motif composition of the brown algal LOX family, we predicted putative motifs using the online tool MEME (Bailey et al., 2009). In total, 15 motifs were identified (Figure 4 and Supplementary Table S2). The presence of conserved motifs in LOX proteins provides additional support for the results of the phylogenetic analysis. Compared with the C2 clade, two distinct motifs occurred only in the C1 clade, that is, motifs 14 and 15 in the N-terminal region of the protein. However, whether these motifs confer a specific function to members of the C1 clade remains to be surveyed. Motifs 1 and 6 contain the representative 38 amino acid residues with the structure His-(X)4-His-(X)4-His-(X)17-His-(X)8-His. In the protein sequence alignment of S. japonica, all of the nine sequences have the conserved motif in this His rich region, except two sequences (SJ03316 and SJ16634) where the first His is replaced by Tyr (Figure 5). Motifs 4 and 5 contain an Asn and Ile residue respectively, which are essential for the binding of iron. Notably, some LOXs lack critical motifs. For instance, SJ05859 lacks the six motifs in C-terminal whereas SJ08291 lacks the eight motifs in N-terminal. One gene in C. okamuranus lacks motif 1 and four motifs in C-terminal.
FIGURE 4. Schematic representation of the conserved motifs in brown algal LOX proteins. Each colored box represents a motif in the protein, with the motif number indicated in the middle of the box. The length of the protein and motif can be estimated using the scale at the bottom.
FIGURE 5. A 38-residue motif among S. japonica LOX sequences. (A) The logo was created with 11 S. japonica LOX sequences. The overall height of each stack indicates the sequence conservation at that position and the height of each residue letter indicates the relative frequency of the corresponding amino acid residue. (B) Sequence alignments of 38-residue motifs in S. japonica LOX members.
Molecular Evolution
The ratio of non-synonymous (dN) to synonymous (dS) nucleotide substitutions ω (dN/dS) provides information about the evolutionary forces operating on a gene or gene region (Yang et al., 1998). ω < 1 suggests negative selection and the removal of adverse mutation. ω = 1 and ω > 1 suggest neutral selection and positive selection, respectively. The LOX family of brown algae was divided into two clades (C1 and C2), and they were further divided into subclades (C1-S and C1-E, C2-S and C2-E). To infer the selection pressure within each clade, we estimated ω values with branch models using the CODEML program in the PAML software package (Table 2) (see Materials and Methods). The log-likelihood values under the one ratio model were lnL = -8197.2974 and -7447.3493, with estimates of ω = 0.10530 and 0.1262 for the C1 and C2 clades LOX, respectively. Under the two-ratio model, the log-likelihood values were lnL = -8193.6592 and -7441.3471 for the C1 and C2 clade LOX genes, respectively. Likelihood ratio tests (LRTs) of the two assumptions indicated that the two-ratio model was significantly more likely (p < 0.01) for both the C1 and C2 LOX clades, suggesting that selective pressure varied among LOX subclades (C1-S vs. C1-E, C2-S vs. C2-E). LOX genes in S. japonica have experienced more relaxed purifying selection than those in E. siliculosus and C. okamuranus. All of the estimated ω values were <1, suggesting that the two LOX clades have been under purifying selection but that the selection constraint on the C2 clade was more relaxed than that on the C1 clade, which is consistent with the lower sequence identity in the C2 clade. However, the ratio averaged over all sites is unlikely to exceed 1, because positive selection is not likely to affect all codon sites of a gene (Yang, 2007). Site model allows the ω ratio to vary among sites over all branches (Roux et al., 2014). Therefore, we employed site models (M0 vs. M3, M7 vs. M8) to test whether episodic positive selection had affected specific amino acid sites in each gene. M0 assumes one same ratio among all sites whereas M3 assumes discrete ratios among sites. Under the M7 model, ω follows a beta distribution between 0 and 1, while no ω > 1 site allowed. The M8 model allows an additional site class with ω > 1 (Yang, 2007). LRT was performed to test which model fits the data best. The site model analysis between M0 and M3 indicated that selection pressure varied among sites both for C1 and C2 clades (p < 0.001). However, there is no significant difference between the lnL values of M7 and M8, indicating that no site was positively selected, which further supported the results from branch model.
Gene Duplications in the LOX Gene Family of Brown Algae
In order to illustrate the expansion of LOXs in S. japonica, we compared the LOX number across different organisms. Proteome size-dependent expansion of metal-binding proteins was observed across all domains of life (Dupont et al., 2010). Thus, the relationship between LOX number and proteome size was taken into consideration (Figure 6). LOXs are rare among most sequenced algae. The green alga Volvox carteri has four members, while other green algae, red algae, and diatoms have maximally one LOX. The two brown algae E. siliculosus and C. okamuranus have four and three LOXs, respectively. Notably, S. japonica, which possesses 11 LOXs, has expanded its LOX family compared with other algae. Furthermore, this expansion goes beyond the mere effect of proteome size, as indicated by the different slopes for brown algae and other algae in the linear regressions in Figure 6A. Moreover, a similar effect was not observed for other metal-binding protein families such as the cytochromes P450 (CYP) family (Figure 6B).
FIGURE 6. (A) The total number of LOXs relative to the total proteome size in a selection species. The following species were included: higher plants – Arabidopsis thaliana, Prunus persica, Populus trichocarpa, Brachypodium distachyon, Medicago truncatula, Oryza sativa, Vitis vinifera, Carica papaya, Cucumis sativus, Pyrus bretschneideri, Malus domestica, Fragaria vesca, Cucumis melo; red algae – Chondrus crispus, Galdieria sulphuraria, Cyanidioschyzon merolae, Porphyridium purpureum; diatoms – Thalassiosira pseudonana, Phaeodactylum tricornutum, Fragilariopsis cylindrus; green algae – Chlamydomonas reinhardtii, Chlorella variabilis, Coccomyxa subellipsoidea, Gonium pectoral, Micromonas pusilla, Micromonas sp. RCC29, Ostreococcus lucimarinus, Volvox carteri. Data for higher plants were taken from Feng et al. (2010); Podolyan et al. (2010); Liu et al. (2011); Li et al. (2014); Zhang et al. (2014); Chen et al. (2015, 2016); Guo et al. (2017). (B) Another family of metal-binding proteins, the cytochromes P450 (CYP) family is displayed as a control set.
According to the physical locations of the LOXs on chromosomes or scaffolds, gene duplication events were constructed. In E. siliculosus, the four LOX genes were assigned to two chromosomes. Within each chromosome, tandem duplication occurred once. However, LOX genes in S. japonica have undergone more rounds of duplication and formed a more expanded LOX gene family compared to E. siliculosus. Tandem duplication plays a major role in the S. japonica LOX expansion history. Three tandem clusters account for 7 out of 11 LOX genes, including SJ05858/SJ05859, SJ21859/SJ21862, and SJ16632/SJ16633/SJ16634. Other genes likely came from segmental duplication after a series of tandem duplications. The dN, dS and dN/dS ratio of these putative paralogous gene pairs were calculated to investigate the divergence fate after duplication (Table 3). All dN/dS ratios were <0.5, suggesting strong purifying selection, which is consistent with the branch/site model results. The largest ratio was observed in the gene pair SJ05859/SJ12567 with the value of 0.4393, indicating the relatively relaxed selection pressure between them. Duplicate paralogs diverge from the ancestral state and accumulate synonymous mutations at similar rates over time (Holland et al., 2017), so the dS value can be used to measure the time since gene divergence. The lowest dS (0.0520) was observed between Ec-20_003610/ Ec-20_003620, suggesting that they were created most recently, whereas the highest dS (1.6936) was found between SJ21859/SJ21862, suggesting an ancient duplication, likely derived from the first round of ancestral tandem duplications. By reconciliation of the physical position of genes, gene tree, and the dS values, we attempted to reconstruct the most parsimonious history of S. japonica LOX gene duplication (Figure 7). For the C1 clade, it seems that one tandem duplication created two ancestral genes (b and c), after which other events, likely segmental duplications, occurred to form two gene pairs; one pair is SJ21859 and SJ21862, the other pair is ancestral d and e. After that, two duplication events occurred sequentially, to create two more pairs of LOXs. On the other hand, in the C2 clade, ancestral gene A duplicated to form SJ08291 and ancestral B, followed by two rounds of tandem duplications to create SJ16632 and ancestral C, D, which further produced SJ16633, SJ16634, and SJ03316.
FIGURE 7. Phylogenetic relationships and hypothetical evolutionary histories of the S. japonica LOX genes in Class 1 and 2. (A) Phylogenetic relationships were reconstructed using MEGA7.0. (B) The duplication histories of Class 1. (C) The duplication histories of Class 2. The letters T and S indicate putative tandem duplication and segmental duplication events, respectively. The letters in the tree and blue, green boxes represent ancestral LOX genes.
Expression of LOX Genes in Brown Algae
Duplicated genes may evolve into different fates, which can be inferred by the divergence of their expression patterns. The expression patterns of LOX gene family members were examined using the publically available microarray data of E. siliculosus and RNA-sequencing data of E. siliculosus and S. japonica. Three of the four LOX genes in E. siliculosus are present in the microarray data, and all the three genes have two contigs/singletons. Hierarchical clustering showed that expression levels under Cu stress were clustered with oxidative stress (Figure 8A), which is consistent with the meta-analysis in Ritter et al. (2014). One gene (Ec-20_003610) was significantly influenced by hyposaline, hypersaline, and oxidative stress (log2-fold change > 1 compared to control). The other two genes were up-regulated by hyposaline, hypersaline, and Cu stress, though the log2-fold change did not reach one. This result indicates that LOX genes in E. siliculosus are responsive to these stress conditions, and may explain the observed accumulation of oxylipins under Cu stress (Ritter et al., 2014). Expression divergence was also observed in different life stages. In E. siliculosus, mature male and female gametophytes have significantly higher LOX mRNA levels than sporophytes or immature gametophytes. Notably, all of the four LOXs were male-biased in immature gametophytes (fold change > 2, adjusted p-value < 0.1; Figure 8B), suggesting their potential roles in sexual dimorphism. For the LOX members in S. japonica, more variation in expression patterns was found in different life stages (Figure 8C and Table 3). Three genes were sex-biased, including one female-biased (SJ16633) and two male-biased (SJ16632 and SJ12567). Two genes in the C1 clade (SJ05858, SJ21859) and two genes in the C2 clade (SJ16632, SJ16633) were ranked as the most highly expressed in sporophytes. However, expression levels of other genes were relatively low, especially in gametophytes. Two genes (SJ21862, SJ03316) were exclusively expressed in sporophytes, whereas no expression was detected in gametophytes, suggesting that they are likely only functional in sporophytes. The other genes also exhibited higher expression levels in sporophytes than in gametophytes, except for one gene (SJ12567), which showed the highest expression in male gametophytes. Another gene (SJ08291) was not expressed in any life stage, and SJ05859 had very a low FPKM of only 0.08. These two genes are likely to be pseudogenes, considering their gene structure and expression level. Notably, the two pairs of LOX genes with lowest expression levels (SJ05859 vs. SJ12567, SJ03316 vs. SJ16634) have higher dN/dS values than other LOX gene pairs (Table 3). Likewise, in E. siliculosus, the lowest overall expression level (Ec-20_003610) and great expression difference was observed in the gene pair with the highest dN/dS value (0.3843 compared to 0.1577 in the other pair of genes). The expression patterns of LOX suggested that functional diversification has occurred between the duplicated LOX gene pairs.
FIGURE 8. Expression profiles of LOX genes in E. siliculosus (A,B) and S. japonica (C) respectively. (A) Log2-transformed fold changes of expression levels compared to the control. Black star indicates the log2 values of fold changes (treatments/control) > 1.0 or < -1.0. (B) Log10-transformed RPKM values. Black star indicates the sex-biased expression compared between male and female gametophytes (fold changes > 2, adjusted p-value < 0.05). (C) Log10-transformed FPKM values. Black star indicates sex-biased or sporophyte-biased expression (fold changes > 2, adjusted p-value < 0.05). sp, sporophytes; fg, female gametophytes; mg, male gametophytes.
Structure Modeling
Three-dimensional protein structure modeling across all eleven LOXs in S. japonica revealed that, in general, these proteins are best modeled by a 15-LOX in the bacterium P. aeruginosa (coverage 79–89%, identity 27.83–41.51%), which further supports the close relationships between them. The simulated structures of S. japonica LOX showed remarkable differences to plant or animal LOXs. For most plant and animal LOXs, a two-domain structure has been identified, but LOXs in brown algae fold into a single domain. The helical lid structure in the P. aeruginosa LOX was suggested to be the functional equivalent of the N-terminal β-barrel in eukaryotic LOXs (Garreta et al., 2013). In LOXs of S. japonica, the N-terminal part seems to also contain the helical lid structure (α2A, α2B), though having shorter helices than the P. aeruginosa LOX (Figure 9 and Supplementary Figure S1). The X-ray structure of 4g33.1 showed that five sites, His377, His382, His555, Asn559, and Ile685 are essential for binding of iron (Garreta et al., 2013; Kalms et al., 2017). Sequence comparison showed that the five sites were conserved in all the S. japonica LOXs, except for protein SJ05859 due to the short sequence length. The predicted 3D structure of SJ16632 showed that the five sites could interact with the iron. The simulated structures of six LOXs in the C1 clade and five LOXs in the C2 clade were superimposed to evaluate the goodness of fit of the overall topologies (Supplementary Figure S2). The result shows that the structure of the α-helix is conserved across all members, whereas structures in the loop regions are variable and less conserved, which might be due to relaxed functional constraints.
FIGURE 9. Modeled 3D structure of the S. japonica LOX protein (SJ16632) and the 3D structure of P. aeruginosa and soybean. Blue sphere: non-heme iron; inset: the predicted active site for binding iron.
Discussion
Lipoxygenases studies have been mainly focused on higher plants. However, with more genomic data becoming available for other photosynthetic lineages (collectively referred to as algae), marine algae have a high potential to provide crucial insights into the evolution of life on the planet (Ulvskov et al., 2013). Our study gives first insights into LOX family evolution in brown algae. LOXs in brown algae have the closest relationship with prokaryotes, suggesting that extant LOXs in both groups derive from the same ancestral protein. They were divided into two clades. In each clade, Ectocarpales and Laminariales are distributed in different subclades, and LOXs in S. japonica have experienced several rounds of duplications. How new genes emerge and gain functionality is a long-standing fundamental question in biology. Gene duplication is a major mechanism for the evolution of novel genes and the development of phenotypic innovation (Lynch and Conery, 2000; Weissenbach et al., 2017). Our study shows that LOX duplication has occurred more frequently in S. japonica compared to Ectocarpales, forming 11 LOX genes. Proteome size-dependent scaling for the number of metal-binding proteins was observed across all of life, that is, the metal-binding proportion of a proteome, such as LOX, increases with the expansion of the proteome (Dupont et al., 2010). The predicted number of genes is 15,891 in E. siliculosus, 15,166 in C. okamuranus, and 18,733 in S. japonica, respectively. To some extent, the expansion of LOXs in S. japonica could therefore be related to the overall gene content expansion. Gene duplications can have one of three outcomes. One copy may simply be silenced (non-functionalization); one copy may acquire a new, beneficial function, with the other copy retaining the original function (neofunctionalization); or each copy may retain a subset of the original function, their total capacity being equal to the single ancestral gene (subfunctionalization) (Ohno, 1970; Lynch and Conery, 2000). Genes produced via gene duplication within a genome are called paralogs. The long-term retention of duplicated genes suggests either subfunctionalization of paralogous genes, or neofunctionalization (Lynch and Walsh, 2007; Weissenbach et al., 2017). Expression divergence is often the first step toward functional diversification. It reduces the chance of non-functionalization of duplicate genes and thereby increases the chance of duplicate genes being retained (Lan et al., 2009; Liu et al., 2015). Clear divergence in expression patterns was observed among brown algal LOXs in all three life stages. Even the most recently diverged paralogs (e.g., Ec-20_003610 vs. Ec-20_003620, SJ05859 vs. SJ12567, SJ16633 vs. SJ16634) differed greatly in their expression, suggesting that those duplicates have a relatively high rate of functional divergence. The E. siliculosus genome has a total of 879 duplicated gene pairs. Of these, there are 28 pairs where both members exhibited sex-biased expression in the same sex (Lipinska et al., 2015). Both pairs of Ectocarpus LOXs were all male-biased, providing a first indication for the potential roles of brown algal LOXs in development. Gene duplications are connected to genome evolution and in turn, morphological complexity (Rensing, 2014). Unlike E. siliculosus, which has an alternative life cycle with two isomorphic generations (both sporophyte and gametophyte are filamentous), S. japonica has two distinctive generations with microscopic gametophytes and large and complex sporophytes, which have initial tissue differentiation, such as holdfast, blade, and stipe (Kawai et al., 2016). S. japonica LOXs exhibit high divergence of expression in different life stages. They are most highly expressed in sporophytes, and four genes are only expressed in sporophytes. Li et al. (2016) demonstrated that 15-LOX activation promotes cell cycle progression via the generation of the lipid mediator 15-hydroxyeicosatetraenoic acid (15-HETE), which is important for fibroblast proliferation (Sutliff, 2017). 15-HETE also exists in kelps (Ritter et al., 2008). In S. japonica LOX genes are expanded and highly expressed in large sporophytes. This could potentially promote the rapid growth and the development of complex morphology via the formation of sufficient quantities of various lipid mediators.
Diversifying selection has been shown to drive expansions or contractions of gene families (Shiu et al., 2006; Gingerich et al., 2007; Lan et al., 2009). Genes involved in stress response, metabolism, and reproduction are frequently reported to be under positive selection, i.e., selection in favor of new genetic variants (Oliver et al., 2010; Koester et al., 2012). Positive selection promotes non-synonymous substitutions, with the ratio of non-synonymous (dN) to synonymous (dS) nucleotide substitutions dN/dS > 1. However, positive selection accounts for only a fraction of the retention of paralogs (Demuth and Hahn, 2009; Lan et al., 2009). Nucleotide substitutions are either under neutral (dN/dS = 1) or purifying selection (dN/dS < 1) if they are deleterious for a population (Teng et al., 2017). Purifying selection has predominated in the LOX family of brown algae, and we found no evidence even for episodic positive selection. Subfunctionalization may well explain why duplicate genes continue to be under purifying selection even after duplication. Selection against deleterious mutations in different functional domains of duplicate genes may promote the conservation of both duplicates and enhance the fixation of functional gene pairs (Lan et al., 2009; Tanaka et al., 2009). This may also apply to the brown algal LOX family. Moreover, LOX-mediated fatty acid transformations contribute to the formation of many brown algal pheromones and structural isomers, which might play a role in the chemical defense of these algae (Pohnert and Boland, 2002). As many as 23 different oxylipins have been found in the brown algal kelp L. digitata (Ritter et al., 2008). The more relaxed selection constraints in giant kelps, such as S. japonica might, to some extent, contribute to the formation of these diverse lipid products. Gene expression levels in yeast are inversely correlated to protein divergence. This relationship is hypothesized to exercise indirect control over mutation rates (Drummond et al., 2005). Positively selected genes are expressed at lower levels than genes subject to neutral or purifying selection and tend to be expressed in restricted conditions or specific tissues (Kosiol et al., 2008; Koester et al., 2012; Teng et al., 2017). Moreover, gene pairs with high dN/dS value tend to be more divergent in expression (Mock et al., 2017). Consistent with those studies, we found that LOX gene pairs with higher dN/dS values have lower expression levels than other LOX gene pairs both in S. japonica and E. siliculosus, though no positive selection was detected. The two pairs of genes with highest dN/dS values were the least expressed in all the three life stages of S. japonica. Three of these four genes were under relaxed purifying selection and exclusively expressed in sporophytes, while the fourth gene (SJ12567) has a higher expression level in male gametophytes. Lipinska et al. (2015) suggested that this sex-biased expression was promoted by an event where gene duplication has freed one member of the pair from selection pressure to evolve a sex-specific function. This gene is also the only one that was highly expressed in gametophytes, indicating a potentially important role in this life stage.
Conclusion
Lipoxygenases constitute an important group of enzymes that metabolize PUFAs and produce various bioactive lipid mediators. LOXs in brown algae have closest relationships with those in prokaryotes. LOX gene expansion occurred in the macroalga S. japonica, mainly by tandem duplication and segmental duplication events. To investigate the evolutionary mechanisms responsible for the retention and subsequent functional divergence of duplicate LOX genes, we conducted molecular evolution, gene duplication, and expression analysis, and reconstructed the evolution history of this gene family in brown algae. LOXs in brown algae have undergone strong purifying selection. Expression of LOXs in E. siliculosus was up regulated under stress conditions. Expression also varied greatly between different LOX genes and life stages in brown algae, suggesting functional divergence since duplication. Our findings give new insights into the evolutionary history and functional diversification of brown algal LOX genes.
Author Contributions
NY planned and designed the work. LT, WH, and XF analyzed and interpretated the data of the work. LT wrote the manuscript. DX, XZ, and SD revised it critically for important intellectual content. All authors approved the version to be published. All authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by China Agriculture Research System (CARS-50); Qingdao Municipal Science and Technology plan project (15-8-2-3-hy; 17-1-1-96-jch); Taishan Scholars Funding; Central Public-interest Scientific Institution Basal Research Fund CAFS (2017HY-YJ01); Special Scientific Research Funds for Central Non-Profit Institutes, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences (20603022016001, 20603022016010); AoShan Talents Program (No. 2015ASTPES03); Director fund Supported by Qingdao National Laboratory for Marine Science and Technology, the Science Fund for Distinguished Young Scholars of Shandong Province (JQ201509); National key Research and Development Plan (2016YFC1402102-2); National Natural Science Foundation of China (41676145). We are grateful to Agnieszka Lipinska for providing LOX expression profiles in different generations of E. siliculosus, and to the ABIMS platform (Roscoff) for providing access to their S. japonica genome browser.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.02018/full#supplementary-material
Footnotes
References
Ackermann, J. A., Hofheinz, K., Zaiss, M. M., and Krönke, G. (2017). The double-edged role of 12/15-lipoxygenase during inflammation and immunity. Biochim. Biophys. Acta 1862, 371–381. doi: 10.1016/j.bbalip.2016.07.014
Bailey, T., Boden, M., Buske, F., Frith, M., Grant, C., Clementi, L., et al. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208. doi: 10.1093/nar/gkp335
Banthiya, S., Pekárová, M., Kuhn, H., and Heydeck, D. (2015). Secreted lipoxygenase from Pseudomonas aeruginosa exhibits biomembrane oxygenase activity and induces hemolysis in human red blood cells. Arch. Biochem. Biophys. 584, 116–124. doi: 10.1016/j.abb.2015.09.003
Barofsky, A., and Pohnert, G. (2007). Biosynthesis of polyunsaturated short chain aldehydes in the diatom Thalassiosira rotula. Org. Lett. 9, 1017–1020. doi: 10.1021/ol063051v
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate - a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300.
Biasini, M., Bienert, S., Waterhouse, A., Arnold, K., Studer, G., Schmidt, T., et al. (2014). SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 42, 252–258. doi: 10.1093/nar/gku340
Blée, E. (2002). Impact of phyto-oxylipins in plant defense. Trends Plant Sci. 7, 315–322. doi: 10.1016/S1360-1385(02)02290-2
Bouckaert, R., Heled, J., Kühnert, D., Vaughan, T., Wu, C.-H., Xie, D., et al. (2014). BEAST 2: a software platform for Bayesian evolutionary analysis. PLOS Comput. Biol. 10:e1003537. doi: 10.1371/journal.pcbi.1003537
Bueno, P., Soto, M. J., Rodríguez-Rosales, M. P., Sanjuan, J., Olivares, J., and Donaire, J. P. (2001). Time-course lipoxygenase, antioxidant enzyme activities and H2O2 accumulation during the early stages of Rhizobium-legume symbiosis. New Phytol. 152, 91–96. doi: 10.1046/j.0028-646x.2001.00246.x
Chen, Z., Chen, D., Chu, W., Zhu, D., Yan, H., and Xiang, Y. (2016). Retention and molecular evolution of lipoxygenase genes in modern rosid plants. Front. Genet. 7:176. doi: 10.3389/fgene.2016.00176
Chen, Z., Chen, X., Yan, H., Li, W., Li, Y., Cai, R., et al. (2015). The lipoxygenase gene family in poplar: identification, classification, and expression in response to MeJA treatment. PLOS ONE 10:e0125526. doi: 10.1371/journal.pone.0125526
Cock, J. M., Lieven, S., Pierre, R., Delphine, S., Allen, A. E., Grigoris, A., et al. (2010). The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature 465, 617–621. doi: 10.1038/nature09016
Crooks, G. E., Hon, G., Chandonia, J. M., and Brenner, S. E. (2004). WebLogo: a sequence logo generator. Genome Res. 14, 1188–1190. doi: 10.1101/gr.849004
Demuth, J. P., and Hahn, M. W. (2009). The life and death of gene families. Bioessays 31, 29–39. doi: 10.1002/bies.080085
d’Ippolito, G., Lamari, N., Montresor, M., Romano, G., Cutignano, A., Gerecht, A., et al. (2009). 15S-Lipoxygenase metabolism in the marine diatom Pseudo-nitzschia delicatissima. New Phytol. 183, 1064–1071. doi: 10.1111/j.1469-8137.2009.02887.x
Dittami, S. M., Scornet, D., Petit, J. L., Ségurens, B., Silva, C. D., Corre, E., et al. (2009). Global expression analysis of the brown alga Ectocarpus siliculosus (Phaeophyceae) reveals large-scale reprogramming of the transcriptome in response to abiotic stress. Genome Biol. 10:R66. doi: 10.1186/gb-2009-10-6-r66
Dobson, L., Reményi, I., and Tusnády, G. E. (2015). CCTOP: a consensus constrained topology prediction web server. Nucleic Acids Res. 43, 408–412. doi: 10.1093/nar/gkv451
Drummond, D. A., Bloom, J. D., Adami, C., Wilke, C. O., and Arnold, F. H. (2005). Why highly expressed proteins evolve slowly. Proc. Natl. Acad. Sci. U.S.A. 102, 14338–14343. doi: 10.1073/pnas.0504070102
Dupont, C. L., Butcher, A., Valas, R. E., Bourne, P. E., and Caetanoanollés, G. (2010). History of biological metal utilization inferred through phylogenomic analysis of protein structures. Proc. Natl. Acad. Sci. U.S.A. 107, 10567–10572. doi: 10.1073/pnas.0912491107
Fatima, N., Ramakrishna, B., Kumar, N., and Reddanna, P. (2017). Cloning, expression, purification of active recombinant human 12R-LOX enzyme and its inhibition by Acalypha indica extracts. PeerJ Preprints 5:e2906v1. doi: 10.7287/peerj.preprints.2906v1
Feng, B., Dong, Z., Xu, Z., An, X., Qin, H., Wu, N., et al. (2010). Molecular analysis of lipoxygenase (LOX) genes in common wheat and phylogenetic investigation of LOX proteins from model and crop plants. J. Cereal Sci. 52, 387–394. doi: 10.1016/j.jcs.2010.06.019
Finn, R. D., Coggill, P., Eberhardt, R. Y., Eddy, S. R., Mistry, J., Mitchell, A. L., et al. (2015). The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 44, D279–D285. doi: 10.1093/nar/gkv1344
Garreta, A., Val-Moraes, S. P., García-Fernández, Q., Busquets, M., Juan, C., Oliver, A., et al. (2013). Structure and interaction with phospholipids of a prokaryotic lipoxygenase from Pseudomonas aeruginosa. FASEB J. 27, k4811–4821. doi: 10.1096/fj.13-235952
Gerwick, W. H., Roberts, M. A., Vulpanovici, A., and Ballantine, D. L. (1999). “Biogenesis and biological function of marine algal oxylipins,” in Lipoxygenases and Their Metabolites, eds S. Nigam and C. R. Pace-Asciak (Boston, MA: Springer), 211–218. doi: 10.1007/978-1-4615-4861-4_20
Gingerich, D. J., Hanada, K., Shiu, S.-H., and Vierstra, R. D. (2007). Large-scale, lineage-specific expansion of a bric-a-brac/tramtrack/broad complex ubiquitin-ligase gene family in rice. Plant Cell 19, 2329–2348. doi: 10.1105/tpc.107.051300
Good, R. T., Gramzow, L., Battlay, P., Sztal, T., Batterham, P., and Robin, C. (2014). The molecular evolution of cytochrome P450 genes within and between Drosophila species. Genome Biol. Evol. 6, 1118–1134. doi: 10.1093/gbe/evu083
Gouet, P., Robert, X., and Courcelle, E. (2003). ESPript/ENDscript: extracting and rendering sequence and 3D information from atomic structures of proteins. Nucleic Acids Res. 31, 3320–3323. doi: 10.1093/nar/gkg556
Gruber, A., Rocap, G., Kroth, P. G., Armbrust, E., and Mock, T. (2015). Plastid proteome prediction for diatoms and other algae with secondary plastids of the red lineage. Plant J. 81, 519–528. doi: 10.1111/tpj.12734
Guo, S., Song, Z., Ma, R., Yang, Y., and Yu, M. (2017). Genome-wide identification and expression analysis of the lipoxygenase gene family during peach fruit ripening under different postharvest treatments. Acta Physiol. Plant. 39, k111–123. doi: 10.1007/s11738-017-2409-6
Hofmann, M., and Eichenberger, W. (1997). Lipid and fatty acid composition of the marine brown alga Dictyopteris membranacea. Plant Cell Physiol. 38, 1046–1052. doi: 10.1093/oxfordjournals.pcp.a029270
Holland, P. W., Marlétaz, F., Maeso, I., Dunwell, T. L., and Paps, J. (2017). New genes from old: asymmetric divergence of gene duplicates and the evolution of development. Philos. Trans. R. Soc. B 372, 20150480. doi: 10.1098/rstb.2015.0480
Hou, Y., Meng, K., Ye, H., Ban, Q., Wang, B., Suo, J., et al. (2015). The persimmon 9-lipoxygenase gene DkLOX3 plays positive roles in both promoting senescence and enhancing tolerance to abiotic stress. Front. Plant Sci. 6:1073. doi: 10.3389/fpls.2015.01073
Hu, B., Jin, J., Guo, A. Y., Zhang, H., Luo, J., and Gao, G. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31, 1296–1297. doi: 10.1093/bioinformatics/btu817
Imbs, A. B., Vologodskaya, A. V., Nevshupova, N. V., Khotimchenko, S. V., and Titlyanov, E. A. (2001). Response of prostaglandin content in the red alga Gracilaria verrucosa to season and solar irradiance. Phytochemistry 58, 1067–1072. doi: 10.1016/S0031-9422(01)00321-1
Jiang, Z. D., Ketchum, S. O., and Gerwick, W. H. (2000). 5-Lipoxygenase-derived oxylipins from the red alga Rhodymenia pertusa. Phytochemistry 53, 129–133. doi: 10.1016/S0031-9422(99)00445-8
Kalms, J., Banthiya, S., Yoga, E. G., Hamberg, M., Holzhutter, H.-G., Kuhn, H., et al. (2017). The crystal structure of Pseudomonas aeruginosa lipoxygenase Ala420Gly mutant explains the improved oxygen affinity and the altered reaction specificity. Biochim. Biophys. Acta 1862, 463–473. doi: 10.1016/j.bbalip.2017.01.003
Kawai, H., Hanyuda, T., and Uwai, S. (2016). Evolution and Biogeography of Laminarialean Kelps. Dordrecht: Springer. doi: 10.1007/978-94-017-7534-2_9
Koester, J. A., Swanson, W. J., and Armbrust, E. V. (2012). Positive selection within a diatom species acts on putative protein interactions and transcriptional regulation. Mol. Biol. Evol. 30, 422–434. doi: 10.1093/molbev/mss242
Kosiol, C., Vinar, T., Da, F. R., Hubisz, M. J., Bustamante, C. D., Nielsen, R., et al. (2008). Patterns of positive selection in six Mammalian genomes. PLOS Genet. 4:e1000144. doi: 10.1371/journal.pgen.1000144
Kousaka, K., Ogi, N., Akazawa, Y., Fujieda, M., Yamamoto, Y., Takada, Y., et al. (2003). Novel oxylipin metabolites from the brown alga Eisenia bicyclis. J. Nat. Prod. 66, 1318–1323. doi: 10.1021/np030049t
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Kumari, P., Reddy, R., and Jha, B. (2014). Quantification of selected endogenous hydroxy-oxylipins from tropical marine Macroalgae. Mar. Biotechnol. 16, k74–87. doi: 10.1007/s10126-013-9533-0
Küpper, F. C., Gaquerel, E., Cosse, A., Adas, F., Peters, A. F., Müller, D. G., et al. (2009). Free fatty acids and methyl jasmonate trigger defense reactions in Laminaria digitata. Plant Cell Physiol. 50, 789–800. doi: 10.1093/pcp/pcp023
Lan, T., Yang, Z.-L., Yang, X., Liu, Y.-J., Wang, X.-R., and Zeng, Q.-Y. (2009). Extensive functional diversification of the Populus glutathione S-transferase supergene family. Plant Cell 21, 3749–3766. doi: 10.1105/tpc.109.070219
Li, M., Li, L., Dunwell, J. M., Qiao, X., Liu, X., and Zhang, S. (2014). Characterization of the lipoxygenase (LOX) gene family in the Chinese white pear (Pyrus bretschneideri) and comparison with other members of the Rosaceae. BMC Genomics 15:444. doi: 10.1186/1471-2164-15-444
Li, Y., Zhang, L., Wang, X., Chen, M., Liu, Y., Xing, Y., et al. (2016). Elk-1-mediated 15-lipoxygenase expression is required for hypoxia-induced pulmonary vascular adventitial fibroblast dynamics. Acta Physiol. 218, 276–289. doi: 10.1111/apha.12711
Liavonchanka, A., and Feussner, I. (2006). Lipoxygenases: occurrence, functions and catalysis. J. Plant Physiol. 163, 348–357. doi: 10.1016/j.jplph.2005.11.006
Lipinska, A., Cormier, A., Luthringer, R., Peters, A. F., Corre, E., Gachon, C. M. M., et al. (2015). Sexual dimorphism and the evolution of sex-biased gene expression in the brown alga Ectocarpus. Mol. Biol. Evol. 32, 1581–1597. doi: 10.1093/molbev/msv049
Liu, H.-J., Tang, Z.-X., Han, X.-M., Yang, Z.-L., Zhang, F.-M., Yang, H.-L., et al. (2015). Divergence in enzymatic activities in the soybean GST supergene family provides new insight into the evolutionary dynamics of whole-genome duplicates. Mol. Biol. Evol. 32, 2844–2859. doi: 10.1093/molbev/msv156
Liu, S. Q., Liu, X. H., and Jiang, L. W. (2011). Genome-wide identification, phylogeny and expression analysis of the lipoxygenase gene family in cucumber. Genet. Mol. Res. 10, 2613–2636. doi: 10.4238/2011.October.25.9
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Lynch, M., and Conery, J. S. (2000). The evolutionary fate and consequences of duplicate genes. Science 290, 1151–1155. doi: 10.1126/science.290.5494.1151
Lynch, M., and Walsh, B. (2007). The Origins of Genome Architecture. Sunderland, MA: Sinauer Associates Sunderland.
Melan, M. A., Dong, X., Endara, M. E., Davis, K. R., Ausubel, F. M., and Peterman, T. K. (1993). An Arabidopsis thaliana lipoxygenase gene can be induced by pathogens, abscisic acid, and methyl jasmonate. Plant Physiol. 101, 441–450. doi: 10.1104/pp.101.2.441
Minor, W., Steczko, J., Stec, B., Otwinowski, Z., Bolin, J. T., Walter, R., et al. (1996). Crystal structure of soybean lipoxygenase L-1 at 1.4 Å resolution. Biochemistry 35, 10687–10701. doi: 10.1021/bi960576u
Mock, T., Otillar, R. P., Strauss, J., McMullan, M., Paajanen, P., Schmutz, J., et al. (2017). Evolutionary genomics of the cold-adapted diatom Fragilariopsis cylindrus. Nature 541, 536–540. doi: 10.1038/nature20803
Molina, A., and Donaire, J. P. (2002). Involvement of endogenous salicylic acid content, lipoxygenase and antioxidant enzyme activities in the response of tomato cell suspension cultures to NaCl. New Phytol. 156, 409–415. doi: 10.1046/j.1469-8137.2002.00527.x
Nishitsuji, K., Arimoto, A., Iwai, K., Sudo, Y., Hisata, K., Fujie, M., et al. (2016). A draft genome of the brown alga, Cladosiphon okamuranus, S-strain: a platform for future studies of ‘mozuku’ biology. DNA Res. 23, 561–570. doi: 10.1093/dnares/dsw039
Oliver, T. A., Garfield, D. A., Manier, M. K., Haygood, R., Wray, G. A., and Palumbi, S. R. (2010). Whole-genome positive selection and habitat-driven evolution in a shallow and a deep-sea Urchin. Genome Biol. Evol. 2, 800–814. doi: 10.1093/gbe/evq063
Petersen, T. N., Brunak, S., Von, H. G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701
Podolyan, A., White, J., Jordan, B., and Winefield, C. (2010). Identification of the lipoxygenase gene family from Vitis vinifera and biochemical characterisation of two 13-lipoxygenases expressed in grape berries of Sauvignon Blanc. Funct. Plant Biol. 37, 767–784. doi: 10.1071/FP09271
Pohnert, G. (2002). Phospholipase A2 activity triggers the wound-activated chemical defense in the diatom Thalassiosira rotula. Plant Physiol. 129, 103–111. doi: 10.1104/pp.010974
Pohnert, G., and Boland, W. (2002). The oxylipin chemistry of attraction and defense in brown algae and diatoms. Nat. Prod. Rep. 19, 108–122. doi: 10.1039/a806888g
Porta, H., and Rocha-Sosa, M. (2001). Lipoxygenase in bacteria: a horizontal transfer event? Microbiology 147, 3199–3200. doi: 10.1099/00221287-147-12-3199
Potin, P., Bouarab, K., Salaün, J.-P., Pohnert, G., and Kloareg, B. (2002). Biotic interactions of marine algae. Curr. Opin. Plant Biol. 5, 308–317. doi: 10.1016/S1369-5266(02)00273-X
Rensing, S. A. (2014). Gene duplication as a driver of plant morphogenetic evolution. Curr. Opin. Plant Biol. 17, 43–48. doi: 10.1016/j.pbi.2013.11.002
Ritter, A., Dittami, S. M., Goulitquer, S., Correa, J. A., Boyen, C., Potin, P., et al. (2014). Transcriptomic and metabolomic analysis of copper stress acclimation in Ectocarpus siliculosus highlights signaling and tolerance mechanisms in brown algae. BMC Plant Biol. 14:116. doi: 10.1186/1471-2229-14-116
Ritter, A., Goulitquer, S., Salaün, J. P., Tonon, T., Correa, J. A., and Potin, P. (2008). Copper stress induces biosynthesis of octadecanoid and eicosanoid oxygenated derivatives in the brown algal kelp Laminaria digitata. New Phytol. 180, 809–821. doi: 10.1111/j.1469-8137.2008.02626.x
Roux, J., Privman, E., Moretti, S., Daub, J. T., Robinson-Rechavi, M., and Keller, L. (2014). Patterns of positive selection in seven ant genomes. Mol. Biol. Evol. 31, 1661–1685. doi: 10.1093/molbev/msu141
Shiu, S.-H., Byrnes, J. K., Pan, R., Zhang, P., and Li, W.-H. (2006). Role of positive selection in the retention of duplicate genes in mammalian genomes. Proc. Natl. Acad. Sci. U.S.A. 103, 2232–2236. doi: 10.1073/pnas.0510388103
Sutliff, R. L. (2017). New insights into mechanisms regulating pulmonary adventitial fibroblast proliferation. Acta Physiol. 219, 17–19. doi: 10.1111/apha.12753
Suyama, M., Torrents, D., and Bork, P. (2006). PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 34, W609–W612. doi: 10.1093/nar/gkl315
Tanaka, K. M., Takahasi, K. R., and Takano-Shimizu, T. (2009). Enhanced fixation and preservation of a newly arisen duplicate gene by masking deleterious loss-of-function mutations. Genet. Res. 91, 267–280. doi: 10.1017/S0016672309000196
Teng, L., Fan, X., Xu, D., Zhang, X., Mock, T., and Ye, N. (2017). Identification of genes under positive selection reveals differences in evolutionary adaptation between brown-algal species. Front. Plant Sci. 8:1429. doi: 10.3389/fpls.2017.01429
Ulvskov, P., Paiva, D. S., Domozych, D., and Harholt, J. (2013). Classification, naming and evolutionary history of glycosyltransferases from sequenced green and red algal genomes. PLOS ONE 8:e76511. doi: 10.1371/journal.pone.0076511
Weissenbach, J., Ilhan, J., Bogumil, D., Hülter, N., Stucken, K., and Dagan, T. (2017). Evolution of chaperonin gene duplication in Stigonematalean cyanobacteria (Subsection V). Genome Biol. Evol. 9, 241–252. doi: 10.1093/gbe/evw287
Wichard, T., Poulet, S. A., Halsband-Lenk, C., Albaina, A., Harris, R., Liu, D., et al. (2005). Survey of the chemical defence potential of diatoms: screening of fifty one species for alpha,beta,gamma,delta-unsaturated aldehydes. J. Chem. Ecol. 31, 949–958. doi: 10.1007/s10886-005-3615-z
Xu, D., Schaum, C. E., Lin, F., Sun, K., Munroe, J. R., Zhang, X. W., et al. (2017). Acclimation of bloom-forming and perennial seaweeds to elevated pCO2 conserved across levels of environmental complexity. Glob. Change Biol. 23, 4828–4839. doi: 10.1111/gcb.13701
Yang, Z. (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586–1591. doi: 10.1093/molbev/msm088
Yang, Z., Nielsen, R., and Hasegawa, M. (1998). Models of amino acid substitution and applications to mitochondrial protein evolution. Mol. Biol. Evol. 15, k1600–1611. doi: 10.1093/oxfordjournals.molbev.a025888
Ye, N., Zhang, X., Miao, M., Fan, X., Zheng, Y., Xu, D., et al. (2015). Saccharina genomes provide novel insight into kelp biology. Nat. Commun. 6:6986. doi: 10.1038/ncomms7986
Zambounis, A., Gaquerel, E., Strittmatter, M., and Salaün, J. P. (2012). Prostaglandin A2 triggers a strong oxidative burst in Laminaria: a novel defense inducer in brown algae? Algae 27, 21–32. doi: 10.4490/algae.2012.27.1.021
Keywords: lipoxygenase, brown algae, gene duplication, selection pressure, functional divergence
Citation: Teng L, Han W, Fan X, Xu D, Zhang X, Dittami SM and Ye N (2017) Evolution and Expansion of the Prokaryote-Like Lipoxygenase Family in the Brown Alga Saccharina japonica. Front. Plant Sci. 8:2018. doi: 10.3389/fpls.2017.02018
Received: 12 August 2017; Accepted: 13 November 2017;
Published: 28 November 2017.
Edited by:
Mathew G. Lewsey, La Trobe University, AustraliaReviewed by:
Manoj Kumar, University of Technology Sydney, AustraliaKui Lin, Beijing Normal University, China
Copyright © 2017 Teng, Han, Fan, Xu, Zhang, Dittami and Ye. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Naihao Ye, yenh@ysfri.ac.cn