Clines, clusters, and the effect of study design on the inference of human population structure

Previously, we observed that without using prior information about individual sampling locations, a clustering algorithm applied to multilocus genotypes from worldwide human populations produced genetic clusters largely coincident with major geographic regions. It has been argued, however, that the degree of clustering is diminished by use of samples with greater uniformity in geographic distribution, and that the clusters we identified were a consequence of uneven sampling along genetic clines. Expanding our earlier dataset from 377 to 993 markers, we systematically examine the influence of several study design variables—sample size, number of loci, number of clusters, assumptions about correlations in allele frequencies across populations, and the geographic dispersion of the sample—on the “clusteredness” of individuals. With all other variables held constant, geographic dispersion is seen to have comparatively little effect on the degree of clustering. Examination of the relationship between genetic and geographic distance supports a view in which the clusters arise not as an artifact of the sampling scheme, but from small discontinuous jumps in genetic distance for most population pairs on opposite sides of geographic barriers, in comparison with genetic distance for pairs on the same side. Thus, analysis of the 993-locus dataset corroborates our earlier results: if enough markers are used with a sufficiently large worldwide sample, individuals can be partitioned into genetic clusters that match major geographic subdivisions of the globe, with some individuals from intermediate geographic locations having mixed membership in the clusters that correspond to neighboring regions.     

Range-wide genetic differentiation of Eugenia dysenterica (Myrtaceae) populations in Brazilian Cerrado

The analysis of geographic patterns in genetic variation has been one of the most important current tools to understand ecological and evolutionary processes underlying population structure. However, inferring such processes from population data may be misleading if biased geographic samples are analyzed. Here we expand previous analyses of Eugenia dysenterica population structure in Brazilian Cerrado, analyzing a larger number of populations distributed throughout a broader geographic region covering most of species' range. We provide new estimates of genetic diversity and population structure based on SSR markers from both neutral and genic regions, using several cluster and ordination techniques. These analyzes reveal a continuous northwestern-southeastern gradient in population differentiation, and not two distinct clusters of populations as suggested in some previous studies. This more comprehensive analysis also reinforces that a simple process of stochastic differentiation do not explain the observed patterns. Moreover, we conclude that explanations for population differentiation may focus on why genetic diversity decreases toward southeastern populations and not necessarily on barriers and interruption of gene flow creating regional patterns of population differentiation.     

Population Structure in the Roundtail Chub (Gila robusta Complex) of the Gila River Basin as Determined by Microsatellites: Evolutionary and Conservation Implications

Ten microsatellite loci were characterized for 34 locations from roundtail chub (Gila robusta complex) to better resolve patterns of genetic variation among local populations in the lower Colorado River basin. This group has had a complex taxonomic history and previous molecular analyses failed to identify species diagnostic molecular markers. Our results supported previous molecular studies based on allozymes and DNA sequences, which found that most genetic variance was explained by differences among local populations. Samples from most localities were so divergent species-level diagnostic markers were not found. Some geographic samples were discordant with current taxonomy due to admixture or misidentification; therefore, additional morphological studies are necessary. Differences in spatial genetic structure were consistent with differences in connectivity of stream habitats, with the typically mainstem species, G. robusta, exhibiting greater genetic connectedness within the Gila River drainage. No species exhibited strong isolation by distance over the entire stream network, but the two species typically found in headwaters, G. nigra and G. intermedia, exhibited greater than expected genetic similarity between geographically proximate populations, and usually clustered with individuals from the same geographic location and/or sub-basin. These results highlight the significance of microevolutionary processes and importance of maintaining local populations to maximize evolutionary potential for this complex. Augmentation stocking as a conservation management strategy should only occur under extreme circumstances, and potential source populations should be geographically proximate stocks of the same species, especially for the headwater forms.     

Genetic structure of Pinus henryi and Pinus tabuliformis: Natural landscapes as significant barriers to gene flow among populations

