Population-specific renal proteomes of marine and freshwater three-spined sticklebacks

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.

Proteomic changes associated with predator-induced morphological defenses in oysters

Inducible prey defenses occur when organisms undergo plastic changes in phenotype to reduce predation risk. When predation pressure varies persistently over space or time, such as when predator and prey co-occur over only part of their biogeographic ranges, prey populations can become locally adapted in their inducible defenses. In California estuaries, native Olympia oyster (Ostrea lurida) populations have evolved disparate phenotypic responses to an invasive predator, the Atlantic oyster drill (Urosalpinx cinerea). In this study, oysters from an estuary with drills, and oysters from an estuary without drills, were reared for two generations in a laboratory common garden, and subsequently exposed to cues from Atlantic drills. Comparative proteomics was then used to investigate molecular mechanisms underlying conserved and divergent aspects of their inducible defenses. Both populations developed smaller, thicker, and harder shells after drill exposure, and these changes in shell phenotype were associated with up-regulation of calcium transport proteins that could influence biomineralization. Inducible defenses evolve in part because defended phenotypes incur fitness costs when predation risk is low. Immune proteins were down-regulated by both oyster populations after exposure to drills, implying a trade-off between biomineralization and immune function. Following drill exposure, oysters from the population that co-occurs with drills grew smaller shells than oysters inhabiting the estuary not yet invaded by the predator. Variation in the response to drills between populations was associated with isoform-specific protein expression. This trend suggests that a stronger inducible defense response evolved in oysters that co-occur with drills through modification of an existing mechanism.

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.

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.

Early-life exposure to the endocrine disruptor 17-α-ethinylestradiol induces delayed effects in adult brain, liver and ovotestis proteomes of a self-fertilizing fish

Early-life represents a critically sensitive window to endocrine disrupting chemicals, potentially leading to long-term repercussions on the phenotype later in life. The mechanisms underlying this phenomenon, referred to as the Developmental Origins of Health and Disease (DOHaD), are still poorly understood. To gain molecular understanding of these effects, we exposed mangrove rivulus (Kryptolebias marmoratus) for 28 days post hatching (dph) to 4 and 120 ng/L 17-α-ethinylestradiol, a model xenoestrogen. After 28 days, fish were raised for 140 days in clean water and we performed quantitative label-free proteomics on brain, liver and ovotestis of 168 dph adults. A total of 820, 888 and 420 proteins were robustly identified in the brain, liver and ovotestis, respectively. Effects of 17-α-ethinylestradiol were tissue and dose-dependent: a total of 31, 51 and 18 proteins were differentially abundant at 4 ng/L in the brain, liver and ovotestis, respectively, compared to 20, 25 and 39 proteins at 120 ng/L. Our results suggest that estrogen-responsive pathways, such as lipid metabolism, inflammation, and the innate immune system were affected months after the exposure. In addition, the potential perturbation of S-adenosylmethionine metabolism encourages future studies to investigate the role of DNA methylation in mediating the long-term effects of early-life exposures. SIGNIFICANCE: The Developmental Origins of Health and Disease (DOHaD) states that early life stages of humans and animals are sensitive to environmental stressors and can develop health issues later in life, even if the stress has ceased. Molecular mechanisms supporting DOHaD are still unclear. The mangrove rivulus is a new fish model species naturally reproducing by self-fertilization, making it possible to use isogenic lineages in which all individuals are highly homozygous. This species therefore permits to strongly reduce the confounding factor of genetic variability in order to investigate the effects of environmental stress on the phenotype. After characterizing the molecular phenotype of brain, liver and ovotestis, we obtained true proteomic reaction norms of these three organs in adults after early life stages have been exposed to the common endocrine disruptor 17-α-ethinylestradiol (EE2). Our study demonstrates long-term effects of early-life endocrine disruption at the proteomic level in diverse estrogen-responsive pathways 5 months after the exposure. The lowest tested and environmentally relevant concentration of 4 ng/L had the highest impact on the proteome in brain and liver, highlighting the potency of endocrine disruptors at low concentrations.

Integrated Proteomic and Transcriptomic Analysis of Differential Expression of Chicken Lung Tissue in Response to NDV Infection during Heat Stress

Newcastle disease virus (NDV) is a devastating worldwide poultry pathogen with major implications for global food security. In this study, two highly inbred and genetically distinct chicken lines, Fayoumis and Leghorns, were exposed to a lentogenic strain of NDV, while under the effects of heat stress, in order to understand the genetic mechanisms of resistance during high ambient temperatures. Fayoumis, which are relatively more resistant to pathogens than Leghorns, had larger numbers of differentially expressed genes (DEGs) during the early stages of infection when compared to Leghorns and subsequently down-regulated their immune response at the latter stages to return to homeostasis. Leghorns had very few DEGs across all observed time points, with the majority of DEGs involved with metabolic and glucose-related functions. Proteomic analysis corroborates findings made within Leghorns, while also identifying interesting candidate genes missed by expression profiling. Poor correlation between changes observed in the proteomic and transcriptomic datasets highlights the potential importance of integrative approaches to understand the mechanisms of disease response. Overall, this study provides novel insights into global protein and expression profiles of these two genetic lines, and provides potential genetic targets involved with NDV resistance during heat stress in poultry.

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.

