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.

Not drowning, (hand)waving? Molecular phylogenetics, biogeography and evolutionary tempo of the ‘Gondwanan’ midge Stictocladius Edwards (Diptera: Chironomidae)

Many insect clades, especially within the Diptera (true flies), have been considered classically ‘Gondwanan’, with an inference that distributions derive from vicariance of the southern continents. Assessing the role that vicariance has played in the evolution of austral taxa requires testing the location and tempo of diversification and speciation against the well-established predictions of fragmentation of the ancient super-continent. Several early (anecdotal) hypotheses that current austral distributions originate from the breakup of Gondwana derive from studies of taxa within the family Chironomidae (non-biting midges). With the advent of molecular phylogenetics and biogeographic analytical software, these studies have been revisited and expanded to test such conclusions better. Here we studied the midge genus Stictocladius Edwards, from the subfamily Orthocladiinae, which contains austral-distributed clades that match vicariance-based expectations. We resolve several issues of systematic relationships among morphological species and reveal cryptic diversity within many taxa. Time-calibrated phylogenetic relationships among taxa accorded partially with the predicted tempo from geology. For these apparently vagile insects, vicariance-dated patterns persist for South America and Australia. However, as often found, divergence time estimates for New Zealand at c. 50 mya post-date separation of Zealandia from Antarctica and the remainder of Gondwana, but predate the proposed Oligocene ‘drowning’ of these islands. We detail other such ‘anomalous’ dates and suggest a single common explanation rather than stochastic processes. This could involve synchronous establishment following recovery from ‘drowning’ and/or deleteriously warming associated with the mid-Eocene climatic optimum (hence ‘waving’, which refers to cycles of drowning events) plus new availability of topography providing of cool running waters, or all these factors in combination. Alternatively a vicariance explanation remains available, given the uncertain duration of connectivity of Zealandia to Australia–Antarctic–South America via the Lord Howe and Norfolk ridges into the Eocene.

Comparison of intraspecific genetic structure among related chironomids (Diptera) from New Zealand and Patagonia : disparity between potential and realized dispersal

Population genetic studies of freshwater invertebrate taxa in New Zealand and South America are currently few despite the geologically and climatically dynamic histories of these regions. The focus of our study was a comparison of the influence on realized dispersal of 2 closely related nonbiting midges (Chironomidae) of population fragmentation on these separated austral land masses. We used a 734-base pair (bp) fragment of cytochrome c oxidase subunit I (COI) to investigate intraspecific genetic structure in Naonella forsythi Boothroyd in New Zealand and Ferringtonia patagonica Edwards in Patagonia. We proposed hypotheses about their potential dispersal and, hence, expected patterns of genetic structure in these 2 species based on published patterns for the closely related Australian taxon Echinocladius martini Cranston. Genetic structure revealed for both N. forsythi and F. patagonica was characterized by several highly divergent (2.0–10.5%) lineages of late Miocene–Pliocene age within each taxon that were not geographically localized. Many were distributed widely. This pattern differed greatly from population structure in E. martini, which was typified by much greater endemicity of divergent genetic lineages. Nevertheless, diversification of lineages in all 3 taxa appeared to be temporally congruent with the onset of late Miocene glaciations in the southern hemisphere that may have driven fragmentation of suitable habitat, promoting isolation of populations and divergence in allopatry. We argue that differences in realized dispersal post-isolation may be the result of differing availability of suitable habitat in interglacial periods.

A continent-wide study reveals clear relationships between regional abiotic conditions and post-dispersal seed predation

