# Phylogenetic Analysis of Nucleotide Sequences from the ITS Region and Biological Characterization of Nematophagous Fungi from Morelos, Mexico

Authors: Gloria Sarahí Castañeda-Ramírez; Pedro Mendoza-de-Gives; Liliana Aguilar-Marcelino; María Eugenia López-Arellano; Jesús Hernández-Romano
Journal: Journal of Mycology
 (2016)

---

## Abstract

We determined the morphological taxonomy of eighteen nematophagous fungi (NF), as well as theirin vitro predatory activity againstHaemonchus contortus infective larvae (L3). Fungi were classified into six genera and three species, the most common of which wereMonacrosporium eudermatum andArthrobotrys oligospora. We then sequenced five NF isolates using ITS4 and ITS5 primers. These sequences showed high identity with sequences from the NCBI database (98-99%). In contrast, alignments among the same genera and species demonstrated 83–97% identity. Polymorphisms observed betweenArthrobotrys andMonacrosporium appear to be associated with differences in biological function, nonspecific mutations, evolutionary processes, feeding behaviour, predatory activity, and microecosystems.


---

## Body

## 1. Introduction

The livestock industry is commonly affected by gastrointestinal parasitic nematodiasis, which causes serious deterioration in the health of both cattle and sheep worldwide [1]. A number of methods of control have been proposed to reduce both the parasitic burden in the animals and the economic consequences for animal production; however, the only control method that is currently available in most countries is the regular administration of chemical anthelmintic drugs as deworming treatments. While such drugs help to reduce the parasitic burden in the animals, disadvantages can include the rising threat of inefficacy due to anthelmintic resistance and residues in animal products that are possibly toxic for human consumption and for the environment. For several years, therefore, alternative methods of control against parasites have been under investigation.Nematophagous fungi are microorganisms in the soil that act as natural antagonists of nematodes and are considered potential biological control agents against agricultural plagues and ruminant parasitic nematodes. Using selected strains of nematophagous fungi for biotechnological purposes requires performing several studies to assess trapping efficiency against target nematodes. Precise taxonomic identification of nematophagous fungi isolates is the first step in selecting possible candidates to assess activity against nematodes in both laboratory and field trials. Traditional taxonomic identification of nematophagous fungi is based on microscopic observation and measurement of different fungal structures, that is, mycelia, spores, conidia, conidiophores, trapping devices, and other typical features [2]. Valuable taxonomic identification guides for nematophagous fungi were published by Drechsler, 1941 [3], Dollfus, 1946 [4], and Duddington, 1955 [5], and a highly descriptive and useful taxonomic key was published by Cooke and Godfrey [2]. More recently, taxonomic revisions concerningDactylaria and theDactylaria-Monacrosporium complex were published by de Hoog [6] and Rubner [7], respectively. These works are considered useful tools for taxonomic identification of nematophagous fungi. Nevertheless, this method of identification can be difficult because structures among fungal species are similar and can easily lead to confusion. For more precise taxonomic identification of nematophagous fungi, comparisons of 28S rDNA, 5.8S rDNA, and β-tubulin gene sequences can be performed [8]. These molecular techniques can be used to identify target strains rapidly and accurately. Additionally, phylogenetic relationships among strains can be compared. Genotyping of nematophagous fungi is considered a valuable and accurate molecular tool for selection of fungi suitable for biotechnological purposes for the livestock and agricultural industries worldwide.The state of Morelos in Mexico is rich in microflora and microfauna, though there is currently only limited information on nematophagous fungi at this location. This study, therefore, aimed to analyze the nucleotide sequences of the ITS regions of nematophagous fungi from Morelos, Mexico, for phylogenetic analysis and also to evaluate the biological activity of isolates against infective larvae of the sheep nematophagous nematodeHaemonchus contortus.


## 2. Materials and Methods

### 2.1. Area of Study for Fungal Isolation

Soil sampling for fungal isolation was carried out at Northern Cuernavaca City in the state of Morelos, Southern Mexico City (Figure1). Morelos possesses significant environmental diversity characterized by a rich variety of vegetation, including conifer and Quercus woodlands. Cuernavaca forms part of the Anahuac valleys and mountains, particularly the Chichinautzin mountain range slope. The northern part of Cuernavaca possesses a subhumid temperate climate with summer rains and high humidity and contains mesophilic woodlands of pine and oak trees. The yearly average temperature is 21.1°C, with yearly average precipitation between 800 and 1,500 mm. The highest temperatures (24–28°C) occur in April and May. The lowest monthly temperature average (below 15°C) occurs in December and January [9]. The soil is humic andosol and ocrico andosol derived from volcano ash that is adequate for forestry use (INAFED) [10, 11]. Such climate conditions promote the development of a wide diversity of floral and faunal species [11]. The high moisture and temperature as well as the soil characteristics that predominate in this woodland area are expected to trigger the proliferation of an important diversity of microorganisms with potential biotechnological uses. Nevertheless, there is no current information on the presence of microorganisms with potential biotechnological uses in this geographic region of Mexico.Figure 1
Map of the Mexican Republic: (a) the State of Morelos, (b) Cuernavaca City.Most samples were collected from the base of different tree species from the Santa María Ahuacatitlán Area. The geographic area was divided into two subareas: Zone 1 was characterized by a mixed woodland with oak trees, holm oak trees, and several types of pine trees, while Zone 2 was a woodland area consisting mainly of pine trees (Figure2; Table 1).Table 1
Sites and source of material collection for fungal isolation from Morelos, Mexico.
Zone 1
Zone 2 (∗)
Identification number
Sampling sites
Identification number
Sampling sites
2.1
Soil from the base of a tree (St Elena street)
3.1
Soil from a particular house
2.2
Soil from a particular house (St Elena street)
3.2
Soil from the base of a pine tree
2.3
Soil from a sidewalk beside the street (St Elena street)
3.3
Compost
2.4
Soil + dry leaves (St Elena street)
3.4
Soil from the base of flowers
2.5
Soil from the Chichinautzin ecological corridor entrance
3.5
Walking Corridor
2.6
Soil from Chichinautzin ecological corridor sidewalk
3.6
Soils collected beside the rail train
2.7
Soil from St. Maria club garden (St. Maria)
3.7
Soils collected beside the rail train
2.8
Soil from a garden at the Unión de las Américas primary school
3.8
Soil from orange colour dots in soil (microbial consortium)
2.9
Soil collected beside the Mexico-Cuernavaca highway
3.9
Soil from gray colour dots in soil (microbial consortium)
2.10
Soil from a particular garden beside the Mexico-Cuernavaca highway
2.11
Soil from the base of a monument (St. Maria)
(
∗) All sampling sites in Zone 2 are situated in Ahuatepec Alarcón community.Figure 2
Map of Cuernavaca City. Triangles showing the area of sampling sites. Photos showing a view of woodlands predominant at Santa María Ahuacactitlán area, Cuernavaca, Morelos State, Mexico


### 2.2. Fungal Isolation

Samples (5 g) of soil, dry leaves, or compost were collected from every sampling site and transported to the lab. A previously described method of sprinkling soil on 2% agar plates for fungal isolation [12] was used. Each sample was sprinkled on the surface of 9 cm diameter water agar plates and incubated at room temperature (18–25°C) and 60% humidity. After 2 days, an undetermined number of mixed stages of free-living nematodes (Panagrellus redivivus) were added to the agar surface of each plate to promote nematophagous fungal development. Every day after the third day, the incubation plate surface was viewed under a stereomicroscope to search for the presence of either nematophagous fungal aerial structures or trapped nematodes. Conidia or conidiophores were transferred with a fine, sterile needle to sterile water agar plates to isolate the observed fungi. This process was repeated as often as necessary until fungal isolates remained as pure culture. Once pure fungal cultures were obtained on water agar plates, fungi were transferred to potato dextrose agar plates and incubated for 20 days at room temperature (18–25°C). After this period, abundant mycelia were produced on the agar surface. Mycelia were recovered from plates by scrapping the agar surface with a sterile spatula. The fungal material was eventually transferred into 1.5 mL microcentrifuge tubes and maintained at −20°C until use.


