Functional genomics of cactus host shifts in Drosophila mojavensis

Enhancing drought stress tolerance in Camelina (Camelina sativa L.) through exogenous application of potassium

The current study was performed under controlled conditions to study the effects of exogenous potassium application on carotenoid contents and drought tolerance in Camelina. Water deficit levels such as 100% FC (control) and 40% FC (drought stress) were imposed after germination of Camelina plants grown to maturity, and different treatments of exogenous K+ were applied at the vegetative stage. We have reported 17 traits of plant growth, physiology, antioxidant enzyme activity, focusing on carotenoids in Camelina to explore their potential yield and yield components. For this purpose, we used multivariate analysis techniques (descriptive statistics, correlation matrix, analysis of variance [ANOVA] and principal components analysis [PCA] to determine the best relation between potassium and studied traits). The results showed a large number of variations in the studied trait under control and water deficit condition. Plant fresh weight (g) was negatively correlated with shoot length and SOD insignificantly correlated with plant fresh weight (g) under water deficit conditions. Potassium loading predicted that foliar application (3 mM K₂ SO₄ ), foliar application (6 mM KNO₃ ), foliar application (12 mM KNO₃ ) and foliar application (12 mM K₂ SO₄ ) are the important doses that contribute the most to enhance the growth, physiological and biochemical activities and carotenoids to improve the Camelina yield under water deficit condition. These doses should be considered in the future to improve the Camelina yield under semi-arid conditions with increased genetic diversity (varietal selection).

Overexpression of NDR1 leads to pathogen resistance at elevated temperatures

Abiotic and biotic environments influence a myriad of plant-related processes, including growth, development, and the establishment and maintenance of interaction(s) with microbes. In the case of the latter, elevated temperature has been shown to be a key factor that underpins host resistance and pathogen virulence. In this study, we elucidate a role for Arabidopsis NON-RACE-SPECIFIC DISEASE RESISTANCE1 (NDR1) by exploiting effector-triggered immunity to define the regulation of plant host immunity in response to both pathogen infection and elevated temperature. We generated time-series RNA sequencing data of WT Col-0, an NDR1 overexpression line, and ndr1 and ics1-2 mutant plants under elevated temperature. Not surprisingly, the NDR1-overexpression line showed genotype-specific gene expression changes related to defense response and immune system function. The results described herein support a role for NDR1 in maintaining cell signaling during simultaneous exposure to elevated temperature and avirulent pathogen stressors.

Acuity and summation strategies differ in vinegar and desert fruit flies

An animal's vision depends on terrain features that limit the amount and distribution of available light. Approximately 10,000 years ago, vinegar flies (Drosophila melanogaster) transitioned from a single plant specialist into a cosmopolitan generalist. Much earlier, desert flies (D. mojavensis) colonized the New World, specializing on rotting cactuses in southwest North America. Their desert habitats are characteristically flat, bright, and barren, implying environmental differences in light availability. Here, we demonstrate differences in eye morphology and visual motion perception under three ambient light levels. Reducing ambient light from 35 to 18 cd/m2 causes sensitivity loss in desert but not vinegar flies. However, at 3 cd/m2, desert flies sacrifice spatial and temporal acuity more severely than vinegar flies to maintain contrast sensitivity. These visual differences help vinegar flies navigate under variably lit habitats around the world and desert flies brave the harsh desert while accommodating their crepuscular lifestyle.

WIND transcription factors orchestrate wound-induced callus formation, vascular reconnection and defense response in Arabidopsis

Wounding triggers de novo organogenesis, vascular reconnection and defense response but how wound stress evoke such a diverse array of physiological responses remains unknown. We previously identified AP2/ERF transcription factors, WOUND INDUCED DEDIFFERENTIATION1 (WIND1) and its homologs, WIND2, WIND3 and WIND4, as key regulators of wound-induced cellular reprogramming in Arabidopsis. To understand how WIND transcription factors promote downstream events, we performed time-course transcriptome analyses after WIND1 induction. We observed a significant overlap between WIND1-induced genes and genes implicated in cellular reprogramming, vascular formation and pathogen response. We demonstrated that WIND transcription factors induce several reprogramming genes to promote callus formation at wound sites. We, in addition, showed that WIND transcription factors promote tracheary element formation, vascular reconnection and resistance to Pseudomonas syringae pv. tomato DC3000. These results indicate that WIND transcription factors function as key regulators of wound-induced responses by promoting dynamic transcriptional alterations. This study provides deeper mechanistic insights into how plants control multiple physiological responses after wounding.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Soft Selective Sweep on Chemosensory Genes Correlates with Ancestral Preference for Toxic Noni in a Specialist Drosophila Population

Understanding how organisms adapt to environmental changes is a major question in evolution and ecology. In particular, the role of ancestral variation in rapid adaptation remains unclear because its trace on genetic variation, known as soft selective sweep, is often hardly recognizable from genome-wide selection scans. Here, we investigate the evolution of chemosensory genes in Drosophila yakuba mayottensis, a specialist subspecies on toxic noni (Morinda citrifolia) fruits on the island of Mayotte. We combine population genomics analyses and behavioral assays to evaluate the level of divergence in chemosensory genes and perception of noni chemicals between specialist and generalist subspecies of D. yakuba. We identify a signal of soft selective sweep on a handful of genes, with the most diverging ones involving a cluster of gustatory receptors expressed in bitter-sensing neurons. Our results highlight the potential role of ancestral genetic variation in promoting host plant specialization in herbivorous insects and identify a number of candidate genes underlying behavioral adaptation.

FRS7 and FRS12 recruit NINJA to regulate expression of glucosinolate biosynthesis genes

