Travelling in time with networks: Revealing present day hybridization versus ancestral polymorphism between two species of brown algae, Fucus vesiculosus and F. spiralis
© Moalic et al; licensee BioMed Central Ltd. 2011
Received: 17 September 2010
Accepted: 31 January 2011
Published: 31 January 2011
Hybridization or divergence between sympatric sister species provides a natural laboratory to study speciation processes. The shared polymorphism in sister species may either be ancestral or derive from hybridization, and the accuracy of analytic methods used thus far to derive convincing evidence for the occurrence of present day hybridization is largely debated.
Here we propose the application of network analysis to test for the occurrence of present day hybridization between the two species of brown algae Fucus spiralis and F. vesiculosus. Individual-centered networks were analyzed on the basis of microsatellite genotypes from North Africa to the Pacific American coast, through the North Atlantic. Two genetic distances integrating different time steps were used, the Rozenfeld (RD; based on alleles divergence) and the Shared Allele (SAD; based on alleles identity) distances. A diagnostic level of genotype divergence and clustering of individuals from each species was obtained through RD while screening for exchanges through putative hybridization was facilitated using SAD. Intermediate individuals linking both clusters on the RD network were those sampled at the limits of the sympatric zone in Northwest Iberia.
These results suggesting rare hybridization were confirmed by simulation of hybrids and F2 with directed backcrosses. Comparison with the Bayesian method STRUCTURE confirmed the usefulness of both approaches and emphasized the reliability of network analysis to unravel and study hybridization
Speciation is a central process in evolution, but its complexity and duration render it difficult or impossible to observe and study as a whole. Studies dealing with the functional and genetic divergence between taxa, particularly when reproductive isolation is incomplete , allow the distinction and analysis of the various stages of this process. The relative influence of the different mechanisms involved in the initiation and maintenance of divergence can then be inferred. Understanding the complexity of evolutionary and ecological mechanisms leading to reproductive isolation and speciation through integration of in situ observations, theoretical models and molecular analysis of genomes, is today one of the major challenges in evolutionary ecology [2–6]. Important efforts towards modeling the processes acting in hybrid zones and understanding the maintenance of divergence have focused on the balance between dispersal and hybrid depression . In these studies, incompatibility resulting from the allopatric divergence between two genomes is considered as a predominant factor causing hybrid depression. During the last decades, the development of the use of molecular markers in a population genetics framework made possible the testing and improvement of these theoretical models in hybrid zones. This screening of genome divergence and incompatibility also allowed testing the schematic models describing speciation as the result of processes spanning from pure vicariance (allopatry) to differential level of gene flow (sympatric speciation, [8–11]).
From a statistical point of view however, analyses at both the scales of populations and genes are still limited by the panel of tools available. At the population level, the analysis of molecular data is limited by the fact that most mathematical models underlying classical population genetics analysis have been developed in an intra-specific framework, and some underlying hypotheses such as random mating or Hardy Weinberg equilibrium do not necessarily stand in real natural populations and are out of scope in hybrid zones. Moreover, the detection and screening of hybrids requires an analysis of the populations and hybrid zones at the individual level whereas most summary statistics deliver estimates at the population level.
Two families of analyses have been proposed recently in order to partially release the underlying assumptions of most classical population genetics data analyses, and take better advantage of the information contained in a given dataset. The first family is built around the coalescent theory and is mostly based on Bayesian computation methods of analysis. It allows individual-centered analyses of genotype clustering [12–16]. The second family is built around the network theory and proposes an exploratory approach centered on the population or individual (also called agents), illustrating the relationship among those agents and inferring their respective roles and importance in the studied system [17–20].
Methods based on the coalescent theory trace the ancestral genealogy of a sample rather than modeling changes of gene frequencies in the population as a whole. These methods have recently been explored in order to improve mathematical models and take advantage of all the information contained in the data [12, 21, 22]. They have proven useful in studying hybrid zones  although they remain complex and time-consuming. On the other hand, this combination between high load of information not being optimally exploited by the summary statistics, and statistical tools available requiring heavy underlying hypotheses, suggests the use of another possible family of methods, coming from science of complex systems based on network theory [24, 25]. Network tools have indeed been developed to take the best advantage of the information contained in a complex data set, while minimizing the assumptions required for the analysis and interpretation of the complex system behavior .
A first step towards the use of the network theory to unravel gene flow has recently been proposed for inter-population analyses by Dyer and Nason . They applied the graph theory developed by Erdös and Renyi [1960, in 27] to describe the complex topology resulting from both history and contemporary genetic interactions among populations of a widely distributed species of cactus. This was further improved to extract from the graph topology hints as the dynamics of information (here gene) flow through the system at the inter-individual  and inter-population levels. Compared to classical methods , these first steps proved useful in illustrating the genetic relationship, and identifying clusters of individuals as well as agents acting as preferential links or sources in the metapopulation system of a threatened seagrass.
In the present study, we explored the usefulness of network analysis to study hybrids identified as agents linking two differentiated clusters of individuals, or gene pools. We tested network analysis in the hybrid zone between our two model species, the two sister species of brown seaweeds Fucus vesiculosus (F_ves) and F. spiralis (F_spi). Both species are ecologically successful and widely distributed (North Atlantic, Channel and North Sea shores). They are characteristic of respectively upper and mid-shore zones on rocky shores. Although distinct genetic entities have been identified within the species F. spiralis [F. spiralis-High and F. spiralis-Low, ] these are nevertheless still one single monophyletic entity, here designated as F. spiralis, distinct from its sister F. vesiculosus . These two sister taxa are a good model system because, despite displaying diagnostic reproductive system (F_ves is dioecious whereas F_spi is hermaphroditic), they may hybridize when encountered in sympatry [31–34], resulting in individuals with intermediate genotypes whose fitness in terms of reproductive investment is not significantly reduced . Whether the occurrence of hybridization is the result of an ongoing and incomplete speciation process or of a secondary contact with ongoing re-homogenization of the two entities is still unclear and a matter of debate [36, 37]. In order to test for the hypothesis of ongoing speciation (i.e., ancestral shared polymorphism) versus present day hybridization due to secondary contact (introgression), we studied the genetic relationships among samples collected both in sympatric and allopatric regions of F_ves and F_spi ranges (see FigureA1 in additional file 1). The global pattern of genetic relationships among individuals was illustrated by networks built with two different metrics integrating different genetic information in term of time and divergence history: the Rozenfeld distance (RD)  and the Shared Allele distance (SAD) . RD, based on loci, helps to resolve ancestral polymorphism through allele length impinged by slow evolutionary processes, while SAD, based on shared alleles, helps to understand recent gene flow characterized by direct allelic exchange free from slow evolutionary process. We aimed at comparing the performance of network and Bayesian coalescent analysis (STRUCTURE, ) to test whether the shared polymorphism and the absence of diagnostic locus was mainly due to retained ancestral polymorphism, or to present day hybridization in sympatric zone. In the first case, we expected shared polymorphism at neutral microsatellites to be distributed evenly across the species distribution range, whereas in the other case a higher proximity among species may be expected in sympatric zones.
Networks of individuals F. spiralis and F. vesiculosus
The individuals of F_spi that are sharing the higher number of links with individuals F_ves were also coming from Northwest Iberia. Despite the global differences denoted for our networks, this last observation is common to both the SAD and RD network topologies. Indeed, both highlight individuals of F_spi and of F_ves from Northwest Iberia relaying gene flow among their two clusters. The combination of our two network analyses with different distances indicates hybridization in our dataset at the Northwest Iberia area.
Hybridization assessment through network analysis
Hybrids datasets used in this study
number of simulated hybrids
number of natural individuals
total number of individuals
Hybridization assessment through a Bayesian clustering approach
When analyzing the datasets where putative natural hybrids were replaced by synthetic hybrids, admixture is mainly recovered. As with network analysis, all individuals of Hybrids F1 are clearly recognized as such by STRUCTURE (Figure 5C). This result is coherent with the network topology, as for backcrosses BC_F_spi, BC_F_ves and BC_F_spi_F_ves, 15/17, 11/17 and 14/17 (9+5) individuals are detected, respectively (Figure 5D, E and 5F). Nevertheless, it should be noticed that individuals BC_F_ves are less detected by STRUCTURE as admixed than BC_F_spi (9 vs. 2 on the totality of backcrosses). This result is different from the network analyses because individuals BC_F_spi tend to be closer to natural individuals of F_spi than BC_F_ves to natural individuals of F_ves.
In this study, we developed a novel application of network analysis to the study of hybridization phenomena between sister species. Networks were constructed based on distinct genetic distances (RD and SAD), differently sensitive to time since divergence (respectively based on allele lengths or on shared alleles), allowing the comparison of the divergence degree at different time-scales.
The first important picture obtained on network built with Rozenfeld Distance (RD) is the straightforward recognition of two well-defined clusters of F. vesiculosus and F. spiralis (Figure 1). This confirms that reproductive isolation and genetic divergence is rather advanced between those two species, despite a significant amount of shared polymorphism. Thus this multi-locus phylogenetic distance may be used to assign individuals to their species of origin. This result also supports the existence of phylogenetic information from microsatellite loci  which is accurately delivered by RD.
The second important result is the occurrence of clusters of individuals exhibiting an intermediate position between both species and maintaining a connection between some F. vesiculosus populations and the F. spiralis cluster as revealed on the SAD network (Figures 2 and 3). This intermediate cluster is formed by tens of individuals, among which those pointed out as intermediate with the RD. In case of the shared polymorphism being mostly ancestral, one may expect to observe some connections anywhere in the distribution range including in allopatric zones, but these F. vesiculosus individuals of intermediate genotypes are only detected in the Northwest Iberia, located specifically at the edge of the locations where both species occur in sympatry. Interestingly, intermediate individuals, although much less numerous, also emerge in two other sympatric areas: Northeastern America and the English Channel, whereas no such individuals are observed in any of the allopatric zones, further supporting the hybridization hypothesis. The SAD distance integrates more recent history, therefore giving it more weight than the RD, which takes into account phylogenetic divergence [19, 38]. It is therefore likely to better reflect present day exchange of genes between the two species. Those intermediate individuals that did cluster with the main F. vesiculosus cluster with RD, with SAD do now connect closer to F. spiralis than to the other con-specifics. This is likely caused by these individuals sharing more alleles with F. spiralis than with the other F. vesiculosus (therefore the branching with SAD distance), but the alleles they do not share with F. spiralis exhibit a strong divergence with this last species and are typical from F. vesiculosus (therefore the clustering of those individuals with F. vesiculosus with RD). The occurrence of all those individuals in the sympatric zone supports these intermediate genotypes as the product of present day hybridization between anciently diverged lineages/species.
The four additional networks built including simulated synthetic hybrids from F1, as well as backcrosses with each of the mother species (Figure 3), confirm this scenario. Indeed, the synthetic hybrids are present at the interface of the two sister-species, although they tend to be more isolated than most of the natural putative hybrids. This isolation seems due to the F1 nature of synthetic hybrids, truly half F. spiralis/half F. vesiculosus, while naturally occurring hybrids may be the product of backcrosses. When the backcrosses are plugged into the network, all possible backcrosses do not lead to a convergence towards the initial topology based on real data only. The admixture of (F1 X F. spiralis) and (F1 X F. vesiculosus) in the same dataset show the most similar network topology to the original one, while backcrosses resulting from (F1 X F. spiralis) or (F1 X F. vesiculosus) plugged alone tend to cluster closer to their species of origin without fitting exactly the same position as the natural intermediates.
The use of network here provides specific advantages, mainly the ability to identify agents (nodes) connecting identified clusters through genetic distance (links). Hybridization indeed usually results in a reticulate tree that cannot be built with classical phylogenetic methods but can be easily grasped on a network. Moreover, the possibility to follow the evolution of network topology through a gradient of genetic distance threshold helps to visualize and understand the attachment preferences to one or the other of the species clusters.
The use of STRUCTURE also illustrates the clear separation between the two species and allows the identification of admixed individuals, although not systematically the same as the intermediate agents appearing on the networks. In order to test for the accuracy of each method in those doubtful cases, the discrepancies between results obtained with both methods can be considered in perspective with geographical locations and ecological conditions in which intermediate agents (network) do not appear as admixed individuals (STRUCTURE), or the other way round. Interestingly, intermediate agents unidentified by STRUCTURE as admixed individuals correspond to samples collected at the edge of the zone where the two species occur in sympatry, Northwest Iberia. On the contrary, the admixed individuals detected by STRUCTURE that do not appear as intermediate agents were collected in the allopatric zone, thereby rendering enigmatic the origin of such admixture.
When analyzing the accuracy of both methods with synthetic hybrids and backcrosses, STRUCTURE shows high reliability in detecting F1 hybrids (100%), whereas only 71% of those are emphasized by an intermediate position on network. The detection of mixed individuals with STRUCTURE logically decreases to 92% for backcrosses with F. spiralis and to 64% for backcrosses with F. vesiculosus. This difference between the two species can be explained by the lower genetic diversity of F_spi, likely due to its reproductive mode, as illustrated by their GDS (see Figure A6 in additional file 7), resulting in an easier detection of the insertion of new alleles in the genetic pool. On the opposite, network analysis, by relying on shared links, shows backcrosses with F_spi having a stronger assimilation to its species of origin than backcrosses with F_ves. The integration of backcross events into the species gene pool seems to be rapid and renders them hard to identify. As a synthesis, both approaches appear as complementary, as STRUCTURE may perform slightly better in the systematic detection of admixed individuals, an advantage however balanced by a lower amount of misleading detection of hybrids (i.e. type I errors) obtained with network analysis. Besides, network analysis allows further screening of links and connections among geographic areas to describe patterns of spatial connectivity.
At the within species scale, finally, network analysis revealed two main genetically distinct clusters across the distributional range of both sister species F. spiralis and F. vesiculosus, (a) a Southern cluster in South Portugal, Morocco, Azores and Canary Islands, and (b) and a Northern cluster in North America, North Sea and Channel (see Figures A2 B and A4 B in additional files 2 and 4). These two regions are common to both species indicating a similar evolutionary history. Indeed, oscillations of climate during the past thousands of years have caused repeated geographic distributional shifts and extinction/recolonization events often experienced by many marine taxa. During the Last Glacial Maximum (LGM, 23-18 ka) in Europe, permafrost extended south at 47° N , leaving temperate species to shift their distribution to potential refugia. Marine species, including intertidal taxa, are also thought to have been displaced to small refugia in the North around the British Isles, Norway or in the Brittany region and more Southerly from the Iberian Peninsula to Mauritania . As the ice melted, species ranges were able to expand back to previous latitudes [42, 43]. On the North American Coast, the LGM may have covered by ice the complete hard substrate available there  and have caused the extinction of rocky shore species in this region. North American rocky shores are thought to have later been recolonized by European populations . For species in the genus Fucus, glacial refugia have been inferred along the Brittany area for F. serratus  or along Northwest Iberia for F. ceranoides  and F. vesiculosus along the American coast shows a genetic signature of a recent recolonization from Europe . In F. spiralis and F. vesiculosus, clustering on our networks is consistent with the presence of two refuge areas in Europe, possibly a Northern refuge (possibly along the Brittany or North Iberian region) and a southern one (possibly located along the southern Iberian Peninsula and/or North Africa). These inferences are thus consistent with a recent mitochondrial phylogeography of these two species for F. spiralis but are contradictory for F. vesiculosus, as the mitochondrial genome of F. vesiculosus supports a single refugial zone . Organelle genomes, with smaller effective population size than nuclear ones, are however more prone to being highly affected by introgressive sweeps, and indeed introgression and massive expansion of F. vesiculosus organelles into other Fucus species has been documented [32, 46]. Isolation into distinct glacial refugia would have been followed by a post-glacial expansion during which the Northern and Southern populations might have converged in a contact region along Northwest Iberia (Northern Portugal/Northwestern Spain) where despite hybridization genetic differentiation is still maintained nowadays (see Figures A2 B and A4 B in additional files 2 and 4) and which now also extends into the Brittany region [29, 30, 32], although in regions not included in our sampling. In the southern region, F. spiralis occurs only on the open coast whereas F. vesiculosus is present only in isolated sheltered areas such as estuaries and coastal lagoons . Consequently, along the broad areas of sympatry of their Northern distribution, the two sister species can hybridize [33, 34, 49], whereas in Southern Europe, their distribution is allopatric and hybridization is highly unlikely due to the limited dispersal capability of their gametes .
A more complex pattern appears on the networks in Northwest Iberia area (Figure 2). The Northwest Iberian close connection to the Northern cluster (North Sea and Channel) indicates a recent secondary expansion either from North Iberia into the Northern region or vice versa, having kept restricted gene flow with the southern clusters.
Another interesting point is the fact that hybrids are mostly localized in the Southern limit of the sympatric zone, where it contacts the allopatric zone ranging from Southern Portugal-Northwestern Africa and Atlantic islands. Hybridization may be favored by unknown factors in this area, such as the rate of hybrid fertility that may change between geographical locations, possibly being higher in Northwest Iberia area as suggest by Billard et al. . Also, reinforcement may be weak due to the recent gene flow from an allopatric zone where it would be lacking, whereas in the middle of the sympatric zone such mechanism would already have developed to maintain species integrity. This is further supported by the peculiar position and clustering of the individuals from the extreme edge of the sympatric zone in Mindelo, where more than half intermediate individuals were detected on the network, and which is the only sympatric population of F. spiralis clustering with the allopatric ones with Structure (k = 3; see Figure A7 in additional file 8). These results are in agreement with a recent multi-gene phylogeny of these species that reveals that in Northwest Iberia F. spiralis from southern origins become extremely introgressed .
The results of the present study suggest three main conclusions: (i) The accuracy of network analysis to unravel hybridization phenomena. The detection of hybrids and introgression is possible and reliable through network analysis, although not systematic as some remain apparently undetected in any analysis attempted, either networks or STRUCTURE. (ii) The putative hybrids detected in the SAD network seem more similar to backcrosses (indifferently with one or the other of the parental species) than to F1 Hybrids. This may be interpreted as a seldom occurrence of first generation hybrids, with a dilution effect on the mosaic genomes produced after backcross events. (iii) The biogeography of both Fucus species addressed through the examination of networks supports the existence of Northern and Southern glacial refuges where both species differentiated in two clusters that came back into contact in Northwest Iberia during post-glacial range changes. Interestingly, this area corresponds to the boundary between the allopatric/sympatric zones, where most putative hybrids were detected, suggesting weaker reinforcement in this area.
Dataset used in this study
A total of 572 individuals of F_spi and F_ves (respectively 329 and 243) were collected throughout the North Hemisphere (see FigureA1 and TableA1 for location distribution in additional files 1 and 9). The analysis covers their entire latitudinal range. A combination of nine microsatellite markers was selected for their amplification consistency and polymorphism in both species. The marker combination includes L20-58-78-38-94 , Fsp 1-2-4  and F90 .
Genetic distances used in this study
Two different genetic distances taking or not phylogenetic information into account were chosen to perform network analysis of our dataset.
First, we used the Goldstein distance for populations  modified for use between individuals, as originally presented by Rozenfeld et al. . As the underlying data are based on multilocus microsatellite genotypes, an individual is characterized by a series of pairs of microsatellite repetitions at k loci with k = 9.
where ai and Ai are the allele length (in number of nucleotides) in both chromosomes at locus i.
which provides a parsimonious (i.e. minimal) representation of the genetic distance, understood as the difference in allele length, between samples A and B.
where the number of shared alleles S is summed over all loci u and nu is the number of loci.
This individual measure can be used to look at population substructure. Bowcock et al.  constructed dendrograms based on this distance calculated from human microsatellite data. Using this technique, a correlation between genetic similarity and geographic location was noted. This distance measure has also proven very successful at placing unknown individuals into the correct subpopulation .
Finally, these distances help to resolve the relationship between individuals at different time-scale. RD, based on loci, helps to resolve individuals' origin at an older time, while SAD, based on shared alleles, helps to understand recent gene flow.
Once we have calculated the matrices of genetic distances between individuals described above containing all the individuals of our dataset, we built networks by considering individuals as nodes and genetic distances between them as links.
as a function of the last distance value removed, thr. N is the total number of nodes not included in the largest cluster and ns is the number of clusters containing s nodes. The resulting curves show a maximum followed by a strong decrease where the percolation threshold is positioned. Once this percolation threshold is identified, we analyzed the network topology and its characteristics at this point.
Estimate the global and local properties of the Network
The different indexes of network theory used to describe and characterize our network are:
The connectivity degree ki of a given node i is the number of other nodes linked to it (i.e., the number of neighbor nodes).
The distribution P(k) gives the proportion of nodes in the network having degree k.
The number of links Ei existing among the neighbours of node i. This quantity takes values between 0 and Ei(max) = ki(ki-1)/2, which is the case of a fully connected neighborhood.
It quantifies how close the node i and its neighbors are to being a clique (complete graph).
The clustering coefficient of the whole network < CC > is defined as the average of all individual clustering coefficients in the system.
Genetic diversity spectrum
The genetic distance between pairs of individuals within all locations was calculated in order to plot the frequency distribution of all pairwise values. This distribution is referred as the genetic diversity spectrum (GDS) as defined by Rozenfeld et al. .
We generated hybrids by random simulation of individuals of which the genotype is half F_ves/half F_spi. To do so, we used a random generator PERL script that selects randomly natural individuals present in locations where hybridization was suspected according to the SD network topology (see results). Then, for each of the nine loci, alleles (one from F_ves and the other from F_spi) were randomly chosen to create hybrids of first generation (F1). Among the natural individuals that were chosen, we excluded the natural putative hybrids, detected at the interface of the clusters of individuals F. spiralis and F. vesiculosus (see Table A2 in additional file 10), in order to avoid backcross hybridization in the same set of simulations. The natural putative hybrids are the individuals F. spiralis and F. vesiculosus connected between them allowing the information flow between the clusters of the species inside the network.
For the backcrosses, we followed the same scheme but we changed one of the natural individual by a hybrid F1. For example, backcrosses F. spiralis are the fusion of natural F. spiralis individuals and the hybrids F1. The detailed methodology used to generate hybrids and insert them into the networks is available in additional file 5.
Admixture proportions of all individuals
In order to test for the ability of the network analysis to detect hybrids, we compared its efficiency with the program STRUCTURE . This method has been previously used in the literature to identify species and hybrids in these species of Fucus, using microsatellite data [29, 34]. STRUCTURE was run on both raw data and on data including synthetic hybrids for two clusters (k = 2) assuming admixture and independent allele frequencies between F_spi and F_ves. Analyses were performed using a burn-in period of 5,000 followed by 10,000 Markov Chain Monte Carlo repetitions. Individuals were considered as intermediate genotypes when they have more than 10% of ancestry coming from one of the 2 species. The confidence interval was computed ('print credible region' parameter) and intermediate genotypes labeled as 'significantly admixed' when the confidence interval was strictly included between 0.1 and 0.9. Results are detailed mentioning n as the number of admixed individuals (admixture proportion included between 0.1 and 0.9) and ns the number of those admixed individuals included in the 90% confidence interval.
The authors thank many international researchers throughout Europe and the USA who collected samples from various regions (namely Arthur Mathieson, Megan Dethier, Karl Gunnarson, Christine Maggs, Christy Henzler, Henning Steen, Carolyn Engel, Claire Daguin, Ana Neto, Candelaria Gil-Rodriguez) and Mirjam van de Vliet for the genotyping work. The authors thank Alejandro Rozenfeld for his methodology of network analysis and the members of the European project EDEN and the networks MARBEF, Marine Genomics Europe and CORONA for helpful discussions. This project was funded by projects of the Portuguese Science Foundation (FCT, and FEDER) and by a FCT (and ESF) fellowship to CP.
- Hewitt GM: Hybrid zones-natural laboratories for evolutionary studies. Trends in Ecology & Evolution. 1988, 3 (7): 158-167.View ArticleGoogle Scholar
- Orr MR, Smith TB: Ecology and speciation. Trends in Ecology & Evolution. 1998, 13 (12): 502-506.View ArticleGoogle Scholar
- Schluter D: Ecology and the origin of species. Trends in Ecology & Evolution. 2001, 16 (7): 372-380.View ArticleGoogle Scholar
- Via S: The ecological genetics of speciation. The American Naturalist. 2002, 159 (Suppl 3): S1-7-S1-7. 10.1086/338368.View ArticlePubMedGoogle Scholar
- Odling-Smee FJ, Laland KN, Feldman MW: Niche construction. 2003, Princeton University PressGoogle Scholar
- Albert V, Jonsson B, Bernatchez L: Natural hybrids in Atlantic eels (Anguilla anguilla, A. rostrata): evidence for successful reproduction and fluctuating abundance in space and time. Molecular ecology. 2006, 15 (7): 1903-1916. 10.1111/j.1365-294X.2006.02917.x.View ArticlePubMedGoogle Scholar
- Barton NH, Hewitt GM: Adaptation, speciation and hybrid zones. Nature. 1989, 341 (6242): 497-503. 10.1038/341497a0.View ArticlePubMedGoogle Scholar
- Jiggins CD, Mallet J: Bimodal hybrid zones and speciation. Trends in Ecology & Evolution. 2000, 15 (6): 250-255.View ArticleGoogle Scholar
- Bierne N, Bonhomme F, David P: Habitat preference and the marine-speciation paradox. Proceedings Biological Sciences/The Royal Society. 2003, 270 (1522): 1399-1406. 10.1098/rspb.2003.2404.View ArticleGoogle Scholar
- Ogden R, Thorpe RS: Molecular evidence for ecological speciation in tropical habitats. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (21): 13612-13615. 10.1073/pnas.212248499.View ArticlePubMedPubMed CentralGoogle Scholar
- Vines TH, Köhler SC, Thiel M, Ghira I, Sands TR, MacCallum CJ, Barton NH, Nürnberger B: The maintenance of reproductive isolation in a mosaic hybrid zone between the fire-bellied toads Bombina bombina and B. variegata. Evolution; International Journal of Organic Evolution. 2003, 57 (8): 1876-1888.View ArticlePubMedGoogle Scholar
- Beerli P, Felsenstein J: Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. Genetics. 1999, 152 (2): 763-773.PubMedPubMed CentralGoogle Scholar
- Chikhi L, Bruford MW, Beaumont MA: Estimation of admixture proportions: a likelihood-based approach using Markov chain Monte Carlo. Genetics. 2001, 158 (3): 1347-1362.PubMedPubMed CentralGoogle Scholar
- Choisy M, Franck P, Cornuet JM: Estimating admixture proportions with microsatellites: comparison of methods based on simulated data. Molecular Ecology. 2004, 13 (4): 955-968. 10.1111/j.1365-294X.2004.02107.x.View ArticlePubMedGoogle Scholar
- Mousset S, Derome N, Veuille M: A test of neutrality and constant population size based on the mismatch distribution. Mol Biol Evol. 2004, 21 (4): 724-731. 10.1093/molbev/msh066.View ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMedPubMed CentralGoogle Scholar
- Dyer RJ, Nason JD: Population Graphs: the graph theoretic shape of genetic structure. Mol Ecol. 2004, 13 (7): 1713-1727. 10.1111/j.1365-294X.2004.02177.x.View ArticlePubMedGoogle Scholar
- Fortuna MA, Albaladejo RG, Fernández L, Aparicio A, Bascompte J: Networks of spatial genetic variation across species. Proceedings of the National Academy of Sciences. 2009, 106 (45): 19044-19049. 10.1073/pnas.0907704106.View ArticleGoogle Scholar
- Rozenfeld AF, Arnaud-Haond S, Hernandez-Garcia E, Eguiluz VM, Matias MA, Serrao E, Duarte CM: Spectrum of genetic diversity and networks of clonal organisms. Journal of the Royal Society, Interface/the Royal Society. 2007, 4 (17): 1093-1102. 10.1098/rsif.2007.0230.View ArticleGoogle Scholar
- Rozenfeld AF, Arnaud-Haond S, Hernandez-Garcia E, Eguiluz VM, Serrao EA, Duarte CM: Network analysis identifies weak and strong links in a metapopulation system. Proc Natl Acad Sci USA. 2008, 105 (48): 18824-18829. 10.1073/pnas.0805571105.View ArticlePubMedPubMed CentralGoogle Scholar
- Griffiths RC, Tavare S: Simulating probability distributions in the coalescent. Theoretical Population Biology. 1994, 46 (2): 131-159. 10.1006/tpbi.1994.1023.View ArticleGoogle Scholar
- Wilson GA, Rannala B: Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003, 163 (3): 1177-1191.PubMedPubMed CentralGoogle Scholar
- Anderson EC, Thompson EA: A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002, 160 (3): 1217-1229.PubMedPubMed CentralGoogle Scholar
- Amaral LA, Scala A, Barthelemy M, Stanley HE: Classes of small-world networks. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (21): 11149-11152. 10.1073/pnas.200327197.View ArticlePubMedPubMed CentralGoogle Scholar
- Strogatz SH: Exploring complex networks. Nature. 2001, 410 (6825): 268-276. 10.1038/35065725.View ArticlePubMedGoogle Scholar
- Proulx SR, Promislow DEL, Phillips PC: Network thinking in ecology and evolution. Trends Ecol Evol. 2005, 20 (6): 345-353. 10.1016/j.tree.2005.04.004.View ArticlePubMedGoogle Scholar
- Barabasi AL, Albert R: Emergence of scaling in random networks. Science (New York, NY. 1999, 286 (5439): 509-512. 10.1126/science.286.5439.509.View ArticleGoogle Scholar
- Arnaud-Haond S, Migliaccio M, Diaz-Almela E, Teixeira S, van de Vliet MS, Alberto F, Procaccini G, Duarte CM, Serrao EA: Vicariance patterns in the Mediterranean Sea: east-west cleavage and low dispersal in the endemic seagrass Posidonia oceanica. J Biogeogr. 2007, 34 (6): 963-976. 10.1111/j.1365-2699.2006.01671.x.View ArticleGoogle Scholar
- Billard E, Serrao E, Pearson G, Destombe C, Valero M: Fucus vesiculosus and spiralis species complex: a nested model of local adaptation at the shore level. Marine Ecology-Progress Series. 2010, 405: 163-174. 10.3354/meps08517.View ArticleGoogle Scholar
- Zardi GI, Nicastro KR, Canovas F, Costa JF, Serrão EA, Pearson GA: Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone. PLoS ONE.
- Coyer JA, Hoarau G, Pearson GA, Serrao EA, Stam WT, Olsen JL: Convergent adaptation to a marginal habitat by homoploid hybrids and polyploid ecads in the seaweed genus Fucus. Biology Letters. 2006, 2 (3): 405-408. 10.1098/rsbl.2006.0489.View ArticlePubMedPubMed CentralGoogle Scholar
- Coyer JA, Hoarau G, Costa JF, Hogerdijk B, Serrão EA, Billard E, Valero M, Pearson GA, Olsen JL: Evolution and diversification within the intertidal brown macroalgae Fucus spiralis/F. vesiculosusspecies complex in the North Atlantic. Molecular Phylogenetics and Evolution.
- Wallace AL, Klein AS, Mathieson AC: Determining the affinities of salt marsh fucoids using microsatellite markers: Evidence of hybridization and introgression between two species of Fucus (Phaeophyta) in a Maine estuary. Journal of Phycology. 2004, 40 (6): 1013-1027. 10.1111/j.1529-8817.2004.04085.x.View ArticleGoogle Scholar
- Engel CR, Daguin C, Serrao EA: Genetic entities and mating system in hermaphroditic Fucus spiralis and its close dioecious relative F. vesiculosus (Fucaceae, Phaeophyceae). Molecular ecology. 2005, 14 (7): 2033-2046. 10.1111/j.1365-294X.2005.02558.x.View ArticlePubMedGoogle Scholar
- Billard E, Daguin C, Pearson G, Serrao E, Engel C, Valero M: Genetic isolation between three closely related taxa: Fucus vesiculosus, F. spiralis, and F. ceranoides (Phaophyceae). Journal of Phycology. 2005, 41 (4): 900-905. 10.1111/j.0022-3646.2005.04221.x.View ArticleGoogle Scholar
- Wallace AL, Klein AS, Mathieson AC: A reply to Engel et al. Molecular ecology. 2006, 15 (4): 1185-1187. 10.1111/j.1365-294X.2005.02759.x.View ArticlePubMedGoogle Scholar
- Engel CR, Daguin C, Serrao EA: When is a hybrid a hybrid? A counter-reply to Wallace et al. Molecular ecology. 2006, 15 (11): 3481-3482. 10.1111/j.1365-294X.2006.02979.x.View ArticlePubMedGoogle Scholar
- Chakraborty R, Jin L: Determination of relatedness between individuals using DNA fingerprinting. Human Biology; an International Record of Research. 1993, 65 (6): 875-895.Google Scholar
- Goldstein DB, Ruiz Linares A, Cavalli-Sforza LL, Feldman MW: An evaluation of genetic distances for use with microsatellite loci. Genetics. 1995, 139 (1): 463-471.PubMedPubMed CentralGoogle Scholar
- Hewitt GM: Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences. 2004, 359 (1442): 183-195. 10.1098/rstb.2003.1388.View ArticleGoogle Scholar
- Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT, Kelly J, Olsen J, Perez KE, Stam W, Vainola R, et al: Evaluating signatures of Glacial Refugia for North Atlantic benthic marine taxa. Ecology. 2008, 89 (11): S108-S122. 10.1890/08-0257.1.View ArticlePubMedGoogle Scholar
- Gabrielsen TM, Brochmann C, Rueness J: The Baltic Sea as a model system for studying postglacial colonization and ecological differentiation, exemplified by the red alga Ceramium tenuicorne. Molecular ecology. 2002, 11 (10): 2083-2095. 10.1046/j.1365-294X.2002.01601.x.View ArticlePubMedGoogle Scholar
- Coyer JA, Peters AF, Stam WT, Olsen JL: Post-ice age recolonization and differentiation of Fucus serratus L. (Phaeophyceae; Fucaceae) populations in Northern Europe. Molecular ecology. 2003, 12 (7): 1817-1829. 10.1046/j.1365-294X.2003.01850.x.View ArticlePubMedGoogle Scholar
- Wares JP, Gaines SD, Cunningham CW: A comparative study of asymmetric migration events across a marine biogeographic boundary. Evolution. 2001, 55 (2): 295-306.View ArticlePubMedGoogle Scholar
- Hoarau G, Coyer JA, Veldsink JH, Stam WT, Olsen JL: Glacial refugia and recolonization pathways in the brown seaweed Fucus serratus. Molecular ecology. 2007, 16 (17): 3606-3616. 10.1111/j.1365-294X.2007.03408.x.View ArticlePubMedGoogle Scholar
- Neiva J, Pearson GA, Valero M, Serrao EA: Surfing the wave on a borrowed board: range expansion and spread of introgressed organellar genomes in the seaweed Fucus ceranoides L. Molecular ecology. 2010, 19 (21): 4812-4822. 10.1111/j.1365-294X.2010.04853.x.View ArticlePubMedGoogle Scholar
- Muhlin JF, Brawley SH: Recent versus relic: discerning the genetic signature of Fucus vesiculosus (Heterokontophyta; Phaeophyceae) in the Northwestern Atlantic. Journal of Phycology. 2009, 45 (4): 828-837. 10.1111/j.1529-8817.2009.00715.x.View ArticlePubMedGoogle Scholar
- Ladah L, Bermudez R, Pearson G, Serrao E: Fertilization success and recruitment of dioecious and hermaphroditic fucoid seaweeds with contrasting distributions near their southern limit. Marine Ecology-Progress Series. 2003, 262: 173-183. 10.3354/meps262173.View ArticleGoogle Scholar
- Billard E, Serrao EA, Pearson GA, Engel CR, Destombe C, Valero M: Analysis of sexual phenotype and prezygotic fertility in natural populations of Fucus spiralis, F. vesiculosus (Fucaceae, Phaeophyceae) and their putative hybrids. European Journal of Phycology. 2005, 40 (4): 397-407. 10.1080/09670260500334354.View ArticleGoogle Scholar
- Pearson GA, Serrao EA: Revisiting synchronous gamete release by fucoid algae in the intertidal zone: fertilization success and beyond?. Integrative and Comparative Biology. 2006, 46 (5): 587-597. 10.1093/icb/icl030.View ArticlePubMedGoogle Scholar
- Perrin C, Daguin C, Van De Vliet M, Engel CR, Pearson GA, Serrao EA: Implications of mating system for genetic diversity of sister algal species: Fucus spiralis and Fucus vesiculosus (Heterokontophyta, Phaeophyceae). European Journal of Phycology. 2007, 42 (3): 219-230. 10.1080/09670260701336554.View ArticleGoogle Scholar
- Goldstein DB, Ruiz Linares A, Cavalli-Sforza LL, Feldman MW: Genetic absolute dating based on microsatellites and the origin of modern humans. Proceedings of the National Academy of Sciences of the United States of America. 1995, 92 (15): 6723-6727. 10.1073/pnas.92.15.6723.View ArticlePubMedPubMed CentralGoogle Scholar
- Bowcock AM, Ruiz-Linares A, Tomfohrde J, Minch E, Kidd JR, Cavalli-Sforza LL: High resolution of human evolutionary trees with polymorphic microsatellites. Nature. 1994, 368 (6470): 455-457. 10.1038/368455a0.View ArticlePubMedGoogle Scholar
- Estoup A, Garnery L, Solignac M, Cornuet JM: Microsatellite variation in honey bee (Apis mellifera L.) populations: hierarchical genetic structure and test of the infinite allele and stepwise mutation models. Genetics. 1995, 140 (2): 679-695.PubMedPubMed CentralGoogle Scholar
- Stauffer D, Aharony A: Introduction to percolation theory. 1994, Taylor & Francis, LondonGoogle Scholar
- Freeman LC: Set of measures of centrality based on betweenness. Sociometry Sociometry. 1977, 40 (1): 35-41.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.