Ecology

Plants and insectsRice (Oryza sativa) cultivar Minghui63 was used in this study. Rice plants were grown in a greenhouse at 27 ± 3 °C with 75 ± 10% RH (relative humidity) and a photoperiod of 16:8 h L:D (light:dark). The cultivation of rice plants followed the same procedure as described previously27. Plants were used for experiments when they were at the tillering stage, which occurred about 44–49 days after sowing.C. suppressalis larvae were reared on an artificial diet as described70. Ten percent honey water solution was provided to supply nutrition for the adults. N. lugens were maintained on a BPH-susceptible rice variety Taichung Native 1 (TN1)38. T. japonicum were obtained from Keyun Industry Co., Ltd (Jiyuan, China). Newly emerged adult wasps were maintained in glass tubes (3.5 cm diameter, 20 cm height) and supplied with 10% honey water solution as a food source and were maintained for at least 6 h to ensure free mating, before females were used for the following experiments. All three species were maintained in climatic chambers at 27 ± 1 °C, 75 ± 5% RH, and a photoperiod of 16:8 h L:D.Performance of caterpillars on insect-infested rice plantsMultiple types of rice plants were prepared: (i) uninfested plants, meaning that potted rice plants remained intact without insect infestation; (ii) SSB-infested plants, each potted rice plant was artificially infested with one 3rd instar SSB larva that had been starved for >3 h for 48 h; (iii) BPH-infested plants, each potted rice plant was artificially infested with a mix of fifteen 3rd and 4th instars BPH nymphs for 48 h; (iv) SSB/BPH-infested plants, each potted rice plant was simultaneously infested with one SSB larva and 15 BPH nymphs for 48 h; (v) SSB → BPH-infested plants, each potted rice plant was artificially infested with one SSB larvae alone for the first 24 h, then 15 BPH nymphs were additionally introduced for another 24 h; (vi) BPH → SSB-infested plants, namely each potted rice plant was artificially infested with 15 BPH nymphs for the first 24 h, then one SSB larvae were additionally introduced for another 24 h. Plant treatments were conducted as described in detail in our previous study27. During herbivory treatment, the uninfested plants were placed in a separate room to avoid possible volatile-mediated interference. During the subsequent bioassays, both SSB caterpillar and BPH nymphs remained in or on the rice plants.Two bioassays were conducted to test the performance of C. suppressalis larvae feeding on differently treated rice plants. The first bioassay included the plant treatments i, ii, iii, and vi, and the second bioassay included the plant treatments i, ii, v, and vi. Three 2-day-old larvae of C. suppressalis were gently introduced onto the middle stem of each rice plant using a soft brush. The infested rice plants were then placed in climatic chambers at 27 ± 1 °C, 75 ± 5% relative humidity, and a photoperiod of 16:8 h L:D. The C. suppressalis larvae were retrieved from the rice plants after 7 days feeding, and they were weighed on a precision balance (CPA2250, Sartorius AG, Germany; readability = 0.01 mg). The mean weight of the three caterpillars on each plant was considered as one biological replicate. The experiment was repeated four times using different batches of plants and herbivores, resulting in a total of 30–46 biological replicates for each treatment.Oviposition-preferences of C. suppressalis females choosing among differently infested rice plantsGreenhouse experimentIn the greenhouse, seven choice tests were conducted with C. suppressalis females including (i) SSB-infested plants versus uninfested plants; (ii) BPH-infested plants versus uninfested plants; (iii) SSB/BPH-infested plants versus uninfested plants; (iv) SSB-infested plants versus BPH-infested plants; (v) SSB-infested plants versus SSB/BPH-infested plants; (vi) BPH-infested plants versus SSB/BPH-infested plants; and (vii) the test in which C. suppressalis females were exposed to all four types of rice plants. The experiments were performed as described in detail by Jiao et al.30. In brief, four potted plants were positioned in the four corners of a cage (80 × 80 × 100 cm) made of 80-mesh nylon nets for each test. For paired comparisons, two potted plants belonging to the same treatment were placed in opposite corners of each age, and in the test with four types of rice plants, each type of plant was positioned in one of the four corners of each cage. Five pairs of freshly emerged moths (less than 1 day) were released in each cage, and a clean Petri dish (9 cm diameter) containing a cotton ball soaked with a 10% honey solution was placed in the center of the cage as food source. After 72 h, the number of individual eggs on each plant were determined. The experiment was conducted in a greenhouse at 27 ± 3 °C, 65 ± 10% RH, and a photoperiod of 16:8 h L:D. Each choice test was repeated with 9–11 times (replicates).Field cage experimentThe oviposition preference of SSB females was further assessed in a field near Langfang City (39.58° N, 116.48° E), China. Four choice tests were conducted: (i) SSB-infested plants versus uninfested plants; (ii) BPH-infested plants versus uninfested plants; (iii) SSB/BPH-infested plants versus uninfested plants; and (iv) SSB/BPH-infested plants versus SSB-infested plants. The treated rice plants were prepared as described above and were transplanted into experimental plots (1.5 × 1.5 m). For each pairwise comparison, six plots of rice plants were covered with a screened cage (8 × 5 × 2.5 m) made of 80-mesh nylon net to prevent moths from entering or escaping. Each of the six plots contained nine rice plants of a particular treatment, with three plots per cage representing the same treatment. Plots were separated by a 1 m buffer and they were alternately distributed in a 3 × 2 grid arrangement in each cage (Supplementary Fig. 4). Approximately 50 mating pairs of newly emerged C. suppressalis adults ( More

Yak and Tibetan sheep thrive under a co-grazing system on the QTP and/or are fed with the same materials; this offers an excellent opportunity to compare the gut microbiota in different host species which share a similar diet. In addition, the grazing systems on the QTP undergo seasonal diets changes in terms of pasture location and forage composition, especially between winter and summer. This presents a good natural “treatment” which helps vary the diets of the yak and Tibetan sheep populations. In the current study, based on a more substantial sample size than the previous study1, we found that diet and environment (represented by seasons winter and summer) superseded host genetics to the family level. That is to say that the gut microbiota of the two animal species showed convergent adaptation to high altitude and harsh environment in QTP, but this convergence had seasonal diets characteristics. These findings may provide a cautionary note for ongoing efforts to link host genetics to gut microbiota composition and function and would provide some food for thought in the breeding of these two livestock groups.The mammalian gut microbiota is acquired from the environment starting at birth, and its assembly and composition is largely shaped by factors such as age, diet, lifestyle, hygiene, and disease state. Researchers subconsciously believe that host species play a greater role than environmental factors when it comes to shaping gut microbiota, especially when there is a large taxonomical difference between the host species. So far, the vast majority of research have focused on the ruminal ecosystem because the rumen is primary site of feed fermentation15,16,17. It is rare to find studies that directly compare the gut microbiota of different species. However, evidence showed that energetically-important microbial products, including VFA (10–13% of total GIT VFA) are produced in the ruminant distal gut3. Hence, it is important to study the composition of distal gut microbiota of ruminants.In this study, at the phylum level, the gut microbiota composition in both groups of livestock was dominated by Bacteroidetes and Firmicutes, which was in agreement with previous reports concerning the yak18. At the same time our result consistently with other study in dairy cows that two dominated phyla Bacteroidetes and Firmicutes found in fecal samples in different seasons were abundant19,20. Firmicutes and Bacteroidetes are responsible for digestion of carbohydrates and proteins, Members of Bacteroidetes having extremely stronger ability to degrade crystalline cellulose. The previous report showed that intestinal microbiome plays an important role in digestion and absorption of the food, and maintaining animals’ health21,22. Intestinal tracts of the ruminants are rich in symbiotic bacteria that helps the body digest plant fibers23,24. Glycans are processed by the distal gut microbiota, generating biologically significant short-chain fatty acids (SCFAs, predominantly acetate, butyrate, and propionate), which serve as the principal energy source for colonocytes25. Fibers may be involved in the regulation of food intake and energy balance via the SCFA-mediated modulation of the secretion of gut hormones26. The higher abundance of Firmicutes and Bacteroidetes in yak may be associated with high-energy consumption at high altitude18.It is worth noting that, at the family level, the dominant genera (Unclassified Ruminococcaceae, Bacteroidaceae, Unclassified BS11, Unclassified Prevotellaceae, Unclassified Christensenellaceae, CF231, Unclassified Mogibacteriaceae and Unclassified Paraprevotellaceae) in the intestines of yak and Tibetan sheep were more greatly influenced by season than genetics (Fig. 5). This has not previously been accurately identified, which may be because there have been few studies into the gut microbial communities of the yak and Tibetan sheep in QTP. So, to improve their husbandry, it is important in the future to study their microbiota profiles using more precise methods such as 16S full-length sequencing or metagenomic sequencing. Ruminococcaceae is a family of autochthonous and benignspecies that primarily inhabit in the caecum and the colon27. It is known that Ruminococcaceae are common in the rumen and hindgut of ruminants, capable of degrading cellulose and starch28. As a member of short chain fatty acid (SCFA) producers, Ruminococcaceae is considered to be the most important fiber and polysaccharides-degrading bacterium in the intestine of herbivores, and produces large amounts of cellulolytic enzymes, including exoglucanases, endoglucanase, glucosidases and hemicellulase29. The microbial community of Yak and Sheep is greatly influenced by alterations in dietary nutrition, Bacteroidaceae have the ability to degrade complex molecules (polysaccharides, proteins) in the intestine18, which can promote the Yak utilizes grasses as its major source of nutrition, due to shortage of grain and other nutrients. Prevotellaceae is responsible for hemicellulose, pectin and high carbohydrate food digestion30. The higher abundance of these microbes may contribute to gaining more energy, and play vital roles in the process of adaption of the hosts to the harsh natural environment15. Bacteroidales BS11 gut group are specialized to active hemicellulose monomeric sugars (e.g., xylose, fucose, mannose and rhamnose) fermentation and short-chain fatty acid (e.g., acetate and butyrate) production that are vital for ruminant energy31. The Bacteroidales BS11 was positively correlated with some metabolites that are involved in amino acid metabolism and biosynthesis, as well as the metabolism of energy sources, such as starch, sucrose, and galactose32.At the genus level, 5-7N15 was most abundant in winter in both animals, on the contrary, the Provotella was predominate. Here, our results indicated that seasonal diets change superseded variations derived from genetic differences between the host species, even though the yak and Tibetan sheep are very different, both taxonomically and in terms of body size. In summer, the forage grass on the Qinghai-Tibet Plateau is dominated by Agropyron cristatum, Elymus nutans, Festuca ovina, Kobresia humilis, Poa pratensis, Stipa aliena, Kobresia pygmaea, Oxytropis biflora, Saussurea hieracioides, Astragalus arnoldii Hemsl. In winter, the main forage was Brachypodium sylvaticum. Carex crebra, Trisetum spicatum and Bupleurum smithii. Stipa has both high palatability and nutritional value, with a high content of crude protein, crude fat, and nitrogen- free extract, and low levels of crude fiber33. The levels of crude protein, crude fat, and nitrogen-free extracts of Brachypodium sylvaticum. Carex crebra, Trisetum spicatum and Bupleurum smithii were lower than that of Stipa, whereas the content of crude fiber was higher than that of Stipa34. Crude protein is the main nutrient of herbage. Crude fat and nitrogen-free extracts provide heat and energy33.Lopes et al. reported that some OTUs known to be functionally relevant for fiber degradation and host development were shared across the entire gastrointestinal tract and present within the feces35. Microbial diversity increases in the distal segments of the gastrointestinal tract. Microbial fermentation appears to be reestablished in the large intestine, with the proportion of acetate, propionate and butyrate being similar to the rumen.Several explanations for this phenomenon are possible. Firstly, both the yak and Tibetan sheep are ruminants. In herbivores, the gut microbiota is dominated by Firmicutes and Bacteroides, the functions of which are related to cellulose digestion36. Therefore, ruminant microbes could possibly be more similar across species than gut microbes from elsewhere.Secondly, the yaks and Tibetan sheep in our study co-grazed from birth to death. As such, the initial gut microbiota source, responsible for populating the remainder of the gut in the months and years after the initial seeding at birth, would necessarily come from the same environment. It has been established that early life events are critical for gut microbiota development and for shaping the adult microbiota. Lifestyle and diet will further influence the composition and function of the gut microbiota. In our study, the investigated animals shared a very similar lifestyle and obtained their diets from the same source. The results revealed that sheep and yaks presented almost identical gut microbiota compositions in the winter, but by the date of collection of the summer samples they were quite different. The reason for this could be that during summer and summer there is pronounced pastoral grass growth, giving the animals more variety and choice in their diets; it is known, after all, that sheep have different diet preferences to yaks37. However, during the winter, the animals have no option but to eat the same food in order to survive until winter.Thirdly, there could be a convergent evolution of gut microbiomes in yaks and Tibetan sheep due to the extremely harsh environment in high-altitude regions1,38. When compared with their low-altitude relatives, cattle (Bos taurus) and ordinary sheep (Ovis aries), metagenomic analyses revealed significant enrichment in rumen microbial genes involving volatile fatty acid-yielding pathways in yaks and Tibetan sheep, whereas methanogenesis pathways were enriched in the cattle metagenome. Analyses of RNA transcriptomes revealed significant upregulation in 36 genes associated with volatile fatty acid transport and absorption in the ruminal epithelium of yaks and Tibetan sheep. This suggests that, aside from host genetics, long-term exposure to harsh environments has allowed the gut microbiome to adapt in order to boost health and survival. In other words, although yaks and Tibetan sheep are very different genetically, their gut microbiota could be similar due to the selection pressures of the high altitude at which they live. Meanwhile, from our data based on functional gene composition (Fig. S3), it is also worth noting that there were no groups clearly distinguished from one another, although the PERMANOVA results indicated both a host and season effect, with the interaction between them being statistically significant. Though factors such as environment and diet (represented by seasons) can trump host genetics, we could not ignore the interplay of these factors as gut microbes are a very complex community.Winter is the harshest period for the survival of yak and Tibetan sheep. To maintain the survival, it’s best to feed the animals with a high protein content. Furthermore, to get more detailed data in different seasons and various dietary habits of yak and sheep, more study should be assessed about intestinal microbiota by collecting feces. More

1.Letko, M., Seifert, S. N., Olival, K. J., Plowright, R. K. & Munster, V. J. Bat-borne virus diversity, spillover and emergence. Nat. Rev. Microbiol. 18, 461–471 (2020). This is a review of the overall diversity of bat-borne coronaviruses and research agenda for enhanced characterization of their zoonotic and pandemic potential.CAS 
PubMed 

Google Scholar 
2.Dobson, B. A. P. et al. Ecology and economics for pandemic prevention. Science 369, 379–381 (2020).CAS 
PubMed 

Google Scholar 
3.Plowright, R. K. et al. Land use-induced spillover: a call to action to safeguard environmental, animal, and human health. Lancet Planet. Heal. 5, e237–e245 (2021).
Google Scholar 
4.Shivaprakash, K. N., Sen, S., Paul, S., Kiesecker, J. M. & Bawa, K. S. Mammals, wildlife trade, and the next global pandemic. Curr. Biol. 31, 3671–3677.e3 (2021).CAS 
PubMed 

Google Scholar 
5.Huong, N. Q. et al. Coronavirus testing indicates transmission risk increases along wildlife supply chains for human consumption in Viet Nam, 2013–2014. PLoS ONE 15, e0237129 (2020). This study provides evidence that coronavirus detection increases along the supply chain of rodents destined for human consumption.CAS 
PubMed 
PubMed Central 

Google Scholar 
6.Xiao, X., Newman, C., Buesching, C. D., Macdonald, D. W. & Zhou, Z.-M. Animal sales from Wuhan wet markets immediately prior to the COVID-19 pandemic. Sci. Rep. 11, 11898 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
7.Chen, L., Liu, B., Yang, J. & Jin, Q. DBatVir: the database of bat-associated viruses. Database 2014, 1–7 (2014).
Google Scholar 
8.Tao, Y. et al. Surveillance of bat coronaviruses in Kenya identifies relatives of human coronaviruses NL63 and 229E and their recombination history. J. Virol. 91, 1–16 (2017).CAS 

Google Scholar 
9.Cui, J., Li, F. & Shi, Z. L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 17, 181–192 (2019).CAS 
PubMed 

Google Scholar 
10.Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273 (2020). This is one of the first studies to discover viruses related to SARS-CoV-2 in wild Rhinolophus spp. bats in China.CAS 
PubMed 
PubMed Central 

Google Scholar 
11.Forni, D., Cagliani, R., Clerici, M. & Sironi, M. Molecular evolution of human coronavirus genomes. Trends Microbiol. 25, 35–48 (2017).CAS 
PubMed 

Google Scholar 
12.Zhou, P. et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature 556, 255–258 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
13.Wang, N. et al. Serological evidence of bat SARS-related coronavirus infection in humans, China. Virol. Sin. 33, 104–107 (2018). This study provides evidence of potentially undetected spillovers of bat-associated coronaviruses in rural human populations in China.PubMed 
PubMed Central 

Google Scholar 
14.Li, H. et al. Human-animal interactions and bat coronavirus spillover potential among rural residents in southern China. Biosaf. Heal. 1, 84–90 (2019).
Google Scholar 
15.Woo, P. C. Y. et al. Discovery of seven novel mammalian and avian coronaviruses in the genus Deltacoronavirus supports bat coronaviruses as the gene source of Alphacoronavirus and Betacoronavirus and avian coronaviruses as the gene source of Gammacoronavirus and Deltacoronavirus. J. Virol. 86, 3995–4008 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
16.Anthony, S. J. et al. Global patterns in coronavirus diversity. Virus Evol. 3, 1–15 (2017). This is a review of ecological patterns of associations between bats and coronaviruses with information up to 2014.
Google Scholar 
17.Li, W. et al. Bats are natural reservoirs of SARS-like coronaviruses. Science 310, 676–679 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
18.Lau, S. K. P. et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc. Natl Acad. Sci. USA 102, 14040–14045 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
19.Anthony, S. et al. Coronaviruses in bats from Mexico. J. Gen. Virol. 94, 1028–1038 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
20.Latinne, A. et al. Origin and cross-species transmission of bat coronaviruses in China. Nat. Commun. 11, 4235 (2020). Using 5 years of surveillance data on coronaviruses in bats in China, the authors show that host switching is common in bat coronaviruses, particularly in Rhinolophus spp.CAS 
PubMed 
PubMed Central 

Google Scholar 
21.Ithete, N. L. et al. Close relative of human Middle East respiratory syndrome coronavirus in bat, South Africa. Emerg. Infect. Dis. 19, 1697–1699 (2013).PubMed 
PubMed Central 

Google Scholar 
22.Råberg, L., Graham, A. L. & Read, A. F. Decomposing health: tolerance and resistance to parasites in animals. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 364, 37–49 (2009).PubMed 

Google Scholar 
23.Schlottau, K. et al. SARS-CoV-2 in fruit bats, ferrets, pigs, and chickens: an experimental transmission study. Lancet Microbe 1, e218–e225 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
24.Munster, V. J. et al. Replication and shedding of MERS-CoV in Jamaican fruit bats (Artibeus jamaicensis). Sci. Rep. 6, 1–10 (2016).
Google Scholar 
25.van Doremalen, N. et al. SARS-like coronavirus WIV1-CoV does not replicate in Egyptian fruit bats (Rousettus aegyptiacus). Viruses 10, 727 (2018).PubMed Central 

Google Scholar 
26.Plowright, R. K. et al. Transmission or within-host dynamics driving pulses of zoonotic viruses in reservoir–host populations. PLoS Negl. Trop. Dis. 10, 1–21 (2016).
Google Scholar 
27.Jeong, J. et al. Persistent infections support maintenance of a coronavirus in a population of Australian bats (Myotis macropus). Epidemiol. Infect. 145, 2053–2061 (2017).CAS 
PubMed 

Google Scholar 
28.Watanabe, S. et al. Bat coronaviruses and experimental infection of bats, the Philippines. Emerg. Infect. Dis. 16, 1217–1223 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
29.Subudhi, S. et al. A persistently infecting coronavirus in hibernating Myotis lucifugus, the North American little brown bat. J. Gen. Virol. 98, 2297–2309 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
30.Widagdo, W. et al. Tissue distribution of the MERS-coronavirus receptor in bats. Sci. Rep. 7, 1–8 (2017).CAS 

Google Scholar 
31.Banerjee, A. et al. Selection of viral variants during persistent infection of insectivorous bat cells with Middle East respiratory syndrome coronavirus. Sci. Rep. 10, 7257 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
32.Wang, M.-N. N. et al. Longitudinal surveillance of SARS-like coronaviruses in bats by quantitative real-time PCR. Virol. Sin. 31, 78–80 (2016).PubMed 
PubMed Central 

Google Scholar 
33.Ge, X. Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535–538 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
34.Hu, B. et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. PLoS Pathog. 13, 1–27 (2017).CAS 

Google Scholar 
35.Smith, C. Australian bat coronaviruses (The University of Queensland, 2015).36.Baldwin, H. J. Epidemiology and ecology of virus and host: bats and coronaviruses in Ghana, West Africa (Macquarie University & Ulm University, 2015).37.Joffrin, L. et al. Bat coronavirus phylogeography in the Western Indian Ocean. Sci. Rep. 10, 1–11 (2020).
Google Scholar 
38.Plowright, R. K. et al. Reproduction and nutritional stress are risk factors for Hendra virus infection in little red flying foxes (Pteropus scapulatus). Proc. R. Soc. B Biol. Sci. 275, 861–869 (2008).
Google Scholar 
39.Peel, A. J. et al. Synchronous shedding of multiple bat paramyxoviruses coincides with peak periods of Hendra virus spillover. Emerg. Microbes Infect. 8, 1314–1323 (2019). This study provides evidence of co-circulation of multiple viruses in single and multispecies roosts of flying foxes, with higher diversity of viruses in mixed-species roosts.PubMed 
PubMed Central 

Google Scholar 
40.Wacharapluesadee, S. et al. Longitudinal study of age-specific pattern of coronavirus infection in Lyle’s flying fox (Pteropus lylei) in Thailand. Virol. J. 15, 1–10 (2018).
Google Scholar 
41.Lau, S. K. P. et al. Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. J. Virol. 84, 2808–2819 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
42.Willoughby, A., Phelps, K. & Olival, K. A comparative analysis of viral richness and viral sharing in cave-roosting bats. Diversity 9, 35 (2017).
Google Scholar 
43.Lloyd-Smith, J. O., Schreiber, S. J., Kopp, P. E. & Getz, W. M. Superspreading and the effect of individual variation on disease emergence. Nature 438, 355–359 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
44.Lee, K. A. Linking immune defenses and life history at the levels of the individual and the species. Integr. Comp. Biol. 46, 1000–1015 (2006).CAS 
PubMed 

Google Scholar 
45.Robinson, D. P. & Klein, S. L. Pregnancy and pregnancy-associated hormones alter immune responses and disease pathogenesis. Horm. Behav. 62, 263–271 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
46.Pauly, M. et al. Novel alphacoronaviruses and paramyxoviruses cocirculate with type 1 and severe acute respiratory system (SARS)-related betacoronaviruses in synanthropic bats of Luxembourg. Appl. Environ. Microbiol. https://doi.org/10.1128/AEM.01326-17 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
47.Drexler, J. F. et al. Amplification of emerging viruses in a bat colony. Emerg. Infect. Dis. 17, 449–456 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
48.Annan, A. et al. Human betacoronavirus 2c EMC/2012-related viruses in bats, Ghana and Europe. Emerg. Infect. Dis. 19, 456–459 (2013).PubMed 
PubMed Central 

Google Scholar 
49.Montecino-Latorre, D. et al. Reproduction of East-African bats may guide risk mitigation for coronavirus spillover. One Heal. Outlook 2, 2 (2020).
Google Scholar 
50.Plowright, R. K., Becker, D. J., McCallum, H. & Manlove, K. R. Sampling to elucidate the dynamics of infections in reservoir hosts. Philos. Trans. R. Soc. B Biol. Sci. https://doi.org/10.1098/rstb.2018.0336 (2019).Article 

Google Scholar 
51.Vanalli, C. et al. Within-host mechanisms of immune regulation explain the contrasting dynamics of two helminth species in both single and dual infections. PLoS Comput. Biol. 16, 1–19 (2020).
Google Scholar 
52.Ge, X. Y. et al. Coexistence of multiple coronaviruses in several bat colonies in an abandoned mineshaft. Virol. Sin. 31, 31–40 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
53.Chu, D. K. W., Peiris, J. S. M., Chen, H., Guan, Y. & Poon, L. L. M. Genomic characterizations of bat coronaviruses (1A, 1B and HKU8) and evidence for co-infections in Miniopterus bats. J. Gen. Virol. 89, 1282–1287 (2008).CAS 
PubMed 

Google Scholar 
54.Drexler, J. F. et al. Genomic characterization of severe acute respiratory syndrome-related coronavirus in European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences. J. Virol. 84, 11336–11349 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
55.Tong, S. et al. Detection of novel SARS-like and other coronaviruses in bats from Kenya. Emerg. Infect. Dis. 15, 482–485 (2009).PubMed 
PubMed Central 

Google Scholar 
56.Wacharapluesadee, S. et al. Diversity of coronavirus in bats from eastern Thailand emerging viruses. Virol. J. 12, 1–7 (2015).
Google Scholar 
57.Valitutto, M. T. et al. Detection of novel coronaviruses in bats in Myanmar. PLoS ONE 15, e0230802 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
58.Anthony, S. J. et al. A strategy to estimate unknown viral diversity in mammals. mBio 4, 289 (2013).
Google Scholar 
59.Prada, D., Boyd, V., Baker, M. L., O’Dea, M. & Jackson, B. Viral diversity of microbats within the south west botanical province of Western Australia. Viruses 11, 1–21 (2019).
Google Scholar 
60.Seltmann, A. et al. Seasonal fluctuations of astrovirus, but not coronavirus shedding in bats inhabiting human-modified tropical forests. Ecohealth 14, 272–284 (2017).PubMed 
PubMed Central 

Google Scholar 
61.Chu, D. K. W., Poon, L. L. M., Guan, Y. & Peiris, J. S. M. Novel astroviruses in insectivorous Bats. J. Virol. 82, 9107–9114 (2008).CAS 
PubMed 
PubMed Central 

Google Scholar 
62.Kemenesi, G. et al. Molecular survey of RNA viruses in Hungarian bats: discovering novel astroviruses, coronaviruses, and caliciviruses. Vector Borne Zoonotic Dis. 14, 846–855 (2014).PubMed 

Google Scholar 
63.Rizzo, F. et al. Coronavirus and paramyxovirus in bats from northwest Italy. BMC Vet. Res. 13, 1–11 (2017).
Google Scholar 
64.Paskey, A. C. et al. The temporal RNA virome patterns of a lesser dawn bat (Eonycteris spelaea) colony revealed by deep sequencing. Virus Evol. 6, 1–14 (2020).
Google Scholar 
65.Davy, C. M. et al. White-nose syndrome is associated with increased replication of a naturally persisting coronaviruses in bats. Sci. Rep. 8, 15508 (2018).PubMed 
PubMed Central 

Google Scholar 
66.Woo, P. C. Y., Lau, S. K. P., Huang, Y. & Yuen, K.-Y. Y. Coronavirus diversity, phylogeny and interspecies jumping. Exp. Biol. Med. 234, 1117–1127 (2009).CAS 

Google Scholar 
67.Fehr, A. R. & Perlman, S. in Coronaviruses. Methods in Molecular Biology Vol 1282 (eds Maier, H., Bickerton, E. & Britton, P.) 1–23 (Humana Press, 2015).68.Duffy, S., Shackelton, L. A. & Holmes, E. C. Rates of evolutionary change in viruses: patterns and determinants. Nat. Rev. Genet. 9, 267–276 (2008).CAS 
PubMed 

Google Scholar 
69.Jenkins, G. M., Rambaut, A., Pybus, O. G. & Holmes, E. C. Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J. Mol. Evol. 54, 156–165 (2002).CAS 
PubMed 

Google Scholar 
70.Eckerle, L. D. et al. Infidelity of SARS-CoV Nsp14-exonuclease mutant virus replication is revealed by complete genome sequencing. PLoS Pathog. 6, 1–15 (2010).
Google Scholar 
71.Ogando, N. S. et al. The curious case of the nidovirus exoribonuclease: its role in RNA synthesis and replication fidelity. Front. Microbiol. 10, 1–17 (2019).
Google Scholar 
72.Nga, P. T. et al. Discovery of the first insect nidovirus, a missing evolutionary link in the emergence of the largest RNA virus genomes. PLoS Pathog. 7, e1002215 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
73.Smith, E., Blanc, H., Vignuzzi, M. & Denison, M. R. Coronaviruses lacking exoribonuclease activity are susceptible to lethal mutagenesis: evidence for proofreading and potential therapeutics. PLoS Pathog. 9, e1003565 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
74.Martin, L. B. et al. Extreme competence: keystone hosts of infections. Trends Ecol. Evol. 34, 303–314 (2019).PubMed 
PubMed Central 

Google Scholar 
75.Graham, R. L. & Baric, R. S. Recombination, reservoirs, and the modular spike: mechanisms of coronavirus cross-species transmission. J. Virol. 84, 3134–3146 (2010).CAS 
PubMed 

Google Scholar 
76.Letko, M. et al. Adaptive evolution of MERS-CoV to species variation in DPP4. Cell Rep. 24, 1730–1737 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
77.Ermonval, M., Baychelier, F. & Tordo, N. What do we know about how hantaviruses interact with their different hosts? Viruses 8, 223 (2016).PubMed Central 

Google Scholar 
78.Su, S. et al. Epidemiology, genetic recombination, and pathogenesis of coronaviruses. Trends Microbiol. 24, 490–502 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
79.Lai, M. M. et al. Recombination between nonsegmented RNA genomes of murine coronaviruses. J. Virol. 56, 449–456 (1985).CAS 
PubMed 
PubMed Central 

Google Scholar 
80.Tian, P.-F. et al. Evidence of recombinant strains of porcine epidemic diarrhea virus, United States, 2013. Emerg. Infect. Dis. 20, 1731–1734 (2014).
Google Scholar 
81.Terada, Y. et al. Emergence of pathogenic coronaviruses in cats by homologous recombination between feline and canine coronaviruses. PLoS ONE 9, e106534 (2014).PubMed 
PubMed Central 

Google Scholar 
82.Decaro, N. et al. Recombinant canine coronaviruses related to transmissible gastroenteritis virus of swine are circulating in dogs. J. Virol. 83, 1532–1537 (2009).CAS 
PubMed 

Google Scholar 
83.Zhang, Y. et al. Genotype shift in human coronavirus OC43 and emergence of a novel genotype by natural recombination. J. Infect. 70, 641–650 (2015).PubMed 

Google Scholar 
84.Pyrc, K., Berkhout, B. & van der hoek, L. in Recent Research in Development of Infection & Immunity 3rd edn 25–48 (Transworld Research Network, 2005).85.Woo, P. C. Y. et al. Phylogenetic and recombination analysis of coronavirus HKU1, a novel coronavirus from patients with pneumonia. Arch. Virol. 150, 2299–2311 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
86.Zhang, X. W., Yap, Y. L. & Danchin, A. Testing the hypothesis of a recombinant origin of the SARS-associated coronavirus. Arch. Virol. 150, 1–20 (2005).CAS 
PubMed 

Google Scholar 
87.Stanhope, M. J., Brown, J. R. & Amrine-Madsen, H. Evidence from the evolutionary analysis of nucleotide sequences for a recombinant history of SARS-CoV. Infect. Genet. Evol. 4, 15–19 (2004).CAS 
PubMed 

Google Scholar 
88.Wang, Y. et al. Origin and possible genetic recombination of the middle east respiratory syndrome coronavirus from the first imported case in China: phylogenetics and coalescence analysis. MBio 6, 1–6 (2015).
Google Scholar 
89.Huang, C. et al. A bat-derived putative cross-family recombinant coronavirus with a reovirus gene. PLoS Pathog. 12, 1–25 (2016). This study provides evidence of cross-family recombination between coronaviruses and reoviruses.
Google Scholar 
90.Drexler, J. F., Corman, V. M. & Drosten, C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. Antivir. Res. 101, 45–56 (2014).CAS 
PubMed 

Google Scholar 
91.Cui, J. et al. Evolutionary relationships between bat coronaviruses and their hosts. Emerg. Infect. Dis. 13, 1526–1532 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
92.Davies, N. et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Sci 372, eabg3055 (2021).CAS 

Google Scholar 
93.Reguera, J., Mudgal, G., Santiago, C. & Casasnovas, J. M. A structural view of coronavirus-receptor interactions. Virus Res. 194, 3–15 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
94.Lim, Y., Ng, Y., Tam, J. & Liu, D. Human coronaviruses: a review of virus–host interactions. Diseases 4, 26 (2016).
Google Scholar 
95.Masters, P. S. & Perlman, S. Coronaviridae. Fields Virol. 1, 825–858 (2013).
Google Scholar 
96.Letko, M., Marzi, A. & Munster, V. Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses. Nat. Microbiol. 5, 562–569 (2020). This study uses a functional viromics platform to rapidly characterize the zoonotic potential of new coronaviruses on the basis of genome sequences.CAS 
PubMed 
PubMed Central 

Google Scholar 
97.van Doremalen, N. et al. Host species restriction of middle east respiratory syndrome coronavirus through its receptor, dipeptidyl peptidase 4. J. Virol. 88, 9220–9232 (2014).PubMed 
PubMed Central 

Google Scholar 
98.Conceicao, C. et al. The SARS-CoV-2 spike protein has a broad tropism for mammalian ACE2 proteins. PLoS Biol. 18, e3001016 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
99.Li, W. et al. Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2. EMBO J. 24, 1634–1643 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
100.Damas, J. et al. Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates. Proc. Natl Acad. Sci. USA 117, 22311–22322 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
101.Shi, J. et al. Susceptibility of ferrets, cats, dogs, and other domesticated animals to SARS–coronavirus 2. Science 1020, 1016–1020 (2020).
Google Scholar 
102.Oude Munnink, B. B. et al. Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans. Science 371, 172–177 (2021). The study provides evidence of spillback of SARS-CoV-2 into mink populations, rapid and widespread transmission, and seeding of mink-associated genetic variants back into humans.CAS 
PubMed 

Google Scholar 
103.Hoffmann, M. et al. Differential sensitivity of bat cells to infection by enveloped RNA viruses: coronaviruses, paramyxoviruses, filoviruses, and influenza viruses. PLoS ONE 8, e72942 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
104.Li, W. et al. Animal origins of the severe acute respiratory syndrome coronavirus: insight from ACE2-S-protein interactions. J. Virol. 80, 4211–4219 (2006).CAS 
PubMed 
PubMed Central 

Google Scholar 
105.Hou, Y. et al. Angiotensin-converting enzyme 2 (ACE2) proteins of different bat species confer variable susceptibility to SARS-CoV entry. Arch. Virol. 155, 1563–1569 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
106.Boni, M. F. et al. Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic. Nat. Microbiol. 5, 1408–1417 (2020). This study provids phylogenetic evidence that suggests that the ancestral lineages from which SARS-CoV-2 may have originated have circulated undetected in bats for decades.CAS 
PubMed 

Google Scholar 
107.Wells, H. L. et al. The evolutionary history of ACE2 usage within the coronavirus subgenus Sarbecovirus. Virus Evol. 7, 1–22 (2021).
Google Scholar 
108.Urbanowicz, R. A. et al. Human adaptation of Ebola virus during the West African outbreak. Cell 167, 1079–1087.e5 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
109.Gupta, A. et al. Extrapulmonary manifestations of COVID-19. Nat. Med. 26, 1017–1032 (2020).CAS 

Google Scholar 
110.Hamming, I. et al. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J. Pathol. 203, 631–637 (2004).CAS 
PubMed 
PubMed Central 

Google Scholar 
111.Lam, T. T. Y. et al. Identifying SARS-CoV-2-related coronaviruses in Malayan pangolins. Nature 583, 282–285 (2020).CAS 
PubMed 

Google Scholar 
112.Martina, B. E. E. et al. SARS virus infection of cats and ferrets. Nature 425, 915 (2003).CAS 
PubMed 
PubMed Central 

Google Scholar 
113.Munster, V. J. et al. Respiratory disease in rhesus macaques inoculated with SARS-CoV-2. Nature 585, 268–272 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
114.Hou, Y. J. et al. SARS-CoV-2 reverse genetics reveals a variable infection gradient in the respiratory tract. Cell 182, 429–446.e14 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
115.Menachery, V. D. et al. Trypsin treatment unlocks barrier for zoonotic bat coronavirus infection. J. Virol. 94, 1–15 (2019).
Google Scholar 
116.Qian, Z., Dominguez, S. R. & Holmes, K. V. Role of the spike glycoprotein of human Middle East respiratory syndrome coronavirus (MERS-CoV) in virus entry and syncytia formation. PLoS ONE 8, 1–12 (2013).
Google Scholar 
117.Belouzard, S., Chu, V. C. & Whittaker, G. R. Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites. Proc. Natl Acad. Sci. USA 106, 5871–5876 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
118.Barlan, A. et al. Receptor variation and susceptibility to middle east respiratory syndrome coronavirus infection. J. Virol. 88, 4953–4961 (2014).PubMed 
PubMed Central 

Google Scholar 
119.Zheng, Y. et al. Lysosomal proteases are a determinant of coronavirus tropism. J. Virol. 92, 1–14 (2018).
Google Scholar 
120.Bertram, S. et al. Cleavage and activation of the severe acute respiratory syndrome coronavirus spike protein by human airway trypsin-like protease. J. Virol. 85, 13363–13372 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
121.Matsuyama, S., Ujike, M., Morikawa, S., Tashiro, M. & Taguchi, F. Protease-mediated enhancement of severe acute respiratory syndrome coronavirus infection. Proc. Natl Acad. Sci. USA 102, 12543–12547 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
122.Yang, Y. et al. Receptor usage and cell entry of bat coronavirus HKU4 provide insight into bat-to-human transmission of MERS coronavirus. Proc. Natl Acad. Sci. USA 111, 12516–12521 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
123.Wacharapluesadee, S. et al. Group C betacoronavirus in bat guano fertilizer, Thailand. Emerg. Infect. Dis. 19, 1349–1351 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
124.Luk, H. K. H., Li, X., Fung, J., Lau, S. K. P. & Woo, P. C. Y. Molecular epidemiology, evolution and phylogeny of SARS coronavirus. Infect. Genet. Evol. 71, 21–30 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
125.Wu, K., Peng, G., Wilken, M., Geraghty, R. J. & Li, F. Mechanisms of host receptor adaptation by severe acute respiratory syndrome coronavirus. J. Biol. Chem. 287, 8904–8911 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
126.Daszak, P., Olival, K. J. & Li, H. A strategy to prevent future epidemics similar to the 2019-nCoV outbreak. Biosaf. Heal. 2, 6–8 (2020).
Google Scholar 
127.Tidemann, C. & Vardon, M. Pests, pestilence, pollen and pot roasts: the need for community based management of flying foxes in Australia. Aust. Biol. 10, 77–83 (1997).
Google Scholar 
128.Mickleburgh, S., Waylen, K. & Racey, P. Bats as bushmeat: a global review. Oryx 43, 217–234 (2009).
Google Scholar 
129.Tuttle, M. D. & Moreno, A. Cave-Dwelling Bats of Northern Mexico: Their Value and Conservation Needs (Bat Conservation International, 2005).130.Rulli, M. C., D’Odorico, P., Galli, N. & Hayman, D. T. S. Land-use change and the livestock revolution increase the risk of zoonotic coronavirus transmission from rhinolophid bats. Nat. Food https://doi.org/10.1038/s43016-021-00285-x (2021).Article 

Google Scholar 
131.Mckee, C. D., Islam, A., Luby, S. P., Salje, H. & Hudson, P. J. The ecology of Nipah virus in Bangladesh: a nexus of land use change and opportunistic feeding behavior in bats. Viruses 13, 169 (2020).
Google Scholar 
132.Kessler, M. K. et al. Changing resource landscapes and spillover of henipaviruses. Ann. N. Y. Acad. Sci. https://doi.org/10.1111/nyas.13910 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
133.Dighe, A., Jombart, T., Van Kerkhove, M. D. & Ferguson, N. A systematic review of MERS-CoV seroprevalence and RNA prevalence in dromedary camels: implications for animal vaccination. Epidemics 29, 100350 (2019).PubMed 
PubMed Central 

Google Scholar 
134.Hegde, S. T. et al. Using healthcare-seeking behaviour to estimate the number of Nipah outbreaks missed by hospital-based surveillance in Bangladesh. Int. J. Epidemiol. 48, 1219–1227 (2019).PubMed 
PubMed Central 

Google Scholar 
135.Glennon, E. E., Jephcott, F. L., Restif, O. & Wood, J. L. N. Estimating undetected Ebola spillovers. PLoS Negl. Trop. Dis. 13, 1–10 (2019).
Google Scholar 
136.Matson, M. J., Chertow, D. S. & Munster, V. J. Delayed recognition of Ebola virus disease is associated with longer and larger outbreaks. Emerg. Microbes Infect. 9, 291–301 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
137.Zheng, B. J. et al. SARS-related virus predating SARS outbreak, Hong Kong. Emerg. Infect. Dis. 10, 176–178 (2004). This study provides serological evidence that populations in Hong Kong sampled in 2001 may have been exposed to SARS-CoV or related viruses in bats or other animals before the first SARS outbreaks.PubMed 
PubMed Central 

Google Scholar 
138.Yu, S. et al. Retrospective serological investigation of severe acute respiratory syndrome coronavirus antibodies in recruits from mainland China. Clin. Diagn. Lab. Immunol. 12, 552–554 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
139.Lloyd-Smith, J. O. et al. Epidemic dynamics at the human-animal interface. Science 326, 1362–1367 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
140.Plowright, R. K. et al. Pathways to zoonotic spillover. Nat. Rev. Microbiol. 15, 502–510 (2017). This study formulates a conceptual model for the multiple layers of ecological and cellular barriers that affect the likelihood of pathogen spillover from animals.CAS 
PubMed 
PubMed Central 

Google Scholar 
141.Klompus, S. et al. Cross-reactive antibodies against human coronaviruses and the animal coronavirome suggest diagnostics for future zoonotic spillovers. Sci. Immunol. 6, eabe9950 (2021).PubMed 

Google Scholar 
142.Field, H. E. Hendra virus ecology and transmission. Curr. Opin. Virol. 16, 120–125 (2016).PubMed 

Google Scholar 
143.Chua, K. B. Nipah virus outbreak in Malaysia. J. Clin. Virol. 26, 265–275 (2003).PubMed 

Google Scholar 
144.Azhar, E. I. et al. Evidence for camel-to-human transmission of MERS coronavirus. N. Engl. J. Med. 370, 2499–2505 (2014).CAS 
PubMed 

Google Scholar 
145.Memish, Z. A. et al. Respiratory tract samples, viral load, and genome fraction yield in patients with middle east respiratory syndrome. J. Infect. Dis. 210, 1590–1594 (2014).CAS 
PubMed 

Google Scholar 
146.Buchholz, U. et al. Contact investigation of a case of human novel coronavirus infection treated in a German hospital, October-November 2012. Eurosurveillance 18, 1–7 (2013).
Google Scholar 
147.Chu, D. K. W. et al. MERS coronaviruses in dromedary camels, Egypt. Emerg. Infect. Dis. 20, 1049–1053 (2014).PubMed 
PubMed Central 

Google Scholar 
148.Zhou, H. et al. Identification of novel bat coronaviruses sheds light on the evolutionary origins of SARS-CoV-2 and related viruses. Cell 184, 4380–4391.e14 (2021). This study reports the discovery of additional coronaviruses related to SARS-CoV-2 in Rhinolophus spp. and use of ecological modelling to highlight areas of southern China and South-East Asia as hotspots of Rhinolophus species diversity.CAS 
PubMed 
PubMed Central 

Google Scholar 
149.Larsen, H. D. et al. Preliminary report of an outbreak of SARS-CoV-2 in mink and mink farmers associated with community spread, Denmark, June to November 2020. Eur. Surveill. 26, 2100009 (2021).CAS 

Google Scholar 
150.Bertzbach, L. D. et al. SARS-CoV-2 infection of Chinese hamsters (Cricetulus griseus) reproduces COVID-19 pneumonia in a well-established small animal model. Transbound. Emerg. Dis. 68, 1075–1079 (2021).CAS 
PubMed 

Google Scholar 
151.Fagre, A. et al. SARS-CoV-2 infection, neuropathogenesis and transmission among deer mice: Implications for spillback to New World rodents. PLoS Pathog. 17, e1009585 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
152.Imai, M. et al. Syrian hamsters as a small animal model for SARS-CoV-2 infection and countermeasure development. Proc. Natl Acad. Sci. USA 117, 16587–16595 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
153.Sia, S. F. et al. Pathogenesis and transmission of SARS-CoV-2 in golden hamsters. Nature 583, 834–838 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
154.Halfmann, P. J. et al. Transmission of SARS-CoV-2 in domestic cats. N. Engl. J. Med. 383, 592–594 (2020).PubMed 

Google Scholar 
155.Griffin, B. D. et al. SARS-CoV-2 infection and transmission in the North American deer mouse. Nat. Commun. 12, 1–10 (2021).
Google Scholar 
156.Palmer, M. V. et al. Susceptibility of white-tailed deer (Odocoileus virginianus) to SARS-CoV-2. J. Virol. 95, e00083–21 (2021).CAS 
PubMed Central 

Google Scholar 
157.Plowright, R. K. & Hudson, P. J. From protein to pandemic: the transdisciplinary approach needed to prevent spillover and the next pandemic. Viruses 13, 1298 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
158.Obameso, J. O. et al. The persistent prevalence and evolution of cross-family recombinant coronavirus GCCDC1 among a bat population: a two-year follow-up. Sci. China Life Sci. 60, 1357–1363 (2017). This study provides evidence of coronavirus evolution in a longitudinally sampled population of bats.PubMed 
PubMed Central 

Google Scholar 
159.Lazov, C. et al. Detection and characterization of distinct alphacoronaviruses in five different bat species in Denmark. Viruses 10, 486 (2018).PubMed Central 

Google Scholar 
160.Hu, D. et al. Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats. Emerg. Microbes Infect. 7, 1–10 (2018).
Google Scholar 
161.Pepin, K. M., Lass, S., Pulliam, J. R. C., Read, A. F. & Lloyd-Smith, J. O. Identifying genetic markers of adaptation for surveillance of viral host jumps. Nat. Rev. Microbiol. 8, 802–813 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
162.Hemida, M. G. et al. Coronavirus infections in horses in Saudi Arabia and Oman. Transbound. Emerg. Dis. 64, 2093–2103 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
163.Zhuang, Q. et al. Surveillance and taxonomic analysis of the coronavirus dominant in pigeons in China. Transbound. Emerg. Dis. 67, 1981–1990 (2020).CAS 

Google Scholar 
164.O’Brien, S. J. et al. Genetic basis for species vulnerability in the cheetah. Science 227, 1428–1434 (1985).PubMed 

Google Scholar 
165.Herrewegh, A. A. P. M., Smeenk, I., Horzinek, M. C., Rottier, P. J. M. & de Groot, R. J. Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double recombination between feline coronavirus type I and canine coronavirus. J. Virol. 72, 4508–4514 (1998).CAS 
PubMed 
PubMed Central 

Google Scholar  More

Effects of temperature on C storage in soils with contrasting stabilisation capacitiesUsing a space-for-time approach, we defined the effect of temperature on C storage as the proportional reduction in C storage for each 10 oC increase in mean annual temperature. In this context, a value of 1 indicates no change in C storage with temperature, values less than 1 indicate C stocks increase with temperature and values greater than 1 indicate C stocks decline with temperature, with, for example, a value of 2 indicating that C stocks halve for every 10 oC increase in temperature. C storage in the top 50 cm of mineral soil declined strongly with increasing temperature, declining by a factor of ~1.4 per 10 oC (Fig. 1b). Critically, the nature of this relationship was modified by soil clay content; C storage in fine-textured soils with greater stabilisation capacities was affected much less by temperature than C storage in coarse-textured soils (Fig. 2a and Supplementary Fig. 1; factors of up to 1.9 per 10 oC for coarse-textured soils versus factors below 1.2 per 10 oC for finer-textured soils). We also demonstrate that the lower effect of temperature on C storage in fine-textured soils was retained after accounting for potentially confounding variation in precipitation, aridity (actual minus potential evapotranspiration), plant productivity, soil pH and cation exchange capacity (CEXC) (Fig. 2b). While we focus on the top 50 cm, due to the potential for vertical profiles of soil C to be affected by temperature22, very similar results were observed for the top 20 cm (Supplementary Fig. 2). In addition, the negative relationship between clay content and the effect of temperature on C storage was observed independently both above and below 15 oC (Fig. 3b, c).Fig. 2: Texture effects on temperature–soil carbon storage relationships.The effect of texture on the relationships between C storage in the top 50 cm of mineral soil and mean annual temperature in the raw data (a), and after accounting for potential confounding variables (b). The y-axes display the proportional reduction in C storage for each 10 oC increase in mean annual temperature, with higher values indicating greater reductions in soil C with temperature. In panel a, the slopes of the relationships (solid line), together with their 95% confidence intervals (dark grey shaded area), are presented for each of the textural categories, with the slope and 95% confidence interval for the full dataset (dotted line and light grey shaded areas) also presented across the graph for comparison. In panel b, the relationship between soil C storage and temperature after accounting for variation in annual precipitation (light blue), gross primary productivity (GPP; dark green), soil pH (purple), aridity (ET/PET; evapotranspiration minus potential evapotranspiration; navy blue), and cation exchange capacity (CEXC, light green) are shown. The slopes of these relationships (solid lines) together with their 95% confidence intervals (shaded area) are presented for each of the textural categories.Full size imageFig. 3: Comparison between soil profile data and JULES model output.The effect of texture on the relationships between C storage in the top 50 cm of mineral soil in the empirical data (solid lines) and JULES output (dashed lines). The slopes of these relationships (solid lines) together with their 95% confidence intervals (shaded area) are presented for each of the textural categories. Results for the full mean annual temperate range (a), as well as for subsets of the data for sites with mean annual temperatures below 15 oC (b, blue) and above 15 oC (c, red) are shown.Full size imageThe lower effect of temperature on soil C storage in fine-textured soils with greater stabilisation capacities was unexpected given the evidence of the high-temperature sensitivity associated with the decomposition of more protected SOM pools13,14,15. However, the findings from our global analysis are in agreement with a recent Europe-wide synthesis23, which, by compiling data from soil physical fractionation studies, demonstrated that mineral-associated C stocks varied less with temperate than freer particulate pools. Therefore, there is growing evidence that the effect of temperature on soil C storage is higher in soils containing a greater proportion of unprotected C.In the literature, there are apparently contradictory conclusions in terms of how C storage varies across fine-scale climate gradients, in which variation in other factors has been minimised. However, these may potentially be resolved by considering differences in the likely extent of SOM stabilisation. For example, on poorly weathered, relatively coarse-textured, silt loam soils in Alaska, mineral soil C stocks declined strongly with temperature24. In contrast, in Hawaiian forests growing on fine-textured soils with high concentrations of Al and Fe oxides, very little change in soil C storage was observed across a gradient of 5 oC in MAT. This was despite the fact that, in these Hawaiian forests, C storage in unprotected pools on the forest floor was found to decline strongly with temperature5. We suggest that differences in the extent of physicochemical protection in the Alaskan versus Hawaiian soils may explain the contrasting results. Thus, apparently contradictory findings may be resolvable within a single framework in which the relative effect of temperature on C storage in mineral soils declines as the soil’s physicochemical stabilisation capacity, and the proportion of C in protected pools, increase.Overall, our analysis identified C stored in high-latitude soils with limited capacities for stabilising organic matter as likely to be most vulnerable to the impacts of climate change. Such stores, therefore, may require particular attention given the high rates of warming taking place in cooler regions. In contrast, the particularly low effect of temperature on C storage in fine-textured soils in warm climates suggests (Fig. 3c) that the C stocks in many tropical soils may be less vulnerable to climate warming. While a soil warming study in a less weathered tropical soil identified the potential for high rates of C release25, our results are consistent with a recent large-scale analysis that concluded that the temperature sensitivity of soil respiration is generally lowest in tropical environments26. However, because C storage in tropical soils has been shown to be potentially vulnerable to drought27, it should not be concluded that C storage in tropical soils will be unaffected by climate change. Our results do, though, suggest that C stocks in coarse-textured soils at high latitudes are likely to be especially vulnerable to warming (Fig. 3b). Finally, while the dataset contains soil profile information for sites across the full mean annual temperature range investigated (0–30 oC), and there were data on a minimum of 500 profiles in every 5 oC temperature increment, increasing the amount of data available for sites with mean annual temperatures below 5 oC and greater than 20 oC would add further confidence to the findings.Because of their greater stabilisation capacities, fine-textured soils store more soil organic matter18. Therefore, fine and coarse-textured soils could contain similar absolute quantities of highly vulnerable C, and the lower effect of temperature in fine-textured soils could reflect the presence of greater quantities of low-vulnerability organic matter4. Therefore, it is likely to be very important to quantify stocks of unprotected pools, such as free particulate C, in soils with contrasting stabilisation capacities, and to investigate how such stocks vary with climate23. This may make it possible to identify whether there are still important stocks of unprotected organic matter that are vulnerable to climate warming in fine-textured soils with high stabilisation capacities2.Predicting and modelling future rates of C releaseAccurately predicting the response of soil C storage to global warming remains a major challenge. While spatial datasets, such as the ones analysed in this paper, add confidence to the prediction that C will be lost overall and help identify the most vulnerable stocks, they provide limited information on the likely rates or dynamics of C release. In this context, long-term surveys can be extremely valuable. For example, a recent study in Chinese grasslands was able to detect warming-induced soil C losses since the 1960s and, consistent with the global analysis presented here, coarser-textured soils lost far greater amounts of SOM28. Experimental soil warming studies also offer opportunities for further determining the factors controlling soil C storage and predicting rates of C release, although recent syntheses have produced conflicting overall findings29,30. Revisiting the networks of warming studies and considering the findings in the context of soil stabilisation capacities and changes in pools of protected and unprotected SOM may allow for a greater understanding of the observed patterns. For example, C losses from subsoils in response to 5 years of whole profile warming were shown to be dominated by the free particulate C pool31. Therefore, understanding the responses of different pools to warming may offer the potential to generate mechanistic understanding, even where changes in total C storage have been difficult to identify. It should though be recognised that there are major challenges in accurately quantifying relatively short-term changes in soil C stocks, and there are many other variables that differ between soil warming studies, including contrasting changes in plant productivity and rates of C input driven by interactions between C and nutrient cycling32. For these reasons, it may not always be possible to determine the role of soil stabilisation capacities in controlling soil C storage responses to experimental warming30, and observations collected across space and time will likely remain important for contextualising experimental results.Soil texture is included as a factor modifying decomposition rates in the terrestrial C cycle modules of a number of Earth systems models (ESMs), including the United Kingdom ESM (UKESM), whose land surface scheme (the Joint UK Land Environment Simulator (JULES)33) is based around the Rothamsted C model34. Therefore, we investigated whether JULES was already able to represent the patterns that we had observed in the empirical data. In direct contrast to the empirical data, JULES predicted very little variation in soil C with temperature in cooler regions (below 15 oC; Fig. 3b), but predicted a strong effect of temperature on C storage across all textural classes above 15 oC (Fig. 3c). The pattern across the full dataset was confounded by the model simulating only a small number of fine-textured soils at high latitudes (Fig. 3a), and the fact that the relationship between temperature and soil C storage differed so strongly above and below 15 oC. However, crucially, JULES failed to reproduce the greater effect of temperature on C storage in coarse-textured soils and overestimated the effect of temperature on C storage in fine-textured soils in warmer regions. These findings question whether JULES is identifying accurately which soil C stocks are most vulnerable to the effects of climate warming. This is important given the considerable geographical variation in (1) rates of climate warming and (2) the amounts of C stored in mineral soil horizons. In recent years, there have been major efforts made towards developing models that include physicochemical stabilisation mechanisms and yet can potentially be run at the global scale35,36,37. Testing whether such models can better simulate the observed spatial patterns of C storage in soils with contrasting stabilisation capacities would increase confidence in projections of future changes in soil C stocks38.Limitations and future perspectivesAs well as influencing rates of key biological processes, climatic variables also control pedogenesis, rates of mineral weathering and therefore influence the reactivity of soil surfaces26,27,39. Directly determining the binding affinity of mineral surfaces is challenging and would require detailed information on the type of clay minerals present, as well as the abundance of key metal oxides35,36,40, but there is, currently, insufficient data to assess these more detailed variables at the global scale35. However, it has been argued that, at broad spatial scales, soil pH may explain an important proportion of variation in mineral-binding affinities35,41. Furthermore, cation exchange capacity (CEXC) varies with the type of clay minerals present and the binding efficiencies of the mineral surfaces42. In global analyses, texture, pH and CEXC tend to be the three edaphic factors that correlate most strongly with soil C storage18,20. For these reasons, we also accounted for variation in both pH and CEXC, and evaluated whether the relationship between soil texture and the effect of temperature on C storage was retained. We found that it was (Fig. 2b). Thus, we conclude, that within this large dataset, clay content remains a strong predictor of soil stabilisation capacities, both overall, and after accounting for factors that potentially control SOM binding affinities.While we consider that our analysis of how SOM stabilisation capacities determine the effects of temperature on soil C storage is robust, it is also high level. Thus, there is considerable opportunity to further investigate different vulnerabilities of specific pools of SOM, contrasting the roles of mineral protection versus occlusion in aggregates7, determining the importance of SOM binding affinities40, and linking protection mechanisms with the sources of the organic matter (e.g. microbial versus plant derived43). A debate has often revolved around whether climatic versus edaphic factors are more important in controlling patterns of soil C storage. Rather, than focusing on which is more important, for predicting future rates of soil C release, we suggest that a key priority should be on identifying how key edaphic factors determine the vulnerability of contrasting soil C stocks to climate warming. In this context, a recent meta-analysis demonstrated the importance of soil properties in controlling the temperature sensitivity of soil respiration, emphasising how responses to global warming will likely vary substantially between different types of soils in contrasting geoclimatic zones26.Using a large global database, we observed declining C storage with temperature in mineral soils, suggesting that there is the potential for strong positive feedback to climate warming. Critically, however, this overall relationship masked differences between soils with contrasting C stabilisation capacities, as indicated by their textural properties. The data suggest that there are stabilised pools of SOM in fine-textured soils that may be relatively insensitive to the impacts of climate change, but that unprotected pools may be substantially more vulnerable to climate warming than currently predicted. Finally, because at least one major ESM was unable to reflect the observed patterns, we argue that ESMs should be evaluated against their ability to simulate the differences in the effects of temperature on C storage in soils with contrasting textural properties in order to reduce uncertainties in projections of the effect of climate change on future soil C storage. More

Tuber yieldPotato tuber yield increased gradually under 0 to 150 kg ha−1 of applied nitrogen (Fig. 1). Compared with the yield in N0, the yields in N60, N120 and N150 were greater by 16.1%, 21.5% and 67.9%, respectively, in 2017 and 18.2%, 27.4% and 44.9%, respectively, in 2018. However, at nitrogen rates of more than 150 kg ha−1, yield did not significantly differ. Furthermore, the fitted parabolic equation of each dataset from 2017 and 2018 showed maximum tuber yields of 19.7 and 20.4 t ha−1, respectively, where the nitrogen rates were 191 and 227 kg ha−1, respectively.Figure 1Potato tuber yield in treatments with different nitrogen fertilizer rates. Different letters indicate significant differences (P  More

1.Milshteyn, A., Schneider, J. S. & Brady, S. F. Mining the metabiome: identifying novel natural products from microbial communities. Chem. Biol. 21, 1211–1223 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
2.Rutledge, P. J. & Challis, G. L. Discovery of microbial natural products by activation of silent biosynthetic gene clusters. Nat. Rev. Microbiol. 13, 509–523 (2015).CAS 
PubMed 

Google Scholar 
3.Dittmann, E., Gugger, M., Sivonen, K. & Fewer, D. P. Natural product biosynthetic diversity and comparative genomics of the cyanobacteria. Trends Microbiol. 23, 642–652 (2015).CAS 
PubMed 

Google Scholar 
4.Cragg, G. M., Kingston, D. G. & Newman, D. J. Anticancer Agents from Natural Products. (CRC press, 2011).5.Cragg, G. M. & Newman, D. J. Natural products: a continuing source of novel drug leads. Biochimica et. Biophysica Acta 1830, 3670–3695 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
6.Singh, R., Kumar, M., Mittal, A. & Mehta, P. K. Microbial metabolites in nutrition, healthcare and agriculture. 3 Biotech 7, 15 (2017).PubMed 
PubMed Central 

Google Scholar 
7.Bohlmann, J. & Keeling, C. I. Terpenoid biomaterials. Plant J. 54, 656–669 (2008).CAS 
PubMed 

Google Scholar 
8.Kang, A. & Lee, T. S. Biotechnology for Biofuel Production and Optimization 35–71 (Elsevier, 2016).9.Nowruzi, B., Sarvari, G. & Blanco, S. The cosmetic application of cyanobacterial secondary metabolites. Algal Res. 49, 101959 (2020).
Google Scholar 
10.Crits-Christoph, A., Diamond, S., Butterfield, C. N., Thomas, B. C. & Banfield, J. F. Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis. Nature 558, 440–444 (2018).CAS 
PubMed 

Google Scholar 
11.Sharrar, A. M. et al. Bacterial secondary metabolite biosynthetic potential in soil varies with phylum, depth, and vegetation type. Mbio. 11, 1–17 (2020).
Google Scholar 
12.Libis, V. et al. Uncovering the biosynthetic potential of rare metagenomic DNA using co-occurrence network analysis of targeted sequences. Nat. Commun. 10, 1–9 (2019).CAS 

Google Scholar 
13.Haro-Moreno, J. M., López-Pérez, M. & Rodriguez-Valera, F. Enhanced recovery of microbial genes and genomes from a marine water column using long-read metagenomics. Front Microbiol. 2410, 1–15 (2021).
Google Scholar 
14.Navarro-Muñoz, J. C. et al. A computational framework to explore large-scale biosynthetic diversity. Nat Chem Biol. 16, 60–68 (2019).PubMed 
PubMed Central 

Google Scholar 
15.Sugimoto, Y. et al. A metagenomic strategy for harnessing the chemical repertoire of the human microbiome. Science 366, 1–11 (2019).
Google Scholar 
16.Katz, M., Hover, B. M. & Brady, S. F. Culture-independent discovery of natural products from soil metagenomes. J. Ind. Microbiol. Biotechnol. 43, 129–141 (2016).CAS 
PubMed 

Google Scholar 
17.Lynch, M. D. & Neufeld, J. D. Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol. 13, 217–229 (2015).CAS 
PubMed 

Google Scholar 
18.Amos, G. C. et al. Comparative transcriptomics as a guide to natural product discovery and biosynthetic gene cluster functionality. Proc. Natl Acad. Sci. USA 114, E11121–E11130 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
19.Rodriguez-Caballero, E. et al. Dryland photoautotrophic soil surface communities endangered by global change. Nat. Geosci. 11, 185–189 (2018).CAS 

Google Scholar 
20.Reddy, B. V. B. et al. Natural product biosynthetic gene diversity in geographically distinct soil microbiomes. Appl. Environ. Microbiol. 78, 3744–3752 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
21.Starkenburg, S. R. et al. Genome of the cyanobacterium Microcoleus vaginatus FGP-2, a photosynthetic ecosystem engineer of arid land soil biocrusts worldwide. J. Bacteriol. 193, 4569–4570 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
22.Belnap, J., Weber, B. & Büdel, B. Biological Soil Crusts: an Organizing Principle in Drylands 3–13 (Springer, 2016).23.Swenson, T. L., Karaoz, U., Swenson, J. M., Bowen, B. P. & Northen, T. R. Linking soil biology and chemistry in biological soil crust using isolate exometabolomics. Nat. Commun. 9, 1–10 (2018).CAS 

Google Scholar 
24.Adamek, M., Spohn, M., Stegmann, E. & Ziemert, N. Antibiotics 23–47 (Springer, 2017).25.Couradeau, E. et al. Bacteria increase arid-land soil surface temperature through the production of sunscreens. Nat. Commun. 7, 1–7 (2016).
Google Scholar 
26.Martins, T. P. et al. Chemistry, bioactivity and biosynthesis of cyanobacterial alkylresorcinols. Nat. Prod. Rep. 36, 1437–1461 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
27.Luesch, H., Moore, R. E., Paul, V. J., Mooberry, S. L. & Corbett, T. H. Isolation of dolastatin 10 from the marine cyanobacterium Symploca species VP642 and total stereochemistry and biological evaluation of its analogue symplostatin 1. J. Nat. Products 64, 907–910 (2001).CAS 

Google Scholar 
28.Chrapusta, E. et al. Microcystins and anatoxin-a in Arctic biocrust cyanobacterial communities. Toxicon 101, 35–40 (2015).CAS 
PubMed 

Google Scholar 
29.Van Goethem, M. W., Swenson, T. L., Trubl, G., Roux, S. & Northen, T. R. Characteristics of wetting-induced bacteriophage blooms in biological soil crust. Mbio. 10, 1–15 (2019).
Google Scholar 
30.Makhalanyane, T. P. et al. Microbial ecology of hot desert edaphic systems. FEMS Microbiol. Rev. 39, 203–221 (2015).CAS 
PubMed 

Google Scholar 
31.Bowker, M. A., Reed, S. C., Maestre, F. T. & Eldridge, D. J. Biocrusts: the living skin of the earth. (Springer, 2018).32.Büdel, B., Dulić, T., Darienko, T., Rybalka, N. & Friedl, T. Biological Soil Crusts: an Organizing Principle in Drylands 55–80 (Springer, 2016).33.Giraldo‐Silva, A., Nelson, C., Barger, N. N. & Garcia‐Pichel, F. Nursing biocrusts: isolation, cultivation, and fitness test of indigenous cyanobacteria. Restor. Ecol. 27, 793–803 (2019).
Google Scholar 
34.Koren, S. et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
35.Kolmogorov, M. et al. metaFlye: scalable long-read metagenome assembly using repeat graphs. Nat. Methods 17, 1103–1110 (2020).CAS 
PubMed 

Google Scholar 
36.Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
37.Mikheenko, A., Saveliev, V. & Gurevich, A. MetaQUAST: evaluation of metagenome assemblies. Bioinformatics 32, 1088–1090 (2016).CAS 
PubMed 

Google Scholar 
38.Howe, A. C. et al. Tackling soil diversity with the assembly of large, complex metagenomes. Proc. Natl Acad. Sci. USA 111, 4904–4909 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
39.Frank, J. A. et al. Improved metagenome assemblies and taxonomic binning using long-read circular consensus sequence data. Sci. Rep. 6, 1–10 (2016).
Google Scholar 
40.Hiraoka, S. et al. Metaepigenomic analysis reveals the unexplored diversity of DNA methylation in an environmental prokaryotic community. Nat. Commun. 10, 1–10 (2019).CAS 

Google Scholar 
41.Blin, K. et al. antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 47, W81–W87 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
42.Garcia-Pichel, F., Loza, V., Marusenko, Y., Mateo, P. & Potrafka, R. M. Temperature drives the continental-scale distribution of key microbes in topsoil communities. Science 340, 1574–1577 (2013).CAS 
PubMed 

Google Scholar 
43.Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410 (1990).CAS 

Google Scholar 
44.Nayfach, S. et al. A genomic catalog of Earth’s microbiomes. Nat. Biotechnol. 499–509 (2020).45.Fuqua, C. & Greenberg, E. P. Listening in on bacteria: acyl-homoserine lactone signalling. Nat. Rev. Mol. Cell Biol. 3, 685–695 (2002).CAS 
PubMed 

Google Scholar 
46.Ciemniecki, J. A. & Newman, D. K. The potential for redox-active metabolites to enhance or unlock anaerobic survival metabolisms in aerobes. J. Bacteriol. 202, 1–14 (2020).
Google Scholar 
47.Rajeev, L. et al. Dynamic cyanobacterial response to hydration and dehydration in a desert biological soil crust. ISME J. 7, 2178–2191 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
48.Nunes da Rocha, U. et al. Isolation of a significant fraction of non-phototroph diversity from a desert biological soil crust. Front. Microbiol. 6, 277 (2015).PubMed 
PubMed Central 

Google Scholar 
49.Felde, V. J. M. N. L., Peth, S., Uteau-Puschmann, D., Drahorad, S. & Felix-Henningsen, P. Soil microstructure as an under-explored feature of biological soil crust hydrological properties: case study from the NW Negev Desert. Biodivers. Conserv. 23, 1687–1708 (2014).
Google Scholar 
50.Bushnell, B. BBMap: a fast, accurate, splice-aware aligner. No. LBNL-7065E. (Lawrence Berkeley National Lab (LBNL), Berkeley, CA (United States), 2014).51.Donia, M. S., Ruffner, D. E., Cao, S. & Schmidt, E. W. Accessing the hidden majority of marine natural products through metagenomics. ChemBioChem. 12, 1230–1236 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
52.Van Der Maaten, L. Accelerating t-SNE using tree-based algorithms. J. Mach. Learn. Res. 15, 3221–3245 (2014).
Google Scholar 
53.Couradeau, E., Giraldo-Silva, A., De Martini, F. & Garcia-Pichel, F. Spatial segregation of the biological soil crust microbiome around its foundational cyanobacterium, Microcoleus vaginatus, and the formation of a nitrogen-fixing cyanosphere. Microbiome 7, 1–12 (2019).
Google Scholar 
54.Magne, F. et al. The firmicutes/bacteroidetes ratio: a relevant marker of gut dysbiosis in obese patients? Nutrients 12, 1474 (2020).CAS 
PubMed Central 

Google Scholar 
55.Baran, R. et al. Exometabolite niche partitioning among sympatric soil bacteria. Nat. Commun. 6, 8289 (2015).CAS 
PubMed 

Google Scholar 
56.Kupriyanova, E. V. et al. Extracellular β-class carbonic anhydrase of the alkaliphilic cyanobacterium Microcoleus chthonoplastes. J. Photochem. Photobiol. B: Biol. 103, 78–86 (2011).CAS 

Google Scholar 
57.Hernandez, M. & Newman, D. Extracellular electron transfer. Cell. Mol. Life Sci. 58, 1562–1571 (2001).CAS 
PubMed 

Google Scholar 
58.Karaoz, U. et al. Large blooms of Bacillales (Firmicutes) underlie the response to wetting of cyanobacterial biocrusts at various stages of maturity. MBio. 9, e01366–01316 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
59.Rodriguez-R, L. M., Gunturu, S., Tiedje, J. M., Cole, J. R. & Konstantinidis, K. T. Nonpareil 3: fast estimation of metagenomic coverage and sequence diversity. MSystems. 3, 1–9 (2018).
Google Scholar 
60.Kopylova, E., Noé, L. & Touzet, H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 28, 3211–3217 (2012).CAS 
PubMed 

Google Scholar 
61.Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 

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

Google Scholar 
63.Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
64.Van der Walt, A. J. et al. Assembling metagenomes, one community at a time. BMC Genom. 18, 1–13 (2017).
Google Scholar 
65.Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinforma. 11, 119 (2010).
Google Scholar 
66.Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).CAS 
PubMed 

Google Scholar 
67.Arkin, A. P. et al. KBase: the United States department of energy systems biology knowledgebase. Nat. Biotechnol. 36, 566 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
68.Kautsar, S. A. et al. MIBiG 2.0: a repository for biosynthetic gene clusters of known function. Nucleic Acids Res. 48, D454–D458 (2020).PubMed 

Google Scholar 
69.von Meijenfeldt, F. B., Arkhipova, K., Cambuy, D. D., Coutinho, F. H. & Dutilh, B. E. Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT. Genome Biol. 20, 217 (2019).
Google Scholar 
70.Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).PubMed 
PubMed Central 

