A new approach to interspecific synchrony in population ecology using tail association

Standard methods for studying the association between two ecologically important variables provide only a small slice of the information content of the association, but statistical approaches are available that provide comprehensive information. In particular, available approaches can reveal tail associations, that is, accentuated or reduced associations between the more extreme values of variables. We here study the nature and causes of tail associations between phenological or population‐density variables of co‐located species, and their ecological importance. We employ a simple method of measuring tail associations which we call the partial Spearman correlation. Using multidecadal, multi‐species spatiotemporal datasets on aphid first flights and marine phytoplankton population densities, we assess the potential for tail association to illuminate two major topics of study in community ecology: the stability or instability of aggregate community measures such as total community biomass and its relationship with the synchronous or compensatory dynamics of the community''s constituent species; and the potential for fluctuations and trends in species phenology to result in trophic mismatches. We find that positively associated fluctuations in the population densities of co‐located species commonly show asymmetric tail associations; that is, it is common for two species’ densities to be more correlated when large than when small, or vice versa. Ordinary measures of association such as correlation do not take this asymmetry into account. Likewise, positively associated fluctuations in the phenology of co‐located species also commonly show asymmetric tail associations. We provide evidence that tail associations between two or more species’ population‐density or phenology time series can be inherited from mutual tail associations of these quantities with an environmental driver. We argue that our understanding of community dynamics and stability, and of phenologies of interacting species, can be meaningfully improved in future work by taking into account tail associations.     

鹞落坪国家级自然保护区吊罐井物种分布格局与种间关联的多尺度分析

为了探讨群落结构形成机制,采用点格局分析方法,对鹞落坪国家级自然保护区吊罐井核心区1hm2样地内36个优势物种在0~50m尺度上的空间分布格局和种间关联性进行了研究。结果显示,36个物种的空间分布格局在小尺度上主要呈聚集分布,在较大尺度上呈随机分布的物种比例最高;对36个物种种间关联性的分析发现,无关联类型是该样地物种间相互作用的主要表现形式,并且在10m尺度上无关联种对的比例最低,正相关和负相关的种对在所有尺度上的比例均未超过20%。群落结构的次生性、生境异质性和扩散限制性是导致该样地物种分布和种间关系格局的主要原因。本研究结果有助于了解物种共存的尺度依赖性特征及其形成机制。     

Dual mycorrhizal colonization of forest-dominating tropical trees and the mycorrhizal status of non-dominant tree and liana species

The contribution of mycorrhizal associations to maintaining tree diversity patterns in tropical rain forests is poorly known. Many tropical monodominant trees form ectomycorrhizal (EM) associations, and there is evidence that the EM mutualism contributes to the maintenance of monodominance. It is assumed that most other tropical tree species form arbuscular mycorrhizal (AM) associations, and while many mycorrhizal surveys have been done, the mycorrhizal status of numerous tropical tree taxa remains undocumented. In this study, we tested the assumption that most tropical trees form AM associations by sampling root vouchers from tree and liana species in monodominant Dicymbe corymbosa forest and an adjacent mixed rain forest in Guyana. Roots were assessed for the presence/absence of AM and EM structures. Of the 142 species of trees and lianas surveyed, three tree species (the monodominant D. corymbosa, the grove-forming D. altsonii, and the non-dominant Aldina insignis) were EM, 137 were exclusively AM, and two were non-mycorrhizal. Both EM and AM structures were observed in D. corymbosa and D. altsonii. These results provide empirical data supporting the assumption that most tropical trees form AM associations for this region in the Guiana Shield and provide the first report of dual EM/AM colonization in Dicymbe species. Dual colonization of the Dicymbe species should be further explored to determine if this ability contributes to the establishment and maintenance of site dominance. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Species importance in a heterospecific foraging association network

There is a growing recognition of the need to integrate non‐trophic interactions into ecological networks for a better understanding of whole‐community organization. To achieve this, the first step is to build networks of individual non‐trophic interactions. In this study, we analyzed a network of interdependencies among bird species that participated in heterospecific foraging associations (flocks) in an evergreen forest site in the Western Ghats, India. We found the flock network to contain a small core of highly important species that other species are strongly dependent on, a pattern seen in many other biological networks. Further, we found that structural importance of species in the network was strongly correlated to functional importance of species at the individual flock level. Finally, comparisons with flock networks from other Asian forests showed that the same taxonomic groups were important in general, suggesting that species importance was an intrinsic trait and not dependent on local ecological conditions. Hence, given a list of species in an area, it may be possible to predict which ones are likely to be important. Our study provides a framework for the investigation of other heterospecific foraging associations and associations among species in other non‐trophic contexts.     

