Abstract

Carissa edulis is a tropical plant belonging to the family Apocynaceae. The species is widely used in the preparation of various herbal medicines. Earlier works in Kenya show that an aqueous extract from the roots of C. edulis has remarkable anti-herpes simplex virus. Due to its medicinal value, the species has been overexploited in its natural range and requires conservation interventions. Studies show that the species has beneficial relationships with arbuscular mycorrhizal fungi (AMF) that can enhance restoration of its population; however, no study has been undertaken to document the diversity of these AMF species. This study evaluated the genetic diversity of AMF associated with the roots of C. edulis within Lake Victoria basin ecosystem of Kenya. A cross-sectional, laboratory-based prospective study was carried out from roots of C. edulis collected from six sites within the ecosystem. Root samples were collected from 6 points (replicates) per site. AMF was assessed through morphological characterization and sequencing of small subunit of ribosomal DNA. Morphological identification identified four genera of AMF (Gigaspora, Acaulospora, Scutellospora, and Glomus) with no significant difference among the sites. Molecular analysis also revealed presence of four genera, but only two (Glomus and Acaulospora) were common for both the analyses with Glomus as the most predominant genera. In all the sites, there were large numbers of spores both in soil and in the roots confirming the association between C. edulis and AMF.

1. Introduction

Carissa edulis is a spiny evergreen shrub or tree having grey bark, with straight woody spines of about 5 cm, often in pairs. Leaves are opposite, leathery, shiny dark green approximately 5 cm with pointed tip and rounded base. The flowers are fragrant in pink-white terminal clusters. The berries are round approximately 1 cm and purple-black when ripe. Carissa edulis is widespread in many parts of Africa [1]. The species grows at forest edges, in forests and woodlands where Euphorbia, Acacia, and Croton commonly occur, especially on rocky hillsides, on clay soils, mainly black cotton soils, in dry and moist low and midlands of 1500–2500 m [1]. The species was once common across Lake Victoria basin ecosystem; however, its distribution is currently dwindling due to overexploitation. East African communities commonly use C. edulis as herbal medicine in different ways. Traditionally, decoctions of C. edulis roots are used as antimalaria and painkiller, and powdered roots are used to treat chest pain, while infusion from roots is used to relieve stomachache [2].

In Kenya, a study by Kenya Medical Research Institute (KEMRI) has shown that an aqueous extract from the roots of C. edulis has remarkable anti-herpes simplex virus (HSV) which is a major opportunistic infection in immunocompromised persons [3]. Herpes simplex virus is a serious threat in HIV/AIDS prone areas in sub-Saharan Africa. As such, C. edulis is highly sought for its medicinal value in East Africa and Africa as a whole. This has led to indiscriminate overexploitation of the species, together with other factors such as deforestation and climate change which has resulted in habitat loss, thereby leaving only patches of C. edulis available across the region. There is therefore the need to develop management strategies for the species conservation including restocking and promotion of on-farm cultivation for sustainable herbal medicine utilization. To achieve this, however, knowledge of both bellow and above ground association with other organisms is imperative. Carissa edulis is believed to have beneficial relationships with some arbuscular mycorrhizal fungi (AMF). However, no studies have been undertaken on the species to document the AMF association within Lake Victoria basin ecosystem to support conservation strategies within the region.

Mycorrhizae are mutualisms formed between fungi and plant roots [4]. The most studied types of mycorrhizae are arbuscular mycorrhizae, ericoid mycorrhizae, arbutoid mycorrhizae, and orchid mycorrhizae. Generally, arbuscular mycorrhizal fungi are very important in land reclamation and functions in sustaining soil fertility and recycling nutrients that increases plant nutrients supply through acquisition of nutrient such as phosphorus and nitrogen that would not be available to plants thus increasing the plant growth vigour [5]. A study in the Coastal region of Eastern Madagascar has linked C. edulis with AMF [6]. Colonization of roots by mycorrhizal fungi also benefits the host plant by improving water uptake. In exchange, the host plant provides the mycorrhizal fungi with carbon from photosynthesis for its energy requirement [7]. It is estimated that 80% of plant species are associated with mycorrhizal fungi, 67% of which is AMF, which may differ in their growth response to a variety of plant species [8]. Thus, the presence or absence of particular AMF species, and overall diversity, can affect the productivity and diversity of host plant communities, both in experimental greenhouse studies and in natural ecosystems [8, 9]. In other instances, species-specific fungal pathogens may promote dominance of one species by limiting recruitment of conspecific species near parent plants thus providing opportunities for only some species to grow [10]. To successfully implement any planting program for the species that has association with AMF, it is important to understand the species growth mechanisms and requirements including below ground associations with beneficial microorganisms. The main objective of the present study was to determine the genetic diversity of AMF found in association with the roots of C. edulis growing within Lake Victoria region of Kenya. To achieve this, both morphological characteristics of AMF and molecular analysis of small subunit of rDNA were employed.

