PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ploscompComputational BiologyView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS Comput Biol. 2010 August; 6(8): e1000876.
Published online 2010 August 5. doi:  10.1371/journal.pcbi.1000876
PMCID: PMC2916855

Connectivity, Cycles, and Persistence Thresholds in Metapopulation Networks

Jordi Bascompte, Editor

Abstract

Synthesising the relationships between complexity, connectivity, and the stability of large biological systems has been a longstanding fundamental quest in theoretical biology and ecology. With the many exciting developments in modern network theory, interest in these issues has recently come to the forefront in a range of multidisciplinary areas. Here we outline a new theoretical analysis specifically relevant for the study of ecological metapopulations focusing primarily on marine systems, where subpopulations are generally connected via larval dispersal. Our work determines the qualitative and quantitative conditions by which dispersal and network structure control the persistence of a set of age-structured patch populations. Mathematical modelling combined with a graph theoretic analysis demonstrates that persistence depends crucially on the topology of cycles in the dispersal network which tend to enhance the effect of larvae “returning home.” Our method clarifies the impact directly due to network structure, but this almost by definition can only be achieved by examining the simplified case in which patches are identical; an assumption that we later relax. The methodology identifies critical migration routes, whose presence are vital to overall stability, and therefore should have high conservation priority. In contrast, “lonely links,” or links in the network that do not participate in a cyclical component, have no impact on persistence and thus have low conservation priority. A number of other intriguing criteria for persistence are derived. Our modelling framework reveals new insights regarding the determinants of persistence, stability, and thresholds in complex metapopulations. In particular, while theoretical arguments have, in the past, suggested that increasing connectivity is a destabilizing feature in complex systems, this is not evident in metapopulation networks where connectivity, cycles, coherency, and heterogeneity all tend to enhance persistence. The results should be of interest for many other scientific contexts that make use of network theory.

Author Summary

Taking advantage of modern network theory, we present a model formulation for determining those factors that control the stability and persistence of complex biological systems. As a case study, we focus on ecological metapopulations, which may be viewed as a set of distinct subpopulations (/sites) that are connected via a dispersal network of arbitrary complexity. Metapopulation persistence is found to depend critically on the topology of cycles, and cyclical components in the connectivity network, because they allow the offspring of the population to eventually “return home” to the sites from which they originated. The methodology identifies critical migration routes, whose presence are vital to overall stability, and are thus of high conservation priority – information that may be of value when designing networks of marine protected areas. In contrast, links that do not participate in a cyclical component have no impact on persistence and thus have low conservation priority. While network theory is highly fashionable in biology, only few studies go deeper than descriptive statistical applications as attempted here. Moreover, the key results are easily extended to other biological contexts (e.g., disease networks), particularly in situations whereby the network controls the dynamics of a complex system.

Introduction

Theoretical biologists and ecologists have long sought to understand the relationships between complexity, connectivity and the stability of large biological systems [1][10]. This interest has only grown in recent years particularly with the new developments in modern network theory and its multidisciplinary applications [11][15]. Here we outline a new synthesis relevant for ecological metapopulations, and provide a framework for untangling the role of dispersal and network structure in maintaining the persistence of a set of patch populations distributed in space. The metapopulation concept has become the theoretical framework that stands behind many modern conservation efforts [2][4], [16][29]. For example, the framework was adopted as part of the EU Habitats Directive Natura 2000 [20], which is the single most important legal tool for biodiversity conservation that has become binding national law of all European Member States. The Directive's goal is the creation of a coherent Europe-wide network of sites to protect important habitats and species. Ecological coherence “is seen in terms of the capacity for individual protected areas to support each other and in the interactions with habitat surrounding protected areas.” Species dispersal between sites provides supportive buffering for impacted habitats (e.g., oil spills) or allows for shifts in species ranges in the face of climate change, thus enhancing overall coherence. Similar metapopulation approaches have been made use of in the development of Marine Protected Areas (MPA's) – areas of the ocean protected from human disturbances – and they are currently strongly advocated as a tactical management tool [20][22].

