Search tips
Search criteria 


Logo of frontmicrobioLink to Publisher's site
Front Microbiol. 2016; 7: 685.
Published online 2016 May 12. doi:  10.3389/fmicb.2016.00685
PMCID: PMC4922216

Concentration of Petroleum-Hydrocarbon Contamination Shapes Fungal Endophytic Community Structure in Plant Roots


Plant-root inhabiting fungi are a universal phenomenon found in all ecosystems where plants are able to grow, even in harsh environments. Interactions between fungi and plant roots can vary widely from mutualism to parasitism depending on many parameters. The role of fungal endophytes in phytoremediation of polluted sites, and characterization of the endophytic diversity and community assemblages in contaminated areas remain largely unexplored. In this study, we investigated the composition of endophytic fungal communities in the roots of two plant species growing spontaneously in petroleum-contaminated sedimentation basins of a former petro-chemical plant. The three adjacent basins showed a highly heterogeneous pattern of pollutant concentrations. We combined a culture-based isolation approach with the pyrosequencing of fungal ITS ribosomal DNA. We selected two species, Eleocharis erythropoda Steud. and Populus balsamifera L., and sampled three individuals of each species from each of three adjacent basins, each with a different concentration of petroleum hydrocarbons. We found that contamination level significantly shaped endophytic fungal diversity and community composition in E. erythropoda, with only 9.9% of these fungal Operational Taxonomic Units (OTUs) retrieved in all three basins. However, fungal community structure associated with P. balsamifera remained unaffected by the contamination level with 28.2% of fungal OTUs shared among all three basins. This could be explained by the smaller differences of pollutant concentrations in the soil around our set of P. balsamifera sampless compared to that around our set of E. erythropoda samples. Our culture-based approach allowed isolation of 11 and 30 fungal endophytic species from surface-sterilized roots of E. erythropoda and P. balsamifera, respectively. These isolates were ribotyped using ITS, and all were found in pyrosequensing datasets. Our results demonstrate that extreme levels of pollution reduce fungal diversity and shape community composition in E. erythropoda. Our findings shed light on the effect of soil petroleum contamination on fungal endophytic communities and could help to develop strategies for improving phytoremediation using fungal endophytes.

Keywords: fungi, endophytes, petroleum-hydrocarbon contamination, pyrosequencing, ribosomal DNA, Eleocharis erythropoda, Populus balsamifera L.


Phytoremediation uses plants to remove pollutants from contaminated soils (Cunningham and Berti, 1993). This decontamination method is cost-effective as well as environmentally friendly biotechnology compared to conventional methods (Salt et al., 1998; Pilon-Smits, 2005). Traditional approaches usually involve excavation and transportation of soils for ex situ chemical extraction of pollutants, an expensive and energy demanding technology, or disposal, a process known as Dig and Dump (Chaudhry et al., 2005; Göhre and Paszkowski, 2006; Kaimi et al., 2007). However, phytoremediation has several downsides, including a generally lower efficiency of pollutant removal, longer treatment and a depth of treated soils limited to the zone colonized by plant roots (Chaudhry et al., 2005; Pilon-Smits, 2005). Thus many attempts have been made to enhance phytoremediation efficiency by selecting high-performing plant species and cultivars, and their associated microorganisms, such as rhizospheric microbes, including symbiotic and endophytic microorganisms (Hassan et al., 2011, 2013; Bell et al., 2014).

Endophyte is defined literally as any organism living inside a plant (from Greek endon = inside and phyton = plant). The term endophyte usually refers to bacteria (Kobayashi and Palumbo, 2000) and fungi (Stone et al., 2000), but algae (Peters, 1991), and insects (Feller, 1995) can also be endophytes. Endophytic bacteria and fungi have been found in all plants studied so far, providing examples of the broad variety of colonizing organisms as well as hosts for endophytic interactions (Schulz and Boyle, 2005; Hyde and Soytong, 2008). There are many different types of interactions between endophytic organisms and plants spanning from mutualism (Sieber, 2002) to parasitism (Kogel et al., 2006; Paszkowski, 2006). Plants may benefit from the presence of endophytes such as increased growth (Petrini, 1991; Varma et al., 1999; Schulz et al., 2006; Rodriguez et al., 2009), and resistance against biotic (Bacon et al., 1977; Carroll, 1988; Melo et al., 2009) and abiotic stresses, e.g., drought tolerance (Schardl et al., 2004; Hamilton and Bauerle, 2012).

The biodiversity of endophytic organisms has been increasingly investigated in recent decades, and methods for studying endophytes within plants have evolved during that time. Endophytic detection was first carried out using plant tissue staining and microscopy (Jumpponen and Trappe, 1998; Jumpponen, 2001; Porras-Alfaro and Bayman, 2011). Isolation of endophytes using surface-sterilized plant tissues on growth media has been intensively used for bacteria and fungi (Schulz and Boyle, 2005; Schulz et al., 2006). These culture-dependent methods have been shown to provide better results than microscopic observation, although they require considerable effort and time, and can also be biased depending on the type of media used or sampling effort (Hyde and Soytong, 2008; Stefani et al., 2015). In recent years, the study of endophytes has shifted towards culture-independent methods using molecular approaches such as PCR, cloning and sequencing, and next-generation sequencing (Ganley et al., 2008). These molecular approaches can be directly applied to any plant tissue and have proven useful to study changes in community structure.

