Differences between nuclear and mitochondrial introgressions of brown trout populations from a restocked main river and its unrestocked tributary

The introgression of domestic genes in wild trout populations as a result of intensive restocking has been observed to be low in numerous rivers of the Mediterranean watershed. With the aim of examining the process of introgression and improving our understanding the reasons for its possible limitations, we analysed the percentage of domestic genes for three protein-coding loci and for mitochondrial haplotypes. Two samples were analysed: a restocked river and a tributary which is no longer stocked. In order to follow the introgression over time, trout were separated into two age groups. The mitochondrial locus showed the same tendencies as those obtained with the nuclear loci, but with generally higher introgression values. A differential introgression among protein-coding loci was also observed. From these results, we discuss the possibility of selective factors acting against the domestic genes in river. The hypothesis of a difference in reproductive success between domestic males and females is also presented to explain the higher introgression of mtDNA.     

Molecular Evidence for Introgression and Loss of Genetic Variability in Salmo (trutta) macrostigma as a Result of Massive Restocking of Apennine Populations (Northern and Central Italy)

We assessed structural gene variation (allozymes and mtDNA) of brown trout to evaluate the genetic variability of Apennine stream populations (Northern and Central Italy) and the possibility of introgression by alien genomes after massive restocking with hatchery strains (Atlantic stocks). Genetic variability within and between Apennine populations was extremely low in our samples. Only two allozyme loci were polymorphic and mean hetero-zygosity was also reduced compared to other brown trout populations. Allelic frequencies determined for both loci were similar to the ones detected in the corresponding hatchery spawners. The reduction or total absence of the Mediterranean nuclear (LDH-5) and mitochondrial (16S rDNA) diagnostic markers suggests the domestic origin of most populations, and the introgression effects carried out by non-native genomes. From a taxonomic point of view, a clear differentiation emerges among basins placed on opposite sides of the Apennine chain (Tyrrhenian and Adriatic regions). In particular, the presence of Mediterranean genotypes and haplotypes characterizing Salmo (trutta) macrostigma is sporadic along the eastern Apennine side, adding additional doubts on the original presence and wide distribution of this salmonid along the Adriatic side of the mountain chain. In spite of conservation programs devoted to preservation of local genetic characteristics of S. t. macrostigma, massive restocking practices with hatchery strains obtained by a few spawners is the major cause of significant `founder effect' and `inbreeding depression' even in Apennine regions.     

Temporal variations of genetic introgression in stocked brown trout populations

Allozymes were used to analyse the genetic impact of hatchery brown trout Salmo trutt a, morpha fario stocked in wild Mediterranean populations to gain a better understanding of the mechanisms of introgression (regulation, elimination or homogenization). Analysis of the genetic structure of populations from the same river drainage basin but subjected to different incidences was performed in space and time (data on two generations and 2 years of sampling). Introgression is associated with high deficits of heterozygotes and linkage disequilibria. Genetic divergence according to age group was observed. These results may indicate selective forces acting against domesticated genes and limiting hybridization between the two forms.     

Comparative analysis of introgression at three marker classes: a case study in a stocked population of brown trout

Comparative analysis of protein loci, microsatellite and mtDNA markers revealed generally comparable estimates for introgression and apparent admixture rates in stocked brown trout populations at two sites in the River Doubs (Rhône dainage, Switzerland), which are 10 km apart and which belong to the same management unit. At one site, a significant deviation between mtDNA and nuclear markers could be explained by stocking of F1 hybrids originating from crosses between hatchery females and males from the local population. Substantial differences between diagnostic protein loci and protein loci having non-fixed private alleles indicated that caution must be exercised when using genetic markers not strictly diagnostic for the distinction of the populations under investigation. Congruent estimates of introgression and apparent admixture rates between diagnostic protein loci and presumed diagnostic microsatellite loci suggest that the latter can be regarded as reliable genetic markers for the estimation of introgression in Mediterranean brown trout populations stocked with trout of Atlantic origin. Significant differences in introgression and apparent admixture rates between the two sites and between age-classes of one study site were observed. Introgression is suggested to depend on environmental factors. Significantly lower introgression rates in age-class 2+ years as compared to juvenile trout might further indicate that introduced Atlantic brown trout and hybrids decrease in proportion between age-classes 1+ and 2+ years.     

THE SELECTION BARRIER BETWEEN POPULATIONS SUBJECT TO STABILIZING SELECTION

