Ecology

1.Mayr, E. The species concept: Semantics versus semantics. Evolution 3, 371–372 (1949).Article 

Google Scholar 
2.Hortal, J. et al. Seven shortfalls that beset large-scale knowledge of biodiversity. Annu. Rev. Ecol. Evol. Syst. 46, 523–549 (2015).Article 

Google Scholar 
3.Mace, G. M. The role of taxonomy in species conservation. Philos. Trans. R. Soc. Lond. B Biol. Sci. 359, 711–719 (2004).PubMed 
PubMed Central 
Article 

Google Scholar 
4.Gustafson, K. D., Kensinger, B. J., Bolek, M. G. & Luttbeg, B. Distinct snail (Physa) morphotypes from different habitats converge in shell shape and size under common garden conditions. Evol. Ecol. Res. 16, 77–89 (2014).
Google Scholar 
5.Aksenova, O. V. et al. Species richness, molecular taxonomy and biogeography of the radicine pond snails (Gastropoda: Lymnaeidae) in the Old World. Sci. Rep. 8, 1–7 (2018).ADS 
CAS 
Article 

Google Scholar 
6.Liu, H. P. & Hershler, R. A new species and range extensions for three other species of pebblesnails (Lithoglyphidae, Fluminicola) from the upper Klamath basin, California-Oregon. ZooKeys 812, 47–67 (2019).Article 

Google Scholar 
7.Alda, P. et al. Systematics and geographical distribution of Galba species, a group of cryptic and worldwide freshwater snails. Mol. Phylogenet. Evol. 157, 107035 (2021).PubMed 
Article 

Google Scholar 
8.Taylor, D. W. Introduction to Physidae (Gastropoda: Hygrophila); biogeography, classification, morphology. Rev. Biol. Trop. 51(Supplement 1), 1–287 (2003).CAS 
PubMed 
PubMed Central 

Google Scholar 
9.Wethington, A. R. & Lydeard, C. A molecular phylogeny of Physidae (Gastropoda: Basommatophora) based on mitochondrial DNA sequences. J. Molluscan Stud. 73, 241–257 (2007).Article 

Google Scholar 
10.Ng, T. H. et al. Molluscs for sale: assessment of freshwater gastropods and bivalves in the ornamental pet trade. PLoS ONE 11, e0161130 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
11.Saito, T., Prozorova, L., Hirano, T., Fukuda, H. & Chiba, S. Endangered freshwater limpets in Japan are actually alien invasive species. Conserv. Genet. 19, 947–958 (2018).Article 

Google Scholar 
12.Lydeard, C., Campbell, D. & Golz, M. Physa acuta Draparnaud, 1805 should be treated as a native of North America, not Europe. Malacologia 59, 347–350 (2016).Article 

Google Scholar 
13.Albrecht, C., Kroll, O., Terrazas, E. M. & Wilke, T. Invasion of ancient Lake Titicaca by the globally invasive Physa acuta (Gastropoda: Pulmonata: Hygrophila). Biol. Invasions 11, 1821–1826 (2009).Article 

Google Scholar 
14.Ng, T. H., Tan, S. K. & Yeo, D. C. Clarifying the identity of the long-established, globally-invasive Physa acuta Draparnaud, 1805 (Gastropoda: Physidae) in Singapore. BioInvasions Rec. 4, 189–194 (2015).Article 

Google Scholar 
15.Collado, G. A. Unraveling cryptic invasion of a freshwater snail in Chile based on molecular and morphological data. Biodivers. Conserv. 26, 567–578 (2017).Article 

Google Scholar 
16.Johnson, P. D. et al. Conservation status of freshwater gastropods of Canada and the United States. Fisheries 38, 247–282 (2013).Article 

Google Scholar 
17.Strong, E. E. & Whelan, N. V. Assessing the diversity of western North American Juga (Semisulcospiridae, Gastropoda). Mol. Phylogenet. Evol. 136, 87–103 (2019).PubMed 
Article 
PubMed Central 

Google Scholar 
18.Hebert, P. D., Ratnasingham, S. & De Waard, J. R. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. R. Soc. Lond. B 270(supplement 1), S96-99 (2003).CAS 

Google Scholar 
19.Stöger, I. & Schrödl, M. Mitogenomics does not resolve deep molluscan relationships (yet?). Mol. Phylogenet. Evol. 69, 376–392 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
20.Cunha, T. J. & Giribet, G. A congruent topology for deep gastropod relationships. Proc. R. Soc. B 286, 20182776 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
21.Varney, R. M. et al. Assessment of mitochondrial genomes for heterobranch gastropod phylogenetics. BMC Ecol. Evol. 21, 6 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
22.Remigio, E. A. & Hebert, P. D. Testing the utility of partial COI sequences for phylogenetic estimates of gastropod relationships. Mol. Phylogenet. Evol. 29, 641–647 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
23.Collins, R. A. & Cruickshank, R. H. The seven deadly sins of DNA barcoding. Mol. Ecol. Resour. 13, 969–975 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
24.Ratnasingham, S. & Hebert, P. D. A DNA-based registry for all animal species: The Barcode Index Number (BIN) system. PLoS ONE 8, e66213 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
25.Whelan, N. V. & Strong, E. E. Morphology, molecules and taxonomy: Extreme incongruence in pleurocerids (Gastropoda, Cerithioidea, Pleuroceridae). Zoolog. Scr. 45, 62–87 (2016).Article 

Google Scholar 
26.Razkin, O., Gómez-Moliner, B. J., Vardinoyannis, K., Martínez-Ortí, A. & Madeira, M. J. Species delimitation for cryptic species complexes: Case study of Pyramidula (Gastropoda, Pulmonata). Zool. Scr. 46, 55–72 (2017).Article 

Google Scholar 
27.Liu, H. P., Hershler, R. & Hovingh, P. Molecular evidence enables further resolution of the western North American Pyrgulopsis kolobensis complex (Caenogastropoda: Hydrobiidae). J. Molluscan Stud. 84, 103–107 (2018).Article 

Google Scholar 
28.Ward, R. D. DNA barcode divergence among species and genera of birds and fishes. Mol. Ecol. Resour. 9, 1077–1085 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
29.Saadi, A. J., Davison, A. & Wade, C. M. Molecular phylogeny of freshwater snails and limpets (Panpulmonata: Hygrophila). Zool. J. Linn. Soc. 190, 518–531 (2020).Article 

Google Scholar 
30.Frest, T. J. & Johannes, E. J. An annotated checklist of Idaho land and freshwater mollusks. J. Idaho Acad. Sci. 36(2), 1–51 (2000).
Google Scholar 
31.Pip, E. & Franck, J. P. Molecular phylogenetics of central Canadian Physidae (Pulmonata: Basommatophora). Can. J. Zool. 86, 10–16 (2008).CAS 
Article 

Google Scholar 
32.Tariel, J., Plénet, S. & Luquet, É. Transgenerational plasticity of inducible defences: Combined effects of grand-parental, parental and current environments. Ecol. Evol. 10, 2367–2376 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
33.Perrin, N. The life history parameters of Physa acuta (Gastropoda, Mollusca) in experimental conditions. Revue Suisse Zoologique 93, 725–736 (1986).Article 

Google Scholar 
34.Taylor, D. W. New species of Physa (Gastropoda: Hygrophila) from the western United States. Malacol. Rev. 21, 43–79 (1988).
Google Scholar 
35.U.S. Fish and Wildlife Service. Determination of endangered or threatened status for five aquatic snails in south central Idaho. Fed. Reg. 57, 59242–59257 (1992).
Google Scholar 
36.Rogers, D. C. & Wethington, A. R. Physa natricina Taylor 1988, junior synonym of Physa acuta Draparnaud, 1805 (Pulmonata: Physidae). Zootaxa 1662, 45–51 (2007).
Google Scholar 
37.Gates, K. K., Kerans, B. L., Keebaugh, J. L., Kalinowski, S. & Vu, N. Taxonomic identity of the endangered Snake River physa, Physa natricina (Pulmonata: Physidae) combining traditional and molecular techniques. Conserv. Genet. 14, 159–169 (2013).Article 

Google Scholar 
38.Moore, A. C., Burch, J. B. & Duda, T. F. Recognition of a highly restricted freshwater snail lineage (Physidae: Physella) in southeastern Oregon: Convergent evolution, historical context, and conservation considerations. Conserv. Genet. 16, 113–123 (2015).Article 

Google Scholar 
39.Dillon, R. T., Robinson, J. D. & Wethington, A. R. Empirical estimates of reproductive isolation among the freshwater pulmonate snails Physa acuta, P. pomilia, and P. hendersoni. Malacologia 49, 283–292 (2007).Article 

Google Scholar 
40.De Queiroz, K. Species concepts and species delimitation. Syst. Biol. 56, 879–886 (2007).PubMed 
Article 
PubMed Central 

Google Scholar 
41.Dyke, A. S., Moore, A. & Robertson, L. Deglaciation of North America. Geological Survey of Canada Open File 1574 (2003).42.Wethington, A. R., Wise, J. & Dillon, R. T. Jr. Genetic and morphological characterization of the Physidae of South Carolina (Gastropoda: Pulmonata: Basommatophora), with description of a new species. Nautilus 123, 282–292 (2009).
Google Scholar 
43.Ebbs, E. T., Loker, E. S. & Brant, S. V. Phylogeography and genetics of the globally invasive snail Physa acuta Draparnaud 1805, and its potential to serve as an intermediate host to larval digenetic trematodes. BMC Evol. Biol. 18, 1–7 (2018).Article 

Google Scholar 
44.Duggan, I. C. The freshwater aquarium trade as a vector for incidental invertebrate fauna. Biol. Invasions 12, 3757–3770 (2010).Article 

Google Scholar 
45.Van Leeuwen, C. H. et al. How did this snail get here? Several dispersal vectors inferred for an aquatic invasive species. Freshw. Biol. 58, 88–99 (2013).Article 

Google Scholar 
46.Coughlan, N. E., Kelly, T. C., Davenport, J. & Jansen, M. A. Up, up and away: Bird-mediated ectozoochorous dispersal between aquatic environments. Freshw. Biol. 62, 631–648 (2017).Article 

Google Scholar 
47.Bony, Y. K. et al. Ecological conditions for spread of the invasive snail Physa marmorata (Pulmonata: Physidae) in the Ivory Coast. Afr. Zool. 43, 53–60 (2008).Article 

Google Scholar 
48.Pierce, K. L. & Morgan, L. A. Is the track of the Yellowstone hotspot driven by a deep mantle plume?—Review of volcanism, faulting, and uplift in light of new data. J. Volcanol. Geotherm. Res. 188, 1–25 (2009).ADS 
CAS 
Article 

Google Scholar 
49.Smith, G. R. et al. Biogeography and timing of evolutionary events among Great Basin fishes. In Great Basin Aquatic Systems History. Smithsonian Contributions to the Earth Sciences Vol. 33 (eds Hershler, R. et al.) 175–234 (Smithsonian Institution Press, 2002).
Google Scholar 
50.Oviatt, C. G. Chronology of Lake Bonneville, 30,000 to 10,000 yr BP. Quatern. Sci. Rev. 110, 166–171 (2015).Article 

Google Scholar 
51.Safran, E. B. et al. Plugs or flood-makers? The unstable landslide dams of eastern Oregon. Geomorphology 248, 237–251 (2015).ADS 
Article 

Google Scholar 
52.Ely, L. L. et al. Owyhee River intracanyon lava flows: Does the river give a dam?. GSA Bull. 124, 1667–1687 (2012).Article 

Google Scholar 
53.Matthews, J. et al. Rapid range expansion of the invasive quagga mussel in relation to zebra mussel presence in the Netherlands and western Europe. Biol. Invasions 16, 23–42 (2014).Article 

Google Scholar 
54.Folmer, O., Black, M., Hoeh, W., Lutz, R. & Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechol. 3, 294–299 (1994).CAS 

Google Scholar 
55.Geller, J., Meyer, C., Parker, M. & Hawk, H. Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Mol. Ecol. Resour. 13, 851–861 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
56.Uit de Weerd, D. R. & Gittenberger, E. Phylogeny of the land snail family Clausiliidae (Gastropoda: Pulmonata). Mol. Phylogenet. Evol. 67, 201–216 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
57.Katoh, K., Rozewicki, J. & Yamada, K. D. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 20, 1160–1166 (2019).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
58.Nixon, K. C. & Wheeler, Q. D. An amplification of the phylogenetic species concept. Cladistics 6, 211–223 (1990).Article 

Google Scholar 
59.Galtier, N. Delineating species in the speciation continuum: A proposal. Evol. Appl. 12, 657–663 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
60.DeSalle, R., Egan, M. G. & Siddall, M. The unholy trinity: Taxonomy, species delimitation and DNA barcoding. Philos. Trans. R. Soc. B Biol. Sci. 360, 1905–1916 (2005).CAS 
Article 

