Proteomics in evolutionary ecology: linking the genotype with the phenotype

Ocean warming alters cellular metabolism and induces mortality in fish early life stages: A proteomic approach

Climate change has pervasive effects on marine ecosystems, altering biodiversity patterns, abundance and distribution of species, biological interactions, phenology, and organisms' physiology, performance and fitness. Fish early life stages have narrow thermal windows and are thus more vulnerable to further changes in water temperature. The aim of this study was to address the sensitivity and underlying molecular changes of larvae of a key fisheries species, the sea bream Sparus aurata, towards ocean warming. Larvae were exposed to three temperatures: 18°C (control), 24°C (warm) and 30°C (heat wave) for seven days. At the end of the assay, i) survival curves were plotted for each temperature treatment and ii) entire larvae were collected for proteomic analysis via 2D gel electrophoresis, image analysis and mass spectrometry. Survival decreased with increasing temperature, with no larvae surviving at 30°C. Therefore, proteomic analysis was only carried out for 18°C and 24°C. Larvae up-regulated protein folding and degradation, cytoskeletal re-organization, transcriptional regulation and the growth hormone while mostly down-regulating cargo transporting and porphyrin metabolism upon exposure to heat stress. No changes were detected in proteins related to energetic metabolism suggesting that larval fish may not have the energetic plasticity needed to sustain cellular protection in the long-term. These results indicate that despite proteome modulation, S. aurata larvae do not seem able to fully acclimate to higher temperatures as shown by the low survival rates. Consequently, elevated temperatures seem to have bottleneck effects during fish early life stages, and future ocean warming can potentially compromise recruitment's success of key fisheries species.

Venomic Analysis of the Poorly Studied Desert Coral Snake, Micrurus tschudii tschudii, Supports the 3FTx/PLA₂ Dichotomy across Micrurus Venoms

The venom proteome of the poorly studied desert coral snake Micrurus tschudii tschudii was unveiled using a venomic approach, which identified ≥38 proteins belonging to only four snake venom protein families. The three-finger toxins (3FTxs) constitute, both in number of isoforms (~30) and total abundance (93.6% of the venom proteome), the major protein family of the desert coral snake venom. Phospholipases A₂ (PLA₂s; seven isoforms, 4.1% of the venom proteome), 1–3 Kunitz-type proteins (1.6%), and 1–2 l-amino acid oxidases (LAO, 0.7%) complete the toxin arsenal of M. t. tschudii. Our results add to the growing evidence that the occurrence of two divergent venom phenotypes, i.e., 3FTx- and PLA₂-predominant venom proteomes, may constitute a general trend across the cladogenesis of Micrurus. The occurrence of a similar pattern of venom phenotypic variability among true sea snake (Hydrophiinae) venoms suggests that the 3FTx/PLA₂ dichotomy may be widely distributed among Elapidae venoms.

Quantitative analysis of oyster larval proteome provides new insights into the effects of multiple climate change stressors

The metamorphosis of planktonic larvae of the Pacific oyster (Crassostrea gigas) underpins their complex life-history strategy by switching on the molecular machinery required for sessile life and building calcite shells. Metamorphosis becomes a survival bottleneck, which will be pressured by different anthropogenically induced climate change-related variables. Therefore, it is important to understand how metamorphosing larvae interact with emerging climate change stressors. To predict how larvae might be affected in a future ocean, we examined changes in the proteome of metamorphosing larvae under multiple stressors: decreased pH (pH 7.4), increased temperature (30 °C), and reduced salinity (15 psu). Quantitative protein expression profiling using iTRAQ-LC-MS/MS identified more than 1300 proteins. Decreased pH had a negative effect on metamorphosis by down-regulating several proteins involved in energy production, metabolism, and protein synthesis. However, warming switched on these down-regulated pathways at pH 7.4. Under multiple stressors, cell signaling, energy production, growth, and developmental pathways were up-regulated, although metamorphosis was still reduced. Despite the lack of lethal effects, significant physiological responses to both individual and interacting climate change related stressors were observed at proteome level. The metamorphosing larvae of the C. gigas population in the Yellow Sea appear to have adequate phenotypic plasticity at the proteome level to survive in future coastal oceans, but with developmental and physiological costs.

Proteomics studies in inner ear disorders: pathophysiology and biomarkers

Although proteomics has been exploited in a wide range of diseases for identification of biomarkers and pathophysiological mechanisms, there are still biomedical disciplines such as otology where proteomics platforms are underused due to technical challenges and/or complex features of the disease. Thus, in the past few years, healthcare and scientific agencies have advocated the development and adoption of proteomic technologies in otological research. However, few studies have been conducted and limited literature is available in this area. Here, we present the state of the art of proteomics in otology, discussing the substantial evidence from recent experimental models and clinical studies in inner-ear conditions. We also delineate a series of critical issues including minute size of the inner ear, delicacy and poor accessibility of tissue that researchers face while undertaking otology proteomics research. Furthermore, we provide perspective to enhance the impact and lead to the clinical implementation of these proteomics-based strategies.

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.

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.

Proteomics in evolutionary ecology

Evolutionary ecologists are traditionally gene-focused, as genes propagate phenotypic traits across generations and mutations and recombination in the DNA generate genetic diversity required for evolutionary processes. As a consequence, the inheritance of changed DNA provides a molecular explanation for the functional changes associated with natural selection. A direct focus on proteins on the other hand, the actual molecular agents responsible for the expression of a phenotypic trait, receives far less interest from ecologists and evolutionary biologists. This is partially due to the central dogma of molecular biology that appears to define proteins as the 'dead-end of molecular information flow' as well as technical limitations in identifying and studying proteins and their diversity in the field and in many of the more exotic genera often favored in ecological studies. Here we provide an overview of a newly forming field of research that we refer to as 'Evolutionary Proteomics'. We point out that the origins of cellular function are related to the properties of polypeptide and RNA and their interactions with the environment, rather than DNA descent, and that the critical role of horizontal gene transfer in evolution is more about coopting new proteins to impact cellular processes than it is about modifying gene function. Furthermore, post-transcriptional and post-translational processes generate a remarkable diversity of mature proteins from a single gene, and the properties of these mature proteins can also influence inheritance through genetic and perhaps epigenetic mechanisms. The influence of post-transcriptional diversification on evolutionary processes could provide a novel mechanistic underpinning for elements of rapid, directed evolutionary changes and adaptations as observed for a variety of evolutionary processes. Modern state-of the art technologies based on mass spectrometry are now available to identify and quantify peptides, proteins, protein modifications and protein interactions of interest with high accuracy and assess protein diversity and function. Therefore, proteomic technologies can be viewed as providing evolutionary biologist with exciting novel opportunities to understand very early events in functional variation of cellular molecular machinery that are acting as part of evolutionary processes.

