Plant proteins of host plants are an important nutrition source used by tephritid flies

Volatile chemicals stimulate chemosensory receptors in tephritid flies when assessing a potential novel host and trying to expand . Therefore, chemosensory-related genes are involved in the initial process of host plant expansion for tephritids. Olfactory-related genes of tephritids are one type of chemosensory gene that includes several gene families of odorant-binding proteins , chemosensory proteins , odorant receptors , ionotropic receptors , and sensory neuron membrane proteins , which are primarily involved in the identification of volatile chemicals, including volatiles of host fruits. After receiving odor chemical signals, these olfactory-related genes are triggered to transduce cascades that send information to specific regions of the brain, which ultimately leads to specific behavioral responses . OBP genes play an important role in the first step of chemosensory identification of insects, including tephritids . OBP genes direct odorant-binding proteins to bind volatile odor molecules specifically by distinct expression to related olfactory receptors that are bound to olfactory receptor neurons in antennae . CSP genes are regarded as playing a similar role as OBP genes involved in the initial process of chemosensory signal transmission to corresponding receptors . OBP and CSP genes are major gene types that lead tephritid flies to respond to different chemosensory chemicals, including volatile chemicals of host plants . Except for these two categories of genes, large plastic growing pots some odor receptor genes also play important roles in host odor recognition of tephritids, such as genes related to odor receptors and ionotropic receptors .

Odorant receptors of insects are composed of at least two proteins: a conserved coreceptor as an ion channel and a specific OR subunit , which determines the ligand specificity and forms structurally ligand-gated ion channels . The OR genes mediate odorant receptors of insects transmitting the odorant molecules they receive into electric signals that are transmitted to a higher-order neural center . IR genes are related to ionotropic glutamate receptors , which are regarded as ion channels . They also play important roles in odor chemical perception . The sensory neuron membrane proteins gene encodes transmembrane domain-containing proteins that belong to a large gene family of CD36 receptors . SNMPs regulates the corresponding proteins to identify chemosensory signals, mainly pheromone chemicals . The GR family is another type of chemosensory protein that is a ligand-gated ion channel broadly expressed in gustatory receptor neurons in taste organs and is mainly involved in taste recognition of CO2 , sugar, and bitterness . When receiving taste signals, GR genes are involved in identifying taste and ingestion. Among tephritid flies, Bactrocera dorsalis and Ceratitis capitata are well-known polyphagous species that have expanded their host plants to more than 250 species . However, Bactrocera minax and Z. cucurbitae are oligophagous species that mainly attack citrus fruits and cucurbit plants, respectively. Bactrocera oleae , Procecidochares utilis , and Carpomya vesuviana are monophagous species infesting olive , crofton weed , and jujube , respectively, and all have limited host plant species . Compared to several major olfactory-related gene families, the two polyphagous species have more genes, with 3 CSPs, 35 OBPs, 74 ORs, and 40 IRs in B. dorsalis and 45 OBPs, 76 ORs, and 70 IRs in C. capitata , than two host-limted species .

A similar situation was observed in the GR family. There are also more GR genes in C. capitata and B. dorsalis than in the host-limited species P. utilis , C. vesuviana , and Z. cucurbitae . The increased numbers of these genes are associated with chemosensory-related gene family expansion via gene duplication and differentiation , which exertimportant roles in tephritid fly adaptation to other hosts and expansion of their host ranges. Obvious chemosensory-related gene expansions were also reported in Tribolium castaneum , Spodoptera frugiperda , and Heliconius melpomene . For example, the pea aphid Acyrthosiphon pisum , with broader host ranges, experienced obvious expansion of the OR, OBP, and GR gene families, with 87 ORs, 18 OBPs, and 78 GRs, compared to the soybean aphid Aphis glycines , with 47 ORs, 10 OBPs, and 61 GRs . Altering gene expression levels also helps tephritids respond to different host plants and realize host expansion. OR13a and OR82 expression are higher in antennae in B. dorsalis in response to 1-octen-3-ol and geranyl acetate, respectively, which are major volatile components of its host fruits, mango and almond fruit . For B. minax, increasing the expression levels of several GR genes regulate the taste process in response to different chemosensory stimuli of hosts .Once a tephritid adult identifies a potential novel host fruit for oviposition or feeding, the plant fruit must be suitable for larval development, which includes overcoming any secondary toxic chemicals in the novel host fruit . Therefore, detoxifcationand other digestion-related genes also play core roles in mediating the host plant expansion of tephritids. Common detoxifcation-related genes of insects include gene families of cytochrome P450s , glutathione S-transferases , UDP-glycosyltransferases , carboxyl/cholinesterases and ATP binding cassettes .

The cytochrome P450 family belonging to phase I enzymes includes various CYP subfamilies for different tephritid species . The GST superfamily consists of phase II enzymes divided into at least seven major subclasses: the delta, epsilon, omega, sigma, theta, zeta, and microsomal classes . The PGE phase II enzymes are a large family that can be divided into 13 clades, including the dietary detoxification class , the hormone/semiochemical processing group , and the neurodevelopmental group . The ABC transporter superfamily belonging to phase III enzymes can be subdivided into eight subfamilies, from ABC-A to ABC-H. The cytochrome p450 gene family of phase I mainly contributes to the catalysis of numerous oxidative reactions during endogenous and exogenous metabolism . The important roles of genes in this family are the metabolism of xenobiotics, plant allelochemicals , and even insecticides. GSTs are multifunctional genes of phase II enzymes that play a crucial role in the detoxification of endogenous and xenobiotic compounds, including plant secondary metabolites and pesticides. CCE families of phase II have been shown to be involved in the detoxification of plant-derived allelochemicals as well as insecticides . The ABC transporter genes of phase III encoding membrane-bound proteins typically function in the ATP-dependent transport of various substrates across biological membranes . The roles of ABC genes are mainly in handling xenobiotics such as plant phytotoxins and insecticides . These genes can participate in regulating detoxification of host plant secondary metabolites of tephritid flies by coding corresponding enzymes, which help to transform toxins entering the insect system into hydrophilic compounds that can be eliminated and in the adaptability of different hosts . The major digestive-related genes include gene families of cysteine proteases, proteases, lipase, glucosidase, and serine proteases . The serine proteases are members of the supergene family, including chymotrypsin, trypsin, thrombin, subtilisin, plasmin, and elastase. subclasses . Various digestive proteases exert important roles in the nutrition digestion of tephritid flies from novel host plants that they try to expand to. However, protease inhibitors of host plants are a widespread defense against herbivores such as tephritids. Therefore, genes coding various proteases react to protease inhibitors by regulating inhibitor-sensitive proteases or expressing proteases that are not targets of the inhibitors . When expanding to other novel hosts, large plastic pots tephritid flies must adapt to different chemical environments from their native hosts. Detoxification-related genes regulate the host expansion of tephritids via gene family expansion similar to chemosensory-related genes. The major gene families of detoxifcation GSTs, P450s, CCEs and ABC transporters are more numerous in polyphagous B. dorsalis and C. capitata than monophagous P. utilis and B. oleae . However, reports about digestive gene family expansion in tephritids are still rare. Overall, detoxifcation and the digestive-related gene family combined with chemosensory-related gene family amplifcation exhibit a close association with host range extension. This gene family expansion is helpful for the host plant expansion of fruit flies. Cases in other insects strengthen this idea. For example, Helicoverpa armigera and Helicoverpa zea are two species of caterpillars that have considerably broader host ranges than any other lepidopterans.

Great expansion of detoxification and digestive gene families was found in the two species. In addition to gene family amplification, detoxification and digestive genes also regulate host expansion of tephritid flies by activating various gene subfamilies, subclasses, or clades. To respond to various toxic environments, including secondary toxic chemicals of different hosts, B. dorsalis primarily triggered the delta subfamily of GSTs, CYP3 and CYP4 subclasses of P450s, A–C clades of CCEs, and ABC-A, ABC-B, and ABC-G subclasses of ABC transporters , C. ceratitis activated the epsilon subfamily of GST, CYP6 and CYP12 of P450s, B clade of CCEs , and P. utilis mainly triggered the delta, epsilon and microsomal subfamilies of GSTs, CYP4, and CYP9 of P450s, C clade of CCEs and ABC-G subclass of ABC transporters , but R. pomonella mainly launched CYP4 and CYP6 of P450s . For the digestive gene family, B. dorsalis and C. capitata primarily triggered aminopeptidase, trypsin and serine peptidase digestive genes, but B. oleae, which is a strictly monophagous species, triggered serine protease and nuclease digestive genes to respond to different host secondary chemical environments . Detoxifcation- or digestion-related genes also facilitate tephritid fly adaptation to different hosts by altering gene expression levels. Rhagoletis zephyria evolved fromRhagoletis pomonella and experienced host expansion from apple to snowberry plants . Increased expression levels were found in some detoxifcation-related genes, including cytochrome P450, glutathione S-transferases, and glycosyltransferase, in R. zephyria facing the apple host environment . Z. cucurbitae is the species that mainly attacks cucurbit plants, and the fly responds to different secondary chemical environments of Mucuna pruriens plants by reducing the expression levels of trypsin and chymotrypsin digestive genes .Although the importance of chemical stimuli is highly emphasized in the host expansion of tephritids, other nonchemical stimuli, such as the color of the novel host fruit, should not be ignored. Many insects locate their host plants primarily by color signals, including beetles, Altica engstroemi , Hylastes ater , and Arhopalus ferus . For tephritid flies, Neoceratitis cyanescens , B. minax , B. dorsalis , and Z. cucurbitae are typical examples of species that appear to select different hosts frst by fruit color rather than chemical signals. Z. cucurbitae realized its host expansion to a novel host, papaya , in Hawaii by strongly relying on the color location of fruits by vision . However, the underlyinggenetic mechanism has not been revealed for Z. cucurbitae. In fact, nonchemical stimuli, such as color, are associated with vision-related genes that allow the identifcation of different hosts . The genes responsible for color discrimination in Diptera are primarily related to opsin proteins in the photoreceptor cells of the eye . Six types of Rh opsin-expressed genes have been identifed as major genes involved in color recognition and photoreception in Diptera insects. The Rh1 and Rh2 opsin genes are associated with motion detection and direction, respectively . Rh3 and Rh4 are UV-sensitive opsin genes, Rh5 is a blue-sensitive gene and Rh6 is a green opsin gene . These opsin genes lead the photoreceptor of eyes to receive various chromophore pigments and then activate a series of visual transduction cascades to launch corresponding color identifcation behavior. In the genome of polyphagous C. capitata, the long wavelengthsensitive genes Rh1, Rh2, and Rh6 and the UV-sensitive genes Rh3 and Rh4 were found, while Rh2-4 and Rh6 were found in the phototransduction pathway of oligophagous B. minax . Moreover, the role of Rh6 in modulating green color discrimination was reported in C. capitata and B. minax . In B. minax, the function of Rh6, which is responsible for green spectral sensitivity, has been identified by knockdown of the gene B. minax in female adults, and B. minax flies significantly reduced their preference for green fruit after cutting Rh6 . Absence of a member of the bluesensitive opsin subfamily was found in both tephritid species C. capitata and B. minax, but Rh5 can be specifically expressed in Drosophila . Reports about vision-related genes directly involved in the host expansion of tephritids are still very few.Tephritid fruit fly hosts expand to other new host plants, and the phenology of the new host is another nonchemical stimulus that affects fly adaptation. The phenology of the host plant fruits includes the time of flowering, fruiting, or maturation .

Bulk density is a measurement of the weight of a volume of soil

The Center’s latest research brief, Farming the College Market: Results of a Consumer Study at UC Santa Cruz, by CASFS social issues associate specialist Jan Perez and social issues specialist Patricia Allen, examines student, staff, and faculty’s food-related concerns, interests, and level of support for specific food criteria, including whether they would be willing to pay more for food produced in an organic and “socially just” manner. Write the authors, “Since the success of farm-to-college programs involves their ability to meet the needs of campus consumers, we undertook a study of our local campus . . . to learn about the perspectives and preferences of campus food consumers.” Perez worked with members of the campus’s Food Systems Working Group, including UCSC Dining Services, Community Agroecology Network, and Students for Organic Solutions, to develop a web-based survey designed to find out what the UCSC community thinks about food system issues. Those responding to the survey included students , staff , and faculty . Based on the online survey, the study found that food issues are important to the UCSC community, particularly regarding concerns for the environment and for people. Key points from this study include: There is significant interest in campus food that is nutritious, safe, supports workers, and is environmentally sound; interest in local food and GE-free food is lower. People are interested in sustainably produced food and a majority of people already purchase food with labels based on these criteria. Many people are willing to pay more for food that meets social justice criteria. A campus community is likely to be receptive to education and discussion about food-system issues. Since people had a great interest in nutrition and food safety, square pots for plants framing discussions of food-system issues in terms of health will meet people’s needs as well as capture their attention for education on other food-system issues, such as working conditions and the environment.