Given the importance of the metapopulation concept for nature conservation, theoretical studies have attempted to gain insights into those factors that make for sustainable populations. Intriguingly, initial research indicated that spatial structure has in fact very little effect, with criteria for metapopulation stability appearing to be identical to the stability conditions for a single patch [23]. Several important studies have shown that variability or heterogeneity in patch connectivity may play a role in enhancing persistence [24], [25], but this work has not been developed further. Interest in these issues has resurfaced in recent years especially in the study of marine ecosystems where connectivity and the role of dispersal in maintaining persistent metapopulations is as controversial as it is enigmatic. Whether marine populations are retentive and recruit back to their native populations, or whether they are open and in the main, disperse with little self-recruitment, is an issue that remains unresolved despite decades of work. As yet, there is little understanding of how the architecture of marine networks control metapopulation persistence, and even fundamental concepts still remain controversial. Hastings and Botsford [3] concluded that multiple criteria are necessary to assess persistence, which explains why they were unable to “obtain a single number [criterion] like the reproductive number for a single population.” Here we derive simple general persistence criteria that are couched in terms of the metapopulation's reproductive potential, network connectivity as well as the topology of cycles in the dispersal network.

Similar to [3], we consider populations that comprise sedentary adults with dispersing juvenile stages, making the model appropriate for marine and terrestrial invertebrates, as well as many plants and fish. However, unlike most other studies, the model includes age-structure and may be generalized further to allow for dispersal across patches between the population's different age-classes. The model is particularly suited to marine systems where questions concerning connectedness, larvae retention and open versus closed populations have become highly topical in the last years.

We begin by considering a population that has m age-classes as described by the vector An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e001.jpg where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e002.jpg is the number of individuals in the k'th age class in year t. For such a population in a single isolated patch, growth may be modelled via the familiar Leslie matrix equations [31]:

equation image
(1)

and the time-step is one year. Here An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e004.jpg is the probability that an individual of age (k-1) years survives to age-k. A fraction An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e005.jpg of juveniles from the population successfully “self-recruit” and return to the population after the dispersal phase, which in marine settings might represent local larvae retention. The parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e006.jpg represents the fertility of age-class-k individuals, in terms of the average number of juveniles produced in the next generation. A graph theoretic interpretation of more complex stage-structured Leslie matrix population models may be found in [32].

The population has a single equilibrium, the extinction state An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e007.jpg whose stability is the basis for understanding population persistence in this model. It is a classical result [2], [31] that the average number of juveniles produced in the lifetime of a typical individual is given by the reproductive number:

equation image
(2)

The condition for population persistence (i.e., a growing population) requires that a typical individual is capable of replacing itself and give rise to at least a single offspring that successfully recruits back to the population. As only a proportion An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e009.jpg of juveniles successfully self-recruit, the net reproductive effort is An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e010.jpg. Thus a growing persisting population requires that the persistence parameter:

equation image
(3)

and the patch is said to be a “source.” Should An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e012.jpg, all age-classes approach a stable extinction state and the patch is said to be a “sink.” These definitions of source and sink follow those used in Armsworth [2] and are based on the original definitions of Pulliam [30], although it should be noted that other definitions are sometimes employed. In words, a sink would correspond to a situation where deaths and emigrations exceed the number of new juveniles [2]. A mathematical analysis of the model's single equilibrium, the extinction state An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e013.jpg, corroborates the above threshold criterion.

Now scaling up, consider a network of n age-structured patch-populations, where juveniles disperse between patches as portrayed in Figure 1. The metapopulation dynamics are given by:

equation image
(4)

In this notation An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e015.jpg is the m-dimensional age-class population vector at patch-i, and each patch has its own associated survival (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e016.jpg) and fertility (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e017.jpg) matrices. More specifically:

equation image
(5)

where the Leslie matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e019.jpg in Eqn.1 is retrieved taking An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e020.jpg. Dispersal processes between the n-patches are defined in the connectivity matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e021.jpg whose elements An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e022.jpg correspond to the proportion of juveniles produced on local population-j that are transported and successfully recruit to local population-i. We note that Eqns. 4 may be viewed as a good approximation to a more complex nonlinear metapopulation model (e.g., possibly with density-dependence), if the latter model is linearized about its extinction state. The asymptotic stability criterion for the extinction state is then the same for both the linearized and full nonlinear model.

Figure 1
The topology of the dispersal pathways between a set of patches as defined by the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e023.jpg.

Results

Connectivity thresholds

A key goal is to examine the effects solely of network architecture on metapopulation persistence. This requires studying the metapopulation with all processes being equal apart from the network structure. We thus initially suppose that local patch populations are identical, and all with the same survival and fertility rates An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e028.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e029.jpg, an assumption that is later relaxed. Based on the properties of the single patch dynamics noted above, it is possible to deduce (see Materials and Methods, as well as Text S1) that metapopulation persistence is controlled by a simple threshold condition: The metapopulation is growing if the fundamental persistence parameter

