Unique arbuscular mycorrhizal fungal communities uncovered in date palm plantations and surrounding desert habitats of Southern Arabia.

The main objective of this study was to shed light on the previously unknown arbuscular mycorrhizal fungal (AMF) communities in Southern Arabia. We explored AMF communities in two date palm (Phoenix dactylifera) plantations and the natural vegetation of their surrounding arid habitats. The plantations were managed traditionally in an oasis and according to conventional guidelines at an experimental station. Based on spore morphotyping, the AMF communities under the date palms appeared to be quite diverse at both plantations and more similar to each other than to the communities under the ruderal plant, Polygala erioptera, growing at the experimental station on the dry strip between the palm trees, and to the communities uncovered under the native vegetation (Zygophyllum hamiense, Salvadora persica, Prosopis cineraria, inter-plant area) of adjacent undisturbed arid habitat. AMF spore abundance and species richness were higher under date palms than under the ruderal and native plants. Sampling in a remote sand dune area under Heliotropium kotschyi yielded only two AMF morphospecies and only after trap culturing. Overall, 25 AMF morphospecies were detected encompassing all study habitats. Eighteen belonged to the genus Glomus including four undescribed species. Glomus sinuosum, a species typically found in undisturbed habitats, was the most frequently occurring morphospecies under the date palms. Using molecular tools, it was also found as a phylogenetic taxon associated with date palm roots. These roots were associated with nine phylogenetic taxa, among them eight from Glomus group A, but the majority could not be assigned to known morphospecies or to environmental sequences in public databases. Some phylogenetic taxa seemed to be site specific. Despite the use of group-specific primers and efficient trapping systems with a bait plant consortium, surprisingly, two of the globally most frequently found species, Glomus intraradices and Glomus mosseae, were not detected neither as phylogenetic taxa in the date palm roots nor as spores under the date palms, the intermediate ruderal plant, or the surrounding natural vegetation. The results highlight the uniqueness of AMF communities inhabiting these diverse habitats exposed to the harsh climatic conditions of Southern Arabia.

Introduction

Geological history and the influence of past climates of Southern Arabia (Preusser et al. 2002) have played a major role in shaping its unique environmental settings. Today, most of the area is hyper-arid (Fisher and Membery 1998), and half of the region has maximum summer temperatures exceeding 48°C (Glennie and Singhvi 2002). Agricultural ecosystems inherited infertile and light structured span class="Chemical">oils, soil salinity, drought (Ministry of Agriculture and Fisheries 1989, 1991), and high phosphorus (P) fixation (Cookson 1996). In the harsh environmental conditions of Southern Arabia, it can be assumed that the mycorrhizal symbiosis between plants and arbuscular mycorrhizal fungi (AMF) is a key factor for plant tolerance to such conditions. Moreover, native AMF communities can be expected to have distinctive physiologies adapted to withstand these conditions. As a multifunctional symbiosis (Newsham et al. 1995), the mycorrhiza is well known to play a crucial role in the P acquisition by plants (Smith and Gianinazzi-Pearson 1988). Through biochemical, biological, and biophysical mechanisms, AMF enhances soil aggregate stability (Rillig and Mummey 2006), a feature of particular relevance for the sandy soils prone to erosion. Under arid conditions, mycorrhizal plants were found also to maintain higher drought tolerance (Augé 2001) and to have better access to P than non-mycorrhizal ones (Neumann and George 2004).

The main crop in Southern Arabia is pan class="Chemical">pan class="Species">date palmpan> (papan>n class="Species">Phoenix dactylifera) (FAO 2009), the possibly most ancient crop in the world (Zaid and de Wet 1999). In addition to its high cultural and economical value, date palm still represents a survival crop (Al-Shahib and Marshall 2003) for the region and also for millions of people in many desert regions of the world. The decline of the world date production has reached an alarming level, and more research and global cooperation were called (FAO 2004). Date palms have already been shown to form AM symbiosis (Sabet 1940; Dreyer et al. 2006). Under greenhouse conditions, AMF promoted the growth of date palm seedlings, especially on nutrient-poor soils (Al-Whaibi and Khaliel 1994). AMF were also proposed as a potential biocontrol agent for bayoud disease (Jaiti et al. 2007; Jaiti et al. 2008), one of the most aggressive diseases threatening date palm (Zaid et al. 1999). It seems highly likely, therefore, that this symbiosis represents a major factor that needs to be considered in the effort to sustain productivity and fitness of date palms. However, a fundamental prerequisite for any management or promotion of this symbiosis is a basic understanding of the occurrence and diversity of AMF in the target ecosystems, because AMF diversity has significant ecological consequences: different members of the AMF community differ in their potential to promote plant fitness and to respond to environmental biotic and abiotic cues (Fitter 2005).

Such information on diversity of AMF associated with pan class="Species">date palm in Southern Arabia or other parts of the world is scarce. To our knowledge, the only study on this topic was conducted in North Africa (Bouamri et al. 2006) and was based exclusively on the classical identification method by spore morphology. In contrast, molecular methods have the potential to detect and identify AMF directly in colonized roots. In particular, by comparing sequences from public data bases, taxa detected in one environment can be compared with those found in other ecosystems (Öpik et al. 2006).

AMF communities associated with pan class="Chemical">pan class="Species">date palmspan> must have been shaped through a complex interaction between two main factors: firstly, the original local environmental settings, including the indigenous AMF community associated with the native plants prior to the introduction of papan>n class="Species">date palms and, second, the impact of agricultural practices including the introduction of a new plant, physical span class="Chemical">oil disturbance, irrigation, and addition of chemical and organic fertilizers. It also includes the probability of accidental introduction of exogenous AMF through the agricultural inputs and machinery. AMF communities in agro-ecosystems were often found to differ from the communities in the surrounding natural habitats. While some agro-ecosystems in a hot and arid area were found to have higher AMF species richness and spore abundance (Li et al. 2007), others in temperate climate zones were found to have lower species richness and infectivity (Helgason et al. 1998; Oehl et al. 2003) compared to surrounding more natural habitats.

