Transcriptomic differentiation underlying marine-to-freshwater transitions in the South American silversides Odontesthes argentinensis and O. bonariensis (Atheriniformes)

Differential expression analyses and detection of SNP loci associated with environmental variables: Are salinity and temperature factors involved in population differentiation and speciation in Odontesthes?

Environmental factors play a key role in individual adaptation to different local conditions. Because of this, studies about the physiological and genetic responses of individuals exposed to different natural environments offer clues about mechanisms involved in population differentiation, and as a subsequent result, speciation. Marine environments are especially suited to survey this kind of phenomena because they commonly harbor species adapted to different local conditions along a geographic continuum. Silversides belonging to Odontesthes are commonly distributed in tropical and temperate regions of South America and exhibit noticeable phenotypic plasticity, which allows them to adapt to contrasting environments. In this study, the genetic expression of O. argentinensis sampled along the Uruguayan Atlantic coast and estuarine adjacent areas was investigated. In addition, the correlation between individual genotypes and environmental variables was also analysed in O. argentinensis and O. bonariensis. Results obtained suggest a differential expression pattern of low magnitude among individuals from the different areas sampled and a correlation between several SNP loci and environmental variables. The analyses carried out did not show a clear differentiation among individuals sampled along different salinity regimens, but enriched GOTerms seem to be driven by water oxygen content. On the other hand, a total of 46 SNPs analysed in O. argentinensis and O. bonariensis showed a correlation with salinity and temperature. Although none of the correlated SNPs and corresponding genes from our both analyses were directly associated with hypoxia, genes related to the cardiovascular system and muscle cell differentiation were found. All these genes are interesting candidates for future studies since they are closely related to the differentially expressed genes. Although salinity was also mentioned as an important parameter limiting introgression between O. argentinensis and O. bonariensis, it was found that salinity does not drive differential expression in O. argentinensis, but rather oxygen levels. Moreover, salinity does not directly affect the structure and genetic divergence of the populations, they appear to be structured based on their degree of isolation and geographical distance between them. Further studies, like genome-wide analyses, could help to elucidate additional genes adapted to the different environments in these silverside species.

Hybridization is strongly constrained by salinity during secondary contact between silverside fishes (Odontesthes, Atheriniformes)

This study investigates a contact zone between two silverside fish species (marine Odontesthes argentinensis and freshwater O. bonariensis) in the estuarine Mar Chiquita lagoon along the Atlantic coast in Argentina (MChL), in which intermediate morphs had been reported. It has been suggested that admixture and introgression occur in MChL between these two species, but direct genetic evidence is lacking. Leveraging samples collected over several years (n = 676), we document the spatial distribution of both species and intermediate morphs within this habitat and collect landmark-based morphometric and multilocus genetic data (9876 loci for n = 110 individuals) to test the hypothesis of hybridization. Our analysis unambiguously characterizes intermediate morphs as F1 or F2 hybrids. We show that the low frequency of hybrid individuals in MChL may be explained by uneven abundance of parental species, which in turn are strongly affected by water salinity, limiting the size of the contact zone. Although hybrids seem to be fertile, their fitness may be reduced by external and intrinsic factors that may limit their success and suggest that this is an unstable hybrid zone. Genetic distinctiveness of both parental species is strongly supported by genome-wide data, explaining a known pattern of mitonuclear discordance as a consequence of hybridization followed by mitochondrial introgression. A clear signature of population genetic structure was detected in O. argentinensis, distinguishing MChL residents from marine populations of this species, that also was supported by distinctive morphometric features among these groups. Previous hypotheses of speciation in these fishes are discussed in the light of the new findings.

Comparative mitogenomics of Clupeoid fish provides insights into the adaptive evolution of mitochondrial oxidative phosphorylation (OXPHOS) genes and codon usage in the heterogeneous habitats

