Simulating the evolution of a clonal trait in plants with sexual and vegetative reproduction

Aims Phenotypic optimality models neglect genetics. However, especially when heterozygous genotypes ire fittest, evolving allele, genotype and phenotype frequencies may not correspond to predicted optima. This was not previously addressed for organisms with complex life histories. Methods Therefore, we modelled the evolution of a fitness-relevant trait of clonal plants, stolon internode length. We explored the likely case of air asymmetric unimodal fitness profile with three model types. In constant selection models (CSMs), which are gametic, but not spatially explicit, evolving allele frequencies in the one-locus and five-loci cases did not correspond to optimum stolon internode length predicted by the spatially explicit, but not gametic, phenotypic model. This deviation was due to the asymmetry of the fitness profile. Gametic, spatially explicit individual-based (SEIB) modeling allowed us relaxing the CSM assumptions of constant selection with exclusively sexual reproduction. Important findings For entirely vegetative or sexual reproduction, predictions. of the gametic SEIB model were close to the ones of spatially explicit CSMs gametic phenotypic models, hut for mixed modes of reproduction they appoximated those of gametic, not spatially explicit CSMs. Thus, in contrast to gametic SEIB models, phenotypic models and, especially for few loci, also CSMs can be very misleading. We conclude that the evolution of trails governed by few quantitative trait loci appears hardly predictable by simple models, that genetic algorithms aiming at technical optimization may actually, miss the optimum and that selection may lead to loci with smaller effects, in derived compared with ancestral lines.

Herbivores limit the population size of big-leaf mahogany trees in an Amazonian forest

The Janzen–Connell hypothesis proposes that specialized herbivores maintain high numbers of tree species in tropical forests by restricting adult recruitment so that host populations remain at low densities. We tested this prediction for the large timber tree species, Swietenia macrophylla, whose seeds and seedlings are preyed upon by small mammals and a host-specific moth caterpillar Steniscadia poliophaea, respectively. At a primary forest site, experimental seed additions to gaps – canopy-disturbed areas that enhance seedling growth into saplings – over three years revealed lower survival and seedling recruitment closer to conspecific trees and in higher basal area neighborhoods, as well as reduced subsequent seedling survival and height growth. When we included these Janzen–Connell effects in a spatially explicit individual-based population model, the caterpillar's impact was critical to limiting Swietenia's adult tree density, with a > 10-fold reduction estimated at 300 years. Our research demonstrates the crucial but oft-ignored linkage between Janzen–Connell effects on offspring and population-level consequences for a long-lived, potentially dominant tree species.

The Effects of Density, Spatial Pattern, and Competitive Symmetry on Size Variation in Simulated Plant Populations

Patterns of size inequality in crowded plant populations are often taken to be indicative of the degree of size asymmetry of competition, but recent research suggests that some of the patterns attributed to size‐asymmetric competition could be due to spatial structure. To investigate the theoretical relationships between plant density, spatial pattern, and competitive size asymmetry in determining size variation in crowded plant populations, we developed a spatially explicit, individual‐based plant competition model based on overlapping zones of influence. The zone of influence of each plant is modeled as a circle, growing in two dimensions, and is allometrically related to plant biomass. The area of the circle represents resources potentially available to the plant, and plants compete for resources in areas in which they overlap. The size asymmetry of competition is reflected in the rules for dividing up the overlapping areas. Theoretical plant populations were grown in random and in perfectly uniform spatial patterns at four densities under size‐asymmetric and size‐symmetric competition. Both spatial pattern and size asymmetry contributed to size variation, but their relative importance varied greatly over density and over time. Early in stand development, spatial pattern was more important than the symmetry of competition in determining the degree of size variation within the population, but after plants grew and competition intensified, the size asymmetry of competition became a much more important source of size variation. Size variability was slightly higher at higher densities when competition was symmetric and plants were distributed nonuniformly in space. In a uniform spatial pattern, size variation increased with density only when competition was size asymmetric. Our results suggest that when competition is size asymmetric and intense, it will be more important in generating size variation than is local variation in density. Our results and the available data are consistent with the hypothesis that high levels of size inequality commonly observed within crowded plant populations are largely due to size‐asymmetric competition, not to variation in local density.

Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

