Genetic structure of the killifish Aphanius fasciatus, Nardo 1827 (Teleostei, Cyprinodontidae), results of mitochondrial DNA analysis

Aphanius fasciatus is a cyprinodont distributed in the salty coastal water of the central and eastern Mediterranean Sea and occasionally in internal fresh water. In this work, the authors have investigated the genetic structure of eight populations of the killifish A. fasciatus from Sardinia and Sicily. The comparison of the mtDNA control region of 237 individuals revealed a total of 49 haplotypes. Several unique haplotypes were present in each population, and no common haplotype was found among Sicilian and Sardinian populations. Almost all Sardinian populations shared a common haplotype, and indeed the four Sicilian populations examined did not share any as determined by the parsimony network analysis. The analysis of molecular variance showed that the percentage of variation among populations is much higher than within each population of A. fasciatus . The overall F _ST value is very high (0·78) and supports an extensive genetic structure of the populations. The observed genetic differentiations of A. fasciatus populations were discussed taking into account the palaeogeographic and palaeoclimatic events that interested the Mediterranean area from Miocenic to Pleistocenic age. The results provide new insight into the knowledge of the pattern of genetic structure and of evolutionary processes occurring in this species.     

Effects of a founder event and supplementary introductions on genetic variation in a captive breeding population of the endangered Spanish killifish

Twelve polymorphic allozyme loci were employed to assess the genetic change in a captive breeding population of the endangered killifish Aphanius baeticus in the Doñana National Park, south‐western Spain. The initial founder event did not significantly reduce the allelic richness or the expected heterozygosity. No genetic bottleneck signature was detected by tests for deviation from mutation‐drift equilibrium. The F _ST between the wild source and captive population, however, was relatively high (0·053 or 0·122 when excluding or including the locus IDHP‐1 * respectively), after just two to three generations in captivity. Two generations after the incorporation of 68 new wild specimens (greater than five generations after founding) decreased the genetic differences and the F _ST(0·041 excluding IDHP‐1 *). The restoration efforts appeared to be helpful and the study of 12 polymorphic loci and a sensitive parameter such as F _ST were useful for monitoring genetic changes in captivity. Nonetheless, future monitoring should include additional highly polymorphic loci (microsatellites) to achieve higher power to detect genetic change. Such restoration and monitoring efforts should help to avoid rapid inbreeding, adaptation to captivity, and to maintain the long‐term evolutionary potential in small isolated populations.     

Does the Mediterranean killifish Aphanius fasciatus (Teleostei: Cyprinodontidae) fit the one-dimensional stepping-stone model of gene flow?

Five samples (30 individuals each) of the Mediterranean cyprinodontid fish, Aphanius fasciatus, were analysed by allozyme electrophoresis for a study concerning genetic structure and gene flow among populations. Forty-three loci, twelve of which were polymorphic at P_0.99 level, were analysed. A. fasciatus showed low levels of genetic polymorphism, with expected heterozygosity values ranging from 0.027 (SE = 0.013) to 0.064 (SE = 0.023). The modified Rogers' genetic distances ranged from 0.053 to 0.202. The observed pattern of genetic differentiation among populations, obtained by multidimensional scaling of the modified Rogers' genetic distances, was consistent with geographic distribution. Coancestry coefficient ( = 0.302, SE = 0.045) indicated a very high degree of genetic subdivision among populations and the estimate of gene flow (the effective number of migrants per generation, Nm = 0.445), computed by Wright's method, showed very restricted gene flow among populations. The regression analysis carried out using log-transformed values of geographic distance against M^ (pairwise values of Nm) gave a regression coefficient, b = –1.03, indicating that a restricted migration rate among populations can occur according to the one-dimensional stepping-stone model. The amount and modality of gene flow in A. fasciatus inferred from the results of the present study are consistent with the absence of dispersal stages and with the fragmented nature of brackish-water habitats along the coasts.     

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.     

Regional differentiation of North American Atlantic salmon at allozyme loci

