Search

Females of many insect species are unreceptive to remating for a period following their first mating. This inhibitory effect may be mediated by either the female or her first mate, or both, and often reflects the complex interplay of reproductive strategies between the sexes. Natural variation in remating inhibition and how this phenotype responds to captive breeding are largely unexplored in insects, including many pest species. We investigated genetic variation in remating propensity in the Queensland fruit fly, Bactrocera tryoni, using strains differing in source locality and degree of domestication. We found up to threefold inherited variation between strains from different localities in the level of intra-strain remating inhibition. The level of inhibition also declined significantly during domestication, which implied the existence of genetic variation for this trait within the starting populations as well. Inter-strain mating and remating trials showed that the strain differences were mainly due to the genotypes of the female and, to a lesser extent, the second male, with little effect of the initial male genotype. Implications for our understanding of fruit fly reproductive biology and population genetics and the design of Sterile Insect Technique pest management programs are discussed.

Female insects commonly have more than one mate during a breeding period ('polyandry'), storing and using sperm from multiple males. In addition to its evolutionary significance, insect polyandry has practical implications for pest management that relies on the sterile insect technique (SIT). The Queensland fruit fly, Bactrocera tryoni (Froggatt), is a major horticultural pest in Australia, and outbreaks are managed by SIT in some regions. The present study provides the first evidence for polyandry in female B. tryoni from field populations from New South Wales (NSW) and Queensland (QLD) through multi-locus genotyping (ten microsatellite markers in four fluorescent multiplexes) of the stored sperm in ovipositing females. Polyandry level was significantly higher in the NSW collection (80.0 %) than the QLD collection (26.1 %), suggesting substantial regional and/or temporal variation. These findings have important implications for the use of SIT to suppress B. tryoni populations and to eradicate outbreaks.

Rectal gland volatiles are key mediators of sexual interactions in tephritid fruit flies. We used solid-phase microextraction (SPME) plus gas chromatography-mass spectrometry (GC-MS) and gas chromatography-flame ionization detection (GC-FID) to substantially expand rectal gland chemical characterisation of the Queensland fruit fly (Bactrocera tryoni (Diptera: Tephritidae); Qfly). The SPME GC-MS analysis identified 24 of the 30 compounds previously recorded from Qfly rectal glands, plus another 21 compounds that had not previously been reported. A few amides and fatty acid esters dominated the chromatograms of males and females respectively, but we also found other esters, alcohols and aldehydes and a ketone. The GC-FID analyses also revealed over 150 others, as yet unidentified, volatiles, generally in lesser amounts. The GC-FID analyses also showed 49 and 12 compounds were male- and female-specific, respectively, both in single sex (virgin) and mixed sex (mostly mated) groups. Another ten compounds were male-specific among virgins but undetected in mixed sex groups, and 29 were undetected in virgins but male-specific in mixed sex groups. The corresponding figures for females were four and zero, respectively. Most short retention time peaks (including a ketone and an ester) were male-specific, whereas most female-biased peaks (including five fatty acid esters) had long retention times. Our results indicate previously unsuspected diversity of rectal gland volatiles that might have pheromone functions in males, but far fewer in females.

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

The Queensland fruit fly, Bactrocera tryoni, is a major pest of Australian horticulture which has expanded its range in association with the spread of horticulture over the last ~ 150 years. Its distribution in northern Australia overlaps that of another fruit fly pest to which some authors accord full species status, Bactrocera aquilonis. We have used reduced representation genome-wide sequencing to genotype 359 individuals taken from 35 populations from across the current range of the two taxa, plus a further 73 individuals from six of those populations collected 15-22 years earlier. We find significant population differentiation along an east-west transect across northern Australia which likely reflects limited but bidirectional gene flow between the two taxa. The southward expansion of B. tryoni has led to relatively little genetic differentiation, and most of it is associated with a move into previously marginal inland habitats. Two disjunct populations elsewhere in Australia and three on Melanesian islands are each clearly differentiated from all others, with data strongly supporting establishment from relatively few founders and significant isolation subsequently. Resequencing of historical samples from one of the disjunct Australian populations shows that its genetic profile has changed little over a 15-year period, while the Melanesian data suggest a succession of 'island hopping' events with progressive reductions in genetic diversity. We discuss our results in relation to the control of B. tryoni and as a model for understanding the genetics of invasion and hybridisation processes.

