Coevolution can explain defensive secondary metabolite diversity in plants

How does plant chemodiversity evolve? Testing five hypotheses in one population genetic model

Plant chemodiversity, the diversity of plant-specialized metabolites, is an important dimension of biodiversity. However, there are so far few mathematical models to test verbal hypotheses on how chemodiversity evolved. Here, we develop such a model to test predictions of five hypotheses: the 'fluctuating selection hypothesis', the 'dominance reversal hypothesis', the interaction diversity hypothesis, the synergy hypothesis, and the screening hypothesis. We build a population genetic model of a plant population attacked by herbivore species whose occurrence fluctuates over time. We study the model using mathematical analysis and individual-based simulations. As predicted by the 'dominance reversal hypothesis', chemodiversity can be maintained if alleles conferring a defense metabolite are dominant with respect to the benefits, but recessive with respect to costs. However, even smaller changes in dominance can maintain polymorphism. Moreover, our results underpin and elaborate predictions of the synergy and interaction diversity hypotheses, and, to the extent that our model can address it, the screening hypotheses. By contrast, we found only partial support for the 'fluctuating selection hypothesis'. In summary, we have developed a flexible model and tested various verbal models for the evolution of chemodiversity. Next, more mechanistic models are needed that explicitly consider the organization of metabolic pathways.

Global patterns and determinants of the initial concentrations of litter carbon components

Plant litter is an important carbon (C) and nutrient pool in terrestrial ecosystems. The C components in plant litter are important because they regulate plant litter decomposition rate, but little is known on the global patterns and determinants of their concentrations in freshly fallen plant litter. Here, we quantified the concentrations of leaf litter C components (i.e., carbohydrate, polyphenol, tannin, and condensed tannin) with 864 measurements from 161 independent publications. We found that (1) the mean concentrations of leaf litter carbohydrate, polyphenol, tannin and condensed tannin were 27.7, 6.08, 8.84 and 5.7 %, respectively; (2) the concentrations of leaf litter C components were affected by taxonomic division, mycorrhizal association, life form, and/or leaf shedding strategy; (3) soil property had similar impacts on the concentrations of the four C compounds, while the influence of mean annual temperature and precipitation varied; and (4) elevation had opposing effects on carbohydrate and polyphenol concentrations, but not on that of tannin and condensed tannin, and only carbohydrate concentration was strongly affected by absolute latitude. In general, our results clearly show the global patterns and drivers of the concentrations of litter C compounds, providing new insights into the role of litter decomposition in global C dynamics.

Nitrogen-mediated volatilisation of defensive metabolites in tomato confers resistance to herbivores

Plants synthesise a vast array of volatile organic compounds (VOCs), which serve as chemical defence and communication agents in their interactions with insect herbivores. Although nitrogen (N) is a critical resource in the production of plant metabolites, its regulatory effects on defensive VOCs remain largely unknown. Here, we investigated the effect of N content in tomato (Solanum lycopersicum) on the tobacco cutworm (Spodoptera litura), a notorious agricultural pest, using biochemical and molecular experiments in combination with insect behavioural and performance analyses. We observed that on tomato leaves with different N contents, S. litura showed distinct feeding preference and growth and developmental performance. Particularly, metabolomics profiling revealed that limited N availability conferred resistance upon tomato plants to S. litura is likely associated with the biosynthesis and emission of the volatile metabolite α-humulene as a repellent. Moreover, exogenous application of α-humulene on tomato leaves elicited a significant repellent response against herbivores. Thus, our findings unravel the key factors involved in N-mediated plant defence against insect herbivores and pave the way for innovation of N management to improve the plant defence responses to facilitate pest control strategies within agroecosystems.

Genomic insights into the evolution of secondary metabolism of Escovopsis and its allies, specialized fungal symbionts of fungus-farming ants

The metabolic intimacy of symbiosis often demands the work of specialists. Natural products and defensive secondary metabolites can drive specificity by ensuring infection and propagation across host generations. But in contrast to bacteria, little is known about the diversity and distribution of natural product biosynthetic pathways among fungi and how they evolve to facilitate symbiosis and adaptation to their host environment. In this study, we define the secondary metabolism of Escovopsis and closely related genera, symbionts in the gardens of fungus-farming ants. We ask how the gain and loss of various biosynthetic pathways correspond to divergent lifestyles. Long-read sequencing allowed us to define the chromosomal features of representative Escovopsis strains, revealing highly reduced genomes composed of seven to eight chromosomes. The genomes are highly syntenic with macrosynteny decreasing with increasing phylogenetic distance, while maintaining a high degree of mesosynteny. An ancestral state reconstruction analysis of biosynthetic pathways revealed that, while many secondary metabolites are shared with non-ant-associated Sordariomycetes, 56 pathways are unique to the symbiotic genera. Reflecting adaptation to diverging ant agricultural systems, we observe that the stepwise acquisition of these pathways mirrors the ecological radiations of attine ants and the dynamic recruitment and replacement of their fungal cultivars. As different clades encode characteristic combinations of biosynthetic gene clusters, these delineating profiles provide important insights into the possible mechanisms underlying specificity between these symbionts and their fungal hosts. Collectively, our findings shed light on the evolutionary dynamic nature of secondary metabolism in Escovopsis and its allies, reflecting adaptation of the symbionts to an ancient agricultural system.IMPORTANCEMicrobial symbionts interact with their hosts and competitors through a remarkable array of secondary metabolites and natural products. Here, we highlight the highly streamlined genomic features of attine-associated fungal symbionts. The genomes of Escovopsis species, as well as species from other symbiont genera, many of which are common with the gardens of fungus-growing ants, are defined by seven chromosomes. Despite a high degree of metabolic conservation, we observe some variation in the symbionts' potential to produce secondary metabolites. As the phylogenetic distribution of the encoding biosynthetic gene clusters coincides with attine transitions in agricultural systems, we highlight the likely role of these metabolites in mediating adaptation by a group of highly specialized symbionts.

Adaptation of Ustilago maydis to phenolic and alkaloid responsive metabolites in maize B73

