Bioengineers and their associated fauna respond differently to the effects of biogeography and upwelling

Temperature and primary production (often linked to nutrient supply) are two of the few factors influencing species diversity and abundances across mesoscale gradients, while at smaller scales the habitat complexity offered by bioengineers is important. Previous studies have illustrated effects of upwelling and biogeography on intertidal bioengineers, but it is not known if these processes influence assemblages associated with those bioengineers in a similar way. We examined the habitat structure offered by two species of mussels and their associated fauna in five regions across 3000 km and three biogeographic provinces of the South African coast, replicating upwelling and non-upwelling areas within each region. Upwelling and region influenced the structure of mussel beds (the density and size of mussels). In contrast, upwelling did not influence mesoscale differences in composition, abundance, and numbers of species of crustaceans, mollusks, or polychaetes in mussel beds. Regardless of trophic level or mode of reproduction, mussel bed fauna were influenced only by region. Regional differences were strongly influenced by biogeography. The associated fauna was, however, also strongly correlated with the structure of the habitats created by mussels. Our results support the importance of upwelling to a critical ecosystem engineer, but show that these effects do not extend directly to the assemblages of associated fauna, which are more influenced by regional-scale effects and biogeography. We suggest that mesoscale patterns in the associated fauna of this bioengineered habitat are driven by the direct effects of biogeography, combined with the influence of biogeography and upwelling on mussel bed structure.     

Habitat structure and complexity as determinants of biodiversity in blue mussel beds on sublittoral rocky shores

Habitat-forming, ecosystem engineer species are common in most marine systems. Still, much uncertainty exists about how individual and population-level traits of these species contribute to ecosystem processes and how engineering species jointly affect biodiversity. In this manipulative field experiment, we examined how biodiversity in marginal blue mussel beds is affected by blue mussel (1) body size, density and patch context and (2) presence of fucoid and algal structures. In the study area, bladder-wrack (Fucus vesiculosus), filamentous algae and blue mussels (Mytilus edulis) coexist at shallow depths in a variety of patch configurations and offer complex habitats with a high variability of resources. We hypothesized that complexity in terms of mussel bed structure and algal presence determines species composition and abundance. Results from the experiment were compared with macrofaunal communities found in natural populations of both engineering species. Results show that the physical structure and blue mussel patch context are important determinants for species composition and abundance. Results further show that the presence of algal structures positively affects diversity in blue mussel habitats due to increased surface availability and complexity that these algae offer. This study shows that blue mussel beds at the very margin of their distribution have an indisputable function for promoting and maintaining biodiversity and suggest that facilitative effects of habitat-modifying species are important on Baltic Sea rocky shores with fundamental importance to community structure.     

Effects of Changes in Biodiversity on Ecosystem Function in Tropical Forests

Changing land use in the tropics has resulted in vast areas of damaged and degraded lands where biodiversity has been reduced. The majority of research on biodiversity has been focused on population and community dynamics and has rarely considered the ecosystem processes that are intimately related. We present a framework for examining the effects of changes in biodiversity on ecosystem function in natural, managed, and damaged tropical forests. Using a whole-ecosystem approach, the framework identifies key nutrient and energy cycling processes and critical junctures or pathways, termed interfaces, where resources are concentrated and transferred between the biotic and abiotic components of the ecosystem. Processes occurring at these interfaces, and the organisms or attributes participating in these processes, exert a strong influence on ecosystem structure. We use examples from Puerto Rico, Southern China, Dominica, and Nicaragua to illustrate how the functional diversity framework can be applied to critically examine the effects of changes in biodiversity on ecosystem function, and the relative success or failure of rehabilitation strategies. The few available data suggest that functional diversity, and not just species richness, is important in maintaining the integrity of nutrient and energy fluxes. High species richness, however, may increase ecosystem resiliency following disturbance by increasing the number of alternative pathways for the flow of resources. We suggest ways in which the framework of functional diversity can be used to design research to examine the effects of changes in biodiversity on ecosystem processes and in the design and evaluation of ecosystem management and land rehabilitation projects in the tropics.     

Local extinction of a foundation species in a hypoxic estuary: integrating individuals to ecosystem