Our goal was to analyze AMF communities associated with pan class="Species">date palm in Southern Arabia and to compare them with the communities associated with the local natural vegetation. Using spore morphotyping, we studied two differently managed pan class="Species">date palm plantations at an oasis and an experimental station, a disturbed and dry area where a ruderal plant was growing at the experimental station, undisturbed habitat adjacent to the experimental station, and a distant natural habitat hosting a desert plant, located in Al-Sharqiya Sands. Additionally, molecular identification methods were applied to investigate the AMF communities colonizing the pan class="Species">date palm roots from the two plantations.

Considering the characteristic environmental settings in the harsh climate of Southern Arabia, we hypothesized that pan class="Chemical">pan class="Species">date palmspan>, ruderal, and natural vegetation are associated with unique AMF communities.

Materials and methods

Study habitats and vegetation

All study habitats (Fig. 1) are located within the territory of the Sultanate of Oman. These habitats are as follows:

Fig. 1

Simplified map of Southern Arabia showing main geographical features of the area (modified from Preusser et al. 2002) including the sampling sites. The four sites are an oasis (filled triangle), an experimental station (empty triangle), an undisturbed habitat (filled circle), and sand dunes (empty circle)

Date palm plantations

The two pan class="Chemical">pan class="Species">date palmpan> plantations were established in a sandy plain. The first was a farm located on the margin of the oasis of Al-Kamel (22°12′56″N, 59°12′9″E). It was cultivated mainly with papan>n class="Species">date palms belonging to different cultivars planted randomly, among them are Khalas, Barmi, Madluki, Helali, Naghal, and Khunaizi. Most of these cultivars originated from the same oasis, but some were brought from areas in the interior of Oman around 300 km northwest. Soil traces from the original land usually are attached to the roots of the imported seedlings. Management followed local traditional farming practices since the establishment of the farm in 1992. The only source of fertilization was the manure produced on the same farm. Irrigation was performed by the traditional flooding system with water flowing from a natural source a few kilometers away through narrow channels called Aflaj (Al-Marshudi 2001; UNESCO 2006). The area between the palm trees was annually tilled using a tractor and then cultivated with Sorghum bicolor. To our knowledge, no synthetic pesticides or chemical fertilizers were used on this farm.

The second plantation was an experimental station (22°14′13″N, 59°11′5″E) belonging to the Ministry of Agriculture of Oman. It is situated 3 km northwest of the oasis, and it was cultivated mainly with pan class="Chemical">pan class="Species">date palmspan> (Fig. 2a). The papan>n class="Species">date palms, all of the same cultivar Khalas, were originally produced by the tissue culture (micropropagation) technique in the tissue culture laboratory of the Ministry of Agriculture of Oman located about 300 km northwest of the study site. Upon transfer to an unsterilized mixture of vermiculite and peat moss, the seedlings passed a period of weaning in a shadow house at the same laboratory. Then they were transferred to the plantation. The management followed conventional guidelines since the establishment of the station in 1990. Synthetic pesticides were used rarely, and chemical fertilizers (mainly NPK) were added according to recommended ecologically sound agricultural practices. A bubbler irrigation system was applied. All the agricultural inputs and managements, including irrigation, were directed to the root zone of the palm trees exclusively, leaving uncultivated and dry inter-cropping strip in these plantations which had lost its original native vegetation, apparently because of the mechanical disturbance but was recolonized by a ruderal plant, pan class="Species">Polygala erioptera (Fig. 2b).

Fig. 2

Plants studied in the different habitats of Southern Arabia: a date palm (Phoenix dactylifera), b P. erioptera, c Z. hamiense, d S. persica, e P. cineraria, and f H. kotschyi

Plants studied in the different habitats of Southern Arabia: a pan class="Chemical">pan class="Species">date palmpan> (papan>n class="Species">Phoenix dactylifera), b P. erioptera, c Z. hamiense, d S. persica, e P. cineraria, and f H. kotschyi

Undisturbed habitat

The third sampled site was an undisturbed habitat adjacent to the pan class="Chemical">pan class="Species">date palmpan> plantation at the experimental station (22°14′11″N, 59°10′53″E). The natural vegetation of this habitat was also representative for the area surrounding the papan>n class="Species">date palm plantation at the oasis. The vegetation consisted mainly of three perennial plant species, namely Zygophyllum hamiense (Fig. 2c), Salvadora persica (Fig. 2d), and Prosopis cineraria (Fig. 2e). The inter-plant area had traces of desiccated grass which seems to grow after rainfall and dry up during the dry season. According to the indigenous knowledge, the areas of both date palm plantations previously looked like this undisturbed natural habitat before agricultural use started.

Sand dunes

Lastly, a different natural habitat, about 50 km southwest of the three main study habitats, was chosen (21°52′39″N, 58°52′44″E). The location was at the bottom of an extinct Holocene lake (Radies et al. 2005) in the current sand sea (Preusser et al. 2002; Radies et al. 2004) that was previously called Wahiba Sands and presently Al-Sharqiya Sands. A native plant, pan class="Chemical">pan class="Species">Heliotropium kotschyipan>, occurring in this habitat (Fig. 2f) was included in our study.

Soil and root sampling

Span class="Chemical">pan class="Chemical">oilpan> samples were collected in August 2006 from each of the four study habitats. Four replicate plots of ca. 200 m2 were randomly chosen at each habitat. From each plot, four plants of each species were selected for obtaining a pooled sample per plant species. In the case of the papan>n class="Species">date palms at the two plantations and the three perennial plant species in the undisturbed habitat, the rhizosphere soil was sampled by first digging a pit of 30-cm diameter and 30-cm depth, close to the plant stems to have access to the roots. Soil was then collected from the border of the pits in order to obtain a sample representing a vertical cross section across the root zones of these plants. Moreover, soil samples were taken below the ruderal plant P. erioptera at the experimental station, at the inter-plant area in the undisturbed habitat, and below H. kotschyi in the sand dunes by sampling the sandy soil around their root systems. This sampling procedure was chosen because the shallow root system did not allow making the soil cross-section. A total of 32 pooled soil samples were obtained as a result of 128 original samples.

