Signatures of Directional and Balancing Selection in the Silverside Basilichthys microlepidotus (Teleostei: Atherinopsidae) Inhabiting a Polluted River

Genome scans reveal signals of selection associated with pollution in fish populations of Basilichthys microlepidotus, an endemic species of Chile

The Maipo River catchment is one of Chile's most polluted basins. In recent decades, discharges of untreated sewage and organic matter have caused eutrophication and water quality degradation. We employed the indigenous silverfish species Basilichthys microlepidotus as a model organism to investigate the process of adaptation and selection on genes influenced by pollution. Using variation at single nucleotide polymorphisms (SNPs), we determined the temporal stability of the population structure patterns previously identified in this species by varying SNPs. We also examined local adaptation to pollution-selected genes. Using the genotypes of 7684 loci in 180 individuals, we identified 429 and 700 loci that may be undergoing selection. We detected these loci using the FSTHET and ARLEQUIN outlier detection software, respectively. Both software packages simultaneously identified a total of 250 loci. B. microlepidotus' population structure did not change over time at contaminated or unpolluted sites. In addition, our analysis found: (i) selection of genes associated with pollution, consistent with observations in other organisms; (ii) identification of candidate genes that are functionally linked to the same biological processes, molecular functions and/or cellular components that previously showed differential expression in the same populations; and (iii) a candidate gene with differential expression and a non-synonymous substitution.

Biomarker selection depends on gene function and organ: the case of the cytochrome P450 family genes in freshwater fish exposed to chronic pollution

Pollution and its effects have been of major concern in recent decades. Many strategies and markers have been developed to assess their effects on biota. Cytochrome P450 (CYP) genes have received significant attention in this context because of their relationship with detoxification and activation of exogenous compounds. While their expression has been identified as a pollution exposure biomarker, in most cases, it has been tested only after acute exposures and for CYP genes associated with exogenous compounds. To elucidate CYP gene expression patterns under chronic pollution exposure, we have used the silverside Basilichthys microlepidotus as a model, which inhabits the Maipo River Basin, a freshwater system with different pollution levels. We performed next-generation RNA sequencing of liver and gill tissues from polluted and non-polluted populations. We found most CYP genes were not dysregulated by pollution, and the seven genes that were present and differentially expressed in liver and gill were mainly downregulated. Three CYP genes associated with exogenous compounds showed differential expression in the gill, while four CYP genes associated with endogenous compounds showed differential expression in the liver. The findings presented here highlight the importance of CYP genes, his family, tissues and his interaction in the context of pollution biomarkers use.

Temporal changes in gene expression and genotype frequency of the ornithine decarboxylase gene in native silverside Basilichthys microlepidotus: Impact of wastewater reduction due to implementation of public policies

Human activity has caused a deterioration in the health and population size of riverine species; thus, public policies have been implemented to mitigate the anthropogenic impacts of water use, watercourse transformation, and pollution. We studied the Maipo River Basin, one of the most polluted with untreated wastewater in Chile, for a period of 12 years (2007-2019). Since the implementation of new public policies, including the operation of a wastewater collector (2012), the Maipo River Basin is currently much less polluted by untreated water than before. To analyze the impact of wastewater reduction in this river basin, we studied the native silverside (Basilichthys microlepidotus), which inhabits both polluted and unpolluted areas of the river basin. Previous studies reported the overexpression of the ornithine decarboxylase (odc) gene, heterozygote deficit, and high frequency of a homozygote odc genotype in silverside populations that inhabit wastewater-polluted sites, suggesting a phenotypic change and genotypic selection in response to pollution. Here, a population affected and another population unaffected by wastewater were studied before and after implementing the wastewater collector. The physicochemical data of water samples, changes in odc expression and microsatellite variability, and odc genotype frequencies were analyzed. The results showed physicochemical changes in the affected site before and after the operation of the wastewater collector. The microsatellite loci showed no changes in either population. The odc expression in the affected site was higher before the operation of the wastewater collector. Significant changes in the genotype frequencies of the odc gene before and after the wastewater collector operation were detected only at the affected site, wherein the homozygous dominant genotype decreased from >59% to <25%. Our results suggest that public policies aimed at mitigating aquatic pollution can indirectly affect both gene expression and genotype frequencies of important functional genes.

