Soil and Rhizosphere Associated Fungi in Gray Mangroves(Avicennia marina)from theRed Sea—A Metagenomic Approach

Marta Filipa Simo˜esAndre´AntunesCristiane A.Ottoni Mohammad Shoaib AminiIntikhab AlamHanin Alzubaidy Noor-Azlin MokhtarJohn A.C.ArcherVladimir B.Bajic*i

1Computational Bioscience Research Center(CBRC),Computer,Electrical and Mathematical Sciences and Engineering Division(CEMSE),King Abdullah University of Science and Technology(KAUST),Thuwal 23955-6900,Saudi Arabia

2Laborato´rio de Biotecnologia Industrial(LBI),Instituto de Pesquisas Tecnolo´gicas do Estado de Sa˜o Paulo(IPT), Sa˜o Paulo 05508-901,Brazil

ORIGINAL RESEARCH

Soil and Rhizosphere Associated Fungi in Gray Mangroves(Avicennia marina)from the
Red Sea—A Metagenomic Approach

Marta Filipa Simo˜es1,a,Andre´Antunes1,b,Cristiane A.Ottoni2,c, Mohammad Shoaib Amini1,d,Intikhab Alam1,e,Hanin Alzubaidy1,f, Noor-Azlin Mokhtar1,g,John A.C.Archer1,h,Vladimir B.Bajic1,*,i

1Computational Bioscience Research Center(CBRC),Computer,Electrical and Mathematical Sciences and Engineering Division(CEMSE),King Abdullah University of Science and Technology(KAUST),Thuwal 23955-6900,Saudi Arabia

2Laborato´rio de Biotecnologia Industrial(LBI),Instituto de Pesquisas Tecnolo´gicas do Estado de Sa˜o Paulo(IPT), Sa˜o Paulo 05508-901,Brazil

Available online 5 November 2015

Handled by Fangqing Zhao

Fungaldiversity; Ascomycota; Basidiomycota; Red Sea; Metagenomics; Bioinformatics

Covering a quarter of the world’s tropical coastlines and being one of the most threatened ecosystems,mangroves are among the major sources of terrestrial organic matter to oceans and harbor a wide microbial diversity.In order to protect,restore,and better understand these ecosystems,researchers have extensively studied their microbiology,yet few surveys have focused on their fungal communities.Our lack of knowledge is even more pronounced for specific fungal populations,such as the ones associated with the rhizosphere.Likewise,the Red Sea gray mangroves(Avicennia marina)remain poorly characterized,and understanding of their fungal communities stillrelies on cultivation-dependent methods.In this study,we analyzed metagenomic datasetsfrom gray mangrove rhizosphere and bulk soil samples collected in the Red Sea coast,to obtain a snapshot of their fungal communities.Our data indicated that Ascomycota was the dominant phylum(76%–85%),while Basidiomycota was less abundant(14%–24%),yet present in higher numbers than usually reported for such environments.Fungalcommunities were more stable within the rhizosphere than within the bulk soil,both at class and genus level.This finding is consistent with the intrinsic patchiness in soilsediments and with the selection of specific microbial communities by plant roots.Our study indicates the presence of several species on this mycobiome that were not previously reported as mangrove-associated.In particular,we detected representatives of several commercially-used fungi,e.g.,producers of secreted cellulases and anaerobic producers of cellulosomes.These results represent additionalinsights into the fungal community of the gray mangroves of the Red Sea,and show that they are significantly richer than previously reported.

Introduction

Mangroves are endangered coastal biotopes that approximately cover a quarter of the world’s tropical coastlines [1–3].They are associated with a wide range of ecologicalbenefits,such as being a major source of terrestrialorganic matter to oceans and are wellrecognized,yet poorly studied,biodiversity hotspots[2,4].Microbes are major components of this biodiversity,with bacteria and fungiconstituting 91%of the total biomass of mangrove ecosystems[5],with the fungal fraction being the least studied.

Fungiare a ubiquitous and very diverse group of organisms currently comprising seven recognized phyla:Basidiomycota, Ascomycota,Glomeromycota,Microsporidia,Blastocladiomycota,Neocallimastigomycota,and Chytridiomycota[6]. Generally,fungi are important soil components as both decomposers and plant symbionts,playing major roles in ecological and biogeochemical processes[7].They contribute significantly to the degradation of mangrove-derived organic matter[8],being its primary mineralizers in mangrove sediments and representing important food source for benthic fauna[2].

