Estimating predation rates from molecular gut content analysis.

Several methods have been published to estimate per capita predation rates from molecular gut content analysis relying on intuitive understanding of predation, but none have been formally derived. We provide a theoretical framework for estimating predation rates to identify an accurate method and lay bare its assumptions. Per capita predation can be estimated by multiplying the prey decay rate and the prey quantity in the predators. This assumes that variation in per capita predation rate is approximately normally distributed, prey decay occurs exponentially, and predation is in steady state. We described several ways to estimate steady state predation, including using only qualitative presence-absence data to estimate the decay rate and in addition, we provided a method for estimating per capita predation rate when predation is not in steady state. We used previously published data on aphid consumption by a ladybird beetle in a feeding trial to calculate the predation rate and compare published methods with this theoretically derived method. The estimated predation rate (3.29 ± 0.27 aphids/h) using our derived method was not significantly different from the actual predation rate, 3.11 aphids/h. In contrast, previously published methods were less accurate, underestimating the predation rate (0.33 ± 0.02 to 1.66 ± 0.8 aphids/h) or overestimating it (3.64 ± 0.30 aphids/h). In summary, we provide methods to estimate predation rates even when variation in predation rates is not exactly normally distributed and not in steady state and demonstrate that the prey decay rate, and not the prey detection period, is required.

DNA High-Throughput Sequencing for Arthropod Gut Content Analysis to Evaluate Effectiveness and Safety of Biological Control Agents.

The search for effective biological control agents without harmful non-target effects has been constrained by the use of impractical (field direct observation) or imprecise (cage experiments) methods. While advances in the DNA sequencing methods, more specifically the development of high-throughput sequencing (HTS), have been quickly incorporated in biodiversity surveys, they have been slow to be adopted to determine arthropod prey range, predation rate and food web structure, and critical information to evaluate the effectiveness and safety of a biological control agent candidate. The lack of knowledge on how HTS methods could be applied by ecological entomologists constitutes part of the problem, although the lack of expertise and the high cost of the analysis also are important limiting factors. In this review, we describe how the latest HTS methods of metabarcoding and Lazaro, a method to identify prey by mapping unassembled shotgun reads, can serve biological control research, showing both their power and limitations. We explain how they work to determine prey range and also how their data can be used to estimate predation rates and subsequently be translated into food webs of natural enemy and prey populations helping to elucidate their role in the community. We present a brief history of prey detection through molecular gut content analysis and also the attempts to develop a more precise formula to estimate predation rates, a problem that still remains. We focused on arthropods in agricultural ecosystems, but most of what is covered here can be applied to natural systems and non-arthropod biological control candidates as well.

Metabarcoding versus mapping unassembled shotgun reads for identification of prey consumed by arthropod epigeal predators.

BACKGROUND: A central challenge of DNA gut content analysis is to identify prey in a highly degraded DNA community. In this study, we evaluated prey detection using metabarcoding and a method of mapping unassembled shotgun reads (Lazaro). RESULTS: In a mock prey community, metabarcoding did not detect any prey, probably owing to primer choice and/or preferential predator DNA amplification, while Lazaro detected prey with accuracy 43-71%. Gut content analysis of field-collected arthropod epigeal predators (3 ants, 1 dermapteran, and 1 carabid) from agricultural habitats in Brazil (27 samples, 46-273 individuals per sample) revealed that 64% of the prey species detections by either method were not confirmed by melting curve analysis and 87% of the true prey were detected in common. We hypothesized that Lazaro would detect fewer true- and false-positive and more false-negative prey with greater taxonomic resolution than metabarcoding but found that the methods were similar in sensitivity, specificity, false discovery rate, false omission rate, and accuracy. There was a positive correlation between the relative prey DNA concentration in the samples and the number of prey reads detected by Lazaro, while this was inconsistent for metabarcoding. CONCLUSIONS: Metabarcoding and Lazaro had similar, but partially complementary, detection of prey in arthropod predator guts. However, while Lazaro was almost 2× more expensive, the number of reads was related to the amount of prey DNA, suggesting that Lazaro may provide quantitative prey information while metabarcoding did not.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution.

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Diverse patterns of constitutive and inducible overexpression of detoxifying enzyme genes among resistant Aphis glycines populations.