2. Materials and Methods

2.1. Study Area

The study was conducted in Homa Bay (Kendu Bay and Gembe hills), Kisumu (Seme and Nyahera), and Siaya (Got Ramogi and Nyabeda) counties along Lake Victoria basin of Kenya where C. edulis is found (Figure 1). Lake Victoria ecosystem is found in the South Western part of Kenya. The ecosystem has warm, dry, and humid climate with mean annual rainfall ranging within 800–1600 mm. It covers an area of about 4100 km2. The region is rich in biodiversity with different plant species, some of which have medicinal value. Carissa edulis is one of the species found in this region which is most sought for because of its medicinal value. The communities that live in this region rely on C. edulis among other traditional medicinal plants for their local ailments, hence the need for increased productivity of the plant and conservation where its population is threatened.

2.2. Field Sample Collection

We collected fine root samples of C. edulis trees selected randomly where the species is found. From each site, five replicates of single standing C. edulis trees separated by at least 100 m were randomly selected and roots collected. Fine roots from each tree were excavated from at least three lateral roots starting from the trunk and working out towards the fine roots. The roots selected for harvesting were located at different directions from the stem, at least 90° from each other. After harvesting, roots were packed in khaki bags and transported to the laboratory. The roots were rinsed in running tap water to remove attached soil and debris. The samples were cut into pieces of about 1 cm and stored in 70% alcohol until subsequent analysis. Composite soil samples of 3 replicates were also collected under the canopy of the tree from where fine roots were collected. The soil samples were collected from two different depths, within 0–15 cm and 15–30 cm. The soil samples were placed in khaki papers, labelled, and transported to the laboratory for fungal spore characterization.

2.3. Preliminary Spore Isolation and Count

Preliminary spore isolation from the soil samples was done by wet sieving as described in [11] and spore count done using a ×50 dissecting microscope. The number of spores was recorded per sample.

2.4. AM Fungal Trap Culture (Bulking)

After preliminary spore count, soil samples were mixed with autoclaved, medium-sized sand in the ratio of 1 : 4. The samples were packed in one-litre plastic pots in a glasshouse. Sorghum was sown as a trap plant at a density of fifteen seeds per pot. After seedling germination, pots were watered once a day to keep 60% of water-holding capacity (WHC). The pot culture experiment ran for 3 months. Two weeks prior to harvesting, the sorghum was subjected to water stress to enhance sporulation of the existing AMF species. The water stress was ensured by watering the sorghum twice a week instead of everyday.

2.5. Spore Isolation and Morphological Characterization

Pot cultures were harvested from the glass house and AMF spores isolated by wet-sieving technique as described in [11]. Withered sorghum was uprooted, soil thoroughly homogenized, and approximately 50 g from each pot collected for spore isolation. The soil samples were soaked in 200 ml of water and sieved through 45- and 65-micron sieves. Debris from the 65 microns sieve was transferred into falcon tubes and mixed with sucrose solution (48%). The mixture was shaken and left to stand overnight. The debris was sieved through a 65-micron sieve, cleaned, and observed under microscope. Spores were counted and characterized according to their sizes, colour, shape, and presence of subtending hyphae. Spores less than 120 µm were classified as small whereas those greater than 120 µm were regarded as large [12]. Some spores were fished using a spore sucker after successful identification for examination under ×50 fluorescent microscope.

2.6. Root Examination for Spore Percentage Infection

Roots were examined for spore percentage infection using ×50 dissecting microscope. The examined roots were subjected to bleaching in 10% () KOH and stained in 0.05% trypan blue [13]. Root colonization was quantified using the gridline-intersect method [14].

