Odoriferous Defensive Stink Gland Transcriptome to Identify Novel Genes Necessary for Quinone Synthesis in the Red Flour Beetle, <i>Tribolium castaneum</i>

<div><p>Chemical defense is one of the most important traits, which endow insects the ability to conquer a most diverse set of ecological environments. Chemical secretions are used for defense against anything from vertebrate or invertebrate predators to prokaryotic or eukaryotic parasites or food competitors. Tenebrionid beetles are especially prolific in this category, producing several varieties of substituted benzoquinone compounds. In order to get a better understanding of the genetic and molecular basis of defensive secretions, we performed RNA sequencing in a newly emerging insect model, the red flour beetle <i>Tribolium castaneum</i> (Coleoptera: Tenebrionidae). To detect genes that are highly and specifically expressed in the odoriferous gland tissues that secret defensive chemical compounds, we compared them to a control tissue, the anterior abdomen. 511 genes were identified in different subtraction groups. Of these, 77 genes were functionally analyzed by RNA interference (RNAi) to recognize induced gland alterations morphologically or changes in gland volatiles by gas chromatography-mass spectrometry. 29 genes (38%) presented strong visible phenotypes, while 67 genes (87%) showed alterations of at least one gland content. Three of these genes showing quinone-less (ql) phenotypes – <i>Tcas-ql VTGl</i>; <i>Tcas-ql ARSB</i>; <i>Tcas-ql MRP</i> – were isolated, molecularly characterized, their expression identified in both types of the secretory glandular cells, and their function determined by quantification of all main components after RNAi. In addition, microbe inhibition assays revealed that a quinone-free status is unable to impede bacterial or fungal growth. Phylogenetic analyses of these three genes indicate that they have evolved independently and specifically for chemical defense in beetles.</p></div

Annotated quinone synthesis-related genes and their relative gland transcriptome expression levels.

<p>A, glucosidases, 19 genes; B, phenol oxidases, 14 genes; C, peroxidases, 18 genes were annotated. In all charts, along the X-axes the different genes are presented, while the Y-axes present log2[fold change of reads in glands against control]. Abbreviations: FC: fold change; s1: sample 1, anterior abdomen as a control; s3: sample 3, male prothoracic glands; s4: sample 4, female prothoracic glands; s5: sample 5, male abdominal glands; s6: sample 6, female abdominal glands; s2: sample 2, <i>tar</i> prothoracic glands. The genes with higher reads in gland samples are marked with red squares.</p

Relative transcriptomic gland expression levels of the <i>Tribolium</i> homologs of the three novel <i>quinone-less</i> genes.

