Microbial community and functions associated with digestion of algal polysaccharides in the visceral tract of Haliotis discus hannai: Insights from metagenome and metatranscriptome analysis

Abstract

Haliotis discus hannai, a species of Pacific abalone, is a highly valuable food source throughout Northeast Asia. As H. discus hannai primarily feed on brown algae and largely extract their energy from algal polysaccharides, understanding the mechanisms by which they digest algal polysaccharides is essential for elucidating their energy metabolism. Gut microbes, as well as the host animal, are involved in the process of polysaccharide degradation. To identify algal polysaccharide-digestion mechanisms and their origin, we analyzed the metagenome and metatranscriptome of abalone visceral extracts. Microbial communities were characterized using the 16S rRNA gene sequences in the metagenome and our results differed significantly from those of previous studies using traditional microbiological methods such as bacterial cultivation and cloning. A greater diversity of bacterial taxa was identified here than was previously identified using cultivation methods. Furthermore, the most abundant bacterial taxa also differed from previous studies, which is not common when comparing the results of bacterial culturing with those of molecular methodologies. Based on the metatranscriptome, overall functions were identified and additional analyses were performed on the coding sequences of algal polysaccharide-digestive enzymes, including alginate lyase. Results of the transcriptomic analyses suggest that alginate lyase in the visceral extracts of H. discus hannai was produced by the host itself, not by visceral bacteria. This is the first next-generation sequencing study performed on abalone to characterize the visceral microbiota and the source of the ability to digest algal polysaccharides by analyzing the metagenome and metatranscriptome together.

About the Authors:

Bo-Hye Nam

Contributed equally to this work with: Bo-Hye Nam, Jisung Jang

Roles Funding acquisition, Project administration, Resources

Affiliation: Biotechnology Research Division, National Institute of Fisheries Science, Busan, Republic of Korea

Jisung Jang

Contributed equally to this work with: Bo-Hye Nam, Jisung Jang

Roles Conceptualization, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing

Affiliations Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea, C&K Genomics, Seoul, Republic of Korea

Kelsey Caetano-Anolles

Roles Writing – review & editing

Affiliation: Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea

Young-Ok Kim

Roles Funding acquisition, Resources

Affiliation: Biotechnology Research Division, National Institute of Fisheries Science, Busan, Republic of Korea

Jung Youn Park

Roles Funding acquisition, Resources

Affiliation: Biotechnology Research Division, National Institute of Fisheries Science, Busan, Republic of Korea

Hawsun Sohn

Roles Funding acquisition, Resources

Affiliation: Cetacean Research Institute, National Institute of Fisheries Science, Nam-gu, Ulsan, Republic of Korea

Sook Hee Yoon

Roles Funding acquisition

Affiliation: Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea

Heebal Kim

Roles Project administration, Resources, Software, Supervision

Affiliations Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Republic of Korea, C&K Genomics, Seoul, Republic of Korea, Department of Agricultural Biotechnology and Research Institute of Agriculture and Life Sciences, Seoul National University, Seoul, Republic of Korea

Woori Kwak

Roles Conceptualization, Formal analysis, Methodology, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: C&K Genomics, Seoul, Republic of Korea

ORCID logo http://orcid.org/0000-0003-0600-067X

Introduction

