Philippa Kaur

Abbott RJ, Brennan AC (2014) Altitudinal gradients, plant hybrid zones and evolutionary novelty. Philos Trans R Soc B Biol Sci 369:6–9Article 

Google Scholar 
Avise JC (2000) Phylogeography: the history and formation of species. Harvard University Press, Cambridge, MABook 

Google Scholar 
Baldassarre DT, White TA, Karubian J, Webster MS (2014) Genomic and morphological analysis of a semipermeable avian hybrid zone suggests asymmetrical introgression of a sexual signal. Evolution 68:2644–2657PubMed 
Article 
PubMed Central 

Google Scholar 
Barr KR, Dharmarajan G, Rhodes OE, Lance R, Leberg PL (2007) Novel microsatellite loci for the study of the black-capped vireo (Vireo atricapillus). Mol Ecol Notes 7:1067–1069CAS 
Article 

Google Scholar 
Barton NH, Gale KS (1993) Hybrid zones and the evolutionary process. In: Harrison RG (ed.) Hybrid Zones and the Evolutionary Process. Oxford University Press, New York, NY
Google Scholar 
Barton NH, Hewitt GM (1989) Adaption, speciation and hybrid zones. Nature 341:497–503CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Billerman SM, Murphy MA, Carling MD (2016) Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement. Ecol Evol 6:7976–7990PubMed 
PubMed Central 
Article 

Google Scholar 
Bell RC, Irian CG (2019) Phenotypic and genetic divergence in reed frogs across a mosaic hybrid zone on São Tomé Island. Biol J Linn Soc 128:672–680Article 

Google Scholar 
Bensch S, Price T, Kohn J (1997) Isolation and characterization of microsatellite loci in a Phylloscopus warbler. Mol Ecol 6:91–92CAS 
PubMed 
Article 

Google Scholar 
Bradbury IR, Bowman S, Borza T, Snelgrove PVR, Hutchings JA, Berg PR et al. (2014) Long distance linkage disequilibrium and limited hybridization suggest cryptic speciation in Atlantic cod. PLoS ONE 9:e106330Article 
CAS 

Google Scholar 
Brelsford A, Irwin DE (2009) Incipient speciation despite little assortative mating: the yellow-rumped warbler hybrid zone. Evolution 63:3050–3060PubMed 
Article 

Google Scholar 
Burg TM, Croxall JP (2004) Global population structure and taxonomy of the wandering albatross species complex. Mol Ecol 13:2345–2355CAS 
PubMed 
Article 

Google Scholar 
Carling MD, Zuckerberg B (2011) Spatio-temporal changes in the genetic structure of the Passerina bunting hybrid zone. Mol Ecol 20:1166–1175PubMed 
Article 

Google Scholar 
Carling MD, Thomassen HA (2012) The role of environmental heterogeneity in maintaining reproductive isolation between hybridizing Passerina (Aves: Cardinalidae) buntings. Int J Ecol 2012:295463Article 

Google Scholar 
Carpenter AM, Graham BA, Spellman GM, Klicka J, Burg TM (2021) Genetic, bioacoustic and morphological analyses reveal cryptic speciation in the warbling vireo complex (Vireo gilvus: Vireonidae: Passeriformes). Zool J Linn Soc zlab036 https://doi.org/10.1093/zoolinnean/zlab036Cicero C, Johnson NK (1998) Molecular phylogeny and ecological diversification in a clade of New World songbirds (genus Vireo). Mol Ecol 7:1359–1370CAS 
PubMed 
Article 

Google Scholar 
Chenuil A, Cahill AE, Délémontey N, Du Salliant du Luc E, Fanton H (2019) Problems and questions posed by cryptic species. A framework to guide future studies. Assessing to conserving biodiversity. History, philosophy and theory of the life sciences, Vol. 24. Springer. Daubenmire, Cham
Google Scholar 
Cheviron ZA, Brumfield RT (2012) Genomic insights into adaptation to high-altitude environments. Heredity 108:354–361CAS 
PubMed 
Article 

Google Scholar 
Coyne JA, Orr HA (2004) Speciation. Sinauer and Associates, Sunderland, Massachusetts
Google Scholar 
Culumber ZW, Shepard DB, Colemans SW, Rosenthal GG, Tobler M (2012) Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus). J Evol Biol 25:1800–1814CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Dubay SG, Witt CC (2014) Differential high-altitude adaptation and restricted gene flow across a mid-elevation hybrid zone in Andean tit-tyrant flycatchers. Mol Ecol 23:3551–3565PubMed 
Article 
PubMed Central 

Google Scholar 
Garroway CJ, Bowman J, Cascaden TJ, Holloway GL, Mahan CG, Malcolm JR et al. (2010) Climate change induced hybridization in flying squirrels. Glob Chang Biol 16:113–121Article 

Google Scholar 
Grabenstein KC, Taylor SA (2018) Breaking barriers: Causes, consequences, and experimental utility of human-mediated hybridization. Trends Ecol Evol 33:198–212PubMed 
Article 

Google Scholar 
Graham BA, Cicero C, Strickland D, Woods JG, Coneybeare H, Dohms KM et al. (2021) Cryptic genetic diversity and cytonuclear discordance characterize contact among Canada jay (Perisoreus canadensis) morphotypes in western North America. Biol J Linn Soc 132:725–740Article 

Google Scholar 
Hammer Ø, Harper DA, Ryan PD (2001) Paleontological statistics software package for education and data analysis. Palaeontol Electron 4:9Haselhorst MSH, Parchman TL, Buerkle CA (2019) Genetic evidence for species cohesion, substructure and hybrids in spruce. Mol Ecol 28:2029–2045PubMed 
Article 

Google Scholar 
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25:1965–1978Article 

Google Scholar 
Hawley DM (2005) Isolation and characterization of eight microsatellite loci from the house finch (Carpodactus mexicanus). Mol Ecol Notes 5:443–445CAS 
Article 

Google Scholar 
Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM (2004) Identification of birds through DNA barcodes. PLoS Biol 2:e312PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Hewitt GM (1988) Hybrid zones-natural laboratories for evolutionary studies. Trends Ecol Evol 3:158–167CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Hewitt GM (2001) Speciation, hybrid zones and phylogeography—or seeing genes in space and time. Mol Ecol 10:537–549CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. Int J Climatol 25:1965–1978Article 

Google Scholar 
Hindley JA, Graham BA, Pulgarin-R PC, Burg TM (2018) The influence of latitude, geographic distance, and habitat discontinuities on genetic variation in a high latitude montane species. Sci Rep. 8:11846CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hinojosa JC, Koubínová D, Szenteczki MA, Pitteloud C, Dincă V, Alvarez N et al. (2019) A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris. Mol Ecol 28:3857–3868PubMed 
Article 
PubMed Central 

Google Scholar 
Hubisz MJ, Falush D, Stephens M, Pritchard JK (2009) Inferring weak population structure with the assistance of sample group information. Mol Ecol Res 9:1322–1332Article 

Google Scholar 
Irwin DE (2020) Assortative mating in hybrid zones is remarkably ineffective in promoting speciation. Evolution 195:E150–E167
Google Scholar 
Johnson NK (1995) Speciation in vireos. I. Macrogeographic patterns of allozymic variation in the Vireo solitarius complex in the contiguous United States. Condor 97:903–919Article 

Google Scholar 
Johnson NK, Cicero C (2004) New mitochondrial DNA data affirm the importance of Pleistocene speciation in North American birds. Evolution 58:1122–1130PubMed 
Article 
PubMed Central 

Google Scholar 
Larson EL, Tinghitella RM, Taylor SA (2019) Insect hybridization and climate change. Front Ecol Evol 7:348Article 

Google Scholar 
Legendre P, Fortin M-J (2010) Comparison of the Mantel test and alternative approaches for detecting complex multivariate relationships in the spatial analysis of genetic data. Mol Ecol Resour 10:831–844PubMed 
Article 

Google Scholar 
Lovell SF, Lein MR, Rogers SM (2021) Cryptic speciation in the warbling vireo (Vireo gilvus). Ornithology 138:ukaa071Article 

Google Scholar 
MacDonald ZG, Dupuis JR, Davis CS, Acorn JH, Nielsen SE, Sperling FAH (2020) Gene flow and climate-associated genetic variation in a vagile habitat specialist. Mol Ecol 29:3889–3906PubMed 
Article 

Google Scholar 
Manthey JD, Klicka J, Spellman GM (2011) Cryptic diversity in a widespread North American songbird: phylogeography of the brown creeper (Certhia americana). Mol Phylogenet Evol 58:502–512PubMed 
Article 

Google Scholar 
Marchetti K, Price T, Richman A (1995) Correlates of wing morphology with foraging behaviour and migration distance in the genus Phylloscopus. J Av Biol 26:177–181Article 

Google Scholar 
Martin H, Touzet P, Dufay M, Gode C, Schmitt E, Lahiani E et al. (2017) Lineages of Silene nutans developed rapid, strong, asymmetric postzygotic reproductive isolation in allopatry. Evolution 71:1519–1531CAS 
PubMed 
Article 

Google Scholar 
Martinez JG, Soler JJ, Soler M, Moller AP, Burke T (1999) Comparative population structure and gene flow of a brood parasite, the great spotted cuckoo (Clamator glandarius) and its primary host, the magpie (Pica pica). Evolution 53:269–278CAS 
PubMed 
PubMed Central 

Google Scholar 
Mettler RD, Spellman GM (2009) A hybrid zone revisited: Molecular and morphological analysis of the maintenance, movement, and evolution of a Great Plains avian (Cardinalidae: Pheucticus) hybrid zone. Mol Ecol 18:3256–3267CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Meirmans PG, Van Tienderen PH (2004) GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. Mol Ecol Notes 4:792–794Article 

Google Scholar 
Nowakowski JK, Szulc J, Remisiewicz M (2014) The further the flight, the longer the wing: relationship between wing length and migratory distance in Old World reed and bush warblers (Acrocephalidae and Locustellidae). Ornis Fennica 91:178–186
Google Scholar 
Pavolova A, Amos JN, Joseph L, Loynes K, Austin JJ, Keogh JS et al. (2013) Perched at the mito-nuclear crossroads: divergent mitochondrial lineages correlate with environment in the face of ongoing nuclear gene flow in an Australian bird. Evol 67:3412–3428Article 
CAS 

Google Scholar 
Piertney SB, Marquiss M, Summers R (1998) Characterization of tetranucleotide microsatellite markers in the Scottish crossbill (Loxia scotica). Mol Ecol 7:1261–1263CAS 
PubMed 
Article 

Google Scholar 
Porras-Hurtado L, Ruiz Y, Santos C, Phillips C, Carracedo A, Lareu MV (2013) An overview of STRUCTURE: Applications, parameter settings, and supporting software. Front Genet 4:98PubMed 
PubMed Central 
Article 
CAS 

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

Google Scholar 
Reding DM, Castañeda-Rico S, Shirazi S, Hofman CA, Cancellare IA, Lance SL et al. (2021) Mitochondrial genomes of the United States distribution of gray fox (Urocyon cinereoargenteus) reveal a major phylogeographic break at the Great Plains suture zone. Front Ecol Evol. https://doi.org/10.3389/fevo.2021.666800.Richardson DS, Jury FL, Dawson DA, Salgueiro P, Komdeur J, Burke T (2003) Fifty Seychelles warbler (Acrocephalus sechellensis) microsatellite loci polymorphic in Sylviidae species and their cross-species amplification in other passerine birds. Mol Ecol 9:2225–2230Article 

Google Scholar 
Riordan EC, Gugger PF, Ortego J, Smith C, Gaddis K, Thompson P et al. (2016) Association of genetic and phenotypic variability with geography and climate in three southern California oaks. Am J Bot 103:73–85PubMed 
Article 

Google Scholar 
Rush AC, Cannings RJ, Irwin DE (2009) Analysis of multilocus DNA reveals hybridization in a contact zone between Empidonax flycatchers. J Avian Biol 40:614–624Article 

Google Scholar 
Sartor CC, Cushman SA, Wan HY, Kretschmer R, Pereira JA, Bou N et al. (2021) The role of the environment in the spatial dynamics of an extensive hybrid zone between two neotropical cats. J Evol Biol 34:614–627PubMed 
Article 

Google Scholar 
Schukman JM, Lira-Noriega A, Townsend Peterson A (2011) Multiscalar ecological characterization of Say’s and eastern phoebes and their zone of contact in the Great Plains. Condor 113:372–384Article 

Google Scholar 
Seehausen O, Takimoto G, Roy D, Jokela J (2008) Speciation reversal and biodiversity dynamics with hybridization in changing environments. Mol Ecol 17:30–44PubMed 
Article 
PubMed Central 

Google Scholar 
Semenchuk GP (1992) The Atlas of Breeding Birds of Alberta. Fed. of Alberta Naturalists, Edmonton, p 243
Google Scholar 
Peakall R, Smouse PE (2012) GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research–an update. Bioinformatics 28:2537–2539CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sorenson MD, Ast JC, Dimcheff DE, Yuri T, Mindell DP (1999) Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates. Mol Phylogent Evol 12:105–114CAS 
Article 

Google Scholar 
Spellman GM, Klicka J (2007) Phylogeography of the white-breasted nuthatch (Sitta carolinensis): diversification in North American pine and oak woodlands. Mol Ecol 16:1729–1740CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Stenzler LM, Fitzpatrick JW (2002) Isolation of microsatellite loci in the Florida scrub jay Aphelocoma coerulescens. Mol Ecol Notes 2:547–550CAS 
Article 

Google Scholar 
Swenson NG (2006) GIS-based niche models reveal unifying climatic mechanisms that maintain location of avian hybrid zones in a North America suture zone. J Evol Biol. 19:717–725CAS 
PubMed 
Article 

Google Scholar 
Swenson NG, Howard DJ (2005) Clustering of contact zones, hybrid zones, and phylogeographic breaks in North America. Am Nat 166:581–591PubMed 
Article 

