The more the merrier: Multi-species experiments in ecology

A major objective in ecology is to find general patterns, and to establish the rules and underlying mechanisms that generate those patterns. Nevertheless, most of our current insights in ecology are based on case studies of a single or few species, whereas multi-species experimental studies remain rare. We underline the power of the multi-species experimental approach for addressing general ecological questions, e. g. on species environmental responses or on patterns of among-and within-species variation. We present simulations that show that the accuracy of estimates of between-group differences is increased by maximizing the number of species rather than the number of populations or individuals per species. Thus, the more species a multi-species experiment includes, the more powerful it is. In addition, we discuss some inevitable methodological challenges of multi-species experiments. While we acknowledge the value of single-or few-species experiments, we strongly advocate the use of multi-species experiments for addressing ecological questions at a more general level.

Effects of two traits of the ecological state equation on our understanding of species coexistence and ecosystem services

Species coexistence has been a fundamental issue to understand ecosystem functioning since the beginnings of ecology as a science. The search of a reliable and all-encompassing explanation for this issue has become a complex goal with several apparently opposing trends. On the other side, seemingly unconnected with species coexistence, an ecological state equation based on the inverse correlation between an indicator of dispersal that fits gamma distribution and species diversity has been recently developed. This article explores two factors, whose effects are inconspicuous in such an equation at the first sight, that are used to develop an alternative general theoretical background in order to provide a better understanding of species coexistence. Our main outcomes are: (i) the fit of dispersal and diversity values to gamma distribution is an important factor that promotes species coexistence mainly due to the right-skewed character of gamma distribution; (ii) the opposite correlation between species diversity and dispersal implies that any increase of diversity is equivalent to a route of “ecological cooling” whose maximum limit should be constrained by the influence of the third law of thermodynamics; this is in agreement with the well-known asymptotic trend of diversity values in space and time; (iii) there are plausible empirical and theoretical ways to apply physical principles to explain important ecological processes; (iv) the gap between theoretical and empirical ecology in those cases where species diversity is paradoxically high could be narrowed by a wave model of species coexistence based on the concurrency of local equilibrium states. In such a model, competitive exclusion has a limited but indispensable role in harmonious coexistence with functional redundancy. We analyze several literature references as well as ecological and evolutionary examples that support our approach, reinforcing the meaning equivalence between important physical and ecological principles.

Distinct lake level lowstand in Lake Prespa (SE Europe) at the time of the 74 (75) ka Toba eruption

Abstract. The 74 (75) ka Toba eruption in Sumatra, Indonesia, is considered to be one of the largest volcanic events during the Quaternary. Tephra from the Toba eruption has been found in many terrestrial and marine sedimentary deposits, and acidity peaks related to the eruption have been used to synchronize ice core records from Greenland and Antarctica. Seismic profiles and sedimentological data from Lake Prespa on the Balkan Peninsula, SE Europe, indicate a lake level lowstand at 73.6±7.7 ka based on ESR dating of shells. Tephrostratigraphy, radiocarbon dating and tuning of the total organic carbon content with the NGRIP isotope record, corroborate that the lake level lowstand was a short-term event superimposed on the general cooling trend at the end of MIS 5, most likely at the onset of the Greenland Stadial GS- 20. Acknowledging that tectonic events or karst processes could have triggered this lake level lowstand, the chronological correspondence between the lowstand and the Toba eruption is intriguing. Therefore a Toba-driven short-term shift to aridity in the Balkan region, leading to lake level changes and triggering spatial expansion events in one of the lake’s most abundant benthic species, the carino mussel Dreissena presbensis, cannot be excluded.

Stratigraphic analysis of lake level fluctuations in Lake Ohrid: an integration of high resolution hydro-acoustic data and sediment cores