Fungalsurveys in mangroves have focused mainly on taxonomic diversity of saprophytic fungi retrieved from intertidal, floating or immersed,pieces of trees and wood debris[9]. Diversity estimates pointed to 625 marine fungispecies associated with mangrove forests,and 269 related to mangrove roots [10].These mangrove fungi are almost exclusively saprophytic and belong primarily to the Ascomycota(e.g.,sac fungi and yeasts)and Basidiomycota(e.g.,mushrooms,rusts,and smuts),which are members of the subkingdom Dikarya[6,11].

A few studies analyzed mangrove-associated fungi [1,8,9,12,13].Highest counts are often found in soil surfaces or in roots and rhizomes,and some studies related their growth peak with higher humidity[14,15].Unfortunately,information on fungal diversity in mangrove rhizospheres,the soil zone located in and around the active roots,is lacking and is mostly based on culture-dependent assessments[16–18].As is well known and widely reported,traditional culturing techniques only succeed in isolating a very limited percentage of microorganisms and failto capture the full microbialdiversity present in the environment[19].Previous reports pointed out thatfrom the total(under)estimated 1.5×106fungal species,only ca. 8%–10%have been identified[6].Culture-independent techniques,e.g.,metagenomics,successfully circumvent such culture-based biases[7,20,21]and are essential for studying the real fungal diversity present in mangroves[1,4,8].

Contrasting with other seas,the Red Sea exhibits an antagonistic salinity-temperature profile:moving from south to north,surface water temperature decreases from 33.8°C down to 21°C;and,salinity increases from 37 to 41 PracticalSalinity Unit(PSU)[22].Such salinities,which are higher than the world average,are further increased in mangrove shallow waters[9,22].The high levels of stress imposed on the mangroves of the Red Sea result in scattered forests,decreased floral diversity,and limited plant height[23–26].Moreover, mangroves in the northern coastline of the Red Sea are mono-specific,and composed exclusively of Acivennia marina (gray mangrove)[9].

Information on fungal diversity in the gray mangroves of the Red Sea is scarce.In one of the very few studies available, Abdel-Wahad et al.[9]used a targeted metagenomic approach to look into fungaldiversity of the soiland rhizosphere in gray mangroves from the Red Sea.They recorded a totalof 29 different fungal species isolated from wood pieces on the mangroves and surrounding beaches,although the rhizosphere remained under-studied.

In order to decrease the paucity of data on fungal communities present in rhizosphere and in the gray mangroves of the Red Sea,we analyzed samples collected from this specific environment.Our results are a valuable addition that further clarifies our understanding of these communities.

Results and discussion

Eukaryotic and fungal representation within the soil and rhizosphere samples

