Reduced Global Genetic Differentiation of Exploited Marine Fish Species

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

Long-distance migrations and seasonal movements of meagre (Argyrosomus regius), a large coastal predator, along the Iberian Peninsula coast

BACKGROUND

The meagre, Argyrosomus regius, is a large coastal predatory fish inhabiting waters from the north-eastern Atlantic and Mediterranean Sea, where it is targeted by commercial and recreational fisheries. Previous genetic studies have found an unexpectedly high population differentiation not only between the Atlantic and the Mediterranean, but also along the Atlantic coast. However, the reasons underpinning this genetic barrier remained unclear. Likewise, even though the species is amongst the world's largest marine teleosts, knowledge about its movement ecology and migratory behaviour remains notably scarce, and primarily reliant on fisheries-dependent data.

METHODS

In this study, we used a combination of acoustic telemetry and pop-up satellite archival tags to investigate the movements of 22 adult meagre (70-143 cm total length) along the Southwestern coast of the Iberian Peninsula.

RESULTS

Our results strongly suggest that the previously reported genetic differentiation is not maintained by limited adult dispersal/movement, as hypothesized. On the contrary, we documented some of the longest individual annual migrations ever recorded for a coastal teleost, up to > 2000 km, with frequent back-and-forth movements between the West and Southern Iberian coasts. Moreover, their detected regional movement patterns support the existence of a marked seasonal behavioural shift, with individuals being less active and moving to deeper waters during winter, and are consistent with spawning philopatry associated to their summer reproductive movements. Finally, we identified putative aggregation areas that may harbour important feeding/overwintering grounds.

CONCLUSIONS

These findings shed new light on the movement and behaviour patterns of meagre that may be of particular importance for the conservation and spatial management of this species throughout its range, and open the door to further research on functional connectivity.

Population genetic structure of Randall's threadfin bream Nemipterus randalli in Indian waters based on mitochondrial and nuclear gene sequences

Nemipterus randalli, commonly known as Randall's threadfin bream, is a commercially important marine finfish. Understanding its genetic structure is critical to effective management and conservation efforts. Previous investigations on population structure in this species were limited by geographic coverage. In this study, we utilized the mitochondrial Cytochrome b gene and nuclear Ribosomal protein gene intron Rp S7 sequences to investigate the population genetic structure, demography and genetic diversity of N. randalli along Indian waters. Our results revealed high haplotype diversity but low nucleotide diversity. AMOVA revealed that the variation among the population was highly significant. Hierarchical AMOVA provided further evidence of significant genetic differentiation between the west and east coasts, which was corroborated by the Bayesian tree and the median-joining network diagram. The mtDNA sequences revealed significant genetic structure between populations based on fixation index analysis following the isolation-by-distance model. Furthermore, the neutrality test and mismatch analysis suggest that N. randalli populations may have experienced a population expansion. However, nuclear marker RpS7, showed a high level of polymorphism, which obscured the population structuring observed with the mitochondrial marker. Consequently, concordant results were not obtained when comparing the mitochondrial and nuclear DNA sequences. The strong genetic differentiation between the east and west coast observed using mitochondrial marker could be attributed to a combination of geographic and environmental factors. These findings lay the groundwork for developing effective conservation and management strategies for N. randalli, considering its genetic structure.

Anti-predator defences are linked with high levels of genetic differentiation in frogs

Predator-prey interactions have been suggested as drivers of diversity in different lineages, and the presence of anti-predator defences in some clades is linked to higher rates of diversification. Warning signals are some of the most widespread defences in the animal world, and there is evidence of higher diversification rates in aposematic lineages. The mechanisms behind such species richness, however, are still unclear. Here, we test whether lineages that use aposematism as anti-predator defence exhibit higher levels of genetic differentiation between populations, leading to increased opportunities for divergence. We collated from the literature more than 3000 pairwise genetic differentiation values across more than 700 populations from over 60 amphibian species. We find evidence that over short geographical distances, populations of species of aposematic lineages exhibit greater genetic divergence relative to species that are not aposematic. Our results support a scenario where the use of warning signals could restrict gene flow, and suggest that anti-predator defences could impact divergence between populations and potentially have effects at a macro-evolutionary scale.

Multidisciplinary estimates of connectivity and population structure suggest the use of multiple units for the conservation and management of meagre, Argyrosomus regius

