Sex-Biased Protein Expression in Threespine Stickleback, Gasterosteus aculeatus

Reporting and reproducibility: Proteomics of fish models in environmental toxicology and ecotoxicology

Environmental toxicology and ecotoxicology research efforts are employing proteomics with fish models as New Approach Methodologies, along with in silico, in vitro and other omics techniques to elucidate hazards of toxicants and toxins. We performed a critical review of toxicology studies with fish models using proteomics and reported fundamental parameters across experimental design, sample preparation, mass spectrometry, and bioinformatics of fish, which represent alternative vertebrate models in environmental toxicology, and routinely studied animals in ecotoxicology. We observed inconsistencies in reporting and methodologies among experimental designs, sample preparations, data acquisitions and bioinformatics, which can affect reproducibility of experimental results. We identified a distinct need to develop reporting guidelines for proteomics use in environmental toxicology and ecotoxicology, increased QA/QC throughout studies, and method optimization with an emphasis on reducing inconsistencies among studies. Several recommendations are offered as logical steps to advance development and application of this emerging research area to understand chemical hazards to public health and the environment.

Proteomics in systems toxicology

Proteins are the ultimate product of gene expression. As they hinge between gene transcription and phenotype, they offer a more realistic perspective of toxicopathic effects, responses and even susceptibility to insult than targeting genes and mRNAs while dodging some inter-individual variability that hinders measuring downstream endpoints like metabolites or enzyme activity. Toxicologists have long focused on proteins as biomarkers but the advent of proteomics shifted risk assessment from narrow single-endpoint analyses to whole-proteome screening, enabling deriving protein-centric adverse outcome pathways (AOPs), which are pivotal for the derivation of Systems Biology informally named Systems Toxicology. Especially if coupled pathology, the identification of molecular initiating events (MIEs) and AOPs allow predictive modeling of toxicological pathways, which now stands as the frontier for the next generation of toxicologists. Advances in mass spectrometry, bioinformatics, protein databases and top-down proteomics create new opportunities for mechanistic and effects-oriented research in all fields, from ecotoxicology to pharmacotoxicology.

Development of a Gill Assay Library for Ecological Proteomics of Threespine Sticklebacks (Gasterosteus aculeatus)

A data-independent acquisition (DIA) assay library for quantitative analyses of proteome dynamics has been developed for gills of threespine sticklebacks (Gasterosteus aculeatus). A raw spectral library was generated by data-dependent acquisition (DDA) and annotation of tryptic peptides to MSMS spectra and protein database identifiers. The assay library was constructed from the raw spectral library by removal of low-quality, ambiguous, and low-signal peptides. Only unique proteins represented by at least two peptides are included in the assay library, which consists of 1506 proteins, 5074 peptides, 5104 precursors, and 25,322 transitions. This assay library was used with DIA data to identify biochemical differences in gill proteomes of four populations representing different eco- and morpho-types of threespine sticklebacks. The assay library revealed unique and reproducible proteome signatures. Warm-adapted, low-plated, brackish-water fish from Laguna de la Bocana del Rosario (Mexico) show elevated HSP47, extracellular matrix, and innate immunity proteins whereas several immunoglobulins, interferon-induced proteins, ubiquitins, proteolytic enzymes, and nucleic acid remodeling proteins are reduced. Fully-plated, brackish-water fish from Westchester Lagoon (Alaska) display elevated ion regulation, GTPase signaling, and contractile cytoskeleton proteins, altered abundances of many ribosomal, calcium signaling and immunity proteins, and depleted transcriptional regulators and metabolic enzymes. Low-plated freshwater fish from Lake Solano (California) have elevated inflammasomes and proteolytic proteins whereas several iron containing and ion regulatory proteins are reduced. Gills of fully-plated, marine fish from Bodega Harbor (California) have elevated oxidative metabolism enzymes and reduced transglutaminase 2, collagens, and clathrin heavy chains. These distinct proteome signatures represent targets for testing ecological and evolutionary influences on molecular mechanisms of gill function in threespine sticklebacks. Furthermore, the gill assay library represents a model for other tissues and paves the way for accurate and reproducible network analyses of environmental context-dependent proteome dynamics in complex organisms.

Reprint of: Environmental toxicology and omics: A question of sex