We integrated across individual, population, community, and ecosystem levels to understand the impact of environmental stress by tracking the foundation species Mytilus edulis in the hypoxic estuary Narragansett Bay, Rhode Island, USA. Our initial surveys revealed that the mussels occurred in nine extensive (2-28 ha) dense (814-9943 individuals/m2) subtidal reefs that attracted a diverse suite of predators (sea stars, crabs, gastropods). Hypoxia occurred in the summer of 2001, and a mussel transplant experiment revealed overall reduced growth rates of individuals, and higher mortality rates among larger mussels. At the population level, large decreases in densities and cover of mussels were correlated with dissolved oxygen concentrations, leading to extinction at one site and reductions of over an order of magnitude at others. Within one year, seven of the eight remaining populations were edged to extinction, and the previously extinct population was recolonized. At the community level, a predator exclusion experiment indicated that predation was an unimportant source of mussel mortality during the hypoxic period, in part due to the emigration of sea stars, as predicted by the Consumer Stress Model. However, mussels were too intolerant to hypoxia to have a net benefit from the predation refuge. The seasonal (summer) occurrence of hypoxia allowed sea stars to return following a lag, as predicted by a stress return time model, and the resumption of predation contributed to the subsequent extinction of mussel populations. At the ecosystem level, the initial filtration rate of the mussel reefs was estimated at 134.6 x 10(6) m3/d, equivalent to filtering the volume of the bay 1.3 times during the 26-d average residence time. That function was reduced by >75% following hypoxia. The effect of hypoxia on each level of organization had consequences at others. For example, size-specific mortality and decreased growth of individuals, and reduced filtration capacity of reefs, indicated a loss of the ability of mussels to entrain planktonic productivity and potential to control future eutrophication and hypoxia. Our study quantified patterns of loss and identified pathways within an integrative framework of feedbacks, summarized in a conceptual model that is applicable to similar foundation species subjected to environmental stress.     

Parameter uncertainties in the modelling of vegetation dynamics—Effects on tree community structure and ecosystem functioning in European forest biomes

Dynamic vegetation models are useful tools for analysing terrestrial ecosystem processes and their interactions with climate through variations in carbon and water exchange. Long-term changes in structure and composition (vegetation dynamics) caused by altered competitive strength between plant functional types (PFTs) are attracting increasing attention as controls on ecosystem functioning and potential feedbacks to climate. Imperfect process knowledge and limited observational data restrict the possibility to parameterise these processes adequately and potentially contribute to uncertainty in model results. This study addresses uncertainty among parameters scaling vegetation dynamic processes in a process-based ecosystem model, LPJ-GUESS, designed for regional-scale studies, with the objective to assess the extent to which this uncertainty propagates to additional uncertainty in the tree community structure (in terms of the tree functional types present and their relative abundance) and thus to ecosystem functioning (carbon storage and fluxes). The results clearly indicate that the uncertainties in parameterisation can lead to a shift in competitive balance, most strikingly among deciduous tree PFTs, with dominance of either shade-tolerant or shade-intolerant PFTs being possible, depending on the choice of plausible parameter values. Despite this uncertainty, our results indicate that the resulting effect on ecosystem functioning is low. Since the vegetation dynamics in LPJ-GUESS are representative for the more complex Earth system models now being applied within ecosystem and climate research, we assume that our findings will be of general relevance. We suggest that, in terms of carbon storage and fluxes, the heavier parameterisation requirement of the processes involved does not widen the overall uncertainty in model predictions.     

Intertidal invertebrates locally enhance primary production

The contribution of autochthonous vs. allochthonous inputs to productivity is an important determinant of ecosystem function across multiple habitats. In coastal marine systems, nutrients are thought to come primarily from the upwelling of deep, nutrient-rich water. Using experimental manipulations of a dominant tide pool animal, the mussel Mytilus californianus, I show that the presence of mussels greatly increases the supply of inorganic nitrogen and phosphorus. Mussels further had a direct effect on productivity: benthic microalgal abundance increased by a factor of 4-8, while the growth of a red alga was four times greater in the presence of mussels. The increase in nitrite and nitrate associated with mussels further suggests nitrifying activity by microbes. These findings have broad implications for coastal marine systems, including that regenerated nutrients may contribute more to productivity than previously recognized and that the presence of animal-generated nutrients sets the stage for numerous positive interactions.     

Dramatic declines in mussel bed community diversity: response to climate change?

