Climate Adaptation and Drift Shape the Genomes of Two Eel-Goby Sister Species Endemic to Contrasting Latitude

Evidence for Selection on Mitochondrial OXPHOS Genes in the Mediterranean Killifish Aphanius fasciatus Valenciennes, 1821

Mitochondrial oxidative phosphorylation (OXPHOS) genes are a system subject to selection under determined environmental constraints despite a neutral evolution model that has long been hypothesized for the mitochondrial genome. In this study, the sequences of ND1, Cytb, and COI OXPHOS genes were analyzed in six populations of the eurythermal and euryhaline killifish A. fasciatus, to detect non-synonymous mutations leading to amino acid changes and to check whether selection acted on them using tests of recombination and selection. The results indicate a high COI and Cytb gene diversity and a high percentage of private haplotypes in all populations. In the Greek population, non-synonymous nucleotide substitutions were observed in the N-terminal region of COI and Cytb. Positively selected sites were also found. The information we obtained from the mitochondrial DNA sequences of A. fasciatus adds to the growing data on selective pressure acting on mitochondrial DNA in non-model species. These results should be explored from the perspective of the local adaptation of eurythermal and euryhaline species and supported using experimental evidence to better understand the interplay between historical climatic events and local adaptation and how each of them contributes to shaping the genetic structure of this species.

Chromosome-level genome assembly and annotation of eel goby (Odontamblyopus rebecca)

The eel gobies fascinate researchers with many important features, including its unique body structure, benthic lifestyle, and degenerated eyes. However, genome assembly and exploration of the unique genomic composition of the eel gobies are still in their infancy. This has severely limited research progress on gobies. In this study, multi-platform sequencing data were generated and used to assemble and annotate the genome of O. rebecca at the chromosome-level. The assembled genome size of O. rebecca is 918.57 Mbp, which is similar to the estimated genome size (903.03 Mbp) using 17-mer. The scaffold N50 is 41.67 Mbp, and 23 chromosomes were assembled using Hi-C technology with a mounting rate of 99.96%. Genome annotation indicates that 53.29% of the genome is repetitive sequences, and 22,999 protein-coding genes are predicted, of which 21,855 have functional annotations. The chromosome-level genome of O. rebecca will not only provide important genomic resources for comparative genomic studies of gobies, but also expand our knowledge of the genetic origin of their unique features fascinating researchers for decades.

Comparative Transcriptome Analyses Provide New Insights into the Evolution of Divergent Thermal Resistance in Two Eel Gobies

Adaptation to thermal conditions in tidal mudflats always involves tolerating frequent fluctuations and often extreme environmental temperatures. Regulation of gene expression plays a fundamental role in the evolution of these thermal adaptations. To identify the key gene regulatory networks associated with the thermal adaptation, we investigated the capability of cold tolerance, as well as the transcriptomic changes under cold stress in two mudflat inhabitants (Odontamblyopus lacepedii and O. rebecca) with contrasting latitude affinity. Our results revealed a remarkable divergent capacity of cold tolerance (CTmin: 0.61 °C vs. 9.57 °C) between the two gobies. Analysis of transcriptomic changes under cold stress unveiled 193 differentially expressed genes exhibiting similar expression profiles across all tissues and species, including several classic metabolic and circadian rhythm molecules such as ACOD and CIART that may represent the core cold response machinery in eel gobies. Meanwhile, some genes show a unique expression spectrum in the more cold-tolerant O. lacepedii suggesting their roles in the enhanced cold tolerance and hence the extreme thermal adaptations. In addition, a weighted gene co-expression network analysis (WGCNA) revealed a subset of metabolic hub genes including MYH11 and LIPT2 showing distinct down-regulation in O. lacepedii when exposed to cold stress which highlights the role of reduced energy consumption in the enhanced cold tolerance of eel gobies. These findings not only provide new insights into how mudflat teleosts could cope with cold stress and their potential evolutionary strategies for adapting to their thermal environment, but also have important implications for sound management and conservation of their fishery resources in a scenario of global climate warming in the marine realm.