Transcriptome analysis of Barbarea vulgaris infested with diamondback moth (Plutella xylostella) larvae

Dynamics of Zea mays transcriptome in response to a polyphagous herbivore, Spodoptera litura

Zea mays defense response is well-crafted according to the physical and chemical weapons utilized by their invaders during the coevolutionary period. Maize plants employ diversified defense strategies and alter the spatiotemporal distribution of several classes of defensive compounds to affect insect herbivore performance. However, only little knowledge is available about the defense orchestration of maize in response to Spodoptera litura, a voracious Noctuidae pest. In order to decipher the defense status of Zea mays (African tall variety) against S. litura, a comparative feeding bioassay was executed, which revealed reduced performance of the herbivore on maize. In order to understand the molecular mechanism behind maize tolerance against S. litura, a microarray-based genome-wide expression analysis was performed. The comparative analysis displayed 792 differentially expressed genes (DEGs), wherein 357 genes were upregulated and 435 genes were downregulated at fold change ≥ 2 and p value ≤ 0.05. The upregulated genes were identified and categorized as defense-related, oxidative stress-related, transcription regulatory genes, protein synthesis genes, phytohormone-related, and primary and secondary metabolism-related. In contrast, downregulated genes were mainly associated with plant growth and development, indicating a balance of growth and defense response and utilization of a highly evolved C-diversion response were noticed. Maize plants showed better tolerance against herbivory and maintained its fitness using a combinatorial strategy. This peculiar response of Zea mays against S. litura offers an excellent possibility of managing polyphagous pests by spicing up the plant's defensive response with tolerance mechanism.

Engineering and optimization of the 2-phenylethylglucosinolate production in Nicotiana benthamiana by combining biosynthetic genes from Barbarea vulgaris and Arabidopsis thaliana

2-Phenylethylglucosinolate (2PE) derived from homophenylalanine is present in plants of the Brassicales order as a defense compound. It is associated with multiple biological properties, including deterrent effects on pests and antimicrobial and health-promoting functions, due to its hydrolysis product 2-phenylethyl isothiocyanate, which confers 2PE as a potential application in agriculture and industry. In this study, we characterized the putative key genes for 2PE biosynthesis from Barbarea vulgaris W.T. Aiton and demonstrated the feasibility of engineering 2PE production in Nicotiana benthamiana Domin. We used different combinations of genes from B. vulgaris and Arabidopsis thaliana (L.) Heynh. to demonstrate that: (i) BvBCAT4 performed more efficiently than AtBCAT4 in biosynthesis of both homophenylalanine and dihomomethionine; (ii) MAM1 enzymes were critical for the chain-elongated profile, while CYP79F enzymes accepted both chain-elongated methionine and homophenylalanine; (iii) aliphatic but not aromatic core structure pathway catalyzed the 2PE biosynthesis; (iv) a chimeric pathway containing BvBCAT4, BvMAM1, AtIPMI and AtIPMDH1 resulted in a two-fold increase in 2PE production compared with the B. vulgaris-specific chain elongation pathway; and (v) profiles of chain-elongated products and glucosinolates partially mirrored the profiles in the gene donor plant, but were wider in N. benthamiana than in the native plants. Our study provides a strategy to produce the important homophenylalanine and 2PE in a heterologous host. Furthermore, chimeric engineering of the complex 2PE biosynthetic pathway enabled detailed understanding of catalytic properties of individual enzymes - a prerequisite for understanding biochemical evolution. The new-to-nature gene combinations have the potential for application in biotechnological and plant breeding.

Transcriptome analysis of maize reveals potential key genes involved in the response to belowground herbivore Holotrichia parallela larvae feeding

The larvae of Holotrichia parallela, a destructive belowground herbivore, causes tremendous damages to maize plants. However, little is known if there are any defense mechanisms in maize roots to defend themselves against this herbivore. In the current research, we carried out RNA-sequencing to investigate the changes in gene transcription level in maize roots after H. parallela larvae infestation. A total of 644 up-regulated genes and 474 down-regulated genes was found. In addition, Gene ontology (GO) annotation analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed. Weighted gene co-expression network analysis (WGCNA) indicated that peroxidase genes may be the hub genes that regulate maize defenses to H. parallela larvae attack. We also found 105 transcription factors, 44 hormone-related genes, and 62 secondary metabolism-related genes within differentially expressed genes (DEGs). Furthermore, the expression profiles of 12 DEGs from the transcriptome analysis were confirmed by quantitative real-time PCR experiments. This transcriptome analysis provides insights into the molecular mechanisms of the underground defense in maize roots to H. parallela larvae attack and will help to select target genes of maize for defense against belowground herbivory.

