Identifying Fishes through DNA Barcodes and Microarrays

Background

International fish trade reached an import value of 62.8 billion Euro in 2006, of which 44.6% are covered by the European Union. Species identification is a key problem throughout the life cycle of fishes: from eggs and larvae to adults in fisheries research and control, as well as processed fish products in consumer protection.

Methodology/Principal Findings

This study aims to evaluate the applicability of the three mitochondrial genes 16S rRNA (16S), cytochrome b (cyt b), and cytochrome oxidase subunit I (COI) for the identification of 50 European marine fish species by combining techniques of “DNA barcoding” and microarrays. In a DNA barcoding approach, neighbour Joining (NJ) phylogenetic trees of 369 16S, 212 cyt b, and 447 COI sequences indicated that cyt b and COI are suitable for unambiguous identification, whereas 16S failed to discriminate closely related flatfish and gurnard species. In course of probe design for DNA microarray development, each of the markers yielded a high number of potentially species-specific probes in silico, although many of them were rejected based on microarray hybridisation experiments. None of the markers provided probes to discriminate the sibling flatfish and gurnard species. However, since 16S-probes were less negatively influenced by the “position of label” effect and showed the lowest rejection rate and the highest mean signal intensity, 16S is more suitable for DNA microarray probe design than cty b and COI. The large portion of rejected COI-probes after hybridisation experiments (>90%) renders the DNA barcoding marker as rather unsuitable for this high-throughput technology.

Conclusions/Significance

Based on these data, a DNA microarray containing 64 functional oligonucleotide probes for the identification of 30 out of the 50 fish species investigated was developed. It represents the next step towards an automated and easy-to-handle method to identify fish, ichthyoplankton, and fish products.     

Application of Equilibrium Models of Solution Hybridization to Microarray Design and Analysis

Background

The probe percent bound value, calculated using multi-state equilibrium models of solution hybridization, is shown to be useful in understanding the hybridization behavior of microarray probes having 50 nucleotides, with and without mismatches. These longer oligonucleotides are in widespread use on microarrays, but there are few controlled studies of their interactions with mismatched targets compared to 25-mer based platforms.

Principal Findings

50-mer oligonucleotides with centrally placed single, double and triple mismatches were spotted on an array. Over a range of target concentrations it was possible to discriminate binding to perfect matches and mismatches, and the type of mismatch could be predicted accurately in the concentration midrange (100 pM to 200 pM) using solution hybridization modeling methods. These results have implications for microarray design, optimization and analysis methods.

Conclusions

Our results highlight the importance of incorporating biophysical factors in both the design and the analysis of microarrays. Use of the probe “percent bound” value predicted by equilibrium models of hybridization is confirmed to be important for predicting and interpreting the behavior of long oligonucleotide arrays, as has been shown for short oligonucleotide arrays.     

Systematic spatial bias in DNA microarray hybridization is caused by probe spot position-dependent variability in lateral diffusion

Background

The hybridization of nucleic acid targets with surface-immobilized probes is a widely used assay for the parallel detection of multiple targets in medical and biological research. Despite its widespread application, DNA microarray technology still suffers from several biases and lack of reproducibility, stemming in part from an incomplete understanding of the processes governing surface hybridization. In particular, non-random spatial variations within individual microarray hybridizations are often observed, but the mechanisms underpinning this positional bias remain incompletely explained.

Methodology/Principal Findings

This study identifies and rationalizes a systematic spatial bias in the intensity of surface hybridization, characterized by markedly increased signal intensity of spots located at the boundaries of the spotted areas of the microarray slide. Combining observations from a simplified single-probe block array format with predictions from a mathematical model, the mechanism responsible for this bias is found to be a position-dependent variation in lateral diffusion of target molecules. Numerical simulations reveal a strong influence of microarray well geometry on the spatial bias.

Conclusions

Reciprocal adjustment of the size of the microarray hybridization chamber to the area of surface-bound probes is a simple and effective measure to minimize or eliminate the diffusion-based bias, resulting in increased uniformity and accuracy of quantitative DNA microarray hybridization.     

Environmental barcoding reveals massive dinoflagellate diversity in marine environments

Background

Dinoflagellates are an ecologically important group of protists with important functions as primary producers, coral symbionts and in toxic red tides. Although widely studied, the natural diversity of dinoflagellates is not well known. DNA barcoding has been utilized successfully for many protist groups. We used this approach to systematically sample known “species”, as a reference to measure the natural diversity in three marine environments.

Methodology/Principal Findings