Despite the increasing number of studies published on endophytes, interactions between fungal endophytic organisms and their hosts are still poorly understood, particularly in the context of phytoremediation (Soleimani et al., 2010). Furthermore, bacterial endophytes have traditionally received more attention, and been more intensively studied in a wider range of contaminant degradations than fungal endophytes (reviewed in Rylott, 2014). While the role of fungi in plant growth-promotion remains nearly unexplored in phytoremediation (Li et al., 2012), some authors have used bacterial and fungal consortia to enhance phytoremediation using poplar (Knoth et al., 2014). It has been shown that the use of endophytic fungi enhances phytoremediation applications due to fungal developmental plasticity and fungal production of a variety of secondary metabolites and enzymes that can degrade or sequester pollutants (Tan and Zou, 2001; Schulz et al., 2002; Hyde and Soytong, 2008). Soleimani et al. (2010) have studied the effect of two-grass species Festuca arundinacea Schreb. and Festuca pratensis Huds. inoculated with endophytic fungi Neotyphodium coenophialum and Neotyphodium uncinatum, respectively, on the degradation of petroleum hydrocarbons in a contaminated soil. Although the authors used plants and endophytic strains that are not native to the contaminated site and the fungi used for inoculation belong to the Balansia group, a group of well-studied endophytic fungi known for their beneficial roles for grasses (Bacon et al., 1977; Schulz et al., 1993), they found that inoculated plants were able to decontaminate soils better than non-inoculated plants (Soleimani et al., 2010).

The aim of this study is to investigate the community structure of endophytic fungi in roots of Eleocharis erythropoda Steud. and Populus balsamifera L., two plant species growing spontaneously in petroleum-contaminated basins of a former petrochemical plant. These plant species have been selected for this study because they were abundant at the sampling site and can be found in three adjacent basins where concentrations of petroleum hydrocarbon pollutants varied from no-contamination to low or extreme contamination. For both species, three individual plants were sampled from each basin. We characterized endophytic fungal communities associated with the roots of the two species using 454 sequencing of the ITS region and evaluated the changes in community composition in response to contaminants concentrations. Fungi were also isolated and culture-dependent community structures were compared to those obtained with amplicon sequencing.

Materials and Methods

Site Description and Sampling

The sampling site, a former petrochemical plant located near Montreal (Varennes, 45°43 N, 73°22 W, Montreal, QC, Canada), consists of three sedimentation basins separated by dikes in which petroleum hydrocarbon wastes were discarded for many decades. A previous study conducted in one of these three basins has shown that the site was characterized by patchy distribution of spontaneous vegetation mainly dominated by Eleocharis obtusa Willd. and Panicum capillare L. (Desjardins et al., 2014). In this study, we chose to sample the roots of E. erythropoda Steud. and P. balsamifera L., because they were found in the three basins although they were not the dominant plant species. We collected three replicates of each plant species in each of the three basins, i.e., 18 root systems sampled. Individual plants were collected with the soil surrounding the root system. For each individual plant, the complete root system was cut and rhizospheric soil was removed and kept. Roots of each individual plant were cleaned and cut in small pieces that were used for DNA extraction and for isolation of endophytes. Composite rhizospheric soil samples were also prepared for each plant species and each basin. Chemical analyses for assessing petroleum hydrocarbon contaminant concentrations were performed on composite samples for each basin by mixing rhizospheric soil of E. erythropoda and P. balsamifera individuals collected within the same basin resulting in three composite samples for each plant species.

Soil Contamination Analysis

Soil analyses were performed by a commercial service provided by Maxxam Analytics (Montreal, QC, Canada). We measured polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons C10–C50 and polychlorinated biphenyl (PCB). Soil analysis data are shown in Table Table11. Analysis revealed concentrations of 3000, 41000, and 91000 mg[center dot]kg-1 of total petroleum hydrocarbon (TPH) in soil where the three individual plants of E. erythropoda were sampled, and concentrations of 280, 380, and 9040 mg[center dot]kg-1 of TPH where the three individual pants of P. balsamifera were sampled.

Table 1
Concentrations of petroleum hydrocarbon pollutants (polycyclic aromatic hydrocarbons (PAHs) (Chaudhry et al., 2005), alkanes [C10–C50] and total polychlorinated biphenyls [PBCs]) in the sediments of the three basins where Eleocharis erythropoda ...

Endophytic Fungal Isolation Procedure

Eleocharis erythropoda and P. balsamifera root samples were used to isolate fungal endophytes. Root samples were processed immediately after their collection from the site. Roots were cleaned using tap water in order to remove soil particles and were put in 50 mL Falcon tubes. After cleaning, roots were surface-sterilized using the following procedure: roots were first washed three times in sterile water, followed by immersion in 3% H2O2 solution for 4 min, 6% NaOCl solution for 7.5 min, 96% ethanol for 2.5 min and 3% H2O2 solution for 3.5 min. Roots were rinsed three times between treatments using autoclaved water and were rinsed again five times at the end of the disinfection procedure.

