Survival and development of potato psyllid (Hemiptera: Triozidae) on Convolvulaceae: Effects of a plant-fungus symbiosis (Periglandula)

Plant species in the family Solanaceae are the usual hosts of potato psyllid, Bactericera cockerelli (Šulc) (Hemiptera: Psylloidea: Triozidae). However, the psyllid has also been shown to develop on some species of Convolvulaceae (bindweeds and morning glories). Developmental success on Convolvulaceae is surprising given the rarity of psyllid species worldwide associated with this plant family. We assayed 14 species of Convolvulaceae across four genera (Convolvulus, Calystegia, Ipomoea, Turbina) to identify species that allow development of potato psyllid. Two populations of psyllids were assayed (Texas, Washington). The Texas population overlaps extensively with native Convolvulaceae, whereas Washington State is noticeably lacking in Convolvulaceae. Results of assays were overlain on a phylogenetic analysis of plant species to examine whether Convolvulaceae distantly related to the typical host (potato) were less likely to allow development than species of Convolvulaceae more closely related. Survival was independent of psyllid population and location of the plant species on our phylogenetic tree. We then examined whether presence of a fungal symbiont of Convolvulaceae (Periglandula spp.) affected psyllid survival. These fungi associate with Convolvulaceae and produce a class of mycotoxins (ergot alkaloids) that may confer protection against plant-feeding arthropods. Periglandula was found in 11 of our 14 species, including in two genera (Convolvulus, Calystegia) not previously known to host the symbiont. Of these 11 species, leaf tissues from five contained large quantities of two classes of ergot alkaloids (clavines, amides of lysergic acid) when evaluated by LC-MS/MS. All five species also harbored Periglandula. No ergot alkaloids were detected in species free of the fungal symbiont. Potato psyllid rapidly died on the five species that harbored Periglandula and contained ergot alkaloids, but survived to adulthood on seven of the nine species in which ergot alkaloids were not detected. These results support the hypothesis that a plant-fungus symbiotic relationship affects the suitability of certain Convolvulaceae to potato psyllid.


