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.

Population genomic and biophysical modeling show different patterns of population connectivity in the spiny lobster Jasus frontalis inhabiting oceanic islands.

Knowledge about connectivity between populations is essential for the fisheries management of commercial species. The lobster Jasus frontalis inhabits two oceanic island groups, the Juan Fernández Archipelago and the Desventuradas Islands, separated by 800 km. Since this species is primarily exploited in the Juan Fernández Archipelago, knowledge of the connectivity patterns among islands is foundational for species management. Here, we used variability at single-nucleotide polymorphisms (SNPs) and individual-based modeling (IBM) to estimate the genetic structure and connectivity between J. frontalis populations in these island groups. The variability at 9090 SNPs suggests two genetic populations, one in the Juan Fernández Archipelago and one in the Desventuradas Islands. Furthermore, IBM suggests an asymmetric connectivity pattern, with particles moving from the Juan Fernández Archipelago to the Desventuradas Islands but not vice versa. Since the IBM analysis suggests asymmetric larval movement between the islands, and the genetic analysis indicates isolation between the Juan Fernández Archipelago and the Desventuradas Islands, larval retention mechanisms such as small-scale oceanographic processes or behavior could hinder larval movement between islands. This study highlights the importance of using more than one methodology to estimate population connectivity.

Spatial and temporal stability in the genetic structure of a marine crab despite a biogeographic break.

Elucidating the processes responsible for maintaining the population connectivity of marine benthic species mediated by larval dispersal remains a fundamental question in marine ecology and fishery management. Understanding these processes becomes particularly important in areas with a biogeographic break and unidirectional water movement along the sides of the break. Based on variability at 4209 single-nucleotide polymorphisms in 234 individuals, we determine the genetic structure, temporal genetic stability, and gene flow among populations of the commercially important mola rock crab Metacarcinus edwardsii in a system in southern Chile with a biogeographic break at latitude 42°S. Specimens were collected at eight sites within its geographic distribution, with collection at four of these sites was performed twice. Using population genetic approaches, we found no evidence of geographic or temporal population differentiation. Similarly, we found no evidence of an effect on gene flow of the biogeographic break caused by the the West Wind Drift Current. Moreover, migration analyses supported gene flow among all sites but at different rates for different pairs of sites. Overall, our findings indicate that M. edwardsii comprises a single large population with high levels of gene flow among sites separated by over 1700 km and demonstrate temporal stability in its genetic structure.

Geographical context outweighs habitat disturbance effects in explaining mistletoe population genetic differentiation at a regional scale.

Genetic differentiation depends on ecological and evolutionary processes that operate at different spatial and temporal scales. While the geographical context is likely to determine large-scale genetic variation patterns, habitat disturbance events will probably influence small-scale genetic diversity and gene flow patterns. Therefore, the genetic diversity patterns that we observe today result from the combination of both processes, but they are rarely assessed simultaneously. We determined the population structure and genetic diversity of a hemiparasitic mistletoe (Tristerix corymbosus) from the temperate rainforests of southern Chile to determine the effects of geographical context and habitat disturbance at a regional scale and if it is affected by the abundance and occurrence of its seed disperser mutualist (the arboreal marsupial Dromiciops gliroides). We genotyped 359 individuals from 12 populations using single nucleotide polymorphisms, across three different geographical contexts and four disturbance conditions. We also used camera traps to estimate the abundance and occurrence of the seed disperser. Our results suggest that genetic differences among populations are related more to geographical context than to habitat disturbance. However, as disturbance increased, D. gliroides abundance and occurrence decreased, and mistletoe inbreeding index (FIS ) increased. We also found highly uneven gene flow among study sites. Despite the high levels of disturbance that these temperate rainforests are facing, our results suggest that mistletoe genetic differentiation at a regional scale was more influenced by historical events. However, habitat disturbance can indirectly affect mistletoe population genetic differentiation via the seed dispersal process, which may increase levels of inbreeding.

High levels of connectivity over large distances in the diadematid sea urchin Centrostephanus sylviae.

Most benthic marine invertebrates with sedentary benthic adult phases have planktonic larvae that permit connectivity between geographically isolated populations. Planktonic larval duration and oceanographic processes are vital to connecting populations of species inhabiting remote and distant islands. In the present study, we analyzed the population genetic structure of the sea urchin Centrostephanus sylviae, which inhabits only the Juan Fernández Archipelago and the Desventuradas islands, separated by more than 800 km. For 92 individuals collected from Robinson Crusoe and Selkirk Islands (Juan Fernández Archipelago) and San Ambrosio Island (Desventuradas Islands), 7,067 single nucleotide polymorphisms (SNPs) were obtained. The results did not show a spatial genetic structure for C. sylviae; relative high migration rates were revealed between the islands. An analysis of the water circulation pattern in the area described a predominant northward water flow with periods of inverted flow, suggesting that larvae could move in both directions. Overall, this evidence suggests that C. sylviae comprises a single large population composed of individuals separated by more than 800 km.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle.

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

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.

Highly Connected Populations and Temporal Stability in Allelic Frequencies of a Harvested Crab from the Southern Pacific Coast.

For marine invertebrates with a benthic adult form and a planktonic larva phase, the connectivity among populations is mainly based on larval dispersal. While an extended larval phase will promote gene flow, other factors such as an intensive fishery and geographical barriers could lead to changes in genetic variability. In this study, the population genetic structure of the commercial crab Metacarcinus edwardsii was analyzed along 700 km of the Chilean coast. The analysis, based on eight microsatellite loci genotyped from megalopae and adult crabs, considered temporal and spatial patterns of genetic variation. The results showed no evidence of spatial patterns in genetic structure, suggesting high connectivity among the sampling sites. The temporal analysis showed no evidence of changes in allele frequencies and no evidence of a recent bottleneck. The lack of spatial structure and allele variation over time could be explained by the interaction of factors such as i) low reproductive variance due to the capability of females to store sperm in the seminal receptacle, which can be used for successive broods, ii) high larval dispersal and iii) high individual reproductive output. Using our data as priors, a genetic modelling approach coincided, predicting this temporal and spatial stability. The same analysis showed that a reduction in population size leads to the loss of genetic variability in populations, as well as of the genetic cohesiveness between populations, pointing out the importance management for species under exploitation, such as M. edwardsii.