Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians

BACKGROUND: Identifying species of organisms by short sequences of DNA has been in the center of ongoing discussions under the terms DNA barcoding or DNA taxonomy. A C-terminal fragment of the mitochondrial gene for cytochrome oxidase subunit I (COI) has been proposed as universal marker for this purpose among animals. RESULTS: Herein we present experimental evidence that the mitochondrial 16S rRNA gene fulfills the requirements for a universal DNA barcoding marker in amphibians. In terms of universality of priming sites and identification of major vertebrate clades the studied 16S fragment is superior to COI. Amplification success was 100% for 16S in a subset of fresh and well-preserved samples of Madagascan frogs, while various combination of COI primers had lower success rates.COI priming sites showed high variability among amphibians both at the level of groups and closely related species, whereas 16S priming sites were highly conserved among vertebrates. Interspecific pairwise 16S divergences in a test group of Madagascan frogs were at a level suitable for assignment of larval stages to species (1-17%), with low degrees of pairwise haplotype divergence within populations (0-1%). CONCLUSION: We strongly advocate the use of 16S rRNA as standard DNA barcoding marker for vertebrates to complement COI, especially if samples a priori could belong to various phylogenetically distant taxa and false negatives would constitute a major problem.     

Phylogeny, species limits, and biogeography of the Brazilian lizards of the genus Eurolophosaurus (Squamata: Tropiduridae) as inferred from mitochondrial DNA sequences

Phylogenetic relationships and divergence times for 10 populations of the three recognized "species" of Brazilian lizards of genus Eurolophosaurus were estimated from 1229bp of cyt b, COI, 12S, and 16S rRNA mitochondrial gene segments. Eurolophosaurus is monophyletic and the basal split within the genus separates E. divaricatus from a clade comprising E. amathites and E. nanuzae. Three populations of E. divaricatus, which occurs along the western bank of Rio S?o Francisco, were consistently grouped together. On the east bank of the river, E. amathites and E. nanuzae from state of Bahia were recovered as the sister group of E. nanuzae populations from state of Minas Gerais. The paraphyly of E. nanuzae and the high divergence levels among populations of E. divaricatus strongly suggest that species limits in Eurolophosaurus should be revised. Even considering an extreme evolutionary rate of 2.8% sequence divergence per million years for the four gene segments analyzed together, E. divaricatus would have separated from the two other species by at least 5.5my ago, and E. amathites from E. nanuzae populations from Bahia and Minas Gerais, respectively, by 1.5 and 3.5my. The paleolacustrine hypothesis and changes in the course of the river potentially explain faunal divergence in the area, but divergences are much older than previously admitted.     

Lack of evidence for reproductive isolation among ecologically specialised lycaenid butterflies

Abstract 1. The evolution of reproductive isolation between recently diverged or incipient species is a critical component of speciation and a major focus of speciation models. In phytophagous insects, host plant fidelity (the habit of mating and ovipositing on a single host species) can contribute to assortative mating and reproductive isolation between populations adapting to alternative hosts. The potential role of host plant fidelity in the evolution of reproductive isolation was examined in a pair of North American blue butterfly species, Lycaeides idas and L. melissa .
2. These species are morphologically distinct and populations of each species utilise different host plants; however they share 410 bp haplotypes of the mitochondrial cytochrome oxidase subunit I (COI) gene, indicating recent divergence.
3. Some populations using native hosts exhibited strong fidelity for their natal host plant over the hosts used by nearby populations. Because these butterflies mate on or near the host plant, the development of strong host fidelity may create reproductive isolation among populations on different hosts and restrict gene flow.
4. Tests of population differentiation using allozyme allele frequency data did not provide convincing evidence of restricted gene flow among populations. Based on morphological differences, observed ecological specialisation, and the sharing of genetic markers, these butterflies appear to be undergoing adaptive radiation driven at least partially by host shifts. Neutral genetic markers may fail to detect the effects of very recent host shifts in these populations due to gene flow and/or the recency of divergence and shared ancestral polymorphism.     

Phylogenetic Relationships of 12 Penaeoidea Shrimp Species Deduced from Mitochondrial DNA Sequences

