Ecology

Mistakes in a landmark paper that reported the first orangutan genomes might have implications for breeding programmes.Credit: Fiona Rogers/Nature Picture Library

Susie the Sumatran orangutan was a genetic pioneer — the first of her species to have her genome fully sequenced. Her genetic library, and that of ten other orangutans, appeared in a landmark paper in Nature in 20111 that has underpinned hundreds of subsequent studies.But in August, researchers revealed that eight of the sequences in this paper had mistakenly been assigned to the wrong orangutans2. Nature issued a correction from the authors of the original paper3.The scale of the errors sparked ire on social media, and some scientists have warned that the mistakes could have repercussions for orangutan breeding programmes. “Well that’s a bit of a f&£k up orang-utan genome researchers — only mildly embarrassing guys and girls”, tweeted Michael Sweet, a molecular ecologist at the University of Derby, UK.
Africa: sequence 100,000 species to safeguard biodiversity
It’s not clear how these swapped identities have affected orangutan research. But researchers involved in the new analysis believe the discovery might highlight how issues in the scientific community — including the pressure to publish and a reliance on peer review to catch mistakes — could allow such errors to slip into the scientific record.“I think there are errors like this in many, many published papers,” says Graham Banes, an evolutionary biologist formerly at the University of Wisconsin–Madison who led the reanalysis of the 2011 paper. “In some ways, we’re lucky that this was just orangutans. What if this was a biomedical paper and people were developing therapies based on published data?”“It’s fairly easy for these things to occur,” adds Robert Fulton, a genomic scientist at Washington University School of Medicine in St Louis, Missouri, who was part of the team behind the original paper and is a co-author on the reanalysis. “What’s important is that that the data are now correct.” Devin Locke, who led the preparation of the 2011 paper and was formerly a colleague of Fulton’s at Washington University, did not respond to questions about the work.Hybrid headacheDetailed ‘reference’ genomes, such as those published in the 2011 Nature paper, are a key tool for biologists. In 2017, Banes and his team were using the genomes to study what happens when different species of orangutan interbreed, a process called hybridization.They noticed that the names given to some of the samples didn’t match the animals’ reported sex. For example, the 2011 paper reported that an orangutan named Dolly was male. But according to the orangutan studbook — a record of orangutans living in zoos — Dolly was female. Even stranger, Banes found that some of the genomes marked as male lacked a Y chromosome. “There was just this series of things that didn’t make sense,” he recalls.
Major wildlife report struggles to tally humanity’s exploitation of species
Banes and his colleagues eventually found that the 2011 paper had misidentified all but two of the orangutan genomes. Some mistakes seem to be the result of typos. In one case, a sample from a male orangutan was given an ID number that actually corresponded to a sample from an African pig in a tissue repository. Other samples seem to have had their identities swapped during laboratory work. The 2011 study helped to pin down when Bornean and Sumatran orangutans split into separate species, and compared their genomes with those of other primates. These conclusions are largely uncompromised by the mix-up. But Banes says that the errors could have implications for other research, including his own.Banes uses genetic data to provide zoos with recommendations about their captive breeding programmes. Zoos try to avoid crossbreeding orangutan species, partly to mimic wild populations and also because hybrids can suffer high rates of miscarriage and birth defects, says Banes. While re-examining the samples from the 2011 paper, the team realized that one of the sequences thought to be Sumatran (Pongo abelii) was actually Tapanuli (Pongo tapanuliensis), a third species of orangutan that was only described in 20174.Unfortunately, the 2011 paper had wrongly assigned the Tapanuli genome to Baldy, a male orangutan, rather than its actual owner, a female orangutan named Bubbles (both are now dead). Banes says that his team came “perilously close” to announcing in a paper that Baldy was Tapanuli.Although Baldy has no living descendants, Bubbles has several offspring at zoos around the world, all of which are Sumatran–Tapanuli hybrids. Zookeepers will now have to decide whether to stop breeding Bubbles’ descendants to avoid further hybridization, says Vincent Nijman, an anthropologist at Oxford Brookes University, UK.‘Bigger concerns’However, Nijman also argues that the errors will have little effect on orangutan conservation as a whole. Zoos often bill their animals as a back-up for endangered species, but conservationists are much more focused on the thousands of orangutans in the wild that are threatened by deforestation. “I think we have bigger concerns than some mixed-up samples,” says Erik Meijaard, a conservation scientist at Borneo Futures, a conservation consultancy company based in Bandar Seri Begawan, Brunei.
Scientists warn deal to save biodiversity is in jeopardy
Michael Krützen, an evolutionary geneticist at the University of Zurich in Switzerland, agrees that although the errors are “annoying”, their impact on downstream research is probably minimal. However, he says that the problems might be an example of how academia’s publish-or-perish environment could lead to “sloppy” work, as researchers race to publish their work in high-tier journals.Banes agrees that this kind of pressure — along with an over-reliance on a peer-review system that does not offer its volunteer reviewers tangible financial or professional benefits — could lead to errors slipping into published manuscripts.A spokesperson for Nature declined to comment on why the errors in the 2011 paper were not caught by peer review, citing concerns about confidentiality. (Nature’s news team is editorially independent of its academic publishing operation). “However, we would like to stress that we take our responsibility to maintain the accuracy of the scientific record very seriously,” they wrote in an e-mail. “If issues are raised about any paper we have published, we will look into them carefully and update the literature where appropriate.”Banes says that it’s important not to blame individual scientists for such errors, not least because it could discourage efforts to correct mistakes in future. “I think any scientist could have made these mistakes,” he says. “But if we all jump out and say, ‘oh my god, how could they have been so stupid?’, no one is ever going to correct anything. That shame is detrimental to science.” More