In this study, we assembled a large cytochrome c oxidase 1 (COI) barcode database from 8 public algal culture collections plus 3 private collections worldwide resulting in 336 individual barcodes linked to specific cultures. We demonstrate that COI can identify to the species level in 15 dinoflagellate genera, generally in agreement with existing species names. Exceptions were found in species belonging to genera that were generally already known to be taxonomically challenging, such as Alexandrium or Symbiodinium. Using this barcode database as a baseline for cultured dinoflagellate diversity, we investigated the natural diversity in three diverse marine environments (Northeast Pacific, Northwest Atlantic, and Caribbean), including an evaluation of single-cell barcoding to identify uncultivated groups. From all three environments, the great majority of barcodes were not represented by any known cultured dinoflagellate, and we also observed an explosion in the diversity of genera that previously contained a modest number of known species, belonging to Kareniaceae. In total, 91.5% of non-identical environmental barcodes represent distinct species, but only 51 out of 603 unique environmental barcodes could be linked to cultured species using a conservative cut-off based on distances between cultured species.

Conclusions/Significance

COI barcoding was successful in identifying species from 70% of cultured genera. When applied to environmental samples, it revealed a massive amount of natural diversity in dinoflagellates. This highlights the extent to which we underestimate microbial diversity in the environment.     

dATP/ATP, a multifunctional nucleotide, stimulates bacterial cell lysis, extracellular DNA release and biofilm development

Background

Signaling by extracellular adenosine 5′-triphosphase (eATP) is very common for cell-to-cell communication in many basic patho-physiological development processes. Rapid release of ATP into the extracellular environment from distressed or injured eukaryotic cells due to pathogens or other etiological factors can serve as a “danger signal”, activating host innate immunity. However, little is known about how or whether pathogenic bacteria respond to this “danger signal”.

Methods and Principal Findings

Here we report that extracellular dATP/ATP can stimulate bacterial adhesion and biofilm formation via increased cell lysis and extracellular DNA (eDNA) release. We demonstrate that extracellular dATP/ATP also stimulates bacterial adherence in vitro to human bronchial epithelial cells.

Conclusions and Significance

These data suggest that bacteria may sense extracellular dATP/ATP as a signal of “danger” and form biofilms to protect them from host innate immunity. This study reveals a very important and unrecognized phenomenon that both bacteria and host cells could respond to a common important signal molecule in a race to adapt to the presence of one another. We propose that extracellular dATP/ATP functions as an “inter-domain” warning signal that serves to induce protective measures in both Bacterial and Eukaryotic cells.     

Growth and asymmetry of soil microfungal colonies from "Evolution Canyon," Lower Nahal Oren, Mount Carmel, Israel

Background

Fluctuating asymmetry is a contentious indicator of stress in populations of animals and plants. Nevertheless, it is a measure of developmental noise, typically obtained by measuring asymmetry across an individual organism''s left-right axis of symmetry. These individual, signed asymmetries are symmetrically distributed around a mean of zero. Fluctuating asymmetry, however, has rarely been studied in microorganisms, and never in fungi.

Objective and Methods

We examined colony growth and random phenotypic variation of five soil microfungal species isolated from the opposing slopes of “Evolution Canyon,” Mount Carmel, Israel. This canyon provides an opportunity to study diverse taxa inhabiting a single microsite, under different kinds and intensities of abiotic and biotic stress. The south-facing “African” slope of “Evolution Canyon” is xeric, warm, and tropical. It is only 200 m, on average, from the north-facing “European” slope, which is mesic, cool, and temperate. Five fungal species inhabiting both the south-facing “African” slope, and the north-facing “European” slope of the canyon were grown under controlled laboratory conditions, where we measured the fluctuating radial asymmetry and sizes of their colonies.

Results

Different species displayed different amounts of radial asymmetry (and colony size). Moreover, there were highly significant slope by species interactions for size, and marginally significant ones for fluctuating asymmetry. There were no universal differences (i.e., across all species) in radial asymmetry and colony size between strains from “African” and “European” slopes, but colonies of Clonostachys rosea from the “African” slope were more asymmetric than those from the “European” slope.

Conclusions and Significance

Our study suggests that fluctuating radial asymmetry has potential as an indicator of random phenotypic variation and stress in soil microfungi. Interaction of slope and species for both growth rate and asymmetry of microfungi in a common environment is evidence of genetic differences between the “African” and “European” slopes of “Evolution Canyon.”     

Why Most Biomedical Findings Echoed by Newspapers Turn Out to be False: The Case of Attention Deficit Hyperactivity Disorder

Context