The authors conclude that, “[Although] It would not be appropriate to extrapolate too much from a study of one campus . . . the results of the UC Santa Cruz study support the idea that colleges and universities are excellent choices for developing farm-to-institution programs and for popular education on food-system issues.” A similar survey was recently distributed to a nationwide audience as part of CASFS research on farm-to-institution programs .Finding a non-fumigant alternative to the soil fumigant methyl bromide has been identified as a top priority by the California Strawberry Commission, and by growers who are facing the phase out of this ozone-depleting pesticide. However, most current state and federal research is focused on alternative fumigants rather than non-fumigant techniques to control soil diseases, weed seeds, and harmful nematodes. Preliminary research conducted at the UCSC Farm has shown promising results with an alternative approach that starves pathogens and weeds of oxygen. Researchers introduce a carbon source such as chopped cover crops, wheat bran, or molasses to the strawberry bed, then irrigate and tarp the beds to create temporary anaerobic conditions. This technique, known as anaerobic soil disinfestation , has been tested for the past several seasons at the Farm and has been shown to control the soil pathogen Verticillium dahliae, a major diseases of strawberries. A new study, funded by a grant from the US Department of Agriculture, will expand this initial work to examine the efficacy of various carbon sources, irrigation techniques, tarp types, and tarping periods to create sufficient anaerobic conditions to control weed seed germination and V. dahliae. The study, conducted by environmental studies professor Carol Shennan, UCSC researcher Joji Muramoto, and colleagues from California and Florida, will also look at using the technique to control diseases, pests, and weeds in Florida cropping systems, which also rely on methyl bromide. The three-year study will take place on conventional farms so that results can be compared to those obtained with standard methyl bromide fumigation.

Following an initial year of field trials in Watsonville and Oxnard, the research team will consult with local growers to determine which ASD options hold promise for larger-scale commercial application.Patricia Allen has been named director of UCSC’s Center for Agroecology and Sustainable Food Systems by Social Sciences Dean Sheldon Kamieniecki. Allen, whose appointment took effect July 1, 2007, had been serving as acting director since January 2007. She joined the campus’s Agroecology Program in 1984. She takes the helm from Carol Shennan, CASFS director for the past 10 years. Following her 2007 sabbatical, Shennan will continue at UCSC as a professor in the Environmental Studies Department. Allen is one of the nation’s leading scholars on the social aspects of sustainable food systems. Her work addresses issues such as labor, gender, and access to food. She is the author of Together at the Table: Sustainability and Sustenance in the American Agrifood System and editor of Food for the Future: Conditions and Contradictions of Sustainability . Allen earned a B.S. in political economy of natural resources from UC Berkeley, an M.S. in international agricultural development from UC Davis, and a Ph.D. in sociology from UCSC in 1998.In the world of sustainable agriculture, it doesn’t get any better than the “Sustie” award, and the UCSC Apprenticeship in Ecological Horticulture took home the top honor at this year’s Ecological Farming Conference. Established in 1988, the “Sustie” award is presented each year by the Ecological Farming Association to “stewards of sustainable agriculture” who have made a significant contribution to the well-being of farming and the planet. Past recipients include chef Alice Waters, publisher Robert Rodale, and several graduates of the apprenticeship itself. UCSC Farm manager Jim Leap and apprenticeship coordinator Diane Nichols accepted the Sustie on behalf of the apprenticeship during the conference’s awards banquet on January 26 at the Asilomar Conference Grounds in Pacific Grove. “There are more than 25 extremely motivated and dedicated individuals who are instrumental in making the training what it is each year,” said Leap. “All of us work collectively to teach and train and run the UCSC Farm, and it is all of us who will be sharing in the acknowledgment that this award represents.” The apprenticeship, which is celebrating its 40th anniversary this year, is the nation’s premier hands-on training program in organic farming and gardening. Widely regarded as one of the most significant influences in the growth of sustainable agriculture, the six-month full-time program has prepared more than 1,200 graduates who have spread their expertise throughout the world. “There’s simply nothing that compares to the apprenticeship for the depth of its program or the breadth of its impact,” said Sheldon Kamieniecki, dean of the Division of Social Sciences at UCSC, who attended the awards banquet. Graduates of the apprenticeship go on to operate commercial farms and market gardens, run community and school gardens, and work at the forefront of international development, food policy, and social justice programs. The impact of the apprenticeship is apparent in the number of graduates who have received Sustie awards, including Cathrine Sneed of The Garden Project in San Francisco; Wendy Johnson, garden manager at Green Gulch Farm in Marin County; Jim Nelson of Camp Joy Gardens in Boulder Creek; Gloria and Steven Decater of Live Power Community Farm in Covelo, CA; Orin Martin, manager of the Alan Chadwick Garden at UCSC; and Kay Thornley, who helped launch UCSC’s Agroecology Program.

Its success is also evident in the number of similar college-based farm-training programs sprouting up at the University of Georgia, Michigan State University, Washington State University, the University of Montana, and other campuses.Perhaps an overview statement might be helpful—let’s define terms. The French intensive, raised bed style of gardening is a handworked system featuring deep cultivation . This technique’s primary effect is on the physical properties of a soil: the aim is to rapidly improve soil structure and fertility. Improved physical properties can positively influence the biological and chemical properties of a soil as well. The main idea is to create a well-drained, well-aerated, large square plant pots fertile soil structure by digging deeply and placing nutrients at specific levels. This gives rise to a profile that enables plant roots to probe/penetrate throughout the bed with ease, especially in a downward direction. Such an arrangement has a continuous system of large and intermediate pore spaces from the surface to the subsoil. Plants’ needs for air, water and nutrients are best met with such a continuous system of pores. The French intensive system is not appropriate in all soils and in all climatic situations. For instance, on deep, improved soils, it’s superfluous, even deleterious. On sandy soils and in hot, windy situations it can “burn up” precious organic matter and cause water losses both through surface evaporation and excessive drainage. As is so often the case in life, there are no panaceas, but we tend to be creatures of habits, creatures of dosages; that is, we want to do the same thing in the same way, with the same amount, repeatedly. The judicious use of deep digging for a few years to develop a soil, followed by lighter, less disruptive surface cultivation and perhaps periodic renewal via deep digging again might be more prudent. Caveat emptor: Digging is a radical act, potentially destructive of soil structure and biological processes Do it skillfully and as infrequently as possible! Conventional wisdom often states that it can take 1,000–2,000 years for 1 foot of topsoil to develop in place. With French intensive it is possible to simulate the creation of 1 foot of topsoil in 3–5 years .Often when people hear French intensive, they automatically think of raised beds. In fact, the beds may be raised slightly or in an exaggerated sense , flat, or even sunken. The degree of “loft” is a function of climate, soil type, and seasonal weather. On a transect from Seattle to Santa Cruz to Santa Fe, the response might be: 1) Seattle, with its high annual rainfall and cool temperatures, can have dark soils with high organic matter and high clay content, and a tendency to remain cold and wet. Thus a raised bed would yield better growth, allowing the soil to warm more quickly. Santa Cruz, with its mild Mediterranean climate, dry summers and wet winters would feature a slightly raised bed during the rainy season and an almost flat bed in summer. Santa Fe might yield a flat or even sunken bed for water catchment, to minimize water loss and afford protection from wind. Permanent beds, be they raised or flat, substantially reduce soil compaction. The bed equals the zone of maximum fertility—you could say “Don’t tread on me,” or only minimally and lightly. The path equals the zone of degradation, with much foot traffic and resultant compaction. Permanent beds foster maintenance of ideal soil structure. While compaction is a primary problem in mechanized agriculture, it can be virtually eliminated in handworked permanent bed systems. In agriculture, it can be said that the back of the tractor is simply undoing the work of the front of the tractor. Some common causes of soil compaction are: ploughing—a “plow pan” develops just below the depth of tillage; 2) machine and foot traffic —the bigger the machine, the greater the number of passes, the greater the compaction; 3) the pounding action of rain drops on open soil, which can destroy surface soil aggregates and lead to crusting and erosion. Natural forces also cause compaction—over time, the fine particles of clay leach downward, accumulate in layers, and create subsurface compaction or a hard pan. Compaction can be measured by an increase in bulk density. It includes pore space as well as solids. It is distinct from particle density, which simply measures the weight of a soil as if there were no pore spaces. Permanent beds also focus efficient placement of fertilizers/nutrients only where plants will be growing.The Apprenticeship in Ecological Horticulture is one step closer to finally building permanent apprentice housing on the UCSC Farm, thanks to long-time Friends’ member Olivia Boyce-Abel. When the campus approved the plans for eight 4-room cabins, the price tag that came with it was $487,000. Boyce-Abel not only pledged $40,000 from her obaboa Foundation for the project but also offered to help inspire other former apprentices and program supporters to give. She has put out a challenge that so far has more than matched her $40,000, with a total of $85,000 raised to date.

The primary techniques used are: single digging and double digging

But the fence removed herders and dogs as well as predators, and over the course of the summer, “the sheep gradually became accustomed to free, unmolested grazing, and forgot the habits learned when herded” . They formed smaller bunches, both while grazing and also to sleep at night, and they moved more lightly over the landscape. “The result was that little or no damage was done to the forage crop in this way. The entire crop was eaten and not wasted” .By weighing twenty lambs “of average size” at the beginning and end of the experiment, and collecting similar data from bands herded on the range nearby, Jardine was able to compare animal performance under the two systems. In 88 days inside the fence, the pastured sheep gained an average of 20 pounds, whereas those from a band on the outside gained only 15 pounds, on average, in 96 days. Overall losses among four outside herds ranged from 1.4 to 3 percent, compared with just 0.5 percent for the fenced herd. Finally, the fenced sheep required only 0.0156 acres per sheep per day. Three open range herds, roughly estimated, required 64-123 percent more acreage, which Jardine attributed to the effects of trampling. “It is safe to conclude that range grazed under the pasturage system will carry 50 per cent more sheep than when grazed under the herding system, where the band is driven to and from camp each day” . In the closing pages of his report for the 1908 season, Jardine took up the economic question: “Will the pasturage system pay?” Here he faced a difficulty, because the coyote-proof pasture had been very expensive to construct: $6,764.31, to be precise, blueberry package including more than $2,000 in materials, $1,000 in transportation costs, and $1,000 to clear heavy timber from portions of the fence line .

This amounted to nearly $850 per mile . Jardine chose not to use these figures, however, arguing that the location was exceptionally remote and heavily timbered. Using more general estimates, he calculated that “the cost on most grazing lands will approach very closely $400 amile” . He then tabulated the financial benefits in increased carrying capacity, heavier sheep, reduced losses to predation, and lower labor costs. Not counting increases in the lamb crop and in the amount and quality of wool , Jardine arrived at an annual return of $746.50, based on a herd of 2,200 sheep in a 2,560-acre pasture for three months. Thus, an initial investment of $3,200 , at 8 percent interest and including maintenance, would pay for itself and begin yielding dividends after six years. He did not include the costs of the hunter in his analysis. The experiment was repeated in 1909 with 2,040 sheep enclosed in the pasture for 99 days, during which time only four perished. The results were nearly identical to those from 1908, although Jardine’s report was more detailed and more emphatic in its declarations of the virtues of “the pasturage system.” The hunter killed one grizzly, two badgers, seven brown bears and seven coyotes , and just one brown bear and three badgers managed to breach the fence . The sheep again displayed a gradual tendency “to depart from their old habits and accommodate themselves to the freedom of the pasture,” so much so that by the end of the season “it was almost impossible to keep them close bunched without a dog” . Losses to poisonous plants, as well as predators, were higher among sheep herded outside the enclosure, and the pastured sheep were again heavier at the end of the season than the herded sheep .

Acreage required per sheep perday was 52-90 percent higher outside the fence, although Jardine attributed some of this to poor quality herders . The issue of herder skill presented a puzzle, which Jardine acknowledged but failed to address directly. “A first-class herder will work all the time” and seldom use a dog, resulting in “quiet, scattered grazing that may approach the pasturage system in efficiency” , whereas “a lazy man…will wear out his dogs, worry the sheep, and destroy the forage” . With respect to weight gain, “there is as much difference in the results obtained by a first-class herder and those obtained by a poor herder as there is between the results under the pasturage system and those secured by the good herder” . And while “range grazed under the pasturage system will carry from 25 to 50 per cent more sheep than when grazed under the herding system,” it was also possible “that an excellent herder can, to a considerable extent, allow his sheep freedom and keep them quiet, thereby increasing the carrying capacity of his range. No doubt there are herders who do this” . All told, “the carrying capacity of the same range utilized by different herders may vary at least 25 per cent” . These were potentially troublesome admissions to make, for both scientific and economic reasons. Herder skill was clearly an important variable in sheep performance, but it was one that Jardine could neither measure nor control. Removing herders might thus be seen as necessary to a properly “scientific”assessment of range grazing. And it could clearly affect any calculation of the economic rationality of building fences and controlling predators. What if better training for herders were a more economical solution?

