Physiology Drives Reworking of Amino Acid δ2H and δ13C in Butterfly Tissues

Studies of animal movement and migration over large geospatial scales have long relied on natural continental-scale hydrogen isotope (δ2H) gradients in precipitation, yet the physiological processes that govern incorporation of δ2H from precipitation into plant and then herbivore tissues remain poorly understood, especially at the molecular level. Establishing a biochemical framework for the propagation of δ2H through food webs would enable us to resolve more complicated regional-scale animal movements and potentially unlock new applications for δ2H data in animal ecology and eco-physiology. Amino acid δ2H analysis offers a promising new avenue by which to establish this framework. We report bulk tissue δ2H, δ13C, and δ15N data as well as amino acid δ2H and δ13C data from three Pipevine swallowtail (Battus philenor) tissues—caterpillars, butterfly bodies, and wings—as well as their obligate plant source: pipevine leaves (Aristolochia macrophylla). Insects are often dominant herbivores in terrestrial food webs and a major food source for many higher-level consumers, so it is particularly important to understand the mechanisms that influence insect tissue δ2H values. Our data reveal extensive δ2H variation within and among individuals of a relatively simple plant-herbivore system that cannot be explained by temporal or geospatial gradients of precipitation δ2H or dietary differences. Variations in essential amino acid δ2H and δ13C indicate that B. philenor acquire these compounds from an additional source that is isotopically distinct from pipevine leaves, potentially gut microbes. We also found multiple isotopic carryover effects associated with metamorphosis. This study emphasizes the strong influence of physiology on consumer-diet δ2H discrimination in a local population of pipevines and swallowtails and provides a template that can be broadly applied to Lepidoptera—the second most diverse insect order—and other holometabolous insects. Understanding these physiological mechanisms is critical to interpreting the large degree of δ2H variation in consumer tissues often observed at a single collection site, which has implications for using δ2H isoscapes to study animal movement. Further investigation into amino acid δ2H holds promise to elucidate how subsets of amino acids may be best utilized to address specific ecological and physiological questions for which bulk tissue δ2H is insufficient.

Studies of animal movement and migration over large geospatial scales have long relied on natural continental-scale hydrogen isotope (δ 2 H) gradients in precipitation, yet the physiological processes that govern incorporation of δ 2 H from precipitation into plant and then herbivore tissues remain poorly understood, especially at the molecular level. Establishing a biochemical framework for the propagation of δ 2 H through food webs would enable us to resolve more complicated regional-scale animal movements and potentially unlock new applications for δ 2 H data in animal ecology and eco-physiology. Amino acid δ 2 H analysis offers a promising new avenue by which to establish this framework. We report bulk tissue δ 2 H, δ 13 C, and δ 15 N data as well as amino acid δ 2 H and δ 13 C data from three Pipevine swallowtail (Battus philenor) tissues-caterpillars, butterfly bodies, and wings-as well as their obligate plant source: pipevine leaves (Aristolochia macrophylla). Insects are often dominant herbivores in terrestrial food webs and a major food source for many higher-level consumers, so it is particularly important to understand the mechanisms that influence insect tissue δ 2 H values. Our data reveal extensive δ 2 H variation within and among individuals of a relatively simple plant-herbivore system that cannot be explained by temporal or geospatial gradients of precipitation δ 2 H or dietary differences. Variations in essential amino acid δ 2 H and δ 13 C indicate that B. philenor acquire these compounds from an additional source that is isotopically distinct from pipevine leaves, potentially gut microbes. We also found multiple isotopic carryover effects associated with metamorphosis. This study emphasizes the strong influence of physiology on consumer-diet δ 2 H discrimination in a local population of pipevines and swallowtails and provides a template that can be broadly applied to Lepidoptera-the second most diverse insect order-and other holometabolous insects. Understanding these physiological mechanisms is critical to interpreting the large degree of δ 2 H variation in consumer tissues often observed at a single collection site, which has implications for using δ 2 H isoscapes to study animal movement. Further investigation into amino acid δ 2 H holds promise to elucidate how subsets of amino acids may be best utilized to address specific ecological and physiological questions for which bulk tissue δ 2 H is insufficient.