Studies of four metagenomic samples from sediments of gray mangrove rhizosphere(RSMgr 01–04)and two samples from bulk soil(CS 01 and CS 02),publicly available under the project name‘‘A.marina rhizosphere”,were retrieved and analyzed at the metagenomics analyzer server(http:// metagenomics.anl.gov).These metagenomic datasets from gray mangroves of the Red Sea revealed that Eukaryota represent a relatively small percentage of allreads.Totalnumber of eukaryotic reads slightly increased from controlsamples(bulk soil)to rhizosphere sediments,(0.6%–0.7%and 2%of total reads,respectively),while fungal abundance was much higher in the rhizosphere sediments(Table 1).Such low abundances occur despite the widely-recognized importance of mangrove fungi,and the factthat they represent the second major ecological group of marine fungi(e.g.,[27,28]).We should note that Kuramae etal.[29]showed thatfungalabundance issignificantlycorrelated with phosphate,while frequent water logging and subsequent episodic anaerobic conditions were proposed as the possible explanations for the low abundance of fungi in some soils.Furthermore,it has been previously reported that the fungal abundance is lower in mangroves with smaller stands and tree-size,as well as less diverse regarding tree flora [30].The mangroves of the northern Red Sea show allof these features.Despite being very rich in carbon(C),mangrove soils are frequently nutrient-poor,with extremely low nutrient availability[31].Mangroves have evolved in tropical oligotrophic tidalenvironments with their soils having characteristically very low contents of total nitrogen(N)and phosphorous(P)[31].Such an effect is even more pronounced in the ultra-oligotrophic environment of the Red Sea[32]. Consequently,mangrove forests in the north of the Red Sea are sparse,with trees displaying decreased height and appearing in patchy and scattered patterns[23–26].

In order to confirm and better represent the fungaldiversity differences between CS and RSMgr samples,we performed principalcomponent analysis(PCA)as described in Materials and methods section with read counts atclass level,which provides enough analysis power.Figure 1 shows the CS and RSMgr samples at the class level in the plane with axes as first and second principal components,respectively.We observed that at the class level,the fungal communities of CS samples were distinct from those of RSMgr samples,with the latter displaying lower intra-group variability.In fact,this is also evident even with the first principal component(PC1)values.

Fungal abundance analysis at phylum level

In contrast to the aforementioned low abundances of Eukaryota in soil and rhizosphere samples,we detected a very high fungalabundance,particularly pronounced for the rhizosphere samples(Table 1).At the phylum level(Figure 2),fungal communities were clearly dominated by Ascomycota(76%–85%) and Basidiomycota(14%–24%).Members of these two phyla are expected to play an important ecological function in the mangroves[1].Ascomycetes from marine environments are an important ecological assembly of mostly saprophytic microbes occurring in different substrata rich in lignin,cellulose,or chitin[33].Other trophic levels are dependent on the lignocellulose-cleaving capability of these fungi that allow this complex substrate to enter the food web[33].Basidiomycetes are also mostly saprophytes[1],yet are mostly excluded from aquatic environments,leading to lower abundances[3].Other previous studies of soils[8,14,34],marine environments [35],and mangroves in general[9,10,12,36]pointed similarly to a predominance of Ascomycota and Basidiomycota. Nevertheless,we have found Basidiomycota to be more frequent here than has been described for other mangrove associated fungal communities(e.g.,[1,8,12,37]).

Figure 1 Principalcomponents analysis ofthe fungalcommunities in the Red Sea gray mangrove samplesAnalysis was performed based on read counts at class level.CS represents bulk soilsamples and RSMgr represents gray mangrove rhizosphere samples.PC1,first principal component,represents 60.2%of the variation in data;PC2,the second principal component,represents 29.9%of the variation in data.

Astudy by Lauber etal.[38]showed thatthe fluctuations in relative abundance of Ascomycota and Basidiomycota in different types of soilwere attributable to variations in C/Nratios and levels of P.Succinctly,P-rich soils contain more Ascomycota,and fewer Basidiomycota,while soils with higher C/N ratios have a higher prevalence of Basidiomycota.These findings are in agreement with our results.RSMgr 01,which possessed the highest P concentration in combination with lowest C/N ratio(considering organic matter and nitrate as proxies for C and N,Table S1),had the highest number of reads for Ascomycota(85%)and the lowest for Basidiomycota(15%), while opposite nutrient distribution was observed for CS 02, which had lower P content and a higher C/N ratio(Figure 2 and Table S1).Interestingly,CS 02 had the highest relative abundance of Basidiomycota(76%Ascomycota and 24% Basidiomycota).

Table 1 Sequencing reads for Eukaryota and fungi in different samples

Figure 2 Eukaryota distribution in different samples from Red Sea gray mangrovesThe read proportion of Eukaryota distribution in different samples.CS,bulk soil;RSMgr,gray mangrove rhizosphere.‘‘Others”include Neocallimastigomycota,Blastocladiomycota,Glomeromycota,and Chytridiomycota.

Fungal abundance analysis at class level

Despite being mostly similar at the phylum level,we noted significant differences at the class level among the samples from the soil(CS)and the rhizosphere(RSMgr),in both Ascomycota(Figure 3)and Basidiomycota(Figure 4).Within the phylum Ascomycota,CS 01 and CS 02 contained comparable percentage of the class Eurotiomycetes(41%and 42%),which is higher than that found in the RSMgr samples(27%–31%). Similarly,more Schizosaccharomycetes were also found in the CS samples(16%and 25%)in comparison to the RSMgr samples(5%–8%).

When comparing the fungalcommunities across allsix samples(Figure 5),we noted that samples can be grouped according to the diversity present.CS samples group apart from the RSMgr samples,while RSMgr 03 and RSMgr 04 show more similar diversities.But the most significant differences were found in the increased percentage of Saccharomycetes(25%–38%vs.16%–18%;P=9.56E-3)and Sordariomycetes (18%–28%vs.11%–15%;P=9.45E-3)in the RSMgr samples when compared to CS samples but decreased content of Eurotiomycetes(26%–31%vs.41%–42%;P=3.90E-4),as shown in Figure 6.

Class Dothideomycetes accounted for 3%–6%of Ascomycota in RSMgr samples(Figure 3).However,high variability was found in the percentage of class Dothideomycetes in the CS samples(15%for CS 01 and 0.6%for CS 02).Even though rhizosphere samples show more class variability,this was not the case for Dothideomycetes.

Within the phylum Basidiomycota(Figure 4),allthe RSMgr samples possessed similar percentage of Ustilaginomycetes (22%–28%).However,RSMgr 01,RSMgr 02,and RSMgr 04 had similar content of Agaricomycetes(43%–46%),Tremellomycetes(23%–24%)and Exobasidiomycetes(6%–8%), whereas RSMgr 03 had much lower percentage of Agaricomycetes(28%)butmore Tremellomycetes(32%)and Exobasidiomycetes(13%).On the other hand,CS samples appeared to have a very different composition of Basidiomycota:75% of Ustilaginomycetes and 24%Agaricomycetes in CS 01; 67%of Tremellomycetes and 32%of Exobasidiomycetes in CS 02).

