Proteomic analysis in kidneys of Atlantic salmon infected with Aeromonas salmonicida by iTRAQ

Proteomics analysis of differentially abundant proteins in the rohu kidney infected with Edwardsiella tarda

Edwardsiella tarda (Et) is a zoonotic gram-negative pathogen with a diverse host range, including fish. However, the in-depth molecular mechanisms underlying the response of Labeo rohita (rohu) kidney to Et are poorly understood. A proteomic and histopathological analysis was performed for the rohu kidney after Et infection. The histopathology of the infected rohu kidney showed vacuolation and necrosis. After LC-MS/MS analysis, ~1240 proteins were identified with ≥2 unique peptides. A total of 96 differentially abundant proteins (DAPs) were observed between the control and Et infected group (ET). Metascape and STRING analysis were used for the gene ontology (GO), and protein-protein interaction network (PPI) for the significant pathways of DAPs. In PPI, low-abundant proteins were mapped to metabolic pathways and oxidative phosphorylation (cox5ab, uqcrfs1). High-abundance proteins were mapped to ribosomes (rplp2), protein process in the ER (hspa8), and immune system (ptgdsb.1, muc2). Our label-free proteomic approach in the rohu kidney revealed abundant enriched proteins involved in vesicle coat (ehd4), complement activation (c3a.1, c9, c7a), phagosome (thbs4, mapk1), metabolic reprogramming (hao1, glud1a), wound healing (vim, alox5), and the immune system (psap) after Et infection. A targeted proteomics approach of multiple reaction monitoring (MRM) validated the DAPs (nprl3, ambp, vmo1a, hspg2, muc2, hao1 and glud1a) between control and ET. Overall, the current analysis of histology and proteome in the rohu kidney provides comprehensive data on pathogenicity and the potential immune proteins against Et.

TMT-based proteomics analysis of sea urchin (Strongylocentrotus intermedius) under high temperature stress

In the context of global warming and continuous high temperatures in the northern part of China in summer, the mortality rate of Strongylocentrotus intermedius through the summer reaches 70-80 %. The protein regulatory mechanism of S. intermedius in response to high temperature stress is still unclear. In order to investigate the protein expression of S. intermedius under high temperature stress, the study was conducted with the high-temperature resistant strain of S. intermedius and the control group of S. intermedius. Tandem Mass Tag (TMT) tagging technique was applied to resolve the protein expression profile of S. intermedius in response to high temperature stress. The results showed that, compared to 15 °C，136 DEPs were screened in high-temperature resistant strain groups of S. intermedius under high temperature stress and 87 DEPs were screened in the control group of S. intermedius. There were 33 common differential proteins in the two groups, such as APOLP, HSP 70, CDC37 and CALM. Further Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses revealed that the up-regulated proteins CALM and HSP70 are significantly enriched in the "Phosphatidylinositol signaling system" and "Protein processing in endoplasmic reticulum" in heat-tolerant S. intermedius strains under high temperature stress. The control group of S. intermedius DEPs were significantly enriched in protein processing in the endoplasmic reticulum. These results provide a theoretical basis for the molecular mechanism of sea urchin heat tolerance and fundamental data for sea urchin selection and breeding for high temperature tolerance.

Plasma proteomic profiling of bacterial cold water disease-resistant and -susceptible rainbow trout lines and biomarker discovery

