Upwelling-level acidification and pH/pCO2 variability moderate effects of ocean acidification on brain gene expression in the temperate surfperch, Embiotoca jacksoni.

Acidification-induced changes in neurological function have been documented in several tropical marine fishes. Here, we investigate whether similar patterns of neurological impacts are observed in a temperate Pacific fish that naturally experiences regular and often large shifts in environmental pH/pCO2 . In two laboratory experiments, we tested the effect of acidification, as well as pH/pCO2 variability, on gene expression in the brain tissue of a common temperate kelp forest/estuarine fish, Embiotoca jacksoni. Experiment 1 employed static pH treatments (target pH = 7.85/7.30), while Experiment 2 incorporated two variable treatments that oscillated around corresponding static treatments with the same mean (target pH = 7.85/7.70) in an eight-day cycle (amplitude ± 0.15). We found that patterns of global gene expression differed across pH level treatments. Additionally, we identified differential expression of specific genes and enrichment of specific gene sets (GSEA) in comparisons of static pH treatments and in comparisons of static and variable pH treatments of the same mean pH. Importantly, we found that pH/pCO2 variability decreased the number of differentially expressed genes detected between high and low pH treatments, and that interindividual variability in gene expression was greater in variable treatments than static treatments. These results provide important confirmation of neurological impacts of acidification in a temperate fish species and, critically, that natural environmental variability may mediate the impacts of ocean acidification.

Evolutionary origin of the Atlantic Cabo Verde nibbler (Girella stuebeli), a member of a primarily Pacific Ocean family of antitropical herbivorous reef fishes.

Nibblers (family Girellidae) are reef fishes that are mostly distributed in the Indo-Pacific, with one exception: Girella stuebeli, which is found in the Cabo Verde Archipelago, in the Atlantic Ocean. We capitalized on this unusual distribution to study the evolutionary history of the girellids, and determine the relationship between G. stuebeli and the remaining nibbler taxa. Based on thousands of genomic markers (RAD sequences), we identified the closest relatives of G. stuebeli as being a clade of three species endemic to the northwestern Pacific, restricted to the Sea of Japan and vicinity. This clade diverged from G. stuebeli approximately 2.2 Mya. Two alternative potential routes of migration may explain this affinity: a western route, from the Tropical Eastern Pacific and the Tropical Western Atlantic, and an eastern route via the Indian Ocean and Southern Africa. The geological history and oceanography of the regions combined with molecular data presented here, suggest that the eastern route of invasion (via the Indian Ocean and Southern Africa) is a more likely scenario.

Tempo and mode of speciation in Holacanthus angelfishes based on RADseq markers.

In this study we estimated the timing of speciation events in a group of angelfishes using 1186 RADseq markers corresponding to 94,880 base pairs. The genus Holacanthus comprises seven species, including two clades of Panama trans-Isthmian geminates, which diverged approximately 3-3.5Mya. These clades diversified within the Tropical Eastern Pacific (TEP, three species) and Tropical Western Atlantic (TWA, two species) which our data suggest to have occurred within the past 1.5My in both ocean basins, but may have proceeded via different mechanisms. In the TEP, speciation is likely to have followed a peripatric pathway, while in the TWA, sister species are currently partially sympatric, thus raising the possibility of sympatric speciation. This study highlights the use of RADseq markers for estimating both divergence times and modes of speciation at a 1-3My timescale.

FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies.