Instead of confronting these issues, Jardine reverted to general claims about labor needs under the pasturage system, superseding by half the estimate that King had given to Coville two years earlier: “It is probable that one energetic man…can properly care for four inclosures similar to the experimental coyote-proof pasture,” meaning “one man would care for from 8,000 to 10,000 head of sheep” . Notwithstanding these problems, Jardine reached the same conclusions in 1909 as he had the year before. “When left unmolested by herders and dogs in an area protected against destructive animals, a band of ewes and lambs will accommodate themselves to the freedom of such a system and will separate into small bunches, coming together occasionally but again separating. With few exceptions they will graze openly and quietly” . Losses will be slight, weight gain and wool clip will improve, carrying capacity will increase by 25-50 percent, and labor costs “will not exceed 25 per cent of the cost under the herding system” . Although he did not present numbers for 1909, he again concluded that “the increase in carrying capacity and decrease in expense of handling in pasture during the lambing season will justifly the cost of construction necessary to inclose [sic] the entire allotment” . The Wallowa experiment was hailed as a remarkable success, and the Forest Service quickly embraced it as guidance for the administration and management of rangelands generally. It appeared to solve numerous problems and satisfly everyone, provided one ignored or excluded the herders. In his “Annual Report to the Forester for the Fiscal Year 1909-1910” for the Branch of Grazing, Potter described the results as “very gratiflying” and summarized them as follows: “The primary objects of the experiment have been accomplished, i.e., it has been demonstrated that the grazing capacity of the Forest lands can be largely increased by improved methods of handling stock, and that the increased cost of such methods, if any, is offset by increases in the number and weight of lambs raised, heavier wool crops, and reduced losses from predatory animals” . Notably, Potter omitted any reference to labor costs in his summation. He suggested that the results be applied “to spring and fall or yearlong ranges” in other national forests. Elsewhere in his report, Potter tabulated the accomplishments of Forest Service personnel assigned to predator control: 269 bears, 129 wolves, 148 wolf pups, 1,155 wildcats, and more than 7,000 coyotes killed in the 11 western states, “an increase of 109 per cent over the number of animals destroyed last year” and representing “a total saving to stockmen of considerably more than one million dollars per year” .A critical physical geography of the Coyote-Proof Pasture Experiment reveals two important insights. First, blueberry packaging the institutional context in which science is practiced may be at least as important as the experiments and findings that the scientists produce. Livingstone has shown the importance of place to the conduct of science, and the geographical particulars—both social and ecological—of the Wallowa site attracted Coville to locate the study there. But the larger context was a national one, and in this case, not only did it shape the questions that were asked and define the terms of success, but the success caused the institutional context itself to change, setting in motion a path dependency for subsequent research and management of rangelands elsewhere.

The particularities of the Wallowa site were abstracted away so that the results could be taken as relevant to rangelands throughout the western US. Second, although the findings of the experiments were presented as scientifically robust, ecological facts—in which predators, livestock, fences, and vegetation interacted in measurable ways—the methods were weak and the ultimate metric for evaluating success was actually an economic one. Costs and returns, and thus profit on investment, determined whether fencing and predator control were worth implementing. In this calculus, the decisive factor was neither fences nor livestock performance but rather the labor of herders. The high cost of fencing could be justified economically only if the fences greatly reduced the need for herders—and in the absence of herders to protect livestock, predators would have to be rendered effectively insignificant. What’s more, even the economic analysis was flawed, with key costs either minimized or excluded in order to reach the desired conclusion. By today’s scientific standards, the Coyote-Proof Pasture Experiment was far from impressive. The methods varied in several ways from one year to the next, and they were reliant on qualitative assessments or herder accounts for some important data. The findings were confounded by sheep breeds, herder practices, vegetation types and other variables, and no direct attempt was made to assess the impacts of the pasturage system on vegetation. Moreover, there were no real controls by which to judge the relative effects of the enclosed pasture, the absence of predators, and the absence of herders or dogs as factors influencing the dependent variables . Estimating the carrying capacity outside the fence was imprecise, since those herds were not confined within fixed boundaries. Finally, the economic calculations excluded the costs of the hunter, and didn’t account for the actual costs of the fence. Virtually every finding looks suspiciously similar to the expectations that Coville carried into the experiment in 1907. Moreover, Jardine’s reports were not subject to peer review, but only to the scrutiny of his superiors in the USDA, who themselves appear to have pre-judged the results. There is almost no chance that the experiment would be recognized as publishable, or even scientific, if it were conducted today. The experiment succeeded not on the basis of its scientific rigor, but instead because it lent authority to ideas that were already viewed favorably within the institutional context that gave rise to it. Coville and Jardine produced a set of knowledge claims that appeared to conform to scientific norms of experimentation: the deliberate manipulation of objects, organisms and people, and the careful recording and interpretation of actions and reactions among them. Most of these primary data were of a broadly ecological character, and thus appeared as apolitical and “objective.” The results were translated into economic terms to assess the practicability of implementing a similar management regime on Western rangelands as a whole, and if the economic analysis was at once flawed and decisive, this contradiction would be resolved by government largesse: the cost of both fencing and predator control would be subsidized by an array of federal agencies over the decades to come. Predator control occurred throughout the national forests and beyond, much of it organized and funded by the Bureau of Biological Survey under federal legislation passed in 1914; fencing was underwritten by the Forest Service, the General Land Office’s Grazing Service, and the Civilian Conservation Corps, whose crews built thousands of miles of fences between 1933 and 1942. Much simpler and cheaper, four-strand barbed wire fences were used—rather than the elaborate Wallowa design—as hunting, trapping, and poisoning reduced predator populations region-wide. If both fencing and predator control were already planned or ongoing, and might well have proceeded without “scientific” support, then the greater significance of the Coyote-Proof Pasture Experiment lies in the ways it altered the institutional context itself, and thereby the trajectory of rangeland administration and research.

The COGS was calculated by dividing the annual operating cost by the annual production

The total cost of the all the UAEs, depicted by the Batch Generic Box in SuperPro, was the most expensive , followed by the belt press filters . For better cost estimates of the UAEs used in the model, a quote was obtained from an industrial manufacturer named REUS® . Pricing forthis equipment was reduced by about 5%, assuming a conversative approach on the expected discount for bulk purchase of equipment .The OPEX, otherwise referred to as annual operating cost, and COGS are shown with and without depreciation, insurance, and local taxes. The annual operating cost and COGS without including depreciation are more representative of the expected cost for production since depreciation is typically spread over years to effort to expense cost over time while simultaneously lowering the value of the asset . Also depreciation is not a cash outlay so does not have a negative impact on profitability . Here, depreciation was calculated using the straight-line method with a depreciation period of 10 years and salvage value of 5% of the DFC. Insurance was estimated to be 1% of the DFC and local taxes were estimated to be 2% of the DFC. As expected, evaluation of the facility’s annual operating costs and COGS including depreciation, insurance and local taxes are higher compared to without.An analysis was performed to investigate the effect of factors such as depreciation, insurance, factory expenses and local taxes on the economics of the model. A breakdown of the annual operating for each case are shown in Figure 3.4 through 3.6. In all cases, nursery grow bag the annual operating costs are composed of the following: raw materials, labor dependence, facility dependence, Laboratory , Quality Assurance , Quality Control , consumables, waste treatment and utilities.

When the first of the three cases were analyzed, it was determined that the largest contributor to the annual operating costs was the facility dependent cost. Here, facility dependent costs include maintenance, depreciation, insurance, taxes, and factory expenses. Maintenance of each equipment was determined using equipment-specific multipliers, default values provided by SuperPro. No pre-existing depreciation of equipment was assumed in the model. Local taxes were assumed to be 2% of the DFC, in alignment with values from municipal tax charts for South Dakota . Percentages for insurance and factory expenses were estimated using a previous SuperPro model performed internally which collaborated with a plant based bio-manufacturing facility for pricing. The second largest contributor to the annual operating cost was the raw materials. A breakdown of the raw materials and their costs can be seen in Figure 3.7. Here, raw material costs amounted to a value of $5.7 million per year. A list of the raw materials used, the quantity used annually, and their total cost can be seen in Table 3.7. Ethanol, or ethyl alcohol as defined in SuperPro, had the biggest yearly expense out of all the raw materials, totaling over $4.1 million. The model estimates 6.1 million kg of ethanol needed to produce 100 MT of resveratrol. It was determined that the cost of ethanol is over two thirds of the total raw material cost, accounting for 71% of the raw material cost. The price used in the model was $2.00 per gallon of ethanol. This value was retrieved using data by Market Insider which tracks the price of commodities like ethanol daily . For the second case, which only excludes depreciation as a factor contributing to facility dependent costs, the largest contributor was raw materials. Raw materials costs are the same in all three cases but now account for 50.9% of the annual operating cost, with the cost of ethanol remaining the largest contributor.

The second largest raw material cost is the knotweed rhizomes used in the process. Roughly 7.3 million kg of knotweed rhizomes are needed in the process simulation to produce 100 MT of resveratrol. The cost for producing a kg of knotweed rhizomes was estimated to be about $0.19, totaling $1.4 million or 24.3% of the raw materials cost. Estimates for a kg of knotweed were calculated to include costs associated with pre-planting, harvesting, postharvesting, farm equipment operating costs, and cash overhead . A breakdown of these calculations is shown in Chapter 2. Following the second case, which provides a clearer estimate of annual operating cost because depreciation is not included, labor was the third highest contributor to the annual operating cost.Here, labor dependent costs refer to the costs associated with labor hours required to effectively operate the downstream processing section. Labor cost for the upstream production of Japanese knotweed rhizome was already included in the overall unit cost for a kg of knotweed. In this model, the downstream operators earn a wage of $22 per hour. The downstream labor cost includes 40% benefits factor, 10% operating supplies factor, 20% supervision factor, and 60% administration factor. An example of this distribution is as follows: for every $20 paid to an operator for an hour of work, there is an additional cost of $8 for benefits, $2 for supplies, $4 for supervision and $12 for administration. Total labor hours amount to 20,398 hours per year with operators devoting most of their time to the blending tank used as an adsorption vessel. This was in alignment with our estimates since the blending tank is one of the two equipment with the highest number of operations needed compared to every other unit in the model. The next factor contributing to the annual operating cost was the Lab, QA, and QC group, which accounted for less than 1 percent. No funds were allocated to on-going research and development whereas, QA and QC were estimated to be 5% of the total labor cost.