Genetic variation for disease resistance is present in salmonid fish; however, the molecular basis is poorly understood, and biomarkers of disease susceptibility/resistance are unavailable. Previously, we selected a line of rainbow trout for high survival following standardized challenge with Flavobacterium psychrophilum (Fp), the causative agent of bacterial cold water disease. The resistant line (ARS-Fp-R) exhibits over 60 percentage points higher survival compared to a reference susceptible line (ARS-Fp-S). To gain insight into the differential host response between genetic lines, we compared the plasma proteomes from day 6 after intramuscular challenge. Pooled plasma from unhandled, PBS-injected, and Fp-injected groups were simultaneously analyzed using a TMT 6-plex label, and the relative abundance of 513 proteins was determined. Data are available via ProteomeXchange, with identifier PXD041308, and the relative protein abundance values were compared to mRNA measured from a prior, whole-body RNA-seq dataset. Our results identified a subset of differentially abundant intracellular proteins was identified, including troponin and myosin, which were not transcriptionally regulated, suggesting that these proteins were released into plasma following pathogen-induced tissue damage. A separate subset of high-abundance, secreted proteins were transcriptionally regulated in infected fish. The highest differentially expressed protein was a C1q family member (designated complement C1q-like protein 3; C1q-LP3) that was upregulated over 20-fold in the infected susceptible line while only modestly upregulated, 1.8-fold, in the infected resistant line. Validation of biomarkers was performed using immunoassays and C1q-LP3, skeletal muscle troponin C, cathelcidin 2, haptoglobin, leptin, and growth and differentiation factor 15 exhibited elevated concentration in susceptible line plasma. Complement factor H-like 1 exhibited higher abundance in the resistant line compared to the susceptible line in both control and challenged fish and thus was a baseline differentiator between lines. C1q-LP3 and STNC were elevated in Atlantic salmon plasma following experimental challenge with Fp. In summary, these findings further the understanding of the differential host response to Fp and identifies salmonid biomarkers that may have use for genetic line evaluation and on-farm health monitoring.

Proteomic analysis of carp seminal plasma provides insights into the immune response to bacterial infection of the male reproductive system

Aeromonas salmonicida is recognized as a significant bacterial pathogen in ulcerative disease of cyprinid fish. However, the mechanism of immunity to these bacteria in common carp is still not well understood, especially the immune regulation in the gonad to bacterial infection. The aims of our study were to analyze changes in the seminal plasma proteome following A. salmonicida infection in carp males. The observed pathological changes in the tissue (liver, spleen, kidney and testis) morphology and upregulation of immune-related genes (tnfa2, il6a) confirmed the successful infection challenge. Using mass spectrometry-based label-free quantitative proteomics, we identified 1402 seminal plasma proteins, and 44 proteins (20 up- and 24 downregulated) were found to be differentially abundant between infected and control males. Most differentially abundant proteins were involved in the immune response mechanisms, such as acute phase response, complement activation and coagulation, inflammation, lipid metabolism, cell-cell and cell-matrix adhesion, creatine-phosphate biosynthesis and germ cell-Sertoli cell junction signaling. Bacterial infection also caused profound changes in expression of selected genes in the testis and hematopoietic organs, which contributed to changes in seminal proteins. The altered seminal proteins and bacterial proteins in seminal plasma may serve as valuable markers of infection in the testis.

Identification of novel biomarkers of inflammation in Atlantic salmon (Salmo salar L.) by a plasma proteomic approach

Monitoring fish welfare has become a central issue for the fast-growing aquaculture industry, and finding proper biomarkers of stress, inflammation and infection is necessary for surveillance and documentation of fish health. In this study, a proteomic approach using mass spectrometry was applied to identify indicators of the acute response in Atlantic salmon blood plasma by comparing Aeromonas salmonicida subsp. salmonicida infected fish and non-infected controls. The antimicrobial proteins cathelicidin (CATH), L-plastin (Plastin-2, LCP1) and soluble toll-like receptor 5 (sTLR5) were uniquely or mainly identified in the plasma of infected fish. In addition, five immune-related proteins showed significantly increased expression in plasma of infected fish: haptoglobin, high affinity immunoglobulin Fc gamma receptor I (FcγR1, CD64), leucine-rich alpha 2 glycoprotein (LRG1), complement C4 (C4) and phospholipase A2 inhibitor 31 kDa subunit-like protein. However, various fibrinogen components, CD209 and CD44 antigen-like molecules decreased in infected fish. Selected biomarkers were further verified by Western blot analysis of plasma and real time PCR of spleen and liver, including CATH1, CATH2 and L-plastin. A significant increase of L-plastin occurred as early as 24 h after infection, and a CATH2 increase was observed from 72 h in plasma of infected fish. Real time PCR of selected genes confirmed increased transcription of CATH1 and CATH2. In addition, serum amyloid A mRNA significantly increased in liver and spleen after bacterial infection. However, transcription of L-plastin was not consistently induced in liver and spleen. The results of the present study reveal novel and promising biomarkers of the acute phase response and inflammation in Atlantic salmon.