Abstract. Ancient Lake Ohrid is a steep-sided, oligotrophic, karst lake that was tectonically formed most likely within the Pliocene and often referred to as a hotspot of endemic biodiversity. This study aims on tracing significant lake level fluctuations at Lake Ohrid using high-resolution acoustic data in combination with lithological, geochemical, and chronological information from two sediment cores recovered from sub-aquatic terrace levels at ca. 32 and 60m water depth. According to our data, significant lake level fluctuations with prominent lowstands of ca. 60 and 35m below the present water level occurred during Marine Isotope Stage (MIS) 6 and MIS 5, respectively. The effect of these lowstands on biodiversity in most coastal parts of the lake is negligible, due to only small changes in lake surface area, coastline, and habitat. In contrast, biodiversity in shallower areas was more severely affected due to disconnection of today sublacustrine springs from the main water body. Multichannel seismic data from deeper parts of the lake clearly image several clinoform structures stacked on top of each other. These stacked clinoforms indicate significantly lower lake levels prior to MIS 6 and a stepwise rise of water level with intermittent stillstands since its existence as water-filled body, which might have caused enhanced expansion of endemic species within Lake Ohrid.

Genomics and the origin of species

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics

Lower-tropospheric humidity: climatology, trends and the relation to the ITCZ

The tropical region is an area of maximum humidity and serves as the major humidity source of the globe. Among other phenomena, it is governed by the so-called Inter-Tropical Convergence Zone (ITCZ) which is commonly defined by converging low-level winds or enhanced precipitation. Given its importance as a humidity source, we investigate the humidity fields in the tropics in different reanalysis data sets, deduce the climatology and variability and assess the relationship to the ITCZ. Therefore, a new analysis method of the specific humidity distribution is introduced which allows detecting the location of the humidity maximum, the strength and the meridional extent. The results show that the humidity maximum in boreal summer is strongly shifted northward over the warm pool/Asia Monsoon area and the Gulf of Mexico. These shifts go along with a peak in the strength in both areas; however, the extent shrinks over the warm pool/Asia Monsoon area, whereas it is wider over the Gulf of Mexico. In winter, such connections between location, strength and extent are not found. Still, a peak in strength is again identified over the Gulf of Mexico in boreal winter. The variability of the three characteristics is dominated by inter-annual signals in both seasons. The results using ERA-interim data suggest a positive trend in the Gulf of Mexico/Atlantic region from 1979 to 2010, showing an increased northward shift in the recent years. Although the trend is only weakly confirmed by the results using MERRA reanalysis data, it is in phase with a trend in hurricane activity�a possible hint of the importance of the new method on hurricanes. Furthermore, the position of the maximum humidity coincides with one of the ITCZ in most areas. One exception is the western and central Pacific, where the area is dominated by the double ITCZ in boreal winter. Nevertheless, the new method enables us to gain more insight into the humidity distribution, its variability and the relationship to ITCZ characteristics.

Impacts of drought stress on water relations and carbon assimilation in grassland species at different altitudes

Grassland is an important ecosystem type which is not only used agriculturally, but also covers sites which cannot be used for other purposes, e.g. in very steep locations or above timberlines. Prolonged summer droughts in the mid-term future, as are predicted for Central Europe, are expected to have a major impact on such ecosystems. To address this topic, rainfall exclusion via shelters was performed on three grassland sites at different altitudes (393, 982 and 1978 m above sea level) in Switzerland. Diurnal drought treatment effects were studied at each study site on a completely sunny day towards the end of an 8–10 week shelter period. Ecophysiological parameters including gas exchange (An, gs and intrinsic WUE) and chlorophyll a fluorescence (Fv/Fm, ΦPSII and NPQ) were considered for several species. The lowland and the Alpine field site were more strongly affected by soil drought than the pre-Alpine site. At all sites, grasses showed different patterns of reductions in stomatal conductance under soil drought compared to legumes and forbs. In addition, grasses were significantly more affected by reductions in assimilation rates at all sites. Time courses of reductions in assimilation rates relative to controls differed between species at the Alpine site, as some species showed reduced assimilation rates at this site in the early morning. Thus, similar rainfall exclusion treatments can trigger different reactions in various species at different sites, which might not become obvious during mere midday measurements. Overall, results suggest strong impacts of prolonged summer drought on grassland net photosynthesis especially at the Alpine site and, within sites, for grasses

Stable two-element control of dTph1 transposition in mutator strains of Petunia by an inactive ACT1 introgression from a wild species