Glucosinolate turnover in Brassicales species to an oxazolidin-2-one, formed via the 2-thione and without formation of thioamide

Glucosinolates are found in plants of the order Brassicales and hydrolyzed to different breakdown products, particularly after tissue damage. In Barbarea vulgaris R.Br. (Brassicaceae), the dominant glucosinolate in the investigated "G-type" is glucobarbarin, (S)-2-hydroxy-2-phenylethylglucosinolate. Formation of the nitrile from glucobarbarin was observed in vitro, while a previously suggested thioamide (synonym thionamide) was not confirmed. Resedine (5-phenyl-1,3-oxazolidin-2-one) was detected after glucobarbarin hydrolysis in crushed B. vulgaris leaves and siliques, but not in intact parts. The abundance increased for several hours after completion of hydrolysis. The corresponding 1,3-oxazolidine-2-thione (OAT), with the common name barbarin, was also formed, and appeared to be the precursor of resedine. Addition of each of two non-endogenous OATs, (S)-5-ethyl-5-methylOAT and (R)-5-vinylOAT (R-goitrin), to a leaf homogenate resulted in formation of the corresponding 1,3-oxazolidin-2-ones (OAOs), confirming the metabolic connection of OAT to OAO. Formation of OAOs was inhibited by prior brief heating of the homogenate, suggesting enzyme involvement. We suggest the conversion of OATs to OAOs to be catalyzed by an enzyme ("oxazolidinethionase") responsible for turnover of OAT formed in intact plants. Resedine had been reported as an alkaloid from another species - Reseda luteola L. (Resedaceae) - naturally containing the glucosinolate glucobarbarin. However, resedine was not detected in intact R. luteola plants, but formed after tissue damage. The formation of resedine in two families suggests a broad distribution of putative OATases in the Brassicales; potentially involved in glucosinolate turnover that needs myrosinase activity as the committed step. In agreement with the proposed function of OATase, several candidate genes for myrosinases in glucosinolate turnover in intact plants were discovered in the B. vulgaris genome. We also suggest that biotechnological conversion of OATs to OAOs might improve the nutritional value of Brassicales protein. HPLC-MS/MS methods for detection of these glucobarbarin products are described.

The early transcriptome response of cassava (Manihot esculenta Crantz) to mealybug (Phenacoccus manihoti) feeding

The mealybug, Phenacoccus manihoti, is a leading pest of cassava (Manihot esculenta Crantz), damaging this crop globally. Although the biological control of this mealybug using natural predators has been established, resistance breeding remains an important means of control. Understanding plant responses to insect herbivory, by determining and identifying differentially expressed genes (DEGs), is a vital step towards the understanding of molecular mechanisms of defence responses in plants and the development of resistant cultivars by gene editing. Morphological and molecular analysis confirmed the mealybug identity as Phenacoccus manihoti (Matile-Ferrero). The transcriptome response of the green mite resistant cassava genotype AR23.1 was compared to P40/1 with no known resistance at 24 and 72 hours of mealybug infestation compared to non-infested mock. A total of 301 and 206 genes were differentially expressed at 24 and 72 of mealybug infestation for AR23.1 and P40/1 genotypes respectively, using a log2 fold change and P-value ≤ 0.05. Gene ontology functional classification revealed an enrichment of genes in the secondary metabolic process category in AR23.1 in comparison with P40/1, while genes in the regulation of molecular function, cellular component biogenesis and electron carrier categories were more significantly enriched in P40/1 than in AR23.1. Biological pathway analysis, based on KEGG, revealed a significant enrichment of plant-pathogen interaction and plant hormonal signal transduction pathways for a cohort of up-regulated and down-regulated DEGs in both genotypes. Defence-related genes such as 2-oxogluterate, gibberellin oxidase and terpene synthase proteins were only induced in genotype AR23.1 and not in P40/1, and subsequently validated by RT-qPCR. The study revealed a difference in response to mealybug infestation in the two genotypes studied, with AR23.1 showing a higher number of differentially expressed transcripts post mealybug infestation at 24 and 72 hours. Candidate defence-related genes that were overexpressed in the AR23.1 genotype post mealybug infestation will be useful in future functional studies towards the control of mealybugs.