The introgression of genes carried by a small group of immigrants is studied. The recipient and the donor populations differ at several autosomal loci subject to weak selection, and two allelic forms of each gene are considered. Fitness variation is determined by additive allelic effects, by dominance effects, and by two-locus additive-by-additive epistatic interaction of the effects of the alleles. The fate of the group of immigrants is quantified by the selection barrier that describes the cumulative mean fitness of the hybrids and hybrid descendants relative to the fitness of the resident population. The monomorphic and the polymorphic loci of the recipient population contribute differently to the selection barrier. If the genetic difference between recipient and donor population is small, then the contribution of the monomorphic loci is dominated by a positive term dependent on the difference in gene frequencies. The contribution of the polymorphic loci depends only on the difference of the leading order in the pairwise linkage disequilibria between the two populations. This contribution may be positive or negative; and, thus, polymorphic loci may either contribute to the barrier or inflate the introgression.     

Genetic structure of five susceptibility gene regions for coronary artery disease: disequilibria within and among regions

We analyzed two-locus disequilibria for 16 polymorphic loci of seven susceptibility genes for coronary artery disease located in five chromosomal regions distributed across four chromosomes. Included were the genes coding for apolipoprotein B (ApoB, chromosome 2, four marker loci), lipoprotein lipase (LPL, chromosome 8, three marker loci), apolipoproteins AI, CIII, AIV (ApoAI–CIII–AIV, chromosome 11, three marker loci), apolipoprotein E (ApoE, chromosome 19, two marker loci), and the low density lipoprotein receptor (LDLR, chromosome 19, four marker loci). Our sample included 540 unrelated individuals from the Rochester, Minn. population. There were no statistically significant deviations of single-locus genotypes from Hardy-Weinberg equilibrium. The strongest associations within genes were for composite diallelic disequilibria; 17/19 were significant (13 at Pr <0.001, 1 at Pr <0.01, 3 at Pr <0.05). These observations suggest marker alleles within genes have a shared evolutionary history reflected by disequilibria that have not been dissipated by recombination. Disequilibrium was not generally concordant with the physical orderings of markers. Only two significant higher-order disequilibria were observed although 12 triallelic disequilibria were at maximum possible values. We observed 19 statistically significant disequilibria (Pr <0.05; 4 composite diallelic, 13 triallelic, and 2 quadriallelic) between 101 pairs of marker loci, where each locus in a pair was from a different unlinked region. These unexpected results are most likely explained by recent historical factors, including worldwide population expansion and amalgamation with continuous admixture, that influence the genetic structure (organization of alleles and non-alleles into genotypes) of a population. We conclude that disequilibria between loci from unlinked regions may be more extensive than is commonly assumed. Our findings also suggest that it is, on average, at least 15 times more likely to not detect significant disequilibrium among unlinked loci when it is really present than to make a false positive inference. Disequilibria between functional loci within or between regions will impact estimates of genetic variance associated with particular functional mutations within a susceptibility gene region. Received: 15 January 1998 / Accepted: 24 June 1998     

Assessment of restocking as a strategy for rehabilitating a native population of brown trout Salmo trutta L. in a fast-flowing mountain stream in the northern French Alps

An attempt was made to extend the area of distribution of a native population of brown trout Salmo trutta belonging to a Mediterranean lineage (ML), which has maintained itself in the Dranse d'Abondance, a fast-flowing alpine stream (Haute-Savoie, France), despite several decades of intensive restocking with brown trout derived from the Atlantic lineage (AL). This was done by releasing an ML component into the predominantly AL population still present on the Ugine, the main tributary of the Dranse d'Abondance. This strategy of rehabilitation restocking was tested using fluoro-marked juveniles produced from a captive breeding stock derived from the wild Dranse d'Abondance ML stock. Samples of 0+ year fish were collected over the period 1995–2003 in order to assess the impact of the restocking. Percentages of fluoro-marked otoliths revealed significant contributions of ML restocking in the 0+ year autumnal standing population, with levels ranging from 34·3 to 61·4%. The change in the genetic characteristics of the 0+ year population produced by natural recruitment was monitored by analysing the unmarked subjects. Frequencies observed at two microsatellite loci revealed a considerable rise (from 0 to 60%) in the level of Mediterranean alleles in the natural 0+ year population since the introduction of restocking using ML individuals.     

Using genotype simulations and Bayesian analyses to identify individuals of hybrid origin in Australian bass: lessons for fisheries management

Simulated genotypes, Bayesian analyses and molecular genetic data were used to detect individuals of hybrid origin and hybrid introgression between the Australian bass ( Macquaria novemaculeata ), a species extensively stocked in Australia, and estuary perch ( Macquaria colonorum ). Based on this analytical framework, 93% of the hybrids up to three generations later could be distinguished from the Australian bass. Individuals of hybrid origin were identified in all three rivers sampled. In addition, this study verified the fertility of hybrids between Australian bass and estuary perch as determined through genomic introgression. This study exemplifies an analytical procedure that has implications for identifying suitable individuals for use in breeding and restocking programmes for other species.     