Google Scholar 
61.Bouchet, P. et al. Revised classification, nomenclator and typification of gastropod and monoplacophoran families. Malacologia 61, 1–526 (2017).Article 

Google Scholar 
62.Wethington, A. R. & Guralnick, R. Are populations of physids from different hot springs distinctive lineages?. Am. Malacol. Bull. 19, 135–144 (2004).
Google Scholar 
63.Huang, Y., Niu, B., Gao, Y., Fu, L. & Li, W. CD-HIT Suite: A web server for clustering and comparing biological sequences. Bioinformatics 26, 680–682 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
64.Nguyen, L. T., Schmidt, H. A., Von Haeseler, A. & Minh, B. Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
Article 

Google Scholar 
65.Puillandre, N., Brouillet, S. & Achaz, G. ASAP: Assemble species by automatic partitioning. Mol. Ecol. Resour. 21, 609–620 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
66.Puillandre, N., Lambert, A., Brouillet, S. & Achaz, G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol. Ecol. 21, 1864–1877 (2012).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
67.Delicado, D., Arconada, B., Aguado, A. & Ramos, M. A. Multilocus phylogeny, species delimitation and biogeography of Iberian valvatiform springsnails (Caenogastropoda: Hydrobiidae), with the description of a new genus. Zool. J. Linn. Soc. 186, 892–914 (2019).Article 

Google Scholar 
68.Kapli, T. et al. Multi-rate Poisson tree processes for single-locus species delimitation under maximum likelihood and Markov chain Monte Carlo. Bioinformatics 33, 1630–1638 (2016).
Google Scholar 
69.Clement, M., Posada, D. C. & Crandall, K. A. TCS: A computer program to estimate gene genealogies. Mol. Ecol. 9, 1657–1659 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
70.Hart, M. W. & Sunday, J. Things fall apart: Biological species form unconnected parsimony networks. Biol. Lett. 3, 509–512 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
71.Minh, B. Q., Nguyen, M. A. & von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 30, 1188–1195 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
72.Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
73.Meier, R., Zhang, G. & Ali, F. The use of mean instead of smallest interspecific distances exaggerates the size of the “barcoding gap” and leads to misidentification. Syst. Biol. 57, 809–813 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
74.Dellicour, S. & Flot, J. F. The hitchhiker’s guide to single-locus species delimitation. Mol. Ecol. Resour. 18, 1234–1246 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
75.Avise, J. C. Phylogeography: The History and Formation of Species (Harvard University Press, 2000).Book 

Google Scholar 
76.Dinapoli, A., Tamer, C., Franssen, S., Naduvilezhath, L. & Klussmann-Kolb, A. Utility of H3-gene sequences for phylogenetic reconstruction—a case study of heterobranch Gastropoda. Bonner Zoologische Beiträge 55(3/4), 191–202 (2006).
Google Scholar 
77.Ayyagari, V. S. & Sreerama, K. Molecular phylogenetic analysis of Pulmonata (Mollusca: Gastropoda) on the basis of histone-3 gene. Beni-Suef Univ. J. Basic Appl. Sci. 8, 1–8 (2019).Article 

Google Scholar  More

Here we show, for the first time, that the toxicity of a pesticide formulation to bees is caused exclusively by a co-formulant (alcohol ethoxylates), rather than the active ingredient. A 0.8 µL acute oral dose of the agricultural fungicide formulation Amistar® caused a range of damage to bees: both lethal, with 23% mortality, and sublethal, with 45% reduced sucrose consumption, 3.8% drop in body weight (whereas the negative control gained 4.8%), and a 302% increase in gut melanisation. For all metrics tested, the Amistar® and alcohol ethoxylates treatments were not statistically different, demonstrating conclusively that the toxicity of the formulation, Amistar®, to bumble bees is driven by the alcohol ethoxylates. These results demonstrate gaps in the regulatory system and highlight the need for a greater research focus on co-formulants.The mortality in the Amistar® treatment, and treatments containing alcohol ethoxylates reached 32% at its highest, which is substantial given that bees are likely to have a high level of exposure to Amistar® and alcohol ethoxylates. The mechanism by which the alcohol ethoxylates cause mortality has not been explicitly isolated, but our results suggest two potential, possibly related, causes. We recorded a 302% increase in the melanised area of bee midguts in the alcohol ethoxylates treatment. A similar effect was observed in Melipona scutellaris exposed to the pure fungicide active ingredient pyraclostrobin alongside a similar reduction in survival37. We suggest that the alcohol ethoxylates are disrupting the structure of the midgut, which the bee immune system is reacting to with melanisation44 (see Fig. 5). In parallel with this gut damage, alcohol ethoxylate treatment drove a 54% reduction in sugar consumption, which persisted throughout the experiment. Supplementary Fig. S3 shows a plot comparing sugar consumption against gut melanisation, with increasing gut melanisation correlated to reduced sugar consumption in the Amistar®, co-formulant mixture and alcohol ethoxylates treatments. Consequently, we propose that mortality was driven by energy depletion due to reduced consumption, which in turn may have been driven by damage to the gut.Figure 5(Left) Bumble bee midgut in the negative control treatment. (Right) Bumble bee midgut in the co-formulant mixture treatment, which contains alcohol ethoxylates. The dark brown patches are areas of melanisation, indicative of damage to the gut. Both bees survived the full 120 h.Full size imageLikely as a consequence of the reduced consumption of sucrose, bumble bees in the alcohol ethoxylates treatment lost 8.4% of their original weight, in stark contrast to the negative control where bees gained 4.8% over the five-day period. This indicates the alcohol ethoxylate treated bees were expending more energy than they were consuming, and thus exhibiting a negative energy balance. This weight loss, while considerable as a percentage of the bee’s total body mass, is also similar in scale to the weight of the sucrose bees consume in one sitting (EA Straw pers. obs.), for which rigorous data do not exist. As such it is possible that a portion of the weight loss is attributable to the reduced sucrose consumption of the bees, meaning they would have less sucrose in their guts at the time of weighing. Sucrose consumption does not, however, explain the failure of alcohol ethoxylate treated bees to gain weight, which was observed in the control treatment. The weight loss, and lack of weight gain, are concerning because they are likely to indicate a reduction in fat reserves, although this has not been experimentally confirmed. Bee fat reserves are important physiologically, in particular in responding to immune threats45,46. Fat reserves allow bees the energetic resources to buffer against challenges, and thus their depletion could expose bees to greater risk from future threats47.The reduced appetite and negative energy balance in alcohol ethoxylates treated bees could have broader effects in the natural environment. Bees pollinate flowers as they forage for nectar and pollen, so a reduction in their appetite could subsequently have effects on ecosystem services. In our experiment, bumble bee appetite was reduced immediately after ingesting a single dose of alcohol ethoxylates or Amistar®. This effect persisted for five days after exposure, indicating a persistent change in consumption behaviour. While nectar-foraging in bumble bees is driven by the needs of the colony48, a reduction in appetite would reduce overall colony nectar consumption, and thus the number of foraging trips made for nectar. Fewer visits to flowers for nectar may lead to reduced pollination, which would be detrimental to crop yields and farm profits. Further studies of how the impacts we have found map onto foraging and pollination are clearly needed. Importantly, the reduction in appetite recorded in our experiment is a sublethal effect, which standard lower tier testing would not detect. When Amistar® is tested on bumble bees for the 2025 renewal of azoxystrobin, this sublethal effect will be missed by regulatory testing, despite the impact it may have on the pollination services such testing is designed to protect. We suggest that a simple modification to the regulatory protocol OECD 247 would be to weigh the sucrose syringes at the start and end of the trials to calculate sucrose consumption, which would allow measurement of this sublethal effect with minimal additional workload.Our results show a slightly, but not significantly, higher level of mortality in the alcohol ethoxylates treatment (30%) than the Amistar® treatment (23%). If this is a real biological difference, one explanation might be that the concentration of alcohol ethoxylates in the Amistar® formulation was lower than that used in the alcohol ethoxylates treatment solution. This is possible because the Amistar® material safety data sheet lists concentrations as a range (10–20% for alcohol ethoxylates), and here we used the upper end of the range. The co-formulant mixture treatment in all metrics was statistically indistinguishable from the alcohol ethoxylates treatment, showing that the toxicity of alcohol ethoxylates is not a result of synergism with other co-formulants.We believe that the implications of our results are not limited to a laboratory setting and a single species, as other published and unpublished research supports our findings. Semi-field flight cage experiments, where Amistar® was applied to a crop, found effects on full bumble bee colonies (Bombus terrestris). Amistar® caused a reduction in average bee weight and a reduction in foraging activity, as our results predict49. This demonstrates that the effects observed in our laboratory testing scale up to effects at a field realistic level. Additionally, in honeybees (Apis mellifera) Amistar® has been found to cause mortality in laboratory experiments at a range of doses50,51, demonstrating the mortality effect found in our experiment is not species specific. However, no mortality was seen in trials on the red mason bee Osmia bicornis (Hellström and Paxton, unpublished data). Additionally, a similar compound, C11 and lower alcohol ethoxylates, has been found in small scale laboratory testing to cause 100% mortality after contact exposure in honeybees31.To measure the exposure of bees to PPP’s, the EU mandates trials that measure chemical residues in pollen and nectar after crops have been sprayed with either active ingredients or formulations34. However, these residue analysis studies only measure active ingredient concentrations, not the co-formulants. As such, we have no systematic data on the exposure of bees to co-formulants7,8,9. This dearth of data means that the exposure of bees to co-formulants is very poorly characterised. To estimate exposure to alcohol ethoxylates, residue data for Amistar®’s active ingredient azoxystrobin could be used as a proxy18,52. However, the chemical properties of alcohol ethoxylates, specifically their surfactant action, make it unlikely that they have an equivalent environmental fate to azoxystrobin, so this would not be appropriate.While we have very little data to quantify bee exposure to alcohol ethoxylates, we know Amistar® can be applied to crops, such as strawberries, during flowering while bees are foraging on them. The Environmental Information Sheet for Amistar® states “[For bees] no risk management is necessary. Amistar® is of low risk to honey bees”53,54,55. In addition, we would note that exposure of bees to alcohol ethoxylates, and related substances, is not exclusively from Amistar®. For example, a cursory search of the Syngenta website56 immediately identified alcohol ethoxylates in five other Syngenta products. Worryingly, the chemical group alcohol ethoxylates sit in, alkoxylated alcohols, are also widely used in adjuvants, which are products which can be added to tank mixtures to modify the action of the agrichemical6. 89 adjuvant products licenced in the UK contain alkoxylated alcohols as the primary ingredient15. To our knowledge, these adjuvants have never been toxicity tested on bees and have no bee exposure mitigation measures in place whatsoever.To complement measures to promote academic research, moving regulatory research beyond its mortality and active ingredient-centric approach to toxicity testing would better reflect the risks pesticides, as used in the field, pose. For regulatory systems to accurately characterise risk they need to estimate the scale of sublethal effects, regardless of initial mortality results33. The results presented here demonstrate that even substances assessed by regulators as ‘bee safe’ can pose a serious hazard to bee health. To reflect potential sublethal differences caused by co-formulation composition, all formulations could undergo a much more rigorous set of lower tier testing or be automatically entered for higher tier testing.In the face of declining bee populations we advocate that a precautionary approach minimising the exposure of bees to potential stressors, where possible, would be prudent. The current legislation allowing application of PPPs directly onto bees and flowering plants does not align with the emerging evidence that co-formulants, adjuvants, herbicides and fungicides can be hazardous to bees8,57. The wealth of untested and undisclosed co-formulants used abundantly in agriculture is a serious and pressing concern for the health of pollinators worldwide. More

1.Cubasch, U. et al. Climate Change 2013: the physical science basis. In Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Vol 1 (eds Stocker, T. F. et al.) 119–158 (Cambridge University Press, 2013).
Google Scholar 
2.Jackson, R. B. et al. Increasing anthropogenic methane emissions arise equally from agricultural and fossil fuel sources. Environ. Res. Lett. https://doi.org/10.1088/1748-9326/ab9ed2 (2020).Article 

Google Scholar 
3.Saunois, M. et al. The global methane budget 2000–2017. Earth Syst. Sci. Data 12, 1561–1623. https://doi.org/10.5194/essd-12-1561-2020 (2020).ADS 
Article 

Google Scholar 
4.Hobson, P. N. & Stewart, C. S. The rumen Microbial Ecosystem (Blackie Academic & Professional, 1997).Book 

Google Scholar 
5.Henderson, G. et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5, 14567. https://doi.org/10.1038/srep14567 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
6.Li, Y. et al. The complete genome sequence of the methanogenic archaeon ISO4-H5 provides insights into the methylotrophic lifestyle of a ruminal representative of the Methanomassiliicoccales. Stand. Genom. Sci. 11, 59. https://doi.org/10.1186/s40793-016-0183-5 (2016).Article 