Abalone are a group of marine snails which constitute a single genus, Haliotis in the Haliotidae family [1]. Within the genus Haliotis, 57 abalone species are classified principally based on their habitat [2]. One characteristic of abalone is that they consume only macroalgae, primarily red and brown algae [3]. The main food of Haliotis discus hannai is Laminaria japonica, a species of brown algae [3]. As polysaccharides are used for energy storage in brown algae, it is crucial to understand the polysaccharide digestion mechanism of abalone to understand their energy metabolism. Brown algae are composed largely of polysaccharides such as alginates, β-glucans (laminarans), cellulose, heteroglycans, and fucoidans [4, 5]. In particular, alginates constitute up to 40% of the dry mass of brown algae [6]. Therefore, digestive enzymes such as alginate lyase, β-1,3-glucanase, cellulase, and agarase are essential for abalone to obtain nutrients from these algal polysaccharides. The ability of abalone visceral extracts to digest algal polysaccharides has been investigated previously using culturing methods on agar plates. For example, 70% of bacteria isolated from the gut of Haliotis gigantea showed an ability to degrade alginate on agar plates [7]. In H. discus hannai, Vibrio halioticoli, a species of gut microbe, produced acetic acid and volatile short chain fatty acids on a growth medium that included sodium alginate due to alginate degradation by multiple enzymes [8]. Likewise, the gut bacteria of abalone can help it obtain energy by supporting the degradation of algal polysaccharides. Gut bacterial enzymes degrade polysaccharides to smaller components that are more easily absorbed by their host [9]. Bacterial culturing methodologies were useful for identifying several species and evaluating their ability to digest algal polysaccharides. However, the results of culturing methods have not accurately reflected the microbial community due to the differing oxygen conditions of the experimental environment and the anaerobic conditions within the actual gut [10]. Additionally, the artificial high-nutrient content found in the growth media used for these cultivation methods differs significantly from the conditions within the visceral environment, which can greatly affect both the relative abundance of bacteria and the dominant species present [10]. To overcome the limitations of culturing methods, we applied high-throughput sequencing to visceral extracts of H. discus hannai. This is the first next-generation sequencing (NGS) study to use 16S rRNA gene sequences to characterize the abalone visceral microbiota. By analyzing the abalone visceral metagenome, we characterized the bacterial community and bacterial abundance in detail.

Until now, both the composition and function of the bacterial community within abalone viscera were poorly understood. In this study, we examined the visceral metatranscriptome of H. discus hannai to determine the functional profile of transcripts from the visceral environment. Additionally, we focused on transcripts related to the energy-yielding pathways that allow H. discus hannai to digest brown algae. By analyzing the phylogeny of related enzymes to identify their origin, we determined the roles of H. discus hannai and its gut bacteria in extracting energy from the indigestible polysaccharides of brown algae.

Materials and methods

Sample collection

Haliotis discus hannai (8-cm shell length) samples used for metagenomic and metatranscriptomic analyses of the abalone digestive system were purchased from Wan Island, Korea in February of 2016. Only the alimentary canal, extracted through viscera dissection, was used in the analysis.

Library construction for bacterial community analysis

The visceral extracts of all samples from alimentary canal were released and pooled into Eppendorf tubes. The mixture of visceral extracts was homogenized with sterile plastic pestles. Genomic DNA was isolated from intestine of abalone using PowerSoil DNA Isoloation kit (Mo Bio). To amplify the variable V3 and V4 region of the 16S rRNA gene, we performed PCR as follows: 10 ng of DNA and 2× Quick Taq HS DyeMix Z (Toyobo) were used to make a 30-μl reaction volume to perform PCR with the following PCR conditions: 95°C for 3 min; 25 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s; followed by a final extension of 72°C for 5 min and a hold at 4°C. The primers used in this analysis were as follows:

1st PCR Forward Primer:

5' -TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3'

1st PCR Reverse Primer:

5'-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAATCC-3'

The PCR products were sequenced from an Illumina MiSeq 300-bp paired-end library, using the standard protocol. Descriptions of the data are provided in S1 Table.

Metatranscriptome library construction

The whole transcriptome sequencing (WTS) library was constructed with the TruSeq Stranded Total RNA Sample Prep Kit (Illumina). Ribosomal RNAs were excluded from the total RNA using an rRNA removal kit, followed by RNA purification with RNAClean XP beads (Beckman Coulter). The remaining RNA was fragmented with EFP Mix under PCR conditions of 94°C for 6 min and a 4°C hold. The fragmented RNAs were used to synthesize first strand cDNA using reverse transcriptase and random primers under PCR conditions of 25°C for 10 min, 42°C for 15 min, 70°C for 15 min, and a 4°C hold. For double stranded cDNA synthesis from the first strand cDNAs, Second Strand Marking Master Mix was used under PCR conditions of 16°C for 1 h and a 4°C hold. Library construction was performed using the Illumina NextSeq 500 reagents following the manufacturer’s protocol. Details are provided in S2 Table.