BACKGROUND: The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a "common garden" approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication. RESULTS: Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait. CONCLUSIONS: Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

BACKGROUND: Bactrocera tryoni and Bactrocera neohumeralis mate asynchronously; the former mates exclusively around dusk while the latter mates during the day. The two species also differ in the colour of the post-pronotal lobe (callus), which is predominantly yellow in B. tryoni and brown in B. neohumeralis. We have examined the genetic relationship between the two characters in hybrids, backcrosses and multigeneration hybrid progeny. RESULTS: Our analysis of the mating time of the parental species revealed that while B. tryoni mate exclusively at dusk, B. neohumeralis females pair with B. neohumeralis males during the day and with B. tryoni males at dusk. We found considerable variance in mating time and callus colour among hybrid backcross individuals of both sexes but there was a strong although not invariant trend for callus colour to co-segregate with mating time in both sexes. To genetically separate these two phenotypes we allowed the interspecific F1 hybrids to propagate for 25 generations (F25) without selection for mating time or callus colour, finding that the advanced hybrid population had moved towards B. tryoni phenotypes for both traits. Selection for day mating in replicate lines at F25 resulted in significant phenotypic shifts in both traits towards B. neohumeralis phenotypes in F26. However, we were unable to completely recover the mating time profile of B. neohumeralis and relaxation of selection for day mating led to a shift back towards dusk mating, but not yellow callus colour, by F35. CONCLUSION: We conclude that the inheritance of the two major species-defining traits is separable but tightly linked and involves more than one gene in each case. It also appears that laboratory conditions select for the B. tryoni phenotypes for mating time. We discuss our findings in relation to speciation theory and the likely effects of domestication during the generation of mass release strains for sterile insect control programmes.

BACKGROUND: Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species. RESULTS: We find disproportionately high rates of gene gains in internal branches in the species' phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species' differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one. CONCLUSION: Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change.

BACKGROUND: Distinguishing orthologous and paralogous relationships between genes across multiple species is essential for comparative genomic analyses. Various computational approaches have been developed to resolve these evolutionary relationships, but strong trade-offs between precision and recall of orthologue prediction remains an ongoing challenge. RESULTS: Here we present Orthonome, an orthologue prediction pipeline, designed to reduce the trade-off between orthologue capture rates (recall) and accuracy of multi-species orthologue prediction. The pipeline compares sequence domains and then forms sequence-similar clusters before using phylogenetic comparisons to identify inparalogues. It then corrects sequence similarity metrics for fragment and gene length bias using a novel scoring metric capturing relationships between full length as well as fragmented genes. The remaining genes are then brought together for the identification of orthologues within a phylogenetic framework. The orthologue predictions are further calibrated along with inparalogues and gene births, using synteny, to identify novel orthologous relationships. We use 12 high quality Drosophila genomes to show that, compared to other orthologue prediction pipelines, Orthonome provides orthogroups with minimal error but high recall. Furthermore, Orthonome is resilient to suboptimal assembly/annotation quality, with the inclusion of draft genomes from eight additional Drosophila species still providing >6500 1:1 orthologues across all twenty species while retaining a better combination of accuracy and recall than other pipelines. Orthonome is implemented as a searchable database and query tool along with multiple-sequence alignment browsers for all sets of orthologues. The underlying documentation and database are accessible at http://www.orthonome.com . CONCLUSION: We demonstrate that Orthonome provides a superior combination of orthologue capture rates and accuracy on complete and draft drosophil

