Harmonia sedecimnotata (F.): Predatory potential, biology, life table, molecular characterization, and field evaluation against Aphis gossypii Glover

The ladybird beetle, Harmonia sedecimnotata (F.) was studied in biology, life table, consumption rates, molecular characterization, and field evaluation. The net reproductive rate (R0), based on the age-stage and two-sex life table, was 43.2 eggs/individual. The female adults lived longer (68.1 d) than the male adults (62.9 d). The rate of consumption increased with progress in each stage of development. Compared to the other larval stages of the predator, the fourth stadium consumed most quantities of Aphis gossypii Glover nymphs (Hemiptera: Aphididae) (200.4). Both female (2214.6) and male (1792.4) consumed more prey (nymphs) than larvae. The net rate of consumption was 1458.92 nymphs of melon aphids. There was no variation in the sequences of the two nucleotides out of 583 bp, H sedecimnotata China (EU392410) and India (MG720024). Our investigations demonstrated that inoculative release of 30 or 40 or 50 adults per 100 m2 attained high reduction of aphids (>90%). Thus, it may be recommended the release rate of 40 adults per 100 m2 to suppress the eggplant aphid population. H. sedecimnotata is therefore one of the most promising biological control agents for cotton aphids that can be achieved for instant control through an inoculative release of adults.

Many of the sucking insect pests cause severe damage to the crop by injecting toxins through their salivary secretions 1-3 and serve as vectors for many plant diseases as well. In both field and protected conditions, they are economically important insect pests 1,4 . The aphid is one of the most important sucking pests and in many ways lower yields of agricultural and horticultural crops. High population densities may damage plants by removing plant nutrients to cause withering and death. Aphids inject toxins from their salivary secretion that cause curling, crinkling and malformations of the plant. They cause secondary plant injury by transmitting many plant diseases. Farmers rely on a variety of insecticides to control these pests 5 . Most of them adversely affect natural enemies such as predators, parasitoids, entomopathogens and have also resulted in residue, resurgence, resistance, and unwanted effects on human and non-target organisms 6,7 . Biological control is a good substitute for toxic insecticides that can protect plants, human beings and the environment [8][9][10][11] . Among the most important beneficial insects that feed on sucking insect pests are coccinellid predators 12 . There are several habitats where the coccinellid beetles are found 13 . They feed on sap-sucking insects like aphids, mealybugs, scales, thrips and mites 14 , as well as other soft-body insects 15 .
More than 5,200 species have been described worldwide in the Coccinellidae family 14 . Around 4,000 predatory ladybird beetle species have been reported from the Indo-Pak subcontinent, including >300 species 16 . The harlequin lady beetle, Harmonia axyridis Pallas (Coleoptera: Coccinellidae) in the genus Harmonia Mulsant is one of the most popular biocontrol agents for greenhouse pest that has become a dangerous invader 17 . Very few studies have been carried out in the other species of the Harmonia genus. The sporadic occurrence of Harmonia sedecimnotata (F.) (Coleoptera: Coccinellidae) from India has recently been reported 18 . Most of the H. sedecimnotata occurs in okra, eggplant, chili, tomato, guava, ornamental, etc. 18 . Polyphagous ladybird predators who feed on various aphid prey species on crop plants exhibit biological fitness variations and are therefore likely to vary in their genotypes 19,20 . An important source for biocontrol agents will be the most effective genotypes of a predator species with higher predation potential 21 . In comparison to other insect species, limited molecular studies were carried out in the members of the Coccinellinae subfamily of the Coccinellidae 20,22 . A partial COI gene was identified and the taxonomic relationship of 16 species was analyzed through four subfamilies 23 . The 5′ region of the COI gene is considered informative and has also been commonly used to detect molecular genetic variation within and between species 24 . For assessing the host-associated variations for insect species, the COI gene has also been effectively employed 3,25 . In this investigation, we used the COI gene for molecular characterization of H. sedecimnotata. The purpose of this study is to select aspects of the consumption rate, biology, life table studies of the H. sedecimnotata fed on Aphis gossypii Glover (Hemiptera: Aphididae), field evaluation and its molecular characterization. Hopefully, this study would help to enhance Integrated Pest Management (IPM) strategies by using this predator as a potential biological control agent for aphid management.