### 2.3. Fungal Identification through Traditional Taxonomy (Morphometric System)

Both the shapes and sizes of different fungal structures (conidia, conidiophores, trapping devices (traps), hyphae, and mycelia) were observed under a microscope. These structures were photographed (at 40x, 75x, and 100x magnification) using a Leica Comax® CRC-3141 colour camera CCD integrated with a microscope. Taxonomic identification of fungal genera and species was achieved following the identification keys published by Cooke and Godfrey [2], de Hoog [6], and Rubner [7].


### 2.4. Molecular Identification

DNA extraction was achieved following the methodology established for entomopathogenic fungi DNA extraction described by Tigano-Milani et al., 1995 [13]. The DNA amplification was achieved using PCR technique with the following primers: ITS5 (GGA AGT AAA AGT CGT AAC AAG G) and ITS4 (TCC TCC GCT TAT TGA TAT GC). These primers were used to amplify regions ITS1, 5.8S, and ITS2, as well as 18S and 28S regions [14]. Amplicons of each fungal strain were obtained using the following mix reaction: DNA 5 μL (100 ng/μL), primers 2 μL each (100 ng/μL), MgCl2 2 μL (2 mM), dNTP mix 1 μL (10 mM), buffer 5x 5 μL, and DNA Taq Polymerase 1 μL (1 U) to obtain 25 μL final volume. PCR conditions were an initial denaturation stage of 94°C for 10 min, a second stage of 35 cycles of denaturation at 94°C for 60 sec, annealing at 60°C for 50 seconds and extension at 72°C for 90 sec, and a final extension stage at 72°C for 10 min. PCR products were visualized by agarose gel electrophoresis. TAE 1x was used as a running buffer. The PCR products were purified using the Gel and PCR Clean-Up system kit (PROMEGA). The purified PCR products were quantified and individual sequencing reactions were prepared, using 100 ng to 120 ng by reaction. Sequencing was performed at the Institute of Biotechnology, UNAM, using an Applied Biosystems Sequencer (7700). Sequence homology analysis was done using the Basic Local Alignment Search Tool (BLAST) [15]. Identity percentages and mismatches count were obtained to confirm the traditional identification based on the morphology of fungal strains. A multiple aligning analysis was performed and analyzed separately to determine the similarity percentages between the selected fungal strains using ITS2 sequence. A total of 36 fungal sequences with different similarity levels were recovered from GenBank. ITS1, ITS2, and 5.8S segments were identified and ITS2 aligning was performed using the CLC Sequence Viewer software [16]. The gap cost and the extension cost were assigned as 10 and 1, respectively. A phenogram was generated using the UPGMA algorithm with the Jukes-Cantor nucleotide distance measure and bootstrap analysis was done with 100 replicates. Additionally, a comparative matrix of transversions, transitions and inversions, was performed in order to compare the selected strains.


### 2.5. Nematodes

#### 2.5.1.Panagrellus redivivus Free-Living Nematodes

AP. redivivus strain belonging to the CENID-Parasitología Veterinaria, INIFAP was used. Nematodes were produced in mass using a medium containing oat flakes and water [17]. Finally, the nematode suspension was transferred using the Baermann funnel technique to obtain a clean nematode suspension [12].


#### 2.5.2.Haemonchus contortus (L3) Infective Larvae

A Mexican strain of the sheep parasitic nematodeHaemonchus contortus from the INIFAP germplasm collection was used. A Pelibuey hair sheep was orally infected with the parasite to act as an “egg donor.” After a 28-day prepatent period, nematode eggs were detected in the sheep’s faeces through the McMaster technique [18]. Sheep faeces were directly collected from the rectum of the animal. Faecal material was processed to elaborate faecal cultures [19]. These cultures were incubated for 7 days at room temperature (17–28°C). Each day, they were mixed with a wooden spoon and hydrated with tap water to promote maximum egg hatching and larva production. Larvae were extracted from faecal cultures using the Baermann funnel technique (for 24 h) to eventually obtain a clean nematode suspension [12]. Residues and detritus were separated and removed from the larval suspension using the sucrose density gradient technique [20]. The larval suspension was washed three times using sterile water to remove sucrose residues. Finally, larvae were resuspended in sterile distilled water and kept at 4°C until use.


### 2.6. Fungal Predatory Activity Assessment

Every fungal strain was cultured in 5 cm Petri dishes containing 2% water agar (n = 10). After 7 days’ incubation at room temperature (25–28°C), approximately five hundredH. contortus infective larvae (L3) were added to the surface of each fungal plate. Additionally, a set of water agar plates (n = 10) without any fungi was used as a control and treated as described above. Later, the agar plate surface was viewed under a microscope (10x and 40x) to visualize the predatory activity of fungi. The total number of larvae contained in each plate of each group was recovered through the Baermann funnel technique to assess fungal predatory activity.The average number of recovered larvae was obtained and compared for both groups to estimate the percentage reduction in larvae due to the predatory action of the nematophagous fungi using the following formula [21]:(1)
A
=
X
X
+
Y
∗
100
,where A is reduction percentage, X is average of recovered larvae from plates with fungi, and Y is average of recovered larvae from plates without fungi.


### 2.7. Data Analysis

Data were analyzed using proofs of normality. A Kruskal-Wallis analysis with a completely random design was used [22]. The average number of recovered larvae from each of the two groups was considered as the variable/response. A multiple average comparison analysis was also performed using SAS [23]. An α value = 0.01 was used.


### 2.8. Database Submission

The nucleotide sequences derived from the nematophagous fungi reported in this work were submitted to GenBank (accession numbers KC663623, KC663625, KC663626, KC663624, and KC663627).


## 2.1. Area of Study for Fungal Isolation

Soil sampling for fungal isolation was carried out at Northern Cuernavaca City in the state of Morelos, Southern Mexico City (Figure1). Morelos possesses significant environmental diversity characterized by a rich variety of vegetation, including conifer and Quercus woodlands. Cuernavaca forms part of the Anahuac valleys and mountains, particularly the Chichinautzin mountain range slope. The northern part of Cuernavaca possesses a subhumid temperate climate with summer rains and high humidity and contains mesophilic woodlands of pine and oak trees. The yearly average temperature is 21.1°C, with yearly average precipitation between 800 and 1,500 mm. The highest temperatures (24–28°C) occur in April and May. The lowest monthly temperature average (below 15°C) occurs in December and January [9]. The soil is humic andosol and ocrico andosol derived from volcano ash that is adequate for forestry use (INAFED) [10, 11]. Such climate conditions promote the development of a wide diversity of floral and faunal species [11]. The high moisture and temperature as well as the soil characteristics that predominate in this woodland area are expected to trigger the proliferation of an important diversity of microorganisms with potential biotechnological uses. Nevertheless, there is no current information on the presence of microorganisms with potential biotechnological uses in this geographic region of Mexico.Figure 1
Map of the Mexican Republic: (a) the State of Morelos, (b) Cuernavaca City.Most samples were collected from the base of different tree species from the Santa María Ahuacatitlán Area. The geographic area was divided into two subareas: Zone 1 was characterized by a mixed woodland with oak trees, holm oak trees, and several types of pine trees, while Zone 2 was a woodland area consisting mainly of pine trees (Figure2; Table 1).Table 1
Sites and source of material collection for fungal isolation from Morelos, Mexico.
Zone 1
Zone 2 (∗)
Identification number
Sampling sites
Identification number
Sampling sites
2.1
Soil from the base of a tree (St Elena street)
3.1
Soil from a particular house
2.2
Soil from a particular house (St Elena street)
3.2
Soil from the base of a pine tree
2.3
Soil from a sidewalk beside the street (St Elena street)
3.3
Compost
2.4
Soil + dry leaves (St Elena street)
3.4
Soil from the base of flowers
2.5
Soil from the Chichinautzin ecological corridor entrance
3.5
Walking Corridor
2.6
Soil from Chichinautzin ecological corridor sidewalk
3.6
Soils collected beside the rail train
2.7
Soil from St. Maria club garden (St. Maria)
3.7
Soils collected beside the rail train
2.8
Soil from a garden at the Unión de las Américas primary school
3.8
Soil from orange colour dots in soil (microbial consortium)
2.9
Soil collected beside the Mexico-Cuernavaca highway
3.9
Soil from gray colour dots in soil (microbial consortium)
2.10
Soil from a particular garden beside the Mexico-Cuernavaca highway
2.11
Soil from the base of a monument (St. Maria)
(
∗) All sampling sites in Zone 2 are situated in Ahuatepec Alarcón community.Figure 2
Map of Cuernavaca City. Triangles showing the area of sampling sites. Photos showing a view of woodlands predominant at Santa María Ahuacactitlán area, Cuernavaca, Morelos State, Mexico


