Efﬁcacy of Entomopathogenic Trichoderma Isolates against Sugarcane Woolly Aphid, Ceratovacuna lanigera Zehntner (Hemiptera: Aphididae)

: Ceratovacuna lanigera Zehntner (Hemiptera: Aphididae) is a destructive insect pest of sugarcane that is responsible for reducing the quality, yield, and sugar content. Ecofriendly methods of pest control, such as the use of biological control, is encouraged since it does not exhibit the harmful effects that are usually seen when using chemical-based pesticides. Here, we evaluated the efﬁcacy of several entomopathogenic fungi isolates against C . lanigera under laboratory conditions. A total of seven isolates were obtained from agricultural soils that were collected from various locations from Sabah, Malaysia, using the insect bait method. Among the seven, four isolates (TMS114, TMS623, TMS628, and TMS707) were positively identiﬁed as T. harzianum (Hypocreales: Hypocreaceae), two isolates (TMS45, TMS321) as T. asperellum (Hypocreales: Hypocreaceae), and one isolate (TMS633) as T. hamatum (Hypocreales: Hypocreaceae) based on a combination of morphological characteristics and molecular data using ITS 1-5.8 S-ITS 2 regions of the rDNA. Four different concentrations (1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 mL − 1 ) of conidia suspensions of the isolates were applied separately on nymphs and adults of C. lanigera. The results showed that mortality in both nymphs and adults was dependent on the conidia concentration and exposure time. Among the four concentrations, the highest mortality of C . lanigera was found to range from 59.99 to 75.70% for nymphs after 72 h and from 57.45 to 72.31% for adults after 6 days of treatments at the concentration of 1 × 10 8 conidia mL − 1 . Trichoderma harzianum (TMS623) showed signiﬁcantly highest mortality of 75.70% for nymphs and 72.31% for adults among the seven isolates. However, the nymph and adult mortality rates for the isolates were 70.00% & 68.00% for TMS114, 66.42% & 63.10% for TMS628, 67.85% & 65.24% for TMS707, 61.42% & 58.80% for TMS45, 59.99% & 57.45% for TMS321 and 63.56% & 60.91% for TMS633, respectively. For nymph mortality, T. harzianum (TMS623) showed a LC 50 (lethal concentration) value of 6.30 × 10 5 conidia mL − 1 and LC 90 value of 3.01 × 10 9 conidia mL − 1 , respectively, and for adult mortality, the values were found of 6.16 × 10 5 conidia mL − 1 for LC 50 and 2.23 × 10 10 conidia mL − 1 for LC 90 , respectively. The lethal time (LT) values for nymph and adult mortality were found of 42.65 h and 3.89 days for LT 50 and 93.32 h and 8.70 days for LT 90 by T. harzianum (TMS623). Therefore, T. harzianum (TMS623) showed more pathogenic potential against C. lanigera and can be further applied to determine its efﬁcacy under ﬁeld conditions. To our best knowledge, this is the ﬁrst report of T. harzianum, T. hamatum , and T. asperellum species that are used as entomopathogenic fungi against sugarcane woolly aphid, C. lanigera .


Introduction
Sugarcane (Saccharum officinarum L.) is an important cash crop and it is the main source of processed sugar. Sugarcane woolly aphid (SWA), Ceratovacuna lanigera Zehntner (Hemiptera: Aphididae), is one of the most destructive insect pests of sugarcane [1]. Adult and nymphs of the sugarcane woolly aphid gather on the underside of leaves and suck the sap from the leaves. They produce honeydew that covers the entire upper surface of the leaf which causes the growth of sooty mold. With continuous infestation and sap sucking from the leaves, the plants will wilt and dry up from the tip downwards. This will eventually lead to the reduction of the length, girth, and weight of the sugarcane, as well as the quality, yield, and sugar content [2,3].
Chemical pesticides can leave poisonous residues in food products and harm nontarget organisms in the environment. Problems with chemical pesticides have led to a global interest in developing alternatives to these chemicals. The use of entomopathogenic fungi (EPFs) as bio-pesticides is a promising alternative to chemical-based pest control. Fungal biocontrol agents can mitigate and reduce the detrimental effects of chemical pesticides in the agriculture industry [4]. Entomopathogenic fungi have been extensively investigated as biological control agents of insect pests in attempts to develop sustainable crop protection [5].
Insect fungal diseases are common and widespread. More than 700 fungal species of 90 genera have the ability to cause diseases [6,7] at all the development phases of a wide range of insect pests in nursery, storage, greenhouses, and field conditions [8][9][10][11]. The pathogenicity of entomopathogenic fungi depends on the fungal species and the virulence of the strain [12]. Entomopathogenic fungi infect their insect host through the cuticle, penetrate them, and spread throughout the body. Once the fungus kills the host, it can grow from the host's cadaver and produce more spores, thus, increasing the chance of others getting infected [13].
The fungal genus Trichoderma (Hypocreales: Hypocreaceae) are usually cosmopolitan, exhibit high genetic diversity, and are commonly found in all geographical climatic areas of different habitats, including agricultural lands, marshes, and deserts [19][20][21]. Trichoderma species can be easily isolated from natural soil, decaying organic plant material, and wood and classified as imperfect fungi belonging to the order Hypocreales of Ascomycota [22,23].
Since different species of insect pests can be controlled by different species of entomopathogenic fungi (EPFs), this suggests that EPFs are useful for the biological control of insect pests. Therefore, it is hypothesized that virulent isolates of EPFs could potentially be used as agents to control and inhibit the growth of sugarcane woolly aphids (SWA). As such, we attempted to identify these virulent isolates from various agricultural soils that were planted with different crops and to determine its efficacy (with different concentrations of conidia suspensions) against sugarcane woolly aphid nymphs and adults under laboratory conditions.