The high copy dTph1 transposon system of Petunia (Solanaceae) is one of the most powerful insertion mutagens in plants, but its activity cannot be controlled in the commonly used mutator strains. We analysed the regulation of dTph1 activity by QTL analysis in recombinant inbred lines of the mutator strain W138 and a wild species (P. integrifolia spp. inflata). Two genetic factors were identified that control dTph1 transposition. One corresponded to the ACT1 locus on chromosome I. A second, previously undescribed locus ACT2 mapped on chromosome V. As a 6-cM introgression in W138, the P. i. inflata act1(S6) allele behaved as a single recessive locus that fully eliminated transposition of all dTph1 elements in all stages of plant development and in a heritable fashion. Weak dTph1 activity was restored in act1(S6)/ACT2(S6) double introgression lines, indicating that the P. i. inflata allele at ACT2 conferred a low level of transposition. Thus, the act1(S6) allele is useful for simple and predictable control of transposition of the entire dTph1 family when introgressed into an ultra-high copy W138 mutator strain. We demonstrate the use of the ACT1(W138)/act1(S6) allele pair in a two-element dTph1 transposition system by producing 10 000 unique and fixed dTph1 insertions in a population of 1250 co-isogenic lines. This Petunia system produces the highest per plant insertion number of any known two-element system, providing a powerful and logistically simple tool for transposon mutagenesis of qualitative as well as quantitative traits.

Plant-plant interactions, environmental gradients and plant diversity: A global synthesis of community-level studies

Previous syntheses on the effects of environmental conditions on the outcome of plant-plant interactions summarize results from pairwise studies. However, the upscaling to the community-level of such studies is problematic because of the existence of multiple species assemblages and species-specific responses to both the environmental conditions and the presence of neighbors. We conducted the first global synthesis of community-level studies from harsh environments, which included data from 71 alpine and 137 dryland communities to: (i) test how important are facilitative interactions as a driver of community structure, (ii) evaluate whether we can predict the frequency of positive plant-plant interactions across differing environmental conditions and habitats, and (iii) assess whether thresholds in the response of plant-plant interactions to environmental gradients exists between ``moderate'' and ``extreme'' environments. We also used those community-level studies performed across gradients of at least three points to evaluate how the average environmental conditions, the length of the gradient studied, and the number of points sampled across such gradient affect the form and strength of the facilitation-environmental conditions relationship. Over 25% of the species present were more spatially associated to nurse plants than expected by chance in both alpine and chyland areas, illustrating the high importance of positive plant-plant interactions for the maintenance of plant diversity in these environments. Facilitative interactions were more frequent, and more related to environmental conditions, in alpine than in dryland areas, perhaps because drylands are generally characterized by a larger variety of environmental stress factors and plant functional traits. The frequency of facilitative interactions in alpine communities peaked at 1000 mm of annual rainfall, and globally decreased with elevation. The frequency of positive interactions in dtyland communities decreased globally with water scarcity or temperature annual range. Positive facilitation-drought stress relationships are more likely in shorter regional gradients, but these relationships are obscured in regions with a greater species turnover or with complex environmental gradients. By showing the different climatic drivers and behaviors of plant-plant interactions in dryland and alpine areas, our results will improve predictions regarding the effect of facilitation on the assembly of plant communities and their response to changes in environmental conditions.

UNDERSTANDING TRENDS ASSOCIATED WITH CLOUDS IN IRRADIATED EXOPLANETS

Unlike previously explored relationships between the properties of hot Jovian atmospheres, the geometric albedo and the incident stellar flux do not exhibit a clear correlation, as revealed by our re-analysis of Q0-Q14 Kepler data. If the albedo is primarily associated with the presence of clouds in these irradiated atmospheres, a holistic modeling approach needs to relate the following properties: the strength of stellar irradiation (and hence the strength and depth of atmospheric circulation), the geometric albedo (which controls both the fraction of starlight absorbed and the pressure level at which it is predominantly absorbed), and the properties of the embedded cloud particles (which determine the albedo). The anticipated diversity in cloud properties renders any correlation between the geometric albedo and the stellar flux weak and characterized by considerable scatter. In the limit of vertically uniform populations of scatterers and absorbers, we use an analytical model and scaling relations to relate the temperature-pressure profile of an irradiated atmosphere and the photon deposition layer and to estimate whether a cloud particle will be lofted by atmospheric circulation. We derive an analytical formula for computing the albedo spectrum in terms of the cloud properties, which we compare to the measured albedo spectrum of HD 189733b by Evans et al. Furthermore, we show that whether an optical phase curve is flat or sinusoidal depends on whether the particles are small or large as defined by the Knudsen number. This may be an explanation for why Kepler-7b exhibits evidence for the longitudinal variation in abundance of condensates, while Kepler-12b shows no evidence for the presence of condensates despite the incident stellar flux being similar for both exoplanets. We include an "observer's cookbook" for deciphering various scenarios associated with the optical phase curve, the peak offset of the infrared phase curve, and the geometric albedo.