Simulating Carbon Stocks and Fluxes of an African Tropical Montane Forest with an Individual-Based Forest Model

Tropical forests are carbon-dense and highly productive ecosystems. Consequently, they play an important role in the global carbon cycle. In the present study we used an individual-based forest model (FORMIND) to analyze the carbon balances of a tropical forest. The main processes of this model are tree growth, mortality, regeneration, and competition. Model parameters were calibrated using forest inventory data from a tropical forest at Mt. Kilimanjaro. The simulation results showed that the model successfully reproduces important characteristics of tropical forests (aboveground biomass, stem size distribution and leaf area index). The estimated aboveground biomass (385 t/ha) is comparable to biomass values in the Amazon and other tropical forests in Africa. The simulated forest reveals a gross primary production of 24 tcha-1yr-1. Modeling above- and belowground carbon stocks, we analyzed the carbon balance of the investigated tropical forest. The simulated carbon balance of this old-growth forest is zero on average. This study provides an example of how forest models can be used in combination with forest inventory data to investigate forest structure and local carbon balances.

Development of a Novel Rabies Simulation Model for Application in a Non-endemic Environment.

Domestic dog rabies is an endemic disease in large parts of the developing world and also epidemic in previously free regions. For example, it continues to spread in eastern Indonesia and currently threatens adjacent rabies-free regions with high densities of free-roaming dogs, including remote northern Australia. Mathematical and simulation disease models are useful tools to provide insights on the most effective control strategies and to inform policy decisions. Existing rabies models typically focus on long-term control programs in endemic countries. However, simulation models describing the dog rabies incursion scenario in regions where rabies is still exotic are lacking. We here describe such a stochastic, spatially explicit rabies simulation model that is based on individual dog information collected in two remote regions in northern Australia. Illustrative simulations produced plausible results with epidemic characteristics expected for rabies outbreaks in disease free regions (mean R0 1.7, epidemic peak 97 days post-incursion, vaccination as the most effective response strategy). Systematic sensitivity analysis identified that model outcomes were most sensitive to seven of the 30 model parameters tested. This model is suitable for exploring rabies spread and control before an incursion in populations of largely free-roaming dogs that live close together with their owners. It can be used for ad-hoc contingency or response planning prior to and shortly after incursion of dog rabies in previously free regions. One challenge that remains is model parameterisation, particularly how dogs' roaming and contacts and biting behaviours change following a rabies incursion in a previously rabies free population.

Colonization history of the Swiss Rhine basin by the bullhead (Cottus gobio): inference under a Bayesian spatially explicit framework

The present distribution of freshwater fish in the Alpine region has been strongly affected by colonization events occurring after the last glacial maximum (LGM), some 20,000 years ago. We use here a spatially explicit simulation framework to model and better understand their colonization dynamics in the Swiss Rhine basin. This approach is applied to the European bullhead (Cottus gobio), which is an ideal model organism to study fish past demographic processes since it has not been managed by humans. The molecular diversity of eight sampled populations is simulated and compared to observed data at six microsatellite loci under an approximate Bayesian computation framework to estimate the parameters of the colonization process. Our demographic estimates fit well with current knowledge about the biology of this species, but they suggest that the Swiss Rhine basin was colonized very recently, after the Younger Dryas some 6600 years ago. We discuss the implication of this result, as well as the strengths and limits of the spatially explicit approach coupled to the approximate Bayesian computation framework.

Eco-evolutionary Dynamics of Individual-Based Food Webs