Mountains as natural barriers often have important effects on intraspecific genetic structure through restraining gene flow and enhancing differentiation among populations. While the Qinling and Daba mountains are considered significant geographic barriers, dividing China into temperate and subtropical regions, little is known about how this barrier influences the genetic patterns of sister species represented in distinct habitats. In this study, we analyzed genetic differentiation and the geographic boundary between Pinus henryi and Pinus tabuliformis using chloroplast microsatellite markers. Our data show high levels of among-population differentiation, consistent with the effects of historical demographic bottlenecks, local adaptation and climate effects. Three main geographic boundaries coinciding with mountain systems indicate natural landscapes, such as large rivers, and habitat loss caused by anthropogenic deforestation, are significant barriers to genetic exchange among populations. The divergence between populations in the eastern and western Qinling Mountains populations may possibly be ascribed to fragmentation driven by climate change. The genetic boundary of P. henryi and P. tabuliformis generally coincides with the previous morphological dividing line based on the unweighted pair group method using arithmetic averages and on spatial analysis of molecular variance.     

Population genetic structure of the predatory,social wasp Vespula pensylvanica in its native and invasive range

Invasive species cause extensive damage to their introduced ranges. Ocean archipelagos are particularly vulnerable to invasive taxa. In this study, we used polymorphic microsatellite markers to investigate the genetic structure of the social wasp Vespula pensylvanica in its native range of North America and its introduced range in the archipelago of Hawaii. Our goal was to gain a better understanding of the invasion dynamics of social species and the processes affecting biological invasions. We found that V. pensylvanica showed no significant genetic isolation by distance and little genetic structure over a span of 2000 km in its native range. This result suggests that V. pensylvanica can successfully disperse across large distances either through natural‐ or human‐mediated mechanisms. In contrast to the genetic patterns observed in the native range, we found substantial genetic structure in the invasive V. pensylvanica range in Hawaii. The strong patterns of genetic differentiation within and between the Hawaiian Islands may reflect the effects of geographic barriers and invasion history on gene flow. We also found some evidence for gene flow between the different islands of Hawaii which was likely mediated through human activity. Overall, this study provides insight on how geographic barriers, invasion history, and human activity can shape population genetic structure of invasive species.     

A Spatial Framework for Understanding Population Structure and Admixture