Only one consumable was defined within the process, as listed in Table 3.8. This consumable was the macroporous resin, NKA-11, which is used within the batch adsorption vessel . The annual cost of the resin amounted to $130,536. Pricing for this resin was estimated using commercial values found online by a large scale supplier Co., Ltd.. Waste treatment costs were also incorporated into the annual operating cost calculations. This waste treatment category incorporates the price to safely dispose of different waste streams generated in the facility. Exiting waste streams are classified under one of four groups: organic waste, aqueous waste, solid waste, or gaseous emissions. Here, the cost of emissions was negligible due to low concentration of nitrogen, oxygen, ethanol, and water vapor being released into the atmosphere. The cost to dispose of each group is as follows: organic waste is $0.01/kg, aqueous waste is $0.001/kg and solid waste is $0.01/kg. The annual amount of waste produced by the process is about 44 million kg, estimated at about $234,000 per year in waste treatment costs. While the cost for waste treatment of each group remained similar, the largest contributor to the cost is attributed to organic waste, which include the disposal of the biomass. Resin disposal costs were negligible. These results are in alignment with the cost and use of raw materials throughout the process. A breakdown of the annual waste treatment contributors is shown in Figure 3.8.Utilities are the last factor which contribute to the annual operating cost. A detailed breakdown by utility type used in the process model annually can be seen in Figure 3.9. Three utility types were used in the model: standard power, steam, and chilled water. Pricing for each utility were set as, $0.10 per kW-h of power, $12.00 per MT of steam, and $0.40 per MT of chilled water. The cost of std power was the highest of all three utilities. Std power was used to operate all major equipment including but not limited to the stirred reactor acting as the enzymatichydrolysis unit, plastic growing bag the blending vessel serving as the adsorption vessel, the ultrasonic assisted extractor unit, and the grinder. Surprisingly, the grinder responsible for crushing and homogenizing the plant material into powder required the largest amount of power to operate. Total utilities costs only amounted to a 1% of the annual operating cost. In this model, the annual operating cost for each case described above were used and divided by 100 MT, resulting in a COGS of $150/kg, $111/kg, and $79.4/kg. A summary of both the process parameters and economic analyses from the base case model is shown in Table 3.9.Throughout the design of the simulation model, certain parameters were initialized to match resveratrol production practices described in both patents and scientific literature. When extrapolating bioprocessing parameters from these sources, a conservative approach was taken so that the lower values from a range of data was used in the model in effort to portray realistic results expected during large scale production. An example of this is including a loss of 5-10% material during a filtration step, instead of expecting an 100% recovery every batch. As acknowledged in Chapter 3, certain bioprocessing parameters used in the model were assumptions, thus leading to some uncertainty whether the process outputs were reliable. For this reason, a sensitivity analysis was performed to assess how a certain variation effects the economics of the model, specifically the CAPEX, OPEX, and COGS. Similarly, it is acknowledged that there are multiple designs which can be implemented to produce 100 MT of resveratrol and this model simply serves as an example of one method. Therefore, certain scenario analysis were performed to assess the effect of design and production amount on key economic parameters.To assess the sensitivity of the base case model, we varied certain process parameters to investigate their impact on the CAPEX, OPEX, and ultimately the COGS. One parameter which was defined in the SuperPro model using a conservative approach was the concentration of resveratrol present in the Japanese knotweed rhizome used for processing . Figure 4.1 demonstrates the relationship between COGS and when the concentration of resveratrol per knotweed rhizome is increased up to 3 mg/g FW. The concentration of polydatinwas also increased proportionally along with resveratrol. Each facility simulation is redesigned for each concentration tested while still reaching 100 MT resveratrol annually.As expected, when a larger concentration of resveratrol was modeled to be present within the Japanese knotweed rhizome entering the process, the model resulted in a decrease in both the CAPEX and COGS. The largest drop in CAPEX occurs directly when the concentration of resveratrol is increased an order of magnitude to 1.0 mg/g. The decrease in CAPEX from 0.5 mg/g to 1.0 mg/g is $13 million, over 3-fold larger than the average drop between increments. Expectedly, the COGS also decreases the largest amount between the first two concentrations. The COGS drops a value of $48/kg to a value of $102/kg, a 32% drop. Values for both economic parameters begin to plateau around a concentration of at 1.5 mg resveratrol/g FW, approximately at values of $27.8 million and $88/kg for CAPEX and COGS , respectively. It was discovered that resveratrol is found in a wide range of concentrations in Japanese knotweed rhizomes. A table listing different resveratrol concentrations found in Japanese knotweed is shown in Chapter 2. As mentioned, an average value of the total resveratrol concentration in knotweed, including polydatin, was about 7 folds higher than just free resveratrol. Using information on resveratrol concentrations for Japanese knotweed rhizomes grown specifically in North America1 , we calculated an average concentration value of 2.6 mg resveratrol/g FW. If a future process used rhizomes under similar conditions , our simulation suggests a cost decrease of about one third for CAPEX and 43% for COGS compared to our base case model operating at 0.5 mg Rsv/g. Notably, the same authors that describe an average concentration of 2.6 mg Rsv/g Japanese knotweed roots also mention certain samples contain concentrations as high as 12 mg Rsv/g FW and 12 mg Polydatin/g FW. Using this information, a simulation operating at the same conditions as the base cased was modeled using concentrations of resveratrol and polydatin at a 1:1 ratio at 12mg/g for each stilbene compounds. The same economic analysis was performed on the 12mg/g case, the resulting CAPEX, OPEX, and COGS values are shown in Table 4.1.In the previous chapter, the cost of ethanol was identified to be the major contributor to the annual operating cost and the largest bottleneck.

Information on plant host taxonomy was gathered on NCBI’s Taxonomy Browser

The marginal ancestral state likelihood estimates of each host for all internal nodes of the ML phylogenetic trees were calculated using the re-rooting method of Yang et al. in the R package phytools , and mapped using the package APE . This method uses the phylogeny of extant taxa to reconstruct ancestral traits of extinct ancestors by analyzing phylogenetic parameters , along with a model of nucleotide substitution, to build posterior probabilities of character states at each interior node by randomly re-rooting the tree at each internal node and calculating the probability of observing the extant distribution of traits over all possibilities of that internal node character identity. The ML estimates at each internal node were calculated based on both the equal rates transition model and the symmetrical rates transition model . The fit of the two models to the data was compared using the Akaike information criterion and can be seen in supplemental table . Scoary was used to test if the pan-genome was correlated with hosts at either the super-order scales or the genus scale by conducting a Fisher’s Exact Test . FET measures the association of each gene in the pangenome to a trait of interest, plastic planters which in this case is plant host. While FET requires no association between datapoints, Scoary uses a phylogeny in order to remove lineage specific interdependencies and corrects the p-value based on those interdependencies.

Significance was evaluated by the “worst pairwise comparison P”, for the phylogenetic corrections, not the naïvep-values from FET. Individual analyses were conducted to test for correlation of gene presence and absence with each of the 29 coded host groups .California is the center of American fruit, vegetable and nut production and is a globally important exporter of plants . As a hub of international plant trade, California has been both a source and a sink for the introduction of novel phytopathogens across the globe. Xylella fastidiosa ssp. fastidiosa is a bacterial pathogen that impacts grapevines and almonds in California. X. fastidiosa subspecies fastidiosa, was likely of Central American origin and possibly introduced via coffee imports to the United States in the 1800s . After its introduction into the U.S., the pathogen has continued to spread, including recent introductions into Spain and Taiwan . In California, outbreaks of X. fastidiosa ssp. fastidiosa have been devastating to the grape industry. In contrast, infections have not caused notable issues in the agriculturally important regions of Central America, where the ssp. has likely existed in the ecosystem for at least thousands of years . This is not a universal trait of X. fastidiosa in C. arabica as opposed to Vitis. In Brazil, X. fastidiosa ssp. pauca infections reduce C. arabica yields, despite the pathogen’s long history in the region . X. fastidiosa is complex in its strain specificity to various host plants, with many conflicting examples of host susceptibility in geographic areas. Nonetheless, documenting host susceptibility to different subspecies and strains of the pathogen is vital. Clarifying the process of host jumps in X. fastidiosa is urgent for global agricultural security, as new outbreaks continue to put plant health at risk. While it is likely that X. fastidiosa ssp. fastidiosa is broadly present in Central America, currently only Costa Rica has reported this pathogen.

Here we refer to “Costa Rican” strains as our representative Central American strains, although there is certainly much larger diversity of X. fastidiosa in the region. X. fastidiosa ssp. fastidiosa frequently infects C. arabica in Costa Rica without substantially impacting production, although subtle leaf curling symptoms correlate with X. fastidiosa infections . In Costa Rica, X. fastidiosa ssp. fastidiosa also infects periwinkle , guava , and avocado . While symptoms have been reported in these hosts, disease severity and progression have not been thoroughly evaluated. Some strains of X. fastidiosa ssp. fastidiosa have been detected in Vitis in Costa Rica, where V. vinifera is not commonly grown but only sequence-type data was collected for these strains. However, a recent study found that related strains from Costa Rica did not cause infections in Vitis . Given these results, infectivity and virulence by the same strain may differ between Vitis and C. arabica. A clear phylogenetic delineation of pathogenicity, which requires extensive crossinoculation experiments, has not yet been determined. This study begins to test the adaptation of the U.S.-introduced clade to Vitis cross-inoculation and computational methods. Within the genus Vitis, there is wide genetic diversity conferring variation in both tolerance and resistance to Pierce’s disease , the disease caused by X. fastidiosa ssp. fastidiosa infections in grapevines. Differences are seen between species such as the naturally tolerant V. arizonica and Muscadinia rotundifolia, and the susceptible V. vinifera . This variation has recently been used to hypothesize that there might have been PD present in the U.S. for longer than has been estimated by all evolutionary data . But disease susceptibility also varies largely within V. vinifera , an intensely bred species that includes a wide and diverse range of cultivars. In contrast to crespera in Costa Rica, PD in the United States is quite virulent, persistent, and economically damaging. Infection with the pathogen causes severe symptoms such as stunting, leaf scorch, shriveled fruit, and eventual plant death in many cases . There are several additional adversely impacted hosts of species within the introduced clade. Of economic importance is almond suffering Almond Leaf Scorch . Other hosts of strains in the introduced clade include maple , plum , ragweed , sweet cherry , and western redbud . C. arabica is a unique host for X. fastidiosa. It has widespread documented infectivity coupled with non-deadly symptoms, and it was likely the source of two economically devastating introduction events: to the United States in grape and to Italy in olive . Nonetheless, C. arabica is generally understudied as a host plant. Data from Brazil show that C. arabica and Citrus strains of X. fastidiosa ssp. pauca do not cross-infect , showing variation in host specificity. However, C. arabica is also a host of at least two of the three major subspecies of X. fastidiosa, showing broad susceptibility. It has been hypothesized that strains of X. fastidiosa do not need to undergo many genetic changes in order to infect C. arabica, potentially making it a host susceptible to low-cost host jumps . While ssp. pauca and ssp. fastidiosa have both been documented to infect C. arabica, slow progression and low virulence have led C. arabica to be described as a “latent carrier” . So far, research has shown that the strains able to infect C. arabica are quite diverse . The infectivity of X. fastidiosa ssp. multiplex has never been tested in C. arabica. Subspecies multiplex is native to the south-eastern United States and is present throughout the U.S., South America, and Europe . However, the advent of coffee production in California coupled with historic documentation of ssp. multiplex in Argentina, Brazil, and Paraguay, call for investigation of whether C. arabica can serve as a host to ssp. multiplex . Upon introduction to the U.S., X. fastidiosa was exposed to new hosts, climate, vectors, and agricultural practices alongside a rapidly growing agricultural industry. While it is likely that the introduced strains were exposed to selective pressures that increased specificity to the conditions in the U.S., mainly climate and host, plastic plant nursery pot that is not the only possibility. It is also plausible that instead of adapting to the specific pressures in the new environment, the introduced population of strains became generally hypervirulent.

This could explain the high level of virulence towards V. vinifera in the U.S. compared to the lower virulence observed in Costa Rica. We investigated both neutral and adaptive changes during the process of naturalization for this pathogen using whole genome sequences.In March 2022, 120 coffee plants were divided into equal groups and inoculated with either a buffer control or one of four CA strains of X. fastidiosa, ssp. fastidiosa: Napa1 , ALS17T5 , Je115 or subspecies multiplex ALS15T2 , or the sterile succinate-citrate-phosphate buffer control . While Napa1 and Je115 are both isolated from V. vinifera, they are from different climates and clades of X. fastidiosa, and might have experienced distinctive selective regimes . ALS17T5 was isolated from a symptomatic almond plant and is in one of several clades of strains infecting only almond trees, however this clade is nested within other clades that infect Vitis . Disease-free C. arabica plants were donated from Frinj coffee company for this experiment. As positive controls, 20 V. vinifera cv. Chardonnay were inoculated as well as 50 Helianthus annus. V. vinifera is known to be highly susceptible to subspecies fastidiosa but not subspecies multiplex, while Helianthus annus is susceptible to both subspecies . Cell suspensions were prepared by suspending week-old cells grown on solid medium in SCP buffer just prior to inoculation. Each suspension was made by scraping 10 streaks of 20µL into 1 mL of buffer. Inoculations were conducted using 2 10µL beads of inoculum and a 00-size entomological pin used to pin prick through the bead of cell suspension several times until the inoculum absorbed into the plant xylem. Inoculations were conducted on small plants with typically 3 full-leaf pairs, and two inoculum beads were placed just above and below the center leaf pair. Plants were not waThered the morning of inoculation to optimize absorption of the inoculum into the xylem vessels. Symptom measurements took place in August, September, October, December, 2022, and January and April 2023, following the March 2022 inoculation. In August and September 2022, all internode lengths along the main stem of the plant were measured along with the heights of each plant. In October, December, January, and April, only the total plant heights were measured to detect stunting. Plants were also visually assessed for foliar scorching.Control and experimental plants were tested for the presence of X. fastidiosa via qPCR or culturing. DNA was extracted using a DNeasy plant mini kit and then quantified using qPCR with a primer pair targeting the gene encoding RecF, RecF1_F+R; for qPCR protocol see Sicard et al. . qPCR was run in duplicate with positive and negative controls on each plate; all samples with Ct values of 37 or higher were considered “undetected” and were considered negative. Culturing was conducted using approximately 0.1 grams of petiole tissue, or the entire petiole and midrib of a C. arabica leaf. Samples were surface sterilized, chopped, ground in a Polytron, and plated on PWG media using the method from Hill and Purcell . On April 3 rd , 2022, petioles that were directly above the inoculations site were sampled from all 50 sunflower plants and X. fastidiosa populations were measured via both qPCR and cell culturing. On June 16th , 2022, V. vinifera plants were sampled from and also visually assessed for symptoms. On July 14th , one leaf from the 2 nd leaf pair above the inoculation point was used for detection via qPCR. On August 1 st, ,15th , and 18th , one leaf from the 2 nd leaf pair above the inoculation point was cultured from each C. arabica plant. On October 4 th , 2022, the opposite leaf was taken for qPCR. While Vitis is not as susceptible infections from ssp. multiplex, typically there is some detectable infection, however those infections have 10-100 fold lower population sizes than infections by ssp. fastidiosa . In February – April 2023, we cultured from the leaves using a different method due to loss of leaves from the ~ 8 cm above the inoculation point. The lowest 8 leaves were collected, and the petioles were pooled for sterilization, tissue grinding, and plating.All analyses were performed using R statistical software version 4.2.3 . Linear mixed models were built using the package lme4 to add a random effect to account for repeated sampling of the plants over the course of the experiment . As analyses of virulence, we tested the effects of variety and the interaction between treatment and sampling date on both internode length, and separately on plant height , each with a random effect of plant ID using the lmer function of lme4. As an analysis of infection persistence via detection assays, we tested the effects of variety and the interaction between treatment and sampling date on the count of plants that tested positive using a generalized linear mixed model with binomial error using the package glmmTMB and the function glmmTMB .