For molecular analysis of AMF communities, roots of pan class="Chemical">pan class="Species">date palmspan> were collected immediately from the sampled spapan>n class="Chemical">oil. They were washed with tap water and kept in centrifuge tubes in a cooling box for transportation. Four pooled samples were obtained from each date palm plantation. DNA was extracted from the roots as soon as they arrived to the laboratory.

In a preliminary experiment, root samples from five different pan class="Chemical">pan class="Species">date palmpan> trees were collected in January 2004 at the experimental station.

Properties of span class="Chemical">pan class="Chemical">oilpan>s from the study habitats are presented in Table 1. Analyzed spapan>n class="Chemical">oils represent pooled samples from the root zones of all plants in each habitat in exception of the experimental station where it represents only the pan class="Species">date palms root zones.

Table 1

Soils chemical and physical properties in the four study habitats

	Pa (mg kg−1)	% N	pH		% Organic matter	% Sand	% Silt	% Clay	Soil texture
			H2O	KCl
Oasis	71	0.18	8.4	7.9	3.7	67.2	24.8	8.1	Sandy loam
Exp. station	93	0.1	8.5	7.9	2.6	60.5	33.6	5.9	Sandy loam
Undisturbed habitat	41.4	<0.1	8.2	8.1	1.5	70.7	25.4	4.3	Sandy loam
Sand dunes	22.6	<0.1	9.1	8.5	0.45	90.7	9.3	0	Sand

aSodium bicarbonate-soluble P, analyzed using the method given by Chapman and Pratt (1961)

Span class="Chemical">pan class="Chemical">oilpan>s chemical and physical properties in the four study habitats

pan class="Chemical">pan class="Chemical">aSodium bicarbonatepan>-soluble P, analyzed using the method given by Chapman and Pratt (1961)

Trap culturing the AMF

For each span class="Chemical">pan class="Chemical">oilpan> sample, two trap culture pots (0.85 L) were established. Pots were filled with 0.45 kg of an autoclaved substrate consisting of Terragreen (American papan>n class="Chemical">aluminum oxide, oil dry US special, type III R, <0.125 mm; Lobbe Umwelttechnik, Iserlohn, Germany). The mycorrhizal inocula (50 g of soil sample) were spread as a layer on the surface of the substrate in the pots and covered with a thin layer of the substrate. Two trap culture methods were used. The first three trap plants (Plantago lanceolata, Hieracium pilosella, and Allium porrum) were randomly planted in every pot. Using more than one trap plant species was expected to increase the recovery of AMF species from the field soils, since different plants might trap different AMF. In the other trap culture system, S. bicolor was used as trap plant. Five 2-week-old AMF-free plantlets of S. bicolor were randomly planted in every pot. The trap cultures were kept in a greenhouse (temperature, 19–23°C; relative humidity, 50%) for 8 months and irrigated with an automatic irrigation system (Tropf-Blumat; Weninger GmbH, Telfs, Austria). The roots were used to assess the AMF colonization levels after 7 months of plant growth as an indicator for the infection potential.

Evaluation of AMF infection potential

The mycorrhizal colonization on pan class="Chemical">pan class="Species">S. bicolorpan> was assessed on 30 root pieces (1-cm segments) from each replication. The root segments were incubated overnight at room temperature in 0.1% cotton blue (w/v) in papan>n class="Chemical">lactic acid, destained in pan class="Chemical">glycerol, and observed under a light microscope for quantification. Infection potential was defined as the percentage of root length colonized by AMF according to Trouvelot et al. (1986).

AMF spore isolation and identification

AMF spores occurring in each field soil sample and those propagated in the trap cultures were extracted by wet sieving and pan class="Chemical">sucrose density gradient centrifugation (Daniels and Skipper 1982). Fifteen grams of air-dried field sn class="Chemical">oil or 30 cm3 of harvested trap culture substrate were passed through 1000-, 500-, 125-, and 32-µm sieves. The 500-µm sieve was checked for large spores, spore clusters, and sporocarps whereas the contents of the 125- and 32-µm sieves were layered onto a water–sucrose solution (70% w/v) gradient and centrifuged at 900×g for 2 min. The resulting supernatant was washed with tap water for 2 min in a 32-µm sieve and transferred to petri dishes. Spores, spore clusters, and sporocarps obtained from all sieves were mounted on slides with polyvinyl-lactic acid-glycerol (Koske and Tessier 1983) or polyvinyl-lactic acid-glycerol mixed 1:1 (v/v) with Melzer’s reagent (Brundrett et al. 1994) and examined under a light microscope (Zeiss; Axioplan) at a magnification of up to ×400. The presence/absence of spores from all AMF species was determined for each sample. Spore abundance was assessed by measuring the total number of AMF spores extracted from the field soil samples while species richness was derived from the total number of AMF species found in the samples of both the field soils and the trap cultures.

Identification was based on original and recent species descriptions and identification manuals (Schenck and Perez 1990; the International Culture Collection of Arbuscular and Vesicular-Arbuscular Endomycorrhizal Fungi, INVAM: http://invam.caf.wvu.edu; Arbuscular Mycorrhizal Fungi (Glomeromycota), Endogone and Complexipes species deposited in the Department of Plant pan class="Chemical">Pathology, University of Agriculture in Szczecin, Poland: http://www.agro.ar.szczecin.pl/~jblaszkowski/).

DNA extraction and polymerase chain reaction

Five centimeters of each pan class="Chemical">pan class="Species">date palmpan> root sample were first ground in liquid papan>n class="Chemical">nitrogen within a 1.5-ml microcentrifuge tube using a pellet pestle. A DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) was used according to the manufacturer’s instructions to extract the total DNA. A volume of 30 µl of elution buffer was used to elute the DNA. DNA extracts were used as template for the first polymerase chain reaction (PCR) reaction. PCR was performed in a nested procedure as previously described (Redecker 2000; Redecker et al. 2003) using Taq polymerase from Amersham (Basel, Switzerland), 2 mM MgCl2, 0.5 µM primers, and 0.13 mM of each desoxynucleotide. The PCR was performed using the universal eukaryote primers NS5 and ITS4 (White et al. 1990) and carried out as follows: 3 min at 94°C, followed by 30 cycles (45 s at 94°C, 50 s at 51°C, and 1 min 30 s at 72°C) and a final extension phase of 10 min at 72°C. For the second PCR, dilutions of 1:100 of the PCR products in TE buffer were used as a template, and a hot start at 61°C was performed to minimize the possibility of nonspecific amplification. Primer sequences are given in Redecker (2000) and Appoloni et al. (2008). Targeted AMF clades were Glomus group A, Glomus group B, Glomus group C, Acaulosporaceae, Archaeospora/Paraglomus, Gigasporaceae, and Paraglomus. The cycling parameters of the second PCR step were the same as those of the first one except that the annealing temperature was raised to 61°C. The second PCR step products were then run on agarose gels (2%:1% NuSieve/SeaKem, Cambrex Bio Science, Rockland, ME, USA) in Tris/Acetate buffer at 120 V for 25 min. When no second PCR product was observed, the amplification process was repeated using combinations of different dilutions of the template DNA and different dilutions of first PCR step products.