Sequencing and assembly of the H. schachtii genomeWe measured (Supplemental Fig. 1), sequenced (BioProject PRJNA722882), and assembled the genome of H. schachtii (population Bonn) using a combination of flow cytometry, Pacific Biosciences sequencing, and Illumina sequencing. H. schachtii has the largest genome (160–170 Mb) of any cyst nematode measured/sequenced to date (Supplementary Table 1). It was sequenced to 192-fold coverage using Pacific Biosciences sequencing (fragment n50 of 16 kb), and 144-fold coverage using Illumina sequencing (150 bp Paired-end reads). The final, polished, contamination-free (Supplemental Fig. 2), assembly (v1.2) included ~179 Mbp contained within 395 scaffolds: 90% of the sequence is contained on scaffolds longer than 281,463 bp (n = 154). The assembly is a largely complete haploid representation of the diploid genome, as evidenced by core eukaryotic genes being largely present, complete and single copy (CEGMA 93.15% complete with an average of 1.12 copies each, and BUSCO (Eukaryota odb9) 79% complete with 8.2% duplicated—Supplementary Table 2). Over three million variants were phased into haplotypes (2029 blocks, N50 239.5 kb, covering 94.7% of the reference) which can be used to predict true protein variants (Supplementary data 1), and 601 larger structural variants were identified (Supplementary data 2).The trans-kingdom, lifestage-specific, transcriptomes of H. schachtii and A. thaliana provide a holistic view of parasitismWe devised a sampling procedure to cover all major life stages/transitions of the parasitic life cycle to generate a simultaneous, chronological, and comprehensive picture of nematode gene expression, and infection-site-specific plant gene expression patterns. We sampled cysts and pre-infective second-stage juveniles (J2s), as well as infected segments of A. thaliana root and uninfected adjacent control segments of root at 10 hours post infection (hpi – migratory J2s, pre-establishment of the feeding site), 48 hpi (post establishment of the feeding site), 12 days post infection females (dpi – virgin), 12 dpi males (differentiated, pre-emergence, most if not all stopped feeding), and 24 dpi females (post mating), each in biological triplicate (Fig. 1A). We generated approximately nine billion pairs of 150 bp strand-specific RNAseq reads (Supplementary data 3) covering each stage in biological triplicate (for the parasite and the host): in the early stages of infection we generated over 400 million reads per replicate, to provide sufficient coverage of each kingdom.Fig. 1: Trans-kingdom, lifestage-specific, transcriptome of H. schachtii and A. thaliana.A Schematic representation of the life cycle of H. schachtii infecting A. thaliana, highlighting the 7 stages sampled in this study. For each stage, the average number of trimmed RNAseq read pairs per replicate is shown, with the proportion of reads mapping to either parasite or host in parentheses. B Principle components 1 and 2 for H. schachtii and A. thaliana expression data are plotted. Arrows indicate progression through the life cycle/real-time. Hours post infection (hpi), days post infection (dpi).Full size imageStrand-specific RNAseq reads originating from host and parasite were deconvoluted by mapping to their respective genome assemblies (H. schachtii v.1.2 and TAIR10). For the parasite, ~500 million Illumina RNAseq read pairs uniquely mapping to the H. schachtii genome were used to generate a set of 26,739 gene annotations (32,624 transcripts – detailed further in the next section), ~77% of which have good evidence of transcription in at least one lifestage (≥10 reads in at least one rep). Similarly for the host, ~2.8 billion Illumina RNAseq read pairs uniquely mapping to the A. thaliana genome show that ~77% of the 32,548 gene models have good evidence of transcription in at least one stage (≥10 reads in at least one rep, even though we only sampled roots). A principal component analysis of the host and parasite gene expression data offers several insights into the parasitic process. Principle component 1 (60% of the variance) and 2 (19% of the variance) of the parasite recapitulate the life cycle in PCA space (Fig. 1B). The 12 dpi female transcriptome is more similar to the 24 dpi female transcriptome than to the 12 dpi male transcriptome. Principle components 1 (75% of the variance) and 2 (10% of the variance) of the host show that the greatest difference between infected and uninfected plant tissue is at the early time points (10 hpi), and that the transcriptomes of infected and uninfected plant material converge over time, possibly due to systemic effects of infection. A 12 dpi male syncytium transcriptome is roughly intermediate between a control root transcriptome and a 12 dpi female syncytium transcriptome. Given that at this stage most if not all of the males will have ceased feeding, this could be due to inadequate formation of the feeding site, or regression of the tissue. In any case, by comparing both principal component analyses, we can see that what is a relatively small difference in the transcriptomes of the feeding sites of males and females is amplified to a relatively large difference in the transcriptomes of the males and females themselves (Fig. 1B).The consequences, and possible causes, of large-scale segmental duplication in the Heterodera lineageTo understand the evolutionary origin(s) of the relatively large number of genes in H. schachtii in particular, and Heterodera spp. in general, we analysed the abundance and categories of gene duplication in the predicted exome. Compared to a related cyst nematode, Globodera pallida (derived using comparable methodology and of comparable contiguity) the exomes of H. schachtii and H. glycines are characterised by a relatively smaller proportion of single-copy genes (as classified by MCSanX toolkit17, and a relatively greater proportion of segmental duplications (at least five co-linear genes with no >25 genes between them), with relatively similar proportions of dispersed duplications (two similar genes with >20 other genes between them), proximal duplications (two similar genes with  +0.5 or  More

Insect rearingThe cotton mealybugs used in this study were originally collected from Rose of Sharon, Hibiscus syriacus L. (Malvales: Malvaceae) in Jinhua, Zhejiang Province, China, in June 2016. They were maintained on fresh tomato plants (cv. Hezuo-903, Shanghai Changzhong Seeds Industry Co., Ltd, China) in a climatically controlled chamber maintained at 27 ± 1 °C, 75% relative humidity (RH), and a photoperiod of 14:10 (L:D). For detailed insect rearing and tomato cultivation methods see ref. 56.Scanning electron microscopy (SEM) of P. solenopsis waxSEM was used to observe changes in wax on the body surface of adult P. solenopsis females according to the methods of Huang et al.57. Briefly, collected insects were taped onto a stub and dried in an ion sputter (Hatachi, Tokyo, Japan) under a vacuum. After gold sputtering, the samples were observed using a TM-1000 SEM (Hatachi, Tokyo, Japan). Photos were scanned from the dorsal part of the third thoracic segment. Thirty insects were used for both RNAi-treated and control groups.Chemical composition analysis of mealybug waxA small soft brush was used to collect wax filaments from the body surface of P. solenopsis females. Prior to use, the brush was washed successively by 70% ethanol, sterile water, and 1× sterile phosphate-buffered saline (PBS, pH 7.4). The wax was collected into a clean chromatography vial for the following experiments. Two vials of wax, each collected from 1000 adult females, were dissolved in 1 ml of methanol and 1 ml of n-hexane, respectively. The vials were stirred gently for 3 min, kept at room temperature for 30 min, and then put into an S06H ultrasonic vibrator (Zealway, Xiamen, China) for 30 min to dissolve the wax sufficiently. The samples were analyzed on a TRACE 1310 (Thermo Scientific, Waltham, USA) gas chromatograph (GC) equipped with an ISQ single quadrupole MS and interfaced with the Chromeleon 7.2 data analysis system (Thermo Scientific, Waltham, USA), with a constant flow of helium at 1 ml/min. For each sample, a splitless injection of 1.0 μl was respectively made into a polar TG-WaxMS (Thermo Scientific, Waltham, USA) and a nonpolar TG-5MS (Thermo Scientific, Waltham, USA) 30 m × 0.25 mm × 0.25 μm capillary column. The temperature program for polar column samples was as follows: 40 °C for 2 min, then 5 °C/min to 240 °C, hold 10 min; the program for nonpolar column samples was: 40 °C for 2 min, then 5 °C/min to 300 °C, hold 5 min. Injector and detector temperatures were, respectively, set at 250 and 230 °C for polar column samples, and at 300 and 300 °C for nonpolar column samples. Mass detection for all samples was run under an EI mode with a 70 eV ionization potential and an effective m/z range of 35–450 at a scan rate of 5 scan/s. Chemical compounds were identified by mapping against the NIST database. The relative content of each compound was calculated by peak area which was determined using the Agilent MassHunter system.RNA extraction and RT-qPCRTotal RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA) following the manufacturer’s instructions, and RNA quality was accessed using agarose gel electrophoresis and a Biodrop μLite. 800 ng of total RNA was used for cDNA synthesis using the HiScript III RT SuperMixfor qPCR (+gDNA wiper) (Vazyme Biotech Co., Ltd., Nanjing, China), according to the manufacturer’s instructions. Quantitative RT-PCR (RT-qPCR) was conducted using an AriaMx real-time PCR system (Agilent Technologies, USA), using a 20 μl reaction containing 2 μl of 10-fold diluted cDNA, 0.8 μl of each primer, and 10 μl ChamQ SYBR Color qPCR Master Mix (Vazyme Biotech Co., Ltd., Nanjing, China). The RT-qPCR thermocycling protocol was 95 °C for 30 s, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The PsActin gene was used as an internal control. At least three biological replicates were used for each experiment. Quantitative variations were evaluated using the relative quantitative method (2−ΔΔCt)58.Transcriptome analysis of integumentary and non-integumentary tissuesTo obtain the integument and other tissues, adult P. solenopsis females were dissected in 1× sterile PBS (pH 7.4) on a sterile Petri dish. Dissected fresh tissues were directly used or frozen in liquid nitrogen and stored at −80 °C for follow-up experiments. We sequenced the transcriptomes of integumentary and non-integumentary tissues (all other tissues without integument) dissected from 150 adult females, with each sample being repeated in triplicate. mRNAs were purified from total RNA via oligo (dT) magnetic beads, and the fragmented mRNAs were then reverse transcribed into cDNA using random primers. Constructed pair-end libraries were sequenced using an Illumina HiSeq X Ten platform in Novogene (Beijing, China). After quality control, the clean RNA-Seq data of the six libraries were aligned with the P. solenopsis genome (http://v2.insect-genome.com/Organism/624) using HISTAT259. Then featureCounts60 and DESeq261 were used for the differential expression analysis of genes. The threshold for differentially expressed genes (DEGs) was defined by log2fold ≥ 1 or ≤−1 and a padj-value  More

Retrieval of E. coli plasmid sequencesAll E. coli sequences were downloaded from the NCBI FTP server in May 2020. To establish an initial collection of plasmids, only complete genomes with an associated plasmid were retained. All genomes were verified for belonging to the species E. coli using kmerfinder (https://cge.cbs.dtu.dk/services/KmerFinder/). Sequence type (ST) was determined via multi-locus sequence typing (MLST) based on the 7-gene Achtman scheme using pubMLST (https:/github.com/tseemann/mlst). Only genomes with exact matches were assigned for each ST and used for subsequent analysis. To ensure our sequences were sufficiently representative of E. coli pathogens expected in nature, a systematic literature search (see description below and Fig. S1) was conducted to establish an expected distribution of STs (Table S1). This information was used to update our initial collection to match the top 4 most prevalent STs (131, 11, 73, and 95). Specifically, to identify supplementary plasmid sequences, genome accession IDs were chosen from EnteroBase based on the following criteria: the strain was matched to the correct ST and had a high-quality genome sequence (based on N50  > 20,000 and the number of contigs  0.1, 2-tailed student t test). For the second method, all kanR plasmids were used, and instead changed the hosts such that DH5αPro cells were in competition with DH5αPro containing a spontaneous rifampicin-resistant mutant (rifR). Any rifR strain was quantified on rifampicin-containing plates, and the second strain was quantified by rifampicin CFU minus CFU obtained on blank plates. We established that rifR exhibited no fitness defects by (1) growth rates between the wild-type (WT) strain (W) and rifR (M) (Fig. S5D), and (2) directly competing the two control strains (Fig. S5E). In both cases, results were statistically indistinguishable (p  > 0.1, two-tailed student t test). KanR/cmR and WT/rifR experiments were each conducted in LB or M9CAG, respectively. In all cases, experiments were repeated with at least three independent biological replicates.Time-kill measurements in the presence of carbenicillinAll strains were grown as previously described. Time-kill experiments entailed hourly measurements of CFU in presence of carbenicillin at either 3.75 μg/mL (3x IC50) or 5 μg/mL (4x IC50) over a span of 2 or 3 h, including time 0. Specifically, overnight cultures were first diluted 1:100 into LB media containing 1 mM IPTG and 50 μg/mL kanamycin and sub-cultured for two hours in a 37 °C incubator with shaking at 250 rpm. Following this, cell density was adjusted as necessary to achieve a starting OD600 of ~0.15 in all cases. Adjusted subcultures were then aliquoted into a 96-well plate and the appropriate carbenicillin treatments were added directly to the well. Plates were sealed with a paper film and placed in a 37 °C incubator with shaking at 250 rpm. Initial collection for time=0 was acquired before carbenicillin treatment. Thereafter, 10 μL of culture was removed from the well every hour, 10-fold serial dilutions were performed and 10 μL was plated on blank LB agar with three technical replicates at each time point. Colonies were counted after plates were grown for 16 h in a 37 °C incubator to determine CFU. This procedure utilized 14 strains of DH5αPro transformed with kanR plasmids of interest – ctrl, katG, lpxM, yfbR, aroH, pld, fdtC, agp, eptC, arcA, argF, mmuM, ahr, and fabG. CFUs were averaged for all technical replicates, and experiments were conducted with at least three independent biological replicates.Oxygen consumption rateOxygen consumption rates (OCR) were obtained with the Resipher device from Lucid Scientific. The selected strains were grown overnight as previously described. Overnight cultures were resuspended in M9CAG media with 1 mM IPTG and 50 μg/mL kanamycin, and placed in 25 °C for one hour to initiate gene expression. Following this, cells were diluted 10x into M9CAG media containing kanamycin and IPTG, and 100 μL was aliquoted per well into a 96-well microtiter plate according to the manufacturer’s instructions. Plates were placed at 30 °C to minimize growth, and oxygen concentration (μM) was measured immediately thereafter. 24 wells were measured consisting of 6 technical replicates for each strain. Given the clear well-well variability (Fig. S8B, C), data shown are for one biological replicate. However, qualitative trends were consistently reproduced in multiple independent experiments.StatisticsIn all cases where t tests and ANOVA’s were used, data was first verified to be normally distributed using Kolmogorov test for normality. Otherwise, Mann-Whitney U-tests were conducted. For panels with multiple tests, Bonferroni correction was used to adjust the p values. To determine whether any metabolic category was significantly dependent on incompatibility groups, we implemented logistic regressions in MATLAB with the function fitglm. Random forest classification was used to establish the relative importance of prevalent metabolic genes and gene categories predicting the presence of antibiotic resistance genes. Chi-square tests were conducted to determine significant co-occurrence of individual antibiotic resistant and metabolism genes. Dissociative relationships were distinguished by the odds ratios from the chi-square tests. To investigate whether the strong associations and disassociations were driven by evolutionary constraints, or simply artifacts of a common ancestor, we re-ran our statistical analysis using Coinfinder [29] to take in our gene presence-absence data, along with the genome phylogeny, and compute the Bonferroni-corrected statistical likelihood of coincidence (either associations or dissociations), thereby accounting for evolutionary relatedness. More

Cline and colleagues analysed spatiotemporal datasets covering 5000 km of popular trout rivers from 1983 to 2017, finding that fishing pressure was four times higher in cold-water sections of rivers than adjacent cool-water sections of rivers, with fisher spending in cold-water sections generating US$500,000 km−1 year−1 and cool-water sections generating US$60,000 km−1 year−1. Overall, 17% and 35% of the current cold-water habitats are projected to be warmer than 18 °C (the threshold for trout thermal extremes) by 2040 and 2080, respectively, with some river sections possibly experiencing habitat losses in excess of 80% by 2080. The combined effects of cold-water habitat loss and increased frequency and severity of drought on fishing pressure could result in 64% declines in fishing river sections by 2040 and 76% declines by 2080. The cumulative impacts of these environmental changes in fishing spending across these rivers could put a total of US$103 million year−1 and US$192 million year−1 at risk by 2040 and 2080, respectively. More

Wurst, S., Sonnemann, I. & Zaller, J. G. Soil Macro-Invertebrates: Their Impact on Plants and Associated Aboveground Communities in Temperate Regions. in Aboveground-Belowground Community Ecology. Ecological Studies 234 (eds. Ohgushi, T., Wurst, S. & Johnson, S. N.) 175–200 (Springer International Publishing, 2018).Wagg, C., Bender, S. F., Widmer, F. & van der Heijden, M. G. A. Soil biodiversity and soil community composition determine ecosystem multifunctionality. PNAS 111, 5266–5270 (2014).ADS 
CAS 
Article 

Google Scholar 
Potapov, A. M. et al. Feeding habits and multifunctional classification of soil-associated consumers from protists to vertebrates. Biol. Rev. 97, 1057–1117 (2022).Article 

Google Scholar 
Potapov, A., Tiunov, A. V. & Scheu, S. Uncovering trophic positions and food resources of soil animals using bulk natural stable isotope composition. Biol. Rev. 94, 37–59 (2019).Article 

Google Scholar 
Voroney, R. P. & Heck, R. J. The Soil Habitat. Soil Microbiology, Ecology and Biochemistry (Elsevier Inc., 2015).Wardle, D. A. The influence of biotic interactions on soil biodiversity. Ecol. Lett. 9, 870–886 (2006).Article 

Google Scholar 
De Deyn, G. B. & Kooistra, L. The role of soils in habitat creation, maintenance and restoration. Philos. Trans. R. Soc. Lond., B, Biol. Sci. 376 (2021).Phillips, H. R. P. et al. Global distribution of earthworm diversity. Science 366, 480–485 (2019).ADS 
CAS 
Article 

Google Scholar 
Geitner, C. et al. Soil and Land Use in the Alps—Challenges and Examples of Soil-Survey and Soil-Data Use to Support Sustainable Development. In Soil Mapping and Process Modelling for Sustainable Land Use Management 221–292 (Elsevier, 2017).FAO. Understanding Mountain Soils: A contribution from mountain areas to the International Year of Soils 2015. (2015).Praeg, N. et al. The role of land management and elevation in shaping soil microbial communities: Insights from the Central European Alps. Soil Biol. Biochem. 150, 107951 (2020).CAS 
Article 

Google Scholar 
Somme, L. Adaptations of terrestrial arthropods to the alpine environment. Biol. Rev. 64, 367–407 (1989).Article 

Google Scholar 
Meyer, E. & Thaler, K. Animal Diversity at High Altitudes in the Austrian Central Alps. In Arctic and Alpine Biodiversity (eds. Chapin, F. S. & Körner, C.) (Springer Berlin Heidelberg, 1995).Guerra, C. A. et al. Blind spots in global soil biodiversity and ecosystem function research. Nat. Commun. 11, 1–13 (2020).Article 

Google Scholar 
Kempson, D., Lloyd, M. & Ghelardi, R. A new extractor for woodland litter. Pedobiologia 3, 1–30 (1963).
Google Scholar 
Tasser, E., Ruffini, F. V. & Tappeiner, U. An integrative approach for analysing landscape dynamics in diverse cultivated and natural mountain areas. Landsc. Ecol. 24, 611–628 (2009).Article 

Google Scholar 
European Environment Agency, (EAA). European Union, Copernicus Land Monitoring Service. https://land.copernicus.eu/pan-european/corine-land-cover/clc2018 (2018).IUSS Working Group WRB. World Reference Base for Soil Resources 2014, update 2015 International soil classification system for naming soils and creating legends for soil maps. World Soil Resources Reports No. 106 (2015).Christian, E. & Zicsi, A. A synoptic key to the earthworms of Austria (Oligochaeta: Lumbricidae). Die Bodenkultur 50, 121–131 (1999).Czusdi, C. & Zicsi, A. Earthworms of Hungary (Annelida, Oligochaeta, Lumbricidae). (Hungarian Natural History Museum, 2003).Schaefer, M. Brohmer – Fauna von Deutschland: Ein Bestimmungsbuch unserer heimischen Tierwelt. (Quelle & Meyer Verlag, 2009).Klausnitzer, B. Exkursionsfauna von Deutschland. (Springer Berlin Heidelberg, 2011).Ellis, D. Taxonomic sufficiency in pollution assessment. Mar. Pollut. Bull. 16, 459 (1985).Article 

Google Scholar 
Caruso, T. & Migliorini, M. Micro-arthropod communities under human disturbance: is taxonomic aggregation a valuable tool for detecting multivariate change? Evidence from Mediterranean soil oribatid coenoses. Acta Oecol. 30, 46–53 (2006).ADS 
Article 

Google Scholar 
Steinwandter, M., Schlick-Steiner, B. C., Seeber, G. U. H., Steiner, F. M. & Seeber, J. Effects of Alpine land-use changes: Soil macrofauna community revisited. Ecol. Evol. 7, 5389–5399 (2017).Article 

Google Scholar 
de Oliveira, S. S. et al. Higher taxa are sufficient to represent biodiversity patterns. Ecol. Indic. 111, 105994 (2020).Article 

Google Scholar 
Lavelle, P. et al. Soil macroinvertebrate communities: A world‐wide assessment. Glob. Ecol. Biogeogr. 31, 1261–1276 (2022).Article 

Google Scholar 
Pik, A. J., Oliver, I. & Beattie, A. J. Taxonomic sufficiency in ecological studies of terrestrial invertebrates. Austral. Ecol. 24, 555–562 (1999).Article 

Google Scholar 
Parisi, V., Menta, C., Gardi, C., Jacomini, C. & Mozzanica, E. Microarthropod communities as a tool to assess soil quality and biodiversity: A new approach in Italy. Agric. Ecosyst. Environ. 105, 323–333 (2005).Article 

Google Scholar 
Ruiz, N. et al. IBQS: A synthetic index of soil quality based on soil macro-invertebrate communities. Soil Biol. Biochem. 43, 2032–2045 (2011).
Google Scholar 
Seeber, J. et al. A 30-years collection of soil macro-invertebrate abundance data from the European Alps. PANGAEA https://doi.org/10.1594/PANGAEA.944405 (2022).de Jong, Y. et al. Fauna Europaea – All European animal species on the web. Biodivers. Data J. 2 (2014).Shannon, C. E. A Mathematical Theory of Communication. Bell Syst. Tech. J. 27, 379–423 (1948).MathSciNet 
Article 

Google Scholar 
Steinberger, K.-H. & Meyer, E. The Spider Fauna of the Nature Reserve Rheindelta (Vorarlberg, Western Austria) (Arachnida: Araneae). Ber. nat.-med. Verein Innsbruck 82, 195–215 (1995).
Google Scholar 
Kopeszki, H. & Meyer, E. Species Composition and Abundance of Collembola in Forest Soils in the Provinces of Bozen and Trient (Italy). Ber. nat.-med. Verein Innsbruck 83, 221–237 (1996).
Google Scholar 
Geitner, C., Mätzler, A., Bou-Vinals, A., Meyer, E. & Tusch, M. Soil characteristics and colonization by earthworms (Lumbricidae) on pastures and hay meadows in the Brixenbach Valley (Kitzbühel Alps, Tyrol). Die Bodenkultur 65, 39–51 (2014).
Google Scholar 
Meyer, E. & Steinberger, K.-H. Über die Bodenfauna in Wäldern Vorarlbergs (Österreich) Bestand und Auswirkungen von Gesteinsmehlapplikationen. Verhandlungen der Gesellschaft für Ökologie 23, 149–164 (1994).
Google Scholar 
Peham, T. & Meyer, E. Kommentierte Artenlisten ausgewählter Bodentiergruppen aus der Erhebung des SoilDiv-Projektes in Südtirol. Gredleriana 14, 227–262 (2014).
Google Scholar 
Steinwandter, M., Rief, A., Scheu, S., Traugott, M. & Seeber, J. Structural and functional characteristics of high alpine soil macro-invertebrate communities. Eur. J. Soil Biol. 86, 72–80 (2018).Article 

Google Scholar 
Steinwandter, M. et al. Does green manuring positively affect the soil macro-invertebrates in vineyards? A case study from Kaltern/Caldaro (South Tyrol, Italy). Gredleriana 18, 17–26 (2018).
Google Scholar 
Steinwandter, M. et al. Raw data from: Does green manuring positively affect the soil macro-invertebrates in vineyards? A case study from Kaltern/ Caldaro (South Tyrol, Italy). PANGAEA https://doi.org/10.1594/PANGAEA.900632 (2019).Damisch, K., Steinwandter, M., Tappeiner, U. & Seeber, J. Soil Macroinvertebrate Distribution along a Subalpine Land Use Transect. Mt. Res. Dev. 40, R1–R10 (2020).Article 

Google Scholar 
Damisch, K., Steinwandter, M., Tappeiner, U. & Seeber, J. Abundance data from soil macro-invertebrates along a subalpine land-use transect. PANGAEA https://doi.org/10.1594/PANGAEA.918958 (2020).Schneider, E., Steinwandter, M. & Seeber, J. A comparison of Alpine soil macro-invertebrate communities from European larch and Swiss pine forests in the LTSER area “Val Mazia/Matschertal”, South Tyrol. Gredleriana 19, 217–228 (2019).
Google Scholar 
Schneider, E., Steinwandter, M. & Seeber, J. Raw data from: A comparison of Alpine soil macro-invertebrate communities from European larch and Swiss pine forests in the LTSER area “Val Mazia/Matschertal”, South Tyrol. PANGAEA https://doi.org/10.1594/PANGAEA.910666 (2020).Seeber, J. et al. Soil invertebrate abundance, diversity, and community composition across steep high elevation snowmelt gradients in the European Alps. Arct. Antarct. Alp. Res. 53, 288–299 (2021).Article 

Google Scholar 
Seeber, J. et al. Abundance data from soil macro- and mesofauna of alpine snowbeds in the European Alps (summer 2017). PANGAEA https://doi.org/10.1594/PANGAEA.935737 (2021).Guariento, E. et al. Management intensification of hay meadows and fruit orchards alters soil macro-invertebrate communities differently. Agronomy 10, 767 (2020).CAS 
Article 

Google Scholar  More

Similar non-native and invasive flora, such as the fever tree (pictured) are found in regions previously occupied by the same European empire.Credit: Alamy

In 1860, a British expedition raided the highland forests of South America, looking for a hot commodity: Cinchona seeds. The bark of these ‘fever’ trees produces the anti-malarial compound quinine, and the British Empire sought a stable source of the drug for its soldiers and civil service in India. After cultivation in the United Kingdom, young Cinchona trees were planted across southern India and what is now Sri Lanka.The British quinine scheme failed — instead, a species introduced to Java, now part of Indonesia, by the Dutch Empire later dominated the global market — but Cinchona trees are still common in parts of India.Such botanical legacies of imperial rule are common, finds a study published on 17 October in Nature Ecology & Evolution1. Regions that were once occupied by the same European colonial power — such as India and Sri Lanka — tend to have similar species of non-native and invasive plants. The longer the regions were occupied, the more their populations of invasive species resemble each other, the research found.Alien floraThe link between European colonialism and invasive species is intuitive, and has been noted by other researchers, says Bernd Lenzner, a macro-ecologist at the University of Vienna who led the study. To test the association, his team turned to the Global Naturalized Alien Flora database, which maps the distribution of nearly 14,000 invasive plant species.
The imperial roots of climate science
Across more than 1,100 regions, including 404 islands, the researchers found that regions once occupied by the British Empire had more similarities in their invasive flora than did ‘artificial’ empires that the team assembled from random regions. This was also the case for regions once part of the Dutch Empire (former Spanish and Portuguese colonies had alien-plant compositions similar to those of the artificial empires).Climate and geography play an important part in explaining the overlap in the diversity of invasive species, modelling by Lenzner’s team found, but so does the length of time regions were occupied by an imperial power. Regions that were central to trade, such as southern India for the British Empire and Indonesia for the Dutch Empire, formed clusters with considerable overlap in invasive-plant composition.The analysis did not look at when individual plant species were introduced or why. But anecdotally, many of the plants that were commonly taken to former empires were once of economic value and their populations were probably established on purpose, says Lenzner.Global trade impactsThe study’s conclusions might be “super obvious”, but they have important implications for conservation, says Nussaïbah Raja, a palaeontologist at Friedrich-Alexander University of Erlangen–Nürnberg in Erlangen, Germany. “We should be taking this history into consideration when we think about management of species.” Appreciating the history of introduced plants — as well as their place in today’s ecosystems — could help conservationists to handle future changes in biodiversity, such as those driven by climate change, Raja adds.Global trade is beginning to overwrite the colonial legacy of introduced plants. For example, the analysis showed similarities between invasive plant populations in Fujian, China, and some parts of Australia. Although both places were once connected by the British Empire, more recent global trade might also be partly responsible for the overlap.“We are still seeing these imprints of the colonial-empire legacies from centuries ago,” Lenzner says. “So what we’re doing and the species we’re redistributing today will be visible far into the future.” More

This study aimed to compare the extent of contrafreeloading in kea to that in Grey parrots, given that the two species exhibit very different levels of play: specifically, kea exhibit complex and frequent play29,30,35,36, whereas Greys exhibit considerably less play than several parrot species29. We found that, at the group level, although the overall amounts of kea classic contrafreeloading were nonsignificant, as a percentage of behaviour, kea generally contrafreeloaded more than Grey parrots in Experiment 1, whereas the opposite was true for Experiment 2. We compare the various behaviour patterns in detail, and propose explanations for our results below.The most interesting comparisons for Smith et al.’s hypothesis are the results from classic contrafreeloading. In Experiment 1, kea performed this behaviour at non-negligible levels, given the supposed rarity of the behaviour5 (two birds at 50%; the others varying between 39 and 47%). In contrast, although one Grey did classically contrafreeload at a statistically significant level, the other three were at ≤ 36%. These data suggest that the kea may have found the task more engaging than did the Greys. However, given that only two kea chose to pop the lid of an empty cup in control trials significantly above chance, whereas three of the four Greys did so significantly above chance and one at chance, we doubt that the kea found the task inherently rewarding. We note that this comparison between both species must be interpreted cautiously due to differences in methodology: For the Greys, the control trials were performed at the end of the study, by which point they may have learnt to associate lid-popping with reward. However, the data from experimental trials in Smith et al.13 are such that their birds would have been primed in the opposite direction: For example, three of those four birds rarely chose the empty lidded cup when free food was available, nor did they classically or super contrafreeload to any significant extent13; an association-driven explanation is therefore unlikely. In contrast, the kea experienced this control condition at the start of the experiment, allowing them 20 trials to become acquainted with the affordances of both options that would be available throughout the study (lid-popping versus not lid-popping). This opportunity was important for kea, as this species has been previously shown to learn about object properties through extensive object manipulation37. That kea popped lids at or above chance in these first 20 control trials suggested two possibilities: (1) After these 20 trials, the task may have been familiar enough to no longer be of much interest (i.e., no longer novel and worthy of consideration) by the time rewarded trials began (recall nonsignificant downward trends for Harley Quinn and Blofeld). (2) They acquired some interest in popping the lids. This latter case seems more likely, as the lid-popping task still likely provided some added value. Kea engaged in non-negligible levels of classic contrafreeloading, such that the chance to pop a lid and eat could be considered more interesting than simply eating an identical but freely available reward. Furthermore, three kea chose a lidded, empty cup over a free, least-preferred reward at least half the time, again suggesting that the activity held some appeal of its own.In Experiment 2 (which corresponds to classic contrafreeloading), all kea preferred freeloading for the walnut without a shell; two Greys, in contrast, nut contrafreeloaded at a statistically significant extent. This variability in behaviour at both the individual and species levels reveals the significance of a task’s proximate and potentially ultimate values in parrots’ choice to contrafreeload. Interestingly, although species like kea are hypothesized to prefer food items requiring high manipulation38,39, nut-cracking—chosen as an activity to provide direct comparison with the Greys13—is not prevalent in kea diet40, and that activity thus may not have been appropriate as an ethologically relevant one for kea. Greys, in contrast, are known to crack nuts in nature41. Future research could use a more ecologically relevant task for the kea, such as working to access food via digging or scraping32.As with Smith et al.’s Greys13, kea in Experiment 1 performed calculated contrafreeloading to a statistically significant extent. All kea did so on over 83% of trials; for the Greys, three birds were close to 90% but one was at only 67%. Kea consistently selected their preferred food out of the two options provided, suggesting that the lid-popping action did not deter kea from selecting their preferred reward. In related trials, where the lid-status of food paired with an empty cup varied, kea, like some Greys13, preferred lidded food over an empty lidless cup, again showing that lid-popping for food was an acceptable task.When examining situations in which food was discarded after contrafreeloading, we found that this choice in Experiment 1 was most common for Bruce. Notably, Bruce lacks a top mandible, making many of the manipulative behaviours more difficult to execute42. Bruce demonstrated consistent food preferences throughout the experiment, however, indicating that the reason some foods were discarded was, indeed, because they were too difficult for him to manipulate. In Experiment 2, Harley Quinn was the most likely to discard the nut, and did so exclusively in trials in which she chose the walnut without the shell (freeloaded). In these occasions, Harley Quinn was observed choosing the nut by tapping on it or the cup.Like the Greys, the kea failed to super contrafreeload to a statistically significant extent. Furthermore, contrafreeloading trials in which a lid was popped but the food underneath was not consumed occurred most often with the least-preferred food. Given kea’s performance on control trials, the super contrafreeloading results are not surprising. Interestingly, when lid-status of food paired with an empty cup varied, some Greys very rarely—and depending on food desirability—preferred to pop the empty cup’s lid rather than consume the free food; as noted earlier, three of eight kea did so on at least half the trials when the food in the lidless cup was their least preferred option (sultanas). Both kea and Greys thus likely placed the appeal of the task along some “value scale” along with that of the available food rewards, the combination influencing their behaviour when the two variables were presented in various permutations. Notably, even in control trials, where no food was involved, no bird of either species found the task aversive, engaging in the behaviour at least 50% of the time. Future research could investigate how a different, more rewarding task would influence this balance and thus contrafreeloading for both species.One possible alternative explanation for kea’s higher rates of contrafreeloading relative to those of Greys could be their natural tendency to probe and manipulate objects, thus causing them to pry off cup lids rather than manipulate lidless (open) cups. Were this action exploratory in nature, we would have observed significant decreases in behaviour as the experiment progressed, but note that we found no significant changes in any bird. Were they consistently drawn to lids and this behaviour were hard-wired, then we should have observed lid-popping appear significantly above chance across all three types of contrafreeloading. However, as discussed previously, kea did not significantly contrafreeload in the classic condition and actively freeloaded in super contrafreeloading conditions, suggesting that they were not simply interacting with lidded cups preferentially, but rather attending to the contents in the two cups and avoiding the additional manipulation of the lid when it led to a less (or, more often than not, equally) preferred food reward.Another potential explanation for the differences observed between kea and Greys might be found in the theoretical overlap between contrafreeloading and play, and how individuals might view the contrafreeloading action as a type of play. As a seemingly nonfunctional, intrinsically motivating behaviour occurring in low-stress environments, incurring a positive mood, varying between conspecifics, and often incomplete and/or repeated14,15, play shares many proximate-level attributes with contrafreeloading13. Our results demonstrate that kea subjects inhabiting a low-stress, captive environment repeatedly chose to engage in classic contrafreeloading to a non-negligible extent and calculated contrafreeloading to a significant extent, varied in their behaviour between individuals, and at times, left the task incomplete (e.g., left food uneaten). Furthermore, evidence for intrinsic motivation to perform a given task is suggested by the kea’s overall differential behaviour between the two experiments, as well as inter-individual differences.Importantly, this study serves only as a first step into determining whether play manifests as a form of contrafreeloading, but cannot ascertain that this is the only possible explanation for the presence or degree of contrafreeloading in the two species. Several alternative explanatory theories regarding the occurrence of contrafreeloading are enumerated in the discussion of Smith et al. (e.g., work ethic; information gathering; relief from boredom)13, and various other potential explanations (beyond playfulness) may reside at the species-level. Grey parrots (Psittacidae) and kea (Strigopidae) are separated by 50–80 million years of evolution43 and differ in their neurobiology (i.e., the size of the shell region related to vocal and possible cognitive abilities44). Differing ecological evolutionary pressures are also likely relevant: an island-based habitat39, a lack of natural predators30,45, and generalist diets40,46,47 are thought to have shaped the playfulness and cognitive abilities of kea30,40,46,47. Greys, in contrast, evolved predominantly on a continent (i.e., although they can be found on islands such as Principe, the Congo Grey is endemic to central Africa48,49), are subject to considerable predation48,50,51,52, and have a relatively less generalist diet (diverse but almost exclusively vegetarian and in which nuts play a significant role; see review in50). Such disparate evolutionary trajectories may offer other potential explanations for the differences in contrafreeloading observed between the two species, and future research could examine differences at genetic and/or neurological levels.The varying rates of contrafreeloading observed between the species could have also been influenced by other factors. For example, although both parrot groups studied here inhabit enriched environments, are habituated to participating in experimental trials, and have access to food ad libitum, their habitats are markedly different. Notably, the Grey subjects live in “man-made” settings (i.e., Griffin and Athena in a lab; Pepper, Franco, and Lucci in private homes), whereas the kea inhabit a naturalistic zoo enclosure. Physical enrichment, although somewhat different in kind, is unlikely to have differed in quantity, as all birds are provided routine naturalistic foraging, and Lucci lives in a free-flight aviary. More likely is the difference in sociality: Relatively more subjects reside together in the kea group (15) compared to the Greys (two groups of two Greys and one Grey living with two birds of differing species), and thus variables such as social stimulation and flock-based foraging techniques could have contributed to the expression of contrafreeloading (note that subadult male kea are known to obtain food through kleptoparasitism32). In order to elucidate the role of habitat on contrafreeloading, future studies could examine the behaviour of species residing in more comparable captive conditions.Future work should aim not only to apply these same methodologies to a broader range of parrot species, but also objectively quantify frequency and complexity of play across a wide range of parrots to allow a direct correlation between play and contrafreeloading over phylogeny in the parrot order. The apparent link between play behaviour and encephalisation in parrots53 offers another possible avenue for cross-species comparisons on contrafreeloading. Future research could also employ cognitive bias tests to quantify the mood of birds before and following contrafreeloading54, directly manipulate subjects’ participation in play behaviours or other control behaviours and observe whether engaging in play can increase contrafreeloading rates at the individual level, or perform behavioural coding of playfulness and/or arousal before and after contrafreeloading. Future research could incorporate more ecologically relevant contrafreeloading tasks to examine this behaviour at both the individual and species level, and approach the phenomenon by using both genetic and neuroscience techniques.In sum, contrafreeloading is, by its very nature, an enigma whose study presents many difficulties. It varies across the diverse contexts within which it is studied, and given that it is rarely exhibited to a statistically significant extent, analyses that require comparing nonsignificant behaviour patterns across individuals and/or species is a challenging undertaking. Many explanations have been proposed, but contrafreeloading is still poorly understood, and its correlation with play is likely only one of several logical rationales. Nevertheless, our findings suggest that interest in play should not be discounted as a contributing factor. More