Trichome density changes genetically and environmentally

When S. ×utahensis plants dehydrated as a result of decreasing substrate volumetric water contents, plants closed their stomata to reduce transpiration and stomatal conductance as a drought acclimation to maintain plant water status and prevent water losses and further dehydration . Although CO2 uptake is limited when stomata are closed and stomatal conductance reduced , plants had a lower proportion of visibly wilted leaves in this study or better aesthetic quality under drought conditions in other reports. These plants were considered drought tolerant in ornamental plant evaluations in semiarid regions in Australia and the United States . Shepherdia ×utahensis reduced its midday stomatal conductance at lower water availability and can be considered as a low water-use landscape plant. Plants have the capacity of regulating stomatal conductance that is related to their habitat aridity. Kjelgren et al. reported that plants native to arid regions, such as Dianella revoluta ‘Breeze’ and Ptilotus nobilis , showed greater reduction of stomatal conductance compared with those from humid areas. Because of restricted transpiration, plants with acclimation capability may reduce leaf size to enhance convective heat loss to mitigate heat stress that causes high leaf-to-air VPD and leaf wilting . The fact that leaf-to-air VPD increased when substrate volumetric water content decreased is likely a direct consequence of increased leaf temperature because leaf vapor pressure is estimated by leaf temperature.

To avoid heat stress, french flower bucket leaf energyis balanced primarily using sensible heat loss under drought . The efficacy of sensible heat loss relates to boundary layer resistance, which is positively correlated to leaf width . Under drought conditions, cell division and leaf expansion are limited , and smaller leaves are beneficial for dissipating heat through convection and conduction to maintain leaf temperature close to air temperature . In this study, S. ×utahensis produced smaller leaves under water stress and leaf size of plants grown at the substrate volumetric water content of 0.05m3 ·m−3 was 51% smaller than those at 0.40m3 ·m−3 . This result is in line with previous studies that consistently reported reductions in leaf size under water stress for drought-tolerant ornamental plants . For instance, Zollinger et al. suggested that small leaves allow Lavandula angustifolia and Penstemon ×mexicali ‘Red Rocks’ to reduce water loss when irrigation intervals were increased from 1week to 4weeks. Toscano et al. also found that leaf size of Viburnum tinus ‘Lucidum’ decreased by 19% to acclimate to drought stress. Shepherdia ×utahensis decreased total leaf area under water stress as a result of reductions in leaf number and size . However, plants with decreased total leaf area have fewer stomata and less light interception, which controls transpiration and leaf temperature, respectively . Reduced total leaf area has been reported as a means of avoiding drought stress in ornamental plants such as Lavandula angustifolia, Pittosporum tobira , and Viburnum tinus ‘Lucidum’ . The root growth of S. ×utahensis was enhanced at low substrate volumetric water content, while shoot growth was inhibited, resulting in a higher root-to-shoot ratio , which helps plants to obtain water more efficiently. Rosa hybrida ‘Ferdy’ and Populus cathayana have been observed to increase root growth to maintain water status under water stress .

Drought-tolerant plants native to the western United States also produce small leaves and deep roots to reduce water demand and loss and increase water uptake . In this study, as substrate volumetric water content decreased, leaves of S. ×utahensis curled as stem water potential became more negative. At the substrate volumetric water content of 0.05m3 ·m−3 , the leaf curling index was 0.17, suggesting that the light interception area was 83% that of flattened leaves. Similarly, Dianella revoluta ‘Breeze’ and Ctenanthe setosa have been shown to minimize sunlight exposure through leaf curling under water deficit . Although light-harvesting efficiency is reduced, leaf curling limits water loss from transpiration and protects plants from overheating to sustain photosystem functions and other biochemical/physiological processes . In addition, as the rooting substrate became drier in this study, specific leaf area decreased, indicating that leaves became thicker , which prevented leaves from overheating. Plants may decrease specific leaf area to acclimate to water stress as reported in Ptilotus nobilis . The trichome density of S. ×utahensis in this study was affected by substrate water availability and plant water status . Water-stressed S. ×utahensis produced densely packed trichomes, resulting in a silvery appearance, while well-watered plants had fewer trichomes to cover epidermal cells and exhibited a greener color . Trichomes promote leaf reflectance , which helps balance energy and reduce heat stress . Positive effects of trichomes on leaf reflectance of visible light have been reported on Verbascum thapsus and Salix commutata . However, because trichomes are broad-spectrum reflectors , the reflectance of PAR, blue, green, and red light is proportional to the trichome density . When substrate volumetric water content decreased, the reflectance of green light did not increase as much as blue light and red light due to the chlorophyll in the epidermal cells . Increased leaf reflectance has been shown to sacrifice the efficacy of light-harvesting pigments and reduce the net assimilation rate when plants are grown in drier conditions. Previous research also suggested that trichomes improved the reflectance of near infrared light . However, in this study, denser trichomes produced in drier substrate did not affect near infrared light reflectance of S. ×utahensis . Slaton et al. reported similar results that near-infrared light reflectance was not affected by increased trichome density in 48 species. More studies are needed to evaluate the effects of trichomes on near-infrared reflectance. Increased trichome density has smaller effects on decreasing gas exchange coed with the effects on leaf reflectance . However, densely packed trichomes covering the stomata of S. ×utahensis may increase resistance to transpiration and reduce water loss . Leaf trichomes also increase leaf roughness and increase the laminar boundary layer to restrict air movement across leaf surfaces to reduce transpiration . Eriogonum corymbosum and S. rotundifolia produce leaf trichomes for better protection from wind and to maintain water status . Densely packed trichomes add an atmospheric boundary layer that imposes additional resistance to water vapor diffusion . However, CO2 influx is also limited by the boundary layer resistance, decreasing the net assimilation rate . Although trichome-induced boundary layer resistance has a smaller effect on transpiration than stomatal conductance , it still provides an advantage for desert plants to survive in dry and hot conditions. The genetic regulation of trichome density of Caragana korshinskii has been reported by Ning et al. . However, it is unclear how xeric plants change their trichome density to acclimate to drought conditions. A negative correlation between leaf trichome density and leaf size or epidermal cell size occurred in this study , which suggests that cell expansion maycontrol trichome density. Low trichome coverage fraction, which was related to greater space between trichomes, showed when epidermal cell density decreased, indicating cell expansion may coordinate trichome density. Ascensão and Pais reported the number of trichomes is determined during leaf lifespan, and leaf cell differentiation does not affect trichome number. Similar results showed in our research that plants had similar total numbers of trichomes per leaf at different substrate volumetric water contents. This may indicate that S. ×utahensis develops trichomes independent of leaf development. In fact, trichomes develop at the early stage of leaf development and often earlier than stomatal development . For instance, trichomes of Inula viscosa are fully developed and reach mature size when leaves are 2mm long; however, a mature leaf is 6–8cm long .

Ocimum basilicum forms trichomes at an early stage of leaf development and trichomes then grow independently . In the same study, trichomes covered young leaves but became more widely spaced when leaf cells started to expand . In our study, the total number of epidermal cells per leaf was similar on plants at different substrate volumetric water contents, bucket flower which indicates cell differentiation might have minor effects on regulating trichome density. In contrast, cell expansion might be the main factor for regulating trichome density because leaf size, epidermal cell size, and the space among trichomes changed along with substrate volumetric water contents and correlated significantly with trichome density of S. ×utahensis . Ehleringer found a negative correlation between leaf size and trichome density of Encelia farinosa, but cell size was not determined. Cell enlargement at high soil moisture levels amplified leaf size and the space among trichomes, reducing the trichome density on the S. ×utahensis leaves in this study. The relationships between trichome density, epidermal cell size and density, and leaf reflectance might indicate changes in cell size predominantly controls trichome density to modify leaf reflectance. Modifying leaf reflectance via the change in cell size helps rapidly acclimate to environmental change without compromising whole leaf function . Cell-expansion driven leaf anatomic change has been widely reported on adjusting stomatal density . For instance, Murphy et al. observed that cell expansion was the predominant factor for coordinating vein and stomata density of eight angiosperm species under sun and shade. Stomatal density decreases and the size of guard cells increases when leaf water potential increases , suggesting cell expansion not only enlarges the distance between epidermal appendages but also increases their size. Environmental factors also promote leaf trichome density. Such factors include increased leaf-to-air VPD and drought , all of which negatively affect plant water status. For instance, high leaf-to-air VPD may increase water loss via transpiration, leading to plant dehydration. Leaf trichome density of Cucumis sativus increased when air humidity decreased from 90 to 20% at 28°C, causing leaf-to-air VPD to increase from 0.4 to 3.0kPa . Shibuya et al. did not investigate cell or leaf expansion of C. sativus, but increased leaf-to-air VPD may promote trichome density because rising leaf-to-air VPDs reduces cell size , making space between trichomes smaller. In this study, higher leaf-to-air VPD and smaller leaves were observed when S. ×utahensis plants grew at the lower θtand the smaller epidermal cell size resulted in greater trichome density. Therefore, because increased leaf-to-air VPD and drought led to a reduction in cell enlargement and denser trichomes in S. ×utahensis, leaf trichome density was regulated using turgor-pressure-driven cell expansion to acclimate to drought conditions.Modern plant trade disturbs historical ecological relationships and creates opportunities for the development of novel pathogenic interactions , often with correlated genetic changes . However, pathogens must be adapted to the environment of the novel host before they meet, or they will not be able to survive and reproduce . That does not mean pathogens necessarily pre-adapted to the exact same host, but either could have adapted to a similar host earlier and retained that adaptation until encountering a novel host. Convergent evolution in diverse pathogen populations can allow for divergent strains to have the ability to infect the same hosts. Three potential mechanisms of genetic change that can accompany host shifts are nucleotide changes leading to different alleles in the core genome of a pathogen , whole gene gain and loss in the pan-genome, leading to unique sets of genes in individual strains, or regulatory/epigenetic changes. Due to the recent increase in whole genome sequencing of plant pathogens, we can now more effectively use phylogenetic analyses to investigate their genetic associations to both novel and historical host plants . Understanding the phylogenetic relationships between specific host and pathogens should allow the development of preemptive plans to protect natural ecosystems as well as agriculture from the emergence of novel pathogens. Xylella fastidiosa is an insect-transmitted, xylem-limited bacterial plant pathogen found across the Americas, and as of recently, globally. X. fastidiosa is considered to be a generalist pathogen, because, as a species; it reportedly infects at least 563 species belonging to 82 botanical families . The lack of host specificity that X. fastidiosa exhibits as a species contrasts with increased plant host specificity in smaller clades and strains . It is still debated whether a pathogen like X. fastidiosa should be considered a generalist species that “leaps” between phylogenetically distant hosts or, alternatively, a crawler at shallower clades . The difference is biological as there are unique implications for either evolutionary path. X. fastidiosa could be repeatedly evolving specialization or it could have biological and genetic traits as a species that make particular hosts of disparate plant taxa suitable. From an applied perspective, there have been recent calls from government agencies for increased focus on understanding the host range of X. fastidiosa. This is because the pathogen has been deemed likely to spread and to be of extremely high risk to crops of agricultural value , 2015. Xylella fastidiosa causes disease in a range of high value crops, including Pierce’s disease of grapevines, citrus variegated chlorosis disease in sweet oranges, almond leaf scorch, leaf scorch of coffee, olive quick decline syndrome , spanning North and South America, Europe, the Middle East, and Taiwan .

Soil data for the Bay Area are obtained from the SSURGO datasets

Dry soils also inhibit reproduction, survival, dispersal, and development of symptoms in host plants . Soil moisture conditions experienced by the pathogen themselves arise from interactions among the precipitation regime, soil depth, drainage, and atmospheric evaporative demand, and thus reflect the interplay of edaphic and climatic conditions. Finally, Pc disease is also often suppressed in rich soils where organic matter content exceeds  5% , probably because of predation by other soil organisms in the diverse microfaunal communities sustained in these soils . Projections of potential future risk therefore require techniques to assess the impact of multiple environmental changes and their interactions on pathogen range and epidemiology. Here, we apply a mechanistic modeling approach to explore how climate change could impact pathogen range and activity. We explore how simultaneous changes in temperature, precipitation, snow-pack extent, and evaporative demand might impact the range of a well-characterized pathogen under different climate scenarios. To do this we couple two existing models that describe controls on the range of the generalist root pathogen Pc in the state of California and surrounding regions in the states of Oregon, Nevada, and Arizona in the southwest USA. Pc occurs in this region but its range is poorly delineated. In other warm climates such as southern Australia and Hawaii, Pc has had a devastating effect on timber production, natural forests and agriculture € .