The sessile lifestyle of plants requires accurate physiology adjustments to be able to thrive in a changing environment. Plants integrate environmental timing signals to control developmental and stress responses. Here, we identified Far1 Related Sequence (FRS) 7 and FRS12, two transcriptional repressors that accumulate in short-day conditions, as regulators of Arabidopsis glucosinolate (GSL) biosynthesis. Loss of function of FRS7 and FRS12 results in plants with increased amplitudes of diurnal expression of GSL pathway genes. Protein interaction analyses revealed that FRS7 and FRS12 recruit the NOVEL INTERACTOR OF JAZ (NINJA) to assemble a transcriptional repressor complex. Genetic and molecular evidence demonstrated that FRS7, FRS12 and NINJA jointly regulate the expression of GSL biosynthetic genes, and thus constitute a molecular mechanism that modulates specialized metabolite accumulation.

Behavioral and Transcriptional Response to Selection for Olfactory Behavior in Drosophila

The detection, discrimination, and behavioral responses to chemical cues in the environment can have marked effects on organismal survival and reproduction, eliciting attractive or aversive behavior. To gain insight into mechanisms mediating this hedonic valence, we applied thirty generations of divergent artificial selection for Drosophila melanogaster olfactory behavior. We independently selected for positive and negative behavioral responses to two ecologically relevant chemical compounds: 2,3-butanedione and cyclohexanone. We also tested the correlated responses to selection by testing behavioral responses to other odorants and life history traits. Measurements of behavioral responses of the selected lines and unselected controls to additional odorants showed that the mechanisms underlying responses to these odorants are, in some cases, differentially affected by selection regime and generalization of the response to other odorants was only detected in the 2,3-butanedione selection lines. Food consumption and lifespan varied with selection regime and, at times, sex. An analysis of gene expression of both selection regimes identified multiple differentially expressed genes. New genes and genes previously identified in mediating olfactory behavior were identified. In particular, we found functional enrichment of several gene ontology terms, including cell-cell adhesion and sulfur compound metabolic process, the latter including genes belonging to the glutathione S-transferase family. These findings highlight a potential role for glutathione S-transferases in the evolution of hedonic valence to ecologically relevant volatile compounds and set the stage for a detailed investigation into mechanisms by which these genes mediate attraction and aversion.

Phylogenetic Relationships and Evolutionary Trends in the Cactus Family

Members of the cactus family are keystone species of arid and semiarid biomes in the Americas, as they provide shelter and resources to support other members of ecosystems. Extraordinary examples are the several species of flies of the genus Drosophila that lay eggs and feed in their rotting stems, which provide a model system for studying evolutionary processes. Although there is significant progress in understanding the evolution of Drosophila species, there are gaps in our knowledge about the cactus lineages hosting them. Here, we review the current knowledge about the evolution of Cactaceae, focusing on phylogenetic relationships and trends revealed by the study of DNA sequence data. During the last several decades, the availability of molecular phylogenies has considerably increased our understanding of the relationships, biogeography, and evolution of traits in the family. Remarkably, although succulent cacti have very low growth rates and long generation times, they underwent some of the fastest diversifications observed in the plant kingdom, possibly fostered by strong ecological interactions. We have a better understanding of the reproductive biology, population structure and speciation mechanisms in different clades. The recent publication of complete genomes for some species has revealed the importance of phenomena such as incomplete lineage sorting. Hybridization and polyploidization are common in the family, and have been studied using a variety of phylogenetic methods. We discuss potential future avenues for research in Cactaceae, emphasizing the need of a concerted effort among scientists in the Americas, together with the analyses of data from novel sequencing techniques.

Host Preference and Olfaction in Drosophila mojavensis

Many organisms live in complex environments that vary geographically in resource availability. This environmental heterogeneity can lead to changes within species in their phenotypic traits. For example, in many herbivorous insects, variation in host plant availability has been shown to influence insect host preference behavior. This behavior can be mediated in part through the insect olfactory system and the odor-evoked responses of olfactory sensory neurons (OSNs), which are in turn mediated by their corresponding odorant receptor genes. The desert dwelling fly Drosophila mojavensis is a model species for understanding the mechanisms underlying host preference in a heterogeneous environment. Depending on geographic region, one to multiple host plant species are available. Here, we conducted electrophysiological studies and found variation in responses of ORNs to host plant volatiles both within and between 2 populations-particularly to the odorant 4-methylphenol. Flies from select localities within each population were found to lack a response to 4-methylphenol. Experiments then assessed the extent to which these electrophysiological differences were associated with differences in several odor-mediated behavioral responses. No association between the presence/absence of these odor-evoked responses and short range olfactory behavior or oviposition behavior was observed. However, differences in odor-induced feeding behavior in response to 4-methylphenol were found. Localities that exhibit an odor-evoked response to the odorant had increased feeding behavior in the presence of the odorant. This study sets the stage for future work examining the functional genetics underlying variation in odor perception.

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.

Functional genomics analysis reveals two novel genes required for littorine biosynthesis

Some medicinal plants of the Solanaceae produce pharmaceutical tropane alkaloids (TAs), such as hyoscyamine and scopolamine. Littorine is a key biosynthetic intermediate in the hyoscyamine and scopolamine biosynthetic pathways. However, the mechanism underlying littorine formation from the precursors phenyllactate and tropine is not completely understood. Here, we report the elucidation of littorine biosynthesis through a functional genomics approach and functional identification of two novel biosynthesis genes that encode phenyllactate UDP-glycosyltransferase (UGT1) and littorine synthase (LS). UGT1 and LS are highly and specifically expressed in Atropa belladonna secondary roots. Suppression of either UGT1 or LS disrupted the biosynthesis of littorine and its TA derivatives (hyoscyamine and scopolamine). Purified His-tagged UGT1 catalysed phenyllactate glycosylation to form phenyllactylglucose. UGT1 and LS co-expression in tobacco leaves led to littorine synthesis if tropine and phenyllactate were added. This identification of UGT1 and LS provides the missing link in littorine biosynthesis. The results pave the way for producing hyoscyamine and scopolamine for medical use by metabolic engineering or synthetic biology.