Collection of Soil Samples
A total of seven soil samples were collected from agricultural soils that were planted with pumpkin, brinjal (two samples from two different locations), okra, sugarcane, maize, and mustard, respectively. The samples were taken from cultivated lands at Tuaran, Menggatal, Penampang, and Papar in Sabah, Malaysia. The soil samples were taken at a depth of approximately 10-15 cm [35] using a sterilized garden shovel after removing debris from the surface. To avoid possible cross-contamination, the garden shovel was sterilized with 70% ethanol immediately after each soil sample was taken. Approximately 500 g of soil were taken randomly from four different sites at each agricultural plot. Next, the four soil samples were mixed together and 2.0 kg of total weight was collected. The samples were placed in labeled polyethylene bags, sealed with an elastic band, and taken to the Microbiology laboratory of the Biotechnology Research Institute of University of Malaysia Sabah (UMS) as soon as possible and stored at 4 • C for subsequent analysis.

Collection of Insect Samples
According to Zimmermann [36], mealworm larvae, Tenebrio molitor L. (Coleoptera: Tenebrionidae) are traditional bait insects for isolation of entomopathogenic fungi by insect bait method as they are easy rearing and susceptible to the fungi. Thus, Tenebrio molitor larvae were obtained from Harvest Fish & Agriculture Sdn. Bhd., Penampang, Sabah, Malaysia and used in the insect bait method. The larvae were then brought to the Microbiology laboratory of BRI, UMS, in clean vials for subsequent assay.

Isolation of Fungi by Insect Bait Method (IBM)
The fungi were isolated from soil samples using the insect bait method (IBM) following by Zimmermann [36]. The soil samples were cleaned from roots and mixed homogeneously, after which an approximately 400 g was placed in a sterilized plastic box (13 cm × 13 cm × 10 cm). A total of four replicates were performed for each soil sample, each labeled according to sample number [37]. Soil samples were moistened with sterile water and then 10 pre-heated mealworm larvae were placed into each box and the larvae were covered with moist soil. The boxes were covered with lids with a mesh cloth and incubated at 22 • C in the dark and inverted every day during the first week for providing maximum contact of spores to larvae. The samples were monitored for the first time after 5 days and this was repeated every 3-4 day for 3 weeks after the initial baiting and the growth of fungi to the larvae was inspected [38].
The larvae that were infected by entomopathogenic fungi were removed from the soil and the fungi were isolated according to Zimmermann [36,38]. The infected larvae were surface-sterilized by soaking 1% (v/v) sodium hypochlorite (NaClO) solution for 3 min, rinsed three times in sterile water, and dried on sterile filter paper.
The surface-sterilized larvae were cut into small pieces (5-7 pieces) and each small piece of larva was placed in a Petri dish (90 mm diameter) containing selective medium, sabouraud dextrose agar + yeast (0.2%) (SDAY) [39]. Streptomycin (0.06%) and tetracycline (0.005%) were used as bacterial inhibitors. The Petri dish was incubated at 28 • C and a humidity of 70%. After 3 days of incubation, the fungal colonies were visible and were then transferred into SDAY medium for fungal culture. After that, the pure cultures were maintained in SDAY medium and subsequently stored at 4 • C for further use.

