Global proteome changes in larvae of Callosobruchus maculatus Coleoptera:Chrysomelidae:Bruchinae) following ingestion of a cysteine proteinase inhibitor

Tandemly duplicated genes encoding polygalacturonase inhibitors are associated with bruchid (Callosobruchus chinensis) resistance in moth bean (Vigna aconitifolia)

Bruchids are stored-grain insect pests responsible for serious seed loss in legume crops. A previous study using an F₂ population (F2OA) derived from a cross between wild moth-bean (Vigna aconitifolia [Jacq.] Maréchal) accession TN67 (resistant) and cultivated moth-bean accession ICPMO056 (susceptible) revealed that resistance to the azuki bean weevil (Callosobruchus chinensis L.) in TN67 was regulated by a single gene located in the major quantitative trait locus-qVacBrc2.1. In this study, qVacBrc2.1 was finely mapped and candidate genes in this locus were identified using F2OA and another large F₂ population (F2NB) derived from the cross mentioned previously. In contrast to the previous study, segregation analysis in the F2NB population revealed that resistance against this pest was controlled by two genes. Furthermore, the addition of novel markers to qVacBrc2.1 and reanalysis of the QTL in the F2OA population demonstrated that qVacBrc2.1 constituted two linked QTLs-qVacBrc2.1-A and qVacBrc2.1-B. The presence of qVacBrc2.1-B was verified using the population F2NB. Comparative genomics using three Vigna spp. strongly suggested the presence of two tandemly duplicated genes, VacPGIP1 and VacPGIP2, which encoded polygalacturonase inhibitors (polygalacturonase-inhibiting proteins) as the candidates for conferring resistance, but only VacPGIP1 could be successfully cloned and sequenced. The alignment of VacPGIP1 coding sequences of TN67 and ICPMO056 revealed eight single nucleotide polymorphisms, three of which altered the amino-acid sequence of the predicted domains of polygalacturonase inhibitors in ICPMO056. Overall, these findings indicate that VacPGIP1 and VacPGIP2 regulated C. chinensis resistance in TN67.

Thirty Years of Mungbean Genome Research: Where Do We Stand and What Have We Learned?

Mungbean is a socioeconomically important legume crop in Asia that is currently in high demand by consumers and industries both as dried beans and in plant-based protein foods. Marker-assisted and genomics-assisted breeding are promising approaches to efficiently and rapidly develop new cultivars with improved yield, quality, and resistance to biotic and abiotic stresses. Although mungbean was at the forefront of research at the dawn of the plant genomics era 30 years ago, the crop is a "slow runner" in genome research due to limited genomic resources, especially DNA markers. Significant progress in mungbean genome research was achieved only within the last 10 years, notably after the release of the VC1973A draft reference genome constructed using next-generation sequencing technology, which enabled fast and efficient DNA marker development, gene mapping, and identification of candidate genes for complex traits. Resistance to biotic stresses has dominated mungbean genome research to date; however, research is on the rise. In this study, we provide an overview of the past progress and current status of mungbean genomics research. We also discuss and evaluate some research results to provide a better understanding of mungbean genomics.

Recruitment of lysosomal cathepsins B, L and D as digestive enzymes in Coleoptera

The recruitment of the lysosomal cathepsins B (CAB), L (CAL) and D (CAD) as luminal digestive enzymes was investigated in 3 species of beetles. Gene expression was determined by RNA-seq in different regions of the midgut and in the carcasses from the transcriptomes of Dermestes maculatus, Tenebrio molitor and Zabrotes subfasciatus. These data together with phylogenetic analyses, allowed us to identify the sequences of the gene coding for digestive and lysosomal CABs, CADs and CALs in T. molitor and Z. subfasciatus and observe the absence of digestive cathepsins in D. maculatus. Comparisons of structures based on the overall similarity of modelled structures were performed and subsite residues in the lysosomal and digestive CALs were identified by molecular docking. The data showed that S2 subsites are very variable, probably as an adaption to a luminal digestive role. The survey of sequences of the gene coding for cathepsins in the genomes of 13 beetle species from different phylogenetic groups showed that expansion of CAL and CAB genes occurred only in the Cucujiformia clade. Several digestive CABs have a reduced occluding loop, probably to act as digestive enzymes. Pollen-feeding was proposed to be the selective pressure to recruit cathepsins as digestive enzymes in Cucujiformia beetles.