Allozyme variation was characterised by starch gel electrophoresis at 23 enzyme coding loci and one regulatory locus in Atlantic salmon from 53 rivers in Eastern Canada, encompassing the majority of the species' North American range. Variation among rivers was highly heterogeneous and eight of the 15 polymorphisms showed regionally restricted distributions. Nearest neighbour joining (NJ) analysis and multi‐dimensional scaling suggest six distinct regional groups; Labrador/Ungava, Gulf of Saint Lawrence, Newfoundland (excluding Gulf rivers), the Atlantic shore/Southern Uplands of Nova Scotia, the inner Bay of Fundy, and the outer Bay of Fundy. Approximately 25% of observed genetic variation was distributed among these regions with a weak though significant overall correlation of genetic and geographic distance (Mantel Test, r  = 0·255, P  = 0·005). Collectively, the rivers showed consistent divergence from European populations with strong bootstrap support for the two clusters across loci in the NJ analysis. Mean heterozygosity was 0·061 for both continental groups, but the European population showed more than twice the variation among populations. F _ST values were 0·076 and 0·176 for North America and Europe, respectively, with an overall F _ST of 0·330.     

Isolation by distance and selection effects on genetic structure of sardines Sardina pilchardus Walbaum

Allozyme data were used to analyse the genetic structure of Sardina pilchardus populations. Fifty samples from 15 locations between the North Sea and Mauritania, including samples from the Azores, Madeira and the Mediterranean Sea, were surveyed. A weak but significant structure was found between all samples ( F _ST= 0.057, P < 0·001). This structure results from a change in the most common allele of SOD* between the North African and the Azores populations separated by the greatest distance. This locus seemed to be under selective pressure according to the test of neutrality, and the variations in allele frequencies may be explained due to isolation by distance (IBD) of coastal populations (from Mauritania to the North Sea) ( r = 0·86, P < 0·001). When SOD* was removed from the analyses, IBD was not observed in coastal populations ( r = 0·236, P > 0·05) but was observed over the whole range ( r = 0·321, P = 0·05). The genetic structure of S. pilchardus is driven by both IBD and selective processes.     

Microsatellite data reveals weak population substructuring in Copadichromis sp.'virginalis kajose', a demersal cichlid from Lake Malawi, Africa

Small but significant differences were found in allele frequencies among five populations (overall F _ST estimate (θ)=0·004, P=0·006; overall R _ST estimate (RHO)=0·019, P <0·00001) of the demersal cichlid Copadichromis sp.'virginalis kajose', collected from five locations in Lake Malawi. Pairwise F _ST estimates revealed significant differences between the most southerly population (Cape Maclear), and the three most northerly populations (Mbamba Bay, Metangula and Chilola). Pairwise R _ST estimates also revealed significant differences between some populations, but no geographical pattern was discernible. There was no evidence of isolation by distance using either the shortest straight-line distance between samples, or the distance around the shoreline following a 50 m depth contour. F _ST estimates were considerably lower than found in previous studies on the mbuna (rock-dwelling species), but higher than those found in a study of three pelagic cichlid species from Lake Malawi. Substructuring in C. sp.'virginalis kajose' appears to be on a similar scale to the Atlantic cod.     

Loss of genetic variation in hatchery-reared Indian major carp, Catla catla

The hypothesis that effective population sizes are low in hatchery-reared catla ( Catla catla ) from Bangladesh, possibly leading to inbreeding and loss of variation, was tested. The study was based on analysis of seven microsatellite loci in three samples of hatchery-reared catla and four samples representing wild populations. Pair-wise estimates of genetic differentiation between samples were low between wild samples (θ ranging from 0·012 to 0·034), but high between hatchery samples (θ ranging from 0·153 to 0·185), suggesting strong genetic drift in hatcheries. Genetic variation, both in terms of expected heterozygosity and allelic richness, was significantly lower in hatchery samples than in samples of wild catla. Application of a method for reconstructing families among offspring without parental genetic data showed that the hatchery samples consisted of very few half- and full-sib families, whereas the wild samples consisted of a high number of families, suggesting that most individuals were unrelated. Finally, estimation of the effective number of parents ( N _b) in the largest sample of hatchery fish confirmed that effective population size was low ( N _b= 14·9 for multiallelic loci and N _b= 10·6 if alleles were pooled into two composite alleles). The results show that low effective population sizes leading to loss of variation and possibly inbreeding depression should be a matter of serious concern in aquaculture production of catla.     

Genetic analysis of fish genomes and populations: allozyme variation within and among Atlantic salmon from Downeast rivers of Maine

