Proteomics in evolutionary ecology: linking the genotype with the phenotype

Selection Across the Three-Dimensional Structure of Venom Proteins from North American Scolopendromorph Centipedes

Gene duplication followed by nucleotide differentiation is one of the simplest mechanisms to develop new functions for genes. However, the evolutionary processes underlying the divergence of multigene families remain controversial. We used multigene families found within the diversity of toxic proteins in centipede venom to test two hypotheses related to venom evolution: the two-speed mode of venom evolution and the rapid accumulation of variation in exposed residues (RAVER) model. The two-speed mode of venom evolution proposes that different types of selection impact ancient and younger venomous lineages with negative selection being the predominant form in ancient lineages and positive selection being the dominant form in younger lineages. The RAVER hypothesis proposes that, instead of different types of selection acting on different ages of venomous lineages, the different types of selection will selectively contribute to amino acid variation based on whether the residue is exposed to the solvent where it can potentially interact directly with toxin targets. This hypothesis parallels the longstanding understanding of protein evolution that suggests that residues found within the structural or active regions of the protein will be under negative or purifying selection, and residues that do not form part of these areas will be more prone to positive selection. To test these two hypotheses, we compared the venom of 26 centipedes from the order Scolopendromorpha from six currently recognized species from across North America using both transcriptomics and proteomics. We first estimated their phylogenetic relationships and uncovered paraphyly among the genus Scolopendra and evidence for cryptic diversity among currently recognized species. Using our phylogeny, we then characterized the diverse venom components from across the identified clades using a combination of transcriptomics and proteomics. We conducted selection-based analyses in the context of predicted three-dimensional properties of the venom proteins and found support for both hypotheses. Consistent with the two-speed hypothesis, we found a prevalence of negative selection across all proteins. Consistent with the RAVER hypothesis, we found evidence of positive selection on solvent-exposed residues, with structural and less-exposed residues showing stronger signal for negative selection. Through the use of phylogenetics, transcriptomics, proteomics, and selection-based analyses, we were able to describe the evolution of venom from an ancient venomous lineage and support principles of protein evolution that directly relate to multigene family evolution.

From Genotype to Phenotype: Raman Spectroscopy and Machine Learning for Label-Free Single-Cell Analysis

Raman spectroscopy has made significant progress in biosensing and clinical research. Here, we describe how surface-enhanced Raman spectroscopy (SERS) assisted with machine learning (ML) can expand its capabilities to enable interpretable insights into the transcriptome, proteome, and metabolome at the single-cell level. We first review how advances in nanophotonics-including plasmonics, metamaterials, and metasurfaces-enhance Raman scattering for rapid, strong label-free spectroscopy. We then discuss ML approaches for precise and interpretable spectral analysis, including neural networks, perturbation and gradient algorithms, and transfer learning. We provide illustrative examples of single-cell Raman phenotyping using nanophotonics and ML, including bacterial antibiotic susceptibility predictions, stem cell expression profiles, cancer diagnostics, and immunotherapy efficacy and toxicity predictions. Lastly, we discuss exciting prospects for the future of single-cell Raman spectroscopy, including Raman instrumentation, self-driving laboratories, Raman data banks, and machine learning for uncovering biological insights.

Omics-based molecular classifications empowering in precision oncology

BACKGROUND

In the past decades, cancer enigmatical heterogeneity at distinct expression levels could interpret disparities in therapeutic response and prognosis. It built hindrances to precision medicine, a tactic to tailor customized treatment informed by the tumors' molecular profile. Single-omics analysis dissected the biological features associated with carcinogenesis to some extent but still failed to revolutionize cancer treatment as expected. Integrated omics analysis incorporated tumor biological networks from diverse layers and deciphered a holistic overview of cancer behaviors, yielding precise molecular classification to facilitate the evolution and refinement of precision medicine.

CONCLUSION

This review outlined the biomarkers at multiple expression layers to tutor molecular classification and pinpoint tumor diagnosis, and explored the paradigm shift in precision therapy: from single- to multi-omics-based subtyping to optimize therapeutic regimens. Ultimately, we firmly believe that by parsing molecular characteristics, omics-based typing will be a powerful assistant for precision oncology.

Proteomic and physiological signatures of altitude adaptation in a Myrsine coriacea population under common garden conditions

Plants exhibit phenotypic plasticity in response to environmental variations, which can lead to stable genetic and physiological adaptations if exposure to specific conditions is prolonged. Myrsine coriacea demonstrates this through its ability to thrive in diverse environments. The objective of the article is to investigate potential differences in protein accumulation and physiological responses of M. coriacea by cultivating plants from seeds collected from four populations at different altitudes in a common garden experiment. Additionally, we aim to evaluate whether these differences exhibit genetic fixation. Through integrated physiological and proteomic analyses, we identified 170 differentially accumulated proteins and observed significant physiological differences among the populations. The high-altitude population (POP1) exhibited a unique proteomic profile with significant down-regulation of proteins involved in carbon fixation and energy metabolism, suggesting a potential reduction in photosynthetic efficiency. Physiological analyses showed lower leaf nitrogen content, net CO2 assimilation rate, specific leaf area, and relative growth rate in stem height for POP1, alongside higher leaf carbon isotopic composition (δ13C) and leaf carbon (C) content. These findings provide insight into the complex interplay between proteomic and physiological adaptations in M. coriacea and underscore the importance of local adaptations. SIGNIFICANCE: We investigate the adaptive responses of M. coriacea, a shrub with a broad phenotypic range, by cultivating plants from seeds collected at four different altitudes in a common garden experiment. These findings provide insight into the complex interplay between proteomic and physiological adaptations in M. coriacea and underscore the importance of local adaptations in the face of climate change. This study contributes to advancing our understanding of the influence of altitude-specific selection pressures on the molecular biology and physiology of plants in natural populations. Our findings provide valuable insights that enhance our ability to predict and comprehend how plants respond to climate change.

Functional convergence in gastric lysozymes of foregut-fermenting rodents, ruminants, and primates is not attributed to convergent molecular evolution

Convergent evolution is a widespread phenomenon. While there are many examples of convergent evolution at the phenotypic scale, convergence at the molecular level has been more difficult to identify. A classic example of convergent evolution across scales is that of the digestive lysozyme found in ruminants and Colobine monkeys. These herbivorous species rely on foregut fermentation, which has evolved to function more optimally under acidic conditions. Here, we explored if rodents with similar dietary strategies and digestive morphologies have convergently evolved a lysozyme with digestive functions. At the phenotypic level, we find that rodents with bilocular stomach morphologies exhibited a lysozyme that maintained higher relative activities at low pH values, similar to the lysozymes of ruminants and Colobine monkeys. Additionally, the lysozyme of Peromyscus leucopus shared a similar predicted protonation state as that observed in previously identified digestive lysozymes. However, we found limited evidence of positive selection acting on the lysozyme gene in foregut-fermenting species and did not identify patterns of convergent molecular evolution in this gene. This study emphasizes that phenotypic convergence need not be the result of convergent genetic modifications, and we encourage further exploration into the mechanisms regulating convergence across biological scales.