Understanding the mechanisms of pyrethroid resistance is essential to the effective management of pesticide resistance in Aphis glycines Matsumura (Hemiptera: Aphididae). We mined putative detoxifying enzyme genes in the draft genome sequence of A. glycines for cytochrome oxidase P450 (CYP), glutathione-S-transferase (GST) and esterases (E4 and carboxylesterases-CES). Aphids from clonal populations resistant to pyrethroids from three sites in Minnesota, USA, were screened against a diagnostic LC99 concentration of either λ-cyhalothrin or bifenthrin and detoxifying enzyme genes expression in survivors was analyzed by qPCR. Their expression profiles were compared relative to a susceptible clonal population. We found 61 CYP (40 full-length), seven GST (all full-length), seven E4 (five full-length) and three CES (two full-length) genes, including 24 possible pseudogenes. The detoxifying enzymes had different expression profiles across resistant aphid populations, possibly reflecting differences in the genetic background and pyrethroid selection pressures as the number of constitutively overexpressed detoxifying enzyme genes was correlated with the level of resistance. Our findings will strengthen the understanding of the pyrethroid resistance mechanisms in A. glycines.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest.

BACKGROUND: Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies. RESULTS: Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications. CONCLUSIONS: Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Systemic and sex-biased regulation of OBP expression under semiochemical stimuli.

Constitutive expression of Odorant-Binding Proteins (OBPs) in antennae and other body parts has been examined mainly to infer their involvement in insect olfaction, while their regulation in response to semiochemical stimuli has remained poorly known. Previous studies of semiochemical response were basically done using electrophysiology, which integrates the response of the set of OBPs present in an antenna or sensillum, without revealing the regulation of OBPs or which ones might be involved. In this study we used boll weevil as a model and mined its OBPs by RNA-Seq to study their simultaneous antennal expression by qPCR under controlled semiochemical stimuli with aggregation pheromone and plant volatiles. In the absence of a semiochemical stimulus, 23 of 24 OBPs were constitutively expressed in the antenna in both sexes. Semiochemicals changed systemically the expression of OBPs in both sexes. There were different patterns of up- and down-regulation in female antennae for each semiochemical stimulus, consistent with female chemical ecology. On the other hand, the only response in males was down-regulation of some OBPs. We suggest that these systemic changes in OBP expression might be related to enhancing detection of the semiochemical stimuli and/or priming the olfactory system to detect other environmental chemicals.

Genome Sequence of Cauliflower Mosaic Virus Identified in Earwigs (Doru luteipes) through a Metagenomic Approach.

Here we report the first complete genome sequence of a cauliflower mosaic virus from Brazil, obtained from the gut content of the predator earwig (Doru luteipes). This virus has a genome of 8,030 nucleotides (nt) and shares 97% genome-wide identity with an isolate from Argentina.

Biochemical aspects of the soybean response to herbivory injury by the brown stink bug Euschistus heros (Hemiptera: Pentatomidae).

Plant defense response is an elaborate biochemical process shown to depend on the plant genetic background and on the biological stressor. This work evaluated the soybean biochemical foliar response to brown stink bug herbivory injury through an analysis of redox metabolism and proteomic 2DE profiles of susceptible (BRS Silvania RR) and resistant (IAC-100) varieties. The activity of lipoxygenase-3, guaiacol peroxidase, catalase and ascorbate peroxidase was monitored every 24 h up to 96 h. In the susceptible variety, injury caused an increase in the activities of lipoxygenase 3 and guaiacol peroxidase, no change in ascorbate peroxidase, and a decrease in catalase. In the resistant variety, injury did not cause an alteration of any of these enzymes. The proteomic profiles were evaluated after 24 h of injury and revealed to have a similar proportion (4-5%) of differential protein expression in both varieties. The differential proteins, identified by mass spectrometry, in the susceptible variety were related to general stress responses, to plant defense, and to fungal infections. However, in the resistant variety, the identified change in protein profile was related to Calvin cycle enzymes. While the susceptible variety showed adaptive changes in redox metabolism and expression of stress-responsive proteins, the resistant showed a defense response to circumvent the biological stressor.

Uptake and transfer of a Bt toxin by a Lepidoptera to its eggs and effects on its offspring.

Research on non-target effects of transgenic crop plants has focused primarily on bitrophic, tritrophic and indirect effects of entomotoxins from Bacillus thuringiensis, but little work has considered intergenerational transfer of Cry proteins. This work reports a lepidopteran (Chlosyne lacinia) taking up a Bt entomotoxin when exposed to sublethal or low concentrations, transferring the entomotoxin to eggs, and having adverse effects on the first filial generation (F1) offspring. Two bioassays were conducted using a sublethal concentration of toxin (100.0 ng/µl Cry1Ac) for adults and a concentration equal to the LC10 (2.0 ng/µl Cry1Ac) for larvae. Cry1Ac is the most common entomotoxin expressed in Bt cotton in Brazil. In the adult diet bioassay there was no adverse effect on the parental generation (P0) adults, but the F1 larvae had higher mortality and longer development time compared to F1 larvae of parents that did not ingest Cry1Ac. For the 3rd instar larvae, there was no measurable effect on the P0 larvae, pupae and adults, but the F1 larvae had higher mortality and longer development time. Using chemiluminescent Western Blot, Cry1Ac was detected in F1 eggs laid by P0 butterflies from both bioassays. Our study indicates that, at least for this species and these experimental conditions, a ∼65 kDa insecticidal protein can be taken up and transferred to descendants where it can increase mortality and development time.