Cloning, expression and characterization of arcelin and its impact on digestive enzymes of the stored product insect pest, Callosobruchus maculatus (F.)

The pulse beetle Callosobruchus maculatus causes potential damage to legume crops by infesting the seeds, leading to a reduction of total protein content. Arcelin found in the wild accessions of the common bean, is an insecticidal protein that has the potency to hamper the metabolism of the bruchid beetle. The arcelin gene from the wild accession of Phaseolus lunatus was isolated and the ORF encoding 158 amino acids was cloned in pET-45b (+) vector. The recombinant clones were transformed in BL21 STAR (DE3) pLysS cells, and the expressed arcelin was purified using Ni-NTA column. The recombinant protein was used in preparing an artificial diet, and the insecticidal activity was elucidated against the bruchid pest C. maculatus. Adult emergence and seed damage were drastically reduced in the treated groups. The response towards ingested diet by digestive enzymes involved in metabolism was elucidated through quantitative gene expression. The highest expression was observed in the aminopeptidase, followed by upregulation of alpha-amylase, glycoside hydrolase family 31 and cathepsin D-like aspartic protease, and downregulation of cathepsin L-like cysteine protease. The recombinant arcelin demonstrates effective insecticidal activity against the bruchid beetle. The changes in digestive enzymes to counteract the anti-nutritional nature of the protein were the strategies of the insect defense mechanism.

The resistance of the cowpea cv. BRS Xiquexique to infestation by cowpea weevil is related to the presence of toxic chitin-binding proteins

The cowpea weevil (Callosobruchus maculatus) is the main pest that attacks cowpea (Vigna unguiculata) seeds during storage, causing nutritional and economic losses in the cowpea crop. Thus, studies aiming to identify resistant cowpea cultivars have been developed. Chitin-binding proteins (CBP), such vicilins and chitinases, have been detected in seeds and related with the toxicity to insects. In this work, we investigated the presence of chitin-binding proteins in the partially resistant cowpea cv. BRS Xiquexique and evaluated their toxicity towards cowpea weevil. The CBP fraction was isolated by chitin affinity chromatography. CBP fraction showed, through 15% SDS PAGE, protein bands with varying molecular masses, mainly below 55 kDa. Proteins present in CBP fraction were identified by Western blotting and mass spectrometry analysis, as vicilins and chitinases. CBP fraction, at 5%, was able to interfere with the development of cowpea weevil, decreasing larval mass and length. A CBV (chitin-binding vicilin) fraction isolated from CBP fraction was toxic, at 2.0%, to C. maculatus, decreasing larval mass and length in 64.3% and 33.23%, respectively. These results suggest that chitin binding proteins, such vicilins and chitinases, may be related to the resistance of cowpea cv. BRS Xiquexique to the infestation by C. maculatus.

hRNAi-mediated knock-down of Sphenophorus levis V-ATPase E in transgenic sugarcane (Saccharum spp interspecific hybrid) affects the insect growth and survival