2.7. DNA Isolation and Polymerase Chain Reaction (PCR)

Fine roots of C. edulis were gathered in Eppendorf tubes, air-dried, and then ground into powder using a mixer mill (MM 300, Retsch, Haan, Germany). DNA was extracted using modified cetyl trimethyl ammonium bromide (CTAB) method described in [15]. Ground plant tissues were transferred to a 2.0 ml microfuge tube containing 500 µl isolation buffer (IB) and mixed by vortexing. The mixture was then centrifuged at 10,000 rpm at 4°C for 3 min. Supernatant was removed and 800 µl of fresh IB was added and then vortexed. CTAB (2% CTAB, 20 mM EDTA, 100 mM Tris-HCl, pH 8.0, 1.4 M NaCl, and 0.2%β-mercaptoethanol) solution (800 µl) was added to the mixture and incubated at 65°C for 40 min and then incubated again at 37°C for another 40 min. Chloroform isoamyl alcohol (CIA) (800 µl) was added and mixed by inversion for 10 min. Upper aqueous was transferred to a new 2.0 µl microfuge tube. An equal volume of CIA was added and mixed by inversion for 10 min. The mixture was again centrifuged at 14000 rpm for 10 min at room temperature. Upper aqueous was transferred to a new 1.5 ml microfuge tube and then 1/10 volume of 3 M NaOAc was added. An equal volume of isopropanol was added and mixed by inversion. The mixture was centrifuged at 15,000 rpm at 4°C for 5 min. The supernatants were discarded and then 800 µl 70% ethanol was added and flipped to wash the DNA. It was centrifuged at 15,000 rpm at 4°C for 5 min to pellet the DNA. Supernatant was discarded and DNA pellet was air-dried. The DNA pellet was dissolved using DNase-free water (100 µl) in readiness for PCR analysis.

2.8. Nested Polymerase Chain Reaction Analysis

Partial small subunit rDNA fragments were amplified using nested PCR [16] with the universal eukaryotic primers NS1 and NS4 and a subsequent amplification round with the glomeromycota-specific primers AML1 and AML2 (Table 1). PCR was carried out using 0.5 µl of 10x dNTPs, 1.0 µM of each primer, 5 U of Taq DNA polymerase, 1.0 µ of 10x buffer, 0.5 µl of 25x Mgcl2, 0.5 µl of 25% PVP, and 0.5 µl of 10 mg/ml BSA in total volume of 10.5 μl. The PCR programme was as follows: initial denaturation at 95°C for 3 min, followed by 35 cycles at 95°C for 30 s, 58°C for 1 min, 72°C for 1 min, followed by a final extension period at 72°C for 10 min. A 2 µl aliquot of each PCR product was analyzed in electrophoresis using a 1.5% agarose gel stained in SYBR safe (0.5 µg ml−1) and photographed under UV light. The first PCR products were used as template DNA in a second PCR. The second reaction was performed using AMF specific primers, AML1 and AML2. The following PCR programme was used: initial denaturation at 95°C for 3 min, followed by 35 cycles of 30 s denaturation at 95°C, 1 min primer annealing at 58°C, and 1 min extension at 72°C, followed by a final extension period of 10 min at 72°C. The nested PCR products were then subjected to restriction fragment length polymorphism (RFLP) using restriction enzymes.

2.9. Restriction Fragment Length Polymorphism (RFLP)

The gel electrophoresis method used was as performed in [19, 20] to confirm the amplification of the target region. Four restriction enzymes used to digest nested PCR products previously amplified by specific primers AML1 and AML2 were HaeIII, MpsI, HinfI, and HhaI. Each of the 21 samples that were amplified for 18 s rDNA was divided into 4 Eppendorf tubes and labelled indicating the name of the enzyme. In each tube, 6 µl PCR products, 2.6 µl H2O, 1 µl buffer, 0.1 µl BSA were added, after which each enzyme was added to its respective labelled tube. The reaction tube was incubated at 37°C for four hours. Amplified fragments were separated using 1.5% agarose gel electrophoresis for 70 mins. The gel was stained by adding 3.5 μl SYBR Safe dye (Invitrogen 10,000x concentrate in DMSO) and molten in a microwave oven for 1 min. The digested products were then mixed with 1x gel loading dye (Thermo Scientific) and 2 μl of the products was loaded into the gel. 100 bp Plus DNA ladder (Thermo Scientific) was used for sizing the fragments. The gels were viewed and documented using UV illumination (ATTA E-Graph).