Tree size and fecundity influence ballistic seed dispersal of two dominant mast-fruiting species in a tropical rain forest

Seed production, seed dispersal, and seedling recruitment are integral to forest dynamics, especially in masting species. Often these are studied separately, yet scarcely ever for species with ballistic dispersal even though this mode of dispersal is common in legume trees of tropical African rain forests. Here, we studied two dominant main-canopy tree species, Microberlinia bisulcata and Tetraberlinia bifoliolata (Caesalpinioideae), in 25 ha of primary rain forest at Korup, Cameroon, during two successive masting events (2007/2010). In the vicinity of c. 100 and 130 trees of each species, 476/580 traps caught dispersed seeds and beneath their crowns c. 57,000 pod valves per species were inspected to estimate tree-level fecundity. Seed production of trees increased non-linearly and asymptotically with increasing stem diameters. It was unequal within the two species’ populations, and differed strongly between years to foster both spatial and temporal patchiness in seed rain. The M. bisulcata trees could begin seeding at 42–44 cm diameter: at a much larger size than could T. bifoliolata (25 cm). Nevertheless, per capita life-time reproductive capacity was c. five times greater in M. bisulcata than T. bifoliolata owing to former’s larger adult stature, lower mortality rate (despite a shorter life-time) and smaller seed mass. The two species displayed strong differences in their dispersal capabilities. Inverse modelling (IM) revealed that dispersal of M. bisulcata was best described by a lognormal kernel. Most seeds landed at 10–15 m from stems, with 1% of them going beyond 80 m (<100 m). The direct estimates of fecundity significantly improved the models fitted. The lognormal also described well the seedling recruitment distribution of this species in 121 ground plots. By contrast, the lower intensity of masting and more limited dispersal of the heavier-seeded T. bifoliolata prevented reliable IM. For this species, seed density as function of distance to traps suggested a maximum dispersal distance of 40–50 m, and a correspondingly more aggregated seedling recruitment pattern ensued than for M. bisulcata. From this integrated field study, we conclude that the reproductive traits of M. bisulcata give it a considerable advantage over T. bifoliolata by better dispersing more seeds per capita to reach more suitable establishment sites, and combined with other key traits they explain its local dominance in the forest. Understanding the linkages between size at onset of maturity, individual fecundity, and dispersal capability can better inform the life-history strategies, and hence management, of co-occurring tree species in tropical forests.

Little evidence for release from herbivores as a driver of plant invasiveness from a multi-species herbivore-removal experiment

Enemy release is frequently posed as a main driver of invasiveness of alien species. However, an experimental multi-species test examining performance and herbivory of invasive alien, non-invasive alien and native plant species in the presence and absence of natural enemies is lacking. In a common garden experiment in Switzerland, we manipulated exposure of seven alien invasive, eight alien non-invasive and fourteen native species from six taxonomic groups to natural enemies (invertebrate herbivores), by applying a pesticide treatment under two different nutrient levels. We assessed biomass production, herbivore damage and the major herbivore taxa on plants. Across all species, plants gained significantly greater biomass under pesticide treatment. However, invasive, non-invasive and native species did not differ in their biomass response to pesticide treatment at either nutrient level. The proportion of leaves damaged on invasive species was significantly lower compared to native species, but not when compared to non-invasive species. However, the difference was lost when plant size was accounted for. There were no differences between invasive, non-invasive and native species in herbivore abundance. Our study offers little support for invertebrate herbivore release as a driver of plant invasiveness, but suggests that future enemy release studies should account for differences in plant size among species.