Modeling the climatic and edaphic limits on its potential range in the US southwest will help determine the risks posed by this pathogen, procona florida container particularly since there is not as yet a detailed understanding of the susceptibility of native species to Pc infection.To address the research questions, we link two existing models that describe: Pc winter survival probabilities and Pc disease severity during the spring growing season . We first model winter soil temperatures. Second, we use the modeled soil temperatures to estimate Pc survival with an existing, validated survival model . Next, we expand a probabilistic soil moisture model to account for snowpack contributions. We use the projected soil moisture fluctuations to drive a stochastic pathogen risk model, previously employed to describe disease range, relative risk on different soil types, and response to different irrigation regimes at sites in Western Australia and Oregon . The output of this model is a metric of the likelihood of disease expansion throughout a hosts’ root system, lying between 0 and 1 . Fourth, we identify regions with high soil organic content from the national State Soil Geographic Database and Soil Survey Geographic Database data sets and exclude these regions from the potential Pc range . Finally, we compute a relative Pc risk as the product of the risk of winter survival, absence of suppressive soil conditions, and moisture-controlled host colonization. We apply this modeling framework to baseline and future climate scenarios across the southwest USA. We repeat the analyses using high resolution downscaled climate and soil data for the San Francisco Bay Area , allowing us to explore the effects of microtopography, orography, and fine-grained changes in soil properties, typical of the California coastline, on the spatial patterns of Pc disease risk and its projected climate sensitivity.We use three different climate datasets: monthly regional historical climate observations interpolated to a 140-km2 grid for the 1950–2000 period ; bias-corrected and spatially downscaled National Center for Atmospheric Research Community Climate System Model 3.0 monthly simulations run for the IPCC A2 and B1 climate scenarios on the same 140 9 140 km grid ; and 30-year average monthly projections from the Geophysical Fluid Dynamics Laboratory climate models for the same scenarios, downscaled to a high resolution grid over the San Francisco Bay area for the 1970–2000 and 2035–2065 periods.

To obtain representative climate characteristics, we average the monthly regional data over 10 year intervals . The climate components obtained for each dataset are summarized in Table 1. We estimate average daily temperature for the Bay Area climate surfaces as a simple mean of monthly maximum and minimum temperature; and estimate the potential evaporation at the southwestern US scale using the Priestly–Taylor equation. One of the limitations of this approach lies in the considerable uncertainty surrounding the use of downscaled climate datasets. Our goal in using these high resolution data sets was to: develop a reasonable baseline drawn from interpolated observations and compare this against plausible climate futures to elucidate the range of potential Pc responses and how they could manifest themselves at regional scales. We do not interpret these future cases as predictions or forecasts, but illustrative scenarios to offer insight into process interactions. Regional soil data, specifically percentage clay, organic matter content and soil depth data are derived from STATSGO . Both data sets derive from soil maps with irregular mapping units ranging from <1 ha in size to >20 000 ha . The soil data were thus mapped to the climate grid scales used for both the Bay Area and the regional scale. We use tabulated hydraulic parameters to estimate the properties of the water retention curve for each soil type on the basis of the percentage clay content, classifying <10% clay soils as sands, 10–15% clay content as loamy sands, 15–20% clay contents as sandy loams, 20–40% clay content as loams and >40% clay contents as clays. Georeferenced locations where Pc has been isolated from natural or agricultural soils in California were obtained from the Phytophthora Database , Forest Phytophthoras of the World , and recent literature sources .The southwestern USA contains a climatically, geologically, and biologically diverse set of landscapes . Climatically, the region encompasses cool, humid, oceanic climates near the coast, arid interior valleys and deserts with continental climates, and montane and nivial uplands in the Sierra Nevada mountain ranges. Precipitation tends to be highly seasonal with the majority falling during a 4-month period from December to March. The diversity of climatic types is reflected in the diversity of terrestrial vegetation, with important vegetation communities including coastal rainforests, dry chapparal, Mediterranean oak savannas, montane forests, and subalpine regions . Within the study area, the California Floristic Province is recognized as a biodiversity hotspot , and the state is home to over 7600 plant species . The San Francisco Bay Area, which is used as a higher resolution case study maintains much of the same complexity of climate, geology, and biology, due to strong ocean-inland gradients in temperature and annual precipitation , and the activity of various tectonic faults in the region . The diversity of climate and soil types in the study area can be directly parameterized within the stochastic soil moisture and pathogen models. Note that we have allowed spatial variations in evaporative demand to be driven by radiation, temperature, and humidity , and have not incorporated an explicit treatment of varying vegetation cover on Emax. The rationale for this lies in the fact that surface conductance loses sensitivity to LAI for low humidity and LAI > 2, suggesting that the first order controls on Emax variation lie in energy limitation rather than vegetation characteristics . The diversity of the vegetation also presents a challenge for parameterization of the pathogen dynamics in the model. It is clearly not feasible to parameterize specific host resistance terms for Pc for all 7600 vascular plant species in the study area; furthermore assessments of susceptibility of the native vegetation to Pc are strikingly incomplete, and remain insufficiently quantitative to allow a direct comparison to be made. The problem is similar to applying this model to Pc in Western Australia, where far more work on Pc susceptibility has been performed. In that comparably biodiverse landscape, 3084 of the 5710 native plant species were found to be susceptible or highly susceptible to Pc infection . This determination, however, required intensive surveys of native vegetation in representative habitats, as well as comparative inoculation greenhouse trials , and still provides only qualitative information about the relative vulnerability of host plants. No comparable investigations have as yet taken place in the US south west, procona London container although numerous common native species are known to be vulnerable to Pc, including manzanita species , chapparal species , MonTherey pine, sycamore, western sword fern, coast live oak , bay laurel and madrone .

In the absence of detailed observations of host–pathogen interactions for the diversity of species in Western Australia, we explored climate limitations on range using a combination of: laboratory estimates of Pc growth rates in sterile media and host resistance parameters taken from a moderately susceptible host and estimated to be 0.1 day 1 . We have repeated this approach for the US southwest with the following rationalizations: moderately susceptible host dynamics ensure that climate effects are clear. For highly susceptible hosts, climatic drivers are relatively unimportant, since any serendipitous infection is likely to lead to mortality. For unsusceptible hosts, climate is equally unimportant. Thus to explore the impacts of climate change, moderately susceptible host assumptions are sensible; Evidence from infections of common tree species suggests that the hosts are indeed moderately susceptible to Pc, with infections most common in wet areas and mortality commonly associated with concurrent pest or drought stress ; for the sake of identifying climate sensitivity specifically, holding host vulnerability constant across the study area provides a useful control, and in the absence of more detailed information, is the most parsimonious modeling approach. As a check on the impacts of this assumption, we have also evaluated the sensitivity of the model results to a doubling or halving of host resistance.The effects of a warming climate on Pc risk vary depending on the risk factors and specific climate scenario being assessed. Warming changes winter survival in a straightforward fashion: under the A2 scenario survival increases dramatically so that the region in which more than half of a Pc population would survive the winter increases from 43% to 72% of the study area. More modest increases in winter-survival probabilities arise under the B1 climate scenario, in which Pc winter survival becomes probable over 65% of the study area . The effects of climate change on soil moisture and Pc spring activity are more complex. At the regional scale, climate change reduces the risk posed by Pc across the majority of the study area. However, the changes are spatially variable. Pc risk declines markedly in the Central Valley. It is largely unchanged in coastal northern California and Oregon, where rainfall levels are projected to remain high. Its range is also unchanged in the south-eastern part of the region, which is significantly water limited under contemporary scenarios and projected to remain so. Pc risk increases in the north-eastern extent of the study area. The increase in Pc risk in this area is greatest in the high emissions scenario. In the lower emissions scenario, comparable increases in Pc extent occur inland in the southern extent of the range. In both cases, these increases indicate an interaction of warmer temperatures with unchanged or slightly enhanced rainfall. The potential complexity of the interactions between changing water and temperature in one of these southern locations is illustrated in Fig. 5 for a site in the southwestern part of the region. In this location, Pc risk increases under the B1 scenario but decreases under the A2 scenario. Figure 5 shows a decomposition of the projected changes into those due to temperature and those due to changes in soil moisture. As shown, increasing temperatures increase Pc risk from the baseline case for both A2 and B1 scenarios, but in the A2 scenario, a decrease in soil moisture more than offsets the effect of warmer spring conditions. Conversely, under the B1 conditions, the slight increase in soil moisture increases pathogen risk at this location, but only when both temperature and soil moisture increase together does the large predicted increase in pathogen risk occur. While these threshold-dynamics are not general across the study range, they illustrate the potential for highly nonlinear pathogen responses to interactions in changing temperature and moisture conditions, and highlight the importance of considering the impacts of synchronous changes in climate on pathogen dynamics. As summarized in Table 2 and illustrated in Fig. 2, water limitation reduces Pc risk over a range of 340 000 km2 across the region for the high emissions scenario, with an average decrease in the Pc risk of 0.28. Pc risk is reduced over 40 000 km2 for the low emissions scenario, by 0.26 on average.

That allele pilU1 is found in the United States is consistent with this view

A second related question concerned the relationship among the recombinant IHR group members. In particular, what could be concluded about the origin of the group given the observation from 2 independent studies that the members appear to form a well-defined cluster of genotypes? Third, we used the sequence data to examine the hypothesis that the introgressed DNA was from X. fastidiosa subsp. fastidiosa, the subspecies that causes Pierce’s disease. X. fastidiosa subsp. fastidiosa is native to Central America, and all known isolates in the United States and northern Mexico can be traced back to a single introduced genotype . IHR would be of limited interest if it simply randomized the genetic differences among the subspecies but had a minimal effect on pathogenesis. For this reason, we were particularly interested in documenting any possible invasion of new plant hosts associated with IHR. The hypothesis is that IHR creates a range of novel genotypes that are far more variable than can arise from a lineage diversifying through point mutations, and this diversity facilitates adaptive evolution of a kind not possible for a clonal lineage. This kind of probabilistic evolutionary hypothesis can rarely be directly proven based on an individual case; however, it makes predictions that, if generally supported, would cause the hypothesis to be accepted. In the case of X. fastidiosa, plastic planter pot compelling evidence supporting the hypothesis would be the invasion of a new native host plant that is uniquely associated with IHR. Our data support this hypothesis: in X. fastidiosa subsp. multiplex, IHR is indeed associated with the invasion of at least 2 new native plant hosts, blueberry and blackberry.

To investigate intersubspecific homologous recombination , we analyzed 31 isolates previously identified as IHR-type and 2 isolates previously identified as intermediate-type X. fastidiosa subsp. multiplex , based on sequence of the 7 housekeeping loci used in the MLST scheme defined by Yuan et al. plus a region of the pilU gene. Together, these 33 isolates made up the recombinant group. Details regarding the isolation and typing of the 33 isolates were provided by Nunney et al. , and a summary of salient features is provided in Table S1 in the supplemental material in that article. All sequences used have previously been published and are available both in GenBank and the MLST website . To detect IHR, we employed a modified version of the introgression test developed by Nunney et al. . In its original form, the test compares a set of target sequences, some of which may have been involved in IHR, to a set of potential donor sequences. Each variable site is classified as F, a fixed difference between the target sequences and the donor sequences, or P, a polymorphic site within the target sequences where at least one variant base is shared with the donor set. In the modified version of the test, the targeted introgression test, the target sequence is known a priori and is compared to two references, the donor group, D , and the ancestral group, A . The minimum number of nucleotide differences between the target and the two references defines a ratio of D to A equivalent to the ratio of F to P and can be tested in the same way . In some cases, there is no breakpoint because the whole locus appears to be an introgressed sequence . Although the signal of introgression across the entire sequenced region may be clear, it is valuable to have a statistical test that documents the strength of the signal. In this case, the null expectation is the ratio that reflects the pairwise differences between the donor and ancestral group versus the pairwise differences within the ancestral group .

We used this ratio to define the expectation of the D/A ratio for a chi-square test of complete introgression. Gene diversity and distance trees were calculated using MEGA5 , and the maximum parsimony tree was created using the PARS program in Phylip . Distance trees and the maximum parsimony tree were used rather than other methods, given the known occurrence of intersubspecific recombination in the data. ClonalFrame was used to provide an independent estimate of the relative importance of recombination versus mutation in the recombinant group.Based on 8 loci sequenced , Nunney et al. identified 9 sequence types belonging to the recombinant group of X. fastidiosa subsp. multiplex. These STs all showed evidence of intersubspecific homologous recombination at one or more of the 8 loci and were characterized by 18 alleles, 10 of which were never found in non-IHR X. fastidiosa subsp. multiplex strains . These 10 alleles were examined for evidence of IHR by comparing them to the previously described non-IHR X. fastidiosa subsp. multiplex alleles and to the known X. fastidiosa subsp. fastidiosa and sandyi alleles . Of these 10, 4 alleles were found to be derived in their entirety from X. fastidiosa subsp. fastidiosa, and 3 were found to be chimeric for X. fastidiosa subsp. multiplex and fastidiosa sequences, with significant evidence of one or more recombination breakpoints. These 7 alleles encompassed 4 loci: leuA, cysG, holC, and pilU. The locus most strongly implicated in IHR wascysG, since all of the 9 recombinant-group STs were characterized at this locus by 1 of 3 cysG alleles unique to the group. The involvement of IHR in the genesis of all 3 of these alleles is illustrated by their close genetic relationship to X. fastidiosa subsp.fastidiosa and sandyi alleles . Allele 12, apart from being found in the recombinant group, is an X. fastidiosa subsp. fastidiosa allele . The other two alleles were found to be chimeric: allele 18 contains a single recombinant region at the 3= end of 342 bp, while allele 6 has two short recombinant regions, one at the 5= end of at least 23 bp and another toward the 5= end of at least 35 bp .

