Patterns of Variation in the Usage of Fatty Acid Chains among Classes of Ester and Ether Neutral Lipids and Phospholipids in the Queensland Fruit Fly.

Modern lipidomics has the power and sensitivity to elucidate the role of insects' lipidomes in their adaptations to the environment at a mechanistic molecular level. However, few lipidomic studies have yet been conducted on insects beyond model species such as Drosophila melanogaster. Here, we present the lipidome of adult males of another higher dipteran frugivore, Bactrocera tryoni. We describe 421 lipids across 15 classes of ester neutral lipids and phospholipids and ether neutral lipids and phospholipids. Most of the lipids are specified in terms of the carbon and double bond contents of each constituent hydrocarbon chain, and more ether lipids are specified to this degree than in any previous insect lipidomic analyses. Class-specific profiles of chain length and (un)saturation are broadly similar to those reported in D. melanogaster, although we found fewer medium-length chains in ether lipids. The high level of chain specification in our dataset also revealed widespread non-random combinations of different chain types in several ester lipid classes, including deficits of combinations involving chains of the same carbon and double bond contents among four phospholipid classes and excesses of combinations of dissimilar chains in several classes. Large differences were also found in the length and double bond profiles of the acyl vs. alkyl or alkenyl chains of the ether lipids. Work on other organisms suggests some of the differences observed will be functionally consequential and mediated, at least in part, by differences in substrate specificity among enzymes in lipid synthesis and remodelling pathways. Interrogation of the B. tryoni genome showed it has comparable levels of diversity overall in these enzymes but with some gene gain/loss differences and considerable sequence divergence from D. melanogaster.

Wide substrate range for a candidate bioremediation enzyme isolated from Nocardioides sp. strain SG-4 G.

Narrow substrate ranges can impact heavily on the range of applications and hence commercial viability of candidate bioremediation enzymes. Here we show that an ester hydrolase from Nocardioides strain SG-4 G has potential as a bioremediation agent against various pollutants that can be detoxified by hydrolytic cleavage of some carboxylester, carbamate, or amide linkages. Previously we showed that a radiation-killed, freeze-dried preparation (ZimA) of this strain can rapidly degrade the benzimidazole fungicide carbendazim due to the activity of a specific ester hydrolase, MheI. Here, we report that ZimA also has substantial hydrolytic activity against phthalate diesters (dimethyl, dibutyl, and dioctyl phthalate), anilide (propanil and monalide), and carbamate ester (chlorpropham) herbicides under laboratory conditions. The reaction products are substantially less toxic, or inactive as herbicides, than the parent compounds. Tests of strain SG-4 G and Escherichia coli expressing MheI found they were also able to hydrolyse dimethyl phthalate, propanil, and chlorpropham, indicating that MheI is principally responsible for the above activities.

Genetic variation for rectal gland volatiles among recently collected isofemale lines and a domesticated strain of Queensland fruit fly, Bactrocera tryoni (Diptera: Tephritidae).

Divergence between populations in mating behaviour can function as a potent premating isolating mechanism and promote speciation. However, very few cases of inherited intraspecific variation in sexual signalling have been reported in tephritid fruit flies, despite them being a highly speciose family. We tested for such variation in one tephritid, the Queensland fruit fly, Bactrocera tryoni (Qfly). Qfly mating behaviour depends on volatiles secreted from male rectal glands but no role for the volatiles from female rectal glands has yet been reported. We previously detected over 100 volatile compounds in male rectal glands and identified over 30 of them. Similar numbers were recorded in females. However, many compounds showed presence/absence differences between the sexes and many others showed quantitative differences between them. Here we report inherited variation among 24 Qfly lines (23 isofemale lines established from recent field collections and one domesticated line) in the abundance of three esters, two alcohols, two amides, an aldehyde and 18 unidentified volatiles in male rectal glands. We did not find any compounds in female rectal glands that varied significantly among the lines, although this may at least partly reflect lower female sample numbers. Most of the 26 male compounds that differed between lines were more abundant in the domesticated line than any of the recently established isofemale lines, which concurs with other evidence for changes in mating behaviour during domestication of this species. There were also large differences in several of the 26 compounds among the isofemale lines, and some of these differences were associated with the regions from which the lines were collected. While some of the variation in different compounds was correlated across lines, much of it was not, implicating involvement of multiple genes. Our findings parallel reports of geographic variation in other Qfly traits and point to inherited differences in reproductive physiology that could provide a basis for evolution of premating isolation between ecotypes.

Diversity and sex differences in rectal gland volatiles of Queensland fruit fly, Bactrocera tryoni (Diptera: Tephritidae).

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.

Population differences and domestication effects on mating and remating frequencies in Queensland fruit fly.

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.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation.

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.

Non-additive gene interactions underpin molecular and phenotypic responses in honey bee larvae exposed to imidacloprid and thymol.