BACKGROUND: Phylogenetic and population genetic studies often deal with multiple sequence alignments that require manipulation or processing steps such as sequence concatenation, sequence renaming, sequence translation or consensus sequence generation. In recent years phylogenetic data sets have expanded from single genes to genome wide markers comprising hundreds to thousands of loci. Processing of these large phylogenomic data sets is impracticable without using automated process pipelines. Currently no stand-alone or pipeline compatible program exists that offers a broad range of manipulation and processing steps for multiple sequence alignments in a single process run. RESULTS: Here we present FASconCAT-G, a system independent editor, which offers various processing options for multiple sequence alignments. The software provides a wide range of possibilities to edit and concatenate multiple nucleotide, amino acid, and structure sequence alignment files for phylogenetic and population genetic purposes. The main options include sequence renaming, file format conversion, sequence translation between nucleotide and amino acid states, consensus generation of specific sequence blocks, sequence concatenation, model selection of amino acid replacement with ProtTest, two types of RY coding as well as site exclusions and extraction of parsimony informative sites. Convieniently, most options can be invoked in combination and performed during a single process run. Additionally, FASconCAT-G prints useful information regarding alignment characteristics and editing processes such as base compositions of single in- and outfiles, sequence areas in a concatenated supermatrix, as well as paired stem and loop regions in secondary structure sequence strings. CONCLUSIONS: FASconCAT-G is a command-line driven Perl program that delivers computationally fast and user-friendly processing of multiple sequence alignments for phylogenetic and population genetic applications and is well suited for incorporation into analysis pipelines.

Strong population differentiation in lingcod (Ophiodon elongatus) is driven by a small portion of the genome.

Delimiting intraspecific genetic variation in harvested species is crucial to the assessment of population status for natural resource management and conservation purposes. Here, we evaluated genetic population structure in lingcod (Ophiodon elongatus), a commercially and recreationally important fishery species along the west coast of North America. We used 16,749 restriction site-associated DNA sequencing (RADseq) markers, in 611 individuals collected from across the bulk of the species range from Southeast Alaska to Baja California, Mexico. In contrast to previous population genetic work on this species, we found strong evidence for two distinct genetic clusters. These groups separated latitudinally with a break near Point Reyes off Northern California, and there was a high frequency of admixed individuals in close proximity to the break. F-statistics corroborate this genetic break between northern and southern sampling sites, although most loci are characterized by low FST values, suggesting high gene flow throughout most of the genome. Outlier analyses identified 182 loci putatively under divergent selection, most of which mapped to a single genomic region. When individuals were grouped by cluster assignment (northern, southern, and admixed), 71 loci were fixed between the northern and southern cluster, all of which were identified in the outlier scans. All individuals identified as admixed exhibited near 50:50 assignment to northern and southern clusters and were heterozygous for most fixed loci. Alignments of RADseq loci to a draft lingcod genome assembly and three other teleost genomes with chromosome-level assemblies suggest that outlier and fixed loci are concentrated on a single chromosome. Similar genomic patterns have been attributed to chromosomal inversions in diverse taxonomic groups. Regardless of the evolutionary mechanism, these results represent novel observations of genetic structure in lingcod and designate clear evolutionary units that could be used to inform fisheries management.

Molecular phylogenetics of snailfishes (Cottoidei: Liparidae) based on MtDNA and RADseq genomic analyses, with comments on selected morphological characters.

Phylogenetic relationships of snailfishes of the family Liparidae were analyzed on the basis of two sets of molecular sequence data: one from the mitochondrial DNA cytochrome c oxidase subunit one gene (COI) and another from restriction-site associated genome-wide sequences (RADseq). The analysis of COI sequence data from at least 122 species of 18 genera from the Pacific, Atlantic, and Southern oceans resulted in a moderately well-resolved phylogeny among the major clades, albeit with significant polytomy among central clades. Nectoliparis was the sister of all other members of the family, followed by Liparis. Liparis, Careproctus, and Paraliparis were paraphyletic. Liparis was recovered in two closely related clades, with L. fucensis sister of all other liparids except Nectoliparis, and both Careproctus and Paraliparis were each recovered among at least three widely separated clades. The RADseq analysis of 26 species of 11 genera from the eastern North Pacific strongly confirmed the overall results of the COI analysis, with the exception of the paraphyly of Liparis due to the absence of L. fucensis. Our results show that the pelvic disc has been independently lost multiple times and the pectoral-fin girdle has been independently reduced in multiple lineages.

