The desert woodrat (Neotoma lepida) induces a diversity of biotransformation genes in response to creosote bush resin

Liver biotransformation enzymes have long been thought to enable animals to feed on diets rich in xenobiotic compounds. However, despite decades of pharmacological research in humans and rodents, little is known about hepatic gene expression in specialized mammalian herbivores feeding on toxic diets. Leveraging a recently identified population of the desert woodrat (Neotoma lepida) found to be highly tolerant to toxic creosote bush (Larrea tridentata), we explored the expression changes of suites of biotransformation genes in response to diets enriched with varying amounts of creosote resin. Analysis of hepatic RNA-seq data indicated a dose-dependent response to these compounds, including the upregulation of several genes encoding transcription factors and numerous phase I, II, and III biotransformation families. Notably, elevated expression of five biotransformation families - carboxylesterases, cytochromes P450, aldo-keto reductases, epoxide hydrolases, and UDP-glucuronosyltransferases - corresponded to species-specific duplication events in the genome, suggesting that these genes play a prominent role in N. lepida's adaptation to creosote bush. Building on pharmaceutical studies in model rodents, we propose a hypothesis for how the differentially expressed genes are involved in the biotransformation of creosote xenobiotics. Our results provide some of the first details about how these processes likely operate in the liver of a specialized mammalian herbivore.

Transcriptome Dynamics of Brassica juncea Leaves in Response to Omnivorous Beet Armyworm (Spodoptera exigua, Hübner)

Cruciferous plants manufacture glucosinolates (GSLs) as special and important defense compounds against insects. However, how insect feeding induces glucosinolates in Brassica to mediate insect resistance, and how plants regulate the strength of anti-insect defense response during insect feeding, remains unclear. Here, mustard (Brassica juncea), a widely cultivated Brassica plant, and beet armyworm (Spodoptera exigua), an economically important polyphagous pest of many crops, were used to analyze the changes in GSLs and transcriptome of Brassica during insect feeding, thereby revealing the plant-insect interaction in Brassica plants. The results showed that the content of GSLs began to significantly increase after 48 h of herbivory by S. exigua, with sinigrin as the main component. Transcriptome analysis showed that a total of 8940 DEGs were identified in mustard challenged with beet armyworm larvae. The functional enrichment results revealed that the pathways related to the biosynthesis of glucosinolate and jasmonic acid were significantly enriched by upregulated DEGs, suggesting that mustard might provide a defense against herbivory by inducing JA biosynthesis and then promoting GSL accumulation. Surprisingly, genes regulating JA catabolism and inactivation were also activated, and both JA signaling repressors (JAZs and JAMs) and activators (MYCs and NACs) were upregulated during herbivory. Taken together, our results indicate that the accumulation of GSLs regulated by JA signaling, and the regulation of active and inactive JA compound conversion, as well as the activation of JA signaling repressors and activators, collectively control the anti-insect defense response and avoid over-stunted growth in mustard during insect feeding.

A persistent major mutation in canonical jasmonate signaling is embedded in an herbivory-elicited gene network

When insect herbivores attack plants, elicitors from oral secretions and regurgitants (OS) enter wounds during feeding, eliciting defense responses. These generally require plant jasmonate (JA) signaling, specifically, a jasmonoyl-L-isoleucine (JA-Ile) burst, for their activation and are well studied in the native tobacco Nicotiana attenuata. We used intraspecific diversity captured in a 26-parent MAGIC population planted in nature and an updated genome assembly to impute natural variation in the OS-elicited JA-Ile burst linked to a mutation in the JA-Ile biosynthetic gene NaJAR4. Experiments revealed that NaJAR4 variants were associated with higher fitness in the absence of herbivores but compromised foliar defenses, with two NaJAR homologues (4 and 6) complementing each other spatially and temporally. From decade-long seed collections of natural populations, we uncovered enzymatically inactive variants occurring at variable frequencies, consistent with a balancing selection regime maintaining variants. Integrative analyses of OS-induced transcriptomes and metabolomes of natural accessions revealed that NaJAR4 is embedded in a nonlinear complex gene coexpression network orchestrating responses to OS, which we tested by silencing four hub genes in two connected coexpressed networks and examining their OS-elicited metabolic responses. Lines silenced in two hub genes (NaGLR and NaFB67) co-occurring in the NaJAR4/6 module showed responses proportional to JA-Ile accumulations; two from an adjacent module (NaERF and NaFB61) had constitutively expressed defenses with high resistance. We infer that mutations with large fitness consequences can persist in natural populations due to compensatory responses from gene networks, which allow for diversification in conserved signaling pathways and are generally consistent with predictions of an omnigene model.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Evolutionary and dietary relationships of wild mammals based on the gut microbiome

Gut microbiome influence the health and evolution of mammals and multiple factors modulate the structure and function of gut microbiome. However, the specific changes of the diets and phylogeny on the gut microbiome were unclear. Here, we compared the gut microbiome of 16 rare wild mammals. All data (>200G 16S rRNA gene sequences) were generated using a high-throughput sequencing platform. Firmicutes and Bacteroidetes were the most predominant phyla in all mammals. However, Proteobacteria was an additionally dominant phylum specifically detected in the microbiome of carnivores and omnivores. Moreover, the dominant phyla in canids were Firmicutes, Bacteroidetes, Proteobacteria, and Fusobacteria. Phylogenetic reconstructions based on the gut microbiome and mitochondrial genome of these mammals were similar. The impact of the host on the microbiome community composition was most evident when considering conspecific and congeneric relationships. Similarity clustering showed that the gut microbiome of herbivores was clustered together, and the other clade comprised both omnivores and carnivores. Collectively, these results revealed that phylogenetic relationships and diet have an important impact on the gut microbiome, and thus the gut microbiome community composition may reflect both the phylogenetic relationships and diets. This study provides valuable basic data to facilitate future efforts related to animal conservation and health.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Rapid PCR-based method for herbivore dietary evaluation using plant-specific primers

Polyphagous pests cause significant economic loss worldwide through feeding damage on various cash crops. However, their diets in agricultural landscapes remain largely unexplored. Pest dietary evaluation in agricultural fields is a challenging task currently approached through visual observation of plant feeding and microscopic identification of semi-digested plant material in pest’s guts. While molecular gut content analysis using metabarcoding approaches using universal primers (e.g., rbcl and trnL) have been successful in evaluating polyphagous pest diet, this method is relatively costly and time-consuming. Hence, there is a need for a rapid, specific, sensitive, and cost-effective method to screen for crops in the gut of pests. This is the first study to develop plant-specific primers that target various regions of their genomes, designed using a whole plant genome sequence. We selected Verticillium wilt disease resistance protein (VE-1) and pathogenesis related protein-coding genes 1–5 (PR-1-5) as our targets and designed species-specific primers for 14 important crops in the agroecosystems. Using amplicon sizes ranging from 115 to 407 bp, we developed two multiplex primer mixes that can separate nine and five plant species per PCR reaction, respectively. These two designed primer mixes provide a rapid, sensitive and specific route for polyphagous pest dietary evaluation in agroecosystems. This work will enable future research to rapidly expand our knowledge on the diet preference and range of crops that pests consume in various agroecosystems, which will help in the redesign and development of new crop rotation regimes to minimize polyphagous pest pressure and damage on crops.

The draft genome of the Asian corn borer yields insights into ecological adaptation of a devastating maize pest

The Asian corn borer (ACB) is the most devastating pest on maize in the western Pacific region of Asia. Despite broad interests in insecticide resistance, seasonal adaptation, and larval color mimicry regarding the ACB system, lacking of reference genomic information and a powerful gene editing approach have hindered the in-depth studies of these aspects. Here we present a 455.7 Mb draft genome of ACB with 98.4% completeness. Comparative genomics analysis showed an evident expansion in gene families of gustatory receptors (105), which is related to polyphagous characteristics. Based on the comparative transcriptome analysis of resistant and susceptible ACB against Bt Cry1Ab toxin, we identified 26 genes related to Cry1Ab resistance. Additionally, transcriptomics of insects exposed to conditions of low temperature and diapause (LT) vs. room temperature and diapause (RT) provided insights into the genetic mechanisms of cold adaptation. We also successfully developed an efficient CRISPR/Cas9-based genome editing system and applied it to explore the role of color pattern genes in the ecological adaptation of ACB. Taken together, our study provides a fully annotated high-quality reference genome and efficient gene editing system to realize the potential of ACB as a study system to address important biological questions such as insecticide resistance, seasonal adaptation, and coloration. These valuable genomic resources will also benefit the development of novel strategies for maize pest management.

