A Bayesian hurdle model for analysis of an insect resistance monitoring database

Motivated by the analysis of the Australian Grain Insect Resistance Database (AGIRD), we develop a Bayesian hurdle modelling approach to assess trends in strong resistance of stored grain insects to phosphine over time. The binary response variable from AGIRD indicating presence or absence of strong resistance is characterized by a majority of absence observations and the hurdle model is a two step approach that is useful when analyzing such a binary response dataset. The proposed hurdle model utilizes Bayesian classification trees to firstly identify covariates and covariate levels pertaining to possible presence or absence of strong resistance. Secondly, generalized additive models (GAMs) with spike and slab priors for variable selection are fitted to the subset of the dataset identified from the Bayesian classification tree indicating possibility of presence of strong resistance. From the GAM we assess trends, biosecurity issues and site specific variables influencing the presence of strong resistance using a variable selection approach. The proposed Bayesian hurdle model is compared to its frequentist counterpart, and also to a naive Bayesian approach which fits a GAM to the entire dataset. The Bayesian hurdle model has the benefit of providing a set of good trees for use in the first step and appears to provide enough flexibility to represent the influence of variables on strong resistance compared to the frequentist model, but also captures the subtle changes in the trend that are missed by the frequentist and naive Bayesian models.

How to sequence and annotate insect mitochondrial genomes for systematic and comparative genomics research

Over the past decade the mitochondrial (mt) genome has become the most widely used genomic resource available for systematic entomology. While the availability of other types of ‘–omics’ data – in particular transcriptomes – is increasing rapidly, mt genomes are still vastly cheaper to sequence and are far less demanding of high quality templates. Furthermore, almost all other ‘–omics’ approaches also sequence the mt genome, and so it can form a bridge between legacy and contemporary datasets. Mitochondrial genomes have now been sequenced for all insect orders, and in many instances representatives of each major lineage within orders (suborders, series or superfamilies depending on the group). They have also been applied to systematic questions at all taxonomic scales from resolving interordinal relationships (e.g. Cameron et al., 2009; Wan et al., 2012; Wang et al., 2012), through many intraordinal (e.g. Dowton et al., 2009; Timmermans et al., 2010; Zhao et al. 2013a) and family-level studies (e.g. Nelson et al., 2012; Zhao et al., 2013b) to population/biogeographic studies (e.g. Ma et al., 2012). Methodological issues around the use of mt genomes in insect phylogenetic analyses and the empirical results found to date have recently been reviewed by Cameron (2014), yet the technical aspects of sequencing and annotating mt genomes were not covered. Most papers which generate new mt genome report their methods in a simplified form which can be difficult to replicate without specific knowledge of the field. Published studies utilize a sufficiently wide range of approaches, usually without justification for the one chosen, that confusion about commonly used jargon such as ‘long PCR’ and ‘primer walking’ could be a serious barrier to entry. Furthermore, sequenced mt genomes have been annotated (gene locations defined) to wildly varying standards and improving data quality through consistent annotation procedures will benefit all downstream users of these datasets. The aims of this review are therefore to: 1. Describe in detail the various sequencing methods used on insect mt genomes; 2. Explore the strengths/weakness of different approaches; 3. Outline the procedures and software used for insect mt genome annotation, and; 4. Highlight quality control steps used for new annotations, and to improve the re-annotation of previously sequenced mt genomes used in systematic or comparative research.

Discovery of a lipid synthesising organ in the auditory system of an insect

Weta possess typical Ensifera ears. Each ear comprises three functional parts: two equally sized tympanal membranes, an underlying system of modified tracheal chambers, and the auditory sensory organ, the crista acustica. This organ sits within an enclosed fluid-filled channel-previously presumed to be hemolymph. The role this channel plays in insect hearing is unknown. We discovered that the fluid within the channel is not actually hemolymph, but a medium composed principally of lipid from a new class. Three-dimensional imaging of this lipid channel revealed a previously undescribed tissue structure within the channel, which we refer to as the olivarius organ. Investigations into the function of the olivarius reveal de novo lipid synthesis indicating that it is producing these lipids in situ from acetate. The auditory role of this lipid channel was investigated using Laser Doppler vibrometry of the tympanal membrane, which shows that the displacement of the membrane is significantly increased when the lipid is removed from the auditory system. Neural sensitivity of the system, however, decreased upon removal of the lipid-a surprising result considering that in a typical auditory system both the mechanical and auditory sensitivity are positively correlated. These two results coupled with 3D modelling of the auditory system lead us to hypothesize a model for weta audition, relying strongly on the presence of the lipid channel. This is the first instance of lipids being associated with an auditory system outside of the Odentocete cetaceans, demonstrating convergence for the use of lipids in hearing.

Opsin clines in butterflies suggest novel roles for insect photopigments