Clupeoid fish can be considered excellent candidates to understand the role of mitochondrial DNA in adaptive evolution, as they have colonized different habitats (marine, brackish, freshwater, tropical and temperate regions) over millions of years. Here, we investigate patterns of tRNA location, codon usage bias, and lineage-specific diversifying selection signals to provide novel insights into how evolutionary improvements of mitochondrial metabolic efficiency have allowed clupeids to adapt to different habitats. Based on whole mitogenome data of 70 Clupeoids with a global distribution we find that purifying selection was the dominant force acting and that the mutational deamination pressure in mtDNA was stronger than the codon/amino acid constraints. The codon usage pattern appears evolved to achieve high translational efficiency (codon/amino acid-related constraints), as indicated by the complementarity of most codons to the GT-saturated tRNA anticodon sites (retained by deamination-induced pressure) and usage of the codons of the tRNA genes situated near to the control region (fixed by deamination pressure) where transcription efficiency was high. The observed shift in codon preference patterns between marine and euryhaline/freshwater Clupeoids indicates possible selection for improved translational efficiency in mitochondrial genes while adapting to low-salinity habitats. This mitogenomic plasticity and enhanced efficiency of the metabolic machinery may have contributed to the evolutionary success and abundance of Clupeoid fish.

Comprehensive analysis of transcriptomics and metabolomics to understand the flesh quality regulation of crucian carp (Carassius auratus) treated with short term micro-flowing water system

The short term micro-flowing purification system (STMFPS) has been shown to improve the flesh quality of freshwater fish. However, few studies have focused on the involved underlying mechanisms. This study explored the effect of STMFPS on the flesh quality of market-size freshwater fish based on the combination of metabolomics and transcriptomics methods. The UPLC-QTOF/MS based metabolomics method was utilized to screen metabolites and predict the possible major metabolic pathways during different STMFPS treatment periods (0 d, 1 d, 5 d and 9 d). Furthermore, the transcriptomic data demonstrated that the differentially expressed genes detected in crucian carp muscle were 2915, 7852 and 7183 after 1 d, 5 d and 9 d STMFPS treatment. Results showed that the TCA cycle, ornithine cycle, purine metabolism and amino acid catabolism play important roles in improving the flesh quality of crucian carp. This study may help to understand the mechanism of improving the flesh quality of aquatic products using STMFPS.

Mitochondrial genome architecture and phylogenetic relationships of Odontesthes argentinensis within Atherinomorpha

Silversides are a widely distributed group across South America, with several species occupying marine, freshwater and estuarine environments. Several authors suggest main transitions among these environments took place during Pleistocene, and were accompanied with rapid speciation events. This scenario produced very limited genetic and morphological differentiation among the species. However, most of these surveys have an incomplete coverage of the intraspecific genetic diversity of the taxa studied. In this work, we reconstructed six mitochondrial genomes of O. argentinensis using transcriptomic data, and used them-in combination with several nuclear markers retrieved from the same transcriptomes-to explore the effect of additional coverage of intraspecific diversity of this species in phylogenetic reconstructions. Unlike previous works, phylogenetic analyses failed to identify O. argentinensis as a monophyletic group in relation with closely related taxa. Our results suggest that several species of the genus, especially those related to O. argentinensis, need further taxonomic revision.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Inferring boundaries among fish species of the new world silversides (Atherinopsidae; genus Odontesthes): new evidences of incipient speciation between marine and brackish populations of Odontesthes argentinensis

