Incongruent estimates of population differentiation among brook charr, Salvelinus fontinalis, from Cape Race, Newfoundland, Canada, based upon allozyme and mitochondrial DNA variation

We determined the amount and temporal stability of genetic differentiation among brook cham sampled from five rivers on Cape Race, Newfoundland, with an electrophoretic analysis of 42 protein coding loci. Fish from four of these rivers were analysed for restriction fragment length polymorphisms in mitochondrial DNA (mtDNA). A single mtDNA clone was observed in all rivers sampled, except one, where 47% offish were from a different and relatively divergent clone (0.31 % sequence divergence). In contrast, Cape Race brook charr show large amounts of genetic differentiation at six enzyme coding loci; Nei's genetic distance ranged between 0,020 and 0.048. This differentiation is relatively stable as no significant differences in allele frequencies were detected between fish sampled from two rivers over two consecutive years. The most divergent population based on protein polymorphism is not that with two mtDNA clonal lineages. In contrast to the commonly held view, mtDNA analyses do not necessarily provide greater resolution of population structure than allozyme analyses.     

Loss of genetic integrity correlates with stocking intensity in brook charr (Salvelinus fontinalis)

Supportive breeding and stocking performed with non‐native or domesticated fish to support sport fishery industry is a common practice throughout the world. Such practices are likely to modify the genetic integrity of natural populations depending on the extent of genetic differences between domesticated and wild fish and on the intensity of stocking. The purpose of this study is to assess the effects of variable stocking intensities on patterns of genetic diversity and population differentiation among nearly 2000 brook charr (Salvelinus fontinalis) from 24 lakes located in two wildlife reserves in Québec, Canada. Our results indicated that the level of genetic diversity was increased in more intensively stocked lakes, mainly due to the introduction of new alleles of domestic origin. As a consequence, the population genetic structure was strongly homogenized by intense stocking. Heavily stocked lakes presented higher admixture levels and lower levels of among lakes genetic differentiation than moderately and un‐stocked lakes. Moreover, the number of stocking events explained the observed pattern of population genetic structure as much as hydrographical connections among lakes in each reserve. We discuss the implications for the conservation of exploited fish populations and the management of stocking practices.     

Quantitative genetic analysis of the physiological stress response in three strains of brook charr Salvelinus fontinalis and their hybrids

Three strains [domestic (D), Laval (L) and Rupert (R)] of brook charr Salvelinus fontinalis and their reciprocal hybrids were submitted to transport stress to measure stress resistance. Primary (cortisol) and secondary (glucose, osmolality and haematocrit) stress responses were measured for each cross. Significant heritabilities were observed for both levels of stress response, with mean ± S.E. heritability (h(2)) = 0.60 ± 0.20 for plasma cortisol and 0.61 ± 0.20 for plasma glucose. There were strain differences whereby the R strain was the least sensitive to stress at the primary and secondary levels. No heterosis was detected, and only one case of outbreeding depression was present. The outbreeding depression was observed in the D(♀) R(♂) hybrid, which had a 27% increase of plasma glucose compared to parental strains. The D(♀) R(♂) and R(♀) L(♂) hybrids had more pronounced variations (increase or decrease) in plasma osmolality than their respective parental strains, but these variations were difficult to relate definitively with the potential secondary stress response. These results indicate a strong potential for genetic improvement in the stress response to transport with the use of purebred crosses while hybridization has little value in this regard.     

Genetic discrimination of wild and hatchery populations of brook charr, Salvelinus fontinalis (Mitchill), in Ontario using mitochondrial DNA analysis