equation image
(6)

while the extinction state is stable if An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e031.jpg. Here An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e032.jpg is the spectral radius, or eigenvalue of largest magnitude, of the connectivity matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e033.jpg. As before, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e034.jpg is the reproductive number of a single patch as defined in Eqn. 2. The special role played by the matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e035.jpg allows us to view the condition as a connectivity threshold. Although previous studies of age-structured metapopulations [2], [4], [21] attempted to deduce persistence criteria, they overlooked the role of the spectral radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e036.jpg.

The metapopulation's complex patch dispersal structure may be visualized in terms of a network or graph, with nodes as patches and edges as dispersal links between patches. The network topology is summarised by the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e037.jpg whose elements An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e038.jpg if there is a direct dispersal route from patch-j to patch-i, while An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e039.jpg in the absence of such a route. In the formulation used here self-loops are excluded from the adjacency matrix, and we can set An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e040.jpg.

Again, in order to specifically elucidate the effects of network structure, it is necessary to ensure that the various processes between patches are kept equal. This motivates a relatively simple but nevertheless useful scheme in which the number of juveniles that immigrate to a patch population is on average a fixed proportion, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e041.jpg, of the source population from which they originate as shown in Figure 1. The average retentivity, or self-recruitment, of each patch is set at An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e042.jpg indicating the proportion of juveniles that complete their life-cycle in the patch. In this scheme highly connected patches recruit better, while less connected patches are disadvantaged with recruits being lost from the metapopulation as a result. Such losses are to be expected in a consistent model of heterogeneous dispersal.

The connectivity matrix may now be written in the particularly simple form: An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e043.jpg where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e044.jpg is the identity matrix. This immediately yields the key relationship whereby the persistence parameter (Eqn. 6) may be written as:

equation image
(7)

We see directly that persistence is controlled by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e046.jpg, the spectral radius of the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e047.jpg. We note that it is not always a simple matter to deduce the value of the spectral radius without resorting to a numerical study of the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e048.jpg adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e049.jpg. Factors that control the spectral radius include the number of nodes (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e050.jpg), the number of edges and the underlying structure of the matrix in terms of the topology of the dispersal routes. However, there are special cases that we will examine, where the particular structure of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e051.jpg makes it possible to solve for the spectral radius analytically.

The expression (Eqn. 7) for the persistence parameter shows clearly the importance of supply-side ecology with new recruits enhancing the possibility of persistence. A metapopulation, in which all n patch populations are sinks (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e052.jpg) requires sufficient subsidy recruitment (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e053.jpg) for persistence; enough that ensures An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e054.jpg. In such a case, the metapopulation is effectively “pulled up by its bootstraps” due to the inflow of larvae circulating through the complex dispersal routes of the network. The additional recruitment here acts as a rescue effect.

An exciting outcome of Eqn. 7 is that it opens the door for investigating the effects of network structure, at both fine and coarse-scale levels. Beginning with coarser-scale features, we focus on network topology. Figure 2 shows three simple network topologies: a) a regular network in which each patch has exactly two connections; b) a random network in which patches have very close to two connections; c) a random heterogeneous network in which patches have large variability in connectivity. For all three networks the average number of connections per patch, or mean degree, is deliberately held the same.

Figure 2
Persistence as a function of heterogeneity in node degree for three example networks.

It is then natural to ask which topology best enhances persistence and how does the topology's heterogeneity affect the threshold An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e058.jpg? This is far more than an academic question given that marine systems can be highly heterogeneous in terms of recruitment, sometimes with orders of magnitude variation to patches within the same metapopulation [33].

The three above cases are now treated separately assuming the networks are undirected; that is assuming a dispersal pathway from patch-i to patch-j implies a converse flow from patch-j to i. However, the results obtained are substantially similar with interpretations that broadly carry over to the case of directed networks [34].

a) Regular networks

In a regular network every patch is connected to An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e059.jpg other patches (see Figure 2a). Again we suppose each patch is recruiting equally from every other, but now with individual migration rates equal to self-recruitment rates An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e060.jpg. Since it is well known for regular networks that the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e061.jpg has spectral radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e062.jpg [34], then metapopulation persistence is ensured if An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e063.jpg. This is illuminating in that it shows that the gain of an individual patch derived from subsidy-recruitment (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e064.jpg) can be considerably larger than self-recruitment (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e065.jpg), and that this “bootstrapping” contribution may be of significance for crossing the persistence threshold.