Fungal abundance analysis at genus and species level

Overall,we found that relative fungal diversity is more stable within the RSMgr than within the CS.Such results are consistent with the intrinsic environmental and biological patchiness in soil sediments,and with the preferential selection of specific microbial communities by plant roots.

Figure 3 Ascomycota distribution in different samples from Red Sea gray mangrovesThe read proportion of Ascomycota distribution in different samples.CS,bulk soil;and RSMgr,gray mangrove rhizosphere.‘‘Others”include Leotiomycetes,Pneumocystidomycetes,Lecanoromycetes,Orbiliomycetes,and Pezizomycetes.

A similar trend was observed at the genus level,with identicalgenera listed in the top ten for RSMgr,and very different profiles for CS(Table 2).It was noticeable that the genus Aspergillus and Schizosaccharomyces clearly dominate in all the samples examined in this study.Both genera are known producers and secretors of a large variety of heterologous proteins[39].For example,Aspergillus oryzae[40]and Schizosaccharomyces pombe[41],the well-known and well-studied examples of heterologous proteins producers,were identified in large amounts and in most of these samples for both bulk and rhizosphere soil(Table S2).Dominance of Aspergillus is in accordance with previous studies[42]reporting that this genus is frequently found in marine sediments.Such dominance also agrees with the isolations made by Thamizhmani and Senthilkumaran[43]from mangrove sediments,where they found several different species of this genus.In addition to Aspergillus,they also identified Emericella and Neurospora in their samples,which were also found in our samples, although the abundance of Emericella was low.

From the metadata hereby presented,a total of 145 different species within 109 different genera were identified (Table S2).Our results bring to light the existence of many additional species on the mycobiome of Red Sea mangrove rhizosphere that were not previously reported as mangroveassociated,since most prior information is based on culturedependent research.Furthermore,we anticipate total fungal diversity in the gray mangroves of the Red Sea to be even higher than shown by our study.It is well known that geochemicalparameters(e.g.,salinity,soilhumidity,and nutrients concentrations)for mangroves vary cyclically,throughout the day,with tides,and with seasons[1].To fully capture totalfungal diversity,additional sampling at different times and locations is essential while taking these variations into account.

Bioprospecting for potentially biotechnological interesting microbes is one of the many utilities that derive from characterizing microbial communities.This task has been largely improved by the use of culture-independent techniques[21]. Metagenomic analysis of fungal diversity of rhizospheres and sediments of mangroves would facilitate the discovery of novel enzymes,bioactivities,and relevant secondary metabolites. Mangroves are a source of cellulosic substrates and are at the transition between terrestrial and aquatic environment. They are a dynamic ecotone subjected to harsh conditions, with fluctuating temperature,salinity,and tides[3].These environmentalcharacteristics lead to fungalcommunity specialized in producing a wide array of hydrolytic enzymes such as cellulases[10].This class of enzymes is under intense study for their involvement in biofuel production from renewable cellulosic substrates[19],as favorite candidates for industrial and biotechnology applications[10].Sahoo et al.[44]noted thatmangrove soil is a good source of microbes able to degrade polythene and plastics.In our samples,we found representatives of commercially-used producers of secreted cellulases: Aspergillus,Fusarium,Phanerochaete,and Penicillium;and also anaerobic fungi producers of cellulosomes,complexes of cellulases with high molecular weight:Orpinomyces, Piromyces,and Neocallimastix[10].

