Consequences of a catadromous life-strategy for levels of mitochondrial DNA differentiation among populations of the Australian bass, Macquaria novemaculeata

The influence of a catadromous life-strategy on levels of spatial genetic structuring in fish is poorly understood. In an effort to gain a better appreciation of how this specialized life-strategy determines population genetic structuring, we assessed variation in the mitochondrial DNA (mtDNA) control region in a catadromous perciform, the Australian bass Macquaria novemaculeata . Nineteen putative haplotypes were resolved using temperature gradient gel electrophoresis from 10 geographically distinct populations. Significant heterogeneity was revealed in haplotype frequencies and their spatial distributions among many locales. Gene partitioning statistics ( AMOVA ) for both raw haplotype frequency data and frequency data with sequence divergences were concordant, indicating that M. novemaculeata populations were moderately genetically structured (Φ_ST = 0.05, 0.06; P < 0.001, respectively). Isolation by distance seems to be a strong structuring force in M. novemaculeata , culminating in no detectable phylogeographic structuring among haplotypes. Low sequence divergences were observed among many haplotypes and it is suggested that these are the result of pruning of maternal lineages by cyclical variations in female reproductive success. This study highlights the importance of life-history patterns and, in particular, spawning locality, in determining spatial structuring of mtDNA variation in catadromous species.     

Population genetic structure of the catadromous Australian bass from throughout its range

One hatchery-derived impoundment and nine wild riverine populations of the catadromous percoid Australian bass Macquaria novemaculeata were examined at six polymorphic allozyme loci to reveal population genetic structure. Subtle genetic differences among many populations indicated that fish did not interact as a panmictic unit, but rather conformed to an isolation by distance pattern of population structure. The hatchery-derived Glenbawn Dam sample differed from the riverine populations in both average heterozygosity and number of polymorphic loci, demonstrating the need for caution in supplementation programs to rehabilitate wild riverine populations.     

Morphological variation in the catadromous Australian bass, from seven geographically distinct riverine drainages

The Australian bass Macquaria novemaculeata is a catadromous species restricted to coastal drainages of south-eastern Australia and recently has undergone reductions in population numbers. Morphometric variation among bass from seven distinct drainages was used to test for the presence of stock structuring. Males and females were treated separately. For both male and female bass, significant heterogeneity in morphology among sample locales was revealed by univariate statistics and discriminant function analysis, with morphological variation found to differ in a clinal pattern. The demonstration of population-specific morphological traits among bass from different river drainages may indicate that inter-populational movement is restricted and that each sample locale constitutes a separate stock of bass.     

Distinguishing between two sympatric Acanthopagrus species from Dapeng Bay, Taiwan, using morphometric and genetic characters

Morphometric and genetic data were used to compare two sympatric and morphologically similar species, Acanthopagrus berda and Acanthopagrus taiwanensis , in Dapeng Bay, South-western Taiwan. A principle component analysis of morphological data indicated a distinction between the two species, with pectoral fin length and eye diameter accounting for 32·27% of the variation. Interspecific sequence divergence, based on mtDNA cytochrome b (0·118 ± 0·01), was larger than intraspecific divergences between haplotypes (0·007 for A. taiwanensis and 0·003 for A. berda ). Individuals of the two species clustered into different groups in three phylogenetic trees with 100% bootstrap support. The mean observed heterozygosity for eight microsatellite loci was 0·471 ± 0·202 for A. taiwanensis and 0·637 ± 0·145 for A. berda . Nei's unbiased measure of interspecific genetic distance ( D _S) was 1·334. F _ST (0·134) and R _ST (0·404) values indicated significant differentiation between species. An unrooted neighbour-joining tree was constructed by allele-sharing distances and the factorial correspondence analysis split all specimens into two distinct clusters. The results of morphometric, mtDNA and microsatellite analyses indicated the presence of two species, A. taiwanensis and A. berda .     

Population structure and phylogeography of Aphyocypris kikuchii (Oshima) based on mitochondrial DNA variation

