Analyses of α-amylase and α-glucosidase in the malaria vector mosquito, Anopheles gambiae, as receptors of Cry11Ba toxin of Bacillus thuringiensis subsp. jegathesan

Diuresis and α-glucosidase inhibition by erythritol in Aedes aegypti (Diptera: Culicidae) and viability for efficacy against mosquitoes

BACKGROUND

Sugar alcohols, such as erythritol, are low-impact candidates for attractive toxic sugar baits (ATSB) to kill mosquitoes. To determine whether erythritol has a viable future in ATSB formulations, a suite of assays was conducted to diagnose toxicity mechanisms and starvation effects on mortality in Aedes aegypti (L.) as a model system.

METHODS

We measured general carbohydrate load, glucosidase levels, and free glucose in intoxicated adult mosquitoes to observe whether sugar digestion was impaired. We assayed the effects of sugar combinations with erythritol on larvae and adults. To measure erythritol effects when mosquitoes were not resource-deprived, additional assays manipulated the prior starvation status.

RESULTS

Up to 50,000 ppm of erythritol in water had no effect on larvae within 72 h, but an ammonia spike indicated diuresis in larvae as early as 4 h (F8,44 = 22.50, P < 0.0001) after sucrose/erythritol combinations were added. Adult consumption of erythritol was diuretic regardless of the sugar pairing, while sucrose and erythritol together generated above 80% mortality (F2,273 = 33.30, P < 0.0001) alongside triple the normal excretion (F5,78 = 26.80, P < 0.0004). Glucose and fructose paired individually with erythritol had less mortality, but still double the fecal excretion. When ingesting erythritol-laced meals, less sugar was detected in mosquitoes as compared to after sucrose meals (χ2 = 12.54, df = 1, P = 0.0004).

CONCLUSIONS

Data showed that erythritol is a linear competitive inhibitor of α-glucosidase, marking it as a novel class of insecticide in the current research climate. However, the efficacy on larvae was null and not persistent in adult mosquitoes when compared across various starvation levels. Despite significant diuresis, the combined effects from erythritol are not acute enough for vector control programs considering ATSB against mosquitoes.

Culex quinquefasciatus membrane-bound alkaline phosphatase is a putative receptor for Lysinibacillus sphaericus Tpp49Aa1 toxin

The binary toxin Cry48Aa1/Tpp49Aa1 produced by Lysinibacillus sphaericus exhibits potent toxicity against Culicidae larvae. Both Cry48Aa1 and Tpp49Aa1 toxins are crucial for binding to the toxin receptor in Culex quinquefasciatus larvae, albeit with different binding sites. Previous studies have identified Glu71, a membrane-bound α-glucosidase, as a putative binding protein for the Cry48Aa1 toxin, involved in the Cry48Aa1/Tpp49Aa1 toxicity. In this study, we employed pulldown assays to identify a group of Tpp49Aa1-binding proteins from C. quinquefasciatus solubilized midgut brush-border membrane proteins (BBMFs). RNA interference assays revealed that the silencing of an alkaline phosphatase gene (referred to as ALP1263) in C. quinquefasciatus resulted in a significant reduction in larval mortality upon exposure to Cry48Aa1/Tpp49Aa1 toxin in vivo. Furthermore, the ALP1263 protein exhibited specific and high-affinity binding to the Tpp49Aa1 toxin, with a dissociation constant (Kd) of approximately 57.3 nM. The dot blot analysis demonstrated that Tpp49Aa1 C-terminal region was essential for its interaction with the ALP1263 protein. In summary, our findings establish ALP1263 as a functional receptor for Tpp49Aa1 and emphasize its role in the toxicity of Cry48Aa1/Tpp49Aa1.

Culex quinquefasciatus alpha-glucosidase serves as a putative receptor of the Cry48Aa toxin from Lysinibacillus sphaericus