Two brood stocks of brook charr, Sulvelinus fontinalis , are currently maintained by the Ontario Ministry of Natural Resources. The Nipigon brood stock originated from Lake Nipigon, in north-central Ontario, while the Hills Lake stock is believed to have been produced by hybridizing several strains (including charr from a Pennsylvania hatchery as well as charr from Ontario) in the past. The mitochondrial DNA (mtDNA) variability of these brood stocks was characterized using 51 hexanucleotide restriction enzymes. Eleven restriction enzymes, Acc I, Am I, Bun I, Bun II, Hind III, Nco I, Nde I, Nhe I, Nsi I, Pst I and Sph I, were polymorphic between the two brood stocks. Eight hatchery mtDNA haplotypes were detected showing a maximum of 0.41 % sequence divergence. Seven haplotypes are present in the Hills Lake strain, and two in the Nipigon strain. These mtDNA haplotypes are useful markers to determine the degree of reproductive success between planted and native fish. In one comparison in southern Ontario, less than 20% of the wild fish sampled from the head water regions of a small drainage entering Lake Erie could have resulted from random introgression with hatchery fish planted further downstream, because most of these fish possessed a unique Acc 1 cut site. This is one of the few examples in stock analysis studies where such a high degree of genetic discrimination is evident between hatchery and native fish.     

Microsatellite variation and genetic structure of brook trout (Salvelinus fontinalis) populations in Labrador and neighboring Atlantic Canada: evidence for ongoing gene flow and dual routes of post-Wisconsinan colonization

In conservation genetics and management, it is important to understand the contribution of historical and contemporary processes to geographic patterns of genetic structure in order to characterize and preserve diversity. As part of a 10-year monitoring program by the Government of Newfoundland and Labrador, Canada, we measured the population genetic structure of the world's most northern native populations of brook trout (Salvelinus fontinalis) in Labrador to gather baseline data to facilitate monitoring of future impacts of the recently opened Trans-Labrador Highway. Six-locus microsatellite profiles were obtained from 1130 fish representing 32 populations from six local regions. Genetic diversity in brook trout populations in Labrador (average H(E)= 0.620) is within the spectrum of variability found in other brook trout across their northeastern range, with limited ongoing gene flow occurring between populations (average pairwise F(ST)= 0.139). Evidence for some contribution of historical processes shaping genetic structure was inferred from an isolation-by-distance analysis, while dual routes of post-Wisconsinan recolonization were indicated by STRUCTURE analysis: K= 2 was the most likely number of genetic groups, revealing a separation between northern and west-central Labrador from all remaining populations. Our results represent the first data from the nuclear genome of brook trout in Labrador and emphasize the usefulness of microsatellite data for revealing the extent to which genetic structure is shaped by both historical and contemporary processes.     

Hydrography and population genetic structure in brook charr (Salvelinus fontinalis, Mitchill) from eastern Canada

Despite the abundance of studies of genetic diversity in freshwater fishes, few have specifically addressed the role of habitat structure in partitioning genetic variance within and among populations. In this study, we analysed the variability of six microsatellite loci among 24 brook charr population samples in order to correlate hydrographic structure with genetic organization. These populations originated from three Canadian National parks (Kouchibouguac, Fundy and Forillon) that showed distinct hydrographic structure. Considering the general characteristics of these habitats, we formulated specific hypotheses in regard to genetic structure, which were principally based on the potential for gene flow and population size associated with each habitat. The hierarchical analysis of molecular variance and the genetic distances computed among populations revealed that habitat structure analyses constitute an important, but insufficient, predictor of genetic structure. We discuss the importance of habitat complexity on genetic structure in the context of management and conservation.     

Pyloric caecal morphology of brook charr,Salvelinus fontinalis,in relation to diet

Synopsis Brook charr, Salvelinus fontinalis, shifts its diet from zoobenthos to pelagic prey when living sympatrically with white sucker, Catostomus commersoni, in lakes of the Laurentian Shield. We tested the hypothesis that this diet difference would have a significant impact on their pyloric caecal morphology in 5 lakes containing allopatric brook charr populations and 6 other lakes containing both brook charr and white sucker. We observed that the mean length of the most posterior caecum of charr was significantly greater in sympatry than in allopatry (X ± 1 SD: 9.91 ± 1.12 mm versus 8.44 ± 0.67 mm). This is equivalent to an increase of 18% of total pyloric caecal mass (dry weight) in sympatric brook charr. These results indicate that this response to differences in diet, well known in birds, also occurs in fish.     

Dynamics of introgressive hybridization assessed by SNP population genomics of coding genes in stocked brook charr (Salvelinus fontinalis)