Mapping the potential mycorrhizal associations of the conterminous United States of America

Mycorrhizal associations are recognized as key symbioses in a changing world, yet our understanding of their geographic distribution and temporal dynamics remains limited. We combined data on mycorrhizal associations and historical dominant vegetation to map the pre-European Settlement mycorrhizal associations of the conterminous United States of America (USA). As a demonstration of the map's utility, we estimated changes in mycorrhizal associations due to urbanization, agriculture and the establishment of non-native species in two regions. We found that the conterminous USA was dominated by vegetation associated with arbuscular mycorrhizas, but that ∼40% of vegetation types included multiple mycorrhizal associations. Shifting land use to agriculture and the introduction of non-native species has disproportionately affected ectomycorrhizas, as did urbanization. These preliminary results set a baseline for mycorrhizal biogeography of the USA and illustrate how synthesis of available data can help us understand the impact of anthropogenic changes on an important mutualism.     

Making the most of population genomic data to understand the importance of chromosomal inversions for adaptation and speciation

Chromosomal inversions are increasingly found to differentiate locally adapted populations. This adaptive role is predictable because reduced recombination protects allelic combinations from gene flow. However, we are far from understanding how frequently inversions contribute to local adaptation and how widespread this phenomenon is across species. In a “From the Cover” article in this issue of Molecular Ecology, Huang, Andrew, Owens, Ostevik, and Rieseberg (2020) provide an important step towards this goal not only by finding adaptive inversions in a sunflower ecotype, but also by reversing the approach used to investigate the link between adaptation and inversions. Most studies compare two phenotypes and uncover divergence at a few regions, of which some can subsequently be identified as inversions. In contrast, Huang et al first catalogue putative inversions and then test genotype–environment associations, which allows them to ask systematically whether inversions may be adaptive and in which ecological contexts. They achieve that by revisiting a previous reduced‐representation sequencing (RAD‐sequencing) data set, demonstrating the suitability of this method to detect inversions in species with limited genomic resources. As such, Huang et al pave the way for a better understanding of the evolutionary role of structural genomic variation and highlight that accounting for inversions in population genomics is now possible, and much needed, in a wider range of organisms.     

Statistical discovery of site inter-dependencies in sub-molecular hierarchical protein structuring

Background

Much progress has been made in understanding the 3D structure of proteins using methods such as NMR and X-ray crystallography. The resulting 3D structures are extremely informative, but do not always reveal which sites and residues within the structure are of special importance. Recently, there are indications that multiple-residue, sub-domain structural relationships within the larger 3D consensus structure of a protein can be inferred from the analysis of the multiple sequence alignment data of a protein family. These intra-dependent clusters of associated sites are used to indicate hierarchical inter-residue relationships within the 3D structure. To reveal the patterns of associations among individual amino acids or sub-domain components within the structure, we apply a k-modes attribute (aligned site) clustering algorithm to the ubiquitin and transthyretin families in order to discover associations among groups of sites within the multiple sequence alignment. We then observe what these associations imply within the 3D structure of these two protein families.

Results

The k-modes site clustering algorithm we developed maximizes the intra-group interdependencies based on a normalized mutual information measure. The clusters formed correspond to sub-structural components or binding and interface locations. Applying this data-directed method to the ubiquitin and transthyretin protein family multiple sequence alignments as a test bed, we located numerous interesting associations of interdependent sites. These clusters were then arranged into cluster tree diagrams which revealed four structural sub-domains within the single domain structure of ubiquitin and a single large sub-domain within transthyretin associated with the interface among transthyretin monomers. In addition, several clusters of mutually interdependent sites were discovered for each protein family, each of which appear to play an important role in the molecular structure and/or function.

Conclusions

Our results demonstrate that the method we present here using a k- modes site clustering algorithm based on interdependency evaluation among sites obtained from a sequence alignment of homologous proteins can provide significant insights into the complex, hierarchical inter-residue structural relationships within the 3D structure of a protein family.

     

Metabolism, function and mass spectrometric analysis of bis(monoacylglycero)phosphate and cardiolipin