The Cry48Aa/Cry49Aa toxin of Lysinibacillus sphaericus shows specific toxicity towards larvae of Culex spp. Individual Cry48Aa and Cry49Aa subunits interact with distinct target sites in the larval midgut and overcome the resistance of Culex to the Bin toxin. However, the toxin-binding proteins have not yet been identified. The present study aimed to identify Cry48Aa-binding proteins in Culex quinquefasciatus. Pulldown assays using C. quinquefasciatus midgut brush-border membrane fractions (BBMFs) identified a class of proteins, including aminopeptidases (APNs), protease m1 zinc metalloproteases, alkaline phosphatases (ALPs), and maltases, that could be potentially involved in the mode of action of this toxin. RNA interference analysis showed that silenced larvae treated with dsRNA of the alpha-glucosidase (named Glu71) gene were more tolerant of the Cry48Aa/Cry49Aa toxin, which induced less than 20% mortality. The amino acid sequence of Glu71 exhibited 42% identity with Cqm1/Cpm1, which acted as a Bin toxin receptor. Toxin binding assays showed that Cry48Aa had a high specific binding capacity for the Glu71 protein, whereas Cry49Aa exhibited no specific binding. Overall, our results showed that Glu71 is a Cry48-binding protein involved in Cry48Aa/Cry49Aa toxicity.

Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots

BACKGROUND

Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood.

RESULTS

Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein.

CONCLUSIONS

Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

Molecular and biological features of Culex quinquefasciatus homozygous larvae for two cqm1 alleles that confer resistance to Lysinibacillus sphaericus larvicides

BACKGROUND

Culex quinquefasciatus resistance to the binary toxin from Lysinibacillus sphaericus larvicides can occur because of mutations in the cqm1 gene that prevents the expression of the toxin receptor, Cqm1 α-glucosidase. In a resistant laboratory-selected colony maintained for more than 250 generations, cqm1_REC and cqm1_REC-2 resistance alleles were identified. The major allele initially found, cqm1_REC , became minor and was replaced by cqm1_REC-2 . This study aimed to investigate the features associated with homozygous larvae for each allele to understand the reasons for the allele replacement and to generate knowledge on resistance to microbial larvicides.

RESULTS

Homozygous larvae for each allele were compared. Both larvae displayed the same level of resistance to the binary toxin (3500-fold); therefore, a change in phenotype was not the reason for the replacement observed. The lack of Cqm1 expression did not reduce the total specific α-glucosidase activity for homozygous cqm1_REC and cqm1_REC-2 larvae, which were statistically similar to the susceptible strain, using artificial or natural substrates. The expression of eight Cqm1 paralog α-glucosidases was demonstrated in resistant and susceptible larvae. Bioassays in which cqm1_REC or cqm1_REC-2 homozygous larvae were reared under stressful conditions showed that most adults produced were cqm1_REC-2 homozygous (69%). Comparatively, in the offspring of a heterozygous sub-colony reared under optimal conditions for 20 generations, the cqm1_REC allele assumed a higher frequency (0.72).

CONCLUSION

Homozygous larvae for each allele exhibited a similar resistant phenotype. However, they presented specific advantages that might favor their selection and can be used in designing resistance management practices. © 2021 Society of Chemical Industry.

Gene duplication and subsequent functional diversification of sucrose hydrolase in Papilio xuthus

Sucrose is the main product of photosynthesis in plants, providing a rich carbon and energy source for the physiological growth and development of insects. In a previous study, we identified a novel sucrose hydrolase (SUH) in the larval midgut of moths. Intriguingly, there are two copies of Suh, namely Suh1 and Suh2, in several species of butterflies. However, the biochemical characteristics of SUHs in butterflies remain unclear. In this study, we found that this duplication and subsequent diversification produced two Suh genes in Papilio xuthus. These two PxSuh genes were significantly divergent in terms of their expression pattern and enzyme properties. PxSuh messenger RNA expression was highest during the larval stage, reduced in the prepupal and pupal stages and, for PxSuh1, slightly increased again in the adult. The observed levels of PxSuh2 were overall below those of PxSuh1 amongst the development stages examined. Compared with PxSUH2, which has maintained the original gene function of maltose hydrolysis, PxSUH1 exhibits substrate specificity for sucrose with an optimum enzyme activity occurring at an alkaline pH. The data show that PxSuh1 is evolutionarily adapted for effective functioning in an alkaline digestive system. Furthermore, we find that functional diversification of Suh facilitates P. xuthus to digestive carbohydrate of host plants. Thus, our findings offer new insights into the ecological and evolutionary adaptation of digestive enzymes in butterflies.

