The cuticle proteome of a planktonic crustacean.

The cuticles of arthropods provide an interface between the organism and its environment. Thus, the cuticle's structure influences how the organism responds to and interacts with its surroundings. Here, we used label-free quantification proteomics to provide a proteome of the moulted cuticle of the aquatic crustacean Daphnia magna, which has long been a prominent subject of studies on ecology, evolution, and developmental biology. We detected a total of 278 high-confidence proteins. Using protein sequence domain and functional enrichment analyses, we identified chitin-binding structural proteins and chitin-modifying enzymes as the most abundant protein groups in the cuticle proteome. Structural cuticular protein families showed a similar distribution to those found in other arthropods and indicated proteins responsible for the soft and flexible structure of the Daphnia cuticle. Finally, cuticle protein genes were also clustered as tandem gene arrays in the D. magna genome. The cuticle proteome presented here will be a valuable resource to the Daphnia research community, informing genome annotations and investigations on diverse topics such as the genetic basis of interactions with predators and parasites.

Interclonal proteomic responses to predator exposure in Daphnia magna may depend on predator composition of habitats.

Phenotypic plasticity, the ability of one genotype to express different phenotypes in response to changing environmental conditions, is one of the most common phenomena characterizing the living world and is not only relevant for the ecology but also for the evolution of species. Daphnia, the water flea, is a textbook example for predator-induced phenotypic plastic defences; however, the analysis of molecular mechanisms underlying these inducible defences is still in its early stages. We exposed Daphnia magna to chemical cues of the predator Triops cancriformis to identify key processes underlying plastic defensive trait formation. To get a more comprehensive idea of this phenomenon, we studied four genotypes with five biological replicates each, originating from habitats characterized by different predator composition, ranging from predator-free habitats to habitats containing T. cancriformis. We analysed the morphologies as well as proteomes of predator-exposed and control animals. Three genotypes showed morphological changes when the predator was present. Using a high-throughput proteomics approach, we found 294 proteins which were significantly altered in their abundance after predator exposure in a general or genotype-dependent manner. Proteins connected to genotype-dependent responses were related to the cuticle, protein synthesis and calcium binding, whereas the yolk protein vitellogenin increased in abundance in all genotypes, indicating their involvement in a more general response. Furthermore, genotype-dependent responses at the proteome level were most distinct for the only genotype that shares its habitat with Triops. Altogether, our study provides new insights concerning genotype-dependent and general molecular processes involved in predator-induced phenotypic plasticity in D. magna.

Stage-specific proteome signatures in early bovine embryo development.

Development of early embryonic stages before activation of the embryonic genome depends on sufficiently stored products of the maternal genome, adequate recruitment and degradation of mRNAs, as well as activation, deactivation, and relocation of proteins. By application of an isobaric tagging for relative and absolute quantification (iTRAQ)-based approach, the proteomes of bovine embryos at the zygote and 2-cell and 4-cell stage with MII oocytes as a reference were quantitatively analyzed. Of 1072 quantified proteins, 87 differed significantly in abundance between the four stages. The proteomes of 2-cell and 4-cell embryos differed most from the reference MII oocyte, and a considerable fraction of proteins continuously increased in abundance during the stages analyzed, despite a strongly attenuated rate of translation reported for this period. Bioinformatic analysis revealed particularly interesting proteins involved in the p53 pathway, lipid metabolism, and mitosis. Verification of iTRAQ results by targeted SRM (selected reaction monitoring) analysis revealed excellent agreement for all five proteins analyzed. By principal component analysis, SRM quantifications comprising a panel of only five proteins were shown to discriminate between all four developmental stages analyzed here. For future experiments, an expanded SRM protein panel will provide the potential to detect developmental disturbances with high sensitivity and enable first insights into the underlying molecular pathways.

Proteomic analysis of Daphnia magna hints at molecular pathways involved in defensive plastic responses.

BACKGROUND: Phenotypic plasticity in defensive traits occurs in many species when facing heterogeneous predator regimes. The waterflea Daphnia is well-known for showing a variety of these so called inducible defences. However, molecular mechanisms underlying this plasticity are poorly understood so far. We performed proteomic analysis on Daphnia magna exposed to chemical cues of the predator Triops cancriformis. D. magna develops an array of morphological changes in the presence of Triops including changes of carapace morphology and cuticle hardening. RESULTS: Using the 2D-DIGE technique, 1500 protein spots could be matched and quantified. We discovered 179 protein spots with altered intensity when comparing Triops exposed animals to a control group, and 69 spots were identified using nano-LC MS/MS. Kairomone exposure increased the intensity of spots containing muscle proteins, cuticle proteins and chitin-modifying enzymes as well as enzymes of carbohydrate and energy metabolism. The yolk precursor protein vitellogenin decreased in abundance in 41 of 43 spots. CONCLUSION: Identified proteins may be either directly involved in carapace stability or reflect changes in energy demand and allocation costs in animals exposed to predator kairomones. Our results present promising candidate proteins involved in the expression of inducible defences in Daphnia and enable further in depth analysis of this phenomenon.