Polyglycerophospholipids (PGPLs) such as bismonoacylglycerophosphate (BMP) and cardiolipin are important membrane phospholipid species for the maintenance of membrane integrity. While BMP serves as membrane curvature regulator in multivesicular bodies for efficient lysosomal enzyme function, cardiolipin stabilizes the electron transfer complex in the inner mitochondrial membrane, which is crucial for physiological ATP production. Beside their membrane modulatory functions PGPLs play an important role in various signaling events. Although a number of disease associations were found for PGPL species, detailed information about their molecular role still remains unknown. This article reviews the known biological functions of PGPLs and the existing mass spectrometric methods. We discuss the different analytical strategies and how ESI–MS/MS can expand our understanding of PGPL homeostasis.     

Habitat niche partitioning by 16 species of Myristicaceae in Amazonian Ecuador

The distribution and spatial pattern of plants in tropical forests have important implications for how species interact with each other and their environments. In this article we use a large-scale permanent census plot to address if the coexistence of 16 co-occurring species of Myristicaceae is aided by topographic and light gradient niche partitioning. We used a second order spatial pattern analysis based on Ripley’s K function to describe species’ distributions and associations among species, and a torus translation procedure to test for associations with three topographically defined habitats. A majority of species displayed spatial aggregation and over half had one or more significant habitat associations. Four species were associated with the ridge habitat, four species with slope habitat, and two with valley habitat. Seven other species showed no habitat association. Within each habitat group, species exhibited a variety of distributions in relation to light availability. However, habitat associations were largely unexplained by differential rates of mortality, growth or recruitment over a 5-year interval. We conclude that although in principle partitioning of the topographic and light environments may double or treble the number of species able to coexist, there is no evidence that partitioning of physical habitats can explain the coexistence of all 16 of these closely related species.     

孔雀河下游断流河道的环境特征及物种间关系

基于孔雀河下游断流河道的环境因子和植被样地数据,采用聚类与CCA排序法,分析了生境的退化特征以及物种间的相互关系,结果表明:1)断流河道退化生境分为绿洲-荒漠过渡类型、轻度荒漠化类型和盐土荒漠化类型。绿洲-荒漠过渡类型地下水位低、盐分含量相对较低,植被盖度相对较高,土壤维持着原砂质壤土,为潜在退化型;轻度荒漠化类型地下水位、土壤质地与含盐量与前者基本相同,土壤未明显退化,但植被盖度低于10%,植物种类与个体数目都较低,属于轻度退化型;盐土荒漠化类型地下水位高、盐分含量高,土壤机械组成中砂粒比重较大、无建群种幸存,属于重度退化型。2)绿洲-荒漠过渡类型总体联结性为显著正联结,正负联结比小于1,生态系统表现为建群种维系物种关系的不稳定状态;轻度荒漠化类型总体联结性为不显著负联结,正负联结比小于1,表现出生态系统进入退化演替的阶段;盐土荒漠化类型总体联结性为显著正联结,正负联结比大于1,表现出重度退化群落的种间平衡状态,物种间以达到稳定共存,其中,真盐生植物对这种平衡的维持起着重要的作用。3)CCA排序表明,绿洲-荒漠过渡类型形成以胡杨为中心的种间正联结,幸存于盐分适中、水分养分相对较高的生境;轻度荒漠化类型,形成以多枝柽柳与刚毛柽柳相互依存的不显著负联结,幸存于土壤养分、水分相对较低的生境;盐土荒漠化类型形成以盐爪爪、盐节木、盐穗木等真盐生植物维系的显著正联结,幸存于土壤贫瘠、地下水位浅、盐分含量高、沙化严重的生境。     

A network approach for inferring species associations from co‐occurrence data

Positive and negative associations between species are a key outcome of community assembly from regional species pools. These associations are difficult to detect and can be caused by a range of processes such as species interactions, local environmental constraints and dispersal. We integrate new ideas around species distribution modeling, covariance matrix estimation, and network analysis to provide an approach to inferring non‐random species associations from local‐ and regional‐scale occurrence data. Specifically, we provide a novel framework for identifying species associations that overcomes three challenges: 1) correcting for indirect effects from other species, 2) avoiding spurious associations driven by regional‐scale distributions, and 3) describing these associations in a multi‐species context. We highlight a range of research questions and analyses that this framework is able to address. We show that the approach is statistically robust using simulated data. In addition, we present an empirical analysis of > 1000 North American tree communities that gives evidence for weak positive associations among small groups of species. Finally, we discuss several possible extensions for identifying drivers of associations, predicting community assembly, and better linking biogeography and community ecology.     