Taxonomic revisions within Embiotocidae (Teleostei, Perciformes) based on molecular phylogenetics.

Embiotocidae, a unique family within the Perciformes that has evolved a complex viviparous natural history, has lacked full resolution and strong support in several interspecific relationships until recently. Here we propose three taxonomic revisions within embiotocid surfperches based on recent molecular phylogenetic analyses that robustly resolve all interspecific relationship in the Eastern Pacific species: Hypsurus caryi (Agassiz, 1853) resurrected to its original name Embiotoca caryi Agassiz, 1853, Rhacochilus vacca (Girard, 1855) shifted into the genus Phanerodon Girard, 1854, and Hyperprosopon anale Agassiz, 1861 separated into the available genus Hypocritichthys Gill, 1862. The proposed changes would leave three previously paraphyletic groups monophyletic (Embiotoca, Hyperprosopon, and Phanerodon) and would maintain the current number of genera at 13.

Altrichthys alelia, a new brooding damselfish (Teleostei, Perciformes, Pomacentridae) from Busuanga Island, Philippines.

A new species of damselfish, Altrichthys aleliasp. n. is described from specimens collected in shallow water (1-8m depth) off Busuanga Island, Palawan Province, Philippines. It differs from the other two species in the genus, A. curatus and A. azurelineatus, in various features including having golden upper body lacking dark edges of dorsal and caudal fins, higher modal number of tubed lateral line scales, as well as differences in two mitochondrial markers, one nuclear marker, and RAD markers.

The complete mitochondrial genome of the black surfperch, Embiotoca jacksoni: Selection and substitution rates among surfperches (Embiotocidae).

The complete 16,515bp nucleotide sequence of the mitochondrial genome was determined for the black surfperch, Embiotoca jacksoni (Perciformes: Embiotocidae). The black surfperch mitochondrial genome contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and the non-coding control region (D-loop), the gene order of which is identical to that observed in most vertebrates. The protein-coding gene sequences of E. jacksoni mitochondrial DNA were compared with two other embiotocid surfperches with available complete mitochondrial genomes, Cymatogaster aggregata and Ditrema temminckii. Across all mitochondrial protein-coding genes in surfperches the weighted average substitution rate was 2.079% per My and average dN/dS ratios for each protein-coding gene ranged from 0.016 in CytB to 0.608 in ND3. Substitution rates and dN/dS ratios were relatively high for ATP8 compared to other protein-coding genes. Although most protein-coding genes showed signals of purifying selection, we found evidence for positive selection in ND3 in E. jacksoni.

Phylogeography of the California sheephead, Semicossyphus pulcher: the role of deep reefs as stepping stones and pathways to antitropicality.

In the past decade, the study of dispersal of marine organisms has shifted from focusing predominantly on the larval stage to a recent interest in adult movement. Antitropical distributions provide a unique system to assess vagility and dispersal. In this study, we have focused on an antitropical wrasse genus, Semicossyphus, which includes the California sheephead, S. pulcher, and Darwin's sheephead, S. darwini. Using a phylogenetic approach based on mitochondrial and nuclear markers, and a population genetic approach based on mitochondrial control region sequences and 10 microsatellite loci, we compared the phylogenetic relationships of these two species, as well as the population genetic characteristics within S. pulcher. While S. pulcher and S. darwini are found in the temperate eastern Pacific regions of the northern and southern hemispheres, respectively, their genetic divergence was very small (estimated to have occurred between 200 and 600 kya). Within S. pulcher, genetic structuring was generally weak, especially along mainland California, but showed weak differentiation between Sea of Cortez and California, and between mainland California and Channel Islands. We highlight the congruence of weak genetic differentiation both within and between species and discuss possible causes for maintenance of high gene flow. In particular, we argue that deep and cooler water refugia are used as stepping stones to connect distant populations, resulting in low levels of genetic differentiation.