The adaptation of pathogenic fungi to plant-specialized metabolites is necessary for their survival and reproduction. The biotrophic fungus Ustilago maydis can cause maize smut and produce tumors in maize (Zea mays), resulting in reduced maize yield and significant economic losses. Qualitative analysis using UPLC-MS/MS revealed that the infection of maize variety B73 with U. maydis resulted in increased levels of phytohormones, phenolics, and alkaloids in maize seedling tissues. However, correlation analysis showed that nearly all compounds in the mechanical damage group were significantly negatively correlated with the shoot growth indexes of maize B73. The correlation coefficients of 2-hydroxy-7-methoxy-1,4-benzoxazin-3-one (HMBOA) and maize B73 shoot length and shoot weight were r = -0.56 (p < 0.01) and r = -0.75 (p < 0.001), respectively. In the inoculation group, these correlations weakened, with the correlation coefficients between HMBOA and maize B73 shoot length and shoot weight being r = 0.02 and r = -0.1, respectively. The correlation coefficients between 6-methoxy-2-benzoxazolinone (MBOA) and the shoot weight were r = -0.73 (p < 0.001) and r = -0.15 in the mechanical damage group and inoculation group, respectively. These findings suggest that increased concentrations of these compounds are more positively associated with mechanical damage than with U. maydis infection. At high concentrations, most of these compounds had an inhibitory effect on U. maydis. This study investigated the ability of U. maydis to regulate various compounds, including phytohormones, phenolic acids, and alkaloids in maize B73, providing evidence that U. maydis has adapted to the specialized metabolites produced by maize B73.

Oviposition strategies in Pieridae butterflies and the role of an egg-killing plant trait therein

Most herbivorous insects are host-plant specialists that evolved detoxification mechanisms to overcome their host plant's toxins. In the evolutionary arms-races between Pieridae butterflies and Brassicaceae plants, some plant species have evolved another defence against the pierids: egg-killing. Underneath the eggs, leaves develop a so-called hypersensitive response (HR)-like cell death. Whether some butterflies have evolved oviposition strategies to counter-adapt against egg-killing remains to be studied. In this study, we assessed the oviposition site location of Pieridae butterflies on their natural host plants. We described the plant tissue on which we located the eggs of the most common Pieridae in the Netherlands: Gonepteryx rhamni, Anthocharis cardamines, Pieris rapae, P. napi, P. brassicae and P. mannii. Additionally, we assessed expression of HR-like cell death in response to the deposited butterfly eggs. We found that both A. cardamines and G. rhamni mainly oviposited on the floral stem and the branch, respectively, and oviposited on host plants from lineages not expected to kill pierid eggs. Accordingly, no HR-like cell death was seen. All Pieris eggs found were located on leaves of their host, the only tissue found to express HR-like cell death. Furthermore, each Pieris species was found to at least occasionally oviposit on Brassica nigra. This was the only plant species in this survey that expressed HR-like cell death in response to the eggs of P. rapae, P. napi and P. brassicae. Our observations demonstrate that HR-like cell death remains an effective defence strategy against these Pieris species and as such did not find evidence for the hypothesized counterstrategies. Surveying certain key species and disentangling the micro-evolution of oviposition strategies within a species would allow us to further investigate potential counter-adaptations that evolved against HR-like cell death. This study provides the basis for further investigation of potential counter-adaptations to egg-killing defences.

Disarming the defenses: Insect detoxification of plant defense-related specialized metabolites

The ability of certain insects to feed on plants containing toxic specialized metabolites may be attributed to detoxification enzymes. Representatives of a few large families of detoxification enzymes are widespread in insect herbivores acting to functionalize toxins and conjugate them with polar substituents to decrease toxicity, increase water solubility and enhance excretion. Insects have also developed specific enzymes for coping with toxins that are activated upon plant damage. Another source of detoxification potential in insects lies in their microbiomes, which are being increasingly recognized for their role in processing plant toxins. The evolution of insect detoxification systems to resist toxic specialized metabolites in plants may in turn have selected for the great diversity of such metabolites found in nature.

Keap1 Negatively Regulates Transcription of Three Counter-Defense Genes and Susceptibility to Plant Toxin Gossypol in Helicoverpa armigera

Expressions of a wide range of cytoprotective counter-defense genes are mainly regulated by the Keap1-Nrf2-ARE signaling pathway in response to oxidative stress from xenobiotics. Gossypol is the major antiherbivore secondary metabolite of cotton, but how the polyphagous pest Helicoverpa armigera copes with this phytochemical to utilize its favorite host plant cotton remains largely elusive. In this study, we first suppressed the Keap1 gene in newly hatched larvae of cotton bollworm by feeding them the siRNA diet for 4 days. All of the larvae were subsequently fed the artificial diet supplied with gossypol or the control diet for 5 days. We identified that the knockdown of the Keap1 gene significantly decreased larval mortality and significantly increased the percentages of larval survival, reaching the fourth instar, compared with ncsiRNA when exposed to a diet containing gossypol. Three counter-defense genes CYP9A17, CYP4L11 and UGT41B3, which were related to the induction or metabolism of gossypol according to the report before, were all significantly up-regulated after the knockdown of the Keap1 gene. The Antioxidant Response Elements (AREs) were also detected in the promoter regions of the three counter-defense genes above. These data indicate that the suppression of the Keap1 gene activates the Keap1-Nrf2-ARE signaling pathway, up-regulates the expressions of counter-defense genes involved in the resistance of oxidative stress and finally contributes to reducing the susceptibility of gossypol. Our results provide more knowledge about the transcriptional regulation mechanisms of counter-defense genes that enable the cotton bollworm to adapt to the diversity of host plants including cotton.

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

The evolution of chemodiversity in plants-From verbal to quantitative models

Plants harbour a great chemodiversity, that is diversity of specialised metabolites (SMs), at different scales. For instance, individuals can produce a large number of SMs, and populations can differ in their metabolite composition. Given the ecological and economic importance of plant chemodiversity, it is important to understand how it arises and is maintained over evolutionary time. For other dimensions of biodiversity, that is species diversity and genetic diversity, quantitative models play an important role in addressing such questions. Here, we provide a synthesis of existing hypotheses and quantitative models, that is mathematical models and computer simulations, for the evolution of plant chemodiversity. We describe each model's ingredients, that is the biological processes that shape chemodiversity, the scales it considers and whether it has been formalized as a quantitative model. Although we identify several quantitative models, not all are dynamic and many influential models have remained verbal. To fill these gaps, we outline our vision for the future of chemodiversity modelling. We identify quantitative models used for genetic variation that may be adapted for chemodiversity, and we present a flexible framework for the creation of individual-based models that address different scales of chemodiversity and combine different ingredients that bring this chemodiversity about.