Information on population structure and connectivity of targeted species is key for proper implementation of spatial conservation measures. We used a combination of genomics, biophysical modelling, and biotelemetry to infer the population structure and connectivity of Atlantic meagre, an important fisheries resource throughout its distribution. Genetic samples from previously identified Atlantic spawning locations (Gironde, Tejo, Guadalquivir, Banc d'Arguin) and two additional regions (Algarve and Senegal) were analysed using genome-wide SNP-genotyping and mitochondrial DNA analyses. Biophysical models were conducted to investigate larval dispersal and connectivity from the known Atlantic spawning locations. Additionally, thirteen fish were double-tagged with biotelemetry transmitters off the Algarve (Portugal) to assess movement patterns and connectivity of adult individuals. This multidisciplinary approach provided a robust overview of meagre population structure and connectivity in the Atlantic. Nuclear SNP-genotyping showed a clear differentiation between the European and African populations, with significant isolation of the few known Atlantic spawning sites. The limited level of connectivity between these subpopulations is potentially driven by adults, capable of wide-ranging movements and connecting sites 500 km apart, as evidenced by tagging studies, whilst larval dispersal inferred by modelling is much more limited (average of 52 km; 95% of connectivity events up to 174 km). Our results show sufficient evidence of population structure, particularly between Africa and Europe but also within Europe, for the meagre to be managed as separate stocks. Additionally, considering the low degree of larvae connectivity, the implementation of marine protected areas in key spawning sites could be crucial towards species sustainability.

Evolutionary origin of genomic structural variations in domestic yaks

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Lack of genetic structure in euryhaline Chelon labrosus from the estuaries under anthropic pressure in the Southern Bay of Biscay to the coastal waters of the Mediterranean Sea

Over the last decade, xenoestrogenic effects have been reported in populations of thicklip grey mullet Chelon labrosus from contaminated estuaries in the Bay of Biscay, resulting in intersex condition. To understand the level of gene flow in individuals of different Basque estuaries microsatellite markers were used to evaluate the population structure and connectivity of C. labrosus from estuaries of the Basque coast. 46 microsatellites were tested and 10 validated for the analysis of 204 individuals collected from 5 selected Basque estuaries and 2 outgroups in the Bay of Cadiz and Thermaic Gulf. The polymorphic microsatellites revealed 74 total alleles, 2-19 alleles per locus. The mean observed heterozygosity (0.49 ± 0.02) was lower than the expected one (0.53 ± 0.01). There was no evidence of genetic differentiation (F_ST = 0.0098, P = 0.0000) among individuals or sites. Bayesian clustering analysis revealed a single population in all sampled locations. The results of this study indicate widespread genetic homogeneity and panmixia of C. labrosus across the current sampling areas spanning the Atlantic and Mediterranean basins. The hypothesis of panmixia could therefore be well supported so individuals inhabiting estuaries with high prevalence of intersex condition should be considered as members of the same single genetic group as those inhabiting adjacent estuaries without incidence of xenoestrogenicity.

Oceanographic connectivity explains the intra-specific diversity of mangrove forests at global scales

The distribution of mangrove intra-specific biodiversity can be structured by historical demographic processes that enhance or limit effective population sizes. Oceanographic connectivity (OC) may further structure intra-specific biodiversity by preserving or diluting the genetic signatures of historical changes. Despite its relevance for biogeography and evolution, the role of oceanographic connectivity in structuring the distribution of mangrove's genetic diversity has not been addressed at global scale. Here we ask whether connectivity mediated by ocean currents explains the intra-specific diversity of mangroves. A comprehensive dataset of population genetic differentiation was compiled from the literature. Multigenerational connectivity and population centrality indices were estimated with biophysical modeling coupled with network analyses. The variability explained in genetic differentiation was tested with competitive regression models built upon classical isolation-by-distance (IBD) models considering geographic distance. We show that oceanographic connectivity can explain the genetic differentiation of mangrove populations regardless of the species, region, and genetic marker (significant regression models in 95% of cases, with an average R-square of 0.44 ± 0.23 and Person's correlation of 0.65 ± 0.17), systematically improving IBD models. Centrality indices, providing information on important stepping-stone sites between biogeographic regions, were also important in explaining differentiation (R-square improvement of 0.06 ± 0.07, up to 0.42). We further show that ocean currents produce skewed dispersal kernels for mangroves, highlighting the role of rare long-distance dispersal events responsible for historical settlements. Overall, we demonstrate the role of oceanographic connectivity in structuring mangrove intra-specific diversity. Our findings are critical for mangroves' biogeography and evolution, but also for management strategies considering climate change and genetic biodiversity conservation.