Delayed HIV diagnosis and initiation of antiretroviral therapy: inequalities by educational level, COHERE in EuroCoord.

OBJECTIVES In Europe and elsewhere, health inequalities among HIV-positive individuals are of concern. We investigated late HIV diagnosis and late initiation of combination antiretroviral therapy (cART) by educational level, a proxy of socioeconomic position. DESIGN AND METHODS We used data from nine HIV cohorts within COHERE in Austria, France, Greece, Italy, Spain and Switzerland, collecting data on level of education in categories of the UNESCO/International Standard Classification of Education standard classification: non-completed basic, basic, secondary and tertiary education. We included individuals diagnosed with HIV between 1996 and 2011, aged at least 16 years, with known educational level and at least one CD4 cell count within 6 months of HIV diagnosis. We examined trends by education level in presentation with advanced HIV disease (AHD) (CD4 <200 cells/μl or AIDS within 6 months) using logistic regression, and distribution of CD4 cell count at cART initiation overall and among presenters without AHD using median regression. RESULTS Among 15 414 individuals, 52, 45,37, and 31% with uncompleted basic, basic, secondary and tertiary education, respectively, presented with AHD (P trend <0.001). Compared to patients with tertiary education, adjusted odds ratios of AHD were 1.72 (95% confidence interval 1.48-2.00) for uncompleted basic, 1.39 (1.24-1.56) for basic and 1.20 (1.08-1.34) for secondary education (P < 0.001). In unadjusted and adjusted analyses, median CD4 cell count at cART initiation was lower with poorer educational level. CONCLUSIONS Socioeconomic inequalities in delayed HIV diagnosis and initiation of cART are present in European countries with universal healthcare systems and individuals with lower educational level do not equally benefit from timely cART initiation.

Association between trends in clinical variables and outcome in intensive care patients with faecal peritonitis: analysis of the GenOSept cohort.

INTRODUCTION Patients admitted to intensive care following surgery for faecal peritonitis present particular challenges in terms of clinical management and risk assessment. Collaborating surgical and intensive care teams need shared perspectives on prognosis. We aimed to determine the relationship between dynamic assessment of trends in selected variables and outcomes. METHODS We analysed trends in physiological and laboratory variables during the first week of intensive care unit (ICU) stay in 977 patients at 102 centres across 16 European countries. The primary outcome was 6-month mortality. Secondary endpoints were ICU, hospital and 28-day mortality. For each trend, Cox proportional hazards (PH) regression analyses, adjusted for age and sex, were performed for each endpoint. RESULTS Trends over the first 7 days of the ICU stay independently associated with 6-month mortality were worsening thrombocytopaenia (mortality: hazard ratio (HR) = 1.02; 95% confidence interval (CI), 1.01 to 1.03; P <0.001) and renal function (total daily urine output: HR =1.02; 95% CI, 1.01 to 1.03; P <0.001; Sequential Organ Failure Assessment (SOFA) renal subscore: HR = 0.87; 95% CI, 0.75 to 0.99; P = 0.047), maximum bilirubin level (HR = 0.99; 95% CI, 0.99 to 0.99; P = 0.02) and Glasgow Coma Scale (GCS) SOFA subscore (HR = 0.81; 95% CI, 0.68 to 0.98; P = 0.028). Changes in renal function (total daily urine output and renal component of the SOFA score), GCS component of the SOFA score, total SOFA score and worsening thrombocytopaenia were also independently associated with secondary outcomes (ICU, hospital and 28-day mortality). We detected the same pattern when we analysed trends on days 2, 3 and 5. Dynamic trends in all other measured laboratory and physiological variables, and in radiological findings, changes inrespiratory support, renal replacement therapy and inotrope and/or vasopressor requirements failed to be retained as independently associated with outcome in multivariate analysis. CONCLUSIONS Only deterioration in renal function, thrombocytopaenia and SOFA score over the first 2, 3, 5 and 7 days of the ICU stay were consistently associated with mortality at all endpoints. These findings may help to inform clinical decision making in patients with this common cause of critical illness.