DNA sequences of an 847 bp fragment of mitochondrial cytochrome oxidase subunit I (COI) gene and a 514 bp fragment of 16s rRNA gene were determined to examine the phylogenetic relationships of 12 Penaeoidea shrimp species (Penaeus chinensis, Penaeus japonicus, Penaeus penicillatus, Penaeus vannamei, Penaeus canaliculatus, Trachypenaeus curvirostris, Metapenaeus affinis, Metapenaeus ensis, Metapenaeopsis barbata, Parapenaeus fissuroides, Parapenaeopsis hardiwickii, Solenocera crassicomis). Both fragments of the swimming crab Portunus trituberculaus chosen as the outgroup were also sequenced. Intraspecific sequence divergence of 0.24-1.2% in the COI gene was found in 5 species, while no intraspecific variation was observed in the 16s rRNA gene. Three phylogenetic trees based on the 1361 bp combined sequences of COI and 16s rRNA were concordant in indicating the following suggestions: (1) phylogenetic relationship of the 11 Penaeidae species based on our result support the opinion of Burkenroad (Burkenroad, M.D. (1983). Crustacean Issues 3:279-290) on the basis of morphological features; (2) it seems more reasonable to class Solenocera crassicorni in the family Penaeidae; (3) the fragment of the COI gene chosen here appears to be a good marker for speciation studies and population analysis in Crustaceans, while the 16s rRNA gene fragment here seems suitable for examining phylogenetic relationships at the species or genus levels in Crustaceans. Our time estimates suggest that Penaeus and Metapenaeus might have separated about 6.38 x 10(6)-7.98 x 10(6) years BP in the post-Miocene, and the species separation within Metapenaeus and Penaeus might occur 0.08 x 10(6)-0.4 x 10(6) years BP in the late Pleistocene.     

Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea,inferred from mitochondria DNA markers

The Asian green mussel Perna viridis is ecologically and economically important in the coastal regions of China. In order to characterize the genetic diversity and population connectivity of P. viridis in South China Sea, a 664 bp region of mitochondrial COI gene and a 293 bp region of 16S rRNA gene were sequenced and analyzed for 78 and 92 individuals from four populations in South China Sea, respectively. A total of 15 haplotypes were defined by 14 variable nucleotide sites in COI gene, and 7 haplotypes by 6 variable nucleotide sites in 16S rRNA gene. High haplotype diversity and low nucleotide diversity were observed in COI gene, while moderate haplotype diversity and low nucleotide diversity were observed in 16S rRNA gene. Pairwise F_ST values of COI gene were all negative and those of 16S rRNA gene ranged from −0.01409 to 0.10289. The results showed that no significant genetic divergence (or shallow genetic structure) and high levels of population connectivity among the four populations of P. viridis in South China Sea.     

Tests of pleistocene speciation in montane grasshoppers (genus Melanoplus) from the sky islands of western North America