Aphyocypris kikuchii is a cyprinid species endemic to northern and eastern Taiwan and is the only primary freshwater fish native east of the Coastal Mountain Range. In total, 92 individuals of A. kikuchii from seven populations in three regions of the island were surveyed for mitochondrial DNA (mtDNA) variation. High haplotype diversity ( h = 0·989) and low nucleotide diversity ( π = 0·009) of mtDNA were detected. Negative values of Tajima's D and unimodal mismatch distributions probably reflect a history of recent demographic expansions from small populations. Three major haplotype clusters displayed geographically non-overlapping distributions, indicating a long-term isolation between regions. Hierarchical analysis of molecular variance showed significant genetic structuring among populations ( Φ _ST= 0·66). Significant haplotype heterogeneity was also detected among populations within regions ( Φ _SC= 0·41, P < 0·001) and among regions ( Φ _CT= 0·43, P < 0·05). Molecular clock estimates of coalescence in the three major mtDNA lineages indicated coalescence in the most recent common ancestor c. 0·11–0·39 million years ago. Haplotypes of cluster B nested as interior nodes in the haplotype network, indicating that migrations from Shueilian (SL) populations to the northern region (cluster A) and to the eastern region (cluster C) may have occurred independently. Lineages A and B + C should be managed as two distinct evolutionarily significant units, while the northern, SL and southern groups should be managed as separate management units.     

Genetic population structure of red snappers (Lutjanus malabaricus Bloch & Schneider, 1801 and Lutjanus erythropterus Bloch, 1790) in central and eastern Indonesia and northern Australia

The genetic population structure of red snapper Lutjanus malabaricus and Lutjanus erythropterus in eastern Indonesia and northern Australia was investigated by allozyme electrophoresis and sequence variation in the control region of mtDNA. Samples were collected from eight sites in Indonesia and four sites in northern Australia for both species. A total of 13 allozyme loci were scored. More variable loci were observed in L. malabaricus than in L. erythropterus . Sequence variation in the control region (left domain) of the mitochondrial genome was assessed by RFLP and direct sequencing. MtDNA haplotype diversity was high ( L. erythropterus , 0·95 and L. malabaricus , 0·97), as was intraspecific sequence divergence, ( L. erythropterus , 0·0–12·5% and L. malabaricus , 0·0–9·5%). The pattern of mtDNA haplotype frequencies grouped both species into two broad fisheries stocks with a genetic boundary either between Kupang and Sape ( L. malabaricus ) or between Kupang and Australian Timor Sea ( L. erythropertus ). The allozyme analyses revealed similar boundaries for L. erythropterus . Seven allozymes stocks compared to two mtDNA stocks of L. malabaricus including Ambon, which was not sampled with mtDNA, however, were reported. Possible reasons for differences in discrimination between the methods include: i) increased power of multiple allozyme loci over the single mtDNA locus, ii) insufficient gene sampling in the mtDNA control region and iii) relative evolutionary dynamics of nuclear (allozyme loci) and mitochondrial DNA in these taxa. Allozyme and haplotype data did not distinguish separate stocks among the four Australian locations nor the central Indonesian (Bali and Sape locations) for both L. malabaricus and L. erythropterus.     

Microsatellite variation in cod Gadus morhua throughout its geographic range