## 2.2. Fungal Isolation

Samples (5 g) of soil, dry leaves, or compost were collected from every sampling site and transported to the lab. A previously described method of sprinkling soil on 2% agar plates for fungal isolation [12] was used. Each sample was sprinkled on the surface of 9 cm diameter water agar plates and incubated at room temperature (18–25°C) and 60% humidity. After 2 days, an undetermined number of mixed stages of free-living nematodes (Panagrellus redivivus) were added to the agar surface of each plate to promote nematophagous fungal development. Every day after the third day, the incubation plate surface was viewed under a stereomicroscope to search for the presence of either nematophagous fungal aerial structures or trapped nematodes. Conidia or conidiophores were transferred with a fine, sterile needle to sterile water agar plates to isolate the observed fungi. This process was repeated as often as necessary until fungal isolates remained as pure culture. Once pure fungal cultures were obtained on water agar plates, fungi were transferred to potato dextrose agar plates and incubated for 20 days at room temperature (18–25°C). After this period, abundant mycelia were produced on the agar surface. Mycelia were recovered from plates by scrapping the agar surface with a sterile spatula. The fungal material was eventually transferred into 1.5 mL microcentrifuge tubes and maintained at −20°C until use.


## 2.3. Fungal Identification through Traditional Taxonomy (Morphometric System)

Both the shapes and sizes of different fungal structures (conidia, conidiophores, trapping devices (traps), hyphae, and mycelia) were observed under a microscope. These structures were photographed (at 40x, 75x, and 100x magnification) using a Leica Comax® CRC-3141 colour camera CCD integrated with a microscope. Taxonomic identification of fungal genera and species was achieved following the identification keys published by Cooke and Godfrey [2], de Hoog [6], and Rubner [7].


## 2.4. Molecular Identification

DNA extraction was achieved following the methodology established for entomopathogenic fungi DNA extraction described by Tigano-Milani et al., 1995 [13]. The DNA amplification was achieved using PCR technique with the following primers: ITS5 (GGA AGT AAA AGT CGT AAC AAG G) and ITS4 (TCC TCC GCT TAT TGA TAT GC). These primers were used to amplify regions ITS1, 5.8S, and ITS2, as well as 18S and 28S regions [14]. Amplicons of each fungal strain were obtained using the following mix reaction: DNA 5 μL (100 ng/μL), primers 2 μL each (100 ng/μL), MgCl2 2 μL (2 mM), dNTP mix 1 μL (10 mM), buffer 5x 5 μL, and DNA Taq Polymerase 1 μL (1 U) to obtain 25 μL final volume. PCR conditions were an initial denaturation stage of 94°C for 10 min, a second stage of 35 cycles of denaturation at 94°C for 60 sec, annealing at 60°C for 50 seconds and extension at 72°C for 90 sec, and a final extension stage at 72°C for 10 min. PCR products were visualized by agarose gel electrophoresis. TAE 1x was used as a running buffer. The PCR products were purified using the Gel and PCR Clean-Up system kit (PROMEGA). The purified PCR products were quantified and individual sequencing reactions were prepared, using 100 ng to 120 ng by reaction. Sequencing was performed at the Institute of Biotechnology, UNAM, using an Applied Biosystems Sequencer (7700). Sequence homology analysis was done using the Basic Local Alignment Search Tool (BLAST) [15]. Identity percentages and mismatches count were obtained to confirm the traditional identification based on the morphology of fungal strains. A multiple aligning analysis was performed and analyzed separately to determine the similarity percentages between the selected fungal strains using ITS2 sequence. A total of 36 fungal sequences with different similarity levels were recovered from GenBank. ITS1, ITS2, and 5.8S segments were identified and ITS2 aligning was performed using the CLC Sequence Viewer software [16]. The gap cost and the extension cost were assigned as 10 and 1, respectively. A phenogram was generated using the UPGMA algorithm with the Jukes-Cantor nucleotide distance measure and bootstrap analysis was done with 100 replicates. Additionally, a comparative matrix of transversions, transitions and inversions, was performed in order to compare the selected strains.


## 2.5. Nematodes

### 2.5.1.Panagrellus redivivus Free-Living Nematodes

AP. redivivus strain belonging to the CENID-Parasitología Veterinaria, INIFAP was used. Nematodes were produced in mass using a medium containing oat flakes and water [17]. Finally, the nematode suspension was transferred using the Baermann funnel technique to obtain a clean nematode suspension [12].


### 2.5.2.Haemonchus contortus (L3) Infective Larvae

A Mexican strain of the sheep parasitic nematodeHaemonchus contortus from the INIFAP germplasm collection was used. A Pelibuey hair sheep was orally infected with the parasite to act as an “egg donor.” After a 28-day prepatent period, nematode eggs were detected in the sheep’s faeces through the McMaster technique [18]. Sheep faeces were directly collected from the rectum of the animal. Faecal material was processed to elaborate faecal cultures [19]. These cultures were incubated for 7 days at room temperature (17–28°C). Each day, they were mixed with a wooden spoon and hydrated with tap water to promote maximum egg hatching and larva production. Larvae were extracted from faecal cultures using the Baermann funnel technique (for 24 h) to eventually obtain a clean nematode suspension [12]. Residues and detritus were separated and removed from the larval suspension using the sucrose density gradient technique [20]. The larval suspension was washed three times using sterile water to remove sucrose residues. Finally, larvae were resuspended in sterile distilled water and kept at 4°C until use.


## 2.5.1.Panagrellus redivivus Free-Living Nematodes

AP. redivivus strain belonging to the CENID-Parasitología Veterinaria, INIFAP was used. Nematodes were produced in mass using a medium containing oat flakes and water [17]. Finally, the nematode suspension was transferred using the Baermann funnel technique to obtain a clean nematode suspension [12].


## 2.5.2.Haemonchus contortus (L3) Infective Larvae

A Mexican strain of the sheep parasitic nematodeHaemonchus contortus from the INIFAP germplasm collection was used. A Pelibuey hair sheep was orally infected with the parasite to act as an “egg donor.” After a 28-day prepatent period, nematode eggs were detected in the sheep’s faeces through the McMaster technique [18]. Sheep faeces were directly collected from the rectum of the animal. Faecal material was processed to elaborate faecal cultures [19]. These cultures were incubated for 7 days at room temperature (17–28°C). Each day, they were mixed with a wooden spoon and hydrated with tap water to promote maximum egg hatching and larva production. Larvae were extracted from faecal cultures using the Baermann funnel technique (for 24 h) to eventually obtain a clean nematode suspension [12]. Residues and detritus were separated and removed from the larval suspension using the sucrose density gradient technique [20]. The larval suspension was washed three times using sterile water to remove sucrose residues. Finally, larvae were resuspended in sterile distilled water and kept at 4°C until use.