Because positive biomedical observations are more often published than those reporting no effect, initial observations are often refuted or attenuated by subsequent studies.

Objective

To determine whether newspapers preferentially report on initial findings and whether they also report on subsequent studies.

Methods

We focused on attention deficit hyperactivity disorder (ADHD). Using Factiva and PubMed databases, we identified 47 scientific publications on ADHD published in the 1990s and soon echoed by 347 newspapers articles. We selected the ten most echoed publications and collected all their relevant subsequent studies until 2011. We checked whether findings reported in each “top 10” publication were consistent with previous and subsequent observations. We also compared the newspaper coverage of the “top 10” publications to that of their related scientific studies.

Results

Seven of the “top 10” publications were initial studies and the conclusions in six of them were either refuted or strongly attenuated subsequently. The seventh was not confirmed or refuted, but its main conclusion appears unlikely. Among the three “top 10” that were not initial studies, two were confirmed subsequently and the third was attenuated. The newspaper coverage of the “top 10” publications (223 articles) was much larger than that of the 67 related studies (57 articles). Moreover, only one of the latter newspaper articles reported that the corresponding “top 10” finding had been attenuated. The average impact factor of the scientific journals publishing studies echoed by newspapers (17.1 n = 56) was higher (p<0.0001) than that corresponding to related publications that were not echoed (6.4 n = 56).

Conclusion

Because newspapers preferentially echo initial ADHD findings appearing in prominent journals, they report on uncertain findings that are often refuted or attenuated by subsequent studies. If this media reporting bias generalizes to health sciences, it represents a major cause of distortion in health science communication.     

Profiling dizziness in older primary care patients: an empirical study

Background

The diagnostic approach to dizzy, older patients is not straightforward as many organ systems can be involved and evidence for diagnostic strategies is lacking. A first differentiation in diagnostic subtypes or profiles may guide the diagnostic process of dizziness and can serve as a classification system in future research. In the literature this has been done, but based on pathophysiological reasoning only.

Objective

To establish a classification of diagnostic profiles of dizziness based on empirical data.

Design

Cross-sectional study.

Participants and Setting

417 consecutive patients of 65 years and older presenting with dizziness to 45 primary care physicians in the Netherlands from July 2006 to January 2008.

Methods

We performed tests, including patient history, and physical and additional examination, previously selected by an international expert panel and based on an earlier systematic review. We used the results of these tests in a principal component analysis for exploration, data-reduction and finally differentiation into diagnostic dizziness profiles.

Results

Demographic data and the results of the tests yielded 221 variables, of which 49 contributed to the classification of dizziness into six diagnostic profiles, that may be named as follows: “frailty”, “psychological”, “cardiovascular”, “presyncope”, “non-specific dizziness” and “ENT”. These explained 32% of the variance.

Conclusions

Empirically identified components classify dizziness into six profiles. This classification takes into account the heterogeneity and multicausality of dizziness and may serve as starting point for research on diagnostic strategies and can be a first step in an evidence based diagnostic approach of dizzy older patients.     

Tetradactyl footprints of an unknown affinity theropod dinosaur from the Upper Jurassic of Morocco

Background

New tetradactyl theropod footprints from Upper Jurassic (Oxfordian-Kimmeridgian) have been found in the Iouaridène syncline (Morocco). The tracksites are at several layers in the intermediate lacustrine unit of Iouaridène Formation. The footprints were named informally in previous works “Eutynichnium atlasipodus”. We consider as nomen nudum.

Methodology/Principal Findings

Boutakioutichnium atlasicus ichnogen. et ichnosp. nov. is mainly characterized by the hallux impression. It is long, strong, directed medially or forward, with two digital pads and with the proximal part of the first pad in lateral position. More than 100 footprints in 15 trackways have been studied with these features. The footprints are large, 38–48 cm in length, and 26–31 cm in width.

Conclusions/Significance

Boutakioutichnium mainly differs from other ichnotaxa with hallux impression in lacking metatarsal marks and in not being a very deep footprint. The distinct morphology of the hallux of the Boutakioutichnium trackmaker –i.e. size and hallux position- are unique in the dinosaur autopodial record to date.     

Genes but not genomes reveal bacterial domestication of Lactococcus lactis

Background

The population structure and diversity of Lactococcus lactis subsp. lactis, a major industrial bacterium involved in milk fermentation, was determined at both gene and genome level. Seventy-six lactococcal isolates of various origins were studied by different genotyping methods and thirty-six strains displaying unique macrorestriction fingerprints were analyzed by a new multilocus sequence typing (MLST) scheme. This gene-based analysis was compared to genomic characteristics determined by pulsed-field gel electrophoresis (PFGE).

