Depth-Dependent Mycoplankton Glycoside Hydrolase Gene Activity in the Open Ocean—Evidence from the Tara Oceans Eukaryote
Total Page:16
File Type:pdf, Size:1020Kb
The ISME Journal (2020) 14:2361–2365 https://doi.org/10.1038/s41396-020-0687-2 BRIEF COMMUNICATION Depth-dependent mycoplankton glycoside hydrolase gene activity in the open ocean—evidence from the Tara Oceans eukaryote metatranscriptomes 1,2 1,3 Nathan Chrismas ● Michael Cunliffe Received: 17 January 2020 / Revised: 20 May 2020 / Accepted: 21 May 2020 / Published online: 3 June 2020 © International Society of Microbial Ecology 2020 Abstract Mycoplankton are widespread components of marine ecosystems, yet the full extent of their functional role remains poorly known. Marine mycoplankton are likely functionally analogous to their terrestrial counterparts, including performing saprotrophy and degrading high-molecular weight organic substrates using carbohydrate-active enzymes (CAZymes). We investigated the prevalence of transcribed oceanic fungal CAZyme genes using the Marine Atlas of Tara Ocean Unigenes database. We revealed an abundance of unique transcribed fungal glycoside hydrolases in the open ocean, including a particularly high number that act upon cellulose in surface waters and the deep chlorophyll maximum (DCM). A variety of 1234567890();,: 1234567890();,: other glycoside hydrolases acting on a range of biogeochemically important polysaccharides including β-glucans and chitin were also found. This analysis demonstrates that mycoplankton are active saprotrophs in the open ocean and paves the way for future research into the depth-dependent roles of marine fungi in oceanic carbon cycling, including the biological carbon pump. Even though our understanding of marine fungal diversity is Many terrestrial fungi occupy key roles as saprotrophs by increasing [1, 2], a comprehensive knowledge of their active decomposing and recycling biogenic matter, making them functional ecology remains limited, especially in the open intrinsic components of healthy functioning ecosystems [3]. ocean [1]. Fungal activity has been detected in corals, deep In coastal waters there is evidence that planktonic fungi sea and coastal sediments, and associated with phyto- (mycoplankton) degrade and utilise phytoplankton-derived plankton blooms, including parasites [1, 2], but the full carbohydrate-rich matter in broadly analogous functional extent that fungi are functionally active throughout the open modes [4]. However, the extent that carbohydrate-based ocean water column is yet to be established. fungal saprotrophy occurs in the open ocean remains largely speculated [1]. Glycoside hydrolases (GHs) are a widespread group of carbohydrate-active enzymes (CAZymes) [5] that degrade Supplementary information The online version of this article (https:// complex polysaccharides and are categorised into substrate- doi.org/10.1038/s41396-020-0687-2) contains supplementary specific families. In terrestrial fungi, secreted CAZymes material, which is available to authorized users. are key to the functional potential of saprotrophs and are the primary mode of degradation of high-molecular weight * Nathan Chrismas [email protected] (HMW) polysaccharides (e.g. cellulose). Coastal sapro- * Michael Cunliffe trophic mycoplankton also employ secreted GHs to degrade [email protected] phytoplankton-derived HMW carbohydrate-based sub- strates [4], but the prevalence and identity of the specific 1 Marine Biological Association of the UK, The Laboratory, Citadel GH families of active open ocean mycoplankton are Hill, Plymouth, UK unknown. 2 School of Geographical Sciences, University of Bristol, University Metagenomes from the Tara Oceans project have been Road, Bristol, UK used to assess mycoplankton diversity [6, 7], but the asso- 3 School of Biological and Marine Sciences, University of ciated metatranscriptomes are yet to be fully explored from Plymouth, Plymouth, UK a fungal perspective. We interrogated the Marine Atlas 2362 N. Chrismas, M. Cunliffe of Tara Ocean Unigenes (MATOU) metatranscriptomic match represents a unique GH (similar genes clustered occurrences database [8] for transcribed fungal GH genes to when similarity < 95% over 90% of the smallest sequence explore broad-scale depth-dependent structuring in the [8]). The database was screened for non-fungal unigenes oceans. The MATOU database consists of all unique using the MATOU taxonomy (Fig. 1a, Supplementary eukaryotic genes assembled from the Tara Oceans meta- Fig. 1). After removal of redundant matches (i.e. where transcriptomes (unigenes), their associated taxonomy and multiple GHs matched to a single unigene), 1,326 unique occurrence within samples (full methods described in [8]). fungal GH unigenes were found (~0.001% of the entire To identify GHs within the MATOU database, refer- unigene catalogue) that occurred 44,386 times in all Tara ence libraries were created for 61 fungal GH families and Oceans samples. clustered using CD-Hit [9]. The MATOU unigenes were The top ten GH families containing the greatest number searched against these libraries using Diamond v.0.9.22 of unique genes were determined by ranking the sum of all [10], yielding a database where each positive unigene unigene occurrences from all samples for each family Fig. 1 Bioinformatics pipeline (a) (b) Total unique and glycoside unigene occurrences 10000 Fungal GH MATOU 8000 hydrolase unigene abundance. 2000 4000 6000 (a) CAZymes Sequences Pipeline describing the steps 0 involved in identifying fungal GH7 CAZymes within the Tara CD-Hit Oceans MATOU database. GH17 GH5 A fungal GH protein sequence Build database GH72 reference database was created GH16 from all the 61 characterised GH GH13 DIAMONDBLAST BLAST GH3 subfamilies found in fungi. The Databases database was consolidated by GH47 clustering sequences at 95 % GH18 identity using CD-Hit before GH43 GH Unigenes GH12 Diamond BLAST databases GH37 were generated for each GH32 subfamily. Unigenes were GH31 MATOU searched against each of these Filter GH11 Taxonomy 61 databases using the following GH135 thresholds: e value > 1e−30, GH45 score > 1, subject Cov > 75%, GH38 GH10 keeping only the best GH20 alignments. Positive matches GH1 were then screened using the GH62 Fungal Non-fungal MATOU taxonomy to Unigenes Unigenes GH30 discriminate between fungal and GH28 non-fungal unigenes. GH51 Occurrences of each unigene GH55 GH63 within the Tara Oceans Fungal MATOU Non-fungal MATOU Metatranscriptomic Metatranscriptomic GH152 transcriptomes were returned. Occurrences Occurrences GH78 (b) Fungal GH groups found in GH93 the MATOU database ranked by GH79 abundance over all Tara Oceans (c) GH81 samples. (c) Total numbers and GH2 250 taxonomy (including GH27 GH133 Ascomycota (green), GH71 Ascomycota Basdiomycota (orange), and 200 GH128 Unassigned (yellow)) of unique Basidiomycota Unassigned Fungal GH65 fungal unigenes from the ten GH115 most abundant GH groups. 150 GH35 GH132 GH75 100 GH36 GH53 GH6 50 GH29 Numbers of unique unigenes GH67 GH84 0 GH49 GH3 GH5 GH7 GH13 GH16 GH17 GH18 GH43 GH47 GH72 GH25 Depth-dependent mycoplankton glycoside hydrolase gene activity in the open ocean—evidence from the Tara. 2363 (a) (b) Surface 100% n=12 n=10 n=12 n=9 n=9 n=4 n=10 TARA_152 90% TARA_145 TARA_11TARA_23TARA_25 TARA_146 TARA_144 TARA_4 TARA_9TARA_22 TARA_26 80% TARA_132 TARA_143 TARA_150 TARA_30 TARA_142 TARA_151 70% TARA_131 TARA_135 TARA_148 TARA_20 TARA_36 TARA_38 TARA_137TARA_139 TARA_7 TARA_18 TARA_136 TARA_149 TARA_39 60% TARA_138 TARA_147 TARA_128 TARA_109 TARA_41 TARA_40 TARA_110 TARA_102 50% TARA_125 TARA_123 TARA_72 TARA_47 TARA_46 TARA_100 TARA_52 TARA_124 TARA_122 TARA_76 TARA_70 TARA_48 TARA_111 40% TARA_78TARA_68TARA_67 TARA_51 TARA_98 TARA_92 TARA_64 30% TARA_93 TARA_80 TARA_66 TARA_65 TARA_82 TARA_81 20% TARA_83 TARA_84 TARA_85 10% 0% [IO] [MS] [NAO][NPO] [SAO] [SO] [SPO] DCM 100% n=9 n=7 n=5 n=7 n=7 n=2 n=6 90% TARA_23 TARA_22 80% TARA_151TARA_9 TARA_25 TARA_132 TARA_135 TARA_143 TARA_36 TARA_150 TARA_7TARA_30 70% TARA_142 TARA_4 TARA_38 TARA_131TARA_137 TARA_18 60% TARA_41TARA_39 TARA_128 TARA_109 TARA_138 TARA_110 50% TARA_72 TARA_47 TARA_100 TARA_102 TARA_76 TARA_52 TARA_111 TARA_51 40% TARA_93 TARA_78 TARA_68 TARA_98 TARA_64 TARA_80 30% TARA_66 TARA_65 TARA_82 TARA_81 20% TARA_85 10% 0% [IO] [MS] [NAO][NPO] [SAO] [SO] [SPO] n=1 n=3 n=1 n=2 Mesopelagic 100% 90% 80% TARA_149 70% TARA_135 60% TARA_138 TARA_109 50% TARA_100 40% TARA_68 TARA_98 30% 20% 10% 0% Key to Regions: [NPO] North Pacific Ocean [NAO][NPO] [SAO] [SPO] [IO] Indian Ocean [SAO] South Atlantic Ocean Key to GH Groups: [MS] Mediterranean Sea [SO] Southern Ocean GH7 GH17 GH5 GH72 GH16 [NAO] North Atlantic Ocean [SPO] South Pacific Ocean GH13 GH3 GH47 GH43 GH18 (c) TARA_100 TARA_109 TARA_135 SUR DCM MES TARA_68 Depth TARA_138 TARA_98 SUR DCM MES 0306090 0306090 0306090 Total unique unigenes/station Fig. 2 Depth distribution of fungal glycoside hydrolases in the station for each of the major oceanic regions sampled. (c) Depth- global oceans. (a) Global map indicating Tara Oceans stations dependent partitioning of fungal GH unigenes in the surface (SUR), searched for fungal GH unigenes in the surface, deep chlorophyll DCM and mesopelagic (MES). maximum (DCM) and mesopelagic. (b) Mean unique unigenes/ 2364 N. Chrismas, M. Cunliffe (Fig. 1b). Overall, the greatest number of unique GHs were in line with phylogenetic studies that show the phylum involved in cellulose/hemicellulose degradation (GH7). dominates open ocean mycoplankton diversity [6]. However, Other substrates of the most abundant GH families included there is a lack of early-diverging taxa (e.g. Chytridiomycota) β-glucans (GH17 and GH72), β-glycans (GH5, GH16,