Google Scholar 
7.Lang, K. et al. New mode of energy metabolism in the seventh order of methanogens as revealed by comparative genome analysis of “Candidatus Methanoplasma termitum”. Appl. Environ. Microbiol. 81, 1338–1352. https://doi.org/10.1128/AEM.03389-14 (2015).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
8.Hoehler, T., Losey, N. A., Gunsalus, R. P. & McInerney, M. J. In Biogenesis of Hydrocarbons (eds Stams, A. & Sousa, D.) 1–26 (Springer, 2018).
Google Scholar 
9.Neill, A. R., Grime, D. W. & Dawson, R. M. C. Conversion of choline methyl groups through trimethylamine into methane in the rumen. Biochem. J. 170, 529–535. https://doi.org/10.1042/bj1700529 (1978).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
10.Erdman, R. A. & Sharma, B. K. Effect of dietary rumen-protected choline in lactating dairy cows. J. Dairy Sci. 74, 1641–1647. https://doi.org/10.3168/jds.S0022-0302(91)78326-4 (1991).CAS 
Article 
PubMed 

Google Scholar 
11.Sharma, B. K. & Erdman, R. A. Effects of dietary and abomasally infused choline on milk production responses of lactating dairy cows. J. Nutr. 119, 248–254. https://doi.org/10.1093/jn/119.2.248 (1989).CAS 
Article 
PubMed 

Google Scholar 
12.Soliva, C. & Hess, H. In Measuring Methane Production from Ruminants: Measuring Methane Emission of Ruminants by In Vitro and In Vivo Techniques (eds Makkar, H. P. & Vercoe, P. E.) 15–31 (Springer, 2007).Chapter 

Google Scholar 
13.Craciun, S. & Balskus, E. P. Microbial conversion of choline to trimethylamine requires a glycyl radical enzyme. PNAS 109, 21307–21312. https://doi.org/10.1073/pnas.1215689109 (2012).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
14.Jameson, E. et al. Anaerobic choline metabolism in microcompartments promotes growth and swarming of Proteus mirabilis. Environ. Microbiol. 18, 2886–2898. https://doi.org/10.1111/1462-2920.13059 (2016).CAS 
Article 
PubMed 

Google Scholar 
15.Herring, T. I., Harris, T. N., Chowdhury, C., Mohanty, S. K. & Bobik, T. A. A bacterial microcompartment is used for choline fermentation by Escherichia coli 536. J. Bacteriol. 200, e00764-e817. https://doi.org/10.1128/JB.00764-17 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
16.EFSA. Scientific Opinion on safety and efficacy of choline chloride as a feed additive for all animal species. EFSA J. 9, 2353 (2011).
Google Scholar 
17.Lewis, D. J. Ammonia toxicity in the ruminant. J. Agric. Sci. 55, 111–117 (1960).CAS 
Article 

Google Scholar 
18.Hogan, J. P. Absorption of ammonia through rumen of sheep. Aust. J. Biol. Sci. 14, 448–450. https://doi.org/10.1071/Bi9610448 (1961).CAS 
Article 

Google Scholar 
19.Sprott, G. D. & Patel, G. B. Ammonia toxicity in pure cultures of methanogenic bacteria. Syst. Appl. Microbiol. 7, 358–363 (1986).CAS 
Article 

Google Scholar 
20.Lewis, D. Ammonia toxicity in the ruminant. J. Agric. Sci. 55(1), 111–117 (1960).CAS 
Article 

Google Scholar 
21.Ungerfeld, E. M., Rust, S. R. & Burnett, R. Increases in microbial nitrogen production and efficiency in vitro with three inhibitors of ruminal methanogenesis. Can. J. Microbiol. 53, 496–503. https://doi.org/10.1139/W07-008 (2007).CAS 
Article 
PubMed 

Google Scholar 
22.Lundgren, B. R., Sarwar, Z., Pinto, A., Ganley, J. G. & Nomura, C. T. Ethanolamine catabolism in Pseudomonas aeruginosa PAO1 is regulated by the enhancer-binding protein EatR (PA4021) and the alternative sigma factor RpoN. J. Bacteriol. 198, 2318–2329. https://doi.org/10.1128/JB.00357-16 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
23.Rychlik, J. L., LaVera, R. & Russell, J. B. Amino acid deamination by ruminal Megasphaera elsdenii strains. Curr. Microbiol. 45, 340–345. https://doi.org/10.1007/s00284-002-3743-4 (2002).CAS 
Article 
PubMed 

Google Scholar 
24.Park, K. & Lee, H. Effects of nitrogen gas flushing in comparison with argon on rumen fermentation characteristics in in vitro studies. J. Anim. Sci. Technol. 62, 52–57. https://doi.org/10.5187/jast.2020.62.1.52 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
25.Hobson, P. N., Summers, R., Postgate, J. R. & Ware, D. A. Nitrogen fixation in the rumen of a living sheep. J. Gen. Microbiol. 77, 225–226. https://doi.org/10.1099/00221287-77-1-225 (1973).CAS 
Article 
PubMed 

Google Scholar 
26.Harada, N., Nishiyama, M. & Matsumoto, S. Inhibition of methanogens increases photo-dependent nitrogenase activities in anoxic paddy soil amended with rice straw. FEMS Microbiol. Ecol. 35, 231–238. https://doi.org/10.1111/j.1574-6941.2001.tb00808.x (2001).CAS 
Article 
PubMed 

Google Scholar 
27.Haaker, H. & Klugkist, J. The bioenergetics of electron transport to nitrogenase. J FEMS Microbiol. Lett. 46, 57–71 (1987).CAS 
Article 

Google Scholar 
28.Edgren, T. & Nordlund, S. The fixABCX genes in Rhodospirillum rubrum encode a putative membrane complex participating in electron transfer to nitrogenase. J. Bacteriol. 186, 2052–2060 (2004).CAS 
Article 

Google Scholar 
29.Igai, K. et al. Nitrogen fixation and nifH diversity in human gut microbiota. Sci. Rep. 6, 31942. https://doi.org/10.1038/srep31942 (2016).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
30.Ungerfeld, E. M. Shifts in metabolic hydrogen sinks in the methanogenesis-inhibited ruminal fermentation: A meta-analysis. Front. Microbiol. 6, 37. https://doi.org/10.3389/fmicb.2015.00037 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
31.Leahy, S. C. et al. The complete genome sequence of Methanobrevibacter sp. AbM4. Stand. Genom. Sci. 8, 215–227. https://doi.org/10.4056/sigs.3977691 (2013).CAS 
Article 

Google Scholar 
32.Hoedt, E. C. et al. Differences down-under: Alcohol-fueled methanogenesis by archaea present in Australian macropodids. ISME J. 10, 2376–2388. https://doi.org/10.1038/ismej.2016.41 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
33.Ungerfeld, E. M. & Kohn, R. A. In Ruminant Physiology: Digestion, Metabolism and Impact of Nutrition on Gene Expression, Immunology and Stress (eds Sejrsen, K. et al.) 55–85 (Wageningen Academic Publishers, 2006).
Google Scholar 
34.van Zijderveld, S. M. et al. Nitrate and sulfate: Effective alternative hydrogen sinks for mitigation of ruminal methane production in sheep. J. Dairy Sci. 93, 5856–5866. https://doi.org/10.3168/jds.2010-3281 (2010).CAS 
Article 
PubMed 

Google Scholar 
35.Lan, W. & Yang, C. Ruminal methane production: Associated microorganisms and the potential of applying hydrogen-utilizing bacteria for mitigation. Sci. Total Environ. 654, 1270–1283. https://doi.org/10.1016/j.scitotenv.2018.11.180 (2019).ADS 
CAS 
Article 
PubMed 

Google Scholar 
36.Loubinoux, J., Bronowicki, J. P., Pereira, I. A., Mougenel, J. L. & Faou, A. E. Sulfate-reducing bacteria in human feces and their association with inflammatory bowel diseases. FEMS Microbiol. Ecol. 40, 107–112. https://doi.org/10.1111/j.1574-6941.2002.tb00942.x (2002).CAS 
Article 
PubMed 

Google Scholar 
37.Gould, D. H., Cummings, B. A. & Hamar, D. W. In vivo indicators of pathologic ruminal sulphide production in steers with diet-induced polioencephalomalacia. J. Vet. Diagn. Invest. 9, 72–76. https://doi.org/10.1177/104063879700900113 (1997).CAS 
Article 
PubMed 

Google Scholar 
38.Anderson, R. C., Rasmussen, M. A., Jensen, N. S. & Allison, M. J. Denitrobacterium detoxificans gen. nov., sp. nov., a ruminal bacterium that respires on nitrocompounds. Int. J. Syst. Evol. Microbiol. 50(Pt 2), 633–638. https://doi.org/10.1099/00207713-50-2-633 (2000).CAS 
Article 
PubMed 

Google Scholar 
39.Anderson, R. C. et al. Ruminal fermentation of anti-methanogenic nitrate- and nitro-containing forages in vitro. Front. Vet. Sci. 3, 62. https://doi.org/10.3389/fvets.2016.00062 (2016).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
40.Zhang, Z. W. et al. Nitrocompounds as potential methanogenic inhibitors in ruminant animals: A review. Anim. Feed Sci. Tech. 236, 107–114. https://doi.org/10.1016/j.anifeedsci.2017.12.010 (2018).CAS 
Article 

Google Scholar 
41.Marounek, M., Fliegrova, K. & Bartos, S. Metabolism and some characteristics of ruminal strains of Megasphaera elsdenii. Appl. Environ. Microbiol. 55, 1570–1573. https://doi.org/10.1128/AEM.55.6.1570-1573.1989 (1989).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
42.Hackmann, T. J., Ngugi, D. K., Firkins, J. L. & Tao, J. Genomes of rumen bacteria encode atypical pathways for fermenting hexoses to short-chain fatty acids. Environ. Microbiol. 19, 4670–4683. https://doi.org/10.1111/1462-2920.13929 (2017).CAS 
Article 
PubMed 

Google Scholar 
43.Janssen, P. H. Influence of hydrogen on rumen methane formation and fermentation balances through microbial growth kinetics and fermentation thermodynamics. Anim. Feed Sci. Technol. 160, 1–22. https://doi.org/10.1016/j.anifeedsci.2010.07.002 (2010).CAS 
Article 

Google Scholar 
44.Greening, C. et al. Diverse hydrogen production and consumption pathways influence methane production in ruminants. ISME J. 13, 2617–2632. https://doi.org/10.1038/s41396-019-0464-2 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
45.Gilmour, M., Flint, H. J. & Mitchell, W. J. Multiple lactate dehydrogenase activities of the rumen bacterium Selenomonas ruminantium. Microbiol. 140(Pt 8), 2077–2084. https://doi.org/10.1099/13500872-140-8-2077 (1994).CAS 
Article 

Google Scholar 
46.Chowdhury, N. P., Kahnt, J. & Buckel, W. Reduction of ferredoxin or oxygen by flavin-based electron bifurcation in Megasphaera elsdenii. FEBS J. 282, 3149–3160. https://doi.org/10.1111/febs.13308 (2015).CAS 
Article 
PubMed 

Google Scholar 
47.Weghoff, M. C., Bertsch, J. & Muller, V. A novel mode of lactate metabolism in strictly anaerobic bacteria. Environ. Microbiol. 17, 670–677. https://doi.org/10.1111/1462-2920.12493 (2015).CAS 
Article 
PubMed 

Google Scholar 
48.Hernandez, J., Benedito, J. L., Abuelo, A. & Castillo, C. Ruminal acidosis in feedlot: from aetiology to prevention. Sci. World J. 2014, 702572. https://doi.org/10.1155/2014/702572 (2014).Article 

Google Scholar 
49.Vuotto, C., Barbanti, F., Mastrantonio, P. & Donelli, G. Lactobacillus brevis CD2 inhibits Prevotella melaninogenica biofilm. Oral Dis. 20, 668–674. https://doi.org/10.1111/odi.12186 (2014).CAS 
Article 
PubMed 

Google Scholar 
50.van Lingen, H. J. et al. Thermodynamic driving force of hydrogen on rumen microbial metabolism: A theoretical investigation. PLoS One 11, e0161362. https://doi.org/10.1371/journal.pone.0161362 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
51.Ungerfeld, E. M., Aedo, M. F., Martinez, E. D. & Saldivia, M. Inhibiting methanogenesis in rumen batch cultures did not increase the recovery of metabolic hydrogen in microbial amino acids. Microorganisms 7, 155. https://doi.org/10.3390/microorganisms7050115 (2019).CAS 
Article 