b) Random networks

The classical Erdős–Rényi (ER) random network assumes there is a probability An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e066.jpg that any two patches -i and -j have a dispersal pathway between them (Figure 2b). As the network has An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e067.jpg patches, the mean number of dispersal pathways per patch is An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e068.jpg. For these random networks, it is well known that the spectral radius of the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e069.jpg is An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e070.jpg [34] and:

equation image
(8)

Thus increasing connectivity simply by the random addition of dispersal routes into a metapopulation increases the mean degree An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e072.jpg, and thus An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e073.jpg, and should therefore be seen as advantageous. The more random pathways, the larger is An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e074.jpg, and the more likely is the possibility of increasing beyond the persistence threshold. This result has important consequences for the ongoing debate as to whether highly connected ecological and biological systems are more stable. Unlike May's [1] prediction where connectivity is seen to be destabilizing in large complex systems, here increasing connectivity clearly enhances metapopulation persistence. It should however be noted that strictly speaking May was referring to species interactions (although see [29]), while our study concerns something quite different, namely network connectivity.

c) Heterogeneous networks

The result can be generalized for heterogeneous random networks having arbitrary degree distribution (Figure 2c). Suppose patch-i has An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e075.jpg dispersal connections to other patches. The variability of the different An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e076.jpg (i.e., the “degree distribution”) is a good index of the heterogeneity in connectivity of the network and is often measured by the squared coefficient of variation An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e077.jpg. Under these conditions,

equation image
(9)

Thus increasing the heterogeneity in connectivity (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e079.jpg) will enhance persistence (even if the mean degree An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e080.jpg remains unchanged) since it increases An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e081.jpg (see Figure 2d). Analogous results were found in [24][26] although not presented in the context of network heterogeneity as here. Curiously, large heterogeneity may be the norm in many metapopulations. Williams and Sale [33] found orders of magnitude variation in recruitment to patches of the same coral species in a single coral lagoon. In some cases substantial variation existed between patches 1km apart. Since then it has been shown that wildly fluctuating recruitment success is not particular to reef fish.

Cycles

Finer-scale network features also play an important role in determining persistence. The present framework allows exploration of how various network motifs [8], [13] or specific subgraphs of the network, might be influential. The following analysis of cyclical motifs indicates the power of this approach. We focus on cycles in directed networks where a dispersal route from patch-A to B (A→B) does not necessarily imply the converse. Typically a closed loop for three patches would have the following structure A→B→C→A, and is termed a directed cycle.

C1) A metapopulation network without any cycles is unable to persist

Firstly, it becomes trivial to show that a metapopulation without any cycles in the dispersal network is unable to persist, as might be expected (see also [2]). In such a metapopulation network, all juveniles or any of their eventual descendants are unable to recruit back to their patch of origin – they never “return home” [3]. When there are no cycles, the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e082.jpg must have spectral radius An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e083.jpg [35] and in the absence of self-recruitment (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e084.jpg), the persistence parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e085.jpg. A stable extinction state is expected.

C2) For metapopulations with self-recruitment but without any other cycles, there is no advantage to dispersal

Self-recruitment An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e086.jpg implies that some proportion of juveniles complete their life-cycle in the patch. As there are no other cycles, all other juveniles fail to return home and the spectral radius remains An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e087.jpg and An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e088.jpg. Thus the criterion for the entire metapopulation to persist, An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e089.jpg, is precisely the same as the criterion for a single self-recruiting patch to persist (Eqn. 3). Consequently there is no advantage for dispersal if there are no cycles and larvae fail to “return home.”

C3) “Lonely links,” those links in a network that are not part of a cycle, have no effect whatsoever on the metapopulation's persistence characteristics

Figure 3a illustrates a complex network with a single simple cycle (red), or cyclic component (see Text S2 as well). Eliminating all lonely links that are not connected to the cyclic component results in the network shown in Figure 3b. As proven for the general case (see Text S2), the persistence parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e090.jpg associated with the two networks is identical and hence both have the same persistence/stability properties. In short, the removal of “lonely links” has no effect on metapopulation persistence making them of low conservation priority.

Figure 3
Example of removing lonely links.