Cloning, RFLP analysis, and sequencing

PCR products were purified using the High Pure Kit from Hoffman LaRoche (Basel, Switzerland) and cloned using the pGEM-T vector (Promega/Catalys, Wallisellen, Switzerland). Thirty-four clones of each PCR product were chosen randomly, and the inserts were re-amplified. Digestion reactions of 3 h were performed using the restriction enzymes HinfI and MboI. Products of the reaction were revealed by pan class="Chemical">pan class="Chemical">agarosepan> gels as described above. Representative clones of restriction types were then re-amplified and purified using the High Pure Kit and sequenced in both directions. The BigDye Terminator Cycle Sequencing Kit (ABI, Foster City, CA, USA) was used for labeling. Samples were run on an ABI 310 capillary sequencer.

Sequence analyses

The glomeromycotan origin of the sequences was initially tested by BLAST search. Using pan class="Chemical">pan class="Chemical">PApan>UP*4b10 (Swofford 2001), sequences from the present study were aligned to sequences in GenBank which showed high similarity to them, to other previously published sequences in addition to yet unpublished sequences from our own database. Phylogenetic trees were obtained by distance analysis (using the neighbor-joining algorithm) in papan>n class="Chemical">PAUP*4b10 using the Kimura two-parameter model and a gamma shape parameter = 0.5. In the resulting phylogenetic trees, sequence phylotypes were defined in a conservative manner as consistently separated monophyletic groups. The sequence phylotypes were named based on the numerical system used to define groups detected by the group of Basel.

Statistical analysis

Significance of differences in AMF spore abundance, species richness, and infection potential between the samples coming from and within different habitats was tested using Fisher's least significant difference at P < 0.05 after one-way analysis of variance. The correlation between the AMF spore abundance and species pan class="Chemical">pan class="Disease">richness and infectionpan> potential was calculated through a simple regression test. A dendrogram was obtained by cluster analysis displaying the similarity of AMF species composition across the host plants in all habitats. The group average clustering method and the squared Euclidean distance metric were used. All the statistical tests were performed using Statgraphics (version 3.1).

Results

AMF morphospecies composition

Altogether, 25 AMF taxa were detected by spore morphotyping in the 32 span class="Chemical">pan class="Chemical">oilpan> samples originating from the two papan>n class="Species">date palm plantations, the disturbed inter-cropping area in the experimental station, the undisturbed habitat, and the sand dunes (Table 2). Out of the detected morphospecies, 18 species (72%) belonged to the genus Glomus, 2 species to each Scutellospora (8%) and Racocetra (8%), and only 1 species (4%) to each of the following genera: Acaulospora, Paraglomus, and Ambispora. Four of the Glomus species (Glomus sp. OMA3, Glomus sp. OMA5, Glomus sp. OMA8, and Glomus sp. OMA9) are undescribed. The most frequent species in the field samples was putatively assigned to Glomus aggregatum, with 26 occurrences in the 32 samples while the least frequent were Glomus microaggregatum, Racocetra gregaria, Acaulospora spinosa, Ambispora gerdemannii, Glomus sp. OMA7, and Glomus sp. OMA8, each with one rare occurrence. However, the most frequent species in the trap cultures was Glomus eburneum with 29 occurrences in the 32 trap culture pots. Most of the species detected in the field samples did not sporulate in the trap cultures.

Table 2

AMF species and their distribution among the different habitats and plant species studied

	Oasis	Exp. station		Undisturbed habitat				Sand dunes
	Date palm	Date palm	P. erioptera	Z. hamiense	S. persica	P. cineraria	Inter-plant	H. kotschyi
Glomus sinuosum	A, B, C, D	A, B, C, D
Glomus sp. OMA6	a, b, c, d	a, b, d
Scutellospora calospora	a, b, d	a, b, c, d
Scutellospora sp. OMA10	a, b, c	a
Glomus microaggregatum		b
Glomus sp. OMA9		a, b, c, d
Racocetra fulgida	a, c
Racocetra gregaria	a
Acaulospora spinosa	b
Glomus etunicatum	ca
Glomus sp. OMA13		ba
Glomus sp. OMA2	a (rare), da	a (rare), ca				da
Glomus constrictum	aa, ba, ca	a, B, c, d	d (rare)
		ca, da
Glomus aggregatum	a, b, c, d	a, b, c, d	a, b, c, d	a, b, c, d	b, c, d	a, b, c, d,	A, b, d	ca
	aa, ba, da	aa, ca				aa, ba, ca, da	aa
Glomus sp. OMA3	a, b, c, d	b, c, d	a, b, c, d	c, d	c, d	a, b, c, d	a, c, d
Glomus eburneum	aa, ba, ca, da	ba, ca, da	aa, ba, da	aa, ba, ca, da	b	A, b, D,	Aa, ba, ca, da	aa, ba, ca
					aa, ba, ca, da	aa, a, ca, da