Mussel beds along the wave-exposed coast of the eastern North Pacific Ocean serve as an important habitat, harboring a high diversity of species. A comparison of California mussel bed community diversity in 2002 to historical data (1960s to 1970s) revealed large declines (mean loss 58.9%), including some declines >141 species (approximately 80% loss). Concurrent work revealed inconsistent changes in mussel populations (biomass and bed thickness) along the California coast, suggesting that diversity declines may be related to large-scale processes rather than local habitat destruction. Potential factors causing declines in mussel community diversity are discussed, with regional climate change associated with the Pacific Decadal Oscillation and climate change induced alterations of ecological interactions and biological processes suggested as likely causes. Although extensive literature has predicted the potential effects of climate change on global diversity, this study is one of the few examples of declines attributed to climate change.     

A note on the sufficiency of low resolution ecosystem models

Analysis of a system of non-linear differential equations illustrates the effects of interactions between biotic and abiotic components of a complex aquatic ecosystem model. A stochastic analysis shows that the variance of the abiotic variables is related in a simple manner to the autocorrelation function of the biotic variables. The results suggest that for oligotrophic and eutrophic conditions, relatively simple ecosystem models may be sufficient for studies of an aquatic environment. Under mesotrophic conditions, the high state variable resolution of a complex model may be necessary.     

Long-Term Ecosystem Dynamics in the Serengeti: Lessons for Conservation

Abstract:  Data from long-term ecological studies further understanding of ecosystem dynamics and can guide evidence-based management. In a quasi-natural experiment we examined long-term monitoring data on different components of the Serengeti-Mara Ecosystem to trace the effects of disturbances and thus to elucidate cause-and-effect connections between them. The long-term data illustrated the role of food limitation in population regulation in mammals, particularly in migratory wildebeest and nonmigratory buffalo. Predation limited populations of smaller resident ungulates and small carnivores. Abiotic events, such as droughts and floods, created disturbances that affected survivorship of ungulates and birds. Such disturbances showed feedbacks between biotic and abiotic realms. Interactions between elephants and their food allowed savanna and grassland communities to co-occur. With increased woodland vegetation, predators' capture of prey increased. Anthropogenic disturbances had direct (hunting) and indirect (transfer of disease to wildlife) effects. Slow and rapid changes and multiple ecosystem states became apparent only over several decades and involved events at different spatial scales. Conservation efforts should accommodate both infrequent and unpredictable events and long-term trends. Management should plan on the time scale of those events and should not aim to maintain the status quo. Systems can be self-regulating through food availability and predator-prey interactions; thus, culling may not be required. Ecosystems can occur in multiple states; thus, there may be no a priori need to maintain one natural state. Finally, conservation efforts outside protected areas must distinguish between natural change and direct human-induced change. Protected areas can act as ecological baselines in which human-induced change is kept to a minimum     

Hierarchical analysis of species distributions and abundance across environmental gradients

Abiotic and biotic processes operate at multiple spatial and temporal scales to shape many ecological processes, including species distributions and demography. Current debate about the relative roles of niche-based and stochastic processes in shaping species distributions and community composition reflects, in part, the challenge of understanding how these processes interact across scales. Traditional statistical models that ignore autocorrelation and spatial hierarchies can result in misidentification of important ecological covariates. Here, we demonstrate the utility of a hierarchical modeling framework for testing hypotheses about the importance of abiotic factors at different spatial scales and local spatial autocorrelation for shaping species distributions and abundances. For the two orchid species studied, understory light availability and soil moisture helped to explain patterns of presence and abundance at a microsite scale (<4 m2), while soil organic content was important at a population scale (<400 m2). The inclusion of spatial autocorrelation is shown to alter the magnitude and certainty of estimated relationships between abundance and abiotic variables, and we suggest that such analysis be used more often to explore the relationships between species life histories and distributions. The hierarchical modeling framework is shown to have great potential for elucidating ecological relationships involving abiotic and biotic processes simultaneously at multiple scales.     

A methodological framework for integrating computational fluid dynamics and ecological models applied to juvenile freshwater mussel dispersal in the Upper Mississippi River