Wolbachia infections have been described in several Drosophila species, but relatively few have been assessed for phenotypic effects. Cytoplasmic incompatibility (CI) is the most common phenotypic effect that has been detected, while some infections cause male killing or feminization, and many Wolbachia infections have few host effects. Here, we describe two new infections in a recently described species, Drosophila pandora, one of which causes near-complete CI and near-perfect maternal transmission (the "CI" strain). The other infection is a male killer (the "MK" strain), which we confirm by observing reinitiation of male production following tetracycline treatment. No incompatibility was detected in crosses between CI strain males and MK strain females, and rare MK males do not cause CI. Molecular analyses indicate that the CI and MK infections are distantly related and the CI infection is closely related to the wRi infection of Drosophila simulans. Two population surveys indicate that all individuals are infected with Wolbachia, but the MK infection is uncommon. Given patterns of incompatibility among the strains, the infection dynamics is expected to be governed by the relative fitness of the females, suggesting that the CI infection should have a higher fitness. This was evidenced by changes in infection frequencies and sex ratios in population cages initiated at different starting frequencies of the infections.

Wolbachia is the most widespread endosymbiotic bacterium of insects and other arthropods that can rapidly invade host populations. Deliberate releases of Wolbachia into natural populations of the dengue fever mosquito, Aedes aegypti, are used as a novel biocontrol strategy for dengue suppression. Invasion of Wolbachia through the host population relies on factors such as high fidelity of the endosymbiont transmission and limited immigration of uninfected individuals, but these factors can be difficult to measure. One way of acquiring relevant information is to consider mitochondrial DNA (mtDNA) variation alongside Wolbachia in field-caught mosquitoes. Here we used diagnostic mtDNA markers to differentiate infection-associated mtDNA haplotypes from those of the uninfected mosquitoes at release sites. Unique haplotypes associated with Wolbachia were found at locations outside Australia. We also performed mathematical and qualitative analyses including modelling the expected dynamics of the Wolbachia and mtDNA variants during and after a release. Our analyses identified key features in haplotype frequency patterns to infer the presence of imperfect maternal transmission of Wolbachia, presence of immigration and possibly incomplete cytoplasmic incompatibility. We demonstrate that ongoing screening of the mtDNA variants should provide information on maternal leakage and immigration, particularly in releases outside Australia. As we demonstrate in a case study, our models to track the Wolbachia dynamics can be successfully applied to temporal studies in natural populations or Wolbachia release programs, as long as there is co-occurring mtDNA variation that differentiates infected and uninfected populations.

Lucilia cuprina is a parasitic fly of major economic importance worldwide. Larvae of this fly invade their animal host, feed on tissues and excretions and progressively cause severe skin disease (myiasis). Here we report the sequence and annotation of the 458-megabase draft genome of Lucilia cuprina. Analyses of this genome and the 14,544 predicted protein-encoding genes provide unique insights into the fly's molecular biology, interactions with the host animal and insecticide resistance. These insights have broad implications for designing new methods for the prevention and control of myiasis.

Mutations in the voltage-sensitive sodium channel gene (Vssc) have been identified in Aedes aegypti and some have been associated with pyrethroid insecticide resistance. Whether these mutations cause resistance, alone or in combination with other alleles, remains unclear, but must be understood if mutations are to become markers for resistance monitoring. We describe High Resolution Melt (HRM) genotyping assays for assessing mutations found in Ae. aegypti in Indonesia (F1565C, V1023G, S996P) and use them to test for associations with pyrethroid resistance in mosquitoes from Yogyakarta, a city where insecticide use is widespread. Such knowledge is important because Yogyakarta is a target area for releases of Wolbachia-infected mosquitoes with virus-blocking traits for dengue suppression. We identify three alleles across Yogyakarta putatively linked to resistance in previous research. By comparing resistant and susceptible mosquitoes from bioassays, we show that the 1023G allele is associated with resistance to type I and type II pyrethroids. In contrast, F1565C homozygotes were rare and there was only a weak association between individuals heterozygous for the mutation and resistance to a type I pyrethroid. As the heterozygote is expected to be incompletely recessive, it is likely that this association was due to a different resistance mechanism being present. A resistance advantage conferred to V1023G homozygotes through addition of the S996P allele in the homozygous form was suggested for the Type II pyrethroid, deltamethrin. Screening of V1023G and S996P should assist resistance monitoring in Ae. aegypti from Yogyakarta, and these mutations should be maintained in Wolbachia strains destined for release in this city to ensure that these virus-blocking strains of mosquitoes are not disadvantaged, relative to resident populations.