Transcriptional Signatures of Immune, Neural, and Endocrine Functions in the Brain and Kidney of Rainbow Trout (Oncorhynchus mykiss) in Response to Aeromonas salmonicida Infection

Rainbow trout (Oncorhynchus mykiss) serves as one of the most important commercial fish with an annual production of around 800,000 tonnes. However, infectious diseases, such as furunculosis caused by Aeromonas salmonicida infection, results in great economic loss in trout culture. The brain and kidney are two important organs associated with “sickness behaviors” and immunomodulation in response to disease. Therefore, we worked with 60 trout and investigated transcriptional responses and enrichment pathways between healthy and infected trout. We observed that furunculosis resulted in the activation of toll-like receptors with neuroinflammation and neural dysfunction in the brain, which might cause the “sickness behaviors” of infected trout including anorexia and lethargy. We also showed the salmonid-specific whole genome duplication contributed to duplicated colony stimulating factor 1 (csf-1) paralogs, which play an important role in modulating brain immunomodulation. Enrichment analyses of kidneys showed up-regulated immunomodulation and down-regulated neural functions, suggesting an immune-neural interaction between the brain and kidney. Moreover, the kidney endocrine network was activated in response to A. salmonicida infection, further convincing the communications between endocrine and immune systems in regulating internal homeostasis. Our study provided a foundation for pathophysiological responses of the brain and kidney in response to furunculosis and potentially offered a reference for generating disease-resistant trout strains.

Genome-wide integrated analysis reveals functions of lncRNA-miRNA-mRNA interactions in Atlantic salmon challenged by Aeromonas salmonicida

Aeromonas salmonicida (A. salmonicida) is a pathogenic bacterium that causes serious problems in the global Atlantic salmon aquaculture industry. In this study, we comprehensively analyzed the profiles of lncRNAs, miRNAs and mRNAs in gills of Atlantic salmon at high-dose A. salmonicida infection (3.06 × 10⁸ CFU/mL), low-dose A. salmonicida infection (3.06 × 10⁵ CFU/mL), and a PBS (100 μL) control. We identified 65 differentially expressed lncRNAs, 41 miRNAs, and 512 mRNAs between the control group and infection groups. Functional analysis showed that these genes were significantly enriched in the p53 signaling pathway, Wnt signaling pathway, mTOR signaling pathway, JAK-STAT signaling pathway, and Toll-like receptor signaling pathway. In addition, we predicted key genes in immune-related pathways and constructed a lncRNA-miRNA-mRNA network based on whole transcriptomic analysis. We further predicted three lncRNA-miRNA-mRNA axes as potential novel biomarkers in regulating the immune response of Atlantic salmon against A. salmonicida infection.

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

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

A transcriptome analysis focusing on splenic immune-related mciroRNAs of rainbow trout upon Aeromonas salmonicida subsp. salmonicida infection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can regulate the immune responses during pathogen infection. Aeromonas salmonicida (A. salmonicida) subsp. salmonicida is the causative agent of furunculosis in salmon and trout. To identify the miRNAs and investigate the specific miRNAs in rainbow trout upon A. salmonicida subsp. salmonicida infection, we performed high throughput sequencing using the spleens of rainbow trout infected with and without an A. salmonicida subsp. salmonicida clinical isolate. A total of 381 known miRNAs and 926 novel miRNAs were identified. Eleven known and 16 novel miRNAs were found to be differentially expressed upon infection. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that the target genes of the differentially expressed miRNAs were closely associated with immune responses and biological regulations. Additionally, over- and suppressed expression of miR-155-5p significantly enhanced and reduced the IL-2 and IL-1β expressions in RTG-2 cells induced by A. salmonicida, respectively. To our knowledge, this is the first experimental study on the miRNAs of rainbow trout upon A. salmonicida infection. The results here might lay a foundation for the further understanding of the roles of miRNAs in the immune responses during A. salmonicida infection in rainbow trout.