The DNA sequence variation defining these patterns is shown in Table 2. The patterns seen in the DNA sequences of the 3 cysGalleles are consistent with the hypothesis of a single IHR that introgressed donor allele 12 into X. fastidiosa subsp. multiplex, followed by subsequent intrasubspecific recombination reintroducing X. fastidiosa subsp. multiplex sequence to create alleles 6 and 18 . There are no inconsistent sites, provided the 5= recombination breakpoint in allele 18 starts after position 71. Introgression of X. fastidiosa subsp. fastidiosa sequence into X. fastidiosa subsp. multiplexwas found in alleles at 3 other loci . In the case of pilU, 7 of the 9 recombinant STs carried either an allele identical to a known X. fastidiosa subsp. fastidiosa allele or 1 bp different from it . Allele 1 is an allele that characterizes most U.S. isolates as well as several STs found in Costa Rica, while allele 9 is unique to the recombinant group. The leuA locus has a single statistically significant recombinant allele, allele 4 . It differed by 2 bp from the X. fastidiosa subsp. fastidiosa allele 9 but by 8 bp from the most similar nonIHR X. fastidiosa subsp. multiplex allele. X. fastidiosa subsp. fastidiosa allele 9 could be the donor for its entirety , although if the recombination region started after site 10 but before position 520 , then only one site would be unexplained. That remaining site carries a base unique to this allele and is probably a novel mutation. If the recombination breakpoint was 3= of position 295 then X. fastidiosa subsp. fastidiosa allele 1 provides as good a match as allele 9 . Similarly, 30 litre plant pots holC allele 7 was also 8 bp different from the most similar non-IHR X. fastidiosa subsp. multiplex allele, providing clear evidence that the 5= end was derived from X. fastidiosa subsp. fastidiosa . The pattern can be explained if X. fastidiosasubsp. fastidiosa allele 19 is the donor of the 5= region ending somewhere between positions 183 and 286, since it leaves no inconsistent bases .Evaluation of the plausibility of a single initial IHR event is complicated by the possibility of subsequent intrasubspecific recombination both within the recombinant group and between the recombinant group and the dominant non-IHR X. fastidiosa subsp. multiplex strains. Plausible sets of recombination events were determined by creating a tree using maximum parsimony applied to the 10 8-locus genotypes . Using allele numbers as characters, there were 2 equally parsimonious trees, each with 14 steps. They differed only in the precise positioning of 22a ; however, assuming a basal introgression of pilU1, only the tree shown in Fig. 3 remained. The hypothetical donor and recipient genotypes were added to root the tree, with the tree dictating gltT3 in the ancestral recipient genotype.

The most parsimonious tree showed that the pattern of introgression was more complex than could result from a single IHR. There are four main events that illustrate this complexity. First, based on this tree, the grouping of STs 27, 28, and 40 is defined by the introgression of holC7, a recombinant allele introduced into the tree far from the basal recombination event. Second, although the mutation of pilU1 could explain the appearance of pilU9, a second introgression of pilU1 would be necessary to account for its appearance in STs 28 and 40. Third, a number of events are necessary to account for the evolution of the cysG locus. While cysG12, an X. fastidiosa subsp. fastidiosa allele introduced in the initial recombination event, could give rise to cysG18 by the introgression of X. fastidiosa subsp. multiplex sequence , this allele appears in two places in the tree, necessitating a lateral transfer within the recombinant group. Despite this complexity, the hypothesis of a single primary IHR event creating the founder of the recombinant group is strongly supported by the pattern seen at the cysGlocus. As noted above, all members of the recombinant group share one of 3 alleles that appear to be derived from a single introgression of donor allele 12. Analysis of theX. fastidiosa subsp. fastidiosa donor.The proposed X. fastidiosa subsp. fastidiosa donor is defined at 4 of the 8 loci: leuA9, cysG12, holC19, and pilU1. Of these 4 alleles, only pilU1 was found in an extensive genetic survey of 86 isolates of X. fastidiosa subsp. fastidiosa within the United States and northern Mexico . The results of this survey, combined with similar genetic data from Costa Rica, led to the conclusion that all isolates of X. fastidiosa subsp. fastidiosa found in North America were derived from a single ancestral strain introduced from Central America . Consistent with this hypothesis was the observation that, in the North American isolates, no allele at the 7 MLST loci orthe pilU locus was more than 1 bp different from the most common allele. Given this background, we can examine the hypothesis that the proposed ancestral donor is consistent with theX. fastidiosa subsp. fastidiosa strains currently found in the United States. Similarly, at leuA there is no inconsistency with U.S. allele 1 if the recombination breakpoint in the recombinant allele 4 was after position 295 . If the breakpoint is before that point, then Costa Rica allele 9 provides a better fit of only 1 bp, a minor difference. In marked contrast, the alleles cysG12 and holC19 have only been found in Costa Rica, not in the United States , and differ markedly from the U.S. alleles. In particular, within the IHR regions, the U.S. X. fastidiosa subsp. fastidiosa alleles cysG1 and holC1 are 5 and 7 bp different, respectively, from the recombinant group sequence, while the Costa Rica alleles precisely match the donor sequence . These large differences require us to reject the hypothesis that the primary X. fastidiosa subsp. fastidiosa donor was derived from the introduced genotype that was the ancestor of all of the North American X. fastidiosa subsp. fastidiosa isolates that have been typed. Estimating recombination rates in the recombinant group of X. fastidiosa subsp. multiplex. The prevalence of recombination over mutation in the evolution of the recombinant group was supported by a ClonalFrame analysis: the estimated ratio of recombination events to mutation was 19,310, with a 95% confidence lower bound of 45.3. Addition of the potential ancestor and/or potential donor genotypes to the analysis maintained high estimates of the lower bound of / . These lower bounds are high for a largely clonal organism, and they illustrate the pervasive involvement of recombination in the genesis of the recombinant group.

Soil drainage was based on depth to water table and hydraulic conductivity

In particular, we highlight the combination of independent and particularly joint effects of climate and soil on trait variation, an interaction that has to date been neglected because few studies include both in a single analysis, at the global scale as we have done here. In doing so, we identify an important gap in knowledge: what is the nature of climate–soil interactions that drive whole-plant trait variation and what distinguishes the majority of climate and soil factors having joint effects on plant traits from those with independent effects? These are the sorts of questions that require answers to increase our capacity to predict plant functional diversity in a changing environment. Such predictive power would contribute to a sound basis for assessing long-term feedbacks between global environmental change and the terrestrial biosphere, helping to constrain parameters of global coupled climate–vegetation models. Humans are currently modifying both climatic and edaphic conditions at the global scale. Climate envelope models used to predict vegetation shifts must be complemented by drivers related to large-scale anthropogenic alterations of soil conditions resulting, for example, from land-use change, atmospheric nitrogen deposition, fertilization, black plastic nursery pots liming and salinization. Our global analysis provides an essential context for finer-scale studies to directly tackle questions of biological processes and mechanisms at landscape and community scales.

Fire is part of the natural disturbance regime of many boreal regions, although recent evidence suggests that anthropogenically induced climate change may be increasing the burned area in North American and Eurasian forests . Because high-latitude ecosystems store approximately 40% of global carbon stocks in biomass and soils , an amount equal to the atmospheric C pool, there has been considerable inTherest in understanding how these systems will respond to climate warming. Combustion of vegetation and forest floor transfers C directly from terrestrial ecosystems to the atmosphere, so increased burned area or fire intensity in the boreal biome could be a strong, positive feedback to atmospheric CO2 concentrations , at least in the early years after fire . At a regional scale, the effect of fire on species composition, soil drainage, and stand age distribution will ultimately regulate whether the CO2 feedback is positive or negative. The response of this long-term signal to the combination of climate change and alThered fire regime is largely unknown for the boreal biome. Patterns of plant species composition, biomass accumulation, and productivity across post-fire succession are important determinants of the amount , structure , residence time , and decomposability of C inputs to these systems. In Interior Alaska, black spruce BSP stands cover approximately 70% of the forested area and occupy landscape positions that range from permafrost-free well-drained soils to permafrost dominated and poorly drained soils . These forests are highly flammable due to their architecture and resin production as well as the thick moss layer on the forest floor , and fire return intervals range from 70 to 100 years . Wildfires tend to be large and high-intensity and while few boreal overstory species survive fire, many understory species re-sprout after fire from buried or protected meristems .

Black spruce are semi-seritonous and release seed after fire, with the majority of trees recruiting in the first 5 years after fire . Stands may or may not go through a deciduous phase where trembling aspen and tall shrubs , which may resprout after fire , dominate over the 15–50 years prior to closure of the black spruce canopy . Moss and lichen expansion across the forest floor follows similar timing, with moss cover reaching its maximum between 30 and 50 years , concurrent with canopy closure by black spruce and a reduction in deciduous litter. The deciduous phase appears to be related to interactive effects of fire severity and site drainage as evidenced by the fact that sites that burn severely , or at a high frequency have the highest abundance of aspen and willow species. At the landscape level, both the severity and frequency of fire appear to be related to soil drainage . Although patterns of species dominance over post-fire succession have been described for Alaskan black spruce stands, there are few published measurements of productivity and biomass after fire. In conjunction with a chronosequence study of soil C dynamics, O‘Neill and others used a mass balance model to conclude that C inputs balanced C losses 7–15 years after fire. Yarie and Billings used forest inventory data from stands across Alaska to estimate generalized biomass accumulation curves for black spruce green timber. They show that biomass accumulation peaked between 75 and 150 years. Simulation modeling of ecosystem C dynamics over post-fire succession suggests that C balance is most sen-sitive to N fixation, moss accumulation, organic layer depth, soil drainage, and fire severity.

Finally, there are several comprehensive studies of post-fire succession in Central and Eastern Canada , but the trees in these sites inhabit different soil drainage and temperature regimes than their Alaskan relatives, potentially resulting in different rates of ecosystem C dynamics . The goals of this study were to describe the changes in community structure and above ground net primary productivity and biomass that occur over post-fire succession in the upland black spruce forests of Interior Alaska. We present measurements that span two different time scales: recovery 1–4 years after fire and recovery over the entire successional cycle. For the former, we followed vegetation recovery for 4 years after the 1999 Donnelly Flats fire near Delta Junction, Alaska. We used a chronosequence approach for the latter by selecting two sequences of sites in the region that varied primarily in time since fire: a mesic sequence on moderately well-drained soil with permafrost and a dry sequence located on well-drained soils without permafrost . These sequences represent transitions in environmental factors that might occur with climate warming, including loss of permafrost and subsequent increases in soil drainage .This study was conducted in the Donnelly Flats area located near Delta Junction in Interior Alaska, in seven upland sites that were previously dominated by black spruce . All sites were located within a 100-km2 area on gently sloped alluvial flats that range from moderately well-drained soils dominated by permafrost to well-drained soils where permafrost was largely absent. Our study include three sites on well-drained soils that burned in stand-killing wildfires in 1999, 1987, and approximately 1921 , hereafter the dry chronosequence, and four sites on moderately well-drained soils that burned in 1999, 1994, 1956, and approximately 1886 , hereafter the mesic chronosequence. Time since last fire was determined by historical record in the younger sites and by tree ring analyses in the older sites. Some or all of these sites have been used to assess the effects of fire on soil C storage and emissions , soil chemistry , hydrogen fluxes , fungal community composition and dynamics , seasonal CO2 and 18O–CO2 fluxes and energy exchange . Within each chronosequence, sites were chosen to have similar state factors other than time . Climate: Micrometerological data collected in the 1999, 1987, and 1921 dry sites and the 1994 , 1999 and 1886 mesic sites support the idea that sites in both chronosequences experienced a similar climatic regime. The regional climate is cold and dry with an annual mean surface air temperature of 2.1C during the 1970–2000 period . Over this same period, mean temperatures in January and July were 20 C and 16.0C, respectively, and mean annual precipitation was 290 mm. Approximately 65% of precipitation fell during June, July, and August. Potential biota: Although all stands were currently or historically dominated by black spruce and were in a close enough proximity that they belong to the same regional pool of potentially colonizing organisms, 30 plant pot the understory vegetation and ground cover varied with soil drainage and stand age .