Transgenic sugarcane expressing V-ATPase subunit E dsRNA affects growth and survival of Sphenophorus levis. Plants being sessile organisms are constantly confronted with several biotic and abiotic stresses. Sugarcane (Saccharum spp) is a major tropical crop widely cultivated for its sugar and other by-products. In Brazil, sugarcane plantations account for significant production losses due to Sphenophorus levis (sugarcane weevil) infestations. With the existing control measures being less effective, there arises a necessity for advanced strategies. Our bioassay injection experiments with V-ATPase E dsRNA in S. levis larvae showed significant mortality and reduction in transcription levels. Furthermore, we down-regulated the V-ATPase E gene of S. levis in transgenic sugarcane using an RNAi approach. The resultant RNAi transgenic lines exhibited reduction in larval growth and survival, without compromising plant performance under controlled environment. Our results illustrate that RNAi-mediated down-regulation of key genes is a promising approach in imparting resistance to sugarcane weevil.

Plants use identical inhibitors to protect their cell wall pectin against microbes and insects

As fundamentally different as phytopathogenic microbes and herbivorous insects are, they enjoy plant-based diets. Hence, they encounter similar challenges to acquire nutrients. Both microbes and beetles possess polygalacturonases (PGs) that hydrolyze the plant cell wall polysaccharide pectin. Countering these threats, plant proteins inhibit PGs of microbes, thereby lowering their infection rate. Whether PG-inhibiting proteins (PGIPs) play a role in defense against herbivorous beetles is unknown. To investigate the significance of PGIPs in insect-plant interactions, feeding assays with the leaf beetle Phaedon cochleariae on Arabidopsis thaliana pgip mutants were performed. Fitness was increased when larvae were fed on mutant plants compared to wild-type plants. Moreover, PG activity was higher, although PG genes were downregulated in larvae fed on PGIP-deficient plants, strongly suggesting that PGIPs impair PG activity. As low PG activity resulted in delayed larval growth, our data provide the first in vivo correlative evidence that PGIPs act as defense against insects.

A second VrPGIP1 allele is associated with bruchid resistance (Callosobruchus spp.) in wild mungbean (Vigna radiata var. sublobata) accession ACC41

Two bruchid species, azuki bean weevil (Callosobruchus chinensis L.) and cowpea weevil (Callosobruchus maculatus F.), are the most important insect pests of mungbean [Vigna radiata (L.) Wilczek] after harvest. Improving bruchid resistance is a major goal for mungbean breeders. Bruchid resistance in mungbean is controlled by a single major locus, Br. The tightly linked VrPGIP1 and VrPGIP2, which encode polygalacturonase-inhibiting proteins (PGIPs), are the candidate genes at the Br locus associated with bruchid resistance. One VrPGIP1 resistance allele and two VrPGIP2 resistance alleles have been identified. In this study, we fine-mapped the bruchid-resistance genes in wild mungbean (V. radiata var. sublobata) accession ACC41 using the F₂ population (574 individuals) derived from the 'Kamphaeng Saen 2' (susceptible) × ACC41 (resistant) cross. A QTL analysis indicated that the resistance to the azuki bean weevil and cowpea weevil in ACC41 is controlled by a major QTL (qBr5.1) and a minor QTL (qBr5.2), which are only 0.3 cM apart. qBr5.1 and qBr5.2 accounted for about 82% and 2% of the resistance variation in the F₂ population, respectively. qBr5.1 was mapped to a 237.35-kb region on mungbean chromosome 5 containing eight annotated genes, including VrPGIP1 and VrPGIP2. An examination of the ACC41 VrPGIP1 and VrPGIP2 sequences revealed a new allele for VrPGIP1 (i.e., VrPGIP1-2). Compared with the wild-type sequence, VrPGIP1-2 has five SNPs, of which four cause amino acid changes (residues 125, 129, 188, and 336). A protein sequence analysis indicated that residues 125 and 129 in VrPGIP1-2 are in a β-sheet B1 region, whereas residues 188 and 336 are in a C₁₀-helix region and at the end of the C-terminal region, respectively. Because the β-sheet B1 region is important for interactions with polygalacturonase (PG), residues 125 and 129 in VrPGIP1-2 likely contribute to bruchid resistance by inhibiting PG. Our results imply that VrPGIP1-2 is associated with the bruchid resistance of wild mungbean accession ACC41. This new resistance allele may be useful for breeding mungbean varieties exhibiting durable bruchid resistance.

