Null model analyses of presence–absence matrices need a definition of independence

Tests for species interactions that involve the comparison of a statistic calculated from observed matrix of species presences and absences with the distribution of the same statistic generated from a null model have been used by ecologists for about 30 years. We argue that the validity of these tests requires a specific definition of independence. In particular, we note that an assumption that is often made is that all presence–absence matrices with the same row and column totals are equally likely if there is no interaction. However, we show using a simple model for species presences and absences without any species interactions that, in general, this assumption should be made with caution. Our model incorporates a definition of independence, allowing the computation of probabilities of different patterns in the null matrices. Other definitions of independence are possible; one of them is outlined using a new generalized linear model approach for carrying out tests applicable to different null models with or without the assumption of keeping row and column totals fixed. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Morphological assembly mechanisms in Neotropical bat assemblages and ensembles within a landscape

Empirical studies on bat assemblages have shown that richness is not appreciably influenced by local processes such as ecological interactions. However, most of these studies have been done in large areas that include high heterogeneity, and they analyse all bat species within such areas, and thus they may be not reflecting local but supra-community conditions. We followed an ecomorphological approach to assess how bat assemblages of species from the families Phyllostomidae and Mormoopidae, and ensembles of frugivorous bats, are assembled in local habitats within a single landscape. We measured the volume of the space defined by wing morphology and quantified the average distance between species within such a volume. Then, we related these measures to local richness. Such relationships were contrasted against relationships with random assemblages to test for statistical differences. At the ensemble level of organization, we found that the frugivorous bat morphological assembly mechanism is different from random patterns, and it corresponds to the volume-increasing model. On the other hand, bat assembly mechanisms may be ubiquitous at the assemblage level, because groups of species coexisting in a local habitat and delimited only by phylogeny include more than one ecological group with no potential to interact. Assembling processes are crucial to an understanding of species diversity in local communities, and ecomorphological analyses are very promising tools that may help in their study.     

Species composition of dabbling duck assemblages: ecomorphological patterns compared with null models

Ecomorphological patterns of breeding dabbling duck (Anas spp.) assemblages were studied in six regions in northern Europe. Observed spacings among species in terms of bill lamellar density and body length were compared with expected spacings based on null models incorporating different levels of constraints (regional species pools, species relative abundances, lake size and habitat requirements of species). Deviations of observed spacings from expected ones were compared with prey abundance and prey size diversity in the lakes. Observed spacings in terms of body length, but not in terms of bill lamellar density, were greater than expected on the basis of null models. The most abundant species were generally relatively more different than less abundant species in terms of body length but not in terms of bill lamellar density. Deviations between observed and expected spacings in terms of body length were more like those predicted by the competition hypothesis in lakes with low food abundance than in lakes with high food abundance. Patterns in bill lamellar spacings were not related to food abundance nor to food size diversity. In general, patterns in body length spacings were consistent with the competition hypothesis whether the null model used in comparisons was constrained or not.     

Bat assemblages on Neotropical land-bridge islands: nested subsets and null model analyses of species co-occurrence patterns