A window into local adaptation.

How organisms adapt to their environment is not only a central topic of evolutionary biology but also a pressing question in the light of recent global change. Unravelling the genetic basis of these local adaptations can help to predict the response of a population to an increase in temperature or the more frequent occurrence of droughts. A popular approach to study the genes that drive local adaptation is the analysis of genotype-environment associations (GEA), testing the correlation of genomic features (typically single-nucleotide polymorphisms, SNPs) and environmental conditions. In this issue of Molecular Ecology Resources, Booker et al. (Molecular Ecology Resources, 2023) present a new approach to GEA, introducing genomic window analysis. They combine the information of neighbouring SNPs instead of analysing each SNP independently, therefore gaining power for detecting genomic signals of environmental adaptation. Using simulations of local adaptation to a heterogeneous environment as well as previously published real data from a natural population of lodgepole pine, they prove the superiority of their method over several established GEA approaches, especially in the case of small sample sizes. Leveraging the information present in closely linked genomic sites, Booker et al. (Molecular Ecology Resources, 2023) take genotype-environment association studies to the next level.

Evolution of Metabolome and Transcriptome Supports a Hierarchical Organization of Adaptive Traits.

Most organismal phenotypes have a polygenic basis, which enables adaptive phenotypic responses on ecological time scales. While adaptive phenotypic changes are highly parallel in replicate populations, this does not apply to the contributing loci. In particular for small populations, the same phenotypic shift can be fueled by different sets of alleles at alternative loci (genetic redundancy). Although this phenomenon is empirically well supported, the molecular basis of the genetic redundancy is not yet understood. To fill this gap, we compared the heterogeneity of the evolutionary transcriptomic and metabolomic response in ten Drosophila simulans populations which evolved parallel high-level phenotypic changes in a novel temperature environment but used different allelic combinations of alternative loci. We showed that the metabolome evolved more parallel than the transcriptome, confirming a hierarchical organization of molecular phenotypes. Different sets of genes responded in each evolved population but led to the enrichment of similar biological functions and a consistent metabolic profile. Since even the metabolomic response was still highly heterogeneous across evolved populations, we propose that selection may operate on pathways/networks.

Effects of larval crowding on the transcriptome of Drosophila simulans.

Larval crowding is one common ecological stressor for many insect species. In Drosophila, high larval density alters multiple widely-studied phenotypes including life-history traits, morphology and behavior. Nevertheless, we still miss a holistic view of the full range of phenotypic changes and the underlying molecular mechanisms. In this study, we analyzed the adult transcriptomes of high and low larval density fly cohorts, and highlighted the molecular basis of the plastic traits. Increased cellular energy metabolism and locomotion, along with reduced reproductive investment, are key responses to high larval density. Moreover, we compared the expression changes among cohorts with different developmental delays caused by larval crowding. The majority of genes induced by larval crowding showed the strongest expression alterations in cohorts with intermediate delay. Furthermore, linear expression changes were observed in genes related to nutrition and detoxification. Comparing different high-density cohorts could provide insights into the varied responses to distinct larval crowding-induced stresses such as space competition, food degradation and waste accumulation.

Detecting selected haplotype blocks in evolve and resequence experiments.

Shifting from the analysis of single nucleotide polymorphisms to the reconstruction of selected haplotypes greatly facilitates the interpretation of evolve and resequence (E&R) experiments. Merging highly correlated hitchhiker SNPs into haplotype blocks reduces thousands of candidates to few selected regions. Current methods of haplotype reconstruction from Pool-seq data need a variety of data-specific parameters that are typically defined ad hoc and require haplotype sequences for validation. Here, we introduce haplovalidate, a tool which detects selected haplotypes in Pool-seq time series data without the need for sequenced haplotypes. Haplovalidate makes data-driven choices of two key parameters for the clustering procedure, the minimum correlation between SNPs constituting a cluster and the window size. Applying haplovalidate to simulated E&R data reliably detects selected haplotype blocks with low false discovery rates. Importantly, our analyses identified a restriction of the haplotype block-based approach to describe the genomic architecture of adaptation. We detected a substantial fraction of haplotypes containing multiple selection targets. These blocks were considered as one region of selection and therefore led to underestimation of the number of selection targets. We demonstrate that the separate analysis of earlier time points can significantly increase the separation of selection targets into individual haplotype blocks. We conclude that the analysis of selected haplotype blocks has great potential for the characterization of the adaptive architecture with E&R experiments.