Interdisciplinary research in hydraulics and ecology for river management and restoration must integrate processes that occur over a wide range of spatial and temporal scales, which presents a challenge to ecohydraulics modelers. Computational fluid dynamics (CFD) models are being more widely used to determine flow fields for ecohydraulics applications. In the Upper Mississippi River (UMR), the mussel dynamics model was developed as a tool for management and conservation of freshwater mussels (Unionidae), which are benthic organisms, imperiled in North America, that are inextricably linked with the hydraulics of river flow. We updated the juvenile dispersal component of the mussel dynamics model by using stochastic Lagrangian particle tracking in a three dimensional flow field output from CFD models of reaches in the UMR. We developed a methodological framework to integrate hydrodynamic data with the mussel dynamics model, and we demonstrate the use of the juvenile dispersal model employed within the methodological framework in two reaches of the UMR. The method was used to test the hypothesis that impoundment affects the relationship of some hydraulic parameters with juvenile settling distribution. Simulation results were consistent with this hypothesis, and the relationships of bed shear stress and Froude number with juvenile settling were altered by impoundment most likely through effects on local hydraulics. The methodological framework is robust, integrates Eulerian and Lagrangian reference frameworks, and incorporates processes over a wide range of temporal and spatial scales, from watershed scale hydrologic processes (decades), to reach scale (km) processes that occur over hours or days, and turbulent processes on spatial scales of meter to millimeter and times scales of seconds. The methods are presently being used to assess the impacts of pre- and early post-settlement processes on mussel distributions, including the effects of bed shear stress, and the sensitivity of the location of the host fish when juveniles excyst, on juvenile settling distribution.     

The land-cover cascade: relationships coupling land and water

We introduce the land-cover cascade (LCC) as a conceptual framework to quantify the transfer of land-cover-disturbance effects to stream biota. We hypothesize that disturbance is propagated through multivariate systems through key variables that transform a disturbance and pass a reorganized disturbance effect to the next hierarchical level where the process repeats until ultimately affecting biota. We measured 31 hydrologic, geomorphic, erosional, and substrate variables and 26 biotic responses that have been associated with land-use disturbance in third- and fourth-order streams in the Blue Ridge physiographic province in western North Carolina (USA). Regression analyses reduced this set of variables to include only those that responded to land cover and/or affected biota. From this reduced variable set, hypotheses were generated that predicted the disturbance pathways affecting each biotic response following the land-cover-cascade design. Cascade pathways began with land cover and ended with biotic responses, passing through at least one intermediate ecosystem abiotic component. Cascade models were tested for predictive ability and goodness-of-fit using path analysis. Biota were influenced by near-stream urban, agricultural, and forest land cover as propagated by hydrologic (e.g., discharge), geomorphic (e.g., stream bank height), erosional (e.g., suspended sediments), and depositional streambed (e.g., substrate size) features occurring along LCC pathways, reflecting abiotic mechanisms mediating land-cover disturbance. Our results suggest that communities are influenced by land-cover change indirectly through a hierarchy of associated abiotic components that propagate disturbance to biota. More generally, the land-cover cascade concept and experimental framework demonstrate an organized approach to the generic study of cascades and the complex relationships between landscapes and streams.     

A box model of carrying capacity for suspended mussel aquaculture in Lagune de la Grande-Entrée,Iles-de-la-Madeleine,Québec

An object-oriented model of environment–mussel aquaculture interactions and mussel carrying-capacity within Lagune de la Grande-Entrée (GEL), Iles-de-la-Madeleine, Québec, was constructed to assist in development of sustainable mussel culture in this region. A multiple box ecosystem model for GEL tied to the output of a hydrodynamic model was constructed using Simile software, which has inherent ability to represent spatial elements and specify water exchange between modelled regions. Mussel growth and other field data were used for model validation. Plackett–Burman sensitivity analysis demonstrated that a variety of bioenergetic parameters of zooplankton and phytoplankton submodels were important in model outcomes. Model results demonstrated that mussel aquaculture can be further developed throughout the lagoon. At present culture densities, phytoplankton depletion is minimal, and there is little food limitation of mussel growth. Results indicated that increased stocking density of mussels in the existing farm will lead to decreased mass per individual mussel. Depending on the location of new farm emplacement within the lagoon, implementation of new aquaculture sites either reduced mussel growth in the existing farm due to depletion of phytoplankton, or exhibited minimum negative impact on the existing farm. With development throughout GEL, an excess of phytoplankton was observed during the year in all modelled regions, even at stocking densities as high as 20 mussels m⁻³. Although mussels cultured at this density do not substantially impact the ecosystem, their growth is controlled by the flux of phytoplankton food and abundance of zooplankton competitors. This model provides an effective tool to examine expansion of shellfish farming to new areas, balancing culture location and density.     

Shift in nutrient and plankton community in eutrophic lake following introduction of a freshwater bivalve