A fundamental goal of ecology is to understand whether ecological communities are structured according to general assembly rules or are essentially dictated by random processes. In the context of fragmentation, understanding assembly patterns and their mechanistic basis also has important implications for conservation. Using distribution data of 20 bat species collected on 11 islands in Gatún Lake, Panama, we tested for non‐randomness in presence–absence matrices with respect to nestedness and negative species co‐occurrence. We examined the causal basis for the observed patterns and conducted separate analyses for the entire assemblage and for various species submatrices reflecting differences in species’ trophic position and mobility. Furthermore, we explored the influence of weighting factors (area, isolation, abundance) on co‐occurrence analyses. Unweighted analyses revealed a significant negative co‐occurrence pattern for the entire assemblage and for phytophagous bats alone. Weighting analyses by isolation retained a pattern of species segregation for the whole assemblage but nullified the non‐random structure for phytophagous bats and suggested negative associations for animalivores and species with low mobility. Area‐ and abundance‐weighted analyses always indicated random structuring. Bat distributions followed a nested subset structure across islands, regardless of whether all species or different submatrices were analysed. Nestedness was in all cases unrelated to island area but weakly correlated with island isolation for incidence matrices of all species, phytophagous bats, and mobile species. Overall, evidence for negative interspecific interactions indicative of competitive effects was weak, corroborating previous studies based on ecomorphological analyses. Our findings indicate that bat assemblages on our study islands are most strongly shaped by isolation effects and species’ differential movement and colonization ability. From a conservation viewpoint this suggests that even in systems with high fragment–matrix contrast, a purely area‐based approach may be inadequate, and structural and functional connectivity among patches are important to consider in reserve planning.     

Predicting Distribution and Density of European Badger (<Emphasis Type="Italic">Meles Meles</Emphasis>) Setts in Denmark

Predictive models of the spatial distribution and abundance of species based on habitat characteristics are finding increasing use in management and conservation. The European badger attracts interest as a model species both for conservation reasons and because of the important role the species is playing in understanding carnivore sociality. We developed a statistical habitat model based on presence/absence data on badger setts. To maximise the utility of the model in management, we limited the choice of model variables to those that had a clear basis in badger ecology and that could be obtained on a nation-wide digital format. We extrapolated the habitat model to a region in Denmark and developed a threshold-independent sett distribution algorithm to estimate sett densities. The habitat model was simpler than previously published models of badger sett habitat selection, but nevertheless had a predictive ability in excess of 80% judged against independent data. The sett distribution algorithm was able to simultaneously reproduce several observed patterns of sett density and distribution over the probability gradient. It thus represents a significant improvement over threshold-dependent methods used to discriminate between suitable and unsuitable habitat predicted by presence/absence regression models. Our approach demonstrates that a model of badger sett habitat suitability with high predictive power can be obtained using easily accessible map-variables and presence/absence data. This is a prerequisite for using habitat models as predictive tools over large areas. The use of a simple sett distribution algorithm circumvents the common problem of subjectively fixing a threshold to discriminate between suitable and unsuitable habitat. In conjunction the models presented here constitute an important contribution to the management of the badger in Denmark and, upon further validation, possibly to similar regions in Northern Europe.     

Fumaric acid: an overlooked form of fixed carbon in Arabidopsis and other plant species

Photoassimilates are used by plants for production of energy, as carbon skeletons and in transport of fixed carbon between different plant organs. Many studies have been devoted to characterizing the factors that regulate photoassimilate concentrations in different plant species. Most studies examining photoassimilate concentrations in C₃ plants have focused on analyzing starch and soluble sugars. However, work presented here demonstrates that a number of C₃ plants, including the popular model organism Arabidopsis thaliana (L.) Heynh., and agriculturally important plants, such as soybean, Glycine max (L.) Merr., contain significant quantities of fumaric acid. In fact, fumaric acid can accumulate to levels of several milligrams per gram fresh weight in Arabidopsis leaves, often exceeding those of starch and soluble sugars. Fumaric acid is a component of the tricarboxylic acid cycle and, like starch and soluble sugars, can be metabolized to yield energy and carbon skeletons for production of other compounds. Fumaric acid concentrations increase with plant age and light intensity in Arabidopsis leaves. Moreover, Arabidopsis phloem exudates contain significant quantities of fumaric acid, raising the possibility that fumaric acid may function in carbon transport. Received: 11 February 2000 / Accepted: 1 April 2000     

Nested communities,invasive species and Holocene extinctions: evaluating the power of a potential conservation tool