Google Scholar 
52.Ng, F. et al. An adhesin from hydrogen-utilizing rumen methanogen Methanobrevibacter ruminantium M1 binds a broad range of hydrogen-producing microorganisms. Environ. Microbiol. 18, 3010–3021. https://doi.org/10.1111/1462-2920.13155 (2016).CAS 
Article 
PubMed 

Google Scholar 
53.Soliva, C. R., Amelchanka, S. L., Duval, S. M. & Kreuzer, M. Ruminal methane inhibition potential of various pure compounds in comparison with garlic oil as determined with a rumen simulation technique (Rusitec). Brit. J. Nutr. 106, 114–122. https://doi.org/10.1017/S0007114510005684 (2011).CAS 
Article 
PubMed 

Google Scholar 
54.Terranova, M. et al. Increasing the proportion of hazel leaves in the diet of dairy cows reduced methane yield and excretion of nitrogen in volatile form, but not milk yield. Anim. Feed Sci. Technol. 276, 114796. https://doi.org/10.1016/j.anifeedsci.2020.114790 (2021).CAS 
Article 

Google Scholar 
55.Ehrlich, G. G., Goerlitz, D. F., Bourell, J. H., Eisen, G. V. & Godsy, E. M. Liquid chromatographic procedure for fermentation product analysis in the identification of anaerobic bacteria. Appl. Environ. Microbiol. 42, 878–885 (1981).ADS 
CAS 
Article 

Google Scholar 
56.Bica, R. et al. Nuclear magnetic resonance to detect rumen metabolites associated with enteric methane emissions from beef cattle. Sci. Rep. 10, 5578. https://doi.org/10.1038/s41598-020-62485-y (2020).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
57.Henderson, G. et al. Effect of DNA extraction methods and sampling techniques on the apparent structure of cow and sheep rumen microbial communities. PLoS One 8, e74787. https://doi.org/10.1371/journal.pone.0074787 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
58.Kittelmann, S. et al. Simultaneous amplicon sequencing to explore co-occurrence patterns of bacterial, archaeal and eukaryotic microorganisms in rumen microbial communities. PLoS One 8, e47879. https://doi.org/10.1371/journal.pone.0047879 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
59.Milanese, A. et al. Microbial abundance, activity and population genomic profiling with mOTUs2. Nat. Comm. 10, 1014. https://doi.org/10.1038/s41467-019-08844-4 (2019).ADS 
CAS 
Article 

Google Scholar 
60.Paoli, L. et al. Uncharted biosynthetic potential of the ocean microbiome. bioRxiv https://doi.org/10.1101/2021.03.24.436479 (2021).Article 

Google Scholar 
61.Bushnell, B. BBMap: A fast, accurate, splice-aware aligner. in 9th Annual Genomics of Energy & Environment Meeting. (Lawrence Berkeley National Lab (LBNL), Berkeley, CA, USA). https://www.osti.gov/servlets/purl/1241166 (2014).62.Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. metaSPAdes: A new versatile metagenomic assembler. Genome Res. 27, 824–834. https://doi.org/10.1101/gr.213959.116 (2017).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
63.Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760. https://doi.org/10.1093/bioinformatics/btp324 (2009).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
64.Kang, D. D. et al. MetaBAT 2: An adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7, e7359. https://doi.org/10.7717/peerj.7359 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
65.Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. https://doi.org/10.1101/gr.186072.114 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
66.Seshadri, R. et al. Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection. Nat. Biotechnol. 36, 359–367. https://doi.org/10.1038/nbt.4110 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
67.Stewart, R. D. et al. Compendium of 4,941 rumen metagenome-assembled genomes for rumen microbiome biology and enzyme discovery. Nat. Biotechnol. 37, 953–961. https://doi.org/10.1038/s41587-019-0202-3 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
68.Hyatt, D. et al. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 11, 119. https://doi.org/10.1186/1471-2105-11-119 (2010).CAS 
Article 

Google Scholar 
69.Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152. https://doi.org/10.1093/bioinformatics/bts565 (2012).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
70.Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acid Res. 28, 27–30. https://doi.org/10.1093/nar/28.1.27 (2000).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
71.Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. https://doi.org/10.1038/nmeth.3176 (2015).CAS 
Article 
PubMed 

Google Scholar 
72.Li, L., Stoeckert, C. J. Jr. & Roos, D. S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189. https://doi.org/10.1101/gr.1224503 (2003).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
73.Allen, B., Drake, M., Harris, N. & Sullivan, T. Using KBase to assemble and annotate prokaryotic genomes. Curr. Protoc. Microbiol. 46, 1E 13 11-11E 13 18. https://doi.org/10.1002/cpmc.37 (2017).Article 

Google Scholar 
74.RStudio Team. R Studio: Integrated development environment for R. Version 1.4.1106 (2021).75.Oksanen, J. et al. The vegan package. Community Ecol. Pack. 10, 631–637 (2007).
Google Scholar 
76.Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. https://doi.org/10.1186/s13059-014-0550-8 (2014).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar  More