INTRODUCTION
For decades, numerous studies have traced terrestrial animal migrations and movements by relying on natural continentalscale gradients in hydrogen isotopes of precipitation and groundwater that intrinsically label animal tissues (Wassenaar and Hobson, 1998;Kelly et al., 2002;Altizer et al., 2015). Because of the logistical and technological challenges associated with tracking, capturing, and sampling highly mobile animals that undergo large migrations, molecular techniques like hydrogen isotope (δ 2 H) analysis are incredibly valuable to studies of movement ecology (Kelly et al., 2002;Hobson and Wassenaar, 2008). Yet, it remains difficult to interpret these isotope data, especially when local site-level isotopic variation limits inferences of geographic affinity. First, bulk tissue (e.g., feather) isotope values are influenced by both the isotope signature at the base of the food web and biochemical modifications that occur within consumers during resource assimilation and tissue synthesis (McMahon and Newsome, 2019). Second, the isotopic composition of a consumer's tissue represents a weighted average of all macromolecules within the sample. Lipids and carbohydrates are known to have vastly different δ 13 C and δ 2 H values than proteins (Estep and Hoering, 1981;Terwilliger et al., 2002;Post et al., 2007); lipids also vary in concentration among different tissues and often fluctuate temporally with body condition and reproductive status, especially in invertebrates (Cherel et al., 2005;Wolf et al., 2009). Moreover, migrations are often associated with significant physiological shifts including changes in metabolic rate, energy expenditure, and fasting, depletion of fat reserves, and variation in hormones-any of which could potentially affect bulk tissue isotope values (Gannes et al., 1998;Gwinner, 2012). Consequently, it can be difficult to determine how baseline isotope values, diet, tissue type, physiology, or metabolism each contribute to consumer bulk tissue isotope data.
Despite these obstacles, bulk tissue δ 2 H data can be successfully implemented in movement and migration studies at large spatial scales if δ 2 H isoscapes are well-constrained and δ 2 H gradients are large enough (Hobson et al., 2003;Bowen et al., 2005). Many of these studies have relied on the feathers of migratory birds because they are metabolically inert once grown, and feather δ 2 H values have repeatedly been shown to correlate with those of the local precipitation where they are synthesized (Chamberlain et al., 1997;Rubenstein et al., 2002;Hobson et al., 2004). However, δ 2 H in animal tissues is influenced by a combination of drinking/food water and diet, and the latter source of hydrogen accounts for the majority used to synthesize proteinaceous tissues like feathers (Hobson, 1999;Wolf et al., 2011Wolf et al., , 2012. For this reason, and because hydrogen in a molecule has many opportunities to exchange with other pools of hydrogen during metabolic reactions, hydrogen cycling is more complex than that of more conventional biomarkers such as carbon or nitrogen isotopes. In fact, nearly all the research into the biochemical processes that govern incorporation of δ 2 H into consumer tissues stems from controlled feeding experiments in the lab (Estep and Dabrowski, 1980;Peters et al., 2012;Fogel et al., 2016;Newsome et al., 2017;Rodriguez Curras et al., 2018), with relatively little being done in natural ecosystems across trophic levels. Establishing a biochemical framework for the propagation of δ 2 H through food webs would expand our ability to resolve smaller regional-scale animal movements and potentially unlock new applications for δ 2 H data in dietary and metabolic studies.
Amino acid hydrogen isotope analysis offers a promising new avenue by which to establish this framework, without the confounding factors associated with bulk tissue analysis. Amino acids are the building blocks of most tissues (e.g., feather, muscle, liver, blood) analyzed by ecologists, making them ideal compounds to study when examining how isotopes propagate through food webs. Although amino acid δ 2 H datasets are currently limited in the literature, a controlled feeding study on E. coli found that the δ 2 H of individual amino acids within a single sample can vary by ∼200 (Fogel et al., 2016)-an extreme amount of intermolecular variation that would be masked by a composite bulk tissue δ 2 H value. Similar patterns among amino acid δ 2 H values were recently reported for muscle tissue of house mice (Mus musculus) reared in controlled laboratory conditions (Newsome et al., 2020), suggesting these δ 2 H values are primarily controlled by biochemical processes governing amino acid synthesis that are shared among prokaryotes and eukaryotes. While this extensive variation can present challenges for ecological studies, it also provides an opportunity to examine the mechanisms driving the transfer of δ 2 H signals up food chains (e.g., metabolic routing) at the molecular level and on a scale that is an order of magnitude larger than the other light isotopes.
Unlike δ 2 H, amino acid δ 13 C analysis is well established in ecology. Amino acids are classified into two groups with respect to δ 13 C-essential and non-essential. Essential amino acids (AA Ess ) are particularly useful tracers in diet and food web studies because their carbon skeletons cannot be synthesized by animals, and so are directly routed into a consumer's tissues from the diet with minimal isotopic alteration (O'Brien et al., 2002;McMahon et al., 2010). Thus, the δ 13 C values of AA Ess in consumer tissues closely reflect those of the primary producers, bacteria, and fungi at the base of their food web, which have unique amino acid δ 13 C patterns or "fingerprints" that likely reflect variation in the biochemical pathways these groups use to synthesize these compounds (Scott et al., 2006;Larsen et al., 2009). In contrast, non-essential amino acids (AA NEss ) can either be directly routed from dietary protein or synthesized de novo by animals from dietary macromolecules (e.g., carbohydrates or lipids) or catabolized from endogenous stores (e.g., muscle or adipose tissue). The AA Ess and AA NEss classification scheme used in δ 13 C studies appears to hold up for δ 2 H, probably because most of the hydrogen (85-90%) in amino acids is bonded to carbon and theoretically does not exchange with body water (Wassenaar and Hobson, 2000;Fogel et al., 2016;Newsome et al., 2017). A controlled feeding study on E. coli found that a higher proportion of environmental water was used to synthesize AA NEss , especially glycolytic amino acids, relative to AA Ess in both protein-rich (tryptone) and protein-free (glucose) media treatments (Fogel et al., 2016). This finding illustrates how amino acid δ 2 H analysis can be employed to identify the sources of hydrogen from both diet (AA Ess ) and drinking or environmental water (AA NEss ) used by consumers to synthesize tissues.
To investigate the biochemical mechanisms by which hydrogen is incorporated into consumer tissues in a natural setting, we rely on a combination of bulk tissue and amino acid isotope data from a simple plant-herbivore trophic transfer in a single watershed. We exploit an obligate caterpillar-host plant association between pipevine plants in the genus Aristolochia and the Pipevine swallowtail (Battus philenor) (Figure 1; Rausher, 1980). Both the caterpillars and adult butterflies rely on flowering C 3 plants in the same geographical region and should retain the isotopic fingerprint of that primary producer group in their own AA Ess . Our approach eliminates the confounding factor of large latitudinal or elevational δ 2 H gradients in meteoric water. Furthermore, because adults of this species feed on nectar primarily composed of sugars (Scott, 1986;Lotts and Naberhaus, 2020), they can be protein-limited. This provides the unique opportunity to investigate isotopic changes associated with metamorphosis and a limited nitrogen pool, with little influence from the adult diet. In theory, therefore, any isotopic offsets in AA Ess or AA NEss among plants, caterpillars, and adults are the result of physiological effects associated with metamorphosis. Butterfly body tissue continues to turn over throughout the adult's short lifetime. However, because little additional nitrogen is consumed during this time, any isotopic offsets in amino acids between butterfly body tissue and wing tissue should either reflect tissue-specific discrimination during metamorphosis or metabolic processing during the adult's lifetime. Butterfly wing chitin is metabolically inert once formed, such that its isotopic signature reflects only larval resources (Wassenaar and Hobson, 1998). Wings also contain metabolically active hemolymph, but it represents a small component of the wing and is unlikely to significantly affect our isotope data or interpretations. See Supporting Information-Text for details.
A complicating influence on insect tissue isotope values may be gut microbiomes that have been identified in several Lepidopteran larvae (Robinson et al., 2010;Hammer et al., 2014) and adults (Paniagua Voirol et al., 2018;Hammer et al., 2020).
Unlike insects, gut microbes are capable of synthesizing AA Ess de novo using all dietary macromolecules, which could impart a unique isotopic signature if microbially-synthesized amino acids are used by the host to build tissues (Christensen and Fogel, 2011;Newsome et al., 2020). Such isotope effects are poorly understood even though insects are often the dominant herbivore group in terrestrial food webs (Price et al., 1980;Chikaraishi et al., 2011). Alternatively, some AA Ess could originate from nectar (Levin et al., 2017). If nectar-sourced amino acids supplied a significant proportion of amino acids for insect tissue synthesis, we would expect that the isotopic compositions of plant and insect tissues to be similar.
Our results reveal a high degree of intra-and inter-individual δ 2 H variation in caterpillars and butterflies despite the absence of geographical variation in the δ 2 H of water and food sources. This variation emphasizes the complexity inherent in natural systems, even ones as simple as an obligate caterpillar-host plant interaction. The influence of metabolism on consumer-diet δ 2 H (and δ 13 C) discrimination in our study provides a template for how δ 2 H might fractionate in response to physiological shifts in metamorphosing insects.