Methodology/Principal Findings

The MLST analysis revealed that L. lactis subsp. lactis is essentially clonal with infrequent intra- and intergenic recombination; also, despite its taxonomical classification as a subspecies, it displays a genetic diversity as substantial as that within several other bacterial species. Genome-based analysis revealed a genome size variability of 20%, a value typical of bacteria inhabiting different ecological niches, and that suggests a large pan-genome for this subspecies. However, the genomic characteristics (macrorestriction pattern, genome or chromosome size, plasmid content) did not correlate to the MLST-based phylogeny, with strains from the same sequence type (ST) differing by up to 230 kb in genome size.

Conclusion/Significance

The gene-based phylogeny was not fully consistent with the traditional classification into dairy and non-dairy strains but supported a new classification based on ecological separation between “environmental” strains, the main contributors to the genetic diversity within the subspecies, and “domesticated” strains, subject to recent genetic bottlenecks. Comparison between gene- and genome-based analyses revealed little relationship between core and dispensable genome phylogenies, indicating that clonal diversification and phenotypic variability of the “domesticated” strains essentially arose through substantial genomic flux within the dispensable genome.     

Imaginal discs--a new source of chromosomes for genome mapping of the yellow fever mosquito Aedes aegypti

Background

The mosquito Aedes aegypti is the primary global vector for dengue and yellow fever viruses. Sequencing of the Ae. aegypti genome has stimulated research in vector biology and insect genomics. However, the current genome assembly is highly fragmented with only ∼31% of the genome being assigned to chromosomes. A lack of a reliable source of chromosomes for physical mapping has been a major impediment to improving the genome assembly of Ae. aegypti.

Methodology/Principal Findings

In this study we demonstrate the utility of mitotic chromosomes from imaginal discs of 4^th instar larva for cytogenetic studies of Ae. aegypti. High numbers of mitotic divisions on each slide preparation, large sizes, and reproducible banding patterns of the individual chromosomes simplify cytogenetic procedures. Based on the banding structure of the chromosomes, we have developed idiograms for each of the three Ae. aegypti chromosomes and placed 10 BAC clones and a 18S rDNA probe to precise chromosomal positions.

Conclusion

The study identified imaginal discs of 4^th instar larva as a superior source of mitotic chromosomes for Ae. aegypti. The proposed approach allows precise mapping of DNA probes to the chromosomal positions and can be utilized for obtaining a high-quality genome assembly of the yellow fever mosquito.     

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

Background

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal “barcoding gene”. An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

Methods and Findings

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a “barcode gap” using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

Conclusions

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a “species”.     

Oligonucleotide Frequencies of Barcoding Loci Can Discriminate Species across Kingdoms

Background

DNA barcoding refers to the use of short DNA sequences for rapid identification of species. Genetic distance or character attributes of a particular barcode locus discriminate the species. We report an efficient approach to analyze short sequence data for discrimination between species.

Methodology and Principal Findings

A new approach, Oligonucleotide Frequency Range (OFR) of barcode loci for species discrimination is proposed. OFR of the loci that discriminates between species was characteristic of a species, i.e., the maxima and minima within a species did not overlap with that of other species. We compared the species resolution ability of different barcode loci using p-distance, Euclidean distance of oligonucleotide frequencies, nucleotide-character based approach and OFR method. The species resolution by OFR was either higher or comparable to the other methods. A short fragment of 126 bp of internal transcribed spacer region in ribosomal RNA gene was sufficient to discriminate a majority of the species using OFR.

Conclusions/Significance

Oligonucleotide frequency range of a barcode locus can discriminate between species. Ability to discriminate species using very short DNA fragments may have wider applications in forensic and conservation studies.     

A gap analysis methodology for collecting crop genepools: a case study with phaseolus beans

Background

The wild relatives of crops represent a major source of valuable traits for crop improvement. These resources are threatened by habitat destruction, land use changes, and other factors, requiring their urgent collection and long-term availability for research and breeding from ex situ collections. We propose a method to identify gaps in ex situ collections (i.e. gap analysis) of crop wild relatives as a means to guide efficient and effective collecting activities.

Methodology/Principal Findings

The methodology prioritizes among taxa based on a combination of sampling, geographic, and environmental gaps. We apply the gap analysis methodology to wild taxa of the Phaseolus genepool. Of 85 taxa, 48 (56.5%) are assigned high priority for collecting due to lack of, or under-representation, in genebanks, 17 taxa are given medium priority for collecting, 15 low priority, and 5 species are assessed as adequately represented in ex situ collections. Gap “hotspots”, representing priority target areas for collecting, are concentrated in central Mexico, although the narrow endemic nature of a suite of priority species adds a number of specific additional regions to spatial collecting priorities.