Geographic patterns of genetic variation within modern populations, produced by complex histories of migration, can be difficult to infer and visually summarize. A general consequence of geographically limited dispersal is that samples from nearby locations tend to be more closely related than samples from distant locations, and so genetic covariance often recapitulates geographic proximity. We use genome-wide polymorphism data to build “geogenetic maps,” which, when applied to stationary populations, produces a map of the geographic positions of the populations, but with distances distorted to reflect historical rates of gene flow. In the underlying model, allele frequency covariance is a decreasing function of geogenetic distance, and nonlocal gene flow such as admixture can be identified as anomalously strong covariance over long distances. This admixture is explicitly co-estimated and depicted as arrows, from the source of admixture to the recipient, on the geogenetic map. We demonstrate the utility of this method on a circum-Tibetan sampling of the greenish warbler (Phylloscopus trochiloides), in which we find evidence for gene flow between the adjacent, terminal populations of the ring species. We also analyze a global sampling of human populations, for which we largely recover the geography of the sampling, with support for significant histories of admixture in many samples. This new tool for understanding and visualizing patterns of population structure is implemented in a Bayesian framework in the program SpaceMix.     

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar (Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad‐scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large‐scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.     

High Differentiation among Eight Villages in a Secluded Area of Sardinia Revealed by Genome-Wide High Density SNPs Analysis

To better design association studies for complex traits in isolated populations it''s important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not “a priori” be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of sub-populations and stratification within genetic isolates.     

Historical and contemporary population genetic connectivity of the European short-snouted seahorse Hippocampus hippocampus and implications for management

This first genetic study of Hippocampus hippocampus covers the species' entire geographic range and employs two mtDNA markers (control region and cytochrome b) to establish patterns of population structuring. A total of 255 specimens from 21 locations were used to obtain 89 concatenated haplotypes. The common haplotype was present in all but one population, however, most haplotypes were unique. The haplotype network had a star-like construction, suggesting expansion from a bottleneck event. F(ST) and AMOVA revealed population subdivision into three geographic regions (English Channel + Bay of Biscay, Mediterranean Sea + Atlantic Ocean Iberian coast + Macaronesian Islands, and West Africa) with barriers to gene flow indentified at Cape Finisterre and the Cape Verde frontal zone. Neutrality tests and nested clade analysis suggest a complex demographic history, with both historic events and contemporary processes shaping patterns of genetic differentiation. The genetic population subdivision detected in this study indicates that H. hippocampus should be managed as three separate units. This is especially pertinent as H. hippocampus populations within the West African region are the only ones known to be specifically targeted for exploitation.     

High interpopulation homogeneity in Central Argentina as assessed by Ancestry Informative Markers (AIMs)

The population of Argentina has already been studied with regard to several genetic markers, but much more data are needed for the appropriate definition of its genetic profile. This study aimed at investigating the admixture patterns and genetic structure in Central Argentina, using biparental markers and comparing the results with those previously obtained by us with mitochondrial DNA (mtDNA) in the same samples. A total of 521 healthy unrelated individuals living in 13 villages of the Córdoba and San Luis provinces were tested. The individuals were genotyped for ten autosomal ancestry informative markers (AIMs). Allele frequencies were compared with those of African, European and Native American populations, chosen to represent parental contributions. The AIM estimates indicated a greater influence of the Native American ancestry as compared to previous studies in the same or other Argentinean regions, but smaller than that observed with the mtDNA tests. These differences can be explained, respectively, by different genetic contributions between rural and urban areas, and asymmetric gene flow occurred in the past. But a most unexpected finding was the marked interpopulation genetic homogeneity found in villages located in diverse geographic environments across a wide territory, suggesting considerable gene flow.     

Multilocus phylogeography of the sea snake Hydrophis curtus reveals historical vicariance and cryptic lineage diversity

The Indo‐Australian archipelago (IAA) supports the world's highest diversity of marine fish, invertebrates and reptiles. Many of the marine fish and invertebrates show congruent phylogeographic patterns, supporting a view that the region's complex geo‐climatic history has played an important role in generating its exceptional biodiversity. Here, we examine population genetic structure of the viviparous sea snake, Hydrophis curtus, to assess how past and present barriers to gene flow in the IAA have contributed to genetic and species diversity in a fully marine reptile. Mitochondrial and anonymous nuclear sequences and ten microsatellite loci were used to identify patterns of historical genetic structure and population expansion, reconstruct dated genealogies and assess levels of recent gene flow. These markers revealed strong concordant geographic structure within H. curtus with a prominent genetic break between populations broadly distributed in the Indian Ocean and the West Pacific. These populations were estimated to have diverged in the late Pliocene or early Pleistocene, and microsatellite admixture analyses suggested limited recent gene flow between them despite the current lack of barriers to dispersal, indicating possible cryptic species. Subsequent divergence in the mid–late Pleistocene was detected within the West Pacific clade among the populations in the Phuket‐Thailand region, South‐East Asia and Australia, and two of these populations also showed genetic signals of recent range expansions. Our results show that climatic fluctuations during the Plio‐Pleistocene generated high levels of cryptic genetic diversity in H. curtus, and add to similar findings for diverse other marine groups in the IAA.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Population genetic structure of orchid bees (Euglossini) in anthropogenically altered landscapes

Habitat degradation and fragmentation are widespread phenomena in tropical regions. Negative effects on the biota are numerous, ranging from interruption of gene flow among populations, to the loss of genetic diversity within populations, to a decline in species richness over time. Orchid bees (Hymenoptera: Apidae: Euglossini) are of major conservation interest due to their function as pollinators of numerous orchid species and other tropical plants. Here, we used microsatellite markers to investigate the effects of geographic distance and habitat fragmentation on gene flow among populations. Populations of Euglossa dilemma in three geographic regions??the Yucat??n peninsula (Mexico), Veracruz (Mexico), and Florida (USA)??were genetically structured predominantly across the regions, with the strength of differentiation among populations being positively correlated with geographic distance. Within geographic regions only little substructure was found, suggesting that dispersal is substantial in the absence of geographic or ecological barriers. In a second study, patterns of genetic differentiation among eight species of Euglossa were not related to habitat fragmentation following deforestation in southern Mexico (Veracruz). Specifically, most bee populations in the 9,800?ha forest remnant of Los Tuxtlas (Volcano San Martin) were neither differentiated from, nor had less genetic diversity than, populations in near-continuous forest separated from Los Tuxtlas by 130?km of agricultural land. Either occasional long distance dispersal across open areas has buffered the expected genetic effects of fragmentation, or the history of fragmentation in southern Mexico is too recent to have caused measurable shifts in allelic composition.     

Major biogeographic barriers in eastern Australia have shaped the population structure of widely distributed Eucalyptus moluccana and its putative subspecies

We have investigated the impact of recognized biogeographic barriers on genetic differentiation of grey box (Eucalyptus moluccana), a common and widespread tree species of the family Myrtaceae in eastern Australian woodlands, and its previously proposed four subspecies moluccana, pedicellata, queenslandica, and crassifolia. A range of phylogeographic analyses were conducted to examine the population genetic differentiation and subspecies genetic structure in E. moluccana in relation to biogeographic barriers. Slow evolving markers uncovering long term processes (chloroplast DNA) were used to generate a haplotype network and infer phylogeographic barriers. Additionally, fast evolving, hypervariable markers (microsatellites) were used to estimate demographic processes and genetic structure among five geographic regions (29 populations) across the entire distribution of E. moluccana. Morphological features of seedlings, such as leaf and stem traits, were assessed to evaluate population clusters and test differentiation of the putative subspecies. Haplotype network analysis revealed twenty chloroplast haplotypes with a main haplotype in a central position shared by individuals belonging to the regions containing the four putative subspecies. Microsatellite analysis detected the genetic structure between Queensland (QLD) and New South Wales (NSW) populations, consistent with the McPherson Range barrier, an east‐west spur of the Great Dividing Range. The substructure was detected within QLD and NSW in line with other barriers in eastern Australia. The morphological analyses supported differentiation between QLD and NSW populations, with no difference within QLD, yet some differentiation within NSW populations. Our molecular and morphological analyses provide evidence that several geographic barriers in eastern Australia, including the Burdekin Gap and the McPherson Range have contributed to the genetic structure of E. moluccana. Genetic differentiation among E. moluccana populations supports the recognition of some but not all the four previously proposed subspecies, with crassifolia being the most differentiated.     

Mitochondrial DNA variability in viviparous and ovoviviparous populations of the urodele Salamandra salamandra

Mitochondrial DNA analysis of 13 populations of S. salamandra along a transect across the North of the Iberian Peninsula showed values of divergence between haplotypes ranging from d = 0.41% to 5.91%. Phenetic and cladistic analysis grouped the isofemale lineages into two main clusters with contrasting phylogeographic patterns. The first group encompasses populations located at each extreme of the Iberian Peninsula. Despite the large geographic distance separating these populations, they exhibit only a minor degree of divergence among haplotypes. In contrast, much higher diversification, in both number of distinct haplotypes, and overall genetic divergence, was observed in the second group of phylogenetically related populations. Surprisingly, this process of radiation and divergence in mtDNA haplotypes occurred in populations in close geographic proximity. All populations sampled in this group are located within a 300 km range, in the central part of our transect across the Northern edge of the Peninsula. Most populations in the central range of our transect exhibit viviparous reproduction — which is derived and highly unusual among urodeles. The genetic distances measured among Asturian (central portion of our transect), viviparous populations are higher than the distances measured between the two main taxonomic clusters. A viviparous population showing an unusual level of mtDNA heterogenetiy is reported and the potential implications of this focus of localized variability are discussed. The dynamics of isofemale lineages among the two reproductive modes was further explored in combination with the previous allozyme data. Several nuclear markers suggest that major mtDNA divergences could be explained by long-term extrinsic barriers to gene flow. Isofemale lineages indicate a narrow secondary contact zone among populations with different reproductive patterns. The existence of viviparous and ovovivparous populations sharing a common haplotype suggests that reproductive transition in S. salamandra could have arisen in absence of genetic mtDNA differentiation. We finally outline a genetic model system where the acquisition of water independence from a primitively aquatic dependent amphibian life cycle can be analyzed from a microevolutionary perspective.     

Estimating divergence times from DNA sequences

The patterns of genetic variation within and among individuals and populations can be used to make inferences about the evolutionary forces that generated those patterns. Numerous population genetic approaches have been developed in order to infer evolutionary history. Here, we present the “Two-Two (TT)” and the “Two-Two-outgroup (TTo)” methods; two closely related approaches for estimating divergence time based in coalescent theory. They rely on sequence data from two haploid genomes (or a single diploid individual) from each of two populations. Under a simple population-divergence model, we derive the probabilities of the possible sample configurations. These probabilities form a set of equations that can be solved to obtain estimates of the model parameters, including population split times, directly from the sequence data. This transparent and computationally efficient approach to infer population divergence time makes it possible to estimate time scaled in generations (assuming a mutation rate), and not as a compound parameter of genetic drift. Using simulations under a range of demographic scenarios, we show that the method is relatively robust to migration and that the TTo method can alleviate biases that can appear from drastic ancestral population size changes. We illustrate the utility of the approaches with some examples, including estimating split times for pairs of human populations as well as providing further evidence for the complex relationship among Neandertals and Denisovans and their ancestors.     

Patterns of Deep-Sea Genetic Connectivity in the New Zealand Region: Implications for Management of Benthic Ecosystems

Patterns of genetic connectivity are increasingly considered in the design of marine protected areas (MPAs) in both shallow and deep water. In the New Zealand Exclusive Economic Zone (EEZ), deep-sea communities at upper bathyal depths (<2000 m) are vulnerable to anthropogenic disturbance from fishing and potential mining operations. Currently, patterns of genetic connectivity among deep-sea populations throughout New Zealand’s EEZ are not well understood. Using the mitochondrial Cytochrome Oxidase I and 16S rRNA genes as genetic markers, this study aimed to elucidate patterns of genetic connectivity among populations of two common benthic invertebrates with contrasting life history strategies. Populations of the squat lobster Munida gracilis and the polychaete Hyalinoecia longibranchiata were sampled from continental slope, seamount, and offshore rise habitats on the Chatham Rise, Hikurangi Margin, and Challenger Plateau. For the polychaete, significant population structure was detected among distinct populations on the Chatham Rise, the Hikurangi Margin, and the Challenger Plateau. Significant genetic differences existed between slope and seamount populations on the Hikurangi Margin, as did evidence of population differentiation between the northeast and southwest parts of the Chatham Rise. In contrast, no significant population structure was detected across the study area for the squat lobster. Patterns of genetic connectivity in Hyalinoecia longibranchiata are likely influenced by a number of factors including current regimes that operate on varying spatial and temporal scales to produce potential barriers to dispersal. The striking difference in population structure between species can be attributed to differences in life history strategies. The results of this study are discussed in the context of existing conservation areas that are intended to manage anthropogenic threats to deep-sea benthic communities in the New Zealand region.     

Sky island bird populations isolated by ancient genetic barriers are characterized by different song traits than those isolated by recent deforestation

Various mechanisms of isolation can structure populations and result in cultural and genetic differentiation. Similar to genetic markers, for songbirds, culturally transmitted sexual signals such as breeding song can be used as a measure of differentiation as songs can also be impacted by geographic isolation resulting in population‐level differences in song structure. Several studies have found differences in song structure either across ancient geographic barriers or across contemporary habitat barriers owing to deforestation. However, very few studies have examined the effect of both ancient barriers and recent deforestation in the same system. In this study, we examined the geographic variation in song structure across six populations of the White‐bellied Shortwing, a threatened and endemic songbird species complex found on isolated mountaintops or “sky islands” of the Western Ghats. While some sky islands in the system are isolated by ancient valleys, others are separated by deforestation. We examined 14 frequency and temporal spectral traits and two syntax traits from 835 songs of 38 individuals across the six populations. We identified three major song clusters based on a discriminant model of spectral traits, degree of similarity of syntax features, as well as responses of birds to opportunistic playback. However, some traits like complex vocal mechanisms (CVM), relating to the use of syrinxes, clearly differentiated both ancient and recently fragmented populations. We suggest that CVMs may have a cultural basis and can be used to identify culturally isolated populations that cannot be differentiated using genetic markers or commonly used frequency‐based song traits. Our results demonstrate the use of bird songs to reconstruct phylogenetic groups and impacts of habitat fragmentation even in complex scenarios of historic and contemporary isolation.     

Desmophyllum dianthus (Esper, 1794) in the Scleractinian Phylogeny and Its Intraspecific Diversity

The cosmopolitan solitary deep-water scleractinian coral Desmophyllum dianthus (Esper, 1794) was selected as a representative model species of the polyphyletic Caryophylliidae family to (1) examine phylogenetic relationships with respect to the principal Scleractinia taxa, (2) check population structure, (3) test the widespread connectivity hypothesis and (4) assess the utility of different nuclear and mitochondrial markers currently in use. To carry out these goals, DNA sequence data from nuclear (ITS and 28S) and mitochondrial (16S and COI) markers were analyzed for several coral species and for Mediterranean populations of D. dianthus. Three phylogenetic methodologies (ML, MP and BI), based on data from the four molecular markers, all supported D. dianthus as clearly belonging to the “robust” clade, in which the species Lophelia pertusa and D. dianthus not only grouped together, but also shared haplotypes for some DNA markers. Molecular results also showed shared haplotypes among D. dianthus populations distributed in regions separated by several thousands of kilometers and by clear geographic barriers. These results could reflect limited molecular and morphological taxonomic resolution rather than real widespread connectivity. Additional studies are needed in order to find molecular markers and morphological features able to disentangle the complex phylogenetic relationship in the Order Scleractinia and to differentiate isolated populations, thus avoiding the homoplasy found in some morphological characters that are still considered in the literature.     

Reconciling patterns of genetic variation with stream structure, earth history and biology in the Australian freshwater fish Craterocephalus stercusmuscarum (Atherinidae)

We examined the consequences of barriers, stream architecture and putative dispersal capability on levels of genetic differentiation among populations of the freshwater fish Craterocephalus stercusmuscarum. Seven polymorphic allozyme loci and sequences of a 498-bp fragment of the ATPase 6 mitochondrial DNA (mtDNA) gene were used to assess patterns of genetic variation among 16 populations from upland and lowland streams of five drainages in northern Queensland, Australia. Concordant patterns at both genetic markers revealed that there were significant levels of genetic subdivision among all populations, while an analysis of molecular variation showed that the distribution of genetic diversity was not consistent with contemporary drainage structure. There were reciprocally monophyletic mtDNA clades and fixed or large frequency differences at allozyme loci either side of instream barriers such as waterfalls. This implied barriers were effective in restricting gene flow between upland and lowland populations separated by waterfalls. However, there were two genetically distinct groups in upland areas, even within the same subcatchment, as well as high levels of genetic subdivision among lowland populations, suggesting barriers alone do not explain the patterns of genetic diversity. The data revealed a complex phylogeographic pattern, which we interpreted to be the result of one or more invasion events of independent lineages to different sections of each drainage, possibly mediated by well documented geomorphological changes. Our results highlight the importance of earth structure and history in shaping population genetic structure in stream organisms where dispersal capability may be limited, and reveal that the contemporary structure of drainages is not necessarily a good indicator of genetic relationships among populations.