Mixtures of Milkweed Cardenolides Protect Monarch Butterflies against Parasites

Plants have evolved a diverse arsenal of defensive secondary metabolites in their evolutionary arms race with insect herbivores. In addition to the bottom-up forces created by plant chemicals, herbivores face top-down pressure from natural enemies, such as predators, parasitoids and parasites. This has led to the evolution of specialist herbivores that do not only tolerate plant secondary metabolites but even use them to fight natural enemies. Monarch butterflies (Danaus plexippus) are known for their use of milkweed chemicals (cardenolides) as protection against vertebrate predators. Recent studies have shown that milkweeds with high cardenolide concentrations can also provide protection against a virulent protozoan parasite. However, whether cardenolides are directly responsible for these effects, and whether individual cardenolides or mixtures of these chemicals are needed to reduce infection, remains unknown. We fed monarch larvae the four most abundant cardenolides found in the anti-parasitic-milkweed Asclepias curassavica at varying concentrations and compositions to determine which provided the highest resistance to parasite infection. Measuring infection rates and infection intensities, we found that resistance is dependent on both concentration and composition of cardenolides, with mixtures of cardenolides performing significantly better than individual compounds, even when mixtures included lower concentrations of individual compounds. These results suggest that cardenolides function synergistically to provide resistance against parasite infection and help explain why only milkweed species that produce diverse cardenolide compounds provide measurable parasite resistance. More broadly, our results suggest that herbivores can benefit from consuming plants with diverse defensive chemical compounds through release from parasitism.

Spatially resolved detection of small molecules from press-dried plant tissue using MALDI imaging

Premise

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) is a chemical imaging method that can visualize spatial distributions of particular molecules. Plant tissue imaging has so far mostly used cryosectioning, which can be impractical for the preparation of large-area imaging samples, such as full flower petals. Imaging unsectioned plant tissue presents its own difficulties in extracting metabolites to the surface due to the waxy cuticle.

Methods

We address this by using established delipidation techniques combined with a solvent vapor extraction prior to applying the matrix with many low-concentration sprays.

Results

Using this procedure, we imaged tissue from three different plant species (two flowers and one carnivorous plant leaf). Material factorization analysis of the resulting data reveals a wide range of plant-specific small molecules with varying degrees of localization to specific portions of the tissue samples, while facilitating detection and removal of signal from background sources.

Conclusions

This work demonstrates applicability of MALDI-MSI to press-dried plant samples without freezing or cryosectioning, setting the stage for spatially resolved molecule identification. Increased mass resolution and inclusion of tandem mass spectrometry are necessary next steps to allow more specific and reliable compound identification.

Isolation and Identification of Herbicidal Active Compounds from Brassica oleracea L. and Exploration of the Binding Sites of Brassicanate A Sulfoxide

Brassica oleracea L. has strong allelopathic effects on weeds. However, the allelochemicals with herbicidal activity in B. oleracea L. are still unknown. In this study, we evaluated the activity of allelochemicals isolated from Brassica oleracea L. based on the germination and growth of model plant Lactuca sativa Linn., grass weed Panicum miliaceum, and broadleaf weed Chenopodium album. Additionally, we employed molecular docking to predict the binding of brassicanate A sulfoxide to herbicide targets. The results of this study showed that eight compounds with herbicidal activity were isolated from B. oleracea L., and the predicted results indicated that brassicanate A sulfoxide was stably bound to dihydroxyacid dehydratase, hydroxymethylpyruvate dioxygenase, acetolactate synthase, PYL family proteins and transport inhibitor response 1. This research provides compound sources and a theoretical foundation for the development of natural herbicides.

Caterpillars Detoxify Diterpenoid from Nepeta stewartiana by the Molting Hormone Gene CYP306A1

Herbivorous insects are well known for detoxifying a broad range of the defense compounds produced by the plants that they feed on, but knowledge of the mechanisms of detoxification is still very limited. Here, we describe a system in which two species of lepidopteran caterpillars metabolize an abietane diterpene from the plants of Nepeta stewartiana Diels to an oxygenated derivative that is less active biologically. We found that this transformation could be catalyzed by a cytochrome P450 enzyme in caterpillars, which are associated with molting. Most interestingly, abietane diterpene targets the molting-associated gene CYP306A1 to alter the content of molting hormones in the insect at specific developmental stages and competitively inhibit molting hormone metabolism. These findings identify the mechanism by which caterpillars are able to detoxify abietane diterpenoid through hydroxylation at the C-19 position, which may be opening up exciting research questions into the mechanisms of interaction between plants and insects.

Potentiating Biosynthesis of Alkaloids and Polyphenolic Substances in Catharanthus roseus Plant Using ĸ-Carrageenan

Catharanthus roseus is a medicinal plant that produces indole alkaloids, which are utilized in anticancer therapy. Vinblastine and vincristine, two commercially important antineoplastic alkaloids, are mostly found in the leaves of Catharanthus roseus. ĸ-carrageenan has been proven as plant growth promoting substance for a number of medicinal and agricultural plants. Considering the importance of ĸ-carrageenan as a promoter of plant growth and phytochemical constituents, especially alkaloids production in Catharanthus roseus, an experiment was carried out to explore the effect of ĸ-carrageenan on the plant growth, phytochemicals content, pigments content, and production of antitumor alkaloids in Catharanthus roseus after planting. Foliar application of ĸ-carrageenan (at 0, 400, 600 and 800 ppm) significantly improved the performance of Catharanthus roseus. Phytochemical analysis involved determining the amount of total phenolics (TP), flavonoids (F), free amino acids (FAA), alkaloids (TAC) and pigments contents by spectrophotometer, minerals by ICP, amino acids, phenolic compounds and alkaloids (Vincamine, Catharanthine, Vincracine (Vincristine), and vinblastine) analysis uses HPLC. The results indicated that all examined ĸ-carrageenan treatments led to a significant (p ≤ 0.05) increase in growth parameters compared to the untreated plants. Phytochemical examination indicates that the spray of ĸ-carrageenan at 800 mg L^-1 increased the yield of alkaloids (Vincamine, Catharanthine and Vincracine (Vincristine)) by 41.85 μg/g DW, total phenolic compounds by 3948.6 μg gallic/g FW, the content of flavonoids 951.3 μg quercetin /g FW and carotenoids content 32.97 mg/g FW as compared to the control. An amount of 400 ppm ĸ-carrageenan treatment gave the best contents of FAA, Chl a, Chl b and anthocyanin. The element content of K, Ca, Cu, Zn and Se increased by treatments. Amino acids constituents and phenolics compounds contents were altered by ĸ-carrageenan.