Morphological Characterization
Morphological characterization of the fungi isolates was carried out following the macroscopic and microscopic characteristics of the culture based on the procedure of Siddiquee [35].

Macroscopic Characterization
To identify the fungi species, several groups of fungi colonies were analyzed by morphological characteristics. The isolates were re-cultured aseptically from the stock slants in SDAY in a Petri dish (90 mm) and incubated at 28 • C, with a humidity of 70% for 5-7 days, and the isolates were monitored every day. The isolates were categorized based on the following macroscopic characteristics: colony growth rate, colony color in both sides of the Petri dish, colony appearance, and the mycelium density [35].

Microscopic Characterization
The slide culture technique was followed for microscopic observation of the isolates according to Siddiquee [35]. The slide culture was made in the Petri dish (90 mm) containing a U-shaped glass rod on the sterile moistened filter paper. A thin film of an agar block (1 cm square) was excised out from the agar plate and transferred to the center of the slide and placed on the glass rod. A total of four edges of the agar block were inoculated with fungal spores or mycelia fragments using a sterile inoculating loop. A sterile cover slip was then placed on the upper surface of the agar cube. Next, the slide was incubated at room temperature for 3 days and then observed under a light microscope (Olympus CX31) with 100× and 400× magnifications. The observations were focused on the size, shape and color of the conidia, the branching pattern of the conidiophores, and the presence of phialides [35]. Fungi genomic DNA was extracted by the cetyltrimethylammonium bromide (CTAB) method according to Cubero [40]. The polymerase chain reaction (PCR) was performed to amplify ITS 1-5.8 S-ITS 2 regions of the rDNA using primers ITS1 (5 -TCC GTA GGT GAA CCT GCG G-3 ) and ITS4 (5 -TCC TCC GCT TAT TGA TAT GC-3 ) [41]. The PCR reactions were performed with an automated thermal cycler (PTC-200 Peltier Thermal Cycler) using 25 µL of reaction mixtures containing 2.5 µL of 10 × PCR buffer, 2.0 µL of 25 mM magnesium chloride (MgCl 2 ), 0.5 µL of 10 mM PCR deoxynucleotides triphosphate mix (dNTPs), 0.25 µL of 5 Units of Taq DNA polymerase, 0.625 µL of each 10 µM primers (ITS1 and ITS4), 16.5 µL of sterilized distilled water (ddH 2 O), and 2 µL of 100 ng/µL template DNA. The ITS 1-5.8 S-ITS 2 regions of the rDNA were amplified according to the following PCR protocols: initial denaturation at 95 • C for 4 min, followed by 39 cycles of denaturation at 95 • C for 45 s, primer annealing at 56 • C for 45 s, and extension at 72 • C for 1 min. The amplification was completed with one additional step of final extension at 72 • C for 10 min.
A total of of 5 µL of the PCR products were electrophoresed in 1% (w/v) of agarose gel in 1 × Tris-borate-EDTA (TBE) buffer (89 mM Tris Base, 89 mM Boric acid, and 2 mM EDTA, pH 8.0) at 70 V for 50 min. The approximate sizes of the amplified fragments were determined by comparing them with a 100 bp DNA ladder (GeneRuler, 100 bp Plus DNA Ladder, 1st BASE). The gel was stained with Ethidium Bromide (EtBr) for 15 min and observed visually under UV light and photographed using the AlphaImager ® HO (Alpha Innotech, TX, United States).

Sequencing and Phylogenetic Analysis
The successfully amplified PCR products were purified using the QIAquick ® PCR Purification Kit (QIAGEN) according to the manufacturer's instructions. DNA sequencing was performed at the Genomic laboratory of the BRI, UMS, using the ABI3130 DNA sequencer. The sequences of the ITS 1-5.8 S-ITS 2 regions were compared with information from public database in NCBI (National Center for Biotechnology Information, Bethesda, MD, USA) (http://www.ncbi.nlm.nih.gov/BLAST/ (accessed on 16 August 2021)). According to the BLAST results, the sequences were taken based on the similarity value ≥97% with known sequences in the NCBI database [42].
Computer-aided alignment of the ITS 1-5.8 S-ITS 2 regions sequence was performed using the BioEdit Sequence Alignment Editor. The alignment was then optimized manually. Single gaps were treated as a fifth nucleotide (A, C, G, T and gaps). All gap positions within the alignment exceeding a single base in length were replaced by question marks. Sequence insertions without homology to any of the other sequences were deleted in the alignment and a single base was left over causing a minimal gap. In this way, the gaps of all sizes were weighted equally corresponding to a hypothetical single evolutionary event. Finally, alignment adjustments were made manually to remove artificial gaps. The aligned sequences were visually inspected and adjustments were made to improve the alignment.
The ITS 1-5.8 S-ITS 2 phylogenetic tree of the rDNA was performed by comparing it with the sequences from NCBI. The trees were constructed according to the unweighted pair-group method based on the arithmetic average (UPGMA) [43] using molecular evolutionary genetics analysis (MEGA X) software version 10 [44]. The UPGMA tree was used as it produced a rooted tree with constant rate of evolution and was suitable for this purpose. The reliability and interior branches and the validity of the trees were tested with 10,000 replications using bootstrap. The distances were computed using the Kimura-2 (K2P) parameter [45].