Study Site and Species
Battus philenor L. is found in warm climates throughout North America, mainly across the southern and eastern portions of the United States and in Mexico (Scott, 1986). This species is non-migratory (Hall, 2020) and females feed on a variety of flowering plants (e.g., thistles, phlox, ironweed) but only oviposit on pipevine species in the genus Aristolochia (Scott, 1986;Lotts and Naberhaus, 2020). In our study system, Battus larvae fed on the foliage of Aristolochia macrophylla Lam., a large liana with ropy stems that ascends into the forest canopy. Large vines can bear hundreds or even thousands of large heartshaped leaves. Little is known about the nectar sources exploited by Battus in the Santeetlah Creek watershed, the amino acid concentrations of those nectars, or the many other protein sources exploited opportunistically. These unknowns generate potential uncertainty when interpreting results and are addressed in the discussion.

Sample Collection
Specimens were collected throughout a narrow elevational band (872-1,191 m above sea level) from the watersheds of Big Santeetlah Creek (5350 ha) and Little Santeetlah Creek, Graham County, North Carolina, United States from 13 to 15 June 2013. These watersheds are heavily forested, dominated by C 3 plants, and situated on the eastern slope of the Unicoi Mountains, a subdivision of the southern Appalachians (Graves and Romanek, 2009). Caterpillars (n = 8) were collected directly from pipevine leaves, which were also collected, along with butterfly samples (n = 8); all insects were stored in cryovials in liquid nitrogen in the field. Specimens were subsequently freeze-dried and stored at room temperature in the lab until isotope analysis. Only a subset of the leaves collected were analyzed (n = 6), each from a different plant along a ∼4 km long transect parallel to Santeetlah Creek. For amino acid analysis, we sampled internal tissues, avoiding the gut. Water samples were also collected from the Santeetlah Creek Watershed and analyzed for bulk δ 2 H (average = -42 ; see Supporting Information-Text for details). Butterflies were sexed in the lab on the basis of dorsal hindwing iridescence, which is brighter in males (Supplementary Figure 1; Hall, 2020). All were males, except G-14. Variation in isotope data were unrelated to butterfly sex and elevation of sample collection (see Supporting Information-Text for details).

Sample Preparation
Approximately 10-15 mg of each dry sample was hydrolyzed for 20 h at 110 • C in 1 mL of 6 N HCl, with a N 2 headspace to remove oxygen. Hydrolysates were evaporated to dryness under a gentle N 2 stream at room temperature then prepared as TFAA derivatives (Silfer et al., 1991). The total free amino acids in the hydrolysates were esterified with acetyl chloride/2-propanol (1:4, v/v) at 110 • C for 1 h then acylated with trifluoroacetic anhydride/dichloromethane (1:1, v/v; DCM) at 110 • C for 10 min. The resultant derivatives were completely dried then brought up in DCM and stored at -25 • C until analysis. A standard mixture of 13 amino acids purchased from Sigma Aldrich (Saint Louis, MO) with known isotope values was derivatized alongside every batch of samples. See Supporting Information for details on bulk tissue analysis and findings.