Insects and bacteriaBean bugs were reared in petri dishes (90 mm in diameter and 20-mm high) at 25 °C under a long-day regimen (16-h light, 8-h dark) and fed with soybean seeds and distilled water containing 0.05% ascorbic acid (DWA). Burkholderia symbiont strain SFA119, a MEP-degrading strain conferring MEP resistant in the bean bug, and its GFP-(green fluorescent protein) labeled derivative, strain SJ586, were used in this study. The symbiont was cultured at 30 °C on YG medium (0.5% yeast extract, 0.4% glucose, and 0.1% NaCl). The GFP-labeled strain was constructed by the Tn7 mini-transposon system, as previously described31.Genome sequencingDNA was extracted from cultured cells of strain SFA1 by the phenol–chloroform extraction as previously described32. The DNA library for Illumina short reads (the mean insert size: 500 bp) was constructed by using the Covaris S2 ultrasonicator (Covaris) and the KAPA HyperPrep Kit (Kapa Biosystems). For the library construction for Nanopore long reads, Native Barcoding Expansion (EXP-NBD104, Oxford Nanopore Technologies) and the Ligation Sequencing Kit (SQK-LSK109, Oxford Nanopore Technologies) were used. The genome sequencing was performed with NextSeq using the 2 × 151-bp protocol (Illumina) and GridION using an R9.4.1 flow cell (Oxford Nanopore Technologies). The Illumina short reads were processed by using Sickle Ver 1.33 (available at https://github.com/najoshi/sickle) for removing the low-quality and shorter reads. After processing the Nanopore long-reads with Porechop Ver 0.2.3 (available at https://github.com/rrwick/Porechop) and Filtlong Ver 0.2.0 (available at https://github.com/rrwick/Filtlong), error correction was performed by using Canu Ver 1.833. These processed short- and long reads were assembled by using Unicycler Ver 0.4.734, resulting in the eight circular replicons (Supplementary Fig. 1). The assembled genome was annotated by DFAST Ver 1.1.035. After the homology searches of the protein sequences by blastp 2.5.0 + 36 against the COG database (PMID: 25428365), circular replicons were visualized with circos v 0.69-837. The chromosomes and plasmids were assigned according to the genome of Caballeronia (Burkholderia) cordobensis strain YI2338.Phylogenetic analysisNucleotide sequences of 16 S rRNA gene of representative Burkholderia spp. and outgroup species were aligned by using SINA v1.2.1139. Protein sequences of MEP-degrading genes (mpd, pnpB, and mhqA) and a plasmid-transfer gene (traH) on plasmid 2 were subjected to the blastp search against the nr database (downloaded in Jul. 2019) and top ~30 hit sequences were retrieved for each gene. Multiple sequencing alignments of each gene were constructed with L-INS-I of mafft v7.40740. Gap-including and ambiguous sites in the alignments were then removed. Unrooted maximum-likelihood (ML) phylogenetic trees were reconstructed with RAxML v8.2.341 using the GTR + Γ model (for 16 S rRNA gene) or the LG + Γ model42 (for other genes). The bootstrap values of 1000 replicates for all internal branches were calculated with a rapid bootstrapping algorithm43.Preparation of SFA1 cultures for RNA-seqBurkholderia symbiont SFA1 was precultured in minimal medium (20 mM phosphate buffer [pH 7.0], 0.01% yeast, 0.1% (NH4)2SO4, 0.02% NaCl, 0.01% MgSO4⋅7H2O, 0.005% CaCl2⋅2H2O, 0.00025% FeSO4⋅7H2O, and 0.00033% EDTA⋅2Na) containing 1.0 mM of MEP on a gyratory shaker (210 rpm) at 30 °C overnight, and subcultured in newly prepared MEP-containing minimal medium under the same conditions for 5 h. As a control, SFA1 was precultured in minimal medium containing 0.1% citrate overnight, and then the overnighter was subcultured in a newly prepared citrate-containing minimal medium under the same conditions for 10 h. The culture was mixed with an equal amount of RNAprotect Bacteria Regent (Qiagen, Valencia, CA, USA), then centrifuged to harvest the cells for the RNA-seq analysis.Preparation of midgut symbiont cells for RNA-seqThe oral administration of the symbiont strain SFA1 was performed as described19,44. The symbiont was inoculated to 2nd instar nymphs, and three days after molting to the 3rd instar, nymphs were transdermally administered with 1 µl of 0.2 µM or 20 µM of MEP (dissolved in acetone). One- or three days after the treatment, insects were dissected and the crypt-bearing symbiotic gut region was subjected to the RNA extraction and RNA-seq analysis. As a control, untreated insects were analyzed.RNA-seq analysisTotal RNA was extracted from triplicate samples from cultures by the hot-phenol method as previously described45 or from the midgut symbiont cells by using RNAiso Plus (Takara Bi, Kusatsu, Shiga, Japan) and the RNeasy mini kit (Qiagen). The extracted total RNA was purified by phenol–chloroform extraction and digestion by DNase (RQ1 RNase-Free DNase, Promega, Fitchburg, WI, USA) and repurified by using a RNeasy Mini Kit. The mRNA in the samples was further enriched by the RiboMinus Transcriptome Isolation Kit bacteria (Thermo Fisher Scientific, Waltham, MA, USA) and the RiboMinus Eukaryote Kit for RNA-Seq (Thermo Fisher Scientific), and purified by using an AMPure XP kit (Beckman Coulter, Brea, CA, USA). The cDNA libraries were constructed from approximately 100 ng of rRNA-depleted RNA samples by the use of a NextUltraRNA library prep kit (New England Biolabs, Ipswich, MA, USA). Size selection of cDNA (200–300 bp) and determination of the size distribution and concentration of the purified cDNA samples were performed as described previously46. In total, 21 cDNA libraries were constructed and sequenced by MiSeq (Illumina, Inc., San Diego, CA, USA). To ensure high sequence quality, the remaining sequencing adapters and the reads with a cutoff Phred score of 15 (for leading and tailing sequences, Phred score of >20) and a length of less than 80 bp in the obtained RNA-seq data were removed by the program Trimmomatic v0.30 using Illumina TruSeq3 adapter sequences for the clipping47. The remaining paired reads were analyzed by FastQC version 0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) for quality control, and Bowtie2 ver. 2.2.248 for mapping on the symbiont genome (DDBJ/EMBL/GenBank accession: AP022305–AP022312). After the conversion of the output BAM files to BED files using the bamtobed program in BEDTools ver. 2.14.349, gene expression levels were calculated in TPM (transcripts per kilobase million) values by using in-house scripts46.Gene deletion and complementationMEP-degrading genes (mpd, pnpA1, and pnpA2) were deleted by the homologous-recombination-based deletion method using pK18mobsacB or pUC18, as previously described50,51. Primers used for the mutagenesis are listed in Supplementary Table 1. For mpd gene deletion, pK18mobsacB was used to construct a markerless mutant. For single deletion of pnpA1 and pnpA2 genes, pUC18 was used to substitute each gene locus with a kanamycin-resistance gene cassette. The double deletion of pnpA1 and pnpA2 genes was performed by substituting pnpA2 gene locus with a tetracycline-resistance gene cassette in the pnpA1-deletion mutant. Gene complementation of mpd was also performed by homologous recombination using plasmid pUC18 with primers listed in Supplementary Table 1. To investigate growth profiles of the wild-type SFA1, the gene-deletion mutants (Δmpd, ΔpnpA1, ΔpnpA2, and ΔpnpA1/ΔpnpA2), and the mph-complement mutant (Δmpd/mpd+) in the MEP-containing minimal medium, the strains were precultured in minimal medium containing 1.0 mM MEP on a gyratory shaker (210 rpm) at 30 °C overnight, and then cultured in newly prepared MEP-containing minimal medium under the same condition. The growth of cultures was estimated by OD600 measurements. To confirm the basic growth abilities of the mutants, these bacterial strains were pre- and subcultured in minimal medium containing 0.1% glucose under the same conditions. These symbiont strains and mutants were inoculated to the bean bug as described above.Quantitative PCRSymbiont titers in the midgut crypts were evaluated by quantitative PCR (qPCR) of bacterial dnaA gene copies. The qPCR was performed by using a KAPA SYBR Fast qPCR Master Mix (Kapa Biosystems) and the LightCycler 96 System (Roche Applied Science) with the following primers: BSdnaA–F (5′-AGC GCG AGA TCA GAC GGT CGT CGA T-3′) and BSdnaA–R (5′-TCC GGC AAG TCG CGC ACG CA-3′).MEP treatment of insectsMEP treatment of R. pedestris was performed as previously described19. Soybean seeds were dipped in 0.2 mM MEP for 5 s and dried at room temperature. In each clean plastic container, 15 individuals of 3rd-instar nymphs were reared on three seeds of the MEP-treated soybean and DWA at 25 °C under the long-day regime, and the number of dead insects was counted 24 h after the treatments. The survival rate of the insects was analyzed under Fisher’s exact test by use of the program R ver. 3.6.3 (available at https://www.R-project.org/). Multiple comparisons were corrected by the Bonferroni method.Bactericidal activities of MEP and its degradation product 3M4NTo measure bactericidal activities of MEP and 3M4N on cultured cells of SFA1, 104 cells of log-phase growing bacteria were mixed with a defined concentration of MEP or 3M4N, and spotted on a YG agar plate. To measure the bactericidal activity against midgut crypt-colonizing cells, the symbiotic organs infected with SFA1 were dissected from 3rd-instar insects, homogenized in PBS, and purified by a 5-µm-size pore Syringe filter to harvest colonizing symbiont cells50. MEP or 3M4N was added to approximately 104 cells of the harvested cells and spotted on a YG agar plate. Bactericidal activities of the chemical compounds were then checked in 24 h after incubation at 30 °C.HPLC detection of in vitro and in vivo MEP-degrading activities of the symbiontTo determine in vitro MEP-degradation activity, cultured cells of SFA1 were prepared as above, and 106 cells were incubated at 25 °C in 200 µl of MEP solution (2 mM MEP in Tris-Hcl [pH 8.5] with 0.1% Triton X-100) in a 1.5-ml microtube. To determine in vivo MEP-degradation activity, the midgut of a 5th-instar insect infected with SFA1 was dissected, the posterior and anterior parts of the crypt-bearing symbiotic region were closed with 0.2-mm polyethylene fishline (Supplementary Fig. 6a), and incubated at 25 °C in 200 µl of the MEP solution. For the in vivo determination, 250 mM of trehalose, known as a major sugar source of insects’ hemolymph52, was added to the MEP solution to keep the tissue fresh. After incubation for different times, the reaction was stopped by adding 400 µl of methanol. After centrifugation, supernatants were subjected to high-performance liquid chromatography (HPLC) analyses to detect MEP and 3M4N, as previously reported21, and precipitated cells and tissues were subjected to DNA extraction and qPCR to estimate symbiont-cell numbers of each reaction.LC–ESI–MS detection of 3M4N in feces from 3M4N-fed insectsAn insect-rearing system for feeding 3M4N and collecting feces is shown in Supplementary Fig. 7. Insects were fed with DW or DW containing 10 mM 3M4N in a plastic container, in which the solution supplier was covered by 0.5-mm mesh, so that insects were able to drink the solution by probing with their proboscis, but did not directly touch the solution by their legs or body. Twenty insects were reared per container and their feces were accumulated on the bottom of the container for five days. The collected feces (DW- or 3M4N-treated) were suspended in 1 ml of MilliQ water, and the water-soluble fractions were extracted by thorough vortexing. Solids and insoluble fractions were removed from the suspension by centrifugation and subsequent filtration using a cellulose-acetate membrane (Φ, 0.20 μm, ADVANTEC, Tokyo, Japan). The resultant fraction was diluted 10-fold by MilliQ water and analyzed by liquid chromatography–electrospray-ionization mass spectrometry (LC–ESI–MS) according to a previous report53,54,55. HPLC was performed using the Nexera X2 system (Shimadzu, Kyoto, Japan) composed of LC-30AD pump, SPD-M30A photodiode-array detector, and SIL-30AC autosampler. Develosil HB ODS-UG column (ID 2.0 mm × L 75 mm, Nomura Chemical Co., Ltd, Aichi, Japan) was employed with a flow rate of 0.2 mL/min. The following gradient system was used for analysis of metabolites: MilliQ water (solvent A) and methanol (solvent B), 90% A and 10% B at 0–5 min, linear gradient from 90% A and 10% B to 20% A and 80% B at 5–15 min, 20% A and 80% B at 15–20 min, and 90% A and 10% B at 20–25 min. Retention time of 3M4N standard reagent was 14.2 min. Electrospray-ionization mass spectrometry (ESI–MS) in positive and negative ion modes was simultaneously performed using amaZon SL (Bruker, Billerica, MA, USA). 3M4N (MW = 153.14) standard showed a clear peak in negative mode at m/z of 151.53.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Nature Video
05 November 2021

Whale-cams reveal how much they really eat

Baleen whales consume twice as much krill as previously estimated.

Sara Reardon

0

Sara Reardon

Sara Reardon is a freelance writer in Bozeman, Montana.

View author publications

You can also search for this author in PubMed
 Google Scholar

Share on Twitter
Share on Twitter

Share on Facebook
Share on Facebook

Share via E-Mail
Share via E-Mail

Tagging whales with cameras and sensors has allowed researchers to calculate how much food these huge creatures are consuming. It’s the most accurate estimate yet and reveals an even more significant impact of whales on ocean ecosystems than was previously known.Read the paper here.

doi: https://doi.org/10.1038/d41586-021-03026-z

Related Articles

Vikings were living in North America exactly a thousand years ago

Strange patterns could be the world’s oldest animal fossil

Divers are replanting reefs one coral at a time

How cuttlefish wear their thoughts on their skin

Swimming shrimp: Causing a stir

Subjects

Zoology

Ecology

Environmental sciences

Latest on:

Zoology

Sponge cells hint at origins of nervous system
News 05 NOV 21

Baby bats try out their ‘sonar’ just after birth
Research Highlight 27 OCT 21

Ivory hunting drives evolution of tuskless elephants
News 21 OCT 21

Ecology

Whales’ gigantic appetites, climate fears — the week in infographics
News 05 NOV 21

COP26 climate pledges: What scientists think so far
News 05 NOV 21

Baleen whale prey consumption based on high-resolution foraging measurements
Article 03 NOV 21

Environmental sciences

Carbon implications of marginal oils from market-derived demand shocks
Article 03 NOV 21

For NGOs, article-processing charges sap conservation funds
Correspondence 02 NOV 21

Embrace open-source sensors for local climate studies
Correspondence 02 NOV 21

Jobs

Postdoctoral Training Fellow

Francis Crick Institute
London, United Kingdom

Head of GeMS

Francis Crick Institute
London, United Kingdom

Postdoctoral scientist (m/f/div) to work on the comparative genomics of gutless marine oligochaetes

Max Planck Institute for Marine Microbiology
Bremen, Germany

Postdoctoral Positions (m/f/div) in Protist Virology

Max Planck Institute for Medical Research (MPIMF)
Heidelberg, Germany More

NEWS
05 November 2021

COP26 climate pledges: What scientists think so far

Nations have promised to end deforestation, curb methane emissions and stop public investment in coal power. Researchers warn that the real work of COP26 is yet to come.

Ehsan Masood

&

Jeff Tollefson

Ehsan Masood

View author publications

You can also search for this author in PubMed
 Google Scholar

Jeff Tollefson

View author publications

You can also search for this author in PubMed
 Google Scholar

Share on Twitter
Share on Twitter

Share on Facebook
Share on Facebook

Share via E-Mail
Share via E-Mail

Methane burns at an oil pit. Among the key pledges so far at COP26 is an agreement to cut methane emissions by 30% by 2030.Credit: Orjan F. Ellingvag/Corbis via Getty

The first few days of the 26th United Nations Climate Change Conference of the Parties (COP26) have seen a flurry of announcements from world leaders promising to tackle climate change — from plans to phase out public finance for coal-fired power, to a pledge to end deforestation. This year, many big names — including US President Joe Biden and Indian Prime Minister Narendra Modi — attended the first two days of the conference to make big announcements.
COP26 climate summit: A scientists’ guide to a momentous meeting
This is different from what has happened at most previous COP summits, says Beth Martin, a specialist in climate negotiation who is part of RINGO (Research and Independent Non-Governmental Organizations), a network of organizations allowed to observe the COP26 negotiations. Usually, the highest-profile figures aren’t present during the first week, but arrive near the end of the meeting to help bridge differences in time for an agreed statement, and for the obligatory ‘UN family photo’.Nature asked researchers what they think of the pledges that have been made so far, as negotiators from some 200 countries prepare to dive into more detailed talks.Methane emissionsOne of the key developments in the first week was an agreement to curb emissions of methane, a powerful greenhouse gas that is second only to carbon dioxide in terms of its impact on the climate. Led by the United States and the European Union, the global methane pledge seeks to curb methane emissions by 30% by 2030, and has been signed by more than 100 countries.
Control methane to slow global warming — fast
“Obviously, as a scientist you’d say, ‘Well, a 50% reduction in the methane emissions by 2030 would be even better,’ but it’s a good start,” says Tim Lenton, who heads the Global Systems Institute at the University of Exeter, UK. “It’s an additional lever that could really help us limit warming.”Research has shown1 that curbing methane emissions using existing technologies could shave up to 0.5 °C off global temperatures by 2100. As with carbon dioxide, however, limiting methane emissions will not happen on its own.With his climate agenda facing challenges in Congress, Biden made methane a centrepiece of his commitments in Glasgow by announcing a new regulation to curb methane emissions from the oil and gas industry. Put forward this week by the US Environmental Protection Agency, the rule would require companies to curb methane emissions from their facilities by 74% over the coming decade, compared with 2005 levels. If implemented as proposed, it could prevent the release of some 37 million tonnes of methane by 2035 — equivalent to more than the annual carbon emissions from the nation’s fleets of passenger vehicles and commercial aircraft.India’s net-zero goalAfter delaying expected updates to India’s climate commitments by more than a year, Modi captured the world’s attention early in the summit by announcing that his country would seek to achieve net-zero emissions by 2070. The deadline is decades after that of many other countries that have made net-zero commitments, and it remains unclear whether India is committing to curbing just carbon dioxide emissions, or the broader category of greenhouse-gas emissions. But scientists say the announcement could mark a significant step forward if India follows through.
Scientists cheer India’s ambitious carbon-zero climate pledge
“We are definitely taken by surprise: this is much more than we were expecting to hear,” says Ulka Kelkar, an economist in Bengaluru who heads the Indian climate programme for the World Resources Institute, an environmental think tank based in Washington DC.Many scientists remain sceptical about mid-century net-zero pledges, in part because it’s easy to make long-term promises but hard to make the difficult short-term decisions that are required to meet those pledges. But India’s commitment includes measurable near-term targets, such as a pledge to provide 50% of the nation’s power through renewable resources and to reduce projected carbon emissions by one billion tonnes of carbon dioxide by 2030.Questions remain about how these targets will be defined and measured, but models indicate that there is a 50% chance such net-zero pledges could limit global warming to 2 °C or less, if fully implemented by all countries.

More than 130 countries have agreed to halt and reverse deforestation by 2030.Credit: Joao Laet/AFP via Getty

Climate cashAmong a cascade of climate-finance announcements this week is a pledge from more than 450 organizations in the financial sector — including banks, fund managers and insurance companies — in 45 countries to move US$130 trillion of funds under their control into investments where the recipient is committed to net-zero emissions by 2050.The pledging institutions, which are part of the Glasgow Financial Alliance for Net Zero, have not yet specified interim targets or timetables to achieve this goal. On 1 November, UN secretary-general António Guterres announced that a group of independent experts would be convened to propose standards for such commitments to net-zero emissions.
The broken $100-billion promise of climate finance – and how to fix it
Governments also announced new investments in clean technologies. And more than 40 countries, including the United Kingdom, Poland, South Korea and Vietnam, have committed to phasing out coal power in the 2030s (for major economies) or 2040s (globally), and to stopping public funding for new coal-fired power plants.“All of this is significant,” says Cristián Samper, an ecologist and president of the Wildlife Conservation Society in New York City. “The involvement of the financial sector and of ministers of finance and energy” in the meeting “is a game-changer”.However, the announcements have been overshadowed by governments’ failure to meet a 2009 pledge to provide $100 billion annually in climate finance for low- and middle-income countries by 2020. Reports suggest that it will take another two years to reach this goal, and that around 70% of the finance will be provided as loans.“We all assumed it would be grant finance. We didn’t pay attention to the fine print or expect that developed countries would hide behind loans,” says climate economist Tariq Banuri, a former director of sustainable development at the UN.Ending deforestationMore than 130 countries have pledged to halt and reverse forest-loss and land degradation by 2030. The signatories, which include Brazil, the Democratic Republic of the Congo and Indonesia, are home to 90% of the world’s forests.It is not the first such commitment: the 2014 New York Declaration on Forests, signed by a broad coalition of nearly 200 countries, regional governments, companies, indigenous groups and others, called for halving deforestation by 2020 and “striving” to end it by 2030.
The United Nations must get its new biodiversity targets right
There is also a long-standing UN pledge to slow down and eventually reverse the loss of biodiversity. But this remains unfulfilled and there is no official monitoring. Researchers say the latest target is unlikely to be met without an enforcement mechanism.Separately, a group of high-income countries has pledged $12 billion in public finance for forest protection between 2021 and 2025, but has not specified how the funding will be provided. A statement from the group, which includes Canada, the United States, the United Kingdom and EU countries, says governments will “work closely with the private sector” to “leverage vital funding from private sources to deliver change at scale”. This suggests that the finance is likely to be dominated by loans. Still, Samper says that there are reasons to be optimistic. Few previous climate COPs discussed nature and forests on the scale now seen in Glasgow. In the past, if biodiversity was mentioned at a climate meeting, “it was like the Martians had landed”, he says, because biodiversity and climate are treated as separate challenges by the UN. “We’ve never seen this much attention. It could be a pivot point.”

doi: https://doi.org/10.1038/d41586-021-03034-z

References1.Ocko, I. B. et al. Environ. Res. Lett. 16, 054042 (2021).Article 

Google Scholar 
Download references

Related Articles

COP26 climate summit: A scientists’ guide to a momentous meeting

Scientists cheer India’s ambitious carbon-zero climate pledge

The broken $100-billion promise of climate finance – and how to fix it

The United Nations must get its new biodiversity targets right

Control methane to slow global warming — fast

Subjects

Funding

Biodiversity

Politics

Climate change

Latest on:

Funding

The African Academy of Sciences is in crisis — responsibility must be shared
Editorial 03 NOV 21

UK research funding to grow slower than hoped
News 28 OCT 21

The high burden of infectious disease
Nature Index 27 OCT 21

Biodiversity

The answer to the biodiversity crisis is not more debt
Editorial 26 OCT 21

Illegal mining in the Amazon hits record high amid Indigenous protests
News 30 SEP 21

Fine-root traits in the global spectrum of plant form and function
Article 29 SEP 21

Politics

All aboard the climate train! Scientists join activists for COP26 trip
News 02 NOV 21

Top climate scientists are sceptical that nations will rein in global warming
News Feature 01 NOV 21

UK research funding to grow slower than hoped
News 28 OCT 21

Jobs

Postdoctoral Training Fellow

Francis Crick Institute
London, United Kingdom

Head of GeMS

Francis Crick Institute
London, United Kingdom

Postdoctoral scientist (m/f/div) to work on the comparative genomics of gutless marine oligochaetes

Max Planck Institute for Marine Microbiology
Bremen, Germany

Postdoctoral Positions (m/f/div) in Protist Virology

Max Planck Institute for Medical Research (MPIMF)
Heidelberg, Germany More

1.Frankham, R., Ballou, J. D. & Briscoe, D. A. Introduction to Conservation Genetics 617 (Cambridge University Press, 2002).Book 

Google Scholar 
2.Spielman, D., Brook, B. W. & Frankham, R. Most species are not driven to extinction before genetic factors impact them. PNAS 101, 15261–15264 (2004).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
3.Ahmadi, M. et al. Evolutionary applications of phylo-genetically-informed ecological niche modelling (ENM) to explore cryptic diversification over cryptic refugia. Mol. Phylogenet. Evol. 127, 712–722 (2018).PubMed 
Article 

Google Scholar 
4.Ashrafzadeh, M. R. et al. Large-scale mitochondrial DNA analysis reveals new light on the phylogeography of Central and Eastern-European Brown hare (Lepus europaeus Pallas, 1778). PLoS ONE 13, e0204653 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
5.Tóth, B. et al. Genetic diversity and structure of common carp (Cyprinus carpio L.) in the Centre of Carpathian Basin: Implications for conservation. Genes 11, 1268 (2020).PubMed Central 
Article 
CAS 
PubMed 

Google Scholar 
6.Broquet, T., Ray, N., Petit, E., Fryxell, J. M. & Burel, F. Genetic isolation by distance and landscape connectivity in the American marten (Martes americana). Landsc. Ecol. 21, 877–889 (2006).Article 

Google Scholar 
7.Cushman, S. A., McKelvey, K. S., Hayden, J. & Schwartz, M. K. Gene flow in complex landscapes: testing multiple hypotheses with causal modeling. Am. Nat. 168, 486–499 (2006).PubMed 
Article 

Google Scholar 
8.Khosravi, R. et al. Effect of landscape features on genetic structure of the goitered gazelle (Gazella subgutturosa) in Central Iran. Conserv. Genet. 19, 323–336 (2018).Article 

Google Scholar 
9.Adavodi, R., Khosravi, R., Cushman, S. A. & Kaboli, M. Topographical features and forest cover influence landscape connectivity and gene flow of the Caucasian pit viper, Gloydius caucasicus (Nikolsky, 1916), Iran. Landsc. Ecol. 34, 2615–2630 (2019).Article 

Google Scholar 
10.Moussy, C. et al. Migration and dispersal patterns of bats and their influence on genetic structure. Mammal Rev. 43, 183–195 (2013).Article 

Google Scholar 
11.Theodoridis, S. et al. Evolutionary history and past climate change shape the distribution of genetic diversity in terrestrial mammals. Nat. Commun. 11, 1–11 (2020).MathSciNet 
Article 
CAS 

Google Scholar 
12.Barilani, M. et al. Detecting introgressive hybridisation in rock partridge populations (Alectoris graeca) in Greece through Bayesian admixture analyses of multilocus genotypes. Conserv. Genet. 8, 343–354 (2007).CAS 
Article 

Google Scholar 
13.Randi, E. Detecting hybridization between wild species and their domesticated relatives. Mol. Ecol. 17, 285–293 (2008).PubMed 
Article 

Google Scholar 
14.Kusza, S., Ashrafzadeh, M. R., Tóth, B. & Jávor, A. Maternal genetic variation in the northeastern Hungarian fallow deer (Dama dama) population. Mamm. Biol. 93, 21–28 (2018).Article 

Google Scholar 
15.Laikre, L., Schwartz, M. K., Waples, R. S., Ryman, N. & GeM Working Group. Compromising genetic diversity in the wild: Unmonitored large-scale release of plants and animals. Trends Ecol. Evol. 25, 520–529 (2010).PubMed 
Article 

Google Scholar 
16.Söderquist, P. et al. Admixture between released and wild game birds: a changing genetic landscape in European mallards (Anas platyrhynchos). Eur. J. Wildl. Res. 63, 98 (2017).Article 

Google Scholar 
17.Robertson, P. A. et al. Pheasant release in Great Britain: Long-term and large-scale changes in the survival of a managed bird. Eur. J. Wildlife Res. 63, 100 (2017).Article 

Google Scholar 
18.Mank, J. E., Carlson, J. E. & Brittingham, M. C. A century of hybridization: Decreasing genetic distance between American black ducks and mallards. Conserv. Genet. 5, 395–403 (2004).CAS 
Article 

Google Scholar 
19.Blanco-Aguiar, J. A. et al. Assessment of game restocking contributions to anthropogenic hybridization: The case of the Iberian red-legged partridge. Anim. Conserv. 11, 535–545 (2008).Article 

Google Scholar 
20.Sanchez-Donoso, I. et al. Are farm-reared quails for game restocking really common quails (Coturnix coturnix)? a genetic approach. PLoS ONE 7, e39031 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
21.Liu, Y. et al. Genome assembly of the common pheasant Phasianus colchicus: A model for speciation and ecological genomics. Genome Biol. Evol. 11, 3326–3331 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
22.Braasch, T., Pes, T., Michel, S. & Jacken, H. The subspecies of the common pheasant Phasianus colchicus in the wild and captivity. Int. J. Galliformes Conserv. 2, 6–13 (2011).
Google Scholar 
23.Robertson, D. H. P. & Hill, D. A. The Pheasant: Ecology, Management and Conservation 281 (Blackwell Scientific Publication, 1988).
Google Scholar 
24.Hill, D. A. & Robertson, P. Hand reared pheasants: how do they compare with wild birds. Game Conserv. Ann. Rep. 17, 76–84 (1986).
Google Scholar 
25.Pfarr, J. True Pheasants: A Noble Quarry 248 (Hancock House Publishers Ltd, 2012).
Google Scholar 
26.BirdLife International. Phasianus colchicus. The IUCN Red List of Threatened Species 2016: e.T45100023A85926819. https://doi.org/10.2305/IUCN.UK.2016-3.RLTS.T45100023A85926819.en27.Ashoori, A. et al. Habitat modeling of the common pheasant Phasianus colchicus (Galliformes: Phasianidae) in a highly modified landscape: application of species distribution models in the study of a poorly documented bird in Iran. Eur. Zool. J. 85, 372–380 (2018).Article 

Google Scholar 
28.Del Hoyo, J., Elliott, A. & Sargatal, J. Handbook of the Birds of the World: New World Vultures to Guineafowl Vol. 2, 638 (Lynx Edicions, 1994).
Google Scholar 
29.Boev, Z. N. Late Pleistocene and Holocene avifauna from three caves in the vicinity of Tran (Pernik District-W Bulgaria). In Proceedings of the First National Conference on Environment and Cultural Heritage in Karst, Sofia, Bulgaria, 10–12 November 2000 (eds Delchev, P. et al.) (Earth and Man National Museum, Association of Environment and Cultural Heritage in Karst, 2001).
Google Scholar 
30.Qu, H. et al. Subspecies boundaries and recent evolution history of the common pheasant (Phasianus colchicus) across China. Biochem. Syst. Ecol. 71, 155–162 (2017).CAS 
Article 

Google Scholar 
31.Kayvanfar, N., Aliabadian, M., Niu, X., Zhang, Z. & Liu, Y. Phylogeography of the common pheasant Phasianus colchicus. Ibis 159, 430–442 (2017).Article 

Google Scholar 
32.Todesco, M. et al. Hybridization and extinction. Evol. Appl. 9, 892–908 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
33.Santilli, F. & Bagliacca, M. Factors influencing pheasant Phasianus colchicus harvesting in Tuscany, Italy. Wildl. Biol. 14, 281–287 (2008).Article 

Google Scholar 
34.Lavadinović, V., Beuković, D. & Popović, Z. Common Pheasant (Phasianus colchicus L.1758) management in Serbia. Contemp. Agric. 68, 71–79 (2019).Article 

Google Scholar 
35.Fenberg, P. B. & Roy, K. Ecological and evolutionary consequences of size-selective harvesting: How much do we know?. Mol. Ecol. 17, 209–220 (2008).PubMed 
Article 

Google Scholar 
36.Allendorf, F. W. & Hard, J. J. Human-induced evolution caused by unnatural selection through harvest of wild animals. PNAS 106, 9987–9994 (2009).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
37.Darimont, C. T. et al. Human predators outpace other agents of trait change in the wild. PNAS 106, 952–954 (2009).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
38.Madden, J. R. & Whiteside, M. A. Selection on behavioural traits during ‘unselective’harvesting means that shy pheasants better survive a hunting season. Anim. Behav. 87, 129–135 (2014).Article 

Google Scholar 
39.Csányi, S. The effect of hand-reared pheasants on the wild population in Hungary: A modelling approach. Hung. Small Game B. 5, 71–82 (2000).
Google Scholar 
40.Queirós, J., Gortázar, C. & Alves, P. C. Deciphering anthropogenic effects on the genetic background of the Red deer in the Iberian Peninsula. Front. Ecol. Evol. 8, 147 (2020).Article 

Google Scholar 
41.Giesel, J. T., Brazeau, D., Koppelman, R. & Shiver, D. Ring-necked pheasant population genetic structure. J. Wildl. Manag. 61, 1332–1338 (1997).Article 

Google Scholar 
42.Qu, J., Liu, N., Bao, X. & Wang, X. Phylogeography of the ring-necked pheasant (Phasianus colchicus) in China. Mol. Phylogenet. Evol. 52, 125–132 (2009).CAS 
PubMed 
Article 

Google Scholar 
43.Liu, Y., Zhan, X., Wang, N., Chang, J. & Zhang, Z. Effect of geological vicariance on mitochondrial DNA differentiation in Common Pheasant populations of the Loess Plateau and eastern China. Mol. Phylogenet. Evol. 55, 409–417 (2010).PubMed 
Article 

Google Scholar 
44.Zhang, L., An, B., Backström, N. & Liu, N. Phylogeography-based delimitation of subspecies boundaries in the Common Pheasant (Phasianus colchicus). Biochem. Genet. 52, 38–51 (2014).PubMed 
Article 
CAS 

Google Scholar 
45.Liu, S. et al. Regional drivers of diversification in the late Quaternary in a widely distributed generalist species, the common pheasant Phasianus colchicus. J. Biogeogr. 47, 2714–2727 (2020).Article 

Google Scholar 
46.Avise, J. C. Phylogeography: The History and Formation of Species 464 (Harvard University Press, 2000).Book 

Google Scholar 
47.Corin, S. E., Lester, P. J., Abbott, K. L. & Ritchie, P. A. Inferring historical introduction pathways with mitochondrial DNA: the case of introduced Argentine ants (Linepithema humile) into New Zealand. Divers. Distrib. 13, 510–518 (2007).Article 

Google Scholar 
48.Oskarsson, M. C. et al. Mitochondrial DNA data indicate an introduction through Mainland Southeast Asia for Australian dingoes and Polynesian domestic dogs. Proc. R. Soc. B. 279, 967–974 (2012).PubMed 
Article 

Google Scholar 
49.Garrett, L. J. et al. Spatio-temporal processes drive fine-scale genetic structure in an otherwise panmictic seabird population. Sci. Rep. 10, 1–12 (2020).Article 
CAS 

Google Scholar 
50.Brito, P. H. Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco). Mol. Ecol. 16, 3423–3437 (2007).CAS 
PubMed 
Article 