2.10. Sequencing Arbuscular Mycorrhizal Fungi DNA

Nested PCR products were used for BigDye Sanger Sequencing where PCR amplicons were electrophoresed on a 2% agarose gel (Sigma-Aldrich Co., USA) stained in 2 μg/ml ethidium bromide (Sigma-Aldrich Co., USA) solution and visualized using the E-box gel documentation system (Vilber Lourmat, France). The amplicons were subsequently purified using Exonuclease I/Shrimp Alkaline Phosphatase (ExoSap-IT) enzyme (Affymetrix, USA) following the manufactures’ guidelines. The amplicons were sequenced on both strands on an automated 3500xL Genetic Analyzer (Applied Biosystems, USA) using the same primers that were used for amplification. Cycle sequencing was performed using the BigDye Terminator Cycle sequencing kit v3.1 (Applied Biosystems, USA).

2.11. Statistical Analysis

The data were subjected to analysis of variance (ANOVA). The means were compared using t-test whereby the mean of the study was compared to the known existing means [21, 22] and was found to be statistically significant with value <0.05. Sequence editing was done using Sequencher V4.2.2 (Gene Codes Corporation, Ann Arbor, MI, USA). All the sequences were aligned using the multiple sequence comparison alignment tool in MAFFT v6 [23]. Sequences of AMF DNA from C. edulis roots were compared with the existing sequences in the National Centre for Biotechnology Information (NCBI) by BLAST. CLUSTAL MEGA was used to perform multiple sequence alignment (MSA), after which phylogenetic tree was generated using MEGA X [24] by maximum likelihood method and Tamura–Nei model [25].

3. Results

3.1. AMF Spore Count

Overall, preliminary AMF spore count in soil samples from all the sampling sites had significantly () lower number of spore count recording an average (41%) compared to the number of spores after bulking (59%).

Both soil sample spore count before and after bulking showed that Homa Bay sampling area had the highest number of spores with a mean count of 105 before bulking and 148 after bulking per 50 g of soil sample. This was followed by Kisumu at 101 spores before bulking and 144 after bulking while Siaya recorded 100 spores before bulking and 142 after bulking (Figure 2). Varying depths for soil sample collection were considered to find out which layer of the soil had more spores. In all the sites, the depth of 0–15 cm had the highest numbers of spores with the mean count of 52.58 before bulking and 78.94 after bulking compared to the 15–30 cm depth that had mean counts of 49.19 before bulking and 65.42 after bulking (Table 2).

Root percentage infection was also done and overall mean of 55.42% was recorded. Kisumu had the highest mean of 58.67% followed by Homa Bay (54%) while Siaya site recorded the least percentage infection of 50.85.

3.2. Morphological Analysis of AMF

The roots of C. edulis from all the six sampling sites that were examined under fluorescent microscope revealed typical structures of AMF infection with arbuscules and vesicles confirming the association between the plant and AMF (Figure 3). Vesicles tend to take the shape of the plant cell while arbuscules appear like round structures in the roots. Both arbuscules and vesicles picked trypan deep blue colour as observed under a fluorescent microscope.

In the soil samples, the study found that some spores appeared singly or in aggregates in an unorganized hyphal matrix with the layers of the spore wall appearing continuous with the wall of subtending hypha. These spores were borne terminally and ranged within 90–400 µm in diameter with varying colours. Some spores were yellow and others were dark brown while others were light brown fitting the description of AMF of genera Glomus. Six hundred and sixty-nine (669) of the spores that fit this description were found in Homa Bay, and 739 were found in Kisumu while 815 were found in Siaya. Another group of spores appeared on or within the neck (subtending hypha) of a sporiferous saccule with yellowish and light brown colouration and the size ranged within 140–260 µm in diameter. This group fits Acaulospora AMF genera. Two hundred and sixty-five (265) were found in Homa Bay, 312 in Kisumu, and 281 in Siaya. The study also found another group of spores that were bilayered, thin wall auxiliary cells with smooth surfaces produced on hypha. These spores also had small bulb on the subtending hyphae. The size ranged within 100–260 µm with yellowish, light brown, and dark brown colouration. These characteristics were found to match those of Scutellospora. Five hundred and forty-eight (548) of these spores were found in Homa Bay, 579 in Kisumu, and 605 in Siaya. Yet another group of spores appeared to be bilayered and had small bulb on the subtending hyphae, and spores had no ornamentations. These spore sizes ranged within 250–360 µm with a greenish yellow colouration. These features fit those of a Gigaspora. One hundred and three (103) of these spores were found in Homa Bay, 134 in Kisumu, and 151 in Siaya.