Previous genetic studies using neutral markers such as allozymes, mtDNA and minisatellite loci have demonstrated varying amounts of population structure in cod Gadus morhua throughout the Atlantic. Microsatellite loci, which are potentially the most informative of presently available neutral genetic markers, have been applied extensively within western and eastern Atlantic areas but not on a range-wide basis. In the present study, six microsatellite DNA loci were used to screen cod samples from nine locations throughout the geographic range from the Scotian Shelf in the West Atlantic to the Barents and Baltic Seas in the east. Overall F _ST value was 0·03 ( P = < 0·001) across all samples. Statistically significant population differences over all loci combined were evident between more geographically distant samples, using either heterogeneity tests or F _ST analysis, with at least one locus showing significant differences between all samples (prior to Bonferroni correction). A significant correlation was observed between genetic and geographical distance, suggesting a higher level of historical and contemporary gene flow between adjacent populations than more distant populations. Samples from either end of the geographic range (Scotian Shelf and Baltic Sea) were particularly distinct when analysed using the STRUCTURE programme and also showed a high level of self-assignment when individuals of either the Scotian Shelf or Baltic Sea were tested against the entire data set. The present microsatellite study demonstrates a high level of geographic population structure between the western Atlantic, middle and eastern Atlantic and Baltic Sea, and thus, the findings should be useful in devising overall management and conservation strategies for the species.     

Hierarchical organization of genetic variation in the Costa Rican livebearing fish Brachyrhaphis rhabdophora (Poeciliidae)

I examined the geographic distribution of genetic variation in the livebearing freshwater fish Brachyrhaphis rhabdophora in northwestern Costa Rica as revealed by allozymes and mitochondrial DNA sequences. Allelic variability at 11 enzyme-coding loci surveyed across 12 localities revealed marked genetic differentiation among populations within drainages (0_P= 0.36) and among drainages within regions (0_D=0.17), but not between northern and southern geographic regions (0_R=– 0.02). Allozyme variation was hierarchically organized such that populations found within stream drainages were more similar to each other than to populations found in adjacent drainages, a result confirmed by cluster analysis. In contrast to the allozyme data, there was extremely little DNA sequence variation among populations in the mitochondrial control region (3 variable nucleotide positions out of 444 bp examined). The difference in genetic divergence between allozyme and mtDNA markers was unexpected and is discussed in terms of biogeographical colonization events and a molecular selective swéep on the mitochondrial genome, both processes that could explain the lack of mitochondrial variability in this highly subdivided species.     

A hybrid zone and bidirectional introgression between two catadromous species: Australian bass Macquaria novemaculeata and estuary perch Macquaria colonorum

The presence and distribution of hybrid individuals and the existence of a hybrid zone between the catadromous Australian bass Macquaria novemaculeata and estuary perch Macquaria colonorum were investigated throughout the range of both species in Australia. Bayesian analyses and genotypic simulations identified 140 putative hybrids (11·5% of the total sample) with varying levels of introgression. Most hybrids were observed in an area extending from the Snowy River to the Albert River suggesting a hybrid zone in the eastern Bass Strait region. Sixteen hybrids, however, were found outside this zone, possibly reflecting the movement of hybrid offspring between estuaries or their inadvertent release during fish stocking programmes. Biparental backcrossing was found to occur suggesting that hybrids were fertile. These results have implications for the management of the extensive stocking programme in M. novemaculeata and for understanding the potential role of habitat degradation and reduced water flow in facilitating hybridization in species with migratory life histories.     

Genetic variation and stock structure of school mackerel and spotted mackerel in northern Australian waters

The total mean sample heterozygosity calculated from eight polymorphic loci was 0·172 (0·047 S.E.) (F_st 0·025) for school mackerel Scomberomorus queenslandicus , and for spotted mackerel S. munroi was 0·110 (0·074 S.E.) (F_st 0·038). There was no evidence of temporal variation for either species as significant genetic differences were not detected between months or year classes within areas. Spatially, school mackerel have a complex stock structure, with stocks being associated with large embayments. In contrast, spotted mackerel appear to comprise a single stock in Australian east coast waters. Both species showed a significant pattern of stock structure between Australian east coast and northern (Arafura Sea) samples.     

Genetic differentiation among populations of the shortfinned eel Anguilla australis from East Australia and New Zealand

