Evolutionary Toxicology and Genomics: Utility of Natural Populations’ Sensitivity and Resistance Mechanisms to Persistent Organic Pollutants in Ecological Risk Assessment and Public Health

Susceptibility to toxic chemical exposures is determined by genotype-environment interactions. Genetically variable natural populations live in spatially, temporally, and chemically complex environments, and their exposures to persistent organic pollutants should be carefully considered when quantifying evolution’s role in mechanisms of sensitivity and resistance. Although most natural populations’ genomes are not sequenced, reasonably robust lower-throughput custom platforms and data imputation methods help infer missing genotypes and improve the analysis of non-model organisms’ sequencing data. Such datasets yield insights about population processes and variance of allelic diversity within and between sensitive and resistant populations. The advantage of this approach is that an a priori choice of biomarkers of an adverse effect is not necessary. Population genomics provide insights into the sequence variations that define differences in gene expression and protein polymorphisms underlying mechanisms of sensitivity and resistance to polluted environment. Complementary “omics” approaches with model-organism and natural populations offer a better understanding of biological effects and should be an integral component of comprehensive ecological and human health risk assessment.


Introduction
The prevalence of today's environmental contaminants is unprecedented. The Toxic Substance Control Act Chemical Substance Inventory lists over 85,000 chemicals in the U.S. market [1]. Adapting to persistent organic pollutants (POPs) exposures is difficult, particularly if adaptive responses require novel genetic variation [2]. An individual's fitness is influenced by the sum of all stressors, which can be additive, antagonistic, and/or synergistic [3][4][5]. Devising effective long-term risk-assessment and mitigation strategies is etc. [6,7]. Due to their general use, POPs have high environmental concentrations; many are hydrophobic and bioaccumulate through the food chain due to their slow rate of metabolism [8][9][10]. Human and animal exposures are mostly through diet (90%), environmental exposures, and during embryo/fetal development via maternal deposition. Chronic POPs exposures can result in acute and chronic toxicity, demonstrated by developmental defects, chronic illnesses, and increased mortality [6]. Some are categorized as carcinogens [11] while many are endocrine disruptors [12][13][14] and may affect nervous, immune, digestive, urinary and reproductive systems [6].
Low-level exposure to POPs during critical fetal development stages can have a lifetime effect [12,15]. POPs concentrations in human serum increase with age and seem to be higher in females than males [16]. Evolutionary responses and adaptations to POPs have been documented since the mid-20th century; initially as pesticide resistance cases among invertebrates and plants [17,18], then among other natural populations, including vertebrates and mammals [19,20]. Although industrial toxic chemicals were provoking rapid evolutionary response among natural populations, their effects were mostly ignored.
To evaluate temporal and spatial exposure risks, government agencies study the most environmentally prevalent POPs by examining bioavailability and dose-response relationships in the laboratory setting [21]. Such assessments are limited within relatively static laboratory settings. For organisms exposed to a mixture of POPs, the adverse effects are assumed to be additive as mixtures of POPs can result in synergistic effects, enhancing the toxicity of each compound [3][4][5]. However, responses to POPs are complex and should be carefully considered when quantifying evolution's role in mechanisms of sensitivity and resistance to complex mixtures of pollutants within and between natural, genetically variable populations [22,23]. Evolutionary changes due to a chronic POP exposure can be quantified in just a few generations [24,25] and include physiological, genetic and epigenetic transgenerational effects influenced by natural selection: studies utilizing natural populations report phenotypic and molecular consequences, including increased mutation rates, receptor desensitization, phenological traits and epigenetic effects [26,27]. While some exposure effects are adaptive, others reduce population size, causing inbreeding depression or genetic drift [28,29] and decreased genetic variation, compromising the ability to adapt to future stressors. Adaptation to POPs can evolve at a cost, such as increased subsequent sensitivity to oxidative stress [30]. Risk assessment of environmental exposures to complex chemical mixtures should not be limited to often oversimplified and outdated controlled laboratory bioassays focused on robust acute and chronic toxicity endpoints such as survival and reproductive viability of model organisms lacking inherent genetic variation.
The post-industrial age of overwhelmingly complex chemical exposures coupled by recent technological advancements in NGS methodology allows for the utility of diverse non-model organisms within an evolutionary context in public health risk assessment.
Scientists started using population genomic approaches to study natural populations lacking robust genomic resources a decade ago.
Given that genome sequences of non-model organisms are accumulating at an unprecedented pace, with 234 animal and 319 contig genome assemblies currently in development [31], technological advancements in high-throughput sequencing from non-model organisms generate critically important data to infer complex genotype-environment interactions in a natural setting. Population genomic studies target the underlying variation found in the DNA among individuals and populations. They also utilize genome-wide sampling of sequence variation under the premise that demography and the evolutionary history of populations affect neutral loci similarly, whereas loci under selection will be affected differently [32].
Data of multiple, although individually analyzed, loci from natural populations enables studies of the genomic evolution, acclimatization, and adaptation to strong selective pressure such as pollution exposure in natural environment. Such data sets can yield insights about population processes and variance of allelic diversity within and between populations. The advantage is that an a priori choice of biomarkers of an adverse effect is not necessary, which is important when organisms are exposed to complex chemical mixtures in their natural environment and biomarkers of exposure are not well established [33]. Moreover, adaptation-related projects looking for genomic signatures of selection within and between natural populations can take advantage of established adaptative phenotypes to environmental gradients. Genome-wide association studies (GWAS) are challenging with non-model organisms because they lack an abundance of readily available, high quality, sequenced genomes, and genome-wide genotyping data. To improve the use of NGS data obtained from non-model organisms, new imputation methods such as Link Imputer [37] help infer missing genotypes and facilitate the analysis of non-model organisms sequencing data.
Susceptibility to toxic substances depends on the organism's genotype and interactions with the polluted environment. Natural populations inhabit spatially and temporally uncontrolled environments where pollutants are present as complex mixtures, so the challenges of "omics" technologies with natural populations are the challenges of basic biological research: which data is biologically meaningful? Considering cost, challenging sample material, lack of biological replicates, and complexity of genotype-environment interactions, choosing an optimal method is difficult. Yet, instead of a single gene, protein, or polymorphism, one can analyze thousands of genes, proteins, and polymorphisms utilizing biologically realistic, genetically variable individuals, extrapolating a population response rather than a strain-specific one. Complementary "omics" approaches offer a better understanding of biological effects at multiple levels: model-organisms are used to study many important questions in biology, population genomics provide insights into the sequence variations that govern differences in gene expression and protein polymorphisms, while transcriptomics provides insights into altered protein levels. A focused approach, using known POP inducers and inhibitors present in such mixtures and analyzing responses of cellular and tissue targets, can provide insights into mechanisms of toxicity seen in natural populations and help devise a more comprehensive risk assessment strategy. Thus, just as the integration of laboratory and field studies improves our knowledge of genes and proteins, so is the integration of laboratory and field studies critical for "omic" approaches and public health risk assessment modeling.