Toshiki Uji
Takuya Kandori
Hiroyuki Mizuta- Laboratory of Aquaculture Genetics and Genomics, Division of Marine Life Science, Faculty of Fisheries Sciences, Hokkaido University, Hakodate, Japan
Saccharina japonica, a significant brown macroalga in the Pacific Ocean, serves as a food source and industrial material. In aquaculture, collecting mature sporophytes for seedling production is essential but challenging due to environmental changes. In this study, transcriptomic analysis of vegetative and sorus tissues was done to identify differentially expressed genes (DEGs) and enhance our understanding of sorus formation regulation in S. japonica. KEGG pathway and Gene Otology (GO) analysis revealed that upregulated DEGs were involved in folate biosynthesis, riboflavin metabolism, and amino acid biosynthesis. In addition, the upregulation of genes associated with cell wall remodeling, such as mannuronan C-5-epimerases, vanadium-dependent haloperoxidases, and NADPH oxidase, was observed in sorus parts. Meanwhile, downregulated DEGs in sorus portions included genes related to chloroplast function. These findings will help us understand the regulatory mechanisms behind sorus formation in S. japonica and extracellular matrix remodeling in brown algae.
1 Introduction
Saccharina japonica, an ecologically significant brown macroalga, is a major component of productive beds on the northwest coast of the Pacific Ocean (Xu et al., 2015). It is one of the most economically important seaweeds in aquaculture, extensively utilized as a food source and raw industrial material (Tseng, 2001). In Japan, S. japonica is an essential ingredient for making dashi (soup stock) as “kombu” in Japanese cuisine.
In S. japonica aquaculture, high-quality wild mature sporophytes are collected from the periphery of aquaculture sites in autumn to produce seedlings. However, due to recent marine environment changes, it has become increasingly difficult to collect mature S. japonica sporophytes that release active zoospores (Yotsukura et al., 2022). Certain hybrid cultivars of S. japonica, despite having excellent agronomic traits, also exhibit late maturation (Liu et al., 2023). This hinders the hybrid cultivars from producing spores at the time of seedling production, limiting their practical application. Thus, understanding sporogenesis regulation is crucial for sustainable aquaculture in S. japonica.
Previous studies investigated the physiological alterations during sorus formation, such as decreased photosynthesis and increased respiration (Nimura and Mizuta, 2001). An exogenous application of abscisic acid (ABA) to sporophyte disks of S. japonica promoted sorus formation, and the amount of ABA in sorus tissue was significantly increased compared to that in vegetative tissue (Nimura and Mizuta, 2002). Sorus also accumulates high levels of phenolics, silicon, and iodoperoxidase (IPO) activity compared to its vegetative parts to increase its defensive capacity (Mizuta and Yasui, 2010; Mizuta and Yasui, 2012) (See review (Uji and Mizuta, 2022)). In addition to physiological studies, systematic analysis identified the genes encoding enzymes for the biosynthesis of cell wall carbohydrates (including alginate, fucoidan, and cellulose) and cytoplasm storage carbohydrates (mannitol, laminarin, and trehalose) during the sporophyte development (Zhang et al., 2020). However, molecular biological findings that complement physiological findings in the sporogenesis of S. japonica are still lacking.
In this study, a comparative transcriptomic analysis was performed on vegetative and sorus tissues of S. japonica sporophytes to investigate the differentially expressed genes (DEGs) during reproduction. The findings of this study will enhance our understanding of the regulatory mechanisms, especially cell wall remodeling involved in sorus formation of S. japonica.
2 Materials and methods
2.1 Algal materials
Mature S. japonica sporophytes cultured in a farming area off the coast of Minamikayabe, Hakodate City, Hokkaido, Japan, were collected in October 2023 and transported in a cool box with refrigerants to our laboratory. After 1 h, the disks (3 cm in diameter) formed sorus (mature) and nonsorus portions that did not contain meristematic regions (vegetative), which were cut by a cork borer from the marginal parts sited at 100–200 cm from the stipe-blade transition along each sporophyte (~3 m long). The release of the zoospores from sorus portions was confirmed. The disks were wiped with a paper towel and washed with sterilized seawater to remove the attached organisms. The disks were cultured in 800 mL of sterile vitamin-free Provasoli’s enriched seawater (Provasoli, 1968) at 10°C under 10–20 μmol photons m−2 s−1 (12 h light/12 h dark cycle) to reduce the effect of dissection on gene expression. After 1 week, the disks were harvested at 10:00 am, immediately frozen with liquid nitrogen, and stored at −80°C until RNA extraction.
2.2 RNA extraction
Total RNA extraction was conducted using a combination of a CTAB-based method and the RNeasy Plant Mini Kit (Qiagen, Hilden, Germany), following the protocol by (Heinrich et al., 2012) with minor modifications. Frozen sporophytes were ground in liquid nitrogen with a mortar and pestle and were transferred to tubes. Then, 1 ml of extraction buffer (composed of 2% CTAB, 1 M NaCl, 100 mM Tris pH 8, 50 mM EDTA, pH 8, 3% Polyvinylpolypyrrolidone) and 40 µl of 1M DTT were added and thoroughly mixed. The resulting mixture was incubated at 45°C for 10 mins. Then, one volume of chloroform: isoamylalcohol (24:1) was vigorously mixed for 10 mins. The tubes were centrifuged for 20 min at 20°C and 12,000 g. Subsequently, 600 µl of the aqueous phase was carefully transferred into a new tube. An additional step involved the gentle addition and mixing of 0.3 volumes of 100% EtOH by inverting the tube. A second chloroform extraction was performed by adding one volume of chloroform: isoamylalcohol (24:1). Following centrifugation, 500 µl of the resulting supernatant was transferred to a new tube. Total RNA extraction was accomplished using an RNeasy Plant Mini Kit, following the manufacturer’s instructions. The extracted RNA was purified using a TURBO DNA-free kit (Invitrogen/Life Technologies, Carlsbad, CA) to obtain DNA-free RNA. RNA samples’ quantity and integrity were assessed using a NanoDrop™ 2000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA) and an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).
2.3 RNA sequencing analysis
Six libraries of complementary DNA (two conditions: vegetative and mature × three replicates) for S. japonica were constructed and subsequently sequenced using an Illumina NovaSeq 6000 instrument at Rhelixa Inc. The obtained reads were trimmed for low-quality reads and adapter sequences using fastp (Chen et al., 2018). After trimming, STAR (Dobin et al., 2013) was used to map high-quality reads to an in-house gene model of S. japonica, which was constructed with a reference genome (ASM882872v1) using BRAKER 2.1.6 (Brůna et al., 2021). The normalized expression of each gene was calculated using RSEM (Li and Dewey, 2011) as transcripts per million. DEGs between vegetative and mature were identified using edgeR (Robinson et al., 2010) on the false discovery rate significance score < 0.05 and |log2 fold change| > 1.
2.4 Functional annotation of RNA-Seq data
The obtained DEGs were annotated as follows. First, gene coding regions of S. japonica were inferred by Braker (Brůna et al., 2021) to use the genome sequences for functional annotation. The eukaryotic protein database OrthoDB v10 (Kriventseva et al., 2019; Zdobnov et al., 2021) was used to predict S. japonica gene coding regions. The protein sequences of the obtained gene coding regions were functionally annotated using EnTAP (Hart et al., 2020) and InterProScan (Jones et al., 2014).
To assess the biological significance of the DEGs, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted. GO terms were assigned to all genes using eggNOG-mapper v2 online (http://eggnog-mapper.embl.de/) (Huerta-Cepas et al., 2019; Cantalapiedra et al., 2021) with the default parameters, except that the min_hit_e-value was set to 0.05. The topGO (Alexa and Rahnenfuhrer, 2023) R package was used for GO enrichment analysis, and GO terms with P < 0.05 were considered significantly enriched in the DEGs. A KOBAS 3.0 software (Bu et al., 2021) was used for KEGG annotation and enrichment analysis based on the KEGG PATHWAY database (https://www.genome.jp/kegg/pathway.html) (Kanehisa and Goto, 2000). Pathways with a corrected p-value (q-value) < 0.05 were defined as significantly enriched pathways for DEGs.
2.5 Quantitative PCR
First-strand cDNA was synthesized from 0.5 µg of total RNA (same RNA used for RNA-seq) using a PrimeScript II 1st strand cDNA Synthesis Kit (TaKaRa Bio, Shiga, Japan). The cDNA was diluted 10-fold for qPCR analysis, and 1.0 µl of the diluted cDNA was used as a template in a 20 μL reaction volume using KOD SYBR® qPCR Mix (TOYOBO, Osaka, Japan), following the manufacturer’s instructions. Real-time PCR was performed with a LightCycler® 96 System (Roche Diagnostics, Basel, Switzerland) under the following conditions: 2 min at 98°C followed by 40 cycles of 10 s at 98°C, 10 s at 55°C, and 30 s at 68°C. The mRNA levels were calculated using the 2−△△Ct method and normalized to levels of the 18S ribosomal RNA gene. The relative expression level was calculated as a ratio of the mRNA level to the transcription level of vegetative samples. All experiments were performed in triplicate. Supplementary Table S1 lists the primers used in this study.
3 Results and discussion
3.1 Identification of DEGs related to reproduction
Using RNA-seq, the transcripts in S. japonica were compared between sorus (mature) and nonsorus portions (vegetative) to identify candidate genes regulating sorus formation. Raw data generated by sequencing ranged from 18.0–24.2 million reads per sample. After filtering, 17.6–23.9 million clean reads were obtained, and the mapping rate was 76.9–86.0%. A summary of the obtained RNA sequencing datasets and mapping rates of clean reads is shown in Table 1.
Table 1 Summary of transcriptome analysis in Saccharina japonica.
A total of 629 DEGs were obtained between mature and vegetative, including 354 upregulated and 275 downregulated genes in sorus portions (Supplementary Table S2). Six DEGs (three each for upregulated and downregulated genes) were selected for qPCR analysis to validate the accuracy of the RNA-seq data. As shown in Figure 1, expression levels of the selected genes were similar in the RT-qPCR and RNA-seq analyses, indicating that the RNA-seq results were reliable. Shown in Tables 2, 3 are the representative genes found to be differentially expressed in mature sporophytes.
Figure 1 Validation of RNA sequencing (RNA-seq) data using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Six representative genes were selected to validate the RNA-seq data by RT-qPCR. The pink bars represent the mean log2-fold change obtained by RT-qPCR, and the blue bars represent the RNA-seq data. Results are presented as relative expression compared to that in vegetative tissues. The primers used for RT-qPCR are listed in Supplementary Table S1.
Table 2 Selected upregulated genes with possible roles in the sporogenesis of Saccharina japonica.
Table 3 Selected downregulated genes with possible roles in the sporogenesis of Saccharina japonica.
3.2 KEGG enrichment and GO analysis
KEGG enrichment analysis indicated that upregulated DEGs in sorus portions could be categorized into several pathways, including folate biosynthesis; riboflavin metabolism; and amino acid biosynthesis, such as serine (Figure 2). Folates are indispensable components of metabolism in all living organisms; they play as donors and acceptors of one-carbon groups in one-carbon transfer reactions that participate in the formation of numerous important biomolecules, such as nucleic acids, pantothenate (vitamin B5), and amino acids (Gorelova et al., 2017). The strong antioxidant properties of folates can be considered a key factor in elucidating their role in enhancing plant tolerance to diverse abiotic stresses and preventing oxidative damage (Alsamadany et al., 2022). Folates also regulate cellular and molecular events that affect plant growth and development, including cell division (Gorelova et al., 2017). In folate metabolism, serine is crucial for the regulation of methyl group transfer by providing tetrahydrofolate metabolism with C1 units (Ros et al., 2014). Previous studies suggest that serine is involved in the biosynthesis of several biomolecules required for cell proliferation (Ros et al., 2014). Riboflavin (vitamin B2) is a vital component required for fundamental metabolism and a precursor of the coenzymes, FAD and FMN (Jiadkong et al., 2024). Riboflavin metabolism includes antioxidant activity, cell signaling, and coenzyme function (Jiadkong et al., 2024). During sorus formation, upregulation of folate and amino acid biosynthesis, and riboflavin metabolism may play an important role in cell proliferation for sporulation and antioxidant system for reactive oxygen species (ROS) homeostasis.
Figure 2 Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the upregulated differentially expressed genes (DEGs). The degree of enrichment increased as the rich factor increased. Larger dots indicate higher numbers of differential genes enriched by the pathway.
Conversely, KEGG enrichment analysis revealed that downregulated DEGs in sorus portions were categorized into pathways, including photosynthesis-antenna proteins, carbon fixation, and secondary metabolite biosynthesis (Figure 3).
Figure 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the downregulated differentially expressed genes (DEGs). The degree of enrichment increased as the rich factor increased. Larger dots indicate higher numbers of differential genes enriched by the pathway.
To elucidate the biological processes, molecular functions, and cellular components associated with DEGs, GO analysis was performed using eggNOG-mapper and topGO. In upregulated DEGs, GO terms related to serine biosynthesis, riboflavin biosynthesis, and cell wall metabolism were enriched (Figure 4). In downregulated DEGs, GO terms related to photosynthesis and chloroplast components were enriched (Figure 5). Consequently, KEGG enrichment analysis and GO analysis revealed that upregulated genes were associated with vitamin and amino acid biosynthesis and cell wall metabolism. In contrast, photosynthesis-related genes were downregulated.
Figure 4 Numbers of enriched Gene Otology (GO) terms for upregulated differentially expressed genes (DEGs). GO terms are presented for three main categories: biological processes, molecular functions, and cellular components. Each GO term is listed in ascending order of p-value.
Figure 5 Numbers of enriched Gene Otology (GO) terms for the downregulated differentially expressed genes (DEGs). GO terms are presented for three main categories: biological processes, molecular functions, and cellular components. Each GO term is listed in ascending order of p-value.
3.3 Upregulation of ECM remodeling
DEG analysis revealed abundant transcripts in the sorus parts, including genes associated with the extracellular matrix (ECM), such as cell wall remodeling (Table 2; Figure 6). The ECM is a complex supramolecular network that imparts both rigidity and flexibility to multicellular tissues (Hynes, 2009). Beyond its structural role, it regulates development and protects cells from biotic and abiotic stresses (Kim et al., 2011). The ECMs of macroalgae, commonly referred to as the cell wall, consist of complex assemblages of cellulose, various hemicelluloses, and unique sulfated polysaccharides (Kloareg et al., 2021). The primary cell wall components in brown algae are anionic polysaccharides, specifically alginates and fucose-containing sulfated polysaccharides (Chi et al., 2018). Alginate plays a more prominent role in supporting cell structure than fucoidan (Kloareg et al., 2021).
Figure 6 Relative expression levels of ECM related genes in S. japonica during sorus formation. RNA samples were prepared from vegetative tissue (VT) and mature tissue formed zoosporangia (MT). Expression levels were assessed using the 18SrRNA gene for normalization. Results are presented as relative expression and compared with that in VT. Data are presented as means ± standard deviations (n = 3). Asterisks indicate significant differences at p < 0.05 between VT and MT (The Mann–Whitney U test).
Alginate is composed of two hexuronic acids: D-mannuronic acid (M) and L-guluronic acid (G). These are arranged in unbranched homopolymeric regions of M and G blocks, interspaced with a random arrangement of both monomers (MG blocks) (Heyraud et al., 1996). Alginate with a low M/G ratio exhibits superior mechanical properties compared to that with a high M/G ratio. Alginate with a high G-content also has a higher antibacterial activity than its low G-content counterpart (Ci et al., 2021). In brown algae, the composition of alginic acid varies with season, region, tissue, and age (Haug et al., 1974; McKee et al., 1992).
Previous research has shown a biosynthetic pathway for alginate in brown algae (Nyvall et al., 2003). The final step involves the epimerization of D-mannuronic residues into L-guluronic residues within the polymer chain, a reaction catalyzed by mannuronan C-5-epimerases (MC5E). Considering both their chemical structure and biological function, which is to provide strength and flexibility to the algal tissue, alginates can be regarded as functional analogs of the pectins found in higher plants. Consequently, mannuronan C-5-epimerases are likely functionally analogous to plant pectin methylesterases in controlling the cell wall matrix texture (Nyvall et al., 2003). In RNA-seq analysis, we found upregulated 8 DEGs encoding enzymes of MC5E (Figure 6), and this result showed similarity to previous results (Zhang et al., 2020). G-rich alginate in brown seaweed serves as a skeleton to provide stiffness or elasticity according to the environment, similar to pectins and cellulose in terrestrial plants (Nyvall et al., 2003). The upregulation of various MC5Es in the sorus parts is likely necessary to specifically customize the relative contents and distributions of G blocks, M blocks, and MG blocks in alginate chains.
In addition to MC5Es, cell wall remodeling-related genes were identified as sorus-preferential genes. DEGs encoding vanadium-dependent haloperoxidases (vHPOs) (g1659), (g277), NADPH oxidase (g3247), and superoxide dismutase (SOD) (g10745) were found to be upregulated in the sorus parts (Table 2). Saccharina japonica has 89 vHPOs, consisting of 21 bromoperoxidases (vBPOs) and 68 IPOs (Liu et al., 2019), emphasizing their role in processes such as chemical defense. vHPOs also play a role in oxidative cross-linking that facilitates the formation of phenolic polymers and their complexation with alginates, contributing to cell wall rigidification (Tarakhovskaya et al., 2015). Previous studies showed that phenolic substances, called phlorotannins, oxidized by vHPOs undergo self-assembly and form a macromolecular cluster with alginate, wherein the phenolic substances are encapsulated within the gel network of alginate (Bitton et al., 2006). The dynamics of vHPO activity during Fucus vesiculosus embryogenesis were highly synchronized with H2O2 content changes. Because H2O2 serves as a vHPO substrate, elevated ROS levels may serve as a prerequisite for an increased enzyme activity (Lemesheva et al., 2020). Genes involved in vHPO activity, including NADPH oxidase and SOD, were also found in upregulated DEGs of mature S. japonica sporophytes. NADPH oxidase, also known as respiratory burst oxidase homolog, is a well-studied enzymatic ROS-producing system (Marino et al., 2012). SODs catalyze the dismutation of O2− produced by NADPH oxidase into O2 and H2O2 (Sagi and Fluhr, 2006). In previous studies, substantial intracellular ROS production was observed in S. japonica sorus, especially zoosporangium and paraphyzes. The high abundance of transcripts encoding NADPH oxidase and SOD in the sorus may contribute to the control of ROS levels for cell wall remodeling in the sorus.
The cell wall undergoes remodeling in a tightly regulated and polarized manner, a process primarily controlled by the cell wall integrity (CWI) signaling pathway (Levin, 2011). CWI signaling activation regulates the production of various carbohydrate cell wall polymers, as well as their polarized delivery to the site of cell wall remodeling. Proteins with a cell wall integrity and stress response component (WSC) domain were first described as cell surface sensors involved in detecting and transmitting cell wall status to the CWI signaling pathway in Saccharomyces cerevisiae (Verna et al., 1997). S. cerevisiae Wsc1 accumulates to sites of enhanced mechanical stress through reduced lateral diffusivity, mediated by the binding of its extracellular WSC domain to cell wall polysaccharides (Neeli-Venkata et al., 2021). RNA-seq analysis showed an upregulated DEG (g17300) containing a WSC domain in the sorus parts. Although there is a lack of information on the WSC domain in macroalgae, g17300 may be an important player in cell wall remodeling during S. japonica sorus formation.
Four upregulated DEGs were found containing a spondin domain as genes related to the ECM. The spondin family, which includes F-spondin and Mindin, comprises molecules attached to ECM (Feinstein et al., 1999). The F-spondin molecule consists of approximately 800 amino acids containing domains homologous to reelin, FS domain, and multiple TSR repeats (Higashijima et al., 1997), while Mindin contains an FS domain and one TSR domain (Umemiya et al., 1997). Structural studies suggest that the FS domain, exhibiting a homologous structure similar to that of the C2 domain, functions as a membrane-targeting module through Ca2+-dependent mechanisms (Tan and Lawler, 2011). There is little information on the spondin family in algae, but four genes that contain only the FS domain (named NySPL1–4) were identified from the red alga Pyropia yezoensis genome (Uji et al., 2022). NySPLs have similar secondary structures to that of FS domains from animals, and their transcripts increased in mature thalli treated with 1-Aminocyclopropane-1-carboxylic acid (ACC), which is a gametogenesis inducer. Similarly, the spondin domain from S. japonica only contains the FS domain (SjSPL1–4), and their transcripts increased in the sorus parts. Thus, FS domain-containing proteins play a crucial role in the formation of reproductive cells in red and brown macroalgae.
ECM plays crucial roles in protecting cells from biotic and abiotic stresses. The present RNA-seq analysis showed an upregulated DEG (g11685) encoding lipoxygenase (LOX) in the sorus parts. In plants and algae, the LOX pathway has been proposed to play a key role in their defense via major oxylipins. Oxylipins derived from the oxidation of polyunsaturated fatty acids induce the establishment of resistance in the kelp Laminaria digitata against infection by its brown algal endophyte Laminariocolax tomentosoides (Küpper et al., 2009). In higher plants, oxylipins from the LOX pathway function in cell wall modifications required for root development and pathogen arrest (Vellosillo et al., 2007). The link between the LOX pathway and cell wall modification in brown algae should be investigated.
3.4 Downregulation of chloroplast function
In RNA-seq analysis, we found downregulated DEGs associated with chloroplast function, such as Tab2 family RNA-binding protein, Tic20, and light-harvesting complex protein (Table 3). In Chlamydomonas reinhardtii, Tab2 was identified as an RNA-binding protein required for the translation of the photosystem I subunit, and Arabidopsis Tab2 is involved in the signaling pathway of light-controlled synthesis of photosystem proteins during early plant development (Dauvillée et al., 2003; Barneche et al., 2006). Tic20 is a crucial component of the protein-conducting channel within the inner membrane preprotein translocon. In S. japonica, Tic20 is specifically positioned in the innermost membrane of the chloroplast, indicating its involvement in facilitating protein transport within the chloroplast (Chen et al., 2019). Light-harvesting complex proteins are known to be involved both in collecting light energy for driving the primary photochemical reactions of photosynthesis (Rochaix and Bassi, 2019). The reduced expression of genes related to chloroplast function is consistent with previous findings that the photosynthetic activity of S. japonica sporophyte was lower in the fertile parts than in the sterile parts (Nimura and Mizuta, 2001). One possible cost of photosynthesis in reproductive structures is the increased damage to DNA by the necessary exposure of reproductive structures to photosynthetically active, and hence (in nature) UV-B, radiation, and the increased potential for the production of ROS (Raven and Griffiths, 2015). Finally, we propose a possible model of the mechanisms of sorus development in S. japonica (Figure 7).
Figure 7 Schematic representation of putative mechanism during S. japonica sorus formation. LOX, lipoxygenase; SOD, superoxide dismutase; vHPO, vanadium-dependent haloperoxidase; SPL, spondin-like, WSC, cell wall integrity and stress response component; MC5E, mannuronan C5 epimerase.
Conclusion
The transcriptomic analysis suggests the upregulation of genes associated with vitamins and amino acid biosynthesis, ECM remodeling, and cell wall metabolism in S. japonica sorus portions, with downregulation of photosynthesis-related genes. Regarding ECM remodeling, DEGs associated with the composition of alginic acid and cell wall cross-linking were upregulated during sorus development. Candidate genes involved in ECM remodeling signaling pathways were also identified. Thus, this study on sorus formation mechanisms in S. japonica is important for a sustainable kelp production and to better understand the regulatory mechanisms of ECM remodeling in brown algae.
Data availability statement
The data presented in the study are deposited in DDBJ under the BioProject Accession number: PRJDB17633.
Author contributions
TU: Writing – review & editing, Writing – original draft. TK: Writing – review & editing, Writing – original draft. HM: Writing – review & editing, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the “Regional University and Local Industrial Vitalization Grant Project” (Hakodate City) adopted by the Cabinet Office (Japan).
Acknowledgments
We sincerely thank Mr. Ikuya Miyajima (Usujiri Fisheries Station, Field Science Center, Hokkaido University) for helping with the cultivation of kelp during this study.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1417582/full#supplementary-material
References
Alexa, A., Rahnenfuhrer, J. (2023). topGO: Enrichment analysis for gene ontology (R package version 2.54.0). doi: 10.18129/B9.bioc.topGO
Alsamadany, H., Mansour, H., Elkelish, A., Ibrahim, M. F. M. (2022). Folic acid confers tolerance against salt stress-induced oxidative damages in snap beans through regulation growth, metabolites, antioxidant machinery and gene expression. Plants 11, 1459. doi: 10.3390/plants11111459
Barneche, F., Winter, V., Crèvecoeur, M., Rochaix, J. D. (2006). ATAB2 is a novel factor in the signalling pathway of light-controlled synthesis of photosystem proteins. EMBO J. 25, 5907–5918. doi: 10.1038/sj.emboj.7601472
Bitton, R., Ben-Yehuda, M., Davidovich, M., Balazs, Y., Potin, P., Delage, L., et al. (2006). Structure of algal-born phenolic polymeric adhesives. Macromol. Biosci. 6, 737–746. doi: 10.1002/mabi.200600073
Brůna, T., Hoff, K. J., Lomsadze, A., Stanke, M., Borodovsky, M. (2021). BRAKER2: Automatic eukaryotic genome annotation with GeneMark-EP+ and AUGUSTUS supported by a protein database. NAR Genomics Bioinf. 3, 1–11. doi: 10.1093/nargab/lqaa108
Bu, D., Luo, H., Huo, P., Wang, Z., Zhang, S., He, Z., et al. (2021). KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 49, 317–325. doi: 10.1093/nar/gkab447
Cantalapiedra, C. P., Hernández-Plaza, A., Letunic, I., Bork, P., Huerta-Cepas, J. (2021). eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829. doi: 10.1093/molbev/msab293
Chen, S., Chen, Y., Gu, J. (2018). Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 34, 884–890. doi: 10.1093/bioinformatics/bty560
Chen, Z. H., Wang, X. L., Li, S., Yao, J. T., Shao, Z. R., Duan, D. L. (2019). Verification of the Saccharina japonica Translocon Tic20 and its Localization in the Chloroplast Membrane in Diatoms. Int. J. Mol. Sci. 20, 4000. doi: 10.3390/ijms20164000
Chi, S., Liu, T., Wang, X. M., Wang, R., Wang, S. S., Wang, G. L., et al. (2018). Functional genomics analysis reveals the biosynthesis pathways of important cellular components (alginate and fucoidan) of Saccharina. Curr. Genet. 64, 259–273. doi: 10.1007/s00294-017-0733-4
Ci, F. F., Jiang, H., Zhang, Z. H., Mao, X. Z. (2021). Properties and potential applications of mannuronan C5-epimerase: A biotechnological tool for modifying alginate. Int. J. Mol. Sci. 168, 663–675. doi: 10.1016/j.ijbiomac.2020.11.123
Dauvillée, D., Stampacchia, O., Girard-Bascou, J., Rochaix, J. D. (2003). Tab2 is a novel conserved RNA binding protein required for translation of the chloroplast psaB mRNA. EMBO J. 22, 6378–6388. doi: 10.1093/emboj/cdg591
Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635
Feinstein, Y., Borrell, V., Garcia, C., Burstyn-Cohen, T., Tzarfaty, V., Frumkin, A., et al. (1999). F-spondin and mindin: Two structurally and functionally related genes expressed in the hippocampus that promote outgrowth of embryonic hippocampal neurons. Development 126, 3637–3648. doi: 10.1242/dev.126.16.3637
Gorelova, V., Ambach, L., Rébeillé, F., Stove, C., van der Straeten, D. (2017). Folates in plants: research advances and progress in crop biofortification. Front. Chem. 5. doi: 10.3389/fchem.2017.00021
Hart, A. J., Ginzburg, S., Xu, M., Fisher, C. R., Rahmatpour, N., Mitton, J. B., et al. (2020). EnTAP: Bringing faster and smarter functional annotation to non-model eukaryotic transcriptomes. Mol. Ecol. Res. 20, 591–604. doi: 10.1111/1755-0998.13106
Haug, A., Larsen, B., Smidsrod, O. (1974). Uronic acid sequence in alginate from different sources. Carbohydr. Res. 32, 217–225. doi: 10.1016/S0008-6215(00)82100-X
Heinrich, S., Valentin, K., Frickenhaus, S., John, U., Wiencke, C. (2012). Transcriptomic analysis of acclimation to temperature and light stress in Saccharina latissima (Phaeophyceae). PLoS One 7 (8), e44342. doi: 10.1371/journal.pone.0044342
Heyraud, A., Gey, C., Leonard, C., Rochas, C., Girond, S., Kloareg, B. (1996). NMR spectroscopy analysis of oligoguluronates and oligomannuronates prepared by acid or enzymatic hydrolysis of homopolymeric blocks of alginic acid. Application to the determination of the substrate specificity of Haliotis tuberculata alginate lyase. Carbohydr. Res. 289, 11–23. doi: 10.1016/0008-6215(96)00060-2
Higashijima, S., Nose, A., Eguchi, G., Hotta, Y., Okamoto, H. (1997). Mindin/F-spondin family: Novel ECM proteins expressed in the zebrafish embryonic axis. Dev. Biol. 192, 211–227. doi: 10.1006/dbio.1997.8760
Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, 309–314. doi: 10.1093/nar/gky1085
Hynes, R. O. (2009). The extracellular matrix: Not just pretty fibrils. Science 326, 1216–1219. doi: 10.1126/science.1176009
Jiadkong, K., Fauzia, A. N., Yamaguchi, N., Ueda, A. (2024). Exogenous riboflavin (vitamin B2) application enhances salinity tolerance through the activation of its biosynthesis in rice seedlings under salinity stress. Plant Sci. 339, 111929. doi: 10.1016/j.plantsci.2023.111929
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: Genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kanehisa, M., Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27
Kim, S. H., Turnbull, J., Guimond, S. (2011). Extracellular matrix and cell signalling: the dynamic cooperation of integrin, proteoglycan and growth factor receptor. J. Endocrinol. 209, 139–151. doi: 10.1530/JOE-10-0377
Kloareg, B., Badis, Y., Cock, J. M., Michel, G. (2021). Role and Evolution of the Extracellular matrix in the acquisition of complex multicellularity in eukaryotes: A macroalgal perspective. Genes 12, 1059. doi: 10.3390/genes12071059
Kriventseva, E. V., Kuznetsov, D., Tegenfeldt, F., Manni, M., Dias, R., Simão, F. A., et al. (2019). OrthoDB v10: Sampling the diversity of animal, plant, fungal, protist, bacterial and viral genomes for evolutionary and functional annotations of orthologs. Nucleic Acids Res. 47, 807–811. doi: 10.1093/nar/gky1053
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
Lemesheva, V., Birkemeyer, C., Garbary, D., Tarakhovskaya, E. (2020). Vanadium-dependent haloperoxidase activity and phlorotannin incorporation into the cell wall during early embryogenesis of Fucus vesiculosus (Phaeophyceae). Eur. J. Phycol. 55, 275–284. doi: 10.1080/09670262.2019.1709131
Levin, D. E. (2011). Regulation of cell wall biogenesis in Saccharomyces cerevisiae: the cell wall integrity signaling pathway. Genetics 189, 1145–1175. doi: 10.1534/genetics.111.128264
Li, B., Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf. 323, 1471–2105. doi: 10.1186/1471-2105-12-323
Liu, T., Wang, X., Wang, G., Jia, S., Liu, G., Shan, G., et al. (2019). Evolution of complex thallus alga: Genome sequencing of saccharina japonica. Front. Genet. 10, 378. doi: 10.3389/fgene.2019.00378
Liu, Y., Liang, Z. R., Zhang, P. Y., Yuan, Y. M., Wu, Y. K., Zhang, D., et al. (2023). Sorus developmental biology of hybrid cultivar in Saccharina japonica: Environmental and endogenous regulation. Aquaculture 565, 739165. doi: 10.1016/j.aquaculture.2022.739165
Marino, D., Dunand, C., Puppo, A., Pauly, N. (2012). A burst of plant NADPH oxidases. Trends Plant Sci. 17, 9–15. doi: 10.1016/j.tplants.2011.10.001
McKee, J. W. A., Kavalieris, L., Brasch, D. J., Brown, M. T., Melton, L. D. (1992). Alginate content and composition of Macrocystis pyrifera from New Zealand. J. Appl. Phycol. 4, 357–369. doi: 10.1007/BF02185794
Mizuta, H., Yasui, H. (2010). Significance of radical oxygen production in sorus development and zoospore germination in Saccharina japonica (Phaeophyceae). Bot. Mar. 53, 409–416. doi: 10.1515/bot.2010.047
Mizuta, H., Yasui, H. (2012). Protective function of silicon deposition in Saccharina japonica sporophytes (Phaeophyceae). J. Appl. Phycol. 24, 1177–1182. doi: 10.1007/s10811-011-9750-8
Neeli-Venkata, R., Diaz, C. M., Celador, R., Sanchez, Y., Minc, N. (2021). Detection of surface forces by the cell-wall mechanosensor Wsc1 in yeast. Dev. Cell. 56, 2856–2870.e7. doi: 10.1016/j.devcel.2021.09.024
Nimura, K., Mizuta, H. (2001). Differences in photosynthesis and nucleic acid content between sterile and fertile parts of the sporophyte of Laminaria japonica (Phaeophyceae). Algae. 16, 151–155.
Nimura, K., Mizuta, H. (2002). Inducible effects of abscisic acid on sporophyte discs from Laminaria japonica Areschoug (Laminariales, Phaeophyceae). J. Appl. Phycol. 14, 159–163. doi: 10.1023/A:1019929712633
Nyvall, P., Corre, E., Boisset, C., Barbeyron, T., Rousvoal, S., Scornet, D., et al. (2003). Characterization of mannuronan C-5-epimerase genes from the brown alga Laminaria digitata. Plant Physiol. 133, 726–735. doi: 10.1104/pp.103.025981
Provasoli, L. (1968). “Media and prospects for the cultivation of marine algae,” in Culture and collections of algae, Proc U S-Japan Conf, Hakone, September 1966. Eds. Watanabe, A., Hattori, A.. Jpn. Soc. Plant. Physiol. (Tokyo), 63–75.
Raven, J. A., Griffiths, H. (2015). Photosynthesis in reproductive structures: costs and benefits. J. Exp. Bot. 66, 1699–1705. doi: 10.1093/jxb/erv009
Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616
Rochaix, J. D., Bassi, R. (2019). LHC-like proteins involved in stress responses and biogenesis/repair of the photosynthetic apparatus. Biochem. J. 476, 581–593. doi: 10.1042/BCJ20180718
Ros, R., Muñoz-Bertomeu, J., Krueger, S. (2014). Serine in plants: biosynthesis, metabolism, and functions. Trends Plant Sci. 19, 564–569. doi: 10.1016/j.tplants.2014.06.003
Sagi, M., Fluhr, R. (2006). Production of reactive oxygen species by plant NADPH oxidases. Plant Physiol. 141, 336–340. doi: 10.1104/pp.106.078089
Tan, K., Lawler, J. (2011). The structure of the Ca2+-binding, glycosylated F-spondin domain of F-spondin—A C2-domain variant in an extracellular matrix protein. BMC Struct. Biol. 11, 22. doi: 10.1186/1472-6807-11-22
Tarakhovskaya, E. R., Bilova, T. E., Maslov, Y. I. (2015). Hydrogen peroxide content and vanadium-dependent haloperoxidase activity in thalli of six species of Fucales (Phaeophyceae). Phycologia 54, 417–424. doi: 10.2216/15-35.1
Tseng, C. K. (2001). Algal biotechnology industries and research activities in China. J. Appl. Phycol. 13, 375–380. doi: 10.1023/A:1017972812576
Uji, T., Mizuta, H. (2022). The role of plant hormones on the reproductive success of red and brown algae. Front. Plant Sci. 13, 1019334. doi: 10.3389/fpls.2022.1019334
Uji, T., Ueda, S., Mizuta, H. (2022). Identification, characterization, and expression analysis of spondin-like and fasciclin-like genes in neopyropia yezoensis, A marine red alga. Phycology 2, 45–59. doi: 10.3390/phycology2010003
Umemiya, T., Takeichi, M., Nose, A. (1997). M-spondin, a novel ECM protein highly homologous to vertebrate F-spondin, is localized at the muscle attachment sites in the Drosophila embryo. Dev. Biol. 186, 165–176. doi: 10.1006/dbio.1997.8591
Vellosillo, T., Martínez, M., López, M. A., Vicente, J., Cascón, T., Dolan, L., et al. (2007). Oxylipins produced by the 9-lipoxygenase pathway in Arabidopsis regulate lateral root development and defense responses through a specific signaling cascade. Plant Cell. 19, 831–846. doi: 10.1105/tpc.106.046052
Verna, J., Lodder, A., Lee, K., Vagts, A., Ballester, R. (1997). A family of genes required for maintenance of cell wall integrity and for the stress response in Saccharomyces cerevisiae. Proc. Natl. Acad. Sci. U.S.A. 94, 13804–13809. doi: 10.1073/pnas.94.25.13804
Xu, D., Wang, D. S., Li, B., Fan, X., Zhang, X. W., Ye, N. H., et al. (2015). Effects of CO2 and seawater acidification on the early stages of Saccharina japonica development. Environ. Sci. Technol. 49, 3548–3556. doi: 10.1021/es5058924
Yotsukura, N., Seki, S., Sasaki, S., Yoshida, M. (2022). Planned seedling production in the distribution area of Saccharina japonica: Sorus formation through the induction of sporophyte maturation and the culture of the seedlings produced. Aquacul. Res. 53, 676–683. doi: 10.1111/are.15612
Zdobnov, E. M., Kuznetsov, D., Tegenfeldt, F., Manni, M., Berkeley, M., Kriventseva, E. V. (2021). OrthoDB in 2020: Evolutionary and functional annotations of orthologs. Nucleic Acids Res. 49, 389–393. doi: 10.1093/nar/gkaa1009
Keywords: brown algae, gene expression, reproduction, Saccharina japonica, sorus
Citation: Uji T, Kandori T and Mizuta H (2024) Identification of differential gene expression related to reproduction in the sporophytes of Saccharina japonica. Front. Plant Sci. 15:1417582. doi: 10.3389/fpls.2024.1417582
Received: 15 April 2024; Accepted: 12 July 2024;
Published: 06 August 2024.
Edited by:
Anna Torelli, University of Parma, ItalyReviewed by:
Rafael R. Robaina, University of Las Palmas de Gran Canaria, SpainVojtech Zarsky, University of Ostrava, Czechia
Copyright © 2024 Uji, Kandori and Mizuta. 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) and the copyright owner(s) 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: Toshiki Uji, t-uji@fish.hokudai.ac.jp