Metagenomic analysis of the bacterial community

The quality of generated reads was checked using FASTQC software [11] and adapter sequences were removed using Trimmomatic software [12]. Paired-end reads were merged using FLASH software [13] to generate a single read including the V3-V4 region of the 16S rRNA gene. Finally, community analysis was conducted using MG-RAST software [14] against the Greengenes database [15].

Metatranscriptome analysis

Quality control processing was conducted using FASTQC and Trimmomatic. The Trinity de novo transcriptome assembly tool was used to identify genes expressed in the visceral extract of abalone [16]. Summary statistics for the total transcriptome obtained from Trinity are shown in S3 Table. Using paired-end read remapping, we filtered out misassembled candidate transcripts displaying low paired-end read mapping from the metatranscript assembly. Candidate coding regions of the assembled transcripts were identified using Transdecoder software [17]. Functional annotation and the origins of coding region sequences were identified by cross-referencing the dataset against the nt database (non-redundant nucleotide database) using BLASTn (e-value < 1e–5, > 50% sequence similarity, and > 50% alignment coverage). Additionally, the overall functions of the metatranscriptome were profiled using the SAMSA2 pipeline [18]. At the preprocessing step, PEAR software was used to merge paired-end reads, and low-quality sequences and adaptor contamination of the merged reads was removed using Trimmomatic [19, 20]. Among the cleaned reads, ribosomal reads were removed using the SortMeRNA program with the SILVA and Rfam rRNA databases [21–23]. Next, without read assembly, a BLAST-like algorithm in the DIAMOND program was applied to annotate the remaining reads using two databases, NCBI’s RefSeq microbial genomes database and the SEED subsystems hierarchical database [24–26]. Finally, annotation results were visualized using scripts included in the SAMSA2 package [18].

Phylogenetic tree construction

Publicly available alginate lyase sequences were collected and used to construct a phylogenetic tree including three sequences obtained from this study. The publicly available sequences included complete bacterial, marine gastropoda, and algal alginate lyase coding sequences. The PRANK program was used for multiple sequence alignment of the dataset [27], and poorly aligned positions were examined and eliminated using the Gblocks program [28]. Sequence identity between aligned sequences was calculated using BioEdit software and some identical sequences were excluded [29]. Phylogenetic trees for the alginate lyase coding nucleotide sequences were constructed using the neighbor joining (NJ) method after replication of 1000 bootstraps using MEGA7 software [30–32]. Maximum composite likelihood was selected as a substitution model and substitution rates among all sites were assumed to be uniform [33].

Results and discussion

Microbial community of abalone visceral extract

We identified and quantified the abalone visceral microbiota at each taxonomic level from the visceral metagenome. Summary statistics of the metagenome for merged paired-end reads are provided in S3 Table. Alpha-diversity was 9.812 and a rarefraction curve is shown in S1 Fig. The bacterial taxa identified within the community and their relative abundance in the viscera of abalone from phylum to genus levels. At the genus level, the evolutionary relationships of the bacteria from the community were visualized along with their relative abundances (Fig 1). A greater diversity of taxa was identified here than in previous studies, which relied on more traditional microbiological methods such as bacterial culturing and cloning. While the total number of bacterial genera identified was less than ten in cultivation studies, paired-end reads from the metagenome in this study were classified into 43 genera [9, 34]. Additionally, differences in the dominant taxa and their relative abundances between previous cultivation studies and the present NGS study were identified at all taxonomic levels. From phylum to genus level, the bacterial abundance at each level in this study was consistent in that the proportion of the two dominant taxa constituted 41% and 33% of the microbiota, respectively. Among microbes identified at the phylum level, Tenericutes and Fusobacteria were the dominant microbes, which constituted 41.4% and 33.1% of the total microbial community, respectively. At the genus level, the majority of Tenericutes (phylum) was Mycoplasma (41.3%) and the majority of Fusobacteria (phylum) was Ilyobacter (33.0%). This suggests low diversity within each of these dominant phyla.