Glomus microcarpum		c (rare)	d (rare)	a, b, d	d	c	a, d
Paraglomus occultum	aa, da					a, b
Glomus sp. OMA5			d	b, c, d	a, b	a	D
Glomus macrocarpum			a, b, c, d		a, b, c, d	c, d	a, b, c, d
Ambispora gerdemannii						c
Glomus sp. OMA7							b
Glomus sp. OMA8							c (rare)
Glomus sp. OMA11				aa, ba	aa
Glomus sp. OMA12			aa			ba

Four replicates were analyzed (a, b, c, d) per plant. Upper case letters show species dominating the spore populations in the corresponding replicate field sample. The spores of Glomus sp.OMA6, Scutellospora sp. OMA10, Glomus sp. OMA2, Glomus sp. OMA3, Glomus sp. OMA9, Glomus sp. OMA5, Glomus sp. OMA7, and Glomus sp. OMA8 (Glomus sp. MEX1; Bashan et al. 2007) resembled spores of Glomus liquidambaris, Scutellospora aurigloba, Glomus versiforme, Glomus invermaium, Glomus glomerulatum, Glomus multicaule, Glomus tortuosum, and Glomus rubiforme, respectively; Glomus sp. OMA3, Glomus sp. OMA5, Glomus sp. OMA8, and Glomus sp. OMA9 are undescribed species; others have not yet been clearly attributable to resembling known species. R. fulgida and R. gregaria are the former Scutellospora fulgida and Scutellospora gregaria, respectively (Oehl et al. 2008)

aIndicates AMF species found in the respective trap cultures

Four replicates were analyzed (a, b, c, d) per plant. Upper case letters show species dominating the spore populations in the corresponding replicate field sample. The spores of pan class="Chemical">pan class="Species">Glomus sppan>.OMA6, papan>n class="Species">Scutellospora sp. OMA10, Glomus sp. OMA2, Glomus sp. OMA3, Glomus sp. OMA9, Glomus sp. OMA5, Glomus sp. OMA7, and Glomus sp. OMA8 (Glomus sp. MEX1; Bashan et al. 2007) resembled spores of Glomus liquidambaris, Scutellospora aurigloba, Glomus versiforme, Glomus invermaium, Glomus glomerulatum, Glomus multicaule, Glomus tortuosum, and Glomus rubiforme, respectively; Glomus sp. OMA3, Glomus sp. OMA5, Glomus sp. OMA8, and Glomus sp. OMA9 are undescribed species; others have not yet been clearly attributable to resembling known species. R. fulgida and R. gregaria are the former Scutellospora fulgida and Scutellospora gregaria, respectively (Oehl et al. 2008)

As shown in Table 2, four AMF species (pan class="Chemical">pan class="Species">Glomus sinuosumpan>, papan>n class="Species">Glomus sp. OMA6, Scutellospora calospora, and Scutellospora sp. OMA10) were exclusively found in both date palm plantations. Other seven AMF species (G. microaggregatum, Glomus sp. OMA9, Racocetra fulgida, R. gregaria, A. spinosa, Glomus etunicatum, and Glomus sp. OMA13) were exclusively found in only one of the date palm plantations. G. aggregatum and G. eburneum were the only two species that were associated with all plants in the four studied habitats including the disturbed area in the experimental station. G. aggregatum, Glomus sp. OMA3, G. eburneum, and Glomus microcarpum were common to the date palm plantations and the undisturbed habitat. In this last location, A. gerdemannii, Glomus sp. OMA7, Glomus sp. OMA8, and Glomus sp. OMA 11 were detected as exclusive taxa. As a consequence of the exclusive presence of different AMF species in different habitats, two main different clusters were formed by the cluster analysis (Fig. 3) suggesting that the AMF communities under the date palms in the two plantations differ from those under the plants in the undisturbed habitat. The AMF community associated with the ruderal plant, P. erioptera, had AMF community composition resembling more to the undisturbed habitat than to that under the date palms. H. kotschyi growing in the sand dunes seems to be closer in its AMF community composition to the vegetation in the undisturbed habitat and to P. erioptera in the disturbed inter-crop area at the experimental station than to the date palms.

Fig. 3

Dendrogram of cluster analysis based on the similarity of AMF morphospecies community composition across the host plants in all studied habitats. Two main clades were formed suggesting that AMF communities associated with date palm in the two plantations have different composition than that of vegetation in the undisturbed habitat (inter-plant area, S. persica, P. cineraria, and Z. hamiense), the disturbed area of the experimental station (P. erioptera), and the sand dunes (H. kotschyi). The group average clustering method and the squared Euclidean distance metric were used

Dendrogram of cluster analysis based on the similarity of AMF morphospecies community composition across the host plants in all studied habitats. Two main clades were formed suggesting that AMF communities associated with pan class="Chemical">pan class="Species">date palmpan> in the two plantations have different composition than that of vegetation in the undisturbed habitat (inter-plant area, papan>n class="Species">S. persica, P. cineraria, and Z. hamiense), the disturbed area of the experimental station (P. erioptera), and the sand dunes (H. kotschyi). The group average clustering method and the squared Euclidean distance metric were used

AMF spore abundance, species richness, and infection potential

The mean spore abundance under the pan class="Chemical">pan class="Species">date palmspan> in both plantations (58.8 and 82.3 spores at the oasis and the experimental station, respectively) was by far higher than that under the ruderal plant, papan>n class="Species">P. erioptera (eight spores) and under natural vegetation in the undisturbed habitat (9.8, 7, 21.5, and 25.8 under Z. hamiense, S. persica, and P. cineraria and at the inter-plant area at the natural site, respectively; Fig. 4a). No AMF spores were found in the field soil coming directly from the root zone of H. kotschyi growing in the sand dunes. The spore abundance was higher under the date palms in the experimental station than under those in the oasis, and it varied among the different host plants in the undisturbed habitat. At this last habitat, the highest mean of spore abundance was found in the inter-plant area. This abundance was significantly higher than that found under Z. hamiense and S. persica but not that under P. cineraria.

Fig. 4