Phenotypic convergence in a natural Daphnia population acclimated to low temperature.

Fluidity of a given membrane decreases at lower ambient temperatures, whereas it rises at increasing temperatures, which is achieved through changes in membrane lipid composition. In consistence with homeoviscous adaptation theory, lower temperatures result in increased tissue concentrations of polyunsaturated fatty acids (PUFAs) in Daphnia magna, suggesting a higher PUFA requirement at lower temperatures. However, so far homeoviscous adaptation has been suggested for single or geographically separated Daphnia genotypes only. Here, we investigated changes in relative fatty acid (FA) tissue concentrations in response to a lower temperature (15°C) within a D. magna population. We determined juvenile growth rates (JGR) and FA patterns of 14 genotypes that were grown on Chlamydomonas klinobasis at 15°C and 20°C. We report significant differences of JGR and the relative body content of various FAs between genotypes at either temperature and between temperatures. Based on slopes of reaction norms, we found genotype-specific changes in FA profiles between temperatures suggesting that genotypes have different strategies to cope with changing temperatures. In a hierarchical clustering analysis, we grouped genotypes according to differences in direction and magnitude of changes in relative FA content, which resulted in three clusters of genotypes following different patterns of changes in FA composition. These patterns suggest a lower importance of the PUFA eicosapentaenoic acid (EPA, C20:5ω3) than previously assumed. We calculated an unsaturation index (UI) as a proxy for membrane fluidity at 15°C, and we neither found significant differences for this UI nor for fitness, measured as JGR, between the three genotype clusters. We conclude that these three genotype clusters represent different physiological solutions to temperature changes by altering the relative share of different FAs, but that their phenotypes converge with respect to membrane fluidity and JGR. These clusters will be subjected to different degrees of PUFA limitation when sharing the same diet.

The genetic architecture of temperature adaptation is shaped by population ancestry and not by selection regime.

BACKGROUND: Understanding the genetic architecture of temperature adaptation is key for characterizing and predicting the effect of climate change on natural populations. One particularly promising approach is Evolve and Resequence, which combines advantages of experimental evolution such as time series, replicate populations, and controlled environmental conditions, with whole genome sequencing. Recent analysis of replicate populations from two different Drosophila simulans founder populations, which were adapting to the same novel hot environment, uncovered very different architectures-either many selection targets with large heterogeneity among replicates or fewer selection targets with a consistent response among replicates. RESULTS: Here, we expose the founder population from Portugal to a cold temperature regime. Although almost no selection targets are shared between the hot and cold selection regime, the adaptive architecture was similar. We identify a moderate number of targets under strong selection (19 selection targets, mean selection coefficient = 0.072) and parallel responses in the cold evolved replicates. This similarity across different environments indicates that the adaptive architecture depends more on the ancestry of the founder population than the specific selection regime. CONCLUSIONS: These observations will have broad implications for the correct interpretation of the genomic responses to a changing climate in natural populations.

A 24 h Age Difference Causes Twice as Much Gene Expression Divergence as 100 Generations of Adaptation to a Novel Environment.

Gene expression profiling is one of the most reliable high-throughput phenotyping methods, allowing researchers to quantify the transcript abundance of expressed genes. Because many biotic and abiotic factors influence gene expression, it is recommended to control them as tightly as possible. Here, we show that a 24 h age difference of Drosophilasimulans females that were subjected to RNA sequencing (RNA-Seq) five and six days after eclosure resulted in more than 2000 differentially expressed genes. This is twice the number of genes that changed expression during 100 generations of evolution in a novel hot laboratory environment. Importantly, most of the genes differing in expression due to age introduce false positives or negatives if an adaptive gene expression analysis is not controlled for age. Our results indicate that tightly controlled experimental conditions, including precise developmental staging, are needed for reliable gene expression analyses, in particular in an evolutionary framework.

Polymorphism-aware protein databases - a prerequisite for an unbiased proteomic analysis of natural populations.