A Review of the Host Plant Location and Recognition Mechanisms of Asian Longhorn Beetle

The Asian longhorn beetle (ALB), Anoplophora glabripennis Motschulsky, is a polyphagous xylophage with dozens of reported host tree species. However, the mechanisms by which individuals locate and recognize host plants are still unknown. We summarize the current knowledge of the host plant list, host kairomones, odorant-binding proteins (OBPs) and microbial symbionts of this beetle and their practical applications, and finally discuss the host localization and recognition mechanisms. A total of 209 species (or cultivars) were reported as ALB host plants, including 101 species of higher sensitivity; host kairomones were preferentially bound to ALB recombinant OBPs, including cis-3-hexen-1-ol, δ-3-carene, nonanal, linalool, and β-caryophyllene. In addition, microbial symbionts may help ALB degrade their host. Complementarity of tree species with different levels of resistance may reduce damage, but trapping effectiveness for adults was limited using a combination of host kairomones and sex pheromones in the field. Therefore, we discuss host location behavior from a new perspective and show that multiple cues are used by ALB to locate and recognize host plants. Further research into host resistance mechanisms and visual signal recognition, and the interaction of sex pheromone synthesis, symbiont microbiota, and host plants may help reveal the host recognition mechanisms of ALBs.

Inducible Gut-Specific Carboxylesterase SlCOE030 in Polyphagous Pests of Spodoptera litura Conferring Tolerance between Nicotine and Cyantraniliprole

Insecticides tolerance in herbivorous arthropods is associated with preadaptation to host plant allelochemicals. However, how plant secondary metabolites activate detoxifying metabolic genes to develop tolerance remains unclear. Herein, the tolerance of Spodoptera litura larvae to cyantraniliprole was increased after nicotine exposure. An S. litura α esterase, SlCOE030, was predominantly expressed in the midgut and induced after exposure to cyantraniliprole, nicotine, and cyantraniliprole plus nicotine. Drosophila melanogaster with ectopically overexpressed SlCOE030 enhanced cyantraniliprole and nicotine tolerance by 4.91- and 2.12-fold, respectively. Compared to UAS-SlCOE030 and Esg-GAL4 lines, the Esg > SlCOE030 line laid more eggs after nicotine exposure. SlCOE030 knockdown decreased the sensitivity of nicotine-treated S. litura larvae to cyantraniliprole. Metabolism assays indicated that recombinant SlCOE030 protein metabolizes cyantraniliprole. Homology modeling and molecular docking analysis demonstrated that SlCOE030 exhibits effective affinities for cyantraniliprole and nicotine. Thus, insect CarEs may result in the development of cross-tolerance between synthetic insecticides and plant secondary metabolites.

Defensive Resistance of Cowpea Vigna unguiculata Control Megalurothrips usitatus Mediated by Jasmonic Acid or Insect Damage

Vigna unguiculata is a vital vegetable crop in Southeast Asia, and Megalurothrips usitatus can cause huge damage to this crop. Enhancing the resistance of V. unguiculata against M. usitatus is a promising way to protect this crop; however, there is limited information regarding the mechanism underlying the resistance of V. unguiculata against M. usitatus. Here, a behavior assay was performed to explore the resistance of V. unguiculata against M. usitatus after insect damage or treatment by jasmonic acid (JA). Furthermore, transcriptome and metabonomics analysis was used to detect the putative mechanism underlying the resistance of V. unguiculata against M. usitatus. The pre-treatment of Vigna unguiculata with JA or infestation with Megalurothrips usitatus alleviated the damage resulting from the pest insect. We further identified differentially expressed genes and different metabolites involved in flavonoid biosynthesis and alpha-linolenic acid metabolism. Genes of chalcone reductase and shikimate O-hydroxycinnamoyltransferase involved in flavonoid biosynthesis, as well as lipoxygenase and acyl-CoA oxidase involved in alpha-linolenic acid metabolism, were upregulated in plants after herbivory or JA supplementation. The upregulation of these genes contributed to the high accumulation of metabolites involved in flavonoid biosynthesis and the alpha-linolenic acid metabolism pathway. These transcriptional and metabolite changes are potentially responsible for plant defense and a putative regulatory model is thus proposed to illustrate the cowpea defense mechanism against insect attack. Our study provides candidate targets for the breeding of varieties with resistance to insect herbivory by molecular technology.

The Distribution of Phytoecdysteroids among Terrestrial Vascular Plants: A Comparison of Two Databases and Discussion of the Implications for Plant/Insect Interactions and Plant Protection

Phytoecdysteroids are a class of plant secondary compounds which are present in a wide diversity of vascular plant species, where they contribute to a reduction in invertebrate predation. Over the past 55 years, a significant body of heterogeneous literature on the presence, identities and/or quantities of ecdysteroids in plant species has accumulated, resulting in the compilation of a first database, the Ecdybase Literature Survey (ELS; 4908 entries, covering 2842 species). A second extensive database on the distribution of ecdysteroids in vascular plants is available as the Exeter Survey (ES; 4540 entries, covering 4155 species), which used standardised extraction and analysis methods to survey seeds/spores. We compare the usefulness of these two databases to provide information on the occurrence of phytoecdysteroids at the order/family levels in relation to the recent molecular classifications of gymnosperms, pteridophytes/lycophytes and angiosperms. The study, in conjunction with the other published literature, provides insights into the distribution of phytoecdysteroids in the plant world, their role in plant protection in nature and their potential future contribution to crop protection. Furthermore, it will assist future investigations in the chemotaxonomy of phytoecdysteroids and other classes of plant secondary compounds.

Phytochemical and Biological Study of Trophic Interaction between Pseudosphinx Tetrio L. Larvae and Allamanda Cathartica L