Comparison of the effect of plant species in different study habitats on a AMF spore abundance in the field, b AMF species richness in the field and the trap cultures, and c the infection potential. Values are reported as means (±SD). Non-significant differences between the means are indicated by similar letters above their error bars

Compan class="Chemical">parison of the effect of plant species in different study habitats on a AMF spore abundance in the field, b AMF species richness in the field and the trap cultures, and c the infection potential. Values are reported as means (±SD). Non-significant differences between the means are indicated by similar letters above their error bars

The AMF species richness shown (Table 2) represents the total number of species found in both the field and trap culture samples. Similar to spore abundance, pan class="Chemical">pan class="Species">date palmspan> in both plantations were found to have significantly higher mean species richness (10.3 and 9.3 species at the oasis and the experimental stations, respectively) than the ruderal plant, papan>n class="Species">P. erioptera (4.3 species) and the vegetation on the undisturbed habitat (4.5, 4, 6, and 5 for Z. hamiense, S. persica, P. cineraria, and inter-plant area, respectively) (Fig. 4b). The lowest mean of species richness (one species) was found for H. kotschyi growing in the sand dunes. Unlike in the case of the spore abundance, the date palms in the two plantations did not differ significantly in species richness. In the undisturbed habitat, the richness also did not vary greatly among different host plants. The highest richness at this habitat was found under P. cineraria where it was significantly higher than under both Z. hamiense and S. persica but not than at the undisturbed inter-plant area. Trap culture revealed two species associated with H. kotschyi growing in the sand dunes, G. aggregatum and G. eburneum (Table 2).

The pan class="Chemical">pan class="Disease">AMF infectionpan> potential was higher in the samples from papan>n class="Species">date palms of the oasis (62.6%) than from the ruderal plant, P. erioptera, of the experimental station (35.7%) and the vegetation in the undisturbed habitat (30.9%, 19%, 20.6%, 21.2% for Z. hamiense, S. persica, P. cineraria, and inter-plant area, respectively; Fig. 4c). However, there was no significant difference in the AMF infection potential between the samples from date palms (45.4%) and the ruderal plant, P. erioptera, in the experimental station and from Z. hamiense in the undisturbed habitat. The date palm samples from the two plantations had no statistically different infection potential, and no significant difference was recorded among samples from the different host plants in the undisturbed habitat. Similar to spore abundance and species richness, the lowest infection potential was found with samples from H. kotschyi growing in the sand dunes (2.7%). It was significantly lower than the infection potential of all the other plants in the other habitats except for that of S. persica growing in the undisturbed habitat.

The results of the simple regression analysis revealed a significant positive correlation between the AMF spore abundance in the field and the pan class="Chemical">pan class="Disease">AMF infectionpan> potential (r = 0.676, P = 0.0001) (Fig. 5a). AMF species richness also had a positive correlation with the papan>n class="Disease">AMF infection potential (r = 0.729, P = 0.0001) (Fig. 5b).

Fig. 5

The relationship between a the infection potential and field spore abundance (r = 0.676, P = 0.0001) and b the infection potential and species richness (r = 0.729, P = 0.0001)

Phylotypes detected in roots of date palms

The primer combinations GLOM1310-GLOM5.8R and GLOM1310-ITS4i were used successfully. These primers were designed to target Glomus group A (Schwarzott et al. 2001). In addition, the primer combination GIGA1313-GIGA5.8R, designed to target members of the Gigasporaceae, was successful in one of the samples originating from the pecies">date palms in the oasis. None of the primer combinations targeting the remaining groups of Glomeromycota yielded PCR products. Clones were selected for sequencing on the basis of their restriction fragment length polymorphism (RFLP) patterns. At least one clone representing each RFLP pattern found in each root sample was sequenced. In the samples of 2006, 263 clones were classified into 19 different RFLP patterns, and 33 representatives were sequenced. In the samples of 2004, 43 clones were classified into ten different RFLP patterns, and 16 samples were sequenced. The phylogenetic tree obtained from the sequences from the pan class="Species">date palms at the two plantations in 2006 and from the experimental station only in 2004 is shown in Fig. 6. Based on this phylogeny, nine phylogenetic taxa colonizing the roots of the date palm were discerned. Eight of these nine taxa belong to the Glomus group A. Defining phylotypes in the Glomeromycota often causes problems because there is no molecular species concept and nuclear-encoded ribosomal RNA genes are polymorphic even within single spores (Sanders et al. 1995; Corradi et al. 2007).

Fig. 6

Phylogenetic tree of sequences obtained in this study from the date palms in the oasis plantation (OA) in 2006 and the experimental station (ES) in 2004 and 2006 and sequences from public databases. Sequences obtained in the present study are shown in boldface and are labeled with the database accession number (e.g., GQ406077), the plantation (e.g., oasis OA), the root sample replicate (e.g., 1), RFLP pattern (e.g., A), and the sampling year. The tree was obtained by neighbor-joining analysis of 308 characters from the 18S rDNA. The numbers above the branches are neighbor-joining bootstrap values from 1,000 replications. The tree was rooted using Paraglomus occultum. Note: S. castanea, S. gregaria, and S. fulgida were all recently renamed R. castanea, R. gregaria, and R. fulgida, respectively (Oehl et al. 2008); Kuklospora colombiana is the former Entrophospora colombiana (Sieverding and Oehl 2006); Ambispora appendicula is the former Archaeospora leptoticha (Spain et al. 2006; Walker 2008)

Phylogenetic tree of sequences obtained in this study from the pecies">date palms in the oasis plantation (OA) in 2006 and the experimental station (ES) in 2004 and 2006 and sequences from public databases. Sequences obtained in the present study are shown in boldface and are labeled with the database accession number (e.g., GQ406077), the plantation (e.g., oasis OA), the root sample replicate (e.g., 1), RFLP pattern (e.g., A), and the sampling year. The tree was obtained by neighbor-joining analysis of 308 characters from the 18S rDNA. The numbers above the branches are neighbor-joining bootstrap values from 1,000 replications. The tree was rooted using pan class="Species">Paraglomus occultum. Note: n class="Species">S. castanea, S. gregaria, and S. fulgida were all recently renamed R. castanea, R. gregaria, and R. fulgida, respectively (Oehl et al. 2008); Kuklospora colombiana is the former Entrophospora colombiana (Sieverding and Oehl 2006); Ambispora appendicula is the former Archaeospora leptoticha (Spain et al. 2006; Walker 2008)

