Introduction
Seaweeds have been used as food and for medicines in East Asian countries for at least 1500 years (Hwang et al. 2019). The rhodophyte genus Pyropia J.Agardh (Bangiales) is a well-known farmed seaweed in East Asian countries such as Korea, Japan, and China. Pyropia is rich in nutrients including polysaccharides, proteins, steroids, carotenoids, long-chain poly-unsaturated fatty acids (PUFAs), and vitamins (Lee et al. 2016). The growing interest in seaweed-based food for the sake of health and nutrition (Van Loo et al. 2017) has increased the market demand for Pyropia, giving it a high potential economic value (Jiang et al. 2018). In addition to Pyropia tenera (Kjellman) N.Kikuchi, M.Miyata, M.S.Hwang et H.G.Choi, P. seriata (Kjellman) N.Kikuchi et M.Miyata, and P. dentata (Kjellman) N.Kikuchi et M.Miyata, P. yezoensis (Ueda) M.S.Hwang et H.G.Choi is one of the most cultivated species in Korea (Hwang et al. 2005).
Many studies have been directed an increasing the production of Pyropia, particularly through the control of algal diseases and using different cultivation technologies (Park et al. 2006,Kim et al. 2014,Kim et al. 2017). Selection, hybridization, and mutation are the basic methods used to develop seaweed strains suitable for mass cultivation, and nine Pyropia strains have been developed using these three methods in Korea from 2012 to 2018 (Park and Hwang 2014,Lee and Choi 2018,Hwang et al. 2019). Physical and chemical mutagens induce genetic variabilities in plants, and enhanced-tolerant mutants have been generated over the past 50 years (Beyaz and Yildiz 2017,Lee and Choi 2018). Several studies have also used mutagens to induce mutations in seaweed (Guimarães et al. 2003,Zayadan et al. 2014,Lee and Choi 2018). Gamma radiation is a physical mutagen generally used to obtain mutants with desirable traits. In this study, a mutant with high growth rate was developed using gamma rays. However, the mechanisms underlying its high growth rate were difficult to understand owing to the limited genetic information. Although some studies have published the draft and plastid genomes of P. yezoensis, information on its genome remains inadequate (Nakamura et al. 2013,Wang et al. 2013).
Transcriptome analysis using high-throughput RNA-sequencing (RNA-seq) offers the opportunity to investigate the molecular components de novo. Several transcriptome-based studies have been performed for non-model organisms such as seaweeds (Bryant et al. 2017,Carvalho et al. 2018,Garcia-Jimenez et al. 2018,Nan et al. 2018). The analysis of expressed transcripts may reveal unknown metabolic pathways or improve our knowledge about specific regulatory mechanisms.
In the present study, we performed de novo RNA-seq analysis to evaluate the mechanism underlying the high growth rate of the mutant P. yezoensis. Computational analysis of the transcriptome analysis result was combined with experiments to reveal the differences between the mutant and wild-type strains. This study aimed to elucidate the metabolic changes in the mutant strain with a high growth rate. The results of this study could improve our understanding about the functions of Pyropia genes and assist in the development of other mutants or recombinant seaweeds. Furthermore, the isolated mutant Pyropia with high growth will be useful in seaweed-based industries.
Materials and methods
Algal culture and development of a mutant with high growth rate
Pyropia yezoensis specimens were obtained from the Seaweed Research Center (Mokpo, Republic of Korea). Gametophytes of PyWT and Py2K were cultured in filter-sterilized Provasoli Enriched Seawater (PES) medium (Provasoli 1963) at 10 °C. A growth chamber set with a 10 h light and 14 h dark cycle (80 μmol photons·m−2·s−1) and filter-sterilized air was used. The medium was replaced every week.
Mutant P. yezoensis with a high growth rate was isolated as previously described (Lee and Choi 2019). Blades (1-2 cm long) were cultured for 4 weeks and exposed to 60Co gamma rays at 2 kGy at a rate of 10 kGy h-1 (Nordion International, Ottawa, Ontario, Canada). The samples were allowed to recover at 10 °C for 24 h in the dark. The surviving cells were selected and cultivated under similar conditions for 30 days. The length and width of blades were measured every week, and blades longer and wider than wild-type blades were isolated. Isolated blades were cultured for 6 weeks under similar conditions. The selected mutant with high growth rate was termed as “Py2K.”
RNA sequencing and de novo transcriptome analysis
The gametophytes were harvested after 4 weeks of cultivation and frozen in liquid nitrogen. Total RNA was prepared and sequenced by Insilicogen Company (Yongin, Republic of Korea) (Baek et al. 2018). The Illumina NextSeq 500 platform was used for pair-end RNA-seq.
Clean reads generated by Trimmomatic (Bolger et al. 2014) were assembled using Trinity (Haas et al. 2013) in an in silico normalization and a de novo assembly mode. The published draft genome of P. yezoensis (Nakamura et al. 2013) was used as the reference indexed with BBsplit in BBtools (Bushnell et al. 2017) to separate more reliable reads as P. yezoensis from clean reads. Putative contaminants were distinguished and eliminated from unmapped clean reads of P. yezoensis draft genome. Putative contaminants were distinguished using Diamond blast (Buchfink et al. 2015) with the NCBI non-redundant (NCBI-nr) database (Pruitt et al. 2005). The most annotated bacterial species were considered as putative contaminants. Mapped reads from P. yezoensis draft genome and unmapped reads from P. yezoensis and filter reads were merged and then de novo assembled using Trinity. The assemblies were separately generated from PyWT and Py2K.
The open-reading frame was identified by TransDecoder (Haas et al. 2013), and duplicates were eliminated using CD-HIT-EST (Li and Godzik 2006) by clustering based on at least 90% similarity. A shared sequence file between PyWT and Py2K was generated by CD-HIT-EST-2D and BBsplit. Transcriptome assembly was qualified using BUSCO v3 with Eukaryota database (Waterhouse et al. 2018).
Functional annotation and differential expressed genes (DEGs) analysis
NCBI-nr, SWISS-PROT (Consortium 2018), Pfam (Sonnhammer et al. 1997), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto 2000) were used for annotation. TransDecoder was used to retrieve amino acid sequences, while BLASTp (Camacho et al. 2009) was employed to descript amino acid sequences from SWISS-PROT. DIAMOND (Buchfink et al. 2015) was used with NCBI-nr to annotate amino acid sequences, while conserved protein domains from Pfam were annotated with hmmscan (Finn et al. 2011). Functional annotation was processed with an e-value ≤ l × 10−5. GO terms were annotated using the top BLASTp matched with sequences in SWISS-PROT as per the Trinotate workflow.
Gene abundance in sequences from PyWT and Py2K was estimated using Salmon (Patro et al. 2017). We used DEseq2 (Love et al. 2014) to identify DEGs between PyWT and Py2K at a false discovery rate (FDR) significant at ≤ 0.05 and an absolute value of log2 ratio ≥ 1. WEGO v2 was used to analyze functional interpretation (Ye et al. 2018).
Real-time quantitative polymerase chain reaction (RT-qPCR) was performed to validate the findings of transcriptome analysis using TB GreenTM Premix Ex TaqTM (Tli RNaseH Plus, TaKaRa) on an Illumina Eco Real-Time system. All reactions were conducted in triplicate. PyGADPH was used as an internal control (Kong et al. 2015). The primer sets specific for selected unigenes were designed using Primer3 (Untergasser, et al. 2012) (On-line Suppl. Tab. 1).
Differential expression level was calculated by the 2(−∆∆Ct) method (Livak and Schmittgen 2001), while relative expression was described by log2 2(−∆∆Ct). All experiments were performed in triplicate.
Antioxidant activity and polyphenol content
The 1,1-diphenyl-2-picryl-hydrazyl (DPPH) assay was performed as described byJao and Ko (2002) andNguyen et al. (2018). Harvested algal samples were powdered in liquid nitrogen using a mortar. Each ground sample (100 mg) was mixed with 80% ethanol at a ratio of 1:10 (w/v) and shaken at 100 rpm at 75 °C for 1 day. The mixture was filtered and condensed, and condensed samples were diluted with dimethyl sulfoxide (DMSO) to 20 mg mL-1. Then, 100 μL DPPH (1.5 × 10-4 M) and a 100 μL sample (500 μg mL-1) were placed in a 96-well plate, mixed, and reacted at 25 °C for 30 min. Absorbance was measured at 517 nm wavelength using a microplate reader (Molecular Devices, CA, San Jose, USA). The ferric ion-reducing antioxidant potential (FRAP) assay was performed with a ferric-reducing antioxidant power (FRAP) assay kit (Sigma-Aldrich, St. Louis, MO, USA) following the manufacturer’s protocol.
Total polyphenolic content was measured as described byMaksimović et al. (2005). Ground algal samples and 80% acetone were mixed at a ratio of 1:100 (w/v) and incubated in a water bath at 20 °C for 6 h. The samples were centrifuged, and the supernatants obtained were incubated with Folin-Ciocalteu reagent (Sigma) at a ratio of 1:2 (v/v) at room temperature for 1 h. The absorbance of the reaction solution was measured at 750 nm wavelength using a microplate reader.
Determination of chlorophyll fluorescence and photosynthetic pigment content
The photochemical efficiency of the photosystem II (PSII; Fv/Fm) was determined for PyWT and Py2K using a microscopy pulse-amplitude modulation apparatus (Walz GmbH, Effeltrich, Germany) as previously described byKim et al. (2006).
Phycoerythrin (PE) and phycocyanin (PC) levels were quantitatively analyzed as per the method described bySampath-Wiley and Neefus (2007). Algal samples were ground in liquid nitrogen using a mortar, and 100 mg of powdered samples were incubated with 800 μL of 0.1 M phosphate buffer (pH 6.8) at 4 °C for 12 h. The incubated samples were centrifuged at 15,000 ×g for 10 min at 4 °C, and the supernatants obtained were mixed with 200 μL phosphate buffer by vortexing. Their absorbance was measured using a microplate reader at 564, 618, and 730 nm (Molecular Devices).
Chlorophyll a (Chl a) was extracted by incubating 100 mg powdered samples with 1.8 mL of 98% acetone for 2 h at 4 °C in the dark. The supernatants obtained after centrifugation at 5000 ×g for 10 min were used to determine spectrophotometrically the content of Chl a, as described byJeffrey and Humphrey (1975). The mean of triplicate ± standard deviation was presented, and the Student t-test was used to determine the significance of the difference between the wild type and the mutant.
Results
De novo assembly and functional annotation of mutant P. yezoensis transcriptome
After exposure to 2 kGy gamma irradiation, a mutant of P. yezoensis with long and wide blades was isolated (Fig. 1A, B). The blades of Py2K were six times longer and three times wider than those of PyWT after 6 weeks of cultivation under similar conditions (Fig. 1C, D). Random amplified polymorphic DNA (RAPD) analysis was conducted to compare genotypic differences between PyWT and Py2K (Al-Zahim et al. 1999). Seven primer sets revealed the different patterns of amplified bands (On-line Suppl. Fig. 1), indicating that the isolated mutant was genetically distinct from the wild-type strain.
RNA was prepared from six samples, and approximately 424 million reads (151 base pairs) were generated by RNA-seq (On-line Suppl. Tab. 2). Elimination of putative contaminating sequences from the transcriptome yielded 197,780 transcripts and 167,165 genes (Tab. 1); N50 contained 672 bp and N10 contained 1,629 bp. While N50 was not long, N10 had sufficient length. The GC content was 52.42%, which was lower than the known GC content of the P. yezoensis genome. The E90 transcript set, indicative of the top 90% expressed transcripts, was 15,299 and the E90N50 was 675 bp. Thus, the assembled transcriptome contained near-complete gene sequence information for 78.2% genes in Eukaryota database (C: 78.2% [S: 43.9%, D: 34.3%], F: 16.5%, M: 5.3%, n: 303; C, complete BUSCOs; S, complete and single-copy BUSCOs; D, complete and duplicated BUSCOs; F, fragmented BUSCOs; M, missing BUSCOs; n, total BUSCO groups searched).
Genes were mostly annotated in NCBI-nr (140,721 genes; 71.65%), followed by Pfam (103,827 genes; 52.86%), GO (98,361 genes; 50.07%), SWISS-PROT (94,390 genes; 48.06%), and KEGG (82,023 genes; 41.76%). In total, 146,654 genes (74.67%) were annotated from at least one database.
DEGs from the wild-type and mutant variants and GO enrichment
We confirmed the high similarity between the heatmaps of samples from each group using Pearson’s correlation coefficients (On-line Suppl. Fig. 2). In all, 15,979 genes were retained at an FDR ≤ 0.05 and absolute value of log2 ratio ≥ 1 (Fig. 2). In Py2K, 9,113 genes were upregulated and 6,866 genes were downregulated as compared to PyWT. A total of 14,658 DEGs (91.73%) were annotated in at least one database, and most DEGs were annotated in NCBI-nr (14,494 DEGs; 90.70%), followed by Pfam (10,703 DEGs; 66.98%), GO (9,531 DEGs; 59.65%), SWISS-PROT (9,251 DEGs; 57.89%), and KEGG (8,042 DEGs; 50.33%).
GO term enrichment analysis revealed the annotation of 17 GO terms (Fig. 3). Four of these GO terms were mostly annotated in cellular process (GO:0009987, Biological Process), catalytic activity (GO:0003824, Molecular Function), cell (GO:0005623, Cellular Component), and cell part (GO:0044464, Cellular Component).
RT-qPCR was conducted to validate the results of transcriptome analysis with 10 randomly selected DEGs. Most genes showed similar expression patterns with DEGs (On-line Suppl. Fig. 3). Thus, RT-qPCR analysis confirmed the validity of the transcriptome analysis.
Upregulated transcripts in the mutant strain
The transcripts that were overexpressed in Py2K were related to antioxidant enzymes, shikimate pathway, nitrogen utilization, ribosome biogenesis, and energy generation. Antioxidant enzymes such as glutathione peroxidase, glutathione S-transferase, glutathione synthetase, peroxiredoxin, catalase, thioredoxin, and aldehyde dehydrogenase were found to be upregulated in the mutant strain (On-line Suppl. Tab. 3). The expression of some genes (chorismite synthase [cs], 3-dehydroquinate dehydratase [DHQase], shikimate kinase [sk], and shikimate dehydrogenase [sdh]), including those involved in the shikimate pathway related to phenol precursors, was upregulated in the mutant (On-line Suppl. Tab. 4, On-line Suppl. Fig. 4). Various transcripts related to the pathways utilizing nitrogen were overexpressed in Py2K (Fig. 4, On-line Suppl. Tab. 5). Moreover, the expression of various transcripts related to energy production such as the Calvin cycle, electron transport chain (ETC), and photosynthesis also increased in Py2K as follows: mitochondrial complex I, II, III, IV, V (genes encoding NADH-ubiquinone oxidoreductase, succinate dehydrogenase, CoQ-cytochrome c reductase, cytochrome c oxidase, and F1F0 ATP synthase, respectively), cytochrome c (Cyt c), isocitrate dehydrogenase (idh), and heme oxygenase (ho) (Fig. 4, On-line Suppl. Tab. 6).
Increased antioxidant capacity and content of photosynthetic pigments in the mutant strain
DPPH and FRAP assays were performed to confirm the increase in the antioxidant capacity of Py2K (Fig. 5A, B). DPPH radical-scavenging activity of Py2K was higher (78.43%) than that of PyWT (69.80%), while FRAP assay revealed a higher level of Fe3+ in Py2K (408.80 nM) than in PyWT (237.06 nM). We analyzed total phenolic content and found that the phenol content of Py2K was higher than that of PyWT (Fig. 5C).
The Fv/Tm value of Py2K was higher than that of PyWT, while Py2K contained more PE and PC than PyWT (Tab. 2). However, the content of Chl a slightly increased in Py2K over that in PyWT.
Discussion
In this study, a mutant P. yezoensis strain with high growth rate, Py2K, was obtained by gamma irradiation. Transcriptome analysis was performed de novo to elucidate differences in phenotypes between mutant Py2K and wild-type strains.
Transcriptome analysis revealed the upregulation in the expression of antioxidant enzymes in the mutant Py2K variant. Normal cellular metabolism generates reactive oxygen species (ROS) as a byproduct (Poljsak et al. 2013). High growth rate may induce the generation of excessive ROS as compared with normal growth conditions, owing to the activation of cellular metabolism. ROS could damage cells and, hence, antioxidants play an important role as scavengers of excessively produced ROS to limit ROS-induced damage. The increase in the antioxidant capacity of mutants developed with gamma rays is common and related to the mechanisms by which mutation is induced, involving generation of ROS or free radicals in cells (Beyaz and Yildiz 2017). The results of DPPH and FRAP assays show that Py2K had higher antioxidant capacity than did PyWT. The FRAP assay measures antioxidant capacity based on reducing power. Higher values in FRAP assay for mutant samples may indicate the additional components of Py2K as sources of reducing power. The content of total phenolic compounds, which act as reducing agents donating a single electron or hydrogen atom (Rabeta and Faraniza 2013), was higher in mutant Py2K than in the wild-type strain. Further, transcriptome analysis showed the activation of the biosynthesis pathway for phenol precursor, probably owing to the increase in the content of phenolics in the mutant strain. Therefore, high expression of antioxidant enzymes and increased levels of phenolic compounds could enhance the ROS-scavenging capacity of mutant Py2K.
The high growth of Py2K could be associated with the increased availability of nitrogen and energy regeneration system. Nitrogen is an essential nutrient and a limiting macro-nutrient in the ocean (Moore et al. 2013,Bristow et al. 2017). Nitrogen is assimilated in the following manner: NO3− uptake → NO3− → NO2− → NH4+ → glutamine → glutamate → amino acids (Lillo 2004). The gene nrt encodes a transporter that is located in the cell membrane and is involved in nitrate uptake. There are two forms of these genes, nrt1 and nrt2. In general, the members of the nrt1 family have low affinity for nitrate, while those of the nrt2 family exhibit high affinity for nitrate (Lillo 2004). In the transcriptome analysis, we failed to observe nrt1 as a DEG; however, nrt2 expression was upregulated in Py2K (On-line Suppl. Tab. 5). Assimilated NO3− is reduced to NH4+ by NR and NiR. The expression of nir is induced by nitrate, which is generated by NR. NH4+ is used to produce glutamine via glutamine synthetase (gs) in the chloroplast, and glutamine is converted to glutamate, which is used for amino acid synthesis. The expression of gs and the genes related to amino acid synthesis (those encoding D-amino-acid N-acetyltransferase, cystathionine gamma-lyase, and arginosuccinate synthase) was upregulated in Py2K. NO3− uptake could be enhanced following increased nrt expression, which resulted in increased synthesis of amino acids that were used as building blocks.
Nucleobases are rich stores of nitrogen in cells (Werner et al. 2013). Purine degradation results in the release of nitrogen, which is re-assimilated into amino acids. Purine ring catabolism via ureide-degrading reactions generates five NH3 molecules from one adenosine molecule. The purine ring catabolism pathway was constructed using ureide-degrading reactions with reference toWerner and Witte (2011),Hafez et al. (2017) and,Melino et al. (2018), although this pathway is controversial in plants. The genes aln and aah (encoding allantoate amidohydrolase) are important for purine ring catabolism, and Arabidopsis mutants deficient in aln and aah had accelerated senescence, rapid transition, and impairment in reproductive growth (Takagi et al. 2018) (Fig. 4). Upregulated expression of aln could contribute to the activation of purine ring catabolism, and activated nitrogen re-assimilation would allow mutant variants to grow at higher rates than the wild-type strain.
Ribosome biogenesis is an important process regulating plant growth and development (Hang et al. 2018). Repressed ribosome biogenesis results in defective vegetable growth and development (Beine-Golovchuk et al. 2018,Maekawa and Yanagisawa 2018, Palm et al. 2019). Repression of ribosome biogenesis is imperative under starvation conditions because it is one of the most energy consuming cellular processes (Kraft et al. 2008,Ohbayashi and Sugiyama 2018). However, transcripts related to ribosome biogenesis were upregulated in Py2K (On-line Suppl. Tab. 5), indicating that their high growth rate may induce the underutilization of cellular nitrogen as building blocks and mediate activation of nitrogen assimilation and re-assimilation; ribosome biogenesis replenished the cellular building blocks that were lacking.
Oxidative phosphorylation generates the bulk of ATP from mitochondria that is essential for the maintenance and growth of cells (Millar et al. 2011). In plants, the electron transport chain (ETC) comprises the mitochondrial complex I, II, III, IV, and V, ubiquinone (UQ), Cyt c, and alternative oxidase (AOX). Enhanced ETC leads to the generation of higher ATP levels, which subsequently support cell maintenance. Cyt c, a small heme protein and an essential component of the ETC in mitochondria, was found to be upregulated in Py2K (Virolainen et al. 2002,Ott et al. 2002,Reape et al. 2008,Matsuura et al. 2016). The presence of bax/bak in the outer membrane and peroxidation of cardiolipin (CL) stimulate the release of Cyt c. The expression of bax/bak is negatively regulated by the mitochondrial complex I (Kim et al. 2003,Westphal et al. 2011). Upregulated expression of mitochondrial complex I may result in the repression of the expression of bax/bak in the mutant variant. Further, CL peroxidation can be induced by ROS generated from oxidative phosphorylation to produce CL hydroperoxide (CL-OOH) and free Cyt c (Perier et al. 2005,de Paepe et al. 2014). Released Cyt c is immediately oxidized to Cyt cox, which activates caspase and combines with caspase-9 and apoptotic protease activating factor 1 (APAF1) to form a complex that promotes apoptosis (McMillin and Dowhan 2002,Vianello et al. 2007,Brown and Borutaite 2008).
Py2K showed upregulated expression of mitochondrial manganese-dependent superoxide dismutase (mnsod), which removes ROS generated from activated ETC (Vanlerberghe 2013) and prevents ROS-mediated CL peroxidation and the release of Cyt c. Cyt c reductants can reduce Cyt cox to Cyt cred (Brown and Borutaite 2008) and may include superoxide, ascorbate, reduced glutathione, cytochrome P450, cytochrome b5, NAD(P)H, and tetramethyl-p-phenylenediamine.
In general, plant mitochondria contain AOX in their ETC systems (Vanlerberghe 2013). AOX is located between mitochondrial complexes II and III, and electrons are distributed to AOX and complex III. AOX is a non-proton pumping oxidase that can reduce the production of ROS generated from complex I and III. In addition, the expression of aox gene is upregulated under abiotic/biotic stress conditions because increased AOX expression is beneficial for maintaining homeostasis and to increase resistance against stress. AOX contributes to redox homeostasis for nitrate assimilation by supplying high reducing power to reduce nitrate to NH4+ (Vanlerberghe 2013). Stress conditions such as deficient cellular nitrogen may upregulate the expression of aox and contribute to nitrate assimilation by supplying high reducing power.
The Calvin cycle provides organic acids, and the ETC oxidizes NADH. The increase in the expression of idh may result in the activation of the Calvin cycle in Py2K mutant and contribute to ETC activation. Therefore, these findings suggest the activation of oxidative phosphorylation in the high growth rate mutant. Upregulated Cyt c expression could induce apoptosis; however, the high level of cytochrome b5 in Py2K may be useful to reduce oxidized Cyt c; thus, Cyt c may not induce apoptosis. Furthermore, the increased expression of aox could provide high reducing power for nitrate assimilation and prevent excessive ROS production (Fig. 4, On-line Suppl. Tab. 6).
Photosynthesis is the process in which energy is derived from light and stored in the form of nutrients. P. yezoensis uses Chl a and phycobiliproteins to obtain light for photosynthesis. Photosynthesis efficiency determines the growth of P. yezoensis, and an increase in photosynthesis efficiency may enhance the growth of P. yezoensis (Gao et al. 1991). Although Chl a is a photosynthetic pigment in plants and algae, the high phycobiliprotein content and phycobiliprotein/Chl a ratio indicate enhanced energy efficiency and rapid growth (Zhang et al. 2014). The increased Fv/Tm value and PE + PC/Chl a ratio in Py2K suggest that the enhanced photosynthetic efficiency may contribute to the high growth rate of the mutant (Tab. 2).
In summary, DEGs analysis revealed the increase in nitrogen uptake owing to the upregulated expression of nrt and that the activation of nitrogen assimilation and re-assimilation may serve as reasons for the high growth rate of the mutant (Fig. 6). In addition, aox expression was upregulated in mitochondria and may have led to nitrogen assimilation by supplying NAD(P)H (Vanlerberghe 2013). Nitrogen availability might have increased and enhanced amino acid synthesis and ribosome biogenesis. Increased contents of amino acids as building blocks and high number of ribosomes could be useful during cell growth. Furthermore, increased nitrogen availability may regulate the synthesis of photosynthetic pigments and generate high levels of PE and PC in Py2K (Gao et al. 2019). The high expression of the genes related to the synthesis of phycobiliprotein precursor was evident in DEG analysis. High levels of PE and PC could enhance photosynthesis in the mutant, thereby increasing its growth rate (Zhang et al. 2014). Further, the differences in the levels of PE and PC between the mutant and PyWT may have led to variations in algal color (Fig. 1A, B). Increased nitrogen availability may activate cellular processes such as photosynthesis and induce the production of ROS. Some genes related to antioxidants were highly expressed in Py2K. Improved antioxidant capability may sufficiently scavenge ROS during the activation of cellular processes and prevent cell death.
In this study, transcriptome analysis was performed to understand the metabolic changes in the newly isolated mutant Pyropia Py2K. Py2K has a high growth rate and would be useful in seaweed cultivation. Furthermore, the transcriptome analysis performed could be useful for an understanding of the expression pattern of genes and the development of new cultivars of seaweeds for which genetic information is missing.