Chemically defended plant tissues present formidable barriers to herbivores. Although mechanisms to resist plant defenses have been identified in ancient herbivorous lineages, adaptations to overcome plant defenses during transitions to herbivory remain relatively unexplored. The fly genus Scaptomyza is nested within the genus Drosophila and includes species that feed on the living tissue of mustard plants (Brassicaceae), yet this lineage is derived from microbe-feeding ancestors. We found that mustard-feeding Scaptomyza species and microbe-feeding Drosophila melanogaster detoxify mustard oils, the primary chemical defenses in the Brassicaceae, using the widely conserved mercapturic acid pathway. This detoxification strategy differs from other specialist herbivores of mustard plants, which possess derived mechanisms to obviate mustard oil formation. To investigate whether mustard feeding is coupled with evolution in the mercapturic acid pathway, we profiled functional and molecular evolutionary changes in the enzyme glutathione S-transferase D1 (GSTD1), which catalyzes the first step of the mercapturic acid pathway and is induced by mustard defense products in Scaptomyza. GSTD1 acquired elevated activity against mustard oils in one mustard-feeding Scaptomyza species in which GstD1 was duplicated. Structural analysis and mutagenesis revealed that substitutions at conserved residues within and near the substrate-binding cleft account for most of this increase in activity against mustard oils. Functional evolution of GSTD1 was coupled with signatures of episodic positive selection in GstD1 after the evolution of herbivory. Overall, we found that preexisting functions of generalized detoxification systems, and their refinement by natural selection, could play a central role in the evolution of herbivory.

SIGNIFICANCE STATEMENTThe origin of land plants >400 million years ago (mya) spurred the diversification of plant-feeding (herbivorous) insects and triggered an ongoing chemical co-evolutionary arms race. Because ancestors of most herbivorous insects first colonized plants >200 mya, the sands of time have buried evidence of how their genomes changed with their diet. We leveraged the serendipitous intersection of two genetic model systems: a close relative of yeast-feeding fruit fly ( Drosophila melanogaste r), the "wasabi fly" ( Scaptomyza flava ), that evolved to consume mustard plants including Arabidopsis thaliana . The yeast-to-mustard dietary transition remodeled the fly's gene repertoire for sensing and detoxifying chemicals. Although many genes were lost, some underwent duplications that encode the most efficient detoxifying enzymes against mustard oils known from animals. Gloss et al. 2019 1ABSTRACT One-quarter of extant Eukaryotic species are herbivorous insects, yet the genomic basis of this extraordinary adaptive radiation is unclear. Recently-derived herbivorous species hold promise for understanding how colonization of living plant tissues shaped the evolution of herbivore genomes. Here, we characterized exceptional patterns of evolution coupled with a recent (<15 mya) transition to herbivory of mustard plants (Brassicaceae, including Arabidopsis thaliana ) in the fly genus Scaptomyza, nested within the paraphyletic genus Drosophila . We discovered a radiation of mustard-specialized Scaptomyza species, comparable in diversity to the Drosophila melanogaster species subgroup. Stable isotope, behavioral, and viability assays revealed these flies are obligate herbivores. Genome sequencing of one species, S. flava, revealed that the evolution of herbivory drove a contraction in gene families involved in chemosensation and xenobiotic metabolism. Against this backdrop of losses, highly targeted gains ("blooms") were found in Phase I and Phase II detoxification gene sub-families, including glutathione S-transferase ( Gst ) and cytochrome P450 ( Cyp450 ) genes. S. flava has more validated paralogs of a single Cyp450 (N=6 for Cyp6g1 ) and Gst (N=5 for GstE5-8 ) than any other drosophilid. Functional studies of the Gst repertoire in S. flava showed that transcription of S. flava GstE5-8 paralogs was differentially regulated by dietary mustard oils, and of 22 heterologously expressed cytosolic S. flava GST enzymes, GSTE5-8 enzymes were exceptionally well-adapted to mustard oil detoxification in vitro . One, GSTE5-8a, was an order of magnitude more efficient at metabolizing mustard oils than GSTs from any other metazoan. The serendipitous intersection of two genetic model organisms, Drosophila and Arabidopsis , helped illuminate how an insect genome was remodeled during the evolutionary transformation to herbivory, identifying mechanisms that facilitated the evolution of the most diverse guild of animal life.