Fungi Conidia Suspensions Preparation
The conidia suspensions were prepared from selected fungal isolates according to the method that was described by Javed [46] and Umaru [47] with some modifications. The fungi isolates were re-cultured in SDAY medium from pure culture stock slants and incubated at 28 • C for 14 days. The conidia were harvested by carefully scraping the surface with a sterile spatula and were suspended in 10 mL of sterile 0.02% (v/v) Tween-80 solution. The suspension of each plate was placed separately in a 100 mL beaker and mixed thoroughly to obtain a homogeneous mixture. The suspension of the conidia was sieved with a sterile gauze. The density of conidia of the suspension was calculated by using a neubauer haemocytometer under the microscope and adjusted to the concentrations of 1 × 10 8 , 1 × 10 7 , 1 × 10 6 , and 1 × 10 5 conidia mL −1 , and stored at 4 • C for the bioassay experiment.

Sugarcane Woolly Aphids Rearing and Bioassay Chamber Preparation
Sugarcane leaves that were attached with nymphs and adults were collected from natural infestation of sugarcane experimental plots that had been previously set up at the net house of plant transgenic facility site of University Malaysia Sabah (UMS), Kota Kinabalu, Sabah, Malaysia, at a latitude of 6 • 02 7.20 N and longitude of 116 • 07 4.20 E. The nymphs and adults were reared separately in plastic containers (25 cm × 16 cm × 10 cm) at a temperature of 25 • C, with a humidity of 70 ± 10%, and a photoperiod ratio of 2:1 (Light:Dark) [46,47]. As a natural diet, fresh sugarcane leaves were used to feed the aphids properly which were replaced after 24-48 h. The nymphs and adults were collected from rearing containers and used for a bioassay [46,48].
The bioassay chamber was prepared using plastic containers (25 cm × 16 cm× 10 cm) that were covered with a mesh cloth to allow for adequate ventilation. To maintain sufficient moisture in the bioassay chamber, sterilized wet filter paper was used at the bottom of the container. Fresh sugarcane leaves (about 20 cm length) were sterilized with 0.005% sodium hypochlorite for 2 min and then rinsed with sterilized water and placed in the center of each container. The sterilized, moist cotton was attached around the opening end of the leaves to provide moisture and protection of the leaves from drying in the chambers. The adults and nymphs that were attached to the sugarcane leaves were collected from the rearing containers and surface sterilized gently for 1 min with 0.001% sodium hypochlorite in water (v/v) and rinsed with sterilized water and then transferred to fresh leaves, which were kept in the bioassay chambers, using soft camel hairbrush.
A total of 50 nymphs per leaf were used and placed in a single bioassay chamber for each treatment. A total of three replicates were performed for each treatment in three separate chambers. For the adults, the same procedure was followed in separate chambers [46,49].

Efficacy Determination Bioassay
The efficacy of the fungal isolates against C. lanigera nymphs and adults was determined according to the methods that were described by Mukherjee [29], Umaru [47], and Shrestha [50]. To determine the efficacy, four different concentrations (1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 conidia mL −1 ) of conidia suspension of each isolate were applied separately to the nymphs and adults, which were placed on fresh sugarcane leaves in the separate bioassay chamber. A total of 1 mL of each concentration was applied separately to 50 nymphs and 50 adults, respectively. The topical application was conducted uniformly by using a micropipette (Eppendorf, Enfield, CT, USA).
A 0.02% Tween 80 solution was applied in the control chamber. A total of three replications were taken for each treatment including the control. All of the fungi isolates were applied in bioassay chambers and then covered with mesh lids and kept in the Microbiology Laboratory of BRI, UMS, at a temperature of 25 • C, humidity of 70 ± 10%, and a photoperiod ratio of 2:1 (Light:Dark) [46,47]. Furthermore, to provide sufficient moisture to the bioassay chamber, the filter paper that was placed at the bottom of the container was sprayed with sterilized distilled water as needed. Fresh new sugarcane leaves were placed near both sides of the old leaves (after 24-48 h) as a food source so that the nymphs and adults could always have enough food. The mortality was counted from 12 to 72 h at 12 h intervals for the nymphs and 1 to 6 days for the adults at 1 day intervals [49]. The nymph and adult mortality rates were counted separately and death was confirmed only for those that appeared to have fungal growth [51]. The mortality of each treatment was calculated and corrected using Abbott's formula [52].

