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.

Effect of broadcast baiting on abundance patterns of red imported fire ants (Hymenoptera: Formicidae) and key local ant genera at long-term monitoring sites in Brisbane, Australia

In 2001, the red imported fire ant (Solenopsis invicta Buren) was identified in Brisbane, Australia. An eradication program involving broadcast bait treatment with two insect growth regulators and a metabolic inhibitor began in September of that year and is currently ongoing. To gauge the impacts of these treatments on local ant populations, we examined long-term monitoring data and quantified abundance patterns of S. invicta and common local ant genera using a linear mixed-effects model. For S. invicta, presence in pitfalls reduced over time to zero on every site. Significantly higher numbers of S. invicta workers were collected on high-density polygyne sites, which took longer to disinfest compared with monogyne and low-density polygyne sites. For local ants, nine genus groups of the 10 most common genera analyzed either increased in abundance or showed no significant trend. Five of these genus groups were significantly less abundant at the start of monitoring on high-density polygyne sites compared with monogyne and low-density polygyne sites. The genus Pheidole significantly reduced in abundance over time, suggesting that it was affected by treatment efforts. These results demonstrate that the treatment regime used at the time successfully removed S. invicta from these sites in Brisbane, and that most local ant genera were not seriously impacted by the treatment. These results have important implications for current and future prophylactic treatment efforts, and suggest that native ants remain in treated areas to provide some biological resistance to S. invicta.

Quality assurance of aerial applications of larvicides for mosquito control : effects of granule and catch tray size on field monitoring programs

Aerial applications of granular insecticides are preferable because they can effectively penetrate vegetation, there is less drift, and no loss of product due to evaporation. We aimed to 1) assess the field efficacy ofVectoBac G to control Aedes vigilax (Skuse) in saltmarsh pools, 2) develop a stochastic-modeling procedure to monitor application quality, and 3) assess the distribution of VectoBac G after an aerial application. Because ground-based studies with Ae. vigilax immatures found that VectoBac G provided effective control below the recommended label rate of 7 kg/ha, we trialed a nominated aerial rate of 5 kg/ha as a case study. Our distribution pattern modeling method indicated that the variability in the number of VectoBac G particles captured in catch-trays was greater than expected for 5 kg/ha and that the widely accepted contour mapping approach to visualize the deposition pattern provided spurious results and therefore was not statistically appropriate. Based on the results of distribution pattern modeling, we calculated the catch tray size required to analyze the distribution of aerially applied granular formulations. The minimum catch tray size for products with large granules was 4 m2 for Altosid pellets and 2 m2 for VectoBac G. In contrast, the minimum catch-tray size for Altosid XRG, Aquabac G, and Altosand, with smaller granule sizes, was 1 m2. Little gain in precision would be made by increasing the catch-tray size further, when the increased workload and infrastructure is considered. Our improved methods for monitoring the distribution pattern of aerially applied granular insecticides can be adapted for use by both public health and agricultural contractors.

Stature estimation of a contemporary Australian sub-population: An evaluation of the Trotter and Gleser method using computed tomography of the femur

Utilising computed tomography scans to allow a virtual analysis of three-dimensional reconstructions of the femur, this project confirms that the traditional 1952 Trotter and Gleser stature estimation equations are inapplicable for a contemporary Queensland population. Therefore, this study introduces modern stature estimation equations for femoral length and fragmentary femoral remains using Bayesian statistics for application in forensic anthropological casework. In addition, it was found that caution needs to be applied when comparing estimated stature to reported stature on the missing persons database due to inaccuracy in Queensland drivers' licences.

Understanding Heliothine (Lepidoptera: Heliothinae) pests: What is a host plant?