Sequence types clustering in the GLOM-A32 group are associated with both a pan class="Chemical">pan class="Species">G. sinuosumpan> and a papan>n class="Species">Glomus clarum sequence at the same time. Some other sequences from semi-arid areas of Namibia also cluster in this group. GLOM-A33, GLOM-A38, GLOM-A39, and the Racocetra group (RACO1) are distantly related to AMF sequences. GLOM-A34, GLOM-A35, GLOM-A36, and GLOM-A37 are not closely related neither to sequences in GenBank nor in our own database. No members from Diversisporaceae, Archaeosporaceae, Acaulosporaceae, or Paraglomeraceae were detected in this study.

The GLOM A-32 phylotype was found exclusively in three out of the four root samples obtained from the pan class="Chemical">pan class="Species">date palmspan> in the experimental station, the GLOM A37 phylotype was found in three of the four root samples from the papan>n class="Species">date palms in the oasis, whereas the GLOM A-39 and RACO 1 groups each were present only in one root sample exclusively in the oasis. GLOM A-33 was found in both plantations but only in one root sample from each plantation.

Discussion

To our knowledge, this is the first study investigating AMF communities in Southern Arabia and uncovering AMF diversity associated with pan class="Chemical">pan class="Species">date palmspan> using morphological and molecular approaches. The detection of 25 morphospecies from a limited number of habitats is remarkable, given that only two morphospecies were known so far from the whole Arabian Peninsula (Khaliel 1989).

Uncovering the AMF community associated with pan class="Chemical">pan class="Species">date palmspan> in plantations and with the vegetation of neighboring undisturbed habitat aimed at elucidating the midterm effect of papan>n class="Species">date palm introduction on the native AMF community. The underlying assumption is that the same AMF community was present in all the three studied habitats before the introduction of the date palms, and therefore, any shift in the AMF community associated with date palms would be due to the land-use change.

pecies">Date palms at both plantations had a different AMF community composition, higher spore abundance, and species richness than the vegetation at the undisturbed habitat. This suggests that there was a high impact of the introduction of pan class="Species">date palms on the native AMF community. In spite of the mechanical disturbance of the uncultivated inter-crop area between the n class="Species">date palm trees in the experimental station, the AMF community under the dominant ruderal plant, P. erioptera, was closer in all the measured parameters to that of the plants in the natural vegetation of the adjacent undisturbed habitat. This is an indication that the propagules of the native AMF in the root zone of this ruderal plant were not affected by the mechanical disturbance, and therefore, the shift in the AMF communities observed in the root zone of date palms is most likely due to the agricultural inputs including irrigation.

Such an unusual trend of increasing AMF abundance and richness upon initiation of agricultural land use might be attributable to the pan class="Chemical">particularity of the studied ecosystems (Ghazanfar and Fisher 1998; Glennie and Singhvi 2002). The natural undisturbed habitats are directly exposed to the harsh climate in Southern Arabia. Drought (Stahl and Christensen 1982; Cui and Nobel 1992) and heat (Koske 1987; Bendavid-Val et al. 1997) are both known for their differential effects on the biology of AMF and consequently on their community composition. Introducing agriculture with irrigation in such a habitat represents a drastic change of the environmental conditions, and thus, a dramatic shift in the AMF community is expected.

It has been reported that agricultural practices (Boddington and Dodd 2000), including soil fertilization (Johnson 1993; Toljander et al. 2008) can affect the species composition of AMF communities. Such practices can lead to a lower AMF species diversity as compared to the corresponding natural sites (Helgason et al. 1998; Oehl et al. 2003). Indeed, in the present study, some of the AMF species occurring in the undisturbed natural habitat were not detected under the pan class="Species">date palms in the adjacent experimental station nor recovered in the respective trap cultures (n class="Species">Glomus sp. OMA5, Glomus macrocarpum, A. gerdemannii, Glomus sp. OMA7, Glomus sp. OMA8, Glomus sp. OMA11, and Glomus sp. OMA12). These species might have been sensitive to the habitat alterations imposed by the plantations. However, the overall species richness was higher under date palms than under the plants growing in the undisturbed habitat. A moderate P fertilization can enhance mycorrhiza formation (Bolan et al. 1984; Koide and Li 1990; Xavier and Germida 1997). Therefore, the sodium bicarbonate-extractable P levels of 71 and 93 ppm in the date palm plantations of the oasis and the experimental station, respectively, might fall within the range that promotes the mycorrhizal association (Schubert and Hayman 1986; Amijee et al. 1989). In an arid and hot ecosystem, similar observations were reported, namely that a low-input agricultural site had higher AMF diversity than an adjacent natural site (Li et al. 2007).

Some AMF species (pan class="Chemical">pan class="Species">G. sinuosumpan>, papan>n class="Species">Glomus sp. OMA6, S. calospora, and Scutellospora sp. OMA10, G. microaggregatum, Glomus sp. OMA9, R. fulgida, R. gregaria, A. spinosa, G. etunicatum, and Glomus sp. OMA13) were found exclusively under date palms, and thus, agricultural inputs might have introduced these species. Soil traces from the tree nurseries is usually attached to the roots of the imported seedlings. Thus, these AMF species might be exotic and might be able to use the new habitat with its low indigenous AMF diversity as an open ecological niche that can be easily colonized, leading to an increase in the AMF species richness.

The lowest AMF spore abundance and species diversity were found in the root zone of H. kotschyi growing in the middle of the sand sea, with only two species recovered in the trap cultures. This apparent rarity might be related to the physical instability of the sandy span class="Chemical">oil at this habitat. Uncovering a rare presence of AMF in this habitat contributes to the efforts (Dutton 1988) to understand and sustain the ecology of this sand sea ecosystem.

