Jonathan Sandoval-Castillo

Permanent Genetic Resources added to the Molecular Ecology Resources Database 1 February 2010-31 March 2010.

This article documents the addition of 142 microsatellite marker loci to the Molecular Ecology Resources database. Loci were developed for the following species: Agriophyllum squarrosum, Amazilia cyanocephala, Batillaria attramentaria, Fungal strain CTeY1 (Ascomycota), Gadopsis marmoratus, Juniperus phoenicea subsp. turbinata, Liriomyza sativae, Lupinus polyphyllus, Metschnikowia reukaufii, Puccinia striiformis and Xylocopa grisescens. These loci were cross‐tested on the following species: Amazilia beryllina, Amazilia candida, Amazilia rutila, Amazilia tzacatl, Amazilia violiceps, Amazilia yucatanensis, Campylopterus curvipennis, Cynanthus sordidus, Hylocharis leucotis, Juniperus brevifolia, Juniperus cedrus, Juniperus osteosperma, Juniperus oxycedrus, Juniperus thurifera, Liriomyza bryoniae, Liriomyza chinensis, Liriomyza huidobrensis and Liriomyza trifolii.

Deep Mitochondrial Divergence in Baja California Populations of an Aquilopelagic Elasmobranch: The Golden Cownose Ray

Assessing the realized effect of dispersal in the genetic makeup of a species has significant evolutionary, ecological, and economical consequences. Here, we investigate the genetic diversity and population differentiation in the aquilopelagic golden cownose ray Rhinoptera steindachneri from the Gulf of California (GC) and the Pacific coast of Baja California (PCBC) using the mitochondrial NADH2 gene. Low levels of genetic diversity were found with only 4 polymerase chain reaction-restriction fragment length polymorphism haplotypes among 76 specimens. Pacific coast organisms were fixed for a unique haplotype not shared with rays from the gulf; 92% of GC rays possessed a single NADH2 haplotype not found in the Pacific. This produced significant differentiation between the GC and the PCBC (Φ(CT) = 0.972, P < 0.001). A pronounced phylogeographic pattern was found in which GC haplotypes were reciprocally monophyletic relative to a very divergent Pacific lineage (d = 10%). Our results indicate that despite high dispersal potential, GC and PCBC golden cownose ray populations are characterized by highly divergent mitochondrial lineages. Although more evidence is needed to corroborate the genetic isolation and systematic status of PCBC and GC golden cownose rays, our results suggest a possible cryptic species in the region.

Strong Population Structure and Shallow Mitochondrial Phylogeny in the Banded Guitarfish, Zapteryx exasperata (Jordan y Gilbert, 1880), from the Northern Mexican Pacific

The northern Mexican Pacific (NMP), the Gulf of California (GC), and Baja California have been recognized as an ecological and evolutionarily dynamic region having experienced significant tectonic and climatic changes leading to the diversification of terrestrial and marine biotas. Zapteryx exasperata is a predominant ray caught in the artisanal fisheries of the NMP. Morphometric and reproductive differences between rays from the GC and the Pacific coast of Baja California (PCBC) regions suggest the presence of distinct populations. We investigate whether this distinction correlates with differences in genetic diversity and differentiation using sequences of the mitochondrial nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene and the noncoding control region (CR) in 63 specimens. Contrary to our expectations, ND2 bore significantly more diversity (h = 0.76) than CR (h = 0.39). Geographic patterns of diversity of CR were opposite to those of ND2, with GC being significantly less (ND2) and more (CR) diverse than PCBC. The diversity of concatenated haplotypes was high (h = 0.84). Low nucleotide diversity suggests the recent coancestry of haplotypes. Marked genetic structure (Φst = 0.23, P < 0.0001) revealed the existence of reproductive isolation and limited matrilineal gene flow between GC and PCBC, which correlates with their phenotypic distinction. These results suggest the influence of factors such as female reproductive philopatry, and ecological or historical vicariant barriers to gene flow. Our results point to the existence of a distinct management unit of banded guitarfish in each region, and add to the increasing evidence attesting to the diversifying nature of this evolutionarily dynamic region.