Heliothine moths (Lepidoptera: Heliothinae) include some of the world's most devastating pest species. Whereas the majority of nonpest heliothinae specialize on a single plant family, genus, or species, pest species are highly polyphagous, with populations often escalating in size as they move from one crop species to another. Here, we examine the current literature on heliothine host-selection behavior with the aim of providing a knowledge base for research scientists and pest managers. We review the host relations of pest heliothines, with a particular focus on Helicoverpa armigera (Hubner), the most economically damaging of all heliothine species. We then consider the important question of what constitutes a host plant in these moths, and some of the problems that arise when trying to determine host plant status from empirical studies on host use. The top six host plant families in the two main Australian pest species (H. armigera and Helicoverpa punctigera Wallengren) are the same and the top three (Asteraceae, Fabaceae, and Malvaceae) are ranked the same (in terms of the number of host species on which eggs or larvae have been identified), suggesting that these species may use similar cues to identify their hosts. In contrast, for the two key pest heliothines in the Americas, the Fabaceae contains approximate to 1/3 of hosts for both. For Helicoverpa zea (Boddie), the remaining hosts are more evenly distributed, with Solanaceae next, followed by Poaceae, Asteraceae, Malvaceae, and Rosaceae. For Heliothis virescens (F.), the next highest five families are Malvaceae, Asteraceae, Solanaceae, Convolvulaceae, and Scrophulariaceae. Again there is considerable overlap in host use at generic and even species level. H. armigera is the most widely distributed and recorded from 68 plant families worldwide, but only 14 families are recorded as a containing a host in all geographic areas. A few crop hosts are used throughout the range as expected, but in some cases there are anomalies, perhaps because host plant relation studies are not comparable. Studies on the attraction of heliothines to plant odors are examined in the context of our current understanding of insect olfaction, with the aim of better understanding the connection between odor perception and host choice. Finally, we discuss research into sustainable management of pest heliothines using knowledge of heliothine behavior and ecology. A coordinated international research effort is needed to advance our knowledge on host relations in widely distributed polyphagous species instead of the localized, piecemeal approaches to understanding these insects that has been the norm to date.

The infectivity and host range of Orgyia anartoides nucleopolyhedrovirus

The painted apple moth (PAM), Teia anartoides (Walker) (Lepidoptera: Lymantriidae) made a recent incursion into New Zealand. A nucleopolyhedrovirus (NPV), Orgyia anartoides NPV (OranNPV), originally isolated from PAM in Australia, was tested for its pathogenicity to PAM and a range of non-target insect species found in New Zealand, to evaluate its suitability as a microbial control for this insect invader. Dosage-mortality tests showed that OranNPV was highly pathogenic to PAM larvae; mean LT50 values for third instars ranged from 17.9 to 8.1 days for doses from 102 to 105 polyhedral inclusion bodies/larva, respectively. The cause of death in infected insects was confirmed as OranNPV. Molecular analysis established that OranNPV can be identified by PCR and restriction digestion, and this process complemented microscopic examination of infected larvae. No lymantriid species occur in New Zealand; however, the virus had no significant effects on species from five other lepidopteran families (Noctuidae, Tortricidae, Geometridae, Nymphalidae and Plutellidae) or on adult honeybees. Thus, all indications from this initial investigation are that OranNPV would be an important tool in the control of PAM in a future incursion of this species into New Zealand.

Effects of artificial shade on attack by the mahogany shoot borer, Hypsipyla robusta (Moore)

Swietenia macrophylla King (Meliaceae: Swietenioideae) provides one of the premier timbers of the world. The mahogany shoot borer Hypsipyla robusta Moore (Lepidoptera: Pyralidae) is an economically important pest of S. macrophylla throughout Asia, Africa and the Pacific. No viable method of controlling this pest is known. Previous observations have suggested that the presence of overhead shade may reduce attack by H. robusta, but this has not been investigated experimentally. This research was therefore designed to assess the influence of light availability on shoot-borer attack on S. macrophylla, by establishing seedlings under three different artificial shade regimes, then using these seedlings to test oviposition preference of adult moths, neonate larval survival and growth and development of shoot borer larvae. Oviposition preference of shoot borer moths was tested on leaves from seedlings grown under artificial shade for 63 weeks. A significant difference in choice was recorded between treatments, with 27.4 ± 1.5 eggs laid under high shade and 87.1 ± 1.8 under low shade. Neonate larval survival on early flushing leaflets of S. macrophylla did not differ significantly between shade treatments. Larval growth rate, estimated by measuring daily frass width, was significantly higher for those larvae fed on seedlings from the high and medium shade treatments (0.1 mm/day), than the low shade treatment (0.06 mm/day). In laboratory-reared larvae, the total mass of frass produced was significantly higher in the high shade treatment (0.4 g) than under the low shade treatment (0.2 g). Longer tunnel lengths were bored by larvae in plants grown under high shade (12.0 ± 2.4 cm) than under low shade (7.07 ± 1.9 cm). However, pupal mass under low shade was 48% higher than that under the high shade treatment, suggesting that plants grown under high shade were of lower nutritional quality for shoot borer larvae. These results indicate that shading of mahogany seedlings may reduce the incidence of shoot borer attack, by influencing both oviposition and larval development. The establishment of mahogany under suitable shade regimes may therefore provide a basis for controlling shoot borer attack using silvicultural approaches.

Author response to : Brough et al. in response to the recently published article : Lottering,N., MacGregor,D.M.,Barry,M.D.,Reynolds,M.S., Gregory,L.S., 2014. Introducing standardized protocols for anthropological measurement of virtual sub-adult crania using computed tomography 2(1),34–38

