Genomic Assessment of the Contribution of the Wolbachia Endosymbiont of Eurosta solidaginis to Gall Induction

We explored the genome of the Wolbachia strain, wEsol, symbiotic with the plant-gall-inducing fly Eurosta solidaginis with the goal of determining if wEsol contributes to gall induction by its insect host. Gall induction by insects has been hypothesized to involve the secretion of the phytohormones cytokinin and auxin and/or proteinaceous effectors to stimulate cell division and growth in the host plant. We sequenced the metagenome of E. solidaginis and wEsol and assembled and annotated the genome of wEsol. The wEsol genome has an assembled length of 1.66 Mbp and contains 1878 protein-coding genes. The wEsol genome is replete with proteins encoded by mobile genetic elements and shows evidence of seven different prophages. We also detected evidence of multiple small insertions of wEsol genes into the genome of the host insect. Our characterization of the genome of wEsol indicates that it is compromised in the synthesis of dimethylallyl pyrophosphate (DMAPP) and S-adenosyl L-methionine (SAM), which are precursors required for the synthesis of cytokinins and methylthiolated cytokinins. wEsol is also incapable of synthesizing tryptophan, and its genome contains no enzymes in any of the known pathways for the synthesis of indole-3-acetic acid (IAA) from tryptophan. wEsol must steal DMAPP and L-methionine from its host and therefore is unlikely to provide cytokinin and auxin to its insect host for use in gall induction. Furthermore, in spite of its large repertoire of predicted Type IV secreted effector proteins, these effectors are more likely to contribute to the acquisition of nutrients and the manipulation of the host's cellular environment to contribute to growth and reproduction of wEsol than to aid E. solidaginis in manipulating its host plant. Combined with earlier work that shows that wEsol is absent from the salivary glands of E. solidaginis, our results suggest that wEsol does not contribute to gall induction by its host.

A multimodal analytical method to simultaneously determine monoacetyldiacylglycerols, medium and long chain triglycerides in biological samples during routine lipidomics

Background

Monoacetyldiglycerides (MAcDG), are acetylated triglycerides (TG) and an emerging class of bioactive or functional lipid with promising nutritional, medical, and industrial applications. A major challenge exists when analyzing MAcDG from other subclasses of TG in biological matrices, limiting knowledge on their applications and metabolism.

Methods

Herein a multimodal analytical method for resolution, identification, and quantitation of MAcDG in biological samples was demonstrated based on thin layer chromatography-flame ionization detection complimentary with C30-reversed phase liquid chromatography-high resolution accurate mass tandem mass spectrometry. This method was then applied to determine the MAcDG molecular species composition and quantity in E. solidaginis larvae. The statistical method for analysis of TG subclass composition and molecular species composition of E. solidaginis larvae was one-way analysis of variance (ANOVA).

Results

The findings suggest that the proposed analytical method could simultaneously provide a fast, accurate, sensitive, high throughput analysis of MAcDG from other TG subclasses, including the fatty acids, isomers, and molecular species composition.

Conclusion

This method would allow for MAcDG to be included during routine lipidomics analysis of biological samples and will have broad interests and applications in the scientific community in areas such as nutrition, climate change, medicine and biofuel innovations.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12944-022-01650-w.

Regulation of an important glycolytic enzyme, pyruvate kinase, through phosphorylation in the larvae of a species of freeze-tolerant insect, Eurosta solidaginis

Larvae of the goldenrod gall fly, Eurosta solidaginis, rely on a freeze tolerance strategy to survive the sub-zero temperatures of Canadian winter. Critical to their survival is the accumulation of polyol cryoprotectants and global metabolic rate depression, both of which require the regulation of glycolysis and reorganization of carbohydrate metabolism. This study explored the role that pyruvate kinase (PK) regulation plays in this metabolic reorganization. PK was purified from control (5 °C-acclimated) and frozen (-15 °C-acclimated) larvae and enzyme kinetic properties, structural stability, and post-translational modifications were examined in both enzyme forms. The K_m phosphoenolpyruvate (PEP) of frozen PK was 20% higher than that of control PK, whereas the V_max of frozen PK was up to 50% lower than that of control PK at the lowest assay temperature, suggesting inhibition of the enzyme during the winter. Additionally, the activity and substrate affinity of both forms of PK decreased significantly at low assay temperatures, and both forms were regulated allosterically by a number of metabolites. Pro-Q™ Diamond phosphoprotein staining and immunoblotting experiments demonstrated significantly higher threonine phosphorylation of PK from frozen animals while acetylation and methylation levels remained constant. Together, these results indicate that PK exists in two structurally distinct forms in E. solidaginis. In response to conditions mimicking the transition to winter, PK appears to be regulated to support metabolic rate depression, the accumulation of polyol cryoprotectants, and the need for extended periods of anaerobic carbohydrate metabolism to allow the animal to survive whole-body freezing.