Understanding the cumulative risk of chemical mixtures at environmentally realistic concentrations is a key challenge in honey bee ecotoxicology. Ecotoxicogenomics, including transcriptomics, measures responses in individual organisms at the molecular level which can provide insights into the mechanisms underlying phenotypic responses induced by one or more stressors and link impacts on individuals to populations. Here, fifth instar honey bee larvae were sampled from a previously reported field experiment exploring the phenotypic impacts of environmentally realistic chronic exposures of the pesticide imidacloprid (5 µg.kg-1 for six weeks) and the acaricide thymol (250 g.kg-1 applied via Apiguard gel in-hive for four weeks), both separately and in combination. RNA-seq was used to discover individual and interactive chemical effects on larval gene expression and to uncover molecular mechanisms linked to reported adult and colony phenotypes. The separate and combined treatments had distinct gene expression profiles which represented differentially affected signaling and metabolic pathways. The molecular signature of the mixture was characterised by additive interactions in canonical stress responses associated with oxidative stress and detoxification, and non-additive interactions in secondary responses including developmental, neurological, and immune pathways. Novel emergent impacts on eye development genes correlated with long-term defects in visual learning performance as adults. This is consistent with these chemicals working through independent modes of action that combine to impact common downstream pathways, and highlights the importance of establishing mechanistic links between molecular and phenotypic responses when predicting effects of chemical mixtures on ecologically relevant population outcomes.

Adaptive responses of free-living and symbiotic microalgae to simulated future ocean conditions.

Marine microalgae are a diverse group of microscopic eukaryotic and prokaryotic organisms capable of photosynthesis. They are important primary producers and carbon sinks but their physiology and persistence are severely affected by global climate change. Powerful experimental evolution technologies are being used to examine the potential of microalgae to respond adaptively to current and predicted future conditions, as well as to develop resources to facilitate species conservation and restoration of ecosystem functions. This review synthesizes findings and insights from experimental evolution studies of marine microalgae in response to elevated temperature and/or pCO2 . Adaptation to these environmental conditions has been observed in many studies of marine dinoflagellates, diatoms and coccolithophores. An enhancement in traits such as growth and photo-physiological performance and an increase in upper thermal limit have been shown to be possible, although the extent and rate of change differ between microalgal taxa. Studies employing multiple monoclonal replicates showed variation in responses among replicates and revealed the stochasticity of mutations. The work to date is already providing valuable information on species' climate sensitivity or resilience to managers and policymakers but extrapolating these insights to ecosystem- and community-level impacts continues to be a challenge. We recommend future work should include in situ experiments, diurnal and seasonal fluctuations, multiple drivers and multiple starting genotypes. Fitness trade-offs, stable versus plastic responses and the genetic bases of the changes also need investigating, and the incorporation of genome resequencing into experimental designs will be invaluable.

Separating two tightly linked species-defining phenotypes in Bactrocera with hybrid recombinant analysis.

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.

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication.

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.

Cuticular Chemistry of the Queensland Fruit Fly Bactrocera tryoni (Froggatt).

The cuticular layer of the insect exoskeleton contains diverse compounds that serve important biological functions, including the maintenance of homeostasis by protecting against water loss, protection from injury, pathogens and insecticides, and communication. Bactrocera tryoni (Froggatt) is the most destructive pest of fruit production in Australia, yet there are no published accounts of this species' cuticular chemistry. We here provide a comprehensive description of B. tryoni cuticular chemistry. We used gas chromatography-mass spectrometry to identify and characterize compounds in hexane extracts of B. tryoni adults reared from larvae in naturally infested fruits. The compounds found included spiroacetals, aliphatic amides, saturated/unsaturated and methyl branched C12 to C20 chain esters and C29 to C33 normal and methyl-branched alkanes. The spiroacetals and esters were found to be specific to mature females, while the amides were found in both sexes. Normal and methyl-branched alkanes were qualitatively the same in all age and sex groups but some of the alkanes differed in amounts (as estimated from internal standard-normalized peak areas) between mature males and females, as well as between mature and immature flies. This study provides essential foundations for studies investigating the functions of cuticular chemistry in this economically important species.

Genome-wide patterns of differentiation over space and time in the Queensland fruit fly.

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.

Cellular and Structural Basis of Synthesis of the Unique Intermediate Dehydro-F420-0 in Mycobacteria.