The key differentially expressed genes and proteins related to immune response in the spleen of pufferfish (Takifugu obscurus) infected by Aeromonas hydrophila

The immune mechanism elicited in pufferfish (Takifugu obscurus) against the invasion of Aeromonas hydrophila is still poorly understood. We examined the spleen of pufferfish at the transcriptome and proteome levels by using Illumina-seq and TMT coupled mass spectrometry after 12 h infection by A. hydrophila, respectively. A total of 2,339 genes (1,512 up-regulated and 827 down-regulated) and 537 (237 up-regulated and 300 down-regulated) proteins were identified. GO and KEGG analyses revealed that the responses to stimulus were the main biological processes, intestinal immune network for IgT production and calcium signaling pathway. Fourteen genes (8 up-regulated and 6 down-regulated) and proteins (5 up-regulated and 9 down-regulated) involved immune responses or signal transduction were validated by qRT-PCR and parallel reaction monitoring to confirm the reliability of the transcriptomic and proteomic analyses, respectively. Moreover, qRT-PCR and flow cytometry were used to detect dynamics of the genes in calcium signaling pathway and changes of concentration of cytoplasm Ca²⁺ in spleen cells within a 72 h challenge. This study provides the findings regarding immune response, especially intestinal immune network for IgT production pathway and calcium signaling pathway at the molecular, protein and cellular in pufferfish after infection by A. hydrophila. These results would provide a new insight and molecular targets into the response to pathogenic infection in pufferfish.

Phosphoproteomic analyses of kidneys of Atlantic salmon infected with Aeromonas salmonicida

Aeromonas salmonicida (A. salmonicida) is a pathogenic bacterium that causes furunculosis and poses a significant global risk, particularly in economic activities such as Atlantic salmon (Salmo salar) farming. In a previous study, we identified proteins that are significantly upregulated in kidneys of Atlantic salmon challenged with A. salmonicida. Phosphoproteomic analyses were conducted to further clarify the dynamic changes in protein phosphorylation patterns triggered by bacterial infection. To our knowledge, this is the first study to characterize phosphorylation events in proteins from A. salmonicida-infected Atlantic salmon. Overall, we identified over 5635 phosphorylation sites in 3112 proteins, and 1502 up-regulated and 77 down-regulated proteins quantified as a 1.5-fold or greater change relative to control levels. Based on the combined data from proteomic and motif analyses, we hypothesize that five prospective novel kinases (VRK3, GAK, HCK, PKCδ and RSK6) with common functions in inflammatory processes and cellular pathways to regulate apoptosis and the cytoskeleton could serve as potential biomarkers against bacterial propagation in fish. Data from STRING-based functional network analyses indicate that fga is the most central protein. Our collective findings provide new insights into protein phosphorylation patterns, which may serve as effective indicators of A. salmonicida infection in Atlantic salmon.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

Anammox bacteria exhibit capacity to withstand long-term starvation stress: A proteomic-based investigation of survival mechanisms