The impact of the freshwater bivalve Corbicula leana on plankton community dynamics was examined during a cyanobacterial bloom period. Nutrient and chlorophyll concentrations, primary productivity, and phytoplankton and zooplankton communities in the experimental enclosures were measured at 2-3 day intervals. The introduction of mussels reduced net primary productivity and phytoplankton and chlorophyll. Chlorophyll decreased immediately following addition of 100 mussels and then increased over time. After 600 mussels were added, chlorophyll decreased continuously from 87to 25 microg l(-1), approaching that in the mussel-free enclosure. Simultaneously, water transparency increased and concentrations of suspended solids and total phosphorus decreased. Mussel addition caused short-term increases in nutrient concentrations, especially following high-density treatment: phytoplankton density decreased, while cell density in the mussel-free enclosure increased. Zooplankton densities in the two enclosures were similar; however, carbon biomass in the mussel enclosure increased, associated with an increase in large zooplankton. The trophic relationship between phytoplankton and zooplankton was positive in the mussel-free enclosure and negative in the mussel-treatment enclosure, possibly reflecting effects of mussels on both consumer and resource control in the plankton community. Thus, filter feeding by Corbicula affects nutrient recycling and plankton community structure in a freshwater ecosystem through direct feeding and competition for food resources.     

Detrimental effects of the isopod, <Emphasis Type="Italic">Edotia doellojuradoi,</Emphasis> on gill morphology and host condition of the mussel, <Emphasis Type="Italic">Mytilus edulis platensis</Emphasis>

There are many reported associations between mussels and other invertebrates, such as pea crabs, polychaetes, turbellarians and copepods, which live in their mantle cavities. The boundary between commensalism and parasitism is often indistinguishable because of insufficient knowledge or because the interaction is variable. Preliminary evidence led to a closer examination of the relationship between the mussel, Mytilus edulis platensis, and an isopod, Edotia doellojuradoi, previously described as commensalism. Monthly intertidal samples of mussels were taken from September 2004 to August 2005 at Caleta Cordova Norte (45°43′S, 67°22′W) in southern Argentina and assessed for the prevalence and abundance levels of isopods. Mussels with and without isopods were measured, examined for evidence of gill damage and their condition (soft tissue dry weight) was determined. The overall isopod prevalence in mussels was 57.9% and infestation varied with mussel length, with maximum occurence at 30.2 mm (medium-sized mussels). Experimental evidence indicated that the position of the isopod inside the mussel depended on the feeding activity of the mussel. Female isopods were observed grasping the ventral food groove of the gill demibranchs and feeding on the mucous food strand produced by the mussel. Juveniles and males were observed clustered together on the dorsal side of the single female in each occupied mussel, suggesting extended maternal care. Gill damage was observed in 58.2% of mussels at the Argentine site and was significantly associated with isopod occurrence. Categorical regression analysis showed that the most important factor associated with the degree of gill lesions was the number of male and juvenile isopods per mussel, followed by the length of female isopods and the developmental stage of juveniles. Conversely, the degree of gill damage decreased with increasing mussel length. Overall, E. doellojuradoi had a significant effect on mussel condition throughout the year, with low flesh weight in mussels with isopods, except during the austral summer and early autumn. In contrast to previous studies, which concluded that the isopod was a commensal, the present study clearly demonstrates that E. doellojuradoi is a parasite of M. edulis platensis. Other symbiotic interactions formerly classified as commensal might not be innocuous on further investigation, especially if samples are taken at multiple sites and at different times of the year.     

Biophysical feedback mediates effects of invasive grasses on coastal dune shape