Variability at seven microsatellite loci was used to survey the genetic population structure of the shortfinned eel Anguilla australis . Samples were collected from six estuaries along the east coast of Australia and from three estuaries around New Zealand. Hierarchical analysis of molecular variance of the five loci with good fit to Hardy–Weinberg genotypic proportions detected highly significant differences among samples ( F _ST= 0·016, P < 0·001). The fixation index between countries ( F _CT= 0·012, P < 0·001) was more than double the index among samples within countries ( F _SC= 0·005, P < 0·05). An unweighted pair-group method with arithmetic mean (UPGMA) tree also supported the separation of Australian and New Zealand populations, as did assignment tests, which correctly assigned 80 and 84% of the individuals to Australia and New Zealand, respectively. Isolation-by-distance appeared among samples overall ( r = 0·807, P < 0·001), but not among samples within countries ( r = 0·027, P > 0·05 in Australia; r = 0·762, P > 0·05 in New Zealand). These findings indicate that populations of A. australis in East Australia and in New Zealand may be reproductively isolated from one another. Genetic differentiation among populations of A. australis was two- to 10-fold higher than that among populations of other temperate eels in the North Atlantic Ocean, suggesting that two group of A. australis may reflect sub-species. Anguilla australis in the two countries have different genetic structures and thus require separate management. Genetic isolation between Australian and New Zealand populations indicates that juveniles recruit independently into these two regions from geographically or temporally isolated spawning areas.     

Population genetic structure and effective population size of ayu (Plecoglossus altivelis), an amphidromous fish

The temporal and spatial population genetic structure of ayu Plecoglossus altivelis (Salmoniformes: Plecoglossidae), an amphidromous fish, was examined using analysis of variation at six microsatellite DNA loci. Intracohort genetic diversities, as measured by the number of alleles and heterozygosity, were similar among six cohorts (2001–2006) within a population (Nezugaseki River), with the mean number of alleles per cohort ranging from 11·0 to 12·5 and the expected heterozygosity ranging from 0·74 to 0·77. Intrapopulational genetic diversities were also similar across the three studied populations along the 50 km coast, with the mean number of alleles and the expected heterozygosity ranging from 11·33 to 11·67 and from 0·75 to 0·76, respectively. The authors observed only one significant difference in pair-wise population differentiation ( F _ST-value) between the cohorts within a population and among three populations. Estimates of the effective population size ( N _e) based on maximum-likelihood method yielded small values (ranging from 94·8 to 135·5), whereas census population size ranged from c. 4800 to 24 000. As a result, the ratio of annual effective population sizes to census population size ( N _e/ N ) ranged from 0·004 to 0·023. These estimates of N _e/ N agree more closely with estimates for marine fishes than that of the larger estimates for freshwater fishes. The present study suggests that ayu which is highly fecund and shows low survival during the early life stages is also characterized by having low value of N _e/ N , similar to marine species with a pelagic life cycle.     

Generalist niche, specialist strategy: the diet of an Australian percichthyid

Dietary analysis revealed that an impoundment population of Australian bass Macquaria novemaculeata holds a generalist niche, but one arising from persistent individual specialization and interindividual variation. This 'individual specialist' strategy appeared adaptive, but the strength of individual specialization was largely independent of variation in diet composition, except during blooms of Daphnia sp. Diet composition and dietary overlap showed only moderate ontogenetic variation, and niche breadth showed no relationship with ontogeny. Macquaria novemaculeata showed an asymmetric predator and prey size distribution, consistent with many aquatic predators, with positive relationships between fish size and average, maximum and minimum prey size. There was no asymmetry in the relative size-based niche breadths of individuals, however, which indicates that the niche is a fixed 'window' of relative prey sizes. The difference in the dietary niche and prey-size relationships of M. novemaculeata at the population and individual levels highlights the necessity of assessing the niche at both these levels.     

Genetic variation, its assessment and implications to the conservation of seagrasses