Identifying the factors governing the maintenance of genetic variation is a central challenge in evolutionary biology. New genomic data, methods and conceptual advances provide increasing evidence that balancing selection, mediated by antagonistic species interactions, maintains functionally-important genetic variation within species and natural populations. Because diverse interactions between plants and herbivorous insects dominate terrestrial communities, they provide excellent systems to address this hypothesis. Population genomic studies of Arabidopsis thaliana and its relatives suggest spatial variation in herbivory maintains adaptive genetic variation controlling defense phenotypes, both within and among populations. Conversely, inter-species variation in plant defenses promotes adaptive genetic variation in herbivores. Emerging genomic model herbivores of Arabidopsis could illuminate how genetic variation in herbivores and plants interact simultaneously.

The serine/threonine protein kinase 3-phosphoinositide-dependent protein kinase 1 (PDK1) is a highly conserved eukaryotic kinase that is a central regulator of many AGC kinase subfamily members. Through its regulation of AGC kinases, PDK1 controls many basic cellular processes, from translation to cell survival. While many of these PDK1-regulated processes are conserved across kingdoms, it is not well understood how PDK1 may have evolved within kingdoms. In order to better understand PDK1 evolution within plants, we have isolated and characterized the PDK1 gene from the moss Physcomitrella patens (PpPDK1), a nonvascular representative of early land plants. PpPDK1 is similar to other plant PDK1s in that it can functionally complement a yeast PDK1 knockout line. However, unlike PDK1 from other plants, the P. patens PDK1 protein does not bind phospholipids due to a lack of the lipid-binding pleckstrin homology domain, which is used for lipid-mediated regulation of PDK1 activity. Sequence analysis of several PDK1 proteins suggests that lipid regulation of PDK1 may not commonly occur in algae and nonvascular land plants. PpPDK1 can phosphorylate AGC kinase substrates from tomato (Solanum lycopersicum) and P. patens at the predicted PDK1 phosphorylation site, indicating that the PpPDK1 substrate phosphorylation site is conserved with higher plants. We have also identified residues within the PpPDK1 kinase domain that affect kinase activity and show that a mutant with highly reduced kinase activity can still confer cell viability in both yeast and P. patens. These studies lay the foundation for further analysis of the evolution of PDK1 within plants.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts with critical biological functions. Despite the wealth of RNA-seq data available, lincRNA identification lags in the plant lineage. In addition, there is a need for a harmonized identification and annotation effort to enable cross-species functional and genomic comparisons. In this study we processed >24 Tbp of RNA-seq data from >16,000 experiments to identify ~130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used Nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify functional motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our dataset with sequence or transcriptional evolutionary conservation, including lincRNAs transcribed adjacent to orthologous genes that display little sequence similarity and likely function as transcriptional regulators. Finally, we used guilt-by-association techniques to further classify these lincRNAs according to putative function. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, short ORFs, and whose expression is modulated by abiotic stress are a few of the annotations that will prioritize and guide future functional analyses.

The tomato AGC kinase Adi3 is phosphorylated by Pdk1 for activation of its cell death suppression activity. The Pdk1 phosphorylation site for activation of Adi3 is at Ser539. However, there is at least one additional Pdk1 phosphorylation site on Adi3 that has an unknown function. Here we identify an Arabidopsis thaliana sequence homologue of Adi3 termed AGC1-3. Two Pdk1 phosphorylation sites were identified on AGC1-3, activation site Ser596 and Ser269, and by homology Ser212 on Adi3 was identified as a second Pdk1 phosphorylation site. While Ser212 is not required for Adi3 autophosphorylation, Ser212 was shown to be required for full phosphorylation of the Adi3 substrate Gal83.

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed &gt;24 Tbp of RNA-seq data from &gt;16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used Nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our dataset with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our dataset. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Protein kinases belonging to the AGC group modulate many diverse cellular processes in all eukaryotes. One important way to regulate AGC kinases is through phosphorylation by the upstream kinase PDK1. PDK1 localization and activity usually depend on interactions with phospholipids, which are mediated by a conserved lipid-binding pleckstrin homology (PH) domain. We recently analyzed putative PDK1 sequences from 17 photosynthetic organisms, finding that PDK1s from vascular and nonvascular species seem to be distinguished by the presence or absence of a PH domain, respectively. The only other reported PDK1 lacking a PH domain is from yeast (Saccharomyces cerevisiae). These observations raise questions about how plant PDK1s and their lipid-binding capabilities have evolved in relation to other eukaryotes, and what this means for PDK1 function. Here we use 100 PDK1 sequences from diverse organisms to discuss possible evolutionary aspects of plant PDK1 structure and lipid binding.