Introgression despite protection: the case of native brown trout in Natura 2000 network in Italy

Genetic structure and relatedness of brown trout (Salmo trutta) populations in the drainage basin of the Ölfusá river, South-Western Iceland

Background

Lake Þingvallavatn in Iceland, a part of the river Ölfusá drainage basin, was presumably populated by brown trout soon after it formed at the end of the last Ice Age. The genetic relatedness of the brown trout in Þingvallavatn to other populations in the Ölfusá drainage basin is unknown. After the building of a dam at the outlet of the lake in 1959 brown trout catches declined, though numbers have now increased. The aim of this study was to assess effects of geographic isolation and potential downstream gene flow on the genetic structure and diversity in brown trout sampled in several locations in the western side of the watershed of River Ölfusá. We hypothesized that brown trout in Lake Þingvallavatn constituted several local spawning populations connected by occasional gene flow before the damming of the lake. We also estimated the effective population size (NE) of some of these populations and tested for signs of a recent population bottleneck in Lake Þingvallavatn.

Methods

We sampled brown trout inhabiting four lakes and 12 rivers within and near the watershed of River Ölfusá by means of electro- and net- fishing. After stringent data filtering, 2,597 polymorphic loci obtained from ddRADseq data from 317 individuals were ascertained as putative neutral markers.

Results

Overall, the genetic relatedness of brown trout in the Ölfusá watershed reflected the connectivity and topography of the waterways. Ancestry proportion analyses and a phylogenetic tree revealed seven distinct clusters, some of which corresponded to small populations with reduced genetic diversity. There was no evidence of downstream gene flow from Lake Þingvallavatn, although gene flow was observed from much smaller mountain populations. Most locations showed low NE values (i.e., ~14.6 on average) while the putative anadromous trout from River Sog and the spawning population from River Öxará, that flows into Lake Þingvallavatn, showed notably higher NE values (i.e., 71.2 and 56.5, respectively). No signals of recent population bottlenecks were detected in the brown trout of Lake Þingvallavatn.

Discussion

This is the first time that the genetic structure and diversity of brown trout in the watershed of River Ölfusá have been assessed. Our results point towards the presence of a metapopulation in the watershed of Lake Þingvallavatn, which has been influenced by restoration efforts and is now dominated by a genetic component originated in River Öxará. Many of the locations studied represent different populations. Those that are isolated in headwater streams and lakes are genetically distinct presenting low genetic diversity, yet they can be important in increasing the genetic variation in downstream populations. These populations should be considered for conservation and direct management.

Exploring Less Invasive Visual Surveys to Assess the Spatial Distribution of Endangered Mediterranean Trout Population in a Small Intermittent Stream

Monitoring the conservation status of endangered freshwater fish using less invasive methods poses challenges for ecologists and conservationists. Visual surveys have been proposed as an alternative to electrofishing, which is a standard methodology that can cause injuries, physiological stress and post-release mortality in organisms. To test the efficacy of visual methods, a study was conducted in an intermittent stream of Sardinia (Italy). Two visual methods were employed: a visual survey from streambanks (VSS) and an underwater visual survey (UVS) using cameras. The aims of this study were (1) to compare the effectiveness of these methods in detecting patch occupancy patterns and (2) to investigate the effect of environmental variables on the detection probability of Mediterranean native trout. Environmental variables characterizing pool habitats were recorded, and generalized linear models (GLMs) were employed to assess the correlation between these variables and trout presence/absence. GLM analysis revealed that UVS had higher detection probability with larger pool volume, whereas submerged macrophytes negatively affected detection probability. Detection from streambanks (VVS) was negatively affected by a high turbulence rate. In conclusion, our study suggests the utility of visual methods to describe patterns of patch occupancy of Mediterranean trout. However, methods can be differently affected by environmental variables. Therefore, monitoring programs using these methods should consider these factors to ensure a reliable description of within-stream trout distribution in intermittent streams.

Mitochondrial diversity and genetic structure of common carp (Cyprinus carpio) in Pearl River and Nandujiang River