The steroid hormone 20-hydroxyecdysone induces phosphorylation and aggregation of stromal interacting molecule 1 for store-operated calcium entry

Oligomerization of stromal interacting molecule 1 (STIM1) promotes store-operated calcium entry (SOCE); however, the mechanism of STIM1 aggregation is unclear. Here, using the lepidopteran insect and agricultural pest cotton bollworm (Helicoverpa armigera) as a model and immunoblotting, RT-qPCR, RNA interference (RNAi), and ChIP assays, we found that the steroid hormone 20-hydroxyecdysone (20E) up-regulates STIM1 expression via G protein-coupled receptors (GPCRs) and the 20E nuclear receptor (EcRB1). We also identified an ecdysone-response element (EcRE) in the 5'-upstream region of the STIM1 gene and also noted that STIM1 is located in the larval midgut during metamorphosis. STIM1 knockdown in larvae delayed pupation time, prevented midgut remodeling, and decreased 20E-induced gene transcription. STIM1 knockdown in a H. armigera epidermal cell line, HaEpi, repressed 20E-induced calcium ion influx and apoptosis. Moreover, 20E-induced STIM1 clustering to puncta and translocation toward the cell membrane. Inhibitors of GPCRs, phospholipase C (PLC), and inositol trisphosphate receptor (IP3R) repressed 20E-induced STIM1 phosphorylation, and we found that two GPCRs are involved in 20E-induced STIM1 phosphorylation. 20E-induced STIM1 phosphorylation on Ser-485 through protein kinase C (PKC), and we observed that Ser-485 phosphorylation is critical for STIM1 clustering, interaction with calcium release-activated calcium channel modulator 1 (Orai1), calcium ion influx, and 20E-induced apoptosis. These results suggest that 20E up-regulates STIM1 phosphorylation for aggregation via GPCRs, followed by interaction with Orai1 to induce SOCE, thereby promoting apoptosis in the midgut during insect metamorphosis.

Pectin Digestion in Herbivorous Beetles: Impact of Pseudoenzymes Exceeds That of Their Active Counterparts

Many protein families harbor pseudoenzymes that have lost the catalytic function of their enzymatically active counterparts. Assigning alternative function and importance to these proteins is challenging. Because the evolution toward pseudoenzymes is driven by gene duplication, they often accumulate in multigene families. Plant cell wall-degrading enzymes (PCWDEs) are prominent examples of expanded gene families. The pectolytic glycoside hydrolase family 28 (GH28) allows herbivorous insects to break down the PCW polysaccharide pectin. GH28 in the Phytophaga clade of beetles contains many active enzymes but also many inactive counterparts. Using functional characterization, gene silencing, global transcriptome analyses, and recordings of life history traits, we found that not only catalytically active but also inactive GH28 proteins are part of the same pectin-digesting pathway. The robustness and plasticity of this pathway and thus its importance for the beetle is supported by extremely high steady-state expression levels and counter-regulatory mechanisms. Unexpectedly, the impact of pseudoenzymes on the pectin-digesting pathway in Phytophaga beetles exceeds even the influence of their active counterparts, such as a lowered efficiency of food-to-energy conversion and a prolongation of the developmental period.

Evolution and functional characterization of CAZymes belonging to subfamily 10 of glycoside hydrolase family 5 (GH5_10) in two species of phytophagous beetles