F420 is a low-potential redox cofactor used by diverse bacteria and archaea. In mycobacteria, this cofactor has multiple roles, including adaptation to redox stress, cell wall biosynthesis, and activation of the clinical antitubercular prodrugs pretomanid and delamanid. A recent biochemical study proposed a revised biosynthesis pathway for F420 in mycobacteria; it was suggested that phosphoenolpyruvate served as a metabolic precursor for this pathway, rather than 2-phospholactate as long proposed, but these findings were subsequently challenged. In this work, we combined metabolomic, genetic, and structural analyses to resolve these discrepancies and determine the basis of F420 biosynthesis in mycobacterial cells. We show that, in whole cells of Mycobacterium smegmatis, phosphoenolpyruvate rather than 2-phospholactate stimulates F420 biosynthesis. Analysis of F420 biosynthesis intermediates present in M. smegmatis cells harboring genetic deletions at each step of the biosynthetic pathway confirmed that phosphoenolpyruvate is then used to produce the novel precursor compound dehydro-F420-0. To determine the structural basis of dehydro-F420-0 production, we solved high-resolution crystal structures of the enzyme responsible (FbiA) in apo-, substrate-, and product-bound forms. These data show the essential role of a single divalent cation in coordinating the catalytic precomplex of this enzyme and demonstrate that dehydro-F420-0 synthesis occurs through a direct substrate transfer mechanism. Together, these findings resolve the biosynthetic pathway of F420 in mycobacteria and have significant implications for understanding the emergence of antitubercular prodrug resistance.IMPORTANCE Mycobacteria are major environmental microorganisms and cause many significant diseases, including tuberculosis. Mycobacteria make an unusual vitamin-like compound, F420, and use it to both persist during stress and resist antibiotic treatment. Understanding how mycobacteria make F420 is important, as this process can be targeted to create new drugs to combat infections like tuberculosis. In this study, we show that mycobacteria make F420 in a way that is different from other bacteria. We studied the molecular machinery that mycobacteria use to make F420, determining the chemical mechanism for this process and identifying a novel chemical intermediate. These findings also have clinical relevance, given that two new prodrugs for tuberculosis treatment are activated by F420.

Oxidative Catabolism of (+)-Pinoresinol Is Initiated by an Unusual Flavocytochrome Encoded by Translationally Coupled Genes within a Cluster of (+)-Pinoresinol-Coinduced Genes in Pseudomonas sp. Strain SG-MS2.

Burkholderia sp. strain SG-MS1 and Pseudomonas sp. strain SG-MS2 have previously been found to mineralize (+)-pinoresinol through a common catabolic pathway. Here, we used comparative genomics, proteomics, protein semipurification, and heterologous expression to identify a flavoprotein from the vanillyl alcohol oxidase/p-cresol methyl hydroxylase (VAO/PCMH) enzyme family in SG-MS2 that carries out the initial hydroxylation of (+)-pinoresinol at the benzylic carbon. The cognate gene is translationally coupled with a downstream cytochrome gene, and the cytochrome is required for activity. The flavoprotein has a unique combination of cofactor binding and cytochrome requirements for the VAO/PCMH family. The heterologously expressed enzyme has a Km of 1.17 µM for (+)-pinoresinol. The enzyme is overexpressed in strain SG-MS2 upon exposure to (+)-pinoresinol, along with 45 other proteins, 22 of which were found to be encoded by genes in an approximately 35.1-kb cluster also containing the flavoprotein and cytochrome genes. Homologs of 18 of these 22 genes, plus the flavoprotein and cytochrome genes, were also found in a 38.7-kb cluster in SG-MS1. The amino acid identities of four of the other proteins within the SG-MS2 cluster suggest they catalyze conversion of hydroxylated pinoresinol to protocatechuate and 2-methoxyhydroquinone. Nine other proteins upregulated in SG-MS2 on exposure to (+)-pinoresinol appear to be homologs of proteins known to comprise the protocatechuate and 2-methoxyhydroquinone catabolic pathways, but only three of the cognate genes lie within the cluster containing the flavoprotein and cytochrome genes.IMPORTANCE (+)-Pinoresinol is an important plant defense compound, a major food lignan for humans and some other animals, and the model compound used to study degradation of the ß-ß' linkages in lignin. We report a gene cluster, in one strain each of Pseudomonas and Burkholderia, that is involved in the oxidative catabolism of (+)-pinoresinol. The flavoprotein component of the α-hydroxylase which heads the pathway belongs to the 4-phenol oxidizing (4PO) subgroup of the vanillyl alcohol oxidase/p-cresol methyl hydroxylase (VAO/PCMH) enzyme family but constitutes a novel combination of cofactor and electron acceptor properties for the family. It is translationally coupled with a cytochrome gene whose product is also required for activity. The work casts new light on the biology of (+)-pinoresinol and its transformation to other bioactive molecules. Potential applications of the findings include new options for deconstructing lignin into useful chemicals and the generation of new phytoestrogenic enterolactones from lignans.

Detoxifying enzyme complements and host use phenotypes in 160 insect species.

We use the genomes of 160 insect species to test the hypothesis that the size of detoxifying enzyme families is greater in species using more chemically diverse food resources. Phylogenetically appropriate contrasts in subsamples of the data generally support the hypothesis. We find relatively high numbers of cytochrome P450, glutathione S-transferase and carboxyl/choline esterase genes in omnivores and herbivores feeding on chemically complex tissues and relatively low numbers of these genes in specialists on relatively simple diets, including plant sap, nectar and pollen, and blood. Among Lepidoptera feeding on green plant tissue and Condylognatha feeding on sap we also find more of these genes in highly polyphagous species, many of which are major agricultural pests. These genomic signatures of food resource use are consistent with the hypothesis that some taxa are preadapted for insecticide resistance evolution.

Genomic changes associated with adaptation to arid environments in cactophilic Drosophila species.

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.

The molecular basis for the neofunctionalization of the juvenile hormone esterase duplication in Drosophila.

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted ß-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.

Detoxification Genes Differ Between Cactus-, Fruit-, and Flower-Feeding Drosophila.

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.