The isolation media used were 0.1 × strength PDA medium (Difco), 0.1 × Malt medium (Difco) and Gel-Gro® medium containing 0.4% Gel-Gro® (ICN Biochemicals Aurora, OH, USA) and 0.03% MgSO4 (Silvani et al., 2008). Streptomycin (100 μg[center dot]μL-1) was added to each medium to limit bacterial growth. The efficiency of surface sterilization was assessed by following up the growth for potential contamination of microorganisms in our media using both last rinsing water and root-imprints. A 150 μL aliquot of the last rinsing water was plated on each medium in triplicate and incubated at 27°C for 4 weeks to assess for the presence of living microorganisms. Root-imprints on plates containing each media type were also used to test for the presence of microorganisms living on the surface of the roots. To isolate endophytes, roots were cut in 5 mm segments, and for each plant sample from each contaminated basin, 10 root segments were plated, with 5 segments per Petri dish on each of the three isolation media. For the Gel-Gro® medium, root segments were placed inside a 150 μL drop of medium as described in Silvani et al. (2008). Plates were incubated at 27°C and checked after 5 days for fungal growth emerging from the root fragment ends, under a dissecting microscope. Fragments showing growing hyphae were transferred on 0.1% of Malt medium and incubated at 27°C with a single fragment per plate as described by Silvani et al. (2008). Fungal isolates were further subcultured in triplicates using the same culture media used for isolation. DNA extraction, PCR amplification, and sequencing of isolated fungal endophytes.

Pure fungal isolates were grown on liquid 0.1% of Malt medium and used for DNA extraction, which was performed using the E.Z.N.A.® Forensic DNA Kit (Omega Bio-Tek, Norcross, GA, USA) according the manufacturer’s instructions with slight modifications: fungal mycelia were first dried out on clean paper towel to remove liquid medium, heated at 50°C in 1.5 mL tubes and then crushed with liquid nitrogen using a micro pestle. No DTT buffer was used in the DNA extraction, the incubation time at step 5 of the protocol was changed to 60 min, and the final elution volume used was 70 μL.

PCR amplification of ITS regions was performed using ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) primers (White et al., 1990). PCR reactions were performed in 50 μL volume containing 2 μL of DNA template, 1 × PCR Buffer, 0.4 μM of each primer, 0.25 mM of dNTPs, 1 μg[center dot]μL-1 of bovine serum albumin, 2.5 mM of MgCl2, 1 U of HotStar Taq polymerase (QIAGEN, Toronto, ON, Canada) or 0.75 U of Kapa Taq polymerase (VWR, Montreal, QC, Canada). Cycling conditions were as follows: 15 min for HotStar Taq or 2 min for Kapa Taq at 95°C, 32 cycles of 94°C for 30 s, 55°C for 1 min and 72°C for 1 min, followed by a final elongation step of 72°C for 10 min. PCR reactions were run on a MaterCycler Pro S thermocycler (Eppendorf, Mississauga, ON, Canada). PCR patterns were revealed with agarose gel electrophoresis using GelRed (Biotium, Hayward, CA, USA). Positive PCR reactions were directly sequenced at the Genome Quebec Innovation Centre (McGill University, Montreal, QC). Geneious Pro v.8.1.8 (Drummond et al., 2009) was used to analyze sequences. Blast nucleotide searches were performed at the NCBI web site1 Nucleotide sequences of ITS were deposited in GenBank under the accession numbers: KU179236-KU179276.

454-Amplicon Sequencing

DNA of surface sterilized roots was extracted using the DNeasy Plant Mini Kit (QIAGEN, Toronto, ON, Canada) according the manufacturer’s instructions. Two hundred micro gram of roots were crushed using a mortar and pestle in liquid nitrogen. The complete ITS region was amplified using the primers ITS1F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) (Gardes and Bruns, 1993) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) (White et al., 1990). Forward primers for each sample were tagged with unique multiplex identifiers (MIDs). PCR amplifications were done in a volume of 20 μL containing 0.275 μM of each primer, 415 ng[center dot]μL-1 of bovine serum albumin, 0.83 μL of Tween 20 1% solution (v/v), 0.83 μL of dimethyl sulfoxide 0.275 mM of dNTPs, final concentration of 2.75 mM of MgCl2, 1 × of Taq polymerase buffer, 0.5 U of Taq polymerase (QIAGEN, Toronto, ON, Canada) and 2 μL of genomic DNA as a template. Amplicons were purified using the QIAquick PCR Purification Kit (QIAGEN, Toronto, ON, Canada). Triplicate samples were pooled together and DNA concentration was quantified using the Quant-iT PicoGreen dsDNA assay kit (Life Technologies, Burlington, ON, Canada). All samples were then pooled together in an equimolar ratio and were sent to sequencing at the Genome Quebec Innovation Centre (McGill University, Montréal, QC, Canada) using the 454 GS-FLX Titanium Lib-L chemistry (Roche, Laval, QC, Canada).

Sequence Analysis

Sequences were pre-processed in Mothur v1.29.2 (Schloss et al., 2009) and then imported in Qiime (Caporaso et al., 2010) for quality filtering. Low quality ends were trimmed using a moving window of 50 bp with a minimum quality score of 25. Reads shorter than 200 bp, longer than 1000 bp, with more than two ambiguities (Ns), with homopolymers longer than 8 bp, or with two or more mismatches in the primer or barcode region were eliminated. Reads were then pruned to a fixed length (400 bp) with Mothur. Usearch v7.0 (Edgar, 2010) was used for chimera checking and clustering into Operational Taxonomic Units (OTUs). Singletons (i.e., OTUs represented by only one read) were eliminated because they are mostly artifacts that can be attributed to sequencing errors (Tedersoo et al., 2010). Finally, reads were mapped back into an OTU table and all global singletons (i.e., OTUs represented by a single sequence in the complete dataset) were eliminated. Taxonomy was assigned to each OTU consensus sequence using the UNITE database (Kõljalg et al., 2013) in Mothur with a minimum bootstrap value of 60%. Raw sequences were deposited in the NCBI Sequence Read Archive and are available under the project number PRJNA272728.