Although anammox bacteria are commonly exposed to long-term starvation during transportation and preservation process, physiological changes in these organisms during long-term starvation are not well understood, nor are the molecular bases of their starvation survival strategies. To reveal survival mechanisms during long-term anaerobic and anoxic starvation (60 days at 20 ± 1 °C), metaproteomic technology was utilized to identify differentially expressed proteins in Candidatus Kuenenia stuttgartiensis. Our results showed that Candidatus Kuenenia stuttgartiensis exhibits a capacity to withstand long-term starvation stress. Although activity decay rates of 0.0129 d^-1 and 0.0049 d^-1 were observed for anammox sludge in anoxic and anaerobic starvation, the relative abundance of Candidatus Kuenenia stuttgartiensis, the shape of anammox granules, and the fraction of viable cells remained constant under both anaerobic and anoxic starvation conditions. Metaproteomics results illustrated that Candidatus Kuenenia stuttgartiensis maintained stable levels of most intracellular proteins, especially enzymes involved in principal metabolic pathways after 60-d of anaerobic or anoxic starvation, thereby allowing cells to regain metabolic activities once substrates became available. Induction of starvation proteins could be a survival strategy employed by Candidatus Kuenenia stuttgartiensis to resist long-term starvation stresses. During anaerobic starvation, 34 proteins were upregulated, five of which were associated with carbohydrate catabolism and oxidation of organic compounds, thereby increasing potential for utilization of endogenous carbon sources to produce energy. During anoxic starvation, only two proteins were upregulated, which may be attributed to insufficient energy for the synthesis of starvation-induced proteins.

Comparative proteomic analysis of hepatic mechanisms of Megalobrama amblycephala infected by Aeromonas hydrophila

Hemorrhage syndrome is one of the most prevalent and epidemic diseases that is mainly caused by Aeromonas hydrophila invasion in Megalobrama amblycephala. Recent studies have uncovered a number of immune enzymes and transcripts that are differently expressed in this disease, but the molecular mechanism elicited still remain largely unknown. Here, we constructed an in vivo A. hydrophila infection to investigate the immune mechanism in M. amblycephala using comparative proteomic approach at the one day after infection. 30 altered protein spots were found to undergo differential expression against A. hydrophila infection in the hepatopancreas of M. amblycephala based on 2-DE and were all successfully identified using MALDI-TOF/TOF, representing 18 unique proteins. These proteins were functionally classified into metabolism, antioxidant, cofactors and vitamins, chaperone and signal transduction. Network interaction and Gene Ontology annotation indicated 13 unique proteins were closely related to immune response and directly regulated by each other. Compared with the control group, A. hydrophila infection significantly decreased the metabolism-related mRNA expressions of ENO3, APOA1, CAT and FASN, but increased the mRNA expressions of MDH, ALDOB and RSP12, which was consistent with the protein expression. Nevertheless, FAH was down-regulated at both levels but had no significant difference in mRNA level, ALDH8a1 was down-regulated at protein level but non-significantly up-regulated at the mRNA level. GSTm was up-regulated at protein level but down-regulated at the mRNA level. Consequently, these results revealed that A. hydrophila infection altered the related antioxidative proteins via complex regulatory mechanisms and reduced the immune ability of M. amblycephala at the one day after infection.

High-throughput proteomic profiling of the fish liver following bacterial infection

Background

High-throughput proteomics was used to determine the role of the fish liver in defense responses to bacterial infection. This was done using a rainbow trout (Oncorhynchus mykiss) model following infection with Aeromonas salmonicida, the causative agent of furunculosis. The vertebrate liver has multifaceted functions in innate immunity, metabolism, and growth; we hypothesize this tissue serves a dual role in supporting host defense in parallel to metabolic adjustments that promote effective immune function. While past studies have reported mRNA responses to A. salmonicida in salmonids, the impact of bacterial infection on the liver proteome remains uncharacterized in fish.

Results

Rainbow trout were injected with A. salmonicida or PBS (control) and liver extracted 48 h later for analysis on a hybrid quadrupole-Orbitrap mass spectrometer. A label-free method was used for protein abundance profiling, which revealed a strong innate immune response along with evidence to support parallel rewiring of metabolic and growth systems. 3076 proteins were initially identified against all proteins (n = 71,293 RefSeq proteins) annotated in a single high-quality rainbow trout reference genome, of which 2433 were maintained for analysis post-quality filtering. Among the 2433 proteins, 109 showed significant differential abundance following A. salmonicida challenge, including many upregulated complement system and acute phase response proteins, in addition to molecules with putative functions that may support metabolic re-adjustments. We also identified novel expansions in the complement system due to gene and whole genome duplication events in salmonid evolutionary history, including eight C3 proteins showing differential changes in abundance.

