Kinship,inbreeding and fine‐scale spatial structure influence gut microbiota in a hindgut‐fermenting tortoise

Herbivorous vertebrates rely on complex communities of mutualistic gut bacteria to facilitate the digestion of celluloses and hemicelluloses. Gut microbes are often convergent based on diet and gut morphology across a phylogenetically diverse group of mammals. However, little is known about microbial communities of herbivorous hindgut‐fermenting reptiles. Here, we investigate how factors at the individual level might constrain the composition of gut microbes in an obligate herbivorous reptile. Using multiplexed 16S rRNA gene sequencing, we characterized the faecal microbial community of a population of gopher tortoises (Gopherus polyphemus) and examined how age, genetic diversity, spatial structure and kinship influence differences among individuals. We recovered phylotypes associated with known cellulolytic function, including candidate phylum Termite Group 3, suggesting their importance for gopher tortoise digestion. Although host genetic structure did not explain variation in microbial composition and community structure, we found that fine‐scale spatial structure, inbreeding, degree of relatedness and possibly ontogeny shaped patterns of diversity in faecal microbiomes of gopher tortoises. Our findings corroborate widespread convergence of faecal‐associated microbes based on gut morphology and diet and demonstrate the role of spatial and demographic structure in driving differentiation of gut microbiota in natural populations.     

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.     

Isolation of tannase‐producing microbiota from the gastrointestinal tracts of some freshwater fish

Tannins are the most abundant among the plant‐derived antinutrients that bind readily with protein and other macromolecules to form indigestible complexes, thereby reducing the nutritional value of the plant feedstuffs. Presence of tannase‐producing gut microbiota in herbivorous animals has been suggested to overcome the antinutritional effects of tannins. However, this topic has been less investigated in herbivorous/omnivorous fish species. The present study was undertaken to evaluate the presence of tannase‐producing autochthonous microbiota in the gastrointestinal (GI) tracts of some culturable freshwater teleosts and to identify most promising tannase‐producing strains by molecular methods. Isolation and enumeration of tannase‐producing autochthonous microbiota have been carried out in the gut of ten culturable freshwater teleosts, namely catla (Catla catla), silver carp (Hypophthalmichthys molitrix), rohu (Labeo rohita), grass carp (Ctenopharyngodon idella), mrigal (Cirrhinus mrigala), common carp (Cyprinus carpio), bata (Labeo bata), kalbasu (Labeo calbasu), tilapia (Oreochromis mossambicus), and Nile tilapia (Oreochromis niloticus). Culturable heterotrophic and tannase‐producing microbial populations evaluated on tryptone soya agar and selective tannic acid agar media, respectively, revealed the maximum in the hindguts of all fish species studied. Out of 72 tannase‐producing colonies, 18 randomly selected isolates were maintained as pure cultures and evaluated quantitatively for tannase production. Among these, four most promising tannase producers were identified by 16S/26S rDNA sequencing following nucleotide blast and deposited in the National Centre for Biotechnology Information (NCBI) GenBank. The strain LR01 isolated from rohu was a bacterium, Enterobacter asburae (GenBank Accession No. GU939631 ). However, the strains CM02, OM01 and LR03 isolated from mrigal, tilapia and rohu were yeasts and identified as Pichia kudriavzevii (GenBank Accession No. GU939629 ), Candida tropicalis (GenBank Accession No. GU911469 ) and Candida parapsilosis (GenBank Accession No. GU939630 ), respectively. To the authors' knowledge, the present study is the first to report tannase‐producing autochthonous microbiota in the gut of freshwater teleosts. Tannin‐degrading microbiota detected in the present study may endow the fish with some ecological advantages by enabling them to overcome the anti‐nutritional effects of plant tannins.     

Comammox Nitrospira are abundant ammonia oxidizers in diverse groundwater‐fed rapid sand filter communities