Recent technological advances have increased the throughput of proteomics, facilitating the characterization of molecular phenotypes on the population level, thus bearing the potential to complement transcriptomic analyses. Reference protein databases are crucial for the analysis and quantification, because only peptides in the protein database can be identified. Any peptide carrying an amino acid variant cannot be identified. Because most proteomic studies, even of natural populations, do not account for polymorphisms, we analysed the influence of variant peptides on quantitative proteomic analyses. We used transcriptomic and proteomic data of two Drosophila melanogaster genotypes and identified genotype-specific variants from RNA-seq data. We introduce a simple pipeline to include these variants in a polymorphism-aware protein database and compared the results to an unmodified reference database. The polymorphism-aware database not only identifies more peptides, but the quantitative values also changed when peptide variants were included. We conclude that proteomic quantification is likely to be biased, in particular for small genes, when polymorphisms are being ignored. Polymorphism-aware databases may be therefore a key step towards improved proteomic data analyses, especially for the analysis of pooled individuals and the comparison of population samples.

Progressive muscle proteome changes in a clinically relevant pig model of Duchenne muscular dystrophy.

Duchenne muscular dystrophy (DMD) is caused by genetic deficiency of dystrophin and characterized by massive structural and functional changes of skeletal muscle tissue, leading to terminal muscle failure. We recently generated a novel genetically engineered pig model reflecting pathological hallmarks of human DMD better than the widely used mdx mouse. To get insight into the hierarchy of molecular derangements during DMD progression, we performed a proteome analysis of biceps femoris muscle samples from 2-day-old and 3-month-old DMD and wild-type (WT) pigs. The extent of proteome changes in DMD vs. WT muscle increased markedly with age, reflecting progression of the pathological changes. In 3-month-old DMD muscle, proteins related to muscle repair such as vimentin, nestin, desmin and tenascin C were found to be increased, whereas a large number of respiratory chain proteins were decreased in abundance in DMD muscle, indicating serious disturbances in aerobic energy production and a reduction of functional muscle tissue. The combination of proteome data for fiber type specific myosin heavy chain proteins and immunohistochemistry showed preferential degeneration of fast-twitch fiber types in DMD muscle. The stage-specific proteome changes detected in this large animal model of clinically severe muscular dystrophy provide novel molecular readouts for future treatment trials.

The influence of simulated microgravity on the proteome of Daphnia magna.

BACKGROUND: The waterflea Daphnia is an interesting candidate for bioregenerative life support systems (BLSS). These animals are particularly promising because of their central role in the limnic food web and its mode of reproduction. However, the response of Daphnia to altered gravity conditions has to be investigated, especially on the molecular level, to evaluate the suitability of Daphnia for BLSS in space. METHODS: In this study, we applied a proteomic approach to identify key proteins and pathways involved in the response of Daphnia to simulated microgravity generated by a two-dimensional (2D) clinostat. We analyzed five biological replicates using 2D-difference gel electrophoresis proteomic analysis. RESULTS: We identified 109 protein spots differing in intensity (P<0.05). Substantial fractions of these proteins are involved in actin microfilament organization, indicating the disruption of cytoskeletal structures during clinorotation. Furthermore, proteins involved in protein folding were identified, suggesting altered gravity induced breakdown of protein structures in general. In addition, simulated microgravity increased the abundance of energy metabolism-related proteins, indicating an enhanced energy demand of Daphnia. CONCLUSIONS: The affected biological processes were also described in other studies using different organisms and systems either aiming to simulate microgravity conditions or providing real microgravity conditions. Moreover, most of the Daphnia protein sequences are well-conserved throughout taxa, indicating that the response to altered gravity conditions in Daphnia follows a general concept. Data are available via ProteomeXchange with identifier PXD002096.

Proteomic analysis of extracellular medium of cryopreserved carp (Cyprinus carpio L.) semen.

During freezing and thawing, spermatozoa are exposed to physical and chemical stressors that result in adverse changes in sperm structures and physiological functions. The present study provides, for the first time, a comprehensive description of protein changes in the extracellular medium of cryopreserved semen. Using 2D-DIGE and a combination of protein fractionation by one-dimensional gel electrophoresis and high performance liquid chromatography electrospray ionization tandem mass spectrometry, 183 proteins released from sperm to an extracellular medium were identified. The majority of released proteins were involved in metabolism and energy production. Moreover, proteins associated with a response to stress, apoptosis, small GTPase mediated signal transduction, transcription, translation, protein folding and turnover, reproduction and DNA repair were identified. The dominant group of released proteins was related to cytoplasm. Moreover, specific proteins associated with the membrane, mitochondria and nucleus were identified. The identification of a high number of proteins released from sperm provides new insight into the mechanism of cryodamage to the particular sperm structure and to specific metabolic pathways, which were affected by cryopreservation. The availability of a catalog of carp sperm proteins altered by cryopreservation provides a crucial tool for the development of novel potential biomarkers of cryoinjuries and for the improvement of a long-term sperm preservation procedure.