In a study of the widespread Australian endemic seagrass Posidonia australis , allozyme analysis identified a wide range in population genetic structure assessed using the multilocus genotype diversity statistic ( D _G). Values of D _G between zero and one were obtained; however, RAPD analysis generally detected higher levels of diversity, where D _G values were all greater than 0.5 ( D _G = 0.67 – 1). Some populations were allozymically monomorphic using allozyme analysis yet were highly polymorphic using RAPD analysis. The differences observed between methods, particularly among allozymically uniform populations, demonstrate the importance of choosing an appropriate method when assessing genotypic diversity. Different methods may reflect different historical aspects of population processes where allozymes reflect broader-scale gene flow and population establishment and DNA fingerprinting methods such as RAPDs may reflect fine-scale local recruitment events and shorter-term population processes. Using either method alone, particularly in genotypically depauperate organisms such as seagrasses and other clonal organisms, will be problematic in assessing their population genetic potential, a parameter being used by conservation managers to decide upon management strategies in rare and endangered organisms. It is recommended that the impact of disturbance assessed using genotypic diversity measures requires more than one technique to provide the most appropriate information for designing subsequent conservation strategies.     

Genetic diversity and effective size of the Atlantic salmon Salmo salar L. inhabiting the River Eo (Spain) following a stock collapse

Historic angling records suggest the occurrence of a drastic decline in the River Eo (Spain) Atlantic salmon population size during the past two decades, as a result of overexploitation and habitat deterioration. In recent years, the population has been apparently recovering, and the present study is aimed to report information on the level of genetic diversity and the effective size of the current population as these may have immediate consequences for its conservation. Eighty-six salmon from two temporal groups (1998–1999 and 2004–2006), representing three generations, were genotyped using a panel of eight microsatellites. Inspite of the recent decline in census numbers and the detection of the signs of a population bottleneck, the population exhibits a high level of genetic diversity, similar to that from other populations, and almost unchanged during the period of study [average allelic richness ( A ) = 14·0 and 13·9, and average heterozygosity ( H _e) = 0·843 and 0·851 in 1998–1999 and 2004–2006, respectively]. The effective population size ( N _e) estimated by two different temporal methods showed a consistent value around 80 salmon, whereas the estimates from the linkage disequilibrium (LD) method provided a value around 165 individuals for either sample. The recent growing number of salmon, as indicated by fisheries records, the relatively large estimates of the ratio N _e/ N (with range 0·23–0·44 for the temporal estimates and 0·31–0·59 for the LD estimates) and the high levels of diversity found suggest that the population has not been greatly affected by the historical census declines and can be expected to recover in the future.     

Indo-West Pacific patterns of genetic differentiation in the high-dispersal starfish Linckia laevigata

Genetic variation in four natural populations of the starfish Linckia laevigata from the Indo-West Pacific was examined using restriction fragment analysis of a portion of the mtDNA including the control region. Digestion with seven restriction enzymes identified 47 haplotypes in a sample of 326 individuals. Samples collected from reef sites within each location were not significantly differentiated based on Φ_ST or spatial distribution of haplotypes, indicating that dispersal is high over short to moderate distances. Evidence of gene flow is further supported by the low divergence among haplotypes and the lack of any clear geographical structuring among different haplotypes in the gene phylogeny. However, analysis of molecular variance ( AMOVA ), Φ_ST and contingency χ² analyses of the spatial distribution of haplotypes demonstrate the presence of significant broad scale population genetic structure among the four widespread locations examined. RFLP data are consistent with high gene flow between the Philippines and Western Australia and moderate gene flow between the Great Barrier Reef (GBR) and Fiji, but only limited gene flow between either the Philippines or Western Australia and either the GBR or Fiji. The presence of mtDNA structure contrasts with previous allozyme data which suggest that dispersal among widely separated locations is equivalent to dispersal among populations within the highly connected GBR studies. This discordance between patterns of gene flow inferred from these two markers cannot be fully accounted for by differences in effective population size for mtDNA. This might suggest that while mtDNA variation may represent contemporary patterns of gene flow, allozyme variation among populations is yet to reach equilibrium between drift and migration over the range surveyed.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Molecular genetic assessment of bipolarity in the anchovy genus Engraulis