Derived esterase activity in Drosophila sechellia contributes to evolved octanoic acid resistance

The dietary specialist fruit fly Drosophila sechellia has evolved resistance to the secondary defence compounds produced by the fruit of its host plant, Morinda citrifolia. The primary chemicals that contribute to lethality of M. citrifolia are the medium-chain fatty acids octanoic acid (OA) and hexanoic acid. At least five genomic regions contribute to this adaptation in D. sechellia and whereas the fine-mapped major effect locus for OA resistance on chromosome 3R has been thoroughly analysed, the remaining four genomic regions that contribute to toxin resistance remain uncharacterized. To begin to identify the genetic basis of toxin resistance in this species, we removed the function of well-known detoxification gene families to determine whether they contribute to toxin resistance. Previous work found that evolution of cytochrome P450 enzymatic activity or expression is not responsible for the OA resistance in D. sechellia. Here, we tested the role of the two other major detoxification gene families in resistance to Morinda fruit toxins - glutathione-S-transferases and esterases - through the use of the pesticide synergists diethyl maleate and tribufos that inhibit the function of these gene families. This work suggests that one or more esterase(s) contribute to evolved OA resistance in D. sechellia.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Genomic analysis of the four ecologically distinct cactus host populations of Drosophila mojavensis

BACKGROUND

Relationships between an organism and its environment can be fundamental in the understanding how populations change over time and species arise. Local ecological conditions can shape variation at multiple levels, among these are the evolutionary history and trajectories of coding genes. This study examines the rate of molecular evolution at protein-coding genes throughout the genome in response to host adaptation in the cactophilic Drosophila mojavensis. These insects are intimately associated with cactus necroses, developing as larvae and feeding as adults in these necrotic tissues. Drosophila mojavensis is composed of four isolated populations across the deserts of western North America and each population has adapted to utilize different cacti that are chemically, nutritionally, and structurally distinct.

RESULTS

High coverage Illumina sequencing was performed on three previously unsequenced populations of D. mojavensis. Genomes were assembled using the previously sequenced genome of D. mojavensis from Santa Catalina Island (USA) as a template. Protein coding genes were aligned across all four populations and rates of protein evolution were determined for all loci using a several approaches.

CONCLUSIONS

Loci that exhibited elevated rates of molecular evolution tend to be shorter, have fewer exons, low expression, be transcriptionally responsive to cactus host use and have fixed expression differences across the four cactus host populations. Fast evolving genes were involved with metabolism, detoxification, chemosensory reception, reproduction and behavior. Results of this study give insight into the process and the genomic consequences of local ecological adaptation.

Pyrethric acid of natural pyrethrin insecticide: complete pathway elucidation and reconstitution in Nicotiana benthamiana

In the natural pesticides known as pyrethrins, which are esters produced in flowers of Tanacetum cinerariifolium (Asteraceae), the monoterpenoid acyl moiety is pyrethric acid or chrysanthemic acid. We show here that pyrethric acid is produced from chrysanthemol in six steps catalyzed by four enzymes, the first five steps occurring in the trichomes covering the ovaries and the last one occurring inside the ovary tissues. Three steps involve the successive oxidation of carbon 10 (C10) to a carboxylic group by TcCHH, a cytochrome P450 oxidoreductase. Two other steps involve the successive oxidation of the hydroxylated carbon 1 to give a carboxylic group by TcADH2 and TcALDH1, the same enzymes that catalyze these reactions in the formation of chrysanthemic acid. The ultimate result of the actions of these three enzymes is the formation of 10-carboxychrysanthemic acid in the trichomes. Finally, the carboxyl group at C10 is methylated by TcCCMT, a member of the SABATH methyltransferase family, to give pyrethric acid. This reaction occurs mostly in the ovaries. Expression in N. benthamiana plants of all four genes encoding aforementioned enzymes, together with TcCDS, a gene that encodes an enzyme that catalyzes the formation of chrysanthemol, led to the production of pyrethric acid.

Host use and host shifts in Drosophila

Over a thousand Drosophila species have radiated onto a wide range of feeding and breeding sites. These radiations involve adaptations for locating, accepting, and growing in hosts with highly differing characteristics. In a number of species, owing to the availability of sequenced genomes, particular steps in host specialization and genes that control them, are being identified. Many cases of specialization involve the ability to detoxify some component of the host. Examples include Drosophila sechellia and the octanoic acid in Morinda citrifolia, alpha-amanitin in mycophagous drosophilids, and the alkaloids in cactophilic species. Owing to the known ecologies of many species for which genomes exist, the Drosophila model system provides an unprecedented opportunity to simultaneously examine the genes underlying HOST LOCATION, HOST ACCEPTANCE and HOST USE, the types of selection acting upon them and any coevolutionary interactions among the genes underlying these steps.

An Improved Genome Assembly for Drosophila navojoa, the Basal Species in the mojavensis Cluster

Three North American cactophilic Drosophila species, D. mojavensis, D. arizonae, and D. navojoa, are of considerable evolutionary interest owing to the shift from breeding in Opuntia cacti to columnar species. The 3 species form the "mojavensis cluster" of Drosophila. The genome of D. mojavensis was sequenced in 2007 and the genomes of D. navojoa and D. arizonae were sequenced together in 2016 using the same technology (Illumina) and assembly software (AllPaths-LG). Yet, unfortunately, the D. navojoa genome was considerably more fragmented and incomplete than its sister species, rendering it less useful for evolutionary genetic studies. The D. navojoa read dataset does not fully meet the strict insert size required by the assembler used (AllPaths-LG) and this incompatibility might explain its assembly problems. Accordingly, when we re-assembled the genome of D. navojoa with the SPAdes assembler, which does not have the strict AllPaths-LG requirements, we obtained a substantial improvement in all quality indicators such as N50 (from 84 kb to 389 kb) and BUSCO coverage (from 77% to 97%). Here we share a new, improved reference assembly for D. navojoa genome, along with a RNAseq transcriptome. Given the basal relationship of the Opuntia breeding D. navojoa to the columnar breeding D. arizonae and D. mojavensis, the improved assembly and annotation will allow researchers to address a range of questions associated with the genomics of host shifts, chromosomal rearrangements and speciation in this group.

