The evolutionary history of the goby Elacatinus puncticulatus in the tropical eastern pacific: Effects of habitat discontinuities and local environmental variability

Genomic assessment reveals signal of adaptive selection in populations of the Spotted rose snapper Lutjanus guttatus from the Tropical Eastern Pacific

Background

The lack of barriers in the marine environment has promoted the idea of panmixia in marine organisms. However, oceanographic conditions and habitat characteristics have recently been linked to genetic structure in marine species. The Tropical Eastern Pacific (TEP) is characterized by dynamic current systems and heterogeneous oceanographic conditions. The Gulf of Panama (part of the equatorial segment for the TEP) is influenced by a complex current system and heterogeneous environment, which has been shown to limit the gene flow for shoreline species. Next Generation Sequencing (NGS) has contributed to detect genetic differences in previously reported panmictic species by the assessment of loci associated with selection and to understand how selection acts affects marine populations. Lutjanus guttatus is a species distributed in the TEP for which previous studies using mitochondrial data recovered a panmictic pattern along its distributional range. In this study, we used SNP data of L. guttatus individuals sampled along its range to evaluate population genetic structure and investigate whether oceanographic factors influence the species’ genetic architecture. Finally, we assessed the role of adaptive selection by evaluating the contribution of outlier and neutral loci to genetic divergence.

Methods

The RADcap method was used to obtain 24 million paired reads for 123 individuals of L. guttatus covering nearly all its distributional area. Genetic variation was assessed using both spatial and non-spatial methods by comparing three different data sets: (i) a Combined Loci (CL dataset = 2003 SNPs); a search for putative loci under selection allowed the evaluation of (ii) Neutral Loci (NL dataset = 1858 SNPs) and (iii) Outlier Loci (OL dataset = 145 SNPs). We used the estimating effective migration surface (EEMS) approach to detect possible barriers to gene flow.

Results

Genetic differences were found in the OL dataset, showing two clusters (Northern and Southern), whereas NL showed no differences. This result may be related to the Selection-Migration balance model. The limit between the Northern and Southern groups was in the Gulf of Panama, which has been previously identified as a barrier to gene flow for other species, mainly due to its heterogeneous oceanographic conditions. The results suggest that selection plays an important role in generating genetic differences in Lutjanus guttatus. A migration corridor was detected that coincides with the Costa Rica Coastal Current that flows from Central America to the Gulf of California, allowing the homogenization of the northern population. In the Southern cluster, a migration corridor was observed with the OL from Panama to Colombia, which could be associated with the currents found in the Gulf of Panama. Genetic variation found in the OL of Lutjanus guttatus highlights the usefulness of NGS data in evaluating the role of selection in population differentiation.

Phylogeography and evolutionary history of the Panamic Clingfish Gobiesox adustus in the Tropical Eastern Pacific

The Panamic Clingfish Gobiesox adustus is widely distributed in the Tropical Eastern Pacific (TEP), from the central Gulf of California, Mexico to Ecuador, including the oceanic Revillagigedo Archipelago, and Isla del Coco. This cryptobenthic species is restricted to very shallow rocky-reef habitats. Here, we used one mitochondrial and three nuclear DNA markers from 155 individuals collected across the distribution range of the species in order to evaluate if geographically structured populations exist and to elucidate its evolutionary history. Phylogenetic analyses recovered a monophyletic group, with four well-supported, allopatric subgroups. Each subgroup corresponded to one of the following well-known biogeographic regions/provinces: 1) the Revillagigedo Archipelago, 2) the Cortez + Mexican provinces (Mexico), 3) the Panamic province (from El Salvador to Ecuador), and 4) Isla del Coco. A molecular-clock analysis showed a mean date for the divergence between clade I (the Revillagigedos and Cortez + Mexican provinces) and clade II (Panamic province and Isla del Coco) in the Pliocene, at ca. 5.33 Mya. Within clade I, the segregation between the Revillagigedos and Cortez + Mexican province populations was dated at ca. 1.18 Mya, during the Pleistocene. Within clade II, the segregation between samples of Isla del Coco and the Panamic province samples was dated at ca. 0.77 Mya, during the Pleistocene. The species tree, Bayesian species delimitation tests (BPP and STACEY), the Φ_ST, AMOVA, and the substantial genetic distances that exist between those four subgroups, indicate that they are independent evolutionary units. These cladogenetic events seem to be related to habitat discontinuities, and oceanographic and geological processes that produce barriers to gene flow for G. adustus, effects of which are enhanced by the intrinsic ecological characteristics of this species.

Smaller fish species in a warm and oxygen-poor Humboldt Current system

Climate change is expected to result in smaller fish size, but the influence of fishing has made it difficult to substantiate the theorized link between size and ocean warming and deoxygenation. We reconstructed the fish community and oceanographic conditions of the most recent global warm period (last interglacial; 130 to 116 thousand years before present) by using sediments from the northern Humboldt Current system off the coast of Peru, a hotspot of small pelagic fish productivity. In contrast to the present-day anchovy-dominated state, the last interglacial was characterized by considerably smaller (mesopelagic and goby-like) fishes and very low anchovy abundance. These small fish species are more difficult to harvest and are less palatable than anchovies, indicating that our rapidly warming world poses a threat to the global fish supply.

Checklist of marine and estuarine fishes from the AlaskaYukon Border, Beaufort Sea, to Cabo San Lucas, Mexico

This paper is a checklist of the fishes that have been documented, through both published and unpublished sources, in marine and estuarine waters, and out 200 miles, from the United States-Canadian border on the Beaufort Sea to Cabo San Lucas, Mexico. A minimum of 241 families and 1,644 species are known within this range, including both native and nonnative species. For each of these species, we include maximum size, geographic and depth ranges, whether it is native or nonnative, as well as a brief mention of any taxonomic issues.