General ecological methods and models that require a minimum amount of information yet are still able to inform conservation planning are particularly valuable. Nested subset analysis has been advocated as such a tool for the prediction of extinction-prone species and populations. However, such advocacy has not been without skepticism and debate, and in the majority of published examples assessing extinction vulnerability, actual extinctions are based on assumptions rather than direct evidence. Here, we empirically test the power of nested subset analysis to predict extinction-prone species, using documented Holocene insular mammal extinctions on three island archipelagos off the west coast of North America. We go on to test whether the introduction of invasive mammals promotes nestedness on islands via extinction. While all three archipelagos were significantly nested before and after the extinction events, nested subset analysis largely failed to predict extinction patterns. We also failed to detect any correlations between the degree of nestedness at the genus-level with area, isolation, or species richness and extinction risk. Biogeography tools, such as nested subset analysis, must be critically evaluated before they are prescribed widely for conservation planning. For these island archipelagos, it appears detailed natural history and taxa-specific ecology may prove critical in predicting patterns of extinction risk.     

Patterns in the co-occurrence of fish species in streams: the role of site suitability,morphology and phylogeny versus species interactions

A number of studies at large scales have pointed out that abiotic factors and recolonization dynamics appear to be more important than biotic interactions in structuring stream-fish assemblages. In contrast, experimental and field studies at small scales show the importance of competition among stream fishes. However, given the highly variable nature of stream systems over time, competition may not be intense enough to generate large-scale complementary distributions via competitive exclusion. Complementary distribution is a recurrent pattern observed in fish communities across stream gradients, though it is not clear which instances of this pattern are due to competitive interactions and which to individual species requirements. In this study, I introduce a series of null models developed to provide a more robust evaluation of species associations by facilitating the distinction between different processes that may shape species distributions and community assembly. These null models were applied to test whether conspicuous patterns in species co-occurrences are more consistent with their differences in habitat use, morphological features and/or phylogenetic constraints, or with species interactions in fish communities in the streams of a watershed in eastern Brazil. I concluded that patterns in species co-occurrences within the studied system are driven by common species-habitat relationships and species interactions may not play a significant role in structuring these communities. I suggest that large-scale studies, where adequate designs and robust analytical tools are applied, can contribute substantially to understanding the importance of different types of processes in structuring stream-fish communities.     

The necessity of functional proteomics: protein species and molecular function elucidation exemplified by in vivo alpha A crystallin N-terminal truncation

Summary. Ten years after the establishment of the term proteome, the science surrounding it has yet to fulfill its potential. While a host of technologies have generated lists of protein names, there are only a few reported studies that have examined the individual proteins at the covalent chemical level defined as protein species in 1997 and their function. In the current study, we demonstrate that this is possible with two-dimensional gel electrophoresis (2-DE) and mass spectrometry by presenting clear evidence of in vivo N-terminal alpha A crystallin truncation and relating this newly detected protein species to alpha crystallin activity regulation by protease cleavage in the healthy young murine lens. We assess the present state of technology and suggest a shift in resources and paradigm for the routine attainment of the protein species level in proteomics.     

Linking limitation to species composition: importance of inter- and intra-specific variation in grazing resistance