## 2.6. Fungal Predatory Activity Assessment

Every fungal strain was cultured in 5 cm Petri dishes containing 2% water agar (n = 10). After 7 days’ incubation at room temperature (25–28°C), approximately five hundredH. contortus infective larvae (L3) were added to the surface of each fungal plate. Additionally, a set of water agar plates (n = 10) without any fungi was used as a control and treated as described above. Later, the agar plate surface was viewed under a microscope (10x and 40x) to visualize the predatory activity of fungi. The total number of larvae contained in each plate of each group was recovered through the Baermann funnel technique to assess fungal predatory activity.The average number of recovered larvae was obtained and compared for both groups to estimate the percentage reduction in larvae due to the predatory action of the nematophagous fungi using the following formula [21]:(1)
A
=
X
X
+
Y
∗
100
,where A is reduction percentage, X is average of recovered larvae from plates with fungi, and Y is average of recovered larvae from plates without fungi.


## 2.7. Data Analysis

Data were analyzed using proofs of normality. A Kruskal-Wallis analysis with a completely random design was used [22]. The average number of recovered larvae from each of the two groups was considered as the variable/response. A multiple average comparison analysis was also performed using SAS [23]. An α value = 0.01 was used.


## 2.8. Database Submission

The nucleotide sequences derived from the nematophagous fungi reported in this work were submitted to GenBank (accession numbers KC663623, KC663625, KC663626, KC663624, and KC663627).


## 3. Results

### 3.1. Traditional Taxonomy

Morphological and morphometric taxonomic features of different fungal structures are shown in Table2. Eighteen fungal isolates were obtained in total from the different geographic areas and from different substrates. Fungi were identified as belonging to three genera and six species of nematophagous fungi. Most of the fungal isolates corresponded to the speciesM. eudermatum (8 isolates); four isolates were identified asA. oligospora, two asA. musiformis, and one asA. kirghizica. Two isolates were determined asTrichothesium globosporum var.globosporum and one isolate asT. g. var.roseum (Table 2). Different morphological features of the nematophagous fungi are shown in Table 3.Table 2
Morphological and morphometric features of the different genera and species of nematophagous fungi isolated from soil samples from Northern Cuernavaca City, Morelos, Mexico.
Genus/species strain
Conidium size (μm)
Conidiophore length (μm)
Trapping devices
Specific features
Length
Width
M. eudermatum (Z-2.1)
27.5–40
15–20
300
Three-dimensional adhesive nets
One single apical conidium, septum situated at the middle conidium, segmented and branched conidiophore
M. eudermatum (Z-2.4)
27.5–40
15–20
450
M. eudermatum (Z-3.1)
25–42.5
12–25
400
M. eudermatum (Z-2.7)
32.5–45
15–22
350
M. eudermatum (Z-3.3)
27–42.5
12–22
370
M. eudermatum (Z-3.5)
27.5–40
15–22
200
M. eudermatum (Z-3.6)
27.5–40
12–20
300
M. eudermatum (Z-3.8)
26–40
13–21
320
A. oligospora (Z-2.2)
17.5–30
12–20
550
Three-dimensional adhesive nets
Conidiophores present conidia clusters; septum situated below the middle conidium
A. oligospora (Z-2.6)
15–25
10–12
300
A. oligospora (Z-2.8)
17.5–25
10–12
180
A. oligospora (Z-2.10)
20–25
7.5–12
340
A. musiformis (Z-2.5)
25–40
7.5–15
350
Three-dimensional adhesive nets
Elongated and segmented conidia; conidia are produced from small denticles
A. musiformis (Z-3.2)
22.5–35
7.5–12
340
A. kirghizica (Z-2.9)
15–22.5
7.5–12
340
Three-dimensional adhesive nets
Ovoid shape conidia; septum is situated slightly below the middle conidium
T. globosporumvar.globosporum (Z-2.3)
15–30
7.5–15
310
Three-dimensional adhesive nets
Conidium distal cell is larger; conidia are frequently produced in clusters
T. g.var.globosporum (Z-2.11)
22.5–30
10–17
350
T. g.var.roseum(Z-3.4)
10–20
10–20
370
Three-dimensional adhesive nets
Apiculated base, distal cell is longerTable 3
Microphotographs showing the aspect of different nematophagous fungi structures (conidia, conidiophore, and trapping devices) isolated from Morelos, Mexico. Scale bars = 20µm.
Genus/species
Identification number
Conidia
Trapping devices
Conidiophore
Monacrosporium eudermatum
Z-2.1
Z-2.4
Z-2.7
Z-3.1
M. eudermatum
Z-3.3
Z-3.5
Z-3.6
Z-3.8
A. oligospora
Z-2.2
Z-2.6
Z-2.8
Z-2.10
A. musiformis
Z-2.5
Z-3.2
A. kirghizica
Z-2.9
Trichothecium globosporumvar. globosporum
Z-2.3
T. g.var. globosporum
Z-2.11
T. g.var. roseum
Z-3.4


### 3.2.In Vitro Predatory Activity of Nematophagous Fungi

Results of the assessment of thein vitro predatory activity of the nematophagous fungi againstH. contortus infective larvae are shown in Table 4. The larval reduction percentages recorded for the different isolates ranged between 32.4 and 96.7%. The lowest predation percentages were observed forM. eudermatum (isolate Z-2.4) at 32.4% andA. kirghizica (Z-2.9) at 36.4%. SpeciesT. globosporum (Z-2.3, Z-2.11),A. oligospora (Z-2.6),A. musiformis (Z-2.5), andM. eudermatum (Z-2.1) showed moderate predation activities (59.1, 46.9, 56.2, 62.2, and 67.6%, resp.). The highest predatory activities were observed withT. globosporum var.roseum (Z-3.4) at 71.2%,M. eudermatum (Z-3.6) at 74.3%,M. eudermatum (Z-3.1) at 78.4%,M. eudermatum (Z-3.8) at 85.4%,M. eudermatum (Z-3.3) at 91.53%,A. oligospora (Z-2.8) at 93.07%,A. oligospora (Z-2.10) at 94.6%,M. eudermatum (Z-2.7) at 95.9%, andA. oligospora (Z-2.2) at 96.7%.Table 4
Mean numberHaemonchus contortus infective larvae and reduction percentages recovered after 7 days’ in vitrointeraction with different nematophagous fungi isolates and reduction percentages.
Fungal isolates
Control groupX
- (±SD)
Treated groupX
- (±SD)
Larvae reduction (%)
Letters
Monacrosporium eudermatum (Z-2.1)
296 ± 92.7
96 ± 50.6
67.6
a, b, c, d, f, g
M. eudermatum (Z-2.4)
296 ± 92.7
200.1 ± 13.14
32.4
b, c, d, e, f, l
M. eudermatum (Z-3.1)
296 ± 92.7
64 ± 42.99
78.4
a, b, c, d, g, h, i, k
M. eudermatum (Z-2.7∗)
296 ± 92.7
12 ± 19.3
95.9
g, i,j, k
M. eudermatum (Z-3.3∗)
520 ± 253.7
44 ± 78.7
91.5
a, c, g, h, i,j, k
M. eudermatum (Z-3.6)
296 ± 92.7
76 ± 60.1
74.3
a, b, c, d, e, g, i, k
M. eudermatum (Z-3.8)
520 ± 253.7
76 ± 74.1
85.4
a, b, c, d, e, g, i, k
Arthrobotrys oligospora (Z-2.2∗)
264 ± 196
8 ± 25.9
96.7
h, i,j, k
A. oligospora (Z-2.6)
520 ± 253.7
228 ± 101.63
56.2
d, e, f, l
A. oligospora (Z-2.8∗)
520 ± 253.7
36 ± 47.9
93.1
a, c, g, h, j,j, k
A. oligospora (Z-2.10∗)
520 ± 253.7
28 ± 53.5
94.6
a, g, h, i,j, k
A. musiformis (Z-2.5)
296 ± 92.7
112 ± 114.39
62.2
a, b, c, d, e, f, g, i
A. kirghizica (Z-2.9)
264 ± 196
168 ± 117.45
36.4
a, b, c, d, e, f, l
Trichothecium globosporum var. (Z-2.3)
264 ± 196
108 ± 86.5
59.1
a, b, c, d, e, f, g, i
Trichothecium globosporum  (Z-2.11)
264 ± 196
140 ± 80.55
46.9
a, c, d, e, f, g, l
T. globosporum  var.roseum (Z-3.4)
264 ± 196
76 ± 36.02
71.2
a, b, c, d, e, g, i, k
n
=
10, SD = standard deviation
∗Isolates selected based on their easy mass production and DNA extraction for molecular analysis.
Different literals in each column indicate statistic differences in the larvae reduction among the different fungal isolates (p
≤
0.01).
a, b, c, d, e, f, and g literals indicate the significant comparison between means.
Statistical significance criteria were based on the Kruskal-Wallis analysis.