The recent discovery of completely nitrifying Nitrospira demands a re‐examination of nitrifying environments to evaluate their contribution to nitrogen cycling. To approach this challenge, tools are needed to detect and quantify comammox Nitrospira. We present primers for the simultaneous quantification and diversity assessement of both comammox Nitrospira clades. The primers cover a wide range of comammox diversity, spanning all available high quality sequences. We applied these primers to 12 groundwater‐fed rapid sand filters, and found comammox Nitrospira to be abundant in all filters. Clade B comammox comprise the majority (～75%) of comammox abundance in all filters. Nitrosomonadaceae were present in all filters, although at low abundance (mean = 1.8%). Ordination suggests that temperature impacts the structure of nitrifying communities, and in particular that increasing temperature favours Nitrospira. The nitrogen content of the filter material, sulfate concentration and surface ammonium loading rates shape the structure of the comammox guild in the filters. This work provides an assay for simultaneous detection and diversity assessment of clades A and B comammox Nitrospira, expands our current knowledge of comammox Nitrospira diversity and demonstrates a key role for comammox Nitrospira in nitrification in groundwater‐fed biofilters.     

Sediment microbial communities in Great Boiling Spring are controlled by temperature and distinct from water communities

Great Boiling Spring is a large, circumneutral, geothermal spring in the US Great Basin. Twelve samples were collected from water and four different sediment sites on four different dates. Microbial community composition and diversity were assessed by PCR amplification of a portion of the small subunit rRNA gene using a universal primer set followed by pyrosequencing of the V8 region. Analysis of 164 178 quality-filtered pyrotags clearly distinguished sediment and water microbial communities. Water communities were extremely uneven and dominated by the bacterium Thermocrinis. Sediment microbial communities grouped according to temperature and sampling location, with a strong, negative, linear relationship between temperature and richness at all taxonomic levels. Two sediment locations, Site A (87–80 °C) and Site B (79 °C), were predominantly composed of single phylotypes of the bacterial lineage GAL35 ($\widehat{p}$=36.1%), Aeropyrum ($\widehat{p}$=16.6%), the archaeal lineage pSL4 ($\widehat{p}$=15.9%), the archaeal lineage NAG1 ($\widehat{p}$=10.6%) and Thermocrinis ($\widehat{p}$=7.6%). The ammonia-oxidizing archaeon ‘Candidatus Nitrosocaldus'' was relatively abundant in all sediment samples <82 °C ($\widehat{p}$=9.51%), delineating the upper temperature limit for chemolithotrophic ammonia oxidation in this spring. This study underscores the distinctness of water and sediment communities in GBS and the importance of temperature in driving microbial diversity, composition and, ultimately, the functioning of biogeochemical cycles.     

Peatland forests are the least diverse tree communities documented in Amazonia,but contribute to high regional beta‐diversity

Western Amazonia is known to harbour some of Earth's most diverse forests, but previous floristic analyses have excluded peatland forests which are extensive in northern Peru and are among the most environmentally extreme ecosystems in the lowland tropics. Understanding patterns of tree species diversity in these ecosystems is important both for quantifying beta‐diversity in this region, and for understanding determinants of diversity more generally in tropical forests. Here we explore patterns of tree diversity and composition in two peatland forest types – palm swamps and peatland pole forests – using 26 forest plots distributed over a large area of northern Peru. We place our results in a regional context by making comparisons with three other major forest types: terra firme forests (29 plots), white‐sand forests (23 plots) and seasonally‐flooded forests (11 plots). Peatland forests had extremely low (within‐plot) alpha‐diversity compared with the other forest types that were sampled. In particular, peatland pole forests had the lowest levels of tree diversity yet recorded in Amazonia (20 species per 500 stems, Fisher's alpha 4.57). However, peatland pole forests and palm swamps were compositionally different from each other as well as from other forest types in the region. Few species appeared to be peatland endemics. Instead, peatland forests were largely characterised by a distinctive combination of generalist species and species previously thought to be specialists of other habitats, especially white‐sand forests. We suggest that the transient nature and extreme environmental conditions of Amazonian peatland ecosystems have shaped their current patterns of tree composition and diversity. Despite their low alpha‐diversity, the unique combination of species found in tree communities in Amazonian peatlands augment regional beta‐diversity. This contribution, alongside their extremely high carbon storage capacity and lack of protection at national level, strengthens their status as a conservation priority.     