Short-term responses of producers highlight that key nutrients (e.g., N, P)—or combinations of these nutrients—limit primary production in aquatic and terrestrial ecosystems. These discoveries continue to provide highly valuable insights, but it remains important to ask whether nutrients always predominantly limit producers despite wide variation in nutrient supply and herbivory among systems. After all, predictions from simple food chain models (derived here) readily predict that limitation by grazers can exceed that by nutrients, given sufficient enrichment. However, shifts in composition of producers and/or increasing dominance of invulnerable stages of a producer can, in theory, reduce grazer limitation and retain primacy of nutrient limitation along nutrient supply gradients. We observed both mechanisms (inter- and intra-species variation in vulnerability to herbivory) working in a two-part mesocosm experiment. We incubated diverse benthic algal assemblages for several months either in the presence or absence of benthic macro-grazers in mesocosms that spread a broad range of nutrient supply. We then conducted short-term assays of nutrient and grazer limitation on these communities. In the “historically grazed” assemblages, we found shifts from more edible, better competitors to more resistant producers over enrichment gradients (as anticipated by the food web model built with a tradeoff in resistance vs. competitive abilities). However, contrary to our expectations, “historically ungrazed” assemblages became dominated by producers with vulnerable juvenile forms but inedible adult forms (long filaments). Consequently, we observed higher resource limitation rather than grazer limitation over this nutrient supply gradient in both “historically grazed” (expected) and “historically ungrazed” (not initially expected). Thus, via multiple, general mechanisms involving resistance to grazing (changes in species composition or variation in stage-structured vulnerability), producer assemblages should remain more strongly or as strongly limited by nutrients than grazers, even over large enrichment gradients. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Effects of environmental heterogeneity on species diversity and composition of terrestrial bryophyte assemblages in tropical montane forests of southern Ecuador

Background: Most studies on tropical bryophytes deal with epiphytic species. This is the first ecological study of tropical forests that focuses specifically on terrestrial bryophytes.

Aim: To investigate the differences between slope and ridge environments in upper montane forests of southern Ecuador in terms of species diversity (richness, abundance), species composition and life forms of terrestrial bryophytes.

Methods: We used Non-metric Multidimensional Scaling (NMDS) to group bryophyte relevés by study location, habitat type and exposure class. Species indicator values were calculated and compared for different habitats.

Results: In total, 140 species were recorded, the majority being liverworts. NMDS analyses and Mantel correlations clearly separated between slope and ridge relevés, and between sunny and shaded microhabitats on ridges. Bryophyte life forms also showed different distribution patterns in slope and in ridge habitats. Mosses were more prominent in sunny than in shaded microhabitats.

Conclusions: Environmental differentiation between ridges and slopes, and small-scale variation in microclimatic conditions caused by differences in exposure, were stronger predictors of species richness and composition than geographical distance between study sites.     

Spatial patterns of coexistence of competing species in patchy habitat

The single-species spatially realistic patch occupancy metapopulation model is, in this study, extended to a metacommunity of many competing species. Competition is assumed to reduce the local carrying capacity (effective patch area), which in turn increases local extinction rates and reduces colonization rates because of smaller population sizes. Each species is described by three parameters: pre-competitive abundance (equilibrium incidence of patch occupancy, which reflects the rate of colonization in relation to extinction rate), the spatial range of migration, and competitive ability. The model ignores spatio–temporal correlations caused by interspecific interactions, because in metacommunities of unequal competitors inhabiting heterogeneous landscapes, correlations in the occurrence of species are driven more by patch heterogeneity than by competition. The model allows the calculation of multispecies equilibria in patchy habitats without simulations. In general, the number of coexisting species in the metacommunity increases with decreasing strength of competition, increasing rate of colonization, and decreasing range of migration. Habitat heterogeneity in the form of spatial variation in patch areas tends to facilitate coexistence. Poor competitors may coexist with superior competitors in the patch network if the former have higher colonization rates (competition–colonization trade-off). When migration distances are short, competition leads to spatial pattern formation: Species tend to have restricted spatial distributions in the network, but contrary to intuitive expectations, often the distributions of many species are nested. Having more dispersive species enhances both local and global diversity, whereas more local migration decreases local but increases global diversity.     

运用聚类分析与Google Maps于大量物种出现记录之研究

物种出现记录包含博物馆动物标本、植物标本、生态调查与物种观察等资料。在台湾生物多样性信息机构(Taiwan Biodiversity Information Facility,TaiBIF)物种出现记录整合平台中,已整合台湾26个数据集,包含超过150万笔物种出现记录,其中约有85%的数据具有地理信息。我们利用数据库中所汇整的鲤科数据,包括11个数据集、超过8,800笔出现记录数据,利用网格式、切割式与密度式3种聚类分析算法分别绘制出不同的空间可视化结果,藉此解决大量物种出现记录于Google Maps上呈现效能与可视化不佳之问题。同时我们也探讨了3种聚类分析法之结果与鲤科的专家意见范围地图(expertopinion range maps)比对的差异。期望透过本研究可快速且有效地呈现物种分布资料,进而帮助研究者挖掘出大量数据所隐含的知识,并为生态保育提供重要参考。     