The physiological differentiation along the midgut of Bombyx mori - inspirations from proteomics and gene expression patterns of the secreted proteins in the ectoperitrophic space

The ectoperitrophic space (EcPS) between the insect midgut epithelial cells and the peritrophic matrix is an unexplored, clean resource for concentrated proteins secreted by the midgut epithelial cells, which offers an ideal opportunity to uncover the midgut functions. In this study, we used Bombyx mori as a model organism and performed comparative proteomic analyses of the secreted proteins in the EcPS at the feeding and wandering stages. A total of 372 proteins were identified from both stages and 70 proteins were predicted to be secreted. Amongst these proteins, 17 secreted digestive proteins were identified and their temporal and spatial transcriptional expression patterns demonstrated that all these proteins were up-regulated at the feeding stage and differentially expressed in different parts of the midgut. Proteins with nutrient reservoir activity and defence activity were found to be up-regulated at the wandering stage. This work is the first to show the presence of digestive enzymes in the EcPS of the insect midgut using a proteomic approach, which provides evidence that suggests a physiological functional differentiation of the insect midgut. It is very clear that the EcPS undergoes dynamic changes in its composition of proteins in response to the changing needs of the insect at different developmental stages.

Effects and mechanisms of Bacillus thuringiensis crystal toxins for mosquito larvae

Bacillus thuringiensis is a Gram-positive aerobic bacterium that produces insecticidal crystalline inclusions during sporulation phases of the mother cell. The virulence factor, known as parasporal crystals, is composed of Cry and Cyt toxins. Most Cry toxins display a common 3-domain topology. Cry toxins exert intoxication through toxin activation, receptor binding and pore formation in a suitable larval gut environment. The mosquitocidal toxins of Bt subsp. israelensis (Bti) were found to be highly active against mosquito larvae and are widely used for vector control. Bt subsp. jegathesan is another strain which possesses high potency against broad range of mosquito larvae. The present review summarizes characterized receptors for Cry toxins in mosquito larvae, and will also discuss the diversity and effects of 3-D mosquitocidal Cry toxin and the ongoing research for Cry toxin mechanisms generated from investigations of lepidopteran and dipteran larvae.

Identification of Cry48Aa/Cry49Aa toxin ligands in the midgut of Culex quinquefasciatus larvae

A binary mosquitocidal toxin composed of a three-domain Cry-like toxin (Cry48Aa) and a binary-like toxin (Cry49Aa) was identified in Lysinibacillus sphaericus. Cry48Aa/Cry49Aa has action on Culex quinquefasciatus larvae, in particular, to those that are resistant to the Bin Binary toxin, which is the major insecticidal factor from L. sphaericus-based biolarvicides, indicating that Cry48Aa/Cry49Aa interacts with distinct target sites in the midgut and can overcome Bin toxin resistance. This study aimed to identify Cry48Aa/Cry49Aa ligands in C. quinquefasciatus midgut through binding assays and mass spectrometry. Several proteins, mostly from 50 to 120 kDa, bound to the Cry48Aa/Cry49Aa toxin were revealed by toxin overlay and pull-down assays. These proteins were identified against the C. quinquefasciatus genome and after analysis a set of 49 proteins were selected which includes midgut bound proteins such as aminopeptidases, amylases, alkaline phosphatases in addition to molecules from other classes that can be potentially involved in this toxin's mode of action. Among these, some proteins are orthologs of Cry receptors previously identified in mosquito larvae, as candidate receptors for Cry48Aa/Cry49Aa toxin. Further investigation is needed to evaluate the specificity of their interactions and their possible role as receptors.