The maternally inherited intracellular bacteria Wolbachia can manipulate host reproduction in various ways that foster frequency increases within and among host populations. Manipulations involving cytoplasmic incompatibility (CI), where matings between infected males and uninfected females produce non-viable embryos, are common in arthropods and produce a reproductive advantage for infected females. CI was associated with the spread of Wolbachia variant wRi in Californian populations of Drosophila simulans, which was interpreted as a bistable wave, in which local infection frequencies tend to increase only once the infection becomes sufficiently common to offset imperfect maternal transmission and infection costs. However, maternally inherited Wolbachia are expected to evolve towards mutualism, and they are known to increase host fitness by protecting against infectious microbes or increasing fecundity. We describe the sequential spread over approximately 20 years in natural populations of D. simulans on the east coast of Australia of two Wolbachia variants (wAu and wRi), only one of which causes significant CI, with wRi displacing wAu since 2004. Wolbachia and mtDNA frequency data and analyses suggest that these dynamics, as well as the earlier spread in California, are best understood as Fisherian waves of favourable variants, in which local spread tends to occur from arbitrarily low frequencies. We discuss implications for Wolbachia-host dynamics and coevolution and for applications of Wolbachia to disease control.

The Mitochondrial Carrier Family (MCF) is a signature group of integral membrane proteins that transport metabolites across the mitochondrial inner membrane in eukaryotes. MCF proteins are characterized by six transmembrane segments that assemble to form a highly-selective channel for metabolite transport. We discovered a novel MCF member, termed Legionellanucleotide carrier Protein (LncP), encoded in the genome of Legionella pneumophila, the causative agent of Legionnaire's disease. LncP was secreted via the bacterial Dot/Icm type IV secretion system into macrophages and assembled in the mitochondrial inner membrane. In a yeast cellular system, LncP induced a dominant-negative phenotype that was rescued by deleting an endogenous ATP carrier. Substrate transport studies on purified LncP reconstituted in liposomes revealed that it catalyzes unidirectional transport and exchange of ATP transport across membranes, thereby supporting a role for LncP as an ATP transporter. A hidden Markov model revealed further MCF proteins in the intracellular pathogens, Legionella longbeachae and Neorickettsia sennetsu, thereby challenging the notion that MCF proteins exist exclusively in eukaryotic organisms.

Clinal variation in quantitative traits is widespread, but its genetic basis awaits identification. Drosophila melanogaster shows adaptive, clinal variation in traits such as body size along latitudinal gradients on multiple continents. To investigate genome wide transcription differentiation between North and South that might contribute to the clinal phenotypic variation, we compared RNA expression patterns during development of D. melanogaster from tropical northern and temperate southern populations using whole genome tiling arrays. We found that genes that were differentially expressed between the cline ends were generally associated with metabolism and growth, and experimental alteration of expression of a sample of them generally resulted in altered body size in the predicted direction, sometimes significantly so. We further identified the serpent (srp) transcription factor binding sites to be enriched near genes up-regulated in expression in the south. Analysis of clinal populations revealed a significant cline in the expression level of srp. Experimental over-expression of srp increased body size, as predicted from its clinal expression pattern, suggesting that it may be involved in regulating adaptive clinal variation in Drosophila. This study identified a handful of genes that contributed to clinal phenotypic variation through altered gene expression level, yet misexpression of individual gene led to modest body size change.