Common carp (Cyprinus carpio) is an important valuable cyprinid in China and has been a popular cultured aquaculture species around the globe. Understanding the genetic diversity of wild native common carp not only provides basic data for the protection and utilisation of common carp resources but also assesses the effect of human activities on the genetic diversity of this species. In this study, genetic diversity and population structure of the common carp from 15 sampling populations in the Pearl River and Nandujiang River were determined using a coalescent mitochondrial locus (MLS), including mitochondrial cytochrome b gene (Cytb) and a control region (D-loop) segment. The haplotype diversity and nucleotide diversity were 0.962 and 0.00628 in the Pearl River and 0.808 and 0.00376 in the Nandujiang River, respectively. Phylogenetic and haplotype network analyses indicated that three sub-species (a) C. c. rubrofuscus, (b) C. c. haematopterus and (c) C. c. carpio all occur in both rivers. AMOVA revealed that the variation within populations (86.2%) was the main source of the total variation. Statistically significant genetic differentiation among different Pearl River populations of C. c. rubrofuscus (F_st  = 0.05-0.25) and relatively high genetic differentiation between the Nandujiang River population and the Pearl River populations (F_st  > 0.238) are apparent. Bayesian clustering analyses detected that global populations consisted of eight genetic clusters and examined that Nandujiang River population included relatively pure genetic clusters. Neutrality tests suggested that native populations experienced recent population expansion, and Extended Bayesian Skyline Plot indicated that the common carp populations likely experienced a historical expansion during 0.125-0.250 MYA. Artificial fish propagation and release, escape from fish farms and Fang Sheng may explain the invasion of non-native sub-species in many river sections, such as Laibin, Rongjiang, Huizhou, Heyuan and Zhaoqing. To conserve the native common carp populations, release station should be established to culture native common carp fry. Overall, the findings can be contributed to complementing scientific knowledge for conservation and management of the wild native common carp.

Size-age population structure of an endangered and anthropogenically introgressed northern Adriatic population of marble trout (Salmo marmoratus Cuv.): insights for its conservation and sustainable exploitation

Salmonid species are main actors in the Italian socio-ecological landscape of inland fisheries. We present novel data on the size-age structure of one of the remnant Italian populations of the critically endangered marble trout Salmo marmoratus, which co-occurs with other stocked non-native salmonids in a large glacial river of the Lake Maggiore basin (Northern Italy-Southern Switzerland). Like other Italian native trout populations, the Toce River marble trout population is affected by anthropogenic introgression with the non-native brown trout S. trutta. Our sample includes 579 individuals, mainly collected in the Toce River main channel. We estimated the length-weight relationship, described the population size-age structure, estimated the age-specific growth trajectories, and fit an exponential mortality model. A subset of the sample was also used to measure numerical and biomass density. The estimated asymptotic maximum length is ~105 cm total length (TL). Mean length at first maturity is ~55 cm TL, and mean length at maximum yield per recruit is ~68 cm TL. Approximately 45-70% of the population are estimated to die annually, along with a fishing annual mortality of ~37%, with an exploitation ratio of ~0.5. The frequency distribution of length classes in a sample collected by angling shows that ~80% of the individuals that could be retained according to the current recreational fishing regulations likely never reproduced, and large fish disproportionally contributing to recruitment are fished and retained. We identify possible overfishing risks posed by present regulations, and propose updated harvest-slot length limits to mitigate such risks. More detailed and long-term datasets on this system are needed to more specifically inform the fishery management and monitor the effects of any change in the management strategy on the size-age structure of the marble trout population of the Toce River.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Identification of Ancestry Informative Markers in Mediterranean Trout Populations of Molise (Italy): A Multi-Methodological Approach with Machine Learning

Brown trout (Salmo trutta), like many other freshwater species, is threated by the release in its natural environment of alien species and the restocking with allochthonous conspecific stocks. Many conservation projects are ongoing and several morphological and genetic tools have been proposed to support activities aimed to restore genetic integrity status of native populations. Nevertheless, due to the complexity of degree of introgression reached up after many generations of crossing, the use of dichotomous key and molecular markers, such as mtDNA, LDH-C1* and microsatellites, are often not sufficient to discriminate native and admixed specimens at individual level. Here we propose a reduced panel of ancestry-informative SNP markers (AIMs) to support on field activities for Mediterranean trout management and conservation purpose. Starting from the genotypes data obtained on specimens sampled in the main two Molise’s rivers (Central-Southern Italy), a 47 AIMs panel was identified and validated on simulated and real hybrid population datasets, mainly through a Machine Learning approach based on Random Forest classifier. The AIMs panel proposed may represent an interesting and cost-effective tool for monitoring the level of introgression between native and allochthonous trout population for conservation purpose and this methodology could be also applied in other species.