Species of new world silversides (Actinopterygii; Atherinopsidae; genus Odontesthes) possess economic relevance, biological interest and ecological importance. In the present paper we: (A) investigate the molecular diversity in marine species of Odontesthes from the South West Atlantic Ocean (SWAO), and analyse their interspecific relationships and divergence by means of DNA Barcoding, including its freshwater congeners, as well. (B) Explore the suitability of DNA Barcoding to analyse the diversity and distribution of haplotypes in Odontesthes argentinensis, the only well documented marine species from the SWAO that exhibit putative estuarine and marine populations. Molecular analysis revealed 100% of agreement between morphological identification and molecular identity. Odontesthes argentinensis, Odontesthes platensis, Odontesthes smitti, Odontesthes nigricans and Odontesthes incisa were assigned to five different barcode index numbers (BINs). Maximum-likelihood analysis showed that all marine species of Odontesthes clustered separately in a unique monophyletic phylogroup, comprising five well defined haplogroups, with genetic divergence between groups ranging from 2.75 to 7.11%. The genetic analysis including freshwater congeners showed that O. incisa clustered alone occupying a basal position. The F_st pairwise comparisons within O. argentinensis support the existence of three population groups: one conformed by Mar Chiquita Lagoon (MCh) specimens, and the others by Mar del Plata/Mar Chiquita coast and San Blas Bay coastal specimens, respectively. The AMOVA showed significant overall differentiation (F_st = 0.238; p = 0.00001) for the entire data set. The previous/present evidence is discussed, and strongly suggests that incipient speciation is occurring in O. argentinensis argentinean populations, and specimens from MCh would be considered at present as the leading candidate of a marine to freshwater incipient speciation event.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Evaluation of Reference Genes to Analyze Gene Expression in Silverside Odontesthes humensis Under Different Environmental Conditions

Some mammalian reference genes, which are widely used to normalize the qRT-PCR, could not be used for this purpose due to its high expression variation. The normalization with false reference genes leads to misinterpretation of results. The silversides (Odontesthes spp.) has been used as models for evolutionary, osmoregulatory and environmental pollution studies but, up to now, there are no studies about reference genes in any Odontesthes species. Furthermore, many studies on silversides have used reference genes without previous validations. Thus, present study aimed to was to clone and sequence potential reference genes, thereby identifying the best ones in Odontesthes humensis considering different tissues, ages and conditions. For this purpose, animals belonging to three ages (adults, juveniles, and immature) were exposed to control, Roundup®, and seawater treatments for 24 h. Blood samples were subjected to flow-cytometry and other collected tissues to RNA extraction; cDNA synthesis; molecular cloning; DNA sequencing; and qRT-PCR. The candidate genes tested included 18s, actb, ef1a, eif3g, gapdh, h3a, atp1a, and tuba. Gene expression results were analyzed using five algorithms that ranked the candidate genes. The flow-cytometry data showed that the environmental challenges could trigger a systemic response in the treated fish. Even during this systemic physiological disorder, the consensus analysis of gene expression revealed h3a to be the most stable gene expression when only the treatments were considered. On the other hand, tuba was the least stable gene in the control and gapdh was the least stable in both Roundup® and seawater groups. In conclusion, the consensus analyses of different tissues, ages, and treatments groups revealed that h3a is the most stable gene whereas gapdh and tuba are the least stable genes, even being considered two constitutive genes.

Transcriptomic differentiation underlying marine-to-freshwater transitions in the South American silversides Odontesthes argentinensis and O. bonariensis (Atheriniformes)

Salinity gradients are critical habitat determinants for freshwater organisms. Silverside fishes in the genus Odontesthes have recently and repeatedly transitioned from marine to freshwater habitats, overcoming a strong ecological barrier. Genomic and transcriptomic changes involved in this kind of transition are only known for a few model species. We present new data and analyses of gene expression and microbiome composition in the gills of two closely related silverside species, marine O. argentinensis and freshwater O. bonariensis and find more than three thousand transcripts differentially expressed, with osmoregulatory/ion transport genes and immune genes showing very different expression patterns across species. Interspecific differences also involve more than one thousand transcripts with nonsynonymous SNPs in the coding sequences, most of which were not differentially expressed. In addition to characterizing gill transcriptomes from wild‐caught marine and freshwater fishes, we test experimentally the response to salinity increases by O. bonariensis collected from freshwater habitats. Patterns of expression in gill transcriptomes of O. bonariensis exposed to high salinity do not resemble O. argentinensis mRNA expression, suggesting lack of plasticity for adaptation to marine conditions in this species. The diversity of functions associated with both the differentially expressed set of transcripts and those with sequence divergence plus marked microbiome differences suggest that multiple abiotic and biotic factors in marine and freshwater habitats are driving transcriptomic differences between these species.