Matrix models for quantifying competitive intransitivity from species abundance data

In a network of competing species, a competitive intransitivity occurs when the ranking of competitive abilities does not follow a linear hierarchy (A > B > C but C > A). A variety of mathematical models suggests that intransitive networks can prevent or slow down competitive exclusion and maintain biodiversity by enhancing species coexistence. However, it has been difficult to assess empirically the relative importance of intransitive competition because a large number of pairwise species competition experiments are needed to construct a competition matrix that is used to parameterize existing models. Here we introduce a statistical framework for evaluating the contribution of intransitivity to community structure using species abundance matrices that are commonly generated from replicated sampling of species assemblages. We provide metrics and analytical methods for using abundance matrices to estimate species competition and patch transition matrices by using reverse-engineering and a colonization-competition model. These matrices provide complementary metrics to estimate the degree of intransitivity in the competition network of the sampled communities. Benchmark tests reveal that the proposed methods could successfully detect intransitive competition networks, even in the absence of direct measures of pairwise competitive strength. To illustrate the approach, we analyzed patterns of abundance and biomass of five species of necrophagous Diptera and eight species of their hymenopteran parasitoids that co-occur in beech forests in Germany. We found evidence for a strong competitive hierarchy within communities of flies and parasitoids. However, for parasitoids, there was a tendency towards increasing intransitivity in higher weight classes, which represented larger resource patches. These tests provide novel methods for empirically estimating the degree of intransitivity in competitive networks from observational datasets. They can be applied to experimental measures of pairwise species interactions, as well as to spatio-temporal samples of assemblages in homogenous environments or environmental gradients.

Biotic and Abiotic Properties Mediating Plant Diversity Effects on Soil Microbial Communities in an Experimental Grassland

Plant diversity drives changes in the soil microbial community which may result in alterations in ecosystem functions. However, the governing factors between the composition of soil microbial communities and plant diversity are not well understood. We investigated the impact of plant diversity (plant species richness and functional group richness) and plant functional group identity on soil microbial biomass and soil microbial community structure in experimental grassland ecosystems. Total microbial biomass and community structure were determined by phospholipid fatty acid (PLFA) analysis. The diversity gradient covered 1, 2, 4, 8, 16 and 60 plant species and 1, 2, 3 and 4 plant functional groups (grasses, legumes, small herbs and tall herbs). In May 2007, soil samples were taken from experimental plots and from nearby fields and meadows. Beside soil texture, plant species richness was the main driver of soil microbial biomass. Structural equation modeling revealed that the positive plant diversity effect was mainly mediated by higher leaf area index resulting in higher soil moisture in the top soil layer. The fungal-to-bacterial biomass ratio was positively affected by plant functional group richness and negatively by the presence of legumes. Bacteria were more closely related to abiotic differences caused by plant diversity, while fungi were more affected by plant-derived organic matter inputs. We found diverse plant communities promoted faster transition of soil microbial communities typical for arable land towards grassland communities. Although some mechanisms underlying the plant diversity effect on soil microorganisms could be identified, future studies have to determine plant traits shaping soil microbial community structure. We suspect differences in root traits among different plant communities, such as root turnover rates and chemical composition of root exudates, to structure soil microbial communities.

Moving forward on facilitation research: response to changing environments and effects on the diversity, functioning and evolution of plant communities

Once seen as anomalous, facilitative interactions among plants and their importance for community structure and functioning are now widely recognized. The growing body of modelling, descriptive and experimental studies on facilitation covers a wide variety of terrestrial and aquatic systems throughout the globe. However, the lack of a general body of theory linking facilitation among different types of organisms and biomes and their responses to environmental changes prevents further advances in our knowledge regarding the evolutionary and ecological implications of facilitation in plant communities. Moreover, insights gathered from alternative lines of inquiry may substantially improve our understanding of facilitation, but these have been largely neglected thus far. Despite over 15 years of research and debate on this topic, there is no consensus on the degree to which plant–plant interactions change predictably along environmental gradients (i.e. the stress-gradient hypothesis), and this hinders our ability to predict how plant–plant interactions may affect the response of plant communities to ongoing global environmental change. The existing controversies regarding the response of plant–plant interactions across environmental gradients can be reconciled when clearly considering and determining the species-specificity of the response, the functional or individual stress type, and the scale of interest (pairwise interactions or community-level response). Here, we introduce a theoretical framework to do this, supported by multiple lines of empirical evidence. We also discuss current gaps in our knowledge regarding how plant–plant interactions change along environmental gradients. These include the existence of thresholds in the amount of species-specific stress that a benefactor can alleviate, the linearity or non-linearity of the response of pairwise interactions across distance from the ecological optimum of the beneficiary, and the need to explore further how frequent interactions among multiple species are and how they change across different environments. We review the latest advances in these topics and provide new approaches to fill current gaps in our knowledge. We also apply our theoretical framework to advance our knowledge on the evolutionary aspects of plant facilitation, and the relative importance of facilitation, in comparison with other ecological processes, for maintaining ecosystem structure, functioning and dynamics. We build links between these topics and related fields, such as ecological restoration, woody encroachment, invasion ecology, ecological modelling and biodiversity–ecosystem-functioning relationships. By identifying commonalities and insights from alternative lines of research, we further advance our understanding of facilitation and provide testable hypotheses regarding the role of (positive) biotic interactions in the maintenance of biodiversity and the response of ecological communities to ongoing environmental changes.

