More stories

  • in

    Grazing pressure on drylands

    Maestre and colleagues collected data using a standardized field survey at 98 sites across 25 countries and 6 continents, fitted linear mixed models to data from all sites and grazing pressure levels, and then applied a multimodel inference procedure to select the set of best-fitting models. The authors found interactions between grazing and biodiversity in almost half of the best-fitting models, where increasing grazing pressure had positive effects on ecosystem services in colder sites with high plant species richness. However, increases in grazing pressure at warmer sites with high rainfall seasonality and low plant species richness interacted with soil properties to either increase or reduce the delivery of multiple ecosystem services. The authors’ findings highlight how increasing herbivore richness could enhance ecosystem service delivery across contrasting environmental and biodiversity conditions, enhancing soil carbon storage and reducing the negative impacts of increased grazing pressure. More

  • in

    Reproductive performance and sex ratio adjustment of the wild boar (Sus scrofa) in South Korea

    Estes, J. A. Predators and ecosystem management. Wildl. Soc. Bull. 24, 390–396 (1996).
    Google Scholar 
    Licht, D. S., Millspaugh, J. J., Kunkel, K. E., Kochanny, C. O. & Peterson, R. O. Using small populations of wolves for ecosystem restoration and stewardship. Bioscience 60, 147–153 (2010).Article 

    Google Scholar 
    Schwartz, C. C., Swenson, J. E. & Miller, S. D. Large carnivores, moose, and humans: A changing paradigm of predator management in the 21st century. Alces J. Devot. Biol. Manag. Moose 39, 41–63 (2003).
    Google Scholar 
    Valente, A. M., Acevedo, P., Figueiredo, A. M., Fonseca, C. & Torres, R. T. Overabundant wild ungulate populations in Europe: Management with consideration of socio-ecological consequences. Mammal Rev. 50, 353–366 (2020).Article 

    Google Scholar 
    Lee, S. D. & Miller-Rushing, A. J. Degradation, urbanization, and restoration: A review of the challenges and future of conservation on the Korean Peninsula. Biol. Cons. 176, 262–276 (2014).Article 

    Google Scholar 
    Kodera, Y. Habitat selection of Japanese wild boar in Iwami district, Shimane Prefecture, western Japan. Wildl. Conserv. Jpn. 6, 119–129 (2001).
    Google Scholar 
    Ohashi, H. et al. Differences in the activity pattern of the wild boar Sus scrofa related to human disturbance. Eur. J. Wildl. Res. 59, 167–177 (2013).Article 

    Google Scholar 
    Ministry of Environment Republic of Korea. Management Plan of Pest Wild Boars. (Sejong, 2010).National Institute of Biological Resources. Analysis of Hunting Effect on Pest Animals and its Monitoring Techniques. (Incheon, 2017).Lee, S. M. & Lee, E. J. Diet of the wild boar (Sus scrofa): Implications for management in forest-agricultural and urban environments in South Korea. PeerJ 7, e7835 (2019).Article 

    Google Scholar 
    Bieber, C. & Ruf, T. Population dynamics in wild boar Sus scrofa: ecology, elasticity of growth rate and implications for the management of pulsed resource consumers. J. Appl. Ecol. 42, 1203–1213 (2005).Article 

    Google Scholar 
    Brogi, R. et al. Capital-income breeding in wild boar: A comparison between two sexes. Sci. Rep. 11, 1–10 (2021).Article 

    Google Scholar 
    Frauendorf, M., Gethöffer, F., Siebert, U. & Keuling, O. The influence of environmental and physiological factors on the litter size of wild boar (Sus scrofa) in an agriculture dominated area in Germany. Sci. Total Environ. 541, 877–882 (2016).Article 
    ADS 
    CAS 

    Google Scholar 
    Massei, G., Genov, P. V. & Staines, B. W. Diet, food availability and reproduction of wild boar in a Mediterranean coastal area. Acta Theriol. 41, 307–320 (1996).Article 

    Google Scholar 
    Sabrina, S., Jean-Michel, G., Carole, T., Serge, B. & Eric, B. Pulsed resources and climate-induced variation in the reproductive traits of wild boar under high hunting pressure. J. Anim. Ecol. 78, 1278–1290 (2009).Article 

    Google Scholar 
    Fonseca, C. et al. Reproduction in the wild boar (Sus scrofa Linnaeus, 1758) populations of Portugal. Galemys 16, 53–65 (2004).
    Google Scholar 
    Morreti, M. Birth distribution, structure and dynamics of a hunted mountain populatin of wild boars (Sus scrofa L.), Ticino, Switzerland. J. Mt. Ecol. 3, 192–196 (1995).
    Google Scholar 
    Rosell, C., Navas, F. & Romero, S. Reproduction of wild boar in a cropland and coastal wetland area: Implications for management. Anim. Biodivers. Conserv. 35, 209–217 (2012).Article 

    Google Scholar 
    Cellina, S. Effects of supplemental feeding on the body condition and reproductive state of wild boar Sus scrofa in Luxembourg PhD Dissertation thesis, University of Sussex, (2008).Gamelon, M. et al. High hunting pressure selects for earlier birth date: Wild boar as a case study. Evol. Int. J. Org. Evol. 65, 3100–3112 (2011).Article 

    Google Scholar 
    Gethöffer, F., Sodeikat, G. & Pohlmeyer, K. Reproductive parameters of wild boar (Sus scrofa) in three different parts of Germany. Eur. J. Wildl. Res. 53, 287–297 (2007).Article 

    Google Scholar 
    Fonseca, C., Da Silva, A., Alves, J., Vingada, J. & Soares, A. Reproductive performance of wild boar females in Portugal. Eur. J. Wildl. Res. 57, 363–371 (2011).Article 

    Google Scholar 
    Gaillard, J.-M., Brandt, S. & Jullien, J.-M. Body weight effect on reproduction of young wild boar (Sus scrofa) females: A comparative analysis. Folia Zool. (Brno) 42, 204–212 (1993).
    Google Scholar 
    Poteaux, C. et al. Socio-genetic structure and mating system of a wild boar population. J. Zool. 278, 116–125 (2009).Article 

    Google Scholar 
    Spitz, F., Valet, G. & Lehr Brisbin, I. Jr. Variation in body mass of wild boars from southern France. J. Mammal. 79, 251–259 (1998).Article 

    Google Scholar 
    Trivers, R. L. & Willard, D. E. Natural selection of parental ability to vary the sex ratio of offspring. Science 179, 90–92 (1973).Article 
    ADS 
    CAS 

    Google Scholar 
    Clutton-Brock, T. H., Albon, S. & Guinness, F. Parental investment in male and female offspring in polygynous mammals. Nature 289, 487–489 (1981).Article 
    ADS 

    Google Scholar 
    Hewison, A. M. & Gaillard, J.-M. Successful sons or advantaged daughters? The Trivers-Willard model and sex-biased maternal investment in ungulates. Trends Ecol. Evol. 14, 229–234 (1999).Article 
    CAS 

    Google Scholar 
    Clutton-Brock, T. H. & Iason, G. R. Sex ratio variation in mammals. Q. Rev. Biol. 61, 339–374 (1986).Article 
    CAS 

    Google Scholar 
    Fernández-Llario, P. & Mateos-Quesada, P. Body size and reproductive parameters in the wild boar Sus scrofa. Acta Theriol. 43, 439–444 (1998).Article 

    Google Scholar 
    Meikle, D. B., Drickamer, L. C., Vessey, S. H., Rosenthal, T. L. & Fitzgerald, K. S. Maternal dominance rank and secondary sex ratio in domestic swine. Anim. Behav. 46, 79–85 (1993).Article 

    Google Scholar 
    Servanty, S., Gaillard, J.-M., Allainé, D., Brandt, S. & Baubet, E. Litter size and fetal sex ratio adjustment in a highly polytocous species: The wild boar. Behav. Ecol. 18, 427–432 (2007).Article 

    Google Scholar 
    Mendl, M., Zanella, A. J., Broom, D. M. & Whittemore, C. T. Maternal social status and birth sex ratio in domestic pigs: An analysis of mechanisms. Anim. Behav. 50, 1361–1370 (1995).Article 

    Google Scholar 
    Cameron, E. Z. Facultative adjustment of mammalian sex ratios in support of the Trivers-Willard hypothesis: Evidence for a mechanism. Proc. R. Soc. Lond. Ser. B Biol. Sci. 271, 1723–1728 (2004).Article 

    Google Scholar 
    Clutton-Brock, T., Albon, S. & Guinness, F. Maternal dominance, breeding success and birth sex ratios in red deer. Nature 308, 358–360 (1984).Article 
    ADS 

    Google Scholar 
    Maillard, D. & Fournier, P. Timing and synchrony of births in the wild boar (Sus scrofa Linnaeus, 1758) in a Mediterranean habitat: The effect of food availability. Galemys 16, 67–74 (2004).
    Google Scholar 
    Bywater, K. A., Apollonio, M., Cappai, N. & Stephens, P. A. Litter size and latitude in a large mammal: the wild boar Sus scrofa. Mammal Rev. 40, 212–220 (2010).
    Google Scholar 
    Orłowska, L., Rembacz, W. & Florek, C. Carcass weight, condition and reproduction of wild boars harvested in north-western Poland. Pest Manag. Sci. 69, 367–370 (2013).Article 

    Google Scholar 
    Carranza, J. Sexual selection for male body mass and the evolution of litter size in mammals. Am. Nat. 148, 81–100 (1996).Article 

    Google Scholar 
    FernáNdez-Llario, P., Carranza, J. & Mateos-Quesada, P. Sex allocation in a polygynous mammal with large litters: The wild boar. Anim. Behav. 58, 1079–1084 (1999).Article 

    Google Scholar 
    McBride, G. The” teat order” and communication in young pigs. Animal Behaviour (1963).McBride, G., James, J. & Hodgens, N. Social behaviour of domestic animals. IV. Growing pigs. Anim. Sci. 6, 129–139 (1964).Article 

    Google Scholar 
    McBride, G., James, J. & Wyeth, G. Social behaviour of domestic animals VII. Variation in weaning weight in pigs. Anim. Sci. 7, 67–74 (1965).Article 

    Google Scholar 
    Geochang County. Geochang Statistical yearbook. (Geochang, 2015).Seoul Metropolitan Government. Seoul Statistical Yearbook. (Seoul, 2017).Animal and Plant Quarantine Agency and Ministry of Food and Drug Safety. Institutional Animal Care and Use Committee Standard Operation Guideline. (Gimcheon, 2020).Magnell, O. & Carter, R. The chronology of tooth development in wild boar – A guide to age determination of linear enamel hypoplasia in prehistoric and medieval pigs. Verterrinarija Ir Zootechnika. T. 40, 43–48 (2007).
    Google Scholar 
    Oroian, T. E., Oroian, R. G., Pasca, I., Oroian, E. & Covrig, I. Methods of age estimation by dentition in Sus scrofa ferus sp. Bull. UASVM Anim. Sci. Biotechnol. 67, 1–2 (2010).
    Google Scholar 
    Vericad, R. Fetal age, conception and birth period estimation on wild boar (Sus scrofa) in West Pyrenees. in Actas del XV Congresso Int. Fauna Cinegética y Silvestre. (Trujillo, 1983).Rosell, C., Fernández-Llario, P. & Herrero, J. The wild boar (Sus scrofa Linnaeus, 1758). Galemys 13, 1–25 (2001).
    Google Scholar 
    R core team. R: A language and environment for statistical computing v. 3.6.0 (Austria, 2019). More

  • in

    Impacts of soil nutrition on floral traits, pollinator attraction, and fitness in cucumbers (Cucumis sativus L.)

    Fichtner, K. & Schulze, E. D. The effect of nitrogen nutrition on growth and biomass partitioning of annual plants originating from habitats of different nitrogen availability. Oecologia 92, 236–241 (1992).Article 
    ADS 
    CAS 

    Google Scholar 
    Rodger, J. G. et al. Widespread vulnerability of flowering plant seed production to pollinator declines. Sci. Adv. 7, eabd3524. https://doi.org/10.1126/sciadv.abd3524 (2021).Article 
    ADS 

    Google Scholar 
    de Groot, C. C., Marcelis, L. F. M., van den Boogaard, R., Kaiser, W. M. & Lambers, H. Interaction of nitrogen and phosphorus nutrition in determining growth. Plant Soil 248, 257–268 (2003).Article 

    Google Scholar 
    Wang, Z. & Li, S. Effects of nitrogen and phosphorus fertilization on plant growth and nitrate accumulation in vegetables. J. Plant Nutr. 27, 539–556 (2004).Article 
    CAS 

    Google Scholar 
    Razaq, M., Zhang, P. & Shen, H. L. Influence of nitrogen and phosphorous on the growth and root morphology of Acer mono. PLoS One 12, e0171321. https://doi.org/10.1371/journal.pone.0171321 (2017).Article 
    CAS 

    Google Scholar 
    Poulton, J. L., Bryla, D., Koide, R. T. & Stephenson, A. G. Mycorrhizal infection and high soil phosphorus improve vegetative growth and the female and male functions in tomato. New Phytol. 154, 255–264 (2002).Article 
    CAS 

    Google Scholar 
    Burkle, L. A. & Irwin, R. E. The effects of nutrient addition on floral characters and pollination in two subalpine plants, Ipomopsis aggregata and Linum lewisii. Plant Ecol. 203, 83–98 (2009).Article 

    Google Scholar 
    Burkle, L. A. & Irwin, R. E. Beyond biomass: measuring the effects of community-level nitrogen enrichment on floral traits, pollinator visitation and plant reproduction. J. Ecol. 98, 705–717 (2010).Article 

    Google Scholar 
    Hoover, S. E. R. et al. Warming, CO2, and nitrogen deposition interactively affect a plant-pollinator mutualism. Ecol. Lett. 15, 227–234 (2012).Article 

    Google Scholar 
    Lau, T. C. & Stephenson, A. G. Effects of soil nitrogen on pollen production, pollen grain size, and pollen performance in Cucurbita pepo (Cucurbitaceae). Am. J. Bot. 80, 763–768 (1993).Article 
    CAS 

    Google Scholar 
    Lau, T. C. & Stephenson, A. Effects of soil phosphorus on pollen production, pollen size, pollen phosphorus content, and the ability to sire seeds in Cucurbita pepo (Cucurbitaceae). Sex. Plant Reprod. 7, 215–220 (1994).Article 

    Google Scholar 
    Atasay, A., Akgül, H., Uçgun, K. & Şan, B. Nitrogen fertilization affected the pollen production and quality in apple cultivars ‘Jerseymac’ and ‘Golden Delicious’. Acta Agric. Scand. Sect. B. Soil Plant Sci. 63, 460–465 (2013).
    Google Scholar 
    Shuel, R. W. Some aspects of the relation between nectar secretion and nitrogen, phosphorus, and potassium nutrition. Can. J. Plant Sci. 37, 220–236 (1957).Article 
    CAS 

    Google Scholar 
    Robacker, D. C., Flottum, P. K., Sammataro, D. & Erickson, E. H. Effects of climatic and edaphic factors on soybean flowers and on the subsequent attractiveness of the plants to honey bees. Field Crops Res. 6, 267–278 (1983).Article 

    Google Scholar 
    Dror, I., Yaron, B. & Berkowitz, B. The human impact on all soil-forming factors during the anthropocene. ACS Environ. Au 2, 11–19 (2022).Article 
    CAS 

    Google Scholar 
    David, T. I., Storkey, J. & Stevens, C. J. Understanding how changing soil nitrogen affects plant–pollinator interactions. Arthropod. Plant Interact. 13, 671–684 (2019).Article 

    Google Scholar 
    Russo, L., Buckley, Y. M., Hamilton, H., Kavanagh, M. & Stout, J. C. Low concentrations of fertilizer and herbicide alter plant growth and interactions with flower-visiting insects. Agric. Ecosyst. Environ. 304, 107141. https://doi.org/10.1016/j.agee.2020.107141 (2020).Article 
    CAS 

    Google Scholar 
    Akter, A. & Klečka, J. Water stress and nitrogen supply affect floral traits and pollination of the white mustard, Sinapis alba (Brassicaceae). PeerJ 10, e13009. https://doi.org/10.7717/peerj.13009 (2022).Article 
    CAS 

    Google Scholar 
    Wu, Y. et al. Soil water and nutrient availability interactively modify pollinator-mediated directional and correlational selection on floral display. New Phytol. https://doi.org/10.1111/nph.18537 (2022).Article 

    Google Scholar 
    Nicolson, S. W. Sweet solutions: nectar chemistry and quality. Philos. Trans. R. Soc. Lond. B Biol. Sci. 377, 2163. https://doi.org/10.1098/rstb.2021.0163 (2022).Article 
    CAS 

    Google Scholar 
    Vaudo, A. D., Tooker, J. F., Grozinger, C. M. & Patch, H. M. Bee nutrition and floral resource restoration. Curr. Opin. Insect Sci. 10, 133–141 (2015).Article 

    Google Scholar 
    Cnaani, J., Thomson, J. D. & Papaj, D. R. Flower choice and learning in foraging bumblebees: effects of variation in nectar volume and concentration. Ethology 112, 278–285 (2006).Article 

    Google Scholar 
    Vaudo, A. D., Patch, H. M., Mortensen, D. A., Tooker, J. F. & Grozinger, C. M. Macronutrient ratios in pollen shape bumble bee (Bombus impatiens) foraging strategies and floral preferences. Proc. Natl. Acad. Sci. U. S. A. 113, E4035–E4042. https://doi.org/10.1073/pnas.1606101113 (2016).Article 
    ADS 
    CAS 

    Google Scholar 
    Vaudo, A. D. et al. Pollen protein: lipid macronutrient ratios may guide broad patterns of bee species floral preferences. Insects 11, 132. https://doi.org/10.3390/insects11020132 (2020).Article 

    Google Scholar 
    Cardoza, Y. J., Harris, G. K. & Grozinger, C. M. Effects of soil quality enhancement on pollinator-plant interactions. Psyche 2012, 581458. https://doi.org/10.1155/2012/581458 (2012).Article 

    Google Scholar 
    Ceulemans, T., Hulsmans, E., Vanden Ende, W. & Honnay, O. Nutrient enrichment is associated with altered nectar and pollen chemical composition in Succisa pratensis Moench and increased larval mortality of its pollinator Bombus terrestris L.. PLoS One 12, e0175160. https://doi.org/10.1371/journal.pone.0175160 (2017).Article 
    CAS 

    Google Scholar 
    Russo, L., Vaudo, A. D., Fisher, C. J., Grozinger, C. M. & Shea, K. Bee community preference for an invasive thistle associated with higher pollen protein content. Oecologia 190, 901–912 (2019).Article 
    ADS 

    Google Scholar 
    Russo, L., Keller, J., Vaudo, A. D., Grozinger, C. M. & Shea, K. Warming increases pollen lipid concentration in an invasive thistle, with minor effects on the associated floral-visitor community. Insects 11, 20. https://doi.org/10.3390/insects11010020 (2019).Article 

    Google Scholar 
    Awmack, C. S. & Leather, S. R. Host plant quality and fecundity in herbivorous insects. Annu. Rev. Entomol. 47, 817–844 (2002).Article 
    CAS 

    Google Scholar 
    Carisey, N. & Bauce, E. Does nutrition-related stress carry over to spruce budworm, Choristoneura fumiferana (Lepidoptera: Tortricidae) progeny?. Bull. Entomol. Res. 92, 101–108 (2002).Article 
    CAS 

    Google Scholar 
    Zhang, G. & Han, X. N: P stoichiometry in Ficus racemosa and its mutualistic pollinator. J. Plant Ecol. 3, 123–130 (2010).Article 

    Google Scholar 
    Visanuvimol, L. & Bertram, S. M. How dietary phosphorus availability during development influences condition and life history traits of the cricket Acheta domesticas. J. Insect Sci. 11, 63. https://doi.org/10.1673/031.011.6301 (2011).Article 

    Google Scholar 
    Dovrat, G., Meron, E., Shachak, M., Golodets, C. & Osem, Y. Plant size is related to biomass partitioning and stress resistance in water-limited annual plant communities. J. Arid Environ. 165, 1–9 (2019).Article 
    ADS 

    Google Scholar 
    Bobbink, R. et al. Global assessment of nitrogen deposition effects on terrestrial plant diversity: a synthesis. Ecol. Appl. 20, 30–59 (2010).Article 
    CAS 

    Google Scholar 
    Tao, L. & Hunter, M. D. Does anthropogenic nitrogen deposition induce phosphorus limitation in herbivorous insects?. Glob. Chang. Biol. 18, 1843–1853 (2012).Article 
    ADS 

    Google Scholar 
    Tognetti, P. M. et al. Negative effects of nitrogen override positive effects of phosphorus on grassland legumes worldwide. Proc. Natl. Acad. Sci. 118(28), e2023718118. https://doi.org/10.1073/pnas.2023718118 (2021).Article 
    CAS 

    Google Scholar 
    Leghari, S. J. et al. Role of nitrogen for plant growth and development: a review. Adv. Environ. Biol. 10, 209–218 (2016).
    Google Scholar 
    Carvalheiro, L. G. et al. Soil eutrophication shaped the composition of pollinator assemblages during the past century. Ecography 43, 209–221 (2020).Article 

    Google Scholar 
    Lefcheck, J. S. Piecewisesem: piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).Article 

    Google Scholar 
    Roulston, T. H., Cane, J. H. & Buchmann, S. L. What governs protein content of pollen: Pollinator preferences, pollen–pistil interactions, or phylogeny?. Ecol. Monogr. 70, 617–643 (2000).
    Google Scholar 
    Pacini, E. & Hesse, M. Pollenkitt—its composition, forms and functions. Flora 200, 399–415 (2005).Article 

    Google Scholar 
    Vaudo, A. D. et al. Bumble bees regulate their intake of essential protein and lipid pollen macronutrients. J. Exp. Biol. 219, 3962–3970 (2016).CAS 

    Google Scholar 
    Vaudo, A. D., Farrell, L. M., Patch, H. M., Grozinger, C. M. & Tooker, J. F. Consistent pollen nutritional intake drives bumble bee (Bombus impatiens) colony growth and reproduction across different habitats. Ecol. Evol. 8, 5765–5776 (2018).Article 

    Google Scholar 
    Treanore, E. D., Vaudo, A. D., Grozinger, C. M. & Fleischer, S. J. Examining the nutritional value and effects of different floral resources in pumpkin agroecosystems on Bombus impatiens worker physiology. Apidologie 50, 542–552 (2019).Article 

    Google Scholar 
    Baker, H. G. & Baker, I. The predictive value of nectar chemistry to the recognition of pollinator types. Israel J. Bot. 39, 157–166 (1990).CAS 

    Google Scholar 
    Thomson, J. D. Pollen transport and deposition by bumble bees in Erythronium: influences of floral nectar and bee grooming. J. Ecol. 74, 329–341 (1986).Article 

    Google Scholar 
    Gonzalez, M. V., Coque, M. & Herrero, M. Influence of pollination systems on fruit set and fruit quality in kiwifruit (Actinidia deliciosa). Ann. Appl. Biol. 132, 349–355 (1998).Article 

    Google Scholar 
    Morandin, L. A., Laverty, T. M. & Kevan, P. G. Effect of bumble bee (Hymenoptera: Apidae) pollination intensity on the quality of greenhouse tomatoes. J. Econ. Entomol. 94, 172–179 (2001).Article 
    CAS 

    Google Scholar 
    Karron, J. D., Mitchell, R. J. & Bell, J. M. Multiple pollinator visits to Mimulus ringens (Phrymaceae) flowers increase mate number and seed set within fruits. Am. J. Bot. 93, 1306–1312 (2006).Article 

    Google Scholar 
    Kiatoko, N., Raina, S. K., Muli, E. & Mueke, J. Enhancement of fruit quality in Capsicum annum through pollination by Hypotrigona gribodoi in Kakamega Western Kenya. Entomol. Sci. 17, 106–110 (2014).Article 

    Google Scholar 
    Abrol, D. P., Gorka, A. K., Ansari, M. J., Al-Ghamdi, A. & Al-Kahtani, S. Impact of insect pollinators on yield and fruit quality of strawberry. Saudi J. Biol. Sci. 26, 524–530 (2019).Article 

    Google Scholar 
    Osman, M. A., Raju, P. S. & Peacock, J. M. The effect of soil temperature, moisture and nitrogen on Striga asiatica (L.) Kuntze seed germination, viability and emergence on sorghum (Sorghum bicolor L. Moench) roots under field conditions. Plant Soil 131, 265–273 (1991).Article 
    CAS 

    Google Scholar 
    Rose, T. J. & Raymond, C. A. Seed phosphorus effects on rice seedling vigour in soils differing in phosphorus status. Agronomy 10(12), 1919. https://doi.org/10.3390/agronomy10121919 (2020).Article 
    CAS 

    Google Scholar 
    Cavatorta, J. et al. ‘Marketmore 97’: a monoecious slicing cucumber inbred with multiple disease and insect resistances. HortScience 42, 707–709 (2007).Article 

    Google Scholar 
    Friedman, J. The evolution of annual and perennial plant life histories: ecological correlates and genetic mechanisms. Annu. Rev. Ecol. Evol. Syst. 51, 461–481 (2020).Article 

    Google Scholar 
    Alzate-Marin, A. L. et al. Warming and elevated CO2 induces changes in the reproductive dynamics of a tropical plant species. Sci. Total Environ. 768, 144899. https://doi.org/10.1016/j.scitotenv.2020.144899 (2021).Article 
    ADS 
    CAS 

    Google Scholar 
    Mu, J. et al. Domesticated honey bees evolutionarily reduce flower nectar volume in a Tibetan lotus. Ecology 95, 3161–3172 (2014).Article 

    Google Scholar 
    Cruden, R. W. Pollen-ovule ratios: a conservative indicator of breeding systems in flowering plants. Evolution 31, 32–46 (1977).
    Google Scholar 
    Costa, C. M. & Yang, S. Counting pollen grains using readily available, free image processing and analysis software. Ann. Bot. 104, 1005–1010 (2009).Article 

    Google Scholar 
    Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).Article 
    CAS 

    Google Scholar 
    Vaudo, A. D., Patch, H. M., Mortensen, D. A., Grozinger, C. M. & Tooker, J. F. Bumble bees exhibit daily behavioral patterns in pollen foraging. Arthropod. Plant. Interact. 8, 273–283 (2014).
    Google Scholar  More

  • in

    Vulcanimicrobium alpinus gen. nov. sp. nov., the first cultivated representative of the candidate phylum “Eremiobacterota”, is a metabolically versatile aerobic anoxygenic phototroph

    Sampling, bacterial isolation, and colony screeningSamples were collected from thermally active sediments as evidenced by temperature and/or emission of steam in November 2010 and 2012 from Harry’s Dream (HD2 and HD3), Warren Cave (WC1, WC2, WC3, WC4, WC7, and WC8), Haggis Hole, (HH), Mammoth Cave (MC), Hut Cave (Hut), and Heroine Cave (HC) (fumarolic ice caves on Mt. Erebus volcano [21]). The samples were divided into four types: Mainly weathered basaltic/phonolitic sand (HD2, HD3, WC3, WC4, WC8 and Hut); pebbles and rock fragments (WC7, HH, and MC); black porous glassy materials that appeared to be solidified lava (WC1, WC2, and HC); and ash sediment (WC10) with little to no organic material (Table S1). The sediments were collected aseptically using sterile 50 mm conical tubes and immediately sealed. Additional information and environmental parameters on the sampling locations are provided in Table S1 and in Tebo et al. [21].To isolate bacteria found in these oligotrophic environments, we employed a culture strategy of long-term incubation in a nutrient-poor medium and screening of slow-growing colonies by direct PCR identification. Reasoner’s 2 A gellan gum medium (10% R2AG) [23] and FS1VG medium [24] were used for bacterial isolation. The 10% R2AG is a 10-fold diluted R2A broth (Nihon Seiyaku, Tokyo, Japan), solidified with 15 g/L gellan gum (Kanto chemical, Tokyo, Japan) with 2 g/L CaCl2. The 10% R2AG and FS1VG media were adjusted to a pH of 4.5 or 6.0, with or without 30 mg/L sodium azide. Samples used for isolation were selected from Warren Cave and Harry’s Dream, where the presence of relatively large numbers of bacteria ( >107/g) was confirmed in a previous report [21], and a total of nine samples, WD1, 2, 3, 4, 7, 8, and 10 and HD2 and 3, were used without any other pretreatment such as drying or dilution. For the sandy and ash samples with fine particles (WC3, 4, 8, 10, and HD2, 3), approximately 50 mg of sample were spread directly on plates. Glassy materials (WC1 and WC2) were embedded directly in plates using 2-3 pieces (approx. 200 mg). Pebbles and rock fragments (WC7) were crushed and approximately 50 mg of debris were spread directly onto the plates. All plates were incubated at 15 °C, 30 °C or 37 °C under dark conditions. New colonies were marked as they appeared and selection of the isolates was performed by picking only colonies that appeared after four weeks of incubation. We selected colonies around the sediment or in the cracks created during spreading using a magnifying glass.For identification, colonies were picked with a sterile toothpick, re-streaked, stabbed on fresh medium, and subsequently suspended in 20 μL sterilized 0.05 M NaOH. Suspensions were heated at 100 °C for 15 min, and supernatants were used as template DNAs for PCR. Partial 16 S rRNA gene sequences were amplified by PCR using commonly used bacterial primer set 27 F (5′-AGATTTGATCCTGGCTCAG-3′) and 1492 R (5′-GGTTACCTTGTTACGACTT-3′) or 536 R (5′-GTA TTA CCG CGG CTG CTG-3′) with TaKaRa ExTaq DNA polymerase (Takara Bio, Shiga, Japan) as previously described [23]. Sequencing was performed at Eurofins Genomics (Louisville, KY, USA), using a 3730xl DNA analyzer (Applied Biosystems, CA, USA). Sequence similarities with closest species were calculated using EZbiocloud’s Identify Service (https://www.ezbiocloud.net/identify). Subsequently, cells in the stab identified as “Ca. Eremiobacterota” by the direct PCR identification were serially diluted and stabbed onto new plates until multiple pure cultures of Eremiobacterota were obtained. The isolate was designated as strain WC8-2.Whole genome sequencing and annotationGenomic DNA was extracted from WC8-2 cells grown in 10% R2A broth (pH6.0) with air/CO2 (90:10, v/v) at 30 °C for 30 days under a 12/12 h light/dark regime with incandescent light (250 μmol m–2s–1), using a Puregene Yeast/Bact. Kit B (Qiagen, Germantown, MD, USA) [25]. Sequencing was performed by Macrogen Japan Corp., on a NovaSeq 6000 (Illumina, Inc., San Diego, CA, USA) and PacBio RSII (Pacific Biosciences of California, Inc. Menlo Park, CA, USA). The gap-free complete genome was assembled de novo using the Unicycler version 0.4.8 hybrid assembly pipeline with default settings [26]. Completeness and contamination levels were estimated using CheckM [27]. The genome was annotated using the DDBJ Fast Annotation and Submission Tool (DFAST) [28] and the BlastKOALA web server version 2.2, and was visualized using CGView Server [29] (http://cgview.ca/).Phylogenetic analysis of “Ca. Eremiobacterota”Identification of strain WC8-2 was performed using the Genome Taxonomy Database Toolkit (GTDB-Tk) (ver. 2.1.0), which produces standardized taxonomic labels that are based on those used in the Genome Taxonomy Database [30]. Terrabacterial genomes including “Ca. Eremiobacterota” MAGs and related genomes were retrieved from the Genome Taxonomy Database (GTDB) (July 2022) and the NCBI RefSeq database (July 2022). Full-length 16 S rRNA gene sequences were retrieved from the WC8-2 genome (WPS_r00030) and the NCBI database (Table S2). Multiple sequences were aligned using SINA (version 1.2.11) [31]. IQ-TREE version 1.6.12 [32] was used to build the phylogeny. ModelFinder [33] was used to determine the optimal evolutionary model for phylogeny building, which selected the TNe+I + G4 model. Branch support was calculated using 1000 ultrafast bootstraps [34]. The pairwise 16 S rRNA gene sequence similarities were determined using SDT software. Phylogenomic analysis based on 400 marker proteins was carried out using PhyloPhlAn v3.0 [35]. Diamond v5.2.32 [36], MAFFT v7.453 [37], and TrimAI were utilized for orthologs searching, multiple sequence alignment within the superphylum Terrabacteria, and gap-trimming, respectively. Gappy sites and sequences with >50% gaps were deleted from the alignments. IQ-TREE version 1.6.12 [32] was used to build the phylogenomic tree. ModelFinder [33] was used to determine the optimal evolutionary model for phylogeny building, which selected the LG + F + R9 model. These analyses were conducted using the “AOBA-B” super- computer (NEC, Tokyo, Japan) with 2CPUs (EPYC7702, AMD, CA, US) and 256GB of RAM. The related similarity of genomes between strain WC8-2 and relatives was estimated using average nucleotide identity (ANI) values, which were calculated using OrthoANI calculator in the EzBio-Cloud web service [38]. The related similarity between strain WC8-2 and its sister phyla with one representative from each class was assessed by pairwise Average Amino acid Identity (AAI) values using the online tools at the Kostas Konstantinidis Lab website Environmental Microbial Genomics Laboratory (http://enve-omics.ce.gatech.edu/aai/). The MAGs used in the tree are listed in Table S3.Phylogeny of photosynthesis- and “atmospheric chemosynthesis”- associated genesWe retrieved phototrophy- and “atmospheric chemosynthesis”-related protein sequences from the WC8-2 genome, “Ca. Eremiobacterota” MAGs, and known phototrophs genomes (Table S3) using the local BLAST server (SequenceServer 1.0.14 [39]) with reference sequences (Table S3) or annotated sequences as queries. Sequences were aligned using MUSCLE and poorly aligned regions were removed using Gblocks version 0.91b [40] or by manual inspection. The alignment sequences were concatenated into a single sequence. The ML tree was constructed using IQ-TREE with the best-hit evolutionary rate model: LG + I + G4 for HhyL, CbbL, BchXYZ and CbbL, LG + F + I + G4 for BchLNB, and LG + F + G4 for PufML, BchI, and BchD. All trees were visualized using iTOL (version 5.0) [41]. All sequences used in the trees are listed in Table S3.Analysis of bacterial communities in the fumarolic ice cavesEnvironmental DNA was extracted from the samples (0.1 g) (WC7, WC8, HD3, MC, Hut, HH, HC) using a DNeasy PowerSoil Kit (Qiagen, Valencia, CA, USA) following the manufacturer’s instructions. The V4 region of the 16 S rRNA gene was PCR amplified using primers with adapter sequences (V3-V4f_MIX) [42]. The PCR cycling was carried out using the following parameters: 94 °C for 2 min followed by 25 cycles at 94 °C for 30 s, 56 °C for 30 s, and 72 °C for 1 min with a final extension at 72 °C for 2 min. Library construction and sequencing were performed at the Bioengineering Lab (Kanagawa, Japan) using MiSeq (Illumina). Briefly, adaptor and primer regions were trimmed using the FASTX-Toolkit v0.0.13 (http://hannonlab.cshl.edu/fastx_toolkit). Read sequences of ≤40 bp with ambiguous bases and low-quality sequences (quality score, ≤Q20), together with their paired-end reads, were filtered out using Sickle v1.33 (https://github.com/najoshi/sickle). High-quality paired-end reads were merged using PEAR v0.9.10 with default settings [43]. Merged sequences of ≤245 and ≥260 bp were discarded using SeqKit v0.8.0 [44]. Operational taxonomic units (OTUs) were classified using QIIME v1.9.1 and the SILVA database (release 132) with 97% identity. To study the phylogeny of the OTUs assigned to the “Ca. Eremiobacterota”, a neighbor-joining (NJ) tree of the OTUs was constructed [45] as described above.Microscopic observationElectron microscopic observations of the cells were performed at Tokai Electron Microscopy (Nagoya, Japan). Briefly, the cells grown in 10% R2A broth (pH6.0) with air/CO2 (90:10, v/v) at 30 °C for 15 days under a 12/12 h light/dark regime with incandescent light (250 μmol m–2s–1) were fixed with 4% paraformaldehyde (PFA) and 4% glutaraldehyde (GA) in 0.1 M phosphate buffer (PB) at pH7.4, and postfixed with 2% OsO4 in 0.1 M PB. Cells were then dehydrated using graded ethanol solutions. The dehydrated cells were polymerized with resin, ultrathin sectioned, stained with 2% uranyl acetate, then secondary-stained with Lead stain solution. A transmission electron microscope (TEM) (JEM-1400Pus; JEOL, Tokyo, Japan) was used to observe the ultrathin sectioned cells at 100 kV acceleration voltage. To observe the negative-stained cells, PFA- and GA-fixed cells were adsorbed on formvar film-coated copper grids and stained with 2% phosphotungstic acid solution (pH 7.0) and observed using a TEM at 100 kV. DAPI (4,6-diamidino2-phenylindole) and Nile-Red staining was performed by incubating 0.1 mL cell suspension with a 1 mL staining solution (1 mg/L DAPI and 1 mg/L Nile Red in PBS buffer) for 10 min. The stained cells were observed under a fluorescence microscope (Olympus AX80T; Olympus Optical; Tokyo, Japan).Growth assayUnless otherwise noted, all cultures were grown in 100 mL butyl stopper- and screw-cap-sealed glass vials containing 50 mL liquid medium (pH6.0) at 30 °C. Growth was monitored by optical density at 600 nm (OD600) using a spectrophotometer (BioSpectrometer Basic; Eppendorf; Tokyo, Japan). Initial cell density was adjusted to 0.005 (OD600). Specific culture conditions are described below. All anaerobic growth tests were conducted with 100% N2 gas in the headspaces and supplemented with a reducing agent (0.3 g/L cysteine-HCl) and a redox indicator (1 mg/L resazurin).Optimal culture conditionsCell growth in different media was examined using 1, 10, 20, 100% (a full strength) R2A broth and Basal_YE with air/CO2 (90:10, v/v) under a 12/12 h light/dark regime with incandescent light (250 μmol m–2s–1) for 25 days. Basal_YE contained (l−1) 0.44 g KH2PO4, 0.1 g (NH4)2SO4, 0.1 g MgSO4.7H2O, 0.3 g yeast extract, and 1 ml trace element SL-8 [refer to DSMZ745]. Cell growth at different temperatures (10, 13, 20, 25, 30, 33 and 37 °C; pH6.0), pH (3.7, 4.5, 6.0, 7.0, 8.0, and 9.0), and NaCl (0, 1, 10, 20, and 30 g [l−1]; pH6.0) was examined using Basal_YE with air/CO2 (90:10, v/v) for 25 days under a 12/12 h light/dark regime with incandescent light (250 μmol m−2s−1). The following pH buffer solutions were used: acetic acid/sodium acetate for pH 4–6, K2HPO4/KH2PO4 for pH 6–8, sodium bicarbonate/sodium carbonate for pH 9–10. To examine colony formation on solid media, WC8-2 was streaked or stabbed on 10% R2AG medium (pH 6.0) and incubated at 30 °C for 30 days under a 100% air atmosphere.Optimal oxygen and carbon dioxide conditionsTo determine the preferred O2 concentration, WC8-2 was grown in Basal_YE (pH 6.0, 30 °C, no NaCl) in butyl stopper-sealed glass bottles with the atmosphere in the headspace adjusted to different N2/O2/CO2 ratios (70:20:10%, 80:10:10%, 89:1:10%, 90:0:10% v/v) after removing oxygen with 100% N2 gas. Cultures were incubated for 25 days under a 12/12 h light/dark regime with incandescent light (250 μmol m–2s–1). To determine the CO2 preference, the gas phase was adjusted to different N2/O2/CO2 ratios (70:20:10%,75:20:5%, 80:20:0% v/v) in sealed bottles or 100% air (plugged with a BIO-SILICO N-38 sponge plug; Shin-Etsu Polymer Co., Ltd, Tokyo, Japan; breathable culture-plug) under the same conditions as the O2 preference test.Photoorganoheterotrophic (or photoorganoautotrophic) and chemoorganoheterotrophic (or chemoorganoautotrophic) growthThe utilization of organic compounds as carbon sources/organic electron donors was tested in Basal medium (Basal_YE without yeast extract, pH 6.0) supplemented with one of the following sources (l-1): 0.3 ml of glycerol, or 0.3 g of sucrose, d-glucose, d-ribose, maltose, l-leucine, l-isoleucine, l-valine, l-serine, l-lysine, taurine, yeast extract or gellan gum, 1 ml of vitamin B12 solution (2 mg/L). Utilization was assessed by measuring growth of the cultures during a 25-day incubation at 30 °C in continuous light (250 μmol m–2s–1) for photoorganoheterotrophy, or continuous dark for chemoorganoheterotrophy under aerobic (air/CO2 [90:10, v/v]) and anaerobic (N2/CO2 [90:10, v/v]) conditions.Photolithoautotrophic and chemolithoautotrophic growthCells were inoculated into amended PSB2 [46] as described below or Basal medium with 5 mM Na2S or Na2S2O3 or 1% H2 (v/v; in the gas phase) as electron donors, and cultivated in continuous light (250 μmol m–2s–1) for photolithoautotrophic growth, or continuous dark for chemolithoautotrophic growth under aerobic (air/CO2 [90:10, v/v]) and anaerobic (N2/CO2 [90:10, v/v]) conditions for 60 days at 30 °C. The amended PBS2 contained (L-1): 0.5 g NH4Cl, 1.0 g KH2PO4, 0.2 g NaCl, 0.4 g MgSO4.7H2O, 0.05 g CaCl2.2H2O, 4.2 g NaHCO3, 1 ml trace element SL-8 [refer to DSMZ745], and 1 ml vitamin B12 solution (2 mg/L), pH6.0.Fermentative or anaerobic growthAnaerobic growth was examined in continuous dark in 20% R2A broth (pH 6.0) with N2/CO2 (90:10, v/v) supplemented with 5 mM Na2SO4, NaNO3, or dimethyl sulfoxide (DMSO) as electron acceptors for 60 days at 30 °C.Pigment assaysCells grown in Basal_YE (pH6.0) with air/CO2 (90:10, v/v) at 30 °C for 14 days (exponential growth phase) under continuous light (250 μmol m–2 s–1) and continuous dark were used for the pigment assays. The absorption spectrum was determined in a cell suspension in 60% (w/v) sucrose and in a 100% methanol extract using a spectrophotometer (V-630; JASCO, Tokyo, Japan) at 350–1100 nm. The BChl a concentration was determined spectroscopically in 100% methanol [47]. Dry cell weight was measured after harvested cells were washed twice with Milli-Q water and dried at 80 °C for 3 days. The extract was also analyzed by HPLC (NEXERA X2; Shimadzu; Kyoto, Japan) equipped with a 4.6 × 250 mm COSMOSIL 5C18-AR (Nakarai Taque; Tokyo, Japan) with isocratic elution of 92.5% (v/v) methanol in water at a flow rate of 1.0 mL/min. BChl a was monitored at 766 nm using a diode-array spectrophotometer detector (SPD-M20A; Shimadzu; Kyoto, Japan).Observation of taxisTo study phototaxis in WC8-2, cells were grown in 20% R2A broth with air/CO2 (90:10, v/v) for 14 days under a 12/12 h light/dark regimen. Cultures were transferred to tissue culture flasks (175 cm2, canted neck, Iwaki, Shizuoka, Japan). Light sources were ultraviolet (UV) at 395 nm (Linkman, Fukui, Japan), blue at 470 nm (CREE, Durham, NC, USA), green at 502 nm (Linkman), red at 653 nm (LENOO, Shinpei, Taiwan), and near-infrared (NIR) at 880 nm (LENOO). The light-emitting device was constructed by assembling LEDs on a breadboard with a power supply. Cultivations were illuminated with each wavelength using a light-emitting device from the underside and incubated at 20 °C for 18 h. Cells aggregating toward a light source was taken to indicate phototaxis. As a control, culture vessels were wrapped in aluminum foil to block light. Images and time-lapse video were captured using an iPhone 6 S camera.Stable carbon isotope ratio mass spectrometry (IRMS)The WC8-2 cells were grown in Basal_YE (pH6.0) with air/13CO2 [90:10, v/v] and air/unlabeled CO2 [90:10, v/v] under continuous light (250 μmol m–2s–1) for photoheterotroph and continuous dark for chemoheterotroph at 30 °C for 14 days (exponential growth phase). Approximately 10 mg culture biomass was collected, washed in HCl overnight, rinsed three times with deionized water, and placed into tin capsules. Stable carbon isotope ratios (δ13C) were analyzed at Shoko Science (Saitama, Japan) using a Delta V Advantage (EA-IRMS; Thermo Fisher Scientific, Bremen, Germany). The standard for C isotope ratio analysis was Vienna PeeDee Belemnite (VPDB). The δ13C values of 13CO2-cultivated cells exceeded the optimum calibration range of the instrument, but were used in this study to provide conclusive evidence that inorganic carbon was incorporated into the biomass.Quantitative reverse transcription PCR (qRT-PCR)Total RNA extraction and cDNA synthesisTotal RNA was extracted from cells grown in Basal_YE (pH6.0) with air/CO2 [90:10, v/v] in continuous light (250 μmol m–2s–1) for photoheterotrophic conditions and in continuous dark for chemoheterotrophic conditions at 30 °C for 14 days (exponential growth phase) using the Total RNA Purification Kit (Norgen, Biotek Corp, Ontario, Canada). DNA was removed from the extracted nucleic acids using an RNase-Free DNase I Kit (Norgen, Biotek Corp) according to the manufacturer’s protocol. The absence of DNA in the RNA samples was confirmed by PCR without reverse transcriptase. cDNA was generated from 500 ng total RNA using a TaKaRa PrimeScript™ 1st strand cDNA Synthesis Kit (TaKaRa Bio) with random hexamers according to the manufacturer’s protocol.Primer design, specificity and efficiencyThe following three photosynthesis- and CO2-fixation-related genes in the WC8-2 genome were selected for qRT-PCR: bchM encoding an enzyme involved in BChl synthesis, pufL encoding the anoxygenic Type II photochemical reaction centers L-subunit, and cbbL encoding the large subunit of type IE RuBisCO. The RNA polymerase subunit beta (rpoB) was used as a housekeeping reference gene. Primers for qRT-PCR were designed with Primer3 (v. 0.4.0) [48] with the following criteria: product size ranging from 80 to 150 bp, optimum Tm of 60 °C and GC content about 50 to 55%. Standard RT-PCR confirmed that each primer set amplified only a single product with expected size (data not shown), and the product was also sequenced using Sanger sequencing at Macrogen Japan Corp to confirm the candidate products. Primer efficiency was calculated for qRT-PCR using the slope of the calibration curve based on a 20-, 40-, 80-, 160-fold dilution series of cDNA samples [49]. In addition, the specificity of the primers was determined by the confirmation of a single peak in the melting curve. All information about the primers is shown in Table S4.qRT-PCRqRT-PCR was performed using a MiniOpticon Real-Time PCR System (Bio-Rad, Marnes la Coquette, France). The reaction mixture contained 10 μL TB Green Premix Ex Taq II (Tli RNaseH Plus, Takara Bio), 0.8 μL 10 mM primer, 2 μL of a 20-, 40-, 80-, 160-fold dilution series of cDNA, and 6.4 μL water. qRT-PCR was performed using the following protocol: denaturation at 95 °C for 30 s; denaturation and amplification at 95 °C for 5 s and 60 °C for 30 s, respectively (40 cycles). Fluorescence was measured at the end of the amplification step, and amplified products were examined by melting curve analysis from 60 to 95 °C. Each reaction was performed in three independent cultivations. Relative gene expression fold change was calculated using the comparative Ct method (2−ΔΔCt) [49]. Normalized expressions were used for reaction in dark. The 2−ΔΔCt values ≤0.5 were defined as downregulated and values ≥2.0 as upregulated, with P  More

  • in

    Building a living shoreline to help combat climate change

    I’m a conservation land manager at the Port of San Diego in California. My team and I aim to manage the tidelands around San Diego Bay, an area of more than 4,850 hectares, three-quarters of which is covered by water at high tide. At least 60% of the bay’s shoreline is ‘hardened’ — that is, it is edged with either a solid seawall or rip rap, piles of artificial boulders.To prevent erosion of the adjacent natural shoreline and restore wetlands, we’re participating in the San Diego Bay Native Oyster Living Shoreline project. As part of that, in December 2021, we placed 360 reef balls — depicted in this photograph from September this year — along 260 metres of shoreline to form the foundation of a native-oyster reef in the south bay. Here, I’m looking for oysters that have settled and are growing on the spheres.The reef balls are made out of ‘baycrete’, a concrete mixture made with local sand and the shells of farmed oysters. These attract wild oysters, which come to live there. We’re targeting the native Olympia oysters (Ostrea lurida), which can filter up to 190 litres of water per day. And sediment should accumulate behind the reef balls, encouraging the growth of eelgrass (Zostera marina). The grass is the foundation of the bay’s food chain.In a couple of years, native oysters will cover the reef balls, forming an artificial reef offshore. This reef will cause storm waves to break farther from the shoreline, protecting the adjacent salt marsh. Just inland from this area is a wetlands habitat refuge for the endangered California least tern (Sternula antillarum browni), and many birds are already hopping onto the reef balls and eating what’s living there.Living shorelines are an important part of sequestering carbon to combat climate change — both eelgrass and oysters store a lot of carbon. The reef balls are win–win–win. I often joke that we’re trying to save the planet one acre (0.4 hectares) at a time. More

  • in

    Ant milk: The mysterious fluid that helps them thrive

    Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain
    the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in
    Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles
    and JavaScript. More

  • in

    Seasonal variation in daily activity patterns of snow leopards and their prey

    Lima, S. L. & Dill, L. M. Behavioral decisions made under the risk of predation: A review and prospectus. Can. J. Zool. 68, 619–640 (1990).Article 

    Google Scholar 
    Ordiz, A., Stoen, O. G., Delibes, M. & Swenson, J. E. Predators or prey? Spatio-temporal discrimination of human-derived risk by brown bears. Oecologia 166, 59–67 (2011).Article 
    ADS 

    Google Scholar 
    Glass, T. W., Breed, G. A., Robards, M. D., Williams, C. T. & Kielland, K. Trade-off between predation risk and behavioural thermoregulation drives resting behaviour in a cold-adapted mesocarnivore. Anim. Behav. 175, 163–174 (2021).Article 

    Google Scholar 
    Daan, S. & Aschoff, J. Circadian rhytms of locomotor activity in captive birds and mammals: Their variation with seasons and latitude. Oecologia 18, 269–316 (1975).Article 
    ADS 

    Google Scholar 
    Kronfeld-Schor, N. & Dayan, T. Partitioning of time as an ecological resource. Annu. Rev. Ecol. Evol. Syst. 34, 153–181 (2003).Article 

    Google Scholar 
    Garcia, R. A., Cabeza, M., Rahbek, C. & Araujo, M. B. Multiple dimensions of climate change and their implications for biodiversity. Science 344, 1247579 (2014).Article 

    Google Scholar 
    Curio, E. The Ethology of Predation (Springer-Verlag, 1976).Book 

    Google Scholar 
    Linkie, M. & Ridout, M. S. Assessing tiger-prey interactions in Sumatran rainforests. J. Zool. 284, 224–229 (2011).Article 

    Google Scholar 
    Heurich, M. et al. Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range. PLoS ONE 9, e114143 (2014).Article 
    ADS 

    Google Scholar 
    Harmsen, B. J., Foster, R. J., Silver, S. C., Ostro, L. E. T. & Doncaster, C. P. Jaguar and puma activity patterns in relation to their main prey. Mamm. Biol. 76, 320–324 (2011).Article 

    Google Scholar 
    Foster, V. C. et al. Jaguar and puma activity patterns and predator-prey interactions in four Brazilian biomes. Biotropica 45, 373–379 (2013).Article 

    Google Scholar 
    Theuerkauf, J. et al. Daily patterns and duration of wolf activity in the Bialowieza forest, Poland. J. Mammal. 84, 243–253 (2003).Article 

    Google Scholar 
    Hebblewhite, M., Merrill, E. H. & McDonald, T. L. Spatial decomposition of predation risk usign resource selection functions: An example in a wolf-elk predator-prey system. Oikos 111, 101–111 (2005).Article 

    Google Scholar 
    Balme, G., Hunter, L. & Slotow, R. Feeding habitat selection by hunting leopards Panthera pardus in a woodland Savanna: Prey catchability versus abundance. Anim. Behav. 74, 589–598 (2007).Article 

    Google Scholar 
    Smith, J. A. et al. Where and when to hunt? Decomposing predation success of an ambush carnivore. Ecology 101, e03172 (2020).Article 

    Google Scholar 
    Hopcraft, J. G. C., Sinclair, A. R. E. & Packer, C. Planning for success: Serengeti lions seek prey accessibility rather than abundance. J. Anim. Ecol. 74, 559–566 (2005).Article 

    Google Scholar 
    Theuerkauf, J. What drives wolves: Fear or hunger? Humans, diet, climate and wolf activity patterns. Ethology 115, 649–657 (2009).Article 

    Google Scholar 
    Funston, P. J., Mills, M. G. & Biggs, H. C. Factors affecting the hunting success of male and female lions in the Kruger National Park. J. Zool. 253, 419–431 (2001).Article 

    Google Scholar 
    Schaller, G. The Serengeti lion (The University of Chicago Press, IL, 1972).
    Google Scholar 
    Bailey, T. N. The African Leopard, Ecology and Behaviour of a Solitary Felid (The Blackburn Press, 1993).Book 

    Google Scholar 
    Jenny, D. & Zuberbühler, K. Hunting behaviour in West African forest leopards. Afr. J. Ecol. 43, 197–200 (2005).Article 

    Google Scholar 
    Packer, C., Swanson, A., Ikanda, D. & Kushnir, H. Fear of darkness, the full moon and the nocturnal ecology of African lions. PLoS ONE 6, e22285 (2011).Article 
    ADS 
    CAS 

    Google Scholar 
    Palmer, M. S., Fieberg, J., Swanson, A., Kosmala, M. & Packer, C. A “dynamic” landscape of fear: Prey responses to spatiotemporal variations in predation risk across the lunar cycle. Ecol. Letters 20, 1364–1373 (2017).Article 
    CAS 

    Google Scholar 
    Steinmetz, R., Seuaturien, N. & Chutipong, W. Tigers, leopards, and dholes in a half-empty forest: Assessing species interactions in a guild of threatened carnivores. Biol. Cons. 163, 68–78 (2013).Article 

    Google Scholar 
    Carter, N., Jasny, M., Gurung, B. & Liu, J. Impacts of people and tigers on leopard spatiotemporal activity patterns in a global biodiversity hotspot. Global Ecol. Conserv. 3, 149–162 (2015).Article 

    Google Scholar 
    George, S. L. & Crooks, K. R. Recreation and large mammal activity in an urban nature reserve. Biol. Cons. 133, 107–117 (2006).Article 

    Google Scholar 
    Beltrán, J. F. & Delibes, M. Environmental determinants of circadian activity of free-ranging Iberian lynxes. J. Mammal. 75, 382–393 (1994).Article 

    Google Scholar 
    McNab, B. K. The standard energetics of mammalian carnivores: Felidae and Hyaenidae. Sikes Can. J. Zool. 78, 2227–2239 (2000).Article 

    Google Scholar 
    Mishra, C. et al. Increasing risks for emerging infectious diseases within a rapidly changing High Asia. Ambio 51, 494–507 (2022).Article 

    Google Scholar 
    Mishra, C., Redpath, S. M. & Suryawanshi, K. R. Livestock predation by snow leopards: Conflicts and the search for solutions. In Snow Leopards (eds McCarthy, T. M. & Mallon, D.) 59–67 (Academic Press, 2016).Chapter 

    Google Scholar 
    Farrington, J. D., and J. Li. 2016. Climate change impacts on snow leopard range. In: McCarthy, T.M., Mallon, D., editors. Snow Leopards. Academic Press.Jackson, R. Home Range, Movements and Habitat use of Snow Leopard in Nepal (Dissertation niversity of London, London, 1996).
    Google Scholar 
    McCarthy, T. M., Fuller, T. K. & Munkhtsog, B. Movements and activities of snow leopards in Southwestern Mongolia. Biol. Cons. 124, 527–537 (2005).Article 

    Google Scholar 
    Salvatori, M. et al. Co-occurrence of snow leopard, wolf and Siberian ibex under livestock encroachment into protected areas across the Mongolian Altai. Biol. Cons. 261, 109294 (2021).Article 

    Google Scholar 
    Rode, J. et al. Population monitoring of snow leopards using camera trapping in Naryn state nature reserve, Kyrgyzstan, between 2016 and 2019. Global Ecol. Conserv. 31, e01850 (2021).Article 

    Google Scholar 
    Sharma, R. K. et al. Spatial variation in population-density of snow leopards in a multiple use landscape in Spiti Valley Trans-Himalay. PLoS ONE 16, e0250900 (2021).Article 
    CAS 

    Google Scholar 
    Kachel, S. M., Karimov, K. & Wirsing, A. J. Predator niche overlap and partitioning and potential interactions in the mountains of Central Asia. J. Mammal. 103, 1019–1029 (2022).Article 

    Google Scholar 
    Johansson, Ö., Simms, A. & McCarthy, T. M. From VHF to satellite GPS collars: Advancements in snow leopard telemetry. In Snow leopards (eds McCarthy, T. M. & Mallon, D.) p355-365 (Academic Press, 2016).Chapter 

    Google Scholar 
    Johansson, Ö. et al. Snow leopard predation in a livestock dominated landscape in Mongolia. Biol. Cons. 184, 251–258 (2015).Article 

    Google Scholar 
    Havmøller, R. W., Jacobsen, N. S., Scharff, N., Rovero, F. & Zimmermann, F. Assessing the activity pattern overlap among leopards (Panthera pardus), potential prey and competitors in a complex landscape in Tanzania. J. Zool. 311, 175–182 (2020).Article 

    Google Scholar 
    Kitchener, A. C., Van Valkenburgh, B. & Yamaguchi, N. Felid form and function. In Biology and Conservation of Wild Felids (eds MacDonald, D. W. & Loveridge, A. J.) 83–106 (Oxford University Press, 2010).
    Google Scholar 
    Fuglesteg, B. N., Haga, Ø. E., Folkow, L. P., Fuglei, E. & Blix, A. S. Seasonal variations in basal metabolic rate, lower critical temperature and responses to temporary starvation in the arctic fox (Alopex lagopus) from Svalbard. Polar Biol. 29, 308–319 (2005).Article 

    Google Scholar 
    Doris, P. A. & Baker, M. A. Effect of dehydration on thermoregulation in cats exposed to high ambient temperatures. J. Appl. Physiol. 51, 46–54 (1981).Article 
    CAS 

    Google Scholar 
    Forrest, J. L. et al. Conservation and climate change: Assessing the vulnerability of snow leopard habitat to treeline shift in the Himalaya. Biol. Cons. 150, 129–135 (2012).Article 

    Google Scholar 
    Sharma, R. K., Bhatnagar, Y. V. & Mishra, C. Does livestock benefit or harm snow leopards?. Biol. Cons. 190, 8–13 (2015).Article 

    Google Scholar 
    Samelius, G. et al. Keeping predators out: Testing fences to reduce livestock depredation at night-time corrals. Oryx 55, 466–472 (2021).Article 

    Google Scholar 
    Hebblewhite, M. & Merrill, E. Modelling wildlife-human relationships for social species with mixed-effects resource selection models. J. Appl. Ecol. 45, 834–844 (2007).Article 

    Google Scholar 
    Johansson, Ö., Malmsten, J., Mishra, C., Lkhagvajav, P. & McCarthy, T. Reversible immobilization of free-ranging snow leopards (Panthera uncia) with a combination of medetomidine and tiletamine-zolazepam. J. Wildl. Dis. 49, 338–346 (2013).Article 

    Google Scholar 
    Johansson, Ö., Kachel, S. & Weckworth, B. Guidelines for telemetry studies on snow leopards. Animals 12, 1663 (2022).Article 

    Google Scholar 
    Bjørneraas, K., Van Moorter, B., Rolandsen, C. M. & Herfindal, I. Screening global positioning system location data for errors using animal movement characteristics. J. Wildl. Manag. 74, 1361–1366 (2010).Article 

    Google Scholar 
    Pålsson O. 2022. Maternal behaviour of the snow leopard (Panthera uncial). MSc thesis. Uppsala University, Uppsala; Sweden https://www.diva-portal.org/smash/get/diva2:1668965/FULLTEXT01.pdf.du Sert, N. P. et al. Reporting animal research: Explanation and elaboration for the ARRIVE guidelines 20. PLoS Biol. 18, e3000411 (2020).Article 

    Google Scholar 
    Nygren, E. 2015. Activity patterns of snow leopards (Panthera uncia) at their kill sites. MSc thesis, Swedish University of Agricultural Sciences, Uppsala, Sweden. https://stud.epsilon.slu.se/8109/1/nygren_e_150625.pdf.Johansson, Ö. et al. Land sharing is essential for snow leopard conservation. Biol. Cons. 203, 1–7 (2016).Article 

    Google Scholar 
    Johansson, Ö. et al. The timing of breeding and independence for snow leopard females and their cubs. Mamm. Biol. 101, 173–180 (2021).Article 

    Google Scholar 
    Nouvellet, P., Rasmussen, G. S. A., Macdonald, D. W., Courchamp, F. & Braae, A. Noisy clocks and silent sunrises: Measurement methods of daily activity pattern. J. Zool. 286, 179–184 (2012).Article 

    Google Scholar 
    Ridout, M. S. & Linkie, M. Estimating overlap of daily activity patterns from camera trap data. J. Agric. Biol. Environ. Stat. 14, 322–337 (2009).Article 
    MathSciNet 
    MATH 

    Google Scholar 
    R Development core team. 2019. R: A language and environment for statistical computing. R foundation for statistical computing Vienna, Austria. www.R-project.org/.Wood, S. N. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J. Roy. Stat. Soc. B 73, 3–36 (2011).Article 
    MathSciNet 
    MATH 

    Google Scholar  More

  • in

    Infection with an acanthocephalan helminth reduces anxiety-like behaviour in crustacean host

    Cézilly, F. & Perrot-Minnot, M. J. Interpreting multidimensionality in parasite-induced phenotypic alterations: Panselectionism versus parsimony. Oikos 119, 1224–1229 (2010).Article 

    Google Scholar 
    Moore, J. Parasites and the Behavior of Animals. (Oxford University Press on Demand, 2002).
    Google Scholar 
    Thomas, F. et al. Do hairworms (Nematomorpha) manipulate the water seeking behaviour of their terrestrial hosts?. J. Evol. Biol. 15, 356–361 (2002).Article 

    Google Scholar 
    Weinersmith, K. L. What’s gotten into you? A review of recent research on parasitoid manipulation of host behavior. Curr. Opin. Insect Sci. 33, 37–42 (2019).Article 

    Google Scholar 
    Hughes, D. P. et al. Behavioral mechanisms and morphological symptoms of zombie ants dying from fungal infection. BMC. Ecol. 11, (2011).Lagrue, C., Kaldonski, N., Perrot-Minnot, M. J., Motreuil, S. & Bollache, L. Modification of hosts’ behavior by a parasite: Field evidence for adaptive manipulation. Ecology 88, 2839–2847 (2007).Article 

    Google Scholar 
    Berdoy, M., Webster, J. P. & Mcdonald, D. W. Fatal attraction in rats infected with Toxoplasma gondii. Proc. R. Soc. B Biol. Sci. 267, 1591–1594 (2000).Article 
    CAS 

    Google Scholar 
    Perrot-Minnot, M. J., Kaldonski, N. & Cézilly, F. Increased susceptibility to predation and altered anti-predator behaviour in an acanthocephalan-infected amphipod. Int. J. Parasitol. 37, 645–651 (2007).Article 

    Google Scholar 
    Cézilly, F. & Perrot-Minnot, M. J. Studying adaptive changes in the behaviour of infected hosts: A long and winding road. Behav. Proc. 68, 223–228 (2005).Article 

    Google Scholar 
    Seppälä, O. & Jokela, J. Host manipulation as a parasite transmission strategy when manipulation is exploited by non-host predators. Biol. Lett. 4, 663–666 (2008).Article 

    Google Scholar 
    Dianne, L. et al. Protection first then facilitation: A manipulative parasite modulates the vulnerability to predation of its intermediate host according to its own developmental stage. Evolution 65, 2692–2698 (2011).Article 

    Google Scholar 
    Iritani, R. & Sato, T. Host-manipulation by trophically transmitted parasites: The switcher-paradigm. Trends Parasitol. 34, 934–944 (2018).Article 

    Google Scholar 
    Poulin, R. & Maure, F. Host manipulation by parasites: A look back before moving forward. Trends Parasitol. 31, 563–570 (2015).Article 

    Google Scholar 
    Herbison, R., Lagrue, C. & Poulin, R. The missing link in parasite manipulation of host behaviour. Parasite Vectors 11, 1–6 (2018).Article 

    Google Scholar 
    Perrot-Minnot, M. J. & Cézilly, F. Investigating candidate neuromodulatory systems underlying parasitic manipulation: Concepts, limitations and prospects. J. Exp. Biol. 216, 134–141 (2013).Article 

    Google Scholar 
    Adamo, S. A. Parasites: Evolution’s neurobiologists. J. Exp. Biol. 216, 3–10 (2013).Article 
    CAS 

    Google Scholar 
    Kaushik, M., Lamberton, P. H. L. & Webster, J. P. The role of parasites and pathogens in influencing generalised anxiety and predation-related fear in the mammalian central nervous system. Horm. Behav. 62, 191–201 (2012).Article 

    Google Scholar 
    Grupe, D. W. & Nitschke, J. B. Uncertainty and anticipation in anxiety: An integrated neurobiological and psychological perspective. Nat. Rev. Neurosci. 14, 488–501 (2013).Article 
    CAS 

    Google Scholar 
    Perry, C. J. & Baciadonna, L. Studying emotion in invertebrates: What has been done, what can be measured and what they can provide. J. Exp. Biol. 220, 3856–3868 (2017).Article 

    Google Scholar 
    Adamec, R. E., Burton, P., Shallow, T. & Budgell, J. NMDA receptors mediate lasting increases in anxiety-like behavior produced by the stress of predator exposure—Implications for anxiety associated with posttraumatic stress disorder. Physiol. Behav. 65, 723–737 (1998).Article 

    Google Scholar 
    Bacqué-Cazenave, J. et al. Serotonin in animal cognition and behavior. Int. J. Mol. Sci. 21, 1–23 (2020).Article 

    Google Scholar 
    Hamilton, T. J., Kwan, G. T., Gallup, J. & Tresguerres, M. Acute fluoxetine exposure alters crab anxiety-like behaviour, but not aggressiveness. Sci. Rep. 6, 4–9 (2016).Article 

    Google Scholar 
    de Bekker, C. et al. Gene expression during zombie ant biting behavior reflects the complexity underlying fungal parasitic behavioral manipulation. BMC Genom. 16, 1–23 (2015).Article 

    Google Scholar 
    Shaw, J. C. et al. Parasite manipulation of brain monoamines in California killifish (Fundulus parvipinnis) by the trematode Euhaplorchis californiensis. Proc. R. Soc. B Biol. Sci. 276, 1137–1146 (2009).Article 
    CAS 

    Google Scholar 
    Fayard, M., Dechaume-Moncharmont, F. X., Wattier, R. & Perrot-Minnot, M. J. Magnitude and direction of parasite-induced phenotypic alterations: A meta-analysis in acanthocephalans. Biol. Rev. 95, 1233–1251 (2020).Article 

    Google Scholar 
    Tain, L., Perrot-Minnot, M. J. & Cézilly, F. Altered host behaviour and brain serotonergic activity caused by acanthocephalans: Evidence for specificity. Proc. R. Soc. B Biol. Sci. 273, 3039–3045 (2006).Article 
    CAS 

    Google Scholar 
    Perrot-Minnot, M. J., Maddaleno, M., Balourdet, A. & Cézilly, F. Host manipulation revisited: No evidence for a causal link between altered photophobia and increased trophic transmission of amphipods infected with acanthocephalans. Funct. Ecol. 26, 1007–1014 (2012).Article 

    Google Scholar 
    Perrot-Minnot, M. J., Sanchez-Thirion, K. & Cézilly, F. Multidimensionality in host manipulation mimicked by serotonin injection. Proc. R. Soc. B Biol. Sci. 281, (2014).Perrot-Minnot, M. J., Banchetry, L. & Cézilly, F. Anxiety-like behaviour increases safety from fish predation in an amphipod crustacea. R. Soc. Open Sci. 4, (2017).Perrot-Minnot, M. J., Balourdet, A. & Musset, O. Optimization of anesthetic procedure in crustaceans: Evidence for sedative and analgesic-like effect of MS-222 using a semi-automated device for exposure to noxious stimulus. Aquat. Toxicol. 240, 105981 (2021).Article 
    CAS 

    Google Scholar 
    Barr, S., Laming, P. R., Dick, J. T. A. & Elwood, R. W. Nociception or pain in a decapod crustacean?. Anim. Behav. 75, 745–751 (2008).Article 

    Google Scholar 
    Fossat, P., Bacqué-Cazenave, J., de Deurwaerdère, P., Delbecque, J. P. & Cattaert, D. Anxiety-like behavior in crayfish is controlled by serotonin. Science 1979(344), 1293–1297 (2014).Article 
    ADS 

    Google Scholar 
    Magee, B. & Elwood, R. W. Shock avoidance by discrimination learning in the shore crab (Carcinus maenas) is consistent with a key criterion for pain. J. Exp. Biol. 216, 353–358 (2013).Article 

    Google Scholar 
    Rakitin, A., Tomsic, D. & Maldonado, H. Habituation and sensitization to an electrical shock in the crab Chasmagnathus. Effect of background illumination. Physiol. Behav. 50, 477–487 (1991).Article 
    CAS 

    Google Scholar 
    Koolhaas, J. M. et al. Stress revisited: A critical evaluation of the stress concept. Neurosci. Biobehav. Rev. 35, 1291–1301 (2011).Article 
    CAS 

    Google Scholar 
    Yuan, T. F. & Hou, G. The effects of stress on glutamatergic transmission in the brain. Mol. Neurobiol. 51, 1139–1143 (2015).Article 
    CAS 

    Google Scholar 
    Fossat, P., Bacqué-Cazenave, J., de Deurwaerdère, P., Cattaert, D. & Delbecque, J. P. Serotonin, but not dopamine, controls the stress response and anxiety-like behavior in the crayfish Procambarus clarkii. J. Exp. Biol. 218, 2745–2752 (2015).
    Google Scholar 
    Benesh, D. P., Valtonen, E. T. & Seppälä, O. Multidimensionality and intra-individual variation in host manipulation by an acanthocephalan. Parasitology 135, 617–626 (2008).Article 
    CAS 

    Google Scholar 
    Kaldonski, N., Perrot-Minnot, M. J. & Cézilly, F. Differential influence of two acanthocephalan parasites on the antipredator behaviour of their common intermediate host. Anim. Behav. 74, 1311–1317 (2007).Article 

    Google Scholar 
    Kaldonski, N., Perrot-Minnot, M. J., Motreuil, S. & Cézilly, F. Infection with acanthocephalans increases the vulnerability of Gammarus pulex (Crustacea, Amphipoda) to non-host invertebrate predators. Parasitology 135, 627–632 (2008).Article 
    CAS 

    Google Scholar 
    Parker, G. A., Ball, M. A., Chubb, J. C., Hammerschmidt, K. & Milinski, M. When should a trophically transmitted parasite manipulate its host?. Evolution 63, 448–458 (2009).Article 

    Google Scholar 
    Paul, E. S. & Mendl, M. T. Animal emotion: Descriptive and prescriptive definitions and their implications for a comparative perspective. Appl. Anim. Behav. Sci. 205, 202–209 (2018).Article 

    Google Scholar 
    Anderson, D. J. & Adolphs, R. A framework for studying emotions across species. Cell 157, 187–200 (2014).Article 
    CAS 

    Google Scholar 
    Weinberger, J. & Klaper, R. Environmental concentrations of the selective serotonin reuptake inhibitor fluoxetine impact specific behaviors involved in reproduction, feeding and predator avoidance in the fish Pimephales promelas (fathead minnow). Aquat. Toxicol. 151, 77–83 (2014).Article 
    CAS 

    Google Scholar 
    Curran, K. P. & Chalasani, S. H. Serotonin circuits and anxiety: What can invertebrates teach us?. Invertebr. Neurosci. 12, 81–92 (2012).CAS 

    Google Scholar 
    Mohammad, F. et al. Ancient anxiety pathways influence Drosophila defense behaviors. Curr. Biol. 26, 981–986 (2016).Article 
    CAS 

    Google Scholar 
    Kavaliers, M. & Colwell, D. D. Decreased predator avoidance in parasitized mice: neuromodulatory correlates. Parasitology 111, 257–263 (1995).Article 

    Google Scholar 
    Chivers, D. P., Wisenden, B. D. & Smith, R. J. F. Damselfly larvae learn to recognize predators from chemical cues in the predator’s diet. Anim. Behav. 52, 315–320 (1996).Article 

    Google Scholar 
    Hazlett, B. A., Acquistapace, P. & Gherardi, F. Differences in memory capabilities in invasive and native crayfish. J. Crustac. Biol. 22, 439–448 (2002).Article 

    Google Scholar 
    R Core Team. R Core Team (2014). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/. R Foundation for Statistical Computing (2014). More