Molecular evolutionary mechanisms driving functional diversification of α-glucosidase in Lepidoptera

The digestive tract of lepidopteran insects is unique given its highly alkaline pH. The adaptive plasticity of digestive enzymes in this environment is crucial to the highly-efficient nutritional absorption in Lepidoptera. However, little is known about the molecular adaptation of digestive enzymes to this environment. Here, we show that lepidopteran α-glucosidase, a pivotal digestive enzyme, diverged into sucrose hydrolase (SUH) and other maltase subfamilies. SUH, which is specific for sucrose, was only detected in Lepidoptera. It suggests that lepidopteran insects have evolved an enhanced ability to hydrolyse sucrose, their major energy source. Gene duplications and exon-shuffling produced multiple copies of α-glucosidase in different microsyntenic regions. Furthermore, SUH showed significant functional divergence (FD) compared with maltase, which was affected by positive selection at specific lineages and codons. Nine sites, which were involved in both FD and positive selection, were located around the ligand-binding groove of SUH. These sites could be responsible for the ligand-binding preference and hydrolytic specificity of SUH for sucrose, and contribute to its conformational stability. Overall, our study demonstrated that positive selection is an important evolutionary force for the adaptive diversification of α-glucosidase, and for the exclusive presence of membrane-associated SUHs in the unique lepidopteran digestive tract.

N-glycosylation influences the catalytic activity of mosquito α-glucosidases associated with susceptibility or refractoriness to Lysinibacillus sphaericus

Cqm1 and Aam1 are α-glucosidases (EC 3.2.1.20) expressed in Culex quinquefasciatus and Aedes aegypti larvae midgut, respectively. These orthologs share high sequence similarity but while Cqm1 acts as a receptor for the Binary (Bin) insecticidal toxin from Lysinibacillus sphaericus, Aam1 does not bind the toxin, rendering Ae. aegypti refractory to this bacterium. Aam1 is heavily glycosylated, contrasting to Cqm1, but little is known regarding how glycosylation impacts on its function. This study aimed to compare the N-glycosylation patterns and the catalytic activities of Aam1 and Cqm1. Mutant proteins were generated where predicted Aam1 N-glycosylation sites (N-PGS) were either inserted into Cqm1 or abrogated in Aam1. The mutants validated four N-PGS which were found to localize externally on the Aam1 structure. These Aam1 and Cqm1 mutants maintained their Bin binding properties, confirming that glycosylation has no role in this interaction. The α-glucosidase activity of both proteins was next investigated, with Aam1 having a remarkably higher catalytic efficiency, influenced by changes in glycosylation. Molecular dynamics showed that glycosylated and nonglycosylated Aam1 models displayed distinct patterns that could influence their catalytic activity. Differential N-glycosylation may then be associated with higher catalytic efficiency in Aam1, enhancing the functional diversity of related orthologs.

Anopheles gambiae Ag55 cell line as a model for Lysinibacillus sphaericus Bin toxin action

Binary toxin (Bin) produced by Lysinibacillus sphaericus is toxic to Culex and Anopheles mosquito larvae. It has been used world-wide for control of mosquitoes that vector disease. The Bin toxin interacts with the glucosidase receptor, Cpm1, in Culex and its orthologue, Agm3, in Anopheles mosquitoes. However, the exact mechanism of its mode of action is not clearly understood. It is essential to understand mode of action of Bin toxin to circumvent the resistance that develops over generations of exposure. A suitable model cell line will facilitate investigations of the molecular action of Bin toxin. Here we report Bin toxin activity on Ag55 cell line that has been derived from an actual target, Anopheles gambiae larvae. The Bin toxin, both in pro and active forms, kills the Ag55 cells within 24h. Bin toxin internalizes in Ag55 cells and also induces vacuolation as tracked by Lysotracker dye. The dose response studies showed that 1.5nM of Bin toxin is sufficient to induce vacuolation and Ag55 cell death. Presence of α-glucosidase gene (Agm3) expression in the Ag55 cells was also confirmed. Thus, Ag55 cells constitute an appropriate model system to decipher the mode of Bin action in mosquito larvae.