Figure 4 Basidiomycota distribution in different samples from Red Sea gray mangrovesThe read proportion of Basidiomycota distribution in different samples.CS,bulk soil;RSMgr,gray mangrove rhizosphere.‘‘Unclassified”are derived from Basidiomycota.

Fungal diversity analyses

Alpha-diversities for the total amount of reads were obtained using the Metagenomic Rapid Annotations with Subsystems Pipeline(MG-RAST)pipeline,as a means to further quantify fungal diversity of the annotated samples.The following evaluation using the species-level annotation distribution showed higher total species diversity in CS samples(Table 3).However,as shown in the previous sections,fungal diversity revealed higher richness in RSMgr samples(Table 3).It is well known that as species richness and evenness(the measure of relative abundance of the different species)increase,so does diversity.Simpson index measures the probability that two individuals randomly selected from a sample belong to the same species,which relates richness with evenness of the population,with higher Simpson index suggesting lower diversity.

As shown in Table 3,RSMgr samples had much lower Simpson indexes,thus demonstrating a higher fungal diversity when compared with CS samples.

Overview and future work

Our study reveals that the diversity of fungal communities in the gray mangroves of the Rea Sea is significantly wider than previously reported.Future studies on fungalcommunity characterizations and bioprospection are recommended for these particular environments(including more thorough sampling efforts),because there is a high probability of very interesting findings.Due to the unique environmental characteristics of the rhizosphere of mangroves,it represents a unique and under-explored source for a poolof uncommon fungiwith particular features of relevance for biotechnology,science,and health research.

Materials and methods

Sample collection

Figure 5 Heat map of the relative abundances of the fungal communities in Red Sea gray mangroves samplesHeat map of the relative abundances of the fungal communities generated with read counts by using STAMP software.CS,bulk soil; RSMgr,gray mangrove rhizosphere.

Sample collection was performed from six differentsites along a 978 m transect of mangrove shore in Thuwal,Saudi Arabia,in December 2011[45].Four samples from sediments ofgray mangrove rhizosphere(RSMgr 01–04),and two samples from bulk soil(CS 01 and CS 02)as controlwere collected.Itis important to note that,at each site,samples were collected from a 10-cm depth aseptically and stored at4°Cpriorto subsequentprocessing within 12 h.Chemical analyses for each sample were performed as follows.Briefly,phosphorous concentrations were measured with microwave-assisted digestion method[45]. Nitrate content was measured directly using Autoanalyzer/ Photometric Analyzer,Aquakem250(Thermo Scientific,Vantaa, Finland).Content of organic matter was calculated with the loss on ignition method[45,46],which was 9.21%-10.12% for RSMgr samples and 2.53%-3.19%for CS samples, respectively.The temperature and salinities were measured with a 5 Star pH/ISE/ORP/DO Conductivity Portable Meter (Thermo Fisher Scientific,Waltham,Massachusetts,USA). The temperature was 21.2°C for allthe samples.The salinities ofthe RSMgr samples were 18.65–23.38 PSU,whereas CS samples had salinities of 8.40–14.23 PSU.

Dataset acquisition

Figure 6 Read percentage of Eurotiomycetes,Sordariomycetes,and Saccharomycetes in CS and RSMgr samplesRead count of the three most statistically relevant classes,Eurotiomycetes,Sordariomycetes,and Saccharomycetes varied in CS and RSMgr samples.A.Comparison of the mean proportions of the three classes.B.Individual bar plot of the three classes.CS represents bulk soiland is marked in blue,whereas RSMgr represents gray mangrove rhizosphere and is marked in orange.The comparisons were made between the group of CS and the group of RSMgr samples(P=3.90E-4 for Eurotiomycetes,9.45E-3 for Sordariomycetes,and 9.56E-3 for Saccharomycetes,respectively).

Metagenomes were obtained through DNA extraction of each sample using the ZR Soil Microbe DNA MidiPrep kit(Zymo Research,Irvine,CA)and pyrosequencing with 454 GS FLX Titanium(Roche Applied Science,Indianapolis,Indiana, USA)in the Bioscience Core Laboratory,King Abdullah University of Science and Technology(Thuwal,SaudiArabia). These data,publicly available on http://metagenomics.anl.gov/, under the project name‘‘A.marina rhizosphere”,detailed on Table S3,were retrieved and analyzed through the MG-RAST [47]at the metagenomics analyzer server.

