The empirical Bayes approach as a tool to identify non-random species associations

The Co-occurrence Matrix and the Correlation Network of Phytophagous Insects Are Driven by Abiotic and Biotic Variables: the Case of Canola

Co-occurrence a correlation profiles are driven by different factors (exogenous and endogenous) and drawing a profile of association between species based on co-occurrence, without assessing how these species vary in terms of ecological niche can lead to wrong conclusions. The objective was to determine the co-occurrence and correlation patterns of phytophagous insects in canola crop and to evaluate how these patterns varied according to the crop stage (phenology-biotic) and sowing times (agricultural practice-abiotic). We found that the patterns of co-occurrence and correlation between species were reflections of population variations due to the phenology and sowing times of canola. Variations in the multi-species abundance matrix were influenced by mean air temperature and accumulated rainfall. The main species associated with canola in southern Brazil, in terms of abundance, were P. xylostella, D. speciosa, and N. viridula. These species were mostly negatively associated. When evaluating their population variations, we found that they explore different temporal niches, whether in terms of phenology or sowing times. Finally, we demonstrate empirically that despite being important, association patterns based on co-occurrence and correlation should be interpreted in light of the understanding of patterns of niche exploitation and temporal variation of species.

Probabilistic approaches for investigating species co-occurrence from presence-absence maps

Background

In this research, we propose probabilistic approaches to identify pairwise patterns of species co-occurrence by using presence-absence maps only. In particular, the two-by-two contingency table constructed from a presence-absence map of two species would be sufficient to compute the test statistics and perform the statistical tests proposed in this article. Some previous studies have investigated species co-occurrence through incidence data of different survey sites. We focus on using presence-absence maps for a specific study plot instead. The proposed methods are assessed by a thorough simulation study.

Methods

A Chi-squared test is used to determine whether the distributions of two species are independent. If the null hypothesis of independence is rejected, the Chi-squared method can not distinguish positive or negative association between two species. We propose six different approaches based on either the binomial or Poisson distribution to obtain p-values for testing the positive (or negative) association between two species. When we test to investigate a positive (or negative) association, if the p-value is below the predetermined level of significance, then we have enough evidence to support that the two species are positively (or negatively) associated.

Results

A simulation study is conducted to demonstrate the type-I errors and the testing powers of our approaches. The probabilistic approach proposed by Veech (2013) is served as a benchmark for comparison. The results show that the type-I error of the Chi-squared test is close to the significance level when the presence rate is between 40% and 80%. For extremely low or high presence rate data, one of our approaches outperforms Veech (2013)'s in terms of the testing power and type-I error rate. The proposed methods are applied to a tree data of Barro Colorado Island in Panama and a tree data of Lansing Woods in USA. Both positive and negative associations are found among some species in these two real data.

A needle in the haystack? Applying species co-occurrence frameworks with fish assemblage data to identify species associations and sharpen ecological hypotheses

Different species can associate or interact in many ways, and methods exist for inferring associations and underlying mechanisms from incidence data (e.g., co-occurrence frameworks). These methods have received criticism despite their recent resurgence in the literature. However, co-occurrence frameworks for identifying nonrandomly associated species pairs (e.g., aggregated or segregated pairs) have value as heuristic tools for sharpening hypotheses concerning fish ecology. This paper provides a case study examining species co-occurrence across 33 stream fish assemblages in southeastern Oklahoma, USA, which were sampled twice (1974 and 2014). This study sought to determine (a) which species were nonrandomly associated, (b) what processes might have driven these associations and (c) how consistent patterns were across time. Associations among most pairs of species (24 species, 276 unique pairs) were not significantly different from random (>80%). Among all significant, nonrandomly associated species pairs (54 unique pairs), 78% (42 pairs) were aggregated and 22% (12 pairs) segregated. Most of these (28 pairs, 52%) were hypothesized to be driven by nonbiotic mechanisms: habitat filtering (20 pairs, 37%), dispersal limitation (two pairs, 0.4%) or both (six pairs, 11%). The remaining 26 nonrandomly associated pairs (48%) had no detectable signal of spatial or environmental factors involved with the association, therefore the potential for biotic interaction was not refuted. Only five species pairs were consistently associated across both sampling periods: stonerollers Campostoma spp. and orangebelly darter Etheostoma radiosum; red shiner Cyprinella lutrensis and bullhead minnow Pimephales vigilax; bluegill sunfish Lepomis macrochirus and redear sunfish Lepomis microlophus; redfin shiner Lythrurus umbratilis and bluntnose minnow Pimephales notatus; and bigeye shiner Notropis boops and golden shiner Notemigonus crysoleucas. Frameworks for identifying nonrandomly associated species pairs can provide insight into broader mechanisms of species assembly and point to potentially interesting species interactions (out of many possible pairs). However, this approach is best applied as a tool for sharpening hypotheses to be investigated further. Rather than a weakness, the heuristic nature is the strength of such methods, and can help guide biologists toward better questions by employing relatively cheap diversity survey data, which are often already in hand, to reduce complex interaction networks down to their nonstochastic parts which warrant further investigation.