A tandem array of UDP-glycosyltransferases from the UGT73C subfamily glycosylate sapogenins, forming a spectrum of mono- and bisdesmosidic saponins

KEY MESSAGE

This study identifies six UGT73Cs all able to glucosylate sapogenins at positions 3 and/or 28 which demonstrates that B. vulgaris has a much richer arsenal of UGTs involved in saponin biosynthesis than initially anticipated. The wild cruciferous plant Barbarea vulgaris is resistant to some insects due to accumulation of two monodesmosidic triterpenoid saponins, oleanolic acid 3-O-β-cellobioside and hederagenin 3-O-β-cellobioside. Insect resistance depends on the structure of the sapogenin aglycone and the glycosylation pattern. The B. vulgaris saponin profile is complex with at least 49 saponin-like metabolites, derived from eight sapogenins and including up to five monosaccharide units. Two B. vulgaris UDP-glycosyltransferases, UGT73C11 and UGT73C13, O-glucosylate sapogenins at positions 3 and 28, forming mainly 3-O-β-D-glucosides. The aim of this study was to identify UGTs responsible for the diverse saponin oligoglycoside moieties observed in B. vulgaris. Twenty UGT genes from the insect resistant genotype were selected and heterologously expressed in Nicotiana benthamiana and/or Escherichia coli. The extracts were screened for their ability to glycosylate sapogenins (oleanolic acid, hederagenin), the hormone 24-epibrassinolide and sapogenin monoglucosides (hederagenin and oleanolic acid 3-O-β-D-glucosides). Six UGTs from the UGT73C subfamily were able to glucosylate both sapogenins and both monoglucosides at positions 3 and/or 28. Some UGTs formed bisdesmosidic saponins efficiently. At least four UGT73C genes were localized in a tandem array with UGT73C11 and possibly UGT73C13. This organization most likely reflects duplication events followed by sub- and neofunctionalization. Indeed, signs of positive selection on several amino acid sites were identified and modelled to be localized on the UGT protein surface. This tandem array is proposed to initiate higher order bisdesmosidic glycosylation of B. vulgaris saponins, leading to the recently discovered saponin structural diversity, however, not directly to known cellobiosidic saponins.

Chemical Compounds and Bioactivity of Aqueous Extracts of Alibertia spp. in the Control of Plutella xylostella L. (Lepidoptera: Plutellidae)

Successive applications of insecticides to control Plutella xylostella L. (Lepidoptera: Plutellidae) have resulted in the emergence of resistant populations of this insect. A novel control measure for this target insect could be the use of botanical insecticides derived from plant tissues. Hence, we experimentally tested aqueous extracts of Alibertia edulis (Rich.), Alibertia intermedia (Mart.), and Alibertia sessilis (Vell.) K. Schum. found in the Brazilian savannah in order to investigate their potential to disrupt the life cycle of P. xylostella. Aqueous extracts of the leaves of A. intermedia and A. sessilis negatively affected the development of P. xylostella in all stages of the life cycle, prolonging the larval stage and causing mortality in the larval or pupal stages. Treatments with A. intermedia and A. sessilis extracts caused the lowest fecundity and the number of hatched larvae. The harmful effects of these aqueous extracts on the life cycle of P. xylostella may be attributable to the flavonoids and other phenolic compounds present in A. intermedia and A. sessilis. These aqueous botanical extracts are low in toxicity when compared to non-aqueous pesticides, and may emerge as an effective approach for control of populations of P. xylostella.

Simulated herbivory in chickpea causes rapid changes in defense pathways and hormonal transcription networks of JA/ethylene/GA/auxin within minutes of wounding