Salmonid fishes rank among species being most severely affected by introgressive hybridization as a result of a long tradition of stocking with hatchery-reared conspecifics. The objectives of this study were (i) to evaluate the genetic consequences of stocking and resulting introgression rates on the genetic integrity of natural populations of brook charr, (ii) to identify genomic regions potentially associated with adaptation to natural and artificial rearing environments, and (iii) to test the null hypothesis that introgression from domesticated brook charr into wild populations is homogeneous among loci. A total of 336 individuals were sampled from nine lakes, which were stocked at different intensities with domestic fish. Individuals were genotyped at 280 SNPs located in transcribed regions and developed by means of next-generation sequencing. As previously reported with microsatellites, we observed a positive relationship between stocking intensity and genetic diversity among stocking groups, and a decrease in population differentiation. Individual admixture proportions also increased with stocking intensity. Moreover, genomic cline analysis revealed 27 SNPs, seven of which were also identified as outliers in a genome scan, which showed an introgression rate either more restricted or enhanced relative to neutral expectations. This indicated that selection, mainly for growth-related biological processes, has favored or hampered the introgression of genomic blocks into the introgressed wild populations. Overall, this study highlights the usefulness of investigating the impact of stocking on the dynamics of introgression of potentially adaptive genetic variation to better understand the consequences of such practice on the genomic integrity of wild populations.     

Phylogeography of the endangered Cathaya argyrophylla (Pinaceae) inferred from sequence variation of mitochondrial and nuclear DNA

Cathaya argyrophylla is an endangered conifer restricted to subtropical mountains of China. To study phylogeographical pattern and demographic history of C. argyrophylla, species-wide genetic variation was investigated using sequences of maternally inherited mtDNA and biparentally inherited nuclear DNA. Of 15 populations sampled from all four distinct regions, only three mitotypes were detected at two loci, without single region having a mixed composition (G(ST) = 1). Average nucleotide diversity (theta(ws) = 0.0024; pi(s) = 0.0029) across eight nuclear loci is significantly lower than those found for other conifers (theta(ws) = 0.003 approximately 0.015; pi(s) = 0.002 approximately 0.012) based on estimates of multiple loci. Because of its highest diversity among the eight nuclear loci and evolving neutrally, one locus (2009) was further used for phylogeographical studies and eight haplotypes resulting from 12 polymorphic sites were obtained from 98 individuals. All the four distinct regions had at least four haplotypes, with the Dalou region (DL) having the highest diversity and the Bamian region (BM) the lowest, paralleling the result of the eight nuclear loci. An AMOVA revealed significant proportion of diversity attributable to differences among regions (13.4%) and among populations within regions (8.9%). F(ST) analysis also indicated significantly high differentiation among populations (F(ST) = 0.22) and between regions (F(ST) = 0.12-0.38). Non-overlapping distribution of mitotypes and high genetic differentiation among the distinct geographical groups suggest the existence of at least four separate glacial refugia. Based on network and mismatch distribution analyses, we do not find evidence of long distance dispersal and population expansion in C. argyrophylla. Ex situ conservation and artificial crossing are recommended for the management of this endangered species.     

Low genetic variability in lake populations of brook trout (Salvelinus fontinalis): A consequence of exploitation?

Previous studies have found lower levels ofgenetic variation in lake than streampopulations of brook trout (Salvelinusfontinalis). We test the generality of thisobservation by examining whether brook troutgenetic variation at 10 allozyme loci differedwithin and among 9 pairs of lake and adjacentstream populations. With one exception, wefound that lake populations had lowerheterozygosity than their adjacent streampopulations. Although the lakes in this studyare small and some have had documented fishmortality events, no association was foundbetween lake size characteristics and thedegree of difference in heterozygosity betweenlakes and their adjacent stream populations. There were, however, negative associationsbetween metrics of fishing mortality and thedifference in heterozygosity between lakes andtheir adjacent stream populations. Thegreater the estimated fishing pressure onlake-dwelling trout, the greater the reductionin heterozygosity in those populationsrelative to their adjacent stream populations. We interpret our findings to suggest thatintensive fishing pressure can significantlyreduce genetic variation. Managers shouldtherefore prevent human-induced mortality atany indication of a large natural mortalityevent to allow populations to increase in sizeas rapidly as possible following a decline.     