The degree of AMF root colonization of pan class="Chemical">pan class="Species">S. bicolorpan> was used as an indirect estimation of the vitality and activity of the AMF community present in the studied habitats. A different response of papan>n class="Species">S. bicolor to the AMF inocula obtained from under the date palms and the vegetation in the undisturbed habitat could indirectly indicate a change on the AMF community structure induced by the land-use change. A higher infection potential was found with the AMF inocula originating from the soil under the date palms than from under plants of the undisturbed habitat (Fig. 4c). This is probably due to the difference of both spore abundance (Fig. 4a) and species richness (Fig. 4b) observed at these habitats. Indeed, a positive correlation was found between AMF spore abundance and species richness and the AMF infection potential (Fig. 5a, b).

Sequences from a pan class="Chemical">pan class="Species">G. sinuosumpan> isolate from Maryland, USA (AJ133706) and from spores obtained from field sites in a semi-dry region of Namibia and identified as papan>n class="Species">G. clarum (AY285854) grouped within the GLOM A-32 clade. Another sequence from a G. clarum isolate in the database (AJ276084) was not closely related to G. sinuosum. Considering that G. sinuosum and G. clarum are easy to distinguish morphologically, this is surprising. If this discrepancy is not due to an error in assigning AY285854 to G. clarum, the most likely explanation would be that a species fitting the morphological description of G. clarum cannot be distinguished from G. sinuosum using 18S DNA sequences. We also found sporocarps of G. sinuosum in the same field but no G. clarum-like spores. Considering this fact and that sequences from cultivated isolates should be more reliable than field-collected spores, all available evidence supports assigning GLOM-A32 to G. sinuosum.

pan class="Chemical">pan class="Species">G. sinuosumpan> was the only named species in our study that could be detected as a morphospecies in the root zone and as a phylospecies in the papan>n class="Species">date palm roots. Although it was the dominant morphospecies in all the samples from date palms of both plantations, it was detected as phylospecies only in the roots of date palms growing in the experimental station. This species seems to be present in various undisturbed ecosystems. G. sinuosum-associated sequences obtained from different regions of the rDNA of AMF were found in unmanaged grassland in Denmark (Rosendahl and Stukenbrock 2004) and in a mountain meadow in Germany (Börstler et al. 2006). Its sporocarps were reported from semi-humid (Tchabi et al. 2008; Tchabi et al. 2009), semi-arid (Muthukumar and Udaiyan 2002), and arid (Shi et al. 2007) undisturbed soils and from low-input agricultural soils (Li et al. 2007). The present study adds date palm with its distinct habitat as a host for G. sinuosum.

The other morphospecies detected in the pan class="Chemical">pan class="Species">date palmpan> rhizosphere might be potential colonizers whose life cycles are controlled by various factors (Azcón-Aguilar et al. 1999) including seasonal fluctuations of their activity (Hijri et al. 2006).

Four taxa were only distantly related to other environmental glomeromycotan phylotypes. GLOM-A33 was related to a sequence obtained from a sandy, alkaline span class="Chemical">pan class="Chemical">oilpan> (Kovács et al. 2007), two characteristics shared with the spapan>n class="Chemical">oil in the present study. The closest relative of GLOM-A38 was an AMF sequence obtained from a Swiss meadow (Sýkorová et al. 2007). GLOM-A39 was another group—in addition to GLOM-A32—that was associated with sequences obtained from roots and root zones of different plants growing in a semi-dry region of Namibia (Uhlmann et al. 2004). The closest relative to RACO 1 group was Racocetra castanea, the former Scutellospora castanea (Oehl et al. 2008).

Four other taxa (GLOM-A34, GLOM-A35, GLOM-A36, GLOM-A37) presently have no closely related sequences in the GenBank and, hence, can neither be related to taxonomically described AMF nor to environmental sequences. We hypothesize that these may be specialists for the ecosystems in Southern Arabia. This hypothesis is in agreement with the observation that there are groups which are specific to certain habitats (Öpik et al. 2006). Recently, an unusual composition of AMF community has been revealed from the distinct ecosystem of the geothermal Yellowstone National pan class="Chemical">Park (USA) (Appoloni et al. 2008).

The exclusive occurrence of GLOM A-32 in the experimental station and GLOM A-37 in the oasis, both in three out of the four root samples, might indicate that these two groups are site-specialist AMF. Some studies showed a community differentiation between different sites, both in temperate (Börstler et al. 2006; Hijri et al. 2006) and in pan class="Chemical">pan class="Disease">dry ecosystemspan> (Wubet et al. 2006). This differentiation might be due to the differences in the agricultural practices that led to different spapan>n class="Chemical">oil chemical properties in the two plantations (Table 1). The two plantations differed also in date palm cultivar diversity, but whether there are cultivar-dependent differences in the association with AMF is not known. The occurrence of what seems to be site-specific taxa in these two plantations might imply that in other different agricultural systems, date palm might be associated with distinct different taxa, but to test this, a larger-scale survey is needed, involving multiple sites of each management type.

In addition to the observed unusual trend of increasing AMF abundance and richness under pan class="Chemical">pan class="Species">date palmspan> supposedly due to agricultural land use, there are two noteworthy points to be mentioned about the composition of the revealed AMF communities. Firstly, neither any of the sequences found nor any morphospecies detected were related to the ubiquitous generalists papan>n class="Species">Glomus intraradices and pan class="Species">Glomus mosseae, the two globally most frequently found AMF (Öpik et al. 2006). We cannot completely rule out that these species might be rare in our study habitats and could potentially be revealed through a more intensive sampling effort and successive trap culturing techniques (Stutz and Morton 1996; Bever et al. 2001) but even such an unusual rarity would be an intriguing finding. Second, we found four undescribed AMF morphospecies and a considerable degree of molecular diversity that either is only distantly related to other known sequences or appears to be new to science. This community composition suggests a striking level of uniqueness most probably resulting from the characteristic environmental settings of Southern Arabia.

Uncovering AMF communities associated with vegetation in Southern Arabia and recognizing their distinctness are important steps toward a better understanding and appreciation of AMF regional and global diversity and ecology.