[Figure omitted. See PDF.]

Fig 1. Microbial community in abalone visceral extract.

Community tree of abalone visceral extract based on 16S rRNA gene sequences at the genus level. The tree displays the evolutionary relationships between bacterial genera. The thickness of the red bar indicates the relative abundance of each genus.

https://doi.org/10.1371/journal.pone.0205594.g001

Most abalone gut microbes, including Mycoplasma and Ilyobacter, have been found to degrade agar when cultured on agar plates [7, 35]. As agar is a mixture of polysaccharides derived from the cell walls of algae, agar degradation is indirect evidence of the ability to digest algal polysaccharides [36]. The agar-degrading bacteria can help the abalone to digest algae within its visceral tract. In the case of H. discus hannai, the digestive enzymes of gut microbes help the host utilize algae containing large amounts of polysaccharides [37].

In addition to Mycoplasma, Vibrio, which also demonstrated agar-degrading ability in culture-dependent studies, has previously been identified as a dominant genus in the abalone gut [7, 38]. However, Vibrio accounted for only 1.63% of the community at the genus level in our results. These conflicting results can be accounted for by temperature and differences in methodologies. First, the room temperature used to culture bacterial isolates in previous studies promoted the growth of Vibrio from abalone gut [39]. In contrast to those cultivation temperatures, our samples were collected at cold winter temperatures. The daily average air temperature of Wan Island ranged from 11.8°C to 14.3°C and the water temperature near Wan Island ranged from 4.3°C to 10.2°C (Feb., 2016, Korea Meteorological Administration) [40, 41]. In addition to temperature, culture media and aerobic conditions were potential causes of overestimation of the proportion of Vibrio. Most of the culture-based studies used alginate-peptone-yeast extract (APY) as a culture medium under aerobic conditions following the method described by Sawabe et al. [42]. This artificial medium and the aerobic conditions differ drastically from the natural conditions within the abalone visceral tract. As a result, Vibrio showed greater abundance in the culture-based study than in the molecular data of this study.

Another novel discovery was the high abundance of the Fusobacteria phylum, including the Ilyobacter genus (also known as Lysobacter). Previous studies of the gut microbiota from abalone, including H. discus hannai, identified Fusobacteria, but its proportion was not dominant [35, 38] as found in our results. This difference was likely due to the obligatory anaerobic and asaccharolytic nature of Fusobacteria and Ilyobacter. In most of the preparation methods used for the gut microbiota it is difficult to maintain the obligatory anaerobic conditions. As a result, the growth of obligate anaerobes, including Fusobacteria, is inhibited [43]. Additionally, the asaccharolytic nature of Fusobacteria inhibits its ability to produce energy under aerobic conditions [44]. Because of these two characteristics related to oxygen, the abundance of Fusobacteria and Ilyobacter would be underestimated when using traditional methods as opposed to a metagenomic approach.

In the first NGS analysis of abalone visceral extract described here, more diverse taxa were identified in detail, which were beyond the detection range of traditional methods, such as bacterial culturing and cloning. Moreover, the relative abundances of microbes could be precisely quantified based on read counts while traditional methods are dependent on ambiguous indicators such as colony forming units (CFU). These two improvements helped us to survey the actual state of the microbial community in abalone viscera.

Functional profile of the abalone visceral metatranscriptome