Effect of Bt genetic engineering on indirect defense in cotton via a tritrophic interaction.

We present a tritrophic analysis of the potential non-intended pleiotropic effects of cry1Ac gene derived from Bacillus thurigiensis (Bt) insertion in cotton (DeltaPine 404 Bt Bollgard® variety) on the emission of herbivore induced volatile compounds and on the attraction of the egg parasitoid Trichogramma pretisoum (Hymenoptera: Trichogrammatidae). Both the herbivore damaged Bt variety and its non-Bt isoline (DeltaPine DP4049 variety) produced volatiles in higher quantity when compared to undamaged plants and significantly attracted the egg parasitoids (T. pretiosum) when compared to undamaged plants. However, Trichogramma pretiosum did not differentiate between the transgenic and nontransgenic varieties, suggesting that the ratios between the compounds released by herbivory damaged -Bt cotton and herbivory damaged-non Bt cotton did not change significantly. Finally, no detrimental effect of the Bt genetic engineering was detected related to the volatile compounds released by Bollgard cotton on the behavior of the natural enemy studied.

Cloning, expression, and purification of the virulence-associated protein D from Xylella fastidiosa.

In this study, an efficient expression system, based on the pET32Xa/LIC vector, for producing a Xylella fastidiosa virulence-associated protein D, found to have a strong similarity to Riemerella anatipestifer and Actinobacillus actinomycetencomitans VapD protein, is presented. The protein has a molecular mass of 17.637 Da and a calculated pI of 5.49. The selected XFa0052 gene was cloned in the pET32Xa/LIC vector and the plasmid was transformed into Escherichia coli BL21 (DE3) strain at 37 degrees C, with an induction time of 2 h and 1 mM IPTG concentration. The protein present in the soluble fraction was purified by immobilized metal affinity chromatography (IMAC), and had its identity determined by mass spectrometry (MALDI-TOF) and N-terminal sequencing. The purified protein was found as a single band on SDS-PAGE and its correct folding was verified by circular dichroism spectroscopy.

Expression and purification of a putative H-NS nucleoid-associated protein from the phytopathogen Xylella fastidiosa.

The H-NS protein is one of the major constituents of the nucleoid structure that has been implicated in the DNA packaging and in the global regulation of gene expression. The study of this transcriptional regulator is an effort to fight Xylella fastidiosa, a citrus pathogen responsible for a range of economically important plant diseases, including the citrus variegated chlorosis (CVC). The putative H-NS ORF was cloned into a pET32-Xa/LIC vector in order to overexpress it coupled with fusion tags in Escherichia coli BL21(DE3). The expressed recombinant protein was purified by immobilized metal affinity chromatography (Ni-NTA resin) and its identity verified by mass spectrometry (MALDI-TOF). Final purification was performed by cation-exchange chromatography (SP Sepharose Fast Flow) and the purified protein was found as a single band on SDS-PAGE. The folding and its DNA binding activity were verified by circular dichroism and fluorescence spectroscopy, respectively.

The influence of self-fertilization performance and copulation behaviour in reproduction by cross-fertilization in groups of Biomphalaria tenagophila (Mollusca, Planorbidae)

The following hypothesis were tested for groups of simultaneous hermaphrodites Biompharia tenagophila: (a) snails that have reproductive success during the process of self-fertilization do not increase their reproductive success after the end of grouping; (b) the copulation behaviour and the presence of one snail whose eggs have a low viability rate influence the partner's reproductive success by cross-fertilization. Groups were constituted by a homozygous pigmented snail and two albinos: one with a viability rate higher than 70 per cent ("good reproducers") and the other less than 10 per cent ("bad reproducers"). All pigmented snails had viability rates higher than 70 per cent. The "good" and "bad" reproducer albine snails had similar copulation behaviour. However, after the end of grouping, the "bad reproducers" continued to have viability rates less than 10 per cent over 30 days. In 100 per cent of the cases that pigmented snails copulated (performing either a male role or simultaneously male and female roles) exclusively with "good" reproducer albinos, they presented high reproductive success (producing, on average of 8.4 pigmented embryos/eggs-mass). However, in 100 per cent of cases that pigmented snails copulated with both partners, the "good" reproducer albine snails produced none or very few embryos (the highest average was 2.2 pigmented embryos/eggs-mass). Therefore, the production of viable embryos by cross-fertilization was more influenced by self-fertilization performance than by copulation behaviour. The presence of a snail whose eggs have a low viability rate could decrease their partner reproductive success.