Quantitative Mass Spectrometry Reveals Partial Translational Regulation for Dosage Compensation in Chicken

There is increasing evidence that dosage compensation is not a ubiquitous feature following sex chromosome evolution, especially not in organisms where females are the heterogametic sex, like in birds. Even when it occurs, compensation can be incomplete and limited to dosage-sensitive genes. However, previous work has mainly studied transcriptional regulation of sex-linked genes, which may not reflect expression at the protein level. Here, we used liquid chromatography-tandem mass spectrometry to detect and quantify expressed levels of more than 2,400 proteins in ten different tissues of male and female chicken embryos. For comparison, transcriptome sequencing was performed in the same individuals, five of each sex. The proteomic analysis revealed that dosage compensation was incomplete, with a mean male-to-female (M:F) expression ratio of Z-linked genes of 1.32 across tissues, similar to that at the RNA level (1.29). The mean Z chromosome-to-autosome expression ratio was close to 1 in males and lower than 1 in females, consistent with partly reduced Z chromosome expression in females. Although our results exclude a general mechanism for chromosome-wide dosage compensation at translation, 30% of all proteins encoded from Z-linked genes showed a significant change in the M:F ratio compared with the corresponding ratio at the RNA level. This resulted in a pattern where some genes showed balanced expression between sexes and some close to 2-fold higher expression in males. This suggests that proteomic analyses will be necessary to reveal a more complete picture of gene regulation and sex chromosome evolution.

Post-transcriptional Mechanisms Contribute Little to Phenotypic Variation in Snake Venoms

Protein expression is a major link in the genotype–phenotype relationship, and processes affecting protein abundances, such as rates of transcription and translation, could contribute to phenotypic evolution if they generate heritable variation. Recent work has suggested that mRNA abundances do not accurately predict final protein abundances, which would imply that post-transcriptional regulatory processes contribute significantly to phenotypes. Post-transcriptional processes also appear to buffer changes in transcriptional patterns as species diverge, suggesting that the transcriptional changes have little or no effect on the phenotypes undergoing study. We tested for concordance between mRNA and protein expression levels in snake venoms by means of mRNA-seq and quantitative mass spectrometry for 11 snakes representing 10 species, six genera, and three families. In contrast to most previous work, we found high correlations between venom gland transcriptomes and venom proteomes for 10 of our 11 comparisons. We tested for protein-level buffering of transcriptional changes during species divergence by comparing the difference between transcript abundance and protein abundance for three pairs of species and one intraspecific pair. We found no evidence for buffering during divergence of our three species pairs but did find evidence for protein-level buffering for our single intraspecific comparison, suggesting that buffering, if present, was a transient phenomenon in venom divergence. Our results demonstrated that post-transcriptional mechanisms did not contribute significantly to phenotypic evolution in venoms and suggest a more prominent and direct role for cis-regulatory evolution in phenotypic variation, particularly for snake venoms.

Protein expression parallels thermal tolerance and ecologic changes in the diversification of a diving beetle species complex

Physiological changes associated with evolutionary and ecological processes such as diversification, range expansion or speciation are still incompletely understood, especially for non-model species. Here we study differences in protein expression in response to temperature in a western Mediterranean diving beetle species complex, using two-dimensional differential gel electrophoresis with one Moroccan and one Iberian population each of Agabus ramblae and Agabus brunneus. We identified proteins with significant expression differences after thermal treatments comparing them with a reference EST library generated from one of the species of the complex (A. ramblae). The colonisation during the Middle Pleistocene of the Iberian peninsula by A. ramblae, where maximum temperatures and seasonality are lower than in the ancestral north African range, was associated with changes in the response to 27 °C in proteins related to energy metabolism. The subsequent speciation of A. brunneus from within populations of Iberian A. ramblae was associated with changes in the expression of several stress-related proteins (mostly chaperons) when exposed to 4 °C. These changes are in agreement with the known tolerance to lower temperatures of A. brunneus, which occupies a larger geographical area with a wider range of climatic conditions. In both cases, protein expression changes paralleled the evolution of thermal tolerance and the climatic conditions experienced by the species. However, although the colonisation of the Iberian peninsula did not result in morphological change, the speciation process of A. brunneus within Iberia involved genetic isolation and substantial differences in male genitalia and body size and shape.

Understanding the cephalopod immune system based on functional and molecular evidence

Cephalopods have the most advanced circulatory and nervous system among the mollusks. Recently, they have been included in the European directive which state that suffering and pain should be minimized in cephalopods used in experimentation. The knowledge about cephalopod welfare is still limited and several gaps are yet to be filled, especially in reference to pathogens, pathologies and immune response of these mollusks. In light of the requirements of the normative, in addition to the ecologic and economic importance of cephalopods, in this review we update the work published to date concerning cephalopod immune system. Significant advances have been reached in relation to the characterization of haemocytes and defensive mechanisms comprising cellular and humoral factors mainly, but not limited, in species of high economic value like Sepia officinalis and Octopus vulgaris. Moreover, the improvement of molecular approaches has helped to discover several immune-related genes/proteins. These immune genes/proteins include antimicrobial peptides, phenoloxidases, antioxidant enzymes, serine protease inhibitor, lipopolysaccharide-induced TNF-α factor, Toll-like receptors, lectins, even clusters of differentiation among others. Most of them have been found in haemocytes but also in gills and digestive gland, and the characterization as well as their precise role in the immune response of cephalopods is still pending to be elucidated. The assessment of immune parameters in cephalopods exposed to contaminants is just starting, but the negative impact of some pollutants on the immune response of the common octopus has been reported. This review summarizes the current status of our knowledge about the cephalopod immune system that seems to be far from simply. On the contrary, the advances gained to date point out a complex innate immunity in cephalopods.