A multi-locus approach to elucidating the evolutionary history of the clingfish Tomicodon petersii (Gobiesocidae) in the Tropical Eastern Pacific

Marine species that are widely distributed in the Tropical Eastern Pacific (TEP) has served as a model for studying biogeographic patterns resulting from the effects of intraregional habitat discontinuities and oceanographic processes on the diversification and evolution of cryptobenthic reef fishes. Tomicodon petersii, a clingfish (Gobiesocidae) endemic to the TEP, is found on very shallow rocky reefs from central Mexico to northern Peru, and in the Cocos and Galapagos islands. We evaluated the effect of likely biogeographic barriers in different parts of the TEP on the diversification process of this species. We used one mitochondrial and three nuclear DNA markers from 112 individuals collected across the distribution range of T. petersii. Our phylogenetic results showed the samples constituted a monophyletic group, with three well-supported, allopatric subgroups: in the Mexican province, the Panamic province (from El Salvador to Ecuador), and the Galapagos Islands. The split between the Mexican and more southerly clades was estimated to occur at the end of the Miocene ca. 5.74 Mya, and the subsequent cladogenetic event separating the Galapagos population from the Panamic population at the junction of the Pliocene and Pleistocene, ca. 2.85 Mya. The species tree, Bayesian species delimitation tests (BPP), STACEY, and substantial genetic distances separating these three populations indicate that these three independent evolutionary units likely include two unnamed species. The cladogenetic events that promoted the formation of those genetically differentiated groups are consistent with disruptive effects on gene flow of habitat discontinuities and oceanographic processes along the mainland shoreline in the TEP and of ocean-island isolation, in conjunction with the species intrinsic life-history characteristics.

Genetic differentiation in the genus Characodon: implications for conservation and taxonomy

The subfamily Goodeinae is a group of fishes endemic to the Mexican highlands. Most of the species are restricted to small and isolated streams or springs. Within this subfamily, the genus Characodon is the earliest diverging lineage of which three species have been described: C. lateralis, C. audax, and C. garmani, with the latter, considered extinct. Characodon lateralis and C. audax are classed as endangered, and have been the subject of taxonomic controversy since their description: previous studies have recognized a genetic differentiation in two groups separated by the El Salto waterfall, but morphological analyses contradict these genetic results. We perform a phylogeographic study using the mitochondrial cytb gene and d-loop region to elucidate the evolutionary history of C. lateralis and C. audax. The results with both markers show the presence of two highly differentiated haplogroups; one distributed north and the other distributed south of the waterfall, with genetic distances of 1.7 and 13.1% with cytb and d-loop respectively, and divergence calculated to have occurred 1.41 Mya. Significant genetic structure was found within each haplogroup and suggests the existence of at least four Evolutionary Significant Units (ESUs) within the examined populations. The possible processes identified as contributing to the formation of differentiated genetic groups are isolation, low population size, recurrent bottlenecks, and the strong sexual selection exhibited by the genus.

Drivers for genetic structure at different geographic scales for Pacific red snapper (Lutjanus peru) and yellow snapper (Lutjanus argentiventris) in the tropical eastern Pacific

The tropical eastern Pacific (TEP) is a highly dynamic region and a model system to study how habitat discontinuities affect the distribution of shorefishes, particularly for species that display ontogenetic habitat shifts, including snappers (Lutjanidae). To evaluate the genetic structure of the Pacific red snapper (Lutjanus peru) and the yellow snapper (Lutjanus argentiventris) throughout their distribution range along the TEP, 13 and 11 microsatellite loci were analysed, respectively. The genetic diversity of L. peru (N = 446) and L. argentiventris (N = 170) was evaluated in 10 and 5 localities, respectively, showing slightly higher but non-significant values in the Gulf of California for both species. The genetic structure analysis identified the presence of significant genetic structure in both species, but the locations of the identified barriers for the gene flow differed between species. The principal driver for the genetic structure at large scales >2500 km was isolation by distance. At smaller scales (<250 km), the habitat discontinuity for juveniles and adults and the environmental differences throughout the distribution range represented potential barriers to gene flow between populations for both species.

Unraveling hierarchical genetic structure in a marine metapopulation: A comparison of three high-throughput genotyping approaches

Marine metapopulations often exhibit subtle population structure that can be difficult to detect. Given recent advances in high-throughput sequencing, an emerging question is whether various genotyping approaches, in concert with improved sampling designs, will substantially improve our understanding of genetic structure in the sea. To address this question, we explored hierarchical patterns of structure in the coral reef fish Elacatinus lori using a high-resolution approach with respect to both genetic and geographic sampling. Previously, we identified three putative E. lori populations within Belize using traditional genetic markers and sparse geographic sampling: barrier reef and Turneffe Atoll; Glover's Atoll; and Lighthouse Atoll. Here, we systematically sampled individuals at ~10 km intervals throughout these reefs (1,129 individuals from 35 sites) and sequenced all individuals at three sets of markers: 2,418 SNPs; 89 microsatellites; and 57 nonrepetitive nuclear loci. At broad spatial scales, the markers were consistent with each other and with previous findings. At finer spatial scales, there was new evidence of genetic substructure, but our three marker sets differed slightly in their ability to detect these patterns. Specifically, we found subtle structure between the barrier reef and Turneffe Atoll, with SNPs resolving this pattern most effectively. We also documented isolation by distance within the barrier reef. Sensitivity analyses revealed that the number of loci (and alleles) had a strong effect on the detection of structure for all three marker sets, particularly at small spatial scales. Taken together, these results illustrate empirically that high-throughput genotyping data can elucidate subtle genetic structure at previously-undetected scales in a dispersive marine fish.