More stories

  • in

    Sedimentary ancient DNA reveals a threat of warming-induced alpine habitat loss to Tibetan Plateau plant diversity

    1.Körner, C. Alpine Plant Life: Functional Plant Ecology of High Mountain Ecosystems; with 47 Tables. (Springer Science & Business Media, 2003).2.Alexander, J. M. et al. Plant invasions into mountains and alpine ecosystems: current status and future challenges. Alp. Bot. 126, 89–103 (2016).Article 

    Google Scholar 
    3.Xu, J. et al. The Melting Himalayas: cascading effects of climate change on water, biodiversity, and livelihoods. Conserv. Biol. 23, 520–530 (2009).CAS 
    PubMed 
    Article 

    Google Scholar 
    4.Liang, Q. et al. Shifts in plant distributions in response to climate warming in a biodiversity hotspot, the Hengduan Mountains. J. Biogeogr. 45, 1334–1344 (2018).Article 

    Google Scholar 
    5.Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    6.Dullinger, S. et al. Extinction debt of high-mountain plants under twenty-first-century climate change. Nat. Clim. Change 2, 619–622 (2012).ADS 
    Article 

    Google Scholar 
    7.Zhang, D. C., Zhang, Y. H., Boufford, D. E. & Sun, H. Elevational patterns of species richness and endemism for some important taxa in the Hengduan Mountains, southwestern China. Biodivers. Conserv. 18, 699–716 (2009).Article 

    Google Scholar 
    8.Tang, Z., Wang, Z., Zheng, C. & Fang, J. Biodiversity in China’s mountains. Front. Ecol. Environ. 4, 347–352 (2006).Article 

    Google Scholar 
    9.Lomolino, Mark V. Elevation gradients of species-density: historical and prospective views. Glob. Ecol. Biogeogr. 10, 3–13 (2001).Article 

    Google Scholar 
    10.Colwell, RobertK. & Lees, D. C. The mid-domain effect: geometric constraints on the geography of species richness. Trends Ecol. Evol. 15, 70–76 (2000).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    11.Dirnböck, T., Essl, F. & Rabitsch, W. Disproportional risk for habitat loss of high-altitude endemic species under climate change. Glob. Change Biol. 17, 990–996 (2011).ADS 
    Article 

    Google Scholar 
    12.Su, X., Han, W., Liu, G., Zhang, Y. & Lu, H. Substantial gaps between the protection of biodiversity hotspots in alpine grasslands and the effectiveness of protected areas on the Qinghai-Tibetan Plateau, China. Agric. Ecosyst. Environ. 278, 15–23 (2019).Article 

    Google Scholar 
    13.Zhang, Y. et al. Spatial phylogenetics of two topographic extremes of the Hengduan Mountains in southwestern China and its implications for biodiversity conservation. Plant Divers. (2020) https://doi.org/10.1016/j.pld.2020.09.001.14.Hopping, K. A., Knapp, A. K., Dorji, T. & Klein, J. A. Warming and land use change concurrently erode ecosystem services in Tibet. Glob. Change Biol. 24, 5534–5548 (2018).ADS 
    Article 

    Google Scholar 
    15.Trivedi, M. R., Berry, P. M., Morecroft, M. D. & Dawson, T. P. Spatial scale affects bioclimate model projections of climate change impacts on mountain plants. Glob. Change Biol. 14, 1089–1103 (2008).ADS 
    Article 

    Google Scholar 
    16.Gavin, D. G. et al. Climate refugia: joint inference from fossil records, species distribution models and phylogeography. N. Phytol. 204, 37–54 (2014).Article 

    Google Scholar 
    17.Birks, H. J. B. et al. Does pollen-assemblage richness reflect floristic richness? A review of recent developments and future challenges. Rev. Palaeobot. Palynol. 228, 1–25 (2016).Article 

    Google Scholar 
    18.Parducci, L. et al. Ancient plant DNA in lake sediments. N. Phytol. 214, 924–942 (2017).CAS 
    Article 

    Google Scholar 
    19.Kramer, A., Herzschuh, U., Mischke, S. & Zhang, C. Late glacial vegetation and climate oscillations on the southeastern Tibetan Plateau inferred from the Lake Naleng pollen profile. Quat. Res. 73, 324–335 (2010).CAS 
    Article 

    Google Scholar 
    20.Kramer, A., Herzschuh, U., Mischke, S. & Zhang, C. Holocene treeline shifts and monsoon variability in the Hengduan Mountains (southeastern Tibetan Plateau), implications from palynological investigations. Palaeogeogr. Palaeoclimatol. Palaeoecol. 286, 23–41 (2010).Article 

    Google Scholar 
    21.Hou, G., Yang, P., Cao, G., Chongyi, E. & Wang, Q. Vegetation evolution and human expansion on the Qinghai–Tibet Plateau since the Last Deglaciation. Quat. Int. 430, 82–93 (2017).Article 

    Google Scholar 
    22.Chen, F. et al. Climate change, vegetation history, and landscape responses on the Tibetan Plateau during the Holocene: a comprehensive review. Quat. Sci. Rev. 243, 106444 (2020).Article 

    Google Scholar 
    23.Singh, U. M., Gupta, V., Rao, V. P., Sengar, R. S. & Yadav, M. K. A review on biological activities and conservation of endangered medicinal herb Nardostachys jatamansi. Int. J. Med. Arom. Plants 3, 113–124 (2013).CAS 

    Google Scholar 
    24.Li, X. H., Zhu, X. X., Niu, Y. & Sun, H. Phylogenetic clustering and overdispersion for alpine plants along elevational gradient in the Hengduan Mountains Region, southwest China: Phylogenetic structure along elevational gradient. J. Syst. Evol. 52, 280–288 (2014).Article 

    Google Scholar 
    25.Kanz, C. et al. The EMBL nucleotide sequence database. Nucleic Acids Res. 33, D29–D33 (2005).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    26.Yu, H. et al. Contrasting floristic diversity of the Hengduan mountains, the Himalayas and the Qinghai-Tibet Plateau Sensu Stricto in China. Front. Ecol. Evol. 8 (2020).27.Scheiner, S. M. et al. The underpinnings of the relationship of species richness with space and time. Ecol. Monogr. 81, 195–213 (2011).Article 

    Google Scholar 
    28.Milner, A. M. et al. Glacier shrinkage driving global changes in downstream systems. Proc. Natl Acad. Sci. USA 114, 9770–9778 (2017).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    29.Opitz, S., Zhang, C., Herzschuh, U. & Mischke, S. Climate variability on the south-eastern Tibetan Plateau since the Lateglacial based on a multiproxy approach from Lake Naleng – comparing pollen and non-pollen signals. Quat. Sci. Rev. 115, 112–122 (2015).ADS 
    Article 

    Google Scholar 
    30.Laliberté, E. et al. How does pedogenesis drive plant diversity? Trends Ecol. Evol. 28, 331–340 (2013).PubMed 
    Article 

    Google Scholar 
    31.Nichols, R. V. et al. Minimizing polymerase biases in metabarcoding. Mol. Ecol. Resour. 18, 927–939 (2018).CAS 
    Article 

    Google Scholar 
    32.Alsos, I. G. et al. Plant DNA metabarcoding of lake sediments: How does it represent the contemporary vegetation. PLoS ONE 13, e0195403 (2018).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    33.Niemeyer, B., Epp, L. S., Stoof-Leichsenring, K. R., Pestryakova, L. A. & Herzschuh, U. A comparison of sedimentary DNA and pollen from lake sediments in recording vegetation composition at the Siberian treeline. Mol. Ecol. Resour. 17, e46–e62 (2017).CAS 
    PubMed 
    Article 

    Google Scholar 
    34.Stein, A., Gerstner, K. & Kreft, H. Environmental heterogeneity as a universal driver of species richness across taxa, biomes and spatial scales. Ecol. Lett. 17, 866–880 (2014).PubMed 
    Article 

    Google Scholar 
    35.Chen, F. H. et al. Agriculture facilitated permanent human occupation of the Tibetan Plateau after 3600 B.P. Science 347, 248–250 (2015).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    36.Group, M. R. I. E. W. et al. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Change 5, 424–430 (2015).ADS 
    Article 

    Google Scholar 
    37.Liang, E. et al. Species interactions slow warming-induced upward shifts of treelines on the Tibetan Plateau. Proc. Natl Acad. Sci. USA 113, 4380–4385 (2016).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    38.Alexander, J. M. et al. Lags in the response of mountain plant communities to climate change. Glob. Change Biol. 24, 563–579 (2018).ADS 
    Article 

    Google Scholar 
    39.Zu, K. et al. Altitudinal biodiversity patterns of seed plants along Gongga Mountain in the southeastern Qinghai–Tibetan Plateau. Ecol. Evol. 9, 9586–9596 (2019).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    40.Sun, H., Zhang, J., Deng, T. & Boufford, D. E. Origins and evolution of plant diversity in the Hengduan Mountains, China. Plant Divers. 39, 161–166 (2017).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    41.Wiens, J. A., Stralberg, D., Jongsomjit, D., Howell, C. A. & Snyder, M. A. Niches, models, and climate change: assessing the assumptions and uncertainties. Proc. Natl Acad. Sci. 106, 19729–19736 (2009).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    42.Kramer, A., Herzschuh, U., Mischke, S. & Zhang, C. Late Quaternary environmental history of the south-eastern Tibetan Plateau inferred from the Lake Naleng non-pollen palynomorph record. Veg. Hist. Archaeobotany 19, 453–468 (2010).Article 

    Google Scholar 
    43.Stuiver, M. & Reimer, P. J. Extended 14C Data Base and Revised CALIB 3.0 14C Age Calibration Program. Radiocarbon 35, 215–230 (1993).Article 

    Google Scholar 
    44.Intcal04 Terrestrial Radiocarbon Age Calibration, 0–26 Cal Kyr BP. Radiocarbon 46, 1029–1058 (2004).45.Taberlet, P. et al. Power and limitations of the chloroplast trnL (UAA) intron for plant DNA barcoding. Nucleic Acids Res. 35, e14–e14 (2007).PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    46.Coissac, E. OligoTag: A Program for Designing Sets of Tags for Next-Generation Sequencing of Multiplexed Samples. in Data Production and Analysis in Population Genomics: Methods and Protocols (eds. Pompanon, F. & Bonin, A.) 13–31 (Humana Press, 2012). https://doi.org/10.1007/978-1-61779-870-2_2.47.De Barba, M. et al. DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet. Mol. Ecol. Resour. 14, 306–323 (2014).PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    48.Boyer, F. et al. obitools: a unix-inspired software package for DNA metabarcoding. Mol. Ecol. Resour. 16, 176–182 (2016).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    49.SøNstebø, J. H. et al. Using next-generation sequencing for molecular reconstruction of past Arctic vegetation and climate: TECHNICAL ADVANCES. Mol. Ecol. Resour. 10, 1009–1018 (2010).PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    50.Willerslev, E. et al. Fifty thousand years of Arctic vegetation and megafaunal diet. Nature 506, 47–51 (2014).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    51.Soininen, E. M. et al. Highly overlapping winter diet in two sympatric lemming species revealed by DNA metabarcoding. PLoS ONE 10, e0115335 (2015).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    52.Ficetola, G. F. et al. Replication levels, false presences and the estimation of the presence/absence from eDNA metabarcoding data. Mol. Ecol. Resour. 15, 543–556 (2015).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    53.R Core Team. R: A Language and Environment for Statistical Computing. (2019).54.Brach, A. R. & Song, H. eFloras: New directions for online floras exemplified by the Flora of China Project. TAXON 55, 188–192 (2006).Article 

    Google Scholar 
    55.Zhao, Y. et al. Evolution of vegetation and climate variability on the Tibetan Plateau over the past 1.74 million years. Sci. Adv. 6, eaay6193 (2020).56.Herzschuh, U. et al. Position and orientation of the westerly jet determined Holocene rainfall patterns in China. Nat. Commun. 10, 2376 (2019).ADS 
    PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    57.Shakun, J. D. et al. Global warming preceded by increasing carbon dioxide concentrations during the last deglaciation. Nature 484, 49–54 (2012).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    58.Marcott, S. A., Shakun, J. D., Clark, P. U. & Mix, A. C. A reconstruction of regional and global temperature for the past 11,300 years. Science 339, 1198–1201 (2013).ADS 
    CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    59.Kessler, M. A., Anderson, R. S. & Stock, G. M. Modeling topographic and climatic control of east-west asymmetry in Sierra Nevada glacier length during the Last Glacial Maximum. J. Geophys. Res. Earth Surf. 111, F02002 (2006).60.Pfeffer, W. T. et al. The Randolph Glacier Inventory: a globally complete inventory of glaciers. J. Glaciol. 60, 537–552 (2014).ADS 
    Article 

    Google Scholar 
    61.Braithwaite, R. J. From Doktor Kurowski’s Schneegrenze to our modern glacier equilibrium line altitude (ELA). Cryosphere 9, 2135–2148 (2015).ADS 
    Article 

    Google Scholar 
    62.Maussion, F. et al. Precipitation seasonality and variability over the Tibetan Plateau as resolved by the high asia reanalysis. J. Clim. 27, 1910–1927 (2014).ADS 
    Article 

    Google Scholar 
    63.Anja, M.-C. et al. GPCC Climatology Version 2011 at 0.25°: Monthly Land-Surface Precipitation Climatology for Every Month and the Total Year from Rain-Gauges built on GTS-based and Historic Data. https://doi.org/10.5676/DWD_GPCC/CLIM_M_V2011_025.64.Yuzhong, Y., Qingbai, W. & Hanbo, Y. Stable isotope variations in the ground ice of Beiluhe Basin on the Qinghai-Tibet Plateau. Quat. Int. 313–314, 85–91 (2013).Article 

    Google Scholar 
    65.Li, X. et al. Near-surface air temperature lapse rates in the mainland China during 1962–2011. J. Geophys. Res. Atmospheres 118, 7505–7515 (2013).ADS 
    Article 

    Google Scholar 
    66.Revelle, W. psych: Procedures for Personality and Psychological Research, Northwestern University, Evanston, Illinois, USA. Northwest. Univ. (2018).67.Zimmermann, H. H. et al. Sedimentary ancient DNA and pollen reveal the composition of plant organic matter in Late Quaternary permafrost sediments of the Buor Khaya Peninsula (north-eastern Siberia). Biogeosciences Discuss. 1–50 (2016) https://doi.org/10.5194/bg-2016-386.68.Oksanen, J. et al. vegan: Community Ecology Package. (2019).69.Hallett, L. M. et al. codyn: An r package of community dynamics metrics. Methods Ecol. Evol. 7, 1146–1151 (2016).Article 

    Google Scholar 
    70.Barbosa, A. M., Brown, J. A., Jimenez-Valverde, A. & Real, R. modEvA: Model Evaluation and Analysis. R package version 1.3.2. (2016).71.Kuhn, M. caret: Classification and Regression Training. R package version 6.0–86. (2020). More

  • in

    Solar geoengineering can alleviate climate change pressures on crop yields

    1.Lawrence, M. G. et al. Evaluating climate geoengineering proposals in the context of the Paris Agreement temperature goals. Nat. Commun. 9, 3734 (2018).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    2.MacMartin, D. G., Ricke, K. L. & Keith, D. W. Solar geoengineering as part of an overall strategy for meeting the 1.5 °C Paris target. Phil. Trans. R. Soc. A 376, 20160454 (2018).ADS 

    Google Scholar 
    3.Crutzen, P. J. Albedo enhancement by stratospheric sulfur injections: a contribution to resolve a policy dilemma? Climatic Change 77, 211–220 (2006).ADS 
    CAS 

    Google Scholar 
    4.Ahlm, L. et al. Marine cloud brightening—as effective without clouds. Atmos. Chem. Phys. 17, 13071–13087 (2017).ADS 
    CAS 

    Google Scholar 
    5.Muri, H. et al. Climate response to aerosol geoengineering: a multimethod comparison. J. Clim. 31, 6319–6340 (2018).ADS 

    Google Scholar 
    6.Kravitz, B. et al. The Geoengineering Model Intercomparison Project (GeoMIP). Atmos. Sci. Lett. 12, 162–167 (2011).ADS 

    Google Scholar 
    7.Robock, A., Oman, L. & Stenchikov, G. L. Regional climate responses to geoengineering with tropical and Arctic SO2 injections. J. Geophys. Res. Atmos. 113, D16101 (2008).ADS 

    Google Scholar 
    8.Tjiputra, J. F., Grini, A. & Lee, H. Impact of idealized future stratospheric aerosol injection on the large-scale ocean and land carbon cycles. J. Geophys. Res. Biogeosci. 121, 2015JG003045 (2016).
    Google Scholar 
    9.Russell, L. M. et al. Ecosystem impacts of geoengineering: a review for developing a science plan. Ambio 41, 350–369 (2012).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    10.Xia, L. et al. Solar radiation management impacts on agriculture in China: a case study in the Geoengineering Model Intercomparison Project (GeoMIP). J. Geophys. Res. Atmos. 119, 8695–8711 (2014).ADS 

    Google Scholar 
    11.Zhan, P., Zhu, W., Zhang, T., Cui, X. & Li, N. Impacts of sulfate geoengineering on rice yield in china: results from a multimodel ensemble. Earth Future 7, 395–410 (2019).ADS 

    Google Scholar 
    12.Parkes, B., Challinor, A. & Nicklin, K. Crop failure rates in a geoengineered climate: impact of climate change and marine cloud brightening. Environ. Res. Lett. 10, 084003 (2015).13.Yang, H. et al. Potential negative consequences of geoengineering on crop production: a study of Indian groundnut. Geophys. Res. Lett. 43, 11786–11795 (2016).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    14.Pongratz, J., Lobell, D. B., Cao, L. & Caldeira, K. Crop yields in a geoengineered climate. Nat. Clim. Change 2, 101–105 (2012).ADS 
    CAS 

    Google Scholar 
    15.Proctor, J., Hsiang, S., Burney, J., Burke, M. & Schlenker, W. Estimating global agricultural effects of geoengineering using volcanic eruptions. Nature 560, 480–483 (2018).ADS 
    CAS 

    Google Scholar 
    16.Tjiputra, J. F. et al. Evaluation of the carbon cycle components in the Norwegian Earth System Model (NorESM). Geosci. Model Dev. 6, 301–325 (2013).ADS 

    Google Scholar 
    17.MacMartin, D. G. & Kravitz, B. Mission-driven research for stratospheric aerosol geoengineering. Proc. Natl Acad. Sci. USA 116, 1089–1094 (2019).ADS 
    CAS 

    Google Scholar 
    18.Lombardozzi, D. L. et al. Simulating agriculture in the community land model version 5. J. Geophys. Res. Biogeosci. 125, e2019JG005529 (2020).ADS 

    Google Scholar 
    19.Popp, A. et al. Land-use futures in the shared socio-economic pathways. Glob. Environ. Change 42, 331–345 (2017).
    Google Scholar 
    20.O’Neill, B. C. et al. The scenario model intercomparison project (ScenarioMIP) for CMIP6. Geosci. Model Dev. 9, 3461–3482 (2016).ADS 

    Google Scholar 
    21.IPCC Climate Change 2013: The Physical Science Basis (eds Stocker, T. F. et al.) (Cambridge Univ. Press, 2013).22.FAOSTAT (FAO, 2019); http://www.fao.org/faostat/en/?#data/QC23.Tai, A. P. K., Martin, M. V. & Heald, C. L. Threat to future global food security from climate change and ozone air pollution. Nat. Clim. Change 4, 817–821 (2014).ADS 
    CAS 

    Google Scholar 
    24.Hsiao, J., Swann, A. L. S. & Kim, S.-H. Maize yield under a changing climate: the hidden role of vapor pressure deficit. Agric. For. Meteorol. 279, 107692 (2019).ADS 

    Google Scholar 
    25.Grossiord, C. et al. Plant responses to rising vapor pressure deficit. New Phytol. https://doi.org/10.1111/nph.16485 (2020).26.Rigden, A. J., Mueller, N. D., Holbrook, N. M., Pillai, N. & Huybers, P. Combined influence of soil moisture and atmospheric evaporative demand is important for accurately predicting US maize yields. Nat. Food 1, 127–133 (2020).
    Google Scholar 
    27.Konings, A. G., Williams, A. P. & Gentine, P. Sensitivity of grassland productivity to aridity controlled by stomatal and xylem regulation. Nat. Geosci. 10, 284–288 (2017).ADS 
    CAS 

    Google Scholar 
    28.Novick, K. A. et al. The increasing importance of atmospheric demand for ecosystem water and carbon fluxes. Nat. Clim. Change 6, 1023–1027 (2016).ADS 
    CAS 

    Google Scholar 
    29.Ainsworth, E. A. & Long, S. P. What have we learned from 15 years of free-air CO2 enrichment (FACE)? A meta-analytic review of the responses of photosynthesis, canopy properties and plant production to rising CO2: Tansley review. New Phytol. 165, 351–372 (2004).
    Google Scholar 
    30.Bishop, K. A., Leakey, A. D. B. & Ainsworth, E. A. How seasonal temperature or water inputs affect the relative response of C3 crops to elevated CO2: a global analysis of open top chamber and free air CO2 enrichment studies. Food Energy Secur. 3, 33–45 (2014).
    Google Scholar 
    31.Ainsworth, E. A. et al. A meta-analysis of elevated CO2 effects on soybean (Glycine max) physiology, growth and yield. Glob. Change Biol. 8, 695–709 (2002).ADS 

    Google Scholar 
    32.Leakey, A. D. B. Rising atmospheric carbon dioxide concentration and the future of C4 crops for food and fuel. Proc. R. Soc. B 276, 2333–2343 (2009).CAS 

    Google Scholar 
    33.Ainsworth, E. A. & Rogers, A. The response of photosynthesis and stomatal conductance to rising CO2: mechanisms and environmental interactions. Plant Cell Environ. 30, 258–270 (2007).CAS 

    Google Scholar 
    34.National Research Council Climate Intervention: Reflecting Sunlight to Cool Earth (National Academies, 2015); https://doi.org/10.17226/1898835.Lutsko, N. J., Seeley, J. T. & Keith, D. W. Estimating impacts and trade-offs in solar geoengineering scenarios with a moist energy balance model. Geophys. Res. Lett. 47, e2020GL087290 (2020).ADS 

    Google Scholar 
    36.Friedlingstein, P. et al. Uncertainties in CMIP5 climate projections due to carbon cycle feedbacks. J. Clim. 27, 511–526 (2014).ADS 

    Google Scholar 
    37.Tilmes, S. et al. The hydrological impact of geoengineering in the geoengineering model intercomparison project (GeoMIP). J. Geophys. Res. Atmos. 118, 11,036–11,058 (2013).
    Google Scholar 
    38.Lawrence, D. M. et al. The community land model version 5: description of new features, benchmarking, and impact of forcing uncertainty. J. Adv. Model. Earth Syst. 11, 4245–4287 (2019).ADS 

    Google Scholar 
    39.Fisher, R. A. et al. Parametric controls on vegetation responses to biogeochemical forcing in the CLM5. J. Adv. Model. Earth Syst. 11, 2879–2895 (2019).ADS 

    Google Scholar 
    40.Bonan, G. B. et al. Model structure and climate data uncertainty in historical simulations of the terrestrial carbon cycle (1850–2014). Glob. Biogeochem. Cycles 33, 1310–1326 (2019).ADS 
    CAS 

    Google Scholar 
    41.Osborne, T., Rose, G. & Wheeler, T. Variation in the global-scale impacts of climate change on crop productivity due to climate model uncertainty and adaptation. Agric. For. Meteorol. 170, 183–194 (2013).ADS 

    Google Scholar 
    42.Peng, B. et al. Improving maize growth processes in the community land model: implementation and evaluation. Agric. For. Meteorol. 250–251, 64–89 (2018).ADS 

    Google Scholar 
    43.Buzan, J. R. & Huber, M. Moist heat stress on a hotter earth. Annu. Rev. Earth Planet. Sci. 48, 623–655 (2020).ADS 
    CAS 

    Google Scholar 
    44.Wieder, W. R. et al. Beyond static benchmarking: using experimental manipulations to evaluate land model assumptions. Glob. Biogeochem. Cycles 33, 1289–1309 (2019).ADS 
    CAS 

    Google Scholar 
    45.Mercado, L. M. et al. Impact of changes in diffuse radiation on the global land carbon sink. Nature 458, 1014–1017 (2009).ADS 
    CAS 

    Google Scholar 
    46.Cheng, S. J. et al. Variations in the influence of diffuse light on gross primary productivity in temperate ecosystems. Agric. For. Meteorol. 201, 98–110 (2015).ADS 

    Google Scholar 
    47.Shao, L. et al. The fertilization effect of global dimming on crop yields is not attributed to an improved light interception. Glob. Change Biol. 26, 1697–1713 (2020).ADS 

    Google Scholar 
    48.Vattioni, S. et al. Exploring accumulation-mode H2SO4 versus SO2 stratospheric sulfate geoengineering in a sectional aerosol–chemistry–climate model. Atmos. Chem. Phys. 19, 4877–4897 (2019).ADS 
    CAS 

    Google Scholar 
    49.Levis, S., Badger, A., Drewniak, B., Nevison, C. & Ren, X. CLMcrop yields and water requirements: avoided impacts by choosing RCP 4.5 over 8.5. Climatic Change 146, 501–515 (2018).ADS 

    Google Scholar 
    50.Fricko, O. et al. The marker quantification of the Shared Socioeconomic Pathway 2: a middle-of-the-road scenario for the 21st century. Glob. Environ. Change 42, 251–267 (2017).
    Google Scholar 
    51.Lauvset, S. K., Tjiputra, J. & Muri, H. Climate engineering and the ocean: effects on biogeochemistry and primary production. Biogeosciences 14, 5675–5691 (2017).ADS 

    Google Scholar 
    52.Hurtt, G. C. et al. Harmonization of land-use scenarios for the period 1500–2100: 600 years of global gridded annual land-use transitions, wood harvest, and resulting secondary lands. Climatic Change 109, 117–161 (2011).ADS 

    Google Scholar 
    53.West, T. O. et al. Cropland carbon fluxes in the United States: increasing geospatial resolution of inventory-based carbon accounting. Ecol. Appl. 20, 1074–1086 (2010).
    Google Scholar 
    54.Lobell, D. B., Schlenker, W. & Costa-Roberts, J. Climate trends and global crop production since 1980. Science 333, 616–620 (2011).ADS 
    CAS 

    Google Scholar 
    55.Farquhar, G., von Caemmerer, Svon & Berry, J. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta 149, 78–90 (1980).CAS 

    Google Scholar 
    56.Collatz, G. J., Ribas-Carbo, M. & Berry, J. A. Coupled photosynthesis-stomatal conductance model for leaves of C4 plants. Funct. Plant Biol. 19, 519–538 (1992).
    Google Scholar 
    57.Medlyn, B. E. et al. Reconciling the optimal and empirical approaches to modelling stomatal conductance. Glob. Change Biol. 17, 2134–2144 (2011).ADS 

    Google Scholar 
    58.Long, S. P., Ainsworth, E. A., Leakey, A. D. B., Nösberger, J. & Ort, D. R. Food for thought: lower-than-expected crop yield stimulation with rising CO2 concentrations. Science 312, 1918–1921 (2006).ADS 
    CAS 

    Google Scholar 
    59.The NCAR Command Language (NCL, Version 6.5.0) (UCAR, NCAR, CISL, TDD, 2018); https://doi.org/10.5065/D6WD3XH5 More

  • in

    A pan-African spatial assessment of human conflicts with lions and elephants

    After preprocessing the data, methods consisted of spatial analyses to map areas at risk of conflict; statistical analyses to identify the most important factors affecting lion and elephant population numbers; economic analyses to estimate the EAA of building and maintaining mitigation fences in areas under severe and high risk of conflict, and fragmentation analyses to assess the impact of fences on migratory mammal species. We describe each step in detail below (see Supplementary Fig. 2 for a flowchart of the analysis). All spatial data were converted to vectors for analysis to reduce commission errors (when a species is mistakenly thought to be present) when converting the species-range maps from vector to raster. Data preprocessing was carried out using the open source database PostgreSQL 11.4 (https://www.postgresql.org/about/) with the GIS extensions of PostGIS 2.5 (https://postgis.net/); conflict mapping and range fragmentation analyses used PostgreSQL 11.4 and PostGIS 2.5, and Python v. 3.7.060; statistical and economic analyses used R v. 3.6.061; sensitivity analyses used PostgreSQL 11.4 and PostGIS 2.5, and Python v. 3.7.060 and R v. 3.6.061.PreprocessingHuman pressuresHuman pressure layers were independently generated from this study. We used Gridded Population of the World Version 4 (GPWv4) as a layer for human population density62. GPWv4 is a minimally modelled data set consisting of estimates of human population (number of persons per raster grid cell) based on non-spatial population data (i.e., tabular counts of population listed by administrative area) and spatially explicit administrative boundary data. Population input data are collected at the most detailed spatial resolution available from the results of the 2010 round of Population and Housing Censuses. The input data are then extrapolated to 2020 using calculated growth rates to produce future population estimates. A proportional allocation gridding algorithm, utilizing ~13.5 million national and subnational administrative units, assigned population counts to 30 arcsecond (~1 km at the equator) grid cells. The population density rasters were created by dividing the population count raster for a given target year by the land-area raster.We used the most recent version of the Gridded Livestock of the World database63, reflecting the compiled and harmonized subnational livestock distribution data for 2010, to extract information on cattle density. The data set provides global population densities of cattle, buffaloes, horses, sheep, goats, pigs, chickens, and ducks in each land pixel at a spatial resolution of 0.083333 decimal degrees (~10 km at the equator). Detailed livestock census statistics are mined from agricultural yearbooks or through direct contacts with ministries or statistical bureaus. The census statistics are usually found in the form of numbers per administrative unit that must be linked to corresponding geographic information system boundaries. Densities are estimated in each census polygon by dividing the number of animals from the census by the surface area of the administrative unit polygon (estimated in an Albert equal-area projection), corrected by a mask excluding unsuitable areas. Livestock densities were then extracted from the subnational census data and were used as the dependent variable in Random Forest models to estimate a density value in each pixel, based on raster predictor variables.We used spatially detailed crop maps available from the Copernicus Global Land Cover map at ~0.001° (~100 m) resolution64. The land-cover map is a discrete map with ten continuous cover fractions (nine base land-cover classes and seasonal water) to provide spatial information about land for a diversity of applications, including biodiversity conservation. Cropland (as percentage of 100 m pixel that is covered by cropland) refers to cultivated and managed agriculture, but does not include perennial woody crops that are classified under the appropriate forest or shrub land-cover type64. Cropland also refers to both irrigated and rainfed agriculture. The land-cover map was generated by compiling the 5-daily PROBA-V multi-spectral image data with a Ground Sampling Distance of ~0.001° as the primary earth observation data and PROBA-V UTM daily multi-spectral image data with a Ground Sampling Distance of ~0.003° (~300 m) as the secondary earth observation data. Next, the 5-daily PROBA-V 100 m and daily 300 m datasets were fused using a Kalman filtering approach. The global overall accuracy of the product for the base year 2015 was calculated through an independent pre-validation and reached 80%.Species-range mapsUpdated range maps showing current distribution for lions and elephants were provided by the International Union for Conservation of Nature (IUCN) Cat and African Elephant Specialist Groups65. In addition to the range maps, the specialist groups provided information on the number of African lions (2018) and elephants (2016) within sites where they are still extant. We also obtained species-range maps for all terrestrial mammal species in orders Cetartiodactyla, Perissodactyla, Primates, and Carnivora occurring in Africa from the IUCN Red List portal (www.iucnredlist.org/). Mammal species in these orders include migratory mammal species (e.g., the common wildebeest Connochaetes taurinus), which might be negatively affected by mitigation fences, e.g., by potentially blocking migratory routes.Protected areasThe data on protected areas were based on the May 2019 release of the World Database on Protected Areas66 (retrieved from http://www.protectedplanet.net). To prevent overestimation of the area coverage of protected areas caused by overlapping designations, we merged polygons into a single layer. We only included in the analysis IUCN categories Ia (Strict Nature Reserve), Ib (Wilderness Area), II (National Park), III (Natural Monument or Feature), and IV (Habitat/Species Management Area), because we wanted to prevent fences from excluding people from protected areas that had been modified by the interaction of nature and people over time (e.g., V, Protected Landscape/Seascape).Mapping potential risk of conflictA database on the spatial distribution of conflict locations between humans and lions and elephants is not available across Africa. We therefore mapped the most prominent factors known to affect conflict: human population density (for both lion and elephant), crop raiding (elephants), and cattle killing (lions)8. Furthermore, spatial modelling of range contractions in carnivores showed that contractions were significantly more likely in regions with high rural human population density, cattle density, and/or cropland4. Therefore, we only retained areas where human, cattle, and crop densities were in the first decile (in our case, the first decile is the decile with the highest human population, crop, and cattle densities) by PostgreSQL/PostGIS. Using only the highest decile likely resulted in a conservative map of spatial conflict.We further classified areas at the highest potential conflict into low, moderate, high, or severe risk of conflict. Specifically, areas at severe risk of conflict are those where the highest human population, crop, and cattle densities all overlap; areas at high risk of conflict are those with overlaps between the highest densities of human population and either crops or cattle; areas at moderate risk of conflict are the areas where the highest crop and cattle densities overlap; and areas at low risk of conflict are those with only one human pressure, i.e., the highest human population, or crop, or cattle density. The remainder was considered as being at no risk of conflict, as it did not meet any of the above criteria, but note the conservative nature of our analysis (see above).The lion and elephant range maps and the protected area layer were intersected to select all protected areas that contain parts of lion and elephant range and/or were adjacent to the species-range maps. The identified protected areas were then merged with the species-range maps to create a new extended range layer (see for an example in Supplementary Fig. 1). These extended range maps were used (i) to identify potential areas where lions and elephants could be restored, and (ii) to avoid interrupting ecological processes (e.g., migrations) and/or causing unintended consequences (e.g., fragment populations) to other biodiversity in neighbouring protected areas.We then identified areas at risk of conflict by intersecting the extended range map layer for lions and elephants with the classified conflict map. In all cases, the intersections were carried out so that the classified conflict areas were either adjacent to, or within a distance of 10 km from, the edge of the extended range map layer. We set this distance to consider the wide-ranging behaviour of both lions and elephants, to account for the fact that conflict decreases at greater distances from protected area boundaries37,38, and to account for the fact that future human pressures will likely increase before conservation actions take place2.We assessed how robust our results were to commission (where human pressure is mistakenly assumed to exist) and omission (where human pressure is mistakenly assumed to be absent) errors in the human pressure maps by carrying out a sensitivity analysis that randomly varied the distances between the extended range maps and the human pressure maps. We first used Latin hypercube sampling, which is a form of sampling used to reduce the number of runs necessary for a Monte Carlo simulation to achieve a reasonably accurate random distribution67, to randomly vary 100 times the distance values between the extended range and human pressure maps. Specifically, we divided the low, moderate, high, and severe conflict lines into 100 m segments, calculated the minimum distance for each segment to human pressure within a 10, 20, and 30 km buffer distance from the edge of the extended range map layer, and then randomly varied that distance 100 times across ±10% of the value. We then averaged the resulting 100 randomly created distance values for each segment and identified which segments fell outside of the analyzed buffer distances of 10, 20, and 30 km. We tested for 20 and 30 km buffer distances, as we wanted to assess the variability of the fencing distance to different buffer sizes. We also estimated the certainty of lion and elephant presence by identifying segments of the perimeter of the range maps of lion and elephant that overlapped with protected areas. We did this as we had information on certain presence of both species from within protected areas, as opposed to areas extending outside of protected areas.Statistical analysesWe used an information theoretic approach68 and Bayesian information criterion to calculate statistical models. We used generalized linear mixed models with a negative‐binomial error distribution to account for over-dispersed count data and a log‐link function to examine factors affecting lion (n = 77) and elephant (n = 191) population sizes in Africa. Generalized linear mixed models were fitted with both random and fixed effects, to capture the data structure. Country was included as a random intercept to represent the hierarchical structure of the data. All variables listed in Supplementary Table 8 were fitted as fixed effects, i.e., with constant regression coefficients across countries. The site-specific variables were calculated only for sites where lions and elephants are currently present and not for the extended ranges. For transboundary sites that stretch across countries, we used the value for Gross Domestic Product, Conservation expenditure, and the Ibrahim Index of African Governance, for the country making the largest area contribution to the site. We compared and ranked models using the Bayesian information criterion68. To avoid multicollinearity among variables, we only selected variables with the strongest effect on population numbers that correlated at r 0.7 was selected as an input into the modelling process. We assessed each model’s relative probability, using Bayesian information criterion weights and the structural goodness-of-fit from the percentage of deviance explained by the model. We determined the magnitude and direction of the coefficients for the independent variables with multi-model averaging implemented in the R package glmulti69. The relative importance of each predictor variable was measured as the sum of the weights over the six top‐ranked models with Bayesian information criterion values closer to that of the best model containing the parameter of interest. Finally, we used a 10-fold cross-validation (a bootstrap resampling procedure using 1000 iterations) to assess the predictive ability of the top-ranked model.Range fragmentation analysesWe assessed how the proposed mitigation fences affected species-range connectivity by calculating the perimeter length-to-area ratio for mammal species in orders Cetartiodactyla, Perissodactyla, Primates, and Carnivora, whose ranges were identified as intersecting with areas at severe risk of conflict. Minimizing the perimeter length-to-area ratio is an important method of optimizing protected area design, resulting in compact reserves with high connectivity that can enhance persistence of the species. The smaller the ratio, the greater the clumping and connectivity of the species ranges. Specifically, we calculated the ratios of perimeter length to area for the ranges of 20 migratory mammalian species (i) under current conditions without fences and (ii) under future conditions where the identified mitigation fences would pass through their ranges. In the latter case, we used a 20 m buffer around the identified fences to account for further habitat clearance due to maintaining clearances around the fences for management purposes.Economic analysesWe used EAA to estimate the return on investment of building and maintaining mitigation fences to reduce cattle loss and crop damage. EAA calculates the constant annual cash flow generated by a project over its lifespan if it were an annuity and the annuity can then be compared to other projects of similar or different lifespan. Therefore, the measure potentially provides an important means for funders/donors to compare different investment opportunities. EAA is calculated by dividing the NPV of a project by the present value of annuity factor39. We started by calculating NPV in countries with areas at severe and high risk of conflict as:$${{NPV}}={sum }_{i}^{n}frac{{R}_{i}}{{left(1+dright)}^{i}}-Z$$
    (1)
    where ({R}_{i}) is net cash flow, (d) is the discount rate specific to each country (Supplementary Table 9), n is the number of time periods, (i) is the cash flow period, and (Z) is the initial investment of building the fences. NPV was calculated over a 10-year investment period. ({R}_{i}) was calculated as:$${R}_{i}=B-C$$
    (2)
    where (B) is the economic benefit derived from mitigation fences and (C) is the cost of maintaining mitigation fences. The economic benefits of mitigation fences for countries with severe risk of conflict refer to the potential reduction in cattle loss (for lions) and crop damage (for elephants) derived from building fences:$$B=L+E$$
    (3)
    where (L) represents the economic benefits of reducing cattle loss and (E) measures the economic benefits of reducing crop damage. For countries with high risk of conflict, the benefit ((B)) is derived from one or the other, i.e., (B) = (L) or (B) = (E).$$L=v * w * P$$
    (4)
    where (v) is the number of cattle that are not lost because of the presence of fences, (w) is the average weight in kg of adult cattle in that country, and (P) is the price of meat per kg paid to producers in that country in 2017 (data can be downloaded from http://www.fao.org/faostat/en/#data/PP). (v) was calculated as the percentage of total cattle present in the 10 km buffer adjacent to severe and high conflict areas, which could potentially be killed, based on published estimates across Africa45. Estimates range from 0.8 to 2.6% of cattle losses, and we decided to use a conservative 1% loss in the analysis (see below for how we accounted for uncertainty in model parameters). (w) was based on the average weight of an adult cow with estimates available at a regional level (west Africa: 262 kg; central Africa: 281 kg; east Africa: 283 kg; and southern Africa: 339 kg)70.$$E=left(d * Aright) * y * P$$
    (5)
    where (d) is the percentage of crop area damaged by elephants; data are taken from published estimates (ranging from 0.2 to 4% and we used a conservative 1% in the analysis)44; (A) is the total area in km2 available as crops in the 10 km buffer adjacent to the areas at severe and high risk of conflict; (y) is the yield (ton/km2) for the crop known to be targeted by elephants (cassava, maize, millet, banana, sorghum, groundnuts)44, which covered the largest area size in that country in 2017 (data calculated from: http://www.fao.org/faostat/en/#data/PP); and (P) is the price per ton paid to producers for that crop in that country in 2017 (data can be downloaded from http://www.fao.org/faostat/en/#data/PP). Although there might be several crops available within the buffer, this information is currently not available at the continental scale. Therefore, we decided to use the most common cultivated crop known to be targeted by elephants in each country.The cost of maintaining mitigation fences ((C)) was calculated as:$$C=f * c$$
    (6)
    where f is the fence length in that country and (c) is the cost for maintaining the fence. We obtained cost estimates of building (Z) and maintaining ((c)) the fences from Pekor et al.33. We used the median estimated current cost of USD 9522 per km for building fences and the median stated annual budget cost of USD 487 per km for adequate fence inspection and maintenance. This is the most up-to-date information validated through peer review on the costs (converted to 2017 USD) across Africa33. Cost estimates varied across surveyed conservation areas because of fence height and materials but included relevant costs of electrification and predator-proof structures33. The data were collected from 29 partially fenced ( More

  • in

    Precision conservation for a changing climate

    The author wishes to thank N. Haan, G. P. Robertson and N. Haddad for their valuable comments. Financial support was provided by USDA/NIFA (awards 2019-67012-29595 and 2015-68007-23133), the US National Science Foundation’s Long-term Ecological Research Program (award 1637653), the US Department of Energy, Office of Science, Office of Biological and Environmental Research (awards DESC0018409 and DE-FC02-07ER64494), and Michigan State University AgBioResearch. More

  • in

    Controversial forestry experiment will be largest-ever in United States

    A clear-cut slope in the Elliott State Forest, Oregon.Credit: Matthew Betts

    Despite lingering tensions among environmentalists and loggers, a plan to launch the largest forestry experiment in the United States — and perhaps the world — last month cleared a major hurdle. Controversially, the study would allow logging in a new research forest, in an attempt to answer a grand question: in a world where wood remains a necessary resource, but biodiversity is declining, what’s the best way to balance timber production with conservation?“We all love wood, and we all need wood,” says Thomas DeLuca, dean of the College of Forestry at Oregon State University (OSU) in Corvallis. “We have to find ways to produce it sustainably, and this project could help us do that.”
    These scientists are setting a forest on fire — and studying it with drones
    If the project — proposed by DeLuca and other researchers at OSU — launches successfully, the newly created Elliott State Research Forest in southwestern Oregon would occupy a roughly 33,000-hectare parcel of land. This would be divided into more than 40 sections, in which scientists would test several forest-management strategies, some including extensive logging. The advisory committee for the project, which comprises environmentalists, hunters, loggers and members of local Indigenous tribes, approved the latest research proposal on 22 April.The plan comes as US President Joe Biden and other international leaders are strengthening commitments to conserve land and biodiversity before a meeting of the United Nations Convention on Biological Diversity later this year. In time, the Elliott research forest could help policymakers to determine how best to define and implement those pledges, says DeLuca.A contested forestFor decades, the land that makes up the Elliott State Forest has been mired in controversy. Logging is big business in the US Pacific Northwest, and this particular state-owned piece of land contains old-growth forest filled with valuable Douglas firs (Pseudotsuga menziesii) and other trees. Other sections have been actively logged and replanted since 1930. It also hosts threatened species such as the spotted owl (Strix occidentalis) and the marbled murrelet (Brachyramphus marmoratus), a seabird that nests in old-growth forests. In 2012, a lawsuit aimed at protecting the marbled murrelet brought commercial logging in the forest to a halt.

    Source: Deanne Carlson, OSU/Oregon Geospatial Enterprise Office; Design: Nature

    The state of Oregon considered various options for the land, before OSU researchers stepped forward with a plan in 2018. Their proposal to convert the property into a research forest would allow logging to resume at a lower level — but in the service of science and conservation, the scientists say. According to the plan, the profit from logging in the Elliott forest — around US$5 million to $7 million annually, says DeLuca — would help to pay for the experiment’s infrastructure and operations.There are dozens of research forests around the globe, including in the United States, and scientists have used them to study everything from ecology and soils to acid rain and the effects of rising atmospheric carbon dioxide levels. But the Elliott research forest would be larger than most of its predecessors, and advocates say that it would provide scientists with the first opportunity to test ecological forestry at such a large scale.A sea change for forestryAs currently designed, the project would leave more than 40% of the forest — a section of old growth that has been regenerating naturally since the area last burnt, a century and a half ago — untouched by logging. In the remaining area, researchers would run a series of replicated experiments, carrying out 4 types of land management across 40 small watersheds. On some plots, selective logging of individual trees would take place across the entire area. On others, clear-cutting would take place on half of the land, with the other half reserved for conservation. Other types of experimental plot would mix these two approaches (see ‘A grand experiment’). To understand the impacts of each management type, scientists would measure a variety of parameters, including levels of carbon in the forest; stream health; and insect, bird and fish diversity.So far, around 20 OSU researchers are involved in the project’s design, but the university hopes eventually to attract more scientists from across the world, who would run their own projects.

    Old-growth areas such as this one would be protected as reserves under a plan to convert the Elliott State Forest into an experimental forest.Credit: Matthew Betts

    The scale and approach of the experiment proposed by OSU would represent a sea change in forestry research, says Sue Baker, a forest ecologist at the University of Tasmania in Hobart, Australia, who is setting up a retrospective study looking at similar questions in Tasmanian forests. “I can’t think of anything similar anywhere in the world where people have been able to manipulate the forest landscape at this scale,” says Baker.More hurdles aheadSince its creation in 1930, the Elliott State Forest has been legally obliged to generate revenue for Oregon’s public schools through logging. Before OSU can take it over, the state must compensate the school fund to the tune of $221 million (the value of the forest); it has so far allocated less than half of that amount.And other hurdles remain. The university must finish a detailed management plan that will lay out rules governing the forest, and it must craft a separate plan for managing threatened and endangered species; this will need to be approved by the US Fish and Wildlife Service.
    How much can forests fight climate change?
    The OSU team has spent the past few years trying to build a broad — and unlikely — coalition for the effort, through public meetings and engagement with local Indigenous tribes, industry, environmentalists and other members of the project’s advisory panel, whose support will be crucial as state leaders weigh their final decision. But tensions haven’t disappeared entirely. Many environmentalists continue to question the logic of clear-cutting forests that absorb and store carbon in the middle of a climate crisis. Rather than perpetuating a long and damaging legacy of clear-cutting, the Elliott forest could be used to pioneer new forestry methods that restore biodiversity and boost carbon storage, says Josh Laughlin, executive director of Cascadia Wildlands, a conservation group based in Eugene, Oregon. “Let’s not make the same mistakes we’ve made over the past 100 years.”Given OSU’s long-standing ties to the timber industry, and controversies surrounding its management of existing research forests, it will also need to overcome scepticism about its role as a land steward, says Bob Van Dyk, a policy director at the Wild Salmon Center, an environmental group based in Portland, Oregon. In 2019, for instance, OSU’s College of Forestry authorized clear-cutting on 6.5 hectares of one of its forests, felling trees that were hundreds of years old.
    When will the Amazon hit a tipping point?
    DeLuca acknowledges that there have been mistakes in the past, but says the university has a solid academic record, and is committed to building a world-class research facility with the Elliott forest. “If we are able to demonstrate practices that accommodate the broadest array of species while still generating timber for meeting human resource needs, we can have a much larger impact,” says DeLuca.Everything will depend on the final management plan, but for now, Van Dyk and other members of the advisory board have unanimously given a provisional green light to the latest proposal. “It’s a good project, and we don’t get many chances to do something that is truly novel and interesting,” he says. More

  • in

    Our radical changes to Earth’s greenery began long ago — with farms, not factories

    A nineteenth-century illustration of a harvest in ancient Greece. Farming intensified around 2000 BC, when the rate of change in Earth’s plant life sped up. Credit: Docutres/Index/Heritage Images/Alamy

    Ecology
    20 May 2021
    Our radical changes to Earth’s greenery began long ago — with farms, not factories

    Humanity’s imprint on plant species and abundance began roughly 4,000 years ago, when agriculture took off.

    Share on Twitter
    Share on Twitter

    Share on Facebook
    Share on Facebook

    Share via E-Mail
    Share via E-Mail

    Human activity began to transform the number and variety of plant species on Earth thousands of years ago, long before the Industrial Revolution, and might have had an even greater impact on vegetation than did the last ice age.Ice entombed much of the planet from roughly 115,000 to some 20,000 years ago. Then, massive glaciers around the world started to retreat and global temperatures rose, resulting in dramatic alterations to Earth’s ecosystems.To investigate how the abundance and composition of global vegetation changed after that thaw, Ondřej Mottl and Suzette Flantua at the University of Bergen in Norway and their colleagues analysed 1,181 fossilized pollen samples from the past 18,000 years. The pollen came from all continents except Antarctica.The researchers found that global vegetation has been transformed, first by the climate changes that accompanied the end of the last glacial period. However, starting about 4,000 years ago, when agriculture intensified, the pace of change in global vegetation accelerated, reaching or exceeding the rate of change at the end of the most recent ice age.

    Science (2021)

    Ecology More

  • in

    Seasonal change is a major driver of soil resistomes at a watershed scale

    1.D’Costa, V. M. et al. Antibiotic resistance is ancient. Nature. 477, 457–461 (2011).PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    2.Allen, H. K. et al. Call of the wild: antibiotic resistance genes in natural environments. Nat. Rev. Microbiol. 8, 251–259 (2010).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    3.Udikovic-Kolic, N., Wichmann, F., Broderick, N. A. & Handelsman, J. Bloom of resident antibiotic-resistant bacteria in soil following manure fertilization. Proc. Natl Acad. Sci. USA. 111, 15202–15207 (2014).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    4.Chen, Q. L. et al. Long-term field application of sewage sludge increases the abundance of antibiotic resistance genes in soil. Environ. Int. 92–93, 1–10 (2016).PubMed 
    Article 
    CAS 
    PubMed Central 

    Google Scholar 
    5.Gillings, M. R. & Stokes, H. W. Are humans increasing bacterial evolvability? Trends Ecol. Evol. 27, 346–352 (2012).PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    6.Zhu, Y. G. et al. Diverse and abundant antibiotic resistance genes in Chinese swine farms. Proc. Natl Acad. Sci. USA. 110, 3435–3440 (2013).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    7.Woods, L. C. et al. Horizontal gene transfer potentiates adaptation by reducing selective constraints on the spread of genetic variation. Proc. Natl Acad. Sci. USA. 117, 26868–26875 (2020).8.World Health Organization. Antimicrobial resistance: global report on surveillance. World Health Organization. (2014).9.Forsberg, K. J. et al. The shared antibiotic resistome of soil bacteria and human pathogens. Science. 337, 1107–1111 (2012).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    10.Zhu, G. et al. Air pollution could drive global dissemination of antibiotic resistance genes. ISME J. 15, 270–281 (2021).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    11.Xiang, Q. et al. Agricultural activities affect the pattern of the resistome within the phyllosphere microbiome in peri-urban environments. J. Hazard Mater. 382, 121068 (2020).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    12.Wang, F. H. et al. High throughput profiling of antibiotic resistance genes in urban park soils with reclaimed water irrigation. Environ. Sci. Technol. 48, 9079–9085 (2014).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    13.Ding, J. et al. Long-term application of organic fertilization causes the accumulation of antibiotic resistome in earthworm gut microbiota. Environ. Int. 124, 145–152 (2019).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    14.Zhou, S. Y. et al. Phyllosphere of staple crops under pig manure fertilization, a reservoir of antibiotic resistance genes. Environ. Pollut. 252, 227–235 (2019).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    15.Wang, F. H., Qiao, M., Chen, Z., Su, J. Q. & Zhu, Y. G. Antibiotic resistance genes in manure-amended soil and vegetables at harvest. J. Hazard Mater. 299, 215–221 (2015).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    16.Marti, R. et al. Impact of manure fertilization on the abundance of antibiotic-resistant bacteria and frequency of detection of antibiotic resistance genes in soil and on vegetables at harvest. Appl. Environ. Microb. 79, 5701–5709 (2013).CAS 
    Article 

    Google Scholar 
    17.Zhu, Y. G. et al. Continental-scale pollution of estuaries with antibiotic resistance genes. Nat. Microbiol. 2, 16270 (2017).CAS 
    PubMed 
    Article 

    Google Scholar 
    18.Du, S. et al. Large-scale patterns of soil antibiotic resistome in Chinese croplands. Sci. Total Environ. 712, 136418 (2020).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    19.Pruden, A., Pei, R. T., Storteboom, H. & Carlson, K. H. Antibiotic resistance genes as emerging contaminants: studies in northern Colorado. Environ. Sci. Technol. 40, 7445–7450 (2006).CAS 
    PubMed 
    Article 

    Google Scholar 
    20.Bahram, M. et al. Structure and function of the global topsoil microbiome. Nature. 560, 233–237 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    21.Hu, H. W. et al. Diversity of herbaceous plants and bacterial communities regulates soil resistome across forest biomes. Environ. Microbiol. 20, 3186–3200 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    22.Han, X. M. et al. Antibiotic resistance genes and associated bacterial communities in agricultural soils amended with different sources of animal manures. Soil Biol. Biochem. 126, 91–102 (2018).CAS 
    Article 

    Google Scholar 
    23.Hu, H. W. et al. Temporal changes of antibiotic-resistance genes and bacterial communities in two contrasting soils treated with cattle manure. FEMS Microbiol. Ecol. 92, fiv169 (2016).24.Zhang, Y. J. et al. Temporal succession of soil antibiotic resistance genes following application of swine, cattle and poultry manures spiked with or without antibiotics. Environ. Pollut. 231, 1621–1632 (2017).CAS 
    PubMed 
    Article 

    Google Scholar 
    25.Zhou, J. et al. Reproducibility and quantitation of amplicon sequencing-based detection. ISME J. 5, 1303–1313 (2011).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    26.Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods. 7, 335–336 (2010).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    27.Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 26, 2460–2461 (2010).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    28.Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microb. 73, 5261–5267 (2007).CAS 
    Article 

    Google Scholar 
    29.Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    30.Su, J. Q. et al. Antibiotic resistome and its association with bacterial communities during sewage sludge composting. Environ. Sci. Technol. 49, 7356–7363 (2015).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    31.Ouyang, W. Y., Huang, F. Y., Zhao, Y., Li, H. & Su, J. Q. Increased levels of antibiotic resistance in urban stream of Jiulongjiang River, China. Appl. Microbiol. Biotechnol. 99, 5697–5707 (2015).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    32.Roberts D. W. labdsv: ordination and multivariate analysis for ecology. R package version 1.8-0. 2016. https://CRAN.R-project.org/package=labdsv.33.Oksanen J. et al. Vegan: community ecology package. R package version 2.2-0. 2014. http://CRAN.R-project.org/package=vegan.34.Jiao, S. et al. Soil microbiomes with distinct assemblies through vertical soil profiles drive the cycling of multiple nutrients in reforested ecosystems. Microbiome. 6, 1–13 (2018).35.Sloan, W. T. et al. Quantifying the roles of immigration and chance in shaping prokaryote community structure. Environ Microbiol. 8, 732–740 (2006).PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    36.Ning, D., Deng, Y., Tiedje, J. M. & Zhou, J. A general framework for quantitatively assessing ecological stochasticity. Proc. Natl Acad. Sci. USA. 116, 16892–16898 (2019).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    37.De Caceres, M. & Legendre, P. Associations between species and groups of sites: indices and statistical inference. Ecology. 90, 3566–3574 (2009).PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    38.Doerks, T., Copley, R. R., Schultz, J., Ponting, C. P. & Bork, P. Systematic identification of novel protein domain families associated with nuclear functions. Genome Res. 12, 47–56 (2002).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    39.Wickham H. ggplot2: elegant graphics for data analysis. (Springer-Verlag, 2009).40.Kassambara A. ggpubr: ‘ggplot2’ based publication ready plots. R package version 0.2. 2018. https://CRAN.R-project.org/package=ggpubr.41.Ahlmann-Eltze C. ggsignif: significance brackets for ‘ggplot2’. R package version 0.4. 0. 2018. https://CRAN.R-project.org/package=ggsignif.42.Zhao, F. K. et al. Soil contamination with antibiotics in a typical peri-urban area in eastern China: seasonal variation, risk assessment, and microbial responses. J. Environ. Sci. (China). 79, 200–212 (2019).Article 

    Google Scholar 
    43.Zhang, Y., Snow, D. D., Parker, D., Zhou, Z. & Li, X. Intracellular and extracellular antimicrobial resistance genes in the sludge of livestock waste management structures. Environ. Sci. Technol. 47, 10206–10213 (2013).CAS 
    PubMed 
    Article 

    Google Scholar 
    44.Mao, D. et al. Persistence of extracellular DNA in river sediment facilitates antibiotic resistance gene propagation. Environ. Sci. Technol. 48, 71–78 (2014).CAS 
    PubMed 
    Article 

    Google Scholar 
    45.Xiang, Q. et al. Spatial and temporal distribution of antibiotic resistomes in a peri-urban area is associated significantly with anthropogenic activities. Environ. Pollut. 235, 525–533 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    46.Forsberg, K. J. et al. Bacterial phylogeny structures soil resistomes across habitats. Nature. 509, 612–616 (2014).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    47.Li, B. et al. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. ISME J. 9, 2490–2502 (2015).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    48.Hu, H. W. et al. Field-based evidence for copper contamination induced changes of antibiotic resistance in agricultural soils. Environ. Microbiol. 18, 3896–3909 (2016).CAS 
    PubMed 
    Article 

    Google Scholar 
    49.Birgander, J., Rousk, J. & Olsson, P. A. Comparison of fertility and seasonal effects on grassland microbial communities. Soil Biol. Biochem. 76, 80–89 (2014).CAS 
    Article 

    Google Scholar 
    50.Fournier, B. et al. Higher spatial than seasonal variation in floodplain soil eukaryotic microbial communities. Soil Biol. Biochem. 147, 107842 (2020).CAS 
    Article 

    Google Scholar 
    51.Zhang, K., Delgado-Baquerizo, M., Zhu, Y. G. & Chu, H. Space is more important than season when shaping soil microbial communities at a large spatial scale. Msystems. 5, e00783–19 (2020).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    52.Ladau, J. & Eloe-Fadrosh, E. A. Spatial, temporal, and phylogenetic scales of microbial ecology. Trends Microbiol. 27, 662–669 (2019).CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar  More

  • in

    Metaplasmidome-encoded functions of Siberian low-centered polygonal tundra soils

    1.Brilli M, Mengoni A, Fondi M, Bazzicalupo M, Liò P, Fani R. Analysis of plasmid genes by phylogenetic profiling and visualization of homology relationships using Blast2Network. BMC Bioinformatics. 2008;9:551.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    2.Dziewit L, Pyzik A, Szuplewska M, Matlakowska R, Mielnicki S, Wibberg D, et al. Diversity and role of plasmids in adaptation of bacteria inhabiting the Lubin copper mine in Poland, an environment rich in heavy metals. Front Microbiol. 2015;6:152.PubMed 
    PubMed Central 

    Google Scholar 
    3.Matyar F, Kaya A, Dinçer S. Antibacterial agents and heavy metal resistance in Gram-negative bacteria isolated from seawater, shrimp and sediment in Iskenderun Bay, Turkey. Sci Total Environ. 2008;407:279–85.CAS 
    PubMed 
    Article 

    Google Scholar 
    4.Dagan T, Artzy-Randrup Y, Martin W. Modular networks and cumulative impact of lateral transfer in prokaryote genome evolution. Proc Natl Acad Sci USA. 2008;105:10039–44.CAS 
    PubMed 
    Article 

    Google Scholar 
    5.Heuer H, Smalla K. Plasmids foster diversification and adaptation of bacterial populations in soil. FEMS Microbiol Rev. 2012;36:1083–104.CAS 
    PubMed 
    Article 

    Google Scholar 
    6.Morozova D, Möhlmann D, Wagner D. Survival of methanogenic archaea from Siberian permafrost under simulated Martian thermal conditions. Orig Life Evol Biosph. 2007;37:189–200.CAS 
    PubMed 
    Article 

    Google Scholar 
    7.Schimel J, Balser TC, Wallenstein M. Microbial stress-response physiology and its implications for ecosystem functioning. Ecology. 2007;88:1386–94.Article 

    Google Scholar 
    8.Margesin R, Miteva V. Diversity and ecology of psychrophilic microorganisms. Res Microbiol. 2011;162:346–61.PubMed 
    Article 

    Google Scholar 
    9.Dutta H, Dutta A. The microbial aspect of climate change. Energy, Ecol Environ. 2016;1:209–32.Article 

    Google Scholar 
    10.Leplae R, Lima-Mendez G, Toussaint A. A first global analysis of plasmid encoded proteins in the ACLAME database. FEMS Microbiol Rev. 2006;30:980–94.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    11.Couturier M, Bex F, Bergquist PL, Maas WK. Identification and classification of bacterial plasmids. Microbiol Rev. 1988;52:375–95.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    12.Carattoli A, Bertini A, Villa L, Falbo V, Hopkins KL, Threlfall EJ. Identification of plasmids by PCR-based replicon typing. J Microbiol Methods. 2005;63:219–28.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    13.Johnson TJ, Wannemuehler YM, Johnson SJ, Logue CM, White DG, Doetkott C, et al. Plasmid replicon typing of commensal and pathogenic Escherichia coli isolates. Appl Environ Microbiol. 2007;73:1976–83.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    14.Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science. 1997;278:631–7.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    15.Smalla K, Jechalke S, Top EM. Plasmid detection, characterization, and ecology. Microbiol Spectr. 2015;3:PLAS-0038-2014.16.Top EM, Holben WE, Forney LJ. Characterization of diverse 2,4-dichlorophenoxyacetic acid-degradative plasmids isolated from soil by complementation. Appl Environ Microbiol. 1995;61:1691–8.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    17.Sayler GS, Hooper SW, Layton AC, King JMH. Catabolic plasmids of environmental and ecological significance. Microb Ecol. 1990;19:1–20.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    18.Elsas JD, Bailey MJ. The ecology of transfer of mobile genetic elements. FEMS Microbiol Ecol. 2006;42:187–97.Article 

    Google Scholar 
    19.Barrón MD, La C, Merlin C, Guilloteau H, Montargès-Pelletier E, Bellanger X. Suspended materials in river waters differentially enrich class 1 integron- and IncP-1 plasmid-carrying bacteria in sediments. Front Microbiol. 2018;9:1443.Article 

    Google Scholar 
    20.Dziewit L, Bartosik D. Plasmids of psychrophilic and psychrotolerant bacteria and their role in adaptation to cold environments. Front Microbiol. 2014;5:596.PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    21.McCann CM, Christgen B, Roberts JA, Su JQ, Arnold KE, Gray ND, et al. Understanding drivers of antibiotic resistance genes in High Arctic soil ecosystems. Environ Int. 2019;125:497–504.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    22.Nesme J, Simonet P. The soil resistome: a critical review on antibiotic resistance origins, ecology and dissemination potential in telluric bacteria. Environ Microbiol. 2015;17:913–30.PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    23.Sjölund M, Bonnedahl J, Hernandez J, Bengtsson S, Cederbrant G, Pinhassi J, et al. Dissemination of multidrug-resistant bacteria into the Arctic. Emerg Infect Dis. 2008;14:70–72.PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    24.Perron GG, Whyte L, Turnbaugh PJ, Goordial J, Hanage WP, Dantas G, et al. Functional characterization of bacteria isolated from ancient Arctic soil exposes diverse resistance mechanisms to modern antibiotics. PLoS One. 2015;10:e0069533.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    25.Hernández J, González-Acuña D. Anthropogenic antibiotic resistance genes mobilization to the polar regions. Infect Ecol Epidemiol. 2016;6:32112.PubMed 
    PubMed Central 

    Google Scholar 
    26.Tan L, Li L, Ashbolt N, Wang X, Cui Y, Zhu X, et al. Arctic antibiotic resistance gene contamination, a result of anthropogenic activities and natural origin. Sci Total Environ. 2018;621:1176–84.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    27.Wang F, Stedtfeld RD, Kim OS, Chai B, Yang L, Stedtfeld TM, et al. Influence of soil characteristics and proximity to antarctic research stations on abundance of antibiotic resistance genes in soils. Environ Sci Technol. 2016;50:12621–9.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    28.Winkel M, Mitzscherling J, Overduin PP, Horn F, Winterfeld M, Rijkers R, et al. Anaerobic methanotrophic communities thrive in deep submarine permafrost. Sci Rep. 2018;8:1291.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    29.Liebner S, Harder J, Wagner D. Bacterial diversity and community structure in polygonal tundra soils from Samoylov Island, Lena Delta, Siberia. Int Microbiol. 2008;11:195–202.CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    30.Taş N, Prestat E, Wang S, Wu Y, Ulrich C, Kneafsey T, et al. Landscape topography structures the soil microbiome in Arctic polygonal tundra. Nat Commun. 2018;9:777.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    31.Reasoner DJ, Geldreich EE. A new medium for the enumeration and subculture of bacteria from potable water. Appl Environ Microbiol. 1985;49:1–7.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    32.Muyzer G, De Waal EC, Uitterlinden AG. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol. 1993;59:695–700.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    33.Mizrahi-Man O, Davenport ER, Gilad Y. Taxonomic classification of bacterial 16S rRNA genes using short sequencing reads: evaluation of effective study designs. PLoS One. 2013;8:e53608.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    34.Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10.Article 

    Google Scholar 
    35.Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    36.Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    37.Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:2584.Article 

    Google Scholar 
    38.Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    39.Hammer DAT, Ryan PD, Hammer Ø, Harper DAT. Past: paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001;4:1–9.
    Google Scholar 
    40.Li D, Luo R, Liu C-M, Leung C-M, Ting H-F, Sadakane K, et al. MEGAHIT v1.0: a fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    41.Krawczyk PS, Lipinski L, Dziembowski A. PlasFlow: predicting plasmid sequences in metagenomic data using genome signatures. Nucleic Acids Res. 2018;46:e35.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    42.Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    43.Kahlke T, Ralph PJ. BASTA—taxonomic classification of sequences and sequence bins using last common ancestor estimations. Methods Ecol Evol. 2019;10:100–3.Article 

    Google Scholar 
    44.Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    45.Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    46.Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Third Int AAAI Conf Weblogs Soc Media. San Jose, California, USA. 2009.47.Jacomy M, Venturini T, Heymann S, Bastian M. ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software. PLoS One. 2014;9:e98679.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    48.McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57:3348–57.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    49.Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DG. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014;42:D737–43.CAS 
    PubMed 
    Article 

    Google Scholar 
    50.Caswell TA, Droettboom M, Hunter J, Lee A, Firing E, Stansby D, et al. matplotlib/matplotlib: REL: v3.1.1. 2019.51.Pham VHT, Kim J. Improvement for isolation of soil bacteria by using common culture media. J Pure Appl Microbiol. 2016;10:49–60.
    Google Scholar 
    52.Romaniuk K, Ciok A, Decewicz P, Uhrynowski W, Budzik K, Nieckarz M, et al. Insight into heavy metal resistome of soil psychrotolerant bacteria originating from King George Island (Antarctica). Polar Biol. 2018;41:1319–33.Article 

    Google Scholar 
    53.Amann RI, Ludwig W, Schleifer KH. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev. 1995;59:143–69.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    54.Tamaki H, Sekiguchi Y, Hanada S, Nakamura K, Nomura N, Matsumura M, et al. Comparative analysis of bacterial diversity in freshwater sediment of a shallow eutrophic lake by molecular and improved cultivation-based techniques. Appl Environ Microbiol. 2005;71:2162–9.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    55.Stewart EJ. Growing unculturable bacteria. J Bacteriol. 2012;194:4151–60.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    56.Ganzert L, Jurgens G, Munster U, Wagner D. Methanogenic communities in permafrost-affected soils of the Laptev Sea coast, Siberian Arctic, characterized by 16S rRNA gene fingerprints. FEMS Microbiol Ecol. 2007;59:476–88.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    57.Bajerski F, Ganzert L, Mangelsdorf K, Padur L, Lipski A, Wagner D. Chryseobacterium frigidisoli sp. nov., a psychrotolerant species of the family Flavobacteriaceae isolated from sandy permafrost from a glacier forefield. Int J Syst Evol Microbiol. 2013;63:2666–71.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    58.Filippidou S, Wunderlin T, Junier T, Jeanneret N, Dorador C, Molina V, et al. A combination of extreme environmental conditions favor the prevalence of endospore-forming Firmicutes. Front Microbiol. 2016;7:1707.PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    59.Kuramae EE, Yergeau E, Wong LC, Pijl AS, Veen JA, Kowalchuk GA. Soil characteristics more strongly influence soil bacterial communities than land-use type. FEMS Microbiol Ecol. 2012;79:12–24.CAS 
    PubMed 
    Article 
    PubMed Central 

    Google Scholar 
    60.Filippidou S, Junier T, Wunderlin T, Lo CC, Li PE, Chain PS, et al. Under-detection of endospore-forming Firmicutes in metagenomic data. Comput Struct Biotechnol J. 2015;13:299–306.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    61.Dziewit L, Cegielski A, Romaniuk K, Uhrynowski W, Szych A, Niesiobedzki P, et al. Plasmid diversity in arctic strains of Psychrobacter spp. Extremophiles. 2013;17:433–44.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    62.Mindlin S, Petrenko A, Kurakov A, Beletsky A, Mardanov A, Petrova M. Resistance of permafrost and modern Acinetobacter lwoffiistrains to heavy metals and arsenic revealed by genome analysis. Biomed Res Int. 2016;2016:3970831.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    63.Moghadam MS, Albersmeier A, Winkler A, Cimmino L, Rise K, Hohmann-Marriott MF, et al. Isolation and genome sequencing of four Arctic marine Psychrobacter strains exhibiting multicopper oxidase activity. BMC Genomics. 2016;17:117.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    64.Pal C, Bengtsson-Palme J, Kristiansson E, Larsson DGJ. Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential. BMC Genomics. 2015;16:964.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    65.Carroll LM, Gaballa A, Guldimann C, Sullivan G, Henderson LO, Wiedmann M. Identification of novel mobilized colistin resistance gene mcr-9 in a multidrug-resistant, colistin-susceptible Salmonella enterica serotype Typhimurium isolate. MBio. 2019;10:e00853–19.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    66.Bleich A, Fox JG. The mammalian microbiome and its importance in laboratory animal research. ILAR J. 2015;56:153–8.CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    67.Grond K, Sandercock BK, Jumpponen A, Zeglin LH. The avian gut microbiota: community, physiology and function in wild birds. J Avian Biol. 2018;49:e01788.Article 

    Google Scholar 
    68.Anganova EV, Savchenkov MF, Stepanenko LA, Savilov ED. Microbiological monitoring of opportunistic Enterobacteriaceae of the Lena river. Gig Sanit. 2016;95:1124–8.69.Tignat-Perrier R, Dommergue A, Thollot A, Keuschnig C, Magand O, Vogel TM, et al. Global airborne microbial communities controlled by surrounding landscapes and wind conditions. Sci Rep. 2019;9:14441.PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    70.Mu C, Zhang F, Chen X, Ge S, Mu M, Jia L, et al. Carbon and mercury export from the Arctic rivers and response to permafrost degradation. Water Res. 2019;161:54–60.CAS 
    PubMed 
    Article 

    Google Scholar 
    71.Petrova M, Kurakov A, Shcherbatova N, Mindlin S. Genetic structure and biological properties of the first ancient multiresistance plasmid pKLH80 isolated from a permafrost bacterium. Microbiology. 2014;160:2253–63.CAS 
    PubMed 
    Article 

    Google Scholar 
    72.Afouda P, Dubourg G, Levasseur A, Fournier P-E, Delerce J, Mediannikov O, et al. Culturing ancient bacteria carrying resistance genes from permafrost and comparative genomics with modern isolates. Microorganisms. 2020;8:1522.CAS 
    PubMed Central 
    Article 
    PubMed 

    Google Scholar  More