Selection of suitable reference genes for gene expression analysis in gills and liver of fish under field pollution conditions

To understand the role of gene expression in adaptive variation, it is necessary to examine expression variation in an ecological context. Quantitative real-time PCR (qPCR) is considered the most accurate and reliable technique to measure gene expression and to validate the data obtained by RNA-seq; however, accurate normalization is crucial. In Chile, the freshwater silverside fish Basilichthys microlepidotus inhabits both polluted and nonpolluted areas, showing differential gene expression related to pollution. In this study, we infer the stability of six potential reference genes (tubulin alpha, hypoxanthine-guanine phosphoribosyltransferase, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, 60S ribosomal protein L13, and 60S ribosomal protein L8) in the gills and liver of silverside individuals inhabiting polluted and nonpolluted areas. To validate the reference genes selected, the most and least stable reference genes were used to normalize two target transcripts, one for each organ. The RefFinder tool was used to analyze and identify the most stably expressed genes. The 60S ribosomal protein L8 gene was ranked as the most stable gene for both organs. Our results show that reference gene selection influences the detection of differences in the expression levels of target genes in different organs and, also highlighting candidate reference genes that could be used in field studies.

Effects of water pollution and river fragmentation on population genetic structure of invasive mosquitofish

We analyzed variation at the GPI-2 locus and eleven microsatellite loci of eastern mosquitofish Gambusia holbrooki populations introduced to the Ebro River (Spain), sampling above and below a dam (Flix Reservoir) where severe chronic pollution has been well documented. Allele frequency changes at the GPI-2 locus in the sites nearest to the polluted sediments agree with previous results from studies in mercury-exposed populations of this highly invasive fish. Genetic distinction of the mosquitofish collected close to the polluted sediments was detected at the GPI locus but also at the presumptive neutral microsatellite loci. Recent migration rates estimated from microsatellites indicated that around 30% of fish collected in a specific location were immigrants from upstream and downstream sources. Such high migration rates probably contribute to the mosquitofish's invasive success and suggest that the consequences on the mosquitofish regional genetic structured of high levels of water toxicants could be mediated by immigration from other sites, but the effect of pollutants on local diversity might be higher than observed here.

Differential gene expression revealed with RNA-Seq and parallel genotype selection of the ornithine decarboxylase gene in fish inhabiting polluted areas

How organisms adapt to unfavorable environmental conditions by means of plasticity or selection of favorable genetic variants is a central issue in evolutionary biology. In the Maipo River basin, the fish Basilichthys microlepidotus inhabits polluted and non-polluted areas. Previous studies have suggested that directional selection drives genomic divergence between these areas in 4% of Amplified Fragment Length Polymorphism (AFLP) loci, but the underlying genes and functions remain unknown. We hypothesized that B. microlepidotus in this basin has plastic and/or genetic responses to these conditions. Using RNA-Seq, we identified differentially expressed genes in individuals from two polluted sites compared with fish inhabiting non-polluted sites. In one polluted site, the main upregulated genes were related to cellular proliferation as well as suppression and progression of tumors, while biological processes and molecular functions involved in apoptotic processes were overrepresented in the upregulated genes of the second polluted site. The ornithine decarboxylase gene (related to tumor promotion and progression), which was overexpressed in both polluted sites, was sequenced, and a parallel pattern of a heterozygote deficiency and increase of the same homozygote genotype in both polluted sites compared with fish inhabiting the non-polluted sites was detected. These results suggest the occurrence of both a plastic response in gene expression and an interplay between phenotypic change and genotypic selection in the face of anthropogenic pollution.