Google Scholar 
51.Suárez, N. M. et al. Phylogeography and genetic struc-ture of the Canarian common chaffinch (Fringilla coelebs) inferred with mtDNA and microsatellite loci. Mol. Phylogenet. Evol. 53, 556–564 (2009).PubMed 
Article 
CAS 

Google Scholar 
52.Piry, S., Luikart, G. & Cornuet, J. M. Computer note. BOTTLENECK: A computer program for detecting recent reductions in the effective size using allele frequency data. J. Hered. 90, 502–503 (1999).Article 

Google Scholar 
53.Rogers, A. R. & Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569 (1992).CAS 
PubMed 

Google Scholar 
54.Slatkin, M. & Hudson, R. R. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 129, 555–562 (1991).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
55.Wang, B. et al. Development and characterization of novel microsatellite markers for the Common Pheasant (Phasianus colchicus) using RAD-seq. Avian Res. 8, 4 (2017).Article 

Google Scholar 
56.Grant, W. A. S. & Bowen, B. W. Shallow population histories in deep evolutionary lineages of marine fishes: Insights from sardines and anchovies and lessons for conservation. J. Hered. 89, 415–426 (1998).Article 

Google Scholar 
57.Barrowclough, G. F., Johnson, N. K. & Zink, R. M. In Current Ornithology Vol. 2 (ed. Johnston, R. F.) 135–154 (Springer, 1985).Chapter 