In this article, we propose to explore the chemical interaction between Pseudosphinx tetrio L. and Allamanda cathartica L. using different analytical methods, including an innovative electrochemical approach (called electrochemical ecology) and multivariate analysis, and we investigate the potential antimicrobial effects (antibacterial and antifungal activities) of this interaction in order to gain a better understanding of their specific interaction. The analytical study presents a similar chemical profile between the leaves of healthy and herbivorous A. cathartica and the excretions of the caterpillars. The similar analytical profile of the leaves of A. cathartica and the excretions of P. tetrio, and the difference with the caterpillar bodies, suggests a selective excretion of compounds by the caterpillar. The measured antimicrobial activities support the physicochemical tests. The natural products found selectively in the excretions (rather than in the body) could explain the ability of P. tetrio to feed on this toxic Apocynaceae species.

Soil-derived bacteria endow Camellia weevil with more ability to resist plant chemical defense

BACKGROUND

Herbivorous insects acquire their gut microbiota from diverse sources, and these microorganisms play significant roles in insect hosts' tolerance to plant secondary defensive compounds. Camellia weevil (Curculio chinensis) (CW) is an obligate seed parasite of Camellia oleifera plants. Our previous study linked the CW's gut microbiome to the tolerance of the tea saponin (TS) in C. oleifera seeds. However, the source of these gut microbiomes, the key bacteria involved in TS tolerance, and the degradation functions of these bacteria remain unresolved.

RESULTS

Our study indicated that CW gut microbiome was more affected by the microbiome from soil than that from fruits. The soil-derived Acinetobacter served as the core bacterial genus, and Acinetobacter sp. was putatively regarded responsible for the saponin-degradation in CW guts. Subsequent experiments using fluorescently labeled cultures verified that the isolate Acinetobacter sp. AS23 can migrate into CW larval guts, and ultimately endow its host with the ability to degrade saponin, thereby allowing CW to subsist as a pest within plant fruits resisting to higher concentration of defensive chemical.

CONCLUSIONS

The systematic studies of the sources of gut microorganisms, the screening of taxa involved in plant secondary metabolite degradation, and the investigation of bacteria responsible for CW toxicity mitigation provide clarified evidence that the intestinal microorganisms can mediate the tolerance of herbivorous insects against plant toxins. Video Abstract.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses

INTRODUCTION

Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown.

OBJECTIVES

We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species.

METHODS

Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues.

RESULTS

Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors.

CONCLUSION

We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.

Interaction diversity explains the maintenance of phytochemical diversity

The production of complex mixtures of secondary metabolites is a ubiquitous feature of plants. Several evolutionary hypotheses seek to explain how phytochemical diversity is maintained, including the synergy hypothesis, the interaction diversity hypothesis, and the screening hypothesis. We experimentally tested a set of predictions derived from these hypotheses by manipulating the richness and structural diversity of phenolic metabolites in the diets of eight plant consumers. Across 3940 total bioassays, there was clear support for the interaction diversity hypothesis over the synergy or screening hypotheses. The number of consumers affected by a particular phenolic composition increased with increasing richness and structural diversity of compounds. Furthermore, the bioactivity of phenolics was consumer-specific. All compounds tested reduced the performance of at least one consumer, but no compounds affected all consumers. These results show how phytochemical diversity may be maintained in nature by a complex selective landscape exerted by diverse communities of plant consumers.

Whitefly hijacks a plant detoxification gene that neutralizes plant toxins

Plants protect themselves with a vast array of toxic secondary metabolites, yet most plants serve as food for insects. The evolutionary processes that allow herbivorous insects to resist plant defenses remain largely unknown. The whitefly Bemisia tabaci is a cosmopolitan, highly polyphagous agricultural pest that vectors several serious plant pathogenic viruses and is an excellent model to probe the molecular mechanisms involved in overcoming plant defenses. Here, we show that, through an exceptional horizontal gene transfer event, the whitefly has acquired the plant-derived phenolic glucoside malonyltransferase gene BtPMaT1. This gene enables whiteflies to neutralize phenolic glucosides. This was confirmed by genetically transforming tomato plants to produce small interfering RNAs that silence BtPMaT1, thus impairing the whiteflies' detoxification ability. These findings reveal an evolutionary scenario whereby herbivores harness the genetic toolkit of their host plants to develop resistance to plant defenses and how this can be exploited for crop protection.

Antioxidant potential of family Cucurbitaceae with special emphasis on Cucurbita genus: A key to alleviate oxidative stress-mediated disorders

Oxidative stress is the imbalance between reactive oxygen species (ROS) production, and accumulation and the ability of a biological system to clear these reactive products. This imbalance leads to cell and tissue damage causing several disorders in human body, such as neurodegeneration, metabolic problems, cardiovascular diseases, and cancer. Cucurbitaceae family consists of about 100 genera and 1,000 species of plants including mostly tropical, annual or perennial, monoecious, and dioecious herbs. The plants from Cucurbita species are rich sources of phytochemicals and act as a rich source of antioxidants. The most important phytochemicals present in the cucurbits are cucurbitacins, saponins, carotenoids, phytosterols, and polyphenols. These bioactive phyto-constituents are responsible for the pharmacological effects including antioxidant, antitumor, antidiabetic, hepatoprotective, antimicrobial, anti-obesity, diuretic, anti-ulcer activity, and antigenotoxic. A wide number of in vitro and in vivo studies have ascribed these health-promoting effects of Cucurbita genus. Results of clinical trials suggest that Cucurbita provides health benefits for diabetic patients, patients with benign prostate hyperplasia, infertile women, postmenopausal women, and stress urinary incontinence in women. The intend of the present review is to focus on the protective role of Cucurbita spp. phytochemicals on oxidative stress-related disorders on the basis of preclinical and human studies. The review will also give insights on the in vitro and in vivo antioxidant potential of the Cucurbitaceae family as a whole.

Metabolomic differences between invasive alien plants from native and invaded habitats