Characterization of the genome and transcriptome of the blue tit Cyanistes caeruleus: polymorphisms, sex-biased expression and selection signals

Decoding genomic sequences and determining their variation within populations has potential to reveal adaptive processes and unravel the genetic basis of ecologically relevant trait variation within a species. The blue tit Cyanistes caeruleus--a long-time ecological model species--has been used to investigate fitness consequences of variation in mating and reproductive behaviour. However, very little is known about the underlying genetic changes due to natural and sexual selection in the genome of this songbird. As a step to bridge this gap, we assembled the first draft genome of a single blue tit, mapped the transcriptome of five females and five males to this reference, identified genomewide variants and performed sex-differential expression analysis in the gonads, brain and other tissues. In the gonads, we found a high number of sex-biased genes, and of those, a similar proportion were sex-limited (genes only expressed in one sex) in males and females. However, in the brain, the proportion of female-limited genes within the female-biased gene category (82%) was substantially higher than the proportion of male-limited genes within the male-biased category (6%). This suggests a predominant on-off switching mechanism for the female-limited genes. In addition, most male-biased genes were located on the Z-chromosome, indicating incomplete dosage compensation for the male-biased genes. We called more than 500,000 SNPs from the RNA-seq data. Heterozygote detection in the single reference individual was highly congruent between DNA-seq and RNA-seq calling. Using information from these polymorphisms, we identified potential selection signals in the genome. We list candidate genes which can be used for further sequencing and detailed selection studies, including genes potentially related to meiotic drive evolution. A public genome browser of the blue tit with the described information is available at http://public-genomes-ngs.molgen.mpg.de.

Phenotypic integration in the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus)

Selection can vary geographically across environments and temporally over the lifetime of an individual. Unlike geographic contexts, where different selective regimes can act on different alleles, age-specific selection is constrained to act on the same genome by altering age-specific expression. Snake venoms are exceptional traits for studying ontogeny because toxin expression variation directly changes the phenotype; relative amounts of venom components determine, in part, venom efficacy. Phenotypic integration is the dependent relationship between different traits that collectively produce a complex phenotype and, in venomous snakes, may include traits as diverse as venom, head shape and fang length. We examined the feeding system of the eastern diamondback rattlesnake (Crotalus adamanteus) across environments and over the lifetime of individuals and used a genotype-phenotype map approach, protein expression data and morphological data to demonstrate that: (i) ontogenetic effects explained more of the variation in toxin expression variation than geographic effects, (ii) both juveniles and adults varied geographically, (iii) toxin expression variation was a result of directional selection and (iv) different venom phenotypes covaried with morphological traits also associated with feeding in temporal (ontogenetic) and geographic (functional) contexts. These data are the first to demonstrate, to our knowledge, phenotypic integration between multiple morphological characters and a biochemical phenotype across populations and age classes. We identified copy number variation as the mechanism driving the difference in the venom phenotype associated with these morphological differences, and the parallel mitochondrial, venom and morphological divergence between northern and southern clades suggests that each clade may warrant classification as a separate evolutionarily significant unit.

Combined venomics, venom gland transcriptomics, bioactivities, and antivenomics of two Bothrops jararaca populations from geographic isolated regions within the Brazilian Atlantic rainforest

Bothrops jararaca is a slender and semi-arboreal medically relevant pit viper species endemic to tropical and subtropical forests in southern Brazil, Paraguay, and northern Argentina (Misiones). Within its geographic range, it is often abundant and is an important cause of snakebite. Although no subspecies are currently recognized, geographic analyses have revealed the existence of two well-supported B. jararaca clades that diverged during the Pliocene ~3.8Mya and currently display a southeastern (SE) and a southern (S) Atlantic rainforest (Mata Atlântica) distribution. The spectrum, geographic variability, and ontogenetic changes of the venom proteomes of snakes from these two B. jararaca phylogroups were investigated applying a combined venom gland transcriptomic and venomic analysis. Comparisons of the venom proteomes and transcriptomes of B. jararaca from the SE and S geographic regions revealed notable interpopulational variability that may be due to the different levels of population-specific transcriptional regulation, including, in the case of the southern population, a marked ontogenetic venom compositional change involving the upregulation of the myotoxic PLA2 homolog, bothropstoxin-I. This population-specific marker can be used to estimate the proportion of venom from the southern population present in the B. jararaca venom pool used for the Brazilian soro antibotrópico (SAB) antivenom production. On the other hand, the southeastern population-specific D49-PLA2 molecules, BinTX-I and BinTX-II, lend support to the notion that the mainland ancestor of Bothrops insularis was originated within the same population that gave rise to the current SE B. jararaca phylogroup, and that this insular species endemic to Queimada Grande Island (Brazil) expresses a pedomorphic venom phenotype. Mirroring their compositional divergence, the two geographic B. jararaca venom pools showed distinct bioactivity profiles. However, the SAB antivenom manufactured in Vital Brazil Institute neutralized the lethal effect of both venoms to a similar extent. In addition, immobilized SAB antivenom immunocaptured most of the venom components of the venoms of both B. jararaca populations, but did not show immunoreactivity against vasoactive peptides. The Costa Rican bothropic-crotalic-lachesic (BCL) antivenom showed the same lack of reactivity against vasoactive peptides but, in addition, was less efficient immunocapturing PI- and PIII-SVMPs from the SE venom, and bothropstoxin-I, a CRISP molecule, and a D49-PLA2 from the venom of the southern B. jararaca phylogroup. The remarkable paraspecificity exhibited by the Brazilian and the Costa Rican antivenoms indicates large immunoreactive epitope conservation across the natural history of Bothrops, a genus that has its roots in the middle Miocene. This article is part of a Special Issue entitled: Omics Evolutionary Ecolog.