Introduction
The potato psyllid, Bactericera cockerelli (Šulc) (Hemiptera: Psylloidea: Triozidae) is a pest of solanaceous crops such as potatoes, tomatoes, and peppers. The psyllid occurs throughout the western and central United States, Canada, Mexico, and Central America [1], and as an introduction in New Zealand and Australia [2,3]. High densities of the psyllid may lead to plant disorders known as "psyllid yellows" [4,5] caused by a toxin that is injected into plants during the psyllid's feeding activities [6]. However, recent crop losses have been caused primarily by a bacterial pathogen, 'Candidatus Liberibacter solanacearum' (Lso), that is transmitted by the psyllid [1]. Difficulties in managing potato psyllid and its associated Liberibacter are in part due to poor understanding of the role that non-crop species have in the biology of the vector. Most species of psyllids are monophagous or oligophagous, limited to development on plants within a single genus or family [7,8]. Potato psyllid is unusual in being able to develop on plants across more than a single family [9,10,11,12]. Non-crop plant species act as reservoirs of the insect during the growing season and may help the psyllid bridge intervals in which crop hosts are unavailable [10,13,14,15,16,17]. It is therefore important to know what noncrop species of plants found in potato or tomato growing regions also support the reproduction and development of potato psyllid.
Although plant species in the family Solanaceae (Solanales) are the typical developmental hosts for potato psyllid, at least some species in the Convolvulaceae (Solanales) also support development [9,10,11,12,14,18]. Observations leading to this conclusion include rearing trials [9,10,11,12] and field records [14,18]. Developmental success on Convolvulaceae is unexpected given that Convolvulaceae is substantially underrepresented among plant families as hosts of Psylloidea. Despite its extensive diversity and widespread distribution [19] Convolvulaceae is listed as a developmental host for only five species of psyllids worldwide, including potato psyllid [20]. Rearing trials with potato psyllid have been limited to two species, Convolvulus arvensis L. (field bindweed) and Ipomoea batatas (L.) Lam. (sweet potato). While potato psyllid is able to complete development on these species, development rates are slow and may be accompanied by nymphal mortality [10,12].
In this study, we examined the development of potato psyllid on species and genera of Convolvulaceae that have not previously been assayed. Our assays targeted species that are native to North America and are thus likely to have an evolutionary history with at least some populations of potato psyllid. Our first objective was to assay a taxonomically broader group of Convolvulaceae than previously done, to determine whether plant suitability extends beyond Co. arvensis and I. batatas. Part of this objective included a comparison of two haplotypes of the psyllid on each plant species. Potato psyllid occurs as a minimum of four unique genetic types or "haplotypes" [21,22] that we now know differ biologically [23,24,25,26,27]. We compared developmental success on Convolvulaceae between two of these haplotypes, the Central haplotype and the Northwestern haplotype. Convolvulaceae is highly diverse in the southern US and Mexico [28] where its presence overlaps extensively with the distribution of the Central haplotype [29]. In contrast, native Convolvulaceae are almost completely absent from the Pacific Northwest region of the US [28] where the Northwestern haplotype of potato psyllid is endemic [21,22,30]. Thus, the Northwestern haplotype is likely to have a much-reduced field history with native Convolvulaceae in comparison to the Central haplotype.
Our second objective was to look for traits that predict whether a given plant species allows psyllid development. We addressed two separate questions in this objective. First, we examined whether suitability is predicted by the location of plant species within a phylogenetic tree. Because of the strong tendency towards host specificity among species of Psylloidea, host switching or dietary expansion by psyllids tends to be phylogenetically conserved [31] such that evolutionary shifts in diets by psyllids are often between closely related plant species [8,31,32]. This specialism prompted us to examine whether plant suitability tracked plant phylogeny. We constructed a phylogenetic tree from DNA-sequence data to examine whether plant species allowing successful development of potato psyllid clustered together in the tree, as would be expected if plant chemistry or other traits affecting psyllid host use also grouped phylogenetically [33].
We then examined whether psyllid development was affected by the presence of a plantfungus mutualism found in Convolvulaceae. The Convolvulaceae is unusual among dicotyledonous plant families in its association with a class of chemicals known as ergot alkaloids [34]. Many species of Convolvulaceae have formed a symbiotic association with clavicipitaceous fungi in the genus Periglandula [35,36,37,38]. This fungus is vertically transmitted, and is present systemically in members of the family Convolvulaceae [39] often forming epiphytic colonies surrounding peltate glandular trichomes on the adaxial leaf surfaces [35,40] where the colonies produce ergot alkaloids [41]. This symbiosis appears to be most common in Ipomoea and related genera, with possibly 450 or more plant species worldwide having the association [34,42]. Similar alkaloids produced by clavicipitaceous fungi in grasses have been shown to have deleterious effects against herbivorous insects [43,44,45]. The defensive properties of ergot alkaloids associated with the Convolvulaceae-Periglandula symbiosis have received almost no attention, although extracts from Ipomoea parasitica (H.B.K.) G. Don, have been found to reduce feeding and digestive efficiency of caterpillars [46]. Our overall goal therefore was to examine whether survival and development of the potato psyllid was correlated with the presence or absence of Periglandula, and to determine whether psyllid development was affected by the types and quantities of fungal alkaloids produced by this symbiosis.