Google Scholar 
Tarr CL, Fleischer RC (1998) Primers for polymorphic GT microsatellites isolated from the Mariana crow, Corvus kubaryi. Mol Ecol 7:253–255CAS 
PubMed 
Article 

Google Scholar 
Tarroso P, Pereira RJ, Martínez-Freiría F, Godinho R, Brito JC (2014) Hybridization at an ecotone: Ecological and genetic barriers between three Iberian vipers. Mol Ecol 23:1108–1123CAS 
PubMed 
Article 

Google Scholar 
Taylor SA, Larson EL, Harrison RG (2015) Hybrid zones: windows on climate change. Trends Ecol Evol 30:398–406PubMed 
PubMed Central 
Article 

Google Scholar 
Toews DPL, Mandic M, Richards JG, Irwin DE (2014) Migration, mitochondria and the yellow-rumped warbler. Evolution 68:241–255CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Toews DPL, Campagna L, Taylor SA, Balakrishnan CN, Baldassarre DT, Deane-Coe PE et al. (2016) Genomic approaches to understanding population divergence and speciation in birds. Auk 133:13–30Article 

Google Scholar 
Toews DPL, Irwin DE (2008) Cryptic speciation in a Holarctic passerine revealed by genetic and bioacoustic analyses. Mol Ecol 17:2691–2705CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
van Els P, Cicero C, Klicka J (2012) High latitudes and high genetic diversity: Phylogeography of a widespread boreal bird, the gray jay (Perisoreus canadensis). Mol Phylogenet Evol 63:456–465PubMed 
Article 
PubMed Central 

Google Scholar 
Voelker G, Rohwer S (1998) Contrasts in scheduling of molt and migration in eastern and western warbling vireos. Auk 155:142–155Article 

Google Scholar 
Walsh J, Billerman SM, Rohwer VG, Butcher BG, Lovette IJ (2020) Genomic and plumage variation across the controversial Baltimore and Bullock’s oriole hybrid zone. Auk 137:1–15Article 

Google Scholar 
Walsh J, Rowe RJ, Olsen BJ, Shriver WG, Kovach AI (2016) Genotype-environment associations support a mosaic hybrid zone between two tidal marsh birds. Ecol Evol 6:279–294PubMed 
Article 
PubMed Central 

Google Scholar 
Walsh P, Metzger D, Higuchi R (1991) Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10:506–513CAS 
PubMed 
PubMed Central 

Google Scholar 
Weir JT, Schluter D (2004) Ice sheets promote speciation in boreal birds. Proc R Soc B 271:1881–1887PubMed 
PubMed Central 
Article 

Google Scholar 
Williams JW (2003) Variations in tree cover in North America since the last glacial maximum. Glob Planet Change 35:1–23Article 

Google Scholar 
Williams DA, Berg EC, Hale AM, Hughes CR (2004) Characterization of microsatellites for parentage studies of white-throated magpie-jays (Calocitta formosa) and brown jays (Cyanocorax morio). Mol Ecol Notes 4:509–511CAS 
Article 

Google Scholar 
Zwartjes PW (2001) Genetic structuring among migratory populations of the black-whiskered vireo, with a comparison to the red-eyed vireo. Condor 103:439–448Article 

Google Scholar  More