The Anopheles-midgut APN1 structure reveals a new malaria transmission-blocking vaccine epitope

Mosquito-based malaria transmission-blocking vaccines (mTBVs) target midgut-surface antigens of the Plasmodium parasite's obligate vector, the Anopheles mosquito. The alanyl aminopeptidase N (AnAPN1) is the leading mTBV immunogen; however, AnAPN1's role in Plasmodium infection of the mosquito and how anti-AnAPN1 antibodies functionally block parasite transmission have remained elusive. Here we present the 2.65-Å crystal structure of AnAPN1 and the immunoreactivity and transmission-blocking profiles of three monoclonal antibodies (mAbs) to AnAPN1, including mAb 4H5B7, which effectively blocks transmission of natural strains of Plasmodium falciparum. Using the AnAPN1 structure, we map the conformation-dependent 4H5B7 neoepitope to a previously uncharacterized region on domain 1 and further demonstrate that nonhuman-primate neoepitope-specific IgG also blocks parasite transmission. We discuss the prospect of a new biological function of AnAPN1 as a receptor for Plasmodium in the mosquito midgut and the implications for redesigning the AnAPN1 mTBV.

A novel sucrose hydrolase from the bombycoid silkworms Bombyx mori, Trilocha varians, and Samia cynthia ricini with a substrate specificity for sucrose

Although membrane-associated sucrase activity has been detected in the midgut of various lepidopteran species, it has not yet been identified and characterized at the molecular level. In the present study, we identified a novel sucrose hydrolase (SUH) gene from the following three bombycoid silkworms: Bombyx mori, Trilocha varians, and Samia cynthia ricini and named them BmSuh, TvSuh, and ScSuh, respectively. The EST dataset showed that BmSuh is one of the major glycoside hydrolase genes in the larval midgut of B. mori. These genes were almost exclusively expressed in the larval midgut in all three species, mainly at the feeding stage. SUHs are classified into the glycoside hydrolase family 13 and show significant homology to insect maltases. Enzymatic assays revealed that recombinant SUHs were distinct from conventional maltases and exhibited substrate specificity for sucrose. The recombinant BmSUH was less sensitive to sugar-mimic alkaloids than TvSUH and ScSUH, which may explain the reason why the sucrase activity in the B. mori midgut was less affected by the sugar-mimic alkaloids derived from mulberry.

Chitosan/DsiRNA nanoparticle targeting identifies AgCad1 cadherin in Anopheles gambiae larvae as an in vivo receptor of Cry11Ba toxin of Bacillus thuringiensis subsp. jegathesan

The Cry11Ba protein of Bacillus thuringiensis subsp. jegathesan crystals has uniquely high toxicity against a spectrum of mosquito species. The high potency of Cry11Ba against Anopheles gambiae is caused by recognition of multiple midgut proteins including glycosyl phosphatidylinositol-anchored alkaline phosphatase AgALP1, aminopeptidase AgAPN2, α-amylase AgAmy1 and α-glucosidase Agm3 that bind Cry11Ba with high affinity and function as putative receptors. The cadherin AgCad2 in An. gambiae larvae also binds Cry11Ba with high affinity (Kd = 12 nM) and is considered a putative receptor, while cadherin AgCad1 bound Cry11Ba with low affinity (Kd = 766 nM), a property not supportive for a Cry11Ba receptor role. Here, we show the in vivo involvement of AgCad1 in Cry11Ba toxicity in An. gambiae larvae using chitosan/DsiRNA nanoparticles to inhibit AgCad expression in larvae. Cry11Ba was significantly less toxic to AgCad1-silenced larvae than to control larvae. Because AgCad1 was co-suppressed by AgCad2 DsRNAi, the involvement of AgCad2 in Cry11Ba toxicity could not be ascertained. The ratio of AgCad1:AgCad2 transcript level is 36:1 for gut tissue in 4th instar larvae. Silencing AgCad expression had no effect on transcript levels of other binding receptors of Cry11Ba. We conclude that AgCad1 and possibly AgCad2 in An. gambiae larvae are functional receptors of Cry11Ba toxin in vivo.

