More stories

  • in

    China economy-wide material flow account database from 1990 to 2020

    China economy-wide material flow identification: system boundary, processes, and materialsThe first step is to define an economy, i.e., the economic (rather than geographical) territory of a country in which the activities and transactions of producer and consumer units are resident. Additionally, the period is a total of thirty-one years, from 1990 to 2020, for the following reasons: (1) statistics before 1990 are of poor quality and are insufficient to allow us to conduct analyses; and (2) so far, statistics have just recently been updated to cover the year of 2020. Furthermore, the analytical framework (hereinafter referred to as China EW-MFA) is developed to explore material utilisation and its environmental consequences within China’s economy.The general structure of China EW-MFA is depicted in Fig. 1, which comprises seven processes. (1) Input of extracted resources: domestic natural resources are extracted from the environment to the economy through human-controlled means. (2) Output of domestic processed materials: after being processed by manufacturers, materials are released from the economy into the environment in the form of by-products and residues, which can be classified by their destinations (i.e., air, land, and water) and pathways (dissipative use and losses). (3) Input and (4) output by cross-border trade: by imports and exports, materials are transported between China’s economy and the economies of the rest of the world. (5) Input and (6) output of balancing items (BI): sometimes, materials identified in the output processes are not considered by inputs, which needs to be balanced. For example, the utilisation of fossil energy materials by combustion causes the emission of carbon dioxide (CO2) into the air, which is identified as system output, but requirements of oxygen (O2) as system input are not counted. (7) Additions to the system: within the economy, materials would have been added to the economy in the form of buildings, infrastructures, durable goods, and household appliances, which are referred to as the net additions to stock (NAS).Fig. 1The general structure of China EW-MFA. To note, white data cells can be obtained directly from official statistics, whereas grey cells are estimated.Full size imageThe last step is to specify the materials concerned in each process. Four types (in blue boxes in Fig. 1) of natural materials are extracted and input into the economy in China, i.e., harvested biomass (33 items), mined metal ores (28 items), quarried non-metallic minerals (155 items), and mined fossil energy materials (6 items in 3 classes). Materials (green boxes) released into the air are greenhouse gases (e.g., CO2, methane (CH4), dinitrogen oxide (N2O)), air pollutants (e.g., particulate matter 10 (PM10), black carbon (BC)), and toxic contaminants of mercury (Hg) in divalent, gaseous elemental, and particulate forms. Those released into the water are inorganic matters (of nitrogen (N), phosphorus (P), Arsenic (As), and four heavy metals of lead (Pb), mercury (Hg), cadmium (Cd), and chromium (Cr)) and organic matters of cyanide, petroleum, and volatile phenol. Materials released into the land are waste disposal in uncontrolled landfills, which are illegal in China. Some materials are dissipated by application, for example, fertilisers, compost, sewage sludge being applied to agricultural land, and pesticides being used to cultivate crops. Some would be unintentionally dissipated from abrasion, corrosion, erosion, and leakages. Materials (in red boxes) are BI, which includes the input of O2 and output of water vapour in the fossil energy material combustion process, the input of O2 and output of water vapour and CO2 in the respiration process of human and cultivated livestock, input and output of water in imported and exported beverages, and the output of water from domestically extracting crops.There are some messages needed to be mentioned: (1) Material of water is not included since its flow volume is more substantial than others, which needs to be independently analysed; (2) Activities of foreign tourists, cross-border transfer of emissions through natural media, etc. are excluded. (3) To be clear, we refer to a data cell as a specific flow process of a specific substance in a specific year, e.g., the number of cereals domestically extracted in 2020.Data acquisition: sources and collectionBased on our China EW-MFA, we first analyse accessibility, reliability, completeness, rules of redistribution, etc., for each data source (yellow boxes in Fig. 1), including China national database, China rural statistical yearbooks, USGS mineral yearbooks, etc. The complete list of data sources and descriptions are presented in Table 1. Then, we store the originally retrieved data source files in a semi- or unstructured format (e.g., CSV, PDF). Next, we manually collect these statistics and reorganise them according to China EW-MFA material types and processes. However, only a tiny part of retrieved statistics can be applied directly, as specified in black colour in Fig. 1.Table 1 Data sources and descriptions.Full size tableData compilation: parameter localisation and data estimationA few inconsistencies in statistics were noticed, which would result in data incompleteness. For example, the domestic extraction of vegetables has been accounted for and published since 1995, before which statistics are unavailable. The domestically harvested timber has been measured in the volume unit of cubic metres, which needs to be converted into the mass unit via density conversion factor. Therefore, acquired statistics have to be estimated, which are specified in grey colour in Fig. 1. The following section elaborates on each data cell’s estimation methods, localised parameters, references, etc. In our uploaded data files, the original statistics, data sources, and compilation methods (using formulas) are all implemented, as explained in the Data Records Section.

    The input of natural resources by domestic extraction

    Vegetables in crops: Statistics of vegetable production (WVegetables)16 during 1990–1994 are unavailable, which is estimated based on the relationship between the production yield (PYield) and areas (AVegetables), as shown in Eq. 1. Here, PYield is assumed to remain constant at 27.04 thousand tonnes per thousand hectares from 1990 to 1995, derived by dividing vegetable production (257,267 thousand tonnes) by areas (9,515 thousand hectares) in 1995.$${W}_{Vegetables}={P}_{Yield}times {A}_{Vegetables}$$
    (1)

    Nuts in crops: One of them is chestnuts. The chestnut production in 2020 is unavailable, which is assumed to be the same as in 2019.

    Crop residues in biomass residues: They are referred to as that harvested production of crops that do not reach the market to be sold but are instead employed as raw materials for commercial purposes such as energy generation and livestock husbandry. This number (Wcrop residues) can be calculated by first determining the number of crop residues available from primary crop production (Wcrop) and the harvest factor (Pharvest factor), and then using the recovery rate (Precovery rate) to determine the number of crop residues used by the economy, as shown in Eq. 2. These parameters have been localized by previous studies17,18, which are adopted in this study, i.e., wheat (1.1 for Pharvest factor and 0.463 for Precovery rate), maize (1.2, 0.463), rice (0.9, 0.463), sugar cane (0.5, 0.9), beetroots (0.7, 0.9), tuber (0.5, 0.463), pulse (1.2, 0.7), cotton (3.4, 0.463), fibre crops (1.8, 0.463), silkworm cocoons (1.8, 0.463), and oil-bearing crops (1.8, 0.463).$${W}_{cropresidues}={W}_{crop}times {P}_{harvestfactor}times {P}_{recoveryrate}$$
    (2)

    Roughage of grazed biomass and fodder crops in biomass residues: In China, the grazed biomass for roughage includes annual forage and perennial forage, whereas fodder crops comprise straw feed, processed straw feed, and all other fodder crops. However, information19 on grazed biomass production is only accessible from 2006 to 2018, whereas fodder crop statistics are only available from 2015 to 2017. Equation 3 and Eq. 4 can be used to estimate unavailable statistics. To note, we assume that China’s domestic roughage supply structure has remained unaltered, which has two meanings. The proportion of total domestic roughage production (WDomestic production) in requirement (WRoughage requirement) has remained constant, while the proportion (PSupply fraction) of grazed biomass and fodder crop in domestic roughage production has been unchanged. The requirement (WRoughage requirement) is determined by the quantity of livestock (QLivestock) and their annual feeding amount (PAnnual intake). PAnnual intake (in tonnes per head per year) has been localised for each type of livestock4, with 4.5 for live cattle and buffaloes, 0.5 for sheep and goats, 3.7 for horses, and 2.2 for mules and asses.$${W}_{Roughagerequirement}={Q}_{Livestock}times {P}_{Annualintake}$$
    (3)
    $${W}_{Domesticproduction}={W}_{Roughagerequirement}times {P}_{Supplyfraction}$$
    (4)

    Timber in wood: As illustrated in Eq. 5, wood production16 is reported in volume units of cubic metres (VTimber), which need to be converted into mass units (WTimber) via density (PDensity). The parameter PDensity is assumed to be 0.58 tonnes per cubic metre, calculated by averaging 0.52 for coniferous types and 0.64 for non-coniferous ones4.$${W}_{Timber}={V}_{Timber}times {P}_{Density}$$
    (5)

    Non-ferrous metals in metal ores: Non-ferrous metal statistics are derived from two sources. China statistics20 are measured in gross ore (WMetal ores in gross ore) but are only available from 1999 to 2017, whereas the USGS statistics21 cover the period of 1990 to 2020 but they are measured in metal or concentrate content (WMetal ores in other units). Therefore, USGS statistics need to be converted with an empirical unit conversion factor (PUnit conversion factor) before being applied to estimate unavailable statistics reported by China, as shown in Eq. 6. Conversion factors are localised for each non-ferrous metal in each year from 2000 to 2017 by using USGS statistics divided by China statistics and then averaged after removing the highest value and the lowest value (i.e., trimmed mean). This factor could capture the general relationship between statistics from two separate sources, which can be used in other long time-series studies on resource management on a particular element in China.$${W}_{Metaloresingrossore}={W}_{Metaloresinotherunits}/{P}_{Unitconversionfactor}$$
    (6)

    Non-metallic minerals: The official China-specific information on non-metallic mineral domestic production is available between 1999 and 201720, the rest of which could be estimated from USGS statistics (1990–2020)21. Also, two differences in reporting standards are observed resulting from the material coverages and reporting units. China statistics contain eighty-eight materials in mineral ores, whereas the USGS only includes twenty in the concentrate unit. Therefore, a conversion factor is developed in this estimation, as shown in Eq. 7. This conversion factor is applied to the total amount of non-metallic mineral production, which is assumed to have been constant from 1990 to 1999 at 11.38% (1999) and 12.56% (2017) from 2017 to 2020.$${W}_{Mineralsingrossore}={W}_{Mineralsinotherunits}/{P}_{Conversionfactor}$$
    (7)

    Coal in fossil energy materials: Coal, mined in China, includes raw coal, peat, stone coal, and oil shale. Except for raw coal, statistics for the rest are only available from 1999 to 201720. The unavailable data (WOther coals) is estimated using Eq. 8 under the assumption that the structure of the coal supply in China barely changes. That is, the proportion (PSupply fraction) of peat, stone coal, and oil shale in raw coal production (WRaw coal) remains constant, so the 1999 proportion is applied to all years before that (earlier years of 1990–1998), while the 2017 proportion is used to the recent years between 2018 and 2020. For example, PSupply fraction for oil shale production was assumed to be 0.014% during 1990–1999, calculated by dividing raw coal production (1,250,000) by oil shale production (179) in 1999. PSupply fraction in the earlier and the recent years are 0.007% and 0.001% for peat, 0.203% and 0.031% for stone coal, and 0.014% and 0.067% for oil shale.

    $${W}_{Othercoals}={W}_{Rawcoal}/{P}_{Supplyfraction}$$
    (8)

    The output of processed materials by release

    Materials released into the air: In China, thirteen materials are released into the air, as shown in Fig. 1. The emission of sulphur dioxide (SO2) is reported in China environmental statistical yearbooks22,23, while the rest is specified in the EDGAR24. However, in EDGAR, statistics for recent years have not yet been updated, which are estimated with the value in the most recent year in our database. For example, nitrous oxide (NOx) records are only available for the years prior to 2016, with 26,365 thousand tonnes in 2015 and 26,837 in 2014. As a result of the observed decreasing trend in NOx emissions, NOx emission data for 2016–2020 is estimated to be 26,000 thousand tonnes. This estimate may be subjective due to constraints, but it would be aligned with European statistics, allowing for international comparisons. Data can be updated after the EDGAR statistics have been updated.

    Materials released into the water: Ten principal materials have been found in China wastewater (both industrial and municipal) that are nitrogen (N), phosphorus (P), organic pollutants of petroleum, volatile phenol and cyanide, heavy metals of mercury (Hg), lead (Pb), cadmium (C·d), and the hexavalent chromium (Cr6+), and arsenic (As). Many statistics22,23 have been of poor quality (e.g., inconsistent material coverages between years). Given that the statistics of pollutants in industrial wastewater cover more periods and contain fewer abnormal observations, the total material emissions can be approximated from those of industrial wastewater. Equations 9 and 10 show the estimation processes. The materials in industrial wastewater (WIndustrial materials) are first identified using material mass concentration (PConcentration) and the weight of industrial wastewater (WIndustrial wastewater), and then the materials in total wastewater (WTotal materials) are identified using the proportion (PContribution) of materials in industrial wastewaters (WIndustrial materials) to the total. The assumption is that PConcentration and PContribution change gradually between years, which enables to use linear interpolation method to estimate unavailable parameters. Consider cyanide: its PConcentration was 23.61 (1‰ ppm) in 2005 and 37.31 in 2002, which was assumed to be 28.18 in 2004 and 32.74 in 2003. PConcentration was assumed to be 100% throughout the years for cyanide because all cyanide emissions in China are driven by industrial wastewater discharges. Later, the total material emissions can be derived by dividing the industrial wastewater mass by PConcentration.$${W}_{Industrialmaterials}={W}_{Industrialwastewater}times {P}_{Concentration}$$
    (9)
    $${W}_{Totalmaterials}={W}_{Industrialmaterials},/,{P}_{Contribution}$$
    (10)

    Materials released to the land: This is zero because uncontrolled landfills are illegal in China.

    Materials dissipated by organic fertiliser use: In China, manure is the primary organic fertiliser, which is excreted by pigs, dairy cows, calves, sheep, horses, asses, mules, camels, chickens, and other animals. As shown in Eq. 11, the manure production (WManure) is estimated through the amounts of raised livestock (QLivestock, heads), the weight of daily manure production (PManure production, kilograms per head per day), the number of days they are raised (PFeeding period, in days per year), and the moisture content of their manure (PDry matter, %) for each type of animal. These parameters are region-specific, which have been localised by Chinese scholars25,26,27 and listed in Table 2.$${W}_{Manure}={Q}_{Livestock}times {P}_{Manureproduction}times {P}_{Feedingperiod}times {P}_{Drymatter}$$
    (11)
    Table 2 Localised parameters for animal manure production.Full size table

    Materials dissipated by mineral fertiliser use: The mineral fertilisers used in China are four types, i.e., nitrogen (N), phosphorus (P), potash (K), and compound. Their usage (WFertiliser usage) is measured in nutrient mass (WNutrient materials), which needs to be converted into the gross mass by dividing their nutrient content (PNutrient content). Equation 12 shows the estimation. This parameter of PNutrient content is localised by the Ministry of Agriculture and Rural Affairs of China28 as 29%, 22%, 35%, and 44% for N- bearing, P- bearing, K-bearing, and compound fertilisers, respectively.$${W}_{Fertiliserusage}={W}_{Nutrientmaterials}/{P}_{Nutrientcontent}$$
    (12)

    Materials dissipated by sewage sludge: Sewage sludge is the residue generated by municipal wastewater treatment. As demonstrated in Eq. 13, its dissipative use (Wss, dissipation) is the untreated amount of production (Wss, production), represented by the parameter of Pss, dissipation rate. Sewage sludge production (Wss, production) statistics are only available for the years 2006–202029, and data for the remaining years can be estimated using Eq. 14 and Eq. 15. In Eq. 14, Pss, production rate represents the relationship between sewage sludge production (Wss, production, 2006–2020) and wastewater treatment (Www, treatment, 2002–2020), and in Eq. 15, Pww, treatment efficiency represents the relationship between the quantity of treated wastewater (Www, treatment, 2002–2020) and the treatment capacity (Www, treatment capacity, 1990–2020). In this estimation, three assumptions are made. The first is to estimate Www, treatment, Pww, treatment efficiency is assumed to be unchanged at 63% during 1990–2001, given it has been increasing from 63% in 2002 to ~80% in recent years. The second is that, in order to estimate Wss, production, Pss, production rate is assumed to be unchanged at 3.5 between 1990 and 2005, suggesting 3.5 tonnes of sewage sludge are generated by processing 10,000 cubic metres of wastewater. This assumption is determined by that Pss, production rate is approximately 3.5 during 2006–2010 while declines sharply and stabilises at around two during 2011–2020. The last is, to estimate the Wss,dissipation, Pss,dissipation rate is assumed to be 5% between 1990 and 2005, given it has been around 5% during 2006–2020.$${W}_{ss,dissipation}={W}_{ss,production}times {P}_{ss,dissipationrate}$$
    (13)
    $${W}_{ss,production}={W}_{ww,treatment}times {P}_{ss,productionrate}$$
    (14)
    $${W}_{ww,treatment}={W}_{ww,treatmentcapacity}times {P}_{ww,treatmentefficiency}$$
    (15)

    Materials dissipated by composting: Composting is a natural process that uses microbes to turn organic materials into other products, which are then used for fertilising and entering the environment. In China, composting has been used to treat two materials: feces and municipal waste, whose quantities (WComposting) were only available from 2003 to 201029. The unavailable data can be estimated using Eq. 16. The dry weight of materials treated by composting (WComposting) is proportionally related to the fresh weight of all treated materials (WTotal), the proportion treated by composting (PComposting rate), and the dry content (PDry matter). Considering that China’s composting capacity has been declining since 2001 due to the implementation of waste incineration power generation technologies30, Pcomposting rate is assumed to be the same as it was in 2003 (9.5%) between 1990 and 2002, and 1.5% in 2010 between 2011 and 2020. The parameter of PDry matter is 50%4.$${W}_{Composting}={W}_{Total}times {P}_{Compostingrate}times {P}_{Drymatter}$$
    (16)

    The input and output by cross-border trade. Statistics of imports and exports have been gathered since 1962 and stored in the UN Comtrade database31. However, the data quality issue of outliers, and missing values, especially in weight, is reportedly identified. In our previous work, we addressed these issues, and an improved database32 is provided. Details about our estimation methods can be found in publications33,34,35. As UN Comtrade lists 5,039 different commodity types (in 6-digit HS0 commodity code), yet only 18 material types are specified in the China EW-MFA, UN Comtrade statistics need to be aligned to the China EW-MFA framework. Therefore, we compared each commodity and each material type between them and established a correspondence table to map UN Comtrade commodity types onto our EW-MFA material types. For example, non-ferrous metal materials of China EW-MFA include commodities, such as copper ores and concentrates (260300 HS0 code), silver powder (710610), manganese, articles thereof, and waste or scrap (811100), etc., whereas biomass residues include cereal straw and husks (121300), lucerne meal and pellets (121410), and other fodder and forage products (121410). This correspondence table between HS0 and EW-MFA classification for imports and exports is provided in Supplementary File 1.

    The input of balancing items

    O2 required for combustion: In BI, requirements for materials can be abstracted as equalling exogenous demands minus intrinsic supplies (Eq. 17). Three parts (two demands and one supply) are considered for O2 requirements by the combustion process: (1) demanding exogenous oxygen to oxidise elements (e.g., carbon, sulphur, nitrogen, etc., except for hydrogen) released into the air, (2) demanding exogenous oxygen to oxidise the hydrogen embedded in fossil energy materials, and (3) providing intrinsic oxygen embedded in fossil energy materials. The first part can be estimated via Eq. 18 by multiplying air emissions (WEmissions) of CO2, N2O, NOx, CO, and SO2 by their oxygen content (POxygen content). For the second (Eq. 19), the oxygen demand is estimated based on the principle of mass balance by converting the hydrogen amount of domestically utilised fossil energy materials (WFossil fuel materials × PHydrogen content) via molar mass conversion factor (PMass conversion factor). PMass conversion factor equals 7.92, derived by the molar mass of one oxygen (16 g/mol) divided by that of two hydrogen atoms (2 × 1.01 g/mol). The last is the intrinsic supplies from fossil fuel materials, which is identified via Eq. 20 by multiplying the domestically utilised amount of fossil fuel materials (WFossil fuel materials) by their oxygen content (POxygen content). The parameters in this estimation are presented in Table 3. As a footnote here, the domestically utilised amount is referred to as the domestic material consumption (DMC), which equals domestic extraction (DE) plus imports (IM) and minus exports (EX).$${W}_{Requirements}={W}_{Demands}-{W}_{Supplies}$$
    (17)
    $${W}_{Demands}={W}_{Emissions}times {P}_{Oxygencontent}$$
    (18)
    $${W}_{Demands}={W}_{Fossilfuelmaterials}times {P}_{Hydrogencontent}times {P}_{Massconversionfactor}$$
    (19)
    $${W}_{Supplies}={W}_{Fossilfuelmaterials}times {P}_{Oxygencontent}$$
    (20)
    Table 3 Parameters related to combustion processes4.Full size table

    O2 required for respiration: O2 is required by the metabolic activities of living organisms, the majority of which are humans and livestock. Bacteria are another sort of organism, which are not included in this estimation because their O2 requirements are too small to be quantified. The respiration-required O2 is related to the total quantity (QOrganisms) and their respiration activity by organism types, as shown in Eq. 21. The respiration activity is represented by the respiration requirement coefficient (PRespiration requirement coefficient), which is the average quantity of O2 that each organism utilises to maintain the metabolic activity, as listed in Table 4.$${W}_{Demands}={Q}_{Organisms}times {P}_{Respirationrequirementcoefficient}$$
    (21)
    Table 4 Parameters related to respiration processes4.Full size table

    Water required for the domestic production of exported beverages: The exported beverages are produced domestically using domestically extracted materials, especially a large amount of water. The weight of water is considered in the output by cross-border trade but is not included in the domestic extraction input. The resulted imbalance can be identified by specifying the water weight in beverages, i.e., multiplying the traded beverage weight (WMaterials) by a parameter of the water content (PWater content), as given in Eq. 22. Fruit and vegetable juices (2009 in HS0 code) and beverages (code 22) are covered in the improved UN Comtrade database32, with PWater content of 85% for the first and 90% for the latter4.

    $${W}_{Water}={W}_{Materials}times {P}_{Watercontent}$$
    (22)

    The output of balancing items.

    Water vapour from combustion: Water vapour emissions by domestically combusting fossil fuel materials are contributed by two paths. The direct evaporation of embedded water is the first path (Eq. 23), which can be derived by multiplying the DMC of fossil fuel materials by their moisture content (PMoisture content). The PMoisture content for each type of fossil fuel material is listed in Table 3. The other is the generation of water vapour during hydrogen oxidation, which can be calculated by converting the oxidised weight of hydrogen to the water weight using the molar mass conversion factor (PMass conversion factor), as given in Eq. 24. PMass conversion factor equals 8.92 by dividing the molar mass of water (18.02 g/mol) by that of two hydrogen atoms (2 × 1.01 g/mol).$${W}_{Water}={W}_{Fossilfuelmaterials}times {P}_{Moisturecontent}$$
    (23)
    $${W}_{Water}={W}_{Fossilfuelmaterials}times {P}_{Hydrogencontent}times {P}_{Massconversionfactor}$$
    (24)

    Water vapour and CO2 from respiration: Respiration activities of organisms will produce water vapour and CO2, whose estimation is similar to that of O2 requirements. As shown in Eq. 25, the respiration-caused gas emissions are related to the number of organisms (QOrganisms) and the respiration activity by organism types. The latter is represented by the parameter of respiration emission coefficient (PRespiration emission coefficient), which is specified in Table 4 for water vapour and CO2 for each type of organism.$${W}_{Emissions}={Q}_{Organisms}times {P}_{Respirationemissioncoefficient}$$
    (25)

    Water from imported beverages: The estimation approach is the same as water by the domestic production of exported beverages, as described in Eq. 16.

    Water in biomass products: Usually, the input of biomass products by domestic extraction16 has been measured in fresh weight, but their corresponding output29 by sewage sludge, composting, etc., are in dry weight, leading to an imbalance in water weight. The water weight in biomass products is calculated by multiplying their domestic extraction amount in fresh weight (WBiomass) by a parameter of moisture content at harvest (PMoisture content), as shown in Eq. 26. The values of PMoisture content by biomass products are presented in Table 5.Table 5 The moisture content at harvest for each biomass product4.Full size table

    $${W}_{Water}={W}_{Biomass}times {P}_{Moisturecontent}$$
    (26)
    Material flow quantificationThe above attempts have quantified material inputs and outputs by flows and presented a detailed profile of material utilisation for each material in China’s economy. In order to depict the economy in a more general way, EW-MFA indicators are assessed by aggregating flows by materials or periods as below.

    Domestic extraction (DE): is referred to as natural materials that are extracted from the domestic environment and are used in the domestic economy, i.e., the total input of natural materials by extraction.

    Domestic processed output (DPO): is referred to as materials that are released to the domestic environment after being processed in the domestic economy, i.e., the total output of processed materials by release.

    Import (IM): is referred to as all goods (in the form of raw materials, semi-finished materials, and final products) that originated from other economies and are further used in the domestic economy. It is calculated as the sum of all imported goods.

    Export (EX): is referred to as all goods that originated from the domestic economy and are transported to other economies to be used. It is calculated as the sum of all exported goods.

    Domestic material input (DMI): is referred to as materials that originated from the domestic environment by extraction and other economies and are available (to be used or to be stored) for the domestic economy. It is calculated as the sum of DE plus IM, as shown in Eq. 27.$$DMI=DE+IM$$
    (27)

    Domestic material consumption (DMC): is referred to as materials that are directly used in the domestic economy after parts of them are exported to other economies. It is calculated as the difference between DMI and EX.

    Physical trade balance (PTB): is referred to as a surplus or deficit of materials for the domestic economy. It is calculated as the difference between IM and EX.

    Net additions to stock (NAS): is referred to as materials that remain in the domestic economy. It is calculated by taking BI items into account, as shown in Eq. 28.

    $$NAS=DMC+B{I}_{in}-DPO-B{I}_{out}$$
    (28) More

  • in

    Nitrogen use aggravates bacterial diversity and network complexity responses to temperature

    Hwang, H. Y. et al. Effect of cover cropping on the net global warming potential of rice paddy soil. Geoderma 292, 49–58 (2017).ADS 
    CAS 
    Article 

    Google Scholar 
    IPCC. Climate change 2013: The physical science basis. The Working Group I contribution to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC, 2013).
    Google Scholar 
    Cardoso, R. M., Soares, P. M. M., Lima, D. C. A. & Miranda, P. M. A. Mean and extreme temperatures in warming climate: EURO CORDEX and WRF regional climate high-resolution projection for Portugal. Clim. Dyn. 52, 129–157 (2019).Article 

    Google Scholar 
    Ding, T., Gao, H. & Li, W. J. Extreme high-temperature event in southern China in 2016 and the possible role of cross-equatorial flows. Int. J. Climatol. 38, 3579–3594 (2018).Article 

    Google Scholar 
    Escalas, A. et al. Functional diversity and redundancy across fish gut, sediment, and water bacterial communities. Environ. Microbiol. 19, 3268–3282 (2017).Article 

    Google Scholar 
    Philippot, L. et al. Loss in microbial diversity affects nitrogen cycling in soil. ISME J. 7, 1609–1619 (2013).CAS 
    Article 

    Google Scholar 
    Li, Y. B. et al. Serratia spp. Are responsible for nitrogen fixation fueled by As(III) oxidation, a novel biogeochemical process identified in mine tailings. Environ. Sci. Technol 56, 2033–2043 (2022).ADS 
    Article 

    Google Scholar 
    Jia, M., Gao, Z. W., Gu, H. J., Zhao, C. Y. & Han, G. D. Effects of precipitation change and nitrogen addition on the composition, diversity, and molecular ecological network of soil bacterial communities in a desert steppe. PLoS ONE 16, e0248194. https://doi.org/10.1371/journal.pone.0248194 (2021).CAS 
    Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Waghmode, T. R. et al. Response of nitrifier and denitrifier abundance and microbial community structure to experimental warming in an agricultural ecosystem. Front. Microbiol. 9, 474. https://doi.org/10.3389/fmicb.2018.00474 (2018).Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Hu, Y. L., Wang, S., Niu, B., Chen, Q. & Zhang, G. Effect of increasing precipitation and warming on microbial community in Tibetan alpine steppe. Environ. Res. 189, 109917. https://doi.org/10.1016/j.envres.2020.109917 (2020).CAS 
    Article 
    PubMed 

    Google Scholar 
    Li, H. et al. Responses of soil bacterial communities to nitrogen deposition and precipitation increment are closely linked with aboveground community variation. Microb. Ecol. 71, 974–989 (2016).CAS 
    Article 

    Google Scholar 
    Wang, H. et al. Experimental warming reduced topsoil carbon content and increased soil bacterial diversity in a subtropical planted forest. Soil Biol. Biochem. 133, 155–164 (2019).CAS 
    Article 

    Google Scholar 
    Haumann, F. A., Gruber, N. & Münnich, M. Sea-Ice Induced Southern Ocean Subsurface Warming and Surface Cooling in a Warming Climate. AGU Advances 1, e2019AV000132. https://doi.org/10.1029/2019AV000132 (2020).ADS 
    Article 

    Google Scholar 
    Ji, F., Wu, Z. H., Huang, J. P. & Chassignet, E. P. Evolution of land surface air temperature trend. Nat. Clim. Chang. 4, 462–466 (2014).ADS 
    Article 

    Google Scholar 
    Sabri, N. S. A., Zakaria, Z., Mohamad, S. E., Jaafar, A. B. & Hara, H. Importance of soil temperature for the growth of temperate crops under a tropical climate and functional role of soil microbial diversity. Microbes Environ. 33, 144–150 (2018).Article 

    Google Scholar 
    McGrady-Steed, J. & Morin, P. T. Biodiversity, density compensation, and the dynamics of populations and functional groups. Ecology 81, 361–373 (2000).Article 

    Google Scholar 
    Jiang, L. Density compensation can cause no effect of biodiversity on ecosystem function. Oikos 116, 324–334 (2007).Article 

    Google Scholar 
    Faust, K. & Raes, J. Microbial interactions: From networks to models. Nat. Rev. Microbiol. 10, 538. https://doi.org/10.1038/nrmicro283 (2012).CAS 
    Article 
    PubMed 

    Google Scholar 
    Gao, X. X. et al. Revegetation significantly increased the bacterial-fungal interactions in different successional stages of alpine grasslands on the Qinghai-Tibetan Plateau. CATENA 205, 105385. https://doi.org/10.1016/j.catena.2021.105385 (2021).CAS 
    Article 

    Google Scholar 
    Morriën, E. et al. Soil networks become more connected and take up more carbon as nature restoration progresses. Nat. Commun. 8, 14349. https://doi.org/10.1038/ncomms14349 (2017).ADS 
    CAS 
    Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Banerjee, S. et al. Agricultural intensification reduces microbial network complexity and the abundance of keystone taxa in roots. ISME J. 13, 1722–1736 (2019).Article 

    Google Scholar 
    Pržulj, N. & Malod-Dognin, N. Network analytics in the age of big data. Science 353, 123–124 (2016).ADS 
    Article 

    Google Scholar 
    Ratzke, C., Barrere, J. M. R. & Gore, J. Strength of species interactions determines biodiversity and stability in microbial communities. Nat. Ecol. Evol. 4, 376–383 (2020).Article 

    Google Scholar 
    Fuhrman, J. A. Microbial community structure and its functional implications. Nature 45, 193–199 (2009).ADS 
    Article 

    Google Scholar 
    Zhao, M. X., Cong, J., Cheng, J. M., Qi, Q. & Zhang, Y. G. Soil microbial community assembly and interactions are constrained by nitrogen and phosphorus in broadleaf forests of southern China. Forest 11, 285. https://doi.org/10.3390/f11030285 (2020).Article 

    Google Scholar 
    Wan, X. L. et al. Biogeographic patterns of microbial association networks in paddy soil within Eastern China. Soil Biol. Biochem. 142, 07696. https://doi.org/10.1016/j.soilbio.2019.107696 (2020).CAS 
    Article 

    Google Scholar 
    Yuan, M. M., Guo, X., Wu, L., Zhang, Y. & Zhou, J. Climate warming enhances microbial network complexity and stability. Nat. Clim. Change 11, 343–348 (2021).ADS 
    Article 

    Google Scholar 
    Lassaletta, L. et al. Food and feed trade as a driver in the global nitrogen cycle: 50-year trends. Biogeochemistry 11, 225–241 (2014).Article 

    Google Scholar 
    Phoenix, G. K. et al. Impacts of atmospheric nitrogen deposition: Responses of multiple plant and soil parameters across contrasting ecosystems in long-term field experiments. Glob. Change Biol. 18, 1197–1215 (2012).ADS 
    Article 

    Google Scholar 
    Nakaji, T., Fukami, M., Dokiya, Y. & Izuta, T. Effects of high nitrogen load on growth, photosynthesis and nutrient status of Cryptomeria japonica and Pinus densiflora seedlings. Trees-Struct. Funct. 15, 453–461 (2001).CAS 
    Article 

    Google Scholar 
    Wang, H. Y. et al. Reduction in nitrogen fertilizer use results in increased rice yields and improved environmental protection. Int. J. Agric. Sustain. 15, 681–692 (2017).Article 

    Google Scholar 
    Zhou, X. G. & Wu, F. Z. Land-use conversion from open field to greenhouse cultivation differently affected the diversities and assembly processes of soil abundant and rare fungal communities. Sci. Total Environ. 788, 147751. https://doi.org/10.1016/j.scitotenv.2021.147751 (2021).ADS 
    CAS 
    Article 
    PubMed 

    Google Scholar 
    Guo, H. et al. Long-term nitrogen & phosphorus additions reduce soil microbial respiration but increase its temperature sensitivity in a Tibetan alpine meadow. Soil Biol. Biochem. 113, 26–34 (2017).CAS 
    Article 

    Google Scholar 
    Zhang, C. et al. Effects of simulated nitrogen deposition on soil respiration components and their temperature sensitivities in a semiarid grassland. Soil Biol. Biochem. 75, 113–123 (2014).CAS 
    Article 

    Google Scholar 
    Zhang, J. J. et al. Different responses of soil respiration and its components to nitrogen and phosphorus addition in a subtropical secondary forest. For. Ecosyst. 8, 37. https://doi.org/10.1186/s40663-021-00313-z (2021).Article 

    Google Scholar 
    Norse, D. & Ju, X. T. Environmental costs of China’s food security. Agric. Ecosyst. Environ. 209, 5–14 (2015).Article 

    Google Scholar 
    Xu, H. F., Du, H., Zeng, F. P., Song, T. Q. & Peng, W. X. Diminished rhizosphere and bulk soil microbial abundance and diversity across succession stages in Karst area, southwest China. Appl. Soil Ecol. 158, 103799. https://doi.org/10.1016/j.apsoil.2020.103799 (2020).Article 

    Google Scholar 
    Li, Y. B. et al. Arsenic and antimony co-contamination influences on soil microbial community composition and functions: Relevance to arsenic resistance and carbon, nitrogen, and sulfur cycling. Environ. Int. 153, 106522. https://doi.org/10.1016/j.envint.2021.106522 (2021).CAS 
    Article 
    PubMed 

    Google Scholar 
    Zhou, J. & Fong, J. J. Strong agricultural management effects on soil microbial community in a non-experimental agroecosystem. Appl. Soil Ecol. 165, 103970. https://doi.org/10.1016/j.apsoil.2021.103970 (2021).Article 

    Google Scholar 
    Bárcenas-Moreno, G., Gómez-Brandón, M., Rousk, J. & Bååth, E. Adaptation of soil microbial communities to temperature: Comparison of fungi and bacteria in a laboratory experiment. Glob. Chang. Biol. 15, 2950–2957 (2009).ADS 
    Article 

    Google Scholar 
    Tan, E. H., Zou, W., Zheng, Z., Yan, X. & Kao, S. J. Warming stimulates sediment denitrification at the expense of anaerobic ammonium oxidation. Nat. Clim. Change 10, 349–355 (2020).ADS 
    CAS 
    Article 

    Google Scholar 
    Supramaniam, Y., Chong, C. W., Silvaraj, S. & Tan, K. P. Effect of short term variation in temperature and water content on the bacterial community in a tropical soil. Appl Soil Ecol. 107, 279–289 (2016).Article 

    Google Scholar 
    Zhu, Y. Z., Li, Y. Y., Zheng, N. G., Chapman, S. J. & Yao, H. Y. Similar but not identical resuscitation trajectories of the soil microbial community based on either DNA or RNA after flooding. Agronomy 10, 502. https://doi.org/10.3390/agronomy10040502 (2020).CAS 
    Article 

    Google Scholar 
    Donhauser, J., Qi, W., Bergk-Pinto, B. & Frey, B. High temperatures enhance the microbial genetic potential to recycle C and N from necromass in high-mountain soils. Glob. Chang. Biol. 27, 1365–1386 (2021).ADS 
    Article 

    Google Scholar 
    Santoyo, G., Hernandez-Pacheco, C., Hernandez-Salmeron, J. & Hernandez-Leon, R. The role of abiotic factors modulating the plant-microbe-soil interactions: Toward sustainable agriculture. A review. Span. J. Agric. Res. 15, e03R01-e11. https://doi.org/10.5424/sjar/2017151-9990 (2017).Article 

    Google Scholar 
    Lefcheck, J. S. et al. Biodiversity enhances ecosystem multifunctionality across trophic levels and habitats. Nat. Commun. 6, 6936. https://doi.org/10.1038/ncomms7936 (2015).ADS 
    CAS 
    Article 
    PubMed 

    Google Scholar 
    Cardinale, B. J. et al. Corrigendum: Biodiversity loss and its impact on humanity. Nature 486, 59–67 (2012).ADS 
    CAS 
    Article 

    Google Scholar 
    Ma, B., Wang, H., Dsouza, M., Lou, J. & Xu, J. Geographic patterns of co-occurrence network topological features for soil microbiota at continental scale in eastern China. ISME J. 10, 1891–1901 (2016).CAS 
    Article 

    Google Scholar 
    Trivedi, C. et al. Losses in microbial functional diversity reduce the rate of key soil processes. Soil Biol. Biochem. 135, 267–274 (2019).CAS 
    Article 

    Google Scholar 
    Melanie, K. et al. Effects of season and experimental warming on the bacterial community in a temperate mountain forest soil assessed by 16S rRNA gene pyrosequencing. FEMS Microbiol. Ecol. 82, 551–562 (2012).Article 

    Google Scholar 
    Zheng, H. F., Liu, Y., Chen, Y., Zhang, J. & Chen, Q. Short-term warming shifts microbial nutrient limitation without changing the bacterial community structure in an alpine timberline of the eastern Tibetan Plateau. Geoderma 360, 113985. https://doi.org/10.1016/j.geoderma.2019.113985 (2020).ADS 
    CAS 
    Article 

    Google Scholar 
    Finlay, B. J. & Cooper, J. L. Microbial diversity and ecosystem function. CEH Integrating Fund second progress report to the Director, Centre for Ecology and Hydrology Nov 1996–Sept (1997).Xing, X. Y. et al. Warming shapes nirS- and nosZ-type denitrifier communities and stimulates N2O emission in acidic paddy soil. Appl. Environ. Microbiol. 87, e02965-e3020. https://doi.org/10.1128/AEM.0296520 (2021).CAS 
    Article 
    PubMed Central 

    Google Scholar 
    Lin, Y. T., Whitman, W. B., Coleman, D. C., Jien, S. H. & Chiu, C. Y. Soil bacterial communities at the treeline in subtropical alpine areas. CATENA 201, 105205. https://doi.org/10.1016/j.catena.2021.105205 (2021).CAS 
    Article 

    Google Scholar 
    Wang, J. C. et al. Impacts of inorganic and organic fertilization treatments on bacterial and fungal communities in a paddy soil. Appl. Soil Ecol. 112, 42–50 (2017).Article 

    Google Scholar 
    Chacón, J. M., Shaw, A. K. & Harcombe, W. R. Increasing growth rate slows adaptation when genotypes compete for diffusing resources. PLoS Comput. Biol. 16, e1007585. https://doi.org/10.1371/journal.pcbi.1007585 (2020).ADS 
    CAS 
    Article 
    PubMed 
    PubMed Central 

    Google Scholar 
    Hartley, I. P., Hopkins, D. W., Garnett, M. H., Sommerkorn, M. & Wookey, P. A. Soil microbial respiration in arctic soil does not acclimate to temperature. Ecol. Lett. 11, 1092–1100 (2008).Article 

    Google Scholar 
    Baath, E. Growth rates of bacterial communities in soils at varying pH: A comparison of the thymidine and leucine incorporation techniques. Microb. Ecol. 36, 316–327 (1998).CAS 
    Article 

    Google Scholar 
    Qin, H. L. et al. Soil moisture and activity of nitrite- and nitrous oxide-reducing microbes enhanced nitrous oxide emissions in fallow paddy soils. Biol. Fertil. Soils 56, 53–67 (2020).CAS 
    Article 

    Google Scholar 
    Chen, Z. et al. Impact of long term fertilization on the composition of denitrifier communities based on nitrite reductase analyses in a paddy soil. Microb. Ecol. 60, 850–861 (2010).CAS 
    Article 

    Google Scholar 
    Wei, G. S. et al. Similar drivers but different effects lead to distinct ecological patterns of soil bacterial and archaeal communities. Soil Biol. Biochem. 144, 107759. https://doi.org/10.1016/j.soilbio.2020.107759 (2020).CAS 
    Article 

    Google Scholar 
    Bastian, F., Bouziri, L., Nicolardot, B. & Ranjard, A. L. Impact of wheat straw decomposition on successional patterns of soil microbial community structure. Soil Biol. Biochem. 41, 262–275 (2009).CAS 
    Article 

    Google Scholar 
    Levins, R. Evolution in Changing Environments: Some Theoretical Explorations (Princeton University Press, 1968).Book 

    Google Scholar  More

  • in

    A dataset of road-killed vertebrates collected via citizen science from 2014–2020

    IPBES. Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services. Zenodo https://doi.org/10.5281/zenodo.5657041 (2019).Laurance, W. F. et al. A global strategy for road building. Nature 513, 229–232 (2014).ADS 
    CAS 
    Article 

    Google Scholar 
    Ibisch, P. L. et al. A global map of roadless areas and their conservation status. Science 354, 1423–1427 (2016).ADS 
    CAS 
    Article 

    Google Scholar 
    Forman, R. T. T., Sperling, D. & Bissonette, J. A. Road Ecology: Science and Solutions. (Island Pr, 2003).van der Ree, R., Smith, D. J. & Grilo, C. Handbook of Road Ecology. (John Wiley & Sons, 2015).Laender. Hunting Statistics. Game casualties 2017/2018: furred game (red deer, roe deer, chamois, moufflon) https://www.statistik.at/web_en/statistics/Economy/agriculture_and_forestry/livestock_animal_production/hunting/index.html (2018).Steiner, W., Leisch, F. & Hacklander, K. A review on the temporal pattern of deer-vehicle accidents: Impact of seasonal, diurnal and lunar effects in cervids. Accident; analysis and prevention 66, (2014).Kioko, J. et al. Driver knowledge and attitudes on animal vehicle collisions in Northern Tanzania. TROPICAL CONSERVATION SCIENCE 8, 352–366 (2015).Article 

    Google Scholar 
    Bíl, M., Andrášik, R. & Janoška, Z. Identification of hazardous road locations of traffic accidents by means of kernel density estimation and cluster significance evaluation. Accident Analysis & Prevention 55, 265–273 (2013).Article 

    Google Scholar 
    Page, Y. A statistical model to compare road mortality in OECD countries. Accident Analysis and Prevention 33, 371–385 (2001).CAS 
    Article 

    Google Scholar 
    Teixeira, F. Z. et al. Are Road-kill Hotspots Coincident among Different Vertebrate Groups? Oecologia Australis 17, 36–47 (2017).Article 

    Google Scholar 
    Canova, L. & Balestrieri, A. Long-term monitoring by roadkill counts of mammal populations living in intensively cultivated landscapes. Biodivers Conserv https://doi.org/10.1007/s10531-018-1638-3 (2018).Brehme, C. S., Hathaway, S. A. & Fisher, R. N. An objective road risk assessment method for multiple species: ranking 166 reptiles and amphibians in California. Landscape Ecol 33, 911–935 (2018).Article 

    Google Scholar 
    Heigl, F. et al. Comparing Road-Kill Datasets from Hunters and Citizen Scientists in a Landscape Context. Remote Sensing 8, (2016).Heigl, F., Horvath, K., Laaha, G. & Zaller, J. G. Amphibian and reptile road-kills on tertiary roads in relation to landscape structure: using a citizen science approach with open-access land cover data. BMC Ecol 17, 24 (2017).Article 

    Google Scholar 
    Dörler, D. & Heigl, F. A decrease in reports on road-killed animals based on citizen science during COVID-19 lockdown. PeerJ 9, e12464 (2021).Article 

    Google Scholar 
    Peer, M. et al. Predicting spring migration of two European amphibian species with plant phenology using citizen science data. Sci Rep 11, 21611 (2021).ADS 
    CAS 
    Article 

    Google Scholar 
    Schwartz, A. L. W. UK Roadkill Records. The Global Biodiversity Information Facility https://doi.org/10.15468/r3xakd (2018).Lin, T. The Taiwan Roadkill Observation Network Data Set. Version 1.3. The Global Biodiversity Information Facility https://doi.org/10.15468/cidkqi (2018).Chandler, M. et al. Contribution of citizen science towards international biodiversity monitoring. Biological Conservation 213, 280–294 (2017).Article 

    Google Scholar 
    Périquet, S., Roxburgh, L., le Roux, A. & Collinson, W. J. Testing the Value of Citizen Science for Roadkill Studies: A Case Study from South Africa. Front. Ecol. Evol. 6, (2018).Abra, F. D., Huijser, M. P., Pereira, C. S. & Ferraz, K. M. P. M. B. How reliable are your data? Verifying species identification of road-killed mammals recorded by road maintenance personnel in São Paulo State, Brazil. Biological Conservation 225, 42–52 (2018).Article 

    Google Scholar 
    Bíl, M., Kubeček, J., Sedoník, J. & Andrášik, R. Srazenazver.cz: A system for evidence of animal-vehicle collisions along transportation networks. Biological Conservation 213, 167–174 (2017). Part A.Article 

    Google Scholar 
    Vercayie, D. & Herremans, M. Citizen science and smartphones take roadkill monitoring to the next level. Nature Conservation 11, 29–40 (2015).Article 

    Google Scholar 
    Waetjen, D. P. & Shilling, F. M. Large Extent Volunteer Roadkill and Wildlife Observation Systems as Sources of Reliable Data. Front. Ecol. Evol. 5, (2017).Shilling, F. M., Perkins, S. E. & Collinson, W. Wildlife/Roadkill Observation and Reporting Systems. in Handbook of Road Ecology 492–501 (John Wiley & Sons, 2015).Eitzel, M. V. et al. Citizen Science Terminology Matters: Exploring Key Terms. Citizen Science: Theory and Practice 2, 1–20 (2017).
    Google Scholar 
    Haklay, M. et al. Contours of citizen science: a vignette study. Royal Society Open Science 8, 202108 (2021).ADS 
    Article 

    Google Scholar 
    Heigl, F., Kieslinger, B., Paul, K. T., Uhlik, J. & Dörler, D. Opinion: Toward an international definition of citizen science. PNAS 116, 8089–8092 (2019).CAS 
    Article 

    Google Scholar 
    Heigl, F. et al. Quality Criteria for Citizen Science Projects on Österreich forscht | Version 1.1. Open Science Framework https://doi.org/10.17605/OSF.IO/48J27 (2018).Heigl, F. et al. Co-Creating and Implementing Quality Criteria for Citizen Science. Citizen Science: Theory and Practice 5, 23 (2020).
    Google Scholar 
    Heigl, F. & Zaller, J. G. Using a Citizen Science Approach in Higher Education: a Case Study reporting Roadkills in Austria. Human Computation 1, (2014).University of Natural Resources and Life Sciences, Vienna. Roadkill, The Global Biodiversity Information Facility, https://doi.org/10.15468/ejb47y (2021).Heigl, F. & Roadkill Community. Roadkill Dataset 2014-2020 Quality level 2, Zenodo, https://doi.org/10.5281/zenodo.5878813 (2022).August, T. A. et al. Citizen meets social science: predicting volunteer involvement in a global freshwater monitoring experiment. Freshwater Science 38, 321–331 (2019).Article 

    Google Scholar 
    IUCN. The IUCN Red List of Threatened Species. Version 2021-3. IUCN Red List of Threatened Species https://www.iucnredlist.org/en (2021). More

  • in

    Disentangling influence over group speed and direction reveals multiple patterns of influence in moving meerkat groups

    Strandburg-Peshkin, A., Papageorgiou, D., Crofoot, M. C. & Farine, D. R. Inferring influence and leadership in moving animal groups. Philos. Trans. R. Soc. Lond. B Biol. Sci. 373(1746), 20170006 (2018).PubMed 
    PubMed Central 

    Google Scholar 
    Garland, J., Berdahl, A. M., Sun, J. & Bollt, E. M. Anatomy of leadership in collective behaviour. Chaos 28(7), 075308 (2018).ADS 
    MathSciNet 
    PubMed 

    Google Scholar 
    King, A. J., Douglas, C. M. S., Huchard, E., Isaac, N. J. B. & Cowlishaw, G. Dominance and affiliation mediate despotism in a social primate. Curr. Biol. 18(23), 1833–1838 (2008).CAS 
    PubMed 

    Google Scholar 
    Lewis, J. S., Wartzok, D. & Heithaus, M. R. Highly dynamic fission–fusion species can exhibit leadership when traveling. Behav. Ecol. Sociobiol. 65(5), 1061–1069 (2011).
    Google Scholar 
    Van Belle, S., Estrada, A. & Garber, P. A. Collective group movement and leadership in wild black howler monkeys (Alouatta pigra). Behav. Ecol. Sociobiol. 67(1), 31–41 (2013).
    Google Scholar 
    Smith, J. E. et al. Collective movements, leadership and consensus costs at reunions in spotted hyaenas. Anim. Behav. 105, 187–200 (2015).
    Google Scholar 
    Kerth, G., Ebert, C. & Schmidtke, C. Group decision making in fission–fusion societies: Evidence from two-field experiments in Bechstein’s bats. Proc. R. Soc. B Biol. Sci. 273(1602), 2785–2790 (2006).
    Google Scholar 
    Nagy, M., Ákos, Z., Biro, D. & Vicsek, T. Hierarchical group dynamics in pigeon flocks. Nature 464(7290), 890–893 (2010).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Giuggioli, L., McKetterick, T. J. & Holderied, M. Delayed response and biosonar perception explain movement coordination in trawling bats. PLoS Comput. Biol. 11(3), e1004089 (2015).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Pettit, B., Ákos, Z., Vicsek, T. & Biro, D. Speed determines leadership and leadership determines learning during pigeon flocking. Curr. Biol. 25(23), 3132–3137 (2015).CAS 
    PubMed 

    Google Scholar 
    Strandburg-Peshkin, A., Farine, D. R., Couzin, I. D. & Crofoot, M. C. Group decisions. Shared decision-making drives collective movement in wild baboons. Science 348(6241), 1358–1361 (2015).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Tokuyama, N. & Furuichi, T. Leadership of old females in collective departures in wild bonobos (Pan paniscus) at Wamba. Behav. Ecol. Sociobiol. 71(3), 55 (2017).
    Google Scholar 
    Montanari, D., O’Hearn, W. J., Hambuckers, J., Fischer, J. & Zinner, D. Coordination during group departures and progressions in the tolerant multi-level society of wild Guinea baboons (Papio papio). Sci. Rep. 11(1), 21938 (2021).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Papageorgiou, D. & Farine, D. R. Shared decision-making allows subordinates to lead when dominants monopolize resources. Sci. Adv. 6(48), 5881 (2020).ADS 

    Google Scholar 
    Bousquet, C. A. H., Sumpter, D. J. T. & Manser, M. B. Moving calls: A vocal mechanism underlying quorum decisions in cohesive groups. Proc. R. Soc. Lond. B Biol. Sci. 278(1711), 1482–1488 (2011).
    Google Scholar 
    Stahl, J., Tolsma, P. H., Loonen, M. J. J. E. & Drent, R. H. Subordinates explore but dominants profit: Resource competition in high Arctic barnacle goose flocks. Anim. Behav. 61(1), 257–264 (2001).PubMed 

    Google Scholar 
    Boinski, S. Social manipulation within and between troops mediate primate group movement. In On the Move: How and Why Animals Travel in Groups (ed. Boinski, S.) (University of Chicago Press, 2000).
    Google Scholar 
    Conradt, L. & Roper, T. J. Consensus decision making in animals. Trends Ecol. Evol. 20(8), 449–456 (2005).PubMed 

    Google Scholar 
    Conradt, L. & Roper, T. J. Conflicts of interest and the evolution of decision sharing. Philos. Trans. R. Soc. Lond. B Biol. Sci. 364(1518), 807–819 (2009).PubMed 

    Google Scholar 
    Byrne, R. W. How monkeys find their way: Leadership, coordination, and cognitive maps of African baboons. In On the Move: How and Why Animals Travel in Groups (eds Boinski, S. & Garber, P. A.) (University of Chicago Press, 2000).
    Google Scholar 
    Conradt, L. & Roper, T. J. Deciding group movements: Where and when to go. Behav. Proc. 84, 675–677 (2010).
    Google Scholar 
    Herbert-Read, J. E. et al. Inferring the rules of interaction of shoaling fish. PNAS 108(46), 18726–18731 (2011).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Katz, Y., Tunstrøm, K., Ioannou, C. C., Huepe, C. & Couzin, I. D. Inferring the structure and dynamics of interactions in schooling fish. PNAS 108(46), 18720–18725 (2011).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Jolles, J. W., Boogert, N. J., Sridhar, V. H., Couzin, I. D. & Manica, A. Consistent individual differences drive collective behavior and group functioning of schooling fish. Curr. Biol. 27(18), 2862–2868 (2017).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Doolan, S. P. & Macdonald, D. W. Breeding and juvenile survival among slender-tailed meerkats (Suricatu suricatta) in the south-western Kalahari: Ecological and social influences. J. Zool. 242(2), 309–327 (1997).
    Google Scholar 
    Clutton-Brock, T. H. & Manser, M. B. Meerkats: Cooperative breeding in the Kalahari. In Cooperative Breeding in Vertebrates (eds Koenig, W. D. & Dickinson, J. L.) (Cambridge University Press, 2016).
    Google Scholar 
    Doolan, S. & Macdonald, D. Diet and foraging behaviour of group living meerkats, Suricata suricatta, in the southern Kalahari. J. Zool. 239, 697–716 (1996).
    Google Scholar 
    Engesser, S. Function of ‘Close’ Calls in a Group Foraging Carnivore, Suricata suricatta (2011).Kranstauber, B., Gall, G. E. C., Vink, T., Clutton-Brock, T. & Manser, M. B. Long-term movements and home-range changes: Rapid territory shifts in meerkats. J. Anim. Ecol. https://doi.org/10.1111/1365-2656.13129 (2019).Article 
    PubMed 

    Google Scholar 
    Manser, M. B. et al. Vocal Complexity in Meerkats and Other Mongoose Species Vol. 46, 281 (Elsevier, 2014).
    Google Scholar 
    Gall, G. E. C. & Manser, M. B. Group cohesion in foraging meerkats: Follow the moving ‘vocal hot spot’. R. Soc. Open Sci. 4, 170004 (2017).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Engesser, S. & Manser, M. B. Collective close calling mediates group cohesion in foraging meerkats via spatially determined differences in call rates. Anim. Behav. 185, 73–82 (2022).
    Google Scholar 
    Gall, G. E. C., Strandburg-Peshkin, A., Clutton-brock, T. & Manser, M. B. As dusk falls: Collective decisions about the return to sleeping sites in meerkats. Anim. Behav. 132, 91–99 (2017).
    Google Scholar 
    Townsend, S. W., Rasmussen, M., Clutton-Brock, T. & Manser, M. B. Flexible alarm calling in meerkats: The role of the social environment and predation urgency. Behav. Ecol. 23(6), 1360–1364 (2012).
    Google Scholar 
    Clutton-Brock, T. H. et al. Contributions to cooperative rearing in meerkats. Anim. Behav. 61(4), 705–710 (2001).
    Google Scholar 
    Griffin, A. S. et al. A genetic analysis of breeding success in the cooperative meerkat (Suricata suricatta). Behav. Ecol. 14(4), 472–480 (2003).
    Google Scholar 
    Thavarajah, N. K., Fenkes, M. & Clutton-Brock, T. H. The determinants of dominance relationships among subordinate females in the cooperatively breeding meerkat. Behaviour 151(1), 89–102 (2014).
    Google Scholar 
    Young, A. J. et al. Stress and the suppression of subordinate reproduction in cooperatively breeding meerkats. Proc. Natl. Acad. Sci. 103(32), 12005–12010 (2006).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Hodge, S. J., Manica, A., Flower, T. P. & Clutton-Brock, T. H. Determinants of reproductive success in dominant female meerkats. J. Anim. Ecol. 77(1), 92–102 (2008).PubMed 

    Google Scholar 
    Bell, M. B. V. et al. Suppressing subordinate reproduction provides benefits to dominants in cooperative societies of meerkats. Nat. Commun. 22(5), 4499 (2014).ADS 

    Google Scholar 
    Bousquet, C. A. H. & Manser, M. B. Resolution of experimentally induced symmetrical conflicts of interest in meerkats. Anim. Behav. 81(6), 1101–1107 (2011).
    Google Scholar 
    Strandburg-Peshkin, A., Clutton-Brock, T. & Manser, M. B. Burrow usage patterns and decision-making in meerkat groups. Behav. Ecol. 31(2), 292–302 (2020).
    Google Scholar 
    Turbé, A. Foraging Decisions and Space Use in a Social Mammal, The Meerkat—Chapter 6: Leadership pby Lactating Female in Meerkats (University of Cambridge, 2006).
    Google Scholar 
    Barelli, C., Reichard, U., Boesch, C. & Heistermann, M. Female white-handed gibbons (Hylobates lar) lead group movements and have priority of access to food resources. Behaviour 145(7), 965–981 (2008).
    Google Scholar 
    Clutton-Brock, T. H. et al. Reproduction and survival of suricates (Suricata suricatta) in the southern Kalahari. Afr. J. Ecol. 37(1), 69–80 (1999).
    Google Scholar 
    Kutsukake, N. & Clutton-Brock, T. H. Do meerkats engage in conflict management following aggression? Reconciliation, submission and avoidance. Anim. Behav. 75(4), 1441–1453 (2008).
    Google Scholar 
    Spong, G. F., Hodge, S. J., Young, A. J. & Clutton-Brock, T. H. Factors affecting the reproductive success of dominant male meerkats: Reproductive success in male meerkats. Mol. Ecol. 17(9), 2287–2299 (2008).PubMed 

    Google Scholar 
    Russell, A. F., Carlson, A. A., McIlrath, G. M., Jordan, N. R. & Clutton-Brock, T. Adaptive size modification in dominant female meerkats. Evolution 58(7), 1600–1607 (2004).PubMed 

    Google Scholar 
    R. Development Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/ (2008).Pinheiro, J. & Bates, D. M. Mixed-Effects Models in S and S-PLUS (Springer-Verlag, 2000).MATH 

    Google Scholar 
    Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biometr. J. 50(3), 346–363 (2008).MathSciNet 
    MATH 

    Google Scholar 
    Makowski, D., Ben-Shachar, M. S., Patil, I. & Lüdecke, D. Methods and algorithms for correlation analysis in R. J. Open Source Softw. 5(51), 2306 (2020).ADS 

    Google Scholar 
    Farine, D. R., Strandburg-Peshkin, A., Couzin, I. D., Berger-Wolf, T. Y. & Crofoot, M. C. Individual variation in local interaction rules can explain emergent patterns of spatial organization in wild baboons. Proc. R. Soc. B 284(1853), 20162243 (2017).PubMed 
    PubMed Central 

    Google Scholar 
    Holekamp, K. E., Boydston, E. E. & Smale, L. Group tarvel in social carnivores. In On the Move (eds Boinski, S. & Garber, P. A.) (University of Chicago Press, 2000).
    Google Scholar 
    Fischhoff, I. R. et al. Social relationships and reproductive state influence leadership roles in movements of plains zebra, Equus burchellii. Anim. Behav. 73(5), 825–831 (2007).
    Google Scholar 
    Furrer, R. D., Kunc, H. P. & Manser, M. B. Variable initiators of group departure in a cooperative breeder: The influence of sex, age, state and foraging success. Anim. Behav. 84(1), 205–212 (2012).
    Google Scholar 
    Clutton-Brock, T. H. et al. Costs of cooperative behaviour in suricates (Suricata suricatta). Proc. R. Soc. B Biol. Sci. 265(1392), 185–190 (1998).CAS 

    Google Scholar 
    MacLeod, K. J. & Clutton-Brock, T. H. Low costs of allonursing in meerkats: Mitigation by behavioral change? Behav. Ecol. 26(3), 697–705 (2015).
    Google Scholar 
    Boinski, S. The coordination of spatial position: A field study of the vocal behaviour of adult female squirrel monkeys. Anim. Behav. 41(1), 89–102 (1991).
    Google Scholar 
    Bode, N. W. F., Franks, D. W. & Wood, A. J. Leading from the front? Social networks in navigating groups. Behav. Ecol. Sociobiol. 66(6), 835–843 (2012).
    Google Scholar 
    Reber, S. A., Townsend, S. W. & Manser, M. B. Social monitoring via close calls in meerkats. Proc. R. Soc. B Biol. Sci. 280(1765), 20131013 (2013).
    Google Scholar 
    Bracken, A. M., Christensen, C., O’Riain, M. J., Fürtbauer, I. & King, A. J. Flexible group cohesion and coordination, but robust leader–follower roles, in a wild social primate using urban space. Proc. R. Soc. B Biol. Sci. 289(1967), 20212141 (2022).
    Google Scholar  More

  • in

    We must get a grip on forest science — before it’s too late

    Climate models need to capture a full spectrum of data from forests such as the Brazilian Amazon.Credit: Florence Goisnard/AFP/Getty

    Humanity’s understanding of how forests are responding to climate change is disconcertingly fragile. Take carbon fertilization, for example — the phenomenon by which plants absorb more carbon dioxide as its concentration in the atmosphere increases. This is one of the principal mechanisms by which nature has so far saved us from the worst of climate change, but there’s little understanding of its future trajectory. In fact, researchers don’t fully understand how climate change interacts with a multitude of forest processes. Complex, unsolved questions include how climate warming affects forest health; how it affects the performance of forests as a carbon sink; and whether it alters the ecosystem services that forests provide. Forests are our life-support system, and we should be more serious about taking their pulse.Six papers in this week’s Nature provide important insights into those questions. They also underline some of the challenges that must be overcome if we are to fully understand forests’ potential in the fight against climate change. These challenges are not only in the science itself, but also relate to how forest scientists collaborate, how they are funded (especially where data collection is concerned) and how they are trained.Forest science is an amalgam of disciplines. Ecologists and plant scientists measure tree growth, soil nutrients and other parameters in thousands of forest plots around the world. Physical scientists monitor factors such as forest height and above-ground forest biomass using remote-sensing data from drones or satellites. Experimental scientists investigate how forests might behave in a warming world by artificially altering factors such as temperature or carbon dioxide levels in experimental plots. Some of the data they generate are absorbed by yet another community: the modellers, who have created dynamic global vegetation models (DGVMs). These simulate how carbon and water cycles change with climate and, in turn, inform broader earth-system and climate models of the type that feed into policymaking.Different DGVMs make different predictions about how long forests will continue to absorb anthropogenic CO2. One reason for these differences is that models are sensitive to assumptions made about the processes in forests. There are many influences — including temperature, moisture, fire and nutrients — that are generally studied in isolation. Yet they interact with each other.Not all DGVMs account for the dampening effect that a lack of soil phosphorus can have on carbon fertilization, for example. Much of central and eastern Amazonia is poor in phosphorus, and research has shown that introducing phosphorus limitation into DGVMs can cut the carbon-fertilization effect1. This week, Hellen Fernanda Viana Cunha at the National Institute for Amazonian Research in Manaus, Brazil, and her colleagues report2 a powerful experimental demonstration of how the soil’s poor phosphorus content limits carbon absorption in an old-growth Amazonian forest.Models simulating the northward spread of boreal forest as temperatures rise are also missing key drivers3, according to Roman Dial at Alaska Pacific University in Anchorage and his colleagues. They report today that a white-spruce population has migrated surprisingly far north into the Arctic tundra. To explain this, it is necessary to take into account winter winds (which facilitate long-distance dispersal) along with the availability of deep snow and soil nutrients (which promote plant growth).Models are often based on a small number of ‘functional tree types’ — for example, ‘evergreen broadleaf’ or ‘evergreen needle leaf’. These are chosen as a proxy for the behaviour of the planet’s more than 60,000 known tree species. Yet ecologists are discovering that the biology of individual species matters when it comes to a tree’s response to climate change.David Bauman at the Environmental Change Institute at the University of Oxford, UK, and his co-workers reported in May that tree mortality on 24 moist tropical plots in northern Australia has doubled in the past 35 years (and life expectancy has halved), apparently owing to the increasing dryness of the air4. But that was an average of the 81 dominant tree species: mortality rates varied substantially between species, a variation that seemed to be related to the density of their wood.Peter Reich at the Institute for Global Change Biology at the University of Michigan in Ann Arbor and his colleagues now report that modest alterations in temperature and rainfall led to varying rates of growth and survival5 for different species in southern boreal-forest trees. The species that prospered were rare.Failure to examine multiple factors simultaneously means that scientists are making findings that challenge the assumptions in models. Spring is coming earlier for temperate forests and most models assume that, by prolonging the growing season, this increases woody-stem biomass. However, observational work carried out in temperate deciduous forests by Kristina Anderson-Teixeira at the Smithsonian Conservation Biology Institute in Front Royal, Virginia, and her colleagues found no sign of this happening6.Modellers are all too aware of the need to incorporate more complexity into their models, and of the potential that increasing amounts of computing power have to assist them in this endeavour. But they need more data.Continuity problemTo obtain comprehensive, valuable data for the models, continuous, long-term observations need to be made, and that depends on the availability of long-term funding. Achieving such continuity is a problem for both remote-sensing and ground-based operations. The former can cost hundreds of millions of dollars, but the value of its long-term data sets is immense, as demonstrated by a team led by Giovanni Forzieri at the University of Florence in Italy. The authors used 20 years of satellite data to show that nearly one-quarter of the world’s intact forests have already reached their critical threshold for abrupt decline7. But even field-based data collection, which costs a pittance by comparison, struggles to achieve financial security.Important ground-based operations include the Forest Global Earth Observatory (ForestGEO), part of the Smithsonian Tropical Research Institute, which is headquartered in Washington DC. This monitors 7.5 million individual trees in plots around the world. The amount of work that goes into this monitoring is formidable. For example, at present, ForestGEO is conducting the eighth five-yearly census of a plot in Peninsular Malaysia. This involves determining the species for each of the 350,000 trees (there are some 800 species growing there) and measuring the circumference of each trunk. It will take 16 skilled people a year to measure all the trees. Delays in the provision of funding to ForestGEO have held up similar censuses at plots in countries including Papua New Guinea, Vietnam, Brunei and Ecuador.

    A ForestGEO researcher making tree measurements at a forest plot in Barro Colorado Island, Panama.Credit: Jorge Aleman, STRI

    The future of the plots in North Queensland, which supplied Bauman with a rare 49 years’ worth of continuous data, is uncertain. They have been monitored since the mid-1970s by the Australian public research-funding agency CSIRO — initially every two years, then, more recently, every five years. In 2019, monitoring of the plots was switched to every 50 years because of funding shortages at CSIRO, leaving scientists searching for new sources of funding.Without continuity of funding, organizations such as ForestGEO can’t equip researchers with the requisite skills or collect data over periods longer than an individual’s time in a specific post or a funder’s cycle. “We have trained people and then lost them due to job insecurity,” says Stuart Davies, who leads ForestGEO.Different groups of forest researchers are trying to address these problems. ForestGEO is coordinating the Alliance for Tropical Forest Science in an effort to make it easier to share data, and to bolster the morale and careers of the skilled technicians and scientists — many of whom live in low- and middle-income countries — who do the bulk of the data collection.But we also need more-imaginative funding mechanisms that lift long-term observational plots out of three- to five-year funding cycles. Space agencies that fund remote-sensing satellites could collaborate with other funding agencies, for example, so that earth-observation missions include a fully funded component for ground-based data collection — which is, after all, crucial for calibrating their results. Journals, too, could do more to value and incentivise the production of long-term data sets.And there is a need for more interdisciplinarity. The US Department of Energy is funding a project called NGEE–Tropics (Next-Generation Ecosystem Experiments–Tropics) in which modellers will work with empirical researchers, both observational and experimental, who study tropical forests to create a full, process-rich model of such forests. This is encouraging, and the idea could be pushed further. What is needed is an initiative that pulls the disciplines together towards a goal of building a better understanding of forest processes. Among other things, such an initiative would encourage researchers in different disciplines to take each other’s data needs into account when planning their projects.For this to work, we need to remember that the edifice of forest science relies on the long-term data that scientists wring from forests over decades. Our chances of overcoming climate change are small, but they will diminish further if we forget the basics of monitoring our home planet. More

  • in

    Metaproteome plasticity sheds light on the ecology of the rumen microbiome and its connection to host traits

    Shotgun sequencing and generation of metagenome-assembled genomesIn our previous study, 78 Holstein Friesian dairy cows were sampled for rumen content, metagenomic shotgun sequencing was carried out, and raw Illumina sequencing reads were assembled into contigs using megahit assembler using default settings [7]. We used a pooled assembly of the original 78 samples to increase the quality of the metagenome-assembled genomes (MAGs) with the syntax: megahit [14] -t 60 -m 0.5 −1 [Illumina R1 files] −2 [Illumina R2 files]. Next, the assembled contigs were indexed using BBMap [15]: bbmap.sh threads = 60 ref = [contigs filename]. Thereafter, reads from each sample were mapped to the assembled contigs using BBTools’ bbwrap.sh script. In order to determine the depth (coverage) of each contig within each sample, the gi_summarize_bam_contig_depths tool was applied with the parameters: gi_summarize_bam_contig_depths –outputDepth depth.txt –pairedContigs paired.txt *.bam –outputDepth depth.txt –pairedContigs paired.txt.Using the depth information, metabat2 [16] was executed to bind genes together into reconstructed genomes, with parameters: metabat2 -t40 -a depth.txt.To evaluate genomic bin quality, we used the CheckM [17] tool, with parameters: checkm lineage_wf [in directory] [out directory] -x faa –genes -t10.Preparing proteomic search libraryWe generated 93 unique high-quality MAGs, and further increased our MAG database by including phyla that were not represented in our set of MAGs. In order to do so, we used the published compendium of 4,941 rumen metagenome-assembled genomes [18] and dereplicated those MAGs using dRep [19]. We then selected MAGs from phylum Spirochaetes, Actinomycetota, Proteobacteria, Firmicutes, Elusimicrobia, Bacillota, Fibrobacteres and Fusobacteria, which had the highest mean coverage in our samples as calculated using BBMap and gi_summarize_bam_contig_depths as described above [15]. This strategy minimized the false discovery rate (FDR), that would have been obtained if larger and unspecific databases would have been employed [20] and allowed the addition of 14 MAGs to our database.In order to create the proteomic search library, genes were identified along the 107 MAGs using the Prodigal tool [21], with parameters: prodigal meta and translated in silico into proteins, using the same tool. Replicates sequences were removed. Protein sequences from the hosting animal (Bos taurus) and common contaminant protein sequences (64,701 in total) were added to the proteomic search library in order to avoid erroneous target protein identification originating from the host or common contaminants. Finally, in order to subsequently assess the percentage of false-positive identifications within the proteomic search [22], the proteomic search library sequences were reversed in order and served as a decoy database.Proteomic analysisThe bacterial fraction from rumen fluid of the 12 selected animals selected from extreme feed efficiency phenotypes, were obtained at the same time as the samples analyzed for metagenomics and stored at −20 °C until extraction. To extract total proteins, a modified protocol from Deusch and Seifert was used [23]. Briefly, cell pellets were resuspended in 100 µl in 50 mM Tris-HCl (pH 7.5; 0.1 mg/ml chloramphenicol; 1 mM phenylmethylsulfonyl fluoride (PMSF)) and incubated for 10 min at 60 °C and 1200 rpm in a thermo-mixer after addition of 150 µl 20 mM Tris-HCl (pH 7.5; 2% sodium dodecyl sulfate (SDS)). After the addition of 500 µl DNAse buffer (20 mM Tris-HCl pH 7.5; 0.1 mg/ml MgCl2, 1 mM PMSF, 1 μg/ml DNAse I), the cells were lysed by ultra-sonication (amplitude 51–60%; cycle 0.5; 4 × 2 min) on ice, incubated in the thermo-mixer (10 min at 37 °C and 1,200 rpm) and centrifuged at 10,000 × g for 10 min at 4 °C. The supernatant was collected and centrifuged again. The proteins in the supernatant were precipitated by adding 20% pre-cooled trichloroacetic acid (TCA; 20% v/v). After centrifugation (12,000 × g; 30 min; 4 °C), the protein pellets were washed twice in pre-cooled (−20 °C) acetone (2 × 10 min; 12,000 × g; 4 °C) and dried by vacuum centrifugation. The protein pellet was resuspended in 2× SDS sample buffer (4% SDS (w/v); 20% glycerin (w/v); 100 mM Tris-HCl pH 6.8; a pinch of bromophenol blue, 3.6% 2‑mercaptoethanol (v/v)) by 5 min sonication bath and vortexing. Samples were incubated for 5 min at 95 °C and separated by 1D SDS-PAGE (Criterion TG 4-20% Precast Midi Gel, BIO-RAD Laboratories, Inc., USA).As previously described, after fixation and staining, each gel line was cut into 10 pieces, destained, desiccated, and rehydrated in trypsin [24]. The in-gel digest was performed by incubation overnight at 37 °C. Peptides were eluted with Aq. dest. by sonication for 15 min The sample volume was reduced in a vacuum centrifuge.Before MS analysis, the tryptic peptide mixture was loaded on an Easy-nLC II or Easy-nLC 1000 (Thermo Fisher Scientific, USA) system equipped with an in-house built 20 cm column (inner diameter 100 µm; outer diameter 360 µm) filled with ReproSil-Pur 120 C18-AQ reversed-phase material (3 µm particles, Dr. Maisch GmbH, Germany). Peptides were eluted with a nonlinear 156 min gradient from 1 to 99% solvent B (95% acetonitrile (v/v); 0.1% acetic acid (v/v)) in solvent A (0.1% acetic acid (v/v)) with a flow rate of 300 ml/min and injected online into an LTQ Orbitrap Velos or Orbitrap Velos Pro (Thermo Fisher Scientific, USA). Overview scan at a resolution of 30,000 in the Orbitrap in a range of 300-2,000 m/z was followed by 20 MS/MS fragment scans of the 20 most abundant precursor ions. Ions without detected charge state as well as singly charged ions were excluded from MS/MS analysis. Original raw spectra files were converted into the common mzXML format, in order to further process it in downstream analysis. The spectra file from each proteomic run of a given sample was searched against the protein search library, using the Comet [25] search engine with default settings.The TPP pipeline (Trans Proteomic Pipeline) [26] was used to further process the Comet [25, 27] search results and produce a protein abundance table for each sample. In detail, PeptideProphet [28] was applied to validate peptide assignments, with filtering criteria set to probability of 0.001, accurate mass binning, non-parametric errors model (decoy model) and decoy hits reporting. In addition, iProphet [28, 29] was applied to refine peptide identifications coming from PeptideProphet. Finally, ProteinProphet [28,29,30] was applied to statistically validate peptide identifications at the protein level. This was carried out using the command: xinteract -N[my_sample_nick].pep.xml -THREADS = 40 -p0.001 -l6 -PPM -OAPd -dREVERSE_ -ip [file1].pep.xml [file2].pep.xml.. [fileN].pep.xml  > xinteract.out 2  > xinteract.err. Then, TPP GUI was used in order to produce a protein table from the resulting ProtXML files (extension ipro.prot.xml).Subsequently, proteins that had an identification probability < 0.9 were also removed as well as proteins supported with less than 2 unique peptides (see Supplementary Table 1).Quantifying metagenomic presence of MAGsA reference database containing all 107 MAGs’ contigs was created (bbmap.sh command, default settings). Then, the paired-end short reads from each sample (FASTQ files) were mapped into the reference database (bbwrap.sh, default settings), producing alignment (SAM) files, which were converted into BAM format. Subsequently, a contig depth (coverage) table was produced using the command jgi_summarize_bam_contig_depths --outputDepth depth.txt --pairedContigs paired.txt *.bam. As each of the MAGs span on more than one contig, MAG depth in each sample was calculated as contig length weighted by the average depth. Finally, to account for unequal sequencing depth, each MAG depth was normalized to the number of short sequencing reads within the given sample.Correlating metagenomic and proteomic structuresIn order to compare metagenomic and proteomic structures, we first calculated the mean coding gene abundance and mean production levels of each of the 1629 detected core proteins over all 12 cows. Both mean gene abundance and mean production level were translated into ranks using the R rank function. The produced proteins were ranked in descending order and the coding genes in the gene abundance vector were reordered accordingly. The two reordered ranked vectors then plotted using the R pheatmap function, and colored using the same color scale.Selection of proteins for downstream analysisAs our goal was to analyze plasticity in microbial protein production in varying environments, e.g., as a function of host state, only MAGs that were identified in all of the 12 proteomic samples were kept for further analysis. Consequently, only proteins that were identified in at least half of the proteomic samples (e.g., in at least six samples) were selected. This last step aimed to reduce spurious correlation results. These filtering steps retained 79 MAGs coding for a total of 1,629 measurable proteins.Feed efficiency state prediction and ordinationIn order to calculate the accuracy in predicting host feed efficiency state based on the different data layers available (16S rRNA (Supplementary Table 2), metagenomics, metaproteomics), the principal component analysis (PCA) axes for all the samples based on the microbial protein production profiles were calculated. Then, twelve cycles of model building and prediction were made. Each time, the two first PCs of each of five cows along with their phenotype (efficiency state) were used to build a Support Vector Machine (SVM) [R caret package] prediction model and one sample was left out. The model was then used to perform subsequent prediction of the left-out animal phenotype (feed efficiency) by feeding the model with that animal’s first two PCs. This leave-one-out methodology was then repeated over all the samples. Finally, the prediction accuracy was determined as the percent of the cases where the correct label was assigned to the left-out sample. For the proteomics data, this procedure was applied on both the raw protein counts, and the protein production normalized based on MAG abundance, which enabled us to compare the prediction accuracies of the microbial protein production to that of the raw protein counts.Identification proteins associated with a specific host stateIn order to split the proteomics dataset into microbial proteins that tend to be produced differently as a function of the host feed efficiency states, each microbial protein profile was correlated to the sample’s host feed efficiency measure (as calculated by RFI) using the Spearman correlation (R function cor), disregarding the p value. Proteins that had a positive correlation to RFI were grouped as inefficiency associated proteins. In contrast, proteins that presented a negative correlation to RFI were grouped as efficiency associated proteins. To test for equal sizes of these two protein groups, a binomial test was performed (R function binom.test) to examine the probability to get a low number of feed efficient proteins from the overall proteins under examination, when the expected probability was set to 0.5.Functional assignment of proteinsProtein functions were assigned based on the KEGG (Kegg Encyclopedia of Genes and Genomes) [31] database. The entire KEGG genes database was compiled into a Diamond [32] search library. Then, the selected microbial proteins were searched against the database using the Diamond search tool. Significant hits (evalue < 5e-5) were further analyzed to identify the corresponding KO (KEGG Ortholog number). Annotations of glycoside hydrolases were performed using dbcan2 [33].Protein level checkerboard distribution across the feed efficiency groupsThe checkerboard distribution in protein production profiles was estimated separately within the feed efficient and inefficient animal groups. To enable the comparison between the two groups’ checkerboardness level, we chose a standardized C-score estimate (Standardized Effect Size C-score - S.E.S C-Score), based on the comparison of the observed C-score to a null-model distribution derived from simulations. The S.E.S C-score was estimated using the oecosimu function from R vegan package with 100,000 simulated null-model communities.Calculating functional redundancyThe functional redundancy within a given group of proteins was measured as the mean number of times a given KO occurred within a given group, while neglecting proteins that have not been assigned a KO level functional annotation.In order to test whether a given group of proteins exhibits more or less functional redundancy than would have been expected, a null distribution for functional redundancy was created, based on the number of proteins in the given group. A random group of proteins was drawn from the entire set, keeping the same sample size as in the tested group, and the process was repeated 100 times. Then, the functional redundancy for each random protein group was calculated. Thereafter, the null distribution was used to obtain a p value to measure the likelihood of obtaining such a value under the null.Examining functional divergenceExamining the functional divergence between the two groups of proteins, e.g. the feed efficiency and inefficiency associated proteins, was done by first counting the amount of shared functional annotations, in terms of KOs between the two groups. Thereafter, a null distribution for the expected count of KOs was built by randomly splitting in an iterative manner the proteins into groups of the same sizes and calculating the number of shared KOs. A p value for the actual count of shared proteins was obtained by ranking the actual count over the null distribution.Calculating average nearest neighbor ratio (ANN ratio)ANN Ratio analysis was carried out independently for each protein function (KO), containing more than 14 proteins with at least 5 proteins within each feed efficiency group. Initially, all proteins assigned to a given KO were split into two sets, in accordance to their feed efficiency affiliation group. Thereafter, proteins within each set were independently projected into two-dimensional space by PCA applied directly to Sequence Matrix [34]. Average nearest neighbor ratio within each set was then calculated within the minimum enclosing rectangle defined by principal component axes PC1 and PC2, as defined by Clark and Evans [35].MAG feed efficiency score calculationMicroorganism feed efficiency score was calculated for each MAG individually by first ranking each protein being produced by the given microbe along the 12 animals, based on the normalized protein production levels. Thereafter, a representative production value for the microbe in each animal was calculated as the average of the ranked (normalized) protein production levels in that animal (using R rank function). This ranking allowed us to alleviate the potential skewing effect of highly expressed proteins. The microorganism’s Feed Efficiency Score was calculated as the difference between its mean representative production value within feed efficient animals to that within feed inefficient animals. Values close to zero will reflect similar distribution between the two animal groups, positive values will indicate higher expression among efficient animals, and negative values will indicate higher expression among inefficient animals. To calculate significance, the actual feed efficiency score was compared to values in a distribution derived from a permutation based null model. Each of the permuted Feed Efficiency Scores (10,000 for each microbe) was obtained by independently shuffling each of the proteins produced by the MAG between the animals, prior to calculating the actual microorganism feed efficiency score. By positioning the absolute score value over its distribution under permuted assumptions (absolute values), we obtained a significance p value.MAG phylogenetic tree construction and phylogenetic signal estimationIn order to assess the link between phylogenetic similarity between the MAGs and their association with feed efficiency, phylogenetic tree estimating evolutionary relationships between the MAGs was constructed using the PhyloPhlAn pipeline [36]. The phylogenetic signal for Microorganism Feed Efficiency Score was estimated by providing the phylogSignal function from R phylosignal [37] package with MAGs phylogenetic tree and respective values. Pagel’s Lambda statistics was chosen for the analysis, owing to its robustness [38].Plot generationAll bar plots, scatter plots and other point plots were generated with R package ggplot2. Heatmaps were produced by either ggplot2 [39] or pheatmap [https://cran.r-project.org/web/packages/pheatmap/index.html] R packages. KEGG map was produced using the online KEGG Mapper tool [40]. Phylocorrelogram was produced with phyloCorrelogram function from R package phylosignal [37].MAG differential production analysisMAGs that contain a minimal number of proteins (50 functions) were selected for differential protein production analysis, in order to have sufficient data to perform statistical tests. For each MAG, the relative production was used in order to calculate the Jaccard pairwise dissimilarity for core protein production between feed efficient and inefficient cows using the R vegan package. Analysis of similarity between efficiency and inefficiency associated proteins for each MAG (ANOSIM) values and p values were then calculated using the same package.Predicting animal feed efficiency state according to GH family countsUsing all GH annotated proteins, a feature table that sums the count of each GH family within each sample was produced. Thereafter a leave-one-out cross-validation (LOOCV) [R caret package] was performed, each time building a Random Forest (RF) prediction model from the GH family counts and efficiency state of 11 samples, leaving one sample outside. Each one of the RF models, in its turn, was applied on the left-out animal to predict its efficiency state. Model accuracy and AUC curve were calculated based on the LOOCV performance. More

  • in

    Healthy foods, COVID rebound — the week in infographics

    Healthy foods are green foodsAn analysis of 57,000 foods reveals which have the best and worst environmental impacts. A team of researchers used an algorithm to estimate how much of each ingredient was in thousands of products sold in major UK supermarket chains. The scientists then gave food items an environmental-impact score out of 100 — with 100 being the worst — by combining the impacts of the ingredients in 100 grams of each product. They considered several factors, including greenhouse-gas emissions and land use.Healthier foods tended to have low environmental impacts, the team found. Products containing lamb and beef — such as ready-made meat pies — had the most serious environmental impact. The lowest-impact foods tended to be made with plants and included bread products, fruits, vegetables, grains and sugar-rich drinks. There were some notable exceptions: both nuts and seafood had a good nutrition score but relatively high environmental impacts.

    Source: M. Clark et al. Proc. Natl Acad. Sci. USA 119, e2120584119 (2022).

    COVID-19 reboundAfter the COVID-19 antiviral Paxlovid began to be used in late 2021, researchers noticed that some people experience a rebound in the virus and symptoms after taking the drug. Two recent studies suggest that it is surprisingly common for SARS-CoV-2 to return in untreated cases of COVID-19 too. To determine the frequency of rebound, Jonathan Li, a physician-scientist at Brigham and Women’s Hospital in Boston, Massachusetts, and his team analysed data from hundreds of people who were randomized to receive a placebo in a large-scale trial of COVID-19 antibody drugs. More than one-quarter of participants who were infected with SARS-CoV-2 reported a rebound in their symptoms, according to the study, which has not yet been peer reviewed. “The main take-home message is that recovery from COVID-19 is not going to be a linear process,” Li says.

    Source: Deo, R. et al. Preprint at medRxiv https://doi.org/10.1101/2022.08.01.22278278 (2022).

    How trees grow in a warmer worldThis chart shows how the early arrival of spring, due to climate change, affects the growth of trees and the amount of carbon they sequester. In a paper in Nature, researchers investigated the consequences of an early start to the growing season in deciduous forests. Leaf emergence is followed by carbon uptake through the process of photosynthesis. Over time, carbon can be captured for long-term sequestration if it contributes to the radial growth of stems or to wood formation. The areas under the curves represent annual growth in terms of: the amount of carbon captured by leaves (top curves, brown); annual radial growth (lower left curve, blue); and increase in woody biomass (lower right curve, red).The authors report that the early arrival of spring, which shifts the margins of the growing season (lighter curves), had little impact on the final annual tree-ring width or the amount of woody biomass produced, whereas high temperatures in summer had a negative effect on radial growth (dotted curve). Other studies (plotted here as dotted curves) indicate that high temperatures and related drought can suppress carbon capture and woody-biomass production.The study provides evidence that warmer springs have advanced the leaf emergence of temperate deciduous forests but have not substantially increased their wood production. This suggests that the extra uptake of carbon dioxide does not contribute to sustainable carbon sequestration in the trunks of long-lived trees, as our News & Views article explains. More

  • in

    The combination of genomic offset and niche modelling provides insights into climate change-driven vulnerability

    Wiens, J. J. Climate-related local extinctions are already widespread among plant and animal species. PLoS Biol. 14, e2001104 (2016).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Scheffers, B. R. et al. The broad footprint of climate change from genes to biomes to people. Science 354, aaf7671 (2016).PubMed 
    Article 
    CAS 

    Google Scholar 
    Parmesan, C. & Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 421, 37–42 (2003).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Urban, M. C. Accelerating extinction risk from climate change. Science 348, 571–573 (2015).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Wingfield, J. C. et al. Organism-environment interactions in a changing world: a mechanistic approach. J. Ornithol. 152, 279–288 (2011).Article 

    Google Scholar 
    Mendoza-Gonzalez, G., Martinez, M. L., Rojas-Soto, O. R., Vazquez, G. & Gallego-Fernandez, J. B. Ecological niche modeling of coastal dune plants and future potential distribution in response to climate change and sea level rise. Glob. Change Biol. 19, 2524–2535 (2013).ADS 
    Article 

    Google Scholar 
    Saunders, S. P. et al. Community science validates climate suitability projections from ecological niche modeling. Ecol. Appl. 30, 17 (2020).Article 

    Google Scholar 
    Peterson, A. T., Cobos, M. E. & Jimenez-Garcia, D. Major challenges for correlational ecological niche model projections to future climate conditions. Ann. N. Y. Acad. Sci. 1429, 66–77 (2018).ADS 
    PubMed 
    Article 

    Google Scholar 
    Mays, H. L. et al. Genomic analysis of demographic history and Ecological niche modeling in the endangered Sumatran Rhinoceros Dicerorhinus sumatrensis. Curr. Biol. 28, 70–76 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    Malcolm, R. J., Liu, C., Neilson, P. R., Hansen, L. & Hannah, L. A. Global warming and extinctions of endemic species from biodiversity hotspots. Conserv. Biol. 20, 538–548 (2005).Article 

    Google Scholar 
    Fitzpatrick, M. C. & Keller, S. R. Ecological genomics meets community-level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation. Ecol. Lett. 18, 1–16 (2015).PubMed 
    Article 

    Google Scholar 
    Gotelli, J. N. & Stanton-Geddes, J. Climate change, genetic markers and species distribution modelling. J. Biogeogr. 42, 1577–1585 (2015).Article 

    Google Scholar 
    Ruegg, K. et al. Ecological genomics predicts climate vulnerability in an endangered southwestern songbird. Ecol. Lett. 21, 1085–1096 (2018).PubMed 
    Article 

    Google Scholar 
    Razgour, O. et al. Considering adaptive genetic variation in climate change vulnerability assessment reduces species range loss projections. Proc. Natl Acad. Sci. USA 116, 10418–10423 (2019).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Valladares, F. et al. The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change. Ecol. Lett. 17, 1351–1364 (2014).PubMed 
    Article 

    Google Scholar 
    Bay, R. A. et al. Genomic signals of selection predict climate-driven population declines in a migratory bird. Science 359, 83–86 (2018).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Rhone, B. et al. Pearl millet genomic vulnerability to climate change in West Africa highlights the need for regional collaboration. Nat. Commun. 11, 5274 (2020).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Rahbek, C. et al. Building mountain biodiversity: geological and evolutionary processes. Science 365, 1114–1119 (2019).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar 
    Fjeldså, J., Bowie, R. C. K. & Rahbek, C. The role of mountain ranges in the diversification of birds. Annu. Rev. Ecol. Evol. Syst. 43, 249–265 (2012).Article 

    Google Scholar 
    Freeman, B. G., Scholer, M. N., Ruiz-Gutierrez, V. & Fitzpatrick, J. W. Climate change causes upslope shifts and mountaintop extirpations in a tropical bird community. Proc. Natl Acad. Sci. USA 115, 11982–11987 (2018).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    He, J. K., Lin, S. L., Li, J. T., Yu, J. H. & Jiang, H. S. Evolutionary history of zoogeographical regions surrounding the Tibetan Plateau. Commun. Biol. 3, 9 (2020).Article 
    CAS 

    Google Scholar 
    Wu, Y. J. et al. Explaining the species richness of birds along a subtropical elevational gradient in the Hengduan Mountains. J. Biogeogr. 40, 2310–2323 (2013).Article 

    Google Scholar 
    del Hoyo, J., Elliott, A., Sargatal, J. & Christie, D. A. Handbook of the Birds of the World (Lynx Edicions, 2013).Qu, Y. et al. Lineage diversification and historical demography of a montane bird Garrulax elliotii – implications for the Pleistocene evolutionary history of the eastern Himalayas. BMC Evolut. Biol. 11, 174 (2011).Article 

    Google Scholar 
    Qu, Y. et al. Long-term isolation and stability explain high genetic diversity in the Eastern Himalaya. Mol. Ecol. 23, 705–720 (2014).PubMed 
    Article 

    Google Scholar 
    Wang, W. J. et al. Glacial expansion and diversification of an East Asian montane bird, the green-backed tit (Parus monticolus). J. Biogeogr. 40, 1156–1169 (2013).Article 

    Google Scholar 
    Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212 (2015).PubMed 
    Article 
    CAS 

    Google Scholar 
    Laine, V. N. et al. Evolutionary signals of selection on cognition from the great tit genome and methylome. Nat. Commun. 7, 9 (2016).Article 
    CAS 

    Google Scholar 
    Ellis, N., Smith, S. J. & Pitcher, C. R. Gradient forests: calculating importance gradients on physical predictors. Ecology 93, 156–168 (2012).PubMed 
    Article 

    Google Scholar 
    Giorgetta, M. A. Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5. J. Adv. Model. Earth Syst. 5, 572–597 (2013).ADS 
    Article 

    Google Scholar 
    Gent, P. R. et al. The community climate system model version 4. J. Clim. 24, 4973–4991 (2011).ADS 
    Article 

    Google Scholar 
    Watanabe, M. et al. Improved climate simulation by MIROC5: mean states, variability, and climate sensitivity. J. Clim. 23, 6312–6335 (2010).ADS 
    Article 

    Google Scholar 
    Voldoire, A. et al. The CNRM-CM5.1 global climate model: description and basic evaluation. Clim. Dyn. 40, 2091–2121 (2013).Article 

    Google Scholar 
    Frichot, E., Schoville, S. D., Bouchard, G. & Francois, O. Testing for associations between loci and environmental gradients using latent factor mixed models. Mol. Biol. Evol. 30, 1687–1699 (2013).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Forester, B. R., Jones, M. R., Joost, S., Landguth, E. L. & Lasky, J. R. Detecting spatial genetic signatures of local adaptation in heterogeneous landscapes. Mol. Ecololgy 25, 104–120 (2016).CAS 
    Article 

    Google Scholar 
    Forester, B. R., Lasky, J. R., Wagner, H. H. & Urban, D. L. Comparing methods for detecting multilocus adaptation with multivariate genotype-environment associations. Mol. Ecol. 27, 2215–2233 (2018).CAS 
    PubMed 
    Article 

    Google Scholar 
    Li, C. et al. Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment. Gigascience 3, 27 (2014).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Pirri, F. et al. Selection-driven adaptation to the extreme Antarctic environment in Emperor penguin. Preprint at bioRxiv https://doi.org/10.1101/2021.12.14.471946 (2021).Wang, L. C. et al. Involvement of the Arabidopsis HIT1/AtVPS53 tethering protein homologuein the acclimation of the plasma membrane to heat stess.J. Exp. Bot. 62, 3609–3620 (2011).CAS 
    PubMed 
    Article 

    Google Scholar 
    Piñol, R. A. et al. Preoptic BRS3 neurons increase body temperature and heart rate via multiple pathways. Cell Metab. 33, 1389–1403 (2021).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Guilherme, A. et al. Neuronal modulation of brown adipose activity through perturbation of white adipocyte lipogenesis. Mol. Metab. 16, 116–125 (2018).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Zhang, Y., Guo, W., zhang, Y., Zhang, H. & Wu, C. Insights into hypoxic adaptation in Tibetan chicken embryos from comparative proteomics. Comp. Biochem. Physiol. Part D. 31, 100602 (2019).CAS 

    Google Scholar 
    Pizzagalli, M. D., Bensimon, A. & Superti-Furga, G. A guide to plasma membrane solute carrier proteins. FEBS J. 288, 2784–2835 (2020).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Qu, Y. et al. Rapid phenotypic evolution with shallow genomic differentiation during early stages of high elevation adaptation in Eurasian Tree Sparrows. Natl Sci. Rev. 7, 113–127 (2020).PubMed 
    Article 

    Google Scholar 
    Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Ferrier, S., Manion, G., Elith, J. & Richardson, K. Using generalized dissimilarity modelling to analyse and predict patterns of beta diversity in regional biodiversity assessment. Diversity Distrib. 13, 252–264 (2007).Article 

    Google Scholar 
    Thuiller, W., Lafourcade, B., Engler, R. & Araújo, M. B. BIOMOD – a platform for ensemble forecasting of species distributions. Ecography 32, 369–373 (2009).Article 

    Google Scholar 
    Chen, Y. et al. Large-scale genome-wide reveals climate adaptive variability in a cosmopolitan pest. Nat. Commun. 12, 7206 (2021).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Clarke, R. T., Rothery, P. & Raybould, A. F. Confidence limits for regression relationships between distance matrices: Estimating gene flow with distance. J. Agric. Biol. Environ. Stat. 7, 361–372 (2002).Article 

    Google Scholar 
    Excoffier, L., Dupanloup, I., Huerta-Sanchez, E., Sousa, V. C. & Foll, M. Robust demographic inference from genomic and SNP data. PLoS Genet. 9, e1003905 (2013).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Foden, W. B. et al. Climate change vulnerability assessment of species. WIREs Clim. Change 10, e551 (2019).Article 

    Google Scholar 
    Smith, T. B. et al. Genomic vulnerablity and soci-economic threats under climate change in an African rainforest bird. Evolut. Appl. 14, 1239–1247 (2021).Article 

    Google Scholar 
    Liu, B., Liang, E. Y., Liu, K. & Camarero, J. J. Species- and elevation-dependent growth responses to climate warming of mountain forests in the Qinling Mountains, central China. Forests 9, 11 (2018).
    Google Scholar 
    Dang, H. S., Zhang, Y. J., Zhang, K. R., Jiang, M. X. & Zhang, Q. F. Climate-growth relationships of subalpine fir (Abies fargesii) across the altitudinal range in the Shennongjia Mountains, central China. Clim. Change 117, 903–917 (2013).ADS 
    Article 

    Google Scholar 
    Lingua, E., Cherubini, P., Motta, R. & Nola, P. Spatial structure along an altitudinal gradient in the Italian central Alps suggests competition and facilitation among coniferous species. J. Veg. Sci. 19, 425–436 (2008).Article 

    Google Scholar 
    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 
    Zhang, R. Z., Zheng, D., Yang, Q. Y. & Liu, Y. H. Physical Geography of Hengduan Mountains (Science Press, 1997).Liu, Y. et al. Sino-Himalayan mountains act as cradles of diversity and immigration centres in the diversification of parrotbills (Paradoxornithidae). J. Biogeogr. 43, 1488–1501 (2016).Bush, A. et al. Incorporating evolutionary adaptation in species distribution modeling reduces projected vulnerability to climate change. Ecol. Lett. 17, 1468–148 (2016).Article 

    Google Scholar 
    Sparks, M. M., Westley, A. A. H., Falke, J. A. & Quinn, T. P. Thermal adaptation and phenotypic plasticity in a warming world: insights from common garden experiments on Alaskan sockeye salmon. Glob. Change Biol. 23, 5203–5217 (2017).ADS 
    Article 

    Google Scholar 
    Merow, C., Wilson, A. M. & Jetz, W. Integrating occurrence data and expert maps for improved species range predictions. Glob. Ecol. Biogeogr. 26, 243–258 (2017).Article 

    Google Scholar 
    Weisenfeld, N. I., Kumar, V., Shah, P., Church, D. M. & Jaffe, D. B. Direct determination of diploid genome sequences. Genome Res. 27, 757–767 (2017).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Luo, R. et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience 1, 18 (2012).PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    She, R., Chu, J. S. C., Wang, K., Pei, J. & Chen, N. GenBlastA: enabling BLAST to identify homologous gene sequences. Genome Res. 19, 143–149 (2009).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Birney, E., Clamp, M. & Durbin, R. GeneWise and genomewise. Genome Res. 14, 988–995 (2004).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    McKenna, A. et al. The genome analysis toolkit: a mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).PubMed 
    PubMed Central 
    Article 
    CAS 

    Google Scholar 
    Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).Article 

    Google Scholar 
    Robinson, J. D., Bunnefeld, L., Hearn, J., Stone, G. N. & Hickerson, M. J. ABC inference of multi-population divergence with admixture from unphased population genomic data. Mol. Ecol. 23, 4458–4471 (2014).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Nazareno, A. G., Bemmels, J. B., Dick, C. W. & Lohmann, L. G. Minimum sample sizes for population genomics: an empirical study from an Amazonian plant species. Mol. Ecol. Resour. 17, 1136–1147 (2017).CAS 
    PubMed 
    Article 

    Google Scholar 
    Willing, E. M., Dreyer, C. & van Oosterhout, C. Estimates of genetic differentiation measured by FST do not necessary require large sample size when using many SNP markers. PLoS One 7, e2649 (2012).Article 
    CAS 

    Google Scholar 
    Keenan, K., Mcginnity, P., Cross, T. F., Crozier, W. W. & Prodöhl, P. A. diveRsity: an Rpackage for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol. Evol. 4, 782–788 (2013).Article 

    Google Scholar 
    Rellstab, C., Gugerli, F., Eckert, I. A., Hancock, M. A. & Holderegger, R. A practical guide to environmental assocaition analysis in landscape genomics. Mol. Ecol. 24, 4348–4370 (2015).PubMed 
    Article 

    Google Scholar 
    Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92 (2012).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Xie, C. et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39, W316–W322 (2011).CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Araújo, M. B. & New, M. Ensemble forecasting of species distributions. Trends Ecol. Evol. 22, 42–47 (2007).PubMed 
    Article 

    Google Scholar 
    Araújo, M. B. et al. Standards for distribution models in biodiversity assessments. Sci. Adv. 5, eaat4858 (2019).ADS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Boria, R. A., Olson, L. E., Goodman, S. M. & Anderson, R. P. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecol. Model. 275, 73–77 (2014).Article 

    Google Scholar 
    Anderson, R. P. & Raza, A. The effect of the extent of the study region on GISmodels of species geographic distributions and estimates of niche evolution: preliminary tests with montane rodents (genus Nephelomys) in Venezuela. J. Biogeogr. 37, 1378–1393 (2010).Article 

    Google Scholar 
    Pearson, R. G., Raxworthy, C., Nakamura, M. & Peterson, A. T. Predicting species distributions from small numbers of occurrence records: a test case using crypticgeckos in Madagascar. J. Biogeogr. 34, 102–117 (2007).Article 

    Google Scholar 
    Heming, N. M., Dambros, C. & Gutiérrez, E. E. ENMwizard: advanced techniques for Ecological Niche Modeling made easy. https://github.com/HemingNM/ENMwizard (2018).Naimi, B., Hamm, N. A. S., Groen, T. A., Skidmore, A. K. & Toxopeus, A. G. Where is positional uncertainty a problem for species distribution modelling. Ecography 37, 191–203 (2014).Article 

    Google Scholar 
    Muscarella, R. et al. ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for MAXENT ecological niche models. Methods Ecol. Evol. 5, 1198–1205 (2014).Article 

    Google Scholar 
    Owens, H. L. et al. Constraints on interpretation of ecological niche models by limited environmental ranges on calibration areas. Ecol. Model. 263, 10–18 (2013).Article 

    Google Scholar 
    Akaike, H. New look at statistical-model identification. IEEE Trans. Autom. Control AC19, 716–723 (1974).ADS 
    MathSciNet 
    MATH 
    Article 

    Google Scholar 
    Phillips, S. J. et al. Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol. Appl. 19, 181–197 (2009).PubMed 
    Article 

    Google Scholar 
    Barbet-Massin, M., Jiguet, F., Albert, C. H. & Thuiller, W. Selecting pseudo-absences for species distribution models: how, where and how many? Methods Ecol. Evol. 3, 327–338 (2012).Article 

    Google Scholar 
    Bellard, C. et al. Will climate change promote future invasions? Glob. Change Biol. 19, 3740–3748 (2013).ADS 
    Article 

    Google Scholar 
    Elith, J., Kearney, M. & Phillips, S. The art of modelling range-shifting species. Methods Ecol. Evol. 1, 330–342 (2010).Article 

    Google Scholar 
    Anantharaman, R., Hall, K., Shah, V. B. & Edelman, A. Circuitscape in Julia: high performance connectivity modelling to support conservation decisions. Proc. JuliaCon Conf. 1, 58 (2020).Article 

    Google Scholar 
    Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).Article 

    Google Scholar 
    Anderson, D. R. & Burnham, K. P. Avoiding pitfalls when using information-theoretic methods. J. Wildl. Manag. 66, 912–918 (2002).Article 

    Google Scholar 
    Van Strien, M. J., Keller, D. & Holderegger, R. A new analytical approach to landscape genetic modelling: least-cost transect analysis and linear mixed models. Mol. Ecol. 21, 4010–4023 (2012).Article 

    Google Scholar 
    Bartoń, K. MuMIn: multi-model inference, R package version 1.9.13 (2013).Zhang, G. et al. Comparative genomics reveal insights into avian genome evolution and adaptation. Science 346, 1311–1320 (2014).ADS 
    CAS 
    PubMed 
    PubMed Central 
    Article 

    Google Scholar 
    Roesti, M., Kueng, B., Moser, D. & Berner, D. The genomics of ecological vicariance in threespine stickleback fish. Nat. Commun. 6, 8767 (2015).ADS 
    CAS 
    PubMed 
    Article 

    Google Scholar  More