Mitochondrial DNA diversity among four sympatric morphs of Arctic charr, Salvelinus alpinus L., from Thingvallavatn, Iceland

Restriction fragment length polymorphisms in mitochondrial DNA (mtDNA) were examined in progeny of four sympatric morphs of Arctic charr ( Salvelinus alpinus L.) from Thingvallavatn, Iceland. The mtDNA analysis with 46 hexanucleotide restriction enzymes indicates that the Thingvallavatn morphs are very closely related. Sequence divergence was less than 0.2% between any of the five clones detected. Although not significant, the topology of the mtDNA UPGMA dendogram was similar to that from a previous allozyme survey. Icelandic Arctic charr show greater affinity to charr from the British Isles than those from North America supporting their current taxonomic distinction.     

Water temperature, light intensity and zooplankton density and the feeding activity of juvenile brook charr (Salvelinus fontinalis)

SUMMARY 1. The objective of this study was to evaluate the effects of zooplankton biomass (as a measure of density), fish biomass, light intensity and water temperature on the attack rate and swimming characteristics (i.e. swimming speed and angle of turn) of juvenile (1+) brook charr (Salvelinus fontinalis) in field enclosures. We used a portable underwater camera system in a series of pelagic enclosures to quantify the feeding behaviour of brook charr over a gradient of natural conditions. 2. In simple linear or non‐linear regression models we found (i) that attack rate and angle of turn were positively related to water temperature, (ii) that attack rate and swimming speed were positively related to zooplankton biomass and light intensity and (iii) that attack rate was positively related to swimming speed. In multiple regression models, fish biomass, light intensity and variance of the angle of turn accounted for 87% of the variation in attack rate. Light intensity and water temperature accounted for 86% of the variation in swimming speed. Fish gut fullness and attack rate accounted for 83% of the variation in the variance of the angle of turn executed by fish. 3. The increase in the number of attacks as zooplankton biomass increases conforms to the general positive functional response observed in other fish species. Our results also support the hypothesis that swimming speed increases with prey biomass. We did not observe a plateau in attack rate as zooplankton biomass increased. As our experiments were performed under natural biotic and abiotic conditions, factors other than zooplankton biomass might affect or limit this response, such as water temperature and light intensity. 4. Because zooplankton biomass was correlated with water temperature and light intensity, it was not possible to evaluate the independent contribution of these factors on the attack rate and swimming characteristics (swimming speeds and angle of turn) of brook charr. However, this study highlighted the impact of these factors on the feeding behaviour of juvenile brook charr when feeding in the pelagic habitat under natural conditions, and their importance in future models of optimal foraging and fish habitat quality.     

Intraspecific variation and phylogeographic patterns of Fagus crenata (Fagaceae) mitochondrial DNA

Mitochondrial (mt) DNA variation in Japanese beech, Fagus crenata (Fagaceae), was studied in 17 populations distributed throughout the species' range. Total genomic DNA of samples from single trees representing each of 12 populations were digested with 18 restriction enzymes and hybridized with three probes containing coxI, coxIII, and atpA gene sequences. Thirty-four of the 54 enzyme/probe combinations showed polymorphisms and all the individuals were subsequently analyzed with six combinations of three probes and two enzymes. Restriction fragment length polymorphisms were evident around all three genes, allowing the identification of eight distinct haplotypes. Haplotype diversity within the populations was found to be very low (HS = 0.031), but population differentiation to be much higher (GST = 0.963). The mtDNA variation was strikingly different from allozyme variation (HS = 0.209; GST = 0.039). Gene flow for maternally inherited mtDNA should be restricted to seed dispersal while nuclear gene flow occurs by both seed and pollen dispersal. Therefore, the difference in the variation between mtDNA and allozymes may be largely a result of the much higher rate of gene flow associated with pollen dispersal than with seed dispersal. The mtDNA variation displayed strong geographic structure, which may reflect the species' distribution in the last glacial maximum and subsequent colonization, and probably also reflects intraspecific phylogeography of the species.     

Heritability of life-history tactics and genetic correlation with body size in a natural population of brook charr (Salvelinus fontinalis)