The Italian Alpine and Subalpine trouts: Taxonomy, Evolution, and Conservation

Simple Summary

In a great part of the world, trout fishing has long inspired human spiritual ideals of immersion into nature and recreation, far removed from the fast-encroaching urbanization. Concurrently, these values and emotions fueled a white-hot business, establishing a florid market of outdoor recreation. Since the 20th century, the trout-culture industry strived to provide anglers with fishing entertainment by stocking massive amounts of non-native trouts in dozens of countries, irrespective of the lakes’ and rivers’ carrying capacity. This had dire consequences on the structural and functional diversity of these ecosystems. “Trout wars” sparked throughout the world between the promoters of stocking activities and the promoters of “wild trout management” and ethics. The “Italian trout war” has been fought on the harsh battleground of trout taxonomy, ecology, distribution, and native vs. non-native interfertile species. Northern Italy, home to the Italian Alpine and subalpine trouts and economic center of the national trout-culture and stocking industry, was particularly affected by this clash. We review here the state of art of this ongoing debate, outlining our scientific view of the taxonomy, evolution, distribution, and sustainable management of the native Italian trouts of northern Italy.

Abstract

During the last 150 years, the trout-culture industry focused on enhancing trout populations by stocking, in response to the growing anglers’ demand and the habitat degradation associated to the rapid urbanization and hydropower development. The industrialized north of Italy, home to the Italian Alpine and subalpine trout populations, is the source of most of the revenues of the national trout-culture industry. Its rapid growth, and the massive introduction of non-native interfertile trouts eroded the genetic diversity of native lineages, leading to harsh confrontations between scientists, institutions, and sportfishing associations. We review here the state of the art of the taxonomy and distribution of the northern Italian native trouts, presenting both scientific results and historical documentation. We think the only native trouts in this region are Salmo marmoratus, widespread in this region, plus small and fragmented populations of S. ghigii, present only in the South-western Alps. We strongly recommend the interruption of stocking of domesticated interfertile non-native trouts in this area, and recommend the adoption of Evolutionary Significant Units for salmonid fishery management. We further propose future research directions for a sustainable approach to the conservation and ecosystem management of the fishery resources and inland waters of northern Italy.

Genetic and Phenotypic Characteristics of the Salmo trutta Complex in Italy

Salmonid fish have become ecological and research models of study in the field of conservation genetics and genomics. Over the last decade, brown trout have received a high level of interest in research and publications. The term Salmo trutta complex is used to indicate the large number of geographic forms present in the species Salmo trutta. In Europe, the S. trutta complex consists (based on mitochondrial DNA control region analysis) of seven major evolutionary lineages: Atlantic (AT), Mediterranean (ME), Adriatic (AD), Danubian (DA), Marmoratus (MA), Duero (DU) and Tigris (TI). In several nations, the difficulty of identifying some lineages derives from their wide phenotypic and geographic plasticity and the presence of mixed lineages (due to introgressive hybridization with domestic AT populations). In Italy, the S. trutta complex populations living in the Tyrrhenian area and on the main islands (Sicily, Sardinia and Corsica) showed high genetic diversity. Currently, on the Italian Red List, the protected (near threatened) populations are the AD and ME lineages. Recent studies based on traditional (mitochondrial and nuclear markers) and NGS (next-generation sequencing) analyses have clarified some genetic differences between the populations of the Tyrrhenian region, Sicily, Sardinia and Corsica. Native populations in Sardinia belong to the AD lineage, while those living in Corsica are mainly characterized by the AD, MA and ME haplotypes. In Sicily, in the area of the Iblei mountains, an AT lineage (North African) exists. According to some authors, the term Salmo macrostigma should only be used for populations in North Africa. The use of genotyping methods based on mtDNA and nuclear markers and the latest generation sequencing techniques can improve the study of populations and evolutionary lineages in areas where there are overlaps and hybridization phenomena.

Conservation Genetics of Mediterranean Brown Trout in Central Italy (Latium): A Multi-Marker Approach