We compared the data to M5NR using the maximum cutoff E-value of 1E-5;we used 60%as the minimum identity cutoff,and the minimum alignment length cutoff of 15,measured in base pairs for RNA databases and in amino acids for proteins.

Statistical analyses

The similarity among the collected samples was analyzed with principal component analysis(PCA),based on relative abundance at class level using domain as the parent level.

Alpha-diversities,the number of distinct species in a given sample,were achieved by the distribution of the species-level annotations(total species from all taxonomic domains) obtained from MG-RAST.

Simpson index for the fungal population was calculated, using the formula:

where nirepresents the total number of organisms of a particular species and N represents the totalnumber of organisms of all species[48].

Heat map was generated using the Statistical Analyses of Metagenomic Profiles(STAMP)software[49]for fungal relative abundances.The associated dendrograms were obtained using the Unweighted Pair Group Method with Arithmetic Mean(UPGMA)with a clustering threshold of 0.75.

Relevant differences in the relative proportions of the classified sequences and mean proportions ofthe most statisticallyrelevant classes found in the samples were detected after analyzing the MG-RAST taxonomic profiles with STAMP. Datasets were analyzed with the two-sided Welch’s test,and we removed all unclassified reads from the analysis.P values of 0.05 were used as a filter to determine the most important taxa,and we only used those categories with more than 2-fold ratio between the proportions and with difference between the proportions of at least 1%.

Authors’contributions

MFS conceived,designed and performed the assays.MFS, AA,and CO analyzed the data.MSA,IA,HA,NAM,and JACA contributed with reagents/materials/analysis tools. MSA performed the statistical analyses.MFS,AA,and CO drafted the manuscript with the help of MSA,IA,and VBB. All authors read and approved the final manuscript.

?

Table 3 Species diversity in the different samples

Competing interests

The authors declare no competing interests.

Acknowledgments

This study was supported by the base research funds to VBB, and the competitive research funding of VBB from King Abdullah University of Science and Technology(KAUST)in Saudi Arabia.

Supplementary material

Supplementary material associated with this article can be found,in the online version,at http://dx.doi.org/10.1016/ j.gpb.2015.07.002.

References

[1]Arfi Y,Bue´e M,Marchand C,Levasseur A,Record E.Multiple markers pyrosequencing reveals highly diverse and host-specific fungal communities on the mangrove trees Avicennia marina and Rhizophora stylosa.FEMS Microbiol Immunol 2012;79:433–44.

[2]Ghizelini AM,Mendonc¸a-Hagler LCS,Macrae A.Microbial diversity in Brazilian mangrove sediments:a mini review.Braz J Microbiol 2012;43:1242–54.

[3]Shearer CA,Descals E,Kohlmeyer B,Kohlmeyer J,Marvanova´ L,Padgett D,et al.Fungal biodiversity in aquatic habitats. Biodivers Conserv 2007;16:49–67.

[4]Mendes R,Garbeva P,Raaijmakers JM.The rhizosphere microbiome:significance of plant beneficial,plant pathogenic, and human pathogenic microorganisms.FEMS Microbiol Rev 2013;37:634–63.

[5]Alongi DM.Bacterial productivity and microbial biomass in tropical mangrove sediments.Microb Ecol 1988;15:59–79.

[6]Simo˜es MF,Pereira L,Santos C,Lima N.Polyphasic identification and preservation of fungal diversity:concepts and applications.In:Malik A,Grohmann E,Alves M,editors.Management of microbial resources in the environment.Dordrecht,Netherlands:Springer;2013.p.91–117.

[7]Liu P,Wang XH,Li JG,Qin W,Xiao CZ,Zhao X,et al. Pyrosequencing reveals fungal communities in the rhizosphere of Xinjiang jujube.Biomed Res Int 2015;2015:972481.

[8]Arfi Y,Marchand C,Wartel M,Record E.Fungal diversity in anoxic-sulfidic sediments in a mangrove soil.Fungal Ecol 2012;5:282–5.

[9]Abdel-Wahab MA,Hodhod MS,Bahkali AH,Jones EB.Marine fungi of Saudi Arabia.Bot Mar 2014;57:323–35.

[10]Alsheikh-Hussain A,Altenaiji EM,Yousef LF.Fungalcellulases from mangrove forests—a short review.J Biochem Technol 2014;5:765–74.

[11]Sarma VV,Hyde KD.A review on frequently occurring fungi in mangroves.Fungal Divers 2001;8:1–34.