Antigonon leptopus invasion is associated with plant community disassembly in a Caribbean island ecosystem

Invasions by non-native plant species are widely recognized as a major driver of biodiversity loss. Globally, (sub-)tropical islands form important components of biodiversity hotspots, while being particularly susceptible to invasions by plants in general and vines in particular. We studied the impact of the invasive vine A. leptopus on the diversity and structure of recipient plant communities on the northern Caribbean island St. Eustatius. We used a paired-plot design to study differences in species richness, evenness and community structure under A. leptopus-invaded and uninvaded conditions. Community structure was studied through species co-occurrence patterns. We found that in plots invaded by A. leptopus, species richness was 40–50% lower, and these plots also exhibited lower evenness. The magnitude of these negative impacts increased with increasing cover of A. leptopus. Invaded plots also showed higher degrees of homogeneity in species composition. Species co-occurrence patterns indicated that plant communities in uninvaded plots were characterized by segregation, whereas recipient plant communities in invaded plots exhibited random co-occurrence patterns. These observations suggest that invasion of A. leptopus is not only associated with reduced species richness and evenness of recipient communities in invaded sites, but also with a community disassembly process that may reduce diversity between sites. Given that A. leptopus is a successful invader of (sub-)tropical islands around the globe, these impacts on plant community structure highlight that this invasive species could be a particular conservation concern for these systems.

Patterns of association of native and exotic boring polychaetes on the southeastern Pacific coast of Chile: the combined importance of negative, positive and random interactions

BACKGROUND

Studies of biological invasions focus on negative interactions between exotic and native biotas, emphasizing niche overlap between species and competitive exclusion. However, the effects of positive interactions and coexistence are poorly known. In this study we evaluate the importance of positive, negative, or random species associations in explaining the coexistence of native and exotic boring polychaetes inhabiting invertebrate hosts, on the southeastern Pacific coast of Chile. We assess three hypotheses to explain the observed patterns: positive species interactions, weak competitive interactions, and competitive intransitivity.

METHODOLOGY

To assess the potential effect of competition between native and exotic polychaetes we analyzed patterns of co-occurrence of species pairs in northern and southern regions, using the metric of the probabilistic model. Since biotic interactions can affect not only native species, we also evaluated correlations between native and exotic polychaete abundance, using reduced major axis regression linear models. To assess the transitivity of competitive hierarchies we used metrics and analytical methods based on abundance matrices to estimate species competition and patch transition matrices.

RESULTS

On average 50% of the species pairs presented significant weak negative associations, all associated with the exotic species Polydora rickettsi; the remaining 50% had random associations, and none showed positive associations. However, the relationship of abundance between native and exotic boring polychates supports a tendency towards coexistence. At local and regional scales, we observed the presence of a transitive network competition structure, where the exotic boring polychaete, P. rickettsi was generally the dominant species.

CONCLUSIONS