Low specificity and nested subset structure characterize mycorrhizal associations in five closely related species of the genus Orchis

Most orchid species rely on mycorrhizae to complete their life cycle. Despite a growing body of literature identifying orchid mycorrhizal associations, the nature and specificity of the association between orchid species and mycorrhizal fungi remains largely an open question. Nonetheless, better insights into these obligate plant–fungus associations are indispensable for understanding the biology and conservation of orchid populations. To investigate orchid mycorrhizal associations in five species of the genus Orchis (O. anthropophora, O. mascula, O. militaris, O. purpurea, and O. simia), we developed internal transcribed spacer‐based DNA arrays from extensive clone library sequence data sets, enabling rapid and simultaneous detection of a wide range of basidiomycetous mycorrhizal fungi. A low degree of specificity was observed, with two orchid species associating with nine different fungal partners. Phylogenetic analysis revealed that the majority of Orchis mycorrhizal fungi are members of the Tulasnellaceae, but in some plants, members of the Thelephoraceae, Cortinariaceae and Ceratobasidiaceae were also found. In all species except one (O. mascula), individual plants associated with more than one fungus simultaneously, and in some cases, associations with ≥3 mycorrhizal fungi at the same time were identified. Nestedness analysis showed that orchid mycorrhizal associations were significantly nested, suggesting asymmetric specialization and a dense core of interactions created by symmetric interactions between generalist species. Our results add support to the growing literature that multiple associations may be common among orchids. Low specificity or preference for a widespread fungal symbiont may partly explain the wide distribution of the investigated species.     

Spatial patterns and associations of dominant woody species in desert–oasis ecotone of South Junggar Basin,NW China

Spatial patterns and associations of plant species are important for revealing how species interact with each other and with the environment, and hence have important implications for the understanding of species interaction and underlying ecological processes with apparent patterns in temperate desert vegetation. In this paper, the function g (r) was used to characterize the spatial patterns and associations of four dominant woody species in three 1-ha desert plant plots in the desert–oasis ecotone of South Junggar Basin, NW China. The complete spatial randomness null model showed four species exhibited significant aggregations at small scales (<20 m). Anabasis aphylla and Nitraria roborowskii, Haloxylon ammodendron and Reaumuria songarica were spatially positive associations at small scales with the independent null model, while A. aphylla and H. ammodendron, A. aphylla and R. songarica, R. songarica and N. roborowskii species pairs exhibited negative associations at small or moderate scales (20–60 m) in our study. The random labeling null model showed that dead standing plants of A. aphylla were largely determined by the combined effects of intra- and interspecific competition. In addition, the results also indicated that the two main factors of habitat heterogeneity and sandy desertification play important roles in determining spatial distribution patterns and associations of woody species in the desert–oasis ecotone of South Junggar Basin. Thus, the differences of species features in spatial patterns and associations should be paid more attention when planning afforestation and developing conservation strategies.     

Spatial patterns and interspecific relations analysis help to better understand species distribution patterns in a Mediterranean high mountain grassland

The aim of this study is to investigate distribution patterns of species occurring at definite spatial scales and to address the main ecological factors that structure a Mediterranean high mountain grassland community. Following the protocols of the GLORIA long-term study, four summits were sampled during the year 2006 and 2007 in Sierra de Guadarrama (Spain). The data recorded on two of those summits were analysed for spatial pattern of species, interspecific associations and environmental relationships were determined from the recorded data. This was done by multivariate analyses (ordinations), SADIE analyses (distance indices) and Pearson χ² test, respectively. Results showed that chamaephyte species were more abundant in less disturbed situations whilst caespitose hemicryptophyte species and mosses were more abundant in disturbed ones. High mountain species were more abundant in less disturbed situations. Higher spatial heterogeneity was detected on northern and eastern slopes, which could be related to the increased environmental severity of northern exposures and to the greater intensity of the disturbances on eastern exposures. The frequent aggregated distribution of a few species may be related to facilitation processes due to their frequent participation in positive interactions. Finally we observed that positive pairwise associations of species were more frequent than negative associations. The relevance of the studied relationships lied in a better understanding of those threats on high mountain biodiversity induced by climate warming.     