Ontogeny of flight in the little brown bat,Myotis lucifugus: behavior,morphology, and muscle histochemistry

Summary Postnatal changes in wing morphology, flight ability, muscle morphology, and histochemistry were investigated in the little brown bat, Myotis lucifugus. The pectoralis major, acromiodeltoideus, and quadriceps femoris muscles were examined using stains for myofibrillar ATPase, succinate dehydrogenase (SDH), and mitochondrial -glycerophosphate dehydrogenase (-GPDH) enzyme reactions. Bats first exhibited spontaneous, drop-evoked flapping behavior at 10 days, short horizontal flight at 17 days, and sustained flight at 24 days of age. Wing loading decreased and aspect ratio increased during postnatal development, each reaching adult range before the onset of sustained flight. Histochemically, fibers from the three muscles were undifferentiated at birth and had lower oxidative and glycolytic capacities compared to other age groups. Cross-sectional areas of fibers from the pectoralis and acromiodeltoideus muscles increased significantly at an age when dropevoked flapping behavior was first observed, suggesting that the neuromuscular mechanism controlling flapping did not develop until this time. Throughout the postnatal growth period, pectoralis and acromiodeltoideus muscle mass and fiber cross-sectional area increased significantly. By day 17 the pectoralis muscle had become differentiated in glycolytic capacity, as indicated by the mosaic staining pattern for -GPDH. By contrast, the quadriceps fibers were relatively large at birth and slowly increased in size during the postnatal period. Fiber differentiation was evident at the time young bats began to fly, as indicated by a mosaic pattern of staining for myosin ATPase. These results indicate that flight muscles (pectoralis and acromiodeltoideus) are less well developed at birth and undergo rapid development just before the onset of flight. By contrast the quadriceps femoris muscle, which is required for postural control, is more developed at birth than the flight muscles and grows more slowly during subsequent development.     

The effect of density-dependent treatment and behavior change on the dynamics of HIV transmission

In this work, we propose a model for heterosexual transmission of HIV/AIDS in a population of varying size with an intervention program in which treatment and/or behavior change of the infecteds occur as an increasing function of the density of the infected class in the population. This assumption has socio-economic implications which is important for public health considerations since density-dependent treatment/behavior change may be more cost-saving than a program where treatment/behavior change occurs linearly with respect to the number of infecteds. We will make use of the conservation law of total sexual contacts which enables us to reduce the two-sex model to a simpler one-sex formulation. Analytical results will be given. Unlike a similar model with linear treatment/behavior change in Hsieh (1996) where conditions were obtained for the eradication of disease, we will show that density-dependent treatment/behavior change cannot eradicate the disease if the disease is able to persist without any treatment/behavior change. This work demonstrates the need to further understand how treatment/behavior change occurs in a society with varying population.     

Mutation,selection, and ancestry in branching models: a variational approach