Shifts in microbial communities do not explain the response of grassland ecosystem function to plant functional composition and rainfall change

Ecosystem functioning in grasslands is regulated by a range of biotic and abiotic factors, and the role of microbial communities in regulating ecosystem function has been the subject of much recent scrutiny. However, there are still knowledge gaps regarding the impacts of rainfall and vegetation change upon microbial communities and the implications of these changes for ecosystem functioning. We investigated this issue using data from an experimental mesotrophic grassland study in south-east England, which had been subjected to four years of rainfall and plant functional composition manipulations. Soil respiration, nitrogen and phosphorus stocks were measured, and the abundance and community structure of soil microbes were characterised using quantitative PCR and multiplex-TRFLP analysis, respectively. Bacterial community structure was strongly related to the plant functional composition treatments, but not the rainfall treatment. However, there was a strong effect of both rainfall change and plant functional group upon bacterial abundance. There was also a weak interactive effect of the two treatments upon fungal community structure, although fungal abundance was not affected by either treatment. Next, we used a statistical approach to assess whether treatment effects on ecosystem function were regulated by the microbial community. Our results revealed that ecosystem function was influenced by the experimental treatments, but was not related to associated changes to the microbial community. Overall, these results indicate that changes in fungal and bacterial community structure and abundance play a relatively minor role in determining grassland ecosystem function responses to precipitation and plant functional composition change, and that direct effects on soil physical and chemical properties and upon plant and microbial physiology may play a more important role.

A conserved pattern in plant-mediated interactions between herbivores

Plant‐mediated interactions between herbivores are important determinants of community structure and plant performance in natural and agricultural systems. Current research suggests that the outcome of the interactions is determined by herbivore and plant identity, which may result in stochastic patterns that impede adaptive evolution and agricultural exploitation. However, few studies have systemically investigated specificity versus general patterns in a given plant system by varying the identity of all involved players. We investigated the influence of herbivore identity and plant genotype on the interaction between leaf‐chewing and root‐feeding herbivores in maize using a partial factorial design. We assessed the influence of leaf induction by oral secretions of six different chewing herbivores on the response of nine different maize genotypes and three different root feeders. Contrary to our expectations, we found a highly conserved pattern across all three dimensions of specificity: The majority of leaf herbivores elicited a negative behavioral response from the different root feeders in the large majority of tested plant genotypes. No facilitation was observed in any of the treatment combinations. However, the oral secretions of one leaf feeder and the responses of two maize genotypes did not elicit a response from a root‐feeding herbivore. Together, these results suggest that plant‐mediated interactions in the investigated system follow a general pattern, but that a degree of specificity is nevertheless present. Our study shows that within a given plant species, plant‐mediated interactions between herbivores of the same feeding guild can be stable. This stability opens up the possibility of adaptations by associated organisms and suggests that plant‐mediated interactions may contribute more strongly to evolutionary dynamics in terrestrial (agro)ecosystems than previously assumed.

Does proximity to conspecific adults influence the establishment of ectomycorrhizal trees in rain forest?