Statistical Analysis

All statistical analyses were conducted in R v.3.0.22 (R Foundation for Statistical Computing). In order to find patterns in fungal endophytes diversity, a rarefaction analysis was computed with the ‘iNEXT’ package (Hsieh et al., 2013). Samples were rarefied at 1000 reads in order to compare richness. We also evaluated total richness with the Chao estimator (Chao, 1984). Sample coverage was computed as suggested by Chao and Jost (2012) for both rarefied and complete samples. ANOVA and Tukey HSD post hoc tests were done using the ‘aov’ and ‘TukeyHSD’ functions of the ‘stats’ package to test for significant differences in richness between contamination level and host plant species. A Principal coordinate analysis (PCoA) was used to explore the patterns of fungal OTU composition across TPH levels and host species. The Hellinger distance was computed using the ‘decostand’ and ‘vegdist’ function in ‘vegan’. This distance metric puts more emphasis on species proportions than on abundances (Ramette, 2007; Anderson et al., 2011) and is therefore appropriate in this case (number of reads not a direct measure of OTU abundance in the environment). A Permutational Multivariate Analysis of Variance was performed with the ‘adonis’ function in ‘vegan’ with 9999 permutations to check the significance of site grouping according to host species and contamination levels (plant and TPH effects). Then we performed post-hoc pairwise comparisons (for the different hosts within the root samples and within the soil samples only) and corrected the p-values for multiple comparisons with in Holm correction using the ‘p.adjust’ function in the ‘stats’ package. Finally, we described endophytic community variations in terms of taxonomy by computing the main number of reads for each major fungal classes.


Culture-Based Isolation of Fungal Endophytes from Sterilized-Plant Roots

A total of 41 endophytic fungal isolates were recovered from E. erythropoda (11 isolates) and P. balsamifera (30 isolates) roots. Nine endophytes were isolated from E. erythropoda roots sampled in soils containing 3000 and 40000 mg[center dot]kg-1 TPH while only two strains were isolated from E. erythropoda roots sampled in the highest contamination level (91000 mg[center dot]kg-1 TPH). Twenty-three fungal strains were recovered from P. balsamifera roots sampled in soil containing 280–380 mg[center dot]kg-1 TPH and seven from the 9040 mg[center dot]kg-1 TPH concentration (Table Table22). Sequencing of the ITS region revealed that Fusarium and Leptosphaeria were the most represented genera with five isolates belonging to each genus, which were both isolated from P. balsamifera roots. All the isolated endophytes belonged to Dikarya (Ascomycota and Basidiomycota), with 39 isolates being Ascomycota and only two isolates being Basidiomycota.

Table 2
List of isolated endophytic fungal isolates from surface-sterilized roots of E. erythropoda (E) and P. balsamifera (P).

Endophyte Fungal Diversity and Community Composition

A total of 217,086 sequences were retrieved from pyrosequencing of the ITS region after quality filtering, which clustered into a total of 592 OTUs. The number of sequences per sample ranged from 647 to 9586, with an observed richness ranging from 10 to 44 OTUs (Table Table33). All rarefaction curves were close to saturation and Good’s coverage was close to 1 (Figure Figure11; Table Table33), indicating that sequencing depth was appropriate. ANOVA did not reveal significant differences in richness between contamination levels for either E. erythropoda or P. balsamifera in rarefied richness or Chao total richness estimators (Figure Figure11).

Table 3
Mean and standard error of the number of reads, richness and coverage of fungal OTUs recovered by 454 pyrosequencing.
Operational Taxonomic Unit richness associated with (A)Eleocharis erythropoda and (B)Populus balsamifera roots. The left hand portion of the figure shows rarefaction curves, solid lines representing interpolation and dashed lines total richness extrapolation ...

Principal coordinate analysis shows that the effect of the host identity on endophytic community structure was larger than contamination levels (Figure Figure22). PERMANOVA indicated a significant segregation of sites according to host species (R2 = 0.38994, p = 0.0001) and to contamination levels (R2 = 0.66753, p = 0.0001). However, the pairwise comparison did not support a significant clustering with respect to contamination levels.

Principal coordinate analysis of fungal endophytes of E. erythropoda and P. balsamifera roots. The fungal communities were binned at phylotype resolution and were ordinated using PCoA based on Hellinger distance. Individual samples were colored based ...

The fungal OTUs retrieved belonged to the phyla Ascomycota, Basidiomycota and Glomeromycota, as well as other unclassified fungal taxa. Endophytic communities appeared to be largely dominated by Ascomycetes in terms of relative read abundance across all contamination levels for both plant species. Taxonomic characterization of endophyte communities displayed distinct patterns in E. erythropoda and P. balsamifera (Figure Figure33). E. erythropoda samples harvested in the more contaminated soils (41000 and 91000 mg[center dot]kg-1 TPH) were dominated by Dothideomycetes, with a large proportion of sequences affiliated with this fungal class (between 44.00 and 97.46%), and most OTUs belonging to the order Pleosporales. Root samples from the lower contamination level of 3000 mg[center dot]kg-1 showed a higher relative abundance of Sordariomycetes (mainly from the order Hypocreales) and of unidentified Ascomycota.