Pupal emergence pattern in cactophilic Drosophila and the effect of host plants

Drosophila buzzatii and D. koepferae are sibling cactophilic species. The former breeds primarily on prickly pears (genus Opuntia) whereas the latter breeds on columnar cacti of the genera Cereus and Trichocereus, although with certain degree of niche overlapping. We examined the interspecific differences in diurnal temporal patterns of adult emergence from puparia and evaluated whether this behavior is affected by rearing in the different cactus hosts available in nature. We detected important host-dependent genetic variation for this trait differentially affecting the emergence schedule of these species. Diurnal pattern of emergence time was directly correlated with developmental time and negatively correlated with adult wing size, suggesting that early emergences are at least indirectly correlated with increased fitness. We discussed our results in terms of their putative effects on fitness and the genetic-metabolic pathways that would be presumably affected by host's nutritional-chemical differences.

Positive selection at sites of chemosensory genes is associated with the recent divergence and local ecological adaptation in cactophilic Drosophila

Background

Adaptation to new hosts in phytophagous insects often involves mechanisms of host recognition by genes of sensory pathways. Most often the molecular evolution of sensory genes has been explained in the context of the birth-and-death model. The role of positive selection is less understood, especially associated with host adaptation and specialization. Here we aim to contribute evidence for this latter hypothesis by considering the case of Drosophila mojavensis, a species with an evolutionary history shaped by multiple host shifts in a relatively short time scale, and its generalist sister species, D. arizonae.

Results

We used a phylogenetic and population genetic analysis framework to test for positive selection in a subset of four chemoreceptor genes, one gustatory receptor (Gr) and three odorant receptors (Or), for which their expression has been previously associated with host shifts. We found strong evidence of positive selection at several amino acid sites in all genes investigated, most of which exhibited changes predicted to cause functional effects in these transmembrane proteins. A significant portion of the sites identified as evolving positively were largely found in the cytoplasmic region, although a few were also present in the extracellular domains.

Conclusions

The pattern of substitution observed suggests that some of these changes likely had an effect on signal transduction as well as odorant recognition and protein-protein interactions. These findings support the role of positive selection in shaping the pattern of variation at chemosensory receptors, both during the specialization onto one or a few related hosts, but as well as during the evolution and adaptation of generalist species into utilizing several hosts.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1250-x) contains supplementary material, which is available to authorized users.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Behavioral evolution accompanying host shifts in cactophilic Drosophila larvae

For plant utilizing insects, the shift to a novel host is generally accompanied by a complex set of phenotypic adaptations. Many such adaptations arise in response to differences in plant chemistry, competitive environment, or abiotic conditions. One less well-understood factor in the evolution of phytophagous insects is the selective environment provided by plant shape and volume. Does the physical structure of a new plant host favor certain phenotypes? Here, we use cactophilic Drosophila, which have colonized the necrotic tissues of cacti with dramatically different shapes and volumes, to examine this question. Specifically, we analyzed two behavioral traits in larvae, pupation height, and activity that we predicted might be related to the ability to utilize variably shaped hosts. We found that populations of D. mojavensis living on lengthy columnar or barrel cactus hosts have greater activity and pupate higher in a laboratory environment than populations living on small and flat prickly pear cactus cladodes. Crosses between the most phenotypically extreme populations suggest that the genetic architectures of these behaviors are distinct. A comparison of activity in additional cactophilic species that are specialized on small and large cactus hosts shows a consistent trend. Thus, we suggest that greater motility and an associated tendency to pupate higher in the laboratory are potential larval adaptations for life on a large plant where space is more abundant and resources may be more sparsely distributed.

Drosophila melanogaster as a Model for Diabetes Type 2 Progression

Drosophila melanogaster has been used as a very versatile and potent model in the past few years for studies in metabolism and metabolic disorders, including diabetes types 1 and 2. Drosophila insulin signaling, despite having seven insulin-like peptides with partially redundant functions, is very similar to the human insulin pathway and has served to study many different aspects of diabetes and the diabetic state. Yet, very few studies have addressed the chronic nature of diabetes, key for understanding the full-blown disease, which most studies normally explore. One of the advantages of having Drosophila mutant viable combinations at different levels of the insulin pathway, with significantly reduced insulin pathway signaling, is that the abnormal metabolic state can be studied from the onset of the life cycle and followed throughout. In this review, we look at the chronic nature of impaired insulin signaling. We also compare these results to the results gleaned from vertebrate model studies.

Moving Speciation Genetics Forward: Modern Techniques Build on Foundational Studies in Drosophila

The question of how new species evolve has been examined at every level, from macroevolutionary patterns of diversification to molecular population genetic analyses of specific genomic regions between species pairs. Drosophila has been at the center of many of these research efforts. Though our understanding of the speciation process has grown considerably over the past few decades, very few genes have been identified that contribute to barriers to reproduction. The development of advanced molecular genetic and genomic methods provides promising avenues for the rapid discovery of more genes that contribute to speciation, particularly those involving prezygotic isolation. The continued expansion of tools and resources, especially for species other than Drosophila melanogaster, will be most effective when coupled with comparative approaches that reveal the genetic basis of reproductive isolation across a range of divergence times. Future research programs in Drosophila have high potential to answer long-standing questions in speciation. These include identifying the selective forces that contribute to divergence between populations and the genetic basis of traits that cause reproductive isolation. The latter can be expanded upon to understand how the genetic basis of reproductive isolation changes over time and whether certain pathways and genes are more commonly involved.