Google Scholar 
71.Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).PubMed 
PubMed Central 

Google Scholar 
72.Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).PubMed 
PubMed Central 

Google Scholar  More

A model of diauxic community assemblyCommunity models studying diauxie should mimic serial dilution cultures instead of chemostats, in order to make their predictions both experimentally and ecologically relevant. Experimentally, microbial community assembly assays frequently utilize serial dilution cultures. Ecologically, diauxic growth is best suited to a “feast and famine” lifestyle, which a serial dilution culture mimics30,31,32. Therefore, throughout this manuscript, we model the assembly of a microbial community undergoing a sequence of growth-dilution cycles (see Fig. 1a). Community assembly occurs gradually through the addition of microbial species from a diverse species pool one at a time. Each species in the pool consumes resources diauxically, i.e., one at a time according to its resource preference.Fig. 1: Model of community assembly with diauxie and serial dilution.a Tables of growth rates and resource preferences of two species α (red) and β (yellow), each capable of consuming all four available resources, R1 to R4. The resource preference sets the sequence in which a microbial species utilizes resources, and the corresponding rates gXi indicate the growth rate while consuming each resource (see “Methods”). b Diauxic growth curve of species α during one serial dilution cycle, which has 4 phases of growth on each individual resource, with rates gα1, gα3, gα2, and gα4, respectively (with a brief lag period between two phases). At the end of each dilution cycle, we dilute the population by a factor D = 100, and supply fresh resources (see “Methods”). c Resource depletion curves corresponding to (b), where each resource is represented by a different color. R1 is exhausted at time T1; then species α consumes R3 which runs out at T3, which is followed by exhaustion of R2 at T2, and so on. d Schematic of serial dilution experiment. During community assembly, new species are added one by one from a species pool. After each successful invasion, the system undergoes several growth-dilution cycles until it reaches a steady state. e Population dynamics corresponding to the assembly process in (d). Panels (b) and (c) correspond to a small section of this process (highlighted in gray), where the community dynamics consist only of species α (red) reaching a steady state.Full size imageWe begin by illustrating the growth of a single species (labeled α) grown in an environment with four resources (Fig. 1a–c). The species first grows on its most preferred resource (R1) with a growth rate gα1 until time T1, when this resource gets exhausted. After a lag period τ, the species switches to growing on its next preferred resource (R3) with growth rate gα3 until time T3, when this resource also gets exhausted. This process of diauxic growth by sequential utilization of resources continues until either all resources are depleted, or the cycle ends at time T. At this point, a fraction 1/D of the medium containing the species is transferred to a fresh medium replete with resources. This corresponds to the dilution of species abundances by a factor D, mimicking serial dilution experiments in the laboratory.After several transfers, species dynamics converge to a steady state, where each species starts a cycle with the same initial abundance as the previous cycle. At this point, we add a small population of a new invader species, chosen randomly from the species pool, to the steady-state community (Fig. 1d, e). (Hence, we assume that species invasions are rare enough such that communities always reach a steady state before the next invasion.) The invader may differ from the resident species in both resource preference order and growth rates on each resource (Fig. 1a). Once introduced, the invader may grow and establish itself in the community in a new steady state (Fig. 1d, e), or it may fail, returning the community to its previous steady state.The growth rates and preference orders completely characterize a species, while the set of resource depletion times (T1, T2, etc.) characterize the current state of the abiotic environment. As we will later show, these resource depletion times are important observables in a community, since they determine the success or failure of an invader.A realistic example of a community captured by our model is the human gut microbiome, specifically the assembly of primary consumers (e.g., Bacteroides species) on the polysaccharides (e.g., starch, cellulose, and mucin) that they consume. Here, there is a significant overlap between the metabolic capabilities of the microbes, but they nevertheless coexist. These species often consume polysaccharides diauxically, and engage in resource competition. Moreover, several of these species have different resource preferences, which others have hypothesized help them coexist26,33.Throughout this paper, we neglect diauxic lag times (τ = 0) for simplicity. We will later show that adding lag times only quantitatively strengthens our main results (see “Discussion” and Fig. 5). We also assume that the supplied resource concentrations are sufficiently large, enabling species to always grow exponentially at their resource-specific growth rates. Further, we assume a balanced supply of resources, i.e., that resources are supplied in equal concentrations (see “Discussion” and Supplementary Text for results in an unbalanced resource supply).We simulated the assembly of 1000 communities, each being colonized from a pool of ~10,000 species (see “Methods”). Species could utilize all 4 supplied resources diauxically. Each species had a random resource preference order and different growth rates on each resource, which were picked randomly from a rectified normal distribution (with mean 0.25 and standard deviation 0.05). We assumed that the growth rate distributions for each of the 4 resources were the same, such that no resource was consistently better than the other. This is a simplifying assumption, but it nevertheless captures a variety of experimental observations showing remarkable growth rate variability of different microbial species on the same carbon sources34,35,36. Community assembly proceeded via introduction of species one at a time, in a random order, until each species had attempted to invade exactly once.Emergent properties of diauxic community assemblyTo study the emergent properties of communities of diauxic species, we followed the assembly process from a species pool via invasion of species one at a time. We used the number of invasion attempts to track time; communities matured over successive invasions. We found that the assembly process became slower over time—successful invasions became rarer as the community matured (Fig. 2a inset). Throughout the assembly process, we recorded four key properties of communities: total resource depletion time, species diversity, complementarity of the community, and prevalence of anomalous species (defined below).Fig. 2: Emergent properties of diauxic microbial communities.In all plots, solid bold lines represent the average over 958 individual community assembly simulations, while gray lines correspond to 100 randomly chosen community assembly simulations. a Total resource depletion time during community assembly (the time taken by the community to deplete all available resources). (Inset) Number of successful invasions during community assembly. b Total species diversity during community assembly (number of surviving species at steady state). c Resource utilization complementarity during community assembly. For each time point, the nth choice complementarity was calculated as a number of unique resources among the n-th preferred choices of all species in the community, divided by the number of unique resources in the environment. For a certain community, the null expectation (complementarity without selection) was defined by the complementarity of a random set of species from the pool that has the same diversity of that community. Colored lines show the average trend of complementarity on each preferred resource choice: top (light blue), second (cyan), third (deep green), and fourth (light green). The red dash-dotted line shows the average trend of null expectation. The gray dash-dotted line at the top corresponds to the perfect complementarity, which is 1. d Frequency of species with anomalous resource preferences during community assembly. The gray dash-dotted line is the expectation of fraction of anomalous species (75%) in the pools.Full size imageResource depletion timeIn each community, resources disappear at specific times and in a well-defined order (Fig. 1c). The total resource depletion time measures how quickly the community consumes all supplied resources. In this way, the total resource depletion time characterizes the overall speed at which a community consumes resources. The total resource depletion time decreased as communities assembled (Fig. 2a, solid line). The rate and degree of this decrease depend on the mean and variance of the growth rate distribution and the number of invasion attempts. In addition, the variability in depletion times between communities reduced over community assembly (Fig. 2a, gray lines; coefficient of variation reduces by 47%, see Fig. S1). Thus the assembly process selects for communities that collectively consume resources quickly.Species diversityThe species diversity was quantified as the number of species coexisting in the steady-state community. In the model, like in other consumer-resource models, the number of coexisting species at steady state is limited by the number of resources, 4 (Fig. 2b, dashed line)10,37. This is a natural consequence of competition for resources in our model (see Supplementary Text, section F for a derivation). Notably, species with the same resource preferences can coexist in the model, as long as the number of species is less than the number of resources (e.g., pairs of E. coli strains can coexist in media with glucose and xylose, see below). We found that the average community diversity increased over time, but the rate slowed as the community matured (Fig. 2b; note the logarithmic x-axis scale). Communities displayed significant variability in the trajectories of increasing diversity (Fig. 2b, gray lines). We discuss the slow increase of diversity, and observed variability, in the next section.Top choice complementarityThe top choice complementarity of a community measured the overlap in the top choice resource of each of the species residing in the community. We defined the top choice complementarity of a community as the number of unique top choice resources among community residents, divided by the number of residents. Thus the top choice complementarity varied between 1, in a maximally complementary community where each resident species had a unique top choice resource (Fig. 2c, right), and 1 divided by the number of coexisting species in the community, where all residents chose the same resource as the top choice (Fig. 2c, left). During community assembly, top choice complementarity stayed close to the maximum value throughout the assembly process (Fig. 2c, blue). This observation was in sharp contrast to the prediction from a null model for the complementarity (Fig. 2c, red). We obtained the null prediction by measuring the complementarity of a group of randomly chosen species from the species pool (group size being the number of coexisting species in the community). This null prediction decreased during the assembly process, due to the increasing community diversity, unlike the top choice complementarity which remained close to the maximum value. We also recorded the complementarity in the second, third, and fourth choice resource of the assembled community (defined similarly to the top choice complementarity). The complementarity of all other choices agreed with the null prediction (Fig. 2c). Together, these observations suggest that communities of coexisting diauxic species exhibit high complementarity on the top-choice resources, in a manner reminiscent of niche partitioning in consumer-resource models.Prevalence of anomalous speciesIntuition gleaned from experiments with E. coli dictates that microbes often grow fastest on their top choice resource (glucose for E. coli)18,20. However, exceptions to this trend also exist, such as Bacteroides species in the human gut that often prefer polysaccharides that they grow slower on22,26,38. Based on this intuition, we defined anomalous microbes as microbes that do not grow fastest on their top choice resource. To investigate which resource preferences might give microbes a competitive advantage during community assembly, we tracked the fraction of anomalous resident species during community assembly. Despite the majority (75%) of species in the pool being anomalous (since growth rates and preferences were randomly picked; see “Methods”), anomalous species were absent in mature communities. The fraction of anomalous resident species decreased rapidly during assembly (Fig. 2d). Thus, anomalous resource preferences are strongly selected against during community assembly. Further investigation revealed a reduced selection pressure against anomalous species if either resource supply was severely imbalanced (i.e., the imbalance has to be comparable to the dilution factor, D = 100), or if the dilution factor was small (see Figs. S4 and S5; also see Supplementary Text, sections C and H). However, microbes with anomalous resource preferences were eventually outcompeted in all conditions.Top choice resources chiefly drive emergent assembly patternsTo understand what factors drove the maintenance of top choice complementarity—despite the steady increase in species diversity, expected to reduce complementarity—we focused on growth on top choice resources. We hypothesized that the reason for the much higher than expected top choice complementarity was the following: diauxic species derived most of their growth, and spent most of their time growing on their top choice resources. Co-utilizing microbes, instead, grow on multiple resources simultaneously, spending roughly equal time on each utilized resource.To test this hypothesis, we first simulated the growth of a single diauxic species in monoculture using our model. We found that indeed, the species derived the overwhelming majority of its growth (measured in generations of growth) and spent most of its time growing on its top choice resource (54%, Fig. 3a, b, left). For a simpler case, where a single species had the same growth rate g while growing on two resources (both supplied at the same concentration), and preferring resource R1 over R2, we derived the ratio of time spent growing on the top choice resource R1 (T1) versus the second choice R2 (T2 − T1). We obtained the following approximate expression for a large dilution factor D (see Supplementary Text, section A):$$frac{{T}_{1}}{{T}_{2}-{T}_{1}}=frac{,{{mbox{log}}}(D/2)}{{{mbox{log}}},(2)},$$
(1)
which explains that the fraction of time spent growing on the top choice resource increases with the dilution factor.Strikingly, the fraction of time spent growing on the top choice resource became even larger if the species grown in monoculture (Fig. 3b, top row) were instead part of a diverse community (i.e., in co-culture with 3 other species, top choice share 70% versus 54% in monoculture, Fig. 3b, bottom row and top row, respectively). This is because of the following reason. In our model, while a species consumes and grows on all available resources in monoculture, in co-culture, it may not have the opportunity to consume all the resources it can grow on because other species might deplete them first. This further skews growth in favor of the top choice resource. Such a phenomenon only occurs in diauxic species, not co-utilizing species (Supplementary Text, section I).Invader successInterestingly, once we understood that the top choice chiefly drove species growth, we could explain the other emergent patterns in diauxic communities. Importantly, the success of an invader depended on the growth rate on their top choice resource. As community assembly proceeded, the top choice growth rate of successful invaders increased consistently (Fig. 3c, blue line), while their growth rates on all other choices remained constant and close to the average growth rate (Fig. 3c, green lines). Selection on the top choice growth rate in diauxic communities is in striking contrast with co-utilizing communities, which we found select for the average growth rate across all resources instead (Supplementary Text, section I). Further, an invader whose top choice resource coincided with the last depleted resource in the community had the highest probability of invasion success (Fig. 3d). Invaders whose top choice resource was not depleted last had lesser time to grow on it, and thus a lower rate of invasion success. By depleting the last resource faster, invaders reduced the total resource depletion time in the community, thus explaining the trend observed in Fig. 2a. In addition, after a successful invasion, the community’s steady state could have a different resource depletion order.Complementarity and diversitySuccessful invasions could be classified into one of two types based on the “invaded resource”, i.e., the invader’s top choice. If the invaded resource was not the top choice of any other resident community member, we called it an invasion of an “unoccupied” resource (Fig. 3e; in our simulations, 33% of cases). If the invaded resource was instead already the top choice of at least one resident, we called it an invasion of an “occupied” resource (Fig. 3e; 67% of cases). Both types of successful invasions had different effects on species diversity, but interestingly, both maintained complementarity (on the top choice, as in Fig. 2c). invasions of unoccupied resources usually increased community diversity by 1 (62% of cases), and were less likely to result in the extinction of one or more other species (38% of cases). This is because, in that case, the invader did not have to compete with other residents for its top choice resource. For communities with a complementarity {T}_{2},$$
(2)
where gα1 and gα2 are the species α’s growth rates on R1 and R2, respectively. The two triangular regions separated by the diagonal define two complementary scenarios: when T1  T2, R1 is depleted second and the species grows on R1 after R2 is depleted.For a given set of initial resource and species concentrations, community dynamics must converge to a steady state lying on the ZNGI of the surviving species (e.g., the bold purple point in Fig. 4a). This point defines the resource depletion times by the resident species at steady state. Changing the resource supply or dilution factor moves this point along the ZNGI.The ZNGI of a species also separates the resource environment space into two regions: a region inside the ZNGI (towards the origin) where that species grows by a factor D. An invader is successful if it is able to grow by a factor ≥D in the community it invades. Geometrically, the invader’s ZNGI must be closer to the origin than the resource environment corresponding to the invaded community (Fig. 4b). In this way, our geometric approach allows easy visualization of invasion criteria.We can also visualize invasion outcomes. A successful invasion of a single-species community leads to either displacement of the resident or coexistence between the invader and resident. For example, in Fig. 4b, because the ZNGI of the invader (blue) lies fully inside the ZNGI of the resident (purple), the invader displaces the resident. This is because the invader reduces the resource depletion times in the environment to a point where the resident can no longer survive, driving it extinct (bold blue point in Fig. 4b). In contrast, in Fig. 4c, the ZNGI of the invader (orange) intersects with the new resident (blue), in a manner that leads to coexistence between both species (albeit at a new set of resource depletion times, i.e., their intersection point in Fig. 4c). In general, whether two species will coexist depends on various factors, such as the supplied resource concentrations, but whenever two species coexist, they will do so at the intersection of their ZNGIs (Supplementary Text, section A). As a corollary, two species whose ZNGIs do not intersect cannot coexist. Notably, the orange and blue species in Fig. 4c coexist stably with each other; a short perturbation to the resource supply is quickly compensated by species growth, and the resource depletion times returned to the coexistence point (see Supplementary Text, section B for details).The geometric approach provides an alternative explanation to why species with complementary top choices are more likely to coexist than species with the same top choice (Fig. 2c). The ZNGIs of species sharing the same top choice are unlikely to intersect with each other (e.g., the blue and purple species in Fig. 4b). This is because of two reasons: (1) their segments in the yellow region are parallel to each other since both species prefer the same resource (R2), and (2) for the slanted segments in the green region to intersect, the blue species would need a higher growth rate on R1 than the purple species. This is as likely as the outcome of a coin toss, since both growth rates derive from the same distribution. Thus, an invasion of an occupied resource often leads to displacement of the resident, not coexistence (Fig. 4b, d) and no change in community diversity, while an invasion of an unoccupied resource often leads to coexistence (Fig. 4c) and an increase in community diversity (Fig. 3e).Fig. 3: Top choice resources chiefly drive community diversity and complementarity.a (top) Table showing the preferences of a diauxic microbial species (purple) for 4 resources, R1 to R4. (bottom) Plots showing the depletion of the 4 resources by the purple species during one serial dilution cycle, when grown alone in our model. b (top) Bar plots showing the time taken by the purple species in (a) to grow on each of the 4 resources. Percentages on each bar represent the fraction of time spent growing on each resource. (bottom) Bar plots showing the number of generations grown, or the number of doublings by the species when growing on each resource. In both cases, the plots on the left show the quantities when the purple species is in monoculture (growing alone), and those on the right show them when the purple species is in a community with 3 other species. c Mean growth rates of successful invaders during community assembly. The blue line corresponds to the invader’s top choice, while the other colors correspond to all other choices. The horizontal dashed line shows the mean growth rate of the species pool. Each quantity represents a moving average from 958 independent community assembly simulations. Error bars represent s.e.m. d Fraction of the successful invasions as a function of the order in which the invader’s top choice resource is depleted, 1 indicating cases where the invader prefers the earliest depleted resource, and 4 where it prefers the last depleted resource. Each bar represents the mean of such a fraction over 958 independent community assembly simulations, and error bars represent s.e.m. e, f Effect of invasions of community diversity and complementarity, based on whether the invader’s top choice was (e) unoccupied or (f) occupied. Cartoons show the typical effect of an invasion. Pie charts show the fraction of invasions that increase, decrease or maintain a community’s species diversity (middle) and complementarity (right). On unoccupied resources, diversity typically increases (62%), but sometimes stays constant (32%) or decreases (6%). On occupied resources, diversity typically stays the same (68%), but sometimes decreases (26%) and rarely increases (6%). In almost all cases complementarity either stays maintained or increases ( >95%), and very rarely decreases ( More

1.IPCC Special Report on Global Warming of 1.5 °C (eds Masson-Delmotte, V. et al.) (WMO, 2018).2.Anderson, C. M. et al. Natural climate solutions are not enough. Science 363, 933–934 (2019).CAS 

Google Scholar 
3.Carton, W., Lund, J. F. & Dooley, K. Undoing equivalence: rethinking carbon accounting for just carbon removal. Front. Clim. 3, 30 (2021).
Google Scholar 
4.Seddon, N. et al. Getting the message right on nature-based solutions to climate change. Glob. Change Biol. 27, 1518–1546 (2021).
Google Scholar 
5.Griscom, B. W. et al. We need both natural and energy solutions to stabilize our climate. Glob. Change Biol. 25, 1889–1890 (2019).
Google Scholar 
6.Fargione, J. et al. Natural climate solutions for the United States. Sci. Adv. 4, eaat1869 (2018).
Google Scholar 
7.Drever, C. R. et al. Natural climate solutions for Canada. Sci. Adv. 7, eabd6034 (2021).CAS 

Google Scholar 
8.Griscom, B. W. et al. Natural climate solutions. Proc. Natl Acad. Sci. USA 114, 11645–11650 (2017).CAS 

Google Scholar 
9.Fuss, S. et al. Negative emissions—Part 2: costs, potentials and side effects. Environ. Res. Lett. 13, 63002 (2018).
Google Scholar 
10.Gregorio, N. et al. in Enhancing Food Security Through Forest Landscape Restoration: Lessons from Burkina Faso, Brazil, Guatemala, Viet Nam, Ghana, Ethiopia and Philippines (eds Kumar, C. et al.) 174–217 (IUCN, 2015).11.Meyer, J. M. Gifford Pinchot, John Muir, and the boundaries of politics in American thought. Polity 30, 267–284 (1997).
Google Scholar 
12.Standard on Biodiversity Offsets (BBOP, 2012).13.Performance Standard 6: Biodiversity Conservation and Sustainable Management of Natural Resources (IFC, 2012).14.Arlidge, W. N. S. et al. A global mitigation hierarchy for nature conservation. Bioscience 68, 336–347 (2018).
Google Scholar 
15.Science-Based Targets for Nature: Initial Guidance for Business (Science Based Targets Network, 2020).16.Seddon, N. et al. Understanding the value and limits of nature-based solutions to climate change and other global challenges. Phil. Trans. R. Soc. B 375, 20190120 (2020).
Google Scholar 
17.Ellis, P. W. et al. Reduced-impact logging for climate change mitigation (RIL-C) can halve selective logging emissions from tropical forests. Ecol. Manag. 438, 255–266 (2019).
Google Scholar 
18.Martin, D. M. Ecological restoration should be redefined for the twenty-first century. Restor. Ecol. 25, 668–673 (2017).
Google Scholar 
19.Veldman, J. W. et al. Where tree planting and forest expansion are bad for biodiversity and ecosystem services. Bioscience 65, 1011–1018 (2015).
Google Scholar 
20.Supporting Canadians and Fighting COVID-19 (Department of Finance Canada, 2020).21.Roe, S. et al. Contribution of the land sector to a 1.5 °C world. Nat. Clim. Change 9, 817–828 (2019).
Google Scholar 
22.Seddon, N. et al. Nature-Based Solutions in Nationally Determined Contributions: Synthesis and Recommendations for Enhancing Climate Ambition and Action by 2020 (IUCN, 2019).23.Carbon Removal Corporate Action Tracker (Institute for Carbon Removal Law and Policy, accessed 6 July 2021); https://research.american.edu/carbonremoval/2020/05/07/carbon-removal-corporate-action-tracker/24.Pendrill, F. et al. Agricultural and forestry trade drives large share of tropical deforestation emissions. Glob. Environ. Change 56, 1–10 (2019).
Google Scholar 
25.Goal 1 Assessment: Striving to End Natural Forest Loss (NYDF Progress Assessment Secretariat, 2020).26.Smith, B. One year later: The path to carbon negative—a progress report on our climate ‘moonshot’. Microsoft Blog (28 January 2021); https://blogs.microsoft.com/blog/2021/01/28/one-year-later-the-path-to-carbon-negative-a-progress-report-on-our-climate-moonshot/27.Ward, C. et al. Smallholder perceptions of land restoration activities: rewetting tropical peatland oil palm areas in Sumatra. Indonesia. Reg. Environ. Change 21, 1 (2020).
Google Scholar 
28.Jacobson, M. & Ham, C. The (un)broken promise of agroforestry: a case study of improved fallows in Zambia. Environ. Dev. Sustain. 22, 8247–8260 (2020).
Google Scholar 
29.West, T. A. P., Börner, J., Sills, E. O. & Kontoleon, A. Overstated carbon emission reductions from voluntary REDD+ projects in the Brazilian Amazon. Proc. Natl Acad. Sci USA 117, 24188–24194 (2020).CAS 

Google Scholar 
30.Cook-Patton, S. C. et al. Lower cost and more feasible options to restore forest cover in the contiguous United States for climate mitigation. One Earth 3, 739–752 (2020).
Google Scholar 
31.Petersen, S. O., Højberg, O., Poulsen, M., Schwab, C. & Eriksen, J. Methanogenic community changes, and emissions of methane and other gases, during storage of acidified and untreated pig slurry. J. Appl. Microbiol. 117, 160–172 (2014).CAS 

Google Scholar 
32.Günther, A. et al. Prompt rewetting of drained peatlands reduces climate warming despite methane emissions. Nat. Commun. 11, 1644 (2020).
Google Scholar 
33.Qin, Z. et al. Delayed impact of natural climate solutions. Glob. Change Biol. 27, 215–217 (2021).
Google Scholar 
34.Anderegg, W. R. L. et al. Climate-driven risks to the climate mitigation potential of forests. Science 368, eaaz7005 (2020).CAS 

Google Scholar 
35.Pagiola, S., Honey-Rosés, J. & Freire-González, J. Assessing the permanence of land-use change induced by payments for environmental services: evidence from Nicaragua. Trop. Conserv. Sci. https://doi.org/10.1177/1940082920922676 (2020).36.Tseng, T.-W. J. et al. Influence of land tenure interventions on human well-being and environmental outcomes. Nat. Sustain. 4, 242–251 (2021).
Google Scholar 
37.Smith, P. et al. Impacts of land-based greenhouse gas removal options on ecosystem services and the United Nations Sustainable Development Goals. Annu. Rev. Environ. Resour. 44, 255–286 (2019).
Google Scholar 
38.Nunez, S., Verboom, J. & Alkemade, R. Assessing land-based mitigation implications for biodiversity. Environ. Sci. Policy 106, 68–76 (2020).
Google Scholar 
39.Chausson, A. et al. Mapping the effectiveness of nature-based solutions for climate change adaptation. Glob. Change Biol. 26, 6134–6155 (2020).
Google Scholar 
40.Infield, M., Entwistle, A., Anthem, H., Mugisha, A. & Phillips, K. Reflections on cultural values approaches to conservation: lessons from 20 years of implementation. Oryx 52, 220–230 (2018).
Google Scholar 
41.Rosenstock, T. S. et al. A planetary health perspective on agroforestry in sub-Saharan Africa. One Earth 1, 330–344 (2019).
Google Scholar 
42.Garrett, H. E. et al. Hardwood silvopasture management in North America. Agrofor. Syst. 61, 21–33 (2004).
Google Scholar 
43.Kroeger, T. et al. Returns on investment in watershed conservation: application of a best practices analytical framework to the Rio Camboriú Water Producer program, Santa Catarina, Brazil. Sci. Total Environ. 657, 1368–1381 (2019).CAS 

Google Scholar 
44.Lamb, D., Erskine, P. D. & Parrotta, J. A. Restoration of degraded tropical forest landscapes. Science 310, 1628–1632 (2005).CAS 

Google Scholar 
45.Ferreira, J. et al. Carbon-focused conservation may fail to protect the most biodiverse tropical forests. Nat. Clim. Change 8, 744–749 (2018).CAS 

Google Scholar 
46.Cook-Patton, S. C. et al. Mapping carbon accumulation potential from global natural forest regrowth. Nature 585, 545–550 (2020).CAS 

Google Scholar 
47.Wilson, S. J., Schelhas, J., Grau, R., Nanni, A. S. & Sloan, S. Forest ecosystem-service transitions: the ecological dimensions of the forest transition. Ecol. Soc. 22, 38 (2017).
Google Scholar 
48.Funk, J. M. et al. Securing the climate benefits of stable forests. Clim. Policy 19, 845–860 (2019).
Google Scholar 
49.Keith, H. et al. Evaluating nature-based solutions for climate mitigation and conservation requires comprehensive carbon accounting. Sci. Total Environ. 769, 144341 (2021).CAS 

Google Scholar 
50.Moomaw, W. R., Masino, S. A. & Faison, E. K. Intact forests in the United States: proforestation mitigates climate change and serves the greatest good. Front. For. Glob. Change 2, 27 (2019).
Google Scholar 
51.Hiraishi, T. et al. 2013 Supplement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories: Wetlands (WMO, 2013).52.Goldstein, A. et al. Protecting irrecoverable carbon in Earth’s ecosystems. Nat. Clim. Change 10, 287–295 (2020).CAS 

Google Scholar 
53.Griscom, B. W. et al. National mitigation potential from natural climate solutions in the tropics. Phil. Trans. R. Soc. B 375, 20190126 (2020).CAS 

Google Scholar 
54.Busch, J. et al. Potential for low-cost carbon dioxide removal through tropical reforestation. Nat. Clim. Change 9, 463–466 (2019).CAS 

Google Scholar 
55.Vargas Zeppetello, L. R. et al. Large scale tropical deforestation drives extreme warming. Environ. Res. Lett. 15, 84012 (2020).
Google Scholar 
56.Spalding, M. D. et al. The role of ecosystems in coastal protection: adapting to climate change and coastal hazards. Ocean Coast. Manag. 90, 50–57 (2014).
Google Scholar 
57.Watson, J. E. M. et al. The exceptional value of intact forest ecosystems. Nat. Ecol. Evol. 2, 599–610 (2018).
Google Scholar 
58.Global Assessment Report on Biodiversity and Ecosystem Services (IPBES, 2019).59.Dobson, A. P. et al. Ecology and economics for pandemic prevention. Science 369, 379–381 (2020).CAS 

Google Scholar 
60.Streck, C. REDD+ and leakage: debunking myths and promoting integrated solutions. Clim. Policy 21, 843–852 (2021).
Google Scholar 
61.Brancalion, P. H. S. et al. The cost of restoring carbon stocks in Brazil’s Atlantic Forest. Land Degrad. Dev. 32, 830–841 (2021).
Google Scholar 
62.Bustamante-Sánchez, M. A. & Armesto, J. J. Seed limitation during early forest succession in a rural landscape on Chiloé Island, Chile: implications for temperate forest restoration. J. Appl. Ecol. 49, 1103–1112 (2012).
Google Scholar 
63.Koch, A., Brierley, C. & Lewis, S. L. Effects of Earth system feedbacks on the potential mitigation of large-scale tropical forest restoration. Biogeosciences 18, 2627–2647 (2021).CAS 

Google Scholar 
64.Zickfeld, K., Azevedo, D., Mathesius, S. & Matthews, H. D. Asymmetry in the climate–carbon cycle response to positive and negative CO2 emissions. Nat. Clim. Change 11, 613–617 (2021).CAS 

Google Scholar 
65.Johnson, K. A. et al. A benefit–cost analysis of floodplain land acquisition for US flood damage reduction. Nat. Sustain. 3, 56–62 (2020).
Google Scholar 
66.Nolte, C. High-resolution land value maps reveal underestimation of conservation costs in the United States. Proc. Natl Acad. Sci. USA 117, 29577–29583 (2020).CAS 

Google Scholar 
67.Reetz, H., Heffer, P. & Bruulsema, T. in Managing Water and Fertilizer for Sustainable Agricultural Intensification (eds Drechsel, P. et al.) 65–86 (IFA, IWMI, IPNI and IPI, 2015).68.Sharma, P. et al. The role of cover crops towards sustainable soil health and agriculture—a review paper. Am. J. Plant Sci. 09, 1935–1951 (2018).CAS 

Google Scholar 
69.Bergeron, M. et al. Reduced soil nutrient leaching following the establishment of tree-based intercropping systems in eastern Canada. Agrofor. Syst. 83, 321–330 (2011).
Google Scholar 
70.Moore, A. A. & Palmer, M. A. Invertebrate biodiveristy in agricultural and urban headwater streams: implications for conservation and management. Ecol. Appl. 15, 1169–1177 (2005).
Google Scholar 
71.Martin, M. P. et al. People plant trees for utility more often than for biodiversity or carbon. Biol. Conserv. 261, 109224 (2021).
Google Scholar 
72.Mendes, T. P., de Assis Montag, L. F., Alvarado, S. T. & Juen, L. Assessing habitat quality on alpha and beta diversity of Odonata larvae (Insect) in logging areas in Amazon forest. Hydrobiologia 848, 1147–1161 (2021).
Google Scholar 
73.Crouzeilles, R. et al. Achieving cost-effective landscape-scale forest restoration through targeted natural regeneration. Conserv. Lett. 13, e12709 (2020).
Google Scholar 
74.Gilroy, J. J. et al. Cheap carbon and biodiversity co-benefits from forest regeneration in a hotspot of endemism. Nat. Clim. Change 4, 503–507 (2014).
Google Scholar 
75.Riahi, K. et al. The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: an overview. Glob. Environ. Change 42, 153–168 (2017).
Google Scholar 
76.Strassburg, B. B. N. et al. Strategic approaches to restoring ecosystems can triple conservation gains and halve costs. Nat. Ecol. Evol. 3, 62–70 (2019).
Google Scholar 
77.Taillardat, P., Thompson, B. S., Garneau, M., Trottier, K. & Friess, D. A. Climate change mitigation potential of wetlands and the cost-effectiveness of their restoration. Interface Focus 10, 20190129 (2020).
Google Scholar 
78.Xu, S., Liu, X., Li, X. & Tian, C. Soil organic carbon changes following wetland restoration: a global meta-analysis. Geoderma 353, 89–96 (2019).CAS 

Google Scholar 
79.Holl, K. D. & Brancalion, P. H. S. Tree planting is not a simple solution. Science 368, 580–581 (2020).CAS 

Google Scholar 
80.Kroeger, T., McDonald, R. I., Boucher, T., Zhang, P. & Wang, L. Where the people are: current trends and future potential targeted investments in urban trees for PM10 and temperature mitigation in 27 U.S. cities. Landsc. Urban Plan. 177, 277–240 (2018).
Google Scholar 
81.McDonald, R. I., Kroeger, T., Zhang, P. & Hamel, P. The value of US urban tree cover for reducing heat-related health impacts and electricity consumption. Ecosystems 23, 137–150 (2020).
Google Scholar 
82.Heris, M. et al. Piloting urban ecosystem accounting for the United States. Ecosyst. Serv. 48, 101226 (2021).
Google Scholar 
83.Ellison, D. et al. Trees, forests and water: cool insights for a hot world. Glob. Environ. Change 43, 51–61 (2017).
Google Scholar 
84.Li, R. et al. Time and space catch up with restoration programs that ignore ecosystem service trade-offs. Sci. Adv. 7, eabf8650 (2021).
Google Scholar 
85.Gaveau, D. L. A. et al. Rise and fall of forest loss and industrial plantations in Borneo (2000–2017). Conserv. Lett. 12, e12622 (2019).
Google Scholar 
86.Griscom, B. W., Goodman, R. C., Burivalova, Z. & Putz, F. E. Carbon and biodiversity impacts of intensive versus extensive tropical forestry. Conserv. Lett. 11, e12362 (2018).
Google Scholar 
87.Gabon’s Proposed National RED+ Forest Reference Level (Gabonese Republic, 2021).88.Umunay, P., Gregoire, T., Gopalakrishna, T., Ellis, P. & Putz, F. Selective logging emissions and potential emission reductions from reduced-impact logging in the Congo Basin. Ecol. Manag. 437, 360–371 (2019).
Google Scholar 
89.Natural Climate Solutions World Atlas (Nature4Climate, accessed 9 December 2020); https://nature4climate.org/n4c-mapper/90.Dave, R. et al. Second Bonn Challenge Progress Report: Application of the Barometer in 2018 (IUCN, 2019); https://doi.org/10.2305/IUCN.CH.2019.06.en91.Kremen, C. & Merenlender, A. M. Landscapes that work for biodiversity and people. Science 362, eaau6020 (2018).
Google Scholar 
92.Seymour, F. Seeing the forests as well as the (trillion) trees in corporate climate strategies. One Earth 2, 390–393 (2020).
Google Scholar 
93.Kronenberg, J. & Mieszkowicz, J. Planting trees for publicity—how much are they worth? Sustainability 3, 1022–1034 (2011).
Google Scholar 
94.Microsoft Carbon Removal: Lessons from an Early Corporate Purchase (Microsoft, 2021).95.Toor, I. A., Smith, E. G., Whalen, J. K. & Naseem, A. Tree-based intercropping in southern Ontario, Canada. Can. J. Agric. Econ. 60, 141–154 (2012).
Google Scholar 
96.Watson, J. E. M., Dudley, N., Segan, D. B. & Hockings, M. The performance and potential of protected areas. Nature 515, 67–73 (2014).CAS 

Google Scholar 
97.zu Ermgassen, S. O. S. E. et al. The ecological outcomes of biodiversity offsets under “no net loss” policies: a global review. Conserv. Lett. 12, e12664 (2019).
Google Scholar 
98.Adoption of the Paris Agreement FCCC/CP/2015/L.9/Rev.1 (UNFCCC, 2015). More