Determination of Lethal Concentrations and Lethal Times
Lethal concentrations (LC 50 & LC 90 ) and lethal times (LT 50 & LT 90 ) of the fungi isolates were determined against the nymphs and adults. The bioassay was conducted with four different conidia concentrations (Described in Section 2.8). The LC 50 & LC 90 values were calculated with different concentrations of 1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 conidia mL −1 and the mortalities were counted at 72 h for the nymphs and 6 days for the adults. The LT 50 & LT 90 values were calculated using a concentration of 1 × 10 8 conidia mL −1 with different times of 24 to 72 h for nymphs and 2 to 6 days for adults.

Statistical Analysis
The laboratory experiments were conducted according to completely randomized design (CRD) with eight treatments, including controls, and three replicates were performed for each treatment. The mortality data for the nymphs and adults were corrected using Abbott's formula [52]. The final data were analyzed by analysis of variance (ANOVA) and the means were compared by the least significant difference (LSD) test at a significance level of 5% using Statistix 10 software. The lethal concentration (LC 50

Fungi Isolates
A total of 26 fungi isolates were found by the insect bait method from seven agricultural soils of Sabah, Malaysia (Table 1), based on colony colors, colony appearance, growth rate, and mycelium density. Out of the 26 isolates, 7 were selected based on the efficacy test against nymphs and adults of sugarcane woolly aphid. Among the seven isolates, four (TMS114, TMS623, TMS628, and TMS707) were identified as Trichoderma harzianum, two isolates (TMS45 and TMS321) as Trichoderma asperellum, and one isolate (TMS633) as Trichoderma hamatum based on a combination of morphological characteristics and gene sequencing analysis of the ITS 1-5.8 S-ITS 2 regions of the rDNA (Table 1).

Morphological Characterization
The isolates were characterized based on the colony appearance, growth pattern, branching pattern in conidiophores, size of conidia, and phialides as shown in Figure 1.

Trichoderma harzianum:
A total of four isolates (TMS623, TMS114, TMS628, and TMS707) were characterized as T. harzianum in accordance with the following description.
Macroscopic characteristics: Trichoderma harzianum was observed to cause intense coloration of the agar from yellow to dark orange pigments in mature colonies after five days of incubation in SDAY at 28 • C ( Figure 1A). An indistinct sweet coconut odor was detected from the isolates. The mycelium was initially smooth, watery white and sparse, until cottony aerial mycelium was produced.
Microscopic characteristics: The conidiophores were spread to the top in pyramidal fashion and highly branched ( Figure 1C). The phialides were flask shaped and typically short and wide in the middle with an average length of 4 to 6 µm. The conidia generally tend to be globose to subglobose in shape with an average length of 2 to 3 µm and pale green in color.
Trichoderma hamatum: A total of one isolate (TMS633) was characterized as T. hamatum according to the following description: Macroscopic characteristics: The colonies of T. hamatum were observed to grow rapidly and covered the entire whole Petri dishes within four to five days after incubation in SDAY at 28 • C ( Figure 1D). It produced a white-color aerial mycelia with yellow conidiation in the centre of the plate and grown densely.

Microscopic characteristics:
The conidiopores were observed to be highly branched in an irregular pattern ( Figure 1F). The phialidides were formed from short lateral branches at the base of the elongation, were generally ellipsoildal to ovoidal with an average length of 5-7 µm, and a width of 3-4 µm at the widest point. The conidia were observed to be green, smooth-walled, and ellipsoidal. The average length of the conidia obtained from 4 to 5 µm.
Trichoderma asperellum A total of two isolates (TMS45 and TMS321) were characterized as T. asperellum based on the following description.
Macroscopic characteristics: The colonies of the T. asperellum isolates that were grown in SDAY at 28 • C for five days formed one or two concentric rings near the inoculation zone with a dense production of conidia ( Figure 1G). Dark green conidia were observed to grow towards the center and the margin.

Trichoderma harzianum:
A total of four isolates (TMS623, TMS114, TMS628, and TMS707) were characterized as T. harzianum in accordance with the following description.