A common dimorphism in life-history tactic in salmonids is the presence of an anadromous pathway involving a migration to sea followed by a freshwater reproduction, along with an entirely freshwater resident tactic. Although common, the genetic and environmental influence on the adoption of a particular life-history tactic has rarely been studied under natural conditions. Here, we used sibship-reconstruction based on microsatellite data and an 'animal model' approach to estimate the additive genetic basis of the life-history tactic adopted (anadromy vs. residency) in a natural population of brook charr, Salvelinus fontinalis. We also assess its genetic correlation with phenotypic correlated traits, body size and body shape. Significant heritability was observed for life-history tactic (varying from 0.52 to 0.56 depending on the pedigree scenario adopted) as well as for body size (from 0.44 to 0.50). There was also a significant genetic correlation between these two traits, whereby anadromous fish were genetically associated with bigger size at age 1 (r(G) = -0.52 and -0.61). Our findings thus indicate that life-history tactics in this population have the potential to evolve in response to selection acting on the tactic itself or indirectly via selection on body size. This study is one of the very few to have successfully used sibship-reconstruction to estimate quantitative genetic parameters under wild conditions.     

Ecological mechanisms favouring behavioural diversification in the absence of morphological diversification: a theoretical examination using brook charr (Salvelinus fontinalis)

1. Behavioural diversification is thought to be an important initial step in the origin of resource polymorphisms. We developed a model for young brook charr (Salvelinus fontinalis Mitchill) to examine four mechanisms that could generate a U-shaped relationship between growth rate (fitness) and the proportion of time spent moving that would favour alternative foraging tactics in the absence of obvious differences in body size and shape. 2. Recently emerged brook charr of similar size and shape inhabit still-water pools along the sides of streams. Some individuals tend to sit and wait for crustacean prey at the pool substrate near the bank, while others tend to search actively for insect prey at the pool surface away from the bank. 3. The ecological mechanisms modelled were (i) the relationship between the rate of prey capture and the proportion of time spent moving is curvilinear, such that net rate of energy gain is maximized at two different levels of activity; (ii) switching between foraging locations and, hence, tactics involves lost opportunity and travel costs; (iii) switching between prey types and, hence, tactics involves a learning cost; and (iv) foraging success is status-dependent with individuals switching between tactics having a lower status than those specializing at a tactic. 4. Singly, no mechanism predicted the U-shaped relationship between growth rate and the proportion of time spent moving. Together, a U-shaped relationship was obtained, indicating that the behavioural diversification and diversifying selection observed in the field may be a consequence of multiple, subtle mechanisms.     

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.     

A comparative mitogenomic analysis of the potential adaptive value of Arctic charr mtDNA introgression in brook charr populations (Salvelinus fontinalis Mitchill)

Wild brook charr populations (Salvelinus fontinalis) completely introgressed with the mitochondrial genome (mtDNA) of arctic charr (Salvelinus alpinus) are found in several lakes of northeastern Québec, Canada. Mitochondrial respiratory enzymes of these populations are thus encoded by their own nuclear DNA and by arctic charr mtDNA. In the present study we performed a comparative sequence analysis of the whole mitochondrial genome of both brook and arctic charr to identify the distribution of mutational differences across these two genomes. This analysis revealed 47 amino acid replacements, 45 of which were confined to subunits of the NADH dehydrogenase complex (Complex I), one in the cox3 gene (Complex IV), and one in the atp8 gene (Complex V). A cladistic approach performed with brook charr, arctic charr, and two other salmonid fishes (rainbow trout [Oncorhynchus mykiss] and Atlantic salmon [Salmo salar]) revealed that only five amino acid replacements were specific to the charr comparison and not shared with the other two salmonids. In addition, five amino acid substitutions localized in the nad2 and nad5 genes denoted negative scores according to the functional properties of amino acids and, therefore, could possibly have an impact on the structure and functional properties of these mitochondrial peptides. The comparison of both brook and arctic charr mtDNA with that of rainbow trout also revealed a relatively constant mutation rate for each specific gene among species, whereas the rate was quite different among genes. This pattern held for both synonymous and nonsynonymous nucleotide positions. These results, therefore, support the hypothesis of selective constraints acting on synonymous codon usage.