Molecular initiating events and downstream transcriptional/proteomic responses provide valuable information for adverse outcome pathways, which can be used predict the effects of chemicals on physiological systems. There has been a paucity of research that addresses sex-specific expression profiling in toxicology and due to cost, time, and logistic considerations, sex as a variable has not been widely considered. In response to this deficiency, federal agencies in the United States, Canada, and Europe have highlighted the importance of including sex as a variable in scientific investigations. Using case studies from both aquatic and mammalian toxicology, we report that there can be less than ~20-25% consensus in how the transcriptome and proteome of each sex responds to chemicals. Chemicals that have been shown to elicit sex-specific responses in the transcriptome or proteome include pharmaceuticals, anti-fouling agents, anticorrosive agents, and fungicides, among others. Sex-specific responses in the transcriptome and proteome are not isolated to whole animals, as investigations demonstrate that primary cell cultures isolated from each sex responds differently to toxicants. This signifies that sex is important, even in cell lines. Sex has significant implications for predictive toxicology, and both male and female data are required to improve robustness of adverse outcome pathways.

BIOLOGICAL SIGNIFICANCE

Clinical toxicology recognizes that sex is an important variable, as pharmacokinetics (ADME; absorption, distribution, metabolism, and excretion) can differ between females and males. However, few studies in toxicology have explored the implication of sex in relation to the transcriptome and proteome of whole organisms. High-throughput molecular approaches are becoming more frequently applied in toxicity screens (e.g. pre-clinical experiments, fish embryos, cell lines, synthetic tissues) and such data are expected to build upon reporter-based cell assays (e.g. receptor activation, enzyme inhibition) used in toxicant screening programs (i.e. Tox21, ToxCast, REACH). Thus, computational models can more accurately predict the diversity of adverse effects that can occur from chemical exposure within the biological system. Our studies and those synthesized from the literature suggest that the transcriptome and proteome of females and males respond quite differentially to chemicals. This has significant implications for predicting adverse effects in one sex when using molecular data generated in the other sex. While molecular initiating events are not expected to differ dramatically between females and males (i.e. an estrogen binds estrogen receptors in both sexes), it is important to acknowledge that the downstream transcriptomic and proteomic responses can differ based upon the presence/absence of co-regulators and inherent sex-specific variability in regulation of transcriptional and translational machinery. Transcriptomic and proteomic studies also reveal that cell processes affected by chemicals can differ due to sex, and this can undoubtedly lead to sex-specific physiological responses.

Comparison of Integrated Responses to Nonlethal and Lethal Hypothermal Stress in Milkfish (Chanos chanos): A Proteomics Study

Milkfish is an important aquaculture species in Taiwan, and its high mortality during cold snaps in winter usually causes huge economic losses. To understand the effect of hypothermal stress and the corresponding compensatory stress response in milkfish, this study aimed to compare liver and gill protein levels between milkfish exposed to nonlethal (18°C), lethal (16°C), and control (28°C) temperatures. Using a proteomics approach based on two-dimensional electrophoresis and nano-LC-MS/MS analysis, this study identified thirty unique protein spots from milkfish livers and gills for which protein abundance was significantly different between nonlethal, lethal, and control temperature groups. Proteins identified in the liver were classified into three different categories according to their cellular function: (1) anti-oxidative stress, (2) apoptotic pathway, and (3) cytoskeleton. Similarly, proteins identified in the gill were sorted in five different functional categories: (1) cytoskeleton, (2) immune response, (3) protein quality control, (4) energy production, and (5) intracellular homeostasis. Based on functional information derived from the identified proteins, we assumed that different levels of hypothermal stress had a different effect and induced a different cellular response. Upon nonlethal hypothermal stress, the identified proteins were involved in anti-oxidative stress and anti-inflammation pathways, suggesting that milkfish had high levels of oxidative stress in the liver and exhibited inflammation response in the gill. Upon lethal hypothermal stress, however, identified proteins were associated with apoptosis in the liver and regulation of intracellular homeostasis in the gill. The present study provided evidence to illustrate different multi-physiological responses to nonlethal and lethal hypothermal stress in milkfish livers and gills.

Deep sexual dimorphism in adult medaka fish liver highlighted by multi-omic approach