Analysis of genetic variation within and among samples of naturally produced Atlantic salmon ( n  = 372) from 7 Maine (U.S.A.) and one Canadian river were conducted based on 54 allozyme loci. Eight of the 54 loci proved polymorphic, and estimated heterozygosities ( H _S) based on all loci ranged from 0·012 to 0·026 (mean = 0·021, s . e . = 0·002). Only one of 56 tests revealed genotypic proportions that deviated significantly from Hardy–Weinberg expectations. Genetic distances ( D ) between samples ranged from 0·002 to 0·022. No obvious association existed between genetic and geographic distances. Cluster analysis of genetic distances revealed the Dennys River sample as the most differentiated when all samples were included in the analysis, though bootstrap support of the cluster analysis was generally weak. G ‐tests revealed significant differences in allele frequencies among samples at five of the polymorphic loci, and the G ‐value summed over all loci also indicated significant differences among samples. F _ST values indicated that 3·4% of the total genetic diversity was due to variability among samples, while 96·6% was due to variability within samples. These results indicate that the Atlantic salmon analyzed in this study had levels of genetic variability and differentiation among samples comparable to native populations from other areas collected across a similar geographic range.     

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.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Allozymic variation and relationships of the endangered cyprinodontid genus Valencia and its implications for conservation

Allozyme analysis was used to determine patterns of genetic variation and relationships within the genus Valencia , a group with two allopatric species, V. letourneuxi inhabiting the Balkan Peninsula, and V. hispanica which is endemic to the Iberian Peninsula. Products of 25 gene loci were analysed with AK *, FH *, LDH -1* and LDH -2* monomorphic when Fundulus heteroclitus macrolepidotus was not taken into account. Nine loci were diagnostic for both species of the genus Valencia . Levels of genetic variation of V. hispanica ( P =0·416, H ^e=0·118) and V. letourneuxi ( P =0·160, H _e=0·040) were higher than values reported previously. High population subdivision ( F _st=0·321) among V. hispanica populations indicates a clear interruption of genetic exchange among populations. High genetic variation and differentiation exhibited by populations of V. hispanica suggest that a recovery program should be carried out with natural stocks from the same localities.     

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.     

Temporal effective size estimates of a managed walleye Sander vitreus population and implications for genetic-based management

The goal of this research was to use the long-term fishery data set and DNA from archived scales of walleye Sander vitreus in Escanaba Lake, WI, U.S.A., to improve the understanding of the underlying mechanism(s) influencing genetic diversity in naturally recruiting populations. The introduced population of S. vitreus in Escanaba Lake has a low mean effective population size ( N _E) between 124·6 and 185·5 despite a mean census size ( N _C) of 4659 ( N _E/ N _C c. 0·04), suggesting an accelerated rate of genetic drift between 1952 and 2002. These values are smaller than the median N _E range of several studies suggesting typical N _E/ N _C ratios of 0·11–0·16 in a wide range of taxa. N _E increased steadily during the past two sampled decades (1992 and 2002) and was consistent with a lowering of the variance in S. vitreus reproductive success, possibly linked to a large, sustained exploitation (mean 28%) rate. Variance in reproductive success is one of the most important factors influencing N _E in species, like S. vitreus , which have a potential for large fecundities and large juvenile mortalities (type III survivorship). The N _B estimates across six sequential cohorts (age classes of S. vitreus , assayed from 1994 to 1999) was consistent with estimates of N _E reported for 1992–2002. These results, coupled with in-depth census and exploitation data, show that the genetic characteristics of Escanaba Lake S. vitreus have changed substantially and that management activities, such as supplemental stocking and harvest practices, have profoundly influenced the genetic dynamics of S. vitreus in this lake.     

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.     

Factors affecting the distribution of genetic variability in the guppy, Poecilia reticulata

Previous studies on Trinidadian guppies have shown an apparent association of genetic variability, expressed as mean heterozygosity (Ho) at allozyme loci, with river order: fish from lowland sites exhibit higher mean Ho than their upland counterparts. Detailed genetic and behavioural observations were undertaken in the present study by repeated sampling throughout two river courses to examine population heterogeneity. Results show that the predation regime mediated mating behaviour of upland and lowland populations in the wild is not as divergent as previous laboratory observations have suggested. Furthermore, genetic variation at the eight loci variable in Trinidadian populations (p₁₀₀=0·13–O·88, H₀=O·03–0·23) appears to bear little association with predation regime, but is related to sub-structuring of populations by habitat features and geographic isolation. We therefore conclude that the observed patterns of genetic variability arise primarily as a consequence of stochastic as opposed to deterministic factors.     

Genetic and morphological variation in a common European cyprinid, Leuciscus cephalus within and across Central European drainages