Major classes of endophytic fungi recovered from roots of (A)E. erythropoda and (B)P. balsamifera Bar length is proportional to relative number of 454 reads.

On the other hand, P. balsamifera samples did not show significant shifts in the taxonomic composition, with unclassified Ascomycota largely dominating all samples. Leotiomycetes accounted between 10 and 20% of the total OTUs of the 280 of TPH samples. While Dothideomycetes dominated in E. erythropoda roots sampled in highly contaminated basins, unclassified Ascomycetes were dominant in P. balsamifera roots.

We compared OTUs shared between samples of each plant species (Figure Figure44). We found 9.9% of OTUs to be shared among plants growing in the different contamination levels in E. erythropoda while this value increased to 28.2% in P. balsamifera The 41000 mg[center dot]kg-1 TPH contamination level showed the highest percentage of unique OTUs in E. erythropoda, while the two other contamination levels (3000 and 91000 mg[center dot]kg-1 TPH) showed similar, but lower percentages of unique OTUs. However, the community composition was completely different between these contamination levels, with only 1.2% of OTUs shared. We also compared fungal isolates of culture-dependent community structure with amplicon sequencing datasets. We found that 454-sequencing did not reveal the same predominant genera of isolated fungal endophytes.

Venn diagram showing the distribution of fungal OTUs among different concentrations of total petroleum hydrocarbon of E. erythropoda and P. balsamifera. Total petroleum hydrocarbons (TPH) concentration is given in mg/Kg.


The community structure of soil microorganisms is influenced by vegetation and environmental conditions such as contamination levels (Ferrari et al., 2011; Bulgarelli et al., 2012; Abbai and Pillay, 2013; Bell et al., 2014; Chen et al., 2014). We investigated the community structure of endophytic fungi in two spontaneous plant species in relation to contamination levels using both high-throughput sequencing and isolation methods. Our results showed that contamination did not lead to significant changes in endophyte diversity measured with the Chao estimator for either plant species. Furthermore PCoA did not show significant variations in endophytic fungal community composition following changes in contamination levels. However, we could observe taxonomic shifts between E. erythropoda and P. balsamifera where E. erythropoda showed a greater abundance of Dothideomycetes in the 41000 and 91000 TPH samples which was not seen in the 3000 TPH samples.

While previous studies generally showed a negative correlation between the contamination levels and soil microbial diversity (Törneman et al., 2008; Ferrari et al., 2011; Bell et al., 2014), our results clearly showed that root endophytic fungi were not similarly affected. The cryptic nature of endophytic fungi, which live inside their plant host roots, compared to soil microorganisms may explain the differences in these results. Indeed, endophytic fungi may be more protected from environmental stresses and thus soil microorganisms might be more affected by stresses such as petroleum contamination. Furthermore, our study investigated endophytic fungal communities associated with spontaneous vegetation while Bell et al. (2014) investigated changes of introduced plant species. Plant species introduction can change soil biochemical properties, such as soil nutrient availability via root exudation, leading to the promotion of certain microorganisms in the rhizosphere (Van Breemen and Finzi, 1998; Reynolds et al., 2003; Bais et al., 2006).

The abundance of Dothideomycetes in E. erythropoda roots growing in high levels of petroleum hydrocarbon corroborates findings of Ferrari et al. (2011) and Bell et al. (2014) who found this group to be dominant in high hydrocarbon conditions. The genus Alternaria accounted for the majority of the Dothideomycetes sequences found in E. erythropoda roots. While generally described as a plant pathogen (Thomma, 2003), some Alternaria species do not exhibit pathogenic behavior and might help plants to deter pathogens through potential antimicrobial activity (Gu, 2009).

Our results suggest that the microbial community structure detected is highly dependent on the approach used. Isolation methods using different types of media and 454-sequencing did not reveal the same predominant genera. Similarly, Stefani et al. (2015) found distinct bacterial and fungal dominant taxa in hydrocarbon-contaminated soils with culture-dependent and culture-independent methods. Roche 454 sequencing offers a powerful molecular toolkit for studying microbial community structures in pooled samples. The main benefit of this approach is that it allows to process large number of samples and provides a large number of sequencing reads with relatively low cost. However, it has its pitfalls as reported by Tedersoo et al. (2010) where the authors compared 454-sequencing approach and the traditional method. Both methods yielded qualitatively similar results, but there were significant differences that affected the taxonomic view of the fungal community because of PCR biases and technical errors, therefore caution should be used for interpretation of results.


The results of our study show that in some cases endophytic fungal community structure in plants is shaped by concentration of petroleum-hydrocarbon contamination. It also appears that culture-dependent and -independent methods do not promote the same strains, as the most abundant reads in 454 sequencing do not match with the most species recovered using isolation methods from samples taken from the same contamination levels. It is therefore important to better understand how plants recruit endophytic fungi and what their role is when inside the plant. This could potentially help develop newer and better strategies to improve phytoremediation using endophytic fungi. Improvement of such techniques could result in a better colonization of the soil by the plants with the help of endophytic fungi.

Author Contributions

GB: Performed the experiments and wrote the paper. AR-B: Analyzed the data and wrote the paper. MS-A: Helped to design the experiments and wrote the paper. MH: Supervised, conceived, designed the experiments, and wrote the paper.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


We thank Petromont Inc. (ConocoPhillips Canada) for allowing us to access to the Varennes field site. We also thank Karen Fisher-Favret for comments on the manuscript.