### 3.3. Molecular Analysis

#### 3.3.1. Homology Analysis

The five isolates of nematophagous fungi that showed the highestin vitro predatory capability from among the 18 isolates, identified by traditional taxonomy asA. oligospora (Z-2.2, Z-2.8, and Z-2.10)M. eudermatum (Z-2.7, Z-3.3), were subject to molecular identification. PCR amplified fragments for Z-2.7 and Z-3.3 isolates were of ≈600 bp length; fragments for Z-2.2, Z-2.8, and Z-2.10 isolates were observed near to 700 bp length. The sequences derived from these PCR products were subject to a homology analysis, and percent identities of 98-99% obtained when the 18S partial gene, ITS1 region, ribosomal 5.8S gene, ITS2 region, and ribosomal 28S partial gene were compared. The morphometric taxonomic identification of the five selected nematophagous fungi was confirmed by the data from the molecular analysis (Table 5).Table 5
Identity percentages with the best hit of sequences containing 18S, ITS1, 5.8S, and ITS2 regions of five nematophagous fungi isolated in Morelos, Mexico, when they are used to do a homology search with BLASTn database. Accession numbers are given in brackets.
Nematophagous fungi sequence identification
Fungi identification
Similarity percentages
Z-2.2/A. oligospora
98% withA. oligospora (AM412778)
Z-2.7/M. eudermatum
98% withM. eudermatum (AF106530)
Z-2.8/A. oligospora
99% withA. oligospora (AM412778)
Z-2.10/A. oligospora
99% withA. oligospora (AM412778)
Z-3.3/M. eudermatum
98% withM. eudermatum (AF106530)


#### 3.3.2. Transitions, Transversions and Insertions

Transitions occur more frequently than transversions; therefore, transitions and transversions were counted, as well as insertions, as a measure of evolutionary distance. The number of transitions between Z-2.2, Z-2.8 and Z-2.10, all classified asA. oligospora, varied from 3 to 5, versus 11 to 16 between this group and Z-2.7 and Z-3.3, classified asM. eudermatum. Only one transversion was identified between Z-2.2, Z-2.8, and Z-2.10, versus 15 to 18 between this group and Z-2.7 and Z-3.3. Three insertions were identified between Z-2.2, Z-2.8, and Z-2.10, versus 18 to 19 between this group and Z-2.7 and Z-3.3. No differences in these three features were detected between Z-2.8 and Z-2.10 and between Z-2.7 and Z-3.3. The lack of differences in this region, which is not subject to selection pressure, suggests that these isolates are the same species or strains (A. oligospora andM. eudermatum, resp.). The correlation between molecular and traditional taxonomy is evident when the frequencies of transitions, transversions, and insertions are compared (Table 6).Table 6
DNA modifications identified in the ITS2 regions of the nematophagous fungal isolates.
Fungal isolates
Kind of change
Z-2.2
A. oligospora
Z-2.7
M. eudermatum
Z-2.8
A. oligospora
Z-2.10
A. oligospora
Z-3.3
M. eudermatum
Z-2.2
A. oligospora
Transitions
0
12
5
3
11
Transversions
0
18
1
1
18
Insertions
0
19
3
3
18
Z-2.7
M. eudermatum
Transitions
0
16
16
0
Transversions
0
15
16
0
Insertions
0
24
24
0
Z-2.8
A. oligospora
Transitions
0
0
16
Transversions
0
0
16
Insertions
0
0
24
Z-2.10
A. oligospora
Transitions
0
16
Transversions
0
16
Insertions
0
24
Z-3.3
M. eudermatum
Transitions
0
Transversions
0
Insertions
0


#### 3.3.3. ITS2 Sequence Analysis

A phenogram was generated using the ITS2 region that included the five selected nematophagous fungi isolates and homologous sequences with different similarity level (Figure3). It shows a clear separation between theArthrobotrys andMonacrosporium genera. The isolates Z-3.3 and Z-2.7 are located in the group ofMonacrosporium (blue in Figure 3), and the isolates Z-2.10, Z-2.8, and Z-2.2 are located within theArthrobotrys group, which is in agreement with the morphological classification. Some organisms classified asArthrobotrys are grouped withMonacrosporium, such as those associated with the accession numbers AY773459, EU3311088, and AY773451, among others; this observation suggests that these organisms could belong to theMonacrosporium genus. Inside of theArthrobotrys group, two main branches are evident, the inferior branch includes the isolates Z-2.10, Z-2.8, and Z-2.2 reported in this work, which share the closest homology with sequences with accession numbers HQ649929 and AY773462, and the organism associated with the last sequence is classified asA. oligospora. Two other sequences also associated withA. oligospora (accession numbers GU799560 and HQ730901) are grouped separately. Similarly, not all the sequences associated withM. eudermatum are grouped together. These observations could indicate that there is a greater diversity than expected after morphological identification.Figure 3
Phenogram generated using the ITS2 region that includes five selected nematophagous fungi isolates and homologous sequences with different similarity level.


## 3.1. Traditional Taxonomy