Hemicelluloses, such as xyloglucan, xylan and mannans, consist of a heterogeneous array of plant-derived polysaccharides that form the plant cell wall. These polysaccharides differ from each other in their structure and physiochemical properties, but they share a β-(1,4)-linked sugar backbone. Hemicelluloses can be hydrolyzed by plant-cell-wall-degrading enzymes (PCWDEs), which are widely distributed in phytopathogenic microbes. Recently, it has become apparent that phytophagous beetles also produce their own PCWDEs. Our previous work identified genes encoding putative mannanases belonging to the subfamily 10 of glycoside hydrolase (GH) family 5 (GH5_10) in the genomes of the leaf beetle, Gastrophysa viridula (Chrysomelidae, Chrysomelinae; one gene), and of the bean beetle, Callosobruchus maculatus (Chrysomelidae, Bruchinae; four genes). In contrast to proteins from other GH5 subfamilies, GH5_10 proteins are patchily distributed within the tree of life and have so far hardly been investigated. We addressed the following questions: Are beetle-derived GH5_10s active PCWDEs? How did they evolve? What is their physiological function? Using heterologous protein expression and enzymatic assays, we show that the G. viridula GH5_10 protein is an endo-β-1,4-mannanase. We also demonstrate that only one out of four C. maculatus GH5_10 proteins is an endo-β-1,4-mannanase, which has additional activity on carboxymethyl cellulose. Unexpectedly, another C. maculatus GH5_10 protein has evolved to use xylan instead of mannans as a substrate. RNAi experiments in G. viridula indicate (i) that the sole GH5_10 protein is responsible for breaking down mannans in the gut and (ii) that this breakdown may rather be accessory and may facilitate access to plant cell content, which is rich in nitrogen and simple sugars. Phylogenetic analyses indicate that coleopteran-derived GH5_10 proteins cluster together with Chelicerata-derived ones. Interestingly, other insect-derived GH5_10 proteins cluster elsewhere, suggesting insects have several independent evolutionary origins.

Novel Alleles of Two Tightly Linked Genes Encoding Polygalacturonase-Inhibiting Proteins (VrPGIP1 and VrPGIP2) Associated with the Br Locus That Confer Bruchid (Callosobruchus spp.) Resistance to Mungbean (Vigna radiata) Accession V2709

Nearly all mungbean cultivars are completely susceptible to seed bruchids (Callosobruchus chinensis and Callosobruchus maculatus). Breeding bruchid-resistant mungbean is a major goal in mungbean breeding programs. Recently, we demonstrated in mungbean (Vigna radiata) accession V2802 that VrPGIP2, which encodes a polygalacturonase inhibiting protein (PGIP), is the Br locus responsible for resistance to C. chinensis and C. maculatus. In this study, mapping in mungbean accession V2709 using a BC₁₁F₂ population of 355 individuals revealed that a single major quantitative trait locus, which controlled resistance to both C. chinensis and C. maculatus, was located in a 237.35 Kb region of mungbean chromosome 5 that contained eight annotated genes, including VrPGIP1 (LOC106760236) and VrPGIP2 (LOC106760237). VrPGIP1 and VrPGIP2 are located next to each other and are only 27.56 Kb apart. Sequencing VrPGIP1 and VrPGIP2 in "V2709" revealed new alleles for both VrPGIP1 and VrPGIP2, named VrPGIP1-1 and VrPGIP2-2, respectively. VrPGIP2-2 has one single nucleotide polymorphism (SNP) at position 554 of wild type VrPGIP2. This SNP is a guanine to cystine substitution and causes a proline to arginine change at residue 185 in the VrPGIP2 of "V2709". VrPGIP1-1 has 43 SNPs compared with wild type and "V2802", and 20 cause amino acid changes in VrPGIP1. One change is threonine to proline at residue 185 in VrPGIP1, which is the same as in VrPGIP2. Sequence alignments of VrPGIP2 and VrPGIP1 from "V2709" with common bean (Phaseolus vulgaris) PGIP2 revealed that residue 185 in VrPGIP2 and VrPGIP1 contributes to the secondary structures of proteins that affect interactions between PGIP and polygalacturonase, and that some amino acid changes in VrPGIP1 also affect interactions between PGIP and polygalacturonase. Thus, tightly linked VrPGIP1 and VrPGIP2 are the likely genes at the Br locus that confer bruchid resistance in mungbean "V2709".