Genetic species identification and population structure of grouper Epinephelus coioides (Hamilton, 1822) collected from fish markets along the Persian Gulf and the Oman Sea

Many ecologically important and valuable fisheries marine species have been misidentified in terms of both the statistical data and market demand. Correct identification at the species level and the population genetic structure of the orange-spotted grouper (Epinephelus coioides), a precious fish in the Persian Gulf and the Oman Sea, was tested using mitochondrial cytochrome oxidase subunit I (DNA barcoding) and D-loop sequencing. The results revealed that the Epinephelus species found in the region, including E. coioides, E. bleekeri, E. polylepis, and E. chlorostigma were all mistakenly grouped together and identified as only E. coioides. Moreover, the analysis of molecular variance (AMOVA) of E. coioides samples using the D-loop showed a significantly unique genetic structure (Φ_ST = 0.068, p < 0.001) within the E. coioides population throughout the Persian Gulf and the Oman Sea, with the pairwise genetic difference between sampling locations in UAE and the Iranian coast. Moreover, D-loop sequences analysis showed two distinct haplotype groups scattered among the sampling locations, which did not correlate with the geographic distance between the sampling locations. These findings indicate that the issue of misidentification should be highlighted in the management and conservation of E. coioides. As this type of misidentification is likely to happen to other threatened marine species as well, the efficacy of using genetic markers for the correct identification, both at the species and the population level, is vital.

MHC class I-α population differentiation in a commercial fish, the European sea bass (Dicentrarchus labrax)

Identifying population structuring in highly fecund marine species with high dispersal rates is challenging, but critical for conservation and stock delimitation for fisheries management. European sea bass (Dicentrarchus labrax) is a commercial species of fisheries and aquaculture relevance whose stocks are declining in the North Atlantic, despite management measures to protect them and identifying their fine population structure is needed for managing their exploitation. As for other marine fishes, neutral genetic markers indicate that eastern Atlantic sea bass form a panmictic population and is currently managed as arbitrarily divided stocks. The genes of the major histocompatibility complex (MHC) are key components of the adaptive immune system and ideal candidates to assess fine structuring arising from local selective pressures. We used Illumina sequencing to characterise allelic composition and signatures of selection at the MHC class I-α region of six D. labrax populations across the Atlantic range. We found high allelic diversity driven by positive selection, corresponding to moderate supertype diversity, with 131 alleles clustering into four to eight supertypes, depending on the Bayesian information criterion threshold applied, and a mean number of 13 alleles per individual. Alleles could not be assigned to particular loci, but private alleles allowed us to detect regional genetic structuring not found previously using neutral markers. Our results suggest that MHC markers can be used to detect cryptic population structuring in marine species where neutral markers fail to identify differentiation. This is particularly critical for fisheries management, and of potential use for selective breeding or identifying escapes from sea farms.

Resolving the origins of invertebrate colonists in the Yangtze River Estuary with molecular markers: Implications for ecological connectivity

Understanding connectivity over different spatial and temporal scales is fundamental for managing of ecological systems. However, controversy exists for wintertime ecological connectivity between the Yangtze River Estuary (YRE) and inner southwestern Yellow Sea. Here, we investigated ecological connectivity between the YRE and inner southwestern Yellow Sea in wintertime by precisely pinpointing the source of the newly colonized populations of a winter-spawning rocky intertidal invertebrate, Littorina brevicula (Philippi, 1844), on artificial structures along the coast of the Yangtze River Delta (YRD) using mitochondrial ND6 sequences and microsatellite data. Clear phylogeographic and genetic differentiation were detected between natural rocky populations south and north of the YRE, which resulted from the lack of hard substrate for rocky invertebrates in the large YRD coast. For the newly colonized populations on the coast of YRD, most individuals (98%) to the south of ~33.5°N were from natural rocky populations to the south of the YRE and most of those (94%) to the north of ~33.5°N were from the northern natural rocky populations, which demonstrated strong ecological connectivity between the inner southwestern Yellow Sea and the YRE in winter time. We presented the first genetic evidence that demonstrated a northward wintertime coastal current in the inner southwestern Yellow Sea, and precisely illustrated the boundary of the coastal current recently proposed by numerical experiment. These results indicated that the YRE serves as an important source of materials and energy for the inner southwestern Yellow Sea in winter, which can be crucial for the function of the Yellow Sea ecosystem.