Bacillus thuringiensis toxins: an overview of their biocidal activity

Bacillus thuringiensis (Bt) is a Gram positive, spore-forming bacterium that synthesizes parasporal crystalline inclusions containing Cry and Cyt proteins, some of which are toxic against a wide range of insect orders, nematodes and human-cancer cells. These toxins have been successfully used as bioinsecticides against caterpillars, beetles, and flies, including mosquitoes and blackflies. Bt also synthesizes insecticidal proteins during the vegetative growth phase, which are subsequently secreted into the growth medium. These proteins are commonly known as vegetative insecticidal proteins (Vips) and hold insecticidal activity against lepidopteran, coleopteran and some homopteran pests. A less well characterized secretory protein with no amino acid similarity to Vip proteins has shown insecticidal activity against coleopteran pests and is termed Sip (secreted insecticidal protein). Bin-like and ETX_MTX2-family proteins (Pfam PF03318), which share amino acid similarities with mosquitocidal binary (Bin) and Mtx2 toxins, respectively, from Lysinibacillus sphaericus, are also produced by some Bt strains. In addition, vast numbers of Bt isolates naturally present in the soil and the phylloplane also synthesize crystal proteins whose biological activity is still unknown. In this review, we provide an updated overview of the known active Bt toxins to date and discuss their activities.

Bt toxin modification for enhanced efficacy

Insect-specific toxins derived from Bacillus thuringiensis (Bt) provide a valuable resource for pest suppression. Here we review the different strategies that have been employed to enhance toxicity against specific target species including those that have evolved resistance to Bt, or to modify the host range of Bt crystal (Cry) and cytolytic (Cyt) toxins. These strategies include toxin truncation, modification of protease cleavage sites, domain swapping, site-directed mutagenesis, peptide addition, and phage display screens for mutated toxins with enhanced activity. Toxin optimization provides a useful approach to extend the utility of these proteins for suppression of pests that exhibit low susceptibility to native Bt toxins, and to overcome field resistance.

Molecular evolution of glycoside hydrolase genes in the Western corn rootworm (Diabrotica virgifera virgifera)

Cellulose is an important nutritional resource for a number of insect herbivores. Digestion of cellulose and other polysaccharides in plant-based diets requires several types of enzymes including a number of glycoside hydrolase (GH) families. In a previous study, we showed that a single GH45 gene is present in the midgut tissue of the western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae). However, the presence of multiple enzymes was also suggested by the lack of a significant biological response when the expression of the gene was silenced by RNA interference. In order to clarify the repertoire of cellulose-degrading enzymes and related GH family proteins in D. v. virgifera, we performed next-generation sequencing and assembled transcriptomes from the tissue of three different developmental stages (eggs, neonates, and third instar larvae). Results of this study revealed the presence of seventy-eight genes that potentially encode GH enzymes belonging to eight families (GH45, GH48, GH28, GH16, GH31, GH27, GH5, and GH1). The numbers of GH45 and GH28 genes identified in D. v. virgifera are among the largest in insects where these genes have been identified. Three GH family genes (GH45, GH48, and GH28) are found almost exclusively in two coleopteran superfamilies (Chrysomeloidea and Curculionoidea) among insects, indicating the possibility of their acquisitions by horizontal gene transfer rather than simple vertical transmission from ancestral lineages of insects. Acquisition of GH genes by horizontal gene transfers and subsequent lineage-specific GH gene expansion appear to have played important roles for phytophagous beetles in specializing on particular groups of host plants and in the case of D. v. virgifera, its close association with maize.