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

Effects of pejus and pessimum zone salinity stress on gill proteome networks and energy homeostasis in Oreochromis mossambicus

Salinity tolerance in fish involves a suite of physiological changes, but a cohesive theory leading to a mechanistic understanding at the organismal level is lacking. To examine the potential of adapting energy homeostasis theory in the context of salinity stress in teleost fish, Oreochromis mossambicus were acclimated to hypersalinity at multiple rates and durations to determine salinity ranges of tolerance and resistance. Over 3000 proteins were quantified simultaneously to analyze molecular phenotypes associated with hypersalinity. A species- and tissue-specific data-independent acquisition (DIA) assay library of MSMS spectra was created. Protein networks representing complex molecular phenotypes associated with salinity acclimation were generated. O. mossambicus has a wide "zone of resistance" from 75 g/kg salinity to 120 g/kg. Crossing into the zone of resistance resulted in marked phenotypic changes including blood osmolality over 400 mOsm/kg, reduced body condition, and cessation of feeding. Protein networks impacted by hypersalinity consist of electron transport chain (ETC) proteins and specific osmoregulatory proteins. Cytoskeletal, cell adhesion, and extracellular matrix proteins are enriched in networks that are sensitive to the critical salinity threshold. These network analyses identify specific proteome changes that are associated with distinct zones described by energy homeostasis theory and distinguish them from general hypersalinity-induced proteome changes.

Tilapia Fish Skin Treatment of Third-Degree Skin Burns in Murine Model

This study explored the feasibility of using fish skin bandages as a therapeutic option for third-degree skin burns. Following the California wildfires, clinical observations of animals with third-degree skin burns demonstrated increased comfort levels and reduced pain when treated with tilapia fish skin. Despite the promises of this therapy, there are few studies explaining the healing mechanisms behind the application of tilapia fish skin. In this study, mice with third-degree burns were treated with either a hydrocolloid adhesive bandage (control) (n = 16) or fish skin (n = 16) 7 days post-burn. Mice were subjected to histologic, hematologic, molecular, and gross evaluation at days 7, 16, and 28 post-burn. The fish skin offered no benefit to overall wound closure compared to hydrocolloids. Additionally, we detected no difference between fish skin and control treatments in regard to hypermetabolism or hematologic values. However, the fish skin groups exhibited 2 times more vascularization and 2 times higher expression of antimicrobial defensin peptide in comparison to controls. Proteomic analysis of the fish skin revealed the presence of antimicrobial peptides. Collectively, these data suggest that fish skin can serve as an innovative and cost-effective therapeutic alternative for burn victims to facilitate vascularization and reduce bacterial infection.

Combination of metabolome and proteome analyses provides insights into the mechanism underlying growth differences in Acipenser dabryanus

To analyze the differences between different-sized Acipenser dabryanus, we randomly selected 600 3-month-old A. dabryanus juveniles. Four months later, the blood and white muscle of these fish were analyzed. The results showed no significant difference in the length-weight relationship (LWR) b value between the large and small A. dabryanus. The levels of serum growth hormone (gh) and insulin-like growth factor 1 (igf1) in the large A. dabryanus were significantly lower than those in the small, whereas the activity levels of Total superoxide dismutase (T-sod) and catalase (cat) were opposite to the results of gh and igf1. A total of 212 and 245 metabolites showed significant changes in the positive and negative polarity mode, respectively. Among 3,308 proteins identified, 69 proteins showed upregulated expression, and 185 proteins showed downregulated expression. These results indicated that the growth advantage of A. dabryanus was closely related to glycolysis, protein synthesis, and antioxidant function.

Transcriptomic and proteomic analyses provide insights into the adaptive responses to the combined impact of salinity and alkalinity in Gymnocypris przewalskii