Genome-wide data delimits multiple climate-determined species ranges in a widespread Australian fish, the golden perch (Macquaria ambigua)

Species range limits often fluctuate in space and time in response to variation in environmental factors and to gradual niche evolution due to changes in adaptive traits. We used genome-wide data to investigate evolutionary divergence and species range limits in a generalist and highly dispersive fish species that shows an unusually wide distribution across arid and semi-arid regions of Australia. We generated ddRAD data (18,979 filtered SNPs and 1.725million bp of sequences) for samples from 27 localities spanning the native range of golden perch, Macquaria ambigua (Teleostei; Percichthyidae). Our analytical framework uses population genomics to assess connectivity and population structure using model-based and model-free approaches, phylogenetics to clarify evolutionary relationships, and a coalescent-based Bayesian species delimitation method to assess statistical support of inferred species boundaries. Addressing uncertainties regarding range limits and taxonomy is particularly relevant for this iconic Australian species because of the intensive stocking activities undertaken to support its recreational fishery and its predicted range shifts associated with ongoing climate change. Strong population genomic, phylogenetic, and coalescent species delimitation support was obtained for three separately evolving metapopulation lineages, each lineage should be considered a distinct cryptic species of golden perch. Their range limits match the climate-determined boundaries of main river basins, despite the ability of golden perch to cross drainage divides. We also identified cases suggestive of anthropogenic hybridization between lineages due to stocking of this recreationally important fish, as well as a potential hybrid zone with a temporally stable pattern of admixture. Our work informs on the consequences of aridification in the evolution of aquatic organisms, a topic poorly represented in the literature. It also shows that genome-scale data can substantially improve and rectify inferences about taxonomy, hybridization and conservation management previously proposed by detailed genetic studies.

Oceanography promotes self-recruitment in a planktonic larval disperser.

The application of high-resolution genetic data has revealed that oceanographic connectivity in marine species with planktonic larvae can be surprisingly limited, even in the absence of major barriers to dispersal. Australia’s southern coast represents a particularly interesting system for studying planktonic larval dispersal, as the hydrodynamic regime of the wide continental shelf has potential to facilitate onshore retention of larvae. We used a seascape genetics approach (the joint analysis of genetic data and oceanographic connectivity simulations) to assess population genetic structure and self-recruitment in a broadcast-spawning marine gastropod that exists as a single meta-population throughout its temperate Australian range. Levels of self-recruitment were surprisingly high, and oceanographic connectivity simulations indicated that this was a result of low-velocity nearshore currents promoting the retention of planktonic larvae in the vicinity of natal sites. Even though the model applied here is comparatively simple and assumes that the dispersal of planktonic larvae is passive, we find that oceanography alone is sufficient to explain the high levels of genetic structure and self-recruitment. Our study contributes to growing evidence that sophisticated larval behaviour is not a prerequisite for larval retention in the nearshore region in planktonic-developing species.

Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc, the greenlip abalone (Haliotis laevigata), along a longitudinal environmental gradient

Populations of broadcast spawning marine organisms often have large sizes and are exposed to reduced genetic drift. Under such scenarios, strong selection associated with spatial environmental heterogeneity is expected to drive localized adaptive divergence, even in the face of connectivity. We tested this hypothesis using a seascape genomics approach in the commercially important greenlip abalone (Haliotis laevigata). We assessed how its population structure has been influenced by environmental heterogeneity along a zonal coastal boundary in southern Australia linked by strong oceanographic connectivity. Our data sets include 9,109 filtered SNPs for 371 abalones from 13 localities and environmental mapping across ~800 km. Genotype–environment association analyses and outlier tests defined 8,786 putatively neutral and 323 candidate adaptive loci. From a neutral perspective, the species is better represented by a metapopulation with very low differentiation (global FST = 0.0081) and weak isolation by distance following a stepping‐stone model. For the candidate adaptive loci, however, model‐based and model‐free approaches indicated five divergent population clusters. After controlling for spatial distance, the distribution of putatively adaptive variation was strongly correlated to selection linked to minimum sea surface temperature and oxygen concentration. Around 80 candidates were annotated to genes with functions related to high temperature and/or low oxygen tolerance, including genes that influence the resilience of abalone species found in other biogeographic regions. Our study includes a documented example about the uptake of genomic information in fisheries management and supports the hypothesis of adaptive divergence due to coastal environmental heterogeneity in a connected metapopulation of a broadcast spawner.