Macroscopic characteristics:
Trichoderma harzianum was observed to cause intense coloration of the agar from yellow to dark orange pigments in mature colonies after five days of incubation in SDAY at 28 °C ( Figure 1A). An indistinct sweet coconut odor was detected from the isolates. The mycelium was initially smooth, watery white and sparse, until cottony aerial mycelium was produced.
Microscopic characteristics: The conidiophores were spread to the top in pyramidal fashion and highly branched ( Figure 1C). The phialides were flask shaped and typically short and wide in the middle

Microscopic characteristics:
The conidiophores are highly branched and symmetrically arranged ( Figure 1I). Phialides are produced at the tips of the primary, secondary, and tertiary branches, rarely straight along the length of the branches, typically in whorls of 2-4 phialides. The ampuliform were slightly enlarged in the middle, an average length of 8 µm. The conidia was observed as a dark green color, globose to subglobose or ovoidal, in the size of 3-4 µm.

DNA Amplification
The genomic DNA of entomopathogenic Trichoderma isolates was amplified using ITS1 and ITS4 primers. An approximately 600 bp of the ITS 1-5.8 S-ITS 2 regions of rDNA were successfully amplified among seven Trichoderma isolates (Figure 2). straight along the length of the branches, typically in whorls of 2-4 phi form were slightly enlarged in the middle, an average length of 8 µm observed as a dark green color, globose to subglobose or ovoidal, in th