Amino Acid Analysis
All samples were analyzed at the EDGE Institute at the University of California Riverside. The δ 13 C and δ 2 H values of individual amino acids were determined by gas chromatographycombustion-isotope ratio mass spectrometry (GC-C-IRMS) using a Trace 1310 gas chromatograph (GC) coupled with a GC-Isolink ll routed through a Conflo IV interface into a Delta V Plus isotope ratio mass spectrometer (Thermo Fisher Scientific, United States). Amino acids were either routed through a combustion reactor at 1,000 • C for δ 13 C or an HTC reactor at 1,420 • C for δ 2 H. Using a Triplus RSH autosampler, samples were injected splitless at 250 • C onto a BPX-5 capillary column (60 m × 0.32 mm inner diameter, 1.0 um film thickness; SGE Analytical Science, United States). The GC oven temperature program was as follows: initial temperature 50 • C for 2 min, followed by a ramp of 15 • C/min to 125 • C, a ramp of 3 • C/min to 160 • C, a ramp of 4 • C/min to 190 • C, a ramp of 6 • C/min to 275 • C, and a ramp of 15 • C/min to 320 • C. Carrier gas (He) flow was 2 mL min −1 .
Stable isotope values are expressed in per mil ( ) using delta notation: δX = [(R sample /R standard ) − 1] × 10 3 where X is 13 C or 2 H, R is the corresponding ratio of 13 C: 12 C or 2 H: 1 H, and R standard is Vienna Pee-Dee Belemnite or VSMOW for δ 13 C and δ 2 H, respectively. We analyzed the δ 13 C and δ 2 H values of 12 individual amino acids, 6 essential-valine (Val), threonine (Thr), leucine (Leu), isoleucine (Ile), phenylalanine (Phe), lysine (Lys)-and 6 non-essential-alanine (Ala), glycine (Gly), serine (Ser), proline (Pro), aspartic acid (Asp), glutamic acid (Glu)-listed in chromatographical order. Because acid hydrolysis converts glutamine (Gln) and asparagine (Asn) into Glu and Asp, respectively, by cleaving the terminal amine group, the measurements we report are combined Gln + Glu and Asn + Asp, hereby referred to as just Glu and Asp, to be consistent with previous studies. Tyr and Arg were excluded from analyses because of missing measurements caused by low concentrations near or below detection limits.
Amino acid standards of known isotopic composition were derivatized alongside and analyzed between unknown samples that were analyzed in triplicate. Reproducibility of replicate sample injections ranged from 0.1 to 0.6 (average < 0.3 ) for δ 13 C and 1.2-5.2 (average 2.5 ) for δ 2 H. Reproducibility of the amino acid standards used for sample corrections averaged < 0.6 for δ 13 C and 2.6 for δ 2 H. Because derivatization reagents add carbon atoms to the amino acids, the standards prepared alongside the samples are used to correct for the reagent carbon added. Similarly, exchangeable hydrogen atoms are removed from the carboxyl and amine groups of the amino acids during derivatization, and an additional hydrogen atom from the amine group is likely exchanged during hydrolysis. So, our measured δ 2 H data are corrected using amino acidspecific mass balance equations that adjust for the percentage of hydrogen atoms removed during derivatization (Fogel et al., 2016). Thus, the amino acid δ 2 H values we report reflect the nonexchangeable (intrinsic) hydrogen atoms only, and not hydrogen contributed by reagents during derivatization.

Statistical Analysis
We used Shapiro-Wilks and Levene's tests to confirm normality and homoscedasticity, respectively, and considered nonparametric equivalents before selecting ANOVA with Tukey HSD post-hoc tests to assess differences between tissue types in each individual amino acid. We further explored patterns in our amino acid δ 2 H and δ 13 C datasets with principal component analysis (PCA). AA Ess and AA NEss were evaluated separately. These statistical analyses were performed in R (version 4.0.3) with RStudio interface (version 1.3.1093; R Development Core Team, 2020) using the stats package.