The oldest dry stand was a lichen woodland , with ground cover dominance split between feather moss and lichens. Vaccinium uliginosum and V. vitis-idaea were the most abundant understory species, with deciduous shrubs and trees, forbs and graminoids present but at low abundance. Many of the same species resprouted or recruited after fire in the 1999 dry site and dominated the understory in the 1987 dry site. Species characteristic of well-drained ecosystems that were present in all dry chronosequence sites and absent from the mesic sites were the grass Festuca altaica and the evergreen shrub Arctostaphylos uva-ursi. These species were present, however, on trails and roadsides around the mesic sites. The oldest mesic stand had continuous feathermoss ground cover and a high abundance of Vaccinium spp. Feathermoss occupied almost the entire ground surface in the 1956 and 1886 mesic sites. In the 1994 mesic site it persisted in patches that appeared to have escaped burning. Vascular nomenclature follows Hulte´n and non-vascular nomenclature follows Vitt and others . Relief: Sites in both chronosequences were within a 100- km2 area with little variation in slope or topography . Parent material: Soils along both chronosequences were mainly derived from the Donnelly moraine and wind blown loess and have been described in detail elsewhere . Differences in drainage between the chronosequences are thought to be related to differences in water table depth and texture . Although great care was taken to control state factors within and between chronosequences, it was difficult to fully constrain the effects of past fires on productivity or biomass pools. In the 1999, 1994, and 1987 sites, fires were stand replacing . In the 1957 mesic site, the relatively small range of tree sizes suggests a single cohort of black spruce. In the mature 1886 mesic and 1921 dry sites where tree sizes are quite variable, however, the number of trees sampled for age was not large enough to determine whether stands are comprised of a single cohort . At the landscape-scale, the severity and frequency of fire are likely to be related to soil drainage . At the site level, however, stochastic factors such as weather conditions, time since last fire, and neighboring vegetation can also affect fire severity. Post-fire vegetation recovery is similarly affected by stochastic processes such as timing of fire in relation to both vegetative and reproductive phenology, proximity of seed source, and/or the effects of past and present climate conditions on demographic processes. Finally, we caution the reader to keep in mind at all times that this is an observational study; we depend on the assumptions of the chronosequence approach to make inferences about time.Above ground biomass of vascular plants, mosses and lichens was measured across all sites by destructive harvest in July 2001 at approximately peak biomass. To more closely examine the dynamics of regrowth in the first several years after fire, biomass was also measured in the 1999 dry site 2 months after the fire as well as mid-summer in 2000–2002; it was also measured in 2000, 2001, and 2002 in the mature dry site for comparison. Trees less than 1.37 m in height that were excluded from the inventory described above were included in these harvests. In harvests of the 1999 dry site, we determined whether each species was a re-sprouter by assessing the presence of charred stems or large rhizomes. We also monitored species or generic richness on an annual basis in these sites by recording the presence of all species within a 144-m2 plot surrounding the 1 m2 harvest blocks. We did not survey species richness in the other sites. In each site, above ground biomass was clipped from either 6 or 12 randomly located 1 m2 quadrats. Organic depth was measured at the four corners of each quadrat and averaged. In the mesic chronosequence sites and the 1987 dry site, above ground biomass of vascular species was clipped from six 1 m2 quadrats randomly located along two 100 m long permanent transects . Mosses and lichens were collected from a 400-cm2 organic soil plug sawed from a randomly selected corner of the 1 m2 quadrat following vascular plant clipping. In the 1987 dry site and the 1994 mesic site, we also harvested tall shrubs in a 4-m2 quadrat surrounding the 1-m2 quadrat to account for their larger stature. Vegetation was harvested similarly in the 1999 and 1921 dry sites, except that 12 quadrats were harvested. Samples were returned to the lab and sorted into species and tissues within 1 day of harvest. Each vascular species was separated into several tissue categories including current year and previous year leaves, current year and previous year stems, and fruits or inflorescences following methods modified from Shaver and Chapin and Chapin and others .

The respective data are too scarce to yet be integrated with global datasets

Individual-level daily egg laying data for three fruit fly species were used for the analysis for the two tephritid fruit fly species including the Mediterranean fruit fly , commonly known as the Medfly , the Mexican fruit fly , commonly known as the Mexfly , and the vinegar fly, Drosophila melanogaster . Husbandry details for each species are described in the above-cited papers. In brief rearing conditions were 25-27° C and 50-75% RH, and 12:12 L:D for the two tephritid species and 10:14 L:D for D. melanogaster. Whereas the two tephritid species were fed a 3:1 sugar-yeast hydrolysate diet, D. melanogaster females were fed a standard agar-gelled Drosophila food medium. Eggs were collected from mesh at one end of the individual cages for the tephritid females and from the food medium for D. melanogaster . These species were chosen primarily for the availability of databases on individual-level lifetime egg laying. However, they were also chosen because they allowed for three levels of comparison to test the robustness of our discoveries—between two dipteran families , between two tephritid genera , and among three different species.Reproductive data were excluded for the first 10 days of adult life for each of the species, a period during which all individuals matured and developed their first eggs. The remaining data for each fly were then parsed into 1-of-2 segment categories: Terminal segment consisting of the sequence of daily reproduction from 11 days before death to and including the day of death. There were 873, 1,071 and 425 terminal segments derived from individual female Medflies, raspberry containers size the Mexflies and D. melanogaster, respectively; and Midlife segments of the same length. The midlife segments were created from contiguous 11-day series of individual flies that started after maturation and ended 11 days before the fly died.

This selection criteria thus excluded individual flies that lived less than 32 days . This segmentation process was repeated until the number of remaining days to the terminalsegment was less than 11. The remaining days were ignored. The 11-day egg-laying sequences were chosen primarily because the results of preliminary investigations revealed that 11-days was the shortest period of egg laying that contained patterns that yielded consistently high performance metrics across all three species. This period was also chosen because longer egg-laying periods would have substantially reduced the number of segments in each of the databases. We computed two metrics from all of the 11-day egg laying sequences for all three species and the terminal and mid-life segments: 1) Total eggs . This metric was chosen based on the observation that most flies laid fewer eggs at the end of life; and egg-laying ratio . This metric was based on the observation that the relative rate of egg laying in most flies decreased at the end of life. This metric was computed as the ratio of the number of eggs produced from 11 to 6 days prior to death to the number of eggs produced from days 5 to 0 prior to death. This ratio specified both a direction and a magnitude of change as an individual fly approached death. Egg ratio was assigned a value of ER=0 if no eggs were laid in either time interval and a value of ER=5.0 if eggs were laid in the first 6 days but none in the last 5 days—a pattern that suggests rapid egg-laying decline to zero but that results in a zero in the denominator . This value was chosen because it was at the mid-range of the highest ER-values when eggs were laid in both first and last segments. This ER value was high enough to serve as a major change metric between the first and last egg-laying subsegments but not so high to be the overriding drivers of the statistical outcomes.Fig. 1 shows the egg laying patterns of selected individual medfly females at 10 different life table deciles.

With the exception of female #5, the egg laying patterns during the terminal segments were consistent with the hypothesis that rate of egg laying near the end of each of their lives was both low and decreasing . However, the egg-laying patterns that characterize the end of female’s life are also sometimes observed at times when they are not approaching death. Indeed, there are also a number of 11-day midlife sequences of some flies that are indistinguishable from these same patterns in the terminal phases. For example, egg laying rates decrease in fly #5 from days 20 to 30 and in fly #8 from days 30 to 40 days. Although these decreasing egg-laying patterns usually predict impending death, both of these flies another 11 days beyond these ages. Similar egg-laying trends are also evident during the midlife of fly #9. Fly #10 produced very few eggs for a 40-day period from ages 20 through 60 days. Thus this visual inspection of egg laying in 10 different individuals reveals the statistical challenge of classifying egg-laying sequences as either terminal or mid-life.Although computing all three of these parameters is straightforward in both chronological or thanatological age, the values for and interpretations of GRR and T differ between the two age categories. Consider the following hypothetical case for clarifying the differences. Suppose that in chronological time, the age-specific egg laying in four female fruit flies was identical for the ages each were alive with reproductive peaks at 50 eggs/day on day 20. However, they each die at different ages 20 days apart—i.e., at ages 20, 40, 60 and 80 days. When their reproductive schedules are considered in thanatological ages, which is to say, relative to their age of death rather than their age of birth , then their reproductive peaks are in thanatological time correspond to ages 0, 20, 40 and 60 days . Although R0 remains the same since every female still produces the same lifetime number of eggs, the values of both GRR and T will in the vast majority of cases be different because the timing of reproduction is relative to age of death rather than to age of birth .

Whereas day 20 was the average age of peak reproduction in chronological time for the hypothetical example, it is day 30 in thanatological time [i.e., /4].A 2×2 contingency table that classifies prediction is shown in Fig. 2 where the numbers along the major diagonal represent the correct decisions made, and the numbers in the other diagonal represent the errors . If the instance is positive and the instance is a positive, it is counted as a true positive. However, if it is classified as a negative it is a false negative. If the instance is negative and the instance is a negative, it is counted as a true negative. However, if it is classified as a positive it is a false positive. This contingency table forms the basis for a number of common metrics given below.Event-history reproductive charts plotted in both chronological and thanatological time for all three fruit fly species are presented in Fig. 3. Several aspects of these graphs merit comment. First, the chronological plots of egg laying patterns for all species reveal the familiar progression starting at eclosion from the pre-reproductive, raspberry plant container maturation period, followed by a periods of high reproduction with relatively high levels of intra-individual and inter-species variability. This period, in turn, is followed by a period of tapering off, the length of which depends largely on an individual’s lifespan. The low levels of egg laying are evident at older ages in all species but is most striking in the oldest D. melanogaster. Second, the event-history plots in thanatological time reveal visually the repositioning of reproduction that occurs when the schedule is normalized with respect death rather than birth . This shift is especially evident in diagonal band of high reproduction that tracks to the left of the cohort survival the curve. In these cases the most advanced ages in thanatological time correspond the very youngest and thus most fecund ages in chronological time. Third, patterns of egg laying near death as seen in all three species and for both time frames differ between short-lived and long-lived individuals. This is outcome of the differences in the underlying “causes” of death at young and old ages. Increasing frailty due to old age is the most likely “cause” of death in the longest lived individuals. This accounts for the progressive decrease in egg production at older ages. However, increasing frailty due to aging is an unlikely “cause” of death for flies that die young and at ages when they are at or near their peak in egg production. Thus no single egg laying pattern or combination of patterns will likely ever apply to all flies regardless of age or cause of death.Comparisons of the average reproductive rates and timing for all three species plotted with respect to both chronological and thanatological ages are given in the different panels shown in Fig. 4. Because the number of eggs laid by the average female in her lifetime is the same regardless of whether the eggs are summed from birth to death or from death to birth, net reproductive rates, R0, will be the same regardless of whether it is considered in chronological or thanatological time. But as noted in Methods section, this not the case for gross reproductive rate computations.

For example, the 50% greater value of GRR for thanatological age relative to the value of this metric for chronological age in D. melanogaster is the result of a subset of extremely long-lived individuals who both matured early and produced many eggs at young ages. When re-plotted these young ages in chronological time represent the “old” ages in thanatological time. Thus individuals who are both long-lived and highly fecund at young chronological ages represent a large fraction of the small number at the tail end of the “death” cohort. Differences in the values of the mean age of parenthood, T, across species for chronological versus thanatological ages revealed that it was 8 days closer to birth than to death in D. melanogaster, but slightly over 4 days closer to death than to birth in the Mexfly and nearly equidistant from birth and death in the Medfly .The means and frequency distributions of the two independent variables we use in the regression model for each of the fly species are shown in the series of plots contained in Fig. 5. These graphs anticipate the outcome of the modelling results by revealing the differences between the metrics in the midlife segments relative to the terminal segments. There were striking differences in the metrics for each of the two categories of segments in D. melanogaster with nearly 5-fold fewer eggs and 3-fold greater average egg ratio in the terminal segments. The signs of the differences in the mean and overall distributions for the two tephritid species were similar but the magnitudes of the differences were much less. We thus anticipate more favorable performance metrics for distinguishing between terminal and mid-life egg laying patterns in D. melanogaster than in the two tephritid species.The logistic regression model yields three general results. First, the model’s overall performance supports the concept that the egg-laying patterns of total eggs and egg ratio as individual flies approached death are, in the majority of cases, distinctly different from those patterns over an 11-day sequence in middle of their lives. This wasevident in the performance metric of fraction correct —i.e., the FC-value for all species exceeded 0.64 using all data and exceeded 0.73 using only the segments in which flies laid 25 eggs or greater . Second, with a minor difference for several of the D. melanogaster metrics, the performance of the regression model was more favorable when applied to the censored data than with use of the uncensored data. The reason for the differences was because there were 11-day egg-laying sequences in which few or no eggs were produced or that egg laying was declining in midlife. These are the patterns for the independent variables that were associated with and thus predictive of the terminal segments. These midlife patterns occurred more in in the two tephritid species than in D. melanogaster and thus helps explain the higher performance levels for the regression model in this species. Third, the performance metrics for D. melanogaster were extraordinarily high relative to those for the two tephritids as well as absolute.