The chemical signatures underlying host plant discrimination by aphids

The diversity of phytophagous insects is largely attributable to speciation involving shifts between host plants. These shifts are mediated by the close interaction between insects and plant metabolites. However, there has been limited progress in understanding the chemical signatures that underlie host preferences. We use the pea aphid (Acyrthosiphon pisum) to address this problem. Host-associated races of pea aphid discriminate between plant species in race-specific ways. We combined metabolomic profiling of multiple plant species with behavioural tests on two A. pisum races, to identify metabolites that explain variation in either acceptance or discrimination. Candidate compounds were identified using tandem mass spectrometry. Our results reveal a small number of compounds that explain a large proportion of variation in the differential acceptability of plants to A. pisum races. Two of these were identified as L-phenylalanine and L-tyrosine but it may be that metabolically-related compounds directly influence insect behaviour. The compounds implicated in differential acceptability were not related to the set correlated with general acceptability of plants to aphids, regardless of host race. Small changes in response to common metabolites may underlie host shifts. This study opens new opportunities for understanding the mechanistic basis of host discrimination and host shifts in insects.

Esterase Profile in Drosophila mercatorum pararepleta (Diptera; Drosophilidae), a Non-cactophilic Species of the repleta Group: Development Patterns and Aspects of Genetic Variability

Luciana Paes de Barros Machado, Natalia Silva Alves, Jaqueline de Oliveira Prestes, Gabriela Ronchi Salomón, Daiane Biegai, Thais Wouk, and Rogério Pincela Mateus (2017) Esterases are a diversified group of isozymes that performs several metabolic functions in Drosophila. In the D. repleta group, this class of enzymes was well described in cactophilic species, existing a lack of studies considering substrate speci city and life cycle expression in the non-cactophilic species. The larvae of cactophilic species of the D. repleta group develop in rotting cacti cladodes, but adults are generalists. Thus, different patterns expression can be found for esterases throughout development. In this work we analyzed esterase pro le and substrate speci city during development, and genetic variability aspects in D. mercatorum pararepleta, a non-cactophilic and generalist species of D. repleta group that was understudied hitherto. Samples of 3rd (F3) and 104th (F104) generations of three D. mercatorum pararepleta strains, obtained after collections in xerophytic enclaves of southeastern Brazil (ITI and SER in São Paulo state and RIP in Paraná state), and of D33 strain (obtained from Cristalina-GO, Midwest of Brazil, and established in the laboratory in 1987) were analyzed. Eight esterase loci, EST-1 to EST-8, were detected. EST-1 and EST-2 were adult exclusive. Only EST-3 and EST-8 were monomorphic; all the others presented between two (EST-6) and six (EST-7) alleles. EST-7 was the only dimeric locus and also the only one that showed to be a preferably β-esterase regarding affinity to α- and β-naphthyl acetates as substrates. The other seven loci were divided into three classes: α-esterase exclusive (EST-2); preferably α-esterase (EST-3, EST-4, EST-5 and EST-8); and α/β-esterase (EST-1 and EST-6). The EST-3, EST-5 and EST-6 loci were not detected in all samples, suggesting that they could have become pseudogenes due to the mutation accumulation after the gene duplication. The allele frequency of EST-7 locus, which showed the highest number of alleles, in adults of D33 and SER-F3 evidenced a higher variability and diversity in the oldest strain (six alleles, Ho = 0.46) than in the youngest ( five alleles, Ho = 0.26). Moreover, the analysis of SER-F104 revealed that this locus became monomorphic. The higher variability in the strain established in the laboratory at least two decades ago, together with the allele fixation in the SER-F104, indicate that the SER strain probably suffered a more severe action of founder effect/bottleneck when it was established in the laboratory and, therefore, even if the maintenance afterwards was performed using a high number of individuals, it did not assured the conservation of the existing genetic variability.

Population differences in host plant preference and the importance of yeast and plant substrate to volatile composition

Divergent selection between environments can result in changes to the behavior of an organism. In many insects, volatile compounds are a primary means by which host plants are recognized and shifts in plant availability can result in changes to host preference. Both the plant substrate and microorganisms can influence this behavior, and host plant choice can have an impact on the performance of the organism. In Drosophila mojavensis, four geographically isolated populations each use different cacti as feeding and oviposition substrates and identify those cacti by the composition of the volatile odorants emitted. Behavioral tests revealed D. mojavensis populations vary in their degree of preference for their natural host plant. Females from the Mojave population show a marked preference for their host plant, barrel cactus, relative to other cactus choices. When flies were given a choice between cacti that were not their host plant, the preference for barrel and organ pipe cactus relative to agria and prickly pear cactus was overall lower for all populations. Volatile headspace composition is influenced by the cactus substrate, microbial community, and substrate-by-microorganism interactions. Differences in viability, developmental time, thorax length, and dry body weight exist among populations and depend on cactus substrate and population-by-cactus interactions. However, no clear association between behavioral preference and performance was observed. This study highlights a complex interplay between the insect, host plant, and microbial community and the factors mediating insect host plant preference behavior.

A core function of EDS1 with PAD4 is to protect the salicylic acid defense sector in Arabidopsis immunity