Q&amp;A: Friends (but sometimes foes) within: the complex evolutionary ecology of symbioses between host and microbes

Over the past decade, there has been a pronounced shift in the study of host–microbe associations, with recognition that many of these associations are beneficial, and often critical, for a diverse array of hosts. There may also be pronounced benefits for the microbes, though this is less well empirically understood. Significant progress has been made in understanding how ecology and evolution shape simple associations between hosts and one or a few microbial species, and this work can serve as a foundation to study the ecology and evolution of host associations with their often complex microbial communities (microbiomes).     

Tropical Plant–Herbivore Networks: Reconstructing Species Interactions Using DNA Barcodes

Plants and their associated insect herbivores, represent more than 50% of all known species on earth. The first step in understanding the mechanisms generating and maintaining this important component of biodiversity is to identify plant-herbivore associations. In this study we determined insect-host plant associations for an entire guild of insect herbivores using plant DNA extracted from insect gut contents. Over two years, in a tropical rain forest in Costa Rica (La Selva Biological Station), we recorded the full diet breadth of rolled-leaf beetles, a group of herbivores that feed on plants in the order Zingiberales. Field observations were used to determine the accuracy of diet identifications using a three-locus DNA barcode (rbcL, trnH-psbA and ITS2). Using extraction techniques for ancient DNA, we obtained high-quality sequences for two of these loci from gut contents (rbcL and ITS2). Sequences were then compared to a comprehensive DNA barcode library of the Zingiberales. The rbcL locus identified host plants to family (success/sequence = 58.8%) and genus (success/sequence = 47%). For all Zingiberales except Heliconiaceae, ITS2 successfully identified host plants to genus (success/sequence = 67.1%) and species (success/sequence = 61.6%). Kindt’s sampling estimates suggest that by collecting ca. four individuals representing each plant-herbivore interaction, 99% of all host associations included in this study can be identified to genus. For plants that amplified ITS2, 99% of the hosts can be identified to species after collecting at least four individuals representing each interaction. Our study demonstrates that host plant identifications at the species-level using DNA barcodes are feasible, cost-effective, and reliable, and that reconstructing plant-herbivore networks with these methods will become the standard for a detailed understanding of these interactions.     

Reptiles Under the Conservation Umbrella of the Greater Sage-Grouse

In conservation paradigms, management actions for umbrella species also benefit co-occurring species because of overlapping ranges and similar habitat associations. The greater sage-grouse (Centrocercus urophasianus) is an umbrella species because it occurs across vast sagebrush ecosystems of western North America and is the recipient of extensive habitat conservation and restoration efforts that might benefit sympatric species. Biologists' understanding of how non-target species might benefit from sage-grouse conservation is, however, limited. Reptiles, in particular, are of interest in this regard because of their relatively high diversity in shrublands and grasslands where sage-grouse are found. Using spatial overlap of species distributions, land cover similarity statistics, and a literature review, we quantified which reptile species may benefit from the protection of intact sage-grouse habitat and which may be affected by recent (since about 1990) habitat restoration actions targeting sage-grouse. Of 190 reptile species in the United States and Canadian provinces where greater sage-grouse occur, 70 (37%) occur within the range of the bird. Of these 70 species, about a third (11 snake and 11 lizard species) have >10% of their distribution area within the sage-grouse range. Land cover similarity indices revealed that 14 of the 22 species (8 snake and 6 lizard species) had relatively similar land cover associations to those of sage-grouse, suggesting greater potential to be protected under the sage-grouse conservation umbrella and greater potential to be affected, either positively or negatively, by habitat management actions intended for sage-grouse. Conversely, the remaining 8 species are less likely to be protected because of less overlap with sage-grouse habitat and thus uncertain effects of sage-grouse habitat management actions. Our analyses of treatment databases indicated that from 1990 to 2014 there were at least 6,400 treatments implemented on public land that covered approximately 4 million ha within the range of the sage-grouse and, of that, >1.5 million ha were intended to at least partially benefit sage-grouse. Whereas our results suggest that conservation of intact sagebrush vegetation communities could benefit ≥14 reptiles, a greater number than previously estimated, additional research on each species' response to habitat restoration actions is needed to assess broader claims of multi-taxa benefits when it comes to manipulative sage-grouse habitat management. Published 2020. This article is a U.S. Government work and is in the public domain in the USA.     