Recent omics advances in hair aging biology and hair biomarkers analysis

Aging is a complex natural process that leads to a decline in physiological functions, which is visible in signs such as hair graying, thinning, and loss. Although hair graying is characterized by a loss of pigment in the hair shaft, the underlying mechanism of age-associated hair graying is not fully understood. Hair graying and loss can have a significant impact on an individual's self-esteem and self-confidence, potentially leading to mental health problems such as depression and anxiety. Omics technologies, which have applications beyond clinical medicine, have led to the discovery of candidate hair biomarkers and may provide insight into the complex biology of hair aging and identify targets for effective therapies. This review provides an up-to-date overview of recent omics discoveries, including age-associated alterations of proteins and metabolites in the hair shaft and follicle, and highlights the significance of hair aging and graying biomarker discoveries. The decline in hair follicle stem cell activity with aging decreased the regeneration capacity of hair follicles. Cellular senescence, oxidative damage and altered extracellular matrix of hair follicle constituents characterized hair follicle and hair shaft aging and graying. The review attempts to correlate the impact of endogenous and exogenous factors on hair aging. We close by discussing the main challenges and limitations of the field, defining major open questions and offering an outlook for future research.

Biomarkers of aging

Aging biomarkers are a combination of biological parameters to (i) assess age-related changes, (ii) track the physiological aging process, and (iii) predict the transition into a pathological status. Although a broad spectrum of aging biomarkers has been developed, their potential uses and limitations remain poorly characterized. An immediate goal of biomarkers is to help us answer the following three fundamental questions in aging research: How old are we? Why do we get old? And how can we age slower? This review aims to address this need. Here, we summarize our current knowledge of biomarkers developed for cellular, organ, and organismal levels of aging, comprising six pillars: physiological characteristics, medical imaging, histological features, cellular alterations, molecular changes, and secretory factors. To fulfill all these requisites, we propose that aging biomarkers should qualify for being specific, systemic, and clinically relevant.

Artificial intelligence-based multi-omics analysis fuels cancer precision medicine

With biotechnological advancements, innovative omics technologies are constantly emerging that have enabled researchers to access multi-layer information from the genome, epigenome, transcriptome, proteome, metabolome, and more. A wealth of omics technologies, including bulk and single-cell omics approaches, have empowered to characterize different molecular layers at unprecedented scale and resolution, providing a holistic view of tumor behavior. Multi-omics analysis allows systematic interrogation of various molecular information at each biological layer while posing tricky challenges regarding how to extract valuable insights from the exponentially increasing amount of multi-omics data. Therefore, efficient algorithms are needed to reduce the dimensionality of the data while simultaneously dissecting the mysteries behind the complex biological processes of cancer. Artificial intelligence has demonstrated the ability to analyze complementary multi-modal data streams within the oncology realm. The coincident development of multi-omics technologies and artificial intelligence algorithms has fuelled the development of cancer precision medicine. Here, we present state-of-the-art omics technologies and outline a roadmap of multi-omics integration analysis using an artificial intelligence strategy. The advances made using artificial intelligence-based multi-omics approaches are described, especially concerning early cancer screening, diagnosis, response assessment, and prognosis prediction. Finally, we discuss the challenges faced in multi-omics analysis, along with tentative future trends in this field. With the increasing application of artificial intelligence in multi-omics analysis, we anticipate a shifting paradigm in precision medicine becoming driven by artificial intelligence-based multi-omics technologies.

The significance of cephalopod beaks as a research tool: An update

The use of cephalopod beaks in ecological and population dynamics studies has allowed major advances of our knowledge on the role of cephalopods in marine ecosystems in the last 60 years. Since the 1960's, with the pioneering research by Malcolm Clarke and colleagues, cephalopod beaks (also named jaws or mandibles) have been described to species level and their measurements have been shown to be related to cephalopod body size and mass, which permitted important information to be obtained on numerous biological and ecological aspects of cephalopods in marine ecosystems. In the last decade, a range of new techniques has been applied to cephalopod beaks, permitting new kinds of insight into cephalopod biology and ecology. The workshop on cephalopod beaks of the Cephalopod International Advisory Council Conference (Sesimbra, Portugal) in 2022 aimed to review the most recent scientific developments in this field and to identify future challenges, particularly in relation to taxonomy, age, growth, chemical composition (i.e., DNA, proteomics, stable isotopes, trace elements) and physical (i.e., structural) analyses. In terms of taxonomy, new techniques (e.g., 3D geometric morphometrics) for identifying cephalopods from their beaks are being developed with promising results, although the need for experts and reference collections of cephalopod beaks will continue. The use of beak microstructure for age and growth studies has been validated. Stable isotope analyses on beaks have proven to be an excellent technique to get valuable information on the ecology of cephalopods (namely habitat and trophic position). Trace element analyses is also possible using beaks, where concentrations are significantly lower than in other tissues (e.g., muscle, digestive gland, gills). Extracting DNA from beaks was only possible in one study so far. Protein analyses can also be made using cephalopod beaks. Future challenges in research using cephalopod beaks are also discussed.

Mechanism exploration of Solanum nigrum L. hyperaccumulating Cd compared to Zn from the perspective of metabolic pathways based on differentially expressed proteins using iTRAQ

It is challenging to determine the mechanism involved in only Cd hyperaccumulation by Solanum nigrum L. owing to the uniqueness of the process. Isobaric tags for relative and absolute quantitation (iTRAQ) were used to explore the mechanism by which S. nigrum hyperaccumulates Cd by comparing the differentially expressed proteins (DEPs) for Cd and Zn accumulation (non-Zn hyperaccumulator). Based on the comparison between the DEPs associated with Cd and Zn accumulation, the relative metabolic pathways reflected by 17 co-intersecting specific proteins associated with Cd and Zn accumulation included phagosome, aminoacyl-tRNA biosynthesis, and carbon metabolism. Apart from the 17 co-intersecting specific proteins, the conjoint metabolic pathways reported by 21 co-intersecting specific proteins associated with Cd accumulation and 30 co-intersecting specific proteins associated with Zn accumulation, the most differentially expressed metabolic pathways might cause Cd TF (Translocation factor)> 1 and Zn TF< 1, including protein export, ribosome, amino sugar, and nucleotide sugar metabolism. The determined DEPs were verified using qRT-PCR with the four key proteins M1CW30, A0A3Q7H652, A0A0V0IFB9, and A0A0V0IAC4. The plasma membrane H⁺-ATPase protein was identified using western blotting. Some physiological indices for protein-related differences indirectly confirmed the above results. These results are crucial to further explore the mechanisms involved in Cd hyperaccumulation.

Genetic characterization of potential venom resistance proteins in California ground squirrels ( Otospermophilus beecheyi ) using transcriptome analyses

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

KATP channel dependent heart multiome atlas