Phylogenetic analysis with multiple markersThe final alignment of 5670 bp length contained 843 variable sites (14.7%). Missing data accounted for 53.5% of the alignment cells and the relative GC content was 39.5%. Our phylogeny suggests that the investigated Holarctic taxa of the niger clade sensu Ref.34 are divided into two major clades with strong statistical support (Fig. 1). The first major clade (L. niger group) consists exclusively of Palearctic species (L. niger, L. platythorax, L. japonicus, L. emarginatus, L. balearicus, L. grandis, L. cinereus, the L. alienus-complex, L. sakagamii, L. productus and L. hayashi), with the exception of an unnamed Nearctic subclade recovered as sister to the rest of the group. This unnamed subclade we describe as a new species below (L. ponderosae sp. nov.). Lasius ponderosae sp. nov. corresponds to what was previously known as the Nearctic form of “L. niger” sensu ref.17, but includes some western Nearctic populations formerly assigned to “L. alienus”17,52 as well. Monophyly of L. ponderosae sp. nov. was fully supported by Bayesian inference (pp = 1) and moderately supported by maximum likelihood (66% bootstrap support, Fig. 1). Lasius ponderosae sp. nov. is distantly related to L. niger; and L.niger is a close relative of L. japonicus and L. platythorax, as well as other Palearctic taxa. The second major clade (L. brunneus group) within the investigated Holarctic members of the L. niger clade contains both Nearctic and Palearctic species not closely related to the taxa of interest (Fig. 1).Figure 1Molecular phylogeny of 26 Holarctic ant taxa belonging to the subgenus Lasius sensu Wilson (1955) and two outgroup taxa (L. pallitarsis and L. mixtus). The phylogeny was calculated under the coalescent model and incorporates data from 9 genes (mtDNA: COI, COII, 16S, nuDNA: Defensin, H3, LR, Wg, Top1 & 28S). Names of species native to the Nearctic are shown in red and those of species native to the Palearctic in blue. Node labels show posterior probability (Bayesian inference) followed by bootstrap support (Maximum likelihood). The scale bar indicates the length of 0.01 substitutions/site.Full size imageDNA-barcodingThe native North American species L. ponderosae sp. nov. contains at least 15 COI-mitotypes (n = 28 sequenced specimens) belonging to four distinct deep lineages, with divergences of up to 5.9%. Haplotype diversity was 0.899 and nucleotide diversity was 0.012. None of the mitotypes of this species was found to be widespread or particularly abundant. In striking contrast, low genetic diversity was found in L. niger across its entire distribution (Fig. 2). No more than 7 different COI-mitotypes were detected in samples from distant localities representing most of the known range (n = 70 specimens from 12 countries), from Spain in the West to the Siberian Baikal-region in the East (Fig. 2). Their maximum pairwise divergence was only 0.6%, with a haplotype diversity of 0.682 and a nucleotide diversity below 0.001. One mitotype of L. niger is highly dominant within the native range, occurring from Western Europe to Central Siberia (mitotype h2 in Fig. 2).Figure 2Mitotype tree and distribution maps for 98 DNA-barcodes belonging to 7 mitotypes of the ant Lasius niger (blue, n = 70) and 15 mitotypes of L. ponderosae sp. nov. (red, n = 28). The red dashed line delimits the expected natural range of L. ponderosae sp. nov.53 Maps have been created using the free R-package “ggmap” v3.0.0 (https://github.com/dkahle/ggmap) in R v4.1.1. Map tiles by Stamen Design, under CC BY 3.0.Full size imageRecent Palearctic L. niger introduction to CanadaPalearctic Lasius niger was introduced to several localities in coastal Canada in recent times, where at least 11 populations were found in two metropolitan areas (Vancouver and Halifax areas, see Table S2 for details). Those populations consist of the most dominant Palearctic mitotype of L. niger (h2). However, in 3 localities in the Vancouver area, 3 specimens with a second mitotype were found (mitotype h4, Fig. 2, Table S2) in syntopy with those carrying the most common mitotype h2. This second Canadian COI-mitotype (h4) was not found among our samples from the Old World, although it only differs by a single nucleotide substitution from mitotypes found there. A review of BOLD data revealed that the Canadian barcoded specimens of L. niger were mostly collected in anthropogenic habitats such as schoolyards (Supplementary Table S2).Description of Lasius ponderosae sp. novLasius ponderosae Schär, Talavera, Rana, Espadaler, Cover, Shattuck and Vila. ZooBank LSID: urn:lsid:zoobank.org:act:22E2743A-2F1C-4870-B318-A1F2DF2B464C Etymology: ponderosae alludes to the ponderosa pine tree (Pinus ponderosa) that is at the centre of occurrence in the ponderosa pine—gambel oak communities in the western Rocky Mountains and northern Arizona.Type material: located at the Museum of Comparative Zoology, Cambridge, USA. Two paratype workers each will be deposited at the collections of University of California Davis (UCDC), the University of Utah (JTLC) and the Natural History Museum of Los Angeles County (LACM).Holotype: worker, Fig. 3a–c. Type locality: USA, Utah: Uintah Co., Uintah Mtns., 2408 m. 18.6 mi N. Jct. Rt. 40 on Rt. 191, 40.66378°N, − 109.47918°E, leg. 15.VII.2013, S. P. Cover; J. D. Rana, collection code SPC 8571. Measurements [mm]: HL: 0.899, HW: 0.823, SL: 0.821, EL: 0.239, EW: 0.189, ProW: 0.56, ML: 1.069, HTL: 0.863, CI: 92, SI: 100.Figure 3Frontal, lateral and dorsal view of the holotype worker (a–c), a paratype gyne (d–f) and a paratype male of Lasius ponderosae sp. nov. (g–i).Full size imageParatypes: 15 workers, two gynes (Fig. 3d–f), two males (Fig. 3g–i) from the same series as the holotype, morphometric data is given in the Appendix, Table S5 and Table S6. CO1 mitotype h17: Genbank Accession no. LT977508.Description of the worker caste: A member of a complex of cryptic species resembling L. niger. Intermediate in overall body size, antennal scape length and eye size and comparable to related species (Table 1). Terminal segment of maxillary palps and torulo-clypeal distance relative to head size shorter than in related Palearctic species (Table 1). Mandibles with 8 or rarely 7 or 9 regular denticles and lacking offset teeth at their basal angle. Penultimate and terminal basal mandibular teeth of subequal size, and the gap in between with subequal area than the basal tooth. Anterior margin of clypeus evenly rounded. Dorsofrontal profile of pronotum slightly angular (Fig. 4a). Propodeal dome short and flat, usually lower than mesonotum (Fig. 4a). Body with abundant and long pilosity, especially lateral propodeum, genae, hind margin and underside of head. Pilosity of tibiae and antennal scapes variable, ranging from almost no setae (“L. alienus”-like phenotype) to very hairy (“L. niger”-like phenotype). Microscopic pubescent hairs on forehead between frontal carinae long and fine. Clypeus typically with only few scattered pubescent hairs (Figs. 3, 4c). Coloration of body dark brown, occasionally yellowish- or reddish-brown or slightly bicolored with head and thorax lighter than abdomen. Femora and antennal scapes brown. Mandibles and distal parts of legs yellowish to dark brown. Specimens of all 3 castes are shown in Fig. 3a–i and morphometric data are summarized in Table 1 and raw measurements are available in Table S5 and S6.Table 1 Morphometric data of Lasius ponderosae sp. nov. and comparison to morphologically similar Palearctic species.Full size tableFigure 4Average thorax profile of Lasius ponderosae sp. nov. (a) and members of the Palearctic L. niger-complex (b). Figures were created by image averaging (L. ponderosae sp. nov n = 35; Palearctic L. niger-complex n = 30 specimens). Frontal view of head and detail of clypeus of the Holotype worker of L. ponderosae sp. nov. (c) and a non-type worker of L. niger (d).Full size imageDiagnosis: Lasius ponderosae sp. nov. workers key out to “L. niger” using Wilson’s 1955 key to the Nearctic Lasius species. However, some populations with reduced pilosity may also be identified as “L. alienus” using this key. Lasius alienus is a Eurasian species not known from North America33. The Nearctic “L. alienus” sensu Wilson (1955) includes both, L. americanus as well as populations of L. ponderosae sp. nov. with sparse setae counts on tibia and/or scapes. Lasius ponderosae sp. nov. can be distinguished from L. americanus by the presence of abundant, long setae surpassing the sides of the head in full face view (nGen  > 5 and nOcc  > 10 vs. nGen  0.8 across models and runs). The strongest predictors were: Annual Mean Temperature (mean variable importance = 0.32), Mean Temperature of Coldest Quarter (0.23), Temperature Annual Range (0.23) and Temperature Seasonality (0.24). The contribution of land cover was low (0.02). The model predicted high probabilities of occurrence of L. niger in the eastern United States and southeastern Canada, including the island of Newfoundland, and small areas of suitable habitat in southwestern Canada and the Aleutians (Fig. 6). The area with high predicted occurrence probability of L. niger in the New World includes the two sites where populations have actually established (which were not used in the modeling): Nova Scotia and Vancouver. Further areas with high occurrence probabilities are New England, Southern Ontario, the Great Lakes-region and the Northern Appalachians. Low occurrence probabilities were found for the central North American prairies as well as arctic, boreal, arid, subtropical and tropical regions (Fig. 6). Considering the highest occurrence probability range (0.8–1 on a 0–1 probability scale), the area of suitable habitats for L. niger is 4,547,537 km2 in Europe and 1,308,920 km2 in North America. For an intermediate to high occurrence probability range (0.5–1) we estimated 5,371,055 km2 in Europe and 3,054,283 km2 in North America, and for the widest probability range (0.2–1) we estimated 6,155,643 km2 of suitable areas in Europe and 6,889,745 km2 in North America (Fig. 6).Figure 6Projected occurrence probability from ecological niche modeling for the Palearctic ant Lasius niger which has been introduced to Canada, based on 19 climatic and one land use variable. The intensity of blue colour indicates the probability of occurrence on a 0–1 scale based on 180 presences (black circles) and 182 absences (white circles) in the native range in the Old World (a). The model was then projected to North America to estimate areas of suitable habitat for this introduced species (b). These maps have been created using the free R-package “ggplot2” v3.3.5 (https://ggplot2.tidyverse.org) in R v4.1.1.Full size image More

Seymour, C. L. et al. Caught on camera: The impacts of urban domestic cats on wild prey in an African city and neighbouring protected areas. Glob. Ecol. Conserv. 23, e01198 (2020).Article 

Google Scholar 
Mori, E. et al. License to Kill? Domestic Cats Affect a Wide Range of Native Fauna in a Highly Biodiverse Mediterranean Country. Front. Ecol. Evol. 7, 477 (2019).Kays, R. et al. The small home ranges and large local ecological impacts of pet cats. Anim. Conserv. 23, 516–523 (2020).Loss, S. R., Will, T. & Marra, P. P. The impact of free-ranging domestic cats on wildlife of the United States. Nat. Commun. 4, 1396 (2013).ADS 
Article 

Google Scholar 
Van Heezik, Y., Smyth, A., Adams, A. & Gordon, J. Do domestic cats impose an unsustain386 able harvest on urban bird populations?. Biol. Conserv. 143, 121–130 (2010).Article 

Google Scholar 
Woods, M., McDonald, R. A. & Harris, S. Predation of wildlife by domestic cats Felis catus in Great Britain. Mammal Rev. 33, 174–188 (2003).Article 

Google Scholar 
Li, Y. et al. Estimates of wildlife killed by free-ranging cats in China. Biol. Conserv. 253, 108929 (2021).Article 

Google Scholar 
Barratt, D. G. Home range size, habitat utilisation and movement patterns of suburban and farm cats Felis catus. Ecography 20, 271–280 (1997).Article 

Google Scholar 
Moseby, K. E., Stott, J. & Crisp, H. Movement patterns of feral predators in an arid environment–implications for control through poison baiting. English. Wildl. Res. 36, 422–435 (2009).Article 

Google Scholar 
Hall, C. M. et al. Factors determining the home ranges of pet cats: A meta-analysis. Biol. Conserv. 203, 313–320 (2016).Article 

Google Scholar 
Castañeda, I. et al. Trophic patterns and home-range size of two generalist urban carnivores: A review. J. Zool. 307, 79–92 (2019).Article 

Google Scholar 
Hebblewhite, M. & Haydon, D. T. Distinguishing technology from biology: A critical review of the use of GPS telemetry data in ecology. Philos. Trans. R. Soc. B Biol. Sci. 365, 2303–2312 (2010).Article 

Google Scholar 
Allen, A. M. et al. Scaling up movements: From individual space use to population patterns. Ecosphere 7, e01524 (2016).
Google Scholar 
Trouwborst, A., McCormack, P. C. & Martínez Camacho, E. Domestic cats and their impacts on biodiversity: A blind spot in the application of nature conservation law. People Nat. 2, 235–250 (2020).Article 

Google Scholar 
Sims, V., Evans, K. L., Newson, S. E., Tratalos, J. A. & Gaston, K. J. Avian assemblage structure and domestic cat densities in urban environments. Divers. Distrib. 14, 387–399 (2008).Article 

Google Scholar 
Lepczyk, C. A., Mertig, A. G. & Liu, J. Landowners and cat predation across rural-to-urban landscapes. Biol. Conserv. 115, 191–201 (2004).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing R Foundation for Statistical Computing (Vienna, Austria, 2021).Heggøy, O. & Shimmings, P. Huskattens predasjon på fugler i Norge. En vurdering basert på en litteraturgjennomgang tech. rep. 36 (2018).Morgan, S. et al. Urban cat (Felis catus) movement and predation activity associated with a wetland reserve in New Zealand. Wildl. Res. 36, 574–580 (2009).Calver, M., Grayson, J., Lilith, M. & Dickman, C. Applying the precautionary principle to the issue of impacts by pet cats on urban wildlife. Biol. Conserv. 144, 1895–1901 (2011).Article 

Google Scholar 
Crowley, S., Cecchetti, M. & Mcdonald, R. Diverse perspectives of cat owners indicate bar riers to and opportunities for managing cat predation of wildlife. Front. Ecol. Environ. 18, 544–549 (2020).Treves, A., Krofel, M., Ohrens, O. & van Eeden, L. M. Predator control needs a standard of unbiased randomized experiments with cross-over design. Front. Ecol. Evol. 7, 462 (2019).Ferreira, G. A., Machado, J. C., Nakano-Oliveira, E., Andriolo, A. & Genaro, G. The effect of castration on home range size and activity patterns of domestic cats living in a natural area in a protected area on a Brazilian island. Appl. Anim. Behav. Sci. 230, 105049 (2020).Bengsen, A. J. et al. Feral cat home-range size varies predictably with landscape productivity and population density. J. Zool. 298, 112–120 (2016).Article 

Google Scholar 
López-Jara, M. J. et al. Free-roaming domestic cats near conservation areas in Chile: Spatial movements, human care and risks for wildlife. Perspect. Ecol. Conserv. 19, 387–398 (2021).Gillies, C. & Clout, M. The prey of domestic cats (Felis catus) in two suburbs of Auckland City, New Zealand. J. Zool. 259, 309–315 (2003).Article 

Google Scholar 
Pirie, T. J., Thomas, R. L. & Fellowes, M. D. E. Pet cats (Felis catus) from urban boundaries use different habitats, have larger home ranges and kill more prey than cats from the suburbs. Landsc. Urban Plan. 220, 104338 (2022).Article 

Google Scholar 
Vucetich, J. A., Hebblewhite, M., Smith, D. W. & Peterson, R. O. Predicting prey population dynamics from kill rate, predation rate and predator-prey ratios in three wolf-ungulate systems. J. Anim. Ecol. 80, 1236–1245 (2011).Article 

Google Scholar 
Kennedy, M., Phillips, B. E. N. L., Legge, S., Murphy, S. A. & Faulkner, R. A. Do dingoes suppress the activity of feral cats in northern Australia?. Austral Ecol. 37, 134–139 (2012).Article 

Google Scholar 
Crooks, K. R. & Soule, M. E. Mesopredator release and avifaunal extinctions in a fragmented system. English. Nature 400, 563–566 (1999).ADS 
CAS 
Article 

Google Scholar 
Ferreira, J. P., Leita, O. I., Santos-Reis, M. & Revilla, E. Human-related factors regulate the spatial ecology of domestic cats in sensitive areas for conservation. PLOS ONE 6, e25970 (2011).ADS 
CAS 
Article 

Google Scholar 
Brook, L. A., Johnson, C. N. & Ritchie, E. G. Effects of predator control on behaviour of an apex predator and indirect consequences for mesopredator suppression. J. Appl. Ecol. 49, 1278–1286 (2012).Article 

Google Scholar 
Laundre, J. W., Hernandez, L. & Altendorf, K. B. Wolves, elk, and bison: Reestablishing the “landscape of fear’’ in Yellowstone National Park, USA. English. Can. J. Zool. 79, 1401–1409 (2001).Article 

Google Scholar 
Ritchie, E. G. & Johnson, C. N. Predator interactions, mesopredator release and biodiversity conservation. English. Ecol. Lett. 12, 9820–998 (2009).Article 

Google Scholar 
Milleret, C. et al. GPS collars have an apparent positive effect on the survival of a large carnivore. Biol. Lett. 17, 20210128 (2021).Cecchetti, M., Crowley, S. L., Goodwin, C. E. D. & McDonald, R. A. Provision of high meat content food and object play reduce predation of wild animals by domestic cats Felis catus. Curr. Biol. 31, 1107-1111.e5 (2021).CAS 
Article 

Google Scholar 
Linklater, W., Farnworth, M., van Heezik, Y., Stafford, K. & Macdonald, E. Prioritizing cat owner behaviors for a campaign to reduce wildlife depredation. Conserv. Sci. Pract. 1, 1:e29 (2019).Selinske, M. J. et al. Identifying and prioritizing human behaviors that benefit biodiversity. Conserv. Sci. Pract. 2, e249 (2020).
Google Scholar 
McDonald, J. L., Maclean, M., Evans, M. R. & Hodgson, D. J. Reconciling actual and perceived rates of predation by domestic cats. Ecol. Evol. 5, 2745–2753 (2015).Article 

Google Scholar 
Bischof, R. et al. Estimating and forecasting spatial population dynamics of apex predators using transnational genetic monitoring. Proc. Natl. Acad. Sci. 117, 30531–30538 (2020).CAS 
Article 

Google Scholar 
Bischof, R., Gjevestad, J. G. O., Ordiz, A., Eldegard, K. & Milleret, C. High frequency GPS bursts and path-level analysis reveal linear feature tracking by red foxes. Sci. Rep. 9, 8849 (2019).ADS 
Article 

Google Scholar 
Gupte, P. R. et al. A guide to pre-processing high-throughput animal tracking data. J. Anim. Ecol. 91, 287–307 (2022).Article 

Google Scholar 
Morris, G. & Conner, L. Assessment of accuracy, fix success rate, and use of estimated horizontal position error (EHPE) to filter inaccurate data collected by a common commercially available GPS logger. PLoS ONE 12, e0189020 (2017).Article 

Google Scholar 
Clapp, J. G., Holbrook, J. D. & Thompson, D. J. GPSeqClus: An R package for sequential clustering of animal location data for model building, model application and field site investigations. Methods Ecol. Evol. 12, 787–793 (2021).Article 

Google Scholar 
Nielson, M., R., Sawyer, H. & McDonald, T. L. BBMM: Brownian Bridge Movement Model R Package Version 3.0 (2013).Horne, J. S., Garton, E. O., Krone, S. M. & Lewis, J. S. Analyzing animal movements using Brownian bridges. Ecology 88, 2354–2363 (2007).Article 

Google Scholar 
Sawyer, H., Kauffman, M. J., Nielson, R. M. & Horne, J. S. Identifying and prioritizing ungulate migration routes for landscape-level conservation. Ecol. Appl. 19, 2016–2025 (2009).Article 

Google Scholar 
Fischer, J. W., Walter, W. D. & Avery, M. L. Brownian bridge movement models to characterize birds’ home ranges. Condor 115, 298–305 (2013).Article 

Google Scholar 
Seidler, R., Long, R., Berger, J., Bergen, S. & Beckmann, J. Identifying impediments to long-distance mammal migrations. Conserv. Biol. 29 (2014).Collins, G. Seasonal distribution and routes of pronghorn in the Northern Great Basin. West. N. Am. Nat. 76, 101–112 (2016).Article 

Google Scholar  More

Paez-Espino D, Eloe-Fadrosh EA, Pavlopoulos GA, Thomas AD, Huntemann M, Mikhailova N, et al. Uncovering Earth’s virome. Nature. 2016;536:425–30.CAS 
PubMed 

Google Scholar 
Anderson PK, Cunningham AA, Patel NG, Morales FJ, Epstein PR, Daszak P. Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers. Trends Ecol Evol. 2004;19:535–44.PubMed 

Google Scholar 
Taylor LH, Latham SM, Woolhouse MEJ. Risk factors for human disease emergence. Philos Trans R Soc B Biol Sci. 2001;356:983–9.CAS 

Google Scholar 
White R, Murray S, Rohweder M. Pilot analysis of global ecosystems: grassland ecosystems. 2000 World Resources Institute. Washington, DC.Zhao Y, Liu Z, Wu J. Grassland ecosystem services: a systematic review of research advances and future directions. Landsc Ecol. 2020;35:793–814.
Google Scholar 
Trubl G, Jang HBin, Roux S, Emerson JB, Solonenko N, Vik DR, et al. Soil viruses are underexplored players in ecosystem carbon processing. mSystems. 2018;3:e00076–18.CAS 
PubMed 
PubMed Central 

Google Scholar 
Emerson JB, Roux S, Brum JR, Bolduc B, Woodcroft BJ, Jang HBin, et al. Host-linked soil viral ecology along a permafrost thaw gradient. Nat Microbiol. 2018;3:870–80.CAS 
PubMed 
PubMed Central 

Google Scholar 
Zablocki O, Adriaenssens EM, Frossard A, Seely M, Ramond J-B, Cowan D. Metaviromes of extracellular soil viruses along a Namib desert aridity gradient. Genome Announc. 2017;5:e01470–16.PubMed 
PubMed Central 

Google Scholar 
Jin M, Guo X, Zhang R, Qu W, Gao B, Zeng R. Diversities and potential biogeochemical impacts of mangrove soil viruses. Microbiome. 2019;7:58.PubMed 
PubMed Central 

Google Scholar 
Adriaenssens EM, Kramer R, Van Goethem MW, Makhalanyane TP, Hogg I, Cowan DA. Environmental drivers of viral community composition in Antarctic soils identified by viromics. Microbiome. 2017;5:83.PubMed 
PubMed Central 

Google Scholar 
Williamson KE, Fuhrmann JJ, Wommack KE, Radosevich M. Viruses in soil ecosystems: an unknown quantity within an unexplored territory. Annu Rev Virol. 2017;4:201–19.CAS 
PubMed 

Google Scholar 
Starr EP, Nuccio EE, Pett-Ridge J, Banfield JF, Firestone MK. Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil. Proc Natl Acad Sci. 2019;116:25900–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
Wu R, Davison MR, Gao Y, Nicora CD, Mcdermott JE, Burnum-Johnson KE, et al. Moisture modulates soil reservoirs of active DNA and RNA viruses. Commun Biol. 2021;4:1–11.
Google Scholar 
Hurwitz BL, Sullivan MB. The Pacific Ocean Virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS One. 2013;8:e57355.CAS 
PubMed 
PubMed Central 

Google Scholar 
Breitbart M, Bonnain C, Malki K, Sawaya NA. Phage puppet masters of the marine microbial realm. Nat Microbiol. 2018;3:754–66.CAS 
PubMed 

Google Scholar 
Wolf YI, Kazlauskas D, Iranzo J, Lucía-Sanz A, Kuhn JH, Krupovic M, et al. Origins and evolution of the Global RNA virome. MBio. 2018;9:e02329–18.PubMed 
PubMed Central 

Google Scholar 
Shi M, Lin XD, Tian JH, Chen LJ, Chen X, Li CX, et al. Redefining the invertebrate RNA virosphere. Nature. 2016;540:539–43.CAS 

Google Scholar 
Callanan J, Stockdale SR, Shkoporov A, Draper LA, Ross RP, Hill C. Expansion of known ssRNA phage genomes: from tens to over a thousand. Sci Adv. 2020;6:eaay5981.CAS 
PubMed 
PubMed Central 

Google Scholar 
Koonin EV, Dolja VV, Krupovic M, Varsani A, Wolf YI, Yutin N, et al. Global organization and proposed megataxonomy of the virus world. Microbiol Mol Biol Rev. 2020;84:e00061-19.PubMed 
PubMed Central 

Google Scholar 
Cobbin JC, Charon J, Harvey E, Holmes EC, Mahar JE. Current challenges to virus discovery by meta-transcriptomics. Curr Opin Virol. 2021;51:48–55.CAS 
PubMed 

Google Scholar 
Trubl G, Hyman P, Roux S, Abedon ST. Coming-of-age characterization of soil viruses: a user’s guide to virus isolation, detection within metagenomes, and viromics. Soil Syst. 2020;4:1–34. MDPI AG.
Google Scholar 
Santos-Medellin C, Zinke LA, ter Horst AM, Gelardi DL, Parikh SJ, Emerson JB. Viromes outperform total metagenomes in revealing the spatiotemporal patterns of agricultural soil viral communities. ISME J. 2021;15:1–15.
Google Scholar 
Adriaenssens EM, Farkas K, Harrison C, Jones DL, Allison HE, McCarthy AJ. Viromic analysis of wastewater input to a river catchment reveals a diverse assemblage of RNA viruses. mSystems. 2018;3:e00025–18.CAS 
PubMed 
PubMed Central 

Google Scholar 
Bibby K, Peccia J. Identification of viral pathogen diversity in sewage sludge by metagenome analysis. Environ Sci Technol. 2013;47:1945–51.CAS 
PubMed 
PubMed Central 

Google Scholar 
Culley A. New insight into the RNA aquatic virosphere via viromics. Virus Res. 2018;244:84–89.CAS 
PubMed 

Google Scholar 
Withers E, Hill PW, Chadwick DR, Jones DL. Use of untargeted metabolomics for assessing soil quality and microbial function. Soil Biol Biochem. 2020;143:107758.CAS 

Google Scholar 
Trubl G, Solonenko N, Chittick L, Solonenko SA, Rich VI, Sullivan MB. Optimization of viral resuspension methods for carbon-rich soils along a permafrost thaw gradient. PeerJ. 2016;4:e1999.PubMed 
PubMed Central 

Google Scholar 
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011;17:10.
Google Scholar 
Joshi N, Fass J. Sickle: a sliding-window, adaptive, quality-based trimming tool for FastQ files. 2011.Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.CAS 
PubMed 
PubMed Central 

Google Scholar 
Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28:3211–7.CAS 
PubMed 

Google Scholar 
Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.CAS 
PubMed 

Google Scholar 
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2014;12:59–60. Nature Publishing Group.PubMed 

Google Scholar 
Huson DH, Beier S, Flade I, Górska A, El-Hadidi M, Mitra S. et al.MEGAN Community Edition – interactive exploration and analysis of large-scale microbiome sequencing data.PLOS Comput Biol. 2016;12:e1004957PubMed 
PubMed Central 

Google Scholar 
Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.
Google Scholar 
Mistry J, Finn RD, Eddy SR, Bateman A, Punta M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013;41:e121–e121.CAS 
PubMed 
PubMed Central 

Google Scholar 
Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
Roux S, Adriaenssens EM, Dutilh BE, Koonin EV, Kropinski AM, Krupovic M, et al. Minimum information about an uncultivated virus genome (MIUViG). Nat Biotechnol. 2018;37:29–37.PubMed 
PubMed Central 

Google Scholar 
Germain P-L, Vitriolo A, Adamo A, Laise P, Das V, Testa G. RNAontheBENCH: computational and empirical resources for benchmarking RNAseq quantification and differential expression methods. Nucleic Acids Res. 2016;44:5054–67.CAS 
PubMed 
PubMed Central 

Google Scholar 
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. 2019.Wickham H. ggplot2: elegant graphics for data analysis. 2016. Springer-Verlag New York.Conway JR, Lex A, Gehlenborg N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics. 2017;33:2938–40.CAS 
PubMed 
PubMed Central 

Google Scholar 
Katoh K. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.CAS 
PubMed 
PubMed Central 

Google Scholar 
Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.PubMed 
PubMed Central 

Google Scholar 
Letunic I, Bork P. Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 2019;47:W256–W259.CAS 
PubMed 
PubMed Central 

Google Scholar 
Roux S, Emerson JB, Eloe-Fadrosh EA, Sullivan MB. Benchmarking viromics: an in silico evaluation of metagenome-enabled estimates of viral community composition and diversity. PeerJ. 2017;5:e3817.PubMed 
PubMed Central 

Google Scholar 
Ayllón MA, Turina M, Xie J, Nerva L, Marzano SYL, Donaire L, et al. ICTV virus taxonomy profile: botourmiaviridae. J Gen Virol. 2020;101:454–5.PubMed 
PubMed Central 

Google Scholar 
Krishnamurthy SR, Janowski AB, Zhao G, Barouch D, Wang D. Hyperexpansion of RNA bacteriophage diversity. PLOS Biol. 2016;14:e1002409.PubMed 
PubMed Central 

Google Scholar 
Hillman BI, Cai G. The family Narnaviridae. Simplest of RNA viruses. Adv Virus Res. 2013;86:149–76.
Google Scholar 
Obbard DJ, Shi M, Roberts KE, Longdon B, Dennis AB. A new lineage of segmented RNA viruses infecting animals. Virus Evol. 2020;6:61.
Google Scholar 
Xu X, Bei J, Xuan Y, Chen J, Chen D, Barker SC, et al. Full-length genome sequence of segmented RNA virus from ticks was obtained using small RNA sequencing data. BMC Genom. 2020;21:1–8.
Google Scholar 
Roossinck MJ. The good viruses: viral mutualistic symbioses. Nat Rev Microbiol. 2011;9:99–108. Nature Publishing Group.CAS 
PubMed 

Google Scholar 
Milgroom MG, Cortesi P. Biological control of chestnut blight with hypovirulence: a critical analysis. Annu Rev Phytopathol. 2004;42:311–38. Annual ReviewsCAS 
PubMed 

Google Scholar 
Zell R, Delwart E, Gorbalenya AE, Hovi T, King AMQ, Knowles NJ, et al. ICTV virus taxonomy profile: Picornaviridae. J Gen Virol. 2017;98:2421–2.CAS 
PubMed 
PubMed Central 

Google Scholar 
Valles SM, Chen Y, Firth AE, Guérin DMA, Hashimoto Y, Herrero S, et al. ICTV virus taxonomy profile: Dicistroviridae. J Gen Virol. 2017;98:355–6.CAS 
PubMed 
PubMed Central 

Google Scholar 
Barrios E. Soil biota, ecosystem services and land productivity. Ecol Econ. 2007;64:269–85.
Google Scholar 
Vainio EJ, Chiba S, Ghabrial SA, Maiss E, Roossinck M, Sabanadzovic S, et al. ICTV virus taxonomy profile: Partitiviridae. J Gen Virol. 2018;99:17–18.CAS 
PubMed 

Google Scholar 
Yong CY, Yeap SK, Omar AR, Tan WS. Advances in the study of nodavirus. PeerJ. 2017;2017:e3841.
Google Scholar 
Schmitt AP, Lamb RA. Escaping from the cell: assembly and budding of negative-strand RNA viruses. In: Kawaoka Y (ed). Biology of negative-strand RNA viruses: the power of reverse genetics. 2004. (Springer Berlin Heidelberg, Berlin, Heidelberg, pp 145–96.Käfer S, Paraskevopoulou S, Zirkel F, Wieseke N, Donath A, Petersen M, et al. Re-assessing the diversity of negative-strand RNA viruses in insects. PLoS Pathog. 2019;15:e1008224.PubMed 
PubMed Central 

Google Scholar 
Bejerman N, Debat H, Dietzgen, RG. The plant negative-sense RNA virosphere: virus discovery through new eyes. Front. Microbiol. 2020;11:588427.PubMed 
PubMed Central 

Google Scholar 
Wolf YI, Silas S, Wang Y, Wu S, Bocek M, Kazlauskas D, et al. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat Microbiol. 2020;5:1262–70.CAS 
PubMed 
PubMed Central 

Google Scholar 
Adriaenssens EM, Kramer R, van Goethem MW, Makhalanyane TP, Hogg I, Cowan DA. Environmental drivers of viral community composition in Antarctic soils identified by viromics. Microbiome. 2017;5:1–14.
Google Scholar 
Mahmoud H, Jose L. Phage and nucleocytoplasmic large viral sequences dominate coral viromes from the Arabian Gulf. Front Microbiol. 2017;8:2063.PubMed 
PubMed Central 

Google Scholar 
Koyama A, Steinweg JM, Haddix ML, Dukes JS, Wallenstein MD. Soil bacterial community responses to altered precipitation and temperature regimes in an old field grassland are mediated by plants. FEMS Microbiol Ecol. 2018;94:fix156.
Google Scholar 
Hurwitz BL, Hallam SJ, Sullivan MB. Metabolic reprogramming by viruses in the sunlit and dark ocean. Genome Biol. 2013;14:R123.PubMed 
PubMed Central 

Google Scholar  More

ColoniesColony setup occurred prior to initiation of the study, between March and May 2017, in Mississippi, USA. Using established methods, queenless colony divisions, obtained from a large commercial beekeeping operation, were equalised to an average calculated population size of ~ 7000 workers112, and housed in 10-frame Langstroth hives (Table S1). After acclimatisation for 24–48 h, they each received an imminently emerging queen cell, containing a queen from one of two stocks, added to the same worker baseline. The stocks used consisted of an Italian ‘Commercial’ stock, propagated from collaborator established breeder queens, and thus representative of the industry standard, and the Varroa-resistant ‘Pol-line’ stock54. To ensure consistency, all queens were reared in the same ‘cell builder’ colonies, based at the USDA Honey Bee Breeding, Genetics and Physiology Laboratory, in Baton Rouge, Louisiana, USA. Colonies from each stock were held in independent apiaries, 80 km apart to maintain physical isolation; and to control genetic fidelity, virgin queens were open mated to drones of the same stock via drone saturation. Fourteen days after queen emergence, colonies were inspected, and mated queens were marked with paint on the thorax, to assist with identification, with white corresponding to Commercial, and blue to Pol-line. Colonies were allowed to acclimatise for six weeks before sampling began, and those that failed to achieve mating success, or had unacceptably high [≥ 3.0 ‘mites per hundred bees’ (MPHB)] Varroa levels, were removed, normalising the average between-stock Varroa difference to  More

The MiDAS global consortium was established in 2018 to coordinate the sampling and collection of metadata from WWTPs across the globe (Supplementary Data 1). Samples were obtained in duplicates from 740 WWTPs in 425 cities, 31 countries on six continents (Fig. 1a). The majority of the WWTPs were configured with the activated sludge process (69.7%) (Fig. 1b), and these were the main focus of the subsequent analyses. Nevertheless, WWTPs based on biofilters, moving bed bioreactors (MBBR), membrane bioreactors (MBR), and granular sludge were also sampled to cover the microbial diversity in other types of WWTPs. The activated sludge plants were designed for carbon removal only (C; 22.1%), carbon removal with nitrification (C,N; 9.5%), carbon removal with nitrification and denitrification (C,N,DN; 40.9%), and carbon removal with nitrogen removal and enhanced biological phosphorus removal, EBPR (C,N,DN,P; 21.7%) (Fig. 1c). The first type represents the simplest design whereas the latter represents the most advanced process type with varying oxic and anoxic stages or compartments.Fig. 1: Sampling of WWTPs across the world.a Geographical distribution of WWTPs included in the study and their process configuration. b Distribution of plant types. MBBR moving bed bioreactor, MBR membrane bioreactor. c Distribution of process types for the activated sludge plants. C carbon removal, C,N carbon removal with nitrification, C,N,DN carbon removal with nitrification and denitrification, C,N,DN,P carbon removal with nitrogen removal and enhanced biological phosphorus removal (EBPR). The values next to the bars are the number of WWTPs in each group.Full size imageMiDAS 4: a global 16S rRNA gene catalogue and taxonomy for WWTPsMicrobial community profiling at high taxonomic resolution (genus- and species-level) using 16S rRNA gene amplicon sequencing requires a reference database with high-identity reference sequences (≥99% sequence identity) for the majority of the bacteria in the samples and a complete seven-rank taxonomy (domain to species) for all reference sequences16,20. To create such a database for bacteria in WWTPs globally, we applied synthetic long-read full-length 16S rRNA gene sequencing20,21 on samples from all WWTPs included in this study.More than 5.2 million full-length 16S rRNA gene sequences were obtained after quality filtering and primer trimming. The sequences were processed with AutoTax20 to yield 80,557 full-length 16S rRNA gene amplicon sequence variant (FL-ASVs). These reference sequences were added to our previous MiDAS 3 database16, providing a combined database (MiDAS 4) with a total of 90,164 unique, chimera-free FL-ASV reference sequences. The absence of detectable chimeric sequences is a unique feature of the database and is achieved due to the attachment of unique molecular identifiers (UMIs) to each end of the original template molecules before any PCR amplification steps21. This allows filtering of true biological sequences from chimera already in the synthetic long-read assembly20,21. The novelty of the FL-ASVs were determined based on the percent identity shared with their closest relatives in the SILVA 138 SSURef NR99 database and the threshold for each taxonomic rank proposed by Yarza et al.22. Out of all FL-ASVs, 88% had relatives above the genus-level threshold (≥94.5% identity) and 56% above the species-level threshold (≥98.7% identity) (Fig. 2 and Table 1).Fig. 2: Novel sequences and de novo taxa defined in the MiDAS 4 reference database.The phylogenetic trees are based on a multiple alignment of all MiDAS 4 reference sequences, which were first aligned against the global SILVA 138 alignment using the SINA aligner, and subsequently pruned according to the ssuref:bacteria positional variability by parsimony filter in ARB to remove hypervariable regions. The eight phyla with most FL-ASVs are highlighted in different colours. Sequence novelty was determined by the percent identity between each FL-ASV and their closest relative in the SILVA_138_SSURef_Nr99 database according to Usearch mapping and the taxonomic thresholds proposed by Yarza et al.22 shown in Table 1. Taxonomy novelty was defined based on the assignment of de novo taxa by AutoTax20.Full size imageTable 1 Novel sequences and de novo taxa observed in the MiDAS 4 reference database.Full size tableMiDAS 4 provides placeholder names for many environmental taxaAlthough only a small percentage of the reference sequences in MiDAS 4 represented new putative taxa at higher ranks (phylum, class, or order) according to the sequence identity thresholds proposed by Yarza et al.22, a large number of sequences lacked lower-rank taxonomic classifications and was assigned de novo placeholder names by AutoTax20 (Fig. 2 and Table 1). In total, de novo taxonomic names were generated by AutoTax for 26 phyla (30.6% of observed), 83 classes (37.2% of observed), 297 orders (46.8% of observed), and more than 8000 genera (86.3% of observed). Without the de novo taxonomy we would not be able to discuss these taxa across studies to unveil their potential role in wastewater treatment systems.Phylum-specific phylogenetic trees were created to determine if the FL-ASV reference sequences that were assigned to de novo phyla were actual phyla or simply artifacts related to the naive sequence identity-based assignment of de novo placeholder taxonomies (Supplementary Fig. 1a). The majority (65 FL-ASVs) created deep branches from within the Alphaproteobacteria together with 16S rRNA gene sequences from mitochondria, suggesting they represented divergent mitochondrial genes rather than true novel phyla. We also observed several FL-ASVs assigned to de novo phyla that branched from the classes Parcubacteria (3 FL-ASVs) and Microgenomatis (22 FL-ASVs) within the Patescibacteria phylum. These two classes were originally proposed as superphyla due to an unusually high rate of evolution of their 16S rRNA genes23,24. It is, therefore, likely that these de novo phyla are also artefacts due to the simple taxonomy assignment approach, which does not take different evolutionary rates into account20. Most of the class- and order-level novelty was found within the Patescibacteria, Proteobacteria, Firmicutes, Planctomycetota and Verrucomicrobiota. (Supplementary Fig. 1b). At the family- and genus-level, we also observed many de novo taxa affiliated to Bacteroidota, Bdellovibrionota and Chloroflexi.MiDAS 4 provides a common taxonomy for the fieldThe performance of the MiDAS 4 database was evaluated based on an independent amplicon dataset from the Global Water Microbiome Consortium (GWMC) project2, which covers ~1200 samples from 269 WWTPs. The raw GWMC amplicon data of the 16S rRNA gene V4 region was resolved into ASVs, and the percent identity to their best hits in MiDAS 4 and other reference databases was calculated (Fig. 3). The MiDAS 4 database had high-identity hits (≥99% identity) for 72.0 ± 9.5% (mean ± SD) of GWMC ASVs with ≥0.01% relative abundance, compared to 57.9 ± 8.5% for the SILVA 138 SSURef NR99 database, which was the best of the universal reference databases (Fig. 3). The relative abundance cutoff selects taxa that likely have a quantitative impact on the ecosystem while filtering out the rare biosphere which includes many bacteria introduced with the influent wastewaters25. Similar analyses of ASVs obtained from the samples included in this study showed, not surprisingly, even better performance with high-identity hits for 90.7 ± 7.9% of V1–V3 ASVs and 90.0 ± 6.6% of V4 ASVs with ≥0.01% relative abundance, compared to 60.6 ± 11.9% and 73.9 ± 10.3% for SILVA (Supplementary Fig. 2a). Although the sampling of WWTPs was focused towards activated sludge plants, the MiDAS 4 database also includes high-identity references for most ASVs in other plant types (granules, biofilters, etc.) (Supplementary Fig. 2b). This suggests that most taxa were shared across plant types, although often present in other relative abundances.Fig. 3: Database evaluation based on amplicon data from the Global Water Microbiome Consortium project.Raw amplicon data from the Global Water Microbiome Consortium project2 was processed to resolve ASVs of the 16S rRNA gene V4 region. The ASVs for each of the samples were filtered based on their relative abundance (only ASVs with ≥0.01% relative abundance were kept) before the analyses. The percentage of the microbial community represented by the remaining ASVs after the filtering was 88.35 ± 2.98% (mean ± SD) across samples. High-identity (≥99%) hits were determined by the stringent mapping of ASVs to each reference database. Classification of ASVs was done using the SINTAX classifier. The violin and box plots represent the distribution of percent of ASVs with high-identity hits or genus/species-level classifications for each database across n = 1165 biologically independent samples. Box plots indicate median (middle line), 25th, 75th percentile (box) and the min and max values after removing outliers based on 1.5x interquartile range (whiskers). Outliers have been removed from the box plots to ease visualisation. Different colours are used to distinguish the different databases.Full size imageUsing MiDAS 4 with the SINTAX classifier, it was possible to obtain genus-level classifications for 75.0 ± 6.9% of the GWMC ASVs with ≥0.01% relative abundance (Fig. 3). In comparison, SILVA 138 SSURef NR99, which was the best of the universal reference databases, could only classify 31.4 ± 4.2% of the ASVs to genus-level. When MiDAS 4 was used to classify amplicons from this study, we obtained genus-level classification for 92.0 ± 4.0% of V1–V3 ASVs and 84.8 ± 3.6% of V4 ASVs (Supplementary Fig. 2a). This is close to the theoretical limit set by the phylogenetic signal provided by each amplicon region analyzed20. Improved classifications were also observed for archaeal V4 ASVs (93.3 ± 10.6% for MiDAS 4 vs 69.3 ± 21.3% for SILVA), although no additional archaeal reference sequences were added to the MiDAS database in this study.MiDAS 4 was also able to assign species-level classifications to 40.8 ± 7.1% of the GWMC ASVs. In contrast, the 16S rRNA gene reference database obtained from GTDB SSU r89, which is the only universal reference database that contains a comprehensive species-level taxonomy, only classified 9.9 ± 2.0% of the ASVs (Fig. 3). For the ASVs created in this study, MiDAS 4 provided a species-level classification for 68.4 ± 6.1% of the V1–V3 and 48.5 ± 6.0% of the V4 ASVs (Supplementary Fig. 2a).Based on the large number of WWTPs sampled, their diversity, and the independent evaluation based on the GWMC dataset2, we expect that the MiDAS 4 reference database essentially covers the large majority of bacteria in WWTPs worldwide. Therefore, the MiDAS 4 taxonomy should act as a shared vocabulary for wastewater treatment microbiologists, providing opportunities for cross-study comparisons and ecological studies at high taxonomic resolution.Comparison of the V1–V3 and V4 primer sets for community profiling of WWTPsBefore investigating what factors shape the activated sludge microbiota, we compared short-read amplicon data created for all activated sludge samples belonging to the four main process types (C; C,N; C,N,DN and C,N,DN,P) collected in the Global MiDAS project using two commonly used primer sets that target the V1–V3 or V4 variable region of the 16S rRNA gene. The V1–V3 primers were chosen because the corresponding region of the 16S rRNA gene provides the highest taxonomic resolution of common short-read amplicons20,26, and these primers have previously shown great correspondence with metagenomic data and quantitative fluorescence in situ hybridisation (FISH) results for wastewater treatment systems17. The V4 region has a lower phylogenetic signal, but the primers used for amplification have better theoretical coverage of the bacterial diversity in the SILVA database20,26.The majority of genera (62%) showed less than twofold difference in relative abundances between the two primer sets, and the rest were preferentially detected with either the V1–V3 or the V4 primer (19% for both) (Fig. 4). We observed that several genera of known importance detected in high abundance by V1–V3 were hardly observed by V4, including Acidovorax, Rhodoferax, Ca. Villigracilis, Sphaerotilus and Leptothrix. Similarly, we observed genera abundant with V4 but strongly underestimated by V1–V3, such as Acinetobacter and Prosthecobacter. A complete list of differentially detected genera (Supplementary Data 2) serves as a valuable tool in combination with in silico primer evaluation for deciding which primer pair to use for targeted studies of specific taxa.Fig. 4: Comparison of relative genus abundance based on V1–V3 and V4 region 16S rRNA gene amplicon data.a Mean relative abundance was calculated based on 709 activated sludge samples. Genera present at ≥0.001% relative abundance in V1–V3 and/or V4 datasets are considered. Genera with less than twofold difference in relative abundance between the two primer sets are shown with gray circles, and those that are overrepresented by at least twofold with one of the primer sets are shown in red (V4) and blue (V1–V3). The twofold difference is an arbitrary choice; however, it relates to the uncertainty we usually encounter in amplicon data. Genus names are shown for all taxa present at a minimum of 0.1% mean relative abundance (excluding those with de novo names). b Heatmaps of the most abundant genera with more than twofold relative abundance difference between the two primer sets.Full size imageBecause the V1–V3 primers provide better classification rates at the genus- and species-level (Supplementary Fig. 2a), we primarily focused on this dataset for the following analyses. It should be noted that the V1–V3 primer set performs poorly on anammox bacteria27,28 and does not target archaea at all. To determine the importance of these groups, we estimated their relative read abundance using the V4 amplicon data. Ca. Brocadia and Ca. Anammoximicrobium were the only anammox genera detected, and the latter was never more than 0.6% abundant. Ca. Brocadia was observed in MBBR reactors and granular sludge in anammox reactors with relative read abundances reaching 29%, but it was below 0.1% relative abundance in all but two of the activated sludge samples investigated. For archaea, the relative read abundance was generally low (median = 0.18%), but for a few WWTPs high (up to 11.7%), so archaea should not be neglected in these cases.Process and environmental factors affecting the activated sludge microbiotaAlpha diversity analysis revealed that the rarefied (10,000 read per sample) richness and diversity in activated sludge plants were most strongly affected by process type, industrial load and continent (Supplementary Fig. 3 and Supplementary Note 1). The richness and diversity increased with the complexity of the treatment process, as found in other studies, reflecting the increased number of niches29. In contrast, it decreased with high industrial loads, presumably because industrial wastewater often is less complex and therefore promotes the growth of fewer specialised species7. The effect of continents is presumably caused by the necessary unbalanced sampling of WWTPs and confounded by the effects of plant types and industrial loads.Distance decay relationship (DDR) analyses were used to determine the effect of geographic distance on the microbial community similarity of activated sludge plants with the four main process types (Supplementary Fig. 4 and Supplementary Note 2). We found that distance decay was only effective within shorter geographical distances (2500 km) at the ASV-level, but higher similarities with OTUs clustered at 97% and even more at the genus-level. This suggests that many ASVs are geographically restricted and functionally redundant in the activated sludge microbiota, so different strains or species from the same genus across the world may provide similar functions.To gain a deeper understanding of the factors that shape the activated sludge microbiota, we examined the genus-level taxonomic beta-diversity using principal coordinate analysis (PCoA) and permutational multivariate analysis of variance (PERMANOVA) analyses (Fig. 5 and Supplementary Note 3). We have chosen taxonomic diversity instead of phylogenetic diversity (UniFrac) because many of the important traits are categorical (yes/no) and only conserved at lower taxonomic ranks (genus/species). The analysis was made at the genus-level due to the high classification rate achieved with MiDAS 4 and because genera were less affected by DDR compared to ASVs. We found that the overall microbial community was most strongly affected by continent and temperature in the WWTPs. However, process type, industrial load and the climate zone also had significant impacts. The percentage of total variation explained by each parameter was generally low, indicating that the global WWTPs microbiota represents a continuous distribution rather than distinct states, as observed for the human gut microbiota30.Fig. 5: Effects of process and environmental factors on the activated sludge microbial community structure. Principal coordinate analyses of Bray–Curtis and Soerensen beta-diversity for genera based on V1–V3 amplicon data. Samples are coloured based on metadata.The fraction of variation in the microbial community explained by each variable in isolation was determined by PERMANOVA (Adonis R2-values). Exact P values 0.1% relative abundance in 80% (strict core), 50% (general core) and 20% (loose core) of all activated sludge plants (Fig. 6a).Fig. 6: Identification of core and conditionally rare or abundant taxa based on V1–V3 amplicon data.a Identification of strict, general and loose core genera based on how often a given genus was observed at a relative abundance above 0.1% in WWTPs. b Identification of conditionally rare or abundant (CRAT) genera based on whether a given genus was observed at a relative abundance above 1% in at least one WWTP. The cumulative genus abundance is based on all ASVs classified at the genus-level. All core genera were removed before identification of the CRAT genera. c, d Number of genera and species, respectively, and their abundance in different process types across the global WWTPs. Values for genera and species are divided into strict core, general core, loose core, CRAT, other taxa and unclassified ASVs. The relative abundance of different groups was calculated based on the mean relative abundance of individual genera or species across samples. C carbon removal, C,N carbon removal with nitrification, C,N,DN carbon removal with nitrification and denitrification, C,N,DN,P carbon removal with nitrogen removal and enhanced biological phosphorus removal (EBPR).Full size imageIn addition to the core taxa, we also identified conditionally rare or abundant taxa (CRAT)32 (Fig. 6b). These are taxa typically present in low abundance but occasionally become prevalent, including taxa related to process disturbances, such as bacteria causing activated sludge foaming or those associated with the degradation of specific residues in industrial wastewater. CRAT have only been studied in a single WWTP treating brewery wastewater, despite their potential effect on performance32,33. CRAT are here defined as taxa which are not part of the core, but present in at least one WWTP with a relative abundance above 1%.Core taxa and CRAT were identified for both the V1–V3 and V4 amplicon data to ensure that critical taxa were not missed due to primer bias. We identified 250 core genera (15 strict, 65 general and 170 loose) and 715 CRAT genera (Supplementary Data 4). The strict core genera (Fig. 7) mainly contained genera with versatile metabolisms found in several environments, including Flavobacterium, Novosphingobium and Haliangium. The general core (Fig. 7) included many known bacteria associated with nitrification (Nitrosomonas and Nitrospira), polyphosphate accumulation (Tetrasphaera, Ca. Accumulibacter) and glycogen accumulation (Ca. Competibacter). The loose core contained well-known filamentous bacteria (Ca. Microthrix, Ca. Promineofilum, Ca. Sarcinithrix, Gordonia, Kouleothrix and Thiothrix), but also Nitrotoga, a less common nitrifier in WWTPs.Fig. 7: Percent relative abundance of strict and general core taxa across process types.The taxonomy for the core genera indicates phylum and genus. For general core species, genus names are also provided. De novo taxa in the core are highlighted in red. C carbon removal, C,N carbon removal with nitrification, C,N,DN carbon removal with nitrification and denitrification, C,N,DN,P carbon removal with nitrogen removal and enhanced biological phosphorus removal (EBPR).Full size imageBecause MiDAS 4 allowed for species-level classification, we also identified core and CRAT species based on the same criteria as for genera (Supplementary Fig. 7 and Supplementary Data 4). This revealed 113 core species (0 strict, 9 general and 104 loose). The general core species (Fig. 7) included Nitrospira defluvii and Tetrasphaera midas_s_5, a common nitrifier and PAO, respectively. Arcobacter midas_s_2255, a potential pathogen commonly abundant in the influent wastewater, was also part of the general core34. The loose core contained additional species associated with nitrification (Nitrosomonas midas_s_139 and Nitrospira nitrosa), polyphosphate accumulation (Ca. Accumulibacter phosphatis, Dechloromonas midas_s_173, Tetrasphaera midas_s_45), as well as known filamentous species (Ca. Microthrix parvicella and midas_s_2 (recently named Ca. M. subdominans35), Ca. Villigracilis midas_s_471 and midas_s_9223, Leptothrix midas_s_884). In addition to the core species, we identified 1417 CRAT species. As CRAT taxa are generally found in low abundance and the current study does not include time series or influent data, we cannot say anything conclusive about their general implications for the ecosystem. However, they may be present due to short-term mass immigration25 or specific operational conditions36 and in both cases, potentially affect the plant operation. They should therefore be considered important target for further investigations together with the core taxa.Many core taxa and CRAT can only be identified with MiDAS 4The core taxa and CRAT included a large proportion of MiDAS 4 de novo taxa. At the genus-level, 106/250 (42%) of the core genera and 500/715 (70%) of the CRAT genera had MiDAS placeholder names. At the species-level, the proportion was even higher. Here placeholder names were assigned to 101/113 (89%) of the core species and 1352/1417 (95%) CRAT species. This highlights the importance of a comprehensive taxonomy that includes the uncultured environmental taxa.The core and CRAT taxa cover the majority of the global activated sludge microbiotaAlthough the core taxa and CRAT represent a small fraction of the total diversity observed in the MiDAS 4 reference database, they accounted for the majority of the observed global activated sludge microbiota (Fig. 6c, d). Accumulated read abundance estimates ranged from 57–68% for the core genera and 11–13% for the CRAT, and combined they accounted for 68–79% of total read abundance in the WWTPs depending on process types. The core taxa represented a larger proportion of the activated sludge microbiota for the more advanced process types, which likely reflects the requirement of more versatile bacteria associated with the alternating redox conditions in these types of WWTPs. The remaining fraction, 21–32%, consisted of 6–8% unclassified genera and genera present in very low abundance, presumably with minor importance for the plant performance. The species-level core taxa and CRAT represented 11–24% and 24–33% accumulated read abundance, respectively. Combined, they accounted for almost 50% of the observed microbiota.Global diversity within important functional guildsThe general change from simple to advanced WWTPs with nutrient removal and the transition to water resource recovery facilities (WRRFs) requires increased knowledge about the bacteria responsible for the removal and recovery of nutrients, so we examined the global diversity of well-described nitrifiers, denitrifiers, PAOs and GAOs (Fig. 8). GAOs were included because they may compete with the PAOs for nutrients and thereby interfere with the biological recovery of phosphorus37. Because MiDAS 4 provided species-level resolution for a large proportion of activated sludge microbiota, we also investigated the species-level diversity within genera affiliated with the functional guilds. A complete overview of species in all genera detected in this global study is provided in the MiDAS field guide (https://www.midasfieldguide.org/guide).Fig. 8: Global diversity of genera belonging to major functional groups.The percent relative abundance represents the mean abundance for each country considering only WWTPs with the relevant process types. Countries are grouped based on continent (shifting colour).Full size imageNitrosomonas and potential comammox Nitrospira were the only abundant (≥0.1% average relative abundance) genera found among ammonia-oxidising bacteria (AOBs), whereas both Nitrospira and Nitrotoga were abundant among the nitrite oxidisers (NOBs), with Nitrospira being the most abundant across all countries (Fig. 8). Nitrobacter was not detected, and Nitrosospira was detected in only a few plants in very low abundance (≤0.01% average relative abundance). At the species-level, each genus had 2–5 abundant species (Supplementary Fig. 8). The most abundant and widespread Nitrosomonas species was midas_s_139. However, midas_s_11707 and midas_s_11733 were dominating in a few countries. For Nitrospira, the most abundant species in nearly all countries was N. defluvii. ASVs classified as the comammox N. nitrosa38,39 was also common in many countries across the world. However, because the comammox trait is not phylogenetically conserved at the 16S rRNA gene level38,39, we cannot conclude that these ASVs represent true comammox bacteria. For Nitrotoga, only two species were detected with notable abundance, midas_s_181 and midas_s_9575. Ammonia-oxidising archaea (AOAs) were not detected with MiDAS 4 due to the lack of reference sequences, and because AOAs are not targeted by the V1–V3 primer pair. However, analyses of our V4 amplicon dataset classified with the SILVA database revealed a considerable relative read abundance of AOAs in Malaysia and the Philippines, but absence or low abundance of AOAs in other countries (Supplementary Fig. 9). Other studies have occasionally found AOAs across the world, but generally in lower abundance than AOBs40,41,42. To ensure detection of AOAs with MiDAS 4, we anticipate adding external reference sequences for AOAs in a future release of the database.Denitrifying bacteria are very common in advanced activated sludge plants, but are generally poorly described. Among the known genera, Rhodoferax, Zoogloea and Thauera were most abundant (Fig. 8). Zoogloea and Thauera are well-known floc formers, sometimes causing unwanted slime formation43. Rhodoferax was the most common denitrifier in Europe, whereas Thauera dominated in Asia. Many denitrifiers could not be classified at the species-level (Supplementary Fig. 10), likely due to highly conserved 16S rRNA genes. An exception was Zoogloea, where midas_s_1080 and Z. caeni and were the most abundant species worldwide.EBPR is performed by PAOs, with three genera recognised as important in full-scale WWTPs: Tetrasphaera, Dechloromonas and Ca. Accumulibacter13. According to relative read abundance, all three were found in EBPR plants globally, with Tetrasphaera as the most prevalent (Fig. 8). Dechloromonas was also abundant in nitrifying and denitrifying plants without EBPR, indicating a more diverse ecology. Four recognised GAOs were found globally: Ca. Competibacter, Defluviicoccus, Propionivibrio and Micropruina, with Ca. Competibacter being the most abundant (Fig. 8). Only a few species (2–6 species) in each genus were dominant across the world for both PAOs (Supplementary Fig. 11) and GAOs (Supplementary Fig. 12), except for Ca. Competibacter, which covered ~20 abundant but country-specific species. Among PAOs, the abundant species were Tetrasphaera midas_s_5, Dechloromonas midas_s_173, (recently named D. phosphorivorans) Ca. Accumulibacter midas_s_315, Ca. A. phosphatis and Ca. A. aalborgensis. Interestingly, some of the most abundant PAOs and GAOs were also abundant in the simple process design with C-removal, indicating more versatile metabolisms.Global diversity of filamentous bacteriaFilamentous bacteria are essential for creating strong activated sludge flocs. However, in large numbers, they can also lead to loose flocs and poor settling properties. This is known as bulking, a major operational problem in many WWTPs. Many can also form foam on top of process tanks due to hydrophobic surfaces. Presently, approximately 20 genera are known to contain filamentous species44, and among those, the most abundant are Ca. Microthrix, Leptothrix, Ca. Villigracilis, Trichococcus and Sphaerotilus (Fig. 9). They are all well-known from studies on mitigation of poor settling properties in WWTPs. Interestingly, Leptothrix, Sphaerotilus and Ca. Villigracilis belong to the genera where abundance-estimation depended strongly on primers, with V4 underestimating their abundance (Fig. 3). Ca. Microthrix and Leptothrix were strongly associated with continents, most common in Europe and less in Asia and North America (Fig. 9).Fig. 9: Global diversity of known filamentous organisms.The percent relative abundance represents the mean abundance for each country across all process types. Countries are grouped based on the continent (shifting colour).Full size imageMany of the filamentous bacteria were linked to specific process types (Supplementary Fig. 13), e.g. Ca. Microthrix were not observed in WWTPs with carbon removal only, and Ca. Amarolinea were only abundant in plants with nutrient removal. The number of abundant species within the genera were generally low, with one species in Trichococcus, two in Ca. Microthrix and approximately five in Leptothrix and Ca. Villigracilis (Supplementary Fig. 14). Only five abundant species were observed for Sphaerotilus. However, a substantial fraction of unclassified ASVs was also observed, demonstrating that certain species within this genus are poorly resolved based on the 16S rRNA gene. Ca. Promineofilum was also poorly resolved at the species-level (Supplementary Fig. 15).Conclusion and perspectivesWe present a worldwide collaborative effort to produce MiDAS 4, an ASV-resolved full-length 16S rRNA gene reference database, which covers more than 31,000 species and enables genus- to species-level resolution in microbial community profiling studies. MiDAS 4 covers the vast majority of WWTP bacteria globally and provides a strongly needed common taxonomy for the field, which provides the foundation for comprehensive linking of microbial taxa in the ecosystem with their functional traits. Presently, hundreds of studies are undertaken to combine engineering and microbial aspects of full-scale WWTPs. However, most ASVs or OTUs in these studies are classified at poor taxonomic resolution (family-level or above) due to the use of incomplete universal reference databases. Because many important functional traits are only conserved at high taxonomic resolution (genus- or species-level), this strongly hampers our ability to transfer taxa-specific knowledge from one study to another. This will change with MiDAS 4, and we expect that reprocessing of data from earlier studies may reveal new perspectives into wastewater treatment microbiology. Our online Global MiDAS Field Guide presents the data generated in this study and summarises present knowledge about all taxa. We encourage researchers within the field to contribute new knowledge to MiDAS using the contact link in the MiDAS website (https://www.midasfieldguide.org/guide/contact).The global microbiota of activated sludge plants has been predicted to harbour a massive diversity with up to one billion species2. However, most of these occur at very low abundance and are of little importance for the treatment process. By focusing only on the abundant taxa, we can see that this number is much smaller, i.e., ~1000 genera and 1500 species. We consider these taxa functionally the most important globally, representing a “most wanted list” for future studies. Some taxa are abundant in most WWTPs (core taxa), and others are occasionally abundant in fewer plants (CRAT). The CRAT have received little attention in the field of wastewater treatment, but they can be of profound importance for WWTP performance. Both groups have a high fraction of poorly characterised species. The high taxonomic resolution provided by MiDAS 4 enables us to identify samples where these important core taxa occur in high abundance. This provides an ideal starting point for obtaining high-quality metagenome-assembled genomes (MAGs), isolation of pure cultures, in addition to targeted culture-independent studies to uncover their physiological and ecological roles.Among the known functional guilds, such as nitrifiers or polyphosphate-accumulating organisms, the same genera were found worldwide, with only a few abundant species in each genus. There were differences in the community structure, and the abundance of dominant species was mainly shaped by process type, temperature, and in some cases, continent. This discovery sends an important message to the field: relatively few species are abundant worldwide, so research or operational results can reliably be transferred from one geographical region to another, stimulating the transition from WWTPs to more sustainable WRRFs.The relatively low number of uncharacterised abundant species also shows that it is within our reach to describe them all in terms of identity, physiology, ecology and dynamics, providing the necessary knowledge for informed process optimisation and management. The number of poorly described genera (i.e. those with only a MiDAS placeholder genus name) was 88 among the 250 core genera (35%) and more than 89% at the species-level, so there is still some work to do to link their identities and function. An important step in this direction is the visualisation of the populations. With the comprehensive set of FL-ASVs, it is possible to design highly specific FISH probes, and to critically evaluate the old probes. In the Danish WWTPs, we have successfully done this for groups in the Acidobacteriota42 based on the MiDAS 3 database18. Our recent retrieval of more than 1000 high-quality MAGs from Danish WWTPs with advanced process design is also an important step to link identity to function43. The HQ-MAGs can be linked directly to MiDAS 4 as they contain complete 16S rRNA genes. They cover 62% (156/250) of the core genera and 61% (69/113) of the core species identified in this study. These MAGs may also form the basis for further studies to link identity and function, e.g. by applying metatranscriptomics44 and other in situ techniques such as FISH combined with Raman45,46, guided by the “most wanted” list provided in this study. We expect that MiDAS 4 will have significant implications for future microbial ecology studies in wastewater treatment systems. More

Zhang, N. et al. Members of the Fusarium solani species complex that cause infections in both humans and plants are common in the environment. J. Clin. Microbiol. 44, 2186–2190 (2006).CAS 
PubMed 
PubMed Central 

Google Scholar 
O’Donnell, K. et al. Molecular Phylogenetic Diversity, Multilocus Haplotype Nomenclature, and In Vitro antifungal resistance within the Fusarium solani species complex. J. Clin. Microbiol. 46, 2477–2490 (2008).PubMed 
PubMed Central 

Google Scholar 
Schroers, H. J. et al. Epitypification of Fusisporium (Fusarium) solani and its assignment to a common phylogenetic species in the Fusarium solani species complex. Mycologia 108, 806–819 (2016).CAS 
PubMed 

Google Scholar 
O’Donnell, K. Molecular phylogeny of the Nectria haematococca-Fusarium solani species complex. Mycologia 92, 919–938 (2000).
Google Scholar 
Gleason, F., Allerstorfer, M. & Lilje, O. Newly emerging diseases of marine turtles, especially sea turtle egg fusariosis (SEFT), caused by species in the Fusarium solani complex (FSSC). Mycology 11, 184–194 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Fernando, N. et al. Fatal Fusarium solani species complex infections in elasmobranchs: the first case report for black spotted stingray (Taeniura melanopsila) and a literature review. Mycoses 58, 422–431 (2015).PubMed 

Google Scholar 
Sarmiento-Ramírez, J. M. et al. Global distribution of two fungal pathogens threatening endangered Sea Turtles. PLoS ONE 9, e85853 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Mayayo, E., Pujol, I. & Guarro, J. Experimental pathogenicity of four opportunist Fusarium species in a murine model. J. Med. Microbiol. 48, 363–366 (1999).CAS 
PubMed 

Google Scholar 
Muhvich, A. G., Reimschuessel, R., Lipsky, M. M. & Bennett, R. O. Fusarium solani isolated from newborn bonnethead sharks, Sphyrna tiburo (L.). J. Fish Dis. 12, 57–62 (1989).
Google Scholar 
Crow, G. L., Brock, J. A. & Kaiser, S. Fusarium solani fungal infection of the lateral line canal system in captive scalloped hammerhead sharks (Sphyrna lewini) in Hawaii. J. Wildl. Dis. 31, 562–565 (1995).CAS 
PubMed 

Google Scholar 
Cabañes, F. J. et al. Cutaneous hyalohyphomycosis caused by Fusarium solani in a loggerhead sea turtle (Caretta caretta L.). J. Clin. Microbiol. 35, 3343–3345 (1997).PubMed 
PubMed Central 

Google Scholar 
Cafarchia, C. et al. Fusarium spp. in Loggerhead Sea Turtles (Caretta caretta): From Colonization to Infection. Vet. Pathol. 57, 139–146 (2019).PubMed 

Google Scholar 
Garcia-Hartmann, M., Hennequin, C., Catteau, S., Béatini, C. & Blanc, V. Cas groupés d’infection à Fusarium solani chez de jeunes tortues marines Caretta caretta nées en captivité. J. Mycol. Med. 28, 113–118 (2017).
Google Scholar 
Orós, J., Delgado, C., Fernández, L. & Jensen, H. E. Pulmonary hyalohyphomycosis caused by Fusarium spp in a Kemp’s ridley sea turtle (Lepidochelys kempi): An immunohistochemical study. N. Z. Vet. J. 52, 150–152 (2004).PubMed 

Google Scholar 
Candan, A. Y., Katılmış, Y. & Ergin, Ç. First report of Fusarium species occurrence in loggerhead sea turtle (Caretta caretta) nests and hatchling success in Iztuzu Beach, Turkey. Biologia (Bratisl). https://doi.org/10.2478/s11756-020-00553-4 (2020).Article 

Google Scholar 
Sarmiento-Ramirez, J. M., van der Voort, M., Raaijmakers, J. M. & Diéguez-Uribeondo, J. Unravelling the Microbiome of eggs of the endangered Sea Turtle Eretmochelys imbricata identifies bacteria with activity against the emerging pathogen Fusarium falciforme. PLoS ONE 9, e95206 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Sarmiento-Ramírez, J. M. et al. Fusarium solani is responsible for mass mortalities in nests of loggerhead sea turtle, Caretta caretta, in Boavista, Cape Verde. FEMS Microbiol. Lett. 312, 192–200 (2010).PubMed 

Google Scholar 
Sarmiento-Ramirez, J. M., Sim, J., Van West, P. & Dieguez-Uribeondo, J. Isolation of fungal pathogens from eggs of the endangered sea turtle species Chelonia mydas in Ascension Island. J. Mar. Biol. Assoc. United Kingdom 97, 661–667 (2017).CAS 

Google Scholar 
Hoh, D., Lin, Y., Liu, W., Sidique, S. & Tsai, I. Nest microbiota and pathogen abundance in sea turtle hatcheries. Fungal Ecol. 47, 100964 (2020).
Google Scholar 
Güçlü, Ö., Bıyık, H. & Şahiner, A. Mycoflora identified from loggerhead turtle (Caretta caretta) egg shells and nest sand at Fethiye beach, Turkey. Afr. J. Microbiol. Res. 4, 408–413 (2010).
Google Scholar 
Gambino, D. et al. First data on microflora of loggerhead sea turtle (Caretta caretta) nests from the coastlines of Sicily. Biol. Open 9, bio045252 (2020).PubMed 
PubMed Central 

Google Scholar 
Bailey, J. B., Lamb, M., Walker, M., Weed, C. & Craven, K. S. Detection of potential fungal pathogens Fusarium falciforme and F. keratoplasticum in unhatched loggerhead turtle eggs using a molecular approach. Endanger. Species Res. 36, 111–119 (2018).
Google Scholar 
Summerbell, R. C. & Schroers, H.-J. Analysis of Phylogenetic Relationship of Cylindrocarpon lichenicola and Acremonium falciforme to the Fusarium solani Species Complex and a Review of similarities in the spectrum of opportunistic infections caused by these fungi. J. Clin. Microbiol. 40, 2866–2875 (2002).CAS 
PubMed 
PubMed Central 

Google Scholar 
Nel, R., Punt, A. E. & Hughes, G. R. Are coastal protected areas always effective in achieving population recovery for nesting sea turtles?. PLoS ONE 8, e63525 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Branch, G. & Branch, M. Living Shores. (Pippa Parker, 2018).Fuller, M. S., Fowles, B. E. & Mclaughlin, D. J. Isolation and pure culture study of marine phycomycetes. Mycologia 56, 745–756 (1964).
Google Scholar 
Greeff, M. R., Christison, K. W. & Macey, B. M. Development and preliminary evaluation of a real-time PCR assay for Halioticida noduliformans in abalone tissues. Dis. Aquat. Organ. 99, 103–117 (2012).CAS 
PubMed 

Google Scholar 
Sandoval-Denis, M., Lombard, L. & Crous, P. W. Back to the roots: a reappraisal of Neocosmospora. Persoonia Mol. Phylogeny Evol. Fungi 43, 90–185 (2019).CAS 

Google Scholar 
O’Donnell, K., Cigelnik, E. & Nirenberg, H. I. Molecular systematics and phylogeography of the Gibberella fujikuroi species complex. Mycologia 90, 465–493 (1998).
Google Scholar 
Geiser, D. M. et al. FUSARIUM-ID v. 1. 0: A DNA sequence database for identifying Fusarium. Eur. J. Plant Pathol. 110, 473–479 (2004).ADS 
CAS 

Google Scholar 
O’Donnell, K. et al. Phylogenetic diversity of insecticolous fusaria inferred from multilocus DNA sequence data and their molecular identification via FUSARIUM-ID and FUSARIUM MLST. Mycologia 104, 427–445 (2012).PubMed 

Google Scholar 
Chehri, K., Salleh, B. & Zakaria, L. Morphological and phylogenetic analysis of Fusarium solani species complex in Malaysia. Microb. Ecol. 69, 457–471 (2015).PubMed 

Google Scholar 
Lanfear, R., Frandsen, P., Wright, A., Senfeld, T. & Calcott, B. PartionFinder 2: new methods for selecting partioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. https://doi.org/10.1093/molbev/msw260 (2016).Article 

Google Scholar 
Ronquist, F. et al. Efficient Bayesian phylogenetic inference and model selection across a large model space. Syst. Biol. 61, 539–542 (2012).PubMed 
PubMed Central 

Google Scholar 
Leslie, J. F. & Summerell, B. A. The Fusarium Laboratory manual (Blackwell Publishing, Hoboken, 2006).
Google Scholar 
Fisher, N. L., Burgess, L. W., Toussoun, T. A. & Nelson, P. E. Carnation leaves as a substrate and for preserving cultures of Fusarium species. Phytopathology 72, 151 (1982).
Google Scholar 
Smyth, C. W. et al. Unraveling the ecology and epidemiology of an emerging fungal disease, sea turtle egg fusariosis (STEF). PLOS Pathog. 15, e1007682 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Rachowicz, L. J. et al. The novel and endemic pathogen hypotheses: Competing explanations for the origin of emerging infectious diseases of wildlife. Conserv. Biol. 19, 1441–1448 (2005).
Google Scholar 
Lombard, L., Sandoval-Denis, M., Cai, L. & Crous, P. W. Changing the game: resolving systematic issues in key Fusarium species complexes. Persoonia Mol. Phylogeny Evol. Fungi 43, i–ii (2019).CAS 

Google Scholar 
Short, D. P. G., Donnell, K. O., Zhang, N., Juba, J. H. & Geiser, D. M. Widespread occurrence of diverse human pathogenic types of the fungus Fusarium detected in plumbing drains. J. Clin. Microbiol. 49, 4264–4272 (2011).PubMed 
PubMed Central 

Google Scholar 
White, T. J., Burns, T., Lee, S. & Taylor, J. Amplification and direct identification of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: a guide to methods and applications (eds Innis, M. A. et al.) 315–322 (Academic Press, San Diego, 1990).
Google Scholar 
Sekimoto, S., Hatai, K. & Honda, D. Molecular phylogeny of an unidentified Haliphthoros-like marine oomycete and Haliphthoros milfordensis inferred from nuclear-encoded small- and large-subunit rRNA genes and mitochondrial-encoded cox2 gene. Mycoscience 48, 212–221 (2007).CAS 

Google Scholar 
Petersen, A. B. & Rosendahl, S. Ø. Phylogeny of the Peronosporomycetes (Oomycota) based on partial sequences of the large ribosomal subunit (LSU rDNA). Mycol. Res. 104, 1295–1303 (2000).CAS 

Google Scholar 
O’Donnell, K. et al. Phylogenetic diversity and microsphere array-based genotyping of human pathogenic fusaria, including isolates from the multistate contact lens-associated U.S. keratitis outbreaks of 2005 and 2006. J. Clin. Microbiol. 45, 2235–2248 (2007).PubMed 
PubMed Central 

Google Scholar 
Migheli, Q. et al. Molecular Phylogenetic diversity of dermatologic and other human pathogenic fusarial isolates from hospitals in Northern and Central Italy. J. Clin. Microbiol. 48, 1076–1084 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar  More

Hyman, L. H. The Invertebrates: Protozoa Through Ctenophora Vol. 1 (McGraw-Hill, 1940).Taylor, M. W., Radax, R., Steger, D. & Wagner, M. Sponge-associated microorganisms: evolution, ecology, and biotechnological potential. Microbiol. Mol. Biol. Rev. 71, 295–347 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
Giles, E. C. et al. Bacterial community profiles in low microbial abundance sponges. FEMS Microbiol. Ecol. 83, 232–241 (2013).CAS 
PubMed 

Google Scholar 
Gloeckner, V. et al. The HMA–LMA dichotomy revisited: an electron microscopical survey of 56 sponge species. Biol. Bull. 227, 78–88 (2014).PubMed 

Google Scholar 
Moitinho-Silva, L. et al. Predicting the HMA–LMA status in marine sponges by machine learning. Front. Microbiol. 8, 752 (2017).PubMed 
PubMed Central 

Google Scholar 
Cárdenas, C. A. et al. High similarity in the microbiota of cold-water sponges of the genus Mycale from two different geographical areas. PeerJ 6, e4935 (2018).PubMed 
PubMed Central 

Google Scholar 
Webster, N. S. & Taylor, M. W. Marine sponges and their microbial symbionts: love and other relationships. Environ. Microbiol. 14, 335–346 (2012).CAS 
PubMed 

Google Scholar 
Freeman, C. J. et al. Microbial symbionts and ecological divergence of Caribbean sponges: a new perspective on an ancient association. ISME J. 14, 1571–1583 (2020).PubMed 
PubMed Central 

Google Scholar 
Bell, J. J. et al. Climate change alterations to ecosystem dominance: how might sponge-dominated reefs function? Ecology 99, 1920–1931 (2018).PubMed 

Google Scholar 
Gardner, T. A., Côté, I. M., Gill, J. A., Grant, A. & Watkinson, A. R. Long-term region-wide declines in Caribbean corals. Science 301, 958–960 (2003).CAS 
PubMed 

Google Scholar 
Lesser, M. P. Benthic–pelagic coupling on coral reefs: feeding and growth of Caribbean sponges. J. Exp. Mar. Biol. Ecol. 328, 277–288 (2006).
Google Scholar 
de Goeij, J. M., Lesser, M. P. & Pawlik, J. R. in Climate Change, Ocean Acidification and Sponges (eds Carballo, J. L. & Bell, J. J.) 373–410 (Springer, 2017); https://doi.org/10.1007/978-3-319-59008-0_8Pita, L., Rix, L., Slaby, B. M., Franke, A. & Hentschel, U. The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome 6, 46 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Slaby, B. M., Hackl, T., Horn, H., Bayer, K. & Hentschel, U. Metagenomic binning of a marine sponge microbiome reveals unity in defense but metabolic specialization. ISME J. 11, 2465–2478 (2017).PubMed 
PubMed Central 

Google Scholar 
Moitinho-Silva, L. et al. Revealing microbial functional activities in the Red Sea sponge Stylissa carteri by metatranscriptomics. Environ. Microbiol. 16, 3683–3698 (2014).CAS 
PubMed 

Google Scholar 
Weisz, J. B., Lindquist, N. & Martens, C. S. Do associated microbial abundances impact marine demosponge pumping rates and tissue densities? Oecologia 155, 367–376 (2008).PubMed 

Google Scholar 
Poppell, E. et al. Sponge heterotrophic capacity and bacterial community structure in high- and low-microbial abundance sponges. Mar. Ecol. 35, 414–424 (2014).
Google Scholar 
McFall-Ngai, M. J. et al. Animals in a bacterial world, a new imperative for the life sciences. Proc. Natl Acad. Sci. USA 110, 3229–3236 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Douglas, A. E. Symbiosis as a general principle in eukaryotic evolution. Cold Spring Harb. Perspect. Biol. 6, a016113 (2014).PubMed 
PubMed Central 

Google Scholar 
Moran, N. A. & Sloan, D. B. The hologenome concept: helpful or hollow? PLoS Biol. 13, e1002311 (2015).PubMed 
PubMed Central 

Google Scholar 
Brooks, A. W., Kohl, K. D., Brucker, R. M., van Opstal, E. J. & Bordenstein, S. R. Phylosymbiosis: relationships and functional effects of microbial communities across host evolutionary history. PLoS Biol. 14, e2000225–e2000229 (2016).PubMed 
PubMed Central 

Google Scholar 
O’Brien, P. A. et al. Diverse coral reef invertebrates exhibit patterns of phylosymbiosis. ISME J. 14, 2211–2222 (2020).PubMed 
PubMed Central 

Google Scholar 
Houwenhuyse, S., Stoks, R., Mukherjee, S. & Decaestecker, E. Locally adapted gut microbiomes mediate host stress tolerance. ISME J. 15, 2401–2414 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moeller, A. H. et al. Experimental evidence for adaptation to species-specific gut microbiota in house mice. mSphere 4, e00387-19 (2019).PubMed 
PubMed Central 

Google Scholar 
van Opstal, E. J. & Bordenstein, S. R. Phylosymbiosis impacts adaptive traits in Nasonia wasps. mBio https://doi.org/10.1128/mBio.00887-19 (2019).Lim, S. J. & Bordenstein, S. R. An introduction to phylosymbiosis. Proc. R. Soc. B https://doi.org/10.1098/rspb.2019.2900 (2020).Pollock, F. J. et al. Coral-associated bacteria demonstrate phylosymbiosis and cophylogeny. Nat. Commun. https://doi.org/10.1038/s41467-018-07275-x (2018).Douglas, A. E. & Werren, J. H. Holes in the hologenome: why host–microbe symbioses are not holobionts. mBio 7, e02099 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hadfield, J. D., Krasnov, B. R., Poulin, R. & Nakagawa, S. A tale of two phylogenies: comparative analyses of ecological interactions. Am. Nat. 183, 174–187 (2014).PubMed 

Google Scholar 
Hill, M. S. et al. Reconstruction of family-level phylogenetic relationships within Demospongiae (Porifera) using nuclear encoded housekeeping genes. PLoS ONE 8, e50437 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Redmond, N. E. et al. Phylogeny and systematics of Demospongiae in light of new small-subunit ribosomal DNA (18S) sequences. Int. Comp. Biol. 53, 388–415 (2013).CAS 

Google Scholar 
Worheide, G. et al. in Advances in Marine Biology: Advances in Sponge Science Vol. 61 (eds Becerro, M. A. et al.) 1–78 (Elsevier, 2012).Schuster, A. et al. Divergence times in demosponges (Porifera): first insights from new mitogenomes and the inclusion of fossils in a birth–death clock model. BMC Evol. Biol. 18, 114 (2018).PubMed 
PubMed Central 

Google Scholar 
Stanley, G. D. & Fautin, D. G. Paleontology and evolution. Orig. Mod. Corals Sci. 291, 1913–1914 (2001).CAS 

Google Scholar 
Brinkmann, C. M., Marker, A. & Kurtböke, D. I. An overview on marine sponge-symbiotic bacteria as unexhausted sources for natural product discovery. Diversity 9, 40 (2017).
Google Scholar 
Rust, M. et al. A multiproducer microbiome generates chemical diversity in the marine sponge Mycale hentscheli. Proc. Natl Acad. Sci. USA 117, 9508–9518 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Faulkner, D. J., Harper, M. K., Haygood, M. G., Salomon, C. E. & Schmidt, E. W. in Drugs from the Sea (ed. Fusetani, N.) 107–119 (Karger, 2000).Loh, T.-L. & Pawlik, J. R. Chemical defenses and resource trade-offs structure sponge communities on Caribbean coral reefs. Proc. Natl Acad. Sci. USA 111, 4151–4156 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Pagel, M. Detecting correlated evolution on phylogenies—a general method for the comparative analysis of discrete characters. Proc. R. Soc. Lond. B 255, 37–45 (1994).
Google Scholar 
Easson, C. G. & Thacker, R. W. Phylogenetic signal in the community structure of host-specific microbiomes of tropical marine sponges. Front. Microbiol. 5, 532 (2014).PubMed 
PubMed Central 

Google Scholar 
Thomas, T. et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat. Commun. 7, 11870 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Schöttner, S. et al. Relationships between host phylogeny, host type and bacterial community diversity in cold-water coral reef sponges. PLoS ONE 8, e55505 (2013).PubMed 
PubMed Central 

Google Scholar 
Robinson, D. R. & Foulds, L. R. Comparison of phylogenetic trees. Math. Biosci. 53, 131–147 (1981).
Google Scholar 
Gloor, G. B., Macklaim, J. M., Pawlowsky-Glahn, V. & Egozcue, J. J. Microbiome datasets are compositional: and this is not optional. Front. Microbiol. 8, 2224 (2017).PubMed 
PubMed Central 

Google Scholar 
Apprill, A. The role of symbioses in the adaptation and stress responses of marine organisms. Annu. Rev. Mar. Sci. 12, 291–314 (2020).
Google Scholar 
Lesser, M. P., Slattery, M. & Mobley, C. Biodiversity and functional ecology of mesophotic coral reefs. Annu. Rev. Ecol. Evol. Syst. 49, 49–71 (2018).
Google Scholar 
Lipps, J. H. & Stanley, G. D. in Coral Reefs at the Crossroads (eds Hubbard, D. K. et al.) 175–196 (Springer, 2016); https://doi.org/10.1007/978-94-017-7567-0_8Macartney, K. J., Slattery, M. & Lesser, M. P. Trophic ecology of Caribbean sponges in the mesophotic zone. Limnol. Oceanogr. 66, 1113–1124 (2021).CAS 

Google Scholar 
McMurray, S. E., Stubler, A. D., Erwin, P. M., Finelli, C. M. & Pawlik, J. R. A test of the sponge-loop hypothesis for emergent Caribbean reef sponges. Mar. Ecol. Prog. Ser. 588, 1–14 (2018).CAS 

Google Scholar 
Olinger, L. K., Strangman, W. K., McMurray, S. E. & Pawlik, J. R. Sponges with microbial symbionts transform dissolved organic matter and take up organohalides. Front. Mar. Sci. 8, 665789 (2021).
Google Scholar 
Haas, A. F. et al. Effects of coral reef benthic primary producers on dissolved organic carbon and microbial activity. PLoS ONE 6, e27973 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sánchez-Baracaldo, P. Origin of marine planktonic cyanobacteria. Sci. Rep. 5, 17418 (2015).PubMed 
PubMed Central 

Google Scholar 
Sanchez-Bracaldo, P., Ridgwell, A. & Raven, J. A. A neoproterozoic transition in the marine nitrogen cycle. Curr. Biol. 24, 652–657 (2014).
Google Scholar 
Falkowski, P. G. et al. The evolution of modern eukaryotic phytoplankton. Science 305, 354–360 (2004).CAS 
PubMed 

Google Scholar 
Wang, D. et al. Coupling of ocean redox and animal evolution during the Ediacaran–Cambrian transition. Nat. Commun. 9, 2575 (2018).PubMed 
PubMed Central 

Google Scholar 
Bellwood, D. R., Goatley, C. H. R. & Bellwood, O. The evolution of fishes and corals on reefs: form, function and interdependence. Biol. Rev. 92, 878–901 (2017).PubMed 

Google Scholar 
Ehrlich, P. R. & Raven, P. H. Butterflies and plants: a study in coevolution. Evolution 18, 586–608 (1964).
Google Scholar 
Després, L., David, J.-P. & Gallet, C. The evolutionary ecology of insect resistance to plant chemicals. Trends Ecol. Evol. 22, 298–307 (2007).PubMed 

Google Scholar 
Richardson, K. L., Gold-Bouchot, G. & Schlenk, D. The characterization of cytosolic glutathione transferase from four species of sea turtles: loggerhead (Caretta caretta), green (Chelonia mydas), olive ridley (Lepidochelys olivacea), and hawksbill (Eretmochelys imbricata). Comp. Biochem. Physiol. C 150, 279–284 (2009).
Google Scholar 
Bayer, K., Jahn, M. T., Slaby, B. M., Moitinho-Silva, L. & Hentschel, U. Marine sponges as Chloroflexi hot spots: genomic insights and high-resolution visualization of an abundant and diverse symbiotic clade. mSystems 3, e00150-18 (2018).PubMed 
PubMed Central 

Google Scholar 
Sachs, J. L., Skophammer, R. G., Bansal, N. & Stajich, J. E. Evolutionary origins and diversification of proteobacterial mutualists. Proc. R Soc. B https://doi.org/10.1098/rspb.2013.2146 (2014).Sachs, J. L., Skophammer, R. G. & Regus, J. U. Evolutionary transitions in bacterial symbiosis. Proc. Natl Acad. Sci. USA 108, 10800–10807 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Seutin, G., White, B. N. & Boag, P. T. Preservation of avian blood and tissue samples for DNA analyses. Can. J. Zool. https://doi.org/10.1139/z91-013 (2011).Sunagawa, S. et al. Generation and analysis of transcriptomic resources for a model system on the rise: the sea anemone Aiptasia pallida and its dinoflagellate endosymbiont. BMC Genomics 10, 258 (2009).PubMed 
PubMed Central 

Google Scholar 
Song, L. & Florea, L. Rcorrector: efficient and accurate error correction for Illumina RNA-seq reads. GigaScience 4, 48 (2015).PubMed 
PubMed Central 

Google Scholar 
Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 29, 644–652 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Chevreux, B., Wetter, T. & Suhai, S. Genome sequence assembly using trace signals and additional sequence information. Comput. Sci. Biol. 99, 45–56 (1999).
Google Scholar 
Li, W. & Godzik, A. CD-HIT: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658–1659 (2006).CAS 

Google Scholar 
Francis, W. R. et al. The genome of the contractile demosponge Tethya wilhelma and the evolution of metazoan neural signalling pathways. Preprint at bioRxiv https://doi.org/10.1101/120998 (2017).Altschul, S. F. A protein alignment scoring system sensitive at all evolutionary distances. J. Mol. Evol. 36, 290–300 (1993).CAS 
PubMed 

Google Scholar 
Srivastava, M. et al. The Amphimedon queenslandica genome and the evolution of animal complexity. Nature 466, 720–726 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Simion, P. et al. A large and consistent phylogenomic dataset supports sponges as the sister group to all other animals. Curr. Biol. https://doi.org/10.1016/j.cub.2017.02.031 (2017).Katoh, K., Misawa, K., Kuma, K.-I. & Miyata, T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066 (2002).CAS 
PubMed 
PubMed Central 

Google Scholar 
Castresana, J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540–552 (2000).CAS 

Google Scholar 
Kalyaanamoorthy, S. et al. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688–2690 (2006).CAS 

Google Scholar 
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 

Google Scholar 
Dohrmann, M. & Wörheide, G. Dating early animal evolution using phylogenomic data. Sci. Rep. 7, 3599 (2017).PubMed 
PubMed Central 

Google Scholar 
Smith, S. A. & O’Meara, B. C. treePL: divergence time estimation using penalized likelihood for large phylogenies. Bioinformatics 28, 2689–2690 (2012).CAS 
PubMed 

Google Scholar 
Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).CAS 
PubMed 

Google Scholar 
Apprill, A., McNally, S., Parsons, R. & Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 75, 129–137 (2015).
Google Scholar 
Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 

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

Google Scholar 
McMurdie, P. J. & Holmes, S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Oksanen, J. et al. vegan: Community Ecology Package. R package version 2.5-5 (2019).Lahti, L. et al. Tools for Microbiome Analysis in R. Microbiome package version 1.17.2 https://github.com/microbiome/microbiome (2017).Harmon, L. J., Weir, J. T., Brock, C. D., Glor, R. E. & Challenger, W. GEIGER: investigating evolutionary radiations. Bioinformatics 24, 129–131 (2008).CAS 
PubMed 

Google Scholar 
Schliep, K. P. phangorn: phylogenetic analysis in R. Bioinformatics 27, 592–593 (2011).CAS 

Google Scholar 
Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223 (2012).
Google Scholar 
Kuhn, M. Building predictive models in R using the caret package. J. Stat. Softw. 28, 1–26 (2008).
Google Scholar 
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Westbrook, A. et al. PALADIN: protein alignment for functional profiling whole metagenome shotgun data. Bioinformatics 33, 1473–1478 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Waddell, B. & Pawlik, J. R. Defenses of Caribbean sponges against invertebrate predators. I. Assays with hermit crabs. Mar. Ecol. Prog. Ser. 195, 125–132 (2000).
Google Scholar 
Paradis, E., Claude, J. & Strimmer, K. APE: analyses of phylogenetics and evolution in R language. FEMS Microbiol. Ecol. 20, 289–290 (2004).CAS 

Google Scholar 
Hadfield, J. D. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33, 1–22 (2010).
Google Scholar 
Nakagawa, S., Johnson, P. C. D. & Schielzeth, H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R. Soc. Interface 14, 20170213 (2017).PubMed 
PubMed Central 

Google Scholar  More