C4) The persistence parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e091.jpg increases with the complexity and richness of the network's cycle structure

Consider first simple cycles, with no sites receiving recruits directly from two distinct cycles. Directed networks composed of only simple cycles, have the property that An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e092.jpg and thus An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e093.jpg, regardless of how many simple cycles are present. If cycles intersect they are no longer simple and in general An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e094.jpg, and thus An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e095.jpg, increases as the number of overlapping cycles increases. Intuitively one notes that maximal “cycle packing” might be considered analogous to a regular network with all n-patches connected to each other giving An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e096.jpg. Thus metapopulation persistence is greatly enhanced when there is intense overlapping of dispersal routes between sites in the network.

C5) Networks may be broken down into nonoverlapping cyclic components

Persistence, is completely determined by the dominant cyclic component. Intriguingly it turns out that valuable information regarding persistence may be obtained by breaking down a network into its nonoverlapping cyclic components. In particular, we are able to show (Text S2) that persistence is completely determined by the dominant cyclic component, that is, the component whose associated eigenvalue is largest in magnitude. The complex network shown in Figure 3c illustrates this concept. Stripping away all lonely links, reveals that the network has two cyclic components (encircled). Extracting the larger eigenvalue associated with the more complex cyclic component, gives us An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e097.jpg and persistence may be checked using [Eqn. 4] (see Text S2). Thus the dominant cyclic component gives complete information concerning the prospects of metapopulation persistence. Note that while “lonely links” have no effect on metapopulation persistence, the example in Figure 4 shows that the existence of a single critical link may in some cases determine the fate of the entire system. Thus removal of the dispersal route (see arrow pointing to critical link in Figure 4) in the dominant cyclic component can reduce An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e098.jpg dramatically, and persistence will be lost. More specifically, removal of the critical link indicated reduces the spectral radius of An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e099.jpg by approximately 33% (from ~1.5 to 1) and thus reduces the persistence parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e100.jpg considerably possibly resulting in extinction. As such this route is of considerable importance from the perspective of conservation.

Figure 4
Persistence is controlled by the network's cyclical components of which there are two here (red patches); one simple component (3 patches, 3 links) and one complex component (6 patches, 9 links) formed by intersecting cycles.

C6) Symmetry

It is revealing to return to the earlier examples of Figure 2a–c that deal exclusively with undirected symmetric networks. As dispersal is often viewed as a function of distance rather than direction, symmetric connectivity is a common assumption in metapopulation models [36]. Yet recent research has revealed that ecological and metapopulation networks are often asymmetric, [7], [37] and according to simulation studies in [7], [36][38], this should have a negative effect on persistence. The framework advanced here provides a simple theoretical explanation, based on the elementary observation that symmetric networks tend to have a propensity of cycles. Even a single connection between two patches in a symmetric network is bidirectional (A→B and B→A) and thus forms a cycle. Increasing asymmetry, tends to remove cycles, and by property C4 reduces An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e109.jpg, thereby having a negative effect on persistence. Text S2 delves into these properties further.

Extensions: Nonidentical patches

So far, persistence has been discussed in the context of networks of identical patch populations. However this work may be extended to obtain persistence/extinction criteria for metapopulations comprised of nonidentical patch populations thus accommodating cases where age structure varies between patches. For convenience, suppose patch-1 of the metapopulation “dominates” in the sense that it has the highest fertility parameters (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e110.jpg for all An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e111.jpg) and highest survival rates (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e112.jpg for all An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e113.jpg) (see Text S3), and thus the highest net reproductive effort (An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e114.jpg). Now construct a metapopulation comprised of n-identical copies of dominating patch-1, retaining the original directed or undirected network structure. Test whether the extinction state of these n-identical patches is stable. If so, it can be shown that the original heterogeneous metapopulation of n-nonidentical patches is also unable to persist (see Text S3).

Along similar lines, suppose patch-1 is the weakest in the sense that it has the lowest fertility parameters and lowest survival rates. If the network of n-identical copies of patch-1 is able to persist, then the original heterogeneous metapopulation must also persist. Thus, the dominant and subordinate patch population of a network, may be used as a guide for determining bounds for the respective extinction or persistence thresholds of the larger heterogeneous metapopulation. These effects are illustrated in Figure 5 and discussed in Text S3.

Figure 5
Metapopulation of four nonidentical patch populations in an “all to all” network.