Transgenic sugarcane overexpressing CaneCPI-1 negatively affects the growth and development of the sugarcane weevil Sphenophorus levis

Transgenic sugarcane expressing CaneCPI-1 exhibits resistance to Sphenophorus levis larvae. Transgenic plants have widely been used to improve resistance against insect attack. Sugarcane is an economically important crop; however, great losses are caused by insect attack. Sphenophorus levis is a sugarcane weevil that digs tunnels in the stem base, leading to the destruction of the crop. This insect is controlled inefficiently by chemical insecticides. Transgenic plants expressing peptidase inhibitors represent an important strategy for impairing insect growth and development. Knowledge of the major peptidase group present in the insect gut is critical when choosing the most effective inhibitor. S. levis larvae use cysteine peptidases as their major digestive enzymes, primarily cathepsin L-like activity. In this study, we developed transgenic sugarcane plants that overexpress sugarcane cysteine peptidase inhibitor 1 (CaneCPI-1) and assessed their potential through feeding bioassays with S. levis larvae. Cystatin overexpression in the transgenic plants was evaluated using semi-quantitative RT-PCR, RT-qPCR, and immunoblot assays. A 50% reduction of the average weight was observed in larvae that fed on transgenic plants in comparison to larvae that fed on non-transgenic plants. In addition, transgenic sugarcane exhibited less damage caused by larval attack than the controls. Our results suggest that the overexpression of CaneCPI-1 in sugarcane is a promising strategy for improving resistance against this insect.

Phytocystatins: Defense Proteins against Phytophagous Insects and Acari

This review deals with phytocystatins, focussing on their potential role as defence proteins against phytophagous arthropods. Information about the evolutionary, molecular and biochemical features and inhibitory properties of phytocystatins are presented. Cystatin ability to inhibit heterologous cysteine protease activities is commented on as well as some approaches of tailoring cystatin specificity to enhance their defence function towards pests. A general landscape on the digestive proteases of phytophagous insects and acari and the remarkable plasticity of their digestive physiology after feeding on cystatins are highlighted. Biotechnological approaches to produce recombinant cystatins to be added to artificial diets or to be sprayed as insecticide-acaricide compounds and the of use cystatins as transgenes are discussed. Multiple examples and applications are included to end with some conclusions and future perspectives.

A gene encoding a polygalacturonase-inhibiting protein (PGIP) is a candidate gene for bruchid (Coleoptera: bruchidae) resistance in mungbean (Vigna radiata)

The Br locus confers bruchid resistance in mungbean; VrPGIP2 (encoding a polygalacturonase inhibitor) is a strong candidate gene for this resistance. The VrPGIP2 sequence differs between resistant and susceptible lines. Azuki bean weevil (Callosobruchus chinensis) and cowpea weevil (Callosobruchus maculatus) are serious insect pests of mungbean during storage. Bruchid resistance in mungbean is controlled by a single dominant locus, Br. Although the Br locus has been located on a genetic map, molecular basis and function of the gene remain unknown. In this study, high-resolution mapping using a BC11F2 population of 418 plants derived from a cross between 'Kamphaeng Saen 1' (KPS1; susceptible) and 'V2802' (resistant) using simple sequence repeat (SSR) markers delimited the Br locus to a genomic region of 38 Kb of chromosome 5 containing two annotated genes. EST-SSR marker DMB-SSR158 co-segregated perfectly with the Br locus. Bioinformatics analyses revealed that DMB-SSR158 corresponds to a gene encoding a polygalacturonase inhibitor (polygalacturonase-inhibiting protein PGIP) and was designated as VrPGIP2. Comparison of VrPGIP2 coding sequences between four bruchid-resistant (V2802, V1128, V2817 and TC1966) and four bruchid-susceptible (KPS1, Sulu-1, CM and an unknown accession) mungbean lines revealed six single nucleotide polymorphisms (SNPs) between the resistant and susceptible groups. Three of the six SNPs resulted in amino acid changes; namely, alanine (A) to serine (S) at position 320, leucine (L) to proline (P) at position 332, and threonine (T) to P at position 335 of the VrPGIP2 sequence in resistant lines, compared with that in susceptible lines. The A to S change at position 320 may affect the interaction between PGIP and polygalacuronase. These results indicate that VrPGIP2 is very likely the gene at the Br locus responsible for bruchid resistance in mungbean.