DNA Amplification
The genomic DNA of entomopathogenic Trichoderma isolates w ITS1 and ITS4 primers. An approximately 600 bp of the ITS 1-5.8 S-ITS were successfully amplified among seven Trichoderma isolates (Figure The sequences of the ITS 1-5.8 S-ITS 2 regions were aligned us quence Alignment Editor and a consensus sequence was generated. quences were compared to the NCBI database using BLAST (a http:/www.ncbi.nlm.nih.gov (accessed on 16th August 2021)) to ide based on the ITS 1-5.8 S-ITS 2 regions of the rDNA. All the isolates sho values between 99-100% with known sequences in the NCBI database the seven isolates, four (TMS114, TMS623, TMS628, and TMS707) w harzianum, two isolates (TMS45 and TMS321) as T. asperellum, and one T. hamatum.

Sequence Analysis of the ITS 1-5.8 S-ITS 2 Regions of the rDNA
The sequences of the ITS 1-5.8 S-ITS 2 regions were aligned using the BioEdit Sequence Alignment Editor and a consensus sequence was generated. All of the gene sequences were compared to the NCBI database using BLAST (available online at http://www.ncbi.nlm. nih.gov (accessed on 16 August 2021)) to identify fungi species based on the ITS 1-5.8 S-ITS 2 regions of the rDNA. All the isolates showed the similarity values between 99-100% with known sequences in the NCBI database. Therefore, among the seven isolates, four (TMS114, TMS623, TMS628, and TMS707) were identified as T. harzianum, two isolates (TMS45 and TMS321) as T. asperellum, and one isolate (TMS633) as T. hamatum.

Phylogenetic Analysis
The phylogenetic tree was constructed using MEGA X software based on the sequence data of the ITS 1-5.8 S-ITS 2 regions of the rDNA. The phylogenetic analysis showed that the seven fungi isolates had the similarities to the previously registered existing sequences on the GenBank database ( Figure 3). According to the tree, three major groups were formed. A total of four isolates (TMS114, TMS623, TMS628, and TMS707) were clustered belong to T. harzianum, two isolates (TMS45 and TMS321) were clustered of T. asperellum, and one isolate (TMS633) was clustered of T. hamatum.

Efficacy Evaluation
The efficacy test was performed for the selection of potential Trichoderma isolates against C. lanigera. The data were obtained from bioassay among isolates against nymphs and adults that were found varying degrees of effectiveness.

Efficacy Evaluation against Nymphs
Trichoderma isolates showed different mortality of C. lanigera nymphs during 72 h of treatment with four different conidia concentrations: in 1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 conidia mL −1 . showed that the seven fungi isolates had the similarities to the previously register isting sequences on the GenBank database ( Figure 3). According to the tree, three groups were formed. A total of four isolates (TMS114, TMS623, TMS628, and TM were clustered belong to T. harzianum, two isolates (TMS45 and TMS321) were clu of T. asperellum, and one isolate (TMS633) was clustered of T. hamatum.

Efficacy Evaluation
The efficacy test was performed for the selection of potential Trichoderma is against C. lanigera. The data were obtained from bioassay among isolates against ny and adults that were found varying degrees of effectiveness.
The time intervals also affected the mortality of the nymphs. After 36 h of treatment, the mortality range of the isolates was between 25.33 to 32%. The mortality was 48.95 to 58.60% after 48 h, 56.66 to 68.54% after 60 h, and 59.99 to 75.70% after 72 h of treatment with the same concentration (1 × 10 8 conidia mL −1 ) of the conidia suspension. T. harzianum (TMS623) caused the highest mortality of 75.70%, followed by T. harzianum (TMS114), which caused 70.00% mortality after 72 h treatment. However, the isolates TMS321, TMS45, TMS633, TMS628, and TMS707 caused 59.33 to 67.85% mortality. Significantly higher nymph mortality was found for T. harzianum (TMS623) compared to other isolates (Figure 4).

Efficacy Evaluation against Adults
Data that were obtained from bioassay experiments that were treated with four different concentrations of the conidia suspension showed that no mortality within one day in the treatments and up to two days in the controls. All of the isolates showed different mortality of adults up to six days of treatment ( Figure 5). Different mortality was observed in different concentrations of the conidia suspension. The highest mortality range of the isolates was found from 57.45 to 72.31% in 1 × 10 8 conidia mL −1 (F = 300.93, p < 0.001), while 33.79 to 42.44% mortality was found in 1 × 10 5 conidia mL −1 (F = 47.62, p < 0.001) within the same six-day treatment period ( Figure 5).
The mortality was also dependent on the duration of the application of conidia suspension. The mortality was 8 to 13.33% after two days, 20.80 to 28.17% after three days, 37.25 to 50.35% after four days, 53.85 to 66.44% after five days, and 57.45 to 72.31% after six days with the same concentration of 1 × 10 8 conidia mL −1 . The highest adult mortality of 72.31% was found after six days of treatment in 1 × 10 8 conidia mL −1 suspension of T. harzianum (TMS623), followed by T. harzianum (TMS114) which caused 68.00% mortality. At the same time, other isolates of TMS321, TMS45, TMS633, TMS628, and TMS707 caused 57.45 to 65.24% mortality. Adult mortality was also significantly higher for T. harzianum (TMS623) than for other isolates ( Figure 5).
The mortality was also dependent on the duration of the application of conidia s pension. The mortality was 8 to 13.33% after two days, 20.80 to 28.17% after three da 37.25 to 50.35% after four days, 53.85 to 66.44% after five days, and 57.45 to 72.31% af six days with the same concentration of 1 × 10 8 conidia mL −1 . The highest adult morta of 72.31% was found after six days of treatment in 1 × 10 8 conidia mL −1 suspension of harzianum (TMS623), followed by T. harzianum (TMS114) which caused 68.00% mortal At the same time, other isolates of TMS321, TMS45, TMS633, TMS628, and TMS707 caus Figure 5. The efficacy of seven Trichoderma isolates against C. lanigera adults that were treated with four different concentrations (1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 conidia mL −1 ) of conidia suspension up to six days treatment. The different letters indicate significant differences (p ≤ 0.05) among the treatments. The columns represent the mean ± SEM (n = 3). SEM = standard error mean.

Discussion
The fungi species of genus Trichoderma is a cosmopolitan fungus and common in all types of soil, agricultural lands, manure, decaying tissue, forest, salt marsh, desert soils, and other form of plant organic matter throughout all climatic zones [54,55]. Entomopathogenic fungi (EPFs) that infect insects are natural pathogens and they can be isolated from the soil and infected or dead insects in the field environment and used as a biological control agent [7].
In this study, T. harzianum, T. hamatum, and T. asperellum were isolated from agricultural soils by the insect bait method. However, according to previous research, the EMP fungi species, T. harzianum and T. asperelum were isolated by insect bait method [29].
In our study, the morphological identification of the Trichoderma species was performed based on the macroscopic and microscopic characteristics of the culture [35]. Molecular identification was performed using the sequence of ITS 1-5.8 S-ITS 2 regions of the rDNA with two primers, ITS 1 and ITS 4. Similarly, in some other previous studies, Trichoderma species were identified using ITS 1-5.8 S-ITS 2 regions of the rDNA with primers ITS 1 and ITS 4 [29,34,56].
Current studies demonstrated the efficacy of different Trichoderma isolates against C. lanigera under laboratory conditions. Bioassays to determine the efficacy of conidia suspension of the isolates showed that all Trichoderma isolates that were tested were pathogenic to both nymphs and adults of C. lanigera, although mortality was different among the isolates (Figures 4 and 5). Similarly, Nawaz [28] reported that Trichoderma species caused 73% mortality of the aphid species, Aphis gossypii (Glover) and 53% mortality of cotton leaf hopper (Jassid), Amrasca bigutulla bigutulla (Hemiptera: Cicadellidae) with the conidia suspension at a concentration of 1 × 10 8 conidia mL −1 after 72 h of treatment under laboratory condition.
Trichoderma harzianum has also been described as an effective biological control agent for the control of the cabbage aphid, Brevicoryne brassicae which caused 70% mortality after 72 h of treatment [34].
The virulence potential of entomopathogenic fungi against target aphid populations varies for different isolates and also varies from strain to strain. Similar aphid species can have varying levels of susceptibility to different types of fungal strain [57].
In the current study, the efficacy results of seven isolates belonged to three species of T. harzianum, T. hamatum, and T. asperellum with four different conidia concentrations (1 × 10 5 , 1 × 10 6 , 1 × 10 7 , and 1 × 10 8 conidia mL −1 ) showed that the mortality of both nymphs and adults are dependent on concentration and time. The mortality increased with increasing conidia concentration and exposure time. The maximum nymph and adult mortality due to Trichoderma isolates was higher with the higher concentrations (1 × 10 8 conidia mL −1 ) of the conidia suspension whereas the mortality was lower at a lower concentration (1 × 10 5 conidia mL −1 ) within the same treatment period of 72 h and 6 days (Figures 4 and 5). This indicates that the mortality is effectively controlled by the density of the conidia suspension. More and more mycelium and conidia grow with a higher density of conidia. This statement is strongly consistent with Nasution [30] who stated that a higher density of conidia increases the presence of conidia in the contact with the insect's body, giving the conidia a better chance to stick, germinate, and penetrate into the insect's body.
As an alternative, the maximum mortality was found for nymphs and adults with a maximum treatment period of 72 h and 6 days with the same conidia concentration (1 × 10 8 conidia mL −1 ) (Figures 4 and 5). This indicates that the longer the EPFs remain in the insect body, the more the mortality increases. The damage to the body tissues of insects from fungal infections will be more severe so the mortality rate will be higher. The previous researcher [30] also reported that the longer the contact time between EPFs and host, the more EPFs can infect the host, leading to increase insect mortality.
Lethal concentration (LC) and lethal time (LT) are a yard measure which can be used to determine the efficacy of a pathogen in relation to the concentration of the suspension of conidia and the time it takes for killing the insect that is exposed at a certain concentration. In our study, LC 50 and LC 90 values as well as LT 50 and LT 90 values were determined for the efficacy of the isolates for the mortality of nymphs and adults of C. lanigera. Among the seven isolates, T. harzianum (TMS623) showed lower LC and LT values compared to other isolates (Tables 2 to 5). Similarly, LC 50 values from 1.15 × 10 6 to 6.76 × 10 6 conidia mL −1 were found for the entomopathogenic fungus, B. bassiana against four aphid species, Brevicoryne brassicae, Schizaphis graminum, Lipaphis erysimi (Hemiptera: Aphididae), and Rhopalosiphum maidis (Hemiptera: Aphididae) [57]. In another previous study, an LC 50 value of 6.75 × 10 6 conidia mL −1 was found for A. flavus (Eurotiales:Trichocomaceae) against seed bug, Elasmolomus pallens (Dallas) (Hemiptera: Rhyparochromidae) [47]. An LT 50 value of 5.54 days was found against Aphis craccivora (Hemiptera: Aphididae) due to entomopathogenic fungi, M. anisopliae [61]. According to Mweke [62], the LT 50 values range of the 23 isolates of M. anisopliae, B. bassiana, and Isaria spp. (Hypocreales: Cordycipitaceae) against Aphis craccivora was found from 3.3 to 6.3 days.
According to LC 50 and LT 50 values of T. harzianum (TMS623) indicated that it has the potential candidate against sugarcane woolly aphid, C. lanigera compared to other isolates that were tested in this study.

Conclusions
Among the seven isolates of different Trichoderma species that were tested in this study, T. harzianum (TMS623) was found to have the potential to be a biological control agent against sugarcane woolly aphid, C. lanigera, based on the mortality of the population that was tested. It is expected to be an alternative to chemical insecticides for the sustainable management of C. lanigera. As far we know, it has been reported for the first time that T. harzianum (TMS623) is a potential candidate for the biological control of C. lanigera and could be used as a biopesticide in agricultural crops.