It should be noted that the effects of age-structure become prominent with more complex migration schemes, for example, between different age-classes from different patches, and representations of ontogenetic shifts in habitat use that are life history dependent [[21]; and see Text S3].

Discussion

Our work has shown the importance of the persistence parameter An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e117.jpg in assessing the fate of the metapopulation. As has been emphasized, a persisting metapopulation is characterized by a state of growth with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e118.jpg, and implies that at least one patch population, say patch-1 of Figure 5, is increasing in numbers with time. This raises the question as to which other patches in the metapopulation persist, and which go extinct? Given that patch-1 increases in time, it must also continuously export larvae to all other patches it connects to, either directly or indirectly. Thus every patch that can be reached by patch-1 via a chain of dispersal pathways will continuously receive new recruits, and must therefore also persist, irrespective of whether it is a sink or source. In contrast, patch populations that fail to receive new recruits because they are disconnected from other patches in the network, cannot persist if they are sinks. This has the important implication that the underlying architecture of the metapopulation network controls the particular set of patch populations that survive (see Figure 4).

Ideally, we also require conditions that establish metapopulation coherency, whereby not only subsets of patches persist, but they are also connected with one another directly or indirectly in a large-scale manner that connects together the majority of the metapopulation network. Coherency offers the metapopulation many advantages. Large scale or global dispersal between sites buffers population fluctuations [39], allows recolonization after local isolated environmental impacts (e.g., oil spills), allows species ranges to shift in the face of climate change [20] and more generally allows “spreading of risks”; properties which are all verified from stochastic metapopulation simulations. Note that it is assumed that dispersal is not strong enough to induce large-scale spatial synchronization, a phenomenon that is prone to enhancing extinction risk [40].

The random Erdős–Rényi (ER) model which exemplifies a broad range of complex networks shows the existence of important coherency thresholds that must be taken into account. Recall that there is a probability p that any two patches-i and j have a dispersal pathway between them. Moreover, there is a critical value An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e119.jpg, such that for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e120.jpg the great majority of the patches in the network become connected to one another by a pathway of edges forming the celebrated “Giant Component” in which there is network-wide connectance [14], [15]. Hence An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e121.jpg is a necessary condition for a coherent (ER) metapopulation. Conversely for An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e122.jpg, the majority of patches are disconnected. In the latter case, metapopulation persistence is impossible in such poorly connected incoherent network. Although only a simplified system, the ER model makes clear the importance of coherency thresholds and serves to further strengthen the philosophy underlying the EU Habitats Directive for the creation of a coherent Europe-wide network of sites to protect important habitats and species [20].

Taken altogether, the above framework has revealed a number of new insights regarding the determinants of persistence, stability and ecological thresholds in complex metapopulations. We note that a number of results derived here are based on several simplifying assumptions. Thus, at least in the first instance, it was necessary to assume that all patches are identical in terms of their structure and dynamics (a condition that was later relaxed). In addition, it was assumed that the rate of dispersal was the same for all pairs of connected patches. In other applications these model limitations may need to be reckoned with. However, for our purposes, only by fixing patches and metapopulation processes equal, is it possible to isolate exactly how different dispersal networks and their topologies govern persistence. This is in fact a major difference between our work here and that of Armsworth [2]. Although Armsworth also found condition C1 above for a general age-structured metapopulation, his analysis did not result in a quantitative formulation connecting the adjacency matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e123.jpg, its spectral radius, and its topological features including especially cycles in the network, and their relation to the persistence threshold.

Our analysis was couched both in terms of coarse and fine-scale network features. With regard to the former, at least for random systems, connectivity and heterogeneity in connectivity were found to be two key factors that enhance metapopulation persistence. In terms of more fine-scale features, the existence of critical dispersal routes have been identified, while cycles have been shown to play a prominent functional role, allowing metapopulations to bootstrap themselves into persistence. All of these factors have conservation applications and would translate most readily into principles that aid in the management of marine populations and in the design of networks such as marine protected areas [20][22]. In addition these outcomes have obvious applications that cross over into many other scientific contexts including the dynamics of disease epidemics that spread via complex large-scale travel networks, internet traffic and computer virus dynamics or metabolic network analyses [10][14].

Materials and Methods

For a metapopulation of n identical patches, the i'th patch population may be notated as An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e124.jpg where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e125.jpg is the number of individuals in the k'th age class in year t. The model of n-identical age-structured patch populations An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e126.jpg may be written using Kronecker product matrix notation. Let An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e127.jpg, then:

equation image
(10)

where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e129.jpg is the An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e130.jpg identity matrix , An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e131.jpg is the lower diagonal matrix with An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e132.jpg and otherwise zero. An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e133.jpg has entries An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e134.jpg and is otherwise zero. (Note that An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e135.jpg.) The above model might also be viewed as a first order approximation to a more general nonlinear metapopulation about the extinction state. After appropriate matrix manipulation, Text S1 shows that the stability of the extinction state depends on the eigenvalues of the stability matrix

equation image
(11)

where An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e137.jpg is the spectral radius of the connectivity matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e138.jpg. Due to the correspondence between the metapopulation stability matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e139.jpg and the Leslie matrix An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e140.jpg above, the metapopulation (Eqn. 11) is persistent and growing if An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e141.jpg, or equivalently An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e142.jpg while the extinction state is stable if An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e143.jpg. While the threshold is reminiscent of results found for unstructured metapopulations [24], [25] however the inclusion of age-structure and/or network structure makes this is a nontrivial result and a challenge set in Refs. [2], (see Text S1.)

For the particular case An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e144.jpg, the spectral radius is given by An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e145.jpg so that An external file that holds a picture, illustration, etc.
Object name is pcbi.1000876.e146.jpg.

Supporting Information

Text S1

The persistence parameter χ. Formulation of age-structured metapopulation model and its mathematical stability criteria.

(0.74 MB PDF)

Text S2

Cycles. Mathematical analysis of the impact of cycles on network stability.

(1.52 MB PDF)

Text S3

Nonidentical patches. Analysis of the metapopulation model for nonidentical patches.

(0.38 MB PDF)

Footnotes

The authors have declared that no competing interests exist.