BACKGROUND: The Australian sheep blowfly Lucilia cuprina (Wiedemann) (Diptera: Calliphoridae) is a destructive pest of the sheep, a model organism for insecticide resistance research, and a valuable tool for medical and forensic professionals. However, genomic information on L. cuprina is still sparse. RESULTS: We report here the construction of an embryonic and 2 larval cDNA libraries for L. cuprina. A total of 29,816 expressed sequence tags (ESTs) were obtained and assembled into 7,464 unique clusters. The sequence collection captures a great diversity of genes, including those related to insecticide resistance (e.g., 12 cytochrome P450s, 2 glutathione S transferases, and 6 esterases). Compared to Drosophila melanogaster, codon preference is different in 13 of the 18 amino acids encoded by redundant codons, reflecting the lower overall GC content in L. cuprina. In addition, we demonstrated that the ESTs could be converted into informative gene markers by capitalizing on the known gene structures in the model organism D. melanogaster. We successfully assigned 41 genes to their respective chromosomes in L. cuprina. The relative locations of these loci revealed high but incomplete chromosomal synteny between L. cuprina and D. melanogaster. CONCLUSIONS: Our results represent the first major transcriptomic undertaking in L. cuprina. These new genetic resources could be useful for the blowfly and insect research community.

BACKGROUND: Progression of development has to be insulated from the damaging impacts of environmental and genetic perturbations to produce highly predictable phenotypes. Molecular chaperones, such as the heat shock proteins (HSPs), are known to buffer various environmental stresses, and are deeply involved in protein homeostasis. These characteristics of HSPs imply that they might affect developmental buffering and canalization. RESULTS: We examined the role of nine Hsp genes using the GAL4/UAS-RNAi system on phenotypic variation of various morphological traits in Drosophila melanogaster. The stability of bristle number, wing size and wing shape was characterized through fluctuating asymmetry (FA) and the coefficient of variation (CV), or among-individual variation. Progeny of the GAL4/Hsp-RNAi crosses tended to have reduced trait means for both wing size and wing shape. Transcriptional knockdown of Hsp67Bc and Hsp22 significantly increased FA of bristle number, while knockdown of Hsp67Ba significantly increased FA and among-individual variation of wing shape but only in males. Suppression of Hsp67Bb expression significantly increased among-individual variation of bristle number. The knockdown of gene expression was confirmed for Hsp67Ba, Hsp67Bc, Hsp22, and Hsp67Bb. Correlation between FA and CV or among-individual variation of each trait is weak and not significant except for the case of male wing shape. CONCLUSION: Four small Hsp genes (Hsp22, Hsp67Ba, Hsp67Bb and Hsp67Bc) showed involvement in the processes of morphogenesis and developmental stability. Due to possible different functions in terms of developmental buffering of these small Hsps, phenotypic stability of an organism is probably maintained by multiple mechanisms triggered by different environmental and genetic stresses on different traits. This novel finding may lead to a better understanding of non-Hsp90 molecular mechanisms controlling variability in morphological traits.

BACKGROUND: Almost all animals, including insects, need to adapt to temperature fluctuations. The molecular basis of thermal adaptation is not well understood, although a number of candidate genes have been proposed. However, a functional link between candidate genes and thermal tolerance has rarely been established. The gene Frost (Fst) was first discovered when Drosophila flies were exposed to cold stress, but the biological function(s) of Fst has so far not been characterized. Because Fst is up-regulated after a cold stress, we tested whether it was essential for chill-coma recovery. METHODOLOGY/PRINCIPAL FINDINGS: A marked increase in Fst expression was detected (by RT-PCR) during recovery from cold stress, peaking at 42-fold after 2 h. The GAL4/UAS system was used to knock down expression of Fst and recovery ability was assessed in transgenic adults following 12 h of chill coma at 0 degrees C. The ability to recover from cold stress (short-, medium- and long-term) was significantly altered in the transgenic adults that had Fst silenced. These findings show that Fst plays an essential role in the recovery from chill coma in both males and females. CONCLUSIONS/SIGNIFICANCE: The Frost gene is essential for cold tolerance in Drosophila melanogaster and may play an important role in thermal adaptation.