Globalization facilitated the spread of invasive alien species (IAS), undermining the stability of the world's ecosystems. We investigated the metabolomic profiles of three IAS species: Chromolaena odorata (Asteraceae) Datura stramonium (Solanaceae), and Xanthium strumarium (Asteraceae), comparing metabolites of individual plants in their native habitats (USA), to their invasive counterparts growing in and around Kruger National Park (South Africa, ZA). Metabolomic samples were collected using RApid Metabolome Extraction and Storage (RAMES) technology, which immobilizes phytochemicals on glass fiber disks, reducing compound degradation, allowing long-term, storage and simplifying biochemical analysis. Metabolomic differences were analyzed using ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) of samples eluted from RAMES disks. Partial Least Squares-Discriminant Analysis (PLS-DA) of metabolomes of individual plants allowed statistical separation of species, native and invasive populations of each species, and some populations on the same continent. Invasive populations of all species were more phytochemically diverse than their native counterparts, and their metabolomic profiles were statistically distinguishable from their native relatives. These data may elucidate the mechanisms of successful invasion and rapid adaptive evolution of IAS. Moreover, RAMES technology combined with PLS-DA statistical analysis may allow taxonomic identification of species and, possibly, populations within each species.

Transcriptional reprogramming strategies and miRNA-mediated regulation networks of Taxus media induced into callus cells from tissues

Background

Taxus cells are a potential sustainable and environment-friendly source of taxol, but they have low survival ratios and slow grow rates. Despite these limitations, Taxus callus cells induced through 6 months of culture contain more taxol than their parent tissues. In this work, we utilized 6-month-old Taxus media calli to investigate their regulatory mechanisms of taxol biosynthesis by applying multiomics technologies. Our results provide insights into the adaptation strategies of T. media by transcriptional reprogramming when induced into calli from parent tissues.

Results

Seven out of 12 known taxol, most of flavonoid and phenylpropanoid biosynthesis genes were significantly upregulated in callus cells relative to that in the parent tissue, thus indicating that secondary metabolism is significantly strengthened. The expression of genes involved in pathways metabolizing biological materials, such as amino acids and sugars, also dramatically increased because all nutrients are supplied from the medium. The expression level of 94.1% genes involved in photosynthesis significantly decreased. These results reveal that callus cells undergo transcriptional reprogramming and transition into heterotrophs. Interestingly, common defense and immune activities, such as “plant–pathogen interaction” and salicylic acid- and jasmonic acid-signaling transduction, were repressed in calli. Thus, it’s an intelligent adaption strategy to use secondary metabolites as a cost-effective defense system. MiRNA- and degradome-sequencing results showed the involvement of a precise regulatory network in the miRNA-mediated transcriptional reprogramming of calli. MiRNAs act as direct regulators to enhance the metabolism of biological substances and repress defense activities. Given that only 17 genes of secondary metabolite biosynthesis were effectively regulated, miRNAs are likely to play intermediate roles in the biosynthesis of secondary metabolites by regulating transcriptional factors (TFs), such as ERF, WRKY, and SPL.

Conclusion

Our results suggest that increasing the biosynthesis of taxol and other secondary metabolites is an active regulatory measure of calli to adapt to heterotrophic culture, and this alteration mainly involved direct and indirect miRNA-induced transcriptional reprogramming. These results expand our understanding of the relationships among the metabolism of biological substances, the biosynthesis of secondary metabolites, and defense systems. They also provide a series of candidate miRNAs and transcription factors for taxol biosynthesis.

Coevolution of Snake Venom Toxic Activities and Diet: Evidence that Ecological Generalism Favours Toxicological Diversity

Snake venom evolution is typically considered to be predominantly driven by diet-related selection pressures. Most evidence for this is based on lethality to prey and non-prey species and on the identification of prey specific toxins. Since the broad toxicological activities (e.g., neurotoxicity, coagulotoxicity, etc.) sit at the interface between molecular toxinology and lethality, these classes of activity may act as a key mediator in coevolutionary interactions between snakes and their prey. Indeed, some recent work has suggested that variation in these functional activities may be related to diet as well, but previous studies have been limited in geographic and/or taxonomic scope. In this paper, we take a phylogenetic comparative approach to investigate relationships between diet and toxicological activity classes on a global scale across caenophidian snakes, using the clinically oriented database at toxinology.com. We generally find little support for specific prey types selecting for particular toxicological effects except that reptile-feeders are more likely to be neurotoxic. We find some support for endothermic prey (with higher metabolic rates) influencing toxic activities, but differently from previous suggestions in the literature. More broadly, we find strong support for a general effect of increased diversity of prey on the diversity of toxicological effects of snake venom. Hence, we provide evidence that selection pressures on the toxicological activities of snake venom has largely been driven by prey diversity rather than specific types of prey. These results complement and extend previous work to suggest that specific matching of venom characteristics to prey may occur at the molecular level and translate into venom lethality, but the functional link between those two is not constrained to a particular toxicological route.

Charting Angiosperm Chemistry: Evolutionary Perspective on Specialized Metabolites Reflected in Chemical Property Space

Plants possess an outstanding chemical diversity of specialized metabolites developed to adapt to environmental niches and increase fitness. The observed diversity is hypothesized to result from various evolutionary mechanisms, such as the continuous branching off and extension of existing biosynthetic pathways or enhanced levels of catalytic promiscuity in certain enzymes. In this study, ChemGPS-NP has been employed to chart the distribution and diversity of physicochemical properties for selected types of specialized metabolites from the angiosperms. Utilizing these charts, it is analyzed how different properties of various types of specialized metabolites change in different plant groups, and the chemical diversity from the volume they occupy in chemical property space is evaluated. In this context, possible underlying evolutionary mechanisms are discussed, which could explain the observed distribution and behavior in chemical property space. Based on these studies, it is demonstrated that evolutionary processes in plant specialized metabolism and the resultant metabolic diversification are reflected in chemical property space.

Plant toxins and parasitoid trophic ecology