Comparative proteomic analysis reveals the role of hydrogen sulfide in the adaptation of the alpine plant Lamiophlomis rotata to altitude gradient in the Northern Tibetan Plateau

We found the novel role of hydrogen sulfide in the adaptation of the alpine plant to altitude gradient in the Northern Tibetan Plateau. Alpine plants have developed strategies to survive the extremely cold conditions prevailing at high altitudes; however, the mechanism underlying the evolution of these strategies remains unknown. Hydrogen sulfide (H2S) is an essential messenger that enhances plant tolerance to environmental stress; however, its role in alpine plant adaptation to environmental stress has not been reported until now. In this work, we conducted a comparative proteomics analysis to investigate the dynamic patterns of protein expression in Lamiophlomis rotata plants grown at three different altitudes. We identified and annotated 83 differentially expressed proteins. We found that the levels and enzyme activities of proteins involved in H2S biosynthesis markedly increased at higher altitudes, and that H2S accumulation increased. Exogenous H2S application increased antioxidant enzyme activity, which reduced ROS (reactive oxygen species) damage, and GSNOR (S-nitrosoglutathione reductase) activity, which reduced RNS (reactive nitrogen species) damage, and activated the downstream defense response, resulting in protein degradation and proline and sugar accumulation. However, such defense responses could be reversed by applying H2S biosynthesis inhibitors. Based on these findings, we conclude that L. rotata uses multiple strategies to adapt to the alpine stress environment and that H2S plays a central role during this process.

Statistical analysis of the individual variability of 1D protein profiles as a tool in ecology: an application to parasitoid venom

Understanding the forces that shape eco-evolutionary patterns often requires linking phenotypes to genotypes, allowing characterization of these patterns at the molecular level. DNA-based markers are less informative in this aim compared to markers associated with gene expression and, more specifically, with protein quantities. The characterization of eco-evolutionary patterns also usually requires the analysis of large sample sizes to accurately estimate interindividual variability. However, the methods used to characterize and compare protein samples are generally expensive and time-consuming, which constrains the size of the produced data sets to few individuals. We present here a method that estimates the interindividual variability of protein quantities based on a global, semi-automatic analysis of 1D electrophoretic profiles, opening the way to rapid analysis and comparison of hundreds of individuals. The main original features of the method are the in silico normalization of sample protein quantities using pictures of electrophoresis gels at different staining levels, as well as a new method of analysis of electrophoretic profiles based on a median profile. We demonstrate that this method can accurately discriminate between species and between geographically distant or close populations, based on interindividual variation in venom protein profiles from three endoparasitoid wasps of two different genera (Psyttalia concolor, Psyttalia lounsburyi and Leptopilina boulardi). Finally, we discuss the experimental designs that would benefit from the use of this method.

Gene-expression novelty in allopolyploid cotton: a proteomic perspective

Allopolyploidization is accompanied by changes in gene expression that are thought to contribute to phenotypic diversification. Here we describe global changes in the single-celled cotton fiber proteome of two natural allopolyploid species (Gossypium hirsutum and G. barbadense) and living models of their diploid parents using two different proteomic approaches. In total, 1323 two-dimensional gel electrophoresis spots and 1652 identified proteins by isobaric tags for relative and absolute quantitation were quantitatively profiled during fiber elongation. Between allopolyploids and their diploid A- and D-genome progenitors, amounts of differential expression ranged from 4.4 to 12.8%. Over 80% of the allopolyploid proteome was additively expressed with respect to progenitor diploids. Interestingly, the fiber proteome of G. hirsutum resembles the parental A-genome more closely, where long, spinable fiber first evolved, than does the fiber proteome of G. barbadense. More protein expression patterns were A-dominant than D-dominant in G. hirsutum, but in G. barbadense, the direction of expression-level dominance switched from the D-genome to the A-genome during fiber development. Comparison of developmental changes between the two allopolyploid species revealed a high level of proteomic differentiation despite their shared ancestry, relatively recent evolutionary divergence, and similar gross morphology. These results suggest that the two allopolyploid species have achieved superficially similar modern fiber phenotypes through different evolutionary routes at the proteome level. We also detected homeolog-specific expression for 1001 proteins and present a novel approach to infer the relationship between homeolog-specific and duplicate expression patterns. Our study provides a proteomic perspective on understanding evolutionary consequences of allopolyploidization, showing how protein expression has been altered by polyploidization and subsequently has diversified among species.

Stressor-induced proteome alterations in zebrafish: a meta-analysis of response patterns

Proteomics approaches are being increasingly applied in ecotoxicology on the premise that the identification of specific protein expression changes in response to a particular chemical would allow elucidation of the underlying molecular pathways leading to an adverse effect. This in turn is expected to promote the development of focused testing strategies for specific groups of toxicants. Although both gel-based and gel-free global characterization techniques provide limited proteome coverage, the conclusions regarding the cellular processes affected are still being drawn based on the few changes detected. To investigate how specific the detected responses are, we analyzed a set of studies that characterized proteome alterations induced by various physiological, chemical and biological stressors in zebrafish, a popular model organism. Our analysis highlights several proteins and protein groups, including heat shock and oxidative stress defense proteins, energy metabolism enzymes and cytoskeletal proteins, to be most frequently identified as responding to diverse stressors. In contrast, other potentially more specifically responding protein groups are detected much less frequently. Thus, zebrafish proteome responses to stress reported by different studies appear to depend mostly on the level of stress rather than on the specific stressor itself. This suggests that the most broadly used current proteomics technologies do not provide sufficient proteome coverage to allow in-depth investigation of specific mechanisms of toxicant action. We suggest that the results of any differential proteomics experiment performed with zebrafish should be interpreted keeping in mind the list of the most frequent responders that we have identified. Similar reservations should apply to any other species where proteome responses are analyzed by global proteomics methods. Careful consideration of the reliability and significance of observed changes is necessary in order not to over-interpret the experimental results and to prevent the proliferation of false positive linkages between the chemical and the cellular functions it perturbs. We further discuss the implications of the identified "top lists" of frequently responding proteins and protein families, and suggest further directions for proteomics research in ecotoxicology. Apart from improving the proteome coverage, further research should focus on defining the significance of the observed stress response patterns for organism phenotypes and on searching for common upstream regulators that can be targeted by specific assays.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Spatial variation in transcript and protein abundance of Atlantic salmon during feeding migration in the Baltic Sea