Morphological and morphometric taxonomic features of different fungal structures are shown in Table2. Eighteen fungal isolates were obtained in total from the different geographic areas and from different substrates. Fungi were identified as belonging to three genera and six species of nematophagous fungi. Most of the fungal isolates corresponded to the speciesM. eudermatum (8 isolates); four isolates were identified asA. oligospora, two asA. musiformis, and one asA. kirghizica. Two isolates were determined asTrichothesium globosporum var.globosporum and one isolate asT. g. var.roseum (Table 2). Different morphological features of the nematophagous fungi are shown in Table 3.Table 2
Morphological and morphometric features of the different genera and species of nematophagous fungi isolated from soil samples from Northern Cuernavaca City, Morelos, Mexico.
Genus/species strain
Conidium size (μm)
Conidiophore length (μm)
Trapping devices
Specific features
Length
Width
M. eudermatum (Z-2.1)
27.5–40
15–20
300
Three-dimensional adhesive nets
One single apical conidium, septum situated at the middle conidium, segmented and branched conidiophore
M. eudermatum (Z-2.4)
27.5–40
15–20
450
M. eudermatum (Z-3.1)
25–42.5
12–25
400
M. eudermatum (Z-2.7)
32.5–45
15–22
350
M. eudermatum (Z-3.3)
27–42.5
12–22
370
M. eudermatum (Z-3.5)
27.5–40
15–22
200
M. eudermatum (Z-3.6)
27.5–40
12–20
300
M. eudermatum (Z-3.8)
26–40
13–21
320
A. oligospora (Z-2.2)
17.5–30
12–20
550
Three-dimensional adhesive nets
Conidiophores present conidia clusters; septum situated below the middle conidium
A. oligospora (Z-2.6)
15–25
10–12
300
A. oligospora (Z-2.8)
17.5–25
10–12
180
A. oligospora (Z-2.10)
20–25
7.5–12
340
A. musiformis (Z-2.5)
25–40
7.5–15
350
Three-dimensional adhesive nets
Elongated and segmented conidia; conidia are produced from small denticles
A. musiformis (Z-3.2)
22.5–35
7.5–12
340
A. kirghizica (Z-2.9)
15–22.5
7.5–12
340
Three-dimensional adhesive nets
Ovoid shape conidia; septum is situated slightly below the middle conidium
T. globosporumvar.globosporum (Z-2.3)
15–30
7.5–15
310
Three-dimensional adhesive nets
Conidium distal cell is larger; conidia are frequently produced in clusters
T. g.var.globosporum (Z-2.11)
22.5–30
10–17
350
T. g.var.roseum(Z-3.4)
10–20
10–20
370
Three-dimensional adhesive nets
Apiculated base, distal cell is longerTable 3
Microphotographs showing the aspect of different nematophagous fungi structures (conidia, conidiophore, and trapping devices) isolated from Morelos, Mexico. Scale bars = 20µm.
Genus/species
Identification number
Conidia
Trapping devices
Conidiophore
Monacrosporium eudermatum
Z-2.1
Z-2.4
Z-2.7
Z-3.1
M. eudermatum
Z-3.3
Z-3.5
Z-3.6
Z-3.8
A. oligospora
Z-2.2
Z-2.6
Z-2.8
Z-2.10
A. musiformis
Z-2.5
Z-3.2
A. kirghizica
Z-2.9
Trichothecium globosporumvar. globosporum
Z-2.3
T. g.var. globosporum
Z-2.11
T. g.var. roseum
Z-3.4


## 3.2.In Vitro Predatory Activity of Nematophagous Fungi

Results of the assessment of thein vitro predatory activity of the nematophagous fungi againstH. contortus infective larvae are shown in Table 4. The larval reduction percentages recorded for the different isolates ranged between 32.4 and 96.7%. The lowest predation percentages were observed forM. eudermatum (isolate Z-2.4) at 32.4% andA. kirghizica (Z-2.9) at 36.4%. SpeciesT. globosporum (Z-2.3, Z-2.11),A. oligospora (Z-2.6),A. musiformis (Z-2.5), andM. eudermatum (Z-2.1) showed moderate predation activities (59.1, 46.9, 56.2, 62.2, and 67.6%, resp.). The highest predatory activities were observed withT. globosporum var.roseum (Z-3.4) at 71.2%,M. eudermatum (Z-3.6) at 74.3%,M. eudermatum (Z-3.1) at 78.4%,M. eudermatum (Z-3.8) at 85.4%,M. eudermatum (Z-3.3) at 91.53%,A. oligospora (Z-2.8) at 93.07%,A. oligospora (Z-2.10) at 94.6%,M. eudermatum (Z-2.7) at 95.9%, andA. oligospora (Z-2.2) at 96.7%.Table 4
Mean numberHaemonchus contortus infective larvae and reduction percentages recovered after 7 days’ in vitrointeraction with different nematophagous fungi isolates and reduction percentages.
Fungal isolates
Control groupX
- (±SD)
Treated groupX
- (±SD)
Larvae reduction (%)
Letters
Monacrosporium eudermatum (Z-2.1)
296 ± 92.7
96 ± 50.6
67.6
a, b, c, d, f, g
M. eudermatum (Z-2.4)
296 ± 92.7
200.1 ± 13.14
32.4
b, c, d, e, f, l
M. eudermatum (Z-3.1)
296 ± 92.7
64 ± 42.99
78.4
a, b, c, d, g, h, i, k
M. eudermatum (Z-2.7∗)
296 ± 92.7
12 ± 19.3
95.9
g, i,j, k
M. eudermatum (Z-3.3∗)
520 ± 253.7
44 ± 78.7
91.5
a, c, g, h, i,j, k
M. eudermatum (Z-3.6)
296 ± 92.7
76 ± 60.1
74.3
a, b, c, d, e, g, i, k
M. eudermatum (Z-3.8)
520 ± 253.7
76 ± 74.1
85.4
a, b, c, d, e, g, i, k
Arthrobotrys oligospora (Z-2.2∗)
264 ± 196
8 ± 25.9
96.7
h, i,j, k
A. oligospora (Z-2.6)
520 ± 253.7
228 ± 101.63
56.2
d, e, f, l
A. oligospora (Z-2.8∗)
520 ± 253.7
36 ± 47.9
93.1
a, c, g, h, j,j, k
A. oligospora (Z-2.10∗)
520 ± 253.7
28 ± 53.5
94.6
a, g, h, i,j, k
A. musiformis (Z-2.5)
296 ± 92.7
112 ± 114.39
62.2
a, b, c, d, e, f, g, i
A. kirghizica (Z-2.9)
264 ± 196
168 ± 117.45
36.4
a, b, c, d, e, f, l
Trichothecium globosporum var. (Z-2.3)
264 ± 196
108 ± 86.5
59.1
a, b, c, d, e, f, g, i
Trichothecium globosporum  (Z-2.11)
264 ± 196
140 ± 80.55
46.9
a, c, d, e, f, g, l
T. globosporum  var.roseum (Z-3.4)
264 ± 196
76 ± 36.02
71.2
a, b, c, d, e, g, i, k
n
=
10, SD = standard deviation
∗Isolates selected based on their easy mass production and DNA extraction for molecular analysis.
Different literals in each column indicate statistic differences in the larvae reduction among the different fungal isolates (p
≤
0.01).
a, b, c, d, e, f, and g literals indicate the significant comparison between means.
Statistical significance criteria were based on the Kruskal-Wallis analysis.


## 3.3. Molecular Analysis

### 3.3.1. Homology Analysis

The five isolates of nematophagous fungi that showed the highestin vitro predatory capability from among the 18 isolates, identified by traditional taxonomy asA. oligospora (Z-2.2, Z-2.8, and Z-2.10)M. eudermatum (Z-2.7, Z-3.3), were subject to molecular identification. PCR amplified fragments for Z-2.7 and Z-3.3 isolates were of ≈600 bp length; fragments for Z-2.2, Z-2.8, and Z-2.10 isolates were observed near to 700 bp length. The sequences derived from these PCR products were subject to a homology analysis, and percent identities of 98-99% obtained when the 18S partial gene, ITS1 region, ribosomal 5.8S gene, ITS2 region, and ribosomal 28S partial gene were compared. The morphometric taxonomic identification of the five selected nematophagous fungi was confirmed by the data from the molecular analysis (Table 5).Table 5
Identity percentages with the best hit of sequences containing 18S, ITS1, 5.8S, and ITS2 regions of five nematophagous fungi isolated in Morelos, Mexico, when they are used to do a homology search with BLASTn database. Accession numbers are given in brackets.
Nematophagous fungi sequence identification
Fungi identification
Similarity percentages
Z-2.2/A. oligospora
98% withA. oligospora (AM412778)
Z-2.7/M. eudermatum
98% withM. eudermatum (AF106530)
Z-2.8/A. oligospora
99% withA. oligospora (AM412778)
Z-2.10/A. oligospora
99% withA. oligospora (AM412778)
Z-3.3/M. eudermatum
98% withM. eudermatum (AF106530)


### 3.3.2. Transitions, Transversions and Insertions