Population genetic response to microsite ecological stress in wild barley, Hordeum spontaneum

Genetic diversity was studied in six subpopulations (a total of 60 individuals) of wild barley, Hordeum spontaneum , the progenitor of cultivated barley, sampled from six stations located along a transect of 300 m across the two opposing slopes of 'Evolution Canyon', a Mediterranean microsite at Lower Nahal Oren, Mt Carmel. The two opposing slopes are separated by between 100 and 400 m and designated SFS (South-Facing Slope) and NFS (North-Facing Slope) with each having three equidistant test stations. The SFS, which receives up to 300% more solar radiation, is drier, ecologically more heterogeneous, fluctuating, and more stressful than the NFS. Analysis of 12 RAPD primers, representing a total of 51 putative loci, revealed a significant inter- and intraslope variation in RAPD band polymorphism. A significantly higher proportion of polymorphic RAPD loci was found amongst the subpopulations on the SFS (mean P = 0.909) than on the NFS (mean P = 0.682), on the basis of the presence and absence of 22 strong bands. Polymorphism generally increased upwards from the bottom to the top of the SFS (0.636, 0.773, 0.955) and NFS (0.409, 0.500, 0.545), respectively. Gametic phase disequilibria estimates, D, revealed SFS and NFS unique predominant combinations which sharply differentiated the two slopes and indicated that there is differential interslope selection favouring slope-specific multilocus combinations of alleles, or blocks of genes over tens to hundreds of meters. This suggests that selection overrides migration. RAPD polymorphism appears to parallel allozyme diversity which is climatically adaptive and driven by natural selection in the same subpopulations at the microsite.     

Range‐wide population structure of European sea bass Dicentrarchus labrax

The euryhaline European sea bass Dicentrarchus labrax L., inhabiting the coasts of the eastern Atlantic Ocean and Mediterranean Sea, has had many opportunities for differentiation throughout its large natural range. However, evidence for this has been incompletely documented geographically and with an insufficient number of markers. Therefore, its full range was sampled at 22 sites and individuals were genotyped with a suite of mapped markers, including 14 microsatellite loci (N = 536) and 46 neutral or gene‐linked single nucleotide polymorphisms (SNPs; N = 644). We confirm that the Atlantic and Mediterranean basins harbour two distinct lineages. Within the Atlantic Ocean no pattern was obvious based on the microsatellite and SNP genotypes, except for a subtle difference between South‐eastern and North‐eastern Atlantic sea bass attributed to limited introgression of alleles of Mediterranean origin. SNP genotypes of the Mediterranean lineage differentiated into three groups, probably under the influence of geographical isolation. The Western Mediterranean group showed genetic homogeneity without evidence for outlier loci. The Adriatic group appeared as a distinct unit. The Eastern Mediterranean group showed a longitudinal gradient of genotypes and most interestingly an outlier locus linked to the somatolactin gene. Overall, the spatial pattern fits those observed with other taxa of between‐basin segregation and within‐basin connectivity, which concurs well with the swimming capabilities of European sea bass. Evidence from a few outlier loci in this and other studies encourages further exploration of its regional connectivity and adaptive evolution.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xii. Linkage Disequilibrium in a Large Local Population

Seven hundred and three second chromosomes were extracted from a Raleigh, North Carolina population of Drosophila melanogaster in 1970. Additionally, four hundred and eighty-nine third chromosomes were extracted from a large cage population founded from the flies in the 1970 Raleigh collection. The alpha glycerol-3-phosphate dehydrogenase-1, malate dehydrogenase-1, alcohol dehydrogenase, and alpha amylase loci were studied from the second chromosomes, and the esterase-6, esterase-C, and octanol dehydrogenase loci were analyzed from the third chromosomes. Inversions, relative viability and fecundity were studied for both classes of chromosomes. The following significant findings were obtained: (1) All loci examined were polymorphic or had at least two alleles at appreciable frequencies. Analysis of the combined data from this experiment with that of Mukai, Mettler and Chigusa (1971) revealed that the frequencies of the genes in the second chromosomes collected in early August were approximately the same over three years. (2) Linkage disequilibria between and among isozyme genes inter se were not detected except in a few cases which can be considered due to non-random sampling. (3) Linkage disequilibria between isozyme genes and polymorphic inversions were detected when the recombination values between the breakage points of the inversions and the genes in question were small. In only a few cases, were second and third order linkage disequilibria including polymorphic inversions detected. (4) Evidence for either variation among genotypes within loci or cumulative effects of heterozygosity was found for viability and fecundity. As a result of these findings, it was tentatively concluded that although selection might be perceptibly operating on some polymorphic isozyme loci, most of the polymorphic isozyme genes are selectively neutral or near-neutral in the populations studied.     