Chickpea (C. arietinum L.) is an important pulse crop in Asian and African countries that suffers significant yield losses due to attacks by insects like H. armigera. To obtain insights into early responses of chickpea to insect attack, a transcriptomic analysis of chickpea leaves just 20 minutes after simulated herbivory was performed, using oral secretions of H. armigera coupled with mechanical wounding. Expression profiles revealed differential regulation of 8.4% of the total leaf transcriptome with 1334 genes up-regulated and 501 down-regulated upon wounding at log₂-fold change (|FC| ≤ −1 and ≥1) and FDR value ≤ 0.05. In silico analysis showed the activation of defenses through up-regulation of genes of the phenylpropanoid pathway, pathogenesis, oxidases and CYTP450 besides differential regulation of kinases, phosphatases and transcription factors of the WRKY, MYB, ERFs, bZIP families. A substantial change in the regulation of hormonal networks was observed with up-regulation of JA and ethylene pathways and suppression of growth associated hormone pathways like GA and auxin within 20 minutes of wounding. Secondary qPCR comparison of selected genes showed that oral secretions often increased differential expression relative to mechanical damage alone. The studies provide new insights into early wound responses in chickpea.

The genome sequence of Barbarea vulgaris facilitates the study of ecological biochemistry

The genus Barbarea has emerged as a model for evolution and ecology of plant defense compounds, due to its unusual glucosinolate profile and production of saponins, unique to the Brassicaceae. One species, B. vulgaris, includes two ‘types’, G-type and P-type that differ in trichome density, and their glucosinolate and saponin profiles. A key difference is the stereochemistry of hydroxylation of their common phenethylglucosinolate backbone, leading to epimeric glucobarbarins. Here we report a draft genome sequence of the G-type, and re-sequencing of the P-type for comparison. This enables us to identify candidate genes underlying glucosinolate diversity, trichome density, and study the genetics of biochemical variation for glucosinolate and saponins. B. vulgaris is resistant to the diamondback moth, and may be exploited for “dead-end” trap cropping where glucosinolates stimulate oviposition and saponins deter larvae to the extent that they die. The B. vulgaris genome will promote the study of mechanisms in ecological biochemistry to benefit crop resistance breeding.

Identification of anthocyanin biosynthesis related microRNAs in a distinctive Chinese radish (Raphanus sativus L.) by high-throughput sequencing

Anthocyanins are widely distributed water-soluble phytochemical pigments belonging to the flavonoid group. To date, limited knowledge is available about the regulatory roles of miRNAs in anthocyanin biosynthesis in plants. To identify the miRNAs associated with anthocyanin biosynthesis in radish, five small RNA (sRNA) libraries constructed from 'Xinlimei' radish roots at 11, 21, 44, 56 and 73 days (d) were examined using high-throughput sequencing technology. A total of 102.02 million (M) clean reads were generated, from which 483 known and 1415 novel miRNAs were identified. Combined with target prediction and annotation, 72 differentially expressed miRNAs (52 known and 20 novel miRNAs) were more likely to participate in anthocyanin biosynthesis. Several target genes for these miRNAs encode a few transcription factors, including Myb domain (MYB), basic helix-loop-helix (bHLH), WD40 repeat, squamosa promoter binding protein like (SPL), auxin response factor (ARF), ethylene insensitive 3 (EIN3), WRKY and MADS-box proteins. Furthermore, the expression patterns of some anthocyanin biosynthesis related miRNAs and their corresponding targets were validated by RT-qPCR. Based on the characterization of anthocyanin biosynthesis related miRNAs and their target genes, a putative miRNA-target module regulating anthocyanin biosynthesis was proposed. This study represents the first genome-wide identification of miRNAs associated with anthocyanin biosynthesis in radish, and provides insights into the molecular mechanisms underlying regulation of anthocyanin biosynthesis in radish and other crops.

Transcriptome Analysis of Ramie (Boehmeria nivea L. Gaud.) in Response to Ramie Moth (Cocytodes coerulea Guenée) Infestation

The ramie moth Cocytodes coerulea Guenée (RM) is an economically important pest that seriously impairs the yield of ramie, an important natural fiber crop. The molecular mechanisms that underlie the ramie-pest interactions are unclear up to date. Therefore, a transcriptome profiling analysis would aid in understanding the ramie defense mechanisms against RM. In this study, we first constructed two cDNA libraries derived from RM-challenged (CH) and unchallenged (CK) ramie leaves. The subsequent sequencing of the CH and CK libraries yielded 40.2 and 62.8 million reads, respectively. Furthermore, de novo assembling of these reads generated 26,759 and 29,988 unigenes, respectively. An integrated assembly of data from these two libraries resulted in 46,533 unigenes, with an average length of 845 bp per unigene. Among these genes, 24,327 (52.28%) were functionally annotated by predicted protein function. A comparative analysis of the CK and CH transcriptome profiles revealed 1,980 differentially expressed genes (DEGs), of which 750 were upregulated and 1,230 were downregulated. A quantitative real-time PCR (qRT-PCR) analysis of 13 random selected genes confirmed the gene expression patterns that were determined by Illumina sequencing. Among the DEGs, the expression patterns of transcription factors, protease inhibitors, and antioxidant enzymes were studied. Overall, these results provide useful insights into the defense mechanism of ramie against RM.