Our results support that native and exotic boring polychaete species coexist through weak competitive interactions. Nevertheless, the large number of random interactions observed indicates that species coexistence can be accounted for by stochastic processes, as proposed by neutral theory. Coexistence may be a frequent result of interactions between native and exotic species, although less apparent than competitive exclusion. Thus, the probabilistic point-of-view used here provides a statistical tool for evaluating coexistence as a result of exotic and native species' interactions, an idea which has been proposed in invasion ecology, but largely lacks empirical support and methodologies for detecting underlying mechanisms. Finally, we found evidence that P. rickettsi is a successful invader by being competitively dominant, but not excluding other species.

Native and Invasive Small Mammals in Urban Habitats along the Commercial Axis Connecting Benin and Niger, West Africa

Based on compiled small mammal trapping data collected over 12 years from Benin and Niger (3701 individual records from 66 sampling sites), located in mainland Africa, we here describe the small mammal community assemblage in urban habitats along the commercial axis connecting the two countries, from the seaport of Cotonou to the Sahelian hinterland, with a particular focus on invasive species. In doing so, we document extant species distributions, which highlight the risks of continuing the range expansion of three synanthropic invasive rodent species, namely black rats (Rattus rattus), brown rats (R. norvegicus), and house mice (Mus musculus). Using various diversity estimates and community ecology approaches, we detect a latitudinal gradient of species richness that significantly decreased Northward. We show that shrews (Crocidura) represent a very important component of micro-mammal fauna in West African towns and villages, especially at lower latitudes. We also demonstrate that invasive and native synanthropic rodents do not distribute randomly in West Africa, which suggests that invasive species dynamics and history differ markedly, and that they involve gradual, as well as human-mediated, long distance dispersal. Patterns of segregation are also observed between native Mastomys natalensis and invasive rats R. rattus and R. norvegicus, suggesting potential native-to-invasive species turn over. Consequences of such processes, especially in terms of public health, are discussed.

Null model analyses of temporal patterns of bird assemblages and their foraging guilds revealed the predominance of positive and random associations

Patterns of species associations have been commonly used to infer interactions among species. If species positively co-occur, they may form predominantly neutral assemblages, and such patterns suggest a relatively weak role for compensatory dynamics. The main objective of this study was to test this prediction on temporal samples of bird assemblages (n = 19, 10-57 years) by the presence/absence and quantitative null models on assemblage and guild levels. These null model outcomes were further analyzed to evaluate the effects of various data set characteristics on the outcomes of the null models. The analysis of two binary null models in combination with three association indices revealed 20% with significant aggregations, 61% with random associations, and only 19% with significant segregations (n = 95 simulations). The results of the quantitative null model simulations detected more none-random associations: 61% aggregations, 6% random associations, and 33% segregations (n = 114 simulations). Similarly, quantitative analyses on guild levels showed 58% aggregations, 20% segregations, and 22% random associations (n = 450 simulations). Bayesian GLMs detected that the outcomes of the binary and quantitative null models applied to the assemblage analyses were significantly related to census plot size, whereas the outcomes of the quantitative analyses were also related to the mean population densities of species in the data matrices. In guild-level analyses, only 9% of the GLMs showed a significant influence of matrix properties (plot size, matrix size, species richness, and mean species population densities) on the null model outcomes. The results did not show the prevalence of negative associations that would have supported compensatory dynamics. Instead, we assume that a similar response of the majority of species to climate-driven and stochastic factors may be responsible for the revealed predominance of positive associations.

Oribatid mites reveal that competition for resources and trophic structure combine to regulate the assembly of diverse soil animal communities