Ecological and Evolutionary Stochasticity Shape Natural Selection

The distribution of biodiversity depends on the combined and interactive effects of ecological and evolutionary processes. The joint contribution of these processes has focused almost exclusively on deterministic effects, even though mechanisms that increase the importance of random ecological processes are expected to also increase the importance of random evolutionary processes. Here we manipulate the sizes of old field fragments to generate correlated sampling effects for a focal population (a gall maker) and its enemy community. Traits and communities were more variable in smaller patches. However, because of the preference of some enemies for some trait values (gall sizes), random variation in population mean trait values exacerbated differences in community composition. The random distribution of traits and interactions created predictable but highly variable patterns of natural selection. Our study highlights how stochastic processes can affect ecological and evolutionary processes structuring the strength and direction of selection locally and at larger scales.

Glucose-6-phosphate dehydrogenase is posttranslationally regulated in the larvae of the freeze-tolerant gall fly, Eurosta solidaginis, in response to freezing

The freeze-tolerant larvae of the goldenrod gall fly (Eurosta solidaginis) undergo substantial alterations to their molecular physiology during the winter including the production of elevated quantities of glycerol and sorbitol, which function as cryoprotectants to survive whole body freezing. Production of these cryoprotectants depends on cytosolic pools of nicotinamide adenine dinucleotide phosphate H (NADPH), a major source being the pentose phosphate pathway (PPP). Glucose-6-phosphate dehydrogenase (G6PDH) mediates the rate-limiting and committed step of the PPP and therefore its molecular properties were explored in larvae sampled from control versus frozen states. G6PDH was purified from control (5°C) and frozen (-15°C) E. solidaginis larvae by a single-step chromatography method utilizing 2',5'-ADP agarose and analyzed to determine its enzymatic parameters. Studies revealed a decrease in K_m for G6P in the frozen animals (to 50% of control values) suggesting an increased flux through the PPP. Immunoblotting of the purified enzyme showed differences in the relative extent of several posttranslational modifications, notably ubiquitination (95% decrease in frozen larvae), cysteine nitrosylation (61% decrease), threonine (4.1 fold increase), and serine phosphorylation (59% decrease). Together these data suggested that the increased flux through the PPP needed to generate NADPH for cryoprotectants synthesis is regulated, at least in part, through posttranslational alterations of G6PDH.

Host-parasitoid evolution in a metacommunity

Patch size and isolation are predicted to alter both species diversity and evolution; yet, there are few empirical examples of eco-evolutionary feedback in metacommunities. We tested three hypotheses about eco-evolutionary feedback in a gall-forming fly, Eurosta solidaginis and two of its natural enemies that select for opposite traits: (i) specialization and poor dispersal ability constrain a subset of natural enemies from occupying small and isolated patches, (ii) this constraint alters selection on the gall fly, causing phenotypic shifts towards traits resistant to generalist and dispersive enemies in small and isolated patches, and (iii) reduced dispersal evolves in small, isolated populations. We sampled patches in a natural metacommunity and found support for all hypotheses; Eurosta's specialist wasp parasitoid attacked fewer galls in small and isolated patches, generating a selection gradient that favoured small galls resistant to predation by a dispersive and generalist bird predator. Phenotype distributions matched this selection gradient, and these phenotypic differences were maintained in a common garden experiment. Finally, we found lower dispersal abilities in small and isolated patches, a phenotypic shift that aids in the maintenance of local adaptation. We suggest that the trophic rank and the species traits of consumers are central to evolution in metacommunities.

HYBRIDIZATION STUDIES ON THE HOST RACES OF EUROSTA SOLIDAGINIS: IMPLICATIONS FOR SYMPATRIC SPECIATION