Genomes of rumen bacteria encode atypical pathways for fermenting hexoses to short‐chain fatty acids

Bacteria have been thought to follow only a few well‐recognized biochemical pathways when fermenting glucose or other hexoses. These pathways have been chiseled in the stone of textbooks for decades, with most sources rendering them as they appear in the classic 1986 text by Gottschalk. Still, it is unclear how broadly these pathways apply, given that they were established and delineated biochemically with only a few model organisms. Here, we show that well‐recognized pathways often cannot explain fermentation products formed by bacteria. In the most extensive analysis of its kind, we reconstructed pathways for glucose fermentation from genomes of 48 species and subspecies of bacteria from one environment (the rumen). In total, 44% of these bacteria had atypical pathways, including several that are completely unprecedented for bacteria or any organism. In detail, 8% of bacteria had an atypical pathway for acetate formation; 21% of bacteria had an atypical pathway for propionate or succinate formation; 6% of bacteria had an atypical pathway for butyrate formation and 33% of bacteria had an atypical or incomplete Embden–Meyerhof–Parnas pathway. This study shows that reconstruction of metabolic pathways – a common goal of omics studies – could be incorrect if well‐recognized pathways are used for reference. Furthermore, it calls for renewed efforts to delineate fermentation pathways biochemically.     

Bacterial translocation from the gastrointestinal tracts of thymectomized mice

The incidence of translocation of viable indigenous bacteria from the gastrointestinal tract to the mesenteric lymph node, spleen, liver, and kidney was compared in neonatally thymectomized mice and sham-thymectomized specific pathogen-free mice. The immunologic responses of the thymectomized mice to sheep erythrocytes were decreased compared to the responses of sham-thymectomized mice. Strictly anaerobic bacteria were isolated from only 1.8% of the organs from thymectomized mice and from none of the organs of shamthymectomized mice. Aerobic or facultatively anaerobic bacteria were cultured from 27.4% of the organs of thymectomized mice. Of the thymectomized mice, 70.7% contained viable aerobic or facultatively anaerobic bacteria in one or more of their organs tested, compared with only 10% of the sham-thymectomized mice.Escherichia coli was the predominant bacterial species isolated from these organs, althoughStaphylococcus aureus, Streptococcus, andCorynebacterium also were present.Bacteroides were the only strictly anaerobic bacteria cultured. Neonatal thymectomy promotes the translocation of certain indigenous bacteria from the gastrointestinal tract to the mesenteric lymph node, spleen, liver, and kidney.     

Evolutionarily distinct amphibians are disproportionately lost from human‐modified ecosystems

Humans continue to alter terrestrial ecosystems, but our understanding of how biodiversity responds is still limited. Anthropogenic habitat conversion has been associated with the loss of evolutionarily distinct bird species at local scales, but whether this evolutionary pattern holds across other clades is unknown. We collate a global dataset on amphibian assemblages in intact forests and nearby human‐modified sites to assess whether evolutionary history influences susceptibility to land conversion. We found that evolutionarily distinct amphibian species are disproportionately lost when forested habitats are converted to alternative land‐uses. We tested the hypothesis that grassland‐associated amphibian lineages have both higher diversification and are pre‐adapted to human landscapes, but found only weak evidence supporting this. The loss of evolutionarily distinct amphibians with land conversion suggests that preserving remnant forests will be vital if we aim to preserve the amphibian tree of life in the face of mounting anthropogenic pressures.     

28S rDNA haplotypes of males are distinct from those of androgenetic hermaphrodites in the clam Corbicula leana

The clam Corbicula leana exists in two forms, hermaphrodites and males. Our previous study on mitochondrial DNA suggested that the male nuclear DNA might have derived from hermaphrodite C. leana relatively recently. To clarify the origin of males in the clam, sequences of the nuclear 28S rDNA divergent domain (which is 441-444 bp long) in androgenetic hermaphrodites and males and dioecious (bisexual) species were analyzed. Unexpectedly, the nuclear 28S rDNA haplotypes of males and hermaphrodites were distinct. Haplotype network analysis indicated that males and hermaphrodites are reproductively isolated from each other without sharing the same nuclear haplotype. These results support a hypothesis that the egg nuclear genome of androgenetic hermaphrodites is replaced by the male sperm genome, and only males develop after fertilization by a male spermatozoon.     