Large-scale identification of extracellular plant miRNAs in mammals implicates their dietary intake

Extracellular microRNAs (miRNAs) have been proposed to function in cross-kingdom gene regulation. Among these, plant-derived miRNAs of dietary origin have been reported to survive the harsh conditions of the human digestive system, enter the circulatory system, and regulate gene expression and metabolic function. However, definitive evidence supporting the presence of plant-derived miRNAs of dietary origin in mammals has been difficult to obtain due to limited sample sizes. We have developed a bioinformatics pipeline (ePmiRNA_finder) that provides strident miRNA classification and applied it to analyze 421 small RNA sequencing data sets from 10 types of human body fluids and tissues and comparative samples from carnivores and herbivores. A total of 35 miRNAs were identified that map to plants typically found in the human diet and these miRNAs were found in at least one human blood sample and their abundance was significantly different when compared to samples from human microbiome or cow. The plant-derived miRNA profiles were body fluid/tissue-specific and highly abundant in the brain and the breast milk samples, indicating selective absorption and/or the ability to be transported across tissue/organ barriers. Our data provide conclusive evidence for the presence of plant-derived miRNAs as a consequence of dietary intake and their cross-kingdom regulatory function within human circulating system.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Reference Gene Selection for Analyzing the Transcription Patterns of Two Fatty Acyl-CoA Reductase Genes From Paracoccus marginatus (Hemiptera: Pseudococcidae)

Paracoccus marginatus (Hemiptera: Pseudococcidae), known as the papaya mealybug, could cause considerable yield loss of several plants. To date, there is no molecular-based study of P. marginatus. Fatty acyl-CoA reductases (FARs) are key enzymes involved in wax synthesis. In the present study, we cloned and characterized coding sequences (CDS) of two FAR genes from P. marginatus. The results showed that PmFAR1 and PmFAR2 CDS were 1,590 and 1,497 bp in length, respectively, and sequence analysis indicated that these two genes both had the conservative motifs belonging to FAR_C superfamily. Furthermore, seven candidate reference genes were analyzed for their expression stability by using common algorithms including comparative ΔCq method, geNorm, NormFinder, BestKeeper, and RefFinder. Eventually, β-actin and GAPDH were the best reference genes in evaluating the expression of those two FAR genes. We found that PmFAR1 and PmFAR2 showed distinct expression patterns in different life stages. Moreover, the transcription of PmFAR1 and PmFAR2 in P. marginatus fed on resistant cassava cultivars was significantly lower compared with those fed on susceptible ones, indicating the potential function of FAR genes in cassava resistance to P. marginatus. The present study might help in better understanding the molecular mechanism of cassava resistance to mealybug.

Adaptive divergence and post-zygotic barriers to gene flow between sympatric populations of a herbivorous mite

Plant-herbivore interactions promote the generation and maintenance of both plant and herbivore biodiversity. The antagonistic interactions between plants and herbivores lead to host race formation: the evolution of herbivore types specializing on different plant species, with restricted gene flow between them. Understanding how ecological specialization promotes host race formation usually depends on artificial approaches, using laboratory experiments on populations associated with agricultural crops. However, evidence on how host races are formed and maintained in a natural setting remains scarce. Here, we take a multidisciplinary approach to understand whether populations of the generalist spider mite Tetranychus urticae form host races in nature. We demonstrate that a host race co-occurs among generalist conspecifics in the dune ecosystem of The Netherlands. Extensive field sampling and genotyping of individuals over three consecutive years showed a clear pattern of host associations. Genome-wide differences between the host race and generalist conspecifics were found using a dense set of SNPs on field-derived iso-female lines and previously sequenced genomes of T. urticae. Hybridization between lines of the host race and sympatric generalist lines is restricted by post-zygotic breakdown, and selection negatively impacts the survival of generalists on the native host of the host race. Our description of a host race among conspecifics with a larger diet breadth shows how ecological and reproductive isolation aid in maintaining intra-specific variation in sympatry, despite the opportunity for homogenization through gene flow. Our findings highlight the importance of explicitly considering the spatial and temporal scale on which plant-herbivore interactions occur in order to identify herbivore populations associated with different plant species in nature. This system can be used to study the underlying genetic architecture and mechanisms that facilitate the use of a large range of host plant taxa by extreme generalist herbivores. In addition, it offers the chance to investigate the prevalence and mechanisms of ecological specialization in nature.

Comparative analysis of mite genomes reveals positive selection for diet adaptation

Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.

Role of digestive protease enzymes and related genes in host plant adaptation of a polyphagous pest, Spodoptera frugiperda

The evolutionary success of phytophagous insects depends on their ability to efficiently exploit plants as a source of energy for survival. Herbivorous insects largely depend on the efficiency, flexibility, and diversity of their digestive physiology and sophistication of their detoxification system to use chemically diverse host plants as food sources. The fall armyworm, Spodoptera frugiperda (J.E. Smith), is a polyphagous pest of many commercially important crops. To elucidate the ability of this insect pest to adapt to host plant mechanisms, we evaluated the impact of primary (corn) and alternate (rice) host plants after 11 generations on gut digestive enzymatic activity and expression profiles of related genes. Results indicated that the total protease and class-specific trypsin- and chymotrypsin-like protease activity of S. frugiperda significantly differed among host plant treatments. The class-specific protease profiles greatly differed in S. frugiperda midguts upon larval exposure to different treatments with inhibitors compared with treatments without inhibitors. Similarly, the single and cumulative effects of the enzyme-specific inhibitors TLCK, TPCK, and E-64 significantly increased larval mortality and reduced larval growth/mass across different plant treatments. Furthermore, the quantitative reverse transcription polymerase chain reaction results revealed increased transcription of two trypsin (SfTry-3, SfTry-7) and one chymotrypsin gene (Sfchym-9), which indicated that they have roles in host plant adaptation. Knockdown of these genes resulted in significantly reduced mRNA expression levels of the trypsin genes. This was related to the increased mortality observed in treatments compared with the dsRED control. This result indicates possible roles of S. frugiperda gut digestive enzymes and related genes in host plant adaptation.

Historic exposure to herbivores, not constitutive traits, explains plant tolerance to herbivory in the case of two Medicago species (Fabaceae)

Mechanisms that allow plants to survive and reproduce after herbivory are considered to play a key role in plant evolution. In this study, we evaluated how tolerance varies in species with different historic exposure to herbivores considering ontogeny. We exposed the range-restricted species Medicago citrina and its closely related and widespread species M. arborea to one and two herbivory simulations (80 % aerial biomass loss). Physiological and growth parameters related to tolerance capacity were assessed to evaluate constitutive values (without herbivory) and induced tolerance after damage. Constitutive traits were not always related to greater tolerance, and each species compensated for herbivory through different traits. Herbivory damage only led to mortality in M. citrina; adults exhibited root biomass loss and increased oxidative stress after damage, but also compensated aerial biomass. Despite seedlings showed a lower death percentage than adults after herbivory in M. citrina, they showed less capacity to recover control values than adults. Moderate tolerance to M. arborea herbivory and low tolerance to M. citrina is found. Thus, although the constitutive characteristics are maintained in the lineage, the tolerance of plants decreases in M. citrina. That represents how plants respond to the lack of pressure from herbivores in their habitat.

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.

MYC2-Activated TRICHOME BIREFRINGENCE-LIKE37 Acetylates Cell Walls and Enhances Herbivore Resistance