The role of niche partitioning in structuring species-rich soil animal communities has been debated for decades and generated the "enigma of soil animal diversity." More recently, resource-based niche partitioning has been hypothesized to play a very limited role in the assembly of soil animal communities. To test this hypothesis, we applied a novel combination of stable isotopes and null models of species co-occurrence to quantify the extent of resource niche partitioning on a diverse oribatid mite community sampled from mature oak woodland.We asked whether species aggregate or segregate spatially and how these patterns correlated with the abundance of estimated trophic guilds. We also estimated the effects of environmental variables on community structure.All measured environmental variables accounted for 12% of variance in community structure, including 8% of pure spatial structure unrelated to measured environmental factors and 2% of pure environmental variance unrelated to spatial variation. Co-occurrence analysis revealed 10 pairs of species that aggregated and six pairs of species that were spatially segregated. Values of δ15N indicated that five out of the 10 pairs of aggregated species occupied the same trophic guild, while values of δ13C indicated that species in these five pairs consumed resources of different quality, supporting a significant role of resource-based niche partitioning. Also, one of the five pairs of segregated species occupied the same trophic guild but had overlapping δ13C values suggesting that these species do not co-occur locally and thus minimize competition for shared resources.Partitioning of resources plays an underestimated role in soil microarthropod communities and different local communities consisted of the same trophic guilds with species identity changing from place to place. The sum of resource partitioning, multi-trophic interactions, and microscale environmental variability in the environment is a viable solution to the enigma of soil animal diversity.

Co-occurrence patterns of rodents at multiple spatial scales: competitive release of generalists following habitat loss?

Theory predicts that habitat generalist species are excluded by specialist species in optimal habitat for specialists, and empirical data commonly show a shift from specialist- to generalist-dominated communities following disturbance. We investigated co-occurrence patterns of habitat generalist and specialist terrestrial rodents at two spatial scales in the Atlantic Forest, aiming at evaluating the following hypotheses: 1) within-patch spatial niche partitioning promotes coexistence of generalists and specialists, leading to checkerboard presence-absence patterns at small (within-patch) rather than large (among-patch) scales; and 2) the decrease in abundance of specialists due to habitat loss promotes a competitive release of generalists, leading to negative covariance in abundance between generalists and specialists among patches. Drawing on a large data set including 363 sites within three patches in continuous forest, and 45 patches within three landscapes, we used C-scores based on presence-absence and abundance data to evaluate spatial segregation. We found consistent segregation between specialists and generalists at the within-patch rather than among-patch scale, but no consistent negative covariance in abundance between generalists and specialists among patches (as covarying species pairs varied across landscapes). Our findings suggest that spatial patterns caused by competition are scale-dependent, and coexistence of generalists and specialists is promoted by within-patch spatial niche partitioning. However, the influence of competitive release on the proliferation of generalists may be outweighed by other factors in fragmented landscapes.

A teoria ecológica prevê que espécies generalistas de habitat são excluídas por espécies especialistas em hábitats ótimos para as especialistas, e dados empíricos comumente mostram uma mudança de dominância das comunidades - de especialistas para generalistas - após distúrbios. Nós investigamos os padrões de coocorrência de roedores terrestres generalistas e especialistas de habitat em duas escalas espaciais na Mata Atlântica, para testar as seguintes hipóteses: 1) a partição espacial do nicho dento de fragmentos promove a coexistência de generalistas e especialistas, levando a padrões de presença-ausência “tabuleiro de damas” em escalas pequenas (dentro de fragmento) mas não em escalas grandes (entre fragmentos); 2) a diminuição da abundância de especialistas devido à perda de habitat promove uma liberação competitiva de generalistas, levando a covariância negativa da abundância de generalistas e especialistas entre fragmentos. A partir de um grande banco de dados - 363 sítios dentro de três fragmentos de floresta contínua, e 45 fragmentos dentro de três paisagens, usamos C-scores baseados em dados de presença/ausência e abundância para avaliar a segregação espacial. Encontramos segregação consistente entre especialistas e generalistas na escala menor (dentro de fragmentos) e não na maior (entre fragmentos), mas não encontramos covariância negativa na abundância de generalistas e especialistas entre fragmentos (dado que os pares de espécies que covariaram mudou entre as paisagens). Nossos resultados sugerem que padrões espaciais causados por competição são dependentes de escala, e que a coexistência de generalistas e especialistas é promovida pela partição espacial de nicho dentro dos fragmentos. No entanto, a influência da liberação competitiva na proliferação de generalistas pode ser superada por outros fatores em paisagens fragmentadas.

Unmasking structural patterns in incidence matrices: an application to ecological data