[12]Alias SA,Zainuddin N,Jones EB.Biodiversity of marine fungiin Malaysian mangroves.Bot Mar 2010;53:545–54.

[13]Schmit JP,Shearer CA.A checklist of mangrove-associated fungi, their geographical distribution and known host plants.Mycotaxon 2003;85:423–77.

[14]Alongi DM.5.Mangrove–Microbe–Soil Relations.In:Kristensen E,Haese RR,Kostka JE,editors.Interactions between macroand microorganisms in marine sediments.American Geophysical Union;2005.p.85–103.

[15]Nambiar GR,Raveendran K.Manglicolous marine fungi on Avicennia and Rhizophora along Kerala coast(India).Middle East J Sci Res 2009;4:48–51.

[16]Lee BK,Baker GE.Fungi associated with the roots of red mangrove.Rhizophora mangle.Mycologia 1973;65:894–906.

[17]Selvakumar V,Panneerselvam A,Vijayakumar N,Savery MA, Thajuddin N.Diversity of endophytic and rhizosphere soil fungi of Avicennia marina in Maravakadu Mangrove Forest.IOSR J Pharm Biol Sci 2014;9:24–8.

[18]Thiripurasundari G,Usharani G,Parthasarathi R,Neelakandan T.Microbialdiversity of the rhizosphere soil of Avicennia marina and Avicennia officinalis collected from mangrove forest.Int J Recent Sci Res 2010;1:039–41.

[19]Kennedy J,Flemer B,Jackson SA,Lejon DP,Morrissey JP, O’gara F,et al.Marine metagenomics:new tools for the study and exploitation of marine microbial metabolism.Mar Drugs 2010;8:608–28.

[20]Cuadros-Orellana S,Leite LR,Smith A,Medeiros JD,Badotti F, Fonseca PLC,et al.Assessment of fungal diversity in the environment using metagenomics:a decade in review.Fungal Genom Biol 2013;3:110.

[21]Guo X,Zhang Q,Zhang X,Zhang J,Gong J.Marine fungal communities in water and surface sediment of a sea cucumber farming system:habitat-differentiated distribution and nutrients driving succession.Fungal Ecol 2015;14:87–98.

[22]Robitzch V,Banguera-Hinestroza E,Sawall Y,Al-Sofyani A, Voolstra CR.Absence of genetic differentiation in the coral Pocillopora verrucosa along environmental gradients of the Saudi Arabian Red Sea.Front Mar Sci2015.http://dx.doi.org/10.3389/ fmars.2015.00005.

[23]Khan M,Kumar A,Muqtadir A.Distribution of mangroves along the Red Sea coast of the Arabian Peninsula:Part 2.The southern coast of western SaudiArabia.e-JournalEarth SciIndia 2010;3:154–62.

[24]Kumar A,Khan MA,Muqtadir A.Distribution of mangroves along the Red Sea coast of the Arabian Peninsula:Part-I:the northern coast of western SaudiArabia.e-JournalEarth SciIndia 2010;3:2842.

[25]PERSGA/GEF.Status of mangroves in the Red Sea and Gulf of Aden.PERSGA Technical Series No.11,PERSGA,Jeddah; 2004.http://www.persga.org/Documents/Mangroves_Status.pdf. [26]Price ARG,Medley PAH,McDowall RJ,Dawson-Shepherd AR, Hogarth PJ,Ormond RFG.Aspects of mangal ecology along the Red Sea coast of Saudi Arabia.J Nat Hist 1987;21:449–64.

[27]Cheng ZS,Pan JH,Tang WC,Chen QJ,Lin YC.Biodiversity and biotechnological potential of mangrove-associated fungi.J For Res 2009;20:63–72.

[28]Sridhar KR.Mangrove fungi in India.Curr Sci 2004;86:1586–7.

[29]Kuramae EE,Yergeau E,Wong LC,Pijl AS,van Veen JA, Kowalchuk GA.Soil characteristics more strongly influence soil bacterial communities than land-use type.FEMS Microbiol Ecol 2012;79:12–24.

[30]Hyde KD,Lee SY.Ecology of mangrove fungi and their role in nutrient cycling:what gaps occur in our knowledge?Hydrobiologia 1995;295:107–18.

[31]Reef R,Feller IC,Lovelock CE.Nutrition of mangroves.Tree Physiol 2010;30:1148–60.

[32]Ngugi DK,Antunes A,Brune A,Stingl U.Biogeography of pelagic bacterioplankton across an antagonistic temperature—salinity gradient in the Red Sea.Mol Ecol 2012;21:388–405.