Plant defenses induced by salicylic acid (SA) are vital for resistance against biotrophic pathogens. In basal and receptor-triggered immunity, SA accumulation is promoted by Enhanced Disease Susceptibility1 with its co-regulator Phytoalexin Deficient4 (EDS1/PAD4). Current models position EDS1/PAD4 upstream of SA but their functional relationship remains unclear. In a genetic and transcriptomic analysis of Arabidopsis autoimmunity caused by constitutive or conditional EDS1/PAD4 overexpression, intrinsic EDS1/PAD4 signaling properties and their relation to SA were uncovered. A core EDS1/PAD4 pathway works in parallel with SA in basal and effector-triggered bacterial immunity. It protects against disabled SA-regulated gene expression and pathogen resistance, and is distinct from a known SA-compensatory route involving MAPK signaling. Results help to explain previously identified EDS1/PAD4 regulated SA-dependent and SA-independent gene expression sectors. Plants have evolved an alternative route for preserving SA-regulated defenses against pathogen or genetic perturbations. In a proposed signaling framework, EDS1 with PAD4, besides promoting SA biosynthesis, maintains important SA-related resistance programs, thereby increasing robustness of the innate immune system.

Genetic basis of octanoic acid resistance in Drosophila sechellia: functional analysis of a fine-mapped region

Drosophila sechellia is a species of fruit fly endemic to the Seychelles islands. Unlike its generalist sister species, D. sechellia has evolved to be a specialist on the host plant Morinda citrifolia. This specialization is interesting because the plant's fruit contains secondary defence compounds, primarily octanoic acid (OA), that are lethal to most other Drosophilids. Although ecological and behavioural adaptations to this toxic fruit are known, the genetic basis for evolutionary changes in OA resistance is not. Prior work showed that a genomic region on chromosome 3R containing 18 genes has the greatest contribution to differences in OA resistance between D. sechellia and D. simulans. To determine which gene(s) in this region might be involved in the evolutionary change in OA resistance, we knocked down expression of each gene in this region in D. melanogaster with RNA interference (RNAi) (i) ubiquitously throughout development, (ii) during only the adult stage and (iii) within specific tissues. We identified three neighbouring genes in the Osiris family, Osiris 6 (Osi6), Osi7 and Osi8, that led to decreased OA resistance when ubiquitously knocked down. Tissue-specific RNAi, however, showed that decreasing expression of Osi6 and Osi7 specifically in the fat body and/or salivary glands increased OA resistance. Gene expression analyses of Osi6 and Osi7 revealed that while standing levels of expression are higher in D. sechellia, Osi6 expression is significantly downregulated in salivary glands in response to OA exposure, suggesting that evolved tissue-specific environmental plasticity of Osi6 expression may be responsible for OA resistance in D. sechellia.

Transcriptome modulation during host shift is driven by secondary metabolites in desert Drosophila

High-throughput transcriptome studies are breaking new ground to investigate the responses that organisms deploy in alternative environments. Nevertheless, much remains to be understood about the genetic basis of host plant adaptation. Here, we investigate genome-wide expression in the fly Drosophila buzzatii raised in different conditions. This species uses decaying tissues of cactus of the genus Opuntia as primary rearing substrate and secondarily, the necrotic tissues of the columnar cactus Trichocereus terscheckii. The latter constitutes a harmful host, rich in mescaline and other related phenylethylamine alkaloids. We assessed the transcriptomic responses of larvae reared in Opuntia sulphurea and T. terscheckii, with and without the addition of alkaloids extracted from the latter. Whole-genome expression profiles were massively modulated by the rearing environment, mainly by the presence of T. terscheckii alkaloids. Differentially expressed genes were mainly related to detoxification, oxidation-reduction and stress response; however, we also found genes involved in development and neurobiological processes. In conclusion, our study contributes new data onto the role of transcriptional plasticity in response to alternative rearing environments.

Nitrogen remobilization and conservation, and underlying senescence-associated gene expression in the perennial switchgrass Panicum virgatum

Improving nitrogen (N) remobilization from aboveground to underground organs during yearly shoot senescence is an important goal for sustainable production of switchgrass (Panicum virgatum) as a biofuel crop. Little is known about the genetic control of senescence and N use efficiency in perennial grasses such as switchgrass, which limits our ability to improve the process. Switchgrass aboveground organs (leaves, stems and inflorescences) and underground organs (crowns and roots) were harvested every month over a 3-yr period. Transcriptome analysis was performed to identify genes differentially expressed in various organs during development. Total N content in aboveground organs increased from spring until the end of summer, then decreased concomitant with senescence, while N content in underground organs exhibited an increase roughly matching the decrease in shoot N during fall. Hundreds of senescence-associated genes were identified in leaves and stems. Functional grouping indicated that regulation of transcription and protein degradation play important roles in shoot senescence. Coexpression networks predict important roles for five switchgrass NAC (NAM, ATAF1,2, CUC2) transcription factors (TFs) and other TF family members in orchestrating metabolism of carbohydrates, N and lipids, protein modification/degradation, and transport processes during senescence. This study establishes a molecular basis for understanding and enhancing N remobilization and conservation in switchgrass.

Variation in life history traits and transcriptome associated with adaptation to diet shifts in the ladybird Cryptolaemus montrouzieri

BACKGROUND

Despite the broad diet range of many predatory ladybirds, the mechanisms involved in their adaptation to diet shifts are not completely understood. Here, we explored how a primarily coccidophagous ladybird Cryptolaemus montrouzieri adapts to feeding on aphids.

RESULTS

Based on the lower survival rate, longer developmental time, and lower adult body weight and reproduction rate of the predator, the aphid Megoura japonica proved being less suitable to support C. montrouzieri as compared with the citrus mealybug Planococcus citri. The results indicated up-regulation of genes related to ribosome and translation in fourth instars, which may be related to their suboptimal development. Also, several genes related to biochemical transport and metabolism, and detoxification were up-regulated as a result of adaptation to the changes in nutritional and non-nutritional (toxic) components of the prey.

CONCLUSION

Our results indicated that C. montrouzieri succeeded in feeding on aphids by regulation of genes related to development, digestion and detoxification. Thus, we argue that these candidate genes are valuable for further studies of the functional evolution of ladybirds led by diet shifts.

Differential responses of the gut transcriptome to plant protein diets in farmed Atlantic salmon

BACKGROUND