Null models have become a crucial tool for understanding structure within incidence matrices across multiple biological contexts. For example, they have been widely used for the study of ecological and biogeographic questions, testing hypotheses regarding patterns of community assembly, species co-occurrence and biodiversity. However, to our knowledge we remain without a general and flexible approach to study the mechanisms explaining such structures. Here, we provide a method for generating 'correlation-informed' null models, which combine the classic concept of null models and tools from community ecology, like joint statistical modelling. Generally, this model allows us to assess whether the information encoded within any given correlation matrix is predictive for explaining structural patterns observed within an incidence matrix. To demonstrate its utility, we apply our approach to two different case studies that represent examples of common scenarios encountered in community ecology. First, we use a phylogenetically informed null model to detect a strong evolutionary fingerprint within empirically observed food webs, reflecting key differences in the impact of shared evolutionary history when shaping the interactions of predators or prey. Second, we use multiple informed null models to identify which factors determine structural patterns of species assemblages, focusing in on the study of nestedness and the influence of site size, isolation, species range and species richness. In addition to offering a versatile way to study the mechanisms shaping the structure of any incidence matrix, including those describing ecological communities, our approach can also be adapted further to test even more sophisticated hypotheses.

Nematodes in a polar desert reveal the relative role of biotic interactions in the coexistence of soil animals

Abiotic factors are major determinants of soil animal distributions and their dominant role is pronounced in extreme ecosystems, with biotic interactions seemingly playing a minor role. We modelled co-occurrence and distribution of the three nematode species that dominate the soil food web of the McMurdo Dry Valleys (Antarctica). Abiotic factors, other biotic groups, and autocorrelation all contributed to structuring nematode species distributions. However, after removing their effects, we found that the presence of the most abundant nematode species greatly, and negatively, affected the probability of detecting one of the other two species. We observed similar patterns in relative abundances for two out of three pairs of species. Harsh abiotic conditions alone are insufficient to explain contemporary nematode distributions whereas the role of negative biotic interactions has been largely underestimated in soil. The future challenge is to understand how the effects of global change on biotic interactions will alter species coexistence.

Tancredi Caruso et al. analyze biodiversity survey data from the McMurdo Dry Valleys, an extreme desert ecosystem in Antarctica in which abiotic factors are thought to determine species distributions. Focusing on three nematode species, they find that abiotic factors alone cannot explain the data and interaction between species have been historically underestimated.

Viral Diversity of House Mice in New York City

The microbiome of wild Mus musculus (house mouse), a globally distributed invasive pest that resides in close contact with humans in urban centers, is largely unexplored. Here, we report analysis of the fecal virome of house mice in residential buildings in New York City, NY. Mice were collected at seven sites in Manhattan, Queens, Brooklyn, and the Bronx over a period of 1 year. Unbiased high-throughput sequencing of feces revealed 36 viruses from 18 families and 21 genera, including at least 6 novel viruses and 3 novel genera. A representative screen of 15 viruses by PCR confirmed the presence of 13 of these viruses in liver. We identified an uneven distribution of diversity, with several viruses being associated with specific locations. Higher mouse weight was associated with an increase in the number of viruses detected per mouse, after adjusting for site, sex, and length. We found neither genetic footprints to known human viral pathogens nor antibodies to lymphocytic choriomeningitis virus.IMPORTANCE Mice carry a wide range of infectious agents with zoonotic potential. Their proximity to humans in the built environment is therefore a concern for public health. Laboratory mice are also the most common experimental model for investigating the pathobiology of infectious diseases. In this survey of mice trapped in multiple locations within New York City over a period of 1 year, we found a diverse collection of viruses that includes some previously not associated with house mice and others that appear to be novel. Although we found no known human pathogens, our findings provide insights into viral ecology and may yield models that have utility for clinical microbiology.

Unraveling Saproxylic Insect Interactions in Tree Hollows from Iberian Mediterranean Forest