Emerging from the ice-fungal communities are diverse and dynamic in earliest soil developmental stages of a receding glacier

We used amplicon sequencing and isolation of fungi from in-growth mesh bags to identify active fungi in three earliest stages of soil development (SSD) at a glacier forefield (0–3, 9–14, 18–25 years after retreat of glacial ice). Soil organic matter and nutrient concentrations were extremely low, but the fungal diversity was high [220 operational taxonomic units (OTUs)/138 cultivated OTUs]. A clear successional trend was observed along SSDs, and species richness increased with time. Distinct changes in fungal community composition occurred with the advent of vascular plants. Fungal communities of recently deglaciated soil are most distinctive and rather similar to communities typical for cryoconite or ice. This indicates melting water as an important inoculum for native soil. Moreover, distinct seasonal differences were detected in fungal communities. Some fungal taxa, especially of the class Microbotryomycetes, showed a clear preference for winter and early SSD. Our results provide insight into new facets regarding the ecology of fungal taxa, for example, by showing that many fungal taxa might have an alternative, saprobial lifestyle in snow-covered, as supposed for a few biotrophic plant pathogens of class Pucciniomycetes. The isolated fungi include a high proportion of unknown species, which can be formally described and used for experimental approaches.     

Phylogenetic analysis reveals positive correlations between adaptations to diverse hosts in a group of pathogen‐like herbivores

A jack of all trades can be master of none—this intuitive idea underlies most theoretical models of host‐use evolution in plant‐feeding insects, yet empirical support for trade‐offs in performance on distinct host plants is weak. Trade‐offs may influence the long‐term evolution of host use while being difficult to detect in extant populations, but host‐use evolution may also be driven by adaptations for generalism. Here we used host‐use data from insect collection records to parameterize a phylogenetic model of host‐use evolution in armored scale insects, a large family of plant‐feeding insects with a simple, pathogen‐like life history. We found that a model incorporating positive correlations between evolutionary changes in host performance best fit the observed patterns of diaspidid presence and absence on nearly all focal host taxa, suggesting that adaptations to particular hosts also enhance performance on other hosts. In contrast to the widely invoked trade‐off model, we advocate a “toolbox” model of host‐use evolution in which armored scale insects accumulate a set of independent genetic tools, each of which is under selection for a single function but may be useful on multiple hosts.     

Statistical challenges of evaluating diversity patterns across environmental gradients in mega‐diverse communities

Ibanez et al. (Journal of Vegetation Science, this issue) applied sample size‐ and coverage‐based rarefaction to analyse the elevational richness pattern in New Caledonian tree communities. We comment on the statistical assumptions behind rarefaction/extrapolation and suggest pooling small plot data to effectively assess/detect the diversity pattern. Broadening the analysis to include abundance‐sensitive diversity measures and phylogenetic information can provide important additional insights.     

Bee communities (Hymenoptera: Apoidea) of small Iowa hill prairies are as diverse and rich as those of large prairie preserves

Natural small, xeric hill prairies in forested landscapes throughout the Midwest United States often contain high diversity and unique species of some organisms because of their unusual landscape context and microclimate. We measured the diversity, richness, and abundance of the bee communities of five hill prairies located in northeastern Iowa and we compared these to values for large prairie preserves in northwestern Iowa, using a Monte Carlo resampling approach to standardize sampling effort between habitat types. We also measured the diversity and richness of the flowering forb communities at the hill prairies and we quantified percentage of the landscape at a 1 km radius in different landscape elements. Bee diversity at the five hill prairies spanned the range of diversity values for large prairies preserves, so although the hill prairies are small (<5 ha), their bee communities are not uniformly depauparate compared to larger western prairie preserves. Bee diversity was significantly related to flowering forb diversity, and may have been influenced by landscape features within 1 km—particularly the percentage of agricultural row crops and open water, which may negatively affect bee diversity at the sites. Iowa’s hill prairie bee communities were largely composed of widespread eastern species, although about 10% of the bee species have more northern or western ranges and appear to be taking advantage of the region’s unique habitat features. Given the dependence of the bee communities on the plant diversity of the sites, management of the plant community to maintain its diversity will also likely benefit bee diversity.     