We annotated metatranscriptome functions using two databases: the SEED subsystems and RefSeq microbial genome databases. Fig 2 shows the annotation results from the SEED subsystems hierarchical database. At the most comprehensive hierarchy level, the ten most abundant major functions accounted for 87.8% of all annotated functions. The most abundant function present was protein metabolism (25.8%), followed by respiration (18.1%), carbohydrates (11.7%), regulation and cell signaling (8.8%), RNA metabolism (6.9%), and stress response (5.1%). A large proportion of the protein metabolism category corresponded to protein biosynthesis and protein degradation. Respiration, the second most predominant function, consisted of four subsystems and three of them (electron-accepting reactions, electron-donating reactions, and ATP synthases) were also prevalent in level 2. The carbohydrates term included 11 carbohydrate metabolisms and the representative function was central carbohydrate metabolism. The regulation and cell signaling group contained four subsystems, including a proteolytic pathway. RNA metabolism consisted of transcription, RNA processing, and RNA modification. The subsystems of stress response were classified based on causes of stress, such as oxygen and heat shock. Likewise, major subsystems of predominant functions at higher levels were also predominant at the lower level. For example, in level 3, ten of the predominant subsystems were included in three categories of level 1 such as ‘protein metabolism’, ‘respiration’, and ‘regulation and cell signaling’. In level 4, the total number of functions was 1957 and the top ten predominant functions accounted for 28.6% of the total.

[Figure omitted. See PDF.]

Fig 2. Hierarchical profiling of functional activities based on the SEED subsystems hierarchical database.

Each pie chart shows the proportion of the top ten functions hierarchically at: (a) level 1, (b) level 2, (c) level 3, and (d) level 4.

https://doi.org/10.1371/journal.pone.0205594.g002

The annotation results from the NCBI RefSeq microbial genomes database are summarized in Table 1. The top ten most abundant functions were similar to the results from the hierarchical level 4 of the SEED database. For example, cytochrome and translation elongation factors were predominant in both sets of annotation results. Furthermore, functions related to respiration, carbohydrate metabolism, and protein metabolism were also predominant. As these three functional categories play fundamental roles in all cells in the visceral environment, their abundances were higher than those of other functions. For example, protein synthesis and degradation, as well as respiration and central carbohydrate metabolism, are essential for all living organisms. Given that this is the first metatranscriptomic study of visceral extracts from marine mollusks, we compared our results to previous studies looking at transcriptomic data from host tissues in general. In comparison with the digestive gland of H. discus discus, transcripts related to cytochrome c oxidase and elongation factors were abundant in the visceral metatranscriptome. Alternatively, transcripts encoding digestive enzymes, such as cellulase and endo-1,4-β-glucanase (laminarase), were more abundant in H. discus discus because the transcriptome was collected from digestive gland [45]. In the case of visceral tissue from a Korean land snail (Koreanohadra kurodana), nucleotide metabolism was less abundant and protein metabolism was more abundant in H. discus hannai [46]. In short, the overall functional abundance of visceral extract was generally similar to that of host tissue except for protein metabolism function. The reason for the difference in functional activity of protein metabolism can ascribed to the active protein metabolism of the visceral microbiota.

[Figure omitted. See PDF.]

Table 1. The top ten enriched functional categories from the metatranscriptome obtained from the RefSeq microbial genomes database.

https://doi.org/10.1371/journal.pone.0205594.t001

Brown algae polysaccharide-degrading enzymes represented in the visceral metatranscriptome