Opsins are ancient molecules that enable animal vision by coupling to a vitamin-derived chromophore to form lightsensitive photopigments. The primary drivers of evolutionary diversification in opsins are thought to be visual tasks related to spectral sensitivity and color vision. Typically, only a few opsin amino acid sites affect photopigment spectral sensitivity. We show that opsin genes of the North American butterfly Limenitis arthemis have diversified along a latitudinal cline, consistent with natural selection due to environmental factors. We sequenced single nucleotide(SNP) polymorphisms in the coding regions of the ultraviolet (UVRh), blue (BRh), and long-wavelength (LWRh) opsin genes from ten butterfly populations along the eastern United States and found that a majority of opsin SNPs showed significant clinal variation. Outlier detection and analysis of molecular variance indicated that many SNPs are under balancing selection and show significant population structure. This contrasts with what we found by analysing SNPs in the wingless and EF-1 alpha loci, and from neutral amplified fragment length polymorphisms, which show no evidence of significant locus-specific or genome-wide structure among populations. Using a combination of functional genetic and physiological approaches, including expression in cell culture, transgenic Drosophila, UV-visible spectroscopy, and optophysiology, we show that key BRh opsin SNPs that vary clinally have almost no effect on spectral sensitivity. Our results suggest that opsin diversification in this butterfly is more consistent with natural selection unrelated to spectral tuning. Some of the clinally varying SNPs may instead play a role in regulating opsin gene expression levels or the thermostability of the opsin protein. Lastly, we discuss the possibility that insect opsins might have important, yet-to-be elucidated, adaptive functions in mediating animal responses to abiotic factors, such as temperature or photoperiod.

Expression, Purification and Characterization of Peanut (Arachis hypogaea) Agglutinin (PNA) from Baculovirus Infected Insect Cells

Peanut (Arachis hypogaea) seed lectin, PNA is widely used to identify tumor specific antigen (T-antigen), Gal beta 1-3GalNAc on the eukaryotic cell surface. The functional amino acid coding region of a cDNA clone, pBSH-PN was PCR amplified and cloned downstream of the polyhedrin promoter in the Autographa californica nucleopolyhedrovirus (AcNPV) based transfer vector pVL1393. Co-transfection of Spodoptera frugiperda cells (Sf9) with the transfer vector, pAcPNA and AcRP6 (a recombinant AcNPV having B-gal downstream of the polyhedrin promoter) DNAs produced a recombinant virus, AcPNA which expresses PNA. Infection of suspension culture of Sf9 cells with plaque purified AcPNA produced as much as 9.8 mg PNA per liter (2.0 x 10(6) cells/ml) of serum-free medium. Intracellularly expressed protein (re-PNA) was purified to apparent homogeneity by affinity chromatography using ECD-Sepharose. Polyclonal antibodies against natural PNA (n-PNA) crossreacted with re-PNA. The subunit molecular weight (30 kDa), hemagglutination activity, and carbohydrate specificity of re-PNA were found to be identical to that of n-PNA, thus confirming the abundant production of a functionally active protein in the baculovirus expression system.

Prospects for predicting insect mortality in relation to changing phosphine concentrations

Typically, in bag-stack or silo fumigations the concentration of phosphine is not constant, and yet most of what is known about phosphine efficacy against grain insects comes from studies with fixed concentrations. Indeed, where changing concentration experiments have been performed, researchers have been unable to explain observed efficacy on the basis of data from fixed concentrations. The ability to predict insect mortality in relation to changing phosphine concentrations would facilitate the development of effective fumigation protocols. In this paper, we explore the prospects for making such predictions. After reviewing published and new results, we conclude that the commonly used concentration x time (Ct) product is unreliable for this purpose. New results, for a strongly resistant strain of Rhyzopertha dominica from Australia, suggest that the relationship Cnt = k may be useful for predicting mortality of this type of insect in changing concentrations. However, in the case of a strain of Sitophilus oryzae with a type of resistance common in Australian S. oryzae, the relationship Cnt = k proved to be less reliable.

A precision apparatus, with solid phase micro-extraction monitoring capability, for incorporation studies of gaseous precursors into insect-derived metabolites

An apparatus is described that facilitates the determination of incorporation levels of isotope labelled, gaseous precursors into volatile insect-derived metabolites. Atmospheres of varying gas compositions can be generated by evacuation of a working chamber followed by admission of the required levels of component gases, using a precision, digitised pressure read-out system. Insects such as fruit-flies are located initially in a small introduction chamber, from which migration can occur downwards into the working chamber. The level of incorporation of labelled precursors is continuously assayed by the Solid Phase Micro Extraction (SPME) technique and GC-MS analyses. Experiments with both Bactrocera species (fruit-flies) and a parasitoid wasp, Megarhyssa nortoni nortoni (Cresson) and oxygen-18 labelled dioxygen illustrate the utility of this system. The isotope effects of oxygen-18 on the carbon-13 NMR spectra of 1,7- dioxaspiro[5,5]undecane are also described.