The potential for alternative plant protein sources to replace limited marine ingredients in fish feeds is important for the future of the fish farming industry. However, plant ingredients in fish feeds contain antinutritional factors (ANFs) that can promote gut inflammation (enteritis) and compromise fish health. It is unknown whether enteritis induced by plant materials with notable differences in secondary metabolism is characterised by common or distinct gene expression patterns, and how using feeds with single vs mixed plant proteins may affect the gut transcriptome and fish performance. We used Atlantic salmon parr to investigate the transcriptome responses of distal gut to varying dietary levels (0-45%) of soy protein concentrate (SPC) and faba bean (Vicia faba) protein concentrate (BPC) following an 8-week feeding trial. Soybean meal (SBM) and fish meal (FM) were used as positive and negative controls for enteritis, respectively. Gene expression profiling was performed using a microarray platform developed and validated for Atlantic salmon.

RESULTS

Different plant protein materials (SPC, BPC and SBM) generated substantially different gut gene expression profiles, with relatively few transcriptomic alterations (genes, pathways and GO terms) common for all plant proteins used. When SPC and BPC were simultaneously included in the diet, they induced less extensive alterations of gut transcriptome than diets with either SPC or BPC singly, probably due to reduced levels of individual ANFs. The mixed plant protein diets were also associated with improved body composition of fish relative to the single plant protein diets, which may provide evidence for a link between the magnitude of changes in gut transcriptome and whole-animal performance.

CONCLUSIONS

Our results indicate that gut transcriptomic profiling provides a useful tool for testing the applicability of alternative protein sources for aquaculture feeds and designing diets with reduced impact of ANFs on fish health. Ultimately, understanding diet-gut interactions and intestinal homeostasis in farmed fish is important to maximise performance and to ensure that aquaculture continues to be a sustainable source of food for a growing world population.

Bacterial RNAs activate innate immunity in Arabidopsis

The common molecular patterns of microbes play a critical role in the regulation of plant innate immunity. However, little is known about the role of nucleic acids in this process in plants. We pre-infiltrated Arabidopsis leaves with total RNAs from Pseudomonas syringae pv. tomato DC3000 (Pto DC3000) and subsequently inoculated these plants with the same bacterial cells. Total Pto DC3000 RNAs pre-infiltrated into Arabidopsis leaves elicited plant immune responses against Pto DC3000. However, sheared RNAs and RNase A application failed to induce immunity, suggesting that intact bacterial RNAs function in plant innate immunity. This notion was supported by the positive regulation of superoxide anion levels, callose deposition, two mitogen-activated protein kinases and defense-related genes observed in bacterial RNA-pre-treated leaves. Intriguingly, the Pto DC3000 population was not compromised in known pattern recognition receptor mutants for chitin, flagellin and elongation factor-Tu (EF-Tu). Plant defense-related mutant analyses further revealed that bacterial RNA-elicited innate immunity was normally required for salicylic and jasmonic acid signaling. Notably, among total RNAs, the abundant bacterial RNA species 16S and 23S ribosomal RNAs were the major determinants of this response. Our findings provide evidence that bacterial RNA serves as a microbe-associated molecular pattern in plants.

Protein domain evolution is associated with reproductive diversification and adaptive radiation in the genus Eucalyptus

Eucalyptus is a pivotal genus within the rosid order Myrtales with distinct geographic history and adaptations. Comparative analysis of protein domain evolution in the newly sequenced Eucalyptus grandis genome and other rosid lineages sheds light on the adaptive mechanisms integral to the success of this genus of woody perennials. We reconstructed the ancestral domain content to elucidate the gain, loss and expansion of protein domains and domain arrangements in Eucalyptus in the context of rosid phylogeny. We used functional gene ontology (GO) annotation of genes to investigate the possible biological and evolutionary consequences of protein domain expansion. We found that protein modulation within the angiosperms occurred primarily on the level of expansion of certain domains and arrangements. Using RNA-Seq data from E. grandis, we showed that domain expansions have contributed to tissue-specific expression of tandemly duplicated genes. Our results indicate that tandem duplication of genes, a key feature of the Eucalyptus genome, has played an important role in the expansion of domains, particularly in proteins related to the specialization of reproduction and biotic and abiotic interactions affecting root and floral biology, and that tissue-specific expression of proteins with expanded domains has facilitated subfunctionalization in domain families.

The floral transcriptome of Eucalyptus grandis

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

A global phylogeny of leafmining Ectoedemia moths (Lepidoptera: Nepticulidae): exploring host plant family shifts and allopatry as drivers of speciation

BACKGROUND

Host association patterns in Ectoedemia (Lepidoptera: Nepticulidae) are also encountered in other insect groups with intimate plant relationships, including a high degree of monophagy, a preference for ecologically dominant plant families (e.g. Fagaceae, Rosaceae, Salicaceae, and Betulaceae) and a tendency for related insect species to feed on related host plant species. The evolutionary processes underlying these patterns are only partly understood, we therefore assessed the role of allopatry and host plant family shifts in speciation within Ectoedemia.

METHODOLOGY

Six nuclear and mitochondrial DNA markers with a total aligned length of 3692 base pairs were used to infer phylogenetic relationships among 92 species belonging to the subgenus Ectoedemia of the genus Ectoedemia, representing a thorough taxon sampling with a global coverage. The results support monophyletic species groups that are congruent with published findings based on morphology. We used the obtained phylogeny to explore host plant family association and geographical distribution to investigate if host shifts and allopatry have been instrumental in the speciation of these leafmining insects.

SIGNIFICANCE

We found that, even though most species within species groups commonly feed on plants from one family, shifts to a distantly related host family have occasionally occurred throughout the phylogeny and such shifts are most commonly observed towards Betulaceae. The largest radiations have occurred within species groups that feed on Fagaceae, Rosaceae, and Salicaceae. Most species are restricted to one of the seven global biogeographic regions, but within species groups representatives are commonly found in different biogeographic regions. Although we find general patterns with regard to host use and biogeography, there are differences between clades that suggest that different drivers of speciation, and perhaps drivers that we did not examine, have shaped diversity patterns in different clades.

