Specialist herbivores have often evolved highly sophisticated mechanisms to counteract defenses mediated by major plant secondary-metabolites. Plant species of the herbivore host range often display high chemical diversity and it is not well understood how specialist herbivores respond to this chemical diversity. Pieris larvae overcome toxic products from glucosinolate hydrolysis, the major chemical defense of their Brassicaceae hosts, by expressing nitrile-specifier proteins (NSP) in their gut. Furthermore, Pieris butterflies possess so-called major allergen (MA) proteins, which are multi-domain variants of a single domain major allergen (SDMA) protein expressed in the guts of Lepidopteran larvae. Here we show that Pieris larvae fine-tune NSP and MA gene expression depending on the glucosinolate profiles of their Brassicaceae hosts. Although the role of MA is not yet fully understood, the expression levels of NSP and MA in larvae that fed on plants whose glucosinolate composition varied was dramatically changed, whereas levels of SDMA expression remained unchanged. In addition, we found a similar regulation pattern among these genes in larvae feeding on Arabidopsis mutants with different glucosinolate profiles. Our results demonstrate that Pieris larvae appear to use different host plant adaptive genes to overcome a wide range of glucosinolate profiles in their host plants.

Adaptive traits that enable organisms to conquer novel niches and experience subsequent diversification are ecologically and evolutionarily important. The larvae of Pieris butterflies express nitrile-specifier proteins (NSPs), a key innovation for overcoming the glucosinolate (GLS)-myrosinase-based defence system of their Brassicales host plants. Nitrile-specifier proteins are a member of the NSP-like gene family, which includes the major allergen (MA) protein, a paralog of NSP with a GLS-disarming function, and a single domain major allergen (SDMA) protein, whose function is unknown. The arms-race between GLS-based defences and the NSP-like gene family is suggested to mediate diversification in both Pierid butterflies and Brassicales plants. Here, we tested whether the expected strong selection on NSP-like gene family correlates with shifts in host plant spectra among Pierid butterflies. We combined feeding experiments using 25 Brassicaceae plants and five Pieris species with larval transcriptome data to investigate the patterns of selection acting on NSP-like gene family members. Although we observed significantly elevated nonsynonymous to synonymous substitution rate ratios in NSPs on branches associated with changes in patterns of host plant usage, no such pattern was observed in MAs or SDMAs. Furthermore, we found evidence for positive selection of NSP at a phylogenetic branch which reflects different host plant spectra. Our data indicate that the NSP-related gene members have evolved differently: NSPs have accumulated more amino acid changes in response to shifting preferences for host plants, whereas MAs and SDMAs appear to be more conserved. Further detailed functional assays of these genes would provide important insights to understand their role in the chemical arms-race between Pieris butterflies and their Brassicales host plants. K E Y W O R D Sarms-race, host plant adaptation, insects, selection | 4959 OKAMURA et Al.

The tremendous diversity of plants and herbivores has arisen from a coevolutionary relationship characterized by plant defense and herbivore counter adaptation. Pierid butterfly species feed on Brassicales plants that produce glucosinolates as a chemical deterrent against herbivory. In turn, the larvae of pierids have nitrile specifier proteins (NSPs) that are expressed in their gut and disarm glucosinolates. Pierid butterflies are known to have diversified in response to glucosinolate diversification in Brassicales. Therefore, each pierid species is expected to have a spectrum of host plants characterized by specific glucosinolate profiles. In this study, we tested whether the larval performance of different Pieris species, a genus in Pieridae (Lepidoptera: Pieridae), was associated with plant defense traits of putative host plants. We conducted feeding assays using larvae of three Pieris species and 10 species of the Brassicaceae family possessing different leaf physical traits and glucosinolate profile measurements. The larvae of Pieris rapae responded differently in the feeding assays compared with the other two Pieris species. This difference was associated with differences in glucosinolate profiles but not with variations in physical traits of the host plants. This result suggests that individual Pieris species are adapted to a subset of glucosinolate profiles within the Brassicaceae. Our results support the idea that the host ranges of Pieris species depend on larval responses to glucosinolate diversification in the host species, supporting the hypothesis of coevolution between butterflies and host plants mediated by the chemical arms race.

20Adaptive traits that enable organisms to conquer novel niches and experience subsequent 21 diversification are ecologically and evolutionarily important. The larvae of Pieris butterflies 22 express nitrile-specifier proteins (NSPs), a key innovation for overcoming the glucosinolate (GLS)-23 myrosinase-based defense system of their Brassicales host-plants. NSPs are a member of the NSP-24 like gene family, which includes the major allergen (MA) protein, a paralog of NSP with a GLS-25 disarming function, and a single domain major allergen (SDMA) protein, whose function is 26 unknown. The arms-race between a highly variable host-plant defense system and members of the 27 NSP-like gene family is suggested to mediate diversification in both Pierid butterflies and 28Brassicales plants. Here, we combined feeding experiments using 25 Brassicaceae plants and five 29Pieris species with larval transcriptome data to investigate the evolutionary forces acting on NSP-30 like gene family members associated with patterns of host-plant usage. Although we observed 31 significantly elevated nonsynonymous to synonymous substitution ratios in NSPs, no such pattern 32 was observed in MAs or SDMAs. Furthermore, we found a signature of positive selection of NSP 33 at a phylogenetic branch which reflects different host-plant preferences. Our data indicate that 34NSPs have evolved in response to shifting preferences for host plants among five Pieris butterflies, 35whereas MAs and SDMAs appear to have more conserved functions. Our results show that the 36 evolution and functional differentiation of key genes used in host-plant adaptation play a crucial 37 role in the chemical arms-race between Pieris butterflies and their Brassicales host-plants. 38 39 40