Sexual dimorphism describes the features that discriminate between the two sexes at various biological levels. Especially, during the reproductive phase, the liver is one of the most sexually dimorphic organs, because of different metabolic demands between the two sexes. The liver is a key organ that plays fundamental roles in various physiological processes, including digestion, energetic metabolism, xenobiotic detoxification, biosynthesis of serum proteins, and also in endocrine or immune response. The sex-dimorphism of the liver is particularly obvious in oviparous animals, as the female liver is the main organ for the synthesis of oocyte constituents. In this work, we are interested in identifying molecular sexual dimorphism in the liver of adult medaka fish and their sex-variation in response to hepatotoxic exposures. By developing an integrative approach combining histology and different high-throughput omic investigations (metabolomics, proteomics and transcriptomics), we were able to globally depict the strong sexual dimorphism that concerns various cellular and molecular processes of hepatocytes comprising protein synthesis, amino acid, lipid and polysaccharide metabolism, along with steroidogenesis and detoxification. The results of this work imply noticeable repercussions on the biology of oviparous organisms environmentally exposed to chemical or toxin issues.

Gender-Specific Toxicological Effects of Chronic Exposure to Pure Microcystin-LR or Complex Microcystis aeruginosa Extracts on Adult Medaka Fish

Cyanobacterial blooms often occur in freshwater lakes and constitute a potential health risk to human populations, as well as to other organisms. However, their overall and specific implications for the health of aquatic organisms that are chronically and environmentally exposed to cyanobacteria producing hepatotoxins, such as microcystins (MCs), together with other bioactive compounds have still not been clearly established and remain difficult to assess. The medaka fish was chosen as the experimental aquatic model for studying the cellular and molecular toxicological effects on the liver after chronic exposures (28 days) to environmentally relevant concentrations of pure MC-LR, complex extracts of MC producing or nonproducing cyanobacterial biomasses, and of a Microcystis aeruginosa natural bloom. Our results showed a higher susceptibility of females to the different treatments compared to males at both the cellular and the molecular levels. Although hepatocyte lysis increased with MC-containing treatments, lysis always appeared more severe in the liver of females compare to males, and the glycogen cellular reserves also appeared to decrease more in the liver of females compared to those in the males. Proteomic investigations reveal divergent responses between males and females exposed to all treatments, especially for proteins involved in metabolic and homeostasis processes. Our observations also highlighted the dysregulation of proteins involved in oogenesis in female livers. These results suggest that fish populations exposed to cyanobacteria blooms may potentially face several ecotoxicological issues.

Subtle tissue and sex-dependent proteome variation in mussel (Mytilus galloprovincialis) populations of the Galician coast (NW Spain) raised in a common environment

The mussel Mytilus galloprovincialis is one of the most important marine resources for aquaculture in Europe, and Galicia (NW Spain) is the EU's leading region for production. Variation in environmental and ecological factors exists in Northern and Southern estuaries of this region, and natural selection could have modulated genetic variation among populations with adaptation to local conditions as the driving force. Results from a previous genetic study using neutral markers suggested subtle genetic differentiation between mussel populations from both estuarine areas. In this new study, mussel samples from Northern and Southern estuaries were brought into a common environment to test for proteome differences due to genetic and permanent non-genetic effects in populations from both estuarine areas, using both foot and mantle border tissues. Because the sex of the mussels was determined through histological tests, sex-specific effects were also examined. Evidence of subtle differences in the foot proteome, dependent on mussel sex, were detected between populations from both estuaries. These differences were more marked for female samples. No evidence of proteome differences was found for the factors estuaries and sex in mantle border tissue. Candidate proteins with a potential role in local adaptation were identified and point to molecular functions that might be involved in responses to different stressors.

Transcriptomic dissection of sexual differences in Bemisia tabaci, an invasive agricultural pest worldwide

Sex difference involving chromosomes and gene expression has been extensively documented. In this study, the gender difference in the sweetpotato whitefly Bemisia tabaci was investigated using Illumina-based transcriptomic analysis. Gender-based RNAseq data produced 27 Gb reads, and subsequent de novo assembly generated 93,948 transcripts with a N50 of 1,853 bp. A total of 1,351 differentially expressed genes were identified between male and female B. tabaci, and majority of them were female-biased. Pathway and GO enrichment experiments exhibited a gender-specific expression, including enriched translation in females, and enhanced structural constituent of cuticle in male whiteflies. In addition, a putative transformer2 gene (tra2) was cloned, and the structural feature and expression profile of tra2 were investigated. Sexually dimorphic transcriptome is an uncharted territory for the agricultural insect pests. Molecular understanding of sex determination in B. tabaci, an emerging invasive insect pest worldwide, will provide potential molecular target(s) for genetic pest control alternatives.

Heterogeneous Genomic Differentiation in marine threespine sticklebacks: adaptation along an environmental gradient

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on-going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F(ST) ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F(ST) ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers--45% of which were identified as outliers--correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.