O-Acetylation of polysaccharides predominantly modifies plant cell walls by changing the physicochemical properties and, consequently, the structure and function of the cell wall. Expression regulation and specific function of cell wall-acetylating enzymes remain to be fully understood. In this report, we cloned a previously identified stunted growth mutant named sucrose uncoupled1 (sun1) in Arabidopsis (Arabidopsis thaliana). SUN1 encodes a member of the TRICHOME BIREFRINGEN-LIKE family, AtTBL37 AtTBL37 is highly expressed in fast-growing plant tissues and encodes a Golgi apparatus-localized protein that regulates secondary cell wall thickening and acetylation. In sun1, jasmonate signaling and expression of downstream chemical defense genes, including VEGETATIVE STORAGE PROTEIN1 and BRANCHED-CHAIN AMINOTRANSFERASE4, are increased but, unexpectedly, sun1 is more susceptible to insect feeding. The central transcription factor in jasmonate signaling, MYC2, binds to and induces AtTBL37 expression. MYC2 also promotes the expression of many other TBLs Moreover, MYC activity enhances cell wall acetylation. Overexpression of AtTBL37 in the myc2-2 background reduces herbivore feeding. Our study highlights the role of O-acetylation in controlling plant cell wall properties, plant development, and herbivore defense.

Comparative gut transcriptome analysis of Diatraea saccharalis in response to the dietary source

The sugarcane borer (Diatraea saccharalis, Fabricius, 1794) is a devastating pest that causes millions of dollars of losses each year to sugarcane producers by reducing sugar and ethanol yields. The control of this pest is difficult due to its endophytic behavior and rapid development. Pest management through biotechnological approaches has emerged in recent years as an alternative to currently applied methods. Genetic information about the target pests is often required to perform biotechnology-based management. The genomic and transcriptomic data for D. saccharalis are very limited. Herein, we report a tissue-specific transcriptome of D. saccharalis larvae and a differential expression analysis highlighting the physiological characteristics of this pest in response to two different diets: sugarcane and an artificial diet. Sequencing was performed on the Illumina HiSeq 2000 platform, and a de novo assembly was generated. A total of 27,626 protein-coding unigenes were identified, among which 1,934 sequences were differentially expressed between treatments. Processes such as defence, digestion, detoxification, signaling, and transport were highly represented among the differentially expressed genes (DEGs). Furthermore, seven aminopeptidase genes were identified as candidates to encode receptors of Cry proteins, which are toxins of Bacillus thuringiensis used to control lepidopteran pests. Since plant-insect interactions have produced a considerable number of adaptive responses in hosts and herbivorous insects, the success of phytophagous insects relies on their ability to overcome challenges such as the response to plant defences and the intake of nutrients. In this study, we identified metabolic pathways and specific genes involved in these processes. Thus, our data strongly contribute to the knowledge advancement of insect transcripts, which can be a source of target genes for pest management.

Tritrophic interactions follow phylogenetic escalation and climatic adaptation

One major goal in plant evolutionary ecology is to address how and why tritrophic interactions mediated by phytochemical plant defences vary across species, space, and time. In this study, we tested three classical hypotheses about plant defences: (i) the resource-availability hypothesis, (ii) the altitudinal/elevational gradient hypothesis and (iii) the defence escalation hypothesis. For this purpose, predatory soil nematodes were challenged to hunt for root herbivores based on volatile cues from damaged or intact roots of 18 Alpine Festuca grass species adapted to distinct climatic niches spanning 2000 meters of elevation. We found that adaptation into harsh, nutrient-limited alpine environments coincided with the production of specific blends of volatiles, highly attractive for nematodes. We also found that recently-diverged taxa exposed to herbivores released higher amounts of volatiles than ancestrally-diverged species. Therefore, our model provides evidence that belowground indirect plant defences associated with tritrophic interactions have evolved under two classical hypotheses in plant ecology. While phylogenetic drivers of volatile emissions point to the defence-escalation hypothesis, plant local adaptation of indirect defences is in line with the resource availability hypothesis.

The evolution and genomic basis of beetle diversity

The order Coleoptera (beetles) is arguably the most speciose group of animals, but the evolutionary history of beetles, including the impacts of plant feeding (herbivory) on beetle diversification, remain poorly understood. We inferred the phylogeny of beetles using 4,818 genes for 146 species, estimated timing and rates of beetle diversification using 89 genes for 521 species representing all major lineages and traced the evolution of beetle genes enabling symbiont-independent digestion of lignocellulose using 154 genomes or transcriptomes. Phylogenomic analyses of these uniquely comprehensive datasets resolved previously controversial beetle relationships, dated the origin of Coleoptera to the Carboniferous, and supported the codiversification of beetles and angiosperms. Moreover, plant cell wall-degrading enzymes (PCWDEs) obtained from bacteria and fungi via horizontal gene transfers may have been key to the Mesozoic diversification of herbivorous beetles-remarkably, both major independent origins of specialized herbivory in beetles coincide with the first appearances of an arsenal of PCWDEs encoded in their genomes. Furthermore, corresponding (Jurassic) diversification rate increases suggest that these novel genes triggered adaptive radiations that resulted in nearly half of all living beetle species. We propose that PCWDEs enabled efficient digestion of plant tissues, including lignocellulose in cell walls, facilitating the evolution of uniquely specialized plant-feeding habits, such as leaf mining and stem and wood boring. Beetle diversity thus appears to have resulted from multiple factors, including low extinction rates over a long evolutionary history, codiversification with angiosperms, and adaptive radiations of specialized herbivorous beetles following convergent horizontal transfers of microbial genes encoding PCWDEs.

An influential meal: host plant dependent transcriptional variation in the beet armyworm, Spodoptera exigua (Lepidoptera: Noctuidae)

BACKGROUND

To understand the genetic mechanisms of insect herbivory, the transcriptional response of insects feeding on different host plant species has to be studied. Here, we generated gene expression data of the generalist herbivore Spodoptera exigua (Hübner) feeding on three selected host plant species and a control (artificial diet). The host plant species used in this study -cabbage (Brassica oleracea), maize (Zea mays) and tobacco (Nicotiana tabacum)- are members of different plant families that each employ specific defence mechanisms and toxins.

RESULTS

Spodoptera exigua larvae had a higher growth rate, indicator for herbivore success, when feeding on Z. mays compared to larvae feeding on B. oleracea or N. tabacum. Larvae feeding on the different host plant species showed divergent transcriptional responses. We identified shared and unique gene expression patterns dependent of the host plant species the larvae fed on. Unique gene expression patterns, containing uniquely upregulated transcripts including specific detoxification genes, were found for larvae feeding on either B. oleracea or N. tabacum. No diet-specific gene cluster was identified for larvae feeding on the host for which larvae showed optimal herbivore success, Z. mays, or artificial diet. In contrast, for larvae feeding on hosts for which they showed low herbivore success, specific diet-dependent gene clusters were identified. Functional annotation of these clusters indicates that S. exigua larvae deploy particular host plant-specific genes for digestion and detoxification.

CONCLUSIONS

The lack of a host plant-specific gene activity for larvae feeding on Z. mays and the artificial diet suggest a general and non-specific gene activity for host plants with optimal herbivore success. Whereas the finding of specific gene clusters containing particular digestion and detoxifying genes expressed in larvae feeding on B. oleracea and N. tabacum, with low herbivore success, imply a host plant-specific gene activity for larvae feeding on host plants with suboptimal herbivore success. This observation leads to the conclusion that a polyphagous herbivore is able to feed on a large variation of host plants due to the flexibility and diversity of genes involved in digestion and detoxification that are deployed in response to particular host plant species.

Genome of the African cassava whitefly Bemisia tabaci and distribution and genetic diversity of cassava-colonizing whiteflies in Africa

The whitefy Bemisia tabaci, a species complex consisting of many morphologically indistinguishable species divided into distinct clades, is one of the most globally important agricultural pests and plant virus vectors. Cassava-colonizing B. tabaci transmits viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). Half of all cassava plants in Africa are affected by these viral diseases, resulting in annual production losses of more than US$ 1 billion. Here we report the draft genome of the cassava whitefly B. tabaci Sub-Saharan Africa - East and Central Africa (SSA-ECA), the super-abundant population that has been associated with the rapid spread of viruses causing the pandemics of CMD and CBSD. The SSA-ECA genome assembled from Illumina short reads has a total size of 513.7 Mb and a scaffold N50 length of 497 kb, and contains 15,084 predicted protein-coding genes. Phylogenetic analysis suggests that SSA-ECA diverged from MEAM1 around 5.26 million years ago. A comprehensive genetic analysis of cassava-colonizing B. tabaci in Africa was also conducted, in which a total of 243 whitefly specimens were collected from 18 countries representing all major cassava-growing regions in the continent and genotyped using NextRAD sequencing. Population genomic analyses confirmed the existence of six major populations linked by gene flow and inferred the distribution patterns of these populations across the African continent. The genome of SSA-ECA and the genetic findings provide valuable resources and guidance to facilitate whitefly research and the development of strategies to control cassava viral diseases spread by whiteflies.