Transitions occur more frequently than transversions; therefore, transitions and transversions were counted, as well as insertions, as a measure of evolutionary distance. The number of transitions between Z-2.2, Z-2.8 and Z-2.10, all classified asA. oligospora, varied from 3 to 5, versus 11 to 16 between this group and Z-2.7 and Z-3.3, classified asM. eudermatum. Only one transversion was identified between Z-2.2, Z-2.8, and Z-2.10, versus 15 to 18 between this group and Z-2.7 and Z-3.3. Three insertions were identified between Z-2.2, Z-2.8, and Z-2.10, versus 18 to 19 between this group and Z-2.7 and Z-3.3. No differences in these three features were detected between Z-2.8 and Z-2.10 and between Z-2.7 and Z-3.3. The lack of differences in this region, which is not subject to selection pressure, suggests that these isolates are the same species or strains (A. oligospora andM. eudermatum, resp.). The correlation between molecular and traditional taxonomy is evident when the frequencies of transitions, transversions, and insertions are compared (Table 6).Table 6
DNA modifications identified in the ITS2 regions of the nematophagous fungal isolates.
Fungal isolates
Kind of change
Z-2.2
A. oligospora
Z-2.7
M. eudermatum
Z-2.8
A. oligospora
Z-2.10
A. oligospora
Z-3.3
M. eudermatum
Z-2.2
A. oligospora
Transitions
0
12
5
3
11
Transversions
0
18
1
1
18
Insertions
0
19
3
3
18
Z-2.7
M. eudermatum
Transitions
0
16
16
0
Transversions
0
15
16
0
Insertions
0
24
24
0
Z-2.8
A. oligospora
Transitions
0
0
16
Transversions
0
0
16
Insertions
0
0
24
Z-2.10
A. oligospora
Transitions
0
16
Transversions
0
16
Insertions
0
24
Z-3.3
M. eudermatum
Transitions
0
Transversions
0
Insertions
0


### 3.3.3. ITS2 Sequence Analysis

A phenogram was generated using the ITS2 region that included the five selected nematophagous fungi isolates and homologous sequences with different similarity level (Figure3). It shows a clear separation between theArthrobotrys andMonacrosporium genera. The isolates Z-3.3 and Z-2.7 are located in the group ofMonacrosporium (blue in Figure 3), and the isolates Z-2.10, Z-2.8, and Z-2.2 are located within theArthrobotrys group, which is in agreement with the morphological classification. Some organisms classified asArthrobotrys are grouped withMonacrosporium, such as those associated with the accession numbers AY773459, EU3311088, and AY773451, among others; this observation suggests that these organisms could belong to theMonacrosporium genus. Inside of theArthrobotrys group, two main branches are evident, the inferior branch includes the isolates Z-2.10, Z-2.8, and Z-2.2 reported in this work, which share the closest homology with sequences with accession numbers HQ649929 and AY773462, and the organism associated with the last sequence is classified asA. oligospora. Two other sequences also associated withA. oligospora (accession numbers GU799560 and HQ730901) are grouped separately. Similarly, not all the sequences associated withM. eudermatum are grouped together. These observations could indicate that there is a greater diversity than expected after morphological identification.Figure 3
Phenogram generated using the ITS2 region that includes five selected nematophagous fungi isolates and homologous sequences with different similarity level.


## 3.3.1. Homology Analysis

The five isolates of nematophagous fungi that showed the highestin vitro predatory capability from among the 18 isolates, identified by traditional taxonomy asA. oligospora (Z-2.2, Z-2.8, and Z-2.10)M. eudermatum (Z-2.7, Z-3.3), were subject to molecular identification. PCR amplified fragments for Z-2.7 and Z-3.3 isolates were of ≈600 bp length; fragments for Z-2.2, Z-2.8, and Z-2.10 isolates were observed near to 700 bp length. The sequences derived from these PCR products were subject to a homology analysis, and percent identities of 98-99% obtained when the 18S partial gene, ITS1 region, ribosomal 5.8S gene, ITS2 region, and ribosomal 28S partial gene were compared. The morphometric taxonomic identification of the five selected nematophagous fungi was confirmed by the data from the molecular analysis (Table 5).Table 5
Identity percentages with the best hit of sequences containing 18S, ITS1, 5.8S, and ITS2 regions of five nematophagous fungi isolated in Morelos, Mexico, when they are used to do a homology search with BLASTn database. Accession numbers are given in brackets.
Nematophagous fungi sequence identification
Fungi identification
Similarity percentages
Z-2.2/A. oligospora
98% withA. oligospora (AM412778)
Z-2.7/M. eudermatum
98% withM. eudermatum (AF106530)
Z-2.8/A. oligospora
99% withA. oligospora (AM412778)
Z-2.10/A. oligospora
99% withA. oligospora (AM412778)
Z-3.3/M. eudermatum
98% withM. eudermatum (AF106530)


## 3.3.2. Transitions, Transversions and Insertions

Transitions occur more frequently than transversions; therefore, transitions and transversions were counted, as well as insertions, as a measure of evolutionary distance. The number of transitions between Z-2.2, Z-2.8 and Z-2.10, all classified asA. oligospora, varied from 3 to 5, versus 11 to 16 between this group and Z-2.7 and Z-3.3, classified asM. eudermatum. Only one transversion was identified between Z-2.2, Z-2.8, and Z-2.10, versus 15 to 18 between this group and Z-2.7 and Z-3.3. Three insertions were identified between Z-2.2, Z-2.8, and Z-2.10, versus 18 to 19 between this group and Z-2.7 and Z-3.3. No differences in these three features were detected between Z-2.8 and Z-2.10 and between Z-2.7 and Z-3.3. The lack of differences in this region, which is not subject to selection pressure, suggests that these isolates are the same species or strains (A. oligospora andM. eudermatum, resp.). The correlation between molecular and traditional taxonomy is evident when the frequencies of transitions, transversions, and insertions are compared (Table 6).Table 6
DNA modifications identified in the ITS2 regions of the nematophagous fungal isolates.
Fungal isolates
Kind of change
Z-2.2
A. oligospora
Z-2.7
M. eudermatum
Z-2.8
A. oligospora
Z-2.10
A. oligospora
Z-3.3
M. eudermatum
Z-2.2
A. oligospora
Transitions
0
12
5
3
11
Transversions
0
18
1
1
18
Insertions
0
19
3
3
18
Z-2.7
M. eudermatum
Transitions
0
16
16
0
Transversions
0
15
16
0
Insertions
0
24
24
0
Z-2.8
A. oligospora
Transitions
0
0
16
Transversions
0
0
16
Insertions
0
0
24
Z-2.10
A. oligospora
Transitions
0
16
Transversions
0
16
Insertions
0
24
Z-3.3
M. eudermatum
Transitions
0
Transversions
0
Insertions
0


## 3.3.3. ITS2 Sequence Analysis

A phenogram was generated using the ITS2 region that included the five selected nematophagous fungi isolates and homologous sequences with different similarity level (Figure3). It shows a clear separation between theArthrobotrys andMonacrosporium genera. The isolates Z-3.3 and Z-2.7 are located in the group ofMonacrosporium (blue in Figure 3), and the isolates Z-2.10, Z-2.8, and Z-2.2 are located within theArthrobotrys group, which is in agreement with the morphological classification. Some organisms classified asArthrobotrys are grouped withMonacrosporium, such as those associated with the accession numbers AY773459, EU3311088, and AY773451, among others; this observation suggests that these organisms could belong to theMonacrosporium genus. Inside of theArthrobotrys group, two main branches are evident, the inferior branch includes the isolates Z-2.10, Z-2.8, and Z-2.2 reported in this work, which share the closest homology with sequences with accession numbers HQ649929 and AY773462, and the organism associated with the last sequence is classified asA. oligospora. Two other sequences also associated withA. oligospora (accession numbers GU799560 and HQ730901) are grouped separately. Similarly, not all the sequences associated withM. eudermatum are grouped together. These observations could indicate that there is a greater diversity than expected after morphological identification.Figure 3
Phenogram generated using the ITS2 region that includes five selected nematophagous fungi isolates and homologous sequences with different similarity level.


## 4. Discussion