Firstly, we would like to thank Ms. Alison Brough and her colleagues for their positive commentary on our published work [1] and their appraisal of our utility of the “off-set plane” protocol for anthropometric analysis. The standardized protocols described in our manuscript have wide applications, ranging from forensic anthropology and paleodemographic research to clinical settings such as paediatric practice and orthopaedic surgical design. We affirm that the use of geometrically based reference tools commonly found in computer aided design (CAD) programs such as Geomagic Design X® are imperative for more automated and precise measurement protocols for quantitative skeletal analysis. Therefore we stand by our recommendation of the use of software such as Amira and Geomagic Design X® in the contexts described in our manuscript...

Temporal Characterization of Ossification of the Crania in Australian Subadults: New standards for Age Estimation using Computed Tomography

After attending this presentation, attendees will gain awareness of the ontogeny of cranial maturation, specifically: (1) the fusion timings of primary ossification centers in the basicranium; and (2) the temporal pattern of closure of the anterior fontanelle, to develop new population-specific age standards for medicolegal death investigation of Australian subadults. This presentation will impact the forensic science community by demonstrating the potential of a contemporary forensic subadult Computed Tomography (CT) database of cranial scans and population data, to recalibrate existing standards for age estimation and quantify growth and development of Australian children. This research welcomes a study design applicable to all countries faced with paucity in skeletal repositories. Accurate assessment of age-at-death of skeletal remains represents a key element in forensic anthropology methodology. In Australian casework, age standards derived from American reference samples are applied in light of scarcity in documented Australian skeletal collections. Currently practitioners rely on antiquated standards, such as the Scheuer and Black1 compilation for age estimation, despite implications of secular trends and population variation. Skeletal maturation standards are population specific and should not be extrapolated from one population to another, while secular changes in skeletal dimensions and accelerated maturation underscore the importance of establishing modern standards to estimate age in modern subadults. Despite CT imaging becoming the gold standard for skeletal analysis in Australia, practitioners caution the application of forensic age standards derived from macroscopic inspection to a CT medium, suggesting a need for revised methodologies. Multi-slice CT scans of subadult crania and cervical vertebrae 1 and 2 were acquired from 350 Australian individuals (males: n=193, females: n=157) aged birth to 12 years. The CT database, projected at 920 individuals upon completion (January 2014), comprises thin-slice DICOM data (resolution: 0.5/0.3mm) of patients scanned since 2010 at major Brisbane Childrens Hospitals. DICOM datasets were subject to manual segmentation, followed by the construction of multi-planar and volume rendering cranial models, for subsequent scoring. The union of primary ossification centers of the occipital bone were scored as open, partially closed or completely closed; while the fontanelles, and vertebrae were scored in accordance with two stages. Transition analysis was applied to elucidate age at transition between union states for each center, and robust age parameters established using Bayesian statistics. In comparison to reported literature, closure of the fontanelles and contiguous sutures in Australian infants occur earlier than reported, with the anterior fontanelle transitioning from open to closed at 16.7±1.1 months. The metopic suture is closed prior to 10 weeks post-partum and completely obliterated by 6 months of age, independent of sex. Utilizing reverse engineering capabilities, an alternate method for infant age estimation based on quantification of fontanelle area and non-linear regression with variance component modeling will be presented. Closure models indicate that the greatest rate of change in anterior fontanelle area occurs prior to 5 months of age. This study complements the work of Scheuer and Black1, providing more specific age intervals for union and temporal maturity of each primary ossification center of the occipital bone. For example, dominant fusion of the sutura intra-occipitalis posterior occurs before 9 months of age, followed by persistence of a hyaline cartilage tongue posterior to the foramen magnum until 2.5 years; with obliteration at 2.9±0.1 years. Recalibrated age parameters for the atlas and axis are presented, with the anterior arch of the atlas appearing at 2.9 months in females and 6.3 months in males; while dentoneural, dentocentral and neurocentral junctions of the axis transitioned from non-union to union at 2.1±0.1 years in females and 3.7±0.1 years in males. These results are an exemplar of significant sexual dimorphism in maturation (p<0.05), with girls exhibiting union earlier than boys, justifying the need for segregated sex standards for age estimation. Studies such as this are imperative for providing updated standards for Australian forensic and pediatric practice and provide an insight into skeletal development of this population. During this presentation, the utility of novel regression models for age estimation of infants will be discussed, with emphasis on three-dimensional modeling capabilities of complex structures such as fontanelles, for the development of new age estimation methods.

Morphometric modelling of ageing in the human pubic symphysis : sexual dimorphism in an Australian population