The decoration of specialized metabolites influences stylar development

Plants produce many different specialized (secondary) metabolites that function in solving ecological challenges; few are known to function in growth or other primary processes. 17-Hydroxygeranylinalool diterpene glycosides (DTGs) are abundant herbivory-induced, structurally diverse and commonly malonylated defense metabolites in Nicotiana attenuata plants. By identifying and silencing a malonyltransferase, NaMaT1, involved in DTG malonylation, we found that DTG malonylation percentages are normally remarkably uniform, but when disrupted, result in DTG-dependent reduced floral style lengths, which in turn result from reduced stylar cell sizes, IAA contents, and YUC activity; phenotypes that could be restored by IAA supplementation or by silencing the DTG pathway. Moreover, the Nicotiana genus-specific JA-deficient short-style phenotype also results from alterations in DTG malonylation patterns. Decorations of plant specialized metabolites can be tuned to remarkably uniform levels, and this regulation plays a central but poorly understood role in controlling the development of specific plant parts, such as floral styles.

Shared Transcriptional Control and Disparate Gain and Loss of Aphid Parasitism Genes

Aphids are a diverse group of taxa that contain agronomically important species, which vary in their host range and ability to infest crop plants. The genome evolution underlying agriculturally important aphid traits is not well understood. We generated draft genome assemblies for two aphid species: Myzus cerasi (black cherry aphid) and the cereal specialist Rhopalosiphum padi. Using a de novo gene prediction pipeline on both these, and three additional aphid genome assemblies (Acyrthosiphon pisum, Diuraphis noxia, and Myzus persicae), we show that aphid genomes consistently encode similar gene numbers. We compare gene content, gene duplication, synteny, and putative effector repertoires between these five species to understand the genome evolution of globally important plant parasites. Aphid genomes show signs of relatively distant gene duplication, and substantial, relatively recent, gene birth. Putative effector repertoires, originating from duplicated and other loci, have an unusual genomic organization and evolutionary history. We identify a highly conserved effector pair that is tightly physically linked in the genomes of all aphid species tested. In R. padi, this effector pair is tightly transcriptionally linked and shares an unknown transcriptional control mechanism with a subset of ∼50 other putative effectors and secretory proteins. This study extends our current knowledge on the evolution of aphid genomes and reveals evidence for an as-of-yet unknown shared control mechanism, which underlies effector expression, and ultimately plant parasitism.

Evidence of an evolutionary hourglass pattern in herbivory-induced transcriptomic responses

Herbivory-induced defenses are specific and activated in plants when elicitors, frequently found in the herbivores' oral secretions, are introduced into wounds during attack. While complex signaling cascades are known to be involved, it remains largely unclear how natural selection has shaped the evolution of these induced defenses. We analyzed herbivory-induced transcriptomic responses in wild tobacco, Nicotiana attenuata, using a phylotranscriptomic approach that measures the origin and sequence divergence of herbivory-induced genes. Highly conserved and evolutionarily ancient genes of primary metabolism were activated at intermediate time points (2-6 h) after elicitation, while less constrained and young genes associated with defense signaling and biosynthesis of specialized metabolites were activated at early (before 2 h) and late (after 6 h) stages of the induced response, respectively - a pattern resembling the evolutionary hourglass pattern observed during embryogenesis in animals and the developmental process in plants and fungi. The hourglass patterns found in herbivory-induced defense responses and developmental process are both likely to be a result of signaling modularization and differential evolutionary constraints on the modules involved in the signaling cascade.

RNA Interference of the Ecdysone Receptor Genes EcR and USP in Grain Aphid (Sitobion avenae F.) Affects Its Survival and Fecundity upon Feeding on Wheat Plants

RNA interference (RNAi) has been widely used in functional genomics of insects and received intensive attention in the development of RNAi-based plants for insect control. Ecdysone receptor (EcR) and ultraspiracle protein (USP) play important roles in molting, metamorphosis, and reproduction of insects. EcR and USP orthologs and their function in grain aphid (Sitobion avenae F.) have not been documented yet. Here, RT-PCR, qRT-PCR, dsRNA feeding assay and aphid bioassay were employed to isolate EcR and USP orthologs in grain aphid, investigate their expression patterns, and evaluate the effect of RNAi on aphid survival and fecundity, and its persistence. The results indicated that SaEcR and SaUSP exhibited similar expression profiles at different developmental stages. Oral administration of dsRNAs of SaEcR and dsSaUSP significantly decreased the survival of aphids due to the down-regulation of these two genes, respectively. The silencing effect was persistent and transgenerational, as demonstrated by the reduced survival and fecundity due to knock-down of SaEcR and SaUSP in both the surviving aphids and their offspring, even after switching to aphid-susceptible wheat plants. Taken together, our results demonstrate that SaEcR and SaUSP are essential genes in aphid growth and development, and could be used as RNAi targets for wheat aphid control.

Horizontal Gene Transfer Contributes to the Evolution of Arthropod Herbivory

Within animals, evolutionary transition toward herbivory is severely limited by the hostile characteristics of plants. Arthropods have nonetheless counteracted many nutritional and defensive barriers imposed by plants and are currently considered as the most successful animal herbivores in terrestrial ecosystems. We gather a body of evidence showing that genomes of various plant feeding insects and mites possess genes whose presence can only be explained by horizontal gene transfer (HGT). HGT is the asexual transmission of genetic information between reproductively isolated species. Although HGT is known to have great adaptive significance in prokaryotes, its impact on eukaryotic evolution remains obscure. Here, we show that laterally transferred genes into arthropods underpin many adaptations to phytophagy, including efficient assimilation and detoxification of plant produced metabolites. Horizontally acquired genes and the traits they encode often functionally diversify within arthropod recipients, enabling the colonization of more host plant species and organs. We demonstrate that HGT can drive metazoan evolution by uncovering its prominent role in the adaptations of arthropods to exploit plants.

Jasmonic acid carboxyl methyltransferase regulates development and herbivory-induced defense response in rice

Jasmonic acid (JA) and related metabolites play a key role in plant defense and growth. JA carboxyl methyltransferase (JMT) may be involved in plant defense and development by methylating JA to methyl jasmonate (MeJA) and thus influencing the concentrations of JA and related metabolites. However, no JMT gene has been well characterized in monocotyledon defense and development at the molecular level. After we cloned a rice JMT gene, OsJMT1, whose encoding protein was localized in the cytosol, we found that the recombinant OsJMT1 protein catalyzed JA to MeJA. OsJMT1 is up-regulated in response to infestation with the brown planthopper (BPH; Nilaparvata lugens). Plants in which OsJMT1 had been overexpressed (oe-JMT plants) showed reduced height and yield. These oe-JMT plants also exhibited increased MeJA levels but reduced levels of herbivore-induced JA and jasmonoyl-isoleucine (JA-Ile). The oe-JMT plants were more attractive to BPH female adults but showed increased resistance to BPH nymphs, probably owing to the different responses of BPH female adults and nymphs to the changes in levels of H2 O2 and MeJA in oe-JMT plants. These results indicate that OsJMT1, by altering levels of JA and related metabolites, plays a role in regulating plant development and herbivore-induced defense responses in rice.

The Arabidopsis immune regulator SRFR1 dampens defences against herbivory by Spodoptera exigua and parasitism by Heterodera schachtii