Aromatic Glucosinolate Biosynthesis Pathway in Barbarea vulgaris and its Response to Plutella xylostella Infestation

The inducibility of the glucosinolate resistance mechanism is an energy-saving strategy for plants, but whether induction would still be triggered by glucosinolate-tolerant Plutella xylostella (diamondback moth, DBM) after a plant had evolved a new resistance mechanism (e.g., saponins in Barbara vulgaris) was unknown. In B. vulgaris, aromatic glucosinolates derived from homo-phenylalanine are the dominant glucosinolates, but their biosynthesis pathway was unclear. In this study, we used G-type (pest-resistant) and P-type (pest-susceptible) B. vulgaris to compare glucosinolate levels and the expression profiles of their biosynthesis genes before and after infestation by DBM larvae. Two different stereoisomers of hydroxylated aromatic glucosinolates are dominant in G- and P-type B. vulgaris, respectively, and are induced by DBM. The transcripts of genes in the glucosinolate biosynthesis pathway and their corresponding transcription factors were identified from an Illumina dataset of G- and P-type B. vulgaris. Many genes involved or potentially involved in glucosinolate biosynthesis were induced in both plant types. The expression patterns of six DBM induced genes were validated by quantitative PCR (qPCR), while six long-fragment genes were validated by molecular cloning. The core structure biosynthetic genes showed high sequence similarities between the two genotypes. In contrast, the sequence identity of two apparent side chain modification genes, the SHO gene in the G-type and the RHO in P-type plants, showed only 77.50% identity in coding DNA sequences and 65.48% identity in deduced amino acid sequences. The homology to GS-OH in Arabidopsis, DBM induction of the transcript and a series of qPCR and glucosinolate analyses of G-type, P-type and F1 plants indicated that these genes control the production of S and R isomers of 2-hydroxy-2-phenylethyl glucosinolate. These glucosinolates were significantly induced by P. xylostella larvae in both the susceptiple P-type and the resistant G-type, even though saponins are the main DBM-resistance causing metabolites in G-type plants. Indol-3-ylmethylglucosinolate was induced in the G-type only. These data will aid our understanding of the biosynthesis and induction of aromatic glucosinolates at the molecular level and also increase our knowledge of the complex mechanisms underpinning defense induction in plants.

Identification and genome organization of saponin pathway genes from a wild crucifer, and their use for transient production of saponins in Nicotiana benthamiana

The ability to evolve novel metabolites has been instrumental for the defence of plants against antagonists. A few species in the Barbarea genus are the only crucifers known to produce saponins, some of which make plants resistant to specialist herbivores, like Plutella xylostella, the diamondback moth. Genetic mapping in Barbarea vulgaris revealed that genes for saponin biosynthesis are not clustered but are located in different linkage groups. Using co-location with quantitative trait loci (QTLs) for resistance, transcriptome and genome sequences, we identified two 2,3-oxidosqualene cyclases that form the major triterpenoid backbones. LUP2 mainly produces lupeol, and is preferentially expressed in insect-susceptible B. vulgaris plants, whereas LUP5 produces β-amyrin and α-amyrin, and is preferentially expressed in resistant plants; β-amyrin is the backbone for the resistance-conferring saponins in Barbarea. Two loci for cytochromes P450, predicted to add functional groups to the saponin backbone, were identified: CYP72As co-localized with insect resistance, whereas CYP716As did not. When B. vulgaris sapogenin biosynthesis genes were transiently expressed by CPMV-HT technology in Nicotiana benthamiana, high levels of hydroxylated and carboxylated triterpenoid structures accumulated, including oleanolic acid, which is a precursor of the major resistance-conferring saponins. When the B. vulgaris gene for sapogenin 3-O-glucosylation was co-expressed, the insect deterrent 3-O-oleanolic acid monoglucoside accumulated, as well as triterpene structures with up to six hexoses, demonstrating that N. benthamiana further decorates the monoglucosides. We argue that saponin biosynthesis in the Barbarea genus evolved by a neofunctionalized glucosyl transferase, whereas the difference between resistant and susceptible B. vulgaris chemotypes evolved by different expression of oxidosqualene cyclases (OSCs).