Estimation of Hardy-Weinberg and pairwise disequilibrium in the apolipoprotein AI-CIII-AIV gene cluster.

Departures from Hardy-Weinberg (HW) equilibria and pairwise disequilibria were estimated in a sample of unrelated healthy individuals typed for six RFLPs in the apo AI-CIII-AIV gene region. The sample was composed of males and females, selected for health, from two populations, those of exclusively French-Canadian (FC) and those of some non-French-Canadian (NFC) ancestry. An approach suggested by Weir and Cockerham, which includes estimates of nonrandom association (disequilibria) between three and four alleles at two loci as well as the traditional associations between two alleles, at two loci was used. The pattern of departures from HW equilibria suggested that the genetic structures of the FC and NFC are different. Departure from HW equilibrium at an RFLP locus could not be predicted from information about other loci in the same gene region. Nonrandom associations were also evident from the pairwise analyses. Two pairs of loci had significant diallelic disequilibria, while two other pairs had significant triallelic disequilibria. All of the RFLP pairs had at least one measure of disequilibrium at its maximum value determined by allele frequencies. Inferences about pairwise disequilibria depended on the statistical approach used. Sizes of the pairwise disequilibria were not correlated with the physical distance between loci. The impact of these disequilibria on RFLP-phenotype association studies is discussed.     

The Genetics of DROSOPHILA SUBOBSCURA Populations. IX. Studies on Linkage Disequilibrium in Four Natural Populations

Gametic frequencies were obtained in four natural populations of D. sub-obscura by extracting wild chromosomes and subsequently analyzing them for inversions and allozymes. The genes Lap and Pept-1, both located within the same inversions of chromosome O, were found in striking nonrandom associations with them of the same kind and degree in all populations studied. On the contrary, the gene Acph, also located within the previously mentioned inversions, was found in linkage disequilibrium with them only in two populations and of opposite directions. This is also the case for the genes Est-9 and Hk, both located within chromosome E inversions. While the gene Est-9 was in strong linkage disequilibrium with the inversions, of the same kind and degree in all populations studied, Hk was found to be in linkage equilibrium. Allele frequencies for the 29 genes studied do not show geographical variation except for the genes Lap, Pept-1 and Est-9, the ones found in linkage disequilibria with the geographically varying gene arrangements. Although mechanical or historical explanations for these equilibria cannot be ruled out, these data cannot be explained satisfactorily by the "middle gene explanation," which states that loci displaying such linkage disequilibria are the ones located near the break points of inversions, while the ones displaying linkage equilibria with them are located in the middle of them. There is no evidence for consistent linkage disequilibria between pairs of loci, except for the closely linked genes of the complex locus, Est-9. This would imply, if it is not a peculiarity of the Est-9 complex, that the linkage disequilibria are found only between very closely linked loci or that, for less closely linked genes, the associations are too weak to be detected by the usual samples sizes.     

Disequilibrium among several linked neutral genes in finite population I. Mean changes in disequilibrium

Formulae are developed for computing changes in expected values in a finite population of linkage disequilibrium among neutral genes from more than two loci, although the exact analysis is taken up to only six loci. An essentially haploid model is used. As with two loci, the three-locus disequilibrium declines exponentially at all generations, but for m > 3 loci a matrix has to be constructed to give joint changes in the m-locus disequilibrium and products of disequilibria with fewer loci, for example of two m/2-locus disequilibria. The asymptotic rates of change in multilocus disequilibria depend on the arrangement of genes on the chromosome as well as its total length, but the initial rate of breakdown of disequilibrium from a line cross base is less dependent on the arrangement. With equally spaced loci the asymptotic rate of breakdown of m locus disequilibrium is roughly proportional to m. Although mutation and interference are excluded from the main analysis, it is shown how they can be incorporated.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Searching the genomes of inbred mouse strains for incompatibilities that reproductively isolate their wild relatives