Funding. This project was supported by the Genome Canada and Genome Quebec funded Genorem project which are greatly acknowledged.


  • Abbai N. S., Pillay B. (2013). Analysis of hydrocarbon-contaminated groundwater metagenomes as revealed by high-throughput sequencing. Mol. Biotechnol. 54 900–912. 10.1007/s12033-012-9639-z [PubMed] [Cross Ref]
  • Anderson M. J., Crist T. O., Chase J. M., Vellend M., Inouye B. D., Freestone A. L., et al. (2011). Navigating the multiple meanings of β diversity: a roadmap for the practicing ecologist. Ecol. Lett. 14 19–28. 10.1111/j.1461-0248.2010.01552.x [PubMed] [Cross Ref]
  • Bacon C. W., Porter J. K., Robbins J. D., Luttrell E. S. (1977). Epichloë typhina from toxic tall fescue grasses. Appl. Environ. Microbiol. 34 576–581. [PMC free article] [PubMed]
  • Bais H. P., Weir T. L., Perry L. G., Gilroy S., Vivanco J. M. (2006). The role of root exudates in rhizosphere interations with plants and other organisms. Annu. Rev. Plant Biol. 57 233–266. [PubMed]
  • Bell T. H., El-Din Hassan S., Lauron-Moreau A., Al-Otaibi F., Hijri M., Yergeau E., et al. (2014). Linkage between bacterial and fungal rhizosphere communities in hydrocarbon-contaminated soils is related to plant phylogeny. ISME J. 8 331–343. 10.1038/ismej.2013.149 [PMC free article] [PubMed] [Cross Ref]
  • Bulgarelli D., Rott M., Schlaeppi K., van Themaat E., Ahmadinejad N., Assenza F., et al. (2012). Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature 488 91–95. 10.1038/nature11336 [PubMed] [Cross Ref]
  • Caporaso J. G., Kuczynski J., Stombaugh J., Bittinger K., Bushman F. D., Costello E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7 335–336. 10.1038/nmeth.f.303 [PMC free article] [PubMed] [Cross Ref]
  • Carroll G. (1988). Fungal endophytes in stems and leaves: from latent pathogen to mutualistic symbiont. Ecology 69 2–9. 10.2307/1943154 [Cross Ref]
  • Chao A. (1984). Nonparametric estimation of the number of classes in a population. Scand. J. Stat. 11 265–270. 10.2307/4615964 [Cross Ref]
  • Chao A., Jost L. (2012). Coverage-based rarefaction and extrapolation: standardizing samples by completeness rather than size. Ecology 93 2533–2547. 10.1890/11-1952.1 [PubMed] [Cross Ref]
  • Chaudhry Q., Blom-Zandstra M., Gupta S., Joner E. J. (2005). Utilising the synergy between plants and rhizosphere microorganisms to enhance breakdown of organic pollutants in the environment. Environ. Sci. Pollut. Res. 12 34–48. 10.1065/espr2004.08.213 [PubMed] [Cross Ref]
  • Chen L., Luo S., Chen J., Wan Y., Li X., Liu C., et al. (2014). A comparative analysis of endophytic bacterial communities associated with hyperaccumulators growing in mine soils. Environ. Sci. Pollut. Res. 21 7538–7547. 10.1007/s11356-014-2670-9 [PubMed] [Cross Ref]
  • Cunningham S., Berti W. (1993). Remediation of contaminated soils with green plants: an overview. In Vitro Cell. Dev. Biol. Plant 29 207–212. 10.1007/bf02632036 [Cross Ref]
  • Desjardins D., Nissim W. G., Pitre F. E., Naud A., Labrecque M. (2014). Distribution patterns of spontaneous vegetation and pollution at a former decantation basin in southern Québec, Canada. Ecol. Eng. 64 385–390. 10.1016/j.ecoleng.2014.01.003 [Cross Ref]
  • Drummond A. J., Ashton B., Buxton S., Cheung M., Cooper A., Duran C., et al. (2009). Geneious version 4.8. Available at:
  • Edgar R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26 2460–2461. 10.1093/bioinformatics/btq461 [PubMed] [Cross Ref]
  • Feller I. C. (1995). Effects of nutrient enrichment on growth and herbivory of dwarf red mangrove (Rhizophora mangle). Ecol. Monogr. 65 477–505. 10.2307/2963499 [Cross Ref]
  • Ferrari B. C., Zhang C., van Dorst J. (2011). Recovering greater fungal diversity from pristine and diesel fuel contaminated sub-antarctic soil through cultivation using both a high and a low nutrient media approach. Front. Microbiol. 2:217 10.3389/fmicb.2011.00217 [PMC free article] [PubMed] [Cross Ref]
  • Ganley R. J., Sniezko R. A., Newcombe G. (2008). Endophyte-mediated resistance against white pine blister rust in Pinus monticola. For. Ecol. Manage. 255 2751–2760. 10.1016/j.foreco.2008.01.052 [Cross Ref]
  • Gardes M., Bruns T. D. (1993). ITS primers with enhanced specificity for basidiomycetes – application to the identification of mycorrhizae and rusts. Mol. Ecol. 2 113–118. 10.1111/j.1365-294X.1993.tb00005.x [PubMed] [Cross Ref]
  • Göhre V., Paszkowski U. (2006). Contribution of the arbuscular mycorrhizal symbiosis to heavy metal phytoremediation. Planta 223 1115–1122. 10.1007/s00425-006-0225-0 [PubMed] [Cross Ref]
  • Gu W. (2009). Bioactive metabolites from Alternaria brassicicola ML-P08 an endophytic fungus residing in Malus halliana. World J. Microbiol. Biotechnol. 25 1677–1683. 10.1007/s11274-009-0062-y [Cross Ref]
  • Hamilton C. E., Bauerle T. L. (2012). A new currency for mutualism? Fungal endophytes alter antioxidant activity in hosts responding to drought. Fungal Diver. 54 39–49. 10.1007/s13225-012-0156-y [Cross Ref]
  • Hassan S. E., Boon E., St-Arnaud M., Hijri M. (2011). Molecular biodiversity of arbuscular mycorrhizal fungi in trace metal-polluted soils. Mol. Ecol. 20 3469–3483. 10.1111/j.1365-294X.2011.05142.x [PubMed] [Cross Ref]
  • Hassan S. E., Hijri M., St-Arnaud M. (2013). Effect of arbuscular mycorrhizal fungi on trace metal uptake by sunflower plants grown on cadmium contaminated soil. N. Biotechnol. 30 780–787. 10.1016/j.nbt.2013.07.002 [PubMed] [Cross Ref]
  • Hsieh T. C., Ma K. H., Chao A. (2013). iNEXT online: interpolation and extrapolation (Version 1.0) [Software]. Available at:
  • Hyde K. D., Soytong K. (2008). The fungal endophyte dilemma. Fungal Diver. 33 163–173.
  • Jumpponen A. (2001). Dark septate endophytes – are they mycorrhizal? Mycorrhiza 11 207–211. 10.1007/s005720100112 [Cross Ref]
  • Jumpponen A., Trappe J. M. (1998). Dark septate endophytes: a review of facultative biotrophic root-colonizing fungi. New Phytol. 140 295–310. 10.1046/j.1469-8137.1998.00265.x [Cross Ref]
  • Kaimi E., Mukaidani T., Tamaki M. (2007). Screening of twelve plant species for phytoremediation of petroleum hydrocarbon-contaminated soil. Plant Product. Sci. 10 211–218. 10.1626/pps.10.211 [Cross Ref]
  • Knoth J. L., Kim S. H., Ettl G. J., Doty S. L. (2014). Biological nitrogen fixation and biomass accumulation within poplar clones as a result of inoculations with diazotrophic endophyte consortia. New Phytol. 201 599–609. 10.1111/nph.12536 [PubMed] [Cross Ref]
  • Kobayashi D. Y., Palumbo J. D. (2000). Bacterial endophytes and their effects on plants and uses in agriculture. Microbial Endophytes 19 199–233.
  • Kogel K.-H., Franken P., Hückelhoven R. (2006). Endophyte or parasite – what decides? Curr. Opin. Plant Biol. 9 358–363. 10.1016/j.pbi.2006.05.001 [PubMed] [Cross Ref]
  • Kõljalg U., Nilsson R. H., Abarenkov K., Tedersoo L., Taylor A. F., Bahram M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22 5271–5277. 10.1111/mec.12481 [PubMed] [Cross Ref]
  • Li H.-Y., Wei D.-Q., Shen M., Zhou Z.-P. (2012). Endophytes and their role in phytoremediation. Fungal Diver. 54 11–18. 10.1007/s13225-012-0165-x [Cross Ref]
  • Melo F. M. P. D., Fiore M. F., Moraes L. A. B. D., Silva-Stenico M. E., Scramin S., Teixeira M. D. A., et al. (2009). Antifungal compound produced by the cassava endophyte Bacillus pumilus MAIIIM4a. Sci. Agric. 66 583–592. 10.1590/S0103-90162009000500002 [Cross Ref]
  • Paszkowski U. (2006). Mutualism and parasitism: the yin and yang of plant symbioses. Curr. Opin. Plant Biol. 9 364–370. 10.1016/j.pbi.2006.05.008 [PubMed] [Cross Ref]
  • Peters A. F. (1991). Field and culture studies of Streblonema macrocystis sp. nov. (Ectocapales, Phaeophyceae) from Chile, a sexual endophyte of giant kelp. Phycologia 30 365–377. 10.2216/i0031-8884-30-4-365.1 [Cross Ref]
  • Petrini O. (1991). “Fungal endophytes of tree leaves,” in Microbial Ecology of Leaves, eds Andrews J., Hirano S., editors. (New York, NY: Springer Verlag; ), 179–197.
  • Pilon-Smits E. (2005). Phytoremediation. Annu. Rev. Plant Biol. 56 15–39. [PubMed]
  • Porras-Alfaro A., Bayman P. (2011). “Hidden fungi, emergent properties: endophytes and microbiomes,” in Annual Review of Phytopathology Vol. 49 eds VanAlfen N. K., Bruening G., Leach J. E., editors. (Palo Alto: Annual Reviews; ), 291–315. [PubMed]
  • Ramette A. (2007). Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 62 142–160. 10.1111/j.1574-6941.2007.00375.x [PMC free article] [PubMed] [Cross Ref]
  • Reynolds H. L., Packer A., Bever J. D., Clay K. (2003). Grassroots ecology: plant-microbe-soil interactions as drivers of plant community structure and dynamics. Ecology 84 2281–2291. 10.1890/02-0298 [Cross Ref]
  • Rodriguez R. J., White J. F., Jr., Arnold A. E., Redman R. S. (2009). Fungal endophytes: diversity and functional roles. New Phytol. 182 314–330. 10.1111/j.1469-8137.2009.02773.x [PubMed] [Cross Ref]
  • Rylott E. L. (2014). Endophyte consortia for xenobiotic phytoremediation: the root to success? Plant Soil 385 389–394. 10.1007/s11104-014-2296-1 [Cross Ref]
  • Salt D. E., Smith R. D., Raskin I. (1998). Phytoremediation. Annu. Rev. Plant Physiol. Plant Mol. Biol. 49 643–668. 10.1146/annurev.arplant.49.1.643 [PubMed] [Cross Ref]
  • Schardl C. L., Leuchtmann A., Spiering M. J. (2004). Symbioses of grasses with seedborne fungal endophytes. Annu. Rev. Plant Biol. 55 315–340. 10.1146/annurev.arplant.55.031903.141735 [PubMed] [Cross Ref]
  • Schloss P. D., Westcott S. L., Ryabin T., Hall J. R., Hartmann M., Hollister E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75 7537–7541. 10.1128/AEM.01541-09 [PMC free article] [PubMed] [Cross Ref]
  • Schulz B., Boyle C. (2005). The endophytic continuum. Mycol. Res. 109 661–686. 10.1017/s095375620500273x [PubMed] [Cross Ref]
  • Schulz B., Boyle C., Draeger S., Rommert A. K., Krohn K. (2002). Endophytic fungi: a source of novel biologically active secondary metabolites. Mycol. Res. 106 996–1004. 10.1017/s0953756202006342 [Cross Ref]
  • Schulz B., Wanke U., Draeger S., Aust H. J. (1993). Endophytes from herbaceous plants and shrubs: effectiveness of surface sterilization methods. Mycol. Res. 97 1447–1450. 10.1016/S0953-7562(09)80215-3 [Cross Ref]
  • Schulz B. J. E., Boyle C. J. C., Sieber T. N. (2006). Microbial Root Endophytes. Berlin: Springer.
  • Sieber T. N. (2002). “Fungal root endophytes,” in Plant Roots – The Hidden Half, eds Kafkafi U., Waisel Y., Eshel A., editors. (Boca Raton: CRC Press; ), 887–917.
  • Silvani V. A., Fracchia S., Fernández L., Pérgola M., Godeas A. (2008). A simple method to obtain endophytic microorganisms from field-collected roots. Soil Biol. Biochem. 40 1259–1263. 10.1016/j.soilbio.2007.11.022 [Cross Ref]
  • Soleimani M., Afyuni M., Hajabbasi M. A., Nourbakhsh F., Sabzalian M. R., Christensen J. H. (2010). Phytoremediation of an aged petroleum contaminated soil using endophyte infected and non-infected grasses. Chemosphere 81 1084–1090. 10.1016/j.chemosphere.2010.09.034 [PubMed] [Cross Ref]
  • Stefani F. O., Bell T. H., Marchand C., de la Providencia I. E., El Yassimi A., St-Arnaud M., et al. (2015). Culture-dependent and -independent methods capture different microbial community fractions in hydrocarbon-contaminated soils. PLoS ONE 10:e0128272 10.1371/journal.pone.0128272 [PMC free article] [PubMed] [Cross Ref]
  • Stone J. K., Bacon C. W., White J. F. (2000). “An overview of endophytic microbes: endophytism defined,” in Microbial Endophytes, eds Bacon C. W., White J. F., editors. (New York, NY: Marcel Dekker, Inc.), 3–29.
  • Tan R. X., Zou W. X. (2001). Endophytes: a rich source of functional metabolites. Nat. Prod. Rep. 18 448–459. 10.1039/b100918o [PubMed] [Cross Ref]
  • Tedersoo L., Nilsson R. H., Abarenkov K., Jairus T., Sadam A., Saar I., et al. (2010). 454 Pyrosequencing and Sanger sequencing of tropical mycorrhizal fungi provide similar results but reveal substantial methodological biases. New Phytol. 188 291–301. 10.1111/j.1469-8137.2010.03373.x [PubMed] [Cross Ref]
  • Thomma B. P. (2003). Alternaria spp.: from general saprophyte to specific parasite. Mol. Plant Pathol. 4 225–236. 10.1046/j.1364-3703.2003.00173.x [PubMed] [Cross Ref]
  • Törneman N., Yang X., Baath E., Bengtsson G. (2008). Spatial covariation of microbial community composition and polycyclic aromatic hydrocarbon concentration in a creosote-polluted soil. Environ. Toxicol. Chem. 27 1039–1046. 10.1897/07-440.1 [PubMed] [Cross Ref]
  • Van Breemen N., Finzi A. C. (1998). Plant-soil interactions: ecological aspects and evolutionary implications. Biogeochemistry 42 1–19. 10.1023/a:1005962124317 [Cross Ref]
  • Varma A., Verma S., Sudha Sahay N., Butehorn B., Franken P. (1999). Piriformospora indica, a cultivable plant-growth-promoting root endophyte. Appl. Environ. Microbiol. 65 2741–2744. [PMC free article] [PubMed]
  • White T., Bruns T., Lee S., Taylor J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols: A Guide to Methods and Applications, eds Innis M., Gelfand D., Shinsky J., White T., editors. (Cambridge: Academic Press; ), 315–322.

Articles from Frontiers in Microbiology are provided here courtesy of Frontiers Media SA