Three ectomycorrhizal legume trees, Microberlinia bisulcata, Tetraberlinia bifoliolata and T. moreliana, form discrete groves in the southern part of Korup National Park, in southwest Cameroon and contribute c. 45–70% of stand basal area locally in a matrix of otherwise species-rich arbuscular mycorrhizal forest. A transplant experiment was performed to assess the importance of ectomycorrhizal infection associated with proximity to parents in seedling establishment of the grove-forming species. Nonectomycorrhizal seedlings of the three species were transplanted into plots of two forest types, one of high (HEM, within-grove) and one of very low (LEM, outside the grove) abundance of all three species as adult trees. For two species (T. moreliana and M. bisulcata) there was no difference in survival over 16 months, but for the third (T. bifoliolata) survival was best in HEM forest, and correlated with the basal area of adult trees of ectomycorrhizal species. Only one species (T. moreliana) increased in biomass over the experimental period; the others declined. There was no effect of forest type on overall growth of any species, but the survivors of two (T. moreliana and M. bisulcata) had heavier stems in the HEM forest. Differences in survival and growth of transplants between the three species were in accord with the ecology of the species as inferred from the frequency distributions of adult tree size in the forest. Seedlings became infected with ectomycorrhizas in both forest types; where there was a difference in extent of infection (T. moreliana) this was not related to survival or growth; and where there was a difference in survival (T. bifoliolata) this was not related to extent of infection. These results confirm that mycorrhizal inoculum associated with conspecific adults is neither a prerequisite nor a guarantee of seedling establishment, but indicates that in some circumstances there might be benefits of being close to parents. Further research is required to unravel the complexities of ectomycorrhizal community structure in this spatially and temporally heterogeneous forest, and to clarify the extent to which the various hosts share ectomycorrhizal partners.

A global seasonal surface ocean climatology of phytoplankton types based on CHEMTAX analysis of HPLC pigments

Much advancement has been made in recent years in field data assimilation, remote sensing and ecosystem modeling, yet our global view of phytoplankton biogeography beyond chlorophyll biomass is still a cursory taxonomic picture with vast areas of the open ocean requiring field validations. High performance liquid chromatography (HPLC) pigment data combined with inverse methods offer an advantage over many other phytoplankton quantification measures by way of providing an immediate perspective of the whole phytoplankton community in a sample as a function of chlorophyll biomass. Historically, such chemotaxonomic analysis has been conducted mainly at local spatial and temporal scales in the ocean. Here, we apply a widely tested inverse approach, CHEMTAX, to a global climatology of pigment observations from HPLC. This study marks the first systematic and objective global application of CHEMTAX, yielding a seasonal climatology comprised of ~1500 1°x1° global grid points of the major phytoplankton pigment types in the ocean characterizing cyanobacteria, haptophytes, chlorophytes, cryptophytes, dinoflagellates, and diatoms, with results validated against prior regional studies where possible. Key findings from this new global view of specific phytoplankton abundances from pigments are a) the large global proportion of marine haptophytes (comprising 32 ± 5% of total chlorophyll), whose biogeochemical functional roles are relatively unknown, and b) the contrasting spatial scales of complexity in global community structure that can be explained in part by regional oceanographic conditions. These publicly accessible results will guide future parameterizations of marine ecosystem models exploring the link between phytoplankton community structure and marine biogeochemical cycles.

Abundance of mesozooplankton in the western part of the Black Sea in summer 1998 during the National Monitoring Programme

The dataset is based on samples collected in the summer of 1998 in the Western Black Sea in front of Bulgaria coast. The whole dataset is composed of 69 samples (from 22 stations of National Monitoring Grid) with data of mesozooplankton species composition abundance and biomass. Samples were collected in discrete layers 0-10, 0-20, 0-50, 10-25, 25-50, 50-100 and from bottom up to the surface at depths depending on water column stratification and the thermocline depth. Zooplankton samples were collected with vertical closing Juday net,diameter - 36cm, mesh size 150 µm. Tows were performed from surface down to bottom meters depths in discrete layers. Samples were preserved by a 4% formaldehyde sea water buffered solution. Sampling volume was estimated by multiplying the mouth area with the wire length. Mesozooplankton abundance: The collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber which is a rectangle form for taxomomic identification and count. Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Oceanology by Lyudmila Kamburska using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972). Taxon-specific abundance: The collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber which is a rectangle form for taxomomic identification and count. Copepods and Cladoceras were identified and enumerated; the other mesozooplankters were identified and enumerated at higher taxonomic level (commonly named as mesozooplankton groups). Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Oceanology by Lyudmila Kamburska using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972).

Abundance of mesozooplankton in the western part of the Black Sea in summer 2001 during the National Monitoring Programme of the Bulgarian Institute of Oceanography