The fitness and reproductive output of fishes can be affected by environmental disturbances. In this study, transcriptomics and label-free proteomics were combined to investigate Atlantic salmon (Salmo salar) sampled from three different field locations within the Baltic Sea (Baltic Main Basin (BMB), Gulf of Finland (GoF), and Bothnian Sea (BS)) during marine migration. The expression of several stress related mRNAs and proteins of xenobiotic metabolism, oxidative stress, DNA damage, and cell death were increased in salmon from GoF compared to salmon from BMB or BS. Respiratory electron chain and ATP synthesis related gene ontology-categories were upregulated in GoF salmon, whereas those associated with RNA processing and synthesis, translation, and protein folding decreased. Differences were seen also in metabolism and immune function related gene expression. Comparisons of the transcriptomic and proteomic profiles between salmon from GoF and salmon from BMB or BS suggest environmental stressors, especially exposure to contaminants, as a main explanation for differences. Salmon feeding in GoF are thus “disturbed by hazardous substances”. The results may also be applied in evaluating the conditions of pelagic ecosystems in the different parts of Baltic Sea.

Proteomics profiling of fiber development and domestication in upland cotton (Gossypium hirsutum L.)

Comparative proteomic analyses were performed to detail the evolutionary consequences of strong directional selection for enhanced fiber traits in modern upland cotton (Gossypium hirsutum L.). Using two complementary proteomic approaches, 2-DE and iTRAQ LC-MS/MS, fiber proteomes were examined for four representative stages of fiber development. Approximately 1,000 protein features were characterized using each strategy, collectively resulting in the identification and functional categorization of 1,223 proteins. Unequal contributions of homoeologous proteins were detected for over a third of the fiber proteome, but overall expression was balanced with respect to the genome-of-origin in the allopolyploid G. hirsutum. About 30% of the proteins were differentially expressed during fiber development within wild and domesticated cotton. Notably, domestication was accompanied by a doubling of protein developmental dynamics for the period between 10 and 20 days following pollination. Expression levels of 240 iTRAQ proteins and 293 2-DE spots were altered by domestication, collectively representing multiple cellular and metabolic processes, including metabolism, energy, protein synthesis and destination, defense and stress response. Analyses of homoeolog-specific expression indicate that duplicated gene products in cotton fibers can be differently regulated in response to selection. These results demonstrate the power of proteomics for the analysis of crop domestication and phenotypic evolution.

The proteome of Atlantic herring (Clupea harengus L.) larvae is resistant to elevated pCO2

Elevated anthropogenic pCO2 can delay growth and impair otolith structure and function in the larvae of some fishes. These effects may concurrently alter the larva's proteome expression pattern. To test this hypothesis, Atlantic herring larvae were exposed to ambient (370 μatm) and elevated (1800 μatm) pCO2 for one-month. The proteome structure of the larvae was examined using a 2-DE and mass spectrometry. The length of herring larvae was marginally less in the elevated pCO2 treatment compared to the control. The proteome structure was also different between the control and treatment, but only slightly: the expression of a small number of proteins was altered by a factor of less than 2-fold at elevated pCO2. This comparative proteome analysis suggests that the proteome of herring larvae is resilient to elevated pCO2. These observations suggest that herring larvae can cope with levels of CO2 projected for near future without significant proteome-wide changes.

Physiological, biochemical and proteomics analysis reveals the adaptation strategies of the alpine plant Potentilla saundersiana at altitude gradient of the Northwestern Tibetan Plateau

This study presents an analysis of leave and rood morphology, biochemical and proteomics approach as adaptation strategies of the alpine plant Potentilla saundersiana in an altitude gradient. Several plant physiological parameter, including root and leaf architecture, leaf photosynthesis capacity, specific leaf area (SLA) and leaf nitrogen concentration, histology and microscopy, anthocyanin and proline contents, antioxidant enzyme activity assay, in-gel enzyme activity staining, H2O2 and O2(-) content, immunoblotting, auxin and strigolactone content and proteomics analysis were evaluated at five different altitudes. P. saundersiana modulated the root architecture and leaf phenotype to enhance adaptation to alpine environmental stress through mechanisms that involved hormone synthesis and signal transduction, particularly the cross-talk between auxin and strigolactone. Furthermore, an increase of antioxidant proteins and primary metabolites as a response to the alpine environment in P. saundersiana was observed. Proteins associated with the epigenetic regulation of DNA stability and post-translational protein degradation was also involved in this process. Based on these findings, P. saundersiana uses multiple strategies to adapt to the high-altitude environment of the Alpine region.

BIOLOGICAL SIGNIFICANCE