Through metatranscriptomic analysis, we identified various expressed genes and related functions. Among them, we focused on genes involved in the digestion of algal polysaccharides and host energy metabolism. The main food and energy source of H. discus hannai is brown algae, Laminaria in particular [3]. Algal polysaccharides, such as alginate, β-glucans (laminarans), cellulose, heteroglycan, and fucoidan are major components of the brown algae [4, 5]. In general, herbivores derive a large amount of energy from dietary polysaccharides and require support from gut their microbiota to supply the digestive enzymes needed to improve digestion [47]. In addition to microbial breakdown, some herbivorous mollusks, including abalone, secrete their own polysaccharide-degrading enzymes [48]. In the case of H. discus hannai, alginate lyase has previously been isolated from the hepatopancreas using chromatography and cDNA cloning [49, 50]. Additionally, several species of gut bacteria from H. discus hannai showed an ability to degrade alginate during culturing experiments on agar plates [8, 42]. Although alginate digestion abilities have been identified previously in both H. discus hannai and their associated gut bacteria, the digestion mechanism in the visceral tract had not yet been explored. We examined the presence and origin of algal polysaccharide-degrading enzymes within the visceral extract metatranscriptome for the first time. After open reading frame (ORF) prediction, we performed a BLAST search of all sequences from the metatranscriptome. Two of the major enzymes required for algal polysaccharide degradation, alginate lyase and laminarinase, were identified. BLAST results suggest that H. discus hannai is the source of both enzymes in the dataset (Table 2 & Table 3). We constructed a phylogenetic tree of the alginate lyase coding sequences to clarify the origin of the alginate lyase coding transcripts identified in this study (Fig 3). The sequences identified in this study grouped with other sequences from Haliotis (Fig 3). Sequences from marine gastropods grouped together and were distinguishable from the groups of bacterial and algal alginate lyases. Although the origin of the alginate lyase-encoding transcript was identified as H. discus hannai, it does not mean that the visceral bacteria did not have any role in degrading algal polysaccharides within the abalone gut. In general, the gut microbiota of aquatic invertebrates supports the energy metabolism of its host by catalyzing the degradation of large polymers into smaller molecules, such as short chain fatty acids [8, 51]. Even though the mechanism has not yet been explained at the molecular level, the ability to degrade algal polysaccharides was identified in many species of visceral bacteria isolated from abalone in culturing experiments [7–9]. Likewise, previous studies have revealed that genes from both abalone and their visceral microbiota play a role in algal digestion. For instance, alginate lyase is essential for degrading alginate, which is the major component of brown algae. The results from this study, including BLAST analysis, a phylogenetic tree, and microbial community assessments, suggest that alginate lyase in the viscera of H. discus hannai was produced by the host itself, not by the visceral bacteria.

[Figure omitted. See PDF.]

Fig 3. Phylogenetic tree of alginate lyase coding sequences created using the neighbor-joining method.

Percentages of bootstrap values over 70% are displayed left of the nodes. Branch length equates to evolutionary distance. The origins of publicly available sequences are given with the corresponding GenBank accession numbers. The different origins of the publicly available sequences are distinguished by different colored boxes to the right of the terminal nodes. Nodes that include sequences from this study are highlighted with a red dot.

https://doi.org/10.1371/journal.pone.0205594.g003

[Figure omitted. See PDF.]

Table 2. BLAST analysis results for alginate lyase coding transcripts using the non-redundant nucleotide database.

https://doi.org/10.1371/journal.pone.0205594.t002

[Figure omitted. See PDF.]

Table 3. BLAST analysis results for laminarase-coding transcripts using the non-redundant nucleotide database.

https://doi.org/10.1371/journal.pone.0205594.t003

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Rarefaction curve of annotated species richness in visceral extract of Haliotis discus hannai.

https://doi.org/10.1371/journal.pone.0205594.s001

(DOCX)

S1 Table. Summary statistics of whole genome sequencing data generated for three Pacific abalone species.

https://doi.org/10.1371/journal.pone.0205594.s002

(DOCX)

S2 Table. Summary statistics for total transcriptome assembly from visceral extract using Trinity.

https://doi.org/10.1371/journal.pone.0205594.s003

(DOCX)

S3 Table. Summary statistics for 16S rRNA library data used in community analysis of visceral extract of Haliotis discus hannai.

https://doi.org/10.1371/journal.pone.0205594.s004

(DOCX)

Citation: Nam B-H, Jang J, Caetano-Anolles K, Kim Y-O, Park JY, Sohn H, et al. (2018) Microbial community and functions associated with digestion of algal polysaccharides in the visceral tract of Haliotis discus hannai: Insights from metagenome and metatranscriptome analysis. PLoS ONE 13(10): e0205594. https://doi.org/10.1371/journal.pone.0205594

References

1. Bouchet P, Gofas S. Haliotis Linnaeus, 1758: World Register of Marine Species; 2017 [cited 2017 07–27]. Available from: http://www.marinespecies.org/aphia.php?p=taxdetails&id=138050.