Google Scholar 
58.Zink, R. M. Phylogeographic studies of North American birds. In Avian Molecular Evolution and Systematics (ed. Mindell, D. P.) 301–324 (Academic Press, 1997).Chapter 

Google Scholar 
59.Payne, R. B. Natal dispersal and population structure in a migratory songbird, the Indigo Bunting. Evolution 45, 49–62 (1991).PubMed 
Article 

Google Scholar 
60.Stenzel, L. E. et al. Long-distance breeding dispersal of snowy plovers in western North America. J. Anim. Ecol. 63, 887–902 (1994).Article 

Google Scholar 
61.Zhu, C. et al. Genetic structure and population dynamics of the silver pheasant (Lophura nycthemera) in southern China. Turk. J. Zool. 44, 31–43 (2020).Article 

Google Scholar 
62.Faragó, S. Élőhelyfejlesztés az Apróvad-Gazdálkodásban: A Fenntartható Apróvad-Gazdálkodás Környezeti Alapjai 341 (Mezőgazda Kiadó, 1997).
Google Scholar 
63.Faragó, S. & Náhlik, A. A Vadállomány Szabályozása: A Fenntartható Vadgazdálkodás Populációökológiai Alapjai 347 (Mezőgazda kiadó, 2011).
Google Scholar 
64.Randi, E. & Lucchini, V. Organization and evolution of the mitochondrial DNA control region in the avian genus Alector-is. J. Mol. Evol. 47, 449–462 (1998).ADS 
CAS 
PubMed 
Article 

Google Scholar 
65.Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
66.Xia, X. DAMBE6: New tools for microbial genomics, phylogenetics, and molecular evolution. J. Hered. 108, 431–437 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
67.Librado, P. & Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).CAS 
PubMed 
Article 

Google Scholar 
68.Posada, D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 25, 1253–1256 (2008).CAS 
PubMed 
Article 

Google Scholar 
69.Nguyen, L. T., Schmidt, H. A., Von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
PubMed 
Article 

Google Scholar 
70.Hasegawa, M., Kishino, H. & Yano, T. A. Dating of the humanape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 160–174 (1985).ADS 
CAS 
PubMed 
Article 

Google Scholar 
71.Bouckaert, R. et al. BEAST 2: A software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 10, e1003537 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
72.Bandelt, H. J., Forster, P. & Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).CAS 
PubMed 
Article 

Google Scholar 
73.Corander, J., Marttinen, P., Sirén, J. & Tang, J. Enhanced Bayesian modelling in BAPS software for learning genetic struc-tures of populations. BMC Bioinform. 9, 1–14 (2008).Article 
CAS 

Google Scholar 
74.Excoffier, L. & Lischer, H. E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).PubMed 
Article 

Google Scholar 
75.Heled, J. & Drummond, A. J. Bayesian inference of population size history from multiple loci. BMC Evol. Biol. 8, 289 (2008).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
76.Team, R.C. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2014).
Google Scholar 
77.Fu, Y. X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147, 915–925 (1997).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
78.Tajima, F. The effect of change in population size on DNA polymorphism. Genetics 123, 597–601 (1989).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
79.Rousset, F. Genepop’007: A complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 8, 103–106 (2008).PubMed 
Article 

Google Scholar 
80.Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 6, 65–70 (1979).MathSciNet 
MATH 

Google Scholar 
81.Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
82.Hubisz, M. J., Falush, D., Stephens, M. & Pritchard, J. K. Inferring weak population structure with the assistance of sample group information. Mol. Ecol. Resour. 9, 1322–1332 (2009).PubMed 
PubMed Central 
Article 

Google Scholar 
83.Earl, D. A. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361 (2012).Article 

Google Scholar 
84.Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).CAS 
PubMed 
Article 

Google Scholar 
85.Langguth, T. et al. Genetic structure and phylogeogra-phy of a European flagship species, the white-tailed sea eagle Haliaeetus albicilla. J. Avian Biol. 44, 263–271 (2013).Article 

Google Scholar 
86.Väli, Ü., Dombrovski, V., Dzmitranok, M., Maciorowski, G. & Meyburg, B. U. High genetic diversity and low differentia-tion retained in the European fragmented and declining Greater Spotted Eagle (Clanga clanga) population. Sci. Rep. 9, 1–11 (2019).Article 
CAS 

Google Scholar 
87.Jakobsson, M. & Rosenberg, N. A. CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23, 1801–1806 (2007).CAS 
PubMed 
Article 

Google Scholar 
88.Rosenberg, N. A. DISTRUCT: A program for the graphical display of population structure. Mol. Ecol. Notes. 4, 137–138 (2004).Article 

Google Scholar 
89.Corander, J., Sirén, J. & Arjas, E. Bayesian spatial modeling of genetic population structure. Comput. Stat. 23, 111–129 (2008).MathSciNet 
Article 

Google Scholar 
90.Galpern, P., Peres-Neto, P. R., Polfus, J. & Manseau, M. MEMGENE: Spatial pattern detection in genetic distance data. Methods Ecol. Evol. 5, 1116–1120 (2014).Article 

Google Scholar 
91.Peakall, R., Ruibal, M. & Lindenmayer, D. B. Spatial autocorrelation analysis offers new insights into gene flow in the Aus-tralian bush rat, Rattus fuscipes. Evolution 57, 1182–1195 (2003).PubMed 
Article 

Google Scholar 
92.Mullins, J. et al. The influence of habitat structure on genetic differentiation in red fox populations in north-eastern Poland. Acta Theriol. 59, 367–376 (2014).Article 

Google Scholar 
93.Oksanen, J. Vegan: R Functions for Vegetation Ecologists. http://cc.oulu.fi/_jarioksa/softhelp/vegan.html. (2005).94.Wright, S. Isolation by distance. Genetics 28, 114 (1943).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
95.Peakall, R. & Smouse, P. E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and re-search: An update. Bioinformatics 28, 2537–2539 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
96.Goudet, J. FSTAT, A Program to Estimate and Test Gene Diversities and Fixation Indices, Version 2.9. 3. http://www2.unil.ch/popgen/softwares/fstat.htm (2001).97.Yeh, F. C. et al. POPGENE version 1.32. Computer Program and Documentation Distributed by the Author. http://www.ualberta.ca/~fyeh/popgene.html (accessed on 23 January 2013).98.Belkhir, K., Borsa, P., Chikhi, L. & Bonhomme, F. Genetix 4.05: WindowsTM Software for Population Genetics (Laboratoire Genome de Populations University of Montpelier II, 1996).
Google Scholar 
99.Luikart, G., Sherwin, W. B., Steele, B. M. & Allendorf, F. W. Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change. Mol. Ecol. 7, 963–974 (1998).CAS 
PubMed 
Article 

Google Scholar 
100.Luikart, G., Allendorf, F. W., Cornuet, J. M. & Sherwin, W. B. Distortion of allele frequency distributions provides a test for recent population bottlenecks. J. Hered. 89, 238–247 (1998).CAS 
PubMed 
Article 

Google Scholar  More

1.Gallo Marin, B. et al. Predictors of COVID-19 severity: A literature review. Rev. Med. Virol. 31, 1–10 (2021).CAS 
PubMed 
Article 

Google Scholar 
2.Worldometers. COVID-19 Coronavirus Pandemic. https://www.worldometers.info/coronavirus/ (Accessed 5 May 2021) (2021).3.The New York Times. Coronavirus World Map: Tracking the Global Outbreak. https://www.nytimes.com/interactive/2020/world/coronavirus-maps.html (Accessed 5 May) (2021).4.Sartorius, B., Lawson, A. & Pullan, R. Modelling and predicting the spatio-temporal spread of COVID-19, associated deaths and impact of key risk factors in England. Sci. Rep. 11, 5378 (2021).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
5.Gamio, L. & Symonds, A. Global Virus Cases Reach New Peak, Driven by India and South America. https://nyti.ms/3xYVO94 (Accessed on 5 May 2021) (2021).6.Samuel, J. et al. COVID-19 public sentiment insights and machine learning for tweets classification. Information 11, 314 (2020).Article 

Google Scholar 
7.Borriello, A., Master, D., Pellegrini, A. & Rose, J. M. Preferences for a COVID-19 vaccine in Australia. Vaccine 39, 473–479 (2021).CAS 
PubMed 
Article 

Google Scholar 
8.Samuel, J. et al. Feeling positive about reopening? New normal scenarios from COVID-19 US reopen sentiment analytics. IEEE Access 8, 142173–142190 (2020).Article 

Google Scholar 
9.Max Roser, E. O.-O., Ritchie, H. & Hasell, J. Coronavirus pandemic (COVID-19). https://ourworldindata.org/coronavirus (Accessed on 5 June 2021) (2020).10.Centers for Disease Control and Prevention. About variants of the virus that causes COVID-19. https://www.cdc.gov/coronavirus/2019-ncov/transmission/variant.html (Accessed on 5 June 2021) (2021).11.Ali, G. G. M. N. et al. Public perceptions of COVID-19 vaccines: Policy implications from US spatiotemporal sentiment analytics. Healthcare 9, 1110 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
12.Al Zobbi, M., Alsinglawi, B., Mubin, O. & Alnajjar, F. Measurement method for evaluating the lockdown policies during the COVID-19 pandemic. Int. J. Environ. Res. Public Health 17, 5574 (2020).CAS 
PubMed Central 
Article 
PubMed 

Google Scholar 
13.Lu, X., Yuan, D., Chen, W. & Fung, J. A machine learning based forecast model for the COVID-19 pandemic and investigation of the impact of government intervention on COVID-19 transmission in China (2020). Preprint on webpage at www.researchsquare.com/article/rs-73671/v1.14.Vinceti, M. et al. Lockdown timing and efficacy in controlling COVID-19 using mobile phone tracking. EClinicalMedicine 25, 100457 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
15.Rahman, M. et al. Machine learning on the COVID-19 pandemic, human mobility and air quality: A review. IEEE Access 9, 72420–72450 (2021).Article 

