Ecology

High diet overlap is assumed to cause competition between the three dominant pelagic planktivorous mesopredators in the Baltic Sea, sprat, herring, and stickleback11,24,25. Despite this assumption, stickleback populations have increased dramatically over the past decades, which raises the question of whether and how resources are partitioned26. While previous studies of fish diet overlap have mainly relied on microscopic identification of gut content, we implemented a DNA metabarcoding approach targeting two different gene regions, the 18S rRNA gene (18S) and the mitochondrial cytochrome c oxidase I gene (COI) to reveal the taxonomic diversity of prey, and a qPCR step to quantify rotifers that are at times abundant in the Baltic Sea. Our study highlights consistency between methods, with DNA metabarcoding resolving the plankton-fish link at the highest taxonomic resolution. Our results suggest a unique niche of stickleback that may enable high population growth in the open water, despite high competition between mesopredators, although this finding needs to be confirmed at larger scale. More than half of the DNA found in herring and sprat stomach contents was assigned to Pseudocalanus, supporting previous observations of high diet overlap between the two clupeids11,12. On the other hand, the diet of stickleback differed substantially from the two clupeids, with rotifers appearing as main prey DNA in spring. The high rotifer biomass in the environment and lack of competition from other predators indicate that this novel niche utilization may support the drastic increase of pelagic stickleback in the Baltic Sea.We find that copepods dominated the gut content of the two clupeids sprat and herring. Pseudocalanus and Temora occupied most of the sequence reads of the clupeid metabarcoding, two species that are often reported as preferred prey in previous studies11,12. Despite high contributions of these two copepods, Pseudocalanus was more than four times as abundant as Temora in clupeid gut contents. A strong preference for this copepod with marine origin can further confirm the increased competition between the clupeids, as Pseudocalanus has decreased due to decreased salinity12 and shares a similar vertical distribution as clupeid during daytime27. Our study using metabarcoding further reveals a large relative quantity (11%) of the ctenophore Mertensia in the gut samples of both clupeids. Similar, Clarke et al.28 reported an important contribution of gelatinous zooplankton to upper trophic levels in the Southern Ocean. Despite high abundances of ctenophores in the Baltic Sea and their assumed importance in marine food webs19, they are not reported as food for planktivorous fish. A possible explanation is the difficulty observing them microscopically, as their digestion rate is faster than crustaceans29, and no hard parts remain in the digestive system. Further, COI detected the presence of cladocerans, which was confirmed by the microscopic survey, but underrepresented with 18S that strongly amplify copepods20. Interestingly, more than twice annelid COI reads, including the benthic macroinvertebrates Bylgides and Marenzellaria, were associated to stickleback (15%) and herring (8%) than to sprat (4%), highlighting their ability to migrate vertically. These interactions suggest that together stickleback and herring contribute to benthic-pelagic coupling when oxygen is not restricting vertical migration in the southern Baltic Sea30.Sprat and herring share a similar feeding niche, which may explain previously observed declines in body mass and stomach fullness, and supports the theory of competition between the two species31. In contrast, stickleback revealed little diet overlap with the other mesopredators. The low relative abundances of Pseudocalanus (1–8%) in metabarcoding analyses indicates that the density-dependent competition may not limit the population growth of stickleback. The copepods that were shared in the diet of stickleback, sprat, and herring were Temora, Acartia, and Centropages have increased over the last decades, as opposed to Pseudocalanus32. Our results show that stickleback are able to feed on a broader spectrum of prey and highlight that stickleback utilizes the rotifer Synchaeta baltica as prey, which is an important component of the plankton community composition in the Baltic Sea18,20. Due to the difference of prey size, we can expect an overrepresentation of copepod to rotifer sequences compared with microscopic count data. High predation rate on S. baltica is supported by both the qPCR assay as well as microscopic counts, although only the eggshells were visible but not the soft-bodied rotifer. Despite the considerably lower carbon content per S. baltica (ca. 6 µg C ind−1) compared to copepods (ca. 20 µg C ind−1)33, the high number of rotifers likely act as a major food source for stickleback. These results propose that stickleback, due to their opportunistic feeding behaviour34 and smaller size35, have a distinct feeding niche from sprat and herring in the open water, as they feed on a smaller size class of zooplankton compared to the clupeids. Thus, we cannot assume the same process of competition between clupeids and stickleback.Rotifers can at times be very abundant in the Baltic Sea, reaching densities up to 25,000 ind m−3, but their natural predators are poorly studied. An increasing trend in biomass of the two main rotifer genera (Synchaeta and Keratella) was observed since the 1990s36. In a recent study, we showed that rotifers might occupy a unique feeding niche, as direct grazers of dinoflagellate spring bloom, as well as in the recycling of organic matter in summer20. The low level of predation on rotifers by clupeid adults ( More

Ecosystem services are the contributions of nature to human well-being — for example, the provision of raw materials, carbon sequestration and recreation. Although relatively new, the study of these essential services has developed rapidly and is now included in many global policies and assessments. However, mapping and modelling these services is restricted by the availability of data that can account for the multidimensional traits of ecosystem services and model coupled social–ecological systems. Traditional datasets, including surveys, interviews, and focus groups, are often not viable on the scale necessary for many ecosystem service assessments. More

We surveyed three islands of Lake Chany: Uzkoredkii (54° 58′ 15′′ N, 77°27′04′′ E), Reden’kii (54° 56′ 05′′ N, 77° 22′ 27′′ 52 E), Korablik (54° 59′ 31′′ N, 77° 40′ 38′′ E). The studied intertidal habitats are rarely reached by humans.Gull nests were counted in colonies by regular surveys over eight years (1993, 1994, 1996–1998, 2001–2003) on the islands of Lake Chany. Colonies were visited daily or sometimes every other day. To minimize the disturbance caused by the investigation, the time spent working, within view of the gulls was restricted to a maximum of forty minutes per study plots. We noted nest content at every visit for the presence of eggs or chicks. In total, there were 1 164 nests under observation. Nests contained 1 (n = 140), 2 (n = 518), 3 (n = 504) or 4 (n = 2) eggs. Modal clutch size of the great black-headed gull is two or three eggs, varying seasonally. The length and width of the eggs were measured using Vernier calipers (division accuracy 0,1 mm) and numbered with a waterproof marker. Egg volumes were estimated using Hoyt’s equation: Volume = 0.51 * Length * Width * Width/100013. We determined the volume of 2117 great black-headed gull eggs.As the laying of eggs has already started by the first visit to the colony, the date of the beginning of egg laying was calculated by subtracting the average length of the incubation period of great black-headed gulls (27 days) from the hatching date of first chick in the nest (n = 559 nests). If the hatching date was not known, the clutch initiation date was determined by subtracting the number of days of incubation from the date that the nest was first discovered (n = 469 nests). The stage of incubation was estimated from the change in position of an incubated egg placed in water14,15. The technique’s accuracy varied throughout incubation and mean prediction error fall between 0–4 days. On average, egg flotation estimated an embryo’s developmental age to within 1.9 ± 1.6 days (mean ± 1 SD)16. Only 47 nests were found during egg laying. Great black-headed gulls usually laid eggs at intervals of two days. Incubation started as soon as the first egg was laid, so eggs hatched asynchronously, one or two days apart.Whenever possible, we determined the within-clutch laying sequence of eggs (1st, 2nd, 3rd, and 4th). A complete laying sequence was established by observation in 47 cases. In about 48% of clutches the position in laying sequence was established on the basis of the sequence of hatching. In other cases, if we could distinguish within-clutch distinct flotation levels of eggs, we numbered eggs according to the stage of incubation. Sometimes this technique for distinguishing egg laying order were used in other seabirds17,18.We recorded the pipping date (i.e. appearance of star-like bursts) and the actual hatching date of the individual eggs. Wet chicks were registered as hatchlings of that day; dry chicks were registered as 1 day old. Chicks older than two days left the nest and moved to a location nearby. Newly hatched gull chicks were captured by hand at nests, ringed, and measured. We determined wing, tarsus, and head length using a ruler with zero-stop and vernier calipers and body weight measured using Pesola spring balances for 747 chicks of great black-headed gulls, and 457 of them hatched from eggs that were measured. More

Turroni F, Milani C, Duranti S, Lugli GA, Bernasconi S, Margolles A, et al. The infant gut microbiome as a microbial organ influencing host well-being. Ital J Pediatr. 2020;46:1–13.Article 

Google Scholar 
Bokulich NA, Chung J, Battaglia T, Henderson N, Jay M, Li H, et al. Antibiotics, birth mode, and diet shape microbiome maturation during early life. Sci Transl Med. 2016;8:343ra82–343ra82.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Akagawa S, Tsuji S, Onuma C, Akagawa Y, Yamaguchi T, Yamagishi M, et al. Effect of delivery mode and nutrition on gut microbiota in neonates. Ann Nutr Metab. 2019;74:132–9.CAS 
PubMed 
Article 

Google Scholar 
Rothschild D, Weissbrod O, Barkan E, Kurilshikov A, Korem T, Zeevi D, et al. Environment dominates over host genetics in shaping human gut microbiota. Nature 2018;555:210–5.Vellend M, Srivastava DS, Anderson KM, Brown CD, Jankowski JE, Kleynhans EJ, et al. Assessing the relative importance of neutral stochasticity in ecological communities. Oikos. 2014;123:1420–30.Article 

Google Scholar 
Fukami T. Historical contingency in community assembly: integrating niches, species pools, and priority effects. Annu Rev Ecol Evol Syst. 2015;46:1–23.Article 

Google Scholar 
Fukami T, Nakajima M. Community assembly: Alternative stable states or alternative transient states? Ecol Lett. 2011;14:973–84.PubMed 
PubMed Central 
Article 

Google Scholar 
Sprockett D, Fukami T, Relman DA. Role of priority effects in the early-life assembly of the gut microbiota. Nat Rev Gastroenterol Hepatol. 2018;15:197–205.PubMed 
PubMed Central 
Article 

Google Scholar 
Debray R, Herbert RA, Jaffe AL, Crits-Christoph A, Power ME, Koskella B. Priority effects in microbiome assembly. Nat Rev Microbiol. 2022;20:109–21.CAS 
PubMed 
Article 

Google Scholar 
Tannock GW, Lawley B, Munro K, Pathmanathan SG, Zhou SJ, Makrides M, et al. Comparison of the compositions of the stool microbiotas of infants fed goat milk formula, cow milk-based formula, or breast milk. Appl Environ Microbiol. 2013;79:3040–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Matsuki T, Yahagi K, Mori H, Matsumoto H, Hara T, Tajima S, et al. A key genetic factor for fucosyllactose utilization affects infant gut microbiota development. Nat Commun. 2016;7:11939.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sakanaka M, Hansen ME, Gotoh A, Katoh T, Yoshida K, Odamaki T, et al. Evolutionary adaptation in fucosyllactose uptake systems supports bifidobacteria-infant symbiosis. Sci Adv. 2019;5:eaaw7696.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Engfer MB, Stahl B, Finke B, Sawatzki G, Daniel H. Human milk oligosaccharides are resistant to enzymatic hydrolysis in the upper gastrointestinal tract. Am J Clin Nutr. 2000;71:1589–96.CAS 
PubMed 
Article 

Google Scholar 
Macrobal A, Sonnenburg JL. Human milk oligosaccharide consumption by intestinal microbiota. Clin Microbiol Infect. 2012;18:12–15.Article 

Google Scholar 
Macrobal A, Barboza M, Froehlich JW, Block DE, German JB, Lebrilla CB, et al. Consumption of human milk oligosaccharides by gut-related microbes. J Agric Food Chem. 2010;58:5334–40.Article 
CAS 

Google Scholar 
Sakanaka M, Gotoh A, Yoshida K, Odamaki T, Koguchi H, Xiao JZ, et al. Varied pathways of infant gut-associated Bifidobacterium to assimilate human milk oligosaccharides: prevalence of the gene set and its correlation with bifidobacteria-rich microbiota formation. Nutrients. 2020;12:71.CAS 
Article 

Google Scholar 
Katayama T. Host-derived glycans serve as selected nutrients for the gut microbe: human milk oligosaccharides and bifidobacteria. Biosci Biotechnol Biochem. 2016;80:621–32.CAS 
PubMed 
Article 

Google Scholar 
Turroni F, Foroni E, Pizzetti P, Giubellini V, Ribbera A, Merusi P, et al. Exploring the diversity of the bifidobacterial population in the human intestinal tract. Appl Environ Microbiol. 2009;75:1534–45.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gore C, Munro K, Lay C, Bibiloni R, Morris J, Woodcock A, et al. Bifidobacterium pseudocatenulatum is associated with atopic eczema: A nested case-control study investigating the fecal microbiota of infants. J Allergy Clin Immunol. 2008;121:135–40.PubMed 
Article 

Google Scholar 
Lewis ZT, Mills DA. Differential establishment of bifidobacteria in the breastfed infant gut. Nestle Nutr Inst Work Ser. 2017;88:149–59.Article 

Google Scholar 
Tannock GW, Lee PS, Wong KH, Lawley B. Why don’t all infants have bifidobacteria in their stool? Front Microbiol. 2016;7:6–10.Article 

Google Scholar 
Reyman M, van Houten MA, van Baarle D, Bosch AATM, Man WH, Chu MLJN, et al. Impact of delivery mode-associated gut microbiota dynamics on health in the first year of life. Nat Commun. 2019;10:1–12.CAS 
Article 

Google Scholar 
Underwood MA, Kalanetra KM, Bokulich NA, Lewis ZT, Mirmiran M, Tancredi DJ, et al. A comparison of two probiotic strains of bifidobacteria in preterm infants. J Pediatr. 2013;163:1585–91.CAS 
PubMed 
Article 

Google Scholar 
Plummer EL, Bulach DM, Murray GL, Jacobs SE, Tabrizi SN, Garland SM. Gut microbiota of preterm infants supplemented with probiotics: sub-study of the ProPrems trial. BMC Microbiol. 2018;18:1–8.Article 
CAS 

Google Scholar 
Kitajima H, Sumida Y, Tanaka R, Yuki N, Takayama H, Fujimura M. Early administration of Bifidobacterium breve to preterm infants: Randomised controlled trial. Arch Dis Child Fetal Neonatal Ed. 1997;76:101–7.Article 

Google Scholar 
Ojima MN, Yoshida K, Sakanaka M, Jiang L, Odamaki T, Katayama T. Ecological and molecular perspectives on responders and non-responders to probiotics and prebiotics. Curr Opin Biotechnol. 2022;73:108–20.CAS 
PubMed 
Article 

Google Scholar 
Suez J, Zmora N, Segal E, Elinav E. The pros, cons, and many unknowns of probiotics. Nat Med. 2019;25:716–29.CAS 
PubMed 
Article 

Google Scholar 
Costeloe K, Hardy P, Juszczak E, Wilks M, Millar MR. Bifidobacterium breve BBG-001 in very preterm infants: A randomised controlled phase 3 trial. Lancet. 2016;387:649–60.PubMed 
Article 

Google Scholar 
Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005;33:5691–702.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Katoh T, Ojima MN, Sakanaka M, Ashida H, Gotoh A, Katayama T. Enzymatic adaptation of Bifidobacterium bifidum to host glycans, viewed from glycoside hydrolyases and carbohydrate-binding modules. Microorganisms. 2020;8:481.CAS 
PubMed Central 
Article 

Google Scholar 
Egan M, Motherway MO, Kilcoyne M, Kane M, Joshi L, Ventura M, et al. Cross-feeding by Bifidobacterium breve UCC2003 during co-cultivation with Bifidobacterium bifidum PRL2010 in a mucin-based medium. BMC Microbiol. 2014;14:1–14.Article 
CAS 

Google Scholar 
Higgins MA, Ryan KS. Generating a fucose permease deletion mutant in Bifidobacterium longum subspecies infantis ATCC 15697. Anaerobe. 2021;68:102320.CAS 
PubMed 
Article 

Google Scholar 
O’Connell Motherway M, Kinsella M, Fitzgerald GF, van Sinderen D. Transcriptional and functional characterization of genetic elements involved in galacto-oligosaccharide utilization by Bifidobacterium breve UCC2003. Micro Biotechnol. 2013;6:67–79.Article 
CAS 

Google Scholar 
Yoshida E, Sakurama H, Kiyohara M, Nakajima M, Kitaoka M, Ashida H, et al. Bifidobacterium longum subsp. infantis uses two different β-galactosidases for selectively degrading type-1 and type-2 human milk oligosaccharides. Glycobiology. 2012;22:361–8.CAS 
PubMed 
Article 

Google Scholar 
Vannette RL, Fukami T. Historical contingency in species interactions: Towards niche-based predictions. Ecol Lett. 2014;17:115–24.PubMed 
Article 

Google Scholar 
Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe. 2015;17:690–703.PubMed 
Article 
CAS 

Google Scholar 
Pu Z, Jiang L. Dispersal among local communities does not reduce historical contingencies during metacommunity assembly. Oikos. 2015;124:1327–36.Article 

Google Scholar 
Chase JM. Community assembly: when should history matter? Oecologia. 2003;136:489–98.PubMed 
Article 

Google Scholar 
Schröder A, Persson L, De Roos AM. Direct experimental evidence for alternative stable states: A review. Oikos. 2005;110:3–19.Article 

Google Scholar 
Asakuma S, Hatakeyama E, Urashima T, Yoshida E, Katayama T, Yamamoto K, et al. Physiology of consumption of human milk oligosaccharides by infant gut-associated bifidobacteria. J Biol Chem. 2011;286:34583–92.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gotoh A, Katoh T, Sakanaka M, Ling Y, Yamada C, Asakuma S, et al. Sharing of human milk oligosaccharides degradants within bifidobacterial communities in faecal cultures supplemented with Bifidobacterium bifidum. Sci Rep. 2018;8:13958.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ashida H, Miyake A, Kiyohara M, Wada J, Yoshida E, Kumagai H, et al. Two distinct α-l-fucosidases from Bifidobacterium bifidum are essential for the utilization of fucosylated milk oligosaccharides and glycoconjugates. Glycobiology. 2009;19:1010–7.CAS 
PubMed 
Article 

Google Scholar 
Roger LC, Costabile A, Holland DT, Hoyles L, McCartney AL. Examination of faecal Bifidobacterium populations in breast- and formula-fed infants during the first 18 months of life. Microbiology. 2010;156:3329–41.CAS 
PubMed 
Article 

Google Scholar 
Avershina E, Storrø O, Øien T, Johnsen R, Wilson R, Egeland T, et al. Bifidobacterial succession and correlation networks in a large unselected cohort of mothers and their children. Appl Environ Microbiol. 2013;79:497–507.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Turroni F, Peano C, Pass DA, Foroni E, Severgnini M, Claesson MJ, et al. Diversity of bifidobacteria within the infant gut microbiota. PLoS One. 2012;7:20–4.Article 
CAS 

Google Scholar 
James K, Bottacini F, Contreras JIS, Vigoureux M, Egan M, Motherway MO, et al. Metabolism of the predominant human milk oligosaccharide fucosyllactose by an infant gut commensal. Sci Rep. 2019;9:1–20.Article 
CAS 

Google Scholar 
Dedon LR, Özcan E, Rani A, Sela DA. Bifidobacterium infantis metabolizes 2′fucosyllactose-derived and free fucose through a common catabolic pathway resulting in 1,2-propanediol secretion. Front Nutr. 2020;7:1–16.Article 
CAS 

Google Scholar 
Sprockett D, Martin M, Costello E, Burns A, Holmes S, Gurven M, et al. Microbiota assembly, structure, and dynamics among tsimane horticulturalists of the Bolivian Amazon. Nat Commun. 2019;11:1–14.Laursen MF, Sakanaka M, von Burg N, Mörbe U, Andersen D, Moll JM, et al. Bifidobacterium species associated with breastfeeding produce aromatic lactic acids in the infant gut. Nat Microbiol. 2021;6:1367–82.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Meng D, Sommella E, Salviati E, Campiglia P, Ganguli K, Djebali K, et al. Indole-3-lactic acid, a metabolite of tryptophan, secreted by Bifidobacterium longum subspecies infantis is anti-inflammatory in the immature intestine. Pediatr Res. 2020;88:209–17.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bunesova V, Lacroix C, Schwab C. Fucosyllactose and l-fucose utilization of infant Bifidobacterium longum and Bifidobacterium kashiwanohense. BMC Microbiol. 2016;16:248.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ruiz-Moyano S, Totten SM, Garrido D, Smilowitz JT, Bruce German J, Lebrilla CB, et al. Variation in consumption of human milk oligosaccharides by infant gut-associated strains of Bifidobacterium breve. Appl Environ Microbiol. 2013;79:6040–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lawson MAE, O’Neill IJ, Kujawska M, Gowrinadh Javvadi S, Wijeyesekera A, Flegg Z, et al. Breast milk-derived human milk oligosaccharides promote Bifidobacterium interactions within a single ecosystem. ISME J. 2020;14:635–48.CAS 
PubMed 
Article 

Google Scholar 
Schwab C, Ruscheweyh HJ, Bunesova V, Pham VT, Beerenwinkel N, Lacroix C. Trophic interactions of infant bifidobacteria and Eubacterium hallii during l-fucose and fucosyllactose degradation. Front Microbiol. 2017;8:1–14.Article 

Google Scholar 
Engels C, Ruscheweyh HJ, Beerenwinkel N, Lacroix C, Schwab C. The common gut microbe Eubacterium hallii also contributes to intestinal propionate formation. Front Microbiol. 2016;7:1–12.Article 

Google Scholar 
Marcobal A, Barboza M, Sonnenburg ED, Pudlo N, Martens EC, Desai P, et al. Bacteroides in the infant gut consume milk oligosaccharides via mucus-utilization pathways. Cell Host Microbe. 2011;10:507–14.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vatanen T, Kostic AD, D’Hennezel E, Siljander H, Franzosa EA, Yassour M, et al. Variation in Microbiome LPS Immunogenicity Contributes to Autoimmunity in Humans. Cell. 2016;165:842–53.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vatanen T, Franzosa EA, Schwager R, Tripathi S, Arthur TD, Vehik K, et al. The human gut microbiome in early-onset type 1 diabetes from the TEDDY study. Nature. 2018;562:589–94.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Moossavi S, Sepehri S, Robertson B, Bode L, Goruk S, Field CJ, et al. Composition and Variation of the Human Milk Microbiota Are Influenced by Maternal and Early-Life Factors. Cell Host Microbe. 2019;25:324–335.e4.CAS 
PubMed 
Article 

Google Scholar 
Martín R, Langa S, Reviriego C, Jiménez E, Marín ML, Xaus J, et al. Human milk is a source of lactic acid bacteria for the infant gut. J Pediatr. 2003;143:754–8.PubMed 
Article 

Google Scholar 
Martín R, Jiménez E, Heilig H, Fernández L, Marín ML, Zoetendal EG, et al. Isolation of bifidobacteria from breast milk and assessment of the bifidobacterial population by PCR-denaturing gradient gel electrophoresis and quantitative real-time PCR. Appl Environ Microbiol. 2009;75:965–9.PubMed 
Article 
CAS 

Google Scholar 
Heikkilä MP, Saris PEJ. Inhibition of Staphylococcus aureus by the commensal bacteria of human milk. J Appl Microbiol. 2003;95:471–8.PubMed 
Article 
CAS 

Google Scholar 
Li Y, Shimizu T, Hosaka A, Kaneko N, Ohtsuka Y, Yamashiro Y. Effects of Bifidobacterium breve supplementation on intestinal flora of low birth weight infants. Pediatr Int. 2004;46:509–15.PubMed 
Article 

Google Scholar 
Nishimoto M, Kitaoka M. Practical preparation of lacto-N-biose I, a candidate for the bifidus factor in human milk. Biosci Biotechnol Biochem. 2007;71:2101–4.CAS 
PubMed 
Article 

Google Scholar 
Duncan SH, Hold GL, Harmsen HJM, Stewart CS, Flint HJ. Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. Int J Syst Evol Microbiol. 2002;52:2141–6.CAS 
PubMed 

Google Scholar 
Browne HP, Forster SC, Anonye BO, Kumar N, Neville BA, Stares MD, et al. Culturing of ‘unculturable’ human microbiota reveals novel taxa and extensive sporulation. Nature. 2016;533:543–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tanizawa Y, Fujisawa T, Kaminuma E, Nakamura Y. Arita M. DFAST and DAGA: Web-based integrated genome annotation tools and resources. Biosci Microbiota, Food Heal. 2016;35:173–84.CAS 
Article 

Google Scholar 
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42:206–14.Article 
CAS 

Google Scholar 
Price MN, Arkin AP. PaperBLAST: Text-mining papers for information about homologs. bioRxiv. 2017;2:1–10.
Google Scholar 
Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42:490–5.Article 
CAS 

Google Scholar 
Saier MH, Reddy VS, Tsu BV, Ahmed MS, Li C, Moreno-Hagelsieb G. The Transporter Classification Database (TCDB): Recent advances. Nucleic Acids Res. 2016;44:D372–9.CAS 
PubMed 
Article 

Google Scholar 
Martínez I, Wallace G, Zhang C, Legge R, Benson AK, Carr TP, et al. Diet-induced metabolic improvements in a hamster model of hypercholesterolemia are strongly linked to alterations of the gut microbiota. Appl Environ Microbiol. 2009;75:4175–84.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Anumula KR. Advances in fluorescence derivatization methods for high-performance liquid chromatographic analysis of glycoprotein carbohydrates. Anal Biochem. 2006;350:1–23.CAS 
PubMed 
Article 

Google Scholar 
Cohenford MA, Abraham A, Abraham J, Dain JA. Colorimetric assay for free and bound l-fucose. Anal Biochem. 1989;177:172–7.CAS 
PubMed 
Article 

Google Scholar 
Kato K, Odamaki T, Mitsuyama E, Sugahara H, Xiao JZ, Osawa R. Age-related changes in the composition of gut Bifidobacterium species. Curr Microbiol. 2017;74:987–95.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.CAS 
PubMed 
PubMed Central 
Article 

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

Google Scholar 
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20:1–13.Article 
CAS 

Google Scholar 
Lu J, Breitwieser FP, Thielen P, Salzberg SL. Bracken: Estimating species abundance in metagenomics data. PeerJ Comput Sci. 2017;2017:1–17.CAS 

Google Scholar 
Milani C, Lugli GA, Fontana F, Mancabelli L, Alessandri G, Longhi G, et al. METAnnotatorX2: A comprehensive tool for deep and shallow metagenomic data set analyses. mSystems. 2021;6:1–15.Article 

Google Scholar 
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. 2019. More

Anderson, P. K. et al. Emerging infectious diseases of plants: Pathogen pollution, climate change and agrotechnology drivers. Trends Ecol. Evol. 19, 535–544 (2004).PubMed 

Google Scholar 
Brasier, C. M. The biosecurity threat to the UK and global environment from international trade in plants. Plant Pathol. 57, 792–808 (2008).
Google Scholar 
Waage, J. K. & Mumford, J. D. Agricultural biosecurity. Philos. Trans. R. Soc. Lond. B Biol. Sci. 363, 863–876 (2008).CAS 
PubMed 

Google Scholar 
IPPC. Surveillance guide—A guide to understand the principal requirements of surveillance programmes for national plant protection organizations. Second edition. http://www.fao.org/documents/card/en/c/cb7139en (2021) https://doi.org/10.4060/cb7139en.Parnell, S., van den Bosch, F., Gottwald, T. & Gilligan, C. A. Surveillance to inform control of emerging plant diseases: An epidemiological perspective. Annu. Rev. Phytopathol. 55, 591–610 (2017).CAS 
PubMed 

Google Scholar 
Cunniffe, N. J., Cobb, R. C., Meentemeyer, R. K., Rizzo, D. M. & Gilligan, C. A. Modeling when, where, and how to manage a forest epidemic, motivated by sudden oak death in California. Proc. Natl. Acad. Sci. 113, 5640–5645 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gottwald, T. R., Dixon, W., Parnell, S. & Riley, T. Huanglongbing: The dragon arrives in the USA. In Huanglongbing-Greening International Workshop, July 14–21 13–14 (2006).Herms, D. A., Stone, A. K. & Chatfield, J. A. Emerald ash borer: The beginning of the end of ash in North America?. Ornam. Plants Annu. Rep. Res. Rev. 2003, 62–71 (2004).
Google Scholar 
Sansford, C. E. Pest Risk Analysis for Hymenoscyphus pseudoalbidus (anamorph Chalara fraxinea) for the UK and the Republic of Ireland. https://webarchive.nationalarchives.gov.uk/ukgwa/20140904094312mp_/http://www.fera.defra.gov.uk/plants/plantHealth/pestsDiseases/documents/hymenoscyphusPseudoalbidusPRA.pdf (2013).Alonso Chavez, V., Parnell, S. & van den Bosch, F. Monitoring invasive pathogens in plant nurseries for early-detection and to minimise the probability of escape. J. Theor. Biol. 407, 290–302 (2016).ADS 
MathSciNet 
PubMed 
MATH 

Google Scholar 
Bourhis, Y., Gottwald, T. R., Lopez-Ruiz, F. J., Patarapuwadol, S. & van den Bosch, F. Sampling for disease absence-deriving informed monitoring from epidemic traits. J. Theor. Biol. 461, 8–16 (2019).ADS 
MathSciNet 
PubMed 
MATH 

Google Scholar 
Mastin, A. J., van den Bosch, F., van den Berg, F. & Parnell, S. Quantifying the hidden costs of imperfect detection for early detection surveillance. Philos. Trans. R. Soc. B Biol. Sci. 374, 20180261 (2019).
Google Scholar 
Mastin, A. J., van den Bosch, F., Gottwald, T. R., Alonso Chavez, V. & Parnell, S. R. A method of determining where to target surveillance efforts in heterogeneous epidemiological systems. PLoS Comput. Biol. 13, e1005712 (2017).PubMed 
PubMed Central 

Google Scholar 
Parnell, S., Gottwald, T. R., Gilks, W. R. & van den Bosch, F. Estimating the incidence of an epidemic when it is first discovered and the design of early detection monitoring. J. Theor. Biol. 305, 30–36 (2012).ADS 
MathSciNet 
CAS 
PubMed 
MATH 

Google Scholar 
Parnell, S., Gottwald, T. R., Cunniffe, N. J., Alonso Chavez, V. & van den Bosch, F. Early detection surveillance for an emerging plant pathogen: A rule of thumb to predict prevalence at first discovery. Proc. R. Soc. B Biol. Sci. 282, 20151478 (2015).
Google Scholar 
Silva, G. et al. Plant pest surveillance: From satellites to molecules. Emerg. Top. Life Sci. 5, 275–287 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Mastin, A. J., Gottwald, T. R., van den Bosch, F., Cunniffe, N. J. & Parnell, S. Optimising risk-based surveillance for early detection of invasive plant pathogens. PLoS Biol. 18, e3000863 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Martelli, G. P., Boscia, D., Porcelli, F. & Saponari, M. The olive quick decline syndrome in south-east Italy: A threatening phytosanitary emergency. Eur. J. Plant Pathol. 144, 235–243 (2015).
Google Scholar 
Saponari, M., Boscia, D., Nigro, F. & Martelli, G. P. Identification of DNA sequences related to Xylella fastidiosa in oleander, almond and olive trees exhibiting leaf scorch symptoms in Apulia (southern Italy). J. Plant Pathol. 95, 668 (2013).
Google Scholar 
Ben Moussa, I. E. et al. Seasonal fluctuations of sap-feeding insect species infected by Xylella fastidiosa in Apulian olive groves of southern Italy. J. Econ. Entomol. 109, 1512–1518 (2016).PubMed 

Google Scholar 
Cornara, D. et al. Transmission of Xylella fastidiosa to grapevine by the meadow spittlebug. Phytopathology 106, 1285–1290 (2016).CAS 
PubMed 

Google Scholar 
Cornara, D. et al. Transmission of Xylella fastidiosa by naturally infected Philaenus spumarius (Hemiptera, Aphrophoridae) to different host plants. J. Appl. Entomol. 141, 80–87 (2017).
Google Scholar 
Saponari, M. et al. Infectivity and transmission of Xylella fastidiosa by Philaenus spumarius (Hemiptera: Aphrophoridae) in Apulia, Italy. J. Econ. Entomol. 107, 1316–1319 (2014).PubMed 

Google Scholar 
European Commission. Commission Implementing Regulation (EU) 2020/1201 of 14 August 2020 as regards measures to prevent the introduction into and the spread within the Union of Xylella fastidiosa (Wells et al.). (2021).EFSA et al. Guidelines for statistically sound and risk-based surveys of Xylella fastidiosa. EFSA. Support. Publ. 17, 1873 (2020).EFSA et al. General guidelines for statistically sound and risk-based surveys of plant pests. EFSA Support. Publ. 17, 1919E (2020).Bourhis, Y., Gottwald, T. & van den Bosch, F. Translating surveillance data into incidence estimates. Philos. Trans. R. Soc. B Biol. Sci. 374, 20180262 (2019).CAS 

Google Scholar 
Cornara, D. et al. Spittlebugs as vectors of Xylella fastidiosa in olive orchards in Italy. J. Pest Sci. 90, 521–530 (2017).
Google Scholar 
Cornara, D., Bosco, D. & Fereres, A. Philaenus spumarius: when an old acquaintance becomes a new threat to European agriculture. J. Pest Sci. 91, 957–972 (2018).
Google Scholar 
Almeida, R. P. P., Blua, M. J., Lopes, J. R. S. & Purcell, A. H. Vector transmission of Xylella fastidiosa: Applying fundamental knowledge to generate disease management strategies. Ann. Entomol. Soc. Am. 98, 775–786 (2005).
Google Scholar 
Purcell, A. H. & Finlay, A. H. Evidence for noncirculative transmission of Pierce’s disease bacterium by sharpshooter leafhoppers. Phytopathology 69, 393–395 (1979).
Google Scholar 
Hill, B. & Purcell, A. H. Acquisition and retention of Xylella fastidiosa by an efficient vector, Graphocephala atropunctata. Phytopathology 85, 209 (1995).
Google Scholar 
Hill, B. L. & Purcell, A. H. Multiplication and movement of Xylella fastidiosa within grapevine and four other plants. Phytopathology 85, 1368 (1995).
Google Scholar 
Huang, Q., Bentz, J. & Sherald, J. L. Fast, easy and efficient DNA extraction and one-step polymerase chain reaction for the detection of Xylella fastidiosa in potential insect vectors. J. Plant Pathol. 88, 77–81 (2006).CAS 

Google Scholar 
Harper, S. J., Ward, L. I. & Clover, G. R. G. Development of LAMP and real-time PCR methods for the rapid detection of Xylella fastidiosa for quarantine and field applications. Phytopathology 100, 1282–1288 (2010).CAS 
PubMed 

Google Scholar 
EFSA et al. Pest survey card on Xylella fastidiosa. EFSA Support. Publ. 16, (2019).Fierro, A., Liccardo, A. & Porcelli, F. A lattice model to manage the vector and the infection of the Xylella fastidiosa on olive trees. Sci. Rep. 9, 8723 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
EPPO. PM 7/24 (4) Xylella fastidiosa. EPPO Bull. 49, 175–227 (2019).Landa, B. B. et al. Emergence of a plant pathogen in Europe associated with multiple intercontinental introductions. Appl. Environ. Microbiol. 86, 1–15 (2019).
Google Scholar 
Castro, C., DiSalvo, B. & Roper, M. C. Xylella fastidiosa: A reemerging plant pathogen that threatens crops globally. PLoS Pathog. 17, e1009813 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Saponari, M., Giampetruzzi, A., Loconsole, G., Boscia, D. & Saldarelli, P. Xylella fastidiosa in olive in Apulia: Where we stand. Phytopathology 109, 175–186 (2019).CAS 
PubMed 

Google Scholar 
Zarco-Tejada, P. J. et al. Previsual symptoms of Xylella fastidiosa infection revealed in spectral plant-trait alterations. Nat. Plants 4, 432–439 (2018).CAS 
PubMed 

Google Scholar 
Gottwald, T. et al. Canine olfactory detection of a vectored phytobacterial pathogen, Liberibacter asiaticus, and integration with disease control. Proc. Natl. Acad. Sci. 117, 3492–3501 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Mendel, J., Furton, K. G. & Mills, D. An Evaluation of scent-discriminating canines for rapid response to agricultural diseases. HortTechnology 28, 102–108 (2018).
Google Scholar 
ECDC. Guidelines for the Surveillance of Invasive Mosquitoes in Europe. (2012).Kading, R. C., Golnar, A. J., Hamer, S. A. & Hamer, G. L. Advanced surveillance and preparedness to meet a new era of invasive vectors and emerging vector-borne diseases. PLoS Negl. Trop. Dis. 12, e0006761 (2018).PubMed 
PubMed Central 

Google Scholar 
Kumagai, L. B. et al. First report of Candidatus Liberibacter asiaticus associated with citrus huanglongbing in California. Plant Dis. 97, 283 (2013).CAS 
PubMed 

Google Scholar 
Ben Moussa, I. E. et al. Evaluation of “Spy Insect” approach for monitoring Xylella fastidiosa in symptomless olive orchards in the Salento peninsula (Southern Italy). IOBC WPRS Bull. 121, 77–84 (2017).
Google Scholar 
Cruaud, A. et al. Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: A case study in Corsica. Sci. Rep. 8, 15628 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
Yaseen, T. et al. On-site detection of Xylella fastidiosa in host plants and in “spy insects” using the real-time loop-mediated isothermal amplification method. Phytopathol. Mediterr. https://doi.org/10.14601/Phytopathol_Mediterr-15250 (2015).Article 

Google Scholar 
López-Mercadal, J. et al. Collection of data and information in Balearic Islands on biology of vectors and potential vectors of Xylella fastidiosa (GP/EFSA/ALPHA/017/01). EFSA Support. Publ. 18, 6925E (2021).
Google Scholar 
Cunty, A. Detection, identification and surveillance of Xylella fastidiosa on vectors in France https://zenodo.org/record/3551122#.XjGqBs77SUl. (2019) https://doi.org/10.5281/zenodo.3551122.Kottelenberg, D., Hemerik, L., Saponari, M. & van der Werf, W. Shape and rate of movement of the invasion front of Xylella fastidiosa spp. pauca in Puglia. Sci. Rep. 11, 1061 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar  More

Sampling, DNA preparation and sequencingStockholmSamples LOW002, LOW003, LOW006, LOW007, LOW008 and PON012 were processed at the Archaeological Research Laboratory at Stockholm University, Sweden, following methods previously described8. In brief, this involved extracting DNA by incubating the bone powder for 24 h at 37 °C in 1.5 ml of digestion buffer (0.45 M EDTA (pH 8.0) and 0.25 mg ml–1 proteinase K), concentrating supernatant on Amicon Ultra-4 (30-kDa molecular weight cut-off (MWCO)) filter columns (MerckMillipore) and purifying on Qiagen MinElute columns. Double-stranded Illumina libraries were prepared using the protocol outlined in ref. 48, with the inclusion of USER enzyme and the modifications described in ref. 49.Samples 367, PDM100, Taimyr-1 and Yana-1 were processed at the Swedish Museum of Natural History in Stockholm, Sweden, following previously described methods8. In brief, this involved extracting DNA using a silica-based method with concentration on Vivaspin filters (Sartorius), according to a protocol optimized for recovery of ancient DNA50. Double-stranded Illumina libraries were prepared using the protocol outlined in ref. 48, with the inclusion of USER enzyme.Samples ALAS_024, VAL_033, ALAS_016, VAL_008, HMNH_007, HMNH_011, VAL_050, VAL_005, DS04, VAL_037, VAL_012, VAL_011, VAL_18A, IN18_016 and IN18_005 were processed at the Swedish Museum of Natural History in Stockholm, Sweden, following previously described methods for permafrost bone and tooth samples51. In brief, this involved DNA extraction using the methodology of ref. 52 and double-stranded Illumina library preparation as described in ref. 48, with dual unique indexes and the inclusion of USER enzyme. Between eight and ten separate PCR reactions with unique indexes were carried out for each sample to maximize library complexity. The libraries were sequenced alongside samples HOV4, AL2242, AL2370, AL2893, AL3272 and AL3284 across three Illumina NovaSeq 6000 lanes with an S4 100-bp paired-end set-up at SciLifeLab in Stockholm.PotsdamSamples JAL48, JAL65, JAL69, JAL358, AH574, AH575 and AH577 were processed at the University of Potsdam. Pre-amplification steps (DNA extraction and library preparation) were conducted in separated laboratory rooms specially equipped for the processing of ancient DNA. Amplification and post-amplification steps were performed in different laboratory rooms. DNA was extracted from bone powder (29–54 mg) following a protocol specially adapted to recover short DNA fragments52. Single-stranded double-indexed libraries were built from 20 µl of DNA extract according to the protocol in ref. 53. The libraries were sequenced on an HiSeq X platform at SciLifeLab in Stockholm.Tübingen/JenaSamples JK2174, JK2175, JK2179, JK2181, JK2183, TU144, TU148, TU839 and TU840 were processed at the University of Tübingen, with DNA extraction and pre-amplification steps undertaken in clean room facilities and post-amplification steps performed in a separate DNA laboratory. Both laboratories fulfil standards for work with ancient DNA54,55. All surfaces of tooth and bone samples were initially UV irradiated for 30 min, to minimize the potential risk of modern DNA contamination. Subsequently, DNA was extracted by applying a well-established guanidine silica-based protocol for ancient samples52. Illumina sequencing libraries were prepared by using 20 µl of DNA extract per library48; afterwards, dual barcodes (indexes) were chemically added to the prime ends of the libraries56. For the samples from Auneau (TU839 and TU840), five sequencing libraries each were prepared; for all other samples processed in Tübingen, three sequencing libraries each were prepared. To detect potential contamination of the chemicals, negative controls were conducted for extraction and library preparation. After preparation of the sequencing libraries, DNA concentration was measured with qPCR (Roche LightCycler) using corresponding primers48. The DNA concentration was given by the copy number of the DNA fragments in 1 µl of the sample.Amplification of the indexed sequencing libraries was performed using Herculase II Fusion under the following conditions: 1× Herculase II buffer, 0.4 µM IS5 primer and 0.4 µM IS6 primer48, Herculase II Fusion DNA polymerase (Agilent Technologies), 0.25 mM dNTPs (100 mM; 25 mM each dNTP) and 0.5–4 µl barcoded library as template in a total reaction volume of 100 µl. The applied amplification thermal profile was processed as follows: initial denaturation for 2 min at 95 °C; denaturation for 30 s at 95 °C, annealing for 30 s at 60 °C and elongation for 30 s at 72 °C for 3 to 20 cycles; and a final elongation step for 5 min at 72 °C. Thereafter, the amplified DNA was purified using a MinElute purification step and DNA was eluted in 20 µl TET. The concentration of the amplified DNA sequencing libraries was measured using a Bioanalyzer (Agilent Technologies) and a DNA1000 lab chip from Agilent Technologies.The sequencing libraries were sequenced on an Illumina HiSeq 4000 platform at the Max Planck Institute for Science of Human History in Jena. The samples from Auneau (TU839 and TU840) were paired-end sequenced applying 2 × 50 + 8 + 8 cycles. All other libraries prepared in Tübingen were single-end sequenced using 75 + 8 + 8 cycles.OxfordSamples AL2657, AL2541, AL2741, AL2744, AL3185, AL2350, CH1109, AL2370, AL3272 and AL3284 were processed at the dedicated ancient DNA facility at the PalaeoBARN laboratory at the University of Oxford, following methods described previously8. In brief, double-stranded libraries were constructed following the protocol in ref. 48. These libraries were sequenced on a HiSeq 2500 (AL2657, AL2541, AL2741, AL2744) or a HiSeq 4000 (AL3185, AL2350, CH1109) instrument at the Danish National Sequencing Center or on a NextSeq 550 instrument (AL2741) at the Natural History Museum of London. For samples AL2370, AL3272 and AL3284, between six and eight separate PCR reactions with unique indexes were carried out on their libraries and they were sequenced alongside samples HOV4, VAL_18A and IN18_016 on an Illumina NovaSeq 6000 lane with an S4 100-bp paired-end set-up at SciLifeLab in Stockholm.CopenhagenSamples CGG13, CGG17, CGG19, CGG20, CGG21, CGG25, CGG26, CGG27, CGG28, CGG34, Tumat1 and IRK were processed at the GLOBE Institute, University of Copenhagen. All pre-PCR work was performed in ancient DNA facilities following ancient DNA guidelines57. The details of extraction, library construction and sequencing for the samples with CGG codes are described in ref. 21, in relation to the publication of mitochondrial data from these specimens. The Tumat1 sample was processed following the exact same protocol. In brief, DNA extraction was performed using a buffer containing urea, EDTA and proteinase K50, double-stranded libraries were prepared with NEBNext DNA Sample Prep Master MixSet 2 (E6070S, New England Biolabs) and Illumina-specific adaptors48, and sequencing was performed on an Illumina HiSeq 2500 platform using 100-bp single-read chemistry. For the IRK sample, DNA was extracted from three subsamples and purified as described in ref. 21. The three DNA extracts and the purified pre-digest of one subsample were incorporated into double-stranded libraries following the BEST protocol58, with the modifications described in ref. 59, and sequenced on a BGISEQ-500 platform using 100-bp single-read chemistry.Santa CruzSamples SC19.MCJ017, SC19.MCJ015, SC19.MCJ010 and SC19.MCJ014 were processed at the UCSC Paleogenomics Lab and were provided by the Yukon Government Paleontology program. All pre-PCR work was performed in a dedicated ancient DNA facility at the University of California, Santa Cruz, following standard ancient DNA methods60. Subsamples (250–350 mg) were sent to the UCI KECK AMS facility for radiocarbon dating, and the remaining amounts were powdered in a Retsch MM400 for extraction. For each sample, ~100 mg of powder was treated with a 0.5% sodium hypochlorite solution before extraction to remove surface contaminants61 and then combined with 1 ml lysis buffer for extraction, following the protocol in ref. 52. Samples were processed in parallel with a negative control. We quantified the extracts using a Qubit 1× dsDNA HS Assay kit (Q33231) before preparing libraries. We prepared single-stranded libraries following the protocol in ref. 62 and amplified the libraries for 9–16 cycles as informed by qPCR. After amplification, we cleaned the libraries using a 1.2× SPRI bead solution and pooled them to an equimolar ratio for in-house shallow quality-control sequencing on a NextSeq 550 paired-end 75-bp run. We then sent the libraries to Fulgent Genetics for deeper sequencing on two paired-end 150-bp lanes on a HiSeq X instrument.ViennaSample HOV4 was processed at the Department of Anthropology, University of Vienna. The sample is a canine tooth, which after sequencing was determined to derive from a dhole (Cuon alpinus). DNA was extracted from its cementum using the methods described in ref. 63 with a modified incubation time of ~18 h. The library was prepared according to the protocol in ref. 48 with the modifications from ref. 64. Five separate PCR reactions with unique indexes were carried out on the library and were sequenced alongside samples VAL_18A, IN18_016, AL2242, AL2370, AL2893, AL3272 and AL3284 on an Illumina NovaSeq 6000 lane with an S4 100-bp paired-end set-up at SciLifeLab in Stockholm.An overview of all samples and their associated metadata is available in Supplementary Data 1.Genome sequence data processingFor paired-end data, read pairs were merged and adaptors were trimmed using SeqPrep (https://github.com/jstjohn/SeqPrep), discarding reads that could not be successfully merged. Reads were mapped to the dog reference genome canFam3.1 using BWA aln (v.0.7.17)65 with permissive parameters, including a disabled seed (-l 16500 -n 0.01 -o 2). Duplicates were removed by keeping only one read from any set of reads that had the same orientation, length and start and end coordinates. For sample Taimyr-1, previously published data13 were merged with newly generated data. Data from samples processed in Copenhagen were processed as described previously66 except that they were also mapped to canFam3.1. Post-mortem damage was quantified using PMDtools (v0.60)67 with the ‘–first’ and ‘–CpG’ arguments.Genotyping and integration with previously published genomesTo construct a comparative dataset for population genetic analyses, we started from a published variant call set compiling 722 modern dog, wolf and other canid genomes from multiple previous studies (NCBI BioProject accession PRJNA448733)40. To this, we added additional modern whole genomes from other studies: 4 African golden wolves and 15 Nigerian village dogs (Genome Sequence Archive (http://gsa.big.ac.cn/), accession PRJCA000335)68, 12 Scandinavian wolves (European Nucleotide Archive accession PRJEB20635)69, 9 North American wolves and coyotes (European Nucleotide Archive accession PRJNA496590)25 and 8 other canids (African hunting dog, dhole, Ethiopian wolf, golden jackal, Middle Eastern grey wolves) (European Nucleotide Archive accession PRJNA494815)22. Reads from these genomes were mapped to the dog reference genome using bwa mem (version 0.7.15)70, marked for duplicates using Picard Tools (v2.21.4) (http://broadinstitute.github.io/picard), genotyped at the sites present in the above dataset using GATK HaplotypeCaller (v3.6)71 with the ‘-gt_mode GENOTYPE_GIVEN_ALLELES’ argument and then merged into the dataset using bcftools merge (http://www.htslib.org/). The following filters were then applied to sites and genotypes across the full dataset: sites with excess heterozygosity (bcftools fill-tags ‘ExcHet’ P value 5.8×. For divergence time analyses, haploid X chromosomes from two different male genomes were combined and the point at which the inferred effective population size for this ‘pseudodiploid’ chromosome increased sharply upwards was taken to correspond to a population divergence. Results were scaled using a mutation rate of 0.4 × 10−8 mutations per site per generation13,87 (with a 25% lower rate for X-chromosome analyses) and a mean generational interval of 3 years13. For effective population size inferences, transition variants were ignored and results were scaled using a transversions-only mutation rate inferred from results on modern genomes. For more details on the MSMC2 analyses, see Supplementary Information section 3.Selection analysesSelection analysis was performed using PLINK (v1.90b5.2)88. This analysis used the 72 ancient wolf genomes and 68 modern wolf genomes (with the latter including a historical Japanese wolf genome73 treated as ancient for analysis purposes, with its age set to 200 bp). A list of the genomes used for this analysis is available in Supplementary Data 2 (“Used for selection scan” column). All SNPs, not only transversions, were used for this analysis. The age of each wolf was set as the phenotype, with values of 0 for modern wolves, and the ‘–linear’ argument was used to test for an association between SNP genotypes and age, also applying the ‘–adjust’ argument to correct P values using genomic control. The application of genomic control34 here aimed to use the magnitude of temporal allele frequency variance observed across the genome to account for what was observed from genetic drift alone given wolf demographic history. Only results for the following sets of sites were retained and included in the Manhattan plot: sites where at least 40 ancient genomes had a genotype call, sites with a minor allele frequency among the ancient wolves of ≥5% and sites that had at least 7 neighbouring sites within a 50-kb window with a P value that was at least 90% as large (on a log10 scale) as the P value of the site itself. The last ‘neighbourhood filter’ aimed to reduce false positives by requiring similar evidence across multiple nearby sites. As a P-value significance cut-off to correct for the genome-wide testing, we used 5 × 10−8, which is commonly used in genome-wide association studies in humans and also in dogs89. We excluded 15 regions where only a single variant reached significance. A detailed table with the 24 detected regions is available in Supplementary Data 3. To test the robustness of this analysis to false positives arising from genetic drift alone, we applied the same analysis to data from neutral coalescent simulations generated using ms90 and found no false positives. For more details, see Supplementary Information section 4.Ancestry modelling with qpAdm and qpWaveWe used the qpAdm and qpWave methods43 from ADMIXTOOLS (v5.0)84 to test ancestry models for wolf and dog targets postdating 23 ka. For the primary analyses, we used the following set of candidate source populations (age estimate in brackets, years bp): Armenia_Hovk1.HOV4 (ancient dhole), Siberia_UlakhanSular.LOW008 (70,772), Germany_Aufhausener.AH575 (57,233), Siberia_BungeToll.CGG29 (48,210), Germany_HohleFels.JK2183 (32,366), Siberia_BelayaGora.IN18_016 (32,020), Yukon_QuartzCreek.SC19.MCJ010 (29,943), Altai_Razboinichya.AL2744 (28,345), Siberia_BelayaGora.IN18_005 (18,148) and Germany_HohleFels.JK2179 (13,229). We used a rotating approach in which, for each target, we tested all possible one-, two- and three-source models that could be enumerated from the above set. Individuals from the set that were not used as a source in a given model served as thereference set (or the ‘right’ population in the qpAdm framework). This means that, in every model, each of the above individuals was always either in the source list or in the reference list. We ranked models on the basis of their P values, but prioritized models with fewer sources using a P-value threshold of 0.01: if a simpler model (meaning a model with fewer sources) had a P value above this threshold, it ranked above a more complex model (meaning a model with more sources) regardless of the P value of the latter. We also failed models with inferred ancestry proportions larger than 1.1 or smaller than −0.1. For single-source models, qpWave was run instead of qpAdm. Both programs were run with the ‘allsnps: YES’ option (without this option, there was very little power to reject models). We describe ancestry assigned to the ancient dhole source (Armenia_Hovk1.HOV4) as ‘unsampled’ ancestry; note that this does not imply that such ancestry is of non-wolf origin, only that it is not represented by (that is, diverged early from and lacks shared genetic drift with) the ancient wolf genomes in the reference set.To test whether any post-23 ka or modern wolf genome available might be a good proxy for the western Eurasian wolf-related ancestry identified in Near Eastern and African dogs, we added the 9,500-year-old Zhokhov dog17 to the rotating set of candidate source populations. Chosen for its high coverage, early date and easterly location, this makes the assumption that the Zhokhov dog is a good representative for the eastern dog ancestry component. Using the African Basenji dog as a target, models involving the Zhokhov dog plus another given wolf thus allowed us to test whether that wolf was a good match for the additional component of ancestry. For more details on the qpAdm and qpWave analyses, see Supplementary Information sections 2 (wolf targets) and  5 (dog targets).Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this paper. More

Misra, V. K. & Draper, D. E. On the role of magnesium ions in RNA stability. Biopolymers 48, 113–135 (1998).CAS 
PubMed 
Article 

Google Scholar 
Apell, H.-J., Hitzler, T. & Schreiber, G. Modulation of the Na, K-ATPase by magnesium ions. Biochemistry 56, 1005–1016 (2017).CAS 
PubMed 
Article 

Google Scholar 
Iseri, L. T. & French, J. H. Magnesium: Nature’s physiologic calcium blocker. Am. Heart J. 108, 188–193 (1984).CAS 
PubMed 
Article 

Google Scholar 
Rubin, H. Central role for magnesium in coordinate control of metabolism and growth in animal cells. Proc. Natl. Acad. Sci. USA 72, 3551–3555 (1975).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
de Baaij, J. H. F., Hoenderop, J. G. J. & Bindels, R. J. M. Magnesium in man: Implications for health and disease. Physiol. Rev. 95, 1–46 (2015).PubMed 
Article 
CAS 

Google Scholar 
Association, A. D. Diagnosis and classification of diabetes mellitus. Diabetes Care 37, S81–S90 (2014).Article 

Google Scholar 
Chatterjee, S., Khunti, K. & Davies, M. J. Type 2 diabetes. Lancet 389, 2239–2251 (2017).CAS 
PubMed 
Article 

Google Scholar 
Gommers, L. M. M., Hoenderop, J. G. J., Bindels, R. J. M. & de Baaij, J. H. F. Hypomagnesemia in type 2 diabetes: A vicious circle?. Diabetes 65, 3–13 (2016).CAS 
PubMed 
Article 

Google Scholar 
Pham, P.-C.T., Pham, P.-M.T., Pham, S. V., Miller, J. M. & Pham, P.-T.T. Hypomagnesemia in patients with type 2 diabetes. Clin. J. Am. Soc. Nephrol. 2, 366–373 (2007).CAS 
PubMed 
Article 

Google Scholar 
Mather, H. M. et al. Hypomagnesaemia in diabetes. Clin. Chim. Acta 95, 235–242 (1979).CAS 
PubMed 
Article 

Google Scholar 
Hubbard, S. R. Crystal structure of the activated insulin receptor tyrosine kinase in complex with peptide substrate and ATP analog. EMBO J. 16, 5572–5581 (1997).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kurstjens, S. et al. Determinants of hypomagnesemia in patients with type 2 diabetes mellitus. Eur. J. Endocrinol. 176, 11–19 (2017).CAS 
PubMed 
Article 

Google Scholar 
Viering, D. H. H. M., de Baaij, J. H. F., Walsh, S. B., Kleta, R. & Bockenhauer, D. Genetic causes of hypomagnesemia, a clinical overview. Pediatr. Nephrol. 32, 1123–1135 (2017).PubMed 
Article 

Google Scholar 
Peacock, J. M. et al. Serum magnesium and risk of sudden cardiac death in the Atherosclerosis Risk in Communities (ARIC) Study. Am. Heart J. 160, 464–470 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Veronese, N. et al. Effect of magnesium supplementation on glucose metabolism in people with or at risk of diabetes: A systematic review and meta-analysis of double-blind randomized controlled trials. Eur. J. Clin. Nutr. 70, 1354–1359 (2016).CAS 
PubMed 
Article 

Google Scholar 
Rodríguez-Morán, M., Simental-Mendía, L. E., Gamboa-Gómez, C. I. & Guerrero-Romero, F. Oral magnesium supplementation and metabolic syndrome: A randomized double-blind placebo-controlled clinical trial. Adv. Chronic Kidney Dis. 25, 261–266 (2018).PubMed 
Article 

Google Scholar 
Grigoryan, R. et al. Multi-collector ICP-mass spectrometry reveals changes in the serum Mg isotopic composition in diabetes type I patients. J. Anal. At. Spectrom. 34, 1514–1521 (2019).CAS 
Article 

Google Scholar 
Bigeleisen, J. & Mayer, M. G. Calculation of equilibrium constants for isotopic exchange reactions. J. Chem. Phys. 15, 261–267 (1947).ADS 
CAS 
Article 

Google Scholar 
Bigeleisen, J. The relative reaction velocities of isotopic molecules. J. Chem. Phys. 17, 675–678 (1949).ADS 
CAS 
Article 

Google Scholar 
McKeegan, K. D. et al. Isotopic compositions of cometary matter returned by stardust. Science 314, 1724–1728 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Jouzel, J. et al. Vostok ice core: A continuous isotope temperature record over the last climatic cycle (160,000 years). Nature 329, 403–408 (1987).ADS 
CAS 
Article 

Google Scholar 
Albarède, F., Télouk, P. & Balter, V. Medical applications of isotope metallomics. Rev. Mineral. Geochem. 82, 851–885 (2017).Article 
CAS 

Google Scholar 
Balter, V. et al. Natural variations of copper and sulfur stable isotopes in blood of hepatocellular carcinoma patients. PNAS 112, 982–985 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Télouk, P. et al. Copper isotope effect in serum of cancer patients. A pilot study. Metallomics 7, 299–308 (2015).PubMed 
Article 
CAS 

Google Scholar 
Lobo, L. et al. Elemental and isotopic analysis of oral squamous cell carcinoma tissues using sector-field and multi-collector ICP-mass spectrometry. Talanta 165, 92–97 (2017).CAS 
PubMed 
Article 

Google Scholar 
Costas-Rodríguez, M. et al. Body distribution of stable copper isotopes during the progression of cholestatic liver disease induced by common bile duct ligation in mice. Metallomics 11, 1093–1103 (2019).PubMed 
Article 
CAS 

Google Scholar 
Lamboux, A. et al. The blood copper isotopic composition is a prognostic indicator of the hepatic injury in Wilson disease. Metallomics 12, 1781–1790 (2020).CAS 
PubMed 
Article 

Google Scholar 
Moynier, F., Creech, J., Dallas, J. & Le Borgne, M. Serum and brain natural copper stable isotopes in a mouse model of Alzheimer’s disease. Sci. Rep. 9, 11894 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Sauzéat, L. et al. Isotopic evidence for disrupted copper metabolism in amyotrophic lateral sclerosis. iScience 6, 264–271 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Krayenbuehl, P.-A., Walczyk, T., Schoenberg, R., von Blanckenburg, F. & Schulthess, G. Hereditary hemochromatosis is reflected in the iron isotope composition of blood. Blood 105, 3812–3816 (2005).CAS 
PubMed 
Article 

Google Scholar 
Anoshkina, Y. et al. Iron isotopic composition of blood serum in anemia of chronic kidney disease. Metallomics 9, 517–524 (2017).CAS 
PubMed 
Article 

Google Scholar 
Morgan, J. L. L. et al. Rapidly assessing changes in bone mineral balance using natural stable calcium isotopes. Proc. Natl. Acad. Sci. USA 109, 9989–9994 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Eisenhauer, A. et al. Calcium isotope ratios in blood and urine: A new biomarker for the diagnosis of osteoporosis. Bone Rep. 10, 100200 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Isaji, Y. et al. Magnesium isotope fractionation during synthesis of chlorophyll a and bacteriochlorophyll a of benthic phototrophs in hypersaline environments. ACS Earth Space Chem. 3, 1073–1079 (2019).ADS 
CAS 
Article 

Google Scholar 
Pokharel, R. et al. Magnesium stable isotope fractionation on a cellular level explored by cyanobacteria and black fungi with implications for higher plants. Environ. Sci. Technol. 52, 12216–12224 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bolou-Bi, E. B., Poszwa, A., Leyval, C. & Vigier, N. Experimental determination of magnesium isotope fractionation during higher plant growth. Geochim. Cosmochim. Acta 74, 2523–2537 (2010).ADS 
CAS 
Article 

Google Scholar 
Wang, Y. et al. Magnesium isotope fractionation reflects plant response to magnesium deficiency in magnesium uptake and allocation: A greenhouse study with wheat. Plant Soil 455, 93–105 (2020).CAS 
Article 

Google Scholar 
Martin, J. E., Vance, D. & Balter, V. Natural variation of magnesium isotopes in mammal bones and teeth from two South African trophic chains. Geochim. Cosmochim. Acta 130, 12–20 (2014).ADS 
CAS 
Article 

Google Scholar 
Martin, J. E., Vance, D. & Balter, V. Magnesium stable isotope ecology using mammal tooth enamel. PNAS 112, 430–435 (2015).ADS 
CAS 
PubMed 
Article 

Google Scholar 
DeFronzo, R. A., Tobin, J. D. & Andres, R. Glucose clamp technique: A method for quantifying insulin secretion and resistance. Am. J. Physiol.-Endocrinol. Metab. 237, E214 (1979).CAS 
Article 

Google Scholar 
Kim, J. K. Hyperinsulinemic-euglycemic clamp to assess insulin sensitivity in vivo. In Type 2 Diabetes: Methods and Protocols, Methods in Molecular Biology (ed. Stocker, C.) 221–238 (Humana Press, 2009).Chapter 

Google Scholar 
DeFronzo, R. A., Hendler, R. & Simonson, D. Insulin resistance is a prominent feature of insulin-dependent diabetes. Diabetes 31, 795–801 (1982).CAS 
PubMed 
Article 

Google Scholar 
Balter, V. et al. Contrasting Cu, Fe, and Zn isotopic patterns in organs and body fluids of mice and sheep, with emphasis on cellular fractionation. Metallomics 5, 1470–1482 (2013).CAS 
PubMed 
Article 

Google Scholar 
Zhang, B., Podolskiy, D. I., Mariotti, M., Seravalli, J. & Gladyshev, V. N. Systematic age-, organ-, and diet-associated ionome remodeling and the development of ionomic aging clocks. Aging Cell 19, e13119 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Morel, J.-D. et al. The mouse metallomic landscape of aging and metabolism. Nat. Commun. 13, 607 (2022).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Grigoryan, R., Costas-Rodríguez, M., Vandenbroucke, R. E. & Vanhaecke, F. High-precision isotopic analysis of Mg and Ca in biological samples using multi-collector ICP-mass spectrometry after their sequential chromatographic isolation—Application to the characterization of the body distribution of Mg and Ca isotopes in mice. Anal. Chim. Acta 1130, 137–145 (2020).CAS 
PubMed 
Article 

Google Scholar 
Goff, S. L., Albalat, E., Dosseto, A., Godin, J.-P. & Balter, V. Determination of magnesium isotopic ratios of biological reference materials via multi-collector inductively coupled plasma mass spectrometry. Rapid Commun. Mass Spectrom. 35, e9074 (2021).PubMed 

Google Scholar 
DeRocher, K. A. et al. Chemical gradients in human enamel crystallites. Nature 583, 66–71 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Johansen, T., Hansen, H. S., Richelsen, B. & Malmlöf, K. The obese Göttingen minipig as a model of the metabolic syndrome: dietary effects on obesity, insulin sensitivity, and growth hormone profile. Comp. Med. 51, 150–155 (2001).CAS 
PubMed 

Google Scholar 
Coelho, P. G. et al. Effect of obesity or metabolic syndrome and diabetes on osseointegration of dental implants in a miniature swine model: A pilot study. J. Oral Maxillofac. Surg. 76, 1677–1687 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Saltiel, A. R. & Kahn, C. R. Insulin signalling and the regulation of glucose and lipid metabolism. Nature 414, 799–806 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Elin, R. J. Assessment of magnesium status. Clin. Chem. 33, 1965–1970 (1987).CAS 
PubMed 
Article 

Google Scholar 
Koopmans, S. J., van der Meulen, J., Dekker, R., Corbijn, H. & Mroz, Z. Diurnal variation in insulin-stimulated systemic glucose and amino acid utilization in pigs fed with identical meals at 12-hour intervals. Horm. Metab. Res. 38, 607–613 (2006).CAS 
PubMed 
Article 

Google Scholar 
Koopmans, S. J., Maassen, J. A., Radder, J. K. & Frölich, M. In vivo insulin responsiveness for glucose uptake and production at eu- and hyperglycemic levels in normal and diabetic rats. Biochimica et Biophysica Acta (BBA) General Subjects 1115, 230–238 (1992).CAS 
Article 

Google Scholar 
Koopmans, S. J. et al. Association of insulin resistance with hyperglycemia in streptozotocin-diabetic pigs: Effects of metformin at isoenergetic feeding in a type 2–like diabetic pig model. Metabolism 55, 960–971 (2006).CAS 
PubMed 
Article 

Google Scholar  More

Mack, R. M. et al. Biotic invasions: Causes, epidemiology, global consequences, and control. Ecol. Appl. 10(3), 689–710. https://doi.org/10.1890/1051-0761(2000)010[0689:BICEGC]2.0.CO;2 (2000).Article 

Google Scholar 
Pyšek, P. et al. A global assessment of invasive plant impacts on resident species, communitiesand ecosystems: The interaction of impact measures, invading species’ traits and environment. Glob. Change Biol. 18, 1725–1737. https://doi.org/10.1111/j.1365-2486.2011.02636.x (2012).ADS 
Article 

Google Scholar 
Wittenberg, R. & Cock, M. J. W. Invasive Alien Species: A Toolkit of Best Prevention and Management Practices (CAB International, 2001).Book 

Google Scholar 
DAISIE. Delivering Alien Invasive Species Inventories for Europe. http://www.europe-aliens.org/speciesFactsheet.do?speciesId=23539# (2018).Hejda, M., Pyšek, P. & Jarošík, V. Impact of invasive plants on the species richness, diversity and composition of invaded communities. J. Ecol. 97, 393–403. https://doi.org/10.1111/j.1365-2745.2009.01480.x (2009).Article 

Google Scholar 
Chmura, D. et al. The influence of invasive Fallopia taxa on resident plant species in two river valleys (southern Poland). Acta Soc. Bot. Pol. 84(1), 23–33. https://doi.org/10.5586/asbp.2015.008 (2015).Article 

Google Scholar 
Stefanowicz, A. M., Stanek, M., Nobis, M. & Zubek, S. Few effects of invasive plants Reynoutria japonica, Rudbeckia laciniata and Solidago gigantea on soil physical and chemical properties. Sci. Total Environ. 574, 938–946. https://doi.org/10.1016/j.scitotenv.2016.09.120 (2017).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Stefanowicz, A. M., Stanek, M., Nobis, M. & Zubek, S. Species-specific effects of plant invasions on activity, biomass and composition of soil microbial communities. Biol. Fertil. Soils 52, 841–852. https://doi.org/10.1007/s00374-016-1122-8 (2016).CAS 
Article 

Google Scholar 
Zubek, S. et al. Invasive plants affect arbuscular mycorrhizal fungi abundance and species richness as well as the performance of native plants grown in invaded soils. Biol. Fertil. Soils 52, 879–893. https://doi.org/10.1007/s00374-016-1127-3 (2016).Article 

Google Scholar 
Krinke, L. et al. Seed bank of an invasive alien, Heracleum mantegazzianum, and its seasonal dynamics. Seed Sci. Res. 15, 239–248. https://doi.org/10.1079/SSR2005214 (2005).Article 

Google Scholar 
Gioria, M. & Osbourne, B. Similarities in the impact of three large invasive plant species on soil seed bank communities. Biol. Invasions 12, 1671–1683. https://doi.org/10.1007/s10530-009-9580-7 (2010).Article 

Google Scholar 
Kundel, D., van Kleunen, M. & Dawson, W. Invasion by Solidago species has limited impacts on soil seed bank communities. Basic Appl. Ecol. 15, 573–580. https://doi.org/10.1016/j.baae.2014.08.009 (2014).Article 

Google Scholar 
Dong, H., Liu, T., Liu, Z. & Song, Z. Fate of the soil seed bank of giant ragweed and its significance in preventing and controlling its invasion in grasslands. Ecol. Evol. https://doi.org/10.1002/ece3.6238 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Harper, J. L. Population Biology of Plants (Academic Press, 1977).
Google Scholar 
Gioria, M. & Pyšek, P. The legacy of plant invasions: Changes in the soil seed bank of invaded plant communities. Bioscience 66(1), 40–53. https://doi.org/10.1093/biosci/biv165 (2015).Article 

Google Scholar 
Gioria, M. & Osborne, B. Resource competition in plant invasions: Emerging patterns and research needs. Front. Plant Sci. 5, 501. https://doi.org/10.3389/fpls.2014.00501 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Holmes, P. M. & Cowling, R. M. Diversity, composition and guild structure relationships between soil-stored seed banks and mature vegetation in alien plant-invaded South African fynbos shrublands. Plant Ecol. 133, 107–122. https://doi.org/10.1023/A:1009734026612 (1997).Article 

Google Scholar 
Gioria, M., Pyšek, P. & Moravcová, L. Soil seed banks in plant invasions: Promoting species invasiveness and long-term impact on plant community dynamics. Preslia 84, 327–350 (2012).
Google Scholar 
Tokarska-Guzik, B. et al. Rośliny Obcego Pochodzenia w Polsce ze Szczególnym Uwzględnieniem Gatunków Inwazyjnych (Generalna Dyrekcja Ochrony Środowiska, 2012).
Google Scholar 
Thompson, K., Bakker, J. P. & Bekker, R. M. The Soil Seed Banks of North West Europe: Methodology, Density and Longevity (Cambridge University Press, 1997).
Google Scholar 
Gioria, M., Le Roux, J. J., Hirsch, H., Moravcová, L. & Pyšek, P. Characteristics of the soil seed bank of invasive and non-invasive plants in their native and alien distribution range. Biol. Invasions 21, 2313–2332 (2019).Article 

Google Scholar 
Pyšek, P. et al. Naturalization of central European plants in North America: Species traits, habitats, propagule pressure, residence time. Ecology 96(3), 762–774. https://doi.org/10.1890/14-1005.1 (2015).Article 
PubMed 

Google Scholar 
Hager, H. A., Rupert, R., Quinn, L. D. & Newman, J. A. Escaped Miscanthus sacchariflorus reduces the richness and diversity of vegetation and the soil seed bank. Biol. Invasions 17, 1833–1847. https://doi.org/10.1007/s10530-014-0839-2 (2015).Article 

Google Scholar 
Robertson, S. G. & Hickman, K. Aboveground plant community and seed bank composition along an invasion gradient. Plant Ecol. 213(9), 1461–1475. https://doi.org/10.1007/s11258-012-0104-7 (2012).Article 

Google Scholar 
Fumanal, B., Gaudot, I. & Bretagnolle, F. Seed-bank dynamics in the invasive plant, Ambrosia artemisiifolia L.. Seed Sci. Res. 18(2), 101–114 (2008).Article 

Google Scholar 
Funk, J. L. et al. Keys to enhancing the value of invasion ecology research for management. Biol. Invasions 22, 2431–2445. https://doi.org/10.1007/s10530-020-02267-9 (2020).Article 

Google Scholar 
Jalas, J. Problems concerning Rudbeckia laciniata (Asteraceae) in Europe Fragmenta Floristica et Geobotanica. Supplementum 2(1), 289–297 (1993).
Google Scholar 
Tokarska-Guzik, B. The Establishment and Spread of Alien Plant Species (Kenophytes) in the Flora of Poland (Prace Naukowe Uniwersytetu Śląskiego w Katowicach, 2005).
Google Scholar 
EPPO. Rudbeckia laciniata (Asteraceae). EPPO Reporting Service—Invsive Plants. European and Mediterranean Plant Protection Organization. https://www.eppo.int/INVASIVE_PLANTS/ias_lists.htm (2009).Zelnik, I. The presence of invasive alien plant species in different habitats: Case study from Slovenia. Acta Biol. Sloven. 55(2), 25–38 (2012).
Google Scholar 
Vojniković, S. Tall cone flower (Rudbeckia laciniata L.)—new invasive species in the flora of Bosnia and Herzegovina. Herbologia 15(1), 39–47. https://doi.org/10.5644/Herb.15.1.05 (2015).Article 

Google Scholar 
Auld, B., Morita, H., Nishida, T., Ito, M. & Michael, P. Shared exotica: Plant invasions of Japan and south eastern Australia. Cunninghamia 8, 147–152 (2003).
Google Scholar 
Akasaka, M., Osawa, T. & Ikegami, M. The role of roads and urban area in occurrence of an ornamental invasive weed: A case of Rudbeckia laciniata L.. Urban Ecosyst. 18, 1021–1030 (2015).Article 

Google Scholar 
GBIF. Global Biodiversity Information Facility. Checklist dataset. https://www.gbif.org/species/3114229 (2021).Francírková, T. Contribution of the invasive ecology of Rudbeckia laciniata in the Czech Republic. In Plant Invasions: Species Ecology and Ecosystem Management (eds Brundu, G. et al.) 89–98 (Backhuys Publishers, 2001).
Google Scholar 
Moravcová, L., Pyšek, P., Jarošík, V., Havlíčková, V. & Zákravský, P. Reproductive characteristics of neophytes in the Czech Republic: Traits of invasive and non-invasive species. Preslia 82, 365–390. https://doi.org/10.1371/journal.pone.0123634 (2010).CAS 
Article 

Google Scholar 
Kościńska-Pająk, M., Musiał, K. & Janiszewska, K. Embryological processes in ovules of Rudbeckia laciniata L. (Asteraceae) from Poland. Mod. Phytomorphol. 5, 19–23 (2014).
Google Scholar 
Urbatsch, L. E. & Cox, P. B. Rudbeckia laciniata in Flora of North America Editorial Committee. http://floranorthamerica.org/Rudbeckia_laciniata (2021).Jankowska-Błaszczuk, M. Zróżnicowanie banków nasion w naturalnych i antropogenicznie przekształconych zbiorowiskach leśnych. Monograph. Bot. 88, 25 (2000).
Google Scholar 
Osawa, T. & Akasaka, M. Management of the invasive perennial herb Rudbeckia laciniata L. (Compositae) using rhizome removal. Jpn. J. Conserv. Ecol. 14(1), 37–43. https://doi.org/10.18960/hozen.14.1_37 (2009).Article 

Google Scholar 
Gleason, H. A. & Cronquist, A. Manual of Vascular Plants of Northeastern United States and Adjacent Canada (The New York Botanical Garden, 1991).Book 

Google Scholar 
Gioria, M. & Osborne, B. The impact of Gunnera tinctoria (Molina) Mirbel invasions on soil seed bank communities. J. Plant Ecol. 2(3), 153–167. https://doi.org/10.1093/jpe/rtp013 (2009).Article 

Google Scholar 
Kleyer, et al. The LEDA Traitbase: A database of life-history traits of Northwest European flora. J. Ecol. 96, 1266–1274. https://doi.org/10.1111/j.1365-2745.2008.01430.x (2008).Article 

Google Scholar 
Ruprecht, E., Fenesi, A. & Nijs, I. Are plasticity in functional traits and constancy in performance traits linked with invasiveness? An experimental test comparing invasive and naturalized plant species. Biol. Invasions 16, 1359–1372. https://doi.org/10.1007/s10530-013-0574-0 (2014).Article 

Google Scholar 
Wróbel, M. Origin and spatial distribution of roadside vegetation within the forest and agricultural areas in Szczecin Lowland (West Poland). Pol. J. Ecol. 54(1), 137–143 (2001).
Google Scholar 
Dajdok, Z. & Pawlaczyk, P. Inwazyjne Gatunki Roślin Mokradłowych Polski (Wydawnictwo Klubu Przyrodnikow, 2009).
Google Scholar 
de Waal, L. C., Child, L. E., Wade, M. & Brock, J. H. Ecology and Management of Invasive Riverside Plants (Wiley, 1994).
Google Scholar 
Pyśek, P. & Prach, K. Plant invasions and the role of riparian habitats: A comparison of four species alien to central Europe. J. Biogeogr. 20, 413–420 (1993).Article 

Google Scholar 
Kucharczyk, M. & Krawczyk, R. Kenophytes as river corridor plants in the vistula and the san river valleys. Teka Komisji Ochrony Kształtowania Środowiska Przyrodniczego 1, 110–115 (2004).
Google Scholar 
Walck, J. L. et al. Defining transient and persistent seed banks in species with pronounced seasonal dormancy and germination patterns. Seed Sci. Res. 15(3), 189–196. https://doi.org/10.1079/SSR2005209 (2005).ADS 
Article 

Google Scholar 
Gioria, M. & Pyšek, P. Early bird catches the worm: Germination as a critical step in plant invasion. Biol. Invasions 19, 1055–1080. https://doi.org/10.1007/s10530-016-1349-1 (2017).Article 

Google Scholar 
Gioria, M., Pyšek, P. & Osborne, B. Timing is everything: Does early and late germination favor invasions by herbaceous alien plants?. J. Plant Ecol. 11(1), 4–16. https://doi.org/10.1093/jpe/rtw105 (2018).Article 

Google Scholar 
Perglová, I. et al. Differences in germination and seedling establishment of alien and native Impatiens species. Preslia 81, 357–375 (2009).
Google Scholar 
Haines, D. F., Larson, D. L. & Larson, J. L. Leafy spurge (Euphorbia esula) affects vegetation more than seed banks in mixed-grass prairies of the Northern Great Plains. Invas. Plant Sci. Manage. 6, 416–432. https://doi.org/10.1614/IPSM-D-12-00076.1 (2013).Article 

Google Scholar 
Gioria, M., Jarosík, V. & Pyšek, P. Impact of invasions by alien plants on soil seed bank communities: Emerging patterns. Perspect. Plant Ecol. Evol. Syst. 16, 132–142. https://doi.org/10.1016/j.ppees.2014.03.003 (2014).Article 

Google Scholar 
Gioria, M. & Osbourne, B. Assessing the impact of plant invasions on soli seed bank communities: Use of univariate and multivariate statistical approaches. J. Veg. Sci. 20, 547–556. https://doi.org/10.1111/j.1654-1103.2009.01054.x (2009).Article 

Google Scholar 
Tokarska-Guzik, B., Bzdega, K., Knapik, D. & Jenczała, G. Changes in plant species richeness in some riparian plant communities as a result of their colonisation by taxa of Reynoutria (Fallopia). Biodivers. Res. Conserv. 1–2, 122–130 (2006).
Google Scholar 
Dölle, M. & Wolfgang, S. The relationship between soil seed bank, above-ground vegetation and disturbance intensity on old-field successional permanent plots. Appl. Veg. Sci. 12, 415–428 (2009).Article 

Google Scholar 
Thompson, K. & Grime, J. P. Seasonal variation in the seed banks of herbaceous species in ten contrasting habitats. J. Ecol. 67, 893–921. https://doi.org/10.2307/2259220 (1979).Article 

Google Scholar 
Czarnecka, J. Microspatial structure of the seed bank of xerothermic grassland—intracommunity differentiation. Acta Soc. Bot. Pol. 73(2), 155–164. https://doi.org/10.5586/asbp.2004.022 (2004).Article 

Google Scholar 
Kalamees, R., Püssa, K., Zobel, K. & Zobel, M. Restoration potential of the persistent soil seed bank in successional calcareous (alvar) grasslands in Estonia. Appl. Veg. Sci. 15, 208–218 (2012).Article 

Google Scholar 
Skowronek, S. et al. Regeneration potential of floodplain forests under the influence of nonnative tree species: Soil seed bank analysis in Northern Italy. Restor. Ecol. 22(1), 22–30. https://doi.org/10.1111/rec.12027 (2014).Article 

Google Scholar  More