Sharing of Diverse Mycorrhizal and Root-Endophytic Fungi among Plant Species in an Oak-Dominated Cool–Temperate Forest

Most terrestrial plants interact with diverse clades of mycorrhizal and root-endophytic fungi in their roots. Through belowground plant–fungal interactions, dominant plants can benefit by interacting with host-specific mutualistic fungi and proliferate in a community based on positive plant–mutualistic fungal feedback. On the other hand, subordinate plant species may persist in the community by sharing other sets (functional groups) of fungal symbionts with each other. Therefore, revealing how diverse clades of root-associated fungi are differentially hosted by dominant and subordinate plant species is essential for understanding plant community structure and dynamics. Based on 454-pyrosequencing, we determined the community composition of root-associated fungi on 36 co-occurring plant species in an oak-dominated forest in northern Japan and statistically evaluated the host preference phenotypes of diverse mycorrhizal and root-endophytic fungi. An analysis of 278 fungal taxa indicated that an ectomycorrhizal basidiomycete fungus in the genus Lactarius and a possibly endophytic ascomycete fungus in the order Helotiales significantly favored the dominant oak (Quercus) species. In contrast, arbuscular mycorrhizal fungi were generally shared among subordinate plant species. Although fungi with host preferences contributed to the compartmentalization of belowground plant–fungal associations, diverse clades of ectomycorrhizal fungi and possible root endophytes were associated not only with the dominant Quercus but also with the remaining plant species. Our findings suggest that dominant-ectomycorrhizal and subordinate plant species can host different subsets of root-associated fungi, and diverse clades of generalist fungi can counterbalance the compartmentalization of plant–fungal associations. Such insights into the overall structure of belowground plant–fungal associations will help us understand the mechanisms that facilitate the coexistence of plant species in natural communities.     

Ecological fitting by phenotypically flexible genotypes: implications for species associations, community assembly and evolution

Ecological fitting is the process whereby organisms colonize and persist in novel environments, use novel resources or form novel associations with other species as a result of the suites of traits that they carry at the time they encounter the novel condition. This paper has four major aims. First, we review the original concept of ecological fitting and relate it to the concept of exaptation and current ideas on the positive role of phenotypic plasticity in evolution. Second, we propose phenotypic plasticity, correlated trait evolution and phylogenetic conservatism as specific mechanisms behind ecological fitting. Third, we attempt to operationalize the concept of ecological fitting by providing explicit definitions for terms. From these definitions, we propose a simple conceptual model of ecological fitting. Using this model, we demonstrate the differences and similarities between ecological fitting and ecological resource tracking and illustrate the process in the context of species colonizing new areas and forming novel associations with other species. Finally, we discuss how ecological fitting can be both a precursor to evolutionary diversity or maintainer of evolutionary stasis, depending on conditions. We conclude that ecological fitting is an important concept for understanding topics ranging from the assembly of ecological communities and species associations, to biological invasions, to the evolution of biodiversity.     

Local similarity analysis reveals unique associations among marine bacterioplankton species and environmental factors

MOTIVATION: Characterizing the diversity of microbial communities and understanding the environmental factors that influence community diversity are central tenets of microbial ecology. The development and application of cultivation independent molecular tools has allowed for rapid surveying of microbial community composition at unprecedented resolutions and frequencies. There is a growing need to discern robust patterns and relationships within these datasets which provide insight into microbial ecology. Pearson correlation coefficient (PCC) analysis is commonly used for identifying the linear relationship between two species, or species and environmental factors. However, this approach may not be able to capture more complex interactions which occur in situ; thus, alternative analyses were explored. RESULTS: In this paper we introduced local similarity analysis (LSA), which is a technique that can identify more complex dependence associations among species as well as associations between species and environmental factors without requiring significant data reduction. To illustrate its capability of identifying relationships that may not otherwise be identified by PCC, we first applied LSA to simulated data. We then applied LSA to a marine microbial observatory dataset and identified unique, significant associations that were not detected by PCC analysis. LSA results, combined with results from PCC analysis were used to construct a theoretical ecological network which allows for easy visualization of the most significant associations. Biological implications of the significant associations detected by LSA were discussed. We also identified additional applications where LSA would be beneficial. AVAILABILITY: The algorithms are implemented in Splus/R and they are available upon request from the corresponding author.