The alpine environment, which is characterized by sharp temperature shifts, high levels of ultraviolet radiation exposure, and low oxygen content, limits plant growth and distribution. Alpine plants have evolved strategies to survive the extremely harsh conditions prevailing at high altitudes; however, the underlying mechanisms remain poorly understood. The alpine plant Potentilla saundersiana is widespread in the Northwestern Tibetan Plateau. Here we adopted a comparative proteomics approach to investigate the mechanisms by which P. saundersiana withstands the alpine environment by examining plants located at five different altitudes. We detected and functionally characterized 118 proteins spots with variable abundance. Proteins involved in antioxidant activity, primary metabolites, epigenetic regulation, and protein post-translational modification play important roles in conferring tolerance to alpine environments. Furthermore, our results indicate that P. saundersiana modulates the root architecture and leaf phenotype to enhance adaptation to alpine environmental stress. These results provide novel insight into the multiple strategies underlying P. saundersiana adaptation to the high-altitude environment of the Northwestern Tibetan Plateau.

Reproducibility and consistency of proteomic experiments on natural populations of a non-model aquatic insect

Population proteomics has a great potential to address evolutionary and ecological questions, but its use in wild populations of non-model organisms is hampered by uncontrolled sources of variation. Here we compare the response to temperature extremes of two geographically distant populations of a diving beetle species (Agabus ramblae) using 2-D DIGE. After one week of acclimation in the laboratory under standard conditions, a third of the specimens of each population were placed at either 4 or 27°C for 12 h, with another third left as a control. We then compared the protein expression level of three replicated samples of 2-3 specimens for each treatment. Within each population, variation between replicated samples of the same treatment was always lower than variation between treatments, except for some control samples that retained a wider range of expression levels. The two populations had a similar response, without significant differences in the number of protein spots over- or under-expressed in the pairwise comparisons between treatments. We identified exemplary proteins among those differently expressed between treatments, which proved to be proteins known to be related to thermal response or stress. Overall, our results indicate that specimens collected in the wild are suitable for proteomic analyses, as the additional sources of variation were not enough to mask the consistency and reproducibility of the response to the temperature treatments.

Comparative transcriptomics of the model mushroom Coprinopsis cinerea reveals tissue-specific armories and a conserved circuitry for sexual development

Background

It is well known that mushrooms produce defense proteins and secondary metabolites against predators and competitors; however, less is known about the correlation between the tissue-specific expression and the target organism (antagonist) specificity of these molecules. In addition, conserved transcriptional circuitries involved in developing sexual organs in fungi are not characterized, despite the growing number of gene expression datasets available from reproductive and vegetative tissue. The aims of this study were: first, to evaluate the tissue specificity of defense gene expression in the model mushroom Coprinopsis cinerea and, second, to assess the degree of conservation in transcriptional regulation during sexual development in basidiomycetes.

Results

In order to characterize the regulation in the expression of defense loci and the transcriptional circuitries controlling sexual reproduction in basidiomycetes, we sequenced the poly (A)-positive transcriptome of stage 1 primordia and vegetative mycelium of C. cinerea A43mutB43mut. Our data show that many genes encoding predicted and already characterized defense proteins are differentially expressed in these tissues. The predicted specificity of these proteins with regard to target organisms suggests that their expression pattern correlates with the type of antagonists these tissues are confronted with. Accordingly, we show that the stage 1 primordium-specific protein CC1G_11805 is toxic to insects and nematodes. Comparison of our data to analogous data from Laccaria bicolor and Schizophyllum commune revealed that the transcriptional regulation of nearly 70 loci is conserved and probably subjected to stabilizing selection. A Velvet domain-containing protein was found to be up-regulated in all three fungi, providing preliminary evidence of a possible role of the Velvet protein family in sexual development of basidiomycetes. The PBS-soluble proteome of C. cinerea primordia and mycelium was analyzed by shotgun LC-MS. This proteome data confirmed the presence of intracellular defense proteins in primordia.

Conclusions

This study shows that the exposure of different tissues in fungi to different types of antagonists shapes the expression pattern of defense loci in a tissue-specific manner. Furthermore, we identify a transcriptional circuitry conserved among basidiomycetes during fruiting body formation that involves, amongst other transcription factors, the up-regulation of a Velvet domain-containing protein.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-492) contains supplementary material, which is available to authorized users.

The salivary secretome of the biting midge, Culicoides sonorensis

Culicoides biting midges (Diptera: Ceratopogonidae) are hematophagous insects with over 1400 species distributed throughout the world. Many of these species are of particular agricultural importance as primary vectors of bluetongue and Schmallenberg viruses, yet little is known about Culicoides genomics and proteomics. Detailed studies of members from other blood-feeding Dipteran families, including those of mosquito (Culicidae) and black fly (Simuliidae), have shown that protein components within the insect's saliva facilitate the blood feeding process. To determine the protein components in Culicoides sonorensis midges, secreted saliva was collected for peptide sequencing by tandem mass spectrometry. Forty-five secreted proteins were identified, including members of the D7 odorant binding protein family, Kunitz-like serine protease inhibitors, maltase, trypsin, and six novel proteins unique to C. sonorensis. Identifying the complex myriad of proteins in saliva from blood-feeding Dipteran species is critical for understanding their role in blood feeding, arbovirus transmission, and possibly the resulting disease pathogenesis.

Litter size manipulation in laboratory mice: an example of how proteomic analysis can uncover new mechanisms underlying the cost of reproduction

Background

Life history theories predict that investment in current reproduction comes at a cost for future reproduction and survival. Oxidative stress is one of the best documented mechanisms underlying costs of reproduction to date. However, other, yet to be described molecular mechanisms that play a short term role during reproduction may explain the negative relationships underlying the cost of reproduction. To identify such new mechanisms, we used a global proteomic determination of liver protein profiles in laboratory adult female mice whose litter size had been either reduced or enlarged after birth. This litter size manipulation was expected to affect females by either raising or decreasing their current reproductive effort. At the same time, global parameters and levels of oxidative stress were also measured in all females.

Results

Based on plasma analyses, females with enlarged litters exhibited increased levels of oxidative stress at the date of weaning compared to females with reduced litters, while no significant difference was found between both the latter groups and control females. None of the liver proteins related to oxidative balance were significantly affected by the experimental treatment. In contrast, the liver protein profiles of females with enlarged and reduced litters suggested that calcium metabolism and cell growth regulation were negatively affected by changes in the number of pup reared.