Resistance to insect growth regulator insecticides in populations of sheep lice as assessed by a moulting disruption assay

Low-volume, backline applications with the benzoylphenyl urea insecticides triflumuron and diflubenzuron represent in excess of 70% of treatments for the control of sheep lice, Bovicola ovis (Schrank) (Phthiraptera: Trichodectidae), in Australia. Reports of reduced effectiveness from 2003 and subsequent controlled treatment trials suggested the emergence of resistance to these compounds in B. ovis populations. A laboratory assay based on the measurement of moulting success in nymphs was developed and used to assess susceptibility to diflubenzuron and triflumuron in louse populations collected from sheep where a control failure had occurred. These tests confirmed the development of resistance to triflumuron and diflubenzuron in at least two instances, with estimated resistance ratios of 67-94X at LC50.

The host range of Isturgia deeraria, an insect considered for the biological control of Acacia nilotica in Australia

The geometrid caterpillar Isturgia deerraria was imported from Kenya into quarantine facilities in Australia as a potential biological control agent for prickly acacia, Acacia nilotica subsp. indica (Benth.) Brenan (family Mimosaceae). The insect was then tested on 30 plant species presented to neonate larvae as a no-choice cut foliage test and 13 plant species presented as a no-choice potted plant test. In these tests the insect was able to complete its life cycle on 13 native Acacia spp. and also on Acacia farnesiana and the exotic ornamental Delonix regia (family Caesalpiniaceae). The tests supported field observations that the insect has a host range spanning many leguminous species and as such the insect could not be considered for release in Australia.

Interrogating an insect society

Insect societies such as those of ants, bees, and wasps consist of 1 or a small number of fertile queens and a large number of sterile or nearly sterile workers. While the queens engage in laying eggs, workers perform all other tasks such as nest building, acquisition and processing of food, and brood care. How do such societies function in a coordinated and efficient manner? What are the rules that individuals follow? How are these rules made and enforced? These questions are of obvious interest to us as fellow social animals but how do we interrogate an insect society and seek answers to these questions? In this article I will describe my research that was designed to see answers from an insect society to a series of questions of obvious interest to us. I have chosen the Indian paper wasp Ropalidia marginata for this purpose, a species that is abundantly distributed in peninsular India and serves as an excellent model system. An important feature of this species is that queens and workers are morphologically identical and physiologically nearly so. How then does an individual become a queen? How does the queen suppress worker reproduction? How does the queen regulate the nonreproductive activities of the workers? What is the function of aggression shown by different individuals? How and when is the queen's heir decided? I will show how such questions can indeed be investigated and will emphasize the need for a whole range of different techniques of observation and experimentation.

Characterising insect plant host relationships facilitates understanding multiple host use

Weed biocontrol relies on host specificity testing, usually carried out under quarantine conditions to predict the future host range of candidate control agents. The predictive power of host testing can be scrutinised directly with Aconophora compressa, previously released against the weed Lantana camara L. (lantana) because its ecology in its new range (Australia) is known and includes the unanticipated use of several host species. Glasshouse based predictions of field host use from experiments designed a posteriori can therefore be compared against known field host use. Adult survival, reproductive output and egg maturation were quantified. Adult survival did not differ statistically across the four verbenaceous hosts used in Australia. Oviposition was significantly highest on fiddlewood (Citharexylum spinosum L.), followed by lantana, on which oviposition was significantly higher than on two varieties of Duranta erecta (‘‘geisha girl’’ and ‘‘Sheena’s gold’’; all Verbenaceae). Oviposition rates across Duranta varieties were not significantly different from each other but were significantly higher than on the two non-verbenaceous hosts (Jacaranda mimosifolia D. Don: Bignoneaceae (jacaranda) and Myoporum acuminatum R. Br.: Myoporaceae (Myoporum)). Production of adult A. compressa was modelled across the hosts tested. The only major discrepancy between model output and their relative abundance across hosts in the field was that densities on lantana in the field were much lower than predicted by the model. The adults may, therefore, not locate lantana under field conditions and/or adults may find lantana but leave after laying relatively few eggs. Fiddlewood is the only primary host plant of A. compressa in Australia, whereas lantana and the others are used secondarily or incidentally. The distinction between primary, secondary and incidental hosts of a herbivore species helps to predict the intensity and regularity of host use by that herbivore. Populations of the primary host plants of a released biological control agent are most likely to be consistently impacted by the herbivore, whereas secondary and incidental host plant species are unlikely to be impacted consistently. As a consequence, potential biocontrol agents should be released only against hosts to which they have been shown to be primarily adapted.