Plants have developed diverse mechanisms to fine tune defence responses to different types of enemy. Cross-regulation between signalling pathways may allow the prioritization of one response over another. Previously, we identified SUPPRESSOR OF rps4-RLD1 (SRFR1) as a negative regulator of ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1)-dependent effector-triggered immunity against the bacterial pathogen Pseudomonas syringae pv. tomato strain DC3000 expressing avrRps4. The use of multiple stresses is a powerful tool to further define gene function. Here, we examined whether SRFR1 also impacts resistance to a herbivorous insect in leaves and to a cyst nematode in roots. Interestingly, srfr1-1 plants showed increased resistance to herbivory by the beet army worm Spodoptera exigua and to parasitism by the cyst nematode Heterodera schachtii compared with the corresponding wild-type Arabidopsis accession RLD. Using quantitative real-time PCR (qRT-PCR) to measure the transcript levels of salicylic acid (SA) and jasmonate/ethylene (JA/ET) pathway genes, we found that enhanced resistance of srfr1-1 plants to S. exigua correlated with specific upregulation of the MYC2 branch of the JA pathway concurrent with suppression of the SA pathway. In contrast, the greater susceptibility of RLD was accompanied by simultaneously increased transcript levels of SA, JA and JA/ET signalling pathway genes. Surprisingly, mutation of either SRFR1 or EDS1 increased resistance to H. schachtii, indicating that the concurrent presence of both wild-type genes promotes susceptibility. This finding suggests a novel form of resistance in Arabidopsis to the biotrophic pathogen H. schachtii or a root-specific regulation of the SA pathway by EDS1, and places SRFR1 at an intersection between multiple defence pathways.

Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp

Adaption to dietary changes is critical in the evolution of host plant ranges in polyphagous insects. We compared three taxa of lepidopteran herbivores from the predominantly generalist genus Spodoptera showing different degrees of polyphagy: Spodoptera littoralis, with a broad host range including both mono- and dicotyledonous plants, and two Spodoptera frugiperda strains [Corn (i.e. maize) (C) and Rice (R)] adapted primarily to different grass species. When feeding on maize we show a lower performance in the broad generalist taxon compared to the grass adapted taxa. Among these taxa, the maize adapted S. frugiperda C-strain generally performed better than the R-strain on maize leaves. On artificial pinto diet, all taxa performed well. Our RNA-Seq analysis of midgut transcriptomes from 3rd instar larvae feeding on maize showed broader transcriptional readjustments in the generalist S. littoralis compared to grass adapted S. frugiperda strains. Substantial alteration in the expression levels of midgut physiological function related transcripts, such as digestive and detoxifying enzymes, transporters, immunity, and peritrophic membrane associated transcripts, existed in all taxa. We found high background expression of UDP-glucosyl transferases, which are known to neutralize maize leaf toxins, in the maize adapted S. frugiperda C-strain, contributing to its fitness on maize compared to the R-strain. Our findings provide evidence for divergent diet specific response of digestive physiology within these Spodoptera taxa. Unexpectedly, the C- and R-strains of S. frugiperda fed on the same diet showed large differences in expression patterns between these two closely related taxa.

Meta-transcriptomics indicates biotic cross-tolerance in willow trees cultivated on petroleum hydrocarbon contaminated soil

BACKGROUND

High concentrations of petroleum hydrocarbon (PHC) pollution can be hazardous to human health and leave soils incapable of supporting agricultural crops. A cheap solution, which can help restore biodiversity and bring land back to productivity, is cultivation of high biomass yielding willow trees. However, the genetic mechanisms which allow these fast-growing trees to tolerate PHCs are as yet unclear.

METHODS

Salix purpurea 'Fish Creek' trees were pot-grown in soil from a former petroleum refinery, either lacking or enriched with C10-C50 PHCs. De novo assembled transcriptomes were compared between tree organs and impartially annotated without a priori constraint to any organism.

RESULTS

Over 45% of differentially expressed genes originated from foreign organisms, the majority from the two-spotted spidermite, Tetranychus urticae. Over 99% of T. urticae transcripts were differentially expressed with greater abundance in non-contaminated trees. Plant transcripts involved in the polypropanoid pathway, including phenylalanine ammonia-lyase (PAL), had greater expression in contaminated trees whereas most resistance genes showed higher expression in non-contaminated trees.

CONCLUSIONS

The impartial approach to annotation of the de novo transcriptomes, allowing for the possibility for multiple species identification, was essential for interpretation of the crop's response treatment. The meta-transcriptomic pattern of expression suggests a cross-tolerance mechanism whereby abiotic stress resistance systems provide improved biotic resistance. These findings highlight a valuable but complex biotic and abiotic stress response to real-world, multidimensional contamination which could, in part, help explain why crops such as willow can produce uniquely high biomass yields on challenging marginal land.

Genetics-based interactions among plants, pathogens, and herbivores define arthropod community structure

Plant resistance to pathogens or insect herbivores is common, but its potential for indirectly influencing plant-associated communities is poorly known. Here, we test whether pathogens' indirect effects on arthropod communities and herbivory depend on plant resistance to pathogens and/or herbivores, and address the overarching interacting foundation species hypothesis that genetics-based interactions among a few highly interactive species can structure a much larger community. In a manipulative field experiment using replicated genotypes of two Populus species and their interspecific hybrids, we found that genetic variation in plant resistance to both pathogens and insect herbivores modulated the strength of pathogens' indirect effects on arthropod communities and insect herbivory. First, due in part to the pathogens' differential impacts on leaf biomass among the two Populus species and the hybrids, the pathogen most strongly impacted arthropod community composition, richness, and abundance on the pathogen-susceptible tree species. Second, we found similar patterns comparing pathogen-susceptible and pathogen-resistant genotypes within species. Third, within a plant species, pathogens caused a fivefold greater reduction in herbivory on insect-herbivore-susceptible plant genotypes than on herbivore-resistant genotypes, demonstrating that the pathogen-herbivore interaction is genotype dependent. We conclude that interactions among plants, pathogens, and herbivores can structure multitrophic communities, supporting the interacting foundation species hypothesis. Because these interactions are genetically based, evolutionary changes in genetic resistance could result in ecological changes in associated communities, which may in turn feed back to affect plant fitness.

The Effect of Host-Plant Phylogenetic Isolation on Species Richness, Composition and Specialization of Insect Herbivores: A Comparison between Native and Exotic Hosts

Understanding the drivers of plant-insect interactions is still a key issue in terrestrial ecology. Here, we used 30 well-defined plant-herbivore assemblages to assess the effects of host plant phylogenetic isolation and origin (native vs. exotic) on the species richness, composition and specialization of the insect herbivore fauna on co-occurring plant species. We also tested for differences in such effects between assemblages composed exclusively of exophagous and endophagous herbivores. We found a consistent negative effect of the phylogenetic isolation of host plants on the richness, similarity and specialization of their insect herbivore faunas. Notably, except for Jaccard dissimilarity, the effect of phylogenetic isolation on the insect herbivore faunas did not vary between native and exotic plants. Our findings show that the phylogenetic isolation of host plants is a key factor that influences the richness, composition and specialization of their local herbivore faunas, regardless of the host plant origin.

Transcriptional responses of invasive and indigenous whiteflies to different host plants reveal their disparate capacity of adaptation

The whitefly Bemisia tabaci contains more than 35 cryptic species. The higher adaptability of Middle East-Asia Minor 1 (MEAM1) cryptic species has been recognized as one important factor for its invasion and displacement of other indigenous species worldwide. Here we compared the performance of the invasive MEAM1 and the indigenous Asia II 3 whitefly species following host plant transfer from a suitable host (cotton) to an unsuitable host (tobacco) and analyzed their transcriptional responses. After transfer to tobacco for 24 h, MEAM1 performed much better than Asia II 3. Transcriptional analysis showed that the patterns of gene regulation were very different with most of the genes up-regulated in MEAM1 but down-regulated in Asia II 3. Whereas carbohydrate and energy metabolisms were repressed in Asia II 3, the gene expression and protein metabolisms were activated in MEAM1. Compared to the constitutive high expression of detoxification genes in MEAM1, most of the detoxification genes were down-regulated in Asia II 3. Enzymatic activities of P450, GST and esterase further verified that the detoxification of MEAM1 was much higher than that of Asia II 3. These results reveal obvious differences in responses of MEAM1 and Asia II 3 to host transfer.