Plasmalemmal ATP sensitive potassium (K_ATP) channels are recognized metabolic sensors, yet their cellular reach is less well understood. Here, transgenic Kir6.2 null hearts devoid of the K_ATP channel pore underwent multiomics surveillance and systems interrogation versus wildtype counterparts. Despite maintained organ performance, the knockout proteome deviated beyond a discrete loss of constitutive K_ATP channel subunits. Multidimensional nano-flow liquid chromatography tandem mass spectrometry resolved 111 differentially expressed proteins and their expanded network neighborhood, dominated by metabolic process engagement. Independent multimodal chemometric gas and liquid chromatography mass spectrometry unveiled differential expression of over one quarter of measured metabolites discriminating the Kir6.2 deficient heart metabolome. Supervised class analogy ranking and unsupervised enrichment analysis prioritized nicotinamide adenine dinucleotide (NAD⁺), affirmed by extensive overrepresentation of NAD⁺ associated circuitry. The remodeled metabolome and proteome revealed functional convergence and an integrated signature of disease susceptibility. Deciphered cardiac patterns were traceable in the corresponding plasma metabolome, with tissue concordant plasma changes offering surrogate metabolite markers of myocardial latent vulnerability. Thus, Kir6.2 deficit precipitates multiome reorganization, mapping a comprehensive atlas of the K_ATP channel dependent landscape.

Integrated Multi-Omics for Novel Aging Biomarkers and Antiaging Targets

Aging is closely related to the occurrence of human diseases; however, its exact biological mechanism is unclear. Advancements in high-throughput technology provide new opportunities for omics research to understand the pathological process of various complex human diseases. However, single-omics technologies only provide limited insights into the biological mechanisms of diseases. DNA, RNA, protein, metabolites, and microorganisms usually play complementary roles and perform certain biological functions together. In this review, we summarize multi-omics methods based on the most relevant biomarkers in single-omics to better understand molecular functions and disease causes. The integration of multi-omics technologies can systematically reveal the interactions among aging molecules from a multidimensional perspective. Our review provides new insights regarding the discovery of aging biomarkers, mechanism of aging, and identification of novel antiaging targets. Overall, data from genomics, transcriptomics, proteomics, metabolomics, integromics, microbiomics, and systems biology contribute to the identification of new candidate biomarkers for aging and novel targets for antiaging interventions.

Shotgun analysis to identify differences in protein expression between granulocytes and hyalinocytes of the European flat oyster Ostrea edulis

Recovery of wild populations of the European flat oyster Ostrea edulis is important for ecosystem health and conservation of this species, because native oyster populations have dramatically declined or disappeared in most European waters. Diseases have contributed to oyster decline and are important constrains for oyster recovery. Understanding oyster immune system should contribute to design effective strategies to fight oyster diseases. Haemocytes play a pivotal role in mollusc immune responses protecting from infection. Two main types of haemocytes, granulocytes and hyalinocytes, are distinguished in O. edulis. A study aiming to explore differential functions between both haemocyte types and, thus, to enrich the knowledge of Ostrea edulis immune system, was performed by comparing the proteome of the two haemolymph cell types, using a shotgun approach through liquid chromatography (LC) coupled to mass spectrometry (MS). Cells from oyster haemolymph were differentially separated by Percoll density gradient centrifugation. Shotgun LC-MS/MS performance allowed the identification of 145 proteins in hyalinocytes and 138 in the proteome of granulocytes. After a comparative analysis, 55 proteins with main roles in defence were identified, from which 28 were representative of granulocytes and 27 of hyalinocytes, plus 11 proteins shared by both cell types. Different proteins involved in signal transduction, apoptosis, oxidative response, processes related with the cytoskeleton and structure, recognition and wound healing were identified as representatives of each haemocyte type. Important signalling pathways in the immune response such as MAPK, Ras and NF-κβ seemed to be more relevant for granulocytes, while the Wnt signalling pathway, particularly relevant for wound healing, more relevant in hyalinocytes. The differences in proteins involved in recognition and in cytoskeleton and structure suggest differential specialisation in processes of phagocytosis and internalisation of pathogens between haemocyte types. Apoptosis seemed more active in granulocytes. The differences in proteins involved in oxidative response also suggest different redox processes in each cell type.

Leveraging Experimental Strategies to Capture Different Dimensions of Microbial Interactions

Microorganisms are a fundamental part of virtually every ecosystem on earth. Understanding how collectively they interact, assemble, and function as communities has become a prevalent topic both in fundamental and applied research. Owing to multiple advances in technology, answering questions at the microbial system or network level is now within our grasp. To map and characterize microbial interaction networks, numerous computational approaches have been developed; however, experimentally validating microbial interactions is no trivial task. Microbial interactions are context-dependent, and their complex nature can result in an array of outcomes, not only in terms of fitness or growth, but also in other relevant functions and phenotypes. Thus, approaches to experimentally capture microbial interactions involve a combination of culture methods and phenotypic or functional characterization methods. Here, through our perspective of food microbiologists, we highlight the breadth of innovative and promising experimental strategies for their potential to capture the different dimensions of microbial interactions and their high-throughput application to answer the question; are microbial interaction patterns or network architecture similar along different contextual scales? We further discuss the experimental approaches used to build various types of networks and study their architecture in the context of cell biology and how they translate at the level of microbial ecosystem.

Host gene expression in wildlife disease: making sense of species-level responses

Emerging infectious diseases are significant threats to wildlife conservation, yet the impacts of pathogen exposure and infection can vary widely among host species. As such, conservation biologists and disease ecologists have increasingly aimed to understand species-specific host susceptibility using molecular methods. In particular, comparative gene expression assays have been used to contrast the transcriptomic responses of disease-resistant and disease-susceptible hosts to pathogen exposure. This work usually assumes that the gene expression responses of disease-resistant species will reveal the activation of molecular pathways contributing to host defence. However, results often show that disease-resistant hosts undergo little gene expression change following pathogen challenge. Here, we discuss the mechanistic implications of these "null" findings and offer methodological suggestions for future molecular studies of wildlife disease. First, we highlight that muted transcriptomic responses with minimal immune system recruitment may indeed be protective for nonsusceptible hosts if they limit immunopathology and promote pathogen tolerance in systems where susceptible hosts suffer from genetic dysregulation. Second, we argue that overly narrow investigation of responses to pathogen exposure may overlook important, constitutively active molecular pathways that underlie species-specific defences. Finally, we outline alternative study designs and approaches that complement interspecific transcriptomic comparisons, including intraspecific gene expression studies and genomic methods to detect signatures of selection. Collectively, these insights will help ecologists extract maximal information from conservation-relevant transcriptomic data sets, leading to a deeper understanding of host defences and, ultimately, the implementation of successful conservation interventions.

Proteomic analysis of F1 hybrids and intermediate variants in a littorina saxatilis hybrid zone