We studied the inheritance of survival ability in host-associated populations of the tephritid fly, Eurosta solidaginis, to test predictions of sympatric speciation models. Eurosta solidaginis induces galls on two species of goldenrod, Solidago altissima and S. gigantea. The host-associated populations have been hypothesized to be host races that originated in sympatry (Craig et al. 1993). We found evidence for disruptive selection for host use, which is a critical assumption of sympatric speciation models. Each host race had higher survival rates on their host plant than on the alternative host. F₁ and backcross hybrids also had lower survival rates than the pure host-race flies on their host plant. Since assortative mating occurs due to host-plant preference (Craig et al. 1993) this would select for divergence in host preference. Low hybrid survival could have been due to strong genetic incompatibilities of the populations or due to host adaptation by each population. Strong genetic incompatibilities would result in poor survival on all host plants, while host adaptation could result in low overall survival with high hybrid survival on some host plants with particularly "benign" environments. High survival of F₁ , F₂ , and backcross hybrids on some plant genotypes in some years supported the host adaptation hypothesis. F₁ flies mated and oviposited normally and produced viable F₂ and backcross hybrids indicating gene flow is possible between the host races. A few flies developed and emerged on the alternative host plant. This demonstrates that genes necessary to utilize the alternative host exist in both host races. This could have facilitated the origin of one of the populations via a host shift from the ancestral host. The inheritance of survival ability appears to be an autosomal trait. We did not find evidence that survival ability was maternally influenced or sex linked.

Expression of miRNAs in response to freezing and anoxia stresses in the freeze tolerant fly Eurosta solidaginis

Insect cold hardiness is associated with substantial metabolic rate suppression, often including developmental diapause as well as metabolic suppression imposed by freezing and freeze-associated oxygen limitation. MicroRNAs, small non-coding transcripts that bind to mRNA, are known modulators of hypometabolism in freeze tolerant insects. To further contribute to the growing signature of stress-responsive miRNAs, this study amplified and quantified changes in the expression levels of four microRNA species, miR-8, miR-9, miR-92b and miR-277, in response to freezing or anoxia exposures of freeze tolerant gall fly larvae, Eurosta solidaginis. MiR-92b levels were significantly elevated by 1.57-fold in frozen E. solidaginis at -15°C as compared with 5°C controls, whereas miR-92b levels were significantly reduced in anoxic E. solidaginis to levels that were 0.77-fold as compared with larvae held under normoxic conditions. The other miRNAs investigated showed no significant changes in stressed larvae. These data demonstrate differential miR-92b expression in frozen/anoxic versus control insect larvae and position this miRNA as a stress responsive marker in this model insect.

The limits of drought-induced rapid cold-hardening: extremely brief, mild desiccation triggers enhanced freeze-tolerance in Eurosta solidaginis larvae

Rapid cold-hardening (RCH) is a highly conserved response in insects that induces physiological changes within minutes to hours of exposure to low temperature and provides protection from chilling injury. Recently, a similar response, termed drought-induced RCH, was described following as little as 6h of desiccation, producing a loss of less than 10% of fresh mass. In this study, we investigated the limits and mechanisms of this response in larvae of the goldenrod gall fly Eurosta solidaginis (Diptera, Tephritidae). The cold-hardiness of larvae increased markedly after as few as 2h of desiccation and a loss of less than 1% fresh mass, as organismal survival increased from 8% to 41% following exposure to -18 °C. Tissue-level effects of desiccation were observed within 1h, as 87% of midgut cells from desiccated larvae remained viable following freezing compared to 57% of controls. We also demonstrated that drought-induced RCH occurs independently of neuroendocrine input, as midgut tissue desiccated ex vivo displayed improved freeze-tolerance relative to control tissue (78-11% survival, respectively). Finally, though there was an increase in hemolymph osmolality beyond the expected effects of the osmo-concentration of solutes during dehydration, we determined that this increase was not due to the synthesis of glycerol, glucose, sorbitol, or trehalose. Our results indicate that E. solidaginis larvae are extremely sensitive to desiccation, which is a triggering mechanism for one or more physiological pathways that confer enhanced freeze-tolerance.

Seasonal accumulation of acetylated triacylglycerols by a freeze-tolerant insect

Most animals store energy as long-chain triacylglycerols (lcTAGs). Trace amounts of acetylated triacylglycerols (acTAGs) have been reported in animals, but are not accumulated, likely because they have lower energy density than lcTAGs. Here we report that acTAGs comprise 36% of the neutral lipid pool of overwintering prepupae of the goldenrod gall fly, Eurosta solidaginis, while only 17% of the neutral lipid pool is made up of typical lcTAGs. These high concentrations of acTAGs, present only during winter, appear to be synthesized by E. solidaginis and are not found in other freeze-tolerant insects, nor in the plant host. The mixture of acTAGs found in E. solidaginis has a significantly lower melting point than equivalent lcTAGs, and thus remains liquid at temperatures at which E. solidaginis is frozen in the field, and depresses the melting point of aqueous solutions in a manner unusual for neutral lipids. We note that accumulation of acTAGs coincides with preparation for overwintering and the seasonal acquisition of freeze tolerance. This is the first observation of accumulation of acTAGs by an animal, and the first evidence of dynamic interconversion between acTAGs and lcTAGs during development and in response to stress.