Parasitoids (parasitic wasps) are ubiquitous components of nearly all communities containing plant-insect herbivore associations. Plant toxin defenses against herbivores may also affect higher trophic levels by directly (e.g., plant toxins encountered in host hemolymph) or indirectly (e.g., plant toxins reduce host size/quality or alter the host's immunity against parasitoids). Yet, whether parasitoids structure plant-herbivore interactions remains relatively understudied. Nevertheless, recent meta-analyses and empirical work emphasize the importance of parasitoids in structuring interactions among lower trophic levels. Two promising areas of research are particularly ripe for future exploration: a) the potential for microbes to alter the interactions among plants, insect herbivores, and parasitoids, and b) the effects of climate change on phenological (mis)matches among trophic levels.

The evolution of variance in sequential defences

The defences used by organisms against predators display a great degree of variability. Defence phenotypes can differ substantially among individuals of the same species, and a single individual can itself deploy a variety of defences. Here, we use a mathematical model that includes mutation and selection to understand the evolutionary origin of this variability in a population of a species that deploys defences sequentially ("first" and "second" defences). Typically, the first defence evolves to have lower variance, i.e. appears more closely accumulated around the ideal phenotype, than the second defence (even when the breaching the first defence incurs more fitness loss than breaching the second defence with the other parameters the same for both defences). However, if the first defence is much less effective in repelling predators, or is much less tolerant of deviation from the ideal phenotype, then the first defence can evolve to have higher variance than the second. Other factors like mutation strength and the losses in the fitness when each defence fails also influence the defence variance. Larger mutation rate incurs larger equilibrium variances, and when the comparative importance in fitness of one defence increases, then the ratio between the variances of this defence and the other defence decreases. Sequentially acting defences are found in many organisms, so we encourage empirical research to test our theoretical predictions.

Gastropods and Insects Prefer Different Solanum dulcamara Chemotypes

Solanum dulcamara (Bittersweet nightshade) shows significant intraspecific variation in glycoalkaloid (GA) composition and concentration. We previously showed that constitutive differences in overall GA levels are correlated with feeding preference of the grey field slug (GFS; Deroceras reticulatum). One particularly preferred accession, ZD11, contained low GA levels, but high levels of previously unknown structurally related uronic acid conjugated compounds (UACs). Here we test whether different slug species as well as insect herbivores show similar feeding preferences among six S. dulcamara accessions with different GA chemotypes. In addition, we investigate whether slug feeding can lead to induced changes in the chemical composition and affect later arriving herbivores. A leaf disc assay using greenhouse-grown plants showed that three slug species similarly preferred accessions with low GA levels. Untargeted metabolomic analyses showed that previous slug feeding consistently increased the levels of N-caffeoyl-putrescine and a structurally related metabolite, but not the levels of GAs and UACs. Slug-induced responses only affected slug preference in one accession. A common garden experiment using the same six accessions revealed that ZD11 received the highest natural gastropod feeding damage, but suffered the lowest damage by specialist flea beetles. The latter preferred to feed on accessions with high GA levels. Our study indicates that different selection pressures imposed by generalist gastropods and specialist insects may explain part of the observed chemical diversity in S. dulcamara.

Diffusible signal factor (DSF)-mediated quorum sensing modulates expression of diverse traits in Xanthomonas citri and responses of citrus plants to promote disease

Background

The gram-negative Xanthomonas genus contains a large group of economically important plant pathogens, which cause severe diseases on many crops worldwide. The diffusible signal factor (DSF) - mediated quorum sensing (QS) system coordinates expression of virulence factors in plant pathogenic Xanthomonas spp. However, the regulatory effects of this system during the Xanthomonas- plant interactions remain unclear from both the pathogen and host aspects.

Results

In this study, we investigated the in planta DSF- mediated QS regulon of X. citri subsp. citri (Xac), the causal agent of citrus canker. We also characterized the transcriptional responses of citrus plants to DSF-mediated Xac infection via comparing the gene expression patterns of citrus trigged by wild type Xac strain 306 with those trigged by its DSF- deficient (∆rpfF) mutant using the dual RNA-seq approach. Comparative global transcript profiles of Xac strain 306 and the ∆rpfF mutant during host infection revealed that DSF- mediated QS specifically modulates bacterial adaptation, nutrition uptake and metabolisms, stress tolerance, virulence, and signal transduction to favor host infection. The transcriptional responses of citrus to DSF-mediated Xac infection are characterized by downregulation of photosynthesis genes and plant defense related genes, suggesting photosynthetically inactive reactions and repression of defense responses. Alterations of phytohormone metabolism and signaling pathways were also triggered by DSF-mediated Xac infection to benefit the pathogen.

Conclusions

Collectively, our findings provide new insight into the DSF- mediated QS regulation during plant-pathogen interactions and advance the understanding of traits used by Xanthomonas to promote infection on host plants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5384-4) contains supplementary material, which is available to authorized users.

Eavesdropping on gall-plant interactions: the importance of the signaling function of induced volatiles

The galling insect manipulates the host plant tissue to its own benefit, building the gall structure where it spends during most of its life cycle. These specialist herbivore insects can induce and manipulate plant structure and metabolism throughout gall development and may affect plant volatile emission. Consequently, volatile emission from altered metabolism contribute to eavesdropping cueing. Eavesdropping can be part of adaptive strategies used by evolution for both galling insects and the entire-associated community in order to cue some interaction response. This is in contrast to some herbivores associated with delayed induced responses, altering plant metabolites during the short time while they feed. Due to the different lifestyles of the galling organism, which are associated with different plant tissues and organs (e.g leaves, flowers or fruits), a distinct diversity of organisms may eavesdrop on induced volatiles interacting with the galls. Furthermore, the eavesdropping cues may be defined according to the phenological coupling between galling organism and host plant, which results from the development of a gall structure. For instance, when plants release volatile-induced defenses after galling insects' activity, another interactor may perceive these volatiles and change its behavior and interactions with host plants and galls. Thus, natural enemies could be attracted by different volatiles emitted by the gall tissues. Considering the duration of the life cycle of the galling organism and the gall, the temporal extent of gall-induced volatiles may include more persistent volatile cues and eavesdropping effects than the volatiles induced by non-galling herbivores. Accordingly, from chemical ecology perspective we expect that galling herbivore-induced volatiles may exhibit robust effects on neighboring-plant interactions including those ones during different plant developmental or phenological periods. Information about multitrophic interactions between insects and plants supports the additional understanding of direct and indirect effects, and allows insight into new hypotheses.

Glycoalkaloid composition explains variation in slug resistance in Solanum dulcamara