Derivatization of isothiocyanates and their reactive adducts for chromatographic analysis

Isothiocyanates form adducts with a multitude of biomolecules, and these adducts need analytical methods. Likewise, analytical methods for hydrophilic isothiocyanates are needed. We considered reaction with ammonia to form thiourea derivatives. The hydrophilic, glycosylated isothiocyanate moringin, 4-(α-L-rhamnopyranosyloxy)benzyl isothiocyanate, was efficiently derivatized to the thiourea derivative by incubation with ammonia. The hydrophobic benzyl isothiocyanate was also efficiently derivatized to the thiourea derivative. The thiourea group provided a UV absorbing chromophore, and the derivatives showed expectable sodium and hydrogen adducts in ion trap mass spectrometry and were suitable for liquid chromatography analysis. Reactive dithiocarbamate adducts constitute the major type of reactive ITC adduct expected in biological matrices. Incubation of a model dithiocarbamate with ammonia likewise resulted in conversion to the corresponding thiourea derivative, suggesting that a variety of matrix-bound reactive isothiocyanate adducts can be determined using this strategy. As an example of the application of the method, recovery of moringin and benzyl isothiocyanate applied to cabbage leaf discs was studied in simulated insect feeding assays. The majority of moringin was recovered as native isothiocyanate, but a major part of benzyl isothiocyanate was converted to reactive adducts.

Glucosinolate hydrolysis products in the crucifer Barbarea vulgaris include a thiazolidine-2-one from a specific phenolic isomer as well as oxazolidine-2-thiones

Two isomeric phenolic glucosinolates, m- and p-hydroxyl derivatives of epiglucobarbarin [(R)-2-hydroxy-2-phenylethylglucosinolate], co-occur in an eastern chemotype (P-type) of the crucifer Barbarea vulgaris along with epiglucobarbarin itself. Levels of the phenolic derivatives in B. vulgaris were low in summer but higher during fall and winter, allowing isolation of all three glucosinolates. Hydrolysis in vitro, catalyzed by Sinapis alba myrosinase at near neutral pH, resulted in expectable oxazolidine-2-thione type hydrolysis products of epiglucobarbarin and its m-hydroxyl derivative. In contrast, a thiazolidine-2-one type product was formed in vitro from p-hydroxy epiglucobarbarin and characterized by UV, IR, MS/MS and 2D NMR. Maceration of leaf material resulted in disappearance of the glucosinolates and formation of the same oxazolidine-2-thione and thiazolidine-2-one products as found in vitro. The detected amounts were comparable to initial amounts of precursor glucosinolates. The corresponding oxazolidine-2-thione type product was also detected quantitatively from glucobarbarin in foliage of a western genotype (G-type). We suggest that p-hydroxy epiglucobarbarin is initially converted into the conventional oxazolidine-2-thione, which would further rearrange to a thiazolidine-2-one due to the activating effect of the p-hydroxyl group. We conclude that a subtle difference between isomeric phenolic glucosinolates results in significantly different natural hydrolysis products.

Expression patterns, molecular markers and genetic diversity of insect-susceptible and resistant Barbarea genotypes by comparative transcriptome analysis

BACKGROUND

Barbarea vulgaris contains two genotypes: the glabrous type (G-type), which confers resistance to the diamondback moth (DBM) and other insect pests, and the pubescent type (P-type), which is susceptible to the DBM. Herein, the transcriptomes of P-type B. vulgaris before and after DBM infestation were subjected to Illumina (Solexa) pyrosequencing and comparative analysis.

RESULTS