Tree hollows are complex microhabitats in which a variety of abiotic and biotic factors shape the community assembly of saproxylic insects. Detecting non-random species co-occurrence patterns is a fundamental goal in ecology in order to understand the assembly mechanisms of communities. We study association patterns of species of Coleoptera and Diptera (Syrphidae), belonging to different trophic guilds, on 72 tree hollows at a local and regional scale in three protected areas in Mediterranean forests using a fixed-fixed null model. Our matrix-level analysis shows a tendency for segregation in species association (species exclusion) at the regional and site levels. However, the high complexity of tree-hollow habitats, offering different resources for a more or less specialized fauna, makes it difficult to prove competition interactions. Indeed, pairwise analysis shows a dominance of non-random aggregation patterns (species coexistence) at the local and regional levels. Both aggregation and segregation of non-random patterns were more common among species from different trophic guilds than within the same guilds, with predators being a common denominator for a high percentage of the inter-guild pairs. Our results suggest that predation and facilitation interactions, together with habitat segregation, are the main factors shaping tree-hollow assemblages, while competition seems to be less important. We conclude that species interactions take an important part of the process of assemblage structuration and special attention should be paid to 'ecosystem engineers' and threatened species in the conservation of tree hollow assemblages.

Factors shaping community assemblages and species co-occurrence of different trophic levels

Species assemblages are the results of various processes, including dispersion and habitat filtering. Disentangling the effects of these different processes is challenging for statistical analysis, especially when biotic interactions should be considered. In this study, we used plants (producers) and leafhoppers (phytophagous) as model organisms, and we investigated the relative importance of abiotic versus biotic factors that shape community assemblages, and we infer on their biotic interactions by applying three-step statistical analysis. We applied a novel statistical analysis, that is, multiblock Redundancy Analysis (mbRA, step 1) and showed that 51.8% and 54.1% of the overall variation in plant and leafhopper assemblages are, respectively, explained by the two multiblock models. The most important blocks of variables to explain the variations in plant and leafhopper assemblages were local topography and biotic factors. Variation partitioning analysis (step 2) showed that pure abiotic filtering and pure biotic processes were relatively less important than their combinations, suggesting that biotic relationships are strongly structured by abiotic conditions. Pairwise co-occurrence analysis (step 3) on generalist leafhoppers and the most common plants identified 40 segregated species pairs (mainly between plant species) and 16 aggregated pairs (mainly between leafhopper species). Pairwise analysis on specialist leafhoppers and potential host plants clearly revealed aggregated patterns. Plant segregation suggests heterogeneous resource availability and competitive interactions, while leafhopper aggregation suggests host feeding differentiation at the local level, different feeding microhabitats on host plants, and similar environmental requirements of the species. Using the novel mbRA, we disentangle for the first time the relative importance of more than five distinct groups of variables shaping local species communities. We highlighted the important role of abiotic processes mediated by bottom-up effects of plants on leafhopper communities. Our results revealed that in-field structure diversification and trophic interactions are the main factors causing the co-occurrence patterns observed.

The distribution of digenean metacercariae within bream (Abramis brama) gill apparatus: preferences, co-occurrence and interactions of parasites

Species-specific microenvironmental preferences and interactions between parasite species have been the focus of many ecological studies. Here, we studied the distribution of ectoparasite species within the gill apparatus of bream (Abramis brama) from Lake Lubāns (Latvia) to establish whether digenean metacercariae: (1) prefer specific patches within the gill apparatus; (2) co-occur in the same patches with monogeneans and copepods within a host individual; and (3) interact with monogeneans and copepods. We recorded all parasites on gill arches of the same host species and used null models to analyse co-occurrences of digenean metacercariae, monogeneans and copepods. Zero-inflated mixture models were used to define the preferred patches of parasites. We found that digenean metacercariae (Bucephalus polymorphus) prefer specific patches of the gill apparatus to encyst, and shared these preferences with monogeneans and copepods, but did not interact with them. We concluded that digenean metacercariae have a species-specific microenvironmental preference to encyst in the gill apparatus and their occurrence (even in high numbers) does not reduce the success of attachment of monogeneans and copepods in the same gill patches.

The Aegean archipelago: a natural laboratory of evolution, ecology and civilisations