Google Scholar 
16.Rahman, M. et al. COVID-19 pandemic severity, lockdown regimes, and people’s mobility: Early evidence from 88 countries. Sustainability 12, 9101 (2020).Article 

Google Scholar 
17.Silva, P. C. et al. COVID-ABS: An agent-based model of COVID-19 epidemic to simulate health and economic effects of social distancing interventions. Chaos Solitons Fractals 139, 110088 (2020).MathSciNet 
PubMed 
PubMed Central 
Article 

Google Scholar 
18.Rahman, M. M. et al. Socioeconomic factors analysis for COVID-19 US reopening sentiment with Twitter and census data. Heliyon 7, e06200 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
19.USAFacts. US COVID-19 cases and deaths by state. https://usafacts.org/visualizations/coronavirus-covid-19-spread-map/ (Accessed 5 June 2021) (2020).20.Wynants, L. et al. Prediction models for diagnosis and prognosis of covid-19: Systematic review and critical appraisal. BMJ. 369(8242), m1328 (2020).21.Swapnarekha, H., Behera, H. S., Nayak, J. & Naik, B. Role of intelligent computing in COVID-19 prognosis: A state-of-the-art review. Chaos Solitons Fractals 138, 109947 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
22.Xiang, Y. et al. COVID-19 epidemic prediction and the impact of public health interventions: A review of COVID-19 epidemic models. Infect. Dis. Model. 6, 324–342 (2021).23.Roy, A. & Kar, B. Characterizing the spread of COVID-19 from human mobility patterns and SocioDemographic indicators. In Proceedings of the 3rd ACM SIGSPATIAL International Workshop on Advances in Resilient and Intelligent Cities, 39–48 (2020).24.Scarpone, C. et al. A multimethod approach for county-scale geospatial analysis of emerging infectious diseases: A cross-sectional case study of COVID-19 incidence in Germany. Int. J. Health Geogr. 19, 32 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
25.Polyzos, S., Samitas, A. & Spyridou, A. E. Tourism demand and the COVID-19 pandemic: An LSTM approach. Tour. Recreat. Res. 46(2), 1777053 (2020).26.Iwendi, C. et al. COVID-19 patient health prediction using boosted random forest algorithm. Front. Public Health 8, 357 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
27.Hou, X. et al. Intra-county modeling of COVID-19 infection with human mobility: assessing spatial heterogeneity with business traffic, age and race (2020). Preprint on webpage at https://doi.org/10.1101/2020.10.04.20206763v1.28.Tang, B. et al. An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov). Infect. Dis. Model. 5, 248–255 (2020).PubMed 
PubMed Central 

Google Scholar 
29.Zhao, S. et al. Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak. Int. J. Infect. Dis. 92, 214–217 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
30.Fanelli, D. & Piazza, F. Analysis and forecast of COVID-19 spreading in China, Italy and France. Chaos Solitons Fractals 134, 109761 (2020).MathSciNet 
PubMed 
PubMed Central 
Article 

Google Scholar 
31.Choi, S. & Ki, M. Estimating the reproductive number and the outbreak size of COVID-19 in Korea. Epidemiol. Health 42, e2020011 (2020).32.Tolles, J. & Luong, T. Modeling epidemics with compartmental models. JAMA 323, 2515–2516 (2020).PubMed 
Article 

Google Scholar 
33.Soures, N. et al. SIRNet: Understanding social distancing measures with hybrid neural network model for COVID-19 infectious spread. arXiv (2020). preprint on webpage at arXiv:2004.10376.34.Chimmula, V. K. R. & Zhang, L. Time series forecasting of COVID-19 transmission in Canada using LSTM networks. Chaos Solitons Fractals 135, 109864 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
35.Wang, D. et al. Agent-based Simulation Model and Deep Learning Techniques to Evaluate and Predict Transportation Trends around COVID-19. arXiv (2020). Preprint on webpage at arXiv:2010.09648.36.Kai, D., Goldstein, G.-P., Morgunov, A., Nangalia, V. & Rotkirch, A. Universal masking is urgent in the covid-19 pandemic: Seir and agent based models, empirical validation, policy recommendations. arXiv (2020). Preprint on webpage at arXiv:2004.13553.37.Panovska-Griffiths, J. et al. Modelling the potential impact of mask use in schools and society on COVID-19 control in the UK. Sci. Rep. 11, 8747 (2021).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
38.Szczepanek, R. Analysis of pedestrian activity before and during COVID-19 lockdown, using webcam time-lapse from cracow and machine learning. PeerJ 8, e10132 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
39.Ahmed, I., Ahmad, M., Rodrigues, J. J., Jeon, G. & Din, S. A deep learning-based social distance monitoring framework for COVID-19. Sustain. Cities Soc. 65, 102571 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
40.Spada, A. et al. Structural equation modeling to shed light on the controversial role of climate on the spread of SARS-CoV-2. Sci. Rep. 11, 8358 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
41.Jiang, F. et al. Artificial intelligence in healthcare: Past, present and future. Stroke Vasc. Neurol. 2, e000101 (2017).Article 

Google Scholar 
42.Khalifa, N. E. M., Taha, M. H. N., Ali, D. E., Slowik, A. & Hassanien, A. E. Artificial intelligence technique for gene expression by tumor RNA-Seq data: A novel optimized deep learning approach. IEEE Access 8, 22874–22883 (2020).Article 

Google Scholar 
43.Dexter, G. P., Grannis, S. J., Dixon, B. E. & Kasthurirathne, S. N. Generalization of machine learning approaches to identify notifiable conditions from a statewide health information exchange. AMIA Summits Transl. Sci. Proc. 2020, 152–161 (2020).PubMed 
PubMed Central 

Google Scholar 
44.Wang, Y., Liao, Z., Mathieu, S., Bin, F. & Tu, X. Prediction and evaluation of plasma arc reforming of naphthalene using a hybrid machine learning model. J. Hazard. Mater. 404, 123965 (2021).CAS 
PubMed 
Article 

Google Scholar 
45.Shao, Y. E., Hou, C.-D. & Chiu, C.-C. Hybrid intelligent modeling schemes for heart disease classification. Appl. Soft Comput. 14, 47–52 (2014).Article 

Google Scholar 
46.Roberts, M. et al. Common pitfalls and recommendations for using machine learning to detect and prognosticate for COVID-19 using chest radiographs and CT scans. Nat. Mach. Intell. 3, 199–217 (2021).Article 

Google Scholar 
47.Sun, J. et al. Forecasting the long-term trend of COVID-19 epidemic using a dynamic model. Sci. Rep. 10, 1–10 (2020).ADS 
CAS 
Article 

Google Scholar 
48.Chen, Y., Qin, R., Zhang, G. & Albanwan, H. Spatial temporal analysis of traffic patterns during the COVID-19 epidemic by vehicle detection using planet remote-sensing satellite images. Remote Sens. 13, 208 (2021).ADS 
Article 

Google Scholar 
49.Briz-Redón, Á. & Serrano-Aroca, Á. A spatio-temporal analysis for exploring the effect of temperature on COVID-19 early evolution in Spain. Sci. Total Environ. 728, 138811 (2020).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
50.Liu, Q. et al. Spatiotemporal patterns of COVID-19 impact on human activities and environment in mainland China using nighttime light and air quality data. Remote Sens. 12, 1576 (2020).ADS 
Article 

Google Scholar 
51.Jarvis, K. F. & Kelley, J. B. Temporal dynamics of viral load and false negative rate influence the levels of testing necessary to combat COVID-19 spread. Sci. Rep. 11, 1–12 (2021).Article 
CAS 

Google Scholar 
52.Ugarte, M. D., Adin, A., Goicoa, T. & Militino, A. F. On fitting spatio-temporal disease mapping models using approximate Bayesian inference. Stat. Methods Med. Res. 23, 507–530 (2014).MathSciNet 
PubMed 
Article 

Google Scholar 
53.Zhang, C. H. & Schwartz, G. G. Spatial disparities in coronavirus incidence and mortality in the United States: An ecological analysis as of May 2020. J. Rural Health 36, 433–445 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
54.Fitzpatrick, K. M., Harris, C. & Drawve, G. Fear of COVID-19 and the mental health consequences in America. Psychol. Trauma Theory Res. Pract. Policy 12, S17–S21 (2020).Article 

Google Scholar 
55.Lalmuanawma, S., Hussain, J. & Chhakchhuak, L. Applications of machine learning and artificial intelligence for covid-19 (SARS-CoV-2) pandemic: A review. Chaos Solitons Fractals 139, 110059 (2020).MathSciNet 
PubMed 
PubMed Central 
Article 

Google Scholar 
56.Adhikari, B., Xu, X., Ramakrishnan, N. & Prakash, B. A. Epideep: Exploiting embeddings for epidemic forecasting. In Proceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining, 577–586 (2019).57.Gautam, Y. Transfer learning for COVID-19 cases and deaths forecast using LSTM network. ISA Transactions (2021).58.Chen, S. et al. Exploring feasibility of multivariate deep learning models in predicting covid-19 epidemic. Front. Public Health 9, 661615 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
59.Venna, S. R. et al. A novel data-driven model for real-time influenza forecasting. IEEE Access 7, 7691–7701 (2018).Article 

Google Scholar 
60.Wu, D. et al. Deepgleam: a hybrid mechanistic and deep learning model for covid-19 forecasting. arXiv (2021). Preprint on webpage at arXiv:2102.06684.61.Zhu, X. et al. Attention-based recurrent neural network for influenza epidemic prediction. BMC Bioinform. 20, 575 (2019).Article 

Google Scholar 
62.Ben Said, A., Erradi, A., Aly, H. & Mohamed, A. Predicting covid-19 cases using bidirectional lstm on multivariate time series. arXiv (2020). Preprint on webpage at arXiv:2009.12325.63.Aktay A., A. et al. Google COVID-19 Community Mobility Reports: anonymization process description (version 1.1). arXiv (2020). Preprint on webpage at arXiv:2004.04145.64.SafeGraph. SafeGraph Places Schema. https://docs.safegraph.com/docs (Accessed 15 Sept 2020) (2020).65.Johns Hopkins University. COVID-19 dashboard by the center for systems science and engineering (CSSE) at Johns Hopkins University (JHU). https://coronavirus.jhu.edu/us-map (Accessed 15 Sept 2020) (2020).66.Ray, E. L. et al. Ensemble forecasts of coronavirus disease 2019 (covid-19) in the us. MedRXiv (2020). Preprint on webpage at https://www.medrxiv.org/content/10.1101/2020.08.19.20177493v1.67.Cartenì, A., Di Francesco, L. & Martino, M. How mobility habits influenced the spread of the COVID-19 pandemic: Results from the Italian case study. Sci. Total Environ. 741, 140489 (2020).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
68.Reich, N. et al. har96, x. Zhang, jinghuichen, G. Espana, X. Xinyue, H. Biegel, L. Castro, Y. Wang, qjhong, E. Lee, A. Baxter, S. Bhatia, E. Ray, and abrennen, and ERDC CV19 Modeling Team (2020). Preprint on webpage at https://github.com/reichlab/covid19-forecast-hub/tree/master/data-processed/COVIDhub-ensemble.69.Hochreiter, S. & Schmidhuber, J. Long short-term memory. Neural Comput. 9, 1735–1780 (1997).CAS 
PubMed 
Article 

Google Scholar 
70.Yao, W., Huang, P. & Jia, Z. Multidimensional lstm networks to predict wind speed. In 2018 37th Chinese Control Conference (CCC), 7493–7497 (IEEE, 2018).71.Cromley, E. K. & McLafferty, S. L. GIS and Public Health (Guilford Press, XXX, 2011).
Google Scholar 
72.Veličković, P. et al. Cross-modal recurrent models for weight objective prediction from multimodal time-series data. In Proceedings of the 12th EAI International Conference on Pervasive Computing Technologies for Healthcare, 178–186 (2018).73.Shi, X. & Yeung, D.-Y. Machine learning for spatiotemporal sequence forecasting: A survey. arXiv (2018). Preprint on webpage at arXiv:1808.06865.74.Bengio, S., Vinyals, O., Jaitly, N. & Shazeer, N. Scheduled sampling for sequence prediction with recurrent neural networks. arXiv (2015). Preprint on webpage at arXiv:1506.03099.75.Taieb, S. B. & Hyndman, R. Boosting multi-step autoregressive forecasts. In International Conference on Machine Learning, Vol. 32, 109–117 (PMLR, 2014). http://proceedings.mlr.press/v32/taieb14.html.76.Veličković, P. et al. Graph attention networks. arXiv (2017). Preprint on webpage at arXiv:1710.10903.77.Hamilton, W. L., Ying, R. & Leskovec, J. Inductive representation learning on large graphs. arXiv (2017). Preprint on webpage at arXiv:1706.02216. More