5.0 gigabase pairs of clean nucleotides were generated. Non-redundant unigenes (33,721) were assembled and 94.1 % of them were annotated. Compared with our previous G-type transcriptome, the expression patterns of many insect responsive genes, including those related to secondary metabolism, phytohormones and transcription factors, which were significantly induced by DBM in G-type plants, were less sensitive to DBM infestation in P-type plants. The genes of the triterpenoid saponin pathway were identified in both G- and P-type plants. The upstream genes of the pathway showed similar expression patterns between the two genotypes. However, gene expression for two downstream enzymes, the glucosyl transferase (UGT73C11) and an oxidosqualene cyclase (OSC), were significantly upregulated in the P-type compared with the G-type plant. The homologous genes from P- and G-type plants were detected by BLAST unigenes with a cutoff level E-value < e(-10). 12,980 gene families containing 26,793 P-type and 36,944 G-type unigenes were shared by the two types of B. vulgaris. 38,397 single nucleotide polymorphisms (SNPs) were found in 9,452 orthologous genes between the P- and G-type plants. We also detected 5,105 simple sequence repeats (SSRs) in the B. vulgaris transcriptome, comprising mono-nucleotide-repeats (2,477; 48.5 %) and triple-nucleotide-repeats (1,590; 31.1 %). Of these, 1,657 SSRs displayed polymorphisms between the P- and G-type. Consequently, 913 SSR primer pairs were designed with a resolution of more than two nucleotides. We randomly chose 30 SSRs to detect the genetic diversity of 32 Barbarea germplasms. The distance tree showed that these accessions were clearly divided into groups, with the G-type grouping with available Western and Central European B. vulgaris accessions in contrast to the P-type accession, B. stricta and B. verna.

CONCLUSIONS

These data represent useful information for pest-resistance gene mining and for the investigation of the molecular basis of plant-pest interactions.

Multiple hydroxyphenethyl glucosinolate isomers and their tandem mass spectrometric distinction in a geographically structured polymorphism in the crucifer Barbarea vulgaris

Two distinct glucosinolate (GSL) chemotypes (P and G-types) of Barbarea vulgaris (Brassicaceae) were known from southern Scandinavia, but whether the types were consistent in a wider geographic area was not known. Populations (26) from Eastern and Central Europe were analyzed for GSLs in order to investigate whether the two types were consistent in this area. Most (21) could be attributed to one of the previously described GSL profile types, the P-type (13 populations) and the G-type (8 populations), based on differences in the stereochemistry of 2-hydroxylation, presence or absence of phenolic glucobarbarin derivatives, and qualitative differences in indole GSL decoration (tested for a subset of 8+6 populations only). The distinction agreed with previous molecular genetic analysis of the same individuals. Geographically, the P-type typically occurred in Eastern Europe while the G-type mainly occurred in Central Europe. Of the remaining five populations, minor deviations were observed in some individuals from two populations genetically assigned to the G-type, and a hybrid population from Finland contained an additional dihydroxyphenethyl GSL isomer attributed to a combinatorial effect of P-type and G-type genes. Major exceptions to the typical GSL profiles were observed in two populations: (1) A G-type population from Slovenia deviated by a high frequency of a known variant in glucobarbarin biosynthesis ('NAS form') co-occurring with usual G-type individuals. (2) A population from Caucasus exhibited a highly deviating GSL profile dominated by p-hydroxyphenethyl GSL that was insignificant in other accessions, as well as two GSLs investigated by NMR, m-hydroxyphenethylGSL and a partially identified m,p disubstituted hydroxy-methoxy derivative of phenethylGSL. Tandem HPLC-MS of seven NMR-identified desulfoGSLs was carried out and interpreted for increased certainty in peak identification and as a tool for partial structure elucidation. The distinct, geographically separated chemotypes and rare variants are discussed in relation to future taxonomic revision and the genetics and ecology of GSLs in B. vulgaris.

Transcriptome analysis of Gossypium hirsutum flower buds infested by cotton boll weevil (Anthonomus grandis) larvae

Background

Cotton is a major fibre crop grown worldwide that suffers extensive damage from chewing insects, including the cotton boll weevil larvae (Anthonomus grandis). Transcriptome analysis was performed to understand the molecular interactions between Gossypium hirsutum L. and cotton boll weevil larvae. The Illumina HiSeq 2000 platform was used to sequence the transcriptome of cotton flower buds infested with boll weevil larvae.

Results