The Aegean archipelago, comprising numerous islands and islets with great heterogeneity in topographic, geological, historical and environmental properties, offers an ideal natural laboratory for ecological and evolutionary research, and has been the stage for a very long interaction between human civilizations and local ecosystems. This work presents insights that have been gained from past and current relevant research in the area, highlighting also the importance of the Aegean archipelago as a useful model to address many major questions in biogeography, ecology and evolutionary processes. Among the most interesting findings from such studies concern the role of habitat heterogeneity as the most important determinant of species richness, the development of a new model (Choros) for the species–area–habitats relationship, the mechanistic aspects of the Small Island Effect, the very high rates of species turnover, the lack of a role for interspecific competition in shaping species co-occurrence patterns in most cases, the importance of non adaptive radiation in diversification of several taxa, the insights into the relative roles of vicariance and dispersal in speciation, the understanding of the interplay between human presence and the establishment of exotic species and extinction of indigenous biotas. Concluding, the Aegean archipelago is an ideal stage for research in evolution, ecology and biogeography, and has the potential to become a model study area at a global level, especially for land-bridge, continental islands.

Enhanced facilitation at the extreme end of the aridity gradient in the Atacama Desert: a community-level approach

Plant facilitation is now recognized as an important process in severe environments. However, there is still no agreement on how facilitation changes as conditions become increasingly severe. The classic stress gradient hypothesis (SGH) predicts a monotonic increase in facilitation, which rises in frequency as conditions approach the extreme end of the environmental gradient. However, few studies have evaluated the validity of the SGH at the community level, the level at which it was formulated. Moreover, few studies have tested the SGH at either extreme of the gradient, and very few have excluded the effect of livestock on community response to stress. In line with the SGH, we hypothesized that several spatial pattern summary statistics would change monotonically from the least to the most arid sites, indicating increasingly aggregated patterns. In this study, we performed an evaluation of the SGH both within communities of shrub species and across a large portion of the Atacama Desert, and we isolated the abiotic component of the SGH. Our environmental gradient covered an extreme aridity gradient (< 20-130 mm annual precipitation). To perform point pattern analysis, we established 13 sites with environmental conditions representing four distinct levels of this gradient. Further, we conducted species co-occurrence analyses at 19 sites along the gradient. Both sets of analyses showed stronger positive spatial associations among plants at the most extreme end of the gradient. This was true regardless of whether we included all individuals, only small individuals located around large ones, or individuals in species pairs. Moreover, species tended to show greater co-occurrence as environmental severity increased. This increase in aggregation in the plant community seems to correlate with an increase in the strength of positive interspecific interactions, rather than greater clustering within each species. These monotonic increases in species co-occurrence and spatial association in more severe environments are consistent with some of the predictions of SGH, and collectively these results suggest that as the climate becomes more arid, positive species pairs interactions tend to be prevalent in the community.

Midgut fungal and bacterial microbiota of Aedes triseriatus and Aedes japonicus shift in response to La Crosse virus infection

Understanding how midgut microbial communities of field-collected mosquitoes interact with pathogens is critical for controlling vector infection and disease. We used 16S rRNA and internal transcribed spacer sequencing to characterize the midgut bacterial and fungal communities of adult females of Aedes triseriatus and Aedes japonicus collected as pupae in tree holes, plastic bins and waste tires and their response to La Crosse virus (LACV) infection. For both mosquito species and across all habitat and virus treatments, a total of 62 bacterial operational taxonomic units (OTUs) from six phyla and 21 fungal OTUs from two phyla were identified. The majority of bacterial (92%) and fungal (71%) OTUs were shared between the mosquito species; however, several OTUs were unique to each species. Bacterial and fungal communities of individuals that took either infectious or noninfectious bloodmeals were less diverse and more homogeneous compared to those of newly emerged adults. Interestingly, LACV-infected A. triseriatus and A. japonicus had higher bacterial richness and lower fungal richness compared to individuals that took a noninfectious bloodmeal, suggesting that viral infection was associated with an increase in bacterial OTUs and a decrease in fungal OTUs. For both mosquito species, several OTUs were identified that had both high fidelity and specificity to mosquito midguts that were infected with LACV. Overall, these findings demonstrate that bacterial and fungal communities that reside in mosquito midguts respond to host diet and viral infection and could play a role in modulating vector susceptibility to LACV.