Phytophagous insects on native and non-native host plants: combining the community approach and the biogeographical approach

During the past centuries, humans have introduced many plant species in areas where they do not naturally occur. Some of these species establish populations and in some cases become invasive, causing economic and ecological damage. Which factors determine the success of non-native plants is still incompletely understood, but the absence of natural enemies in the invaded area (Enemy Release Hypothesis; ERH) is one of the most popular explanations. One of the predictions of the ERH, a reduced herbivore load on non-native plants compared with native ones, has been repeatedly tested. However, many studies have either used a community approach (sampling from native and non-native species in the same community) or a biogeographical approach (sampling from the same plant species in areas where it is native and where it is non-native). Either method can sometimes lead to inconclusive results. To resolve this, we here add to the small number of studies that combine both approaches. We do so in a single study of insect herbivory on 47 woody plant species (trees, shrubs, and vines) in the Netherlands and Japan. We find higher herbivore diversity, higher herbivore load and more herbivory on native plants than on non-native plants, generating support for the enemy release hypothesis.

Environmental RNAi in herbivorous insects

Environmental RNAi (eRNAi) is a sequence-specific regulation of endogenous gene expression in a receptive organism by exogenous double-stranded RNA (dsRNA). Although demonstrated under artificial dietary conditions and via transgenic plant presentations in several herbivorous insects, the magnitude and consequence of exogenous dsRNA uptake and the role of eRNAi remains unknown under natural insect living conditions. Our analysis of coleopteran insects sensitive to eRNAi fed on wild-type plants revealed uptake of plant endogenous long dsRNAs, but not small RNAs. Subsequently, the dsRNAs were processed into 21 nt siRNAs by insects and accumulated in high quantities in insect cells. No accumulation of host plant-derived siRNAs was observed in lepidopteran larvae that are recalcitrant to eRNAi. Stability of ingested dsRNA in coleopteran larval gut followed by uptake and transport from the gut to distal tissues appeared to be enabling factors for eRNAi. Although a relatively large number of distinct coleopteran insect-processed plant-derived siRNAs had sequence complementarity to insect transcripts, the vast majority of the siRNAs were present in relatively low abundance, and RNA-seq analysis did not detect a significant effect of plant-derived siRNAs on insect transcriptome. In summary, we observed a broad genome-wide uptake of plant endogenous dsRNA and subsequent processing of ingested dsRNA into 21 nt siRNAs in eRNAi-sensitive insects under natural feeding conditions. In addition to dsRNA stability in gut lumen and uptake, dosage of siRNAs targeting a given insect transcript is likely an important factor in order to achieve measurable eRNAi-based regulation in eRNAi-competent insects that lack an apparent silencing amplification mechanism.

No evidence of geographical structure of salicinoid chemotypes within Populus tremula

Salicinoids are well-known defense compounds in salicaceous trees and careful screening at the population level is warranted to fully understand their diversity and function. European aspen, Populus tremula, is a foundation species in Eurasia and highly polymorphic in Sweden. We exhaustively surveyed 102 replicated genotypes from the Swedish Aspen collection (SwAsp) for foliar salicinoids using UHPLC-ESI-TOF/MS and identified nine novel compounds, bringing the total to 19 for this species. Salicinoid structure followed a modular architecture of a salicin skeleton with added side groups, alone or in combination. Two main moieties, 2'-cinnamoyl and 2'-acetyl, grouped the SwAsp population into four distinct chemotypes, and the relative allocation of salicinoids was remarkably constant between different environments, implying a highly channeled biosynthesis of these compounds. Slightly more than half of the SwAsp genotypes belonged to the cinnamoyl chemotype. A fraction synthesized the acetyl moiety alone (∼7%) or in combination with cinnamoyl (∼2%), and close to forty percent lacked either of the two characteristic moieties, and thus resemble P. tremuloides in their salicinoid profile. The two most abundant chemotypes were evenly distributed throughout Sweden, unlike geographical patterns reported for SwAsp phenology traits, plant defense genes, and herbivore community associations. Here we present the salicinoid characterization of the SwAsp collection as a resource for future studies of aspen chemical ecology, salicinoid biosynthesis, and genetics.

Herbivore body condition response in altered environments: mule deer and habitat management

The relationships between habitat, body condition, life history characteristics, and fitness components of ungulates are interwoven and of interest to researchers as they strive to understand the impacts of a changing environment. With the increased availability of portable ultrasound machines and the refinement of hormonal assays, assessment of ungulate body condition has become an accessible monitoring strategy. We employed body condition scoring, estimation of % ingesta-free body fat (%IFBF), assessment of free thyroid hormones (FT4 and FT3), and assessment of pregnancy, as metrics to determine if landscape-level habitat treatments affected body condition of adult (≥1.5 years old) female mule deer (Odocoileus hemionus). All body condition related metrics were measured on 2 neighboring study areas — a reference area that had received no habitat treatments and a treatment study area that had received mechanical removal of pinyon pine (Pinyus edulis) - Utah juniper (Juniperus osteosperma) forest, chemical control of weeds, and reseeding with preferred mule deer browse species. A consistent trend of higher %IFBF was observed in the treatment study area than in the reference study area , although variation of estimates was larger than hypothesized. A similar pattern was observed with higher thyroid hormones concentrations being observed in the treatment study area, but large amounts of variation within concentration estimates were also observed. The consistent pattern of higher body condition related estimates in our treatment study area provides evidence that large mammalian species are sensitive to landscape change, although variation within estimates underlie the challenge in detecting population level impacts stemming from environmental change.

Horizontal gene transfer and functional diversification of plant cell wall degrading polygalacturonases: Key events in the evolution of herbivory in beetles

Plant cell walls are the largest reservoir of organic carbon on earth. To breach and utilize this carbohydrate-rich protective barrier, microbes secrete plant cell wall degrading enzymes (PCWDEs) targeting pectin, cellulose and hemicelluloses. There is a growing body of evidence that genomes of some herbivorous insects also encode PCWDEs, raising questions about their evolutionary origins and functions. Among herbivorous beetles, pectin-degrading polygalacturonases (PGs) are found in the diverse superfamilies Chrysomeloidea (leaf beetles, long-horn beetles) and Curculionoidea (weevils). Here our aim was to test whether these arose from a common ancestor of beetles or via horizontal gene transfer (HGT), and whether PGs kept their ancestral function in degrading pectin or evolved novel functions. Transcriptome data derived from 10 beetle species were screened for PG-encoding sequences and used for phylogenetic comparisons with their bacterial, fungal and plant counterparts. These analyses revealed a large family of PG-encoding genes of Chrysomeloidea and Curculionoidea sharing a common ancestor, most similar to PG genes of ascomycete fungi. In addition, 50 PGs from beetle digestive systems were heterologously expressed and functionally characterized, showing a set of lineage-specific consecutively pectin-degrading enzymes, as well as conserved but enzymatically inactive PG proteins. The evidence indicates that a PG gene was horizontally transferred ∼200 million years ago from an ascomycete fungus to a common ancestor of Chrysomeloidea and Curculionoidea. This has been followed by independent duplications in these two lineages, as well as independent replacement in two sublineages of Chrysomeloidea by two other subsequent HGTs. This origin, leading to subsequent functional diversification of the PG gene family within its new hosts, was a key event promoting the evolution of herbivory in these beetles.

Sex- and tissue-specific profiles of chemosensory gene expression in a herbivorous gall-inducing fly (Diptera: Cecidomyiidae)

BACKGROUND

The chemical senses of insects mediate behaviors that are closely linked to survival and reproduction. The order Diptera contains two model organisms, the vinegar fly Drosophila melanogaster and the mosquito Anopheles gambiae, whose chemosensory genes have been extensively studied. Representing a third dipteran lineage with an interesting phylogenetic position, and being ecologically distinct by feeding on plants, the Hessian fly (Mayetiola destructor Say, Diptera: Cecidomyiidae) genome sequence has recently become available. Among plant-feeding insects, the Hessian fly is unusual in 'reprogramming' the plant to create a superior food and in being the target of plant resistance genes, a feature shared by plant pathogens. Chemoreception is essential for reproductive success, including detection of sex pheromone and plant-produced chemicals by males and females, respectively.