The analysis generated a total of 327,489,418 sequence reads that were aligned to the G. hirsutum reference transcriptome. The total number of expressed genes was over 21,697 per sample with an average length of 1,063 bp. The DEGseq analysis identified 443 differentially expressed genes (DEG) in cotton flower buds infected with boll weevil larvae. Among them, 402 (90.7%) were up-regulated, 41 (9.3%) were down-regulated and 432 (97.5%) were identified as orthologues of A. thaliana genes using Blastx. Mapman analysis of DEG indicated that many genes were involved in the biotic stress response spanning a range of functions, from a gene encoding a receptor-like kinase to genes involved in triggering defensive responses such as MAPK, transcription factors (WRKY and ERF) and signalling by ethylene (ET) and jasmonic acid (JA) hormones. Furthermore, the spatial expression pattern of 32 of the genes responsive to boll weevil larvae feeding was determined by “in situ” qPCR analysis from RNA isolated from two flower structures, the stamen and the carpel, by laser microdissection (LMD).

Conclusion

A large number of cotton transcripts were significantly altered upon infestation by larvae. Among the changes in gene expression, we highlighted the transcription of receptors/sensors that recognise chitin or insect oral secretions; the altered regulation of transcripts encoding enzymes related to kinase cascades, transcription factors, Ca²⁺ influxes, and reactive oxygen species; and the modulation of transcripts encoding enzymes from phytohormone signalling pathways. These data will aid in the selection of target genes to genetically engineer cotton to control the cotton boll weevil.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-854) contains supplementary material, which is available to authorized users.

High-throughput analysis and characterization of Astragalus membranaceus transcriptome using 454 GS FLX

Astragalus membranaceus (Fisch.) Bge (AR), one of the most important medicinal plants in Asia, was found to exhibit various bioactivities. Due to limited genomic and transcriptomic data, the biosynthetic pathway of the major bioactive compound in AR, is currently unclear. In this study, 454 GS FLX technology was employed to produce a substantial expressed sequence tag (EST) dataset from the AR. In all, 742721 high-quality reads from the AR were produced using Roche GS FLX Titanium. A total of 9893 unique sequences were obtained and annotated by a similarity search against the public databases, and involved in the secondary metabolic pathway, which would facilitate deciphering the molecular mechanism of secondary metabolism in AR. The assembled sequences were annotated with gene names and Gene Ontology (GO) terms. GO revealed the unique sequences that could be assigned to 34 vocabularies. In the KEGG mapping, unique sequences were established as associated with 46 biochemical pathways. These results provided the largest EST collections in AR and will contribute to biosynthetic and biochemical studies that lead to drug improvement. With respect to the genes related to metabolism and biosynthesis pathway were also found. Our work demonstrated the utility of 454 GS FLX as a method for the rapid and cost-effective identification of AR transcriptome, and this EST dataset will be a powerful resource for further studies such as taxonomy, molecular breeding, and secondary metabolism in AR.

Acylated glucosinolates with diverse acyl groups investigated by high resolution mass spectrometry and infrared multiphoton dissociation

With the aim of developing a procedure for detecting and identifying intact acylated glucosinolates (a-GLSs) found in trace quantities in natural plant samples, extracts of Barbarea vulgaris seeds were analyzed by reversed-phase liquid chromatography coupled with electrospray ionization and Fourier-transform ion cyclotron resonance mass spectrometry (RPLC-ESI FTICR MS). After a preliminary optimization of fragmentation conditions, based on a non-acylated parent glucosinolate (glucobarbarin) and three previously identified a-GLSs (the 6'-isoferuloyl esters of glucobarbarin, gluconasturtiin and glucobrassicin), infrared multiphoton dissociation (IRMPD) was employed for a tandem MS-based elucidation of the molecular structures of novel a-GLSs. As a result, three acylated derivatives of glucobarbarin, esterified at the thioglucose moiety with a coumaric acid isomer, sinapic acid or an isomer and a dimethoxycinnamic acid isomer, were identified. In addition, a further acylated glucosinolate was tentatively identified as the isoferuloyl ester of an unidentified hydroxylic derivative of glucobarbarin. This is the first demonstration of diversity in the acyl moieties of thioglucose-acylated glucosinolates, which may reflect the substrate specificity of the endogenous acyl transferase. As expected, 6'-isoferuloyl-glucobarbarin was detected as the main acylated GLS in extracts of B. vulgaris seeds. A quantitative estimate suggested that non-isoferuloyl substituted glucobarbarins correspond to ca. 0.026% of the level of 6'-isoferuloyl glucobarbarin. The formation of an uncommon distonic radical anion, most likely generated in the gas phase upon methyl radical (CH3·) loss from the isoferuloyl anion, is demonstrated.