Results
Consumption rate. The various larval stages, as well as the adult stage, had significant effects on per-day, percentage and total consumption of A. gossypii. The age-specific consumption rate, k x (mean number of aphids fed by H. sedecimnotata age x) and age-specific net consumption rate, q x (weighted number of prey fed by H. sedecimnotata age x) were formulated by considering sex and stage differentiation (Fig. 1). The net consumption rate (C 0 ), the transformation rate (Q p ) and the total consumption of the cohort were 1458.92, 33.77 and 36473 nymphs of melon aphids respectively ( Table 1).
The rate of consumption per day was significantly different (F = 75.519; df = 5, 54; P < 0.001) (Fig. 2a). In the fourth stadium it was found highest compared to other stadia. The daily consumption rate among adult beetles The k x (age-specific consumption rate), q x (the age-specific net consumption rate) and C y (the cumulative net consumption rate) of Harmonia sedecimnotata fed on Aphis gossypii. The highest values of k x and q x occur at age 86 d and 19 d, respectively. Data in interaction analyzed with least squares means (LSD) and means separated with a standard error of the mean at P ≤ 0.01. In a column means followed by the same letter are not significantly different at P ≤ 0.001 as determined by Tukey's HSD test.

Stage and statistics
was the highest in female adults compared to male adults. Percent consumption was significantly varied (Fig. 2b). Fourth stadium percent consumption was found to be higher than three other stadia. The second stadium was the least percentage of aphids consumption. In adult beetles, aphid intake was found to be the highest in female adults compared with male adults. During the 4th day and the 6th day the first stadium consumed more and least prey respectively (Fig. 3a). At 3rd day the second stadium killed more prey and at 2nd day the lowest. During the 5th day the third stadium consumed most aphids and during the 1st day the least. The fourth stadium preyed most aphids at 9th day and less at 11th day. Male adult fed more aphids on the 69th day and less on 62nd and 68th day (Fig. 3b). Female adult ate most aphids at 69th and 90th day and it was least at 68th day. The most daily consumption rate in adults was higher than in grubs.
The total consumption rate (F = 77.429; df = 5, 54; P < 0.001) was significantly different (Fig. 4). During the first to fourth stadium, a drastic increase in the total consumption rate was found (Fig. 4a). The total number of aphids consumed per larva (±SE) during the larval stage was 230.88 ± 30.98 (Table 1). Compared to the second stadium, the 1st stadium killed more aphids but less than the 3rd and 4th stadium (Fig. 4a). The 2nd stadium had fewer aphids than the 1st, 3rd, and 4th stadia. The 3rd stadium consumed more aphids than 2nd and 1st stadia, but less than the 4th stadium. The 4th stadium killed more prey than any of the preceding stages. Most aphids were consumed by female adults compared to male adults (Fig. 4b). Because they lived longer, both females and males killed more prey than larvae.
Life tables and population parameters. All 85 H. sedecimnotata eggs hatched. Table 2 presents an average developmental period for egg, larva, pupa, longevity of adults and fecundity. It was observed that the predator lay eggs in groups and the duration was around 5.0 d. The average overall development period (±SE) was 29.6 ± 0.62 d. The average number of eggs laid per female was 75.7.  The probability of age-stage survival rate (S xj ) for eggs was 1.0. The probability of S xj was differed between larval stadia and also differed between adults. The overlap between different stages showed that more variability was found among individuals in the developmental rates (Fig. 5). The age-specific fecundity (m x ), age -specific survival (l x ), and age -specific maternity (l x m x ) ranged from 0.04 to 0.88, 1.62 to 10.71 and 0.84 to 6.00, respectively (     (Table 3). During two seasons, a gradual decrease in aphid infestation was observed in all treatments. A percentage reduction in an aphid population in both seasons influenced by release rates. The aphid infestation in both seasons decreased at the release rate of 40 and 50 adults. After 2 releases of 30 or 40 or 50 adults in both seasons, percentage reduction in an aphid population reached >90%. After second release in release rate 40 and 50 adults were the highest percentage reduction in the aphid population in both season 1 and 2.   Only bootstrap values larger than 60% are shown. The sequences in diamond was generated in this study and deposited in the NCBI database.

Discussion
By feeding on A. gossypii alone, the coccinellid predator (H. sedecimnotata) could reproduce. Both larvae and adults were voracious feeders on all stages of cotton aphids. The prey consumption values were significantly higher with the advancement in each larval stage except the second stadium, which consumed fewer aphids. An increase in feeding potential with age of the larvae was demonstrated in Cryptolaemus montrouzieri Mulsant (Coleoptera: Coccinellidae) 26 which is in accordance with our findings. Compared to the other larval stages of the predator, the fourth stadium consumed most aphids. These results are consistent with earlier findings that reported more than 653.25 Cervaphis quercus (Takahashi) (Hemiptera: Aphididae) 27 , 88.75 Brevicoryne brassicae (L.) (Hemiptera: Aphididae) 28 and 92.13 Tuberculatus nervatus Chakrabarti and Raychaudhuri (Hemiptera: Aphididae) 29 were consumed in the fourth stadium of Harmonia dimidiata (F.) (Coleoptera: Coccinellidae). In the later larval stadia, the consumption rate increased and continued in the adult stage, resulting in one individual being able to kill 4338.6 aphids throughout life.
Most aphids were consumed by the females (F = 77.429; df = 5, 54; P < 0.001) compared to male and larval stages. This is in line with the earlier finding that Cheilomenes sexmaculata (F.) (Coleoptera: Coccinellidae) female killed a maximum of 1624, Aphis craccivora Koch (Hemiptera: Aphididae) than larval stages and male 30 . The consumption rate of adults in the current investigation was higher than the larval stages. Similarly, when fed on A. craccivora, the feeding effectiveness of C. sexmaculata adults was higher than that of larvae 30 . Using the life tables due to non-predatory stages (eggs and pupae), two gaps in consumption rate were depicted. This will determine the life tables' significance and also helped to determine the release time in a biocontrol program.
An overall adult consumption rate was higher than the larva feeding rate because adults were able to live longer in the ecosystems and search for their prey. In addition, more eggs were laid by the female beetles and emerged larvae could search for their prey. C. sexmaculata adults had a higher consumption rate than larvae fed on A. craccivora 30 . The release of coccinellid adults (@ 10 adults/tree or 1500 adults/ha) could suppress the populations of sucking pests 31 . Thus, adults have the potential to be used to suppress sucking pests in biocontrol programs.  www.nature.com/scientificreports www.nature.com/scientificreports/ The female adults lived longer than the male adults in our investigation. A longer female longevity (81.0 d) and male longevity (73.0 d) of H. dimidiata was demonstrated when reared on B. brassicae 28 than our results indicate shorter lifespans fed on A. gossypii by H. sedecimnotata. Our study showed that H. sedecimnotata could reproduce a lifelong fecundity of 1080 eggs. The age-specific fecundity (m x ) and age-specific maternity (l x m x ) showed irregular peaks in female fecundity similar to those shown by Lemnia biplagiata (Swartz) (Coleoptera: Coccinellidae), another coccinellid 32 . The first emerged adult's survival rate (l a ) was 0.65. Since its most important hypothesis, this will help in further statistical analysis.
The development rate among individuals and the role of male individuals in a population is ignored by traditional female age-specific life tables 21 . The stage survival rates and overlaps would not have been found if traditional life tables were used in this study. The life table of Tamarixia radiata Waterson (Hymenoptera: Eulophidae) was recorded using the female age-specific life table 33 . Since the female age-specific life table avoids the male population and could not adequately explain the stage difference, it is possible to obtain 'adult age' in measuring population parameters and no clear F and R 0 relationship. Problems were addressed in adapting female life tables to insect populations 34 . It was possible to observe the stage differentiation of H. sedecimnotata by using the age-stage, two-sex life table and the emergence and termination of various stages. Nevertheless, this gives a single age-specific survival rate when the survival rates from different stages are combined. At an older age, the survival rate can only be equal or lower than at an earlier age (x). It is also possible to observe certain monotonous properties of the l x curve in all studies focused on age-specific life table.
Different methods and parameters were used to estimate population growth potential and insect mass rearing efficiencies. For example, Adalia bipunctata L. (Coleoptera: Coccinellidae) can be grown massively on Encallipterus tillae (L.) and Euceraphis betulae (L.) (Hemiptera: Aphididae) as they were more suitable prey for A. bipunctata based on their rates of development and survival and more abundance in ecosystems 35 . Based on data from the life table, Scymnus apetzi Mulsant (Coleoptera: Coccinellidae), S. subvillosus (Goeze), and Exochomus nigromaculatus (Goeze) (Coleoptera: Coccinellidae) have the potential to be used in the IPM program to manage Hyalopterus pruni (Geoffer) (Hemiptera: Aphididae) 36 .
The mitochondrial genes and the 16 S rDNA gene were used to investigate the insect relationships between the genus level and the level of the species 37 . The mitochondrial genes provide support throughout the topology below the genus level and this will help confirm the species' current positions in coccinellid evolution. In this investigation, 100 percent similarity was found between the sequences of the two nucleotides, H. sedecimnotata China (EU392410) and India (MG720024) out of 583 bp. These patterns might be controlled by the function of the gene and might also be caused by mitochondrial genes 38 . The 16 S rDNA was used to investigate higher taxonomic levels of Anisoptera and Odonata groups 39 . Most important to know about the entire phylogeny of ladybird beetles is the incorporation of morphological characteristics, ecological characteristics and DNA sequence.
Harmonia sedecimnotata is suitable for mass rear in the laboratory, which is an essential requirement for inoculative releases. The most important considerations for inoculation are the number of beetles to be released and the timing of releases. Inoculative releases of H. sedecimnotata showed a higher reduction in the aphid population in the eggplant. Earlier, Rodolia iceryae Janson (Coleoptera: Coccinellidae) was reintroduced as an inoculative release into a plantation where the beetle was absent, and within one year the population of lcerya pattersoni Newstead (Hemiptera: Monophlebidae) was reduced to below the economic threshold 40 . Ladybird beetle predation can slow population growth and reduce aphid densities in the early growing season and the growth of a crop 41 . After 3 releases of 30 Hippodamia varigata (Goeze) (Coleoptera: Coccinellidae) larvae on rose leaves and flower buds, more than 95% aphid reduction was achieved 42 . In this study, the release rates of 30 or 40 or 50 adults on eggplant achieved a high reduction in aphid (>90%). Thus, to suppress the aphid population on eggplant, we recommend the release rate of 40 adults per 100 m 2 .
The life tables and consumption rates were the most important parameters for predator-prey dynamics models, but most pest management based biocontrol programs do not use these life tables. We determined the development rate, survival rate, life tables, and biology of H. sedecimnotata fed on A. gossypii in our investigations, and this predator also had higher consumption rates. In order to know the complete phylogeny of the ladybird beetles, research on ecological character, morphological character and character in the DNA sequence is essential. H. sedecimnotata is a good biocontrol agent that can be used in an inoculative release (40 beetles/100 m 2 ) to control A. gossypii. This biological information revealed that H. sedecimnotata has the potential to be exploited on various crops in India to control sucking insects.

Methods
Aphis gossypii source. An A. gossypii culture was maintained at a temperature of 25 ± 5 °C and 75 ± 5% RH on potted okra plants in a polyhouse of 20 mesh cm −2 . In the experimental field of the ICAR Research Complex for NEH Region, Mizoram Centre, Kolasib, Mizoram, India, Okra (cv. Parbhani Kranti) was grown on a plot sized 10 × 10 m. Okra plants were exposed to natural aphids infestation and on okra plants the aphids began colonizing.
Harmonia sedecimnotata source. Adult beetles were obtained from various host plants and grown on aphids at 23 ± 5 °C temperature and 70 ± 5% RH in the Biological Control Laboratory, ICAR Research Complex for the NEH Region, Mizoram Centre, Kolasib, Mizoram, India. The beetles were fed daily with aphids in plastic containers (7 cm in diameter and 10 cm height) covered by a fine nylon netting for ventilation.
Studies of the consumption rate, biology and life tables. The 10 pairs of males and females that emerged freshly were kept in a plastic container. The plastic containers were secured for ventilation by fine nylon netting and kept at a temperature of 23 ± 1 °C at 70 ± 1% RH. These pairs fed nymphs of melon aphids along with okra leave as food for the nymphs. For the biology and life tables, 85 freshly laid eggs were collected and held in a plastic container under the same conditions. Newly emerging larvae were kept for individual rearing in a new Scientific RepoRtS | (2020) 10:3079 | https://doi.org/10.1038/s41598-020-59809-3 www.nature.com/scientificreports www.nature.com/scientificreports/ plastic container. The first stadium, second stadium, third stadium, fourth stadium, and adults were provided with 30,40,50,75, and 100 nymphs of melon aphids per cup daily. The observations on prey consumption were recorded daily in all larval stages until the stadium changes and in the case of adult beetles until death. Small plastic pipette tubes were kept for oviposition (1.2 cm in diameter, 3 cm in length). The time taken from egg laying to hatching was recorded as well as the duration of larval and pupal stages. Also calculated were the total life cycle, pre-oviposition, oviposition, post-oviposition, fecundity, and adult longevity. The experiment was repeated twice (November to February 2013 and 2014) under laboratory conditions and the mean of both experiments was used to determine life table and consumption rate of H. sedecimnotata.
Analysis of life tables. Based on raw data, the age-stage and two-sex life tables 43,44 were shown using the TWOSEX-MSChart computer program 45 . It calculated the age-stage-specific survival rate (S xj ), the age-stage-specific fecundity (f xj ), the age-specific survival rate (l x ), the age-specific fecundity (m x ), the net reproductive rate (R 0 ), the intrinsic rate of increase (r), the finite rate of increase (λ), and the mean generation time (T) 43 .
Analysis of consumption rate. Total, per-day and percentage of aphid consumption was calculated for each stage of life of predatory beetle. The daily consumption rates were calculated using the CONSUME-MSChart computer program 46 . It calculated the age-specific consumption rate (k x ), the age-stage-specific consumption (C xj ), the age-specific net consumption rate (q x ), the net consumption rate (C 0 ), the cumulative consumption rate (C y ) and the transformation rate (Q p ) 31,47 . Molecular characterization. The DNA was extracted with a slight modification using the cetyl trimethylammonium bromide (CTAB) method 48  Sequences for H. sedecimnotata Indian isolate was deposited in the NCBI database and the accession number was MG720024. Homology searches conducted using the Basic Local Alignment Search Tool (BLAST) and differences in H. sedecimnotata COX-1 sequences were determined using the BioEdit version 7.0.5.3 sequence alignment editor.
The Clustal W program (a MEGA 4.0 sub-program) has been used to align the sequence and where necessary the alignment has been refined 50 . The statistics of molecular character, pairwise nucleotide sequence divergence, and base frequencies for each gene fragment were estimated using MEGA 4.0 software. Neighbor-Joining (NJ) and Minimum-Evolution (ME) for phylogenetic construction implemented in MEGA 4.0. The predator's two releases were made at 30 days apart. Each release was released at five rates of 10, 20, 30, 40, and 50 adults (1 F: 1 M) per plot (100 m 2 ). The first release was made during the vegetative stage when the aphid population was above the economic threshold (10-29% affected plants). The second release covered new leaves and shoots and increased populations in newly emerging aphids nymphs and adults, making it 30 days after the first release. In the same plants/plots, both releases were made. During the study, the control was free to release coccinellids. In order to avoid the movement of predators from one treatment to another, each treatment was separated by a distance of 25 m and no insecticide spray was given throughout the experimental period.

Field evaluation of
Nymph and adult aphid densities were estimated 24 h before each release, and observations of post-treatment were at intervals of fortnight. Densities of melon aphids were determined by counting individuals on the upper and under sides of three leaves from the bottom, middle, and top. Thirty leaf samples (each "leaf sample" consisted of three leaves from bottom, middle, and top) by selecting ten plants from each replication. The plants were randomly examined on each sample date for nymph and adult stages of melon aphids in each plot. Mean aphid population per leaf was estimated by an average of three leaves (bottom, middle and top). The aphid population Scientific RepoRtS | (2020) 10:3079 | https://doi.org/10.1038/s41598-020-59809-3 www.nature.com/scientificreports www.nature.com/scientificreports/ was recorded with the aid of a magnifying glass (MAG, Generic Brand, India). Data were estimated weekly in the canopy area between 8.0 and 10.0 h. For each release, the percentage reduction in aphids over control was calculated.
Statistical analysis. Data were analyzed using SAS Software Version 9.3 51 for biology, and consumption rate with completely randomized design (CRD) and field assessment with RCBD. Data on total, per-day, and percentage consumption for larvae or adults were analyzed using analysis of variance (ANOVA) and the ANOVA was performed on the original values and the means were analyzed by least squares means (LSD) and the mean separated using Tukey's HSD test at P ≤ 0.001. Field assessment data were transformed into square root and arc sine wherever required, as described by Poisson for statistical analysis 52 . ANOVA was done on the transformed values of field evaluation data and the same values were analyzed by least squares means (LSD) and separated using Tukey's HSD test at P ≤ 0.001.