Fetal mummification in silky shark (Carcharhinus falciformis) from the Gulf of California, Mexico

The aim of this work was the description of fetal mummification in silky shark. During July of 2001 in Santa Maria fishing camp, Gulf of California Mexico, fourteen pregnant female of silky sharks were captured. In all of them the embryos were found to be dehydrated. However, total observations were possible in only two litters. Six embryos in each analyzed litters were found. All the embryos were mummified, complete and without maceration. Because the litters were totally mummificated in several females, it was suggested that fetal mummification was produced by the environmental condition.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray‐Darling Basin (MDB), Australia, using a genome‐wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype–environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Genome-wide SNPs resolve a key conflict between sequence and allozyme data to confirm another threatened candidate species of river blackfishes (Teleostei: Percichthyidae: Gadopsis)

Conflicting results from different molecular datasets have long confounded our ability to characterise species boundaries. Here we use genome-wide SNP data and an expanded allozyme dataset to resolve conflicting systematic hypotheses on an enigmatic group of fishes (Gadopsis, river blackfishes, Percichthyidae) restricted to southeastern Australia. Previous work based on three sets of molecular markers: mtDNA, nuclear intron DNA and 51 allozyme loci was unable to clearly resolve the status of a putative fifth candidate species (SWV) within Gadopsis marmoratus. Resolving the taxonomic status of candidate species SWV is particularly critical as based on IUCN criteria this taxon would be considered Critically Endangered. After all filtering steps we retained a subset of 10,862 putatively unlinked SNP loci for population genetic and phylogenomic analyses. Analyses of SNP loci based on maximum likelihood, fastSTRUCTURE and DAPC were all consistent with the previous and updated allozyme results supporting the validity of the candidate Gadopsis species SWV. Immediate conservation actions should focus on preventing take by anglers, protection of water resources to sustain perennial reaches and drought refuge pools, and aquatic and riparian habitat protection and improvement. In addition, a formal morphological taxonomic review of the genus Gadopsis is urgently required.

Mitochondrial DNA is unsuitable to test for isolation by distance

Tests for isolation by distance (IBD) are the most commonly used method of assessing spatial genetic structure. Many studies have exclusively used mitochondrial DNA (mtDNA) sequences to test for IBD, but this marker is often in conflict with multilocus markers. Here, we report a review of the literature on IBD, with the aims of determining (a) whether significant IBD is primarily a result of lumping spatially discrete populations, and (b) whether microsatellite datasets are more likely to detect IBD when mtDNA does not. We also provide empirical data from four species in which mtDNA failed to detect IBD by comparing these with microsatellite and SNP data. Our results confirm that IBD is mostly found when distinct regional populations are pooled, and this trend disappears when each is analysed separately. Discrepancies between markers were found in almost half of the studies reviewed, and microsatellites were more likely to detect IBD when mtDNA did not. Our empirical data rejected the lack of IBD in the four species studied, and support for IBD was particularly strong for the SNP data. We conclude that mtDNA sequence data are often not suitable to test for IBD, and can be misleading about species’ true dispersal potential. The observed failure of mtDNA to reliably detect IBD, in addition to being a single-locus marker, is likely a result of a selection-driven reduction in genetic diversity obscuring spatial genetic differentiation.