RESULTS

We identified genes encoding 122 odorant receptors (OR), 28 gustatory receptors (GR), 39 ionotropic receptors (IR), 32 odorant binding proteins, and 7 sensory neuron membrane proteins in the Hessian fly genome. We then mapped Illumina-sequenced transcriptome reads to the genome to explore gene expression in male and female antennae and terminal abdominal segments. Our results reveal that a large number of chemosensory genes have up-regulated expression in the antennae, and the expression is in many cases sex-specific. Sex-specific expression is particularly evident among the Or genes, consistent with the sex-divergent olfactory-mediated behaviors of the adults. In addition, the large number of Ors in the genome but the reduced set of Grs and divergent Irs suggest that the short-lived adults rely more on long-range olfaction than on short-range gustation. We also report up-regulated expression of some genes from all chemosensory gene families in the terminal segments of the abdomen, which play important roles in reproduction.

CONCLUSIONS

We show that a large number of the chemosensory genes in the Hessian fly genome have sex- and tissue-specific expression profiles. Our findings provide the first insights into the molecular basis of chemoreception in plant-feeding flies, representing an important advance toward a more complete understanding of olfaction in Diptera and its links to ecological specialization.

Differential activation of sporamin expression in response to abiotic mechanical wounding and biotic herbivore attack in the sweet potato

BACKGROUND

Plants respond differently to mechanical wounding and herbivore attack, using distinct pathways for defense. The versatile sweet potato sporamin possesses multiple biological functions in response to stress. However, the regulation of sporamin gene expression that is activated upon mechanical damage or herbivore attack has not been well studied.

RESULTS

Biochemical analysis revealed that different patterns of Reactive oxygen species (ROS) and antioxidant mechanism exist between mechanical wounding (MW) and herbivore attack (HA) in the sweet potato leaf. Using LC-ESI-MS (Liquid chromatography electrospray ionization mass spectrometry analysis), only the endogenous JA (jasmonic acid) level was found to increase dramatically after MW in a time-dependent manner, whereas both endogenous JA and SA (salicylic acid) increase in parallel after HA. Through yeast one-hybrid screening, two transcription factors IbNAC1 (no apical meristem (NAM), Arabidopsis transcription activation factor (ATAF), and cup-shaped cotyledon (CUC)) and IbWRKY1 were isolated, which interact with the sporamin promoter fragment of SWRE (sporamin wounding-responsive element) regulatory sequences. Exogenous application of MeJA (methyl jasmonate), SA and DIECA (diethyldithiocarbamic acid, JAs biosynthesis inhibitor) on sweet potato leaves was employed, and the results revealed that IbNAC1 mediated the expression of sporamin through a JA-dependent signaling pathway upon MW, whereas both IbNAC1 and IbWRKY1 coordinately regulated sporamin expression through JA- and SA-dependent pathways upon HA. Transcriptome analysis identified MYC2/4 and JAZ2/TIFY10A (jasmonate ZIM/tify-domain), the repressor and activator of JA and SA signaling among others, as the genes that play an intermediate role in the JA and SA pathways, and these results were further validated by qRT-PCR (quantitative real-time polymerase chain reaction).

CONCLUSION

This work has improved our understanding of the differential regulatory mechanism of sporamin expression. Our study illustrates that sweet potato sporamin expression is differentially induced upon abiotic MW and biotic HA that involves IbNAC1 and IbWRKY1 and is dependent on the JA and SA signaling pathways. Thus, we established a model to address the plant-wounding response upon physical and biotic damage.

Plant resistance to aphid feeding: behavioral, physiological, genetic and molecular cues regulate aphid host selection and feeding

Aphids damage major world food and fiber crops through direct feeding and transmission of plant viruses. Fortunately, the development of many aphid-resistant crop plants has provided both ecological and economic benefits to food production. Plant characters governing aphid host selection often dictate eventual plant resistance or susceptibility to aphid herbivory, and these phenotypic characters have been successfully used to map aphid resistance genes. Aphid resistance is often inherited as a dominant trait, but is also polygenic and inherited as recessive or incompletely dominant traits. Most aphid-resistant cultivars exhibit constitutively expressed defenses, but some cultivars exhibit dramatic aphid-induced responses, resulting in the overexpression of large ensembles of putative aphid resistance genes. Two aphid resistance genes have been cloned. Mi-1.2, an NBS-LRR gene from wild tomato, confers resistance to potato aphid and three Meloidogyne root-knot nematode species, and Vat, an NBS-LRR gene from melon, controls resistance to the cotton/melon aphid and to some viruses. Virulence to aphid resistance genes of plants occurs in 17 aphid species--more than half of all arthropod biotypes demonstrating virulence. The continual appearance of aphid virulence underscores the need to identify new sources of resistance of diverse sequence and function in order to delay or prevent biotype development.

Genome-wide analysis of Nilaparvata lugens nymphal responses to high-density and low-quality rice hosts

The brown planthopper (BPH) Nilaparvata lugens is an economically important pest on rice plants. In this study, the higher population density and yellow-ripe stage of rice plants were used to construct adverse survival conditions (ASC) against BPH nymphs. Simultaneously, the low population density and tillering stage of rice plants were used to establish a suitable survival condition (SSC) as a control. Solexa/Illumina sequencing was used to identify genes of BPH nymphs responding to ASC. Significantly longer duration development of BPH nymphs and significantly lower brachypterous ratio of BPH adults were observed by ASC compared with SSC. A total of 2 544 differentially expressed genes (DEGs) were obtained and analyzed by BLASTx, Gene Ontology and KEGG Orthology. Gene ontology analysis revealed that the DEGs were mainly involved in categories of cell, cell part, cellular process, binding, catalytic, organelle and metabolic processes. 1 138 DEGs having enzyme commission numbers were assigned to different metabolic pathways. The largest clusters were neurodegenerative diseases (137, 12.0%), followed by carbohydrate metabolism (113, 9.9%), amino acid metabolism (94, 8.3%), nucleotide metabolism (76, 6.7%), energy metabolism (64, 5.6%), translation (60, 5.3%), lipid metabolism (58, 5.1%), and folding, sorting and degradation (52, 4.6%). Expressing profile of 11 DEGs during eight nymphal developmental stages of BPH were analyzed by quantitative real-time polymerase chain reaction. The 11 genes exhibited differential expression between ASC and SSC during at least one developmental stage. The DEGs identified in this study provide molecular proof of how BPH reconfigures its gene expression profile to adapt to overcrowding and low-quality hosts.

Testing local host adaptation and phenotypic plasticity in a herbivore when alternative related host plants occur sympatrically

Host race formation in phytophagous insects can be an early stage of adaptive speciation. However, the evolution of phenotypic plasticity in host use is another possible outcome. Using a reciprocal transplant experiment we tested the hypothesis of local adaptation in the aphid Brevicoryne brassicae. Aphid genotypes derived from two sympatric host plants, Brassica oleracea and B. campestris, were assessed in order to measure the extent of phenotypic plasticity in morphological and life history traits in relation to the host plants. We obtained an index of phenotypic plasticity for each genotype. Morphological variation of aphids was summarized by principal components analysis. Significant effects of recipient host on morphological variation and life history traits (establishment, age at first reproduction, number of nymphs, and intrinsic growth rate) were detected. We did not detected genotype × host plant interaction; in general the genotypes developed better on B. campestris, independent of the host plant species from which they were collected. Therefore, there was no evidence to suggest local adaptation. Regarding plasticity, significant differences among genotypes in the index of plasticity were detected. Furthermore, significant selection on PC1 (general aphid body size) on B. campestris, and on PC1 and PC2 (body length relative to body size) on B. oleracea was detected. The elevation of the reaction norm of PC1 and the slope of the reaction norm for PC2 (i.e., plasticity) were under directional selection. Thus, host plant species constitute distinct selective environments for B. brassicae. Aphid genotypes expressed different phenotypes in response to the host plant with low or nil fitness costs. Phenotypic plasticity and gene flow limits natural selection for host specialization promoting the maintenance of genetic variation in host exploitation.