Gymnocypris przewalskii is the only high-land endemic teleost living in Qinghai Lake, the largest saline-alkaline lake in China. Its osmoregulatory physiology remains elusive due to a lack of precise identification of the response proteins. In the present study, DIA/SWATH was used to identify differentially expressed proteins (DEPs) under alkaline (pH = 10.1, carbonate buffer), saline (12‰, sodium chloride), and saline-alkaline [carbonate buffer (pH = 10.1) plus 11‰ sodium chloride] stresses. A total of 66,056 unique peptides representing 7,150 proteins and 230 DEPs [the false discovery rate (FDR) ≤ 0.05, fold change (FC) ≥ 1.5] were identified under different stresses. Comparative analyses of the proteome and transcriptome indicated that over 86% of DEPs did not show consistent trends with mRNA. In addition to consistent enrichment results under different stresses, the specific DEPs involved in saline-alkaline adaptation were primarily enriched in functions of homeostasis, hormone synthesis and reactions of defense response, complement activation and reproductive development. Meanwhile, a protein-protein interaction (PPI) network analysis of these specific DEPs indicated that the hub genes were ITGAX, MMP9, C3, F2, CD74, BTK, ANXA1, NCKAP1L, and CASP8. This study accurately isolated the genes that respond to stress, and the results could be helpful for understanding the physiological regulation mechanisms regarding salinity, alkalinity, and salinity-alkalinity interactions.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

Integrated Analysis of Liver Transcriptome, miRNA, and Proteome of Chinese Indigenous Breed Ningxiang Pig in Three Developmental Stages Uncovers Significant miRNA-mRNA-Protein Networks in Lipid Metabolism

Liver is an important metabolic organ of mammals. During each transitional period of life, liver metabolism is programmed by a complex molecular regulatory system for multiple physiological functions, many pathways of which are regulated by hormones and cytokines, nuclear receptors, and transcription factors. To gain a comprehensive and unbiased molecular understanding of liver growth and development in Ningxiang pigs, we analyzed the mRNA, microRNA (miRNA), and proteomes of the livers of Ningxiang pigs during lactation, nursery, and fattening periods. A total of 22,411 genes (19,653 known mRNAs and 2758 novel mRNAs), 1122 miRNAs (384 known miRNAs and 738 novel miRNAs), and 1123 unique proteins with medium and high abundance were identified by high-throughput sequencing and mass spectrometry. We show that the differences in transcriptional, post-transcriptional, or protein levels were readily identified by comparing different time periods, providing evidence that functional changes that may occur during liver development are widespread. In addition, we found many overlapping differentially expressed genes (DEGs)/differentially expressed miRNAs (DEMs)/differentially expressed proteins (DEPs) related to glycolipid metabolism in any group comparison. These overlapping DEGs/DEMs/DGPs may play an important role in functional transformation during liver development. Short Time-series Expression Miner (STEM) analysis revealed multiple expression patterns of mRNA, miRNA, and protein in the liver. Furthermore, several diverse key Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including immune defense, glycolipid metabolism, protein transport and uptake, and cell proliferation and development, were identified by combined analysis of DEGs and DGPs. A number of predicted miRNA–mRNA–protein pairs were found and validated by qRT-PCR and parallel reaction monitoring (PRM) assays. The results provide new and important information about the genetic breeding of Ningxiang pigs, which represents a foundation for further understanding the molecular regulatory mechanisms of dynamic development of liver tissue, functional transformation, and lipid metabolism.

Nonlinear effects of environmental salinity on the gill transcriptome versus proteome of Oreochromis niloticus modulate epithelial cell turnover

A data-independent acquisition (DIA) assay library for targeted quantitation of thousands of Oreochromis niloticus gill proteins using a label- and gel-free workflow was generated and used to compare protein and mRNA abundances. This approach generated complimentary rather than redundant data for 1899 unique genes in gills of tilapia acclimated to freshwater and brackish water. Functional enrichment analyses identified mitochondrial energy metabolism, serine protease and immunity-related functions, and cytoskeleton/ extracellular matrix organization as major processes controlled by salinity in O. niloticus gills. Non-linearity in salinity-dependent transcriptome versus proteome regulation was revealed for specific functional groups of genes. The relationship was more linear for other molecular functions/ cellular processes, suggesting that the salinity-dependent regulation of O. niloticus gill function relies on post-transcriptional mechanisms for some functions/ processes more than others. This integrative systems biology approach can be adopted for other tissues and organisms to study cellular dynamics for many changing ecological contexts.