2. Abbot R, Dance S. Compendium of seashells. A full-colour guide to more than 4,200 of the world’s marine shells. California: Odyssey Publishing; 2000.

3. Bansemer MS, Qin JG, Harris JO, Howarth GS, Stone DA. Nutritional requirements and use of macroalgae as ingredients in abalone feed. Reviews in Aquaculture. 2016;8(2):121–35.

4. Rupérez P, Toledano G. Indigestible fraction of edible marine seaweeds. Journal of the Science of Food and Agriculture. 2003;83(12):1267–72.

5. Wells ML, Potin P, Craigie JS, Raven JA, Merchant SS, Helliwell KE, et al. Algae as nutritional and functional food sources: revisiting our understanding. Journal of applied phycology. 2017;29(2):949–82. pmid:28458464

6. Draget KI, Smidsrød O, Skjåk‐Bræk G. Alginates from algae. Biopolymers Online: Biology• Chemistry• Biotechnology• Applications. 2005;6.

7. Tanaka R, Shibata T, Miyake H, Mori T, Tamaru Y, Ueda M, et al. Temporal fluctuation in the abundance of alginate‐degrading bacteria in the gut of abalone Haliotis gigantea over 1 year. Aquaculture Research. 2016;47(9):2899–908.

8. Sawabe T, Setoguchi N, Inoue S, Tanaka R, Ootsubo M, Yoshimizu M, et al. Acetic acid production of Vibrio halioticoli from alginate: a possible role for establishment of abalone–V. halioticoli association. Aquaculture. 2003;219(1–4):671–9.

9. Erasmus JH, Cook PA, Coyne VE. The role of bacteria in the digestion of seaweed by the abalone Haliotis midae. Aquaculture. 1997;155(1–4):377–86.

10. Hugenholtz P. Exploring prokaryotic diversity in the genomic era. Genome biology. 2002;3(2):reviews0003. 1. pmid:11864374

11. Andrews S. FastQC: A quality control tool for high throughput sequence data. Reference Source. 2010.

12. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014:btu170.

13. Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63. pmid:21903629

14. Glass EM, Wilkening J, Wilke A, Antonopoulos D, Meyer F. Using the metagenomics RAST server (MG-RAST) for analyzing shotgun metagenomes. Cold Spring Harbor Protocols. 2010;2010(1):pdb. prot5368.

15. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Applied and environmental microbiology. 2006;72(7):5069–72. pmid:16820507

16. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nature biotechnology. 2011;29(7):644–52. pmid:21572440

17. Haas B, Papanicolaou A. TransDecoder (find coding regions within transcripts). 2016.

18. Westreich ST, Treiber ML, Mills DA, Korf I, Lemay DG. SAMSA2: A standalone metatranscriptome analysis pipeline. bioRxiv. 2017:195826.

19. Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics. 2013;30(5):614–20. pmid:24142950

20. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. pmid:24695404

21. Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28(24):3211–7. pmid:23071270

22. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic acids research. 2007;35(21):7188–96. pmid:17947321

23. Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR. Rfam: an RNA family database. Nucleic acids research. 2003;31(1):439–41. pmid:12520045

24. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nature methods. 2015;12(1):59–60. pmid:25402007

25. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic acids research. 2013;42(D1):D206–D14.

26. Tatusova T, Ciufo S, Fedorov B, O’Neill K, Tolstoy I. RefSeq microbial genomes database: new representation and annotation strategy. Nucleic acids research. 2013;42(D1):D553–D9.

27. Löytynoja A, Goldman N. Phylogeny-aware gap placement prevents errors in sequence alignment and evolutionary analysis. Science. 2008;320(5883):1632–5. pmid:18566285

28. Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Molecular biology and evolution. 2000;17(4):540–52. pmid:10742046

29. Hall TA, editor BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium series; 1999: [London]: Information Retrieval Ltd., c1979–c2000.

30. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molecular biology and evolution. 2016;33(7):1870–4. pmid:27004904

31. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular biology and evolution. 1987;4(4):406–25. pmid:3447015

32. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39(4):783–91. pmid:28561359

33. Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proceedings of the National Academy of Sciences of the United States of America. 2004;101(30):11030–5. pmid:15258291

34. Zhao J, Shi B, Jiang Q-r, Ke C-h. Changes in gut-associated flora and bacterial digestive enzymes during the development stages of abalone (Haliotis diversicolor). Aquaculture. 2012;338:147–53.

35. Tanaka R, Ootsubo M, Sawabe T, Ezura Y, Tajima K. Biodiversity and in situ abundance of gut microflora of abalone (Haliotis discus hannai) determined by culture-independent techniques. Aquaculture. 2004;241(1):453–63.

36. Chi W-J, Chang Y-K, Hong S-K. Agar degradation by microorganisms and agar-degrading enzymes. Applied microbiology and biotechnology. 2012;94(4):917–30. pmid:22526785

37. Sim S-J, Baik KS, Park SC, Choe HN, Seong CN, Shin T-S, et al. Characterization of alginate lyase gene using a metagenomic library constructed from the gut microflora of abalone. Journal of industrial microbiology & biotechnology. 2012;39(4):585–93.

38. Huang ZB, Guo F, Zhao J, Li WD, Ke CH. Molecular analysis of the intestinal bacterial flora in cage‐cultured adult small abalone, Haliotis diversicolor. Aquaculture Research. 2010;41(11):e760–e9.

39. Lee M-J, Lee J-J, Chung HY, Choi SH, Kim B-S. Analysis of microbiota on abalone (Haliotis discus hannai) in South Korea for improved product management. International Journal of Food Microbiology. 2016;234:45–52. pmid:27371902

40. Korea Meteorological Administration K. Analysis geometry of South Korea from December 2015 to February 2016. 2016.

41. Korea Meteorological Administration K. marine buoy cosmos 2016. Available from: http://www.kma.go.kr/mini/marine/marine_buoy_cosmos.jsp.

42. Sawabe T, Oda Y, Shiomi Y, Ezura Y. Alginate degradation by bacteria isolated from the gut of sea urchins and abalones. Microbial Ecology. 1995;30(2):193–202. pmid:24185485

43. Schink B, Janssen P, Brune A. Genus III. Ilyobacter Stieb and Schink 1985, 375 VP (Effective publication: Stieb and Schink 1984, 145). Bergey’s manual of systematic bacteriology. 2011;4:759–61.

44. Bennett K, Eley A. Fusobacteria: new taxonomy and related diseases. Journal of medical microbiology. 1993;39(4):246–54. pmid:8411084

45. Munasinghe H, Kang HS, Lee J. Analysis of digestive gland expressed sequence tag library from the disk abalone, Haliotis discus discus. Journal of the World Aquaculture Society. 2006;37(1):96–106.

46. Kang SW, Patnaik BB, Hwang H-J, Park SY, Chung JM, Song DK, et al. Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers. Molecular Genetics and Genomics. 2016;291(5):1999–2014. pmid:27507702

47. Flint HJ, Bayer EA, Rincon MT, Lamed R, White BA. Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nature reviews Microbiology. 2008;6(2):121. pmid:18180751

48. Ojima T. Polysaccharide-degrading enzymes from herbivorous marine invertebrates. Marine Enzymes for Biocatalysis: Elsevier; 2013. p. 333–71.

49. Shimizu E, Ojima T, Nishita K. cDNA cloning of an alginate lyase from abalone, Haliotis discus hannai. Carbohydrate research. 2003;338(24):2841–52. pmid:14667705

50. Suzuki H, Suzuki K-i, Inoue A, Ojima T. A novel oligoalginate lyase from abalone, Haliotis discus hannai, that releases disaccharide from alginate polymer in an exolytic manner. Carbohydrate research. 2006;341(11):1809–19. pmid:16697989

51. Harris JM. The presence, nature, and role of gut microflora in aquatic invertebrates: a synthesis. Microbial ecology. 1993;25(3):195–231. pmid:24189919