Lewi Stone is grateful for the support of the Israel Ministry of Sciences and the European Union Future Emerging Technologies contract no. 231807 (EPIWORK). Yael Artzy-Randrup is grateful for support from the Howard Hughes Medical Institute, Rothschild Post-Doctorial Fellowship and the Sara Lee Schupf Post-Doctorial Award funded by the Clore Foundation and S. Donald Sussman. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. May RM. Will a large complex system be stable? Nature. 1972;283:413–414. [PubMed]
2. Armsworth PR. Recruitment limitation, population regulation, and larval connectivity in reef fish metapopulations. Ecology. 2002;83:1092–1104.
3. Hastings A, Botsford LW. Persistence of spatial populations depends on returning home. Proc Natl Acad Sci USA. 2006;103:6067–6072. [PubMed]
4. Fahrig L, Merriam G. Habitat patch connectivity and population survival. Ecology. 1985;66:1762–1768.
5. Holland MD, Hastings A. Strong effect of dispersal network structure on ecological dynamics. Nature. 2008;456:792–794. [PubMed]
6. Bastolla U, Fortuna MA, Pascual-García A, Ferrera A, Luque B, et al. The architecture of mutualistic networks minimizes competition and increases biodiversity. Nature. 2009;458:1018–1020. [PubMed]
7. Bascompte J, Jordano P, Olesen JM. Asymmetric coevolutionary networks facilitate biodiversity maintenance. Science. 2006;312:431–433. [PubMed]
8. Bascompte J, Stouffer DB. The assembly and disassembly of ecological networks. Phil Trans R Soc B. 2009;364:1781–1787. [PMC free article] [PubMed]
9. Montoya JM, Pimm SL, Ricard V, Sole RV. Ecological networks and their fragility. Nature. 2006;442:259–264. [PubMed]
10. Saavedra S, Reed-Tsochas F, Uzzi B. A simple model of bipartite cooperation for ecological and organizational networks. Nature. 2009;457:463–466. [PubMed]
11. Colizza V, Vespignani A. Invasion threshold in heterogeneous metapopulation networks. Phys Rev Lett. 2007;99(14) Article Number:148701. [PubMed]
12. Clauset A, Moore C, Newman MEJ. Hierarchical structure and the prediction of missing links in networks. Nature. 2008;453:98–101. [PubMed]
13. Artzy-Randrup Y, Fleishman SJ, Ben-Tal N, Stone L. Comment on “Network Motifs: Simple Building Blocks of Complex Networks.” Science. 2004;305:1107c. [PubMed]
14. Barabási AL. 2002. Linked: How Everything Is Connected to Everything Else, ISBN 0-452-28439-2.
15. Newman MEJ. The structure and function of complex networks. SIAM Review. 2003;45:167–256.
16. Caley MJ, Carr MH, Hixon MA, Hughes TP, Jones GP, et al. Recruitment and the local dynamics of open marine populations. Annu Rev Ecol Syst. 1996;27:477–500.
17. Cowen RK, Sponaugle S. Larval Dispersal and Marine Population Connectivity. Annu Rev Mar Sci. 2009;1:443–66. [PubMed]
18. Cowen RK, Gawarkiewicz G, Pineda J, Thorrold SR, Werner FE. Population connectivity in marine systems: An overview. Oceanography. 2007;20:14–20.
19. Gaines SD, Gaylord B, Gerber LR, Hastings A, Kinlan BP. Connecting places: The ecological consequences of dispersal in the sea. Oceanography. 2007;20:90–99.
20. Johnson MP, Crowe TP, McAllen R, Allcock AL. Characterizing the marine Natura 2000 network for the Atlantic region. Aquatic Conserv Mar Freshw Ecosyst. 2008;18:86–97.
21. Gerber LR, Botsford LW, Hastings A, Possingham HP, Gaines SD, et al. Population models for marine reserve design: a retrospective and prospective synthesis. Ecol App. 2003;13:S47–S64.
22. Fogarty MJ, Botsford LW. Population connectivity and spatial management of marine fisheries. Oceanography. 2007;20:112–123.
23. Rohani P, May RM, Hassell MP. Metapopulations and Equilibrium Stability: The Effects of Spatial Structure. J Theor Biol. 1996;181:97–109. [PubMed]
24. Adler FR, Neurnberger B. Persistence in patchy irregular landscapes. Theor Pop Biol. 1994;45:41–75.
25. Ovaskainen O, Hanski I. Spatially structured metapopulation models: Global and local assessment of metapopulation capacity. Theor Pop Biol. 2001;60:281–302. [PubMed]
26. Ovaskainen O, Hanski I. Extinction threshold in metapopulation models. Annales Zoologici Fennici. 2003;40:81–97.
27. Hanski I, Ovaskainen O. Metapopulation theory for fragmented landscapes. Theor Pop Biol. 2003;64:119–127. [PubMed]
28. Ovaskainen O, Sato K, Bascompte J, Hanski I. Metapopulation models for extinction threshold in spatially correlated landscapes. J Theor Biol. 2002;215:95–108. [PubMed]
29. Gardner MR, Ashby WR. Connectance of large dynamic (cybernetic) systems: Critical values for stability. Nature. 1970;228:784–784. [PubMed]
30. Pulliam HR. Sources, sinks and population regulation. Am Nat. 1988;132:652–661.
31. Caswell H. Matrix population models: Construction, analysis, and interpretation. 1989. Sinauer, Sunderland, Mass.
32. de-Camino-Beck T, Lewis MA. A new method for calculating net reproductive rate from graph reduction with applications to the control of invasive species. Bull Math Biol. 2007;69:1341–1354. [PubMed]
33. Williams D, Sale PF. Spatial and temporal patterns of recruitment of juvenile coral reef fishes to coral habitats within One Tree Lagoon, Great Barrier Reef. Mar Biol. 1981;65:245–253.
34. Restrepo JG, Ott E, Hunt BR. Approximating the largest eigenvalue of network adjacency matrices. Phys Rev E. 2007;76:056119. [PubMed]
35. McKay BD. Acyclic digraphs and eigenvalues of (0,1) matrices. Jnl. of Integer Sequences. 2004;7:1–5.
36. Bode M, Burrage K, Possingham HP. Using complex network metrics to predict the persistence of metapopulations with asymmetric connectivity patterns. Ecol Modell. 2004;214:201–209.
37. Vuilleumier S, Possingham HP. Does colonization asymmetry matter in metapopulations? Proc R Soc B. 2006;273:1637–1642. [PMC free article] [PubMed]
38. Salomon Y, Connolly S, Bode L. Effects of asymmetric dispersal on the coexistence of competing species. Ecol Lett. 2010;13:409–541. [PubMed]
39. Eckert GL. Effects of the planktonic period on marine population fluctuations. Ecology. 2003;84:372–383.
40. Earn DJD, Levin SA, Rohani P. Coherence and conservation. Science. 2000;290:1360–1364. [PubMed]

Articles from PLoS Computational Biology are provided here courtesy of Public Library of Science