Aim Large-scale patterns linking energy availability, biological productivity and diversity form a central focus of ecology. Despite evidence that the activity and abundance of animals may be limited by climatic variables associated with regional biological productivity (e.g. mean annual precipitation and annual actual evapotranspiration), it is unclear whether plant–granivore interactions are themselves influenced by these climatic factors across broad spatial extents. We evaluated whether climatic conditions that are known to alter the abundance and activity of granivorous animals also affect rates of seed removal. Location Eleven sites across temperate North America. Methods We used a common protocol to assess the removal of the same seed species (Avena sativa) over a 2-day period. Model selection via the Akaike information criterion was used to determine a set of candidate binomial generalized linear mixed models that evaluated the relationship between local climatic data and post-dispersal seed predation. Results Annual actual evapotranspiration was the single best predictor of the proportion of seeds removed. Annual actual evapotranspiration and mean annual precipitation were both positively related to mean seed removal and were included in four and three of the top five models, respectively. Annual temperature range was also positively related to seed removal and was an explanatory variable in three of the top four models. Main conclusions Our work provides the first evidence that energy and precipitation, which are known to affect consumer abundance and activity, also translate to strong, predictable patterns of seed predation across a continent. More generally, these findings suggest that future changes in temperature and precipitation could have widespread consequences for plant species composition in grasslands, through impacts on plant recruitment.

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.

Temporal and spatial patterns of seed dispersal of Musa acuminata by Cynopterus sphinx

The foraging behavior of greater short-nosed fruit bats (Cynopterus sphinx) on wild banana (Musa acuminata) and subsequent dispersal of seeds were studied in the Tropical Rainforest Conservation Area, Xishuangbanna Tropical Botanical Garden, Yunnan province, by direct observation of marked fruits, mist netting, and seed collection. The mean number (± SE) of individual C. sphinx captured by mist net were 2.2 ± 0.33/day and 1.4 ± 0.32/day in the rainy season (September to October) and dry season (November to December), respectively; the difference was not significant. The number of seed pellets expelled was 9.0 ± 1.12/day and 7.2 ± 1.37/day in the rainy and dry seasons respectively; again the difference was not significant. The removal curves for marked fruit were similar for 10 focal trees. Fruits were consumed heavily within two weeks after ripening and all the marked fruit were removed within one month. The difference in seed dispersal was significant between different feeding roosts indicating that patterns of seed dispersal may not be uniform. We found the seeds of M. acuminata can be dispersed by C. sphinx to a distance of about 200 m, and C. sphinx can be considered as an effective seed disperser of M. acuminata.

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.

Trapping to monitor tephritid movement: Results, best practice, and assessment of alternatives

Movement of tephritid flies underpins their survival, reproduction, and ability to establish in new areas and is thus of importance when designing effective management strategies. Much of the knowledge currently available on tephritid movement throughout landscapes comes from the use of direct or indirect methods that rely on the trapping of individuals. Here, we review published experimental designs and methods from mark-release-recapture (MRR) studies, as well as other methods, that have been used to estimate movement of the four major tephritid pest genera (Bactrocera, Ceratitis, Anastrepha, and Rhagoletis). In doing so, we aim to illustrate the theoretical and practical considerations needed to study tephritid movement. MRR studies make use of traps to directly estimate the distance that tephritid species can move within a generation and to evaluate the ecological and physiological factors that influence dispersal patterns. MRR studies, however, require careful planning to ensure that the results obtained are not biased by the methods employed, including marking methods, trap properties, trap spacing, and spatial extent of the trapping array. Despite these obstacles, MRR remains a powerful tool for determining tephritid movement, with data particularly required for understudied species that affect developing countries. To ensure that future MRR studies are successful, we suggest that site selection be carefully considered and sufficient resources be allocated to achieve optimal spacing and placement of traps in line with the stated aims of each study. An alternative to MRR is to make use of indirect methods for determining movement, or more correctly, gene flow, which have become widely available with the development of molecular tools. Key to these methods is the trapping and sequencing of a suitable number of individuals to represent the genetic diversity of the sampled population and investigate population structuring using nuclear genomic markers or non-recombinant mitochondrial DNA markers. Microsatellites are currently the preferred marker for detecting recent population displacement and provide genetic information that may be used in assignment tests for the direct determination of contemporary movement. Neither MRR nor molecular methods, however, are able to monitor fine-scale movements of individual flies. Recent developments in the miniaturization of electronics offer the tantalising possibility to track individual movements of insects using harmonic radar. Computer vision and radio frequency identification tags may also permit the tracking of fine-scale movements by tephritid flies by automated resampling, although these methods come with the same problems as traditional traps used in MRR studies. Although all methods described in this chapter have limitations, a better understanding of tephritid movement far outweighs the drawbacks of the individual methods because of the need for this information to manage tephritid populations.