Salicylic acid binding of mitochondrial alpha-ketoglutarate dehydrogenase E2 affects mitochondrial oxidative phosphorylation and electron transport chain components and plays a role in basal defense against tobacco mosaic virus in tomato

Salicylic acid (SA) plays a critical role in plant defense against pathogen invasion. SA-induced viral defense in plants is distinct from the pathways mediating bacterial and fungal defense and involves a specific pathway mediated by mitochondria; however, the underlying mechanisms remain largely unknown. The SA-binding activity of the recombinant tomato (Solanum lycopersicum) alpha-ketoglutarate dehydrogenase (Slα-kGDH) E2 subunit of the tricarboxylic acid (TCA) cycle was characterized. The biological role of this binding in plant defenses against tobacco mosaic virus (TMV) was further investigated via Slα-kGDH E2 silencing and transient overexpression in plants. Slα-kGDH E2 was found to bind SA in two independent assays. SA treatment, as well as Slα-kGDH E2 silencing, increased resistance to TMV. SA did not further enhance TMV defense in Slα-kGDH E2-silenced tomato plants but did reduce TMV susceptibility in Nicotiana benthamiana plants transiently overexpressing Slα-kGDH E2. Furthermore, Slα-kGDH E2-silencing-induced TMV resistance was fully blocked by bongkrekic acid application and alternative oxidase 1a silencing. These results indicated that binding by Slα-kGDH E2 of SA acts upstream of and affects the mitochondrial electron transport chain, which plays an important role in basal defense against TMV. The findings of this study help to elucidate the mechanisms of SA-induced viral defense.

Genomics of ecological adaptation in cactophilic Drosophila

Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

Identification of differentially expressed genes in female Drosophila antonietae and Drosophila meridionalis in response to host cactus odor

Background

Studies of insect-plant interactions have provided critical insights into the ecology and evolution of adaptive processes within and among species. Cactophilic Drosophila species have received much attention because larval development occurs in the necrotic tissues of cacti, and both larvae and adults feed on these tissues. Such Drosophila-cactus interactions include effects of the host plant on the physiology and behavior of the flies, especially so their nutritional status, mating condition and reproduction. The aim of this work was to compare the transcriptional responses of two species, Drosophila antonietae and Drosophila meridionalis, and identify genes potentially related to responses to odors released by their host cactus, Cereus hildmannianus. The two fly species are sympatric in most of their populations and use this same host cactus in nature.

Results

We obtained 47 unique sequences (USs) for D. antonietae in a suppression subtractive hybridization screen, 30 of these USs had matches with genes predicted for other Drosophila species. For D. meridionalis we obtained 81 USs, 46 of which were orthologous with genes from other Drosophila species. Functional information (Gene Ontology) revealed that these differentially expressed genes are related to metabolic processes, detoxification mechanisms, signaling, response to stimuli, and reproduction. The expression of 13 genes from D. meridionalis and 12 from D. antonietae were further analyzed by quantitative real time-PCR, showing that four genes were significantly overexpressed in D. antonietae and six in D. meridionalis.

Conclusions

Our results revealed the differential expression of genes related to responses to odor stimuli by a cactus, in two associated fly species. Although the majority of activated genes were similar between the two species, we also observed that certain metabolic pathways were specifically activated, especially those related to signaling pathways and detoxification mechanisms. The activation of these genes may reflect different metabolic pathways used by these flies in their interaction with this host cactus. Our findings provide insight into how the use of C. hildmannianus may have arisen independently in the two fly species, through genetic differentiation in metabolic pathways to effectively explore this cactus as a host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0191-2) contains supplementary material, which is available to authorized users.

Evidence for functional convergence in genes upregulated by herbivores ingesting plant secondary compounds

BACKGROUND

Nearly 40 years ago, Freeland and Janzen predicted that liver biotransformation enzymes dictated diet selection by herbivores. Despite decades of research on model species and humans, little is known about the biotransformation mechanisms used by mammalian herbivores to metabolize plant secondary compounds (PSCs). We investigated the independent evolution of PSC biotransformation mechanisms by capitalizing on a dramatic diet change event-the dietary inclusion of creosote bush (Larrea tridentata)-that occurred in the recent evolutionary history of two species of woodrats (Neotoma lepida and N. bryanti).

RESULTS

By comparing gene expression profiles of two populations of woodrats with evolutionary experience to creosote and one population naïve to creosote, we identified genes either induced by a diet containing creosote PSCs or constitutively higher in populations with evolutionary experience of creosote. Although only one detoxification gene (an aldo-keto reductase) was induced by both experienced populations, these populations converged upon functionally equivalent strategies to biotransform the PSCs of creosote bush by constitutively expressing aldehyde and alcohol dehydrogenases, Cytochromes P450s, methyltransferases, glutathione S-transferases and sulfotransferases. The response of the naïve woodrat population to creosote bush was indicative of extreme physiological stress.

CONCLUSIONS

The hepatic detoxification system of mammals is notoriously complex, with hundreds of known biotransformation enzymes. The comparison herein of woodrat taxa that differ in evolutionary and ecological experience with toxins in creosote bush reveals convergence in the overall strategies used by independent species after a historical shift in diet. In addition, remarkably few genes seemed to be important in this dietary shift. The research lays the requisite groundwork for future studies of specific biotransformation pathways used by woodrats to metabolize the toxins in creosote and the evolution of diet switching in woodrats. On a larger level, this work advances our understanding of the mechanisms used by mammalian herbivores to process toxic diets and illustrates the importance of the selective relationship of PSCs in shaping herbivore diversity.