The dataset is based on samples collected in the summer of 2001 in the Western Black Sea in front of Bulgaria coast (transects at c. Kaliakra and c. Galata). The whole dataset is composed of 26 samples (from 10 stations of National Monitoring Grid) with data of mesozooplankton species composition abundance and biomass. Samples were collected in discrete layers 0-10, 10-20, 10-25, 25-50, 50-75, 75-90. Zooplankton samples were collected with vertical closing Juday net,diameter - 36cm, mesh size 150 µm. Tows were performed from surface down to bottom meters depths in discrete layers. Samples were preserved by a 4% formaldehyde sea water buffered solution. Sampling volume was estimated by multiplying the mouth area with the wire length. Mesozooplankton abundance: The collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber which is a rectangle form for taxomomic identification and count. Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Oceanology by Lyudmila Kamburska and Kremena Stefanova using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972). Taxon-specific abundance: The collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber which is a rectangle form for taxomomic identification and count. Copepods and Cladoceras were identified and enumerated; the other mesozooplankters were identified and enumerated at higher taxonomic level (commonly named as mesozooplankton groups). Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Oceanology by Lyudmila Kamburska and Kremena Stefanova using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972).

(Data S1) Palynological analysis of IODP and ODP sediment cores around the Antarctic margin

The circum-Antarctic Southern Ocean is an important region for global marine food webs and carbon cycling because of sea-ice formation and its unique plankton ecosystem. However, the mechanisms underlying the installation of this distinct ecosystem and the geological timing of its development remain unknown. Here, we show, on the basis of fossil marine dinoflagellate cyst records, that a major restructuring of the Southern Ocean plankton ecosystem occurred abruptly and concomitant with the first major Antarctic glaciation in the earliest Oligocene (~33.6 million years ago). This turnover marks a regime shift in zooplankton-phytoplankton interactions and community structure, which indicates the appearance of eutrophic and seasonally productive environments on the Antarctic margin. We conclude that earliest Oligocene cooling, ice-sheet expansion, and subsequent sea-ice formation were important drivers of biotic evolution in the Southern Ocean.

Soil characteristics and Collembola and Acari abundance in control and warming plots at Abisco Research Station

Extreme weather events can have negative impacts on species survival and community structure when surpassing lethal thresholds. Extreme winter warming events in the Arctic rapidly melt snow and expose ecosystems to unseasonably warm air (2-10 °C for 2-14 days), but returning to cold winter climate exposes the ecosystem to lower temperatures by the loss of insulating snow. Soil animals, which play an integral part in soil processes, may be very susceptible to such events depending on the intensity of soil warming and low temperatures following these events. We simulated week-long extreme winter warming events - using infrared heating lamps, alone or with soil warming cables - for two consecutive years in a sub-Arctic dwarf shrub heathland. Minimum temperatures were lower and freeze-thaw cycles were 2-11 times more frequent in treatment plots compared with control plots. Following the second event, Acari populations decreased by 39%; primarily driven by declines of Prostigmata (69%) and the Mesostigmatic nymphs (74%). A community-weighted vertical stratification shift occurred from smaller soil dwelling (eu-edaphic) Collembola species dominance to larger litter dwelling (hemi-edaphic) species dominance in the canopy-with-soil warming plots compared with controls. The most susceptible groups to these winter warming events were the smallest individuals (Prostigmata and eu-edaphic Collembola). This was not apparent from abundance data at the Collembola taxon level, indicating that life forms and species traits play a major role in community assembly following extreme events. The observed shift in soil community can cascade down to the micro-flora affecting plant productivity and mineralization rates. Short-term extreme weather events have the potential to shift community composition through trait composition with potentially large consequences for ecosystem development.

Diversity of zooplankton assemblages at Station L4, Western English Channel, measured using next generation DNA sequencing

Background: Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel Next Generation Sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify the richness and diversity of a mixed zooplankton assemblage from a productive monitoring site in the Western English Channel. Methodology/Principle Findings: Plankton WP2 replicate net hauls (200 µm) were taken at the Western Channel Observatory long-term monitoring station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,042 sequences were obtained for all samples. The sequences clustered in to 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 138 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 75 taxonomic groups. Conclusions: The percentage of OTUs assigned to major eukaryotic taxonomic groups broadly aligns between the metagenetic and morphological analysis and are dominated by Copepoda. However, the metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for estimating diversity and species richness of zooplankton communities.