<p>A, <i>Tcas-ql VTGl</i> (<i>GT39</i>) homologs. B, <i>Tcas-ql ARSB</i> (<i>GT62</i>) homologs. C, <i>Tcas-ql MRP</i> (<i>GT63</i>) homologs. In all charts, along the X-axes the different genes are presented (using the same GI numbers and expression level square codes as in <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003596#pgen-1003596-g007" target="_blank">Figure 7</a>), while Y-axes present log2[fold change of reads in glands against control]. For abbreviations see <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003596#pgen-1003596-g003" target="_blank">Figure 3</a>.</p

Quantification of main volatile glandular chemicals by GC-MS in wild-type and novel <i>quinone-less</i> gene RNAi-knock-downs.

<p>Comparisons in female thoracic glands (ft) (A), female abdominal glands (fa) (B), male thoracic glands (mt) (C), and male abdominal glands (ma) (D). Y-axis: amount in nanogram; X-axis: wild-type and different RNAi-knock-downs. Abbreviations: wt: wild-type; E: dsEGFP-injected control; GT39: <i>Tcas-ql VTGl</i>; GT62: <i>Tcas-ql ARSB</i>; GT63: <i>Tcas-ql MRP</i>; MBQ-sum: methyl-1,4-benzoquinone; EBQ-sum: ethyl-1,4-benzoquinone; 15ene: 1-pentadecene; 17diene: 1,8-heptadecadiene; 17ene: 1-heptadecene. The error bars indicate standard deviations at N = 15–30.</p

Quantification of the main volatiles in wild-type odoriferous stink glands.

<p>The amounts are indicated as mean±standard deviation. N = 15–30. Abbreviations: fem.: female; thr: thoracic glands; abd: abdominal glands.</p

Phenotype classifications of 77 highly gland-specifically expressed genes by RNAi.

<p>A, Description of phenotype strengths: 1, Very strong: at least one type of chemical was undetectable or less than 5% left in thoracic or/and abdominal glands; 2, Strong: at least one type of chemical was 75%–95% reduced or increased by more than 75% in thoracic or/and abdominal glands; 3, Strong or very strong but occasionally: phenotype was similar to 1 or 2 but not observed in all the injected beetles; 4, Some changes: at least one type of chemical was 25%–75% reduced or increased in thoracic or/and abdominal glands; 5, Neglectable changes: less than 25% reduction or increase in any type of chemical; 6, uncloned or untested. B, Strong and very strong gland phenotypes (strengths 1–3 in panel A) in details, some genes had more than one phenotype. Except for the last two classifications, examples are provided in <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003596#pgen-1003596-g004" target="_blank">Figure 4</a>: localized brown or melanized dots (4B1-B5 and G2-G5); condensed (4C1-C2); colorless (4E1-F2); empty-looking (4D3).</p

Phenol oxidase (PO) activity assays of wild-type and novel <i>quinone-less</i> gene RNAi knock-downs in females.

<p>The Y-axis indicates the square root of PO Vmax, red boxes are boxplots, green lines represent the mean value, the gray line represents the grand mean, while the X-axis presents wild-type, control injections, and different RNAi-knock-downs (N = 12–25). Buffer: buffer-injection control; EGFP: dsEGFP-injection control; GT39: <i>Tcas-ql VTGl</i>; GT62: <i>Tcas-ql ARSB</i>; GT63: <i>Tcas-ql MRP</i>. The asterisks (*) mark the t-test results comparing to wild-type: ***: p<0.001; **: 0.001</p

Microbe growth inhibition assays of wild-type and RNAi-knock-down glands.

<p>The fungus <i>A. niger</i> (A) and the bacterium <i>A. globiformis</i> (B) were analyzed for gland-mediated growth inhibition. Y-axes indicate the areas of respective inhibition zones (cm²). X-axes: sex-specific wild-type (wt) and different RNAi-knock-downs (m: male; fm: female; GT12: gene causing alkene-less phenotype; GT63: <i>Tcas-ql MRP</i>). Non-parametric comparisons were made between wild-type and knock-downs using Wilcoxon method, ***, P<0.001; **, 0.001</p

GO annotation of odoriferous glands transcriptome data.

<p>A, analyses of the genes abundant in control (Ctl; 1451 genes abundant in anterior abdomen transcriptome), wild-type glands (Glands; 1206 genes), and the genes identified in all ten subtraction groups (<a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003596#pgen-1003596-g001" target="_blank">Figure 1</a>) together (G1–G10; 511 genes). B, analyses of individual Group 1 to Group 10, respectively and the control (Ctl; 290 genes possessing at least 64 times higher reads in anterior abdomen than the wild-type gland samples). X-axis: different GO terms (level 2); Y-axis: percentage of the genes classified in each group.</p

Visible morphological gland phenotypes after RNAi.

<p>A1–A4, wild-type: A1, thr (darker regions on two sides); A2, abd (darker regions on the tip); A3, thrg reservoir filled with yellowish oily fluid ; A4, abdg; B1–B5, <i>GT20</i> knock-down: B1, thr, visible melanized dots on both sides; B2, abd, darker gland regions; B3, thrg, left one, small localized brown dots, right one, tiny localized brown dots accompanied by colorless oily secretions; B4, thrg, localized melanized dots; B5, abdg, localized melanized (left) or brown (right) dots; C1–C2, <i>GT47</i> knock-down: C1, thr; C2, thrg, condensed secretions, i.e. not as semitransparent as in wild-type and close to solid; D1–D4, <i>GT62</i> knock-down: D1–D2, thr and abd with the gland regions almost invisible through cuticle; D3, thrg, very few colorless secretions, in some cases, empty-looking (not shown); D4, abdg, colorless secretions; E1–E2, <i>GT39</i> knock-down, thrg (E1) and abdg (E2), colorless secretions; F1–F2, <i>GT63</i> knock-down, thrg (F1) and abdg (F2, ventral view), colorless secretions; G1–G5, <i>GT02</i> knock-down, at stage A10 localized melanized dots (G1, abd; G2, abdg), at stage 24 localized brown dots to different extent (G3 and G4, thrg, G5, abdg); H1–H4, <i>GT25</i> knock-down, various phenotypes: H1 and H2, thrg, tiny (H1) or big (H2) localized brown dots; H3, left half abdg, oily colorless secretions; H4, abdg, localized brown dots. The beetles used for G3–G5 were from stage A24, all the others were at stage A10. The statistics of the visible phenotypes are in <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003596#pgen-1003596-g005" target="_blank">Figure 5B</a>. Abbreviations used in this legend: stage A10, adults ten days; A24, twenty four days after eclosure; abd, abdominal tip, ventral view; abdg, dissected abdominal tip with abdominal glands, dorsal view; thr, prothorax, ventral view; thrg, dissected prothoracic gland(s), with the opening on top, the secretory cells are attached to the reservoir.</p