There has a been a resurgence of debate on whether the Pleistocene glaciations inhibited speciation. This study tests a model of Pleistocene speciation, estimating the phylogenetic relationships and divergence times of 10 species of montane grasshoppers, genus Melanoplus, using 1300 bp of the mitochondrial gene cytochrome oxidase I (COI). Based on average pairwise distances (corrected for multiple substitutions using Kimura's two-parameter model), all species appear to have originated within the Pleistocene. Sequence divergences between species are less than 4%, corresponding to divergence times less than 1.7 million years ago. Branching patterns among the species suggest that speciation was associated with more than one glacial-interglacial cycle. A likelihood-ratio test rejected a model of simultaneous species origins, the predicted branching pattern if species arose from the fragmentation of a widespread ancestor. These grasshoppers live in an area that was previously glaciated and, as inhabitants of the northern Rocky Mountain sky islands, underwent latitudinal and probably altitudinal shifts in distribution in response to climatic fluctuations. Given the repeated distributional shifts and range overlap of the taxa, there most likely has been ample opportunity for population mixing. However, despite periodic glacial cycles, with more than 10 major glaciations over the past million years and climatic fluctuations over as short a time scale as 10(3) to 10(4) years, the dynamic history of the Pleistocene did not preclude speciation. Although relationships among some taxa remain unresolved, these grasshopper species, even with their recent origins, exhibit genetic coherence and monophyletic or paraphyletic gene trees. The frequency of glacial cycles suggests that the speciation process must have been extremely rapid. These species of grasshoppers are morphologically very similar, differing primarily in the shape of the male genitalia. These characters are posited to be under sexual selection, may play an important role in reproductive isolation, and are known to diverge rapidly. This suggests the rapidity of evolution of reproductive isolation may determine whether species divergences occurred during the Pleistocene glaciations.     

Weak premating isolation between Clitarchus stick insect species despite divergent male and female genital morphology

Documenting natural hybrid systems builds our understanding of mate choice, reproductive isolation and speciation. The stick insect species Clitarchus hookeri and C. tepaki differ in their genital morphology and hybridize along a narrow peninsula in northern New Zealand. We utilize three lines of evidence to understand the role of premating isolation and species boundaries: (a) genetic differentiation using microsatellites and mitochondrial DNA; (b) variation in 3D surface topology of male claspers and 2D morphometrics of female opercular organs; and (c) behavioural reproductive isolation among parental and hybrid populations through mating crosses. The genetic data show introgression between the parental species and formation of a genetically variable hybrid swarm. Similarly, the male and female morphometric data show genital divergence between the parental species as well as increased variation within the hybrid populations. This genital divergence has not resulted in reproductive isolation between species, instead weak perimating isolation has enabled the formation of a hybrid swarm. Behavioural analysis demonstrates that the entire mating process influences the degree of reproductive isolation between species undergoing secondary contact. Mechanical isolation may appear strong, whereas perimating isolation is weak.     

COI is better than 16S rRNA for DNA barcoding Asiatic salamanders (Amphibia: Caudata: Hynobiidae)

The 5' region of the mitochondrial DNA (mtDNA) gene cytochrome c oxidase I (COI) is the standard marker for DNA barcoding. However, because COI tends to be highly variable in amphibians, sequencing is often challenging. Consequently, another mtDNA gene, 16S rRNA gene, is often advocated for amphibian barcoding. Herein, we directly compare the usefulness of COI and 16S in discriminating species of hynobiid salamanders using 130 individuals. Species identification and classification of these animals, which are endemic to Asia, are often based on morphology only. Analysis of Kimura 2-parameter genetic distances (K2P) documents the mean intraspecific variation for COI and 16S rRNA genes to be 1.4% and 0.3%, respectively. Whereas COI can always identify species, sometimes 16S cannot. Intra- and interspecific genetic divergences occasionally overlap in both markers, thus reducing the value of a barcoding gap to identify genera. Regardless, COI is the better DNA barcoding marker for hynobiids. In addition to the comparison of two potential markers, high levels of intraspecific divergence in COI (>5%) suggest that both Onychodactylus fischeri and Salamandrella keyserlingii might be composites of cryptic species.     

Extensive loss of Wnt genes in Tardigrada

Evolution of reproductive isolation is an important process, generating biodiversity and driving speciation. To better understand this process, it is necessary to investigate factors underlying reproductive isolation through various approaches but also in various taxa. Previous studies, mainly focusing on diploid animals, supported the prevalent view that reproductive barriers evolve gradually as a by-product of genetic changes accumulated by natural selection by showing a positive relationship between the degree of reproductive isolation and genetic distance. Haplodiploid animals are expected to generate additional insight into speciation, but few studies investigated the prevalent view in haplodiploid animals. In this study, we investigate whether the relationship also holds in a haplodiploid spider mite, Amphitetranychus viennensis (Zacher). We sampled seven populations of the mite in the Palaearctic region, measured their genetic distance (mtDNA) and carried out cross experiments with all combinations. We analyzed how lack of fertilization rate (as measure of prezygotic isolation) as well as hybrid inviability and hybrid sterility (as measures of postzygotic isolation) varies with genetic distance. We found that the degree of reproductive isolation varies among cross combinations, and that all three measures of reproductive isolation have a positive relationship with genetic distance. Based on the mtDNA marker, lack of fertilization rate, hybrid female inviability and hybrid female sterility were estimated to be nearly complete (99.0–99.9% barrier) at genetic distances of 0.475–0.657, 0.150–0.209 and 0.145–0.210, respectively. Besides, we found asymmetries in reproductive isolation. The prevalent view on the evolution of reproductive barriers is supported in the haplodiploid spider mite we studied here. According to the estimated minimum genetic distance for total reproductive isolation in parent population crosses in this study and previous work, a genetic distance of 0.15–0.21 in mtDNA (COI) appears required for speciation in spider mites. Variations and asymmetries in the degree of reproductive isolation highlight the importance of reinforcement of prezygotic reproductive isolation through incompatibility and the importance of cytonuclear interactions for reproductive isolation in haplodiploid spider mites.     

Variation in reproductive isolation across a species range

Reproductive isolation is often variable within species, a phenomenon that while largely ignored by speciation studies, can be leveraged to gain insight into the potential mechanisms driving the evolution of genetic incompatibilities. We used experimental greenhouse crosses to characterize patterns of reproductive isolation among three divergent genetic lineages of Campanulastrum americanum that occur in close geographic proximity in the Appalachian Mountains. Substantial, asymmetrical reproductive isolation for survival due to cytonuclear incompatibility was found among the lineages (up to 94% reduction). Moderate reductions in pollen viability, as well as cytoplasmic male sterility, were also found between some Mountain populations. We then compared these results to previously established patterns of reproductive isolation between these Mountain lineages and a fourth, widespread Western lineage to fully characterize reproductive isolation across the complete geographic and genetic range of C. americanum. Reproductive isolation for survival and pollen viability was consistent across studies, indicating the evolution of the underlying genetic incompatibilities is primarily determined by intrinsic factors. In contrast, reproductive isolation for germination was only found when crossing Mountain populations with the Western lineage, suggesting the underlying genetic incompatibility is likely influenced by environmental or demographic differences between the two lineages. Cytoplasmic male sterility was also limited in occurrence, being restricted to a handful of Mountain populations in a narrow geographic range. These findings illustrate the complexity of speciation by demonstrating multiple, independent genetic incompatibilities that lead to a mosaic of genetic divergence and reproductive isolation across a species range.     

Molecular phylogeny, taxonomy and evolution of the land snail genus Iberus (Pulmonata: Helicidae)

Partial DNA sequences of two mitochondrial genes [cytochrome oxidase subunit I (COI) and 16S rRNA] from 59 specimens of Iberus were used to test the validity of the described morphospecies of this genus, and examine genetic divergences within and between main phylogenetic groups. Both gene fragments showed phylogenetic concordance. The COI gene was found to be faster evolving than the 16S gene and was fully protein-coding with no insertions or deletions. 16S rRNA was more informative than COI for resolving basal nodes. Both individual and combined analyses of the two gene fragments revealed five main phylogroups. These five groups are genetically unique lineages that are allopatrically distributed and considered to have full species status. Further subdivisions were also considered. Shell morphology was suitable for delimiting species boundaries, but several incongruences between morphology and mtDNA phylogeny were observed. These incongruences were considered consequence of hybridization between Iberus cobosi and Iberus marmoratus , and the result of shell shape polymorphism in Iberus rositai . According to spatial patterns of sequence divergence, life habits and shell morphology may be concluded that the keeled-flat shelled snails independently originated several times within Iberus and they could represent cases of similar shell adaptation to a karstic arid environment.     

Biogeography in a continental island: population structure of the relict endemic centipede Craterostigmus tasmanianus (Chilopoda, Craterostigmomorpha) in Tasmania using 16S rRNA and COI

We used 16S ribosomal RNA (rRNA) and cytochrome c oxidase subunit I (COI) sequence data to investigate the population structure in the centipede Craterostigmus tasmanianus Pocock, 1902 (Chilopoda: Craterostigmomorpha: Craterostigmidae) and to look for possible barriers to gene flow on the island of Tasmania, where C. tasmanianus is a widespread endemic. We first confirmed a molecular diagnostic character in 28S rRNA separating Tasmanian Craterostigmus from its sister species Craterostigmus crabilli (Edgecombe and Giribet 2008) in New Zealand and found no shared polymorphism in this marker for the 2 species. In Tasmania, analysis of molecular variance analysis showed little variation at the 16S rRNA and COI loci within populations (6% and 13%, respectively), but substantial variation (56% and 48%, respectively) among populations divided geographically into groups. We found no clear evidence of isolation by distance using a Mantel test. Bayesian clustering and gene network analysis both group the C. tasmanianus populations in patterns which are broadly concordant with previously known biogeographical divisions within Tasmania, but we did not find that genetic distance varied in a simple way across cluster boundaries. The coarse-scale geographical sampling on which this study was based should be followed in the future by sampling at a finer spatial scale and to investigate genetic structure within clusters and across cluster boundaries.     

A molecular phylogenetic approach to the phylum Cycliophora provides further evidence for cryptic speciation in Symbion americanus

Cryptic speciation is a phenomenon that has been reported in a wide range of marine invertebrates, including the American cycliophoran Symbion americanus . Although cryptic speciation is often proposed based on phylogeographical analyses, here we demonstrate the utility of multilocus molecular phylogenetic approaches in identifying cryptic lineages within the phylum Cycliophora. Thirteen individuals from three putative cryptic lineages of S. americanus , as well as two individuals each from the European species S. pandora and a new species of Symbion living on Homarus gammarus , were sequenced for up to 4.8 kb of genomic DNA over four loci (COI, 16S rRNA, 18S rRNA and 28S rRNA). Phylogenetic analyses of individual loci as well as combined data clearly support a division between American and European cycliophorans. Moreover, the American cycliophorans consistently form three well-supported clades, which is congruent with the presence of three putative reproductively isolated lineages in a species complex. Further studies are necessary to more precisely describe the evolution of reproductive isolation within S. americanus .     

DNA barcoding of Scandinavian birds reveals divergent lineages in trans-Atlantic species

Birds are a taxonomically well-described group of animals, yet DNA barcoding, i.e., the molecular characterization of species using a standardized genetic marker, has revealed unexpected patterns of genetic divergences among North American birds. We performed a comprehensive COI (cytochrome c oxidase subunit I) barcode survey of 296 species of Scandinavian birds, and compared genetic divergences among 78 trans-Atlantic species whose breeding ranges include both Scandinavia and North America. Ninety-four percent of the Scandinavian species showed unique barcode clusters; the remaining 6% had overlapping barcodes with one or more congeneric species, which may reflect incomplete lineage sorting or a single gene pool. Four species showed large intra-specific divergences within Scandinavia, despite no apparent morphological differentiation or indications of reproductive isolation. These cases may reflect admixture of previously isolated lineages, and may thus warrant more comprehensive phylogeographic analyses. Nineteen (24%) of 78 trans-Atlantic species exhibited divergent genetic clusters which correspond with regional subspecies. Three of these trans-Atlantic divergences were paraphyletic. Our study demonstrates the effectiveness of COI barcodes for identifying Scandinavian birds and highlights taxa for taxonomic review. The standardized DNA barcoding approach amplified the power of our regional studies by enabling independently obtained datasets to be merged with the established avian barcode library.     

Lack of taxonomic differentiation in an apparently widespread freshwater isopod morphotype (Phreatoicidea: Mesamphisopidae: Mesamphisopus) from South Africa

The unambiguous identification of phreatoicidean isopods occurring in the mountainous southwestern region of South Africa is problematic, as the most recent key is based on morphological characters showing continuous variation among two species: Mesamphisopus abbreviatus and M. depressus. This study uses variation at 12 allozyme loci, phylogenetic analyses of 600 bp of a COI (cytochrome c oxidase subunit I) mtDNA fragment and morphometric comparisons to determine whether 15 populations are conspecific, and, if not, to elucidate their evolutionary relationships. Molecular evidence suggested that the most easterly population, collected from the Tsitsikamma Forest, was representative of a yet undescribed species. Patterns of differentiation and evolutionary relationships among the remaining populations were unrelated to geographic proximity or drainage system. Patterns of isolation by distance were also absent. An apparent disparity among the extent of genetic differentiation was also revealed by the two molecular marker sets. Mitochondrial sequence divergences among individuals were comparable to currently recognized intraspecific divergences. Surprisingly, nuclear markers revealed more extensive differentiation, more characteristic of interspecific divergences. This disparity and the mosaic pattern of differentiation may be driven by stochastic population crashes and genetic bottlenecks (caused by seasonal habitat fluctuations), coupled with genetic drift.     

Does reproductive isolation evolve faster in larger populations via sexually antagonistic coevolution?

Sexual conflict over reproductive investment can lead to sexually antagonistic coevolution and reproductive isolation. It has been suggested that, unlike most models of allopatric speciation, the evolution of reproductive isolation through sexually antagonistic coevolution will occur faster in large populations as these harbour greater levels of standing genetic variation, receive larger numbers of mutations and experience more intense sexual selection. We tested this in bruchid beetle populations (Callosobruchus maculatus) by manipulating population size and standing genetic variability in replicated lines derived from founders that had been released from sexual conflict for 90 generations. We found that after 19 generations of reintroduced sexual conflict, none of our treatments had evolved significant overall reproductive isolation among replicate lines. However, as predicted, measures of reproductive isolation tended to be greater among larger populations. We discuss our methodology, arguing that reproductive isolation is best examined by performing a matrix of allopatric and sympatric crosses whereas measurement of divergence requires crosses with a tester line.     

Long- term effects of previous experience determine nutrient discrimination abilities in birds

Background

The evolution of reproductive traits, such as hybrid incompatibility (postzygotic isolation) and species recognition (prezygotic isolation), have shown their key role in speciation. Theoretical modeling has recently predicted that close linkage between genes controlling pre- and postzygotic reproductive isolation could accelerate the conditions for speciation. Postzygotic isolation could develop during the sympatric speciation process contributing to the divergence of populations. Using hybrid fitness as a measure of postzygotic reproductive isolation, we empirically studied population divergence in perch (Perca fluviatilis L.) from two genetically divergent populations within a lake.

Results

During spawning time of perch we artificially created parental offspring and F₁ hybrids of the two populations and studied fertilization rate and hatching success under laboratory conditions. The combined fitness measure (product of fertilization rate and hatching success) of F₁ hybrids was significantly reduced compared to offspring from within population crosses.

Conclusion

Our results suggest intrinsic genetic incompatibility between the two populations and indicate that population divergence between two populations of perch inhabiting the same lake may indeed be promoted by postzygotic isolation.     

Polymorphism of mtDNA regions in brown hare (Lepus europaeus) populations from Vojvodina (Serbia and Montenegro)

It is well known that there is heterogeneity of mitochondrial DNA (mtDNA) within and among natural populations of same species. The polymorphism level of particular regions of mtDNA gives valuable results in detection of population genetic structure. The aim of this paper was to detect polymorphism of three mtDNA regions: cytochrome oxidase I (COI), Control region, and 12S/16S rRNA, by the mtDNA RFLP-PCR method, in three Lepus europaeus populations from Vojvodina province (Serbia and Montenegro). Polymorphism was detected within the two regions, COI and Control region, while 12S/16S rRNA region was monomorphic in all 77 individuals. Eight haplotypes were detected in the brown hare population in Vojvodina, and three were unique for the Srem brown hare population.     

Gametic isolation in guppies (Poecilia reticulata)

Post-mating reproductive isolating mechanisms may be among the earliest reproductive barriers to emerge among incipient species. Trinidadian guppy, Poecilia reticulata, populations in the Caroni and Oropouche drainages in Northern Trinidad exhibit marked genetic divergence and provide an ideal system in which to search for these barriers. We inseminated virgin females with equal amounts of sperm from two males, a 'native' male from the female's own population and a 'foreign' male from the other drainage. Artificial insemination ensured that mating order and mate choice did not affect the outcome. Paternities were assigned to the resulting broods using microsatellite markers. As predicted, sperm from native males had precedence over foreign sperm. Moreover, this effect was symmetrical for both drainages. In contrast, we detected no native sperm precedence in controls, in which females received sperm from the same and another population within the same drainage. Our results show that gametic isolation can arise between geographically proximate, though genetically divergent, populations of a single species and highlight the potential role of this process in speciation.     

Fine-scale genetic structure, estuarine colonization and incipient speciation in the marine silverside fish Odontesthes argentinensis

The identification of incipient ecological species represents an opportunity to investigate current evolutionary process where adaptive divergence and reproductive isolation are associated. In this study we analysed the genetic structure of marine and estuarine populations of the silverside fish Odontesthes argentinensis using nine microsatellite loci and 396 bp of the mitochondrial DNA (mtDNA) control region. Our main objective was to investigate the relationship among estuarine colonization, divergent selection and speciation in silversides. Significant genetic structure was detected among all marine and estuarine populations. Despite the low phylogeographic structure in mtDNA haplotypes, there was clear signal of local radiations of haplotypes in more ancient populations. Divergence among marine populations was interpreted as a combined result of homing behaviour, isolation by distance and drift. On the other hand, ecological shifts due to the colonization of estuarine habitats seem to have promoted rapid adaptive divergence and reproductive isolation in estuarine populations, which were considered as incipient ecological species. This conclusion is supported by the existence of a set of environmental factors required for successful reproduction of estuarine ecotypes. The pattern of genetic structure indicates that phenotypic and reproductive divergence evolved in the face of potential gene flow between populations. We suggest that the 'divergence-with-gene-flow' model of speciation may account for the diversification of estuarine populations. The approach used can potentially identify 'incipient estuarine species', being relevant to the investigation of the evolutionary relationships of silversides in several coastal regions of the world.