Complex relationships between species niches and environmental heterogeneity affect species co-occurrence patterns in modelled and real communities

Species co-occurrence analysis is commonly used to assess how interspecific interactions dictate community assembly. Non-random co-occurrences, however, may also emerge from niche differences as well as environmental heterogeneity. The relationships between species co-occurrence patterns, environmental heterogeneity and species niches are not fully understood, due to complex interactions among them. To analyse the relationships among these patterns and processes, I developed synthetic community models and analysed a large dataset of tree species across the conterminous United States. Niche overlap and environmental heterogeneity had significant and contrasting effects on species co-occurrence patterns, in both modelled and real communities. Niche breadth, in turn, affected the effect sizes of both variables on species co-occurrence patterns. The effect of niche breadth on the relationship between co-occurrence and niche overlap was markedly consistent between modelled and real communities, while its effect on the relationship between co-occurrence and environmental heterogeneity was mostly consistent between real and modelled data. The results of this analysis highlight the complex and interactive effects of species niche overlap, niche breadth and environmental heterogeneity on species co-occurrence patterns. Therefore, inferring ecological processes from co-occurrence patterns without accounting for these fundamental characteristics of species and environments may lead to biased conclusions.

Holocene shifts in the assembly of plant and animal communities implicate human impacts

Understanding how ecological communities are organized and how they change through time is critical to predicting the effects of climate change. Recent work documenting the co-occurrence structure of modern communities found that most significant species pairs co-occur less frequently than would be expected by chance. However, little is known about how co-occurrence structure changes through time. Here we evaluate changes in plant and animal community organization over geological time by quantifying the co-occurrence structure of 359,896 unique taxon pairs in 80 assemblages spanning the past 300 million years. Co-occurrences of most taxon pairs were statistically random, but a significant fraction were spatially aggregated or segregated. Aggregated pairs dominated from the Carboniferous period (307 million years ago) to the early Holocene epoch (11,700 years before present), when there was a pronounced shift to more segregated pairs, a trend that continues in modern assemblages. The shift began during the Holocene and coincided with increasing human population size and the spread of agriculture in North America. Before the shift, an average of 64% of significant pairs were aggregated; after the shift, the average dropped to 37%. The organization of modern and late Holocene plant and animal assemblages differs fundamentally from that of assemblages over the past 300 million years that predate the large-scale impacts of humans. Our results suggest that the rules governing the assembly of communities have recently been changed by human activity.

Temporal Overlap and Co-Occurrence in a Guild of Sub-Tropical Tephritid Fruit Flies

Studies of community assembly have emphasized snapshot comparisons of spatially replicated samples from "natural" assemblages. Agro-ecosystems are characterized by relatively little habitat heterogeneity and no dispersal barriers for actively flying insects. Therefore, dynamic patterns of species segregation and aggregation are more likely to reflect the direct or indirect effects of species interactions. We studied the temporal organization of a guild of 21 congeneric species of Anastrepha that colonized fruit orchards in Monte Alegre do Sul, São Paulo, Brazil. This assemblage also included the introduced Mediterranean fruit fly Ceratitis capitata. One hundred six consecutive weekly censuses (11 Jan 2002-16 Jan 2004) of flies in guava, loquat, and peach orchards revealed a pattern of minimum abundance during the coldest months of each year (June and July) and a maximum abundance during periods of flowering and fruit ripening. Overall, phenological overlap was greater than expected by chance. However, conditioned on the pattern of seasonal abundances, temporal occurrence and abundance matrices exhibited patterns of significant species segregation and anti-nestedness. In each year, the 3 orchards contained a small number of species pairs that exhibited statistically significant temporal segregation or aggregation. Most aggregated and segregated pairs reflected seasonal shifts in species presences that were not related to variation in air temperature. Most of the significant pairwise associations involved C. capitata: 8 of the 11 segregated pairs and 2 of the 7 aggregated pairs. These results suggest that species interactions between introduced and native species can be an important determinant of species associations in agro-ecosystems.