7-Ketocholesterol overcomes drug resistance in chronic myeloid leukemia cell lines beyond MDR1 mechanism

Chronic myeloid leukemia (CML) is a myeloproliferative disease with a characteristic BCR-ABL tyrosine kinase (TK) fusion protein. Despite the clinical efficacy accomplished by TKIs therapies, disease progression may affect patient response rate to these inhibitors due to a multitude of factors that could lead to development of a mechanism known as multidrug resistance (MDR). 7-Ketocholesterol (7KC) is an oxidized cholesterol derivative that has been extensively reported to cause cell death in a variety of cancer models. In this study, we showed the in vitro efficacy of 7KC against MDR leukemia cell line, Lucena. 7KC treatment induced reduction in cell viability, together with apoptosis-mediated cell death. Moreover, downregulation of MDR protein caused intracellular drug accumulation and 7KC co-incubation with either Daunorubicin or Vincristine reduced cell viability compared to the use of each drug alone. Additionally, quantitative label-free mass spectrometry-based protein quantification showed alteration of different molecular pathways involved in cell cycle arrest, induction of apoptosis and misfolded protein response. Conclusively, this study highlights the effect of 7KC as a sensitizing agent of multidrug resistance CML and elucidates its molecular mechanisms.

SIGNIFICANCE

CML patients treated with tyrosine kinase inhibitors (TKIs) have showed a 5-year estimated overall survival of 89%, with cumulative complete cytogenetic response of 87%. However, development of drug resistance is a common feature of the disease progression. This study aimed at showing the effect of 7KC as a cytotoxic and sensitizing agent of multidrug resistance CML cell lines. The cellular and molecular basis of this compound were elucidated using a comprehensive strategy based on quantitative proteomic and cell biology assays. We showed that 7KC induced cell death and overcomes drug resistance in CML through mechanisms that go beyond the classical MDR1 pathways.

Positive selection of digestive Cys proteases in herbivorous Coleoptera

Positive selection is thought to contribute to the functional diversification of insect-inducible protease inhibitors in plants in response to selective pressures exerted by the digestive proteases of their herbivorous enemies. Here we assessed whether a reciprocal evolutionary process takes place on the insect side, and whether ingestion of a positively selected plant inhibitor may translate into a measurable rebalancing of midgut proteases in vivo. Midgut Cys proteases of herbivorous Coleoptera, including the major pest Colorado potato beetle (Leptinotarsa decemlineata), were first compared using a codon-based evolutionary model to look for the occurrence of hypervariable, positively selected amino acid sites among the tested sequences. Hypervariable sites were found, distributed within -or close to- amino acid regions interacting with Cys-type inhibitors of the plant cystatin protein family. A close examination of L. decemlineata sequences indicated a link between their assignment to protease functional families and amino acid identity at positively selected sites. A function-diversifying role for positive selection was further suggested empirically by in vitro protease assays and a shotgun proteomic analysis of L. decemlineata Cys proteases showing a differential rebalancing of protease functional family complements in larvae fed single variants of a model cystatin mutated at positively selected amino acid sites. These data confirm overall the occurrence of hypervariable, positively selected amino acid sites in herbivorous Coleoptera digestive Cys proteases. They also support the idea of an adaptive role for positive selection, useful to generate functionally diverse proteases in insect herbivores ingesting functionally diverse, rapidly evolving dietary cystatins.