Effects of 1-oleoyl-2-acetyl glycerol are distinct from those of phorbol ester in rat pancreatic acini

1-oleoyl-2-acetyl glycerol (OAG), a potent activator of protein kinase C, inhibited the binding of 125I-labelled epidermal growth factor (EGF) in isolated rat pancreatic acini. Unlike cholecystokinin-octapeptide (CCK8) and the C-kinase activator 12-O-tetradecanoyl phorbol-13-acetate (TPA), two inhibitors of 125I-EGF endocytosis in the pancreas, OAG had no effect on the distribution of bound ligand between the cell surface and intracellular compartments. Unlike TPA, OAG failed to potentiate the inhibitory effects of the calcium ionophore A23187 on 125I-EGF cell-associated radioactivity and had no effect on either basal or carbachol-stimulated amylase release in acini. These data suggest that the actions of the synthetic diacyl-glycerol OAG are not fully equivalent with the action of other known activators of protein kinase C in the pancreatic acinar cell.     

BFA‐induced compartments from the Golgi apparatus and trans‐Golgi network/early endosome are distinct in plant cells

Brefeldin A (BFA) is a useful tool for studying protein trafficking and identifying organelles in the plant secretory and endocytic pathways. At low concentrations (5–10 μg ml^?1), BFA caused both the Golgi apparatus and trans‐Golgi network (TGN), an early endosome (EE) equivalent in plant cells, to form visible aggregates in transgenic tobacco BY‐2 cells. Here we show that these BFA‐induced aggregates from the Golgi apparatus and TGN are morphologically and functionally distinct in plant cells. Confocal immunofluorescent and immunogold electron microscope (EM) studies demonstrated that BFA‐induced Golgi‐ and TGN‐derived aggregates are physically distinct from each other. In addition, the internalized endosomal marker FM4‐64 co‐localized with the TGN‐derived aggregates but not with the Golgi aggregates. In the presence of the endocytosis inhibitor tyrphostin A23, which acts in a dose‐ and time‐dependent manner, SCAMP1 (secretory carrier membrane protein 1) and FM4‐64 are mostly excluded from the SYP61‐positive BFA‐induced TGN aggregates, indicating that homotypic fusion of the TGN rather than de novo endocytic trafficking is important for the formation of TGN/EE‐derived BFA‐induced aggregates. As the TGN also serves as an EE, continuously receiving materials from the plasma membrane, our data support the notion that the secretory Golgi organelle is distinct from the endocytic TGN/EE in terms of its response to BFA treatment in plant cells. Thus, the Golgi and TGN are probably functionally distinct organelles in plants.     

I(Ca(TTX)) channels are distinct from those generating the classical cardiac Na(+) current

The Na(+) current component I(Ca(TTX)) is functionally distinct from the main body of Na(+) current, I(Na). It was proposed that I(Ca(TTX)) channels are I(Na) channels that were altered by bathing media containing Ca(2+), but no, or very little, Na(+). It is known that Na(+)-free conditions are not required to demonstrate I(Ca(TTX).) We show here that Ca(2+) is also not required. Whole-cell, tetrodotoxin-blockable currents from fresh adult rat ventricular cells in 65 mm Cs(+) and no Ca(2+) were compared to those in 3 mM Ca(2+) and no Cs(+) (i.e., I(Ca(TTX))). I(Ca(TTX)) parameters were shifted to more positive voltages than those for Cs(+). The Cs(+) conductance-voltage curve slope factor (mean, -4.68 mV; range, -3.63 to -5.72 mV, eight cells) is indistinguishable from that reported for I(Ca(TTX)) (mean, -4.49 mV; range, -3.95 to -5.49 mV). Cs(+) current and I(Ca(TTX)) time courses were superimposable after accounting for the voltage shift. Inactivation time constants as functions of potential for the Cs(+) current and I(Ca(TTX)) also superimposed after voltage shifting, as did the inactivation curves. Neither of the proposed conditions for conversion of I(Na) into I(Ca(TTX)) channels is required to demonstrate I(Ca(TTX)). Moreover, we find that cardiac Na(+) (H1) channels expressed heterologously in HEK 293 cells are not converted to I(Ca(TTX)) channels by Na(+)-free, Ca(2+)-containing bathing media. The gating properties of the Na(+) current through H1 and those of Ca(2+) current through H1 are identical. All observations are consistent with two non-interconvertable Na(+) channel populations: a larger that expresses little Ca(2+) permeability and a smaller that is appreciably Ca(2+)-permeable.     