Conclusions

This study provides the first high-throughput proteomic examination of the fish liver in response to bacterial challenge, revealing novel markers for the host defense response, and evidence of metabolic remodeling in conjunction with activation of innate immunity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5092-0) contains supplementary material, which is available to authorized users.

iTRAQ-based quantitative proteomic analysis of the earthworm Eisenia fetida response to Escherichia coli O157:H7

Soil environment contaminated by Escherichia coli O157:H7 which come from the waste of infected animals. Earthworms can live in the pathogens-polluted soil by their innate immunity. How the proteins of earthworms E. fetida will response to E. coli O157:H7-contaminated-soil still unclear? To identify the defense proteins under E. coli O157:H7 stress, we performed a proteomic analysis of earthworm under E. coli O157:H7 exposure through an iTRAQ technology. In total, we found 283 non-redundant proteins, including fibrinolytic protease 1, lombricine kinase, lysozyme, gelsolin, coelomic cytolytic factor-1, antimicrobial peptide lumbricin-l, lysenin, and et al. The proteins participate in metabolic processes, transcription, defense response to bacterium, translation, response to stress, and transport. The study will contribute to understand why earthworm can live in the pathogens-polluted environment.

Proteome Profiles of Head Kidney and Spleen of Rainbow Trout (Oncorhynchus Mykiss)

The head kidney and spleen are major lymphoid organs of the teleost fish. The authors identify proteome profiles of head kidney and spleen of rainbow trout (Oncorhynchus mykiss) using a shotgun proteomic approach. Gene ontology annotation of proteins is predicted using bioinformatic tools. This study represents detailed proteome profiles of head kidney and spleen of rainbow trout, with a total of 3241 and 2542 proteins identified, respectively. It is found that lymphoid organs are equipped with a variety of functional proteins related to defense, receptor, signal transduction, antioxidant, cytoskeleton, transport, binding, and metabolic processes. The identified proteome profiles will serve as a template for understanding lymphoid organ functions in salmonids and will increase the amount of spectra information of rainbow trout proteins in the public data repository PRIDE. This data can be accessed via ProteomeXchange with identifiers PXD008473 and PXD008478.

Comparative proteomic analysis of eggplant (Solanum melongena L.) heterostylous pistil development

Heterostyly is a common floral polymorphism, but the proteomic basis of this trait is still largely unexplored. In this study, self- and cross-pollination of L-morph and S-morph flowers and comparison of embryo sac development in eggplant (Solanum melongena L.) suggested that lower fruit set from S-morph flowers results from stigma-pollen incompatibility. To explore the molecular mechanism underlying heterostyly development, we conducted isobaric tags for relative and absolute quantification (iTRAQ) proteomic analysis of eggplant pistils for L- and S-morph flowers. A total of 5,259 distinct proteins were identified during heterostyly development. Compared S-morph flowers with L-morph, we discovered 57 and 184 differentially expressed proteins (DEPs) during flower development and maturity, respectively. Quantitative real time polymerase chain reactions were used for nine genes to verify DEPs from the iTRAQ approach. During flower development, DEPs were mainly involved in morphogenesis, biosynthetic processes, and metabolic pathways. At flower maturity, DEPs primarily participated in biosynthetic processes, metabolic pathways, and the formation of ribosomes and proteasomes. Additionally, some proteins associated with senescence and programmed cell death were found to be upregulated in S-morph pistils, which may lead to the lower fruit set in S-morph flowers. Although the exact roles of these related proteins are not yet known, this was the first attempt to use an iTRAQ approach to analyze proteomes of heterostylous eggplant flowers, and these results will provide insights into biochemical events taking place during the development of heterostyly.