Conclusions

Plasma oxidative stress levels in reproductive mice revealed that the degree of investment in reproduction can actually incur a cost in terms of plasmatic oxidative stress, their initial investment in reproduction being close to maximum and remaining at a same level when the energy demand of lactation is reduced. Liver proteomic profiles in reproductive females show that hepatic oxidative stress is unlikely to be involved in the cost of reproduction. Reproductive females rather exhibited liver protein profiles similar to those previously described in laboratory ageing mice, thus suggesting that hepatic cell pro-ageing processes may be involved in the cost of reproduction. Overall, our data illustrate how a proteomic approach can unravel new mechanisms sustaining life-history trade-offs, and reproduction costs in particular.

Proteomics to study adaptations in marine organisms to environmental stress

Comparisons of proteomic responses of closely related congeners and populations have shown which cellular processes are critical to adapt to environmental stress. For example, several proteomic species comparisons showed that increasing abundances of oxidative stress proteins indicate that reactive oxygen species (ROS) represent a ubiquitous signal and possible co-stressor of warm and cold temperature, acute hyposaline and low pH stress, possibly causing a shift from pro-oxidant NADH-producing to anti-oxidant NADPH-producing and -consuming metabolic pathways. Changes in cytoskeletal and actin-binding proteins in response to several stressors, including ROS, suggest that both are important structural and functional elements in responding to stress. Disruption of protein homeostasis, e.g., increased abundance of molecular chaperones, was severe in response to acute heat stress, inducing proteolysis, but was also observed in response to chronic heat and cold stress and was concentrated to the endoplasmic reticulum during hyposaline stress. Small GTPases affecting vesicle formation and transport, Ca(2+)-signaling and ion transport responded to salinity stress in species- and population-specific ways. Aerobic energy metabolism was in general down-regulated in response to temperature, hypoxia, hyposalinity and low pH stress, but other metabolic pathways were activated to respond to increased oxidative stress or to switch metabolic fuels. Thus, comparative proteomics is a powerful approach to identify functionally adaptive variation. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Annual reversible plasticity of feeding structures: cyclical changes of jaw allometry in a sea urchin

A wide variety of organisms show morphologically plastic responses to environmental stressors but in general these changes are not reversible. Though less common, reversible morphological structures are shown by a range of species in response to changes in predators, competitors or food. Theoretical analysis indicates that reversible plasticity increases fitness if organisms are long-lived relative to the frequency of changes in the stressor and morphological changes are rapid. Many sea urchin species show differences in the sizes of jaws (demi-pyramids) of the feeding apparatus, Aristotle's lantern, relative to overall body size, and these differences have been correlated with available food. The question addressed here is whether reversible changes of relative jaw size occur in the field as available food changes with season. Monthly samples of the North American Pacific coast sea urchin Strongylocentrotus purpuratus were collected from Gregory Point on the Oregon (USA) coast and showed an annual cycle of relative jaw size together with a linear trend from 2007 to 2009. Strongylocentrotus purpuratus is a long-lived species and under field conditions individuals experience multiple episodes of changes in food resources both seasonally and from year to year. Their rapid and reversible jaw plasticity fits well with theoretical expectations.

Looking for protein expression signatures in European eel peripheral blood mononuclear cells after in vivo exposure to perfluorooctane sulfonate and a real world field study