Source of plants and insects
Insect bioassays included a screening of 11 species of native Convolvulaceae distributed across three plant genera (Convolvulus, Ipomoea, and Turbina), and three introduced species in Convolvulus and Calystegia including the widespread pest field bindweed, Co. arvensis (Table 1). All species except Calystegia silvatica overlap geographically with psyllids of the Central haplotype ( Table 1). The Northwestern haplotype overlaps geographically with Co. arvensis, possibly with Ca. silvatica, and is likely to have some overlap with the four species of Ipomoea that are grown extensively as summer ornamentals (Table 1). It is unlikely that these ornamentals are able to survive the winter conditions of the Pacific Northwest.
Test plants were examined in side-by-side comparisons with potato, Solanum tuberosum L. ('Russet Burbank') (Solanaceae), a typical and highly suitable host for potato psyllid. Plants were grown either from seeds or from stem cuttings (sources listed in Table 1). Seeds were scarified using sandpaper and soaked in gibberellic acid (1000 ppm in water) for 24 h prior to planting. Plants were grown in 10-cm pots (volume~473. Potato psyllids to be used in assays were obtained from colonies maintained at the USDA-ARS facility in Wapato, WA. The parental insects for colonies were collected from potato fields near Weslaco, TX in March 2017 (Central haplotype, APHIS permit P526P-17-00366) and from solanaceous weeds growing near Prosser, WA in the summer and autumn of 2016 (Northwestern haplotype). The colonies were maintained on potato ('Russet Burbank') at 22˚C and a 16:8 h light: dark cycle. Colonies were assayed preceding the study using high resolution melting analysis to confirm haplotype status [21]. Colonies were checked periodically for Lso infection using PCR detection methods [50]. Colonies were Lso-free for the duration of the study.

Suitability of Convolvulaceae to potato psyllid
Our primary objective was to determine whether the potato psyllid is able to complete development on targeted plant species. Given the large size of the experimental design (two psyllid haplotypes x 15 plant species), we limited our measures of psyllid performance to two traits: egg-to-adult survival (as a yes/no variable), and egg-to-adult development time (in days). Ten adults (unsexed) from both haplotypes were collected from their respective colony cages. The ten psyllids of a given haplotype were confined for egg-laying on a test plant kept individually in a 7.5 L plastic container (Cambro 1 , Huntington Beach, CA) modified to allow ventilation at 22˚C and a 16:8 h light: dark cycle. Once 20 or more eggs were present on a test plant, the For plant species on which psyllids developed successfully, we recorded the number of days required to develop from egg deposition to production of the first adult (i.e., the minimum time required to complete development). We did not attempt to compare survival rates among plant species due to difficulties in obtaining accurate estimates of egg numbers on plants without also damaging the plant. Once new adults were seen in a container, that plant and container was dismantled. A leaf was collected from the plant for DNA extraction and biochemical analysis (described below).
On species which failed to support development, mortality almost invariably occurred as first instar nymphs often within 48 h of hatch. When monitoring showed that all nymphs on a given plant were dead, the assay for that plant and container was dismantled, and leaf samples were collected for DNA extraction and ergot alkaloid quantification. We had five replicates per plant species per psyllid haplotype combination. The large number of treatments, combined with uneven germination of seed, did not allow us to conduct the five replicates simultaneously. Thus, each replicate was initiated on a separate date, with date of assay included in the statistical analyses as a blocking factor (see Statistical analyses).

Phylogenetic mapping of Convolvulaceae
DNA was extracted using a cetyltrimethylammonium bromide (CTAB) precipitation method [51]. Two different universal plant barcoding primer sets were used. The first primer set targeted approximately 500 bp of the internal transcribed spacer region (ITS): ITS2F (ATGCGA TACTTGGTGTGAAT) and ITS3R (GACGCTTCTCCAGACTACAAT) [52]. The second primer set targeted approximately 684 bp region of the chloroplast maturase K gene (matK): matK 472-F (CCCRTYCATCTGGAAATCTTGGTT) and matK 1248-R (GCTRTRATAATGAGAAAGA TTTCTGC) [53]. PCR conditions used for both primer sets were similar, consisting of an initial denaturation step of 94˚C for 5 min followed by 35 cycles of 94˚C for 30 s, 56˚C for 30 s, and 72˚C for 42 s, followed by a final extension at 72˚C for 10 min. Each 20μl reaction contained Amplitaq Gold 360 PCR Master Mix (Invitrogen, Carsbad, CA), 500nM of each primer, and DNA template (10-20 ng). Upon amplification, bands were excised from agarose gels, purified using GenElute minus ethidium bromide spin columns (Sigma, St. Louis, MO), and were cloned using a TOPO TA cloning kit with TOP10 E. coli chemical competent cells (Invitrogen, Carlsbad, CA). The QIAprep spin mini prep kit (Qiagen, Valencia, CA) was used to prepare plasmid DNA for sequencing by MC Laboratories (MC Lab, San Francisco, CA). Sequences were deposited into GenBank (Table 2).
DNA sequences were aligned and consensus sequences were made using Geneious R10 software (North America Biomatters Inc, Newark, NJ). The phylogenetic tree was constructed using a Tamura-Nei model and neighbor-Joining method with the Tree Builder function of Geneious R10 [54]. Phylogenetic distances for tree construction were estimated based upon concatenated sequences of ITS and matK regions. Potato was treated as an outgroup.

Detection of the Convolvulaceae-Periglandula association
Because Periglandula is not always readily visible on plants even when the fungus is present, extraction and analysis of DNA-sequences is often used to confirm infestation. Presence or absence of the dmaW gene, encoding 4-(γ,γ -dimethylallyl) tryptophan synthase and required for the determinant step of ergot alkaloid synthesis, was evaluated using PCR [35,55]. Plant DNA extracted using CTAB method was used to amplify approximately 1050 bp region using dmaWF5 (GACCGTAAACGAGTCAGGAA) and dmaWR2 (AAATACACCTGGGGCTCG) primers.
PCR conditions consisted of an initial denaturation step of 95˚C for 5 min followed by 40 cycles of 95˚C for 1 min, 52˚C for 1 min, and 72˚C for 45 s, followed by a final extension at 72˚C for 5 min. Each 20μl reaction contained Amplitaq Gold 360 PCR Master Mix (Invitrogen, Carsbad, CA), 500nM of each primer, and DNA template (10-20 ng). Upon amplification, bands were excised from agarose gels, purified using GenElute minus ethidium bromide spin columns, and were cloned (methods described in previously). Sequencing again was done by MC Laboratories. Sequences were deposited into GenBank ( Table 2).
Fully expanded leaves were collected from assayed plants at the end of suitability tests and subjected to air drying at the room temperature~22-25˚C for 3-5d. Dried tissue was ground using either a mortar and pestle or a cyclone sample mill with a 0.5 mm screen (UDY Corporation, Fort Collins CO). Extraction solution (79:20:1 ACN:water:acetic acid) was added to ground sample at a ratio of 4 mL/g and turned for 90 min in the dark [56]. The sample was then centrifuged for 2 min at 1462 x g. Dilution solution (250 μL 20:79:1 ACN:water:acetic acid) was added to 250 μL supernatant, vortexed for 10 sec, then placed in an amber HPLC vial for ergot alkaloid analysis by LC-MS/MS.
An ABI/SCIEX 3200 QTRAP LC-MS/MS system (Applied Biosystems, Foster City, CA USA) was used to monitor for ergoline and ergopeptide compounds via positive electrospray ionization, with separation performed using a Perkin Elmer (Waltham, MA USA) Series 200 autosampler and HPLC connected to a Gemini C18 column (150 x 4.6 mm, 5 μ, Phenomenex (Torrance, CA USA)) with a 4 x 3 mm security guard cartridge of similar packing [56]. Mobile phases consisted of 5 mM ammonium acetate and methanol:water:acetic acid in a ratio of 10:89:1(v/v/v) (A) or 97:2:1 (B) and were run in a gradient program at 1 mL/min. Multiple reaction monitoring (MRM) of two transitions (quantitative and qualitative) per compound was used to detect the ergot alkaloids ergonovine, ergotamine, ergocornine, α-ergocryptine, ergocristine, ergovaline, ergine, ergosine and their epimers, as well as agroclavine, chanoclavine, lysergol, lysergic acid, oxidized luol, dihydrolysergol, chanoclavine, dihydroergosine, dihydroergotamine, festuclavine, fumigaclavine and elymoclavine. The presence of a mycotoxin was confirmed when the signal was equal to or greater than a signal-to-noise (S/N) ratio of 3:1 (limit of detection (LOD)), and both quantitative and qualitative transitions were present. Samples were quantitated blind as to sample identity against a standard curve using Analyst 1.6.2 and MultiQuant 3.0.1 (Applied Biosystems). The limit of quantitation (LOQ) was defined as the concentration at which the analyte had a precision and accuracy that did not exceed greater than 20% of the coefficient of variation [57]. LOD and LOQ for detected mycotoxins were as follows: ergotamine, ergocornine, ergosine ergocristine and agroclavine (1, 1 ng/mL); ergonovine (1, 2 ng/mL); lysergol (2, 2 ng/mL); lysergic acid (20 and 50 ng/mL). No commercial standards were available for chanoclavine, festuclavine, elymoclavine, elymoclavine fructoside, ergine and dihydrolysergol. These compounds were compared on a scale of present ("+" indicating low, "++" indicating high) or not present ("-") amongst the plant species extracted based on relative peak area.

Statistical analyses
Effects of plant species and psyllid haplotype on mean psyllid development time were examined using a generalized linear mixed model (Proc GLIMMIX) [58]. Plant species, psyllid haplotype, and the species x haplotype interaction were included as fixed effects. The analysis was limited to plant species on which psyllids completed development to the adult stage. We specified an underlying gamma distribution using a DIST = gamma statement. This distribution is useful for modeling time-to-occurrence data [59]. The ILINK function was used to back-transform means into the original units (number of days to first adult). The CONTRAST statement was used to examine a priori defined comparisons among plant species following a significant plant species effect in the overall model (see results). The survival data (yes/no) were not analyzed statistically, as the two haplotypes showed identical results as to what plant species supported development (see results).
Genetic distances between species in a tree were calculated automatically by the Geneious R10 software using the Tamura-Nei model. This approach expresses distance as nucleotide substitutions per site. We used these distances to determine whether plant suitability for psyllids decreased as genetic distance from potato increased. We conducted a two-sample t-test [58] to determine whether mean genetic distance from potato differed between plant species allowing development and plant species not allowing development. If suitability was affected by genetic distance from the typical host (potato), we expected mean distance to be smaller for plants on which psyllids survived than plants on which psyllids failed to survive.

Psyllid developmental success and plant phylogeny
Eggs were present within 24-48 h of adding egg-laying psyllids on all plant species except for Ipomoea pandurata and Turbina corymbosa which generally required 72 h before eggs were present. Egglayers were noticeably reluctant to settle on these two plant species, and instead were seen to spend considerable time wandering on the sides of the cages. Psyllids of both haplotypes failed to complete development in all five replicates on Convolvulus equitans, Calystegia silvatica, Ipomoea imperati, Ipomoea leptophylla, Ipomoea pandurata, Ipomoea tricolor, and Turbina corymbosa. Nymphs invariably died within a week of hatch on these species. It was not clear from our observations whether nymphs fed to any extent prior to death. On some species, notably I. imperati, I. pandurata, and T. corymbosa, mortality occurred within 24-48 h of egg hatch. Psyllids of both haplotypes completed development on potato, Convolvulus arvensis, Convolvulus tricolor, Ipomoea alba, Ipomoea cordatotriloba, Ipomoea hederacea, Ipomoea ternifolia, and Ipomoea nil.
A tree generated from ITS and matK sequences resolved the 14 species of Convolvulaceae into two major groups (Fig 1: Convolvuleae, Ipomoeeae), consistent with subfamilial groupings shown elsewhere in substantially more detailed taxonomic work [60]. The assay data were overlain on the tree to search for evidence that plant phylogeny predicted survival. Within   Fig 1. Phylogeny of assayed Convolvulaceae based on ITS and matK sequences. Node confidence was calculated using Neighbor Joining tree (Bootstrap replicates = 100). Species in red font followed by an insect kill icon failed to allow survival to adult stage; species in green font followed by a psyllid adult icon allowed egg-to-adult development. https://doi.org/10.1371/journal.pone.0201506.g001 Tribe Convolvuleae, psyllids of both haplotypes developed successfully on Co. arvensis and Co. tricolor, but failed to survive on Ca. silvatica and Co. equitans, despite their phylogenetic closeness to Co. arvensis (Fig 1). Within Tribe Ipomoeeae, psyllids of both haplotypes developed successfully on five species of Ipomoea, but failed to develop on four other Ipomoea or on T. corymbosa. Phylogenetic distance from the control host plant (potato) was calculated for each species of Convolvulaceae in the phylogenetic tree (distances calculated in Geneious 1 , S1 Table). A two-sample t-test demonstrated that mean phylogenetic distance from potato was statistically identical between plant species that allowed psyllid survival versus species on which psyllids failed to survive (P = 0.4563; Fig 2), confirming observations in the phylogenetic tree that phylogenetic nearness of a species to potato did not predict whether the psyllid would complete development on the plant (Fig 1).
We did not record actual rates of survival, so the assays cannot tell us whether percent survival on plant species that allowed egg-to-adult development was similar to survival on potato.
To determine if developmental rates on Convolvulaceae were similar to rates on potato, we compared mean number of days from oviposition to production of the first adult (N = 5 replicates per plant species and psyllid haplotype) between psyllids on potato and on those Convolvulaceae allowing survival to the adult stage. Development times varied between~20-35 days depending upon psyllid haplotype and plant species (Fig 3). Mean number of days between egg deposition and emergence of the first adult differed statistically between psyllid haplotypes (F = 15.5; df = 1, 57.0; P <0.001) and among plant species (F = 2.8; df = 7, 57.1; P = 0.013); the haplotype x plant species interaction was not significant (F = 1.4; df = 7, 57.1; P = 0.22) indicating that the effects of plant species on psyllid development time was similar between the two haplotypes. The Central haplotype developed more rapidly (mean = 24.7 ± 1.1 d) than the Northwestern haplotype (mean = 29.4 ± 1.3 d), when averaged across host plant. We extracted contrasts to examine two a priori defined comparisons of interest. A test of mean development time on potato vs. Convolvulaceae was significant (F = 18.01; df = 1, 57.02; P < 0.001), and showed that mean development time on potato was statistically shorter than development time on Convolvulaceae, averaged over haplotype (Fig 3). A second set of contrasts was extracted to examine whether there was evidence for plant effects within the Convolvulaceae, ignoring potato. Averaged over the two haplotypes, there was no evidence that development time of psyllids varied among species of Convolvulaceae (Fig 3:

Psyllid developmental success and a plant-fungus symbiosis
Visible evidence for the presence of Periglandula was most pronounced in two species, T. corymbosa and I. leptophylla (Fig 4A and 4B). The fungal colonies were found on the adaxial surfaces of younger leaves. Because visible evidence for presence of fungal colonies was rare, we used a molecular approach for detection of the fungus. Analysis of DNA-sequences led to detection of the dmaW gene in 11 of 14 plant species (Fig 4C and 4D;  [34,42]. Ergot alkaloids are categorized into three classes (clavines, simple amides of lysergic acid, and ergopeptines) based on their structural complexity and occurrence in the biochemical pathway [61,62]. Compounds from two classes (clavines, amides of lysergic acid) were detected in leaf tissues of plant species in which the dmaW gene (indicating presence of Periglandula) was also detected (Table 3). Compounds included eight clavines and two lysergic acid amides (Table 3). No ergopeptines were detected. Additionally, no ergot alkaloids were detected in species not shown to host Periglandula (Co. arvensis, I. nil, I. alba). However, the presence of Periglandula did not always lead to detection of alkaloids in plant tissues. Alkaloid content may vary with plant age or organ, with higher concentrations typically occurring in seeds and seedlings over vegetative parts [46,63]. This variation, combined with the possibility that ergot alkaloid concentrations can fall below detection limits, may lead to a failure in confirming presence of ergot alkaloids despite detection of Periglandula by molecular methods [55]. We observed often striking differences in alkaloid profiles between plant species that allowed psyllid development and species on which the psyllid failed to develop (Table 3). Plants in which clavines and amides of lysergic acid were readily detected were invariably fatal to nymphal psyllids (Table 3). Mortality was quite rapid on these species. Nymphs always died as first instars generally within 24-48 h following egg hatch (Kaur and Horton pers. observation). With two exceptions (Ca. silvatica, Co. equitans), plant species in which alkaloids were not detected allowed egg-to-adult development (Table 3). Psyllids failed to develop successfully on these two species despite a failure to detect alkaloids and despite detection of Periglandula in host tissues (Fig 4C, Table 3). Whether ergot alkaloids were actually present, but not detected, is not known. Lack of survival on Co. equitans may have been caused in part by the plant's extreme hairiness, as the pubescence was found to interfere with the ability of psyllids to feed and settle (from visual observations). Psyllids did successfully develop on four other species in which the dmaW gene was detected (Co. tricolor, I. cordatotriloba, I. hederacea, I. ternifolia). However, no ergot alkaloids were detected in leaf tissues from these four species, despite presence of the fungus (Table 3).

Discussion
This study adds to the list of Convolvulaceae that support egg-to-adult development of potato psyllid, and shows conclusively that the psyllid is able to develop on Convolvulaceae other than the two species (Co. arvensis and I. batatas) previously listed in literature accounts [9,11,12]. These additional taxa included an ornamental species of Convolvulus (Co. tricolor) likely of Mediterranean origin [64] and five species of New World Ipomoea. The Ipomoea comprised a mix of species that are grown as ornamentals (I. alba, I. nil, I. hederacea), and two species (I. cordatotriloba, I. ternifolia) that are present naturally in regions of Central America, Mexico, and the southwestern U.S. [28,65,66,67]. Previous accounts of association between potato psyllid and Convolvulaceae include rearing trials and field observations. Some care must be taken in interpretation of the field records, as field observations can lead to inflated ideas of true host range of psyllids due to the willingness of these insects to colonize and feed upon plant species that nonetheless fail to support nymphal development [7,11,68]. In this study, we followed the strict published guidelines of Burckhardt et al. [7] in defining psyllid "host plant" as a species that allows egg-to-adult development. A failure to appreciate this distinction has led to confusion about the host range of potato psyllid [11]. We obtained egg-laying and egg hatch (presence of nymphs) on all 14 species of Convolvulaceae that were assayed in this study, but development to the adult stage was limited to seven of these species.
Psyllids of the Central and Northwestern haplotypes were identical with respect to what plant species allowed successful development. The haplotypes did differ in development rates on species allowing development, with psyllids of the Central haplotype developing more rapidly than psyllids of the Northwestern haplotype. Other studies have shown that haplotypes of potato psyllid differ in biological traits, including settling and oviposition behavior [69], development rates [25], body size [24,25], and composition of endosymbiont communities [27]. The Central haplotype developed more rapidly on cultivated and weedy Solanaceae than psyllids of the Northwestern haplotype [25], which is consistent with our observations. It is likely that differences in development times were partly or largely due to differences between haplotypes in body size. Psyllids of the Northwestern haplotype are conspicuously larger than psyllids of the Central haplotype [24,25] and it seems likely that the size differences translated into differences in development times between the haplotypes.
We examined whether survival of psyllids on a given plant species could be predicted by location of the species in a phylogenetic tree. The Psylloidea have shown the ability to track phylogenetic diversification of plants within lineages, and host switching or dietary expansion in evolutionary or ecological time by psyllids appear to occur most often between phylogenetically related plants species [8,31,32]. One outcome of this sort of phylogenetic tracking is the expectation that dietary breadth for a given psyllid species would likely encompass phylogenetically related plant species rather than a mixture of less-related species. The phylogenetic tree developed from our sequencing work is consistent with trees constructed by earlier phylogenetic work for the Convolvulaceae [60]. Our sequences resolved the fourteen assayed species into two clades which fall respectively into two major tribes [60], the Convolvuleae and Ipomoeeae. The Ipomoeeae was further resolved into two clades [60]: the Argyreiinae, which includes one of our assayed species (Turbina corymbosa); and the Astripomoeinae clade, which contains the remaining Ipomoeeae (all species of Ipomoea) that were assayed. Our data failed to show that developmental success of psyllids was affected by location of plants in our phylogenetic tree. Plant species that allowed development were represented in both Tribes of Convolvulaceae that were assayed here, as were species that failed to allow development (Fig 1).
Observations in the literature indicate that species of Convolvulaceae may often harbor a class of alkaloids (ergot alkaloids) known in grasses to confer resistance to insect herbivory [43,45,61]. These compounds are produced in grasses by fungal species in the family Clavicipitaceae (genus Epichloё) which have formed a mutualistic relationship with grasses. A similar mutualistic association between Convolvulaceae and clavicipitaceous fungi in a different genus (Periglandula) has been shown to explain the presence of ergot alkaloids in Convolvulaceae [35,36,37,38,39]. The visual presence of fungal colonies on at least some of our targeted species (Fig 4A and 4B), combined with extensive literature confirming the presence of ergot alkaloids in Convolvulaceae, prompted us to examine whether psyllid development or lack of development was correlated with the presence or absence of ergot alkaloids.
We detected Periglandula in a surprisingly large proportion of assayed plants (11 of 14 species), including in two genera (Convolvulus, Calystegia) not previously known to host this fungal symbiont. Previous surveys have suggested that the occurrence of ergot alkaloids (and thus this mutualistic association) was limited to the tribe Ipomoeeae and two clades (Argyreiinae and Astripomoeinae) within this tribe (data based on analyses of 46 species) [34,42]. It has now been estimated that approximately 50% of Ipomoeeae species, or upwards of 450 species worldwide, could contain ergot alkaloids. These observations understandably have led researchers to focus on the tribe Ipomoeeae in efforts to document presence of ergot alkaloids [34,42,37] and it is possible that this focus has led workers to substantially underestimate the taxonomic diversity of Convolvulaceae actually harboring ergot alkaloids. The few reports in the literature suggesting that ergot alkaloids in the Convolvulaceae occur outside of the Ipomoeeae, including in Calystegia and Convolvulus, have been categorized as "unverified" [34]. Our results are the first to demonstrate that the presence of Periglandula in Convolvulaceae does indeed extend outside of Ipomoeeae.
Ergot alkaloids representing two classes (clavines and amides of lysergic acid) were detected in five of our assayed plant species (Table 3). Previous literature accounts summarized in Eich (2008) report these same two classes of alkaloids in four of these five species (failing to list only I. pandurata). These same accounts identified many of the same specific compounds that were identified in this study [34,37]. All species in this study which showed presence of ergot alkaloids were shown (with PCR) to also host Periglandula. No ergot alkaloids were detected in the three species in which we failed to detect Periglandula (I. nil, I. alba, Co. arvensis). This result is consistent with other studies of these three species [34].
We invariably observed 100% mortality of psyllid nymphs on species in which ergot alkaloids were detected (Table 3). Mortality occurred very rapidly following egg hatch, generally within 24-48 h of hatch. Assuming that nymphal mortality was due to the presence of these alkaloids, the next logical question is what mode of action explains our results? Absence of development could have been caused by direct toxicity of the alkaloids or because the compounds deter feeding. At this time, we cannot separate these effects. Insecticidal activity of this class of alkaloids could arise from their capacity to act as agonists or antagonists to neurotransmitter receptors and subsequent malfunctioning of the central nervous system [61]. However, it is also possible that the compounds deterred feeding enough that newly hatched nymphs rapidly desiccated and died. An evaluation of these competing effects will require additional assays, likely including assays that allow measurement of feeding rates (e.g., production of honeydew). Studies in which synthetic analogues of targeted compounds are assayed would also be useful, as use of synthesized compounds would allow insect responses to be examined relative to specific concentrations of compounds or to mixtures of compounds [45,70].
We failed to detect ergot alkaloids in six species that nonetheless were shown by PCR to harbor Periglandula (Table 3). It is unclear if the alkaloids were actually present but were not at detectable levels, if ergot alkaloids were present but were different compounds than targeted by our biochemistry work, or if indeed alkaloids were not present at all. Efforts to detect ergot alkaloids in Convolvulaceae can lead to inconsistent results, even in assays of plant species known from previous studies to harbor the chemicals [34,55]. These inconsistencies may be the consequence of any of a number of factors, including sensitivity of the analytical approach chosen to look for alkaloids, age of the plant seed or conditions under which the seed was stored, age of the plant, which plant structures are examined, and incorrect taxonomic work leading to mistakes in species identification [34,71]. Alkaloid levels within a single plant may vary with plant structure. Levels in vegetative tissues, as were targeted here, may be lower than levels in other plant parts, such as seed or newly expanded cotyledons [46,63]. It may be that analysis and extraction of plant structures other than those that were targeted here (the fully expanded leaf) would have led to detection of ergot alkaloids in those species found to harbor Periglandula but in which we failed to detect the chemicals. Potato psyllid successfully completed development on five species in which Periglandula was present but in which ergot alkaloids were not detected. If ergot alkaloids do have psyllicidal effects, as suggested by our results in Fig 4C and 4D and Table 3, then successful development by psyllids on those five Periglandula-positive species from which we failed to detect alkaloids may indicate that alkaloids were indeed not present, or that they were at levels low enough to allow psyllid development and to escape biochemical detection.
Symbiotic association between plants and clavicipitaceous fungi is best known for monocotyledonous plants (Poaceae, Cyperaceae and Junaceae), where (as with Convolvulaceae) the symbioses lead to production of ergot alkaloids [41,72,73]. These associations may lead to any of several benefits for the plant, notably protection against herbivores, but including also nondefense type functions such as enhanced growth rates of the plant or increased ability to withstand drought or other abiotic stresses [74,75,76,77]. Observations that benefits to plants may include multiple types of effects, combined with observations showing that these effects are not always predictable across studies, species, or environments, have led to a large body of literature debating the actual evolutionary processes leading to these associations [77,78,79,80]. Our results provide correlative evidence that presence of ergot alkaloids in Convolvulaceae prevents development of psyllid nymphs, suggesting that the Periglandula-Convolvulaceae symbiosis does lead to protection of plants against insect herbivores. Our results also showed, however, that presence of the fungus does not necessarily indicate that psyllids would not survive on the plant host, as species in which Periglandula was present but from which alkaloids were not detected did allow egg-to-adult development by psyllids. Future studies will include screening of a larger diversity of Convolvulaceae than assayed here, comprising both Periglandula-positive and Periglandula-negative species, and we believe that this larger study will shed additional light on the role of this fungal symbiosis in affecting fitness of phloem-feeding insects.
Supporting information S1 Table. Genetic distance from potato as calculated by Geneious1. (DOCX)