Cytokinins in the ball gall of Solidago altissima and in the gall forming larvae of Eurosta solidaginis

•  The identities and concentrations of cytokinins in goldenrod ball galls, normal goldenrod stems, and in the gall-forming larvae of the tephritid fly Eurosta solidaginis were determined in order to gain insight into the mechanism of ball gall formation on Solidago altissima. •  Zeatin, zeatin riboside, isopentenyladenine, and isopentenyladenosine were identified and quantified by gas chromatography-mass spectrometry in ball galls, goldenrod stems, and in larvae of E. solidaginis. •  Concentrations of the four cytokinins were higher in gall tissues than in stem tissues when expressed on a weight per stem length basis but not when expressed on a weight per weight basis. Isopentenyladenine was the most abundant cytokinin in larvae from developing galls and in larvae from fully formed galls exhibiting 'green islands'. In first instar larvae, isopentenyladenine was present at a much higher concentration (c. 50 fold) than in control stem tissues. •  The presence of an E. solidaginis larva results in higher cytokinin amounts in a given length of stem than are found in its absence. Larvae of E. solidaginis may act as point sources of isopentenyladenine in developing ball galls. ABBREVIATIONS: BHT, butylated hydroxytoluene; DMSO, dimethylsulphoxide; DMSO^- , dimethylsulphoxide anion; GC-MS, gas chromatography mass spectrometry; HPLC, high pressure liquid chromatography; iP, isopentenyladenine; iPA, isopentenyladenosine; Me-dHZ, methyl-dihydrozeatin; Me-dHZR, methyl-dihydrozeatin riboside; Me-iP, methyl-isopentenyladenine; Me-iPA, methyl-isopentenyladenosine; Me-Z, methyl-zeatin; Me-ZOG, methyl-zeatin-O-glucoside; Me-ZR, methyl-zeatin riboside; Me-ZROG, methyl-zeatin riboside-O-glucoside; m/z, mass to charge ratio; PVPP, polyvinylpolypyrrolidone; SIM, selected ion monitoring; TEAB, triethylammonium bicarbonate; Z, zeatin; ZR, zeatin riboside.

Indole-3-acetic acid and ball gall development on Solidago altissima

•  Indole-3-acetic acid (IAA) concentrations are reported in ball galls, stem tissue, and in the gall-forming larvae of the tephritid fly Eurosta solidaginis, to gain insight into the mechanism of ball gall formation on the goldenrod, Solidago altissima. •  IAA was identified and quantified by gas chromatography-mass spectrometry, with [¹³ C₆ ]-IAA as an internal standard. •  Concentrations of IAA were higher in gall tissues than in stem tissues on a weight per stem length basis, but not on a weight per weight basis. First instar larvae of E. solidaginis contained IAA at 9 µg g^-1 f. wt or 1.21 ng per larva, a concentration 33-times higher than in control stem tissues. Gall development was unaffected by the early removal of the apical bud and leaves from galls, and lateral bud release was inhibited by galls in terminal positions in the absence of the apical bud, indicating that ball galls act as sources of auxin. •  The presence of an E. solidaginis larva results in higher amounts of IAA in a given length of stem than are found in its absence. E. solidaginis larvae may act as point sources of IAA in developing ball galls.

Variation in selection pressures on the goldenrod gall fly and the competitive interactions of its natural enemies

Larvae of the tephritid fly Eurosta solidaginis induce ball-shaped galls on the stem of tall goldenrod, Solidago altissima. Survival probability depends on gall size; in small galls the larva is vulnerable to parasitoid oviposition, whereas larvae in large galls are more frequently eaten by avian predators. Fly populations from 20 natural old fields in central Pennsylvania were monitored in 1983 and 1984 to examine the distribution of the selection intensity imposed by natural enemies, the parasitoids Eurytoma gigantea and E. obtusiventris, the inquiline Mordellistena unicolor, and the predatory birds Dendrocopus pubescens and Parus atricapillus. Mordellistena and E. obtusiventris are able to attack galls of all diameters while E. gigantea and the predatory birds preferentially assaulted small and large diameter galls, respectively. Eurosta in intermediate sized galls had the highest survivorship, hence selection had a stabilizing component. However, parasitoid attack was more frequent than bird attack, and the two did not exactly balance, thus there was also a directional component. The mean directional selection intensity on gall size was 0.21 standard deviations of the mean, indicating that larger gall size was favored. Interactions among the insect members of the Eurosta natural enemy guild are complex and frequent.