In natural environments, plants have to deal with a wide range of different herbivores whose communities vary in time and space. It is believed that the chemical diversity within plant species has mainly arisen from selection pressures exerted by herbivores. So far, the effects of chemical diversity on plant resistance have mostly been assessed for arthropod herbivores. However, also gastropods, such as slugs, can cause extensive damage to plants. Here we investigate to what extent individual Solanum dulcamara plants differ in their resistance to slug herbivory and whether this variation can be explained by differences in secondary metabolites. We performed a series of preference assays using the grey field slug (Deroceras reticulatum) and S. dulcamara accessions from eight geographically distinct populations from the Netherlands. Significant and consistent variation in slug preference was found for individual accessions within and among populations. Metabolomic analyses showed that variation in steroidal glycoalkaloids (GAs) correlated with slug preference; accessions with high GA levels were consistently less damaged by slugs. One, strongly preferred, accession with particularly low GA levels contained high levels of structurally related steroidal compounds. These were conjugated with uronic acid instead of the glycoside moieties common for Solanum GAs. Our results illustrate how intraspecific variation in steroidal glycoside profiles affects resistance to slug feeding. This suggests that also slugs should be considered as important drivers in the co-evolution between plants and herbivores.

Origin and maintenance of chemical diversity in a species-rich tropical tree lineage

Plant secondary metabolites play important ecological and evolutionary roles, most notably in the deterrence of natural enemies. The classical theory explaining the evolution of plant chemical diversity is that new defences arise through a pairwise co-evolutionary arms race between plants and their specialized natural enemies. However, plant species are bombarded by dozens of different herbivore taxa from disparate phylogenetic lineages that span a wide range of feeding strategies and have distinctive physiological constraints that interact differently with particular plant metabolites. How do plant defence chemicals evolve under such multiple and potentially contrasting selective pressures imposed by diverse herbivore communities? To tackle this question, we exhaustively characterized the chemical diversity and insect herbivore fauna from 31 sympatric species of Amazonian Protieae (Burseraceae) trees. Using a combination of phylogenetic, metabolomic and statistical learning tools, we show that secondary metabolites that were associated with repelling herbivores (1) were more frequent across the Protieae phylogeny and (2) were found in average higher abundance than other compounds. Our findings suggest that generalist herbivores can play an important role in shaping plant chemical diversity and support the hypothesis that chemical diversity can also arise from the cumulative outcome of multiple diffuse interactions.

The Hardy Rubber Tree Genome Provides Insights into the Evolution of Polyisoprene Biosynthesis

Eucommia ulmoides, also called hardy rubber tree, is an economically important tree; however, the lack of its genome sequence restricts the fundamental biological research and applied studies of this plant species. Here, we present a high-quality assembly of its ∼1.2-Gb genome (scaffold N50 = 1.88 Mb) with at least 26 723 predicted genes for E. ulmoides, the first sequenced genome of the order Garryales, which was obtained using an integrated strategy combining Illumina sequencing, PacBio sequencing, and BioNano mapping. As a sister taxon to lamiids and campanulids, E. ulmoides underwent an ancient genome triplication shared by core eudicots but no further whole-genome duplication in the last ∼125 million years. E. ulmoides exhibits high expression levels and/or gene number expansion for multiple genes involved in stress responses and the biosynthesis of secondary metabolites, which may account for its considerable environmental adaptability. In contrast to the rubber tree (Hevea brasiliensis), which produces cis-polyisoprene, E. ulmoides has evolved to synthesize long-chain trans-polyisoprene via farnesyl diphosphate synthases (FPSs). Moreover, FPS and rubber elongation factor/small rubber particle protein gene families were expanded independently from the H. brasiliensis lineage. These results provide new insights into the biology of E. ulmoides and the origin of polyisoprene biosynthesis.

Does coevolution with a shared parasite drive hosts to partition their defences among species?

When mimicry imposes costs on models, selection may drive the model's phenotype to evolve away from its mimic. For example, brood parasitism often drives hosts to diversify in egg appearance among females within a species, making mimetic parasitic eggs easier to detect. However, when a single parasite species exploits multiple host species, parasitism could also drive host egg evolution away from other co-occurring hosts, to escape susceptibility to their respective mimics. This hypothesis predicts that sympatric hosts of the same parasite should partition egg phenotypic space (defined by egg colour, luminance and pattern) among species to avoid one another. We show that eggs of warbler species parasitized by the cuckoo finch Anomalospiza imberbis in Zambia partition phenotypic space much more distinctly than do eggs of sympatric but unparasitized warblers. Correspondingly, cuckoo finch host-races better match their own specialist host than other local host species. In the weaver family, parasitized by the diederik cuckoo Chrysococcyx caprius, by contrast, parasitized species were more closely related and overlapped extensively in phenotypic space; correspondingly, cuckoos did not match their own host better than others. These results suggest that coevolutionary arms races between hosts and parasites may be shaped by the wider community context in which they unfold.

Papain-like cysteine proteases in Carica papaya: lineage-specific gene duplication and expansion

Background

Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases structurally related to papain, play important roles in plant development, senescence, and defense responses. Papain, the first cysteine protease whose structure was determined by X-ray crystallography, plays a crucial role in protecting papaya from herbivorous insects. Except the four major PLCPs purified and characterized in papaya latex, the rest of the PLCPs in papaya genome are largely unknown.

Results

We identified 33 PLCP genes in papaya genome. Phylogenetic analysis clearly separated plant PLCP genes into nine subfamilies. PLCP genes are not equally distributed among the nine subfamilies and the number of PLCPs in each subfamily does not increase or decrease proportionally among the seven selected plant species. Papaya showed clear lineage-specific gene expansion in the subfamily III. Interestingly, all four major PLCPs purified from papaya latex, including papain, chymopapain, glycyl endopeptidase and caricain, were grouped into the lineage-specific expansion branch in the subfamily III. Mapping PLCP genes on chromosomes of five plant species revealed that lineage-specific expansions of PLCP genes were mostly derived from tandem duplications. We estimated divergence time of papaya PLCP genes of subfamily III. The major duplication events leading to lineage-specific expansion of papaya PLCP genes in subfamily III were estimated at 48 MYA, 34 MYA, and 16 MYA. The gene expression patterns of the papaya PLCP genes in different tissues were assessed by transcriptome sequencing and qRT-PCR. Most of the papaya PLCP genes of subfamily III expressed at high levels in leaf and green fruit tissues.

Conclusions

Tandem duplications played the dominant role in affecting copy number of PLCPs in plants. Significant variations in size of the PLCP subfamilies among species may reflect genetic adaptation of plant species to different environments. The lineage-specific expansion of papaya PLCPs of subfamily III might have been promoted by the continuous reciprocal selective effects of herbivore attack and plant defense.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4394-y) contains supplementary material, which is available to authorized users.