Identification of the genes that underlie reproductive isolation provides important insights into the process of speciation. According to the Dobzhansky-Muller model, these genes suffer disrupted interactions in hybrids due to independent divergence in separate populations. In hybrid populations, natural selection acts to remove the deleterious heterospecific combinations that cause these functional disruptions. When selection is strong, this process can maintain multilocus associations, primarily between conspecific alleles, providing a signature that can be used to locate incompatibilities. We applied this logic to populations of house mice that were formed by hybridization involving two species that show partial reproductive isolation, Mus domesticus and Mus musculus. Using molecular markers likely to be informative about species ancestry, we scanned the genomes of 1) classical inbred strains and 2) recombinant inbred lines for pairs of loci that showed extreme linkage disequilibria. By using the same set of markers, we identified a list of locus pairs that displayed similar patterns in both scans. These genomic regions may contain genes that contribute to reproductive isolation between M. domesticus and M. musculus. This hypothesis can now be tested using laboratory crosses and surveys of introgression in the wild.     

Disequilibrium among several linked neutral genes in finite population. II. Variances and covariances of disequilibria

A method is derived for computing the variances and covariances of linkage disequilibria between neutral genes in finite populations, which is based on a linear transformation of results given previously for the mean values of disequilibria. The formulae obtained are limited to moments of sixth order or less, such as the variance of the three-locus disequilibrium. It is shown that there is no covariance between any pair of disequilibria in populations starting equilibrium. The pattern of change with time in variance of the three-locus disequilibrium from populations initially in equilibrium is similar to that for two loci, except that the highest values are achieved rather earlier and are smaller.     

The Genetical Response to Natural Selection by Varied Environments. IV. Gametic Disequilibrium in Spatially Varied Environments

Gametic disequilibria between allozyme loci were related to spatial variation of the environment in caged populations of Drosophila melanogaster . Two experiments, one with flies collected at "Chateau Tahbilk," South Australia, and the other with flies from "Groningen," The Netherlands, were sampled at generations 16 and 32. Spatial variation of the environment was stimulated using three food media. Eight polymorphic allozyme loci were used to estimate gametic disequilibria from digenic combinations of allotypes. All populations were duplicated within an environment and maintained at about 2500 adults. Standardized gametic disequilibria were compared by a weighted least squares analysis of the z-transformed statistical correlation of allele frequencies. Gametic disequilibria were strongly dependent upon food niche and food-niche interactions. The effects also varied with sampling time and were similar in duplicate populations. Gametic disequilibria were most often detected in the "Groningen"-derived populations and their strength was not strongly associated with recombination fraction. Many of the disequilibria concerned unlinked loci. The strength of selection was probably considerable and populations were evolving genetic architectures which reflected niche selection by the different foods without marked genetic isolation between foods; gene frequencies did not vary between niches within a population cage.     

Patterns of introgression across an expanding hybrid zone: analysing historical patterns of gene flow using nonequilibrium approaches

Previous studies suggest that the hybrid zone between two taxa in the Piriqueta caroliniana complex in central Florida (south-eastern North America) has recently expanded with hybrids replacing parental genotypes across a broad region of the Florida peninsula. Here I use patterns of genetic disequilibria and levels of differentiation among populations to infer historical patterns of introgression across this broad hybrid zone. There were strong positive associations among taxon-specific alleles at the southern extent of hybridization, but disequilibria values were close to zero across the central and northern portions of the hybrid zone. Levels of among-population differentiation in the central portion of the hybrid zone were relatively low, and increased towards the northern extent of hybridization. The high levels of disequilibria are coincident with the sharp clines at the southern end of the hybrid zone, suggesting that there is a tension zone in this region that is maintained by selection against hybrids and dispersal from parental regions. The levels of disequilibria within populations and patterns of differentiation among populations are consistent with historical introgression and northward expansion of this hybrid zone, which may have slowed or ceased in recent generations.     

Disequilibrium among several linked neutral genes in finite population 1. mean changes in disequilibrium

Formulae are developed for computing changes in expected values in a finite population of linkage disequilibrium among neutral genes from more than two loci, although the exact analysis is taken up to only six loci. An essentially haploid model is used. As with two loci, the three-locus disequilibrium declines exponentially at all generations, but for m > 3 loci a matrix has to be constructed to give joint changes in the m-locus disequilibrium and products of disequilibria with fewer loci, for example of two $\text{m}\text{2}\text{-}\text{locus}$ disequilibria. The asymptotic rates of change in multilocus disequilibria depend on the arrangement of genes on the chromosome as well as its total length, but the initial rate of breakdown of disequilibrium from a line cross base is less dependent on the arrangement. With equally spaced loci the asymptotic rate of breakdown of m locus disequilibrium is roughly proportional to m. Although mutation and interference are excluded from the main analysis, it is shown how they can be incorporated.