Despite the prominent use of the pubic symphysis for age estimation in forensic anthropology, little has been documented regarding the quantitative morphological and micro-architectural changes of this surface. Specifically, utilising post-mortem computed tomography data from a large, contemporary Australian adult population, this study aimed to evaluate sexual dimorphism in the morphology and bone composition of the symphyseal surface; and temporal characterisation of the pubic symphysis in individuals of advancing age. The sample consisted of multi-slice computed tomography (MSCT) scans of the pubic symphysis(slice thickness: 0.5 mm, overlap: 0.1 mm) of 200 individuals of Caucasian ancestry aged 15–70 years, obtained in 2011. Surface rendering reconstruction of the symphyseal surface was conducted in OsiriX1 (v.4.1) and quantitative analyses in Rapidform XOSTM and OsteomeasureTM. Morphometric variables including inter-pubic distance, surface area, circumference, maximum height and width of the symphyseal surface and micro-architectural assessment of cortical and trabecular bone compositions were quantified using novel automated engineering software capabilities. The major results of this study are correlated with the macroscopic ossification and degeneration pattern of the symphyseal surface, demonstrating significant age-related changes in the morphometric and bone tissue variables between 15 and 70 years. Regardless of sex, the overall dimensions of the symphyseal surface increased with age, coupled with a decrease in bone mass in the trabecular and cortical bone compartments. Significant differences between the ventral, dorsal and medial cortical surfaces were observed, which may be correlated to bone formation activity dependent on muscle activity and ligamentous attachments. Our study demonstrates significant sexual dimorphism at this site, with males exhibiting greater surface dimensions than females. These baseline results provide a detailed insight into the changes in the structure of the pubic symphysis with ageing and sexually dimorphic features associated with the cortical and trabecular bone profiles.

Human identification from video using advanced gait recognition techniques

The solutions proposed in this thesis contribute to improve gait recognition performance in practical scenarios that further enable the adoption of gait recognition into real world security and forensic applications that require identifying humans at a distance. Pioneering work has been conducted on frontal gait recognition using depth images to allow gait to be integrated with biometric walkthrough portals. The effects of gait challenging conditions including clothing, carrying goods, and viewpoint have been explored. Enhanced approaches are proposed on segmentation, feature extraction, feature optimisation and classification elements, and state-of-the-art recognition performance has been achieved. A frontal depth gait database has been developed and made available to the research community for further investigation. Solutions are explored in 2D and 3D domains using multiple images sources, and both domain-specific and independent modality gait features are proposed.

One and the same: Integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis

The invasive fruit fly Bactrocera invadens Drew, Tsuruta & White, and the Oriental fruit fly Bactrocera dorsalis (Hendel) are highly destructive horticultural pests of global significance. Bactrocera invadens originates from the Indian subcontinent and has recently invaded all of sub-Saharan Africa, while B. dorsalis principally occurs from the Indian subcontinent towards southern China and South-east Asia. High morphological and genetic similarity has cast doubt over whether B. invadens is a distinct species from B. dorsalis. Addressing this issue within an integrative taxonomic framework, we sampled from across the geographic distribution of both taxa and: (i) analysed morphological variation, including those characters considered diagnostic (scutum colour, length of aedeagus, width of postsutural lateral vittae, wing size, and wing shape); (ii) sequenced four loci (ITS1, ITS2, cox1 and nad4) for phylogenetic inference, and; (iii) generated a cox1 haplotype network to examine population structure. Molecular analyses included the closely related species, Bactrocera kandiensis Drew & Hancock. Scutum colour varies from red-brown to fully black for individuals from Africa and the Indian subcontinent. All individuals east of the Indian subcontinent are black except for a few red-brown individuals from China. The postsutural lateral vittae width of B. invadens is narrower than B. dorsalis from eastern Asia, but the variation is clinal, with subcontinent B. dorsalis populations intermediate in size. Aedeagus length, wing shape and wing size cannot discriminate between the two taxa. Phylogenetic analyses failed to resolve B. invadens from B. dorsalis, but did resolve B. kandiensis. Bactrocera dorsalis and B. invadens shared cox1 haplotypes, yet the haplotype network pattern does not reflect current taxonomy or patterns in thoracic colour. Some individuals of B. dorsalis/B. invadens possessed haplotypes more closely related to B. kandiensis than to conspecifics, suggestive of mitochondrial introgression between these species. The combined evidence fails to support the delimitation of B. dorsalis and B. invadens as separate biological species. Consequently, existing biological data for B. dorsalis may be applied to the invasive population in Africa. Our recommendation, in line with other recent publications, is that B. invadens be synonymized with B. dorsalis.