The spores were counted per 50 g of soil sample and were found to belong to four genera of AMF: Gigaspora, Acaulospora, Scutellospora, and Glomus. These genera were identified using morphological characteristics of the spores according to INVAM (http://invam.caf.wvu.edu/). Glomus was the most abundant of all the AMF with the mean spore count of 61.31 followed by Scutellospora (47.67), Acaulospora (23.56), and Gigaspora (11.06) (Table 3).

The AMF spores were characterized according to their shapes, sizes, colour, and the presence of subtending hyphae (Figure 4) after bulking. In terms of shape, the majority of spores appeared either spherical or oval. The spherical spores were the highest with the mean spore count of 71.33 while oval ones had mean spore count of 62.58.

In terms of size, there were more large spores recording mean count of 70 than small spores at 62.11. During the characterization, some spores appeared yellow, dark brown, and dark while others light brown with mean counts of 40.89, 37.47, 31.81, and 21.36, respectively. The presence of subtending hyphae was also considered and the majority of the spores counted at the time of examination were lacking the subtending hyphae (mean count of 80.75) compared to spores that had the subtending hyphae at 49.97.

3.3. Molecular Analysis of AMF

Molecular analysis of partial small subunit rDNA from the roots of C. edulis using nested PCR yielded the expected results of 1100 bp by NS1 and NS4 primer pairs and 800 bp by AML1 and AML2 primer pairs according to [16].

Four different restriction enzymes, HaeIII, MpsI, HinfI, and HhaI were used in restriction fragment length polymorphism (RFLP) and revealed molecular diversity (Supplementary Figures 1 and 2). They generated different bands of different sizes against the 100 bp molecular marker. HaeIII generated two prominent sets of bands for AMF. One set generated at 900 bp while another set was generated at 350 bp. HhaI generated two bands in some samples and three bands in others at 700 bp, 400 bp, and 300 bp. In HinfI enzyme, some bands were generated in some samples at 600 bp, 500 bp, and 400 bp and MspI generated fewer bands of all the four restriction enzymes at 800 bp and 400 bp. It was not easy to tell which of the four enzymes was the most effective because each enzyme generated different bands per sample. The bands were assumed to represent different AMF species.

The bands generated by RFLP could not reveal the species or genera of the AMF isolated; therefore, AMF DNA sequencing was performed. Nineteen AMF DNA amplicons that generated clear bands were selected and sequenced. The sequences were phylogenetically analyzed. The obtained phylogenetic tree (Figure 5) generated 14 species. These include seven, three, two, and 1 species of Glomus, Acaulospora, uncultured Glomeromycotina, and uncultured Rhizophagus, respectively. Sequencing of AMF DNA revealed that Glomus was the predominant genera infecting the roots of C. edulis plant within Lake Victoria ecosystem. Even though Kisumu area seems to be having larger number of Glomus, the distribution of the genus is fairly uniform across the region. Acaulospora was evenly distributed across the three sampling sites though slightly high in Homa Bay. The distribution of the uncultured Glomeromycotina clone was also uniformly distributed across the region. Uncultured Rhizophagus was only found in one sampling site in Siaya area. DNA sequencing did not reveal anything on Scutellospora and Gigaspora as did in morphological characterization while it discovered a species of Rhizophagus that was not identified during the morphological analysis.

4. Discussions

Mycorrhizae symbiosis is a key component in helping plants cope with adverse environmental conditions through regulation of their growth and development [26, 27] and had been reported about 400 million years ago [28]. As reported by Smith and Read [29], AMF have mutualistic relationships with more than 80% of terrestrial plant species and the present study shows that AMF successfully develops symbiotic mycorrhizae with C. edulis within the study area.

The main aim for setting up pot trap culture (bulking) using Sorghum bicolor in this study was to propagate AMF spores in good numbers for clearer microscopic analysis and morphological characterization. As described in [30], the procedure is widely used in AMF ecological studies and usually provides sufficient quantities of viable spores with almost all morphological traits that can aid accurate identification and detection of nonsporulating AMF species. A similar strategy has been used in other studies [31] and led to 18% increase in spore count after bulking compared to the preliminary spore count. Besides increased sporulation of AMF in trap cultures, some studies have shown species richness after bulking [32]; however, in the present study, no increase in species richness was observed, similar to what was reported in [33].

Morphological assessment (spore size, colour) of AMF spores from soil samples collected under C. edulis plants suggested the presence of four genera (Gigaspora, Acaulospora, Scutellospora, and Glomus) with variations within the genera. We found that the size of Glomus spores ranged within 90–400 µm with varying colours: yellow, dark brown, and light brown. Acaulospora ranged within 140–260 µm in diameter with yellowish and light brown colouration. Scutellospora had its spore diameter ranging from 100 to 260 µm with yellowish, light brown, and dark brown colouration while Gigaspora recoded spores with diameter ranging from 250 to 360 µm with greenish yellow colouration. Such diversity may reflect the ability to survive or functionality of the species [34]. It was not within the scope of this study to determine exactly what causes variations in the colour and the size of the AMF spores of the same species; however, soil chemistry was suspected to play a key role in the diversity observed. Koch et al. [35] reported similar variations; however, more genetic diversity studies should be conducted within or among AMF species for better understanding of their functional roles in these ecosystems. Another study by [36] found that elevated temperatures reduced the AMF spore diameter by almost 10% indicating that environmental condition could be playing a role in spore size variation. However, [37] found that it was not easy to classify AMF species based on the spore size because of great spore size overlap among AMF species. A similar result was reported in the woodlands of Northern Ethiopia by [38]. The sampling sites in the present study are generally known for rampant livestock grazing, which may have had an impact on the species richness; however, this was not examined. Furthermore, AMF belonging to the Glomus genus are known to be tolerant if not resistant to soil and ecosystem disturbances [39].

The percentage of C. edulis root colonization obtained in this study (54.42%) was comparatively high and similar to those reported by other workers such as [21] for Jatropha curcas in Guantánamo, Cuba and [22] for the same species in Viçosa, Brazil. These findings suggest that terrestrial plants majorly rely on mycorrhizal symbiosis and percentage infection with AMF of about 50% to 60% is enough to thrive the plants in diverse climatic conditions [40]. Among the four genera identified in the present study, Glomus was the most abundant AMF infecting roots of C. edulis within the Lake Victoria ecosystem probably because it is generally sporageous and resilient to disturbance. This finding corroborates the earlier reports suggesting that Glomus is the major AMF genera infecting the roots of most terrestrial plants across the world [41]. Probably, high adaptability of this genus to the variations of temperature and soils characteristics could explain its abundance. In addition, its ability to survive in wider soil pH range and resilience to soil disturbances may be making Glomus dominant among terrestrial plants [42]. In addition, works by [43] suggested that the difference in development process might also explain the reason why Acaulospora and Glomus species dominate roots of most terrestrial plants. This is because they require less time to produce spores than Gigaspora and Scutellospora species. Members of Gigasporaceae typically establish an extensive mycelium in soil and produce fewer spores than those of Acaulosporaceae and Glomaceae. The authors in [38, 44] established that Glomus and Acaulospora are resistant to soil disturbances and changes in ecosystem. This makes Glomus a good alternative for production of inoculum to be used even in the period of formation of plant seedlings including those of C. edulis.

The variation in AMF spore abundance among the sites could not be attributed to any factor and only viewed as a mere chance even though it was suspected that the existence of Rhus natalensis in close proximity with C. edulis had something to do with the large number of spores in some areas. This is an argument that requires further investigation. In general, AMF spore densities are affected by soil erosion and mechanical disturbance due to livestock grazing [38]. In the same way, [44] found that soil disturbance reduces AM fungi spore densities; particularly AM fungi Entrophospora is sensitive to human disturbance and soil erosion.

Our study found that the uppermost soil horizon (depth of 0–15 cm) had more spores than that of deeper soil horizon (15–30 cm). This was attributed to the fact that C. edulis just like many other plants has both tap and lateral roots and the fine roots projecting from lateral roots tend to spread just beneath the earth surface that fall within 0–15 cm where interaction between the roots and AMF spores occurs most. This was not far from the findings reported by Ezeokoli et al. [45] where they collected soil samples randomly at depths of ≤20 cm (topsoil) and >20 cm (subsoil) for AMF spore density assessment and found that topsoil of ≤20 cm had higher number of spores than the subsoil. Our findings also concur with those of [46] showing that the function of AMF in efficient and intensive extraction of soluble nutrients is concentrated in the uppermost soil horizon where organic matter is found and most root growth occurs. This also explains why disturbance of soil especially by grazing livestock would reduce the spore density according to [44].

Restriction fragment length polymorphism (RFLP) analysis revealed the presence of genetic diversity for AMF within the samples collected. However, the restriction by the enzymes could not categorically determine the genus and species of the AMF; therefore DNA sequencing was done to help to assign samples to a genus or a species group. Molecular analysis through sequencing just like morphological analysis revealed that Glomus was the predominant AMF genera sporulating with C. edulis; however, there was inconsistency between molecular and morphological analyses in that only Glomus and Acaulospora were common in both analyses. In molecular analysis, besides Glomus and Acaulospora, uncultured Glomeromycotina, and uncultured Rhizophagus were also identified while in the morphological analysis Scutellospora and Gigaspora were reported. The primers AML1 and AML2 allow robust amplification of AMF sequences directly from field roots, comparable to that obtained with the widely used NS31–AM1 combination [47]. From phylogenetic analysis, all the sequences fell within the phylum Glomeromycota. Lee et al. [47] found that AML1 and AML2 primer pairs had much better specificity and coverage of all known AMF groups. Owing to better specificity of the primer pairs, Gigaspora and Scutellospora as revealed by morphological analysis in this study would have been picked and amplified. Therefore, to miss the two genera by AML1 and AML2 could mean they were sporulating with other plants within the proximity other than the target plant. Čatská et al. [48] also found dissimilar results between molecular and morphological analyses and attributed it to the reduced number of spores of some species found in the field or the dilution of spore in the sample preparation. The gene SSU rDNA that was amplified in this study had been studied to develop PCR-based methodology to determine the presence of AMF in soil and roots using the same set of primers in [49] of Cardiospermum halicacabum L. and Ricinus communis L. This gene evolves relatively slow and may not provide sufficient information to adequately characterize the AMF guild within a community. A more appropriate approach could rely on the sequence variability of the internal transcribed spacer regions (ITS) between the small and large subunits of the rRNA genes [50, 51].

5. Conclusion

In conclusion, we have assessed for the first time within Lake Victoria ecosystem the diversity of mycorrhizal fungi in the soil and roots of C. edulis. A wide variety of arbuscular mycorrhizal spores were identified in the rhizosphere of this plant by morphological investigations, while the microscopic observations of intraradicular hyphae indicated a high colonization rate of C. edulis. The results obtained in this study constitute a solid scientific basis and another pathway to preserve and maybe to enhance restoration of C. edulis populations within the region.

Data Availability

Data are available on request from the authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank Kenya Forestry Research Institute (KEFRI) for technical support and the entire molecular laboratory staff at KEFRI for the facility. They also appreciate Prof. Wallace Bulimo of University of Nairobi and the head of Walter Reed Project (WRP) Influenza Reference Laboratory, Nairobi, together with his staff for offering the DNA sequencing facilities.

Supplementary Materials

Supplementary data consists of electropherograms of PCR products of the target gene area and the patterns of restricted fragments of the PCR products after digestion with restriction enzymes. Figure S1: agarose gel showing PCR products from 1 to 21 while M in both ends represent 100 bp marker. Figure S2: example of a HaeIII restriction fragment length polymorphism (RFLP) electropherogram. (Supplementary Materials)