A data-independent acquisition (DIA) assay library for quantitation of environmental effects on the kidney proteome of Oreochromis niloticus

Interactions of organisms with their environment are complex and environmental regulation at different levels of biological organization is often nonlinear. Therefore, the genotype to phenotype continuum requires study at multiple levels of organization. While studies of transcriptome regulation are now common for many species, quantitative studies of environmental effects on proteomes are needed. Here we report the generation of a data-independent acquisition (DIA) assay library that enables simultaneous targeted proteomics of thousands of Oreochromis niloticus kidney proteins using a label- and gel-free workflow that is well suited for ecologically relevant field samples. We demonstrate the usefulness of this DIA assay library by discerning environmental effects on the kidney proteome of O. niloticus. Moreover, we demonstrate that the DIA assay library approach generates data that are complimentary rather than redundant to transcriptomic data. Transcript and protein abundance differences in kidneys of tilapia acclimated to freshwater and brackish water (25 g/kg) were correlated for 2114 unique genes. A high degree of non-linearity in salinity-dependent regulation of transcriptomes and proteomes was revealed suggesting that the regulation of O. niloticus renal function by environmental salinity relies heavily on post-transcriptional mechanisms. The application of functional enrichment analyses using STRING and KEGG to DIA assay data sets is demonstrated by identifying myo-inositol metabolism, antioxidant and xenobiotic functions, and signalling mechanisms as key elements controlled by salinity in tilapia kidneys. The DIA assay library resource presented here can be adopted for other tissues and other organisms to study proteome dynamics during changing ecological contexts.

NOD-like receptor signaling pathway activation: A potential mechanism underlying negative effects of benzo(α)pyrene on zebrafish

Benzo(α)pyrene (BaP) is one of typical polycyclic aromatic hydrocarbons (PAHs) in aquatic environments and has been shown to cause toxic effects to aquatic animals. Although the negative effects of BaP have been investigated, the potential toxic mechanisms remain uncharacterized. To explore the potential mechanisms mediating the toxic effects of BaP, zebrafish (Danio rerio) were exposed to BaP for 15 days and the toxic effects of BaP in zebrafish liver were investigated using physiological and transcriptomic analyses. After 15-day BaP exposure, zebrafish liver exhibited abnormalities including increased cytoplasmic vacuolation, inflammatory cell infiltration, swelled nuclei and irregular pigmentation. BaP exposure also induced oxidative stress to the liver of zebrafish. Transcriptomic profiles revealed 5129 differentially expressed genes (DEGs) after 15-days of BaP exposure, and the vast majority of DEGs were up-regulated under BaP treatment. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggest that genes related to immune response were significantly dysregulated. Furthermore, the nucleotide-binding, oligomerization domain (NOD)-like receptor signaling pathway was significantly enriched and most of the genes in this pathway exhibited enhanced expression after BaP exposure. These results partially explained the mechanisms underlying the toxic effects of BaP on zebrafish liver. In conclusion, BaP has the potential to induce physiological responses in zebrafish liver through altering associated genes.

Proteome analysis of rainbow trout (Oncorhynchus mykiss) liver responses to chronic heat stress using DIA/SWATH

Aquaculture of rainbow trout (Oncorhynchus mykiss) is severely hampered by high temperatures in summer, and understanding the regulatory mechanisms controlling responses to chronic heat stress may assist the development of measures to relieve heat stress. In the present study, biochemical parameters revealed a strong stress response in rainbow trout at 24 °C, including activation of stress defence and immune systems. Liver proteome analysis under heat stress (24 °C) and control (18 °C) conditions using DIA/SWATH identified precursors (90,827), peptides (67,028), proteins (6770) and protein groups (5124), among which 460 differentially abundant proteins (DAPs; q-value < 0.05, fold change >1.5), 201 and 259 were up- and down-regulated, respectively. Many were related to heat shock proteins (HSPs), metabolism and immunity. Gene Ontology (GO) analysis showed that some DAPs induced at high temperature were involved in regulating cell homeostasis, metabolism, adaptive stress and stimulation. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified metabolic pathways, protein processing in endoplasmic reticulum, PPAR signalling, and complement and coagulation cascades. Protein-protein interaction (PPI) network analysis indicated that HSP90b1 and C3 may cooperative to affect cell membrane integrity under heat stress. Our findings assist the development of strategies to relieve heat stress in rainbow trout.