RESULTS
Amino Acid δ 13 C The AA Ess δ 13 C of the four tissue types were distinct from one another, a finding very different from O'Brien et al. (2002O'Brien et al. ( , 2003O'Brien et al. ( , 2005's laboratory studies (Figures 2A,B). PCA of AA Ess δ 13 C clustered the samples according to their tissue type fairly well, with caterpillars the least variable and butterfly bodies the most variable ( Figure 3A). Caterpillars separated from pipevines and butterfly wings along the first (PC1) and second (PC2) principal component axes. Differences between pipevines and caterpillars were driven by Thr, Leu, and Ile, while differences between caterpillars and adult butterflies (especially wings) were driven by Val, Phe, and Lys. Higher Thr and Leu also distinguished butterfly bodies from pipevines. In the ANOVA, all AA Ess showed significant differences between some of the tissue types (Supplementary Table 5). Thr differed significantly between butterfly bodies and each of the other three tissues. Val differed for every pair except pipevines and butterfly bodies. Leu differed between pipevines and each of the other tissues. Ile was significantly different for half of the pairings (pipevines-caterpillars, pipevines-butterfly bodies, caterpillarsbutterfly wings). Phe and Lys differed between butterfly wings and every other tissue type.
With PCA for AA NEss δ 13 C, we found very tight clustering and excellent separation among pipevines, caterpillars, and butterfly wings. However, the two ellipses for butterfly tissues overlap ( Figure 3B) because bodies show a high degree of variation, producing an ellipse that stretches across the entire PC1 axis. Consequently, half of the butterfly body samples plot near the wings, while the other half plot below pipevines and caterpillars. PC1, and the separation between butterfly wings on one end and caterpillars and pipevines on the other end, is strongly influenced by Ala, Pro, and Asp. PC2 is strongly influenced by Gly and Ser, and accounts for even more of the separation FIGURE 2 | δ 13 C (A,B) and δ 2 H (C,D) data for essential and non-essential amino acids. Our average δ 2 H value (-42 ) for local water at the time of sample collection is indicated by the dashed gray line. among pipevines, caterpillars, and butterfly wings. Interestingly, the half of the butterfly body samples that plotted most closely with pipevines and caterpillars along PC1 plot the farthest from them along PC2. ANOVA indicated that every AA NEss except Glu showed significant differences among some of the tissue types (Supplementary Table 5). Ala and Gly differed between caterpillars and each of the butterfly tissues. Ser differed between butterfly bodies and the other tissues. Pro and Asp behaved similarly and only differed significantly between two pairings--pipevines-caterpillars, caterpillars-butterfly bodies.

Amino Acid δ 2 H
Similar to the δ 13 C data, AA Ess δ 2 H were distinct from one another, supporting the above finding that essential amino acids were not simply routed from plants to caterpillars to adult butterflies (Figures 2C,D, 4). PCA of AA Ess δ 2 H shows the tightest clustering of samples by their tissue types, with four nonoverlapping ellipses (Figure 3C). The two butterfly tissues show the greatest variability and the greatest similarity to each other. PC1 provides excellent separation of pipevines and caterpillars from both butterfly tissues. PC2 similarly provides substantial separation between pipevines and caterpillars. Differences between pipevines and adult butterfly tissues were mainly due to Lys, Leu, and Phe; Thr was most responsible for separating adult butterfly tissues from the caterpillars and plant samples, and Val and Ile drove the separation between caterpillars and butterfly bodies. ANOVA revealed that Thr differed significantly between all but two pairs-pipevines-caterpillars and the two butterfly tissues (Supplementary Table 5). Leu was significantly different between all pipevine pairs and between caterpillars and butterfly wings. For Ile, caterpillars differed from both butterfly tissues. Phe was significantly different for all pairings involving butterfly tissues. Lys was significantly different for all pipevine and butterfly wing pairs. Like AA Ess δ 2 H, all PCA ellipses for AA NEss δ 2 H are nonoverlapping ( Figure 3D). In contrast to δ 13 C, AA NEss δ 2 H varies most among pipevines and least among butterfly wings. PC1 mainly separates butterfly wings from pipevines and caterpillars and is controlled by Pro, Gly, and Glu. PC2 is governed primarily by Asp and Ser, and to a lesser extent, Ala. These three amino acids create the separation of caterpillars from the other three tissue types. Asp and Ser are also responsible for the separation between the two butterfly tissues. The pipevine ellipse is oblong, such that it stretches across most of PC2. Within this ellipse space, the pipevine samples form two clusters-four samples group at the top of PC2 to the left of the butterfly body samples while the other two samples are centrally located between caterpillars and butterfly wings. This suggests high variability among pipevine Asp and Ser. In ANOVA, Ala was significantly different in half of the pairings-pipevines-butterfly bodies, caterpillarsbutterfly bodies, and caterpillars-butterfly wings (Supplementary Table 5). Gly differed for every butterfly wing pairing. Ser was significantly different for all six pairs. Pro differed for every pairing except that of the two butterfly tissues. Like Ala, Asp differed significantly in half of the pairs-pipevines-caterpillars, caterpillars-butterfly bodies, and the two butterfly tissues. Glu only differed between caterpillars-butterfly wings, and the two butterfly tissues.

DISCUSSION
We found extensive isotopic variation within and between individuals and among tissue types despite an obligate caterpillarhost plant relationship, a C 3 herbivorous diet for both larval and adult B. philenor, and sampling from a single watershed (Figures 2, 3 and Supplementary Table 4). Our data reveal the extent to which physiological processes-particularly those associated with metamorphosis-contribute to isotopic variation and persist into the adult lifestage as carryover effects. Additionally, our results suggest that the gut microbiome contributes to the amino acid pool used by B. philenor for tissue synthesis throughout their life cycle. Understanding what drives δ 2 H variation among primary producers and herbivores is critical for isoscape development and δ 2 H-based animal migration studies.
Essential Amino Acid δ 13 C Because AA Ess cannot be synthesized by B. philenor, they should be obtained directly from dietary amino acids. Direct routing from dietary protein is assumed to result in minimal isotopic fractionation such that isotope values of AA Ess in caterpillars/butterflies and pipevines would be similar. In contrast, the use of AA Ess synthesized by microbes from hydrogen or carbon in dietary biochemicals would yield a mismatch in the isotopic composition of AA Ess between consumer and dietary protein. Adult B. philenor consume very little-if any-protein, which means their AA Ess primarily derive from larval resources (Scott, 1986;O'Brien et al., 2005;Lotts and Naberhaus, 2020). Thus, we expected the AA Ess δ 13 C values of B. philenor tissues to closely resemble each other and that of pipevines. However, our data show significant variation in AA Ess δ 13 C values among pipevines, caterpillars, and butterfly tissues (Figure 2A and Supplementary Table 3), which must be driven by foraging behavior and/or gut microbial contributions. One such foraging behavior is compensatory feeding, wherein insects increase food consumption-and frass production-to compensate for a low-protein diet (Lavoie and Oberhauser, 2004;Le Gall and Behmer, 2014). This foraging strategy could potentially lessen the severity of protein-limitation in B. philenor caterpillars but the extent to which it might have a carryover effect on adult fitness is unknown.
Alternatively, some nectarivorous butterflies combat protein limitation by supplementing their adult diet with amino acids, which can increase their fecundity or even lifespan (O'Brien et al., 2005;Otis et al., 2006;Cahenzli and Erhardt, 2013).
These supplements can potentially be acquired from nectar (Cahenzli and Erhardt, 2013;Levin et al., 2017), pollen (Erhardt and Baker, 1990;O'Brien et al., 2003), or mud-puddling, a behavior of aggregating at moist soil, tears, animal scat, or carcasses that is typically employed by males (Otis et al., 2006). However, amino acids are only present in trace amounts in nectar (Baker and Baker, 1973;Baker, 1975), B. philenor does not consume pollen (Konagaya et al., 2015), and mud-puddling is likely used to acquire sodium rather than protein (Otis et al., 2006). Moreover, even if amino acid supplements were contributing to the AA Ess pool in butterfly body tissues, this would not explain differences in δ 13 C values among pipevines, caterpillars, and butterfly wings. In four out of six AA Ess , it is either the caterpillar or the wing tissue that deviates most strongly from the δ 13 C fingerprint of pipevine leaves (Figure 2A). Thus, while B. philenor butterflies may opportunistically consume some amino acids from nectar, our data suggest that amino acid supplements do not substantially contribute to their AA Ess budgets.
A more likely explanation for observed differences in the δ 13 C of AA Ess between plant and consumer is subsidization by the gut microbiome. Although it remains unknown whether or not most Lepidopteran species host a gut microbiome, our data suggest that it contributes to the AA Ess pool used by both caterpillars and butterflies for tissue synthesis. How the microbial signature deviates from that of dietary AA Ess depends on whether gut microbes use carbohydrates or lipids-which are 13 C-enriched or depleted, respectively, relative to dietary proteins-as a carbon source to synthesize AA Ess . For example, gut microbes appear to synthesize Val, Leu, and Ile in caterpillars and Thr in butterflies from carbohydrates as the δ 13 C of these particular amino acids is substantially more positive (∼10 ) than those in leaf tissue (Figure 2A). Conversely, Val, Phe, and Lys in wing tissues seem to reflect larval gut microbial contributions from lipid stores during metamorphosis because these amino acids have more negative δ 13 C (∼7 ) than those in leaf tissue. We also found evidence of direct routing of AA Ess between pipevines and caterpillars and during metamorphosis. δ 13 C values of Leu are identical among the insect tissues and only ∼5 lower in pipevines (Supplementary Table 3).
Non-essential Amino Acid δ 13 C AA NEss δ 13 C values are known to reflect a combination of direct routing from dietary protein and de novo synthesis by consumers from dietary carbohydrates and lipids (Jim et al., 2006;Newsome et al., 2014). Because pipevine leaves and nectar have less protein content than animal tissue, we expected that caterpillars and butterflies would use other dietary macromolecules (e.g., nectar carbohydrates) and limited endogenous nitrogen stores to construct some of the AA NEss they need for metabolism and reproduction. In support of this prediction, we did not find a common AA NEss δ 13 C fingerprint among tissue types, indicating extensive metabolic remodeling of amino acids during each insect life history stage (Figure 2B and Supplementary Figure 4B). Glu was the only AA NEss that had a similar δ 13 C composition among tissue types (Supplementary Table 3); we hypothesize that this pattern results from Glu being the most abundant amino acid in animals, so its isotopic composition is relatively stable in comparison to amino acids that have smaller pool sizes. In contrast, Pro and Asp in wings are 13 C-depleted relative to caterpillars, suggesting that pupae, or their gut microbes, draw from lipid stores to synthesize these AA NEss for wing synthesis. Lipids are the most abundant energy reserve available during metamorphosis (Arrese and Soulages, 2010), and lipid carbon is easily introduced as acetyl CoA into the TCA cycle, from which Pro and Asp are derived (Horton et al., 2002).
Because adult B. philenor consume carbohydrates almost exclusively, it follows that amino acids synthesized from intermediaries in glycolysis (Ala, Gly, Ser) should reflect the strongest carbohydrate δ 13 C signature of all the tissue types. In support of this prediction, Ser δ 13 C is on average ∼20 higher in butterfly bodies than in the other three tissues. Of all the AA NEss measured, Ser has the most variable δ 13 C values (Figure 2B and Supplementary Table 3), which range > 50 among pipevines and butterflies and almost 40 in butterfly bodies alone (Supplementary Figure 3B). This is a substantial amount of δ 13 C variation for a single plantherbivore system and indicates extensive de novo synthesis, particularly by the nectarivorous adult butterflies. Ala δ 13 C values are significantly higher by ∼9 in caterpillars relative to pipevine leaves (Supplementary Table 5), also indicative of synthesis from carbohydrates. Interestingly, caterpillars have uniquely low Gly δ 13 C values, suggesting either a mass balance effect associated with metamorphosis or synthesis of this AA NEss with lipid carbon. Overall, isotope data for glycolytic amino acids demonstrate that the magnitude of 13 C-enrichment imparted by carbohydrate carbon is not universal. Rather, it likely depends on the interaction between the pool sizes and fluxes of each amino acid dictated by supply and demand of these metabolites, as well as the individual fractionation factors of the enzymes used to synthesize them. Essential Amino Acid δ 2 H AA Ess δ 2 H data further show that the hydrogen isotope composition of individual amino acids is systematically modified according to the biochemical pathways used for synthesis and that eukaryotes and many prokaryotes utilize the same basic metabolic machinery to do so (Rodriguez Curras et al., 2018;Newsome et al., 2020). As was the case with E. coli (Fogel et al., 2016) and house mice (Newsome et al., 2020) reared in controlled feeding experiments, Ile has the lowest δ 2 H values of all amino acids in our dataset. Ile also has the most uniform δ 2 H values, suggesting it was mainly routed directly from the larval diet into the insect tissues. Furthermore, the extreme variation we found in AA Ess δ 2 H within individuals and among tissue types ( Figure 2C and Supplementary Table 4) mirrors that seen previously in E. coli and mice. AA Ess δ 2 H varied over 300 in a single caterpillar ( Figure 4C) and variation was similar in the other tissues measured, even within a single pipevine leaf (∼230 ).
In every AA Ess other than Ile, δ 2 H was highest in either butterfly wings (Phe, Lys) or both wings and butterfly bodies (Thr, Val, Leu) relative to pipevines and caterpillars. This 2 Henrichment in adult insect tissues suggests that gut microbes, rather than direct incorporation from nectar, are contributing AA Ess to their hosts using some combination of hydrogen sourced from the metabolic breakdown of carbohydrates, and possibly from drinking water. Environmental water typically has higher δ 2 H values than all dietary hydrogen sources because of the large negative discrimination during the synthesis of organic compounds in photosynthesis (Estep and Hoering, 1981). The fact that AA Ess in wings are similar to or more enriched than the butterfly bodies points to body water as the predominant hydrogen source over carbohydrates because although stored glycogen can be accessed for glycolytic fuel during metamorphosis, without frequent feeding these carbohydrates are far less abundant than fat reserves and are quickly depleted (Arrese and Soulages, 2010). Although the adult gut microbiome may also play a role in butterfly body δ 2 H, carryover effects from the larval stage appear to persist throughout the butterfly's lifetime and are explored further below.

Non-essential Amino Acid δ 2 H
In agreement with the analogous δ 13 C data, AA NEss δ 2 H data suggest considerable de novo synthesis in B. philenor, particularly in Ser that ranges > 300 (-187 -137 ) among tissue types ( Figure 2D and Supplementary Table 4). Variation in Ser δ 2 H is not surprising given the essential function of serine proteases in insect metabolism-such as aiding digestion and development (Cao and Jiang, 2018)-which involves serine donating a hydrogen then gaining one from water (Radisky et al., 2006). Generally, structurally simple amino acids like Ser are more likely to exchange hydrogen with body water because a larger proportion of their hydrogen is bound to oxygen and nitrogen relative to more complex amino acids (Newsome et al., 2020). In controlled feeding studies on E. coli (Fogel et al., 2016) and mice (Newsome et al., 2020), drinking water appeared to have a particularly strong effect on the δ 2 H of glycolytic amino acidsthis was especially true when organisms were fed diets containing no or low protein contents. Our data indicate that the same is true of our protein-limited herbivores.
The influence of body water on tissue δ 2 H is likely amplified in a holometabolous insect like B. philenor. Inside its chrysalis, which the caterpillar must construct from protein, the pupa undergoes a complete morphological remodel. While wings are an iconic part of this transformation, the pupa must also shorten the gut, grow limbs, expand the tracheal system, reshape the mouthparts, and develop a hard exoskeleton (Lowe et al., 2013). This complex process requires energy, expends water, and generates waste, and it is unclear how differential routing of amino acids might affect fractionation of carbon and hydrogen isotopes among newly synthesized tissues. What is clear, though, is that there are myriad opportunities during metamorphosis for isotopic discrimination to occur.
Evaluation of over 50 Lepidopteran species revealed that pupae typically lose a large percentage (≤80%) of their weight during metamorphosis (Molleman et al., 2011). Most of the weight lost is water (Wappner and Quesada-Allué, 1996;Molleman et al., 2011). This substantial water loss likely results in 2 H-enrichment of the remaining body water pool, which is then used to synthesize 2 H-enriched amino acids for adult butterfly tissues. δ 2 H patterns in Ala, Pro, and many microbially synthesized AA Ess follow this prediction (Figures 2C,D). If wings are synthesized later in metamorphosis than adult body tissues, they may acquire higher δ 2 H values as a consequence of this gradual process, particularly if they are strongly influenced by hydrogen in body water. For instance, δ 2 H values of Gly, Ser, and Asp in butterfly wings are 40-50 more positive than in butterfly bodies (Supplementary Figure 4A and Supplementary  Table 4). Furthermore, prior studies have identified bulk tissue 15 N-and 2 H-enrichments following metamorphosis in multiple species, including Lepidoptera, and primarily ascribed them to the excretion of 15 N-depleted waste products and the inclusion of 2 H-enriched body water for tissue synthesis (Doi et al., 2007;Tibbets et al., 2008;Peters et al., 2012; see Supporting Information-Text for additional discussion of Pro).

Inter-Individual Variation in Amino Acid Isotope Values
We focus this discussion on δ 2 H, though we lean on δ 13 C for context, because understanding inter-individual δ 2 H variation among primary producers and herbivores is critical for isoscape development and further elucidates the influence of physiology on tissue δ 2 H values. Also see Supporting Information-Text for a brief discussion of mycoheterotrophy as a potential explanation for inter-individual variation in plant isotope values. We found that bulk tissue δ 2 H values of caterpillars are often less variable than those of pipevines (Figures 4C,D). Relative to all other amino acids, Asp has the most variable δ 2 H values in both pipevine and caterpillar samples; Asp δ 13 C values are also variable in caterpillars (Supplementary Figure 3C). In plants, aspartate is a precursor for asparagine, which is used to transport and store nitrogen, and a precursor for Thr, Ile, Lys, and methionine synthesis (Azevedo et al., 2006). In caterpillars, Asp is abundant in a storage protein (Telfer and Kunkel, 1991) that is catabolized in adult Lepidoptera and used as a critical substrate for reproduction (O'Brien et al., 2002). These key metabolic functions may explain the unusually high inter-individual variability observed in Asp. Despite these exceptions, the similarity in AA Ess δ 13 C and δ 2 H fingerprints among individual caterpillars implies relatively consistent fractionation during the plant-herbivore trophic transfer (Figure 4C and Supplementary Figure 3C). This uniformity strongly contrasts with the high degree of inter-individual variation among adult butterflies (Figures 4A,B and Supplementary Figures 3A,B). It is not surprising that caterpillars would be less variable than butterfly body tissue, given the generalist diet of the adults and the obligate host plant relationship maintained by the larvae. However, the fact that caterpillars also tend to be less variable than wings-which are largely metabolically inert-emphasizes the dominant role played by physiology during metamorphosis in determining δ 2 H values of adult tissues.
Concurrently assessing both δ 13 C and δ 2 H fingerprints for individual butterfly bodies reveals an additional interesting feature of our dataset that is best illustrated by the highly variable Ser and Thr: the relative rankings of individuals' 13 C-enrichment are fairly consistent independent of amino acid (Supplementary Figure 3B). That is, the three butterflies with the highest Ser δ 13 C also have the highest Thr δ 13 C. The same is not true for δ 2 H ( Figure 4B). This is likely a reflection of the increased complexity of hydrogen cycling relative to that of carbon, but perhaps it also signifies a promising opportunity. δ 2 H analysis of specific subsets of amino acids may be best suited to answering different kinds of ecological or physiological questions. For example, glycolytic amino acids have previously been proposed as the most faithful proxies for δ 2 H of drinking water (Newsome et al., 2020) and so may be most beneficial for assigning geographic origin to migratory animals, particularly when bulk tissue δ 2 H is insufficient. Because caterpillars are an important food source for many higher trophic level organisms, their AA Ess δ 2 H (and δ 13 C) values could be useful dietary tracers. Moreover, AA Ess δ 2 H and δ 13 C provide a way to assess the presence and role of more complex metabolic processes like gut microbial input in systems where the isotopic composition of diet is well characterized. Several microbially synthesized AA Ess as well as Ala, Pro, and Asp may be particularly valuable for studying physiology and biochemistry during metamorphosis, including the pupa's use of a 2 H-enriched body water pool to synthesize AA NEss or AA Ess (via its microbiome), amino acid routing, and the order in which adult tissues are synthesized.

CONCLUSION
The use of δ 2 H data in movement and migration studies is built on the premise that bulk animal tissue δ 2 H is primarily determined by precipitation δ 2 H. The presumption being that δ 2 H signatures are transferred from meteoric water to primary producers to animal consumers in a systematic fashion such that δ 2 H is a faithful intrinsic record of tissue synthesis location. Yet, our data show that bulk tissue δ 2 H values can range over 30 within a single population of pipevine-feeding caterpillars and amino acid δ 2 H values can vary by 300 within and among these herbivorous insects collected from a single watershed over a narrow elevation range. Given the absence of geographical, dietary, or temporal variation among our study organisms, this is a clear indication that physiology drives δ 2 H and δ 13 C fractionation during the trophic transfer from primary producers to caterpillars and especially during metamorphosis. Insects are a major herbivore group in terrestrial food webs and a common food source for consumers, including birds that are frequently of interest in δ 2 H-based studies of avian migration. Therefore, understanding how δ 2 H and δ 13 C values of insect tissues are biochemically mediated and vary among individuals is particularly critical for elucidating how δ 2 H propagates through food webs. Bulk tissue δ 2 H can be an invaluable tool for answering questions regarding breeding and overwintering grounds that are separated by thousands of kilometers along steep and well-established gradients in precipitation δ 2 H. However, further investigation into amino acid δ 2 H may yield a complementary tool that can be more routinely applied to studies of bird and insect movement across smaller geospatial scales or to better understand the physiological processes associated with metamorphosis in insects.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
SN, GG, and MF formulated the project idea and collected samples. KM performed stable isotope and statistical analyses and led the writing of the manuscript. All authors engaged in intellectual exchange and manuscript editing and gave final approval for manuscript submission.

FUNDING
Funding for this project was provided by an NSF award to SN (DEB-1343015) and MF (DEB-1437845).