[33]Velez P,Gonza´lez MC,Rosique-Gil E,Cifuentes J,del Rocı´o Reyes-Montes M,Capello-Garcı´a S,et al.Community structure and diversity of marine ascomycetes from coastal beaches of the southern Gulf of Mexico.Fungal Ecol 2013;6:513–21.

[34]Mandeel QA.Microfungal community associated with rhizosphere soilof Zygophyllum qatarense in arid habitats of Bahrain.J Arid Environ 2002;50:665–81.

[35]Edgcomb VP,Beaudoin D,Gast R,Biddle JF,Teske A.Marine subsurface eukaryotes:the fungal majority.Environ Microbiol 2011;13:172–83.

[36]Abdel-Azeem AM,Salem FM.Biodiversity of laccase producing fungi in Egypt.Mycosphere 2012;3:900–20.

[37]Richards TA,Jones MD,Leonard G,Bass D.Marine fungi:their ecology and molecular diversity.Ann Rev Mar Sci 2012;4:495–522.

[38]Lauber CL,Strickland MS,Bradford MA,Fierer N.The influence of soil properties on the structure of bacterial and fungal communities across land-use types.Soil Biol Biochem 2008;40:2407–15.

[39]Petrescu-Da˘nila˘E,Voicu M,Sta˘nescu R,Stoica B,Rusu M. Fission yeast Schizosaccharomyces pombe as a producer and secretor of heterologous proteins.Rom Biotechnol Lett 2009;14:4201–10.

[40]Oda K,Kakizono D,Yamada O,Iefuji H,Akita O,Iwashita K. Proteomic analysis of extracellular proteins from Aspergillus oryzae grown under submerged and solid-state culture conditions. Appl Environ Microbiol 2006;72:3448–57.

[41]Takegawa K,Tohda H,Sasaki M,Idiris A,Ohashi T, Mukaiyama H,et al.Production of heterologous proteins using the fission-yeast(Schizosaccharomyces pombe)expression system. Biotechnol Appl Biochem 2009;53:227–35.

[42]Jones EBG,Sakayaroj J,Suetrong S,Somrithipol S,Pang KL. Classification of marine Ascomycota,anamorphic taxa and Basidiomycota.Fungal Divers 2009;35:1–187.

[43]Thamizhmani R,Senthilkumaran R.Diversity of fungiin selected mangroves along the east coast of India.Int J Curr Microbiol Appl Sci 2012;1:29–33.

[44]Sahoo K,Dhal NK.Potential microbial diversity in mangrove ecosystems:a review.Indian J Marine Sci 2009;38:249–56.

[45]Alzubaidy H,Essack M,Malas TB,Bokhari A,Motwalli O, Kamanu FK,et al.Rhizosphere microbiome metagenomics of gray mangroves(Avicennia marina)in the Red Sea.Gene 2015. http://dx.doi.org/10.1016/j.gene.2015.10.032.

[46]Dean WE.Determination of carbonate and organic matter in calcareous sediments and sedimentary rocks by loss on ignition: comparison with other methods.J Sed Petrol 1974;44:242–8.

[47]Meyer F,Paarmann D,D’Souza M,Olson R,Glass EM,Kubal M,et al.The metagenomics RAST serve:a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 2008;9:386–94.

[48]Zak JC,Willig MR.Fungalbiodiversity patterns.In:Mueller G, Foster M,Bills G,editors.Biodiversity of fungi:inventory and monitoring methods.Amsterdam:Academic Press;2004.p.59–75.

[49]Parks DH,Tyson GW,Hugenholtz P,Beiko RG.STAMP: statistical analysis of taxonomic and functional profiles.Bioinformatics 2014;30:3123–4.

13 May 2015;revised 8 July 2015;accepted 22 July 2015

.

E-mail:vladimir.bajic@kaust.edu.sa(Bajic VB).aORCID:0000-0002-8767-9487.bORCID:0000-0001-7668-9842.cORCID:0000-0003-4069-509X.dORCID:0000-0001-7486-7190.eORCID:0000-0001-5306-847X.fORCID:0000-0002-3388-4039.gORCID:0000-0002-0863-7597.hORCID:0000-0002-3302-3933.iORCID:0000-0001-5435-4750.

Peer review under responsibility of Beijing Institute of Genomics, Chinese Academy of Sciences and Genetics Society of China.