Fibroin silk proteins from the nonmulberry silkworm Philosamia ricini are biochemically and immunochemically distinct from those of the mulberry silkworm Bombyx mori

Silk proteins were isolated from the cocoons of the nonmulberry silkworm, Philosamia ricini. Three polypeptides of 97, 66, and 45 kDa were identified. The 66-kDa molecule represented sericin, whereas the 97-kDa and the 45-kDa polypeptides linked together through a disulfide bond constituted the fibroin protein. Antibodies raised against the 97-kDa P. ricini fibroin heavy chain reacted specifically with this molecule and did not recognize fibroin heavy chain from another nonmulberry silkworm, Antheraea assama or from the mulberry silkworm, Bombyx mori, suggesting the presence of P. ricini species-specific determinants in this heavy chain. Antibodies generated against fibroin light chain of P. ricini also showed similar reactivity pattern. Immunoblot analysis with proteins isolated from the silk glands of P. ricini at different stages of larval development showed that the expression of fibroin heavy chain was developmentally and spatially regulated. The protein was most abundant in the 5th instar larva, and could be detected in the middle and the posterior but not the anterior silk glands. The amino acid composition of the 97-kDa fibroin protein showed abundance of glutamic acid and did not contain (Gly-Ala)(n) motifs, a characteristic feature of B. mori fibroin heavy chain. Our study reveals significant differences between the nonmulberry silkworm P. ricini and the mulberry silkworm B. mori in the biochemical composition and immunochemical characteristics of fibroin heavy chain. These differences might be responsible for the differences seen in the quality of silk produced by these two silkworms.     

Responses of belowground communities to large aboveground herbivores: Meta‐analysis reveals biome‐dependent patterns and critical research gaps

The importance of herbivore–plant and soil biota–plant interactions in terrestrial ecosystems is amply recognized, but the effects of aboveground herbivores on soil biota remain challenging to predict. To find global patterns in belowground responses to vertebrate herbivores, we performed a meta‐analysis of studies that had measured abundance or activity of soil organisms inside and outside field exclosures (areas that excluded herbivores). Responses were often controlled by climate, ecosystem type, and dominant herbivore identity. Soil microfauna and especially root‐feeding nematodes were negatively affected by herbivores in subarctic sites. In arid ecosystems, herbivore presence tended to reduce microbial biomass and nitrogen mineralization. Herbivores decreased soil respiration in subarctic ecosystems and increased it in temperate ecosystems, but had no net effect on microbial biomass or nitrogen mineralization in those ecosystems. Responses of soil fauna, microbial biomass, and nitrogen mineralization shifted from neutral to negative with increasing herbivore body size. Responses of animal decomposers tended to switch from negative to positive with increasing precipitation, but also differed among taxa, for instance Oribatida responded negatively to herbivores, whereas Collembola did not. Our findings imply that losses and gains of aboveground herbivores will interact with climate and land use changes, inducing functional shifts in soil communities. To conceptualize the mechanisms behind our findings and link them with previous theoretical frameworks, we propose two complementary approaches to predict soil biological responses to vertebrate herbivores, one focused on an herbivore body size gradient, and the other on a climate severity gradient. Major research gaps were revealed, with tropical biomes, protists, and soil macrofauna being especially overlooked.