Philippa Kaur

Study areasThe study took place in the cities of Basel, Lausanne, Lugano and Zurich, in Switzerland. Basel, Lausanne and Zurich are located north of the Alps, in the geographical region of the Central Plateau (Supplementary Fig. S1). This region stretches from Lake Geneva in the southwest to Lake Constance in the northeast and is the most densely populated region in Switzerland. Zurich is the largest city of Switzerland and encompasses 88 km2 with a total human resident population of 420,21741. Lausanne and Basel are smaller than Zurich, with a surface of 41 and 24 km2 and a total population of 139,408 and 173,232, respectively41. The climate in these three cities is moderately continental, with cold winters often reaching freezing temperatures in January, and warm summers. Lugano is located in Ticino, south of the Alps (Supplementary Fig. S1), where the climate is strongly affected by the Mediterranean Sea, with mild winters and summers warm and humid, sometimes hot. Lugano is the smallest of the four cities with 50,603 residents in 26 km241.Aedes albopictus is well established in Lugano since 2009 and an integrated vector management is constantly implemented to contain the numbers of the mosquito at a manageable level. This consists of an intensive surveillance, with oviposition traps distributed according to a grid system, several control interventions, such as the removal of breeding sites and the systematic application of larvicides in public areas, mainly in catch basins, and extensive public information campaigns24,26. In Basel, two populations of Ae. albopictus are established since 2018: a first population in an area adjacent to the motorway toll on the border with France and a second population in an area near the border with Germany27. The mosquito has also been recorded repeatedly at various locations in the city of Basel and the surveillance indicates that the mosquito is spreading42. Control actions are taken exclusively within the perimeter of repeated detections of the mosquito and include regular treatment of catch basins with larvicides, distribution of flyers and door-to-door information campaigns42. In Zurich, Ae. albopictus was first detected in 2016 in a bus station for international coach services located in the centre of the city, near the main train station. Thanks to immediate surveillance and control actions (i.e., treatment of catch basins in the area with larvicides), to date there is no established population within the perimeter of the bus station despite continuous repeated introductions40. A small population was also detected in 2018 in a suburban neighbourhood in the Wollishofen district of Zurich, approximately 5 km southwest from the international bus station. Also in this case, immediate surveillance and control actions, including larval control and door-to-door information, were taken with success and no adults, eggs or aquatic stages have been found in 2020 and 202140. In Lausanne, no tiger mosquito has been reported to date (Swiss Mosquito Network, http://www.mosquitoes-switzerland.ch (accessed on 17 February 2022)).Microclimate dataBased on a previous investigation we conducted in Ticino, Basel and Zurich20, we focused the microclimate monitoring on ordinary stormwater catch basins positioned on the side of public roads. In each city, we monitored ten catch basins located either in urban context (defined as areas with high-density development, consisting of apartment blocks, commercial or industrial units) or in residential areas consisting mainly of houses with private gardens located in peri-urban area (Supplementary Table S1, Supplementary Fig. S2). The catch basins were usually homogeneous in dimension, in the same city, although we recorded variations in depth. In Basel, we included catch basins located in the urban area near the border with France, in which Ae. albopictus is established. In Zurich, we included catch basins located in the international bus station, where Ae. albopictus was recorded in summer, and in the residential area of Wollishofen, where a small population of Ae. albopictus was detected and then likely eradicated. In Lausanne, some catch basins were selected in potential points of introduction of the mosquito (e.g., near a campsite, the main train station, etc.). In Lugano, Ae. albopictus was established in all the locations selected.A sensor device was installed in each selected catch basin. The sensor devices were built in house. The development of the devices and the Wireless Sensor Network (WSN) has been described in detail by Strigaro et al.29. Briefly, the device consisted of a waterproof plastic box containing a LoPy Micro-Controller Unit (Pycom, Guildford, United Kingdom), a waterproof temperature probe (accuracy of ± 0.5 °C), a light sensor (measuring illuminance arriving at the sensor device, in lux), an SD card, the rechargeable batteries and other parts. The main box, with the light sensor, was hung on the inside wall of the catch basin. The temperature probe was attached to the wall at a depth ranging from 0.3 to 0.5 m, depending on the depth of the catch basin and the level of the water in the catch basin. The probe was placed in direct contact with the inside wall of the catch basin, in order to measure the microclimatic conditions where the mosquito eggs are potentially laid. The data collected was transmitted to a data warehouse based on istSOS, an open-source Python based implementation of the Sensor Observation Service standard (SOS) of the Open Geospatial Consortium (OGC)43. The data was transmitted through the Swisscom Low Power Network (LPN) LoRaWAN (Swisscom Ltd, Ittigen, Switzerland): the data sent by the sensor devices was received by a Swisscom Gateway and then sent to the data warehouse29.In addition to the sensor devices installed in the catch basins, four devices were installed outside four catch basins in each city, except in Lugano, where three devices were installed. These external devices were placed in vegetation representing potential resting habitats for Ae. albopictus adults in the reproductive season, at 1–2 m above the ground and analyzed to confirm the close similarity between measured external temperatures and MeteoSwiss gridded temperature data. However, since the main goal of the data collection was to model the differences between MeteoSwiss gridded temperature data and catch basins’ temperatures, only a small number of external sensors were deployed. Microclimate data were collected from beginning of December 2019 to end of February 2020, a period defined as cold season, with acquisition interval set at one hour. In Lugano, data collection started on the 12th or 13th of December 2019.Local climate dataWe used two types of local climate data. The first type is the momentary hourly free-air temperatures recorded at 2 m above ground level by permanent weather stations. The weather stations belong to SwissMetNet, the automatic monitoring network of MeteoSwiss. For each city, we selected the weather station closest to the study area (Supplementary Table S1, Supplementary Fig. S2) and temperature data were retrieved from https://gate.meteoswiss.ch/idaweb (source: MeteoSwiss, Zurich-Airport, Switzerland; accessed on 12 August 2021).The second type of local climate data is the MeteoSwiss spatial climate daily datasets (source: MeteoSwiss). These temperature datasets are constructed through interpolation of daily minimum, maximum, and mean temperatures from a network of approximately 90 SwissMetNet permanent weather stations to a 1 km resolution grid in the Swiss coordinate system CH190344,45. This results in three temperature datasets describing the km-scale distribution of day-to-day temperature variations in Switzerland. We referred to them as gridded temperature data. Each monitored catch basin and external device was assigned, based on its geographical position, to the corresponding 1 km × 1 km cell of the climate grid. Each cell was identified with its MeteoSwiss (MS) number (Supplementary Table S1).Data analysisThe hourly temperatures were used to compute daily mean, maximum and minimum temperatures and daily temperature ranges, which were calculated as the difference between the maximum and minimum daily temperature. Temperatures of catch basins and external habitats were compared to temperatures of permanent weather stations and to the gridded temperatures both graphically and using the nonparametric Mann–Whitney U-test, for which a P value of  More

Ethical statementThis research project was ethically reviewed and approved by the Ethics Committee of the Instituto Gulbenkian de Ciência (license reference: A009.2018), and by the Portuguese National Entity that regulates the use of laboratory animals (DGAV – Direção Geral de Alimentação e Veterinária (license reference: 008958). All experiments conducted on animals followed the Portuguese (Decreto-Lei n° 113/2013) and European (Directive 2010/63/EU) legislations, concerning housing, husbandry and animal welfare.Escherichia coli clonesThe ancestral invader E. coli strain expresses a Yellow Fluorescent Protein (YFP), and carries streptomycin and ampicillin resistance markers for easiness of isolation from the mouse feces [galK::amp (pZ12)::PLlacO−1-YFP, strR (rpsl150), ΔlacIZYA::scar]. An E. coli strain used for the in vivo competition experiments is isogenic to the ancestral invader but expresses a Cyan Fluorescent Protein (CFP) and carries streptomycin and chloramphenicol resistance markers [galK::chlor (pZ12)::PLlacO−1-CFP, strR (rpsl150), ΔlacIZYA::scar]. The resident E. coli lineage was isolated from the feces along time using McConkey + 0.4% lactose medium, as previously described9. All the resident clones sampled from each mouse belong to E.coli phylogenetic group B9. The invader E. coli strains (YFP and CFP) derive from the K-12 MG1655 strain (DM08) and exhibit a gat negative phenotype, gatZ::IS112. The resident E. coli clone used for the competition experiments in the mouse gut expresses a mCherry fluorescent protein and a chloramphenicol resistance marker, allowing to distinguish the invader and resident strains in the mice feces.E. coli clones were grown at 37 °C under aeration in liquid media Luria broth (LB) from SIGMA — or McConkey and LB agar plates. Media were supplemented with antibiotics streptomycin (100 µg/mL), ampicillin (100 µg/mL) or chloramphenicol (30 µg/mL) when specified.Serial plating of 1X PBS dilutions of feces in LB agar plates supplemented with the appropriate antibiotics were incubated overnight and YFP, CFP or mCherry-labeled bacterial numbers were assessed by counting the fluorescent colonies using a fluorescent stereoscope (SteREO Lumar, Carl Zeiss). The detection limit for bacterial plating was ~300 CFU/g of feces9.In vivo evolution and competition experimentsAll mice (Mus musculus) used in this study were supplied by the Rodent Facility at Instituto Gulbenkian de Ciência (IGC) and were given ad libitum access to food (Rat and Mouse No.3 Breeding (Special Diets Services) and water. Mice were kept at 20-24 °C and 40-60% humidity with a 12-h light-dark cycle. For the in vivo evolution experiment we used the gut colonization model previously established9. Briefly, mice drank water with streptomycin (5 g/L) only for 24 h before a 4 h starvation period of food and water. The animals were then inoculated by gavage with 100 µL of an E. coli bacterial suspension of ~108 colony-forming units (CFUs). Mice A2, B2, D2, E2, G2, H2 and I2 were successfully colonized with the invader E. coli, while mice C2 and F2 failed to be colonized. Six- to eight-week-old C57BL/6 J non-littermate female mice were kept in individually ventilated cages under specified pathogen-free (SPF) barrier conditions at the IGC animal facility. Fecal pellets were collected during more than one year ( >400 days) and stored in 15% glycerol at −80 °C for later analysis. In the competition experiments between the invader ancestral E. coli and evolved populations, we colonized the mice using a 1:1 ratio of each genotype, with bacterial loads being assessed and frozen on a daily basis after gavage.In vivo competition experiments in which the two modes of selection (directional and diversifying) were acting for a longer time period were performed using evolved invader E. coli populations colonizing mice D2, B2 and A2, H2. Here we used both male (n = 8) and female (n = 8) C57BL/6 J mice aged six- to eight-week-old treated with streptomycin during 3 days before gavage. E. coli populations evolving for short time periods do not allow for strong conclusions on which mode of selection is taking place. Evolved invader populations such as I2 or G2 were therefore not used for in vivo fitness assays. To assess the impact of the mouse resident E. coli in the competitive fitness of dgoR we performed one-to-one competitions between the invader ancestral and dgoR KO clones. We first homogenized the mice microbiotas by co-housing the animals during seven days. The animals (n = 6, female C57BL/6 J mice aged six- to eight-week-old) were then maintained under co-housing and given streptomycin-supplemented (5 g/L) water during seven days to break colonization resistance and eradicate their resident E. coli. At this point, the co-housed mice were removed from the antibiotic-supplemented water for two days. The following day, one group of mice was gavaged with an mCherry-expressing resident E. coli (n = 3 mice) while the other group (n = 3) was not, with all animals being individually caged from this point on and receiving normal water without antibiotic. The day after gavage, all mice were colonized with a mix (1:1) of the invader ancestral and the dgoR KO clones, and the bacterial loads were assessed and frozen on a daily basis.Microbiota analysisFecal DNA was extracted with a QIAamp DNA Stool MiniKit (Qiagen), according to the manufacturer’s instructions and with an additional step of mechanical disruption32. 16 S rRNA gene amplification and sequencing was carried out at the Gene Expression Unit from Instituto Gulbenkian de Ciência, following the service protocol. For each sample, the V4 region of the 16 S rRNA gene was amplified in triplicate, using the primer pair F515/R806, under the following PCR cycling conditions: 94 °C for 3 min, 35 cycles of 94 °C for 60 s, 50 °C for 60 s, and 72 °C for 105 s, with an extension step of 72 °C for 10 min. Samples were then pair-end sequenced on an Illumina MiSeq Benchtop Sequencer, following Illumina recommendations. Sampling for microbiota analysis was performed until the microbiota composition stabilized (~1 year after the antibiotic perturbation).QIIME2 version 2017.1133 was used to analyze the 16 S rRNA sequences by following the authors’ online tutorials (https://docs.qiime2.org/2017.11/tutorials/). Briefly, the demultiplexed sequences were filtered using the “denoise-single” command of DADA2 version 1.1434, and forward and reverse sequences were trimmed in the position in which the 25th percentile’s quality score got below 20. Diversity analysis was performed following the QIIME2 tutorial35. Beta diversity distances were calculated through Unweighted Unifrac36. For taxonomic analysis, OTU were picked by assigning operational taxonomic units at 97% similarity against the Greengenes database version 13 (Greengenes 13_8 99% OTUs (250 bp, V4 region 515 F/806 R))37.Whole-genome sequencing and analysis pipelineDNA was extracted38 from E. coli populations (mixture of  > 1000 clones) or a single clone growing in LB plates supplemented with antibiotic to avoid contamination. DNA concentration and purity were quantified using Qubit and NanoDrop, respectively. The DNA library construction and sequencing were carried out by the IGC genomics facility using the Illumina Miseq platform. Processing of raw reads and variants analysis was based on the previous work39. Briefly, sequencing adapters were removed using fastp version 0.20.040 and raw reads were trimmed bidirectionally by 4 bp window sizes across which an average base quality of 20 was required to be retained. Further retention of reads required a minimum length of 100 bps per read containing at least 50% base pairs with phred scores at or above 20. BBsplit (part of BBMap version 38.9)41 was used to remove likely contaminating reads as explained previously39. Separate reference genomes were used for the alignment of invader (K-12 (substrain MG1655; Accession Number: NC_000913.2)) and resident (Accession Number: SAMN15163749) E. coli genomes. Alignments were performed via three alignment approaches: BWA-sampe version 0.7.1742, MOSAIK version 2.743, and Breseq version 0.35.144,45. Final average alignment depths for invader and resident populations across time points equalled 302 (median = 236) and 253 (median = 235), respectively. While Breseq provides variant analysis in addition to alignment, other variant calling approaches were used to identify putative variation in the sequenced genomes, and to verify data from Breseq. A naïve pipeline39 using the mpileup utility within SAMtools version 1.946 and a custom script written in python was employed. Only reads with a minimum mapping quality of 20 were considered for analysis, and variant calling was limited to bases with call qualities of at least 30. At these positions, a minimum of 5 quality reads had to support a putative variant on both strands (with strand bias, pos. strand / neg. strand, above 0.2 or below 5) for further consideration. Finally, mutations were retained if detected in more than one of the alignment approaches, and if they reached a minimum frequency of 5% at a minimum of one time point sampled. Further simple and complex small variants were considered from freebayes version 0.9.2147 with similar thresholds, while insertion sequence movements and other mobile element activity was inferred via is mapper version 248 and panISa version 0.1.649, as well as Breseq, as previously described39. All putative variants were verified manually in IGV version 2.750,51. Raw sequencing reads were deposited in the sequence read archive under bioproject PRJNA666769. Population dynamics of lineage-specific dynamics and the resulting Muller plots were inferred manually and are meant strictly as a means of presenting the data. In order to generate these plots, mutations were sorted by frequency (descending for each time point at which the population was sampled). The largest frequency mutations were considered major lineages within which minor frequency mutations occurred. Assuming that a mutation, which arises subsequent to a preexisting mutation (an already differentiated lineage) cannot exceed the frequency of that preexisting mutation at any point, and will fluctuate in frequency with the preexisting one, we assigned mutations to the lineages within each population. While this resolved the majority of high frequency and medium frequency mutations, low-frequency mutations within the Muller plots cannot be placed with high confidence, and are only included for completeness.Prophage induction rateTo calculate the maximum prophage induction rate we grew E. coli lysogenic clones, starting with the same initial OD600 values: ~0.1 (Bioscreen C system, Oy Growth Curves Ab Ltd), with agitation at 37 °C in LB medium in the presence or absence of mitomycin C along time (5 µg/mL)9. The OD600 values were normalized by dividing the ones in the presence of mitomycin C by the ones in the absence of mitomycin C (sampling interval: 30 min). The LN of this ratios along time originates a lysis curve, where the maximum slope corresponds to the maximal prophage induction rate for each clone analyzed. We tested evolved clones from mouse A2, H2 and G2 against the ancestral clone which only carries the Nef and the KingRac prophages. We also tested clones of the resident strain that had evolved in the presence of the invader for more than 400 days (these clones were sampled from mouse A2).
E. coli growth rate, growth curves, cell aggregation, biofilm and motility capacityTo calculate the maximum bacterial growth rate, we grew E. coli lysogenic clones, starting with the same initial OD600 values: ~0.1 (Bioscreen C system, Oy Growth Curves Ab Ltd), with agitation at 37 °C in LB medium along time using reading intervals of 30 min. The LN of the OD600 values along time originates a growth curve, where the maximum slope corresponds to the maximum bacterial growth rate for each clone analyzed.To test for metabolic differences of the psuK/fruA mutation, growth curves of evolved lysogenic E. coli clones, bearing the Nef and KingRac prophages, with or without the psuK/fruA mutation were performed with the same initial OD600 value (~0.03) for each clone. The clones were grown in glucose (0.4%) minimal medium (MM9-SIGMA) with or without pseudouridine (80 μM) and absorbance values were obtained using the Bioscreen C apparatus during 12 h.Frozen stocks of E. coli clones were used to seed tubes with 5 mL of liquid LB. These were incubated overnight at 37 °C under static conditions to assess the formation of cell flocks/clumps, observable to the naked eye, in order to evaluate the formation of cell aggregates. Biofilm was tested according a previously published protocol52 and to evaluate the motility capacity we adapted the protocol from Croze and colleagues53. Briefly, overnight E. coli clonal cultures grown with agitation at 37 °C in 5 mL LB medium supplemented with streptomycin (100 ug/mL) were adjusted to the same absorbance and a 3uL volume was dropped on top of soft agar (0.25%). Plates were incubated at 37 °C and photos were taken at day 1, 2 and 5 post-inoculation to assess swarming motility phenotype.Number of E. coli generations during mouse gut colonizationTo estimate the number of generations of E. coli in the mouse gut, we used a previously described protocol to measure the fluorescent intensity of a probe specific to E. coli 23 S rRNA (as a measure of ribosomal content) that correlates with the growth rate of the bacterial cells54. We measured the number of generations of the ancestral E. coli clone while colonizing the gut of 2 mice, treated during 24 h with streptomycin (5 g/L) before gavage, during 25 days.Plasmid DNA extraction and PCR detection of ~69Kb (repA) and ~109Kb (repB) plasmidsPlasmid DNA was extracted from overnight cultures using a Plasmid Mini Kit (Qiagen), according to the manufacturer’s guidelines. Specific primers for the amplification of repA and repB genes, were used to determine the frequency of the 68935 bp (~69 Kb) and 108557 bp (~109 Kb) plasmids, respectively, in the invader E. coli population.The primers used for repA gene were:repA-Forward: 5’-CAGTCCCCTAAAGAATCGCCCC-3’ and repA-Reverse: 5’-TGACCAGGAGCGGCACAATCGC-3’.For repB the primer sequences were:repB-Forward: 5’-GTGGATAAGTCGTCCGGTGAGC-3’ and repB-Reverse: 5’-GTTCAAACAGGCGGGGATCGGC3’.PCR amplification of plasmid-specific genes was performed in 12 isolated random clones from mouse A2 at days 104 and 493. PCR reactions were performed in a total volume of 25 μL, containing 1 μL of plasmid DNA, 1X Taq polymerase buffer, 200 μM dNTPs, 0.2 μM of each primer and 1.25 U Taq polymerase. PCR reaction conditions: 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 65 °C for 30 s and 72 °C for 30 s, finalizing with 5 min at 72 °C. DNA was visualized on a 2% agarose gel stained with GelRed and run at 160 V for 60 min.Construction of the dgoR KO mutantP1 transduction was used to construct a ΔdgoR mutant (dgoR KO). This KO strain was created by replacing the wild-type dgoR in the invader ancestral YFP-expressing genetic background by the respective knock-out from the KEIO collection, strain JW562755, in which the dgoR sequence is replaced by a kanamycin resistance cassette. The presence of the cassette was confirmed by PCR using primers dgoK-F: GCGATGTAGCGAGCTGTC, and yidX-R: GGGAATAAACCGGCAGCC. PCR reactions were performed in a total volume of 25 μL, containing 1 μL of DNA, 1X Taq polymerase buffer, 200 μM dNTPs, 0.2 μM of each primer and 1.25 U Taq polymerase. PCR reaction conditions: 95 °C for 3 min, followed by 35 cycles of 95 °C for 30 s, 65 °C for 30 s and 72 °C for 30 s, finalizing with 5 min at 72 °C. DNA was visualized in a 2% agarose gel stained with GelRed and run at 160 V for 60 min.RNA extraction, DNAse treatment, RT-PCR and qPCRThe Qiagen RNeasy Mini Kit was used for RNA extraction. RNA concentration and quality were evaluated in the Nanodrop 2000 and by gel-electrophoresis. DNase treatment was performed with the RQ1 DNase (Promega) by adding 0.5 μl of DNase to 1 μg of RNA and 1 μl buffer in a final volume of 15 ul, followed by incubation 30 min at 37 °C. Afterwards, 1 ul of stop solution was added and incubation for 15 min at 65 °C was performed to inactivate the DNase. As a control for complete DNA digest a PCR was performed on the reactions including positive controls. Reverse transcription was performed with M-MLV RT[-H] (Promega) by mixing 1 μg of RNA with 0.5 μl random primers (Promega) and nuclease free water to a volume of 15 μl, incubation at 70 °C for 5 min and a quick cool down on ice. Afterwards the reverse transcription was accomplished by adding 5 μl of RT buffer, 0.5 μl RT enzyme and 2 μl dNTP mix, followed by incubation for 10 min at 25 °C, 50 min at 50 °C and 10 min at 70 °C. The resulting cDNA was diluted 100-fold in nuclease free water before changes in gene expression were detected using the The QuantStudio 7Flex (Applied Biosystems) with iTaq Universal SYBR Green Supermix (BioRad) and the following cycling protocol: Hold stage: 2 min at 50 °C, 10 min at 95 °C. PCR stage (40 cycles): 15 s at 95 °C, 30 s at 58 °C, 30 s at 60 °C. Melt curve stage: 15 s at 95 °C, 1 min at 50 °C then increments of 0.05 °C/s until 95 °C. Melt curve analysis was performed to verify product homogeneity. All reactions included six biological and three technical replicates for each sample. A relative quantification method of analysis with normalization against the endogenous control rrsA and employing the primer specific efficiencies was used according to the Pfaffl method (add reference). The primers used were designed with PrimerQuest (idt). The used primer sequences were: psuK – TGCGTTAGCAGCGATTGA, AATTTACGCCTGGTGGAGTAG; arcA – GATTCATGGTACGGGACAGTAG, CCGTGACAACGAAGTCGATAA; yjtD – CGCACATGGATCTGGTGATA, GGCGTGGCGTAGTAATGATA and rrsR – GTCAGCTCGTGTTGTGAAATG, CCCACCTTCCTCCAGTTTATC.Statistics and reproducibilityCorrelation between microbiota diversity measures and E. coli loads (CFU) or persistence (1-presence or 0-absence) was performed in R using the statistical package rmcorr (version 0.5.2)56 and lme4 (version 1.1-10)57, respectively. The rate of accumulation of new ISs in vivo was compared using Wilcoxon paired signed ranked test for expected and observed insertions, while the rate of selective sweeps correlation was performed using the Spearman Correlation test. Selective sweeps were taken to be mutations or HGT events that reached  > 95% frequency in the population and kept high frequency until the end of the colonization. Statistical analysis of prophage induction as well as biofilm levels was performed using the Mann-Whitney test in GraphPad Prism (version 8.4.3). A single sample T-Test was used test if the growth rate of evolved invader clones deviates from the mean of the ancestral. A Wilcoxon rank sum test with continuity correction was used to compare the relative expression levels of the evolved clones with the ancestral. P values of x0 are of order of their inverse selection coefficient (up to logarithmic corrections):$${{{{{rm{G}}}}}}left(xright)approx 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{s}.$$
(1)

Clonal interference under uniform directional selection. This mode occurs in asexual populations when adaptive mutations become frequent enough to interfere with one another59,60,61. Only a fraction of the established adaptive mutations reaches fixation; sojourn times to intermediate frequencies are set by a global coalescence rate (widetilde{sigma }) that is higher than the typical selection coefficient of individual mutations62:

$${{{{{rm{G}}}}}}left(xright) < 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{widetilde{sigma}}.$$ (2) Details of these dynamics depend on the spectrum of selection coefficients and on the overall mutation rate, which set the strength of clonal interference. For moderate interference, where a few concurrent beneficial mutations compete for fixation, we expect a roughly exponential drop of the frequency propagator, (Gleft(xright)sim {{exp }}left(-lambda xright)), reflecting the probability that a trajectory reaches frequency x without interference by a stronger competing clade. Moderate interference generates an effective neutrality for weaker beneficial mutations and at higher frequencies63. This regime has been mapped for influenza64. In the asymptotic regime of a travelling fitness wave, where many beneficial mutations are simultaneously present, the fate of a mutation is settled in the range of small frequencies; that is, at the tip of the wave65. In this regime, emergent neutrality affects the vast majority of beneficial mutations and most of the frequency regime66. Hence, the frequency propagator rapidly drops to its asymptotic value (Gleft(x=1right)ll 1.) Adaptation under diversifying selection. More complex selection scenarios involve selection within and between ecotypes, i.e., subpopulations occupying distinct ecological niches67,68. An important factor generating niches and ecotypes is the differential use of food and other environmental resources. In this mode, ecotype-specific, conditionally beneficial mutations reach intermediate frequencies after a time given by their within-ecotype selection coefficients, but fixation can be slowed down or suppressed by diversifying (negative frequency-dependent) cross-ecotype selection18, $${{{{{rm{G}}}}}}left(xright)approx 1,{{{{{rm{T}}}}}}left(xright)sim frac{1}{s},left(xlesssim ,frac{1}{2}right)$$ (3) $${{{{{rm{G}}}}}}left(xright) < 1,{{{{{rm{T}}}}}}left(xright)gg frac{1}{s}left(xto 1right).$$ (4) The details depend on the details of the eco-evolutionary model (synergistic vs. antagonistic interactions, carrying capacities, amount of resource competition vs. explicitly frequency-dependent selection). In a model with directional selection within ecotypes, conditionally beneficial mutations rapidly fix within ecotypes, but lead only to finite shifts of the ecotype frequencies. In the simplest case, the resulting dynamics of ecotype frequencies is diffusive, resulting in an effectively neutral turnover of ecotypes18. Given negative frequency-dependent selection between ecotypes, fixations become even rarer and can be completely suppressed; that is, ecotypes can become stable on the time scales of observation. The separation of time and selection scales between intra- and cross-ecotype frequency changes is expected to be a robust feature of ecotype-dependent selection: sojourn of adaptive alleles to intermediate frequencies is fast, fixation is slower and rarer. In other words, ecotype-dependent selection is characterized by two regimes of coalescence times T(x).Frequency propagators and the coalescence time spectra expected under these evolutionary modes are qualitatively sketched in Supplementary Fig. 11. For periodic sweeps under directional selection (dark green, left column), G(x) depends weakly on x and T(x) is set by rapid sweeps for all x. For clonal interference under directional selection (green, center column), G(x) decreases substantially with increasing x and T(x) becomes uniformly shorter. Under negative frequency-dependent selection (brown, right column), G(x) decreases substantially with increasing x, while T(x) substantially increases for large x and diverges in case of strong frequency-dependent selection generating stable ecotypes (dashed lines). (see Supplementary Fig. 11 for the results of simulations assuming a model of direction selection or assuming a resource competition model where ecotype formation occurs31.The ({{{{{boldsymbol{p}}}}}})-({{{{{boldsymbol{tau }}}}}}) selection testThis test is based on qualitative characteristics of the functions G(x), T(x) and does not depend on details of the evolutionary process. We evaluate G(x) and T(x) for host-specific families of frequency trajectories; sojourn times are counted from an initial frequency x0=0.01. Origination times at this frequency are inferred by backward extrapolation of the first observed trajectory segment; the reported results are robust under variations of the threshold x0 and the extrapolation procedure. We then compute two summary statistics: the probability (p) that a mutation established at an intermediate frequency xm reaches near-fixation at a frequency xf,$$p=frac{{{{{{rm{G}}}}}}({x}_{f})}{{{{{{rm{G}}}}}}({x}_{m})},$$ (5) and the corresponding fraction of sojourn times,$$tau=,frac{{{{{{rm{T}}}}}}({x}_{f})}{{{{{{rm{T}}}}}}({x}_{m})}.$$ (6) Here we use xm=0.3 and xf=0.95 to limit the uncertainties of empirical trajectories at low and high frequency; however, the selection test is robust under variation of these frequencies. We find evidence for different modes of evolution: The long-term frequency trajectories of mice B2, D2 and E2 are consistent with predominantly frequency-dependent selection (Fig. 2, Fig. 4a–c). The propagator G(x) is a strongly decreasing function of x, resulting in fixation probabilities (p) 0.6, as measured by time ratios τ  > 3.

The trajectories of mice A2, G2, and I2 show a signature of recurrent selective sweeps and clonal interference under uniform directional selection (Fig. 4a–c). The propagator G(x) is a decreasing function of x, resulting in fixation probabilities (p=0.2-0.8), depending on the strength of clonal interference. Fixation times are short, giving time ratios (tau lesssim 2).

The shorter trajectory of mouse H2 signals periodic sweeps under uniform directional selection (Fig. 3, Fig. 4a–c). The origination rate of mutations is lower than in the longer trajectories, and G(x) shows a weak decrease with (p=1.) Sojourn times T(x) are short and grow uniformly with x, resulting in a time ratio τ=2.25. This pattern is expected under directional selection in the low mutation regime: (Tleft(xright)={{log }}left[x/(1-x)right]/{s}) for individual mutations with a uniform selection coefficient s, leading to τ=2.0 for xm=0.3 and xf=0.95 (this value is marked as a dashed line in Fig. 4c).

The trajectories of non-mutator lines in the long-term in vitro evolution experiment of Good et al1, evaluated over the first 7500 generations, show an overall signal of clonal interference under uniform directional selection (Fig. 4c, Supplementary Fig. 12). The frequency propagators G(x) are strongly decreasing functions of x and sojourn times T(x) grow uniformly with x. We find (p=0.2-0.8) and (tau lesssim 2), similar to the pattern in mice A2, G2, and I2.

Code for Selection testsThe code for selection tests from the mutation frequency trajectories can be found in the Supplementary Information file.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Study site and animalsWe studied n = 64 flocks of Japanese tits in mixed deciduous-coniferous forests in Nagano and Gumma (36°17-31’N, 138°26-39’E), Japan. Although most of the birds had not been individually colour-ringed, all the experimental trials were conducted at least 400 m apart; previous observations on colour-ringed individuals showed that this distance was enough to ensure the collection of data from different individuals30. In this site, one of the major predators of small birds is the bull-headed shrike, which is often mobbed by small birds including Japanese tits.Playback stimulusTo test whether Japanese tits recognize an alert-recruitment call sequence as a single unit, we prepared four treatments: (i) one-speaker playback of alert-recruitment call sequences, (ii) two-speaker playback of alert-recruitment call sequences with alert and recruitment calls played from different speakers, (iii) one-speaker playback of recruitment-alert call sequences, (iv) two-speaker playback of recruitment-alert call sequences with recruitment and alert calls played from different speakers (Fig. 3). We created sound files for these treatments using the software program Audacity 2.1.3 (http://www.audacityteam.org). For one-speaker treatments, we composed mono sound files where call sequences were repeated onto a single channel, whereas for two-speaker treatments, we composed stereo sound files where either alert or recruitment calls were repeated onto the right or left channels, respectively. All the files contained an equal number of alert calls (30 calls) and recruitment calls (30 calls) at the same rate (one call every 3 s), resulting in 90-s of stimuli (Fig. 3), which corresponds to the range of the natural calling rate of alert-recruitment sequences during mobbing by Japanese tits10. For all stimuli, within-call-sequence intervals between alert and recruitment calls were constant (0.1 s), which is within the range of intervals of these calls in natural call sequences11,17. In contrast, between-call-sequence intervals varied from 1.50 to 1.81 (median = 1.68) due to the difference in call length, but were constant across playback stimuli within the same “block” where the four treatments were created using the same call exemplars (see below). While alert calls are composed of three distinct note types, recruitment calls are strings of the same note type that vary in repetition number. Since the repetition number can vary depending on predator type10, we conducted predator exposure experiments to Japanese tit flocks (n = 12) and recorded call sequences towards a bull-headed shrike life-like specimen. In response to a shrike specimen, tits produced alert-recruitment call sequences with a recruitment note repetition number ranging from 5 to 15. Since the interquartile range of repetition number was 6.75 to 10, we used recruitment calls with 7–10 notes as playback stimuli in this study. In consideration for the possible influence of sound editing procedure, we created all the stimuli in the same manner; we copied alert and recruitment call parts separately from recording files, and pasted them onto background noise files to produce all four types of stimuli. Playback amplitudes were constant across treatments, 70 dB at 1.0 m measured using a sound level meter (SM-325, AS ONE Corporation). Therefore, the differences between treatments only depend on whether these calls are produced as sequences from the same source and how the calls are ordered.We carefully designed experiments to control for the possibility that individual-based acoustic features in alert and recruitment calls might influence tits’ responses. First, we prepared 16 unique sets of alert and recruitment calls using either calls from the same bird (n = 8 source individuals, n = 8 unique call sets) or from two different birds (n = 16 source individuals, n = 8 unique call sets). Then, we created the four types of treatments (i.e., alert-recruitment call sequences from the same speaker, from different speakers, and in reversed order from the same speaker and from different speakers) from each of the alert-recruitment call sets, resulting in 16 blocks of playback stimuli (Supplementary Table 3). This allows us to test the possible influence of individual-based acoustic variation on receivers’ responses.We were also careful to avoid the possible influence of population-level signatures of acoustic features: we only used Japanese tits’ call sequences that had been previously recorded from the same study population. We saved the sound files in .wav format (16-bit accuracy, 48-kHz sampling rate) onto a playback device (iPhone 8, Apple Inc.). We used the default Music app (Apple Inc.) to playback the sound files.ExperimentWe (TNS and YKM) conducted experimental trials from 26 October to 4 December 2020 and during the period of 0800 and 1600 h (Japan Standard Time). We did not conduct trials under wet and windy weather conditions, since these may influence behavioural patterns of forest birds31. First, we searched for and located a flock of Japanese tits. Upon finding a flock, we fixed a taxidermic specimen of bull-headed shrike in a perching posture on the branch at 1.8 ± 0.2 m (mean ± s.d., n = 64) above the ground. Then, we placed either one or two Bluetooth speakers (SoundLink Micro, BOSE) on tree branches at 1.6 ± 0.2 m (mean ± s.d., n = 96) above the ground, and oriented them upwards to control for the possible influence of directionality. We set the distance between the shrike specimen and the speaker(s) at 5 m. For trials with two speakers, we set the distance between speakers at 10 m, mimicking the situation in which two birds are calling (Fig. 3). The shrike specimen was first covered with a black cloth and was exposed by removing the cloth just before each trial.We began playbacks when at least two Japanese tits were present within 15 m from the shrike specimen. During 90-s of playbacks, we recorded (i) whether birds approached within 2-m of the shrike specimen during the playback and (ii) whether birds exhibited wing flicking displays12,13. We counted the number of individuals within 15 m from the shrike during 90-s of playbacks and considered it as flock size. During trials, we sat on the ground at ca. 10 m from the shrike specimen to decrease the influence of the observers’ presence on bird behaviour. To account for the inter-observer reliability32, we calculated intra-class correlation coefficient (ICC; icc function in the R package irr) between us. The lowest ICC was 0.998, indicating high degree of inter-observer reliability for the two behavioural measurements. We also video-recorded the responses of tits using a digital video camera (FDR-AX60, SONY). After completion of each trial, we checked the video recording and made an on-the-spot confirmation of the exact location at which each bird made the closest approach to the shrike specimen during the 90-s of playbacks. Then, using a tape measure, we recorded the minimum approach distance of birds to the shrike specimen. Thus, our final data set consisted of the most reliable observations confirmed by two experimenters and video evidence.The order of trials was randomized within each block (n = 16 blocks), each of which is composed of a unique alert-recruitment call set but includes four treatments differing in the number of speakers and call order. Therefore, responses to all four treatments were observed under largely similar conditions. In a few trials, the first bird to approach the shrike specimen was from a heterospecific species, such as a varied tit (n = 1) or a long-tailed tit (n = 1). To account for the possibility that these birds evoke mobbing behaviour in Japanese tits, we only used the data from instances where the first individual to approach the shrike was a Japanese tit. Otherwise, we repeated the same treatment at a different site.We used 64 unique playbacks created from 16 unique sets of alert-recruitment calls for 64 trials in order to avoid pseudoreplication33. We prepared two specimens of male bull-headed shrikes and used each of them for the equal number of trials. We did not use specimens of female shrikes since females migrate from the study site in late summer and only males were observed during the study period.Statistical analysisWe analyzed the effect of playback treatments on the mobbing behaviours of Japanese tits using generalized linear mixed models in R34,35. We used the proportions of Japanese tits in flocks that (i) approached within 2-m of the shrike specimen and (ii) exhibited wing flicking displays. For the analysis of predator approach, we prepared two vectors (i.e., the number of Japanese tits that approached the shrike specimen and the number of Japanese tits that did not approach the shrike specimen). Then, we created a single response variable by binding together these two vectors using cbind function. Similarly, for the analysis of wing flicking displays, we created a single response variable by binding two vectors (i.e., the number of tits that exhibit wing flicking and the number of tits that did not exhibit wing flicking). We fitted playback treatments as a fixed term, and flock size (maximum number of Japanese tits observed during 90-s of playback) and the way of creating playback stimuli (whether the two call types were recorded from a single individual or two individuals) as covariates. We also included identity of alert-recruitment call sets that were used for creating playback stimuli (i.e., call sets from either one or two source individuals) and identity of shrike specimens as random terms. We used a binomial error distribution and logit-link function (glmer in the R package lme4) for these models. Statistical significance was calculated by log-likelihood ratio tests using anova in the R package stats. We further conducted post-hoc pairwise comparisons between treatments by using estimated marginal means (emmeans in the R package emmeans). When making pairwise comparisons, we adjusted p-values by applying a false discovery rate control for multiple testing36. All tests were two-sided and the significance level was set at α = 0.05. Exact p-values are reported when p ≥ 0.0001.Ethics statementAll protocols were approved by the ethics committee of Kyoto University, the Ministry of the Environment, and the Forestry Agency of Japan, and adhered to Guidelines for the Use of Animals of the Association for the Study of Animal Behaviour/Animal Behavior Society37.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

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).ADS 
PubMed 
Article 

Google Scholar 
Bruno, J. F. & Selig, E. R. Regional decline of coral cover in the Indo-Pacific: Timing, extent, and subregional comparisons. PLoS ONE 2, e711 (2007).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
De’Ath, G., Fabricius, K. E., Sweatman, H. & Puotinen, M. The 27-year decline of coral cover on the Great Barrier Reef and its causes. Proc. Natl. Acad. Sci. U. S. A. 109, 17995–17999 (2012).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hughes, T. P., Graham, N. A. J., Jackson, J. B. C., Mumby, P. J. & Steneck, R. S. Rising to the challenge of sustaining coral reef resilience. Trends Ecol. Evol. 25, 633–642 (2010).PubMed 
Article 

Google Scholar 
Pandolfi, J. M., Connolly, S. R., Marshall, D. J. & Cohen, A. L. Projecting coral reef futures under global warming and ocean acidification. Science 333, 418–422 (2011).ADS 
PubMed 
Article 

Google Scholar 
Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377 (2017).ADS 
PubMed 
Article 

Google Scholar 
Bindoff, N. L. et al. Chapter 5: Changing Ocean, Marine Ecosystems, and Dependent Communities. IPCC Special Report on the Ocean and Cryosphere in a Changing Climate (2019).Richmond, R. H. Reproduction and recruitment in corals: Critical links in the persistence of reefs. In Life and Death of Coral Reefs (ed. Birkeland, C. E.) 175–197 (Springer, 1997).Chapter 

Google Scholar 
Trapon, M. L., Pratchett, M. S., Hoey, A. S. & Baird, A. H. Influence of fish grazing and sedimentation on the early post-settlement survival of the tabular coral Acropora cytherea. Coral Reefs 32, 1051–1059 (2013).ADS 
Article 

Google Scholar 
Gallagher, C. & Doropoulos, C. Spatial refugia mediate juvenile coral survival during coral–predator interactions. Coral Reefs 36, 51–61 (2017).ADS 
Article 

Google Scholar 
Vermeij, M. J. A. & Sandin, S. A. Density-dependent settlement and mortality structure the earliest life phases of a coral population. Ecology 89, 1994–2004 (2008).PubMed 
Article 

Google Scholar 
Vermeij, M. J. A., Smith, J. E., Smith, C. M., Vega Thurber, R. & Sandin, S. A. Survival and settlement success of coral planulae: Independent and synergistic effects of macroalgae and microbes. Oecologia 159, 325–336 (2009).ADS 
PubMed 
Article 

Google Scholar 
Ricardo, G. F., Jones, R. J., Nordborg, M. & Negri, A. P. Settlement patterns of the coral Acropora millepora on sediment-laden surfaces. Sci. Total Environ. 609, 277–288 (2017).ADS 
PubMed 
Article 

Google Scholar 
Brunner, C. A., Uthicke, S., Ricardo, G. F., Hoogenboom, M. O. & Negri, A. P. Climate change doubles sedimentation-induced coral recruit mortality. Sci. Total Environ. 768, 143897 (2021).ADS 
PubMed 
Article 

Google Scholar 
Birrell, C. L., McCook, L. J., Willis, B. L. & Diaz-Pulido, G. A. Effects of benthic algae on the replenishment of corals and the implications for the resilience of coral reefs. In Oceanography and Marine Biology: An Annual Review 25–63 (CRC Press, 2008).Chapter 

Google Scholar 
Karcher, D. B. et al. Nitrogen eutrophication particularly promotes turf algae in coral reefs of the central Red Sea. PeerJ 2020, 1–25 (2020).
Google Scholar 
Kirschner, C. M. & Brennan, A. B. Bio-inspired antifouling strategies. Annu. Rev. Mater. Res. 42, 211–229 (2012).ADS 
Article 

Google Scholar 
Webster, N. S. et al. Metamorphosis of a scleractinian coral in response to microbial biofilms. Appl. Environ. Microbiol. 70, 1213–1221 (2004).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Heyward, A. J. & Negri, A. P. Natural inducers for coral larval metamorphosis. Coral Reefs 18, 273–279 (1999).Article 

Google Scholar 
Negri, A. P., Webster, N. S., Hill, R. T. & Heyward, A. J. Metamorphosis of broadcast spawning corals in response to bacteria isolated from crustose algae. Mar. Ecol. Prog. Ser. 223, 121–131 (2001).ADS 
Article 

Google Scholar 
Tebben, J. et al. Induction of larval metamorphosis of the coral Acropora millepora by tetrabromopyrrole isolated from a Pseudoalteromonas bacterium. PLoS ONE 6, 1–8 (2011).Article 

Google Scholar 
Sneed, J. M., Sharp, K. H., Ritchie, K. B. & Paul, V. J. The chemical cue tetrabromopyrrole from a biofilm bacterium induces settlement of multiple Caribbean corals. Proc. R. Soc. B Biol. Sci. 281, 1–9 (2014).
Google Scholar 
Tebben, J. et al. Chemical mediation of coral larval settlement by crustose coralline algae. Sci. Rep. 5, 1–11 (2015).Article 

Google Scholar 
Carpenter, R. C. & Edmunds, P. J. Local and regional scale recovery of Diadema promotes recruitment of scleractinian corals. Ecol. Lett. 9, 268–277 (2006).Article 

Google Scholar 
Box, S. J. & Mumby, P. J. Effect of macroalgal competition on growth and survival of juvenile Caribbean corals. Mar. Ecol. Prog. Ser. 342, 139–149 (2007).ADS 
Article 

Google Scholar 
Linares, C., Cebrian, E. & Coma, R. Effects of turf algae on recruitment and juvenile survival of gorgonian corals. Mar. Ecol. Prog. Ser. 452, 81–88 (2012).ADS 
Article 

Google Scholar 
McCook, L. J., Jompa, J. & Diaz-Pulido, G. Competition between corals and algae on coral reefs: A review of evidence and mechanisms. Coral Reefs 19, 400–417 (2001).ADS 
Article 

Google Scholar 
Nugues, M. M., Smith, G. W., Van Hooidonk, R. J., Seabra, M. I. & Bak, R. P. M. Algal contact as a trigger for coral disease. Ecol. Lett. 7, 919–923 (2004).Article 

Google Scholar 
Fong, J. et al. Allelopathic effects of macroalgae on Pocillopora acuta coral larvae. Mar. Environ. Res. 151, 104745. https://doi.org/10.1016/j.marenvres.2019.06.007 (2019).Article 
PubMed 

Google Scholar 
Hauri, C., Fabricius, K. E., Schaffelke, B. & Humphrey, C. Chemical and physical environmental conditions underneath mat- and canopy-forming macroalgae, and their effects on understorey corals. PLoS ONE 5, 1–9 (2010).Article 

Google Scholar 
Bay, L. K. et al. Reef Restoration and Adaptation Program : Intervention Technical Summary. A report provided to the Australian Government by the Reef Restoration and Adaptation Program. (2019).Anthony, K. R. N. et al. Interventions to help coral reefs under global change—A complex decision challenge. PLoS ONE 15, 1–14 (2020).Article 

Google Scholar 
Vardi, T. et al. Six priorities to advance the science and practice of coral reef restoration worldwide. Restor. Ecol. 29, 1–7 (2021).Article 

Google Scholar 
Heyward, A. J., Rees, M. & Smith, L. D. Coral spawning slicks harnessed for large-scale coral culture. Progr. Abstr. Int. Conf. Sci. Asp. Coral Reef Assess. Monit. Restor. 104, 188–189 (1999).
Google Scholar 
Harrison, P., Villanueva, R. & De la Cruz, D. Coral Reef Restoration using Mass Coral Larval Reseeding (Southern Cross University, 2016).
Google Scholar 
de la Cruz, D. W. & Harrison, P. L. Enhancing coral recruitment through assisted mass settlement of cultured coral larvae. PLoS ONE 15, e0242847. https://doi.org/10.1371/journal.pone.0242847 (2020).Article 

Google Scholar 
Chamberland, V. F., Snowden, S., Marhaver, K. L., Petersen, D. & Vermeij, M. J. A. The reproductive biology and early life ecology of a common Caribbean brain coral, Diploria labyrinthiformis (Scleractinia: Faviinae). Coral Reefs 36, 83–94 (2017).ADS 
Article 

Google Scholar 
Randall, C. J. et al. Sexual production of corals for reef restoration in the Anthropocene. Mar. Ecol. Prog. Ser. 635, 203–232 (2020).ADS 
Article 

Google Scholar 
Miller, M. W. et al. Settlement yields in large-scale in situ culture of Caribbean coral larvae for restoration. Restor. Ecol. https://doi.org/10.1111/rec.13512 (2021).Article 

Google Scholar 
Baria-Rodriguez, M. V., de la Cruz, D. W., Dizon, R. M., Yap, H. T. & Villanueva, R. D. Performance and cost-effectiveness of sexually produced Acropora granulosa juveniles compared with asexually generated coral fragments in restoring degraded reef areas. Aquat. Conserv. Mar. Freshw. Ecosyst. 29, 891–900 (2019).Article 

Google Scholar 
Doropoulos, C., Elzinga, J., ter Hofstede, R., van Koningsveld, M. & Babcock, R. C. Optimizing industrial-scale coral reef restoration: Comparing harvesting wild coral spawn slicks and transplanting gravid adult colonies. Restor. Ecol. 27, 758–767 (2019).Article 

Google Scholar 
Kuffner, I. B., Andersson, A. J., Jokiel, P. L., Rodgers, K. S. & MacKenzie, F. T. Decreased abundance of crustose coralline algae due to ocean acidification. Nat. Geosci. 1, 114–117 (2008).ADS 
Article 

Google Scholar 
Webster, N. S., Uthicke, S., Botté, E. S., Flores, F. & Negri, A. P. Ocean acidification reduces induction of coral settlement by crustose coralline algae. Glob. Change Biol. 19, 303–315 (2013).ADS 
Article 

Google Scholar 
Randall, C. J., Giuliano, C., Heyward, A. J. & Negri, A. P. Enhancing coral survival on deployment devices with microrefugia. Front. Mar. Sci. 8, 662263. https://doi.org/10.3389/fmars.2021.662263 (2021).Article 

Google Scholar 
Kuffner, I. B. et al. Inhibition of coral recruitment by macroalgae and cyanobacteria. Mar. Ecol. Prog. Ser. 323, 107–117 (2006).ADS 
Article 

Google Scholar 
Arnold, S. N., Steneck, R. S. & Mumby, P. J. Running the gauntlet: Inhibitory effects of algal turfs on the processes of coral recruitment. Mar. Ecol. Prog. Ser. 414, 91–105 (2010).ADS 
Article 

Google Scholar 
Speare, K. E., Duran, A., Miller, M. W. & Burkepile, D. E. Sediment associated with algal turfs inhibits the settlement of two endangered coral species. Mar. Pollut. Bull. 144, 189–195 (2019).PubMed 
Article 

Google Scholar 
Tebben, J., Guest, J. R., Sin, T. M., Steinberg, P. D. & Harder, T. Corals like it waxed: Paraffin-based antifouling technology enhances coral spat survival. PLoS ONE 9, 1–8 (2014).Article 

Google Scholar 
Almeida, E., Diamantino, T. C. & de Sousa, O. Marine paints: The particular case of antifouling paints. Prog. Org. Coat. 59, 2–20 (2007).Article 

Google Scholar 
Negri, A. P., Smith, L. D., Webster, N. S. & Heyward, A. J. Understanding ship-grounding impacts on a coral reef: Potential effects of anti-foulant paint contamination on coral recruitment. Mar. Pollut. Bull. 44, 111–117 (2002).PubMed 
Article 

Google Scholar 
Smith, L. D., Negri, A. P., Philipp, E., Webster, N. S. & Heyward, A. J. The effects of antifoulant-paint-contaminated sediments on coral recruits and branchlets. Mar. Biol. 143, 651–657 (2003).Article 

Google Scholar 
Jacobson, A. H. & Willingham, G. L. Sea-nine antifoulant: An environmentally acceptable alternative to organotin antifoulants. Sci. Total Environ. 258, 103–110 (2000).ADS 
PubMed 
Article 

Google Scholar 
Silva, V. et al. Isothiazolinone biocides: Chemistry, biological, and toxicity profiles. Molecules 25, 991. https://doi.org/10.3390/molecules25040991 (2020).Article 
PubMed Central 

Google Scholar 
da Silva, A. R., da Guerreiro, A. S., Martins, S. E. & Sandrini, J. Z. DCOIT unbalances the antioxidant defense system in juvenile and adults of the marine bivalve Amarilladesma mactroides (Mollusca: Bivalvia). Comp. Biochem. Physiol. Part C 250, 109169 (2021).
Google Scholar 
Cima, F. et al. Preliminary evaluation of the toxic effects of the antifouling biocide Sea-Nine 211TM in the soft coral Sarcophyton cf. glaucum (Octocorallia, Alcyonacea) based on PAM fluorometry and biomarkers. Mar. Environ. Res. 83, 16–22 (2013).PubMed 
Article 

Google Scholar 
Wendt, I., Backhaus, T., Blanck, H. & Arrhenius, Å. The toxicity of the three antifouling biocides DCOIT, TPBP and medetomidine to the marine pelagic copepod Acartia tonsa. Ecotoxicology 25, 871–879 (2016).PubMed 
Article 

Google Scholar 
Chen, L. et al. Identification of molecular targets for 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) in teleosts: New insight into mechanism of toxicity. Environ. Sci. Technol. 51, 1840–1847 (2017).ADS 
PubMed 
Article 

Google Scholar 
Martins, S. E., Fillmann, G., Lillicrap, A. & Thomas, K. V. Review: Ecotoxicity of organic and organo-metallic antifouling co-biocides and implications for environmental hazard and risk assessments in aquatic ecosystems. Biofouling 34, 34–52 (2018).PubMed 
Article 

Google Scholar 
Moon, Y. S., Kim, M., Hong, C. P., Kang, J. H. & Jung, J. H. Overlapping and unique toxic effects of three alternative antifouling biocides (Diuron, Irgarol 1051 ®, Sea-Nine 211 ® ) on non-target marine fish. Ecotoxicol. Environ. Saf. 180, 23–32 (2019).PubMed 
Article 

Google Scholar 
Su, Y. et al. Toxicity of 4,5-dichloro-2-n-octyl-4-isothiazolin-3-one (DCOIT) in the marine decapod Litopenaeus vannamei. Environ. Pollut. 251, 708–716 (2019).PubMed 
Article 

Google Scholar 
Fonseca, V. B., da Guerreiro, A. S., Vargas, M. A. & Sandrini, J. Z. Effects of DCOIT (4,5-dichloro-2-octyl-4-isothiazolin-3-one) to the haemocytes of mussels Perna perna. Comp. Biochem. Physiol Part C 232, 108737. https://doi.org/10.1016/j.cbpc.2020.108737 (2020).Article 

Google Scholar 
Ferreira, V. et al. Effects of nanostructure antifouling biocides towards a coral species in the context of global changes. Sci. Total Environ. 799, 149324 (2021).ADS 
PubMed 
Article 

Google Scholar 
de Campos, B. G. et al. A preliminary study on multi-level biomarkers response of the tropical oyster Crassostrea brasiliana to exposure to the antifouling biocide DCOIT. Mar. Pollut. Bull. 174, 112141 (2022).Article 

Google Scholar 
Maia, F. et al. Incorporation of biocides in nanocapsules for protective coatings used in maritime applications. Chem. Eng. J. 270, 150–157 (2015).Article 

Google Scholar 
Santos, J. V. N. et al. Can encapsulation of the biocide DCOIT affect the anti-fouling efficacy and toxicity on tropical bivalves?. Appl. Sci. 10, 1–12 (2020).Article 

Google Scholar 
Detty, M. R., Ciriminna, R., Bright, F. V. & Pagliaro, M. Environmentally benign sol-gel antifouling and foul-releasing coatings. Acc. Chem. Res. 47, 678–687 (2014).PubMed 
Article 

Google Scholar 
Korschelt, K., Tahir, M. N. & Tremel, W. A Step into the future: Applications of nanoparticle enzyme mimics. Chemistry 24, 9703–9713 (2018).PubMed 
Article 

Google Scholar 
Herget, K. et al. Haloperoxidase mimicry by CeO2-x nanorods combats biofouling. Adv. Mater. 29, 1–8 (2017).Article 

Google Scholar 
Korschelt, K. et al. CeO2-: X nanorods with intrinsic urease-like activity. Nanoscale 10, 13074–13082 (2018).PubMed 
Article 

Google Scholar 
Herget, K., Frerichs, H., Pfitzner, F., Tahir, M. N. & Tremel, W. Functional enzyme mimics for oxidative halogenation reactions that combat biofilm formation. Adv. Mater. 30, 1–28 (2018).Article 

Google Scholar 
Doropoulos, C., Ward, S., Marshell, A., Diaz-Pulido, G. & Mumby, P. J. Interactions among chronic and acute impacts on coral recruits: The importance of size-escape thresholds. Ecology 93, 2131–2138 (2012).PubMed 
Article 

Google Scholar 
Ji, Z. et al. Designed synthesis of CeO2 nanorods and nanowires for studying toxicological effects of high aspect ratio nanomaterials. ACS Nano 6, 5366–5380 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Herget, K. et al. Supporting Information: Haloperoxidase mimicry by CeO2-x nanorods combats biofouling. Adv. Mater. 29, 1603823 (2017).Article 

Google Scholar 
Sokolova, A. et al. Spontaneous multiscale phase separation within fluorinated xerogel coatings for fouling-release surfaces. Biofouling 28, 143–157 (2012).PubMed 
Article 

Google Scholar 
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
ImageJ Release Notes. https://imagej.nih.gov/ij/notes.html.Arganda-Carreras, I. et al. Trainable Weka Segmentation: A machine learning tool for microscopy pixel classification. Bioinformatics 33, 2424–2426 (2017).PubMed 
Article 

Google Scholar 
Arganda-Carreras, I. et al. Supplementary Data: Trainable Weka Segmentation: A Machine Learning Tool for Microscopy Pixel Classification: Trainable Weka Segmentation User Manualhttps://doi.org/10.1093/bioinformatics/btx180 (2017).Vyas, N., Sammons, R. L., Addison, O., Dehghani, H. & Walmsley, A. D. A quantitative method to measure biofilm removal efficiency from complex biomaterial surfaces using SEM and image analysis. Sci. Rep. 6, 2–11 (2016).Article 

Google Scholar 
Carbone, D. A., Gargano, I., Pinto, G., De Natale, A. & Pollio, A. Evaluating microalgae attachment to surfaces: A first approach towards a laboratory integrated assessment. Chem. Eng. Trans. 57, 73–78 (2017).
Google Scholar 
Moreno Osorio, J. H. et al. Early colonization stages of fabric carriers by two Chlorella strains. J. Appl. Phycol. 32, 3631–3644 (2020).Article 

Google Scholar 
Ricardo, G. F. et al. Impacts of water quality on Acropora coral settlement: The relative importance of substrate quality and light. Sci. Total Environ. 777, 146079. https://doi.org/10.1016/j.scitotenv.2021.146079 (2021).ADS 
Article 
PubMed 

Google Scholar 
Macadam, A., Nowell, C. J. & Quigley, K. Machine learning for the fast and accurate assessment of fitness in coral early life history. Remote Sens. 13, 1–17 (2021).Article 

Google Scholar 
Negri, A. P. & Heyward, A. J. Inhibition of Fertilization and Larval Metamorphosis of the Coral Acropora millepora (Ehrenberg, 1834) by Petroleum Products. Mar. Pollut. Bull. 41, 420–427 (2000).Article 

Google Scholar 
Nordborg, F. M., Flores, F., Brinkman, D. L., Agustí, S. & Negri, A. P. Phototoxic effects of two common marine fuels on the settlement success of the coral Acropora tenuis. Sci. Rep. 8, 1–12 (2018).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. (2021).Wickham, H. et al. Welcome to the Tidyverse. J. Open Source Softw. 4, 1686. https://doi.org/10.21105/joss.01686 (2019).ADS 
Article 

Google Scholar 
Pinheiro, J., Bates, D., Debroy, S., Sarkar, D. & R Core Team. Linear and nonlinear mixed effects models contact. Linear nonlinear Mix. Eff. Model. 3, 103–135 (2021).
Google Scholar 
Fox, J. & Weisberg, S. An R Companion to Applied Regression (Sage Publications, 2019).
Google Scholar 
Lenth, R. V. Emmeans: Estimated Marginal Means. https://cran.r-project.org/package=emmeans (2021).Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 9, 378–400 (2017).Article 

Google Scholar 
Dafforn, K. A., Lewis, J. A. & Johnston, E. L. Antifouling strategies: History and regulation, ecological impacts and mitigation. Mar. Pollut. Bull. 62, 453–465 (2011).PubMed 
Article 

Google Scholar 
Wu, R. et al. Room temperature synthesis of defective cerium oxide for efficient marine anti-biofouling. Adv. Compos. Hybrid Mater. https://doi.org/10.1007/s42114-021-00256-7 (2021).Article 

Google Scholar 
Hu, M. et al. Nanozymes in nanofibrous mats with haloperoxidase-like activity to combat biofouling. ACS Appl. Mater. Interfaces 10, 44722–44730 (2018).PubMed 
Article 

Google Scholar 
He, X. et al. Haloperoxidase mimicry by CeO2-x nanorods of different aspect ratios for antibacterial performance. ACS Sustain. Chem. Eng. 8, 6744–6752 (2020).Article 

Google Scholar 
Saxena, P. & Harish,. Nanoecotoxicological reports of engineered metal oxide nanoparticles on algae. Curr. Pollut. Rep. 4, 128–142 (2018).Article 

Google Scholar 
Xu, Y. et al. Effects of cerium oxide nanoparticles on bacterial growth and behaviors: Induction of biofilm formation and stress response. Environ. Sci. Pollut. Res. 26, 9293–9304 (2019).Article 

Google Scholar 
Xu, Y. et al. Mechanistic understanding of cerium oxide nanoparticle-mediated biofilm formation in Pseudomonas aeruginosa. Environ. Sci. Pollut. Res. 25, 34765–34776 (2018).Article 

Google Scholar 
Tang, Y. et al. Hybrid xerogel films as novel coatings for antifouling and fouling release. Biofouling 21, 59–71 (2005).PubMed 
Article 

Google Scholar 
Gunari, N. et al. The control of marine biofouling on xerogel surfaces with nanometer-scale topography. Biofouling 27, 137–149 (2011).PubMed 
Article 

Google Scholar 
Maia, F. et al. Silica nanocontainers for active corrosion protection. Nanoscale 4, 1287–1298 (2012).ADS 
PubMed 
Article 

Google Scholar 
Martins, R. et al. Effects of a novel anticorrosion engineered nanomaterial on the bivalve: Ruditapes philippinarum. Environ. Sci. Nano 4, 1064–1076 (2017).Article 

Google Scholar 
Gutner-Hoch, E. et al. Antimacrofouling efficacy of innovative inorganic nanomaterials loaded with booster biocides. J. Mar. Sci. Eng. 6, 15. https://doi.org/10.3390/jmse6010006 (2018).Article 

Google Scholar 
Negri, A. P. & Heyward, A. J. Inhibition of coral fertilisation and larval metamorphosis by tributyltin and copper. Mar. Environ. Res. 51, 17–27 (2001).PubMed 
Article 

Google Scholar 
Morse, D. E., Hooker, N., Morse, A. N. C. & Jensen, R. A. Control of larval metamorphosis and recruitment in sympatric agariciid corals. J. Exp. Mar. Biol. Ecol. 116, 193–217 (1988).Article 

Google Scholar 
Harrington, L., Fabricius, K., De’ath, G. & Negri, A. Recognition and selection of settlement substrata determine post-settlement survival in corals. Ecology 85, 3428–3437 (2004).Article 

Google Scholar 
Jorissen, H., Baumgartner, C., Steneck, R. S. & Nugues, M. M. Contrasting effects of crustose coralline algae from exposed and subcryptic habitats on coral recruits. Coral Reefs 39, 1767–1778 (2020).Article 

Google Scholar 
Figueiredo, J. et al. Toxicity of innovative anti-fouling nano-based solutions to marine species. Environ. Sci. Nano 6, 1418–1429 (2019).Article 

Google Scholar 
Shafir, S., Abady, S. & Rinkevich, B. Improved sustainable maintenance for mid-water coral nursery by the application of an anti-fouling agent. J. Exp. Mar. Biol. Ecol. 368, 124–128 (2009).Article 

Google Scholar  More

Bi, X. et al. Response of grassland productivity to climate change and anthropogenic activities in arid regions of Central Asia. Peer J. 8, e9797–e9797 (2020).PubMed 
PubMed Central 

Google Scholar 
Zhou, W. et al. Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 83, 303–313 (2017).
Google Scholar 
Liu, Y. Y. et al. Assessing the effects of climate variation and human activities on grassland degradation and restoration across the globe. Ecol. Indic. 106, 105504–105504 (2019).
Google Scholar 
Zhang, Y. et al. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012. Sci. Total Environ. 563–564, 210–220 (2016).ADS 
PubMed 

Google Scholar 
Wang, Z. et al. Quantitative assess the driving forces on the grassland degradation in the Qinghai-Tibet Plateau, China. Ecol. Inf. 33, 32–44 (2016).CAS 

Google Scholar 
He, C. Y., Tian, J., Gao, B. & Zhao, Y. Y. Differentiating climate- and human-induced drivers of grassland degradation in the Liao River Basin, China. Environ. Monit. Assess. 187(1), 4199 (2015).PubMed 

Google Scholar 
Liu, Y. Y. et al. Grassland dynamics in responses to climate variation and human activities in China from 2000 to 2013. Sci. Total Environ. 690, 27–39 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Jiang, L. L., Jiapaer, G., Bao, A. M., Guo, H. & Ndayisaba, F. Vegetation dynamics and responses to climate change and human activities in Central Asia. Sci. Total Environ. 599–600, 967–980 (2017).ADS 
PubMed 

Google Scholar 
Chen, T. et al. Disentangling the relative impacts of climate change and human activities on arid and semiarid grasslands in Central Asia during 1982–2015. Sci. Total Environ. 653, 1311–1325 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Gang, C. et al. The impacts of land conversion and management measures on the grassland net primary productivity over the Loess Plateau, Northern China. Sci. Total Environ. 645, 827–836 (2018).ADS 
CAS 
PubMed 

Google Scholar 
Guo, D. & Wang, H. Simulation of permafrost and seasonally frozen ground conditions on the Tibetan Plateau, 1981–2010. J. Geophys. Res. Atmos. 118, 5216–5230 (2013).ADS 

Google Scholar 
Yang, Y. et al. Comparative assessment of grassland degradation dynamics in response to climate variation and human activities in China, Mongolia, Pakistan and Uzbekistan from 2000 to 2013. J. Arid Environ. 135, 164–172 (2016).ADS 

Google Scholar 
Li, C. X., Jong, R., Schmid, B., Wulf, H. & Michael, E. S. Changes in grassland cover and in its spatial heterogeneity indicate degradation on the Qinghai-Tibetan Plateau. Ecol. Indic. 119, 106641 (2020).
Google Scholar 
Li, F., Chen, W., Zeng, Y., Zhao, Q. J. & Wu, B. F. Improving estimates of grassland fractional vegetation cover based on a pixel dichotomy model: A case study in Inner Mongolia, China. Remote Sens. 6, 4705–4722 (2014).ADS 

Google Scholar 
Wang, J., Brown, D. G. & Chen, J. Q. Drivers of the dynamics in net primary productivity across ecological zones on the Mongolian plateau. Landsc. Ecol. 28(4), 725–739 (2014).
Google Scholar 
Han, D. M. et al. Evaluation of semiarid grassland degradation in north China from multiple perspectives. Ecol. Eng. 112, 41–50 (2018).
Google Scholar 
Liu, H. X. et al. Response of vegetation productivity to climate change and human activities in the Shaanxi–Gansu–Ningxia region, China. J. Indian Soc. Remote Sens. 46(7), 1081–1092 (2018).
Google Scholar 
Zheng, K. et al. Impacts of climate change and human activities on grassland vegetation variation in the Chinese Loess Plateau. Sci. Total Environ. 660, 236–244 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Yan, Y. C., Liu, X. P., Wen, Y. Y. & Ou, J. P. Quantitative analysis of the contributions of climatic and human factors to grassland productivity in northern China. Ecol. Indic. 103, 542–553 (2019).
Google Scholar 
Wang, H. et al. Impacts of climate change on net primary productivity in arid and semiarid regions of China. Chin. Geogra. Sci. 26, 35–47 (2016).CAS 

Google Scholar 
Thomas, M. et al. Human land-use practices lead to global long-term increases in photosynthetic capacity. Remote Sens. 6(6), 5717–5731 (2014).
Google Scholar 
Becerril-Pina, R., Mastachi-Loza, C. A., Gonzalez-Sosa, E., Diaz-Delgado, C. & Ba, K. M. Assessing desertification risk in the semi-arid highlands of central Mexico. J. Arid Environ. 120, 4–13 (2015).ADS 

Google Scholar 
Evans, J. & Geerken, R. Discrimination between climate and human-induced dryland degradation. J. Arid Environ. 57(4), 535–554 (2004).ADS 

Google Scholar 
Meng, M. et al. Vegetation change in response to climate factors and human activities on the Mongolian Plateau. Peer J. 7, e7735 (2019).PubMed 
PubMed Central 

Google Scholar 
Burrell, A. L., Evans, J. P. & Liu, Y. Detecting dryland degradation using time series segmentation and residual trend analysis (TSS-RESTREND). Remote Sens Environ. 197, 43–57 (2017).ADS 

Google Scholar 
Gedefaw, M. G., Geli, H. M. E. & Abera, T. A. Assessment of rangeland degradation in New Mexico using time series segmentation and residual trend analysis (TSS-RESTREND). Remote Sens. 13(9), 1618–1618 (2021).ADS 

Google Scholar 
Zhang, F. Changes of Grassland Net Primary Productivity in the Qinghai Tibet Plateau During the Past 34 Years and Analysis of Its Local Degradation Characteristics (Lanzhou University, 2021).
Google Scholar 
Li, L. H. et al. Current challenges in distinguishing climatic and anthropogenic contributions to alpine grassland variation on the Tibetan Plateau. Ecol. Evol. 8(11), 5949–5963 (2018).PubMed 
PubMed Central 

Google Scholar 
Zhu, Z. C. et al. Greening of the earth and its drivers. Nat. Clim. Change. 6, 791–795 (2016).ADS 
CAS 

Google Scholar 
Song, L. C., Ma, W. W., Li, G., Liu, S. N. & Lu, G. Effect of temperature changes on nitrogen mineralization in soils with different degradation gradients in Gahai Wetland. Acta Pratacul. Sin. 30(09), 27–37 (2021).
Google Scholar 
Dai, L. C. et al. Effect of grazing management strategies on alpine grassland on the northeastern Qinghai-Tibet Plateau. Ecol. Eng. 173, 106418 (2021).
Google Scholar 
Liu, Y. Y. et al. Evaluating the dynamics of grassland net primary productivity in response to climate change in China. Glob. Ecol. Conserv. 28, e01574 (2021).
Google Scholar 
Bestelmeyer, B. T., Duniway, M. C., James, D. K., Burkett, L. M. & Havstad, K. M. A test of critical thresholds and their indicators in a desertification-prone ecosystem: More resilience than we thought. Ecol. Lett. 16, 339–345 (2013).PubMed 

Google Scholar 
Kéfi, S. et al. Early warning signals of ecological transitions: Methods for spatial patterns. PLoS ONE 9(3), e92097 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Li, J. Z. et al. IKONOS image-based extraction of the distribution area of Stellera chamaejasme L. in Qilian County of Qinghai Province, China. Remote Sens. 8(2), 148 (2016).ADS 

Google Scholar 
Liu, Y. Q. & Lu, C. H. Quantifying grass coverage trends to identify the hot plots of grassland degradation in the Tibetan Plateau during 2000–2019. Int. J. Environ. Res. Public Health. 18(2), 416 (2021).MathSciNet 
PubMed Central 

Google Scholar 
Kendall, M. G. Rank Correlation Methods (Griffin, 1948).MATH 

Google Scholar 
Mann, H. B. Nonparametric tests against trend. Econometrica 13, 245–259 (1945).MathSciNet 
MATH 

Google Scholar 
Zhang, Z. M. & Lu, C. H. Clustering analysis of soybean production to understand its spatiotemporal dynamics in the North China Plain. Sustainability. 12(15), 6178 (2020).
Google Scholar 
Pei, T. T. et al. The sensitivity of vegetation phenology to extreme climate indices in the Loess Plateau, China. Sustainability. 13(14), 7623–7623 (2021).
Google Scholar 
Lu, B. B., Charlton, M., Harris, P. & Fotheringham, A. S. Geographically weighted regression with a non-Euclidean distance metric: A case study using hedonic house price data. Int. J. Geogr. Inf. Sci. 28(4), 660–681 (2014).
Google Scholar 
Sun, L. Q., Zhang, F. H., Yang, S. W., Qiu, A. G. & Zhang, X. L. The method of selecting geographically and temporally weight regression variable based on stepwise regression. Sci. Surv. Mapp. 44(01), 73–78+97 (2019).
Google Scholar 
Jiang, W. G. et al. Spatio-temporal analysis of vegetation variation in the Yellow River basin. Ecol. Indic. 51, 117–126 (2015).
Google Scholar 
Ndayisaba, F. et al. Understanding the spatial temporal vegetation dynamics in Rwanda. Remote Sens. 8(2), 129 (2016).ADS 

Google Scholar 
Kéfi, S. et al. Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems. Nature 449(7159), 213–217 (2007).ADS 
PubMed 

Google Scholar 
Chen, J. J., Yi, S. H. & Qin, Y. The contribution of plateau pika disturbance and erosion on patchy alpine grassland soil on the Qinghai-Tibetan Plateau: Implications for grassland restoration. Geoderma 297, 1–9 (2017).ADS 
CAS 

Google Scholar 
Cai, H. Y., Yang, X. H. & Xu, X. L. Human-induced grassland degradation/restoration in the central Tibetan Plateau: The effects of ecological protection and restoration projects. Ecol. Eng. 83, 112–119 (2015).
Google Scholar 
Wang, P., Lassoie, J. P., Morreale, S. J. & Dong, S. K. A critical review of socioeconomic and natural factors in ecological degradation on the Qinghai-Tibetan Plateau. China. Rangel. J. 37(1), 1–9 (2015).
Google Scholar 
Lu, C. B. & Hou, L. F. Cause analysis and Control Countermeasures of grassland degradation in Qilian County, Qinghai Province. Today Anim. Husb. Vet. Med. 34(02), 62 (2018).
Google Scholar 
Guo, X. W. et al. Light grazing significantly reduces soil water storage in Alpine Grasslands on the Qinghai-Tibet Plateau. Sustainability. 12(6), 2523–2523 (2020).
Google Scholar 
Bai, Y. F. et al. Climate warming benefits alpine vegetation growth in Three-River Headwater Region, China. Sci. Total Environ. 742, 140574–140574 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Chen, T. et al. Unraveling the relative impacts of climate change and human activities on grassland productivity in Central Asia over last three decades. Sci. Total Environ. 743, 140649 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Li, A., Wu, J. G. & Huang, J. H. Distinguishing between human-induced and climate-driven vegetation changes: A critical application of RESTREND in inner Mongolia. Landsc. Ecol. 27(7), 969–982 (2012).CAS 

Google Scholar 
Wu, J. S. et al. Disentangling climatic and anthropogenic contributions to nonlinear dynamics of alpine grassland productivity on the Qinghai-Tibetan Plateau. J. Environ. Manag. 281, 111875–111875 (2020).
Google Scholar 
Gang, C. et al. Comparative assessment of grassland NPP dynamics in response to climate change in China, North America, Europe and Australia from 1981 to 2010. J. Agron. Crop Sci. 201(1), 57–68 (2015).
Google Scholar 
Gang, C. C. et al. Quantitative assessment of the contributions of climate change and human activities on global grassland degradation. Environ. Earth Sci. 72(11), 4273–4282 (2014).
Google Scholar 
Chen, Y. Z. et al. Grassland carbon sequestration ability in China: A new perspective from terrestrial aridity zones. Rangeland Ecol. Manag. 69(1), 84–94 (2016).
Google Scholar 
Mowll, W. et al. Climatic controls of aboveground net primary production in semi-arid grasslands along a latitudinal gradient portend low sensitivity to warming. Oecologia 177(4), 959–969 (2015).ADS 
PubMed 

Google Scholar 
Zhou, Y. et al. Climate contributions to vegetation variations in central Asian Drylands: Pre- and post-USSR collapse. Remote Sens. 7(3), 2449–2470 (2015).ADS 

Google Scholar 
Ji, Y. et al. Variation of net primary productivity and its drivers in China’s forests during 2000–2018. For. Ecosyst. 7(1), 1–11 (2020).CAS 

Google Scholar 
Zeng, B. & Yang, T. B. Impacts of climate warming on vegetation in Qaidam Area from 1990 to 2003. Environ. Monit. Assess. 144(1–3), 403–417 (2008).PubMed 

Google Scholar 
Duan, A. M. & Xiao, Z. X. Does the climate warming hiatus exist over the Tibetan Plateau?. Sci. Rep. 5(1), 13711 (2015).ADS 
PubMed 
PubMed Central 

Google Scholar 
Fang, J. Y. et al. Precipitation patterns alter growth of temperate vegetation. Geophys. Res. Lett. 32(21), L21411 (2005).ADS 

Google Scholar 
Zhao, X., Tan, K., Zhao, S. & Fang, J. Changing climate affects vegetation growth in the arid region of the northwestern China. J. Arid Environ. 75(10), 946–952 (2011).ADS 

Google Scholar 
Ukkola, A. M. et al. Reduced streamflow in water-stressed climates consistent with CO2 effects on vegetation. Nat. Clim. Change. 6(1), 75–78 (2016).ADS 

Google Scholar 
Dong, S. K., Shang, Z. H., Gao, J. X. & Boone, R. B. Enhancing sustainability of grassland ecosystems through ecological restoration and grazing management in an era of climate change on Qinghai-Tibetan Plateau. Agric. Ecosyst. Environ. 287(C), 106684 (2019).
Google Scholar 
Xu, H. P. et al. Responses of plant productivity and soil nutrient concentrations to different alpine grassland degradation levels. Environ Monit Assess. 191(11), 678 (2019).CAS 
PubMed 

Google Scholar 
Wen, W. Y. et al. Research on soil net nitrogen mineralization in Stipa grandis grassland with different stages of degradation. Geosci J. 20(4), 485–494 (2016).ADS 
CAS 

Google Scholar 
She, Y. et al. Vegetation attributes and soil properties of alpine grassland in different degradation stages on the Qinghai-Tibet Plateau, China: A meta-analysis. Arab J Geosci. 15, 193 (2022).
Google Scholar 
Xu, G. P. Study on the Change of Vegetation and Soil Nutrients of Alpine Meadow Under Different Degradation Degrees in Eastern Qilian Mountains (Gansu Agricultural University, 2006).
Google Scholar 
Anderson, K. et al. Vegetation expansion in the subnival Hindu Kush Himalaya. Glob. Chang. Biol. 26(3), 1608–1625 (2020).ADS 
PubMed 
PubMed Central 

Google Scholar 
Chen, B. X. et al. The impact of climate change and anthropogenic activities on alpine grassland over the Qinghai-Tibet Plateau. Agric. For. Meteorol. 189–190, 11–18 (2014).ADS 

Google Scholar 
Zhang, X. W., Li, G., Dong, K. H. & Zhao, X. Effects of grazing and enclosure on community characteristics and biodiversity in Leymus chinensis grassland. J. Grassl. Forage Sci. 4, 22–27 (2019).
Google Scholar 
Huang, K. et al. The influences of climate change and human activities on vegetation dynamics in the Qinghai-Tibet Plateau. Remote Sens. 8(10), 876 (2016).ADS 

Google Scholar 
Duan, Q. T., Luo, L. H., Zhao, W. Z., Zhuang, Y. L. & Liu, F. Mapping and evaluating human pressure changes in the Qilian mountains. Remote Sens. 13(12), 2400–2400 (2021).ADS 

Google Scholar 
Wang, Y. et al. Performance and obstacle tracking to natural forest resource protection project: A rangers’ case of Qilian mountain, China. Int. J. Environ. Res. Public Health. 17(16), 5672 (2020).PubMed Central 

Google Scholar 
Li, Z. Y. et al. Changes in nutrient balance, environmental effects, and green development after returning farmland to forests: A case study in Ningxia, China. Sci. Total Environ. 735, 139370 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Li, C. X., de Jong, R., Schmid, B., Wulf, H. & Schaepman, M. E. Spatial variation of human influences on grassland biomass on the Qinghai-Tibetan plateau. Sci. Total Environ. 665, 678–689 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Li, X. L. et al. Rangeland degradation on the Qinghai-Tibet Plateau: Implications for rehabilitation. Land Degrad. Dev. 24, 72–80 (2011).
Google Scholar 
Li, C. B. et al. Regional vegetation dynamics and its response to climate change—a case study in the Tao River Basin in Northwestern China. Environ. Res. Lett. 9(12), 125003–125003 (2014).ADS 

Google Scholar 
Liu, Y. Y. et al. Untangling the effects of management measures, climate and land use cover change on grassland dynamics in the Qinghai-Tibet Plateau, China. Land Degrad. Dev. 32(17), 4974–4987 (2021).
Google Scholar 
Hou, X. Chinese Grassland Science (Science Press, 2013) (In Chinese).
Google Scholar  More

Royse DJ. A global perspective on the high five: Agaricus, Pleurotus, Lentinula, Auricularia and Flammulina. In: Singh M, editor. Proceedings of the 8th International Conference on Mushroom Biology and Mushroom Products. New Delhi; 2014. p. 1–6.Vos AM, Heijboer A, Boschker HTS, Bonnet B, Lugones LG, Wosten HAB. Microbial biomass in compost during colonization of Agaricus bisporus. AMB Express. 2017; 7:12.Jurak E, Punt AM, Arts W, Kabel MA, Gruppen H. Fate of carbohydrates and lignin during composting and mycelium growth of Agaricus bisporus on wheat straw based compost. PLoS ONE. 2015;10:e0138909.PubMed 
PubMed Central 
Article 

Google Scholar 
Beyer DM. Basic procedures for Agaricus mushroom growing PennState Extension: the Pennsylvania State University. 2003. https://extension.psu.edu/basic-procedures-for-agaricus-mushroom-growing.Wang L, Mao J, Zhao H, Li M, Wei Q, Zhou Y, et al. Comparison of characterization and microbial communities in rice straw- and wheat straw-based compost for Agaricus bisporus production. J Ind Microbiol Biotechnol. 2016;43:1249–60.CAS 
PubMed 
Article 

Google Scholar 
Adams JDW, Frostick LE. Investigating microbial activities in compost using mushroom (Agaricus bisporus) cultivation as an experimental system. Bioresour Technol. 2008;99:1097–102.CAS 
PubMed 
Article 

Google Scholar 
Liu L, Wang S, Guo X, Zhao T, Zhang B. Succession and diversity of microorganisms and their association with physicochemical properties during green waste thermophilic composting. Waste Manage. 2018;73:101–12.CAS 
Article 

Google Scholar 
Reyes-Torres M, Oviedo-Ocana ER, Dominguez I, Komilis D, Sanchez A. A systematic review on the composting of green waste: feedstock quality and optimization strategies. Waste Manage. 2018;77:486–99.CAS 
Article 

Google Scholar 
Pardo‐Giménez A, González JEP, Zied DC. Casing materials and techniques in Agaricus bisporus cultivation. In: Zied DC, Pardo‐Giménez A, editors. Edible and medicinal mushrooms technology and applications. Chichester, UK: Wiley; 2017. p. 149–74.Baars JJP, Scholtmeijer K, Sonnenberg ASM, van Peer A. Critical factors involved in primordia building in Agaricus bisporus: a review. Molecules. 2020;25:2984.Vieira FR, Pecchia JA. Bacterial community patterns in the Agaricus bisporus cultivation system, from compost raw materials to mushroom caps. Microb Ecol. 2021;84:20–32.PubMed 
Article 

Google Scholar 
Kristensen JB, Thygesen LG, Felby C, Jorgensen H, Elder T. Cell-wall structural changes in wheat straw pretreated for bioethanol production. Biotechnol Biofuels. 2008;1:1–9.Article 

Google Scholar 
Jurak E, Patyshakuliyeva A, de Vries RP, Gruppen H, Kabel MA. Compost grown Agaricus bisporus lacks the ability to degrade and consume highly substituted xylan fragments. PLoS ONE. 2015;10:e0134169.PubMed 
PubMed Central 
Article 

Google Scholar 
Ryckeboer J, Mergaert J, Vaes K, Klammer S, De Clercq D, Coosemans J, et al. A survey of bacteria and fungi occurring during composting and self-heating processes. Ann Microbiol. 2003;53:349–410.
Google Scholar 
Kutzner HJ. Microbiology of composting. In: Rehm H-J, Reed G, editors. Biotechnology. 11c. 2nd ed. Verlag: Wiley-VCH; 2000. p. 35–100.Carrasco J, Garcia-Delgado C, Lavega R, Tello ML, De Toro M, Barba-Vicente V, et al. Holistic assessment of the microbiome dynamics in the substrates used for commercial champignon (Agaricus bisporus) cultivation. Microb Biotechnol. 2020;13:1933–47.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vieira FR, Pecchia JA. Bacterial community patterns in the Agaricus bisporus cultivation system, from compost raw materials to mushroom caps. Microb Ecol. 2021;82. https://doi.org/10.1007/s00248-021-1833-5.Vieira FR, Pecchia JA. An exploration into the bacterial community under different pasteurization conditions during substrate preparation (composting–Phase II) for Agaricus bisporus cultivation. Microb Ecol. 2018;75:318–30.CAS 
PubMed 
Article 

Google Scholar 
Cao GT, Song TT, Shen YY, Jin QL, Feng WL, Fan LJ, et al. Diversity of bacterial and fungal communities in wheat straw compost for Agaricus bisporus cultivation. Hortscience. 2019;54:100–9.CAS 
Article 

Google Scholar 
Wiegant WM. Growth characteristics of the thermophilic fungus Scytalidium thermophilum in relation to production of mushroom compost. Appl Environ Microbiol. 1992;58:1301–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fermor T, Randle P, Smith J. Compost as a substrate and its preparation. In: Flegg PB, Spencer DM, Wood D, editors. The biology and technology of the cultivated mushroom. Chichester, UK: John Wiley & Sons, Ltd; 1985. p. 81–109.Straatsma G, Samson RA, Olijnsma TW, Op den Camp HJM, Gerrits JPG, Griensven LJLDV. Ecology of thermophilic fungi in mushroom compost, with emphasis on Scytalidium thermophilum and growth stimulation of Agaricus bisporus mycelium. Appl Environ Microbiol. 1994;60:454–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ross RC, Harris PJ. An investigation into the selective nature of mushroom compost. Sci Hortic. 1983;19:55–64.Article 

Google Scholar 
Coello-Castillo MM, Sanchez JE, Royse DJ. Production of Agaricus bisporus on substrates pre-colonized by Scytalidium thermophilum and supplemented at casing with protein-rich supplements. Bioresour Technol. 2009;100:4488–92.CAS 
PubMed 
Article 

Google Scholar 
Szekely A, Sipos R, Berta B, Vajna B, Hajdu C, Marialigeti K. DGGE and T-RFLP analysis of bacterial succession during mushroom compost production and sequence-aided T-RFLP profile of mature compost. Microb Ecol. 2009;57:522–33.PubMed 
Article 

Google Scholar 
Kertesz M, Safianowicz K, Bell TL. New insights into the microbial communities and biological activities that define mushroom compost. Sci Cultiv Edible Fungi. 2016;19:161–5.
Google Scholar 
McGee CF, Byrne H, Irvine A, Wilson J. Diversity and dynamics of the DNA and cDNA-derived bacterial compost communities throughout the Agaricus bisporus mushroom cropping process. Ann Microbiol. 2017;67:751–61.CAS 
Article 

Google Scholar 
McGee CF, Byrne H, Irvine A, Wilson J. Diversity and dynamics of the DNA- and cDNA-derived compost fungal communities throughout the commercial cultivation process for Agaricus bisporus. Mycologia. 2017;109:475–84.CAS 
PubMed 
Article 

Google Scholar 
Yeates C, Gillings MR. Rapid purification of DNA from soil for molecular biodiversity analysis. Lett Appl Microbiol. 1998;27:49–53.CAS 
Article 

Google Scholar 
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
White TJ, Bruns T, Lee S, Taylor JW. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, editors. PCR protocols: a guide to methods and applications. New York: Academic Press; 1990. p. 315–22.
Google Scholar 
Lever MA, Torti A, Eickenbusch P, Michaud AB, Santl-Temkiv T, Jorgensen BB. A modular method for the extraction of DNA and RNA, and the separation of DNA pools from diverse environmental sample types. Front Microbiol. 2015;6:476.Muyzer G, Waal ECD, Uitterlinden AG. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol. 1993;59:695–700.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108:4516–22.CAS 
PubMed 
Article 

Google Scholar 
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation For Statistical Computing; 2019.Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina aplicon data. Nat Meth. 2016;13:581–3.CAS 
Article 

Google Scholar 
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucl Acids Res. 2012;41:D590–6.PubMed 
PubMed Central 
Article 

Google Scholar 
Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics. 2010;27:592–3.PubMed 
PubMed Central 
Article 

Google Scholar 
Wright ES. Using DECIPHER v2.0 to analyze big biological sequence data in R. R J. 2016;8:352–9.Article 

Google Scholar 
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dixon P. VEGAN, a package of R functions for community ecology. J Veget Sci. 2003;14:927–30.Article 

Google Scholar 
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2016.Sharma HS, Kilpatrick M. Mushroom (Agaricus bisporus) compost quality factors for predicting potential yield of fruiting bodies. Can J Microbiol. 2000;46:515–9.CAS 
PubMed 
Article 

Google Scholar 
Seaby DA. Mushroom (Agaricus bisporus) yield modelling for the bag method of mushroom production using commercial yields and from micro plots. Sci Cultiv Edible Fungi. 1995;14:409–16.
Google Scholar 
O’Donoghue DC. Relationship between some compost factors and their effects on yield of Agaricus. Mushroom Sci. 1965;6:245–54.
Google Scholar 
Andersen B, Sorensen JL, Nielsen KF, van den Ende BG, de Hoog S. A polyphasic approach to the taxonomy of the Alternaria infectoria species-group. Fungal Genet Biol. 2009;46:642–56.CAS 
PubMed 
Article 

Google Scholar 
van den Brink J, Samson RA, Hagen F, Boekhout T, de Vries RP. Phylogeny of the industrial relevant, thermophilic genera Myceliophthora and Corynascus. Fungal Divers. 2012;52:197–207.Article 

Google Scholar 
Souza TP, Marques SC, Santos D, Dias ES. Analysis of thermophilic fungal populations during phase II of composting for the cultivation of Agaricus subrufescens. World J Microbiol Biotechnol. 2014;30:2419–25.PubMed 
Article 

Google Scholar 
Vajna B, Szili D, Nagy A, Márialigeti K. An improved sequence-aided T-RFLP analysis of bacterial succession during oyster mushroom substrate preparation. Microb Ecol. 2012;64:702–13.CAS 
PubMed 
Article 

Google Scholar 
Du R, Yan J, Li S, Zhang L, Zhang S, Li J, et al. Cellulosic ethanol production by natural bacterial consortia is enhanced by Pseudoxanthomonas taiwanensis. Biotechnol Biofuels. 2015;8:10.Kato S, Haruta S, Cui ZJ, Ishii M, Igarashi Y. Stable coexistence of five bacterial strains as a cellulose-degrading community. Appl Environ Microbiol. 2005;71:7099–106.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Haruta S, Cui Z, Huang Z, Li M, Ishii M, Igarashi Y. Construction of a stable microbial community with high cellulose-degradation ability. Appl Microbiol Biotechnol. 2002;59:529–34.CAS 
PubMed 
Article 

Google Scholar 
Vajna B, Adrienn N, Sajben-Nagy E, Manczinger L, Szijártó N, Kádár Z, et al. Microbial community structure changes during oyster mushroom substrate preparation. Appl Microbiol Biotechnol. 2010;86:367–75.CAS 
PubMed 
Article 

Google Scholar 
Karadag D, Özkaya B, Ölmez E, Nissilä ME, Çakmakçı M, Yıldız Ş, et al. Profiling of bacterial community in a full-scale aerobic composting plant. Int Biodeter Biodeg. 2013;77:85–90.CAS 
Article 

Google Scholar 
Rathinam NK, Gorky, Bibra M, Salem DR, Sani RK. Bioelectrochemical approach for enhancing lignocellulose degradation and biofilm formation in Geobacillus strain WSUCF1. Bioresour Technol. 2020;295:122271.Song TT, Shen YY, Jin QL, Feng WL, Fan LJ, Cao GT, et al. Bacterial community diversity, lignocellulose components, and histological changes in composting using agricultural straws for Agaricus bisporus production. PeerJ. 2021;9:e10452.PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang X, Zhong Y, Yang S, Zhang W, Xu M, Ma A, et al. Diversity and dynamics of the microbial community on decomposing wheat straw during mushroom compost production. Bioresour Technol. 2014;170:183–95.CAS 
PubMed 
Article 

Google Scholar 
Goodfellow M, Maldonado LA, Quintana ET. Reclassification of Nonomuraea flexuosa (Meyer 1989) Zhang et al. 1998 as Thermopolyspora flexuosa gen. nov., comb. nov., nom. rev. Int J Syst Evol Microbiol. 2005;55:1979–83.CAS 
PubMed 
Article 

Google Scholar 
Lin SB, Stutzenberger FJ. Purification and characterization of the major beta-1,4-endoglucanase from Thermomonospora curvata. J Appl Bacteriol. 1995;79:447–53.CAS 
PubMed 
Article 

Google Scholar 
Kukolya J, Nagy I, Láday M, Tóth E, Oravecz O, Márialigeti K, et al. Thermobifida cellulolytica sp. nov., a novel lignocellulose-decomposing actinomycete. Int J Syst Evol Microbiol. 2002;52:1193–9.CAS 
PubMed 

Google Scholar 
Weon H-Y, Lee S-Y, Kim B-Y, Noh H-J, Schumann P, Kim J-S, et al. Ureibacillus composti sp. nov. and Ureibacillus thermophilus sp. nov., isolated from livestock-manure composts. Int J Syst Evol Microbiol. 2007;57:2908–11.CAS 
PubMed 
Article 

Google Scholar 
Poli A, Laezza G, Gul-Guven R, Orlando P, Nicolaus B. Geobacillus galactosidasius sp. nov., a new thermophilic galactosidase-producing bacterium isolated from compost. Syst Appl Microbiol. 2011;34:419–23.CAS 
PubMed 
Article 

Google Scholar 
Gavande PV, Basak A, Sen S, Lepcha K, Murmu N, Rai V, et al. Functional characterization of thermotolerant microbial consortium for lignocellulolytic enzymes with central role of Firmicutes in rice straw depolymerization. Sci Rep. 2021;11:3032.Xu JQ, Lu YY, Shan GC, He XS, Huang JH, Li QL. Inoculation with compost-born thermophilic complex microbial consortium induced organic matters degradation while reduced nitrogen loss during co-composting of dairy manure and sugarcane leaves. Waste Biomass Valor. 2019;10:2467–77.CAS 
Article 

Google Scholar 
Yoon JH, Kang SJ, Im WT, Lee ST, Oh TK. Chelatococcus daeguensis sp nov., isolated from wastewater of a textile dye works, and emended description of the genus Chelatococcus. Int J Syst Evol Microbiol. 2008;58:2224–8.CAS 
PubMed 
Article 

Google Scholar 
Zhou C, Liu Z, Huang Z-L, Dong M, Yu X-L, Ning P. A new strategy for co-composting dairy manure with rice straw: addition of different inocula at three stages of composting. Waste Manage. 2015;40:38–43.CAS 
Article 

Google Scholar 
Gómez A. New technology in Agaricus bisporus cultivation. In: Zied DC, Pardo-Giménez A, editors. Edible and medicinal mushrooms. Chichester, UK: John Wiley & Sons; 2017. p. 211–20.von Minnigerode HF, editor. Method for controlling and regulating the composting process. Proceedings of the Eleventh International Scientific Congress on the Cultivation of Edible Fungi. Sydney, Australia: The International Society for Mushroom Science; 1981.Jurak E, Gruppen H, Kabel MA, Eggink G, Meyer AS, van der Maarel MJEC, et al. How mushrooms feed on compost: conversion of carbohydrates and lignin in industrial wheat straw based compost enabling the growth of Agaricus bisporus. Wageningen University—Graduate School VLAG; 2015.Miller FC, Macauley BJ, Harper ER. Investigation of various gases, pH and redox potential in mushroom composting Phase-I stacks. Aust J Exper Agric. 1991;31:415–25.Article 

Google Scholar 
Miller FC, Harper ER, Macauley BJ, Gulliver A. Composting based on moderately thermophilic and aerobic conditions for the production of commercial growing compost. Aust J Exper Agric. 1990;30:287–96.Article 

Google Scholar 
Carrasco J, Preston GM. Growing edible mushrooms: a conversation between bacteria and fungi. Environ Microbiol. 2020;22:858–72.PubMed 
Article 

Google Scholar  More

Jones, O. R. et al. Diversity of ageing across the tree of life. Nature 505, 169–173 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Nussey, D. H., Froy, H., Lemaitre, J.-F., Gaillard, J.-M. & Austad, S. N. Senescence in natural populations of animals: Widespread evidence and its implications for bio-gerontology. Ageing Res. Rev. 12, 214–225 (2013).PubMed 
Article 

Google Scholar 
Sergio, F. et al. Variation in age-structured vital rates of a long-lived raptor: Implications for population growth. Basic Appl. Ecol. 12, 107–115 (2011).Article 

Google Scholar 
Sydeman, W. J., Huber, H. R., Emslie, S. D., Ribic, C. A. & Nur, N. Age-specific weaning success of northern elephant seals in relation to previous breeding experience. Ecology 76, 2204–2217 (1991).Article 

Google Scholar 
Gaillard, J.-M., Festa-Bianchet, M., Yoccoz, N. G., Loison, A. & Toïgo, C. Temporal variation in fitness components and population dynamics of large herbivores. Annu. Rev. Ecol. Syst. 31, 367–393 (2000).Article 

Google Scholar 
Caswell, H. Matrix Population Models. Online Library (Wiley, 2001).
Google Scholar 
Colchero, F. et al. The diversity of population responses to environmental change. Ecol. Lett. 22, 342–353 (2019).PubMed 
Article 

Google Scholar 
Morris, W. F. & Doak, D. F. Quantitative Conservation Biology (Sinauer, 2002).
Google Scholar 
Stott, I., Townley, S. & Hodgson, D. J. A framework for studying transient dynamics of population projection matrix models: A synthesis of transient demography. Ecol. Lett. 14, 959–970 (2011).PubMed 
Article 

Google Scholar 
Vitousek, P. M. Beyond global warming: Ecology and global change. Ecology 75, 1861–1876 (1994).Article 

Google Scholar 
Gamelon, M. et al. Influence of life-history tactics on transient dynamics: A comparative analysis across mammalian populations. Am. Nat. 184, 673–683 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Jackson, J., Mar, K. U., Htut, W., Childs, D. Z. & Lummaa, V. Changes in age?structure over four decades were a key determinant of population growth rate in a long-lived mammal. J. Anim. Ecol. 89, 2268–2278 (2020).PubMed 
Article 

Google Scholar 
Lebreton, J.-D., Burnham, K. P., Clobert, J. & Anderson, D. R. Modeling survival and testing biological hypotheses using marked animals: A unified approach with case studies. Ecol. Monogr. 62, 67–118 (1992).Article 

Google Scholar 
Oro, D. & Pradel, R. Recruitment of audouin’s gull to the ebro delta colony at metapopulation level in the western mediterranean. Mar. Ecol. Prog. Ser. 180, 267–273 (1999).ADS 
Article 

Google Scholar 
Payo-Payo, A. et al. Predator arrival elicits differential dispersal, change in age structure and reproductive performance in a prey population. Sci. Rep. 8, 1–7 (2018).CAS 
Article 

Google Scholar 
Dunshea, G. et al. Telomeres as age markers in vertebrate molecular ecology. Mol. Ecol. Resour. 11, 225–235 (2011).CAS 
PubMed 
Article 

Google Scholar 
De Paoli-Iseppi, R. et al. Age estimation in a long-lived seabird (Ardenna tenuirostris) using dna methylation-based biomarkers. Mol. Ecol. Resour. 19, 411–425 (2019).PubMed 
Article 

Google Scholar 
Horvath, S. Dna methylation age of human tissues and cell types. Genome Biol. 14, 1–20 (2013).Article 

Google Scholar 
Xia, X., Chen, W., McDermott, J. & Han, J.-D.J. Molecular and phenotypic biomarkers of aging. F1000Research 6, 860 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Parrott, B. B. & Bertucci, E. M. Epigenetic aging clocks in ecology and evolution. Trends Ecol. Evol. 34, 767–770 (2019).PubMed 
Article 

Google Scholar 
Monaghan, P., Maklakov, A. A. & Metcalfe, N. B. Intergenerational transfer of ageing: Parental age and offspring lifespan. Trends Ecol. Evol. 35, 927–937 (2020).PubMed 
Article 

Google Scholar 
Moorad, J. A. & Nussey, D. H. Evolution of maternal effect senescence. Proc. Natl. Acad. Sci. 113, 362–367 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Ivimey-Cook, E. & Moorad, J. The diversity of maternal-age effects upon pre-adult survival across animal species. Proc. R. Soc. B 287, 20200972 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Krist, M. Egg size and offspring quality: A meta-analysis in birds. Biol. Rev. 86, 692–716 (2011).PubMed 
Article 

Google Scholar 
Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).MATH 
Article 

Google Scholar 
Cutler, D. R. et al. Random forests for classification in ecology. Ecology 88, 2783–2792 (2007).PubMed 
Article 

Google Scholar 
Oro, D., Hernández, N., Jover, L. & Genovart, M. From recruitment to senescence: Food shapes the age-dependent pattern of breeding performance in a long-lived bird. Ecology 95, 446–457 (2014).PubMed 
Article 

Google Scholar 
Lemaitre, J.-F. & Gaillard, J.-M. Reproductive senescence: New perspectives in the wild: Reproductive senescence in the wild. Biol. Rev. 92, 2182–2199 (2017).PubMed 
Article 

Google Scholar 
Forslund, P. & Part, T. Age and reproduction in birds-hypotheses and tests. Trends Ecol. Evol. 10, 374–378 (1995).CAS 
PubMed 
Article 

Google Scholar 
Reid, J. M., Bignal, E. M., Bignal, S., McCracken, D. I. & Monaghan, P. Age-specific reproductive performance in red-billed choughs Pyrrhocorax pyrrhocorax: Patterns and processes in a natural population. J. Anim. Ecol. 72, 765–776 (2003).Article 

Google Scholar 
Cam, E., Aubry, L. M. & Authier, M. The conundrum of heterogeneities in life history studies. Trends Ecol. Evol. 31, 872–886 (2016).PubMed 
Article 

Google Scholar 
Hamel, S., Côté, S. D., Gaillard, J.-M. & Festa-Bianchet, M. Individual variation in reproductive costs of reproduction: High-quality females always do better. J. Anim. Ecol. 78, 143–151 (2009).PubMed 
Article 

Google Scholar 
Hamel, S., Gaillard, J.-M. & Yoccoz, N. G. Introduction to: Individual heterogeneity—the causes and consequences of a fundamental biological process (2018).Gimenez, O., Cam, E. & Gaillard, J.-M. Individual heterogeneity and capture-recapture models: What, why and how?. Oikos 127, 664–686 (2018).Article 

Google Scholar 
Räsänen, K. & Kruuk, L. E. B. Maternal effects and evolution at ecological time-scales. Funct. Ecol. 21, 408–421 (2007).Article 

Google Scholar 
Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: The role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010).PubMed 
Article 

Google Scholar 
Mills, J. A. et al. Archiving primary data: Solutions for long-term studies. Trends Ecol. Evol. 30, 581–589 (2015).PubMed 
Article 

Google Scholar 
Conde, D. A. et al. Data gaps and opportunities for comparative and conservation biology. Proc. Natl. Acad. Sci. 116, 9658–9664 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Beamonte-Barrientos, R., Velando, A., Drummond, H. & Torres, R. Senescence of maternal effects: Aging influences egg quality and rearing capacities of a long-lived bird. Am. Nat. 175, 469–480 (2010).PubMed 
Article 

Google Scholar 
Descamps, S. et al. Detecting population heterogeneity in effects of north Atlantic oscillations on seabird body condition: Get into the rhythm. Oikos 119, 1526–1536 (2010).Article 

Google Scholar 
Froy, H., Phillips, R. A., Wood, A. G., Nussey, D. H. & Lewis, S. Age-related variation in reproductive traits in the wandering albatross: Evidence for terminal improvement following senescence. Ecol. Lett. 16, 642–649 (2013).PubMed 
Article 

Google Scholar 
Giron, D. & Casas, J. Mothers reduce egg provisioning with age. Ecol. Lett. 6, 273–277 (2003).Article 

Google Scholar 
Lecomte, V. J. et al. Patterns of aging in the long-lived wandering albatross. Proc. Natl. Acad. Sci. 107, 6370–6375 (2010).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Massot, M. et al. An integrative study of ageing in a wild population of common lizards: Integrative study of ageing. Funct. Ecol. 25, 848–858 (2011).Article 

Google Scholar 
Nussey, D. H. et al. Inter and intrasexual variation in aging patterns across reproductive traits in a wild red deer population. Am. Nat. 174, 342–357 (2009).PubMed 
Article 

Google Scholar 
Vega-Trejo, R., Kruuk, L. E. B., Jennions, M. D. & Head, M. L. What happens to offspring when parents are inbred, old or had a poor start in life? Evidence for sex-specific parental effects. J. Evol. Biol. 31, 1138–1151 (2018).PubMed 
Article 

Google Scholar 
Attard, M. R. G. et al. A new, three-dimensional geometric morphometric approach to assess egg shape. PeerJ 6, e5052 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Holveck, M.-J., Guerreiro, R., Perret, P., Doutrelant, C. & Gregoire, A. Eggshell coloration indicates female condition during egg-laying: A field experiment in blue tits. Biol. J. Lin. Soc. 128, 181–200 (2019).Article 

Google Scholar 
Siefferman, L., Navara, K. J. & Hill, G. E. Egg coloration is correlated with female condition in eastern bluebirds (Sialia sialis). Behav. Ecol. Sociobiol. 59, 651–656 (2006).Article 

Google Scholar 
Valcu, C.-M. et al. Life history shapes variation in egg composition in the blue tit Cyanistes caeruleus. Commun. Biol. 2, 1–14 (2019).Article 

Google Scholar 
Genovart, M., Oro, D. & Tenan, S. Immature survival, fertility, and density dependence drive global population dynamics in a long-lived species. Ecology 99, 2823–2832 (2018).CAS 
PubMed 
Article 

Google Scholar 
Fernandez-Chacon, A. et al. When to stay, when to disperse and where to go: Survival and dispersal patterns in a spatially structured seabird population. Ecography 36, 1117–1126 (2013).Article 

Google Scholar 
Tavecchia, G., Pradel, R., Genovart, M. & Oro, D. Density-dependent parameters and demographic equilibrium in open populations. Oikos 116, 1481–1492 (2007).Article 

Google Scholar 
Oro, D. Perturbation, Behavioural Feedbacks, and Population Dynamics in Social Animals: When to Leave and Where to Go (Oxford Univ Press, 2020).Book 

Google Scholar 
Ferrer, M. & Penteriani, V. A process of pair formation leading to assortative mating: Passive age-assortative mating by habitat heterogeneity. Anim. Behav. 66, 137–143 (2003).Article 

Google Scholar 
Ludwig, S. C. & Becker, P. H. Supply and demand: Causes and consequences of assortative mating in common terns sterna hirundo. Behav. Ecol. Sociobiol. 62, 1601–1611 (2008).Article 

Google Scholar 
Oro, D., Jover, L. & Ruiz, X. Influence of trawling activity on the breeding ecology of a threatened seabird, audouin’s gull Larus audouinii. Mar. Ecol. Prog. Ser. 139, 19–29 (1996).ADS 
Article 

Google Scholar 
Hoyt, D. F. Practical methods of estimating volume and fresh weight of bird eggs. Auk 96, 73–77 (1979).
Google Scholar 
Oro, D., Ruiz, X., Jover, L., Pedrocchi, V. & González-Solís, J. Diet and adult time budgets of audouin’s gull Larus audouinii in response to changes in commercial fisheries. Ibis 139, 631–637 (1997).Article 

Google Scholar 
Hurrell, J. W., Kushnir, Y., Ottersen, G. & Visbeck, M. An overview of the north atlantic oscillation. In Geophysical Monograph Series Vol. 134 (eds Hurrell, J. W. et al.) (American Geophysical Union, 2003).R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2013).RStudio Team. RStudio: Integrated Development Environment for R (RStudio, PBC, 2022).Breiman, L. (ed.) Classification and Regression Trees repr. (Chapman and Hall, 1998).
Google Scholar 
Liaw, A. & Wiener, M. Classification and regression by randomforest. R News 2, 18–22 (2002).
Google Scholar 
Wickham, H., François, R., Henry, L. & Müller, K. dplyr: A Grammar of Data Manipulation (2018). R package version 0.7.6.Kuhn, M. The caret package. J. Stat. Softw. 28, 25 (2012).
Google Scholar  More

Bonan, G. B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449. https://doi.org/10.1126/science.1155121 (2008).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Le Quere, C. et al. Global carbon budget 2016. Earth Syst. Sci. Data 8, 605–649. https://doi.org/10.5194/essd-8-605-2016 (2016).ADS 
Article 

Google Scholar 
DavidMorales-Hidalgo, D., Oswalt, S. N. & Somanathan, E. Status and trends in global primary forest, protected areas, and areas designated for conservation of biodiversity from the Global Forest Resources Assessment 2015. Forest Ecol. Manag. 352, 68–77. https://doi.org/10.1016/j.foreco.2015.06.011 (2015).Article 

Google Scholar 
Kauppi, P. E., Sandstrom, V. & Lipponen, A. Forest resources of nations in relation to human well-being. PLoS One 13, e0196248. https://doi.org/10.1371/journal.pone.0196248 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Anderegg, W. R. L. et al. Climate-driven risks to the climate mitigation potential of forests. Science 368, 1327. https://doi.org/10.1126/science.aaz7005 (2020).CAS 
Article 

Google Scholar 
Wilson, M. C. et al. Habitat fragmentation and biodiversity conservation: Key findings and future challenges. Landsc. Ecol. 31, 219–227. https://doi.org/10.1007/s10980-015-0312-3 (2016).Article 

Google Scholar 
Haddad, N. M. et al. Habitat fragmentation and its lasting impact on earth’s ecosystems. Sci. Adv. 1, e1500052. https://doi.org/10.1126/sciadv.1500052 (2015).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Holl, K. D. Restoring tropical forests from the bottom up. Science 355, 455–456. https://doi.org/10.1126/science.aam5432 (2017).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Audino, L. D., Murphy, S. J., Zambaldi, L., Louzada, J. & Comita, L. S. Drivers of community assembly in tropical forest restoration sites: Role of local environment, landscape, and space. Ecol. Appl. 27, 1731–1745. https://doi.org/10.1002/eap.1562 (2017).Article 
PubMed 

Google Scholar 
Temperton, V. M., Hobbs, R. J., Nuttle, T. & Halle, S. in Assembly Rules and Restoration Ecology: Bridging the Gap Between Theory and Practice [Science and Practice of Ecological Restoration]. i–xv, 1–439 (2004).Young, T. P., Chase, J. M. & Huddleston, R. T. Community succession and assembly: Comparing, contrasting and combining paradigms in the context of ecological restoration. Ecol. Restor. 19, 5–18 (2001).Article 

Google Scholar 
Vellend, M. The Theory of Ecological Communities (Princeton University Press, 2016).
Google Scholar 
HilleRisLambers, J., Adler, P. B., Harpole, W. S., Levine, J. M. & Mayfield, M. M. Rethinking community assembly through the lens of coexistence theory. Annu. Rev. Ecol. Evol. Syst. 43(43), 227–248. https://doi.org/10.1146/annurev-ecolsys-110411-160411 (2012).Article 

Google Scholar 
Connell, J. H. On the role of natural enemies in preventing competitive exclusion in some marine animals and in rain forest trees. In Dynamics of Populations (eds Den Boer, P. J. & Gradwell, G. R.) (Centre for Agricultural Publishing and Documentation, 1971).
Google Scholar 
Janzen, D. H. Herbivores and the number of tree species in tropical forests. Am. Nat. 104, 501. https://doi.org/10.1086/282687 (1970).Article 

Google Scholar 
Schmid, J. S., Taubert, F., Wiegand, T., Sun, I. F. & Huth, A. Network science applied to forest megaplots: Tropical tree species coexist in small-world networks. Sci. Rep. https://doi.org/10.1038/s41598-020-70052-8 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Wang, H. X. et al. Prevalence of inter-tree competition and its role in shaping the community structure of a natural Mongolian scots pine (Pinus sylvestris var. mongolica) forest. Forests https://doi.org/10.3390/f8030084 (2017).Article 

Google Scholar 
Hubbell, S. P. et al. Light-gap disturbances, recruitment limitation, and tree diversity in a neotropical forest. Science 283, 554–557. https://doi.org/10.1126/science.283.5401.554 (1999).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Janik, D. et al. Breaking through beech: A three-decade rise of sycamore in old-growth European forest. Forest Ecol. Manag. 366, 106–117. https://doi.org/10.1016/j.foreco.2016.02.003 (2016).Article 

Google Scholar 
Svatek, M., Rejzek, M., Kvasnica, J., Repka, R. & Matula, R. Frequent fires control tree spatial pattern, mortality and regeneration in argentine open woodlands. Forest Ecol. Manag. 408, 129–136. https://doi.org/10.1016/j.foreco.2017.10.048 (2018).Article 

Google Scholar 
Giammarchi, F. et al. Effects of the lack of forest management on spatiotemporal dynamics of a subalpine Pinus cembra forest. Scand. J. Forest Res. 32, 142–153. https://doi.org/10.1080/02827581.2016.1207802 (2017).Article 

Google Scholar 
Janik, D. et al. Patterns of Fraxinus angustifolia in an alluvial old-growth forest after declines in flooding events. Eur. J. Forest Res. 135, 215–228. https://doi.org/10.1007/s10342-015-0925-8 (2016).Article 

Google Scholar 
Bagchi, R. et al. Defaunation increases the spatial clustering of lowland western amazonian tree communities. J. Ecol. 106, 1470–1482. https://doi.org/10.1111/1365-2745.12929 (2018).Article 

Google Scholar 
Zhang, L. Y., Dong, L. B., Liu, Q. & Liu, Z. G. Spatial patterns and interspecific associations during natural regeneration in three types of secondary forest in the central part of the greater Khingan mountains, Heilongjiang province, China. Forests https://doi.org/10.3390/f11020152 (2020).Article 

Google Scholar 
Obiang, N. L. E. et al. Determinants of spatial patterns of canopy tree species in a tropical evergreen forest in Gabon. J. Veg. Sci. 30, 929–939. https://doi.org/10.1111/jvs.12778 (2019).Article 

Google Scholar 
Wiegand, T. et al. Spatially explicit metrics of species diversity, functional diversity, and phylogenetic diversity: Insights into plant community assembly processes. Annu. Rev. Ecol. Evol. Syst. 48(48), 329–351. https://doi.org/10.1146/annurev-ecolsys-110316-022936 (2017).Article 

Google Scholar 
Gabriel, E. Spatial point patterns: Methodology and applications with R. Math. Geosci. 49, 815–817. https://doi.org/10.1007/s11004-016-9670-x (2017).CAS 
Article 
MATH 

Google Scholar 
Baddeley, A., Rubak, R. & Turner, R. Spatial Point Patterns, Methodology and Applications with R (CRC Press, 2016).MATH 

Google Scholar 
Wiegand, T. & Moloney, K. A. Rings, circles, and null-models for point pattern analysis in ecology. Oikos 104, 209–229. https://doi.org/10.1111/j.0030-1299.2004.12497.x (2004).Article 

Google Scholar 
Plotkin, J. B., Chave, J. M. & Ashton, P. S. Cluster analysis of spatial patterns in Malaysian tree species. Am. Nat. 160, 629–644. https://doi.org/10.1086/342823 (2002).Article 
PubMed 

Google Scholar 
Ripley, B. D. Modeling spatial patterns. J. R. Stat. Soc. B 39, 172–212 (1977).
Google Scholar 
He, F. L. & Gaston, K. J. Estimating species abundance from occurrence. Am. Nat. 156, 553–559. https://doi.org/10.1086/303403 (2000).Article 
PubMed 

Google Scholar 
Diggle, P. Statistical Analysis of Spatial Point Patterns (Academic Press, 1983).MATH 

Google Scholar 
Pielou, E. C. The use of point-to-plant distances in the study of the pattern of plant-populations. J. Ecol. 47, 607–613. https://doi.org/10.2307/2257293 (1959).Article 

Google Scholar 
Losapio, G., Montesinos-Navarro, A. & Saiz, H. Perspectives for ecological networks in plant ecology. Plant Ecol. Divers. 12, 87–102. https://doi.org/10.1080/17550874.2019.1626509 (2019).Article 

Google Scholar 
Fuller, M. M., Wagner, A. & Enquist, B. J. Using network analysis to characterize forest structure. Nat. Resour. Model. 21, 225–247. https://doi.org/10.1111/j.1939-7445.2008.00004.x (2008).MathSciNet 
Article 

Google Scholar 
Montoya, J. M., Pimm, S. L. & Sole, R. V. Ecological networks and their fragility. Nature 442, 259–264. https://doi.org/10.1038/nature04927 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Proulx, S. R., Promislow, D. E. L. & Phillips, P. C. Network thinking in ecology and evolution. Trends Ecol. Evol. 20, 345–353. https://doi.org/10.1016/j.tree.2005.04.004 (2005).Article 
PubMed 

Google Scholar 
Nakagawa, Y., Yokozaw, M. & Hara, T. Complex network analysis reveals novel essential properties of competition among individuals in an even-aged plant population. Ecol. Complex 26, 95–116. https://doi.org/10.1016/j.ecocom.2016.03.005 (2016).Article 

Google Scholar 
Wiegand, T. & Moloney, K. A. Handbook of Spatial Point Pattern Analysis in Ecology (CRC Press, 2013).Book 

Google Scholar 
Barthelemy, M. Spatial networks. Phys. Rep. Rev. Sect. Phys. Lett. 499, 1–101. https://doi.org/10.1016/j.physrep.2010.11.002 (2011).MathSciNet 
CAS 
Article 

Google Scholar 
Keren, S. Modeling tree species count data in the understory and canopy layer of two mixed old-growth forests in the Dinaric region. Forests https://doi.org/10.3390/f11050531 (2020).Article 

Google Scholar 
Podlaski, R. Models of the fine-scale spatial distributions of trees in managed and unmanaged forest patches with Abies alba Mill. and Fagus sylvatica L. Forest Ecol. Manag. 439, 1–8 (2019).Article 

Google Scholar 
Levin, S. A. Theoretical ecology—Principles and applications, 3rd edition. Science 316, 1699–1700. https://doi.org/10.1126/science.1141870 (2007).CAS 
Article 

Google Scholar 
Martinez-Lopez, V., Garcia, C., Zapata, V., Robledano, F. & De la Rua, P. Intercontinental long-distance seed dispersal across the Mediterranean basin explains population genetic structure of a bird-dispersed shrub. Mol. Ecol. 29, 1408–1420. https://doi.org/10.1111/mec.15413 (2020).Article 
PubMed 

Google Scholar 
Dale, M. R. T. & Fortin, M. J. From graphs to spatial graphs. Annu. Rev. Ecol. Evol. Syst. 41, 21–38. https://doi.org/10.1146/annurev-ecolsys-102209-144718 (2010).Article 

Google Scholar 
Silva, C. A. et al. Treetop: A shiny-based application and R package for extracting forest information from LiDAR data for ecologists and conservationists. Methods Ecol. Evol. 13, 1164–1176. https://doi.org/10.1111/2041-210x.13830 (2022).Article 

Google Scholar 
Tatsumi, S., Yamaguchi, K. & Furuya, N. Forestscanner: A mobile application for measuring and mapping trees with LiDAR-equipped iPhone and iPad. Methods Ecol. Evol. https://doi.org/10.1111/2041-210x.13900 (2022).Article 

Google Scholar 
Ferraz, A., Saatchi, S. S., Longo, M. & Clark, D. B. Tropical tree size-frequency distributions from airborne LiDAR. Ecol. Appl. 30, e02154. https://doi.org/10.1002/eap.2154 (2020).Article 
PubMed 

Google Scholar 
Bianchi, E., Bugmann, H., Hobi, M. L. & Bigler, C. Spatial patterns of living and dead small trees in subalpine Norway spruce forest reserves in Switzerland. Forest Ecol. Manag. 494, 119315. https://doi.org/10.1016/j.foreco.2021.119315 (2021).Article 

Google Scholar 
Tatsumi, S., Owari, T., Yin, M. F. & Ning, L. Z. Neighborhood analysis of underplanted Korean pine demography in larch plantations: Implications for uneven-aged management in northeast china. Forest Ecol. Manag. 322, 10–18. https://doi.org/10.1016/j.foreco.2014.03.022 (2014).Article 

Google Scholar 
Cornett, M. W., Reich, P. B. & Puettmann, K. J. Canopy feedbacks and microtopography regulate conifer seedling distribution in two Minnesota conifer-deciduous forests. Ecoscience 4, 353–364. https://doi.org/10.1080/11956860.1997.11682414 (1997).Article 

Google Scholar 
Wang, X. F., Zheng, G., Yun, Z. X. & Moskal, L. M. Characterizing tree spatial distribution patterns using discrete aerial LiDAR data. Remote Sens. Basel 12, 712. https://doi.org/10.3390/rs12040712 (2020).ADS 
Article 

Google Scholar 
Matérn, B. Spatial variation: Stochastic models and their application to some problems in forest surveys and other sampling investigations. Meddelanden från Statens Skogsforskningsinstitut 49, 1–144 (1960).MathSciNet 

Google Scholar 
Matérn, B. Spatial Variation. Lecture Notes in Statistics Vol. 36 (Springer, 1986).Book 

Google Scholar 
Thomas, M. A generalisation of Poisson’s binomial limit for use in ecology. Biometrika 36, 18–25 (1949).MathSciNet 
CAS 
Article 

Google Scholar 
Lotwick, H. W. Simulation of some spatial hard core models, and the complete packing problem. J. Stat. Comput. Simul. 15, 295–314 (1982).MathSciNet 
Article 

Google Scholar 
Strauss, D. J. A model for clustering. Biometrika 62, 467–475 (1975).MathSciNet 
Article 

Google Scholar 
Cressie Noel, A. C. Statistics for Spatial Data (Wiley-Interscience, 1993).Book 

Google Scholar 
Besag, J. E. Contribution to the discussion of the paper by Ripley. J. R. Stat. Soc. 39, 193–195 (1977).MathSciNet 

Google Scholar  More