To test three hypotheses accounting for bipolar distributions in Engraulis , seven of eight taxa (except E. eurystole ) were surveyed with allozymes (34 loci) and 521 bp of the mitochondrial DNA cytochrome b gene. Both allozymes and mtDNA sequences revealed deep separations between New World and Old World anchovies with a mean allozyme genetic distance D  = 1·26 and net mtDNA sequence divergence d  = 15%. These values reflected separations of 5 to 10 million years. Contrary to previous phylogenetic hypotheses, which place north‐east Pacific E. mordax and south‐east Pacific E. ringens as sister taxa, the south‐west Atlantic E. anchoita and E. ringens are most closely related to each other. The north‐east Pacific E. mordax is the closest lineage to the Old World Engraulis , a group of taxa showing low divergences typical of population‐level separations (mean D  = 0·06; mean d  = 0·87%). Bipolarities of sister taxa in the east Atlantic and west Pacific reflect recent dispersals. Bipolarities in the east Pacific and west Atlantic represent paraphyletic taxa in lineages isolated since the Miocene. None of the bipolarities can be attributed to tectonic separations or competitive displacements from the tropics, but the latter situation should be re‐evaluated with comparisons to tropical anchovies of the New World.     

Morphological and genetic differentiation in anadromous smelt Osmerus mordax (Mitchill): disentangling the effects of geography and morphology on gene flow

Morphological analyses were combined with genetic analyses at nine microsatellite loci to examine the determinants of gene flow at 21 spawning locations of rainbow smelt Osmerus mordax along the east coast of Canada. Associations between morphology, geography and gene flow were examined using a computational geometric approach and partial Mantel tests. Significant barriers to gene flow and discontinuities in morphology were observed between Newfoundland and mainland Canada, as well as within Newfoundland samples. On regional scales, contrasting patterns were present with restricted gene flow between Newfoundland populations ( F _ST= c . 0·11) and high gene flow between mainland populations ( F _ST= c . 0·017). Within Newfoundland populations, geographic distance was significantly associated with gene flow ( r = 0·85, P < 0·001) contrasting mainland samples where gene flow was most associated with phenotypic divergence ( r = 0·33, P < 0·001). At large spatial scales, weak ( r = 0·19, P = 0·02) associations between gene flow and geographic distance were observed, and moderate associations were also observed between gene flow and morphology ( r = 0·28, P < 0·001). The presence of significant genetic isolation by distance in Newfoundland samples and the clear discontinuity associated with the Cabot Strait suggest geography may be the primary determinant of gene flow. Interestingly, the association between genetic and morphological divergence within mainland samples and overall, supports the hypothesis that gene flow may be moderated by morphological divergence at larger spatial scales even in high gene flow environments.     

Mitochondrial DNA phylogeography of Glyptothorax fokiensis and Glyptothorax hainanensis in Asia

Whole mitochondrial DNA cytochrome b sequences in 62 fish from 13 locations in Southeast China identified two major clades corresponding to two allopatric taxa, Glyptothorax fokiensis fokiensis and Glyptothorax fokiensis hainanensis . Reciprocal monophyly and a molecular clock separation between these two taxa of 2·3 million years indicate these taxa should be elevated to species. Mismatch distributions and Fu's F _S statistic suggest that both G. fokiensis and G. hainanensis have experienced recent population expansions. Analysis of molecular variance indicates that most of the genetic variation resides among populations within both species, with Φ _ST= 0·645 for G. fokiensis and 0·801 for G. hainanensis , suggesting restricted gene flow among populations. Significant correlations between the geographic and the genetic distances provide support for the importance of geographic isolations between populations. Nested clade analysis also confirms low levels of genetic exchanges between the two major groups and between populations within each group. The phylogeographical pattern among populations of Glyptothorax in East Asia can be attributed to historical fragmentations, demographic expansions and occasional long-distance dispersals stimulated by tectonic activity and Ice Age climate changes.