Nematophagous fungi are considered important tools for the control of parasitic diseases in both the agricultural and livestock industries worldwide. The present research shows evidence about the rich mycobiota of nematophagous fungi in soil from the studied area at Morelos, Mexico. The obtained fungal isolates show significant potential for use againstHaemonchus contortus, which is one of the most pathogenic parasitic nematodes of ruminants around the world. The fungal species most commonly found wasM. eudermatum. Some species of theMonacrosporium genus have been considered potential agents for control of animal parasitic diseases [24–26]. Another species we identified wasA. oligospora which is the most intensively studied nematophagous species worldwide and a number of important research works are published every year about its use against parasitic nematodes of either plants [27, 28] or animals [29, 30]. A drastic structural change in the nematode cuticle when exposed toA. oligospora has been attributed to proteolytic enzyme activity produced by this species [31]. In addition, secondary metabolites with antibacterial activity produced byA. oligospora are currently being identified [14], giving this species additional biotechnological value. Knowledge of metabolic and chemical aspects of nematophagous fungi has led to important advances, and vanguard proposals on the genomics and proteomics ofA. oligospora are currently being generated [32].Another fungal species we identified wasT. globosporum. In an extensive search on this species, only limited information could be found. Filtrates of different nematophagous fungi, includingT. globosporum, were assessed to identify a possible hatching inhibition effect against the eggs ofMeloidogyne incognita. More than 77% ofM. incognita eggs were inhibited from hatching by theT. globosporum filtrate [20]. In the case ofA. musiformis, several studies have been performed with this species and other genera/species. High capture and destruction efficiency against phytonematodes has been recorded, particularly against the causal agent of root gall, the nematodeM. incognita, in tomato plants [33]. Similarly, this species has been assessed with other fungal species, and a high predatory activity (>97%) againstH. contortus has been recently recorded [34]. In general, the predatory activity of the nematophagous fungi obtained in the present study was variable. SomeA. oligospora isolates (Z-2.2 and Z-2.10), as well asM. eudermatum (Z-2.7) isolates, had high predatory activity (>94%). Variable results have been reported with differentA. oligospora isolates by other researchers. For example, a 93.3% rate ofin vitro predatory activity was found againstH. contortus infective larvae for anotherA. oligospora isolate [34]; however, very low predation values (37%) were obtained against the same parasitic nematode with anotherA. oligospora isolate [35]. In contrast, there is little information about the predatory activity ofM. eudermatum, especially against ruminant parasitic nematodes.In anin vitro trial,M. eudermatum,A. oligospora, andA. robusta conidia were mixed with sheep faeces containingH. contortus eggs; 95.7%, 98.3%, and 10.1% larval reduction was recorded after 15 days’ interaction [29]. Other nematophagous fungal species belonging to theMonacrosporium genus have been evaluated against sheep parasitic nematodes. Thein vitro predatory activity ofM. gephyropagum andA. robusta againstStrongyloides papillosus resulted in 93.1% and 32.3% larval reductions, respectively [26]. Regarding the activity ofA. musiformis, the results of our present research indicate a moderatein vitro larval reduction efficiency (62.1%). In another work, this species showed a low predatory activity (40%) against the same parasite/stage [35]. In contrast, this species has been previously reported to have high predatory activity (97%) against the same parasite/stage [34]. On the other hand,A. kirghizica showed a lowin vitro predatory activity (36.4%) against the assessed nematode in the present study. This result is similar (35%) to the one obtained in another study with the same species [35]. The predatory activity ofA. kirghizica species has been evaluated as a biological control agent against agricultural nematode pests. This species was assessed with other genera/species of nematophagous fungi, that is,Trichothecium pravicoli andA. doliformis under field conditions. The three fungi showed high predatory activity against cucumber parasitic nematodes and root-gall nematodes were substantially reduced (90%) [33].The present research shows a large variability in predatory capability not only among nematophagous fungal species but also among different isolates from the same genera/species. Such a phenomenon could be due to nematophagous fungi possessing a natural “dual” feeding habit in which fungi can develop in two ways: either by feeding on vegetation residues where they obtain their Nitrogen source (N2) from both plants and detritus in soil [12] or by obtaining their nutritional requirements from live nematodes. In the latter case, they become predator fungi [2]. The biological environments of nematophagous fungal isolates could differ in many respects, particularly the presence or absence of nematodes and the nematode population size. Such facts could determine the predatory or saprophytic behaviour of nematophagous fungi. Likewise, the five nematophagous fungi that were selected based on their highin vitro predatory activity (A. oligospora, isolates Z-2.2, Z-2.8, and Z-2.10 andM. eudermatum isolates Z-2.7 and Z-3.3) were sequenced and showed 100% identity toA. oligospora andM. eudermatum sequences from the GenBank-NCBI database. Phylogenetic and comparative analyses showed that the Z-2.10 and Z-2.8 isolates are virtually identical; aside from one transition and one insertion in the ITS1 region, no other change was found in either of the other two analyzed regions. However, these two sequences differ from the Z-2.2 isolate that came from the same site. It is important to highlight such differences because the three isolates are apparently identical from the morphological point of view and in agreement with the consulted taxonomic identification keys [2, 6, 7]. One possible explanation for the fact that some species of fungi are morphologically identical but show some molecular differences is that variants such as subspecies or varieties of fungi are under evolutionary processes influenced by environmental factors. The phylogram shows that two other sequences obtained from GenBank (http://www.ncbi.nlm.nih.gov/genbank/) coincide with the isolates selected in the present study. TheA. oligospora Spanish strain (accession number HQ649929.1) is identical to the Mexican isolates (Z-2.10 and Z-2.8) from the present study; however, the Chinese strain (accession number AY773462.1) shows differences that confirm the presence of genetic diversity. A comparison of morphological features between these two isolates could be interesting to confirm the morphological variability. When conidia and conidiophores from Z-2.10 and Z-2.2 were compared under the microscope variability in length and width was observed, being Z-2.2 conidia wider and longer than Z-2.10, which confirm the morphological variation between the same fungal species. A more detailed molecular analysis must be carried out to achieve a more accurate classification of isolates from different parts of the world. The twoM. eudermatum isolates were classified in a different group and were separated from other reported sequences for the same species (accession numbers AF106530 and AF106528) isolated from Germany. Similar toA. oligospora, there are differences in the DNA sequences that allow a more accurate classification of fungi. The results suggest that the current classification system (based on morphological characteristics) could mask a greater diversity among isolates belonging to theArthrobotrys andMonacrosporium genera. Additionally, the assignation of isolates from different parts of the world to various species based only on morphological characteristics could produce inaccurate results. Similar studies on fungal predatory activity against ectoparasites includingBeauveria,Verticillium,Paecilomyces,Tolypocladium, andMetarhizium suggest that there is intrageneric variation that has been associated with biological processes due to mitochondrial DNA (mtDNA) polymorphisms [36]. The polymorphism that we identified inArthrobotrys andMonacrosporium in the present study could be attributed to differences in biological functions, unspecific mutations, evolutionary processes, predatory behaviour, and so forth. An example of such variability can be observed in the differences in trapping devices developed by the different nematophagous fungal strains [32, 37]. During the evolution ofA. musiformis, a point mutation in the β-tubulin gene produces an altered tubulin protein, whose major function is determining fungal cell structure. This mutation inA. musiformis affects fungal viability when exposed to anthelmintic compounds [37] and this mutation has been associated withA. musiformis polymorphism. In contrast, this mutation has not been found in other species of the genusArthrobotrys. This fact suggests that diversity is present in the genus.The present study identified a number of fungal isolates from Morelos, Mexico, with potential as biological control candidates for the livestock industry. The precise identification of nematophagous fungal isolates using traditional taxonomic methods complemented with molecular tools and characterization of biological activity against livestock parasites and ability to adapt to diverse laboratory conditions are vital features that must be considered in identifying elite strains to be used as potential candidates for the control of parasites in livestock.



---

*Source: 8502629-2016-05-09.xml*