The past decade has seen the rise of high resolution datasets. One of the main surprises of analysing such data has been the discovery of a large genetic, phenotypic and behavioural variation and heterogeneous metabolic rates among individuals within natural populations. A parallel discovery from theory and experiments has shown a strong temporal convergence between evolutionary and ecological dynamics, but a general framework to analyse from individual-level processes the convergence between ecological and evolutionary dynamics and its implications for patterns of biodiversity in food webs has been particularly lacking. Here, as a first approximation to take into account intraspecific variability and the convergence between the ecological and evolutionary dynamics in large food webs, we develop a model from population genomics and microevolutionary processes that uses sexual reproduction, genetic-distance-based speciation and trophic interactions. We confront the model with the prey consumption per individual predator, species-level connectance and prey–predator diversity in several environmental situations using a large food web with approximately 25,000 sampled prey and predator individuals. We show higher than expected diversity of abundant species in heterogeneous environmental conditions and strong deviations from the observed distribution of individual prey consumption (i.e. individual connectivity per predator) in all the environmental conditions. The observed large variance in individual prey consumption regardless of the environmental variability collapsed species-level connectance after small increases in sampling effort. These results suggest (1) intraspecific variance in prey–predator interactions has a strong effect on the macroscopic properties of food webs and (2) intraspecific variance is a potential driver regulating the speed of the convergence between ecological and evolutionary dynamics in species-rich food webs. These results also suggest that genetic–ecological drift driven by sexual reproduction, equal feeding rate among predator individuals, mutations and genetic-distance-based speciation can be used as a neutral food web dynamics test to detect the ecological and microevolutionary processes underlying the observed patterns of individual and species-based food webs at local and macroecological scales.

Intermittent spike-wave dynamics in a heterogeneous, spatially extended neural mass model

Generalised epileptic seizures are frequently accompanied by sudden, reversible transitions from low amplitude, irregular background activity to high amplitude, regular spike-wave discharges (SWD) in the EEG. The underlying mechanisms responsible for SWD generation and for the apparently spontaneous transitions to SWD and back again are still not fully understood. Specifically, the role of spatial cortico-cortical interactions in ictogenesis is not well studied. We present a macroscopic, neural mass model of a cortical column which includes two distinct time scales of inhibition. This model can produce both an oscillatory background and a pathological SWD rhythm. We demonstrate that coupling two of these cortical columns can lead to a bistability between out-of-phase, low amplitude background dynamics and in-phase, high amplitude SWD activity. Stimuli can cause state-dependent transitions from background into SWD. In an extended local area of cortex, spatial heterogeneities in a model parameter can lead to spontaneous reversible transitions from a desynchronised background to synchronous SWD due to intermittency. The deterministic model is therefore capable of producing absence seizure-like events without any time dependent adjustment of model parameters. The emergence of such mechanisms due to spatial coupling demonstrates the importance of spatial interactions in modelling ictal dynamics, and in the study of ictogenesis.

Transmission of Chlamydia trachomatis through sexual partnerships: a comparison between three individual-based models and empirical data

Chlamydia trachomatis is the most common bacterial sexually transmitted infection (STI) in many developed countries. The highest prevalence rates are found among young adults who have frequent partner change rates. Three published individual-based models have incorporated a detailed description of age-specific sexual behaviour in order to quantify the transmission of C. trachomatis in the population and to assess the impact of screening interventions. Owing to varying assumptions about sexual partnership formation and dissolution and the great uncertainty about critical parameters, such models show conflicting results about the impact of preventive interventions. Here, we perform a detailed evaluation of these models by comparing the partnership formation and dissolution dynamics with data from Natsal 2000, a population-based probability sample survey of sexual attitudes and lifestyles in Britain. The data also allow us to describe the dispersion of C. trachomatis infections as a function of sexual behaviour, using the Gini coefficient. We suggest that the Gini coefficient is a useful measure for calibrating infectious disease models that include risk structure and highlight the need to estimate this measure for other STIs.