We consider the evolution of populations under the joint action of mutation and differential reproduction, or selection. The population is modelled as a finite-type Markov branching process in continuous time, and the associated genealogical tree is viewed both in the forward and the backward direction of time. The stationary type distribution of the reversed process, the so-called ancestral distribution, turns out as a key for the study of mutation–selection balance. This balance can be expressed in the form of a variational principle that quantifies the respective roles of reproduction and mutation for any possible type distribution. It shows that the mean growth rate of the population results from a competition for a maximal long-term growth rate, as given by the difference between the current mean reproduction rate, and an asymptotic decay rate related to the mutation process; this tradeoff is won by the ancestral distribution. We then focus on the case when the type is determined by a sequence of letters (like nucleotides or matches/mismatches relative to a reference sequence), and we ask how much of the above competition can still be seen by observing only the letter composition (as given by the frequencies of the various letters within the sequence). If mutation and reproduction rates can be approximated in a smooth way, the fitness of letter compositions resulting from the interplay of reproduction and mutation is determined in the limit as the number of sequence sites tends to infinity. Our main application is the quasispecies model of sequence evolution with mutation coupled to reproduction but independent across sites, and a fitness function that is invariant under permutation of sites. In this model, the fitness of letter compositions is worked out explicitly. In certain cases, their competition leads to a phase transition.      

Fish species richness and incidence patterns in isolated and connected stream pools: effects of pool volume and spatial position

I tested the effects of pool size and spatial position (upstream or downstream) on fish assemblage attributes in isolated and connected pools in an upland Oklahoma stream, United States. I hypothesized that there would be fundamental differences between assemblages in these two pool types due to the presence or absence of colonization opportunities. Analyses were carried out at three ecological scales: (1) the species richness of pool assemblages, (2) the species composition of pool assemblages, and (3) the responses of individual species. There were significant species-volume relationships for isolated and connected pools. However, the relationship was weaker and there were fewer species, on average, in isolated pools. For both pool types, species incidences were significantly nested such that species-poor pools tended to be subsets of species-rich pools, a common pattern that ultimately results from species-specific differences in colonization ability and/or extinction susceptibility. To examine the potential importance of these two processes in nestedness patterns in both pool types, I made the following two assumptions: (1) probability of extinction should decline with increasing pool size, and (2) probability of immigration should decline in an upstream direction (increasing isolation). When ordered by pool volume, only isolated pools were significantly nested suggesting that these assemblages were extinction-driven. When ordered by spatial position, only connected pools were significantly nested (more species downstream) suggesting that differences in species-specific dispersal abilities were important in structuring these assemblages. At the individual-species level, volume was a significant predictor of occurrence for three species in isolated pools. In connected pools, two species showed significant position effects, one species showed a pool volume effect, and one species showed pool volume and position effects. These results demonstrate that pool size and position within a watershed are important determinants of fish species assemblage structure, but their importance varies with the colonization potential of the pools. Isolated pool assemblages are similar to the presumed relaxed faunas of montane forest fragments and land bridge islands, but at much smaller space and time scales. Received: 6 December 1996 / Accepted: 10 December 1996     

Breeding dispersal in a heterogeneous landscape: the influence of habitat and nesting success in greater snow geese

Despite numerous studies on breeding dispersal, it is still unclear how habitat heterogeneity and previous nesting success interact to determine nest-site fidelity at various spatial scales. In this context, we investigated factors affecting breeding dispersal in greater snow geese (Anser caerulescens atlanticus), an Arctic breeding species nesting in two contrasting habitats (wetlands and mesic tundra) with variable pattern of snowmelt at the time of settlement in spring. From 1994 to 2005, we monitored the nesting success and breeding dispersal of individually marked females. We found that snow geese showed a moderate amount of nest-site fidelity and considerable individual variability in dispersal distance over consecutive nesting attempts. This variability can be partly accounted for by the annual timing of snowmelt. Despite this environmental constraint, habitat differences at the colony level consistently affected nesting success and settlement patterns. Females nesting in wetlands had higher nesting success than those nesting in mesic tundra. Moreover, geese responded adaptively to spatial heterogeneity by showing fidelity to their nesting habitat, independently of snowmelt pattern. From year to year, geese were more likely to move from mesic to high-quality wetland habitat, regardless of previous nesting success and without cost on their subsequent nesting performance. The unpredictability of snowmelt and the low cost of changing site apparently favour breeding-site dispersal although habitat quality promotes fidelity at the scale of habitat patches.