Brown trout is considered a complex of incipient species, including several phylogenetic lineages, whose natural distribution in the Mediterranean area has been altered, since the beginning of the 1900s, by massive introductions of domestic strains of Atlantic origin to support fisheries. Introduced trout naturalize in new suitable environments and extensively hybridize with native populations. Here, we characterized putatively neutral and adaptive genetic variability and population structure of Mediterranean brown trout from six river catchments in central peninsular Italy, as revealed by both mitochondrial (Control Region) and nuclear (microsatellites, LDH-C1, major histocompatibility complex) markers. We quantified the admixture of wild populations with hatchery strains and evaluated the effects of domestic trout introductions on shaping population genetics. Our analyses indicated: (1) a composite picture of genetic variability in the area, with the presence of all native Mediterranean trout mitochondrial lineages (“Adriatic”, “Mediterranean”, “marmoratus”), various frequencies of allochthonous genotypes and different rates of introgression among sampling sites; (2) asymmetric mito-nuclear introgression; (3) increasing nuclear marker diversity with increasing levels of admixture across populations; (4) strong population structure coupled with relatively low effective population size. Data allowed the identification of five management units and we propose specific actions to support ongoing and future conservation strategies within the examined area.

ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout

Background

Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples.

Results

Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections.

Conclusions

Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00698-7.

Genomic Hatchery Introgression in Brown Trout (Salmo trutta L.): Development of a Diagnostic SNP Panel for Monitoring the Impacted Mediterranean Rivers

Brown trout (Salmo trutta L.) populations have been restocked during recent decades to satisfy angling demand and counterbalance the decline of wild populations. Millions of fertile brown trout individuals were released into Mediterranean and Atlantic rivers from hatcheries with homogeneous central European stocks. Consequently, many native gene pools have become endangered by introgressive hybridization with those hatchery stocks. Different genetic tools have been used to identify and evaluate the degree of introgression starting from pure native and restocking reference populations (e.g., LDH-C* locus, microsatellites). However, due to the high genetic structuring of brown trout, the definition of the "native pool" is hard to achieve. Additionally, although the LDH-C* locus is useful for determining the introgression degree at the population level, its consistency at individual level is far from being accurate, especially after several generations were since releases. Accordingly, the development of a more powerful and cost-effective tool is essential for an appropriate monitoring to recover brown-trout-native gene pools. Here, we used the 2b restriction site-associated DNA sequencing (2b-RADseq) and Stacks 2 with a reference genome to identify single-nucleotide polymorphisms (SNPs) diagnostic for hatchery-native fish discrimination in the Atlantic and Mediterranean drainages of the Iberian Peninsula. A final set of 20 SNPs was validated in a MassARRAY® System genotyping by contrasting data with the whole SNP dataset using samples with different degree of introgression from those previously recorded. Heterogeneous introgression impact was confirmed among and within river basins, and was the highest in the Mediterranean Slope. The SNP tool reported here should be assessed in a broader sample scenario in Southern Europe considering its potential for monitoring recovery plans.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Genotyping of Two Mediterranean Trout Populations in Central-Southern Italy for Conservation Purposes Using a Rainbow-Trout-Derived SNP Array

Mediterranean trout is a freshwater fish of particular interest with economic significance for fishery management, aquaculture and conservation biology. Unfortunately, native trout populations' abundance is significantly threatened by anthropogenic disturbance. The introduction of commercial hatchery strains for recreation activities has compromised the genetic integrity status of native populations. This work assessed the fine-scale genetic structure of Mediterranean trout in the two main rivers of Molise region (Italy) to support conservation actions. In total, 288 specimens were caught in 28 different sites (14 per basins) and genotyped using the Affymetrix 57 K rainbow-trout-derived SNP array. Population differentiation was analyzed using pairwise weighted F_ST and overall F-statistic estimated by locus-by-locus analysis of molecular variance. Furthermore, an SNP data set was processed through principal coordinates analysis, discriminant analysis of principal components and admixture Bayesian clustering analysis. Firstly, our results demonstrated that rainbow trout SNP array can be successfully used for Mediterranean trout genotyping. In fact, despite an overwhelming number of loci that resulted as monomorphic in our populations, it must be emphasized that the resulted number of polymorphic loci (i.e., ~900 SNPs) has been sufficient to reveal a fine-scale genetic structure in the investigated populations, which is useful in supporting conservation and management actions. In particular, our findings allowed us to select candidate sites for the collection of adults, needed for the production of genetically pure juvenile trout, and sites to carry out the eradication of alien trout and successive re-introduction of native trout.