Vegetation at the aquatic-terrestrial interface can alter landscape features through its growth and interactions with sediment and fluids. Even similar species may impart different effects due to variation in their interactions and feedbacks with the environment. Consequently, replacement of one engineering species by another can cause significant change in the physical environment. Here we investigate the species-specific ecological mechanisms influencing the geomorphology of U.S. Pacific Northwest coastal dunes. Over the last century, this system changed from open, shifting sand dunes with sparse vegetation (including native beach grass, Elymus mollis), to densely vegetated continuous foredune ridges resulting from the introduction and subsequent invasions of two nonnative grass species (Ammophila arenaria and Ammophila breviligulata), each of which is associated with different dune shapes and sediment supply rates along the coast. Here we propose a biophysical feedback responsible for differences in dune shape, and we investigate two, non-mutually exclusive ecological mechanisms for these differences: (1) species differ in their ability to capture sand and (2) species differ in their growth habit in response to sand deposition. To investigate sand capture, we used a moveable bed wind tunnel experiment and found that increasing tiller density increased sand capture efficiency and that, under different experimental densities, the native grass had higher sand capture efficiency compared to the Ammophila congeners. However, the greater densities of nonnative grasses under field conditions suggest that they have greater potential to capture more sand overall. We used a mesocosm experiment to look at plant growth responses to sand deposition and found that, in response to increasing sand supply rates, A. arenaria produced higher-density vertical tillers (characteristic of higher sand capture efficiency), while A. breviligulata and E. mollis responded with lower-density lateral tiller growth (characteristic of lower sand capture efficiency). Combined, these experiments provide evidence for a species-specific effect on coastal dune shape. Understanding how dominant ecosystem engineers, especially nonnative ones, differ in their interactions with abiotic factors is necessary to better parameterize coastal vulnerability models and inform management practices related to both coastal protection ecosystem services and ecosystem restoration.     

An ecosystem model for estimating potential shellfish culture production in sheltered coastal waters

A generic ecosystem model has been developed for estimating the potential production of shellfish culture and the effect of that cultivation on the pelagic ecosystem in sheltered coastal waters. The model describes the dynamics of a simple food web, nutrient cycling and growth of shellfish. The design of the model is closely tied to the temporal and spatial scales that are important in determining the sustainable production level for a particular embayment. The pelagic ecosystem, mussel energetics, population dynamics and hydrodynamics are coupled to allow fully dynamic predictions of the effect of the shellfish density. When applied to Beatrix Bay, an intensive culture embayment in the Pelorus Sound of New Zealand, the model successfully captured main features of the observed system behaviour. The hydrodynamic regime of the bay controls mussel growth and production. Although high fluxes of water into the bay suppress nutrient and carbon cycling signals in the system, the model simulations demonstrated that the mussel cultivation can have considerable effects on the ecosystem of the bay including food depletion and nutrient cycling. One of the most obvious effects is nutrient enhancement through mussel excretion at low cultivation densities, which promotes primary production particularly during the N-limitation period in summer. The sensitivity analysis identified uncertainty in some parameters and indicated areas for which experimental studies could lead to model improvement. The modelling exercise has established a primary predictive tool for managing mussel aquaculture of a coastal embayment to estimate relationships between the stock level and the growth rate of mussels, and the potentially achievable harvest and stocking density.     

Specificity between Neotropical tree seedlings and their fungal mutualists leads to plant-soil feedback

A growing body of evidence obtained largely from temperate grassland studies suggests that feedbacks occurring between plants and their associated soil biota are important to plant community assemblage. However, few studies have examined the importance of soil organisms in driving plant-soil feedbacks in forested systems. In a tropical forest in central Panama, we examined whether interactions between tree seedlings and their associated arbuscular mycorrhizal fungi (AMF) lead to plant-soil feedback. Specifically, do tropical seedlings modify their own AMF communities in a manner that either favors or inhibits the next cohort of conspecific seedlings (i.e., positive or negative feedback, respectively)? Seedlings of two shade-tolerant tree species (Eugenia nesiotica, Virola surinamensis) and two pioneer tree species (Luehea seemannii, Apeiba aspera) were grown in pots containing identical AMF communities composed of equal amounts of inoculum of six co-occurring AMF species. The different AMF-host combinations were all exposed to two light levels. Under low light (2% PAR), only two of the six AMF species sporulated, and we found that host identity did not influence composition of AMF spore communities. However, relative abundances of three of the four AMF species that produced spores were influenced by host identity when grown under high light (20% PAR). Furthermore, spores of one of the AMF species, Glomus geosporum, were common in soils of Luehea and Eugenia but absent in soils of Apeiba and Virola. We then conducted a reciprocal experiment to test whether AMF communities previously modified by Luehea and Apeiba differentially affected the growth of conspecific and heterospecific seedlings. Luehea seedling growth did not differ between soils containing AMF communities modified by Luehea and Apeiba. However, Apeiba seedlings were significantly larger when grown with Apeiba-modified AMF communities, as compared to Apeiba seedlings grown with Luehea-modifed AMF communities. Our experiments suggest that interactions between tropical trees and their associated AMF are species-specific and that these interactions may shape both tree and AMF communities through plant-soil feedback.