Comparative analysis of weighted gene co-expression networks in human and mouse

The application of complex network modeling to analyze large co-expression data sets has gained traction during the last decade. In particular, the use of the weighted gene co-expression network analysis framework has allowed an unbiased and systems-level investigation of genotype-phenotype relationships in a wide range of systems. Since mouse is an important model organism for biomedical research on human disease, it is of great interest to identify similarities and differences in the functional roles of human and mouse orthologous genes. Here, we develop a novel network comparison approach which we demonstrate by comparing two gene-expression data sets from a large number of human and mouse tissues. The method uses weighted topological overlap alongside the recently developed network-decomposition method of s-core analysis, which is suitable for making gene-centrality rankings for weighted networks. The aim is to identify globally central genes separately in the human and mouse networks. By comparing the ranked gene lists, we identify genes that display conserved or diverged centrality-characteristics across the networks. This framework only assumes a single threshold value that is chosen from a statistical analysis, and it may be applied to arbitrary network structures and edge-weight distributions, also outside the context of biology. When conducting the comparative network analysis, both within and across the two species, we find a clear pattern of enrichment of transcription factors, for the homeobox domain in particular, among the globally central genes. We also perform gene-ontology term enrichment analysis and look at disease-related genes for the separate networks as well as the network comparisons. We find that gene ontology terms related to regulation and development are generally enriched across the networks. In particular, the genes FOXE3, RHO, RUNX2, ALX3 and RARA, which are disease genes in either human or mouse, are on the top-10 list of globally central genes in the human and mouse networks.

A composite network of conserved and tissue specific gene interactions reveals possible genetic interactions in glioma

Differential co-expression network analyses have recently become an important step in the investigation of cellular differentiation and dysfunctional gene-regulation in cell and tissue disease-states. The resulting networks have been analyzed to identify and understand pathways associated with disorders, or to infer molecular interactions. However, existing methods for differential co-expression network analysis are unable to distinguish between various forms of differential co-expression. To close this gap, here we define the three different kinds (conserved, specific, and differentiated) of differential co-expression and present a systematic framework, CSD, for differential co-expression network analysis that incorporates these interactions on an equal footing. In addition, our method includes a subsampling strategy to estimate the variance of co-expressions. Our framework is applicable to a wide variety of cases, such as the study of differential co-expression networks between healthy and disease states, before and after treatments, or between species. Applying the CSD approach to a published gene-expression data set of cerebral cortex and basal ganglia samples from healthy individuals, we find that the resulting CSD network is enriched in genes associated with cognitive function, signaling pathways involving compounds with well-known roles in the central nervous system, as well as certain neurological diseases. From the CSD analysis, we identify a set of prominent hubs of differential co-expression, whose neighborhood contains a substantial number of genes associated with glioblastoma. The resulting gene-sets identified by our CSD analysis also contain many genes that so far have not been recognized as having a role in glioblastoma, but are good candidates for further studies. CSD may thus aid in hypothesis-generation for functional disease-associations.

Comparative proteomics analysis of teleost intermuscular bones and ribs provides insight into their development

BACKGROUND

Intermuscular bones (IBs) and ribs both are a part of skeletal system in teleosts, but with different developing process. The chemical composition of fish IBs and ribs as well as the underlying mechanism about their development have not been investigated. In the present study, histological structures showed that one bone cavity containing osteoclasts were existed in ribs, but not in IBs of Megalobrama amblycephala. We constructed the first proteomics map for fish bones including IBs and ribs, and identified the differentially expressed proteins between IBs and ribs through iTRAQ LC-MS/MS proteomic analysis.

RESULTS

The proteins extracted from IBs and ribs at 1- to 2-year old M. amblycephala were quantified 2,342 proteins, with 1,451 proteins annotated with GO annotation in biological processes, molecular function and cellular component. A number of bone related proteins as well as pathways were identified in the study. A total of 93 and 154 differently expressed proteins were identified in comparison groups of 1-IB-vs-1-Rib and 2-IB-vs-2-Rib, which indicated the obvious differences of chemical composition between these two bone tissues. The two proteins (vitronectin b precursor and matrix metalloproteinase-2) related to osteoclasts differentiation were significantly up-regulated in ribs compared with IBs (P < 0.05), which was in accordance with the results from histological structures. In comparison groups of 1-IB-vs-2-IB and 1-Rib-vs-2-Rib, 33 and 51 differently expressed proteins were identified and the function annotation results showed that these proteins were involved in regulating bone development and differentiation. Subsequently, 11 and 13 candidate proteins in comparison group of 1-IB-vs-1-Rib and 1-IB-vs-2-IB related to bone development were validated by MRM assays.

CONCLUSIONS

Our present study suggested the different key proteins involved in the composition of fish ribs and IBs as well as their growth development. These findings could provide important clues towards further understanding of fish skeletal system and the roles of proteins playing in regulating diverse biological processes in fish.