Overcompensation in response to herbivory in Arabidopsis thaliana: the role of glucose-6-phosphate dehydrogenase and the oxidative pentose-phosphate pathway

That some plants benefit from being eaten is counterintuitive, yet there is now considerable evidence demonstrating enhanced fitness following herbivory (i.e., plants can overcompensate). Although there is evidence that genetic variation for compensation exists, little is known about the genetic mechanisms leading to enhanced growth and reproduction following herbivory. We took advantage of the compensatory variation in recombinant inbred lines of Arabidopsis thaliana, combined with microarray and QTL analyses to assess the molecular basis of overcompensation. We found three QTL explaining 11.4, 10.1, and 26.7% of the variation in fitness compensation, respectively, and 109 differentially expressed genes between clipped and unclipped plants of the overcompensating ecotype Columbia. From the QTL/microarray screen we uncovered one gene that plays a significant role in overcompensation: glucose-6-phosphate-1-dehydrogenase (G6PDH1). Knockout studies of Transfer-DNA (T-DNA) insertion lines and complementation studies of G6PDH1 verify its role in compensation. G6PDH1 is a key enzyme in the oxidative pentose-phosphate pathway that plays a central role in plant metabolism. We propose that plants capable of overcompensating reprogram their transcriptional activity by up-regulating defensive genes and genes involved in energy metabolism and by increasing DNA content (via endoreduplication) with the increase in DNA content feeding back on pathways involved in defense and metabolism through increased gene expression.

Forward genetics by genome sequencing reveals that rapid cyanide release deters insect herbivory of Sorghum bicolor

Whole genome sequencing has allowed rapid progress in the application of forward genetics in model species. In this study, we demonstrated an application of next-generation sequencing for forward genetics in a complex crop genome. We sequenced an ethyl methanesulfonate-induced mutant of Sorghum bicolor defective in hydrogen cyanide release and identified the causal mutation. A workflow identified the causal polymorphism relative to the reference BTx623 genome by integrating data from single nucleotide polymorphism identification, prior information about candidate gene(s) implicated in cyanogenesis, mutation spectra, and polymorphisms likely to affect phenotypic changes. A point mutation resulting in a premature stop codon in the coding sequence of dhurrinase2, which encodes a protein involved in the dhurrin catabolic pathway, was responsible for the acyanogenic phenotype. Cyanogenic glucosides are not cyanogenic compounds but their cyanohydrins derivatives do release cyanide. The mutant accumulated the glucoside, dhurrin, but failed to efficiently release cyanide upon tissue disruption. Thus, we tested the effects of cyanide release on insect herbivory in a genetic background in which accumulation of cyanogenic glucoside is unchanged. Insect preference choice experiments and herbivory measurements demonstrate a deterrent effect of cyanide release capacity, even in the presence of wild-type levels of cyanogenic glucoside accumulation. Our gene cloning method substantiates the value of (1) a sequenced genome, (2) a strongly penetrant and easily measurable phenotype, and (3) a workflow to pinpoint a causal mutation in crop genomes and accelerate in the discovery of gene function in the postgenomic era.

Wide host ranges of herbivorous beetles? Insights from DNA bar coding

There are very few studies that have investigated host-specificity among tropical herbivorous insects. Indeed, most of the trophic interactions of herbivorous insects in Southeast Asian tropical rainforests remain unknown, and whether polyphagous feeding is common in the herbivores of this ecosystem has not been determined. The present study employed DNA bar coding to reveal the trophic associations of adult leaf-chewing chrysomelid beetles in a Bornean rainforest. Plant material ingested by the adults was retrieved from the bodies of the insects, and a portion of the chloroplast rbcL sequence was then amplified from this material. The plants were identified at the family level using an existing reference database of chloroplast DNA. Our DNA-based diet analysis of eleven chrysomelid species successfully identified their host plant families and indicated that five beetle species fed on more than two families within the angiosperms, and four species fed on several families of gymnosperms and/or ferns together with multiple angiosperm families. These findings suggest that generalist chrysomelid beetles associated with ecologically and taxonomically distant plants constitute a part of the plant-insect network of the Bornean rainforest.

A porin-like protein from oral secretions of Spodoptera littoralis larvae induces defense-related early events in plant leaves

Insect herbivory on plants is a complex incident consisting of at least two different aspects, namely mechanical damage and chemical challenge, as feeding insects introduce oral secretions (OS) into the wounded tissue of the attacked plant. Mechanical wounding alone is sufficient to induce a set of defense-related reactions in host plants, but some early events such as membrane potential (Vm) changes and cytosolic Ca²⁺-elevations can be triggered only by herbivores suggesting that OS-derived molecules are involved in those processes. Following an assay-guided purification based on planar lipid bilayer membrane technique in combination with proteomic analysis, a porin-like protein (PLP) of most likely bacterial origin was determined from collected OS of Spodoptera littoralis larvae. PLP exhibited channel-forming activity. Further, early defense-related events in plant-insect interaction were evaluated by using a purified fraction and α-hemolysin (α-HL) as a commercial pore-forming compound. Both up-regulated the calmodulin-like CML42 in Arabidopsis thaliana, which only responds to oral secretion and not to wounding. An elevation of in vivo [Ca²⁺](cyt) was not observed. Because membrane channel formation is a widespread phenomenon in plant-insect interactions, this PLP might represent an example for microbial compounds from the insect gut which are initially involved in plant-insect interactions.

Understanding plant defence responses against herbivore attacks: an essential first step towards the development of sustainable resistance against pests

Plant-herbivore relationships are complex interactions encompassing elaborate networks of molecules, signals and strategies used to overcome defences developed by each other. Herbivores use multiple feeding strategies to obtain nutrients from host plants. In turn, plants respond by triggering defence mechanisms to inhibit, block or modify the metabolism of the pest. As part of these defences, herbivore-challenged plants emit volatiles to attract natural enemies and warn neighbouring plants of the imminent threat. In response, herbivores develop a variety of strategies to suppress plant-induced protection. Our understanding of the plant-herbivore interphase is limited, although recent molecular approaches have revealed the participation of a battery of genes, proteins and volatile metabolites in attack-defence processes. This review describes the intricate and dynamic defence systems governing plant-herbivore interactions by examining the diverse strategies plants employ to deny phytophagous arthropods the ability to breach newly developed mechanisms of plant resistance. A cornerstone of this understanding is the use of transgenic tools to unravel the complex networks that control these interactions.

Comparative genomic analysis of the microbiome [corrected] of herbivorous insects reveals eco-environmental adaptations: biotechnology applications

Metagenome analysis of the gut symbionts of three different insects was conducted as a means of comparing taxonomic and metabolic diversity of gut microbiomes to diet and life history of the insect hosts. A second goal was the discovery of novel biocatalysts for biorefinery applications. Grasshopper and cutworm gut symbionts were sequenced and compared with the previously identified metagenome of termite gut microbiota. These insect hosts represent three different insect orders and specialize on different food types. The comparative analysis revealed dramatic differences among the three insect species in the abundance and taxonomic composition of the symbiont populations present in the gut. The composition and abundance of symbionts was correlated with their previously identified capacity to degrade and utilize the different types of food consumed by their hosts. The metabolic reconstruction revealed that the gut metabolome of cutworms and grasshoppers was more enriched for genes involved in carbohydrate metabolism and transport than wood-feeding termite, whereas the termite gut metabolome was enriched for glycosyl hydrolase (GH) enzymes relevant to lignocellulosic biomass degradation. Moreover, termite gut metabolome was more enriched with nitrogen fixation genes than those of grasshopper and cutworm gut, presumably due to the termite's adaptation to the high fiber and less nutritious food types. In order to evaluate and exploit the insect symbionts for biotechnology applications, we cloned and further characterized four biomass-degrading enzymes including one endoglucanase and one xylanase from both the grasshopper and cutworm gut symbionts. The results indicated that the grasshopper symbiont enzymes were generally more efficient in biomass degradation than the homologous enzymes from cutworm symbionts. Together, these results demonstrated a correlation between the composition and putative metabolic functionality of the gut microbiome and host diet, and suggested that this relationship could be exploited for the discovery of symbionts and biocatalysts useful for biorefinery applications.