Conclusions/Significance

Results of the gap analysis method mostly align very well with expert opinion of gaps in ex situ collections, with only a few exceptions. A more detailed prioritization of taxa and geographic areas for collection can be achieved by including in the analysis predictive threat factors, such as climate change or habitat destruction, or by adding additional prioritization filters, such as the degree of relatedness to cultivated species (i.e. ease of use in crop breeding). Furthermore, results for multiple crop genepools may be overlaid, which would allow a global analysis of gaps in ex situ collections of the world''s plant genetic resources.     

Genomic signatures of strain selection and enhancement in Bacillus atrophaeus var. globigii, a historical biowarfare simulant

Background

Despite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS).

Results

Archival strains and two “present day” type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the “military” isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of “military” isolates on lactate-containing media, and showed that the “military” strains exhibited a hypersporulating phenotype.

Conclusions

Our analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.     

The Effectiveness of Three Regions in Mitochondrial Genome for Aphid DNA Barcoding: A Case in Lachininae

Background

The mitochondrial gene COI has been widely used by taxonomists as a standard DNA barcode sequence for the identification of many animal species. However, the COI region is of limited use for identifying certain species and is not efficiently amplified by PCR in all animal taxa. To evaluate the utility of COI as a DNA barcode and to identify other barcode genes, we chose the aphid subfamily Lachninae (Hemiptera: Aphididae) as the focus of our study. We compared the results obtained using COI with two other mitochondrial genes, COII and Cytb. In addition, we propose a new method to improve the efficiency of species identification using DNA barcoding.

Methodology/Principal Findings

Three mitochondrial genes (COI, COII and Cytb) were sequenced and were used in the identification of over 80 species of Lachninae. The COI and COII genes demonstrated a greater PCR amplification efficiency than Cytb. Species identification using COII sequences had a higher frequency of success (96.9% in “best match” and 90.8% in “best close match”) and yielded lower intra- and higher interspecific genetic divergence values than the other two markers. The use of “tag barcodes” is a new approach that involves attaching a species-specific tag to the standard DNA barcode. With this method, the “barcoding overlap” can be nearly eliminated. As a result, we were able to increase the identification success rate from 83.9% to 95.2% by using COI and the “best close match” technique.

Conclusions/Significance

A COII-based identification system should be more effective in identifying lachnine species than COI or Cytb. However, the Cytb gene is an effective marker for the study of aphid population genetics due to its high sequence diversity. Furthermore, the use of “tag barcodes” can improve the accuracy of DNA barcoding identification by reducing or removing the overlap between intra- and inter-specific genetic divergence values.     

Shared and unique patterns of embryo development in extremophile poeciliids

Background

Closely related lineages of livebearing fishes have independently adapted to two extreme environmental factors: toxic hydrogen sulphide (H₂S) and perpetual darkness. Previous work has demonstrated in adult specimens that fish from these extreme habitats convergently evolved drastically increased head and offspring size, while cave fish are further characterized by reduced pigmentation and eye size. Here, we traced the development of these (and other) divergent traits in embryos of Poecilia mexicana from benign surface habitats (“surface mollies”) and a sulphidic cave (“cave mollies”), as well as in embryos of the sister taxon, Poecilia sulphuraria from a sulphidic surface spring (“sulphur mollies”). We asked at which points during development changes in the timing of the involved processes (i.e., heterochrony) would be detectible.

Methods and Results

Data were extracted from digital photographs taken of representative embryos for each stage of development and each type of molly. Embryo mass decreased in convergent fashion, but we found patterns of embryonic fat content and ovum/embryo diameter to be divergent among all three types of mollies. The intensity of yellow colouration of the yolk (a proxy for carotenoid content) was significantly lower in cave mollies throughout development. Moreover, while relative head size decreased through development in surface mollies, it increased in both types of extremophile mollies, and eye growth was arrested in mid-stage embryos of cave mollies but not in surface or sulphur mollies.

Conclusion

Our results clearly demonstrate that even among sister taxa convergence in phenotypic traits is not always achieved by the same processes during embryo development. Furthermore, teleost development is crucially dependent on sufficient carotenoid stores in the yolk, and so we discuss how the apparent ability of cave mollies to overcome this carotenoid-dependency may represent another potential mechanism explaining the lack of gene flow between surface and cave mollies.