- Research article
- Open access
- Published:
Comparative genomics of four strains of the edible brown alga, Cladosiphon okamuranus
BMC Genomics volume 21, Article number: 422 (2020)
Abstract
Background
The brown alga, Cladosiphon okamuranus (Okinawa mozuku), is one of the most important edible seaweeds, and it is cultivated for market primarily in Okinawa, Japan. Four strains, denominated S, K, O, and C, with distinctively different morphologies, have been cultivated commercially since the early 2000s. We previously reported a draft genome of the S-strain. To facilitate studies of seaweed biology for future aquaculture, we here decoded and analyzed genomes of the other three strains (K, O, and C).
Results
Here we improved the genome of the S-strain (ver. 2, 130 Mbp, 12,999 genes), and decoded the K-strain (135 Mbp, 12,511 genes), the O-strain (140 Mbp, 12,548 genes), and the C-strain (143 Mbp, 12,182 genes). Molecular phylogenies, using mitochondrial and nuclear genes, showed that the S-strain diverged first, followed by the K-strain, and most recently the C- and O-strains. Comparisons of genome architecture among the four strains document the frequent occurrence of inversions. In addition to gene acquisitions and losses, the S-, K-, O-, and C-strains possess 457, 344, 367, and 262 gene families unique to each strain, respectively. Comprehensive Blast searches showed that most genes have no sequence similarity to any entries in the non-redundant protein sequence database, although GO annotation suggested that they likely function in relation to molecular and biological processes and cellular components.
Conclusions
Our study compares the genomes of four strains of C. okamuranus and examines their phylogenetic relationships. Due to global environmental changes, including temperature increases, acidification, and pollution, brown algal aquaculture is facing critical challenges. Genomic and phylogenetic information reported by the present research provides useful tools for isolation of novel strains.
Background
Brown algae are not only significant primary producers of marine ecosystems, but also have been used as a food resource since ancient times. Recently, they have been cultivated commercially for this purpose. In Japan, the majority of edible brown algae (the class Phaeophyceae) include members of the order Laminariales, Saccharina japonica (“kombu” in Japanese) and Undaria pinnatifida (“wakame”), and the order Chordariales, Cladosiphon okamuranus (“Okinawa mozuku”) and Nemacystus decipiens (“ito-mozuku”). Especially in Okinawa (the southwestern prefecture of Japan), C. okamuranus and N. decipiens have been farmed since the 1980s and 1990s, respectively. Approximately 17,000 and 800 tons of these two species were produced in fiscal year 2017. In addition, C. okamuranus and N. decipiens are sources of fucoidan [1], a sulfated polysaccharide that has anti-coagulant, anti-thrombin-like, and tumor-suppressant activities [2].
For the last three or four decades, global environmental conditions have changed drastically, mainly due to human activities. Greenhouse gas emissions are warming and acidifying the oceans. Pollution from agriculture and sewage is degrading seawater quality, and such problems pose a greater threat to coral reefs in transparent seawater. Aquaculture of “Okinawa mozuku” and “ito-mozuku” in Okinawa has been carried out along the coast, close to coral reefs. Accordingly, brown-alga aquaculture has also been threatened in recent years.
Given the environmental threats facing phaeophyte aquaculture, it is essential to identify and maintain strains with different physiological features. As mentioned above, culturing of C. okamuranus commenced in the 1980s and the predominant strain is the S-strain (‘Shikenjo-kabu’, registered as Inou-no-Megumi in Japanese). C. okamuranus has frond-like sporophytes, the main axes of which are 2 ~ 5 mm in diameter (Fig. 1a). The S-strain exhibits comparatively long lateral branches, and the body is not tough or fibrous.
Given that algae share little similarity with other organisms, the limited amount of genomic data available makes it difficult to determine functions of unknown algal genes [3]. Although nucleic acid extraction from algae is complicated because of the abundance of polysaccharides in cell walls and in the extracellular matrix, we decoded a draft genome of the S-strain in 2016 [4]. During the past 35 years, three other strains have been cultivated, including the K-strain, originally from the Katsuren coast, the O-strain from the Onna coast, and the C-strain from the Chinen coast (Fig. 1b; Supplementary Figure S1). The K-strain comprises thicker, tougher lateral branches (Fig. 1a). The O-strain is composed of smaller, denser lateral branches (Fig. 1a), and the C-strain is intermediate in size, with thinner lateral branches (Fig. 1a). Genetic characterization of the strains is essential for future aquaculture. To help support this industry, we decoded draft genomes of the S- (ver. 2), K-, O- and C-strains.
Results
Genome constituents of four strains
The draft genome of Cladosiphon okamuranus (S-strain, ver. 1) has been reported previously [4]. The approximately 140-Mbp genome was estimated to contain 13,460 protein-coding genes. The S-strain genome assembly and annotation was improved bioinformatically in this study to approximately 130 Mbp with 12,999 predicted genes (S-strain ver. 2). Nearly 93% of gene models were supported by corresponding mRNAs. Repetitive sequences comprised 11.2% of the genome (Table 1). Quality scores of this newly assembled genome are comparable to those of Nemacystus decipiens [5], Saccharina japonica [6] and Ectocarpus siliculosus [7, 8], although it is difficult to compare these assemblies directly due to the different methods used for each species.
Using Illumina platforms, we sequenced and assembled draft genomes of the K-, O-, and C-strains (Supplementary Table S1), which are summarized in Table 1 and Supplementary Figure S2. Scaffold N50s of the four genomes ranged from 416 ~ 1051 kb (Table 1). The GC content of the four strains was identical, amounting to 54% of the DNA sequences (Table 1). The genome size of the K-strain was 135 Mbp, with an estimated 12,511 protein-coding genes. That of the O-strain was 140 Mbp with an estimated 12,548 genes, and that of the C-strain was 143 Mbp with 12,182 estimated genes (Table 1). CEGMA [9] completeness and partial scores were 86 and 94% (Table 1), indicating that genome assemblies of the four strains were adequate for comparative analyses of their genomic and gene constituents.
Repetitive sequences were estimated to constitute 11.2, 9.9, 11.5, and 12.6% of the S-, K-, O-, and C-strain genomes, respectively. This suggests that differences in genome size are not always associated with the number of repetitive sequences, since the smaller S-strain genome (130 Mbp) contained 11.2% repetitive sequences in contrast to the mid-sized K-strain genome (135 Mbp) genome with 9.9% repetitive sequences (Table 1).
DNA transposons comprised 0.4–0.6% of total genome sequences, and RNA transposons comprised an additional 0.6–0.8% (Table 2). Simple repeats and unclassified repeats accounted for 2.2–2.5% and 4.7–5.6% of the genome sequences, respectively. Since the total number of repetitive sequences was larger in the O- and C-strains than in the K-strain, the proportion of each class of repetitive sequences was also larger. No increase or decrease of repetitive sequences was specific to a given strain.
Molecular phylogeny of the four strains
All four strains, which were isolated in the early 2000s, have been maintained continuously at the Okinawa Prefectural Fisheries Research and Extension Center (OPFREC). Although the four strains each show unique morphology, the origins of these differences and phylogenetic relationships of the strains are unclear. In order to examine their evolutionary history, we investigated molecular phylogeny using sequences of 32 mitochondrial genes of 41 brown algae and randomly selected 200 single-copy orthologous nuclear genes of six brown algae.
As was evident in the resulting tree, Nemacystus and Cladosiphon form a distinct clade, corresponding to the order Chordariales (Fig. 2a, b). In addition, Cladosiphon okamuranus constitutes its own distinct clade. The S-strain diverged first, followed by the K-strain, and finally the C- and O-strains. All nodes within the order Chordariales were supported by 100% bootstrap values. The branch length or divergence time between the S-strain and the K/C/O-strains is longer than that between the K-strain and the C/O-strains. The branch length between the C- and O-strains was very short. This suggests that the S-strain is the likely ancestor of the four strains, and that the K-strain probably blanched out first from the K/C/O-strains ancestor. The C- and O-strains are likely the most recently developed.
Conservation of genome architecture among the four strains
The results of AliTV [10] (Fig. 3a) and Dot-plot analysis using D-Genies [11] (Fig. 3b-f) suggested that the four strains of C. okamuranus exhibit similarities exceeding 90% for each comparison. Genome-wide sequence similarity was not so high between the C. okamuranus S-strain and N. decipiens (compare the upper lane with next lane), and it was low between the C. okamuranus C-strain and E. siliculosus.
A similar profile of genome-wide sequence resemblance was evident in the Dot-plot analysis (Fig. 3b-f). Linearity of dot plots was evident between the S- and K-strain (Fig. 3c), between the K- and O-strain (Fig. 3d), and between the C- and O-strain (Fig. 3e). On the other hand, a weaker linear correlation was evident in sequence comparisons between the S-strain and N. decipiens (Fig. 3b), and almost no relationship exists between C-strain and E. siliculosus (Fig. 3f). The overall similarity was ~ 50% between the C. okamuranus S-strain and N. decipiens (Fig. 3a, b), and less than 7% between the C. okamuranus C-strain and E. siliculosus (Fig. 3a, f). C. okamuranus belongs to the family Chordariaceae of the order Chordariales, whereas N. decipiens pertains to the family Spermatochnaceae in the same order. E. siliculosus belongs to a different order, the Ectocarpales. Differences in genome-wide sequence similarity may depend on time since divergence, during which neutral DNA-sequence changes likely occurred.
Inversions appear to have occurred during diversification of the four strains, especially between the O- and C-strains (Fig. 3a). Although results from AliTV and Dot-plot analysis cannot be directly compared, the dot-plot profiles of O- and C-strains suggested similar frequencies of inversions.
Synteny analysis of the four strain genomes
Next we examined synteny of genes in the four C. okamuranus genomes. Analyses using i-ADHoRe [12] identified 550 genomic regions showing shared synteny among the four genomes (Fig. 4). For example, scaffold #276 of the S-strain manifests a syntenic region comprising six genes (gene IDs, g8865, g8866, g8867, g8868, g8869 and g8870) (Fig. 4a). The K-, O-, and C-strains each retained scaffolds corresponding to S-#276, (K-#485, O-#136, and C-#002). In this synteny, g8865 was present only in the S-strain, suggesting a gene loss in the ancestor of the K/O/C strains after their divergence from the S-strain. This provides support for the phylogenetic gap between the S-strain and K/O/C-strains, discussed earlier.
Scaffold #228 of the S-strain offers a second example comprising 13 genes (gene IDs, g8477, g8476, g8475 g8473, g8472, g8471, g8470, g8469, g8468, g8467, g8466, g8465, and g8464) (Fig. 4b). All genes occur in exactly the same order in corresponding scaffolds K-#485 and O-#136, although two differences were evident in the C-strain. One is an inversion of g17337, corresponding to g8473 of the S-strain, moving beyond g17336 and next to 17,339 (Fig. 4b). The other is a lack of g8465 (Fig. 4b). This result provides further support for the supposition that the C-strain genome was uniquely modified after its divergence.
Analysis of orthologous gene families
The genome project documented approximately 12,500 genes in genomes of the four strains of Cladosiphon ([4], the present study), 15,156 genes in the Nemacystus genome [5] and 17,418 genes in the Ectocarpus genome [6] (Table 1). Orthologous analysis of numbers of gene families in each genome (Fig. 5a) demonstrated that 8367 gene families were shared or conserved by the six genomes (Fig. 5b). On the other hand, 4489 families were unique to Ectocarpus, 2532 to Nemacystus, and 405 to Cladosiphon, respectively. In addition, this analysis demonstrated the presence of unique families in each strain: 187 families unique to the S-strain, 210 to the K-strain, 225 to the O-strain, and 155 to the C-strain (Fig. 5b). There were many patterns depending on how families were shared by combinations of the six genomes (Fig. 5b). For example, different numbers of gene families were shared among different combinations of strains: 123 (S/O/K), 60 (S/K/C), 59 (S/O), 59 (S/K), 55 (S/C), 53 (O/K), 42 (S/O/C), 34 (O/K/C), 26 (K/C), and 18 (O/C), providing more information about diversification of the four strains.
We further compared gene families in the four Cladosiphon strains using OrthoFinder [13]. The four strains shared 9544 gene families (Fig. 5c), but the S-, K-, O-, and C-strains include 457, 344, 367, and 262 unique families each, constituting 3.5, 2.7, 2.9, and 2.2% of the orthologous genes in those strains, respectively. These gene families may be involved in the evident morphological diversification of the four strains, and may also support different physiology as well, although that remains to be explored.
It would be highly desirable to know the functions of these unique gene families; however, extensive Blast [14] searches failed to identify orthologies for most of them, meaning that their functions are novel. Therefore, we performed GO annotation analysis of these genes (Supplementary Table S2), based on “molecular function” (Fig. 6a), “biological process” (Fig. 6b), and “cellular component” (Fig. 6c). In general, the four strains displayed similar GO annotation profiles. All four contained similar numbers of genes related to “catalytic activity” and “binding” in the category “molecular function” (Fig. 6a), and to “cellular process” and “metabolic process”, subcategories of “biological process” (Fig. 6b), and to “cell and organelle”, under the heading of “cellular component” (Fig. 6c). Some unknown genes were specific to only one or two strains (Fig. 6).
Discussion
Brown algae have served as a food resource since ancient times. In Japan, especially in Okinawa, Cladosiphon okamuranus (“Okinawa mozuku”) has been commercially farmed since the 1980s, yielding approximately 15 kt per year. Four strains of C. okamuranus with different morphologies, the S-, K-, O-, and C-strains, have been maintained at the Okinawa Prefecture Fisheries Research and Extension Center since the early 2000s. However, algal aquaculture in Okinawa now faces various threats, mostly due to surface seawater temperature increases and declining seawater quality. In such circumstances, it is desirable to understand genetic properties and the evolutionary trajectories of the strains. Such genic and genomic information may help improve aquaculture methods and/or production of new strains with greater ability to withstand environmental stresses. To this end, we decoded the 130 ~ 143 Mbp nuclear genomes of the four C. okamuranus strains. The S-strain genome was first assembled in a previous study [4], but was improved in this study, and the K-, O- and C-strain genomes were assembled for the first time in this study. The quality of all four genomes is comparable to those of other algal genomes.
Molecular phylogeny using mitochondrial and nuclear DNA sequences allowed us to infer the evolutionary trajectories of the four strains (Fig. 2). The S-strain diverged first, then the K-strain, and finally the O- and C-strains. Both analyses, one based on mitochondrial DNA sequences (Fig. 2a) and the other on nuclear DNA sequences (Fig. 2b), resulted in identical tree topologies. In addition, all nodes had 100% bootstrap support, indicating that this represents the most probable history of the four strains. Only one difference was noticed between the two trees. That is, the branch that connects Nemacystus and the Cladosiphon S-strain was much shorter in the mitochondrial tree than in the nuclear tree. This suggests an occurrence of introgression of mitochondrial genome from Nemacystus to the S-strain or common ancestor of C. okamuranus. This question should be examined in future studies.
Each strain exhibits different morphology (Fig. 1a). Together with results of molecular phylogeny, we speculate that the morphology of the Cladosiphon okamuranus ancestor was likely close to that of the extant S-strain, with comparatively long lateral branches and a soft body, because the closely-related Nemacystus decipiens exhibits morphology similar to that of the S-strain. From this ancestor, the K-strain became a middle-sized sporophyte with thicker and tougher lateral branches, and the O-strain diverged further with smaller and more dense lateral branches. The C-strain is intermediate in size with thinner lateral branches (Fig. 1a), respectively. This information appears useful for further isolation of C. okamuranus strains, which is ongoing at OPFREC.
Genome architecture of each Cladosiphon strain varied from 130 to 143 Mbp. It is generally accepted that there is some degree of correlation between the genome size, the proportion of repetitive sequence, and the number of genes [15, 16]. In Cladosiphon however such assumption is not so simple (Table 1). Instead, each of the four genomes has likely undergone its own modification. Although long-read sequence technology might reveal the relationship between the genome size and the repetitive sequence, technical issues like the presence of lots of polysaccharide in DNA samples have impeded to introduce such methods. It remains to establish new methods, by which we make these questions clear in near future.
The most interesting finding was that the four strains contain unique gene families that equate to 2 to 3.5% for all genes (Supplementary Table S2). Such proportions are quite higher than that of seven strains of Arabidopsis thaliana (0.79%) [17]. This might suggest that the genetic distances among the four strains of Cladosiphon might reach the grade, in which the four strains are distinguished as “related species or subspecies.” Although we conducted BLAST searches of these genes, most of them showed no sequence similarity to any genes in the NCBI database; hence, their functions are unknown. GO annotations for these genes indicated that most of them are associated with biological functions, such as metabolic processes and cellular components. It is likely that these gene families are involved in development of the distinctive morphology of the four Cladosiphon strains, although details of molecular mechanisms remain to be explored. It is entirely possible that they also underlie variations in morphology. Our investigation also revealed the occurrence of many inversions in the genomes during their diversification (Figs. 3a; 5b). Although frequent inversion of genomic regions is known in land plants [18], this may be the first report showing this phenomenon in brown algal genomes.
Land plant genomes have been analyzed to support basic biology and applied use as well. Comparative genomics of the deciphered genomes show the expansion of genes that may be required for their survival [19]. Diatoms, like brown algae, are classified as Stramenopiles, and various diatom genomes have been decoded due to their relatively smaller genome sizes. Analyses of diatom genomes suggest the occurrence of horizontal gene transfers from bacteria, and documents a lack of photosynthesis-related genes as well [20, 21]. On the other hand, comparative genomics of brown algae are still limited. Analyses have been restricted to genes related to fucoidan biosynthesis, extracelluar components, stress responses, and vanadium-dependent haloperoxidases [4, 5, 22]. In addition to taxon-level genomics, comparative genomics of strains / subspecies also contribute to our understanding of macroalgal evolution.
Conclusions
Genome analyses of four C. okamuranus strains with distinctively different morphologies revealed that they possess unique gene families, and phylogenetic analysis indicated that the S-strain is closest to the ancestral strain. C. okamuranus has been cultivated in Okinawa since the 1980s. Approximately 17,000 tons of C. okamuranus were produced in fiscal year 2017, representing more than 90% of the entire Japanese harvest. Due to global environmental changes, including temperature increases, ocean acidification, and pollution, brown algal aquaculture is facing critical challenges. Efforts to maintain and improve culture methods are required, and heat stress experiments are under way at OPFREC to identify strains that tolerate warmer ocean temperatures. The genomic information reported in this study may help to develop and characterize new strains.
Methods
Brown algal strains
Originally, four strains of Cladosiphon okamuranus were isolated from four local fisheries (Fig. 1b). The S-strain came from Iheya in 2008, the O-strain from Onna in 2008, the K-strain from Katsuren in 2009, and the C-strain was collected in Chinen in 2012 (Fig. 1b). From 2n clonal cycles, 2n zoospores are isolated from sporophytes, and maintained as 2n germlings on plates. Each strain exhibits characteristic morphology, distinguishing it from the others. Each strain has been maintained as a stock culture at OPFREC, Okinawa, Japan. Cladosiphon is cultivated at 22.5 °C with a 12-h light-dark cycle in seawater containing 0.5% KW21 (Daiichi Seimo Co., Ltd., Kumamoto, Japan). The life cycle of C. okamuranus includes both haploid (n) and diploid (2n) generations (Fig. 1c). 2n germlings mature into sporophytes that are harvested for market.
DNA extraction, genome sequencing, and assembly
For DNA extraction, 2n germlings of each strain were frozen in liquid nitrogen and stored at − 80 °C until use. They were crushed to powder with a frozen-cell-crusher, Cryo-Press (Microtec Co., Ltd., Chiba, Japan) and genomic DNA was extracted from the powder using a DNA extraction kit, DNA-Suisui-VS (Rizo Co., Ltd., Ibaraki, Japan). The strains have been maintained as protonemas without contamination from other eukaryotes, allowing us to extract strain-specific, genomic DNA [4].
Illumina MiSeq, HiSeq Rapid, and HiSeq 4000 platforms were used for sequencing [23]. Libraries were prepared with slight protocol modifications recommended by the manufacturer. Fragmented genomic DNA was further purified using Blue Pippin (Sage Science, Beverly, MA, USA). A paired-end library consisting of 700-bp clones was prepared for the MiSeq using a TruSeq DNA PCR-Free LT Sample Prep Kit (Illumina, San Diego, CA, USA), and mate-pair libraries with insert lengths of 2-, 3-, 4-, 5-, 6-, 7-, 8-, 10-, 12-, and 14-kbp were prepared for the HiSeq 4000 using a Nextera Mate Pair Sample Prep Kit (Illumina) (Supplementary Table S1).
K-mer counting and estimation of genome sizes were accomplished with JELLYFISH 2.2.0 [24, 25] and GenomeScope [26]. Adapter sequences were trimmed from all reads using Trimmomatic-0.30 [27]. High-quality, paired-end reads (quality > 20) were assembled de novo using Newbler 2.9 [28] and Platanus 1.2.4 [29] to create contigs. Better assembled sequences were used for downstream analysis. Subsequent scaffolding of the Newbler or Platanus output was performed using SSPACE 3.0 [30], based on Illumina mate-pair information [31]. Gaps inside scaffolds were closed using GapCloser 1.12 [32].
Diploid sequences were merged using HaploMerger2 v3.4 [33] and BLASTN (1e− 50) aligned by more than 50%. CEGMA 2.5 was used to evaluate genome assembly. Sequences that likely originated from bacteria and other microbiota were removed from the assembled genome with Maxbin version 2.2 [34] and RNAmmer 1.2 [35]. Paired-end genomic DNA reads that were not used in each of the strain-genome assembly processes were collected with kneaddata v0.6.1 (https://bitbucket.org/biobakery/kneaddata/wiki/Home). Those reads were assembled with novoPlasty (version2.7.2) [36] for the chloroplast and mitochondrial genomes of C. okamuranus.
Transcriptome analyses
RNA was isolated from sporophytes (2–5 cm) and 2n germlings (Fig. 1c). Total RNA was extracted according to manufacturer instructions, using DNase and RNeasy Plant mini kits (QIAGEN, Hilden, Germany). Transcriptome libraries were prepared using a TruSeq Stranded mRNA Library Prep kit (Illumina). RNA was sequenced as per manufacturer instructions for the Illumina HiSeq Rapid and HiSeq 4000. Only sequences of high quality (quality > 20) were assembled, using Velvet 1.2.10 [37] and Oases 0.2.08 [38].
Transposable elements and repetitive sequences
Tandem repeats were detected and classified using RepeatModeler-1.1.8 (http://www.repeatmasker.org/RepeatModeler/). A de novo repeat library was generated with RepeatScout (version 1.0.5) [39]. Transposons and SINE in the scaffold were identified using RepeatMasker (ver. 4.0.7, http://www.repeatmasker.org/RMDownload.html) with the Repbase (version 21.01) [40].
Gene model prediction
For each strain genome, a set of gene model predictions (C. okamuranus Gene Model ver. 1) was generated with AUGUSTUS 3.2.1 [41], which was trained on transcriptome contigs recommended by PASA 2.2.0 [42]. Gene models were predicted by running AUGUSTUS, by mapping RNA-seq and transcriptome data to the genome, and by using results of repetitive sequence analysis to greatly improve gene prediction accuracy. Predicted gene models were refined with PASA.
Gene functional annotation
In order to identify putative C. okamuranus orthologous genes, reciprocal BLAST analysis was performed. This was carried out using mutual best hits with genes from N. decipiens, Ectocarpus siliculosus, and the non-redundant protein sequences database from NCBI against C. okamuranus gene models (BLASTP) or their assembly (TBLASTN). A second approach used for encoded proteins with one or more specific protein domains was to screen the models using HMMER (hmmer3) [43] against the Pfam database (Pfam-A.hmm, release 24.0, http://pfam.sanger.ac.uk) [44], which contains approximately 11,000 conserved domains. Encoded proteins were also analyzed using InterProScan 5.25–64.0 [45] and the SwissProt database for gene ontology and functional annotation.
Identification of orthologous gene groups and synteny analysis
Protein sequences of C. okamuranus, N. decipiens, and E. siliculosus were analyzed with OrthoFinder version 2.0.0, using default parameters to identify orthologous gene groups. Synteny of six brown algal genomes was analyzed with i-ADHoRe 3.0 using default parameters.
Gene collection for phylogenetic tree analysis
Sets of related sequences were subjected to phylogenetic analyses to more precisely determine orthologous relationships among C. okamuranus strains. Mitochondrial genome sequences of 37 brown algae were downloaded from the NCBI database or our genome browsers. Mitochondrial genomes were annotated using GeSeq. cDNA sequences of Atp6, Atp8, Atp9, Cox1, Cox3, Cob, Nad1, Nad2, Nad3, Nad4, Nad4l, Nad5, Nad6, Nad7, Nad9, Rpl2, Rpl5, Rpl14, Rpl16, Rpl31, Rps2, Rps3, Rps4, Rps7, Rps8, Rps10, Rps11, Rps12, Rps13, Rps14, Rps19, and Tatc genes from brown algae were collected. 200 randomly selected, single-copy, orthologous nuclear genes of six brown algae were collected using OrthoFinder output (Supplementary Table S3). 32 mitochondrial and 200 nuclear gene sequences were independently aligned using MAFFT [46] with default options. Spurious sequences or poorly aligned regions were filtered using trimAl [47]. Then filtered sequences were concatenated. Phylogenetic trees were constructed using the maximum likelihood method (GTR-gamma model) with RAxML version 8.2.11 [48] with partition analysis excluding third codon positions and 1000 bootstrap replicates.
Genome-wide comparative analysis among six brown algal genomes
AliTV was used to examine the overall similarity of six genomes, including four strains of C. okamuranus, with N. decipiens, and E. siliculosus as outgroups. Using D-Genies, Dot-plot analysis was also used to compare sequence similarities between them (15 pairwise comparisons of the six genomes, using a 10-kb sliding window).
Availability of data and materials
Sequence data were deposited in DDBJ as BioProject ID, PRJDB9346. A genome browser has been established for the assembled genome sequences using the JavaScript-based Genome Browser (JBrowse) 1.11.6 [49]. Assembled sequence and gene models are accessible at http://marinegenomics.oist.jp/gallery/.
Abbreviations
- GO:
-
Gene Ontology
- C.o. :
-
Cladosiphon okamuranus
- N.d. :
-
Nemacystus decipiens
- E.s. :
-
Ectocarpus siliculosus
- OPFREC:
-
Okinawa Prefectural Fisheries Research and Extension Center
References
Tako M, Nakada T, Hongou F. Chemical characterization of Fucoidan from commercially cultured Nemacystus decipiens (Itomozuku). Biosci Biotechnol Biochem. 1999;63(10):1813–5.
Baba M, Snoeck R, Pauwels R, de Clercq E. Sulfated polysaccharides are potent and selective inhibitors of various enveloped viruses, including herpes simplex virus, cytomegalovirus, vesicular stomatitis virus, and human immunodeficiency virus. Antimicrob Agents Chemother. 1988;32(11):1742–5.
Chan CX, Ho CL, Phang SM. Trends in seaweed research. Trends Plant Sci. 2006;11(4):165–6.
Nishitsuji K, Arimoto A, Iwai K, Sudo Y, Hisata K, Fujie M, Arakaki N, Kushiro T, Konishi T, Shinzato C, et al. A draft genome of the brown alga, Cladosiphon okamuranus, S-strain: a platform for future studies of ‘mozuku’ biology. DNA Res. 2016;23(6):561–70.
Nishitsuji K, Arimoto A, Higa Y, Mekaru M, Kawamitsu M, Satoh N, Shoguchi E. Draft genome of the brown alga, Nemacystus decipiens, Onna-1 strain: fusion of genes involved in the sulfated fucan biosynthesis pathway. Sci Rep. 2019;9(1):4607.
Cormier A, Avia K, Sterck L, Derrien T, Wucher V, Andres G, Monsoor M, Godfroy O, Lipinska A, Perrineau MM, et al. Re-annotation, improved large-scale assembly and establishment of a catalogue of noncoding loci for the genome of the model brown alga Ectocarpus. New Phytol. 2017;214(1):219–32.
Ye N, Zhang X, Miao M, Fan X, Zheng Y, Xu D, Wang J, Zhou L, Wang D, Gao Y, et al. Saccharina genomes provide novel insight into kelp biology. Nat Commun. 2015;6:6986.
Cock JM, Sterck L, Rouze P, Scornet D, Allen AE, Amoutzias G, Anthouard V, Artiguenave F, Aury JM, Badger JH, et al. The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010;465(7298):617–21.
Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23(9):1061–7.
Ankenbrand MJ, Hohlfeld S, Hackl T, Forster F. AliTV-interactive visualization of whole genome comparisons. Peerj Comput Sci. 2017;3:e116.
Cabanettes F, Klopp C. D-GENIES: dot plot large genomes in an interactive, efficient and simple way. PeerJ. 2018;6:e4958.
Proost S, Fostier J, De Witte D, Dhoedt B, Demeester P, Van de Peer Y, Vandepoele K. i-ADHoRe 3.0--fast and sensitive detection of genomic homology in extremely large data sets. Nucleic Acids Res. 2012;40(2):e11.
Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16:157.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421.
Hou Y, Lin S. Distinct gene number-genome size relationships for eukaryotes and non-eukaryotes: gene content estimation for dinoflagellate genomes. PLoS One. 2009;4(9):e6978.
Feng R, Wang X, Tao M, Du G, Wang Q. Genome size and identification of abundant repetitive sequences in Vallisneria spinulosa. PeerJ. 2017;5:e3982.
Jiao WB, Schneeberger K. Chromosome-level assemblies of multiple Arabidopsis genomes reveal hotspots of rearrangements with altered evolutionary dynamics. Nat Commun. 2020;11(1):989.
Fransz P, Linc G, Lee CR, Aflitos SA, Lasky JR, Toomajian C, Ali H, Peters J, van Dam P, Ji X, et al. Molecular, genetic and evolutionary analysis of a paracentric inversion in Arabidopsis thaliana. Plant J. 2016;88(2):159–78.
Zhang J, Fu XX, Li RQ, Zhao X, Liu Y, Li MH, Zwaenepoel A, Ma H, Goffinet B, Guan YL, et al. The hornwort genome and early land plant evolution. Nat Plants. 2020;6(2):107–18.
Bowler C, Allen AE, Badger JH, Grimwood J, Jabbari K, Kuo A, Maheswari U, Martens C, Maumus F, Otillar RP, et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature. 2008;456(7219):239–44.
Nakayama T, Kamikawa R, Tanifuji G, Kashiyama Y, Ohkouchi N, Archibald JM, Inagaki Y. Complete genome of a nonphotosynthetic cyanobacterium in a diatom reveals recent adaptations to an intracellular lifestyle. Proc Natl Acad Sci U S A. 2014;111(31):11407–12.
Liu T, Wang X, Wang G, Jia S, Liu G, Shan G, Chi S, Zhang J, Yu Y, Xue T, et al. Evolution of complex Thallus alga: genome sequencing of Saccharina japonica. Front Genet. 2019;10:378.
Bentley DR. Whole-genome re-sequencing. Curr Opin Genet Dev. 2006;16(6):545–52.
Marcais G, Kingsford C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011;27(6):764–70.
Hirakawa H, Shirasawa K, Kosugi S, Tashiro K, Nakayama S, Yamada M, Kohara M, Watanabe A, Kishida Y, Fujishiro T, et al. Dissection of the Octoploid strawberry genome by deep sequencing of the genomes of Fragaria species. DNA Res. 2014;21(2):169–81.
Vurture GW, Sedlazeck FJ, Nattestad M, Underwood CJ, Fang H, Gurtowski J, Schatz MC. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics. 2017;33(14):2202–4.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Silva GG, Dutilh BE, Matthews TD, Elkins K, Schmieder R, Dinsdale EA, Edwards RA. Combining de novo and reference-guided assembly with scaffold_builder. Source Code Biol Med. 2013;8(1):23.
Kajitani R, Toshimoto K, Noguchi H, Toyoda A, Ogura Y, Okuno M, Yabana M, Harada M, Nagayasu E, Maruyama H, et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 2014;24(8):1384–95.
Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27(4):578–9.
Shoguchi E, Beedessee G, Tada I, Hisata K, Kawashima T, Takeuchi T, Arakaki N, Fujie M, Koyanagi R, Roy MC, et al. Two divergent Symbiodinium genomes reveal conservation of a gene cluster for sunscreen biosynthesis and recently lost genes. BMC Genomics. 2018;19(1):458.
Li R, Fan W, Tian G, Zhu H, He L, Cai J, Huang Q, Cai Q, Li B, Bai Y, et al. The sequence and de novo assembly of the giant panda genome. Nature. 2010;463(7279):311–7.
Huang S, Kang M, Xu A. HaploMerger2: rebuilding both haploid sub-assemblies from high-heterozygosity diploid genome assembly. Bioinformatics. 2017;33(16):2577–9.
Wu YW, Tang YH, Tringe SG, Simmons BA, Singer SW. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome. 2014;2:26.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.
Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017;45(4):e18.
Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18(5):821–9.
Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012;28(8):1086–92.
Price AL, Jones NC, Pevzner PA. De novo identification of repeat families in large genomes. Bioinformatics. 2005;21(Suppl 1):i351–8.
Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005;110(1–4):462–7.
Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008;24(5):637–44.
Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK Jr, Hannick LI, Maiti R, Ronning CM, Rusch DB, Town CD, et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003;31(19):5654–66.
Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14(9):755–63.
Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, et al. Pfam: clans, web tools and services. Nucleic Acids Res. 2006;34(Database issue):D247–51.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.
Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30(14):3059–66.
Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25(15):1972–3.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
Skinner ME, Uzilov AV, Stein LD, Mungall CJ, Holmes IH. JBrowse: a next-generation genome browser. Genome Res. 2009;19(9):1630–8.
Acknowledgements
We thank Ms. Haruhi Narisoko and Ms. Wakana Chinen for culturing Cladosiphon okamuranus, and Dr. Steven D. Aird for editing the manuscript.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
KN conceived the project, performed the wet experiments, analyzed the data, and co-authored the manuscript. AS and KH performed computational analyses. YY cultured Cladosiphon germlings. MF and MK created mate-pair libraries and performed sequencing. ES assisted with wet experiments. NS co-conducted the project and co-authored the manuscript. The author(s) read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Figure S1.
Location of Okinawa in Japan Okinawa is the southernmost prefecture of Japan. Figure S2. A summary of genome size estimates for three strains of Cladosiphon okamuranus using GenomeScope.
Additional file 2: Table S1.
A summary of DNA and RNA sequence data from three strains of Cladosiphon okamuranus.Table S2. Numbers of GO annotated genes in Cladosiphon okamuranus strain-specific gene groups. Table S3. Gene IDs of randomly selected 200 single-copy orthologous nuclear genes.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Nishitsuji, K., Arimoto, A., Yonashiro, Y. et al. Comparative genomics of four strains of the edible brown alga, Cladosiphon okamuranus. BMC Genomics 21, 422 (2020). https://doi.org/10.1186/s12864-020-06792-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-020-06792-8