The decline of European eel population can be attributed to many factors such as pollution by xenobiotics present in domestic and industrial effluents. Perfluorooctane sulfonate (PFOS) is a ubiquitous compound of a particular concern in Europe. PFOS can reach high concentrations in tissues of organisms and many toxic effects have been reported in fish. This study aimed at evaluating the toxicological effects of PFOS in European eel peripheral blood mononuclear cells (PBMCs) at the protein expression level. To identify proteins whose expression was modified by PFOS, we performed a proteomic analysis on the post-nuclear fraction of PBMCs after a chronic exposure (28 days) of yellow eels to zero, 1 or 10 μg/L PFOS. This in vivo study was completed by a proteomic field study on eels sampled in Belgian rivers presenting different PFOS pollution degrees. Proteins were separated by two-dimensional in-gel electrophoresis (2D-DIGE) to compare the post-nuclear fraction of PBMCs from the reference group with cells from fish exposed to the pollutant of interest. On the 28 spots that were significantly (p < 0.05; ANOVA followed by a Dunnett post-hoc test) affected by PFOS in the in vivo experiment, a total of 17 different proteins were identified using nano-LC ESI-MS/MS and the Peptide and Protein Prophet of Scaffold software. In the field experiment, 18 significantly (p < 0.05; ANOVA followed by Dunnett's test) affected spots conducted to the identification of 16 different proteins. Interestingly, only three proteins were found in common between in vivo and in situ experiments: plastin-2, alpha-enolase and glyceraldehyde 3-phosphate dehydrogenase. Comparing the results with a previous study, plastin-2 and alpha-enolase were also been found to be affected after in vitro exposure of PBMCs during 48 h to either 10 μg or 1 mg PFOS/L. Potential use of these proteins as biomarkers of PFOS exposure in European eel could indicate early warning signals.

Integrating phenotypic plasticity within an Ecological Genomics framework: recent insights from the genomics, evolution, ecology, and fitness of plasticity

E.B. Ford's 1964 book Ecological Genetics was a call for biologists to engage in multidisciplinary work in order to elucidate the link between genotype, phenotype, and fitness for ecologically relevant traits. In this review, we argue that the integration of an ecological genomics framework in studies of phenotypic plasticity is a promising approach to elucidate the causal links between genes and the environment, particularly during colonization of novel environments, environmental change, and speciation. This review highlights some of the questions and hypotheses generated from a mechanistic, evolutionary, and ecological perspective, in order to direct the continued and future use of genomic tools in the study of phenotypic plasticity.

Ecological genomics of host behavior manipulation by parasites

Among the vast array of niche exploitation strategies exhibited by millions of different species on Earth, parasitic lifestyles are characterized by extremely successful evolutionary outcomes. Some parasites even seem to have the ability to 'control' their host's behavior to fulfill their own vital needs. Research efforts in the past decades have focused on surveying the phylogenetic diversity and ecological nature of these host-parasite interactions, and trying to understand their evolutionary significance. However, to understand the proximal and ultimate causes of these behavioral alterations triggered by parasitic infections, the underlying molecular mechanisms governing them must be uncovered. Studies using ecological genomics approaches have identified key candidate molecules involved in host-parasite molecular cross-talk, but also molecules not expected to alter behavior. These studies have shown the importance of following up with functional analyses, using a comparative approach and including a time-series analysis. High-throughput methods surveying different levels of biological information, such as the transcriptome and the epigenome, suggest that specific biologically-relevant processes are affected by infection, that sex-specific effects at the level of behavior are recapitulated at the level of transcription, and that epigenetic control represents a key factor in managing life cycle stages of the parasite through temporal regulation of gene expression. Post-translational processes, such as protein-protein interactions (interactome) and post translational modifications (e.g. protein phosphorylation, phosphorylome), and processes modifying gene expression and translation, such as interactions with microRNAs (microRNAome), are examples of promising avenues to explore to obtain crucial insights into the proximal and ultimate causes of these fascinating and complex inter-specific interactions.

Proteomic characterization of the hemolymph of Octopus vulgaris infected by the protozoan parasite Aggregata octopiana

The immune system of cephalopods is poorly known to date. The lack of genomic information makes difficult to understand vital processes like immune defense mechanisms and their interaction with pathogens at molecular level. The common octopus Octopus vulgaris has a high economic relevance and potential for aquaculture. However, disease outbreaks provoke serious reductions in production with potentially severe economic losses. In this study, a proteomic approach is used to analyze the immune response of O. vulgaris against the coccidia Aggregata octopiana, a gastrointestinal parasite which impairs the cephalopod nutritional status. The hemocytes and plasma proteomes were compared by 2-DE between sick and healthy octopus. The identities of 12 differentially expressed spots and other 27 spots without significant alteration from hemocytes, and 5 spots from plasma, were determined by mass spectrometry analysis aided by a six reading-frame translation of an octopus hemocyte RNA-seq database and also public databases. Principal component analysis pointed to 7 proteins from hemocytes as the major contributors to the overall difference between levels of infection and so could be considered as potential biomarkers. Particularly, filamin, fascin and peroxiredoxin are highlighted because of their implication in octopus immune defense activity. From the octopus plasma, hemocyanin was identified. This work represents a first step forward in order to characterize the protein profile of O. vulgaris hemolymph, providing important information for subsequent studies of the octopus immune system at molecular level and also to the understanding of the basis of octopus tolerance-resistance to A. octopiana.

BIOLOGICAL SIGNIFICANCE

The immune system of cephalopods is poorly known to date. The lack of genomic information makes difficult to understand vital processes like immune defense mechanisms and their interaction with pathogens at molecular level. The study herein presented is focused to the comprehension of the octopus immune defense against a parasite infection. Particularly, it is centered in the host-parasite relationship developed between the octopus and the protozoan A. octopiana, which induces severe gastrointestinal injuries in octopus that produce a malabsorption syndrome. The common octopus is a commercially important species with a high potential for aquaculture in semi-open systems, and this pathology reduces the condition of the octopus populations on-growing in open-water systems resulting in important economical loses. This is the first proteomic approach developed on this host-parasite relationship, and therefore, the contribution of this work goes from i) ecological, since this particular relationship is tending to be established as a model of host-parasite interaction in natural populations; ii) evolutionary, due to the characterization of immune molecules that could contribute to understand the functioning of the immune defense in these highly evolved mollusks; and iii) to economical view. The results of this study provide an overview of the octopus hemolymph proteome. Furthermore, proteins influenced by the level of infection and implicated in the octopus cellular response are also showed. Consequently, a set of biomarkers for disease resistance is suggested for further research that could be valuable for the improvement of the octopus culture, taken into account their high economical value, the declining of landings and the need for the diversification of reared species in order to ensure the growth of the aquaculture activity. Although cephalopods are model species for biomedical studies and possess potential in aquaculture, their genomes have not been sequenced yet, which limits the application of genomic data to research important biological processes. Similarly, the octopus proteome, like other non-model organisms, is poorly represented in public databases. Most of the proteins were identified from an octopus' hemocyte RNA-seq database that we have performed, which will be the object of another manuscript in preparation. Therefore, the need to increase molecular data from non-model organisms is herein highlighted. Particularly, here is encouraged to expand the knowledge of the genomic of cephalopods in order to increase successful protein identifications. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Proteomic profiling of developing cotton fibers from wild and domesticated Gossypium barbadense

Pima cotton (Gossypium barbadense) is widely cultivated because of its long, strong seed trichomes ('fibers') used for premium textiles. These agronomically advanced fibers were derived following domestication and thousands of years of human-mediated crop improvement. To gain an insight into fiber development and evolution, we conducted comparative proteomic and transcriptomic profiling of developing fiber from an elite cultivar and a wild accession. Analyses using isobaric tag for relative and absolute quantification (iTRAQ) LC-MS/MS technology identified 1317 proteins in fiber. Of these, 205 were differentially expressed across developmental stages, and 190 showed differential expression between wild and cultivated forms, 14.4% of the proteome sampled. Human selection may have shifted the timing of developmental modules, such that some occur earlier in domesticated than in wild cotton. A novel approach was used to detect possible biased expression of homoeologous copies of proteins. Results indicate a significant partitioning of duplicate gene expression at the protein level, but an approximately equal degree of bias for each of the two constituent genomes of allopolyploid cotton. Our results demonstrate the power of complementary transcriptomic and proteomic approaches for the study of the domestication process. They also provide a rich database for mining for functional analyses of cotton improvement or evolution.