Identification of key proteins involved in stickleback environmental adaptation with system-level analysis

Using abundance measurements of 1,490 proteins from four separate populations of three-spined sticklebacks, we implemented a system-level approach to correlate proteome dynamics with environmental salinity and temperature and the fish's population and morphotype. We identified robust and accurate fingerprints that classify environmental salinity, temperature, morphotype, and the population sample origin, observing that proteins with specific functions are enriched in these fingerprints. Highly apparent functions represented in all fingerprints include ion transport, proteostasis, growth, and immunity, suggesting that these functions are most diversified in populations inhabiting different environments. Applying a differential network approach, we analyzed the network of protein interactions that differs between populations. Looking at specific population combinations of differential interaction, we identify sets of connected proteins. We find that these sets and their corresponding enriched functions reflect key processes that have diverged between the four populations. Moreover, the extent of divergence, i.e., the number of enriched functions that differ between populations, is highest when all three environmental parameters are different between two populations. Key nodes in the differential interaction network signify functions that are also inherent in the fingerprints, most prominently proteostasis-related functions. However, the differential interaction network also reveals additional functions that have diverged between populations, notably cytoskeletal organization and morphogenesis. The strength of these analyses is that the results are purely data driven. With such an unbiased approach applied on a large proteomic data set, we find the strongest signals given by the data, making it possible to develop more discriminatory and complex biomarkers for specific contexts of interest.

Proteomics of Osmoregulatory Responses in Threespine Stickleback Gills

The gill proteome of threespine sticklebacks (Gasterosteus aculeatus) differs greatly in populations that inhabit diverse environments characterized by different temperature, salinity, food availability, parasites, and other parameters. To assess the contribution of a specific environmental parameter to such differences it is necessary to isolate its effects from those of other parameters. In this study the effect of environmental salinity on the gill proteome of G. aculeatus was isolated in controlled mesocosm experiments. Salinity-dependent changes in the gill proteome were analyzed by Liquid chromatography/Tandem mass spectrometry data-independent acquisition (DIA) and Skyline. Relative abundances of 1691 proteins representing the molecular phenotype of stickleback gills were quantified using previously developed MSMS spectral and assay libraries in combination with DIA quantitative proteomics. Non-directional stress responses were distinguished from osmoregulatory protein abundance changes by their consistent occurrence during both hypo- and hyper-osmotic salinity stress in six separate mesocosm experiments. If the abundance of a protein was consistently regulated in opposite directions by hyper- versus hypo-osmotic salinity stress, then it was considered an osmoregulatory protein. In contrast, if protein abundance was consistently increased irrespective of whether salinity was increased or decreased, then it was considered a non-directional response protein. KEGG pathway analysis revealed that the salivary secretion, inositol phosphate metabolism, valine, leucine, and isoleucine degradation, citrate cycle, oxidative phosphorylation, and corresponding endocrine and extracellular signaling pathways contain most of the osmoregulatory gill proteins whose abundance is directly proportional to environmental salinity. Most proteins that were inversely correlated with salinity map to KEGG pathways that represent proteostasis, immunity, and related intracellular signaling processes. Non-directional stress response proteins represent fatty and amino acid degradation, purine metabolism, focal adhesion, mRNA surveillance, phagosome, endocytosis, and associated intracellular signaling KEGG pathways. These results demonstrate that G. aculeatus responds to salinity changes by adjusting osmoregulatory mechanisms that are distinct from transient non-directional stress responses to control compatible osmolyte synthesis, transepithelial ion transport, and oxidative energy metabolism. Furthermore, this study establishes salinity as a key factor for causing the regulation of numerous proteins and KEGG pathways with established functions in proteostasis, immunity, and tissue remodeling. We conclude that the corresponding osmoregulatory gill proteins and KEGG pathways represent molecular phenotypes that promote transepithelial ion transport, cellular osmoregulation, and gill epithelial remodeling to adjust gill function to environmental salinity.