Genetic variability and differentiation of the common European cyprinid Leusciscus cephalus was investigated within and across the drainages of the rivers Rhine, Elbe and Danube using starch gel electrophoresis and morphological characters. Eleven of 24 presumptive enzyme loci showed some allozymic variation. Compared to other freshwater fish, genetic variability within all sites was high ( H _e=0·074–0·113). This was explained by the wide geographical and ecological range of the chub. In contrast, genetic and morphological differentiation between sites and drainages was low (F_ST=0·15), which may be due to the high dispersal ability of the chub and translocations of spawn by waterfowl. The phylogeographic structure of Leuciscus cephalus in north-eastern Bavaria showed low congruence with Cottus gobio populations sampled in the same area.     

Low genetic variability of the koala Phascolarctos cinereus in south-eastern Australia following a severe population bottleneck

Genotyping of koalas at CA-repeat microsatellite loci has revealed significant differences in the levels of allelic diversity ( A ) and expected heterozygosity ( H¯ _E) between populations from north-eastern and south-eastern Australia. In the 10 populations studied, allelic diversity ranged from 8.0 in the Nowendoc population to 1.7 in the Kangaroo Is. population, and values of H¯ _E ranged from 0.831 in the Nowendoc population to 0.331 in the Kangaroo Is. population. Data from pooled populations revealed koalas from the north-eastern region had significantly higher levels of allelic diversity ( A = 11.5 ± 1.4) than those from south-eastern Australia ( A = 5.3 ± 1.0). Furthermore significantly higher heterozygosity levels were found in the north-eastern ( H¯ _E= 0.851) vs. the south-eastern ( H¯ _E= 0.436) regions of Australia. Following a near-extinction bottleneck in the 1920s, mainland Victorian and Kangaroo Is. koalas have been involved in an extensive program of relocations. The source populations of the relocated animals were islands in Westernport Bay, which were founded by very few individuals in the late 1800s and early 1900s. The significantly lower levels of variation between south-eastern Australian populations suggests that human intervention has had a severe effect on levels of genetic diversity in this region, and this may have long-term genetic consequences.     

Genetic differentiation in pikeperch (Sander lucioperca): the relative importance of gene flow, drift and common history

To investigate the relative importance of homogenizing factors, such as gene flow, and diversifying factors, such as drift, genetic variation in pikeperch ( Sander lucioperca ) in two Fennoscandian regions (North and South) was analysed with microsatellites. Allelic richness and the degree of differentiation were significantly higher in the North ( F _ST= 0·20) than in the South ( F _ST= 0·064). In northern areas, assignments of genotypes were almost exclusively to the population of origin, but in southern areas, the proportion of correct assignments was significantly lower. Most samples exhibited significant heterozygote deficits, and the level of relatedness was higher than expected from randomness. These combined results suggest that there has been more gene flow between populations in southern areas than in northern areas, where the importance of genetic drift has been greater. Effective population sizes were small ( c. 100) and did not differ between areas. The effect of a common history appears minor, and thus processes such as genetic drift and gene flow have been more influential in shaping the patterns of genetic diversity in this species.     

Conservation Genetics of an Endangered Herb, Hanabusaya asiatica (Campanulaceae)

Abstract: Hanabusaya asiatica (Nakai) Nakai (Campanulaceae), a bee- pollinated, perennial herb, is restricted to the mountainous regions of the eastern-central Korean peninsula. Allozyme analyses for 348 individuals assessed the levels of genetic diversity for five populations. Spatial autocorrelation statistics were also used to examine the spatial distribution of allozyme polymorphisms. The species maintains high levels of allozyme diversity ( H _eS = 0.217) and it exhibits low allozyme differentiation among populations ( G _ST = 0.132) compared with other endemics (mean H _e = 0.096, G _ST = 0.248). There is an apparent pattern of isolation by distance among populations. These results suggest that H. asiatica is at a genetic equilibrium. A considerable deficit in numbers of heterozygotes suggests mating among relatives in populations. At least three populations of H. asiatica should be sampled or conserved to capture or maintain > 99 % of the genetic diversity in the species as a whole. Within local populations, individuals are distributed in a structured, isolation by distance, manner. Approximate genetic patch width in the populations of H. asiatica examined is 5 - 8 m. For conservation purposes, it is suggested that, in general, the sampling of H. asiatica should be conducted at intervals in order to efficiently sample the genetic diversity across an entire population.