Pug-Headedness Anomaly in a Wild and Isolated Population of Native Mediterranean Trout Salmo trutta L., 1758 Complex (Osteichthyes: Salmonidae)

Skeletal anomalies are commonplace among farmed fish. The pug-headedness anomaly is an osteological condition that results in the deformation of the maxilla, pre-maxilla, and infraorbital bones. Here, we report the first record of pug-headedness in an isolated population of the critically endangered native Mediterranean trout Salmo trutta L., 1758 complex from Sardinia, Italy. Fin clips were collected for the molecular analyses (D-loop, LDH-C1* locus. and 11 microsatellites). A jaw index (JI) was used to classify jaw deformities. Ratios between the values of morphometric measurements of the head and body length were calculated and plotted against values of body length to identify the ratios that best discriminated between malformed and normal trout. Haplotypes belonging to the AD lineage and the genotype LDH-C1*100/100 were observed in all samples, suggesting high genetic integrity of the population. The analysis of 11 microsatellites revealed that observed heterozygosity was similar to the expected one, suggesting the absence of inbreeding or outbreeding depression. The frequency of occurrence of pug-headedness was 12.5% (two out of 16). One specimen had a strongly blunted forehead and an abnormally short upper jaw, while another had a slightly anomaly asymmetrical jaw. Although sample size was limited, variation in environmental factors during larval development seemed to be the most likely factors to trigger the deformities.

Population status and ecology of the Salmo trutta complex in an Italian river basin under multiple anthropogenic pressures

Salmonids inhabiting Mediterranean rivers are of particular concern for biodiversity conservation, as they are threatened by various stressors, including habitat alterations, overfishing, climate change, and introgressive hybridization with alien species. In the Tiber River basin (Central Italy), genetic introgression phenomena of the native Salmo cettii with the non-native Salmo trutta hinder the separate analysis of the two species, which are both included in the S. trutta complex. Little is known about the factors currently limiting the trout populations in this area, particularly with respect to climate change. With the intention of filling this gap, the aims of the current study were to (a) quantify changes in the climate and (b) analyze the distribution, status, and ecology of trout populations, in the context of changing abiotic conditions over the last decades. Fish stock assessments were carried out by electrofishing during three census periods (1998-2004, 2005-2011, and 2012-2018) at 129 sites. The trend over time of meteorological parameters provided evidence for increased air temperature and decreased rainfall. Multivariate analysis of trout densities and environmental data highlighted the close direct correlation of trout abundance with water quality, altitude, and current speed. Climate-induced effects observed over time in the sites where trout were sampled have not yet led to local extinctions or distribution shifts, indicating a marked resilience of trout, probably due to the buffering effect of intrinsic population dynamics. Decreasing body conditions over time and unbalanced age structures support the hypothesis that variations in hydraulic regime and water temperature could overcome these compensatory effects, which may lead to a severe decline in trout populations in the near future. In a climate change context, habitat availability plays a key role in the distribution of cold-water species, which often do not have the possibility to move upstream to reach their thermal optimum because of water scarcity in the upper river stretches.

Native bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix populations in the Pearl River are threatened by Yangtze River introductions as revealed by mitochondrial DNA

Bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix have been two economically important aquaculture species in China for centuries. In the past decades, bighead and silver carp have been introduced from the Yangtze River to many river systems in China, including the Pearl River, in annual, large-scale, stocking activities to enhance wild fisheries. Nonetheless, few studies have assessed the ecological or genetic impacts of such introductions on native conspecific fish populations. An mtDNA D-loop segment of 978 bp from 213 bighead carp samples from 9 populations and a 975 bp segment from 204 silver carp samples from 10 populations were obtained to evaluate genetic diversity and population integrity. Results from a haplotype network analysis revealed that most haplotypes of the Pearl River clustered with those of Yangtze River origin and only a small proportion were distinct, suggesting that both the native Pearl River bighead and silver carp populations are currently dominated by genetic material from the Yangtze River. The genetic diversity of Pearl River populations is high in both species because of this inter-population gene flow, but the diversity of native Pearl River populations is low. To preserve the native genetic diversity, stocking of non-native fingerlings should cease immediately and native Pearl River bighead and silver carp fish farms should be established. This research demonstrates the danger to native biodiversity across China because of the substantial, ongoing stock-enhancement activities without prior genetic assessment.