Proteomic characterization and comparison of ram (Ovis aries) and buck (Capra hircus) spermatozoa proteome using a data independent acquisition mass spectometry (DIA-MS) approach

Fresh semen is most commonly used in an artificial insemination of small ruminants, because of low fertility rates of frozen sperm. Generally, when developing and applying assisted reproductive technologies, sheep and goats are classified as one species. In order to optimize sperm cryopreservation protocols in sheep and goat, differences in sperm proteomes between ram and buck are necessary to investigate, which may contribute to differences in function and fertility of spermatozoa. In the current work, a data-independent acquisition-mass spectrometry proteomic approach was used to characterize and make a comparison of ram (Ovis aries) and buck (Capra hircus) sperm proteomes. A total of 2,109 proteins were identified in ram and buck spermatozoa, with 238 differentially abundant proteins. Proteins identified in ram and buck spermatozoa are mainly involved in metabolic pathways for generation of energy and diminishing oxidative stress. Specifically, there are greater abundance of spermatozoa proteins related to the immune protective and capacity activities in ram, while protein that inhibit sperm capacitation shows greater abundance in buck. Our results not only provide novel insights into the characteristics and potential activities of spermatozoa proteins, but also expand the potential direction for sperm cryopreservation in ram and buck.

Proteomics for understanding pathogenesis, immune modulation and host pathogen interactions in aquaculture

Proteomic analyses techniques are considered strong tools for identifying and quantifying the protein contents in different organisms, organs and secretions. In fish biotechnology, the proteomic analyses have been used for wide range of applications such as identification of immune related proteins during infections and stresses. The proteomic approach has a significant role in understanding pathogen surviving strategies, host defence responses and subsequently, the fish pathogen interactions. Proteomic analyses were employed to highlight the virulence related proteins secreted by the pathogens to invade the fish host's defence barriers and to monitor the kinetics of protein contents of different fish organs in response to infections. The immune related proteins of fish and the virulence related proteins of pathogens are up or down regulated according to their functions in defence or pathogenesis. Therefore, the proteomic analyses are useful in understanding the virulence mechanisms of microorganisms and the fish pathogen interactions thereby supporting the development of new effective therapies. In this review, we focus and summarise the recent proteomic profiling studies exploring pathogen virulence activities and fish immune responses to stressors and infections.

Integrated Analysis of Transcriptomic, miRNA and Proteomic Changes of a Novel Hybrid Yellow Catfish Uncovers Key Roles for miRNAs in Heterosis

Heterosis is a complex biological phenomenon in which hybridization produces offspring that exhibit superior phenotypic characteristics compared with the parents. Heterosis is widely utilized in agriculture, for example in fish farming; however, its underlying molecular basis remains elusive. To gain a comprehensive and unbiased molecular understanding of fish heterosis, we analyzed the mRNA, miRNA, and proteomes of the livers of three catfish species, Pelteobagrus fulvidraco, P. vachelli, and their hybrid, the hybrid yellow catfish "Huangyou-1" (P. fulvidraco ♀ × P. vachelli ♂). Using next-generation sequencing and mass spectrometry, we show that the nonadditive, homoeolog expression bias and expression level dominance pattern were readily identified at the transcriptional, post-transcriptional, or protein levels, providing the evidence for the widespread presence of dominant models during hybridization. A number of predicted miRNA-mRNA-protein pairs were found and validated by qRT-PCR and PRM assays. Furthermore, several diverse key pathways were identified, including immune defense, metabolism, digestion and absorption, and cell proliferation and development, suggesting the vital mechanisms involved in the generation of the heterosis phenotype in progenies. We propose that the high parental expression of genes/proteins (growth, nutrition, feeding, and disease resistance) coupled with low parental miRNAs of the offspring, are inherited from the mother or father, thus indicating that the offspring were enriched with the advantages of the father or mother. We provide new and important information about the molecular mechanisms of heterosis, which represents a significant step toward a more complete elucidation of this phenomenon.