Proteomic analysis was carried out on the Crab (upper-shore) and Wave (lower-shore) ecotypes of Littorina saxatilis from a hybrid zone at Silleiro Cape, Spain. Proteome profiles of individual snails were obtained. Protein expression in F₁ hybrid snails bred in the laboratory and snails with intermediate shell phenotypes collected from the mid-shore were compared with Crab and Wave ecotypes using analytical approaches used to study dominance. Multivariate analysis over many protein spots showed that the F₁ snails are distinct from both ecotypes but closer to the Wave ecotype. The intermediate snails are highly variable, some closer to the Crab and others to the Wave ecotype. Considered on a protein by protein basis, some proteins are significantly closer in expression to the Crab and others to the Wave ecotype for both F₁ and intermediate snails. Furthermore, a significant majority of proteins were closer in expression to the Wave ecotype for the F₁, consistent with the multivariate analysis. No such significant majority toward either the Crab or Wave ecotype was observed for the intermediate snails. The closer similarity of F₁ and Wave ecotype expression patterns could be the result of similar selective pressures in the similar mid-shore and low-shore environments. For a significantly larger number of proteins, intermediate snails were closer in expression to the ecotype having the lower expression, for both Crab and Wave ecotypes. This is somewhat unexpected as lower expression might be expected to be an indication of impairment of function and lower fitness. Proteomic analysis could be important for the identification of candidate proteins useful for gaining improved understanding of adaptation and barriers to gene flow in hybrid zones.

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.

Kawasaki Disease Patient Stratification and Pathway Analysis Based on Host Transcriptomic and Proteomic Profiles

The aetiology of Kawasaki disease (KD), an acute inflammatory disorder of childhood, remains unknown despite various triggers of KD having been proposed. Host 'omic profiles offer insights into the host response to infection and inflammation, with the interrogation of multiple 'omic levels in parallel providing a more comprehensive picture. We used differential abundance analysis, pathway analysis, clustering, and classification techniques to explore whether the host response in KD is more similar to the response to bacterial or viral infections at the transcriptomic and proteomic levels through comparison of 'omic profiles from children with KD to those with bacterial and viral infections. Pathways activated in patients with KD included those involved in anti-viral and anti-bacterial responses. Unsupervised clustering showed that the majority of KD patients clustered with bacterial patients on both 'omic levels, whilst application of diagnostic signatures specific for bacterial and viral infections revealed that many transcriptomic KD samples had low probabilities of having bacterial or viral infections, suggesting that KD may be triggered by a different process not typical of either common bacterial or viral infections. Clustering based on the transcriptomic and proteomic responses during KD revealed three clusters of KD patients on both 'omic levels, suggesting heterogeneity within the inflammatory response during KD. The observed heterogeneity may reflect differences in the host response to a common trigger, or variation dependent on different triggers of the condition.

Transcriptome Analysis of the Japanese Pine Sawyer Beetle, Monochamus alternatus, Infected with the Entomopathogenic Fungus Metarhizium anisopliae JEF-197

The Japanese pine sawyer (JPS) beetle, Monochamus alternatus Hope (Coleoptera: Cerambycidae), damages pine trees and transmits the pine wilt nematode, Bursaphelenchus xylophilus Nickle. Chemical agents have been used to control JPS beetle, but due to various issues, efforts are being made to replace these chemical agents with entomopathogenic fungi. We investigated the expression of immune-related genes in JPS beetle in response to infection with JEF-197, a Metarhizium anisopliae isolate, using RNA-seq. RNA samples were obtained from JEF-197, JPS adults treated with JEF-197, and non-treated JPS adults on the 8th day after fungal treatment, and RNA-seq was performed using Illumina sequencing. JPS beetle transcriptome was assembled de novo and differentially expressed gene (DEG) analysis was performed. There were 719 and 1953 up- and downregulated unigenes upon JEF-197 infection, respectively. Upregulated contigs included genes involved in RNA transport, ribosome biogenesis in eukaryotes, spliceosome-related genes, and genes involved in immune-related signaling pathways such as the Toll and Imd pathways. Forty-two fungal DEGs related to energy and protein metabolism were upregulated, and genes involved in the stress response were also upregulated in the infected JPS beetles. Together, our results indicate that infection of JPS beetles by JEF-197 induces the expression of immune-related genes.

Genome-wide DNA methylation of the liver reveals delayed effects of early-life exposure to 17-α-ethinylestradiol in the self-fertilizing mangrove rivulus

Organisms exposed to endocrine disruptors in early life can show altered phenotype later in adulthood. Although the mechanisms underlying these long-term effects remain poorly understood, an increasing body of evidence points towards the potential role of epigenetic processes. In the present study, we exposed hatchlings of an isogenic lineage of the self-fertilizing fish mangrove rivulus for 28 days to 4 and 120 ng/L of 17-α-ethinylestradiol. After a recovery period of 140 days, reduced representation bisulphite sequencing (RRBS) was performed on the liver in order to assess the hepatic genome-wide methylation landscape. Across all treatment comparisons, a total of 146 differentially methylated fragments (DMFs) were reported, mostly for the group exposed to 4 ng/L, suggesting a non-monotonic effect of EE2 exposure. Gene ontology analysis revealed networks involved in lipid metabolism, cellular processes, connective tissue function, molecular transport and inflammation. The highest effect was reported for nipped-B-like protein B (NIPBL) promoter region after exposure to 4 ng/L EE2 (+ 21.9%), suggesting that NIPBL could be an important regulator for long-term effects of EE2. Our results also suggest a significant role of DNA methylation in intergenic regions and potentially in transposable elements. These results support the ability of early exposure to endocrine disruptors of inducing epigenetic alterations during adulthood, providing plausible mechanistic explanations for long-term phenotypic alteration. Additionally, this work demonstrates the usefulness of isogenic lineages of the self-fertilizing mangrove rivulus to better understand the biological significance of long-term alterations of DNA methylation by diminishing the confounding factor of genetic variability.

The physiological and molecular mechanisms of N transfer in Eucalyptus and Dalbergia odorifera intercropping systems using root proteomics

BACKGROUND

The mixing of Eucalyptus with N₂-fixing trees species (NFTs) is a frequently successful and sustainable cropping practice. In this study, we evaluated nitrogen (N) transfer and conducted a proteomic analysis of the seedlings of Eucalyptus urophylla × E. grandis (Eucalyptus) and an NFT, Dalbergia (D.) odorifera, from intercropping and monocropping systems to elucidate the physiological effects and molecular mechanisms of N transfer in mixed Eucalyptus and D. odorifera systems.

RESULTS

N transfer occurred from D. odorifera to Eucalyptus at a rate of 14.61% in the intercropping system, which increased N uptake and growth in Eucalyptus but inhibited growth in D. odorifera. There were 285 and 288 differentially expressed proteins by greater than 1.5-fold in Eucalyptus and D. odorifera roots with intercropping vs monoculture, respectively. Introduction of D. odorifera increased the stress resistance ability of Eucalyptus, while D. odorifera stress resistance was increased by increasing levels of jasmonic acid (JA). Additionally, the differentially expressed proteins of N metabolism, such as glutamine synthetase nodule isozyme (GS), were upregulated to enhance N competition in Eucalyptus. Importantly, more proteins were involved in synthetic pathways than in metabolic pathways in Eucalyptus because of the benefit of N transfer, and the two groups of N compound transporters were found in Eucalyptus; however, more functional proteins were involved in metabolic degradation in D. odorifera; specifically, the molecular mechanism of the transfer of N from D. odorifera to Eucalyptus was explained by proteomics.

CONCLUSIONS

Our study suggests that N transfer occurred from D. odorifera to Eucalyptus and was affected by the variations in the differentially expressed proteins. We anticipate that these results can be verified in field experiments for the sustainable development of Eucalyptus plantations.

A Primer and Guidelines for Shotgun Proteomic Analysis in Non-model Organisms

During the last decade, we have witnessed outstanding advances in proteomics led mostly by great technological improvements in mass spectrometry field allowing high-throughput production of high-quality data used for massive protein identification and quantification. From a practical viewpoint, these advances have been mainly exploited in research projects involving model organisms with abundant genomic and proteomic information available in public databases. However, there is a growing number of organisms of high interest in different disciplines, such as ecological, biotechnological, and evolutionary research, yet poorly represented in these databases. Important advances in massive parallel sequencing technology and easy accessibility of this technology to many research laboratories have made nowadays possible to produce customized genomic and proteomic databases of any organism. Along this line, the use of proteogenomic approaches by combining in the same analysis the data obtained from different omic levels has emerged as a very useful and powerful strategy to run shotgun proteomic experiments specially focused on non-model organisms. In this chapter, we provide detailed procedures to undertake shotgun quantitative proteomic experiments following either a label-free or an isobaric labeling approach in non-model organisms, emphasizing also a few key aspects related to experimental design and data analysis.

Proteogenomic insight into the basis of the insecticide tolerance/resistance of the pollen beetle Brassicogethes (Meligethes) aeneus

The pollen beetle is a major pest of oilseed rape. Although various resistance mechanisms have been identified, such as kdr (mutation in the sodium channel) and metabolic resistance (CYP overexpression), other "hidden" factors also exist. Some studies have stressed the importance of epistasis as a genetic background. The combination of kdr and metabolic resistance appears to be unfavorable under field conditions in the absence of pesticide selection. The regulation of detoxification enzymes can play an important role, but we highlight different detoxification markers compared to those emphasized in other studies. We also stress the importance of studying the role of markers identified as pathogenesis-related protein 5-like (PR5; upregulated by insecticides) and highlight the role of RNA (DEAD-box) helicases (downregulated by insecticides). Thus, we suggest the importance of epigenetic drivers of resistance/tolerance to pesticides. The key results are similar to those of our previous study, in which deltamethrin treatment of the pollen beetle was also investigated by a proteogenomic approach. Indeed, the mechanism leading to resistance of the pollen beetle may be an innate mechanism that the pollen beetle can also employ in natural habitats, but under field conditions (pesticide exposure), this mechanism is used to survive in response to insecticides. SIGNIFICANCE: Pesticide resistance is a serious problem that hampers the successful production of crops. Understanding the mechanisms of insecticide resistance is highly important for successful pest control, especially when considering integrated pest management. Here, using a proteogenomic approach, we identified novel markers for understanding pollen beetle resistance to pesticides. In addition, future studies will reveal the role of these markers in the multiresistance of pollen beetle populations. We highlight that the proteins identified as PR5, which are known to occur in beetles and are similar to those in plants, may be responsible for tolerance to multiple stresses. In addition, our results indicate that the RNA helicases that exhibited changes in expression may be the epigenetic drivers of multiresistance. The nature of these changes remains an open question, and their relevance in different situations (responses to different stresses) in natural habitats in the absence of pesticides can be proposed.

Advances in venomics: Modern separation techniques and mass spectrometry

Snake venoms are complex chemical mixtures of biologically active proteins and non-protein components. Toxins have a wide range of targets and effects to include ion channels and membrane receptors, and platelet aggregation and platelet plug formation. Toxins target these effectors and effects at high affinity and selectivity. From a pharmacological perspective, snake venom compounds are a valuable resource for drug discovery and development. However, a major challenge to drug discovery using snake venoms is isolating and analyzing the bioactive proteins and peptides in these complex mixtures. Getting molecular information from complex mixtures such as snake venoms requires proteomic analyses, generally combined with transcriptomic analyses of venom glands. The present review summarizes current knowledge and highlights important recent advances in venomics with special emphasis on contemporary separation techniques and bioinformatics that have begun to elaborate the complexity of snake venoms. Several analytical techniques such as two-dimensional gel electrophoresis, RP-HPLC, size exclusion chromatography, ion exchange chromatography, MALDI-TOF-MS, and LC-ESI-QTOF-MS have been employed in this regard. The improvement of separation approaches such as multidimensional-HPLC, 2D-electrophoresis coupled to soft-ionization (MALDI and ESI) mass spectrometry has been critical to obtain an accurate picture of the startling complexity of venoms. In the case of bioinformatics, a variety of software tools such as PEAKS also has been used successfully. Such information gleaned from venomics is important to both predicting and resolving the biological activity of the active components of venoms, which in turn is key for the development of new drugs based on these venom components.

Promising biomarkers of human aging: In search of a multi-omics panel to understand the aging process from a multidimensional perspective

The aging process has been linked to the occurrence of chronic diseases and functional impairments, including cancer, sarcopenia, frailty, metabolic, cardiovascular, and neurodegenerative diseases. Nonetheless, aging is highly variable and heterogeneous and represents a challenge for its characterization. In this sense, intrinsic capacity (IC) stands as a novel perspective by the World Health Organization, which integrates the individual wellbeing, environment, and risk factors to understand aging. However, there is a lack of quantitative and qualitative attributes to define it objectively. Therefore, in this review we attempt to summarize the most relevant and promising biomarkers described in clinical studies at date over different molecular levels, including epigenomics, transcriptomics, proteomics, metabolomics, and the microbiome. To aid gerontologists, geriatricians, and biomedical researchers to understand the aging process through the IC. Aging biomarkers reflect the physiological state of individuals and the underlying mechanisms related to homeostatic changes throughout an individual lifespan; they demonstrated that aging could be measured independently of time (that may explain its heterogeneity) and to be helpful to predict age-related syndromes and mortality. In summary, we highlight the areas of opportunity and gaps of knowledge that must be addressed to fully integrate biomedical findings into clinically useful tools and interventions.

Human and Nonhuman Primate Lineage-Specific Footprints in the Salivary Proteome

Proteins in saliva are needed for preprocessing food in the mouth, maintenance of tooth mineralization, and protection from microbial pathogens. Novel insights into human lineage-specific functions of salivary proteins and clues to their involvement in human disease can be gained through evolutionary studies, as recently shown for salivary amylase AMY1 and salivary agglutinin DMBT1/gp340. However, the entirety of proteins in saliva, the salivary proteome, has not yet been investigated from an evolutionary perspective. Here, we compared the proteomes of human saliva and the saliva of our closest extant evolutionary relatives, chimpanzees and gorillas, using macaques as an outgroup, with the aim to uncover features in saliva protein composition that are unique to each species. We found that humans produce a waterier saliva, containing less than half total protein than great apes and Old World monkeys. For all major salivary proteins in humans, we could identify counterparts in chimpanzee and gorilla saliva. However, we discovered unique protein profiles in saliva of humans that were distinct from those of nonhuman primates. These findings open up the possibility that dietary differences and pathogenic pressures may have shaped a distinct salivary proteome in the human lineage.

Proteomics of Osmoregulatory Responses in Threespine Stickleback Gills

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

Proteomic similarity of the Littorinid snails in the evolutionary context

BACKGROUND

The introduction of DNA-based molecular markers made a revolution in biological systematics. However, in cases of very recent divergence events, the neutral divergence may be too slow, and the analysis of adaptive part of the genome is more informative to reconstruct the recent evolutionary history of young species. The advantage of proteomics is its ability to reflect the biochemical machinery of life. It may help both to identify rapidly evolving genes and to interpret their functions.

METHODS

Here we applied a comparative gel-based proteomic analysis to several species from the gastropod family Littorinidae. Proteomes were clustered to assess differences related to species, geographic location, sex and body part, using data on presence/absence of proteins in samples and data on protein occurrence frequency in samples of different species. Cluster support was assessed using multiscale bootstrap resampling and the stability of clustering-using cluster-wise index of cluster stability. Taxon-specific protein markers were derived using IndVal method. Proteomic trees were compared to consensus phylogenetic tree (based on neutral genetic markers) using estimates of the Robinson-Foulds distance, the Fowlkes-Mallows index and cophenetic correlation.

RESULTS

Overall, the DNA-based phylogenetic tree and the proteomic similarity tree had consistent topologies. Further, we observed some interesting deviations of the proteomic littorinid tree from the neutral expectations. (1) There were signs of molecular parallelism in two Littoraria species that phylogenetically are quite distant, but live in similar habitats. (2) Proteome divergence was unexpectedly high between very closely related Littorina fabalis and L. obtusata, possibly reflecting their ecology-driven divergence. (3) Conservative house-keeping proteins were usually identified as markers for cryptic species groups ("saxatilis" and "obtusata" groups in the Littorina genus) and for genera (Littoraria and Echinolittorina species pairs), while metabolic enzymes and stress-related proteins (both potentially adaptively important) were often identified as markers supporting species branches. (4) In all five Littorina species British populations were separated from the European mainland populations, possibly reflecting their recent phylogeographic history. Altogether our study shows that proteomic data, when interpreted in the context of DNA-based phylogeny, can bring additional information on the evolutionary history of species.

Seasonal proteome variation in intertidal shrimps under a natural setting: Connecting molecular networks with environmental fluctuations

The ability of intertidal organisms to maintain their performance via molecular and physiological adjustments under low tide, seasonal fluctuations and extreme events ultimately determines population viability. Analyzing this capacity in the wild is extremely relevant since intertidal communities are under increased climate variability owing to global changes. We addressed the seasonal proteome signatures of a key intertidal species, the shrimp Palaemon elegans, in a natural setting. Shrimps were collected during spring and summer seasons at low tides and were euthanized in situ. Environmental variability was also assessed using hand-held devices and data loggers. Muscle samples were taken for 2D gel electrophoresis and protein identification through mass spectrometry. Proteome data revealed that 55 proteins (10.6% of the proteome) significantly changed between spring and summer collected shrimps, 24 of which were identified. These proteins were mostly involved in cytoskeleton remodelling, energy metabolism and transcription regulation. Overall, shrimps modulate gene expression leading to metabolic and structural adjustments related to seasonal differences in the wild (i.e. abiotic variation and possibly intrinsic cycles of reproduction and growth). This potentially promotes performance and fitness as suggested by the higher condition index in summer-collected shrimps. However, inter-individual variation (% coefficient of variation) in protein levels was quite low (min-max ranges were 0.6-8.3% in spring and 1.2-4.8% in summer), possibly suggesting reduced genetic diversity or physiological canalization. Protein plasticity is relevant to cope with present and upcoming environmental variation related to anthropogenic forcing (e.g. global change, pollution) but low inter-individual variation may limit evolutionary potential of shrimp populations.

Chronic kidney disease: a review of proteomic and metabolomic approaches to membranous glomerulonephritis, focal segmental glomerulosclerosis, and IgA nephropathy biomarkers

Chronic Kidney Disease (CKD) is a global health problem annually affecting millions of people around the world. It is a comprehensive syndrome, and various factors may contribute to its occurrence. In this study, it was attempted to provide an accurate definition of chronic kidney disease; followed by focusing and discussing on molecular pathogenesis, novel diagnosis approaches based on biomarkers, recent effective antigens and new therapeutic procedures related to high-risk chronic kidney disease such as membranous glomerulonephritis, focal segmental glomerulosclerosis, and IgA nephropathy, which may lead to end-stage renal diseases. Additionally, a considerable number of metabolites and proteins that have previously been discovered and recommended as potential biomarkers of various CKD_s using ‘-omics-’ technologies, proteomics, and metabolomics were reviewed.

Comparative proteomics of stenotopic caddisfly Crunoecia irrorata identifies acclimation strategies to warming

Species' ecological preferences are often deduced from habitat characteristics thought to represent more or less optimal conditions for physiological functioning. Evolution has led to stenotopic and eurytopic species, the former having decreased niche breadths and lower tolerances to environmental variability. Species inhabiting freshwater springs are often described as being stenotopic specialists, adapted to the stable thermal conditions found in these habitats. Whether due to past local adaptation these species have evolved or have lost intra-generational adaptive mechanisms to cope with increasing thermal variability has, to our knowledge, never been investigated. By studying how the proteome of a stenotopic species changes as a result of increasing temperatures, we investigate if the absence or attenuation of molecular mechanisms is indicative of local adaptation to freshwater springs. An understanding of compensatory mechanisms is especially relevant as spring specialists will experience thermal conditions beyond their physiological limits due to climate change. In this study, the stenotopic species Crunoecia irrorata (Trichoptera: Lepidostomatidae, Curtis 1834) was acclimated to 10, 15 and 20°C for 168 hr. We constructed a homology-based database and via liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based shotgun proteomics identified 1,358 proteins. Differentially abundant proteins and protein norms of reaction revealed candidate proteins and molecular mechanisms facilitating compensatory responses such as trehalose metabolism, tracheal system alteration and heat-shock protein regulation. A species-specific understanding of compensatory physiologies challenges the characterization of species as having narrow tolerances to environmental variability if that characterization is based on occurrences and habitat characteristics alone.

Chemical Responses to the Biotic and Abiotic Environment by Early Diverging Metazoans Revealed in the Post-Genomic Age

For many years methodological constraints limited insights on the molecular biology of non-model organisms. However, the development of various sequencing platforms has led to an explosion of transcriptomic and genomic data on non-model systems. As a consequence the molecular drivers of organismal phenotypes are becoming clearer and the chemicals that animals use to detect and respond to their environments are increasingly being revealed-this latter area inspired our symposium theme. The papers in this volume broadly address this theme by their more specific focus in one of the following general areas: 1) sensory biology and the molecular basis of perception, 2) chemicals deployed to deal with the biotic and abiotic environment, and 3) chemical interactions along the parasite-mutualist continuum. Here we outline and synthesize the content of these papers-an exercise which demonstrates that sophisticated gene repertoires enable early diverging metazoans to encode many of the signaling, sensory, defensive, and offensive capacities typically associated with animals that have complex nervous systems. We then consider opportunities and associated challenges that may delay progress in comparative functional biochemistry, a reinvigorated field that can be expected to rapidly expand with new 'omics data. Future knowledge of chemical adaptations should afford new perspectives on the comparative evolution of chemical mediators.

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.

RNA-seq coupled to proteomic analysis reveals high sperm proteome variation between two closely related marine mussel species

Speciation mechanisms in marine organisms have attracted great interest because of the apparent lack of substantial barriers to genetic exchange in marine ecosystems. Marine mussels of the Mytilus edulis species complex provide a good model to study mechanisms underlying species formation. They hybridise extensively at many localities and both pre- and postzygotic isolating mechanisms may be operating. Mussels have external fertilisation and sperm cells should show specific adaptations for survival and successful fertilisation. Sperm thus represent key targets in investigations of the molecular mechanisms underlying reproductive isolation. We undertook a deep transcriptome sequencing (RNA-seq) of mature male gonads and a 2DE/MS-based proteome analysis of sperm from Mytilus edulis and M. galloprovincialis raised in a common environment. We provide evidence of extensive expression differences between the two mussel species, and general agreement between the transcriptomic and proteomic results in the direction of expression differences between species. Differential expression is marked for mitochondrial genes and for those involved in spermatogenesis, sperm motility, sperm-egg interactions, the acrosome reaction, sperm capacitation, ATP reserves and ROS production. Proteins and their corresponding genes might thus be good targets in further genomic analysis of reproductive barriers between these closely related species. SIGNIFICANCE: Model systems for the study of fertilization include marine invertebrates with external fertilisation, such as abalones, sea urchins and mussels, because of the ease with which large quantities of gametes released into seawater can be collected after induced spawning. Unlike abalones and sea urchins, hybridisation has been reported between mussels of different Mytilus spp., which thus makes them very appealing for the study of reproductive isolation at both pre- and postzygotic levels. There is a lack of empirical proteomic studies on sperm samples comparing different Mytilus species, which could help to advance this study. A comparative analysis of sperm proteomes across different taxa may provide important insights into the fundamental molecular processes and mechanisms involved in reproductive isolation. It might also contribute to a better understanding of sperm function and of the adaptive evolution of sperm proteins in different taxa. There is now growing evidence from genomics studies that multiple protein complexes and many individual proteins might have important functions in sperm biology and the fertilisation process. From an applied perspective, the identification of sperm-specific proteins could also contribute to the improved understanding of fertility problems and as targets for fertility control.

Limited proteomic response in the marine snail Melarhaphe neritoides after long-term emersion

Rocky intertidal organisms are commonly exposed to environmental gradients, promoting adaptations to these conditions. Emersion time varies along the intertidal range and in the supralittoral zone is frequently larger than a single tidal cycle, even lasting for weeks. The planktonic-dispersing gastropod Melarhaphe neritoides is a common species of the high shore, adapted to reduce water loss in order to survive during long-term emersion. In this study, we investigated the molecular response, at the proteome level, of M. neritoides collected in high-shore tide pools to a series of emersion periods, from 8 to 24 days, in laboratory conditions. We compared this response to individuals maintained submerged during this period, because this was their original habitat. We also included a reversion treatment in the study, in which emersed individuals were returned to the submerged conditions. Although we detected an increase in overall protein concentration with longer emersion periods, contrary to general expectation, the two dimensional electrophoresis (2DE)-based proteomic analysis did not show significant differences between the treatments at the level of individual protein spots, even after an emersion period of 24 days. Our results suggest that the metabolism remains unaltered independent of the treatment carried out or the changes are very subtle and therefore difficult to detect with our experimental design. We conclude that M. neritoides could be equally adapted to emersion and submersion without drastic physiological changes.

Advancing behavioural genomics by considering timescale

Animal behavioural traits often covary with gene expression, pointing towards a genomic constraint on organismal responses to environmental cues. This pattern highlights a gap in our understanding of the time course of environmentally responsive gene expression, and moreover, how these dynamics are regulated. Advances in behavioural genomics explore how gene expression dynamics are correlated with behavioural traits that range from stable to highly labile. We consider the idea that certain genomic regulatory mechanisms may predict the timescale of an environmental effect on behaviour. This temporally minded approach could inform both organismal and evolutionary questions ranging from the remediation of early life social trauma to understanding the evolution of trait plasticity.

Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a QST - FST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Eco-Metabolomics and Metabolic Modeling: Making the Leap From Model Systems in the Lab to Native Populations in the Field

Experimental high-throughput analysis of molecular networks is a central approach to characterize the adaptation of plant metabolism to the environment. However, recent studies have demonstrated that it is hardly possible to predict in situ metabolic phenotypes from experiments under controlled conditions, such as growth chambers or greenhouses. This is particularly due to the high molecular variance of in situ samples induced by environmental fluctuations. An approach of functional metabolome interpretation of field samples would be desirable in order to be able to identify and trace back the impact of environmental changes on plant metabolism. To test the applicability of metabolomics studies for a characterization of plant populations in the field, we have identified and analyzed in situ samples of nearby grown natural populations of Arabidopsis thaliana in Austria. A. thaliana is the primary molecular biological model system in plant biology with one of the best functionally annotated genomes representing a reference system for all other plant genome projects. The genomes of these novel natural populations were sequenced and phylogenetically compared to a comprehensive genome database of A. thaliana ecotypes. Experimental results on primary and secondary metabolite profiling and genotypic variation were functionally integrated by a data mining strategy, which combines statistical output of metabolomics data with genome-derived biochemical pathway reconstruction and metabolic modeling. Correlations of biochemical model predictions and population-specific genetic variation indicated varying strategies of metabolic regulation on a population level which enabled the direct comparison, differentiation, and prediction of metabolic adaptation of the same species to different habitats. These differences were most pronounced at organic and amino acid metabolism as well as at the interface of primary and secondary metabolism and allowed for the direct classification of population-specific metabolic phenotypes within geographically contiguous sampling sites.

Biological information systems: Evolution as cognition-based information management

An alternative biological synthesis is presented that conceptualizes evolutionary biology as an epiphenomenon of integrated self-referential information management. Since all biological information has inherent ambiguity, the systematic assessment of information is required by living organisms to maintain self-identity and homeostatic equipoise in confrontation with environmental challenges. Through their self-referential attachment to information space, cells are the cornerstone of biological action. That individualized assessment of information space permits self-referential, self-organizing niche construction. That deployment of information and its subsequent selection enacted the dominant stable unicellular informational architectures whose biological expressions are the prokaryotic, archaeal, and eukaryotic unicellular forms. Multicellularity represents the collective appraisal of equivocal environmental information through a shared information space. This concerted action can be viewed as systematized information management to improve information quality for the maintenance of preferred homeostatic boundaries among the varied participants. When reiterated in successive scales, this same collaborative exchange of information yields macroscopic organisms as obligatory multicellular holobionts. Cognition-Based Evolution (CBE) upholds that assessment of information precedes biological action, and the deployment of information through integrative self-referential niche construction and natural cellular engineering antecedes selection. Therefore, evolutionary biology can be framed as a complex reciprocating interactome that consists of the assessment, communication, deployment and management of information by self-referential organisms at multiple scales in continuous confrontation with environmental stresses.

Differences in protein expression among five species of stream stonefly (Plecoptera) along a latitudinal gradient in Japan

Proteome variation among natural populations along an environmental gradient may provide insights into how the biological functions of species are related to their local adaptation. We investigated protein expression in five stream stonefly species from four geographic regions along a latitudinal gradient in Japan with varying climatic conditions. The extracted proteins were separated by two-dimensional gel electrophoresis and identified by matrix-assisted laser desorption/ionization of time-of-flight (MALDI TOF/TOF), yielding 446 proteins. Low interspecies variation in the proteome profiles was observed among five species within geographical regions, presumably due to the co-occurring species sharing the environments. However, large spatial variations in protein expression were found among four geographic regions, suggesting strong regulation of protein expression in heterogeneous environments, where the spatial variations were positively correlated with water temperature. We identified 21 unique proteins expressed specifically in a geographical region and six common proteins expressed throughout all regions. In warmer regions, metabolic proteins were upregulated, whereas proteins related to cold stress, the photoperiod, and mating were downregulated. Oxygen-related and energy-production proteins were upregulated in colder regions with higher altitudes. Thus, our proteomic approach is useful for identifying and understanding important biological functions related to local adaptations by populations of stoneflies.

Molecular Plasticity under Ocean Warming: Proteomics and Fitness Data Provides Clues for a Better Understanding of the Thermal Tolerance in Fish

Ocean warming is known to alter the performance and fitness of marine organisms albeit the proteome underpinnings of species thermal tolerance are still largely unknown. In this 1-month experiment we assessed the vulnerability of the gilt-head sea bream Sparus aurata, taken here as a biological model for some key fisheries species, to ocean warming (control 18°C, nursery ground temperature 24°C and heat wave 30°C). Survival was impaired after 28 days, mainly at 30°C although fishes' condition was unaltered. Muscle proteome modulation was assessed at 14 and 21 days, showing that protein expression profiles were similar between fish exposed to 18 and 24°C, differing from fish exposed to 30°C. Fish subjected to 24°C showed an enhanced glycolytic potential and decreased glycogenolysis mainly at 14 days of exposure. Fish subjected to 30°C also showed enhanced glycolytic potential and up-regulated proteins related to gene expression, cellular stress response (CSR), and homeostasis (mostly cytoskeletal dynamics, acid-base balance, chaperoning). However, inflammatory processes were elicited at 21 days along with a down-regulation of the tricarboxylic acid cycle. Thus, juvenile fish seem able to acclimate to 24°C but possibly not to 30°C, which is the predicted temperature for estuaries during heat waves by the year 2100. This may be related with increasing constraints on organism physiology associated with metabolic scope available for performance and fitness at higher temperatures. Consequently, recruitment of commercial sea breams may be in jeopardy, highlighting the need for improved management plans for fish stocks.

Absolute venomics: Absolute quantification of intact venom proteins through elemental mass spectrometry

We report the application of a hybrid element and molecular MS configuration for the parallel absolute quantification of μHPLC-separated intact sulfur-containing venom proteins, via ICP triple quadrupole MS and ³²S/³⁴S isotope dilution analysis, and identification by ESI-QToF-MS of the toxins of the medically important African black-necked spitting cobra, Naja nigricollis (Tanzania); New Guinea small-eyed snake, Micropechis ikaheka; and Papuan black snake, Pseudechis papuanus. The main advantage of this approach is that only one generic sulfur-containing standard is required to quantify each and all intact Cys- and/or Met-containing toxins of the venom proteome. The results of absolute quantification are in reasonably good agreement with previously reported relative quantification of the most abundant protein families. However, both datasets depart in the quantification of the minor ones, showing a tendency for this set of proteins to be underestimated in standard peptide-centric venomics approaches. The molecular identity, specific toxic activity, and concentration in the venom, are the pillars on which the toxicovenomics-aimed discovery of the most medically-relevant venom toxins, e.g. those that need to be neutralized by an effective therapeutic antivenom, should be based. The pioneering venom proteome-wide absolute quantification shown in this paper represents thus a significant advance towards this goal. The potential of ICP triple quadrupole MS in proteomics in general, and venomics in particular, is critically discussed.

BIOLOGICAL SIGNIFICANCE

Animal venoms provide excellent model systems for investigating interactions between predators and prey, and the molecular mechanisms that contribute to adaptive protein evolution. On the other hand, numerous cases of snake bites occur yearly by encounters of humans and snakes in their shared natural environment. Snakebite envenoming is a serious global public health issue that affects the most impoverished and geopolitically disadvantaged rural communities in many tropical and subtropical countries. Unveiling the temporal and spatial patterns of venom variability is of fundamental importance to understand the molecular basis of envenoming, a prerequisite for developing therapeutic strategies against snakebite envenoming. Research on venoms has been continuously enhanced by advances in technology. The combined application of next-generation transcriptomic and venomic workflows has demonstrated unparalleled capabilities for venom characterization in unprecedented detail. However, mass spectrometry is not inherently quantitative, and this analytical limitation has sparked the development of methods to determine absolute abundance of proteins in biological samples. Here we show the potential of a hybrid element and molecular MS configuration for the parallel ESI-QToF-MS and ICP-QQQ detection and absolute quantification of intact sulfur-containing venom proteins via ³²S/³⁴S isotope dilution analysis. This configuration has been applied to quantify the toxins of the medically important African snake Naja nigricollis (Tanzania), and the Papuan species Micropechis ikaheka and Pseudechis papuanus.

Evolutionary Transcriptomics and Proteomics: Insight into Plant Adaptation

Comparative transcriptomics and proteomics (T&P) have brought biological insight into development, gene function, and physiological stress responses. However, RNA-seq and high-throughput proteomics remain underutilised in studies of plant adaptation. These methodologies have created discovery tools with the potential to significantly advance our understanding of adaptive diversification. We outline experimental recommendations for their application. We discuss analysis models and approaches that accelerate the identification of adaptive gene sets and integrate transcriptome, proteome, phenotypic, and environmental data. Finally, we encourage widespread uptake and future developments in T&P that will advance our understanding of evolution and adaptation.

Protein-species quantitative venomics: looking through a crystal ball

In this paper we discuss recent significant developments in the field of venom research, specifically the emergence of top-down proteomic applications that allow achieving compositional resolution at the level of the protein species present in the venom, and the absolute quantification of the venom proteins (the term “protein species” is used here to refer to all the different molecular forms in which a protein can be found. Please consult the special issue of Jornal of Proteomics “Towards deciphering proteomes via the proteoform, protein speciation, moonlighting and protein code concepts” published in 2016, vol. 134, pages 1-202). Challenges remain to be solved in order to achieve a compact and automated platform with which to routinely carry out comprehensive quantitative analysis of all toxins present in a venom. This short essay reflects the authors’ view of the immediate future in this direction for the proteomic analysis of venoms, particularly of snakes.