key: cord-0966088-19zh9evl
authors: Hu, Hao‐Liang; Kang, Yu; Zeng, Yong; Zhang, Ming; Liao, Qiong; Rong, Ming‐Qiang; Zhang, Qin; Lai, Ren
title: Region‐resolved proteomics profiling of monkey heart
date: 2019-01-15
journal: J Cell Physiol
DOI: 10.1002/jcp.28052
sha: a8a2882316256ca578ec31863059897cb66a4039
doc_id: 966088
cord_uid: 19zh9evl

Nonhuman primates (NHPs) play an indispensable role in biomedical research because of their similarities in genetics, physiological, and neurological function to humans. Proteomics profiling of monkey heart could reveal significant cardiac biomarkers and help us to gain a better understanding of the pathogenesis of heart disease. However, the proteomic study of monkey heart is relatively lacking. Here, we performed the proteomics profiling of the normal monkey heart by measuring three major anatomical regions (vessels, valves, and chambers) based on iTRAQ‐coupled LC‐MS/MS analysis. Over 3,200 proteins were identified and quantified from three heart tissue samples. Furthermore, multiple bioinformatics analyses such as gene ontology analysis, protein–protein interaction analysis, and gene‐diseases association were used to investigate biological network of those proteins from each area. More than 60 genes in three heart regions are implicated with heart diseases such as hypertrophic cardiomyopathy, heart failure, and myocardial infarction. These genes associated with heart disease are mainly enriched in citrate cycle, amino acid degradation, and glycolysis pathway. At the anatomical level, the revelation of molecular characteristics of the healthy monkey heart would be an important starting point to investigate heart disease. As a unique resource, this study can serve as a reference map for future in‐depth research on cardiac disease‐related NHP model and novel biomarkers of cardiac injury.

especially mitral valve disease, could contribute to cardiac dysfunction, and ultimately lead to heart failure (Levine et al., 2015) . Atrial and ventricular hypertrophy can also cause pathological conditions, including hypertension and valvular insufficiency (Angeli, Reboldi, & Verdecchia, 2014; Bekeredjian & Grayburn, 2005; Cuspidi, Negri, Tadic, Sala, & Parati, 2014) .

So far, molecular mechanisms in several cardiac diseases are relatively unknown. Proteomic technology, which aimed to reveal intrinsic relationships between genotype and phenotype, provides useful tools for the diagnosis and treatment of heart-related diseases.

A draft of the human proteome assembled from human tissues, cell lines, and body fluids helps us to enhance the understanding of the human body at the molecular level Wilhelm et al., 2014) . Recently, the protein composition of the dissimilar regions and cell types of the healthy human heart has been revealed, and offer a valuable reference for future research aimed to uncover biomarkers indicative of cardiac compartments deterioration (Doll et al., 2017) .

Transcriptomic (Ounzain et al., 2015) , proteomic (Lau et al., 2016; Lu, Sinha, Sharma, Kislinger, & Gramolini, 2014; Xiang et al., 2013) , and phosphoproteomic (Lundby & Olsen, 2013; Peng et al., 2014) studies have also been applied to the analysis of mammalian hearts, and provide attractive candidate biomarkers for the clinical study.

Nonhuman primates (NHPs) play an important role in biomedical research because of their similarities in genetics, physiological, and neurological function to humans. The differential proteins of serum and heart tissue in NHP models associated with cardiac disorder provide significant insight to continue the cardiac biomarker discovery and understand the cellular mechanism of drug-induced heart injury (Liu et al., 2013; Song et al., 2014) . At the anatomical level, the revelation of molecular characteristics of the monkey heart, particularly the healthy state monkey heart, would be an important starting point to investigate heart disease. But, the study about protein composition in dissimilar regions of the monkey heart is relatively less.

Here, we conducted a systematic quantitative proteomics analysis on monkey heart proteins from three compartments, including vessels, valves, and chambers. The high-sensitivity sample preparation, peptide fractionation, and an advanced iTRAQ-coupled LC-MS workflow were used to quantify a total of more than 3,200 proteins. Our results establish proteomic and functional differences between normal monkey heart compartments and pinpoint potential cardiac diseases-related proteins in different regions.

Three non-naive rhesus macaques (male) with ages ranging from 4 to 6 years were assigned to the study. The selected animal was killed under deep anesthesia with pentobarbital. Then the heart with proximal aortic artery and pulmonary artery connected was excised. The study was approved by the Ethics Committee on Animal Experimentation of the Hunan Normal University (Protocol 259/2012) and the experiments were carried out according to the ethical guidelines for the Care and Use of Laboratory Animals (Ministry of Science and Technology of China, 2006) .

The heart anatomy of rhesus monkey was previously reported to be similar to that of human (Sleeper, Gaughan, Gleason, & Burkett, 2008) .

The harvest of vessels, valves, and myocardial tissues was done by an experienced cardiologist under the assistance of a veterinary. At first, discreetly isolated coronary vessels ran at the surface of heart from myocardium and then harvested them. After that, opened the left atrium (LA) and the left ventricular (LV) wall was dissected starting at the anterolateral commissure to harvest LA and LV myocardial tissues, and anterior, and posterior mitral leaflets. The similar process was performed for the right heart to collect the right atrium and ventricular myocardial tissues. The chordae tendineae consisted in atrioventricular valves (mitral and tricuspid) was not to be taken as the valve samples, as its distinctive organizational structure and cell types differed from semilunar aortic and pulmonary valves. Finally, the rest of valves were obtained after opening the aorta and pulmonary artery. Besides, the heart was identified with normal size and structure and without heart disease, that is, no signs of hypertrophy, atrophy or scar of myocardium, thickening, shortening, or stiffening of valves and atherosclerosis of great vessels, by using visual inspection of the gross appearance of explanted tissue by the cardiologist during the tissue harvest process. All tissues were snap frozen in liquid nitrogen and stored at −80°C until analysis.

Tissue samples were washed in phosphate-buffered saline (PBS) for three times and ground in liquid nitrogen. The powdered tissues were resuspended dissolved in lysis buffer (7 M urea, 2 M thiourea, 65 mM dithiothreitol, and 0.1 mM phenylmethylsulfonyl fluoride) and ultrasonicated for 3 min followed by dissolving for 1 hr at 4°C. After centrifugation at 14,000 × g for 1 hr at 4°C, the supernatant was collected and the protein concentration was determined by the 2D Quantification Kit (Amersham Biosciences, Freiburg, Baden-Württemberg, Germany).

Trypsin digestion and iTRAQ labeling were performed according to the manufacturer's protocol (Sicex, Framingham, MA) . In brief, 100 μg proteins of each sample were reduced, alkylated and precipitated using acetone and resuspended using 50 mM triethylammonium bicarbonate (TEAB), pH 8.5. Then the protein samples were digested by trypsin (Promega, Madison, WI) overnight at 37°C. The samples were subsequently labeled with iTRAQ reagents as following; Vessel (MO1) with 115 tags; Vessel (MO2) with 116 tags; Vessel (MO3) with 116 tags; Mix (MO1-3) with 114 tags. And, the same tag method was used in the other two tissues (valve and chamber). A total of six sets of iTRAQ analysis (per regional sample was repeated two times) were constructed in each monkey. Each set of iTRAQ labeled peptides was mixed, dried, and desalted for further analysis.

The labeled peptides were diluted with the loading buffer (5 mM ammonium formate containing 2% acetonitrile, pH 10) and fractionated HU ET AL. | 13721 by using high pH reversed phase liquid chromatography on a UPLC system (Waters, Milford, MA). A C-18 column (2.1 × 250 mm X Bridge BEH300) was used for the separation. The gradient elution was performed by 0-25% B (5 mM ammonium formate containing 98% acetonitrile, pH 10, 5-35 min) and 25-45% B (35-48 min) on high pH RPLC column (Waters, Xbridge C18 3.5 µm, 150 × 2.1 mm, Milford, MA) with a flow rate of 300 µl/min. A total of 20 fractions were collected which were then combined to 10 fractions.

Each fraction was dried, dissolved in buffer A (0.1% FA, H 2 O) and analyzed on a LC-MS/MS platform which was equipped with an Orbitrap Fusion™ mass spectrometry (Thermo Fisher Scientific, San Jose, CA) and a nano Easy-Spray (Thermo Fisher Scientific, San Jose, CA). Peptides were separated with a linear gradient from 5% B (98% ACN with 0.1% formic acid) to 40% B in 90 min within a total analysis of 120 min. Orbitrap Fusion™ mass spectrometer was operated in the data-dependent mode to switch automatically between MS and MS/MS acquisition. Survey fullscan MS spectra (m/z 300-1200) were acquired with a mass resolution of 70 K, followed by 10 sequential high energy collisional dissociation (HCD) MS/MS scans (100-2000) with a resolution of 17.5 K. The dynamic exclusion time was set as 30 s. For MS/MS, precursor ions were activated using 27% normalized collision energy.

All the raw data were submitted to Proteome Discoverer software (version 1.3, Thermo scientific, Sunnyvale, CA) using the MASCOT search engine for protein identification and quantitation.

The data was search against rhesus Uniprot database with 161,948 entries. Parameters searching were set as following: Quantitation: iTRAQ-8plex; enzyme: Trypsin; variable modification: Oxidation (Met); fixed modification: Methyl methanethiosulfonate (Cys). MS1 tolerance: 20 ppm; MS/MS tolerance: 20 mmu; with one missed cleavage on trypsin was allowed. The false discovery rate (FDR) of 1% was used to control protein level identification based on the target-decoy strategy. The peptides which were unique to the given proteins and with a confidence higher than 95% were used for the calculation of the protein ratios.

Based on cytoscape platform, the plugin ClueGO in conjunction with CluePedia created a dynamic and visualized network to analyze biological functions and signaling pathways of monkey cardiac proteins. All the significant pathways or terms with a statistical threshold (p < 0.05) were functionally grouped and highlighted with different colors (Bindea et al., 2009; Bindea, Galon, & Mlecnik, 2013 

The Reactome database can be obtained from http://www.reactome. org/, which comprises a set of core modules called the reaction (Fabregat et al., 2018) . Entities (nucleic acids, proteins, complexes, etc.) participating in reactions form a network of biological interactions and are grouped into several pathways like classical intermediary metabolism, signaling, transcriptional regulation, apoptosis, and disease (Matthews et al., 2009) . To illuminate the correlation between monkey heart genes and human pathways, we imported the healthy monkey cardiac gene list into Reactome online tools for pathway analysis.

The STRING (Search Tool for the Retrieval of Interacting Genes/ Proteins) online engine, which aimed to assemble, evaluate, and disseminate protein-protein association information, was used to estimate the interactions between proteins in each area (Matthews et al., 2009) . To narrow down the lists and enhance the readability of the result, we merely selected the core proteins which were contained in principal biological processes. Be worth what carry is, the STRING data lists were imported into Cytoscape for network reconstruction and visualization, and were emphasized by the level of the combined score. Furthermore, the interaction information between proteins in dissimilar heart regions was presented in Figure 5 and Supporting Information S3-4.

To better understand the association between cardiac genes and diseases, the core gene lists from protein-protein association analysis in each heart region above were put through the DisGeNET database (http://www.disgenet.org/web/DisGeNET /menu/home) for gene-disease association analysis. The associations are all supported by experimental evidence or previous publications. The associations from the DisGeNET are all supported by expert curated repositories, GWAS catalogs, animal models, and the scientific literature (Pinero et al., 2017) . Then, we mapped our lists to the knowledge base (more than 500 thousand gene-disease associations [GDAs] ), so that all the GDAs could be extracted; one gene may be involved in several diseases, and one disease may include numerous genes. Finally, the cardiovascular disorder-related genes in each heart region were further enriched by signaling pathways involved.

To ensure the reliability of the reported data acquired from LC-MS/MS, we performed the following steps for data quality control. First, 1%

FDR was performed to control protein identification and quantification before mass spectra data retrieval. Second, unused ≥ 1.3 was set to ensure the confidence of their proteins higher than 95% after database searching. Third, we in turn removed the results identified by reverse database, the proteins with extremely high or low ratios, and the proteins with abnormal quantification between biological repetition and technical repetition. Statistical significance of DAVID bioinformatics enrichment analysis was assumed when p-values were < 0.05.

Meanwhile, multiple statistical values were estimated for the purpose of enrichment evaluation in each term, such as count and percentage of enriched genes, p-value, fold enrichment score, Bonferroni value, and false discovery rate (FDR).

Three anatomically compartments from the normal monkey heart were separated for MS analysis: The vessels, valves, and chambers. To obtain accurate proteome information, all the F I G U R E 1 The schematic and protein distributions of the monkey heart proteome. (a) Schematic depiction of cardiac proteomic analysis of the normal monkey heart tissues, including experimental design (upper panel), protein sample collection, and bioinformatics analysis (lower panel); (b) Bar plot comparing the total number of cardiac proteins in dissimilar heart regions; (c) Venn diagram showing the count of respective and common proteins in each cardiac region; (d) The subcellular localization of total heart proteins in normal monkey; (e) The top 10 terms illustrating the subcellular distribution of cardiac proteins in vessels [Color figure can be viewed at wileyonlinelibrary.com] F I G U R E 2 The bioinformatics analysis of the healthy monkey heart genes. (a) Enrichment network of functional terms from monkey cardiac proteome (p < 0.05). Each dot represents one functional term, and the same color dots stands for terms which have similar biological functions; (b) The pathway annotation network in total heart proteins was classified into 25 signaling pathways (p < 0.05), which mainly contain three major pathway categories including energy metabolism, signal transduction, and disease regulation; (c) Main groups and hierarchical relationship of Reactome pathways. Based on the similarity of their physiological functions, Reactome pathways were grouped into several categories [Color figure can be viewed at wileyonlinelibrary.com] sample preparation was performed by reducing sample loss and increasing quantification accuracy after tissue homogenization in liquid nitrogen. The workflow of this study was presented in Figure   1a 

To gain a better understanding of proteins from normal monkey heart, we analyzed their biological characteristics and signaling pathways enrichment (Supporting Information Supplementary File 2). For an overall assessment of subcellular distribution of total cardiac proteins, the top 10 terms with subcellular localization annotations were used by F I G U R E 3 GO analysis of unique genes from dissimilar heart regions (p < 0.05). (a) The count of vessel-specific genes in monkey heart was highlighted in red; (b) The enrichment of molecular functions in vessel-specific genes; (c) Histogram plot displaying the terms and number of vessel-specific genes in monkey heart; (d) The enriched biological processes network of heart vessel-specific genes; (e-h) The number (e), molecular function categories (f) and terms (g), and biological processes (h) of heart valves-specific proteins were showed, respectively; (i-l) Similar GO visualized network specific for monkey heart chambers-related genes were also displayed including the number (i), molecular function categories (j) and terms (k), and biological processes (l). GO: Gene Ontology [Color figure can be viewed at wileyonlinelibrary.com] the counts of enriched genes, which represent the majority of protein distribution in monkey heart. Among all heart proteins, 36.7% were extracellular proteins, 18.2% were related to the cytoplasm, 12.8% located in mitochondrion, 9.6% functioned in the membrane, 8.9% were distributed in the nucleoplasm, and 7.1% were found in the cytosol F I G U R E 5 Regional signaling pathway analysis of the genes associated with heart diseases. (a) The genes related to heart diseases in cardiac vessels such as CS, SDHA, IDH1, and DLAT worked together to regulate TCA cycle pathway; (b) The disease-related genes in valves, like ACAA1, ALDH9A1, and HADHB, performed their functions in amino acid and fatty acid degradation pathway; (c) The genes in chambers were mainly involved in multiple pathways including glycolysis and propanoate metabolism pathway. TCA: citrate cycle [Color figure can be viewed at wileyonlinelibrary.com] term). Meanwhile, for each term, a series of statistical values such as the count of enriched genes, p-value and Bonferroni value were taken for the purpose of enrichment evaluation. More than half of all biological terms were focused on metabolism-related processes, including ATP metabolic process, organonitrogen compound metabolic process, and small molecule metabolic process. Myocardial contraction-related processes, such as heart contraction, myofibril assembly, and fibril organization, are also thought as an essential part in all biological terms. Besides, other processes like transcription processes also played a vital role in regulating heart function.

To investigate the signaling pathways involved in the corresponding proteins of monkey heart, all identified proteins were submitted to Cytoscape platform for KEGG pathway analysis and functional annotation clustering (FAC). The pathway annotation network in total heart proteins was classified into 25 signaling pathways. Three major types of pathways involving energy metabolism, signal transduction, and disease regulation were considered as the important parts ( Figure 2b and Supporting Information S2). Monkey cardiac proteins were partially enriched into multiple metabolic pathways such as carbon metabolism, fatty acid degradation, and amino acid degradation. Twenty-two cardiac proteins were clustered into the citrate cycle (TCA cycle) pathway, which account for 65.71% of the total proteins in this signaling pathway.

The genetic information associated with each pathway was provided in Supporting Information Table S1 .

To further elucidate the correlation between monkey heart proteins and human pathways, Reactome pathways analysis was used on the basis of the orthologous correspondences informed from the ENSEMBL Compara database. According to physiological similarity and hierarchy, human biological pathways were classified into different groups including immune system, signal transduction, metabolism, disease, and cell cycle and hemostasis ( Figure 2c ). As shown in Figure 2c , 1,330 matched genes were enriched in 1,535 specific pathways, which were accompanied by several statistical values including Entities p value and Entities FDR. Based on the degree of Entities FDR, these genes were mainly involved in several pathway categories such as extracellular matrix organization, muscle contraction, and programmed cell death (Table 1) . More narrowly, most of these genes are enriched in the citric acid cycle, respiratory electron transport, Complex I biogenesis, and other pathways (Supporting Information Table S2 ).

Gene Ontology (GO) analysis of unique proteins in each region was used to elucidate their characteristic functions of heart tissue (Supporting Information Supplementary File 3) . The contribution of main organelles to the amount of cellular protein in each region is similar to the cellular component in total heart proteins. Similar to the cellular component in total heart proteins, more than half of all spatial-specific cardiac proteins were expressed in extracellular exosome, mitochondrion, plasma membrane, cytosol, and endoplasmic reticulum (Figure 1e and Supporting Information S1a,b).

Proteomics analysis of heart vessels mainly detected the genes among aorta and pulmonary artery that does not distinguish between the concave and convex areas. In heart vessels, the region-specific proteins were enriched in multiple molecular functions categories such as oxidoreductases activity, RNA binding, cytoskeletal protein binding, and cell adhesion molecule binding, and were distributed on 12 specific terms (Figure 3a-c) . Further, the BP enrichment network of heart vessel-specific proteins was displayed in Figure 3d . These proteins were implicated in numerous biological functions including aerobic respiration, cellular component assembly, generation of precursor metabolite and energy, and organic acid metabolic process (Supporting Information Table S3 ). For example, up to 41 vesselexpressed genes are involved in regulation of oxidation-reduction process in vivo, including UQCRC2, EIF6, ALDOC, and ADSSL1. In addition, ALDOC, EIF6, and UQCRC2 have the potential to play a significant role in multiple metabolic processes involving nucleoside monophosphate, ribonucleoside triphosphate, and ATP metabolic process.

As shown in Figure 3e -g, the unique proteins in cardiac valves were focused on multiple MF categories involving oxidoreductases activity, RNA binding, cytoskeletal protein binding, and cell adhesion molecule

binding. Our BP analysis shows that most valves-specific proteins contributed to metabolic processes (Supporting Information Table S4 ).

Cellular component assembly, protein activation cascade and organelle organization processes were also prominent in heart valves (Figure 3h) .

At the genetic level, over 100 genes specifically expressed in heart valves are associated with regulation of organelle organization, 75 genes with proteolysis process, and 132 genes with metabolic processes.

In atria and ventricles, the molecular functions of specific proteins mainly were involved in activate various enzymes, and bind related molecules and complexes (Figure 3i-k) . Furthermore, BP analysis of the chambers-specific proteins highlighted the characteristic processes about multiple generating energy and metabolism, response to oxidative stress, and protein complex subunit organization ( Figure 3l and Supporting Information Table S5 ). Among these, up to 20 metabolism-related biological processes and 83 chambersspecific genes are enriched in the Group 7.

The forming protein-protein interactions network is rather intricate because of the counts of unique proteins in each region varied from 400 to 600. To narrow down the list and enhance the readability of the result, the specific proteins of each region were categorized by diverse biological processes in which these proteins participate. The target gene lists, derived from diverse biological processes in each heart region, were imported into STRING and Cytoscape for protein-protein association network analysis, reconstruction, and visualization.

As showed in Figure 4a -e and Supporting Information S3-4, over 20 specific protein interaction networks with diverse biological processes were generated. There is a closely relationship between proteins in the same biological process. For example, the relativity of vessel-specific proteins enriched in aerobic respiration was showed in 

To gain deeper insights into the association between cardiac genes and diseases, we intentionally investigated the crucial gene from protein-protein interaction networks above which are most likely to reflect the biological functions in each heart region. A total of 5,522

GDAs were obtained by matching the core gene lists to the entire GDA data set (Supporting Information Supplementary File 4) . Among these associations, many diseases such as cardiovascular diseases, nutritional and metabolic diseases and endocrine system diseases were highly associated with the physiology and pathology of different heart regions.

From pathophysiological perspective, it would be desirable to further elucidate the relevance between monkey heart genes and cardiac-related illnesses. As displayed as Table 2 , more than 66 genes related to heart diseases (17, 26, and 23 genes in vessels, valves and chambers, respectively) were described in three dissimilar heart regions, which are involved in several cardiac-related illnesses including hypertrophic cardiomyopathy, heart failure, and myocardial infarction.

To illuminate potential differences of the heart diseases-related genes in different heart regions, the gene lists mentioned above were enriched into diverse signaling pathways. The enrichment genes such as CS, SDHA, IDH1, and DLAT in cardiac vessels worked together to regulate the citrate cycle (TCA cycle) pathway ( Figure 5a) . And, the genes in heart valves, like ACAA1, ALDH9A1, and HADHB, were mainly functional in amino acid and fatty acid degradation pathway ( Figure 5b) . However, the genes in atria and ventricles were mainly involved in glycolysis and propanoate metabolism pathway ( Figure 5c ).

As the important part of studying the heart function, the regionresolved proteomics profiling of monkey heart is comparatively lacking. In this study, we established the first spatial-resolved proteomic map of the healthy monkey heart. A total of 3,226 proteins were cataloged in monkey heart, and conveyed multiple information for these genes, including biological characteristics, protein-protein interaction network and genotype-phenotype relationships. Among these genes, more than 60 regional-specific genes are implicated with heart diseases such as hypertrophic cardiomyopathy, heart failure, and myocardial infarction. These genes are mainly enriched in citrate cycle, amino acid degradation, and glycolysis pathway.

Although the orthologous genes among mammals are not fully equal to the genes in human, the homology models in mammals (i.e., mouse, rat, and NHP) are widely used to deduce the function of human genes (Alessio et al., 2017; Gharib & Robinson-Rechavi, 2011 diseases such as severe acute respiratory syndrome (SARS), AIDS disease, and transthyretin (TTR) amyloidosis. The common marmoset with SARS-associated coronavirus infection is suited for understanding SARS pathogenesis, and will allow for efficient testing of vaccines and various treatment modalities (Greenough et al., 2005) .

(SHIV)-infected macaques have become a well-established primate model to study AIDS and pulmonary arterial disease pathogenesis (George et al., 2013; Letvin et al., 1985; Voelkel, Cool, & Flores, 2008) . Aged vervet monkeys have been seen as an ideal animal model of TTR amyloidosis, and have the TTR Ile122 allele, which is well known as a frequent mutation causing the leptomeningeal phenotype of human TTR amyloidosis (Ueda et al., 2012) .

Besides, multiple nonhuman primate models of heart disease afford an opportunity to understand the pathology of human heart disease. Monkey model of myocardial infarction is generated by ischemia followed by reperfusion to study the extent of remuscularization of the infarcted heart by transplanting human cardiomyocyte (Chong et al., 2014) . Acute myocardial infarction is also induced in monkey by left anterior descending artery ligation (Yang et al., 2011; Yoshioka et al., 2005) . Adhesion molecule P-selectin of coronary endothelium is overexpressed in a monkey heart ischemia reperfusion model, which allow stable adherence with leukocytes that provide possibility for detecting myocardial ischemia by molecular imaging of P-selectin with targeted ultrasound contrast agents (Kaufmann et al., 2010; Kaufmann, Lewis, Xie, Mirza-Mohd, & Lindner, 2007; Thomas et al., 2010) . Sensitive proarrhythmia model establishing by a chronic atrioventricular block in cynomolgus monkey possesses essentially the same pathophysiological adaptations with human (Sugiyama, 2008) . CCR5 inhibition observably impacts myocardial viral load and prevents cardiomyocyte functional impairment in SIV-infected macaque model of diastolic dysfunction (Kelly et al., 2014) . CCR5 blockade also attenuates postsurgical stress responses and favorably modulates pathogenic alloimmunity in nonhuman primate cardiac allograft model (Schroder et al., 2007) .

Monkey models of heart disease also serve as a suitable large-animal model to examine diabetic cardiomyopathy (Gong et al., 2013; Qian et al., 2015) , postmenopausal atherosclerosis (Melendez et al., 2015) , and heart failure (Qiu et al., 2008) . Interestingly, the major pathological process of heart disease is different in humans and chimpanzees, although heart disease is common in the two species (Varki et al., 2009) . The commonest cause of heart disease in human arises predominantly from coronary artery atherosclerosis, whereas sudden death of cardiac etiology in chimpanzees is generally related to myocardial fibrosis (Varki et al., 2009) .

Dissecting monkey heart proteins may provide a fresh perspective into understanding heart-related biological functions. In our current study, there is a greater proportion of monkey heart proteins from extracellular exosome (36.7%), plasma membrane (9.6%), and nucleoplasm (8.9%) compared with human heart proteome (Doll et al., 2017) . But, the count of monkey heart proteins and the detected monkey heart proteins in myofibril, mitochondrion, and endoplasmic reticulum less than in healthy human heart (Doll et The region-specific genes of monkey heart enriched in protein-protein interactions above were imported to a gene-disease association analysis. ** These genes associated with heart diseases were listed and involved in one or more cardiac-related illnesses. *** 2017). In addition, the biological processes of identified proteins in normal monkey heart were classified into metabolism and myocardial contraction-related processes. The metabolism processes detected in monkey heart consists mainly of multiple metabolic processes, which regulate the chemical reactions and pathways involving ATP, any low molecular weight molecule, and organonitrogen compound. Heart contraction processes are primarily comprised of regulation of heart contraction, heart rate, and blood circulation. Moreover, fibril organization is mainly implicated with the assembly, arrangement, and disassembly of several fibrils. These major biological processes play an equally important role in the regulation of heart function.

Mitochondrial dysfunction would lower the rate of oxidative phosphorylation and cause energy deficiency, which have been implicated with many diseases, including heart disease and neurological disorder (Fillmore, Mori, & Lopaschuk, 2014; Ryan, Hoek, Fon, & Wade-Martins, 2015; Subramaniam & Chesselet, 2013) . Previous studies showed that several proteins, like ATP5B, CYC1, and COX5A, involved in the oxidative phosphorylation pathway are significantly downregulated in myocardial infarct-like monkey heart (Song et al., 2014) . In addition to these proteins, more than 100 proteins, which are associated with oxidative phosphorylation and neurological disorders pathway, were identified in normal monkey heart tissue. Similarly, Reactome pathways analysis showed that most of these identified proteins are enriched in energy metabolic pathway including TCA cycle, respiratory electron transport and mitochondrial protein import pathway.

The disturbance in the energy metabolism is inextricably linked with cardiac disorders (Lopaschuk, Ussher, Folmes, Jaswal, & Stanley, 2010) .

In the GDA study, various genes from heart regions were implicated with several heart illnesses including hypertrophic cardiomyopathy, heart failure, and myocardial infarction. Among these, vessel-expressed genes like CS, SDHA, IDH1, and DLAT are enriched in TCA cycle pathway. Citrate synthase (CS), as the rate-limiting enzyme of TCA cycle, reflects mitochondrial oxidative capacity in cardiac tissue (de Castro Bras et al., 2014; Jaenisch, Bertagnolli, Borghi-Silva, Arena, & Lago, 2017) . ACAA, ALDH9A1, and HADHB expressed in heart valves can perform these functions in metabolism pathways including amino acid and fatty acid degradation. Besides, the genes in chambers were mostly involved in metabolism pathways including glycolysis and propanoate metabolism pathway.

In our study, we identified some proteins implicated in the signaling pathway of heart disorders such as hypertrophic cardiomyopathy (HCM), arrhythmogenic right ventricular cardiomyopathy (ARVC), and dilated cardiomyopathy (DCM) pathway. Among these proteins, ACTC1, MYBPC, MYL, TNC, TNI, and TPM belonged to sarcomeric or intrasarcomeric protein which enhances Ca 2+ sensitivity of ATPase activity and force production (Gangadharan et al., 2017; Gomez et al., 2014; Sequeira et al., 2013) , PRKCG, and RYR consisted in mitochondria and sarcoplasmic reticulum, respectively, which increases cytosolic free Ca 2+ level, take mutual effort in left ventricle hypertrophy and can potentially lead to heart failure (Rouet-Benzineb et al., 1996; Uenoyama et al., 2010) . In terms of ARVC mechanisms, DSG and DSP, which induce loss of electrical coupling between cardiac myocytes and lead to myocyte cell death (Azaouagh, Churzidse, Konorza, & Erbel, 2011; Pilichou et al., 2006) , and RYR2, which disrupts normal Ca 2+ cycling (Kannankeril et al., 2006; Tang, Tian, Wang, Fill, & Chen, 2012; Yano, Yamamoto, Kobayashi, & Matsuzaki, 2009) , indirectly participate in symptomatic right ventricular arrhythmias mediation. In addition, protein kinase A (PKA) and sarcoplasmic protein RYR and PLN destroy intracellular Ca 2+ homeostasis resulting in further left ventricle chamber dilation and systolic dysfunction (Gramolini et al., 2008; Sun et al., 2010; Yaniv et al., 2011) .

In summary, we have applied an iTRAQ-coupled LC-MS/MS method to identify cardiac region-specific mapping from healthy monkeys. Results from multiple bioinformatics analysis enable us to better understand the relationship between biological function, molecule mechanism, and disease development in cardiac proteins, which may help other researchers to develop comprehensive cognition on physiology and etiology studies regarding cardiac proteins. Furthermore, our study provided a reference map for future in-depth research on cardiac disease-related NHP model and novel biomarkers of cardiac injury. 

The authors declare no conflicts of interest.

http://orcid.org/0000-0003-4936-3411

Ming-Qiang Rong http://orcid.org/0000-0002-4865-1035

Misidentified human gene functions with mouse models: The case of the retinoblastoma gene family in senescence

Hypertension, inflammation and atrial fibrillation

Arrhythmogenic right ventricular cardiomyopathy/dysplasia: A review and update

Valvular heart disease: Aortic regurgitation

Dynamics of cell generation and turnover in the human heart

CluePedia cytoscape plugin: Pathway insights using integrated experimental and in silico data

ClueGO: A cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks

Mechanoresponsive networks controlling vascular inflammation

Citrate synthase is a novel in vivo matrix metalloproteinase-9 substrate that regulates mitochondrial function in the postmyocardial infarction left ventricle. Antioxidants & Redox Signaling

Human embryonic-stem-cell-derived cardiomyocytes regenerate non-human primate hearts

Left atrial enlargement and right ventricular hypertrophy in essential hypertension

Region and cell-type resolved quantitative proteomic map of the human heart

The Reactome pathway knowledgebase

Mitochondrial fatty acid oxidation alterations in heart failure, ischaemic heart disease and diabetic cardiomyopathy

Molecular mechanisms and structural features of cardiomyopathy-causing troponin T mutants in the tropomyosin overlap region

Physiologic changes in a nonhuman primate model of HIV-associated pulmonary arterial hypertension

When orthologs diverge between human and mouse

Endothelial cell dysfunction and the pathobiology of atherosclerosis

Mutation analysis of the main hypertrophic cardiomyopathy genes using multiplex amplification and semiconductor next-generation sequencing

Comparison of the clinical manifestations of type 2 diabetes mellitus between rhesus monkey (Macaca mulatta lasiotis) and human being

Comparative proteomics profiling of a phospholamban mutant mouse model of dilated cardiomyopathy reveals progressive intracellular stress responses

Pneumonitis and multi-organ system disease in common marmosets (Callithrix jacchus) infected with the severe acute respiratory syndrome-associated coronavirus

Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists

Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources

Respiratory muscle training improves diaphragm citrate synthase activity and hemodynamic function in rats with heart failure

Mice with the R176Q cardiac ryanodine receptor mutation exhibit catecholamine-induced ventricular tachycardia and cardiomyopathy

Molecular imaging of the initial inflammatory response in atherosclerosis: Implications for early detection of disease

Detection of recent myocardial ischaemia by molecular imaging of P-selectin with targeted contrast echocardiography

CCR5 inhibition prevents cardiac dysfunction in the SIV/macaque model of HIV

A draft map of the human proteome

In vivo modulation of endothelial polarization by Apelin receptor signalling

A large dataset of protein dynamics in the mammalian heart proteome

Induction of AIDS-like disease in macaque monkeys with T-cell tropic retrovirus STLV-III

Mitral valve disease--morphology and mechanisms

The emerging role of valve interstitial cell phenotypes in regulating heart valve pathobiology

Serum biomarkers reveal long-term cardiac injury in isoproterenol-treated African green monkeys

Myocardial fatty acid metabolism in health and disease

Proteomic analysis of human fetal atria and ventricle

Phosphoproteomics taken to heart

Reactome knowledgebase of human biological pathways and processes

Circadian influence on metabolism and inflammation in atherosclerosis

Beneficial effects of soy supplementation on postmenopausal atherosclerosis are dependent on pretreatment stage of plaque progression

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs

Top-down proteomics reveals concerted reductions in myofilament and Z-disc protein phosphorylation after acute myocardial infarction

Mutations in desmoglein-2 gene are associated with arrhythmogenic right ventricular cardiomyopathy

DisGeNET: A comprehensive platform integrating information on human disease-associated genes and variants

Diastolic dysfunction in spontaneous type 2 diabetes rhesus monkeys: A study using echocardiography and magnetic resonance imaging

Characterization of a novel cardiac isoform of the cell cycle-related kinase that is regulated during heart failure

Protein kinase C isoform expression in normal and failing rabbit hearts

Mitochondrial dysfunction and mitophagy in Parkinson's: From familial to sporadic disease

CCR5 blockade modulates inflammation and alloimmunity in primates

Perturbed length-dependent activation in human hypertrophic cardiomyopathy with missense sarcomeric gene mutations

Echocardiographic reference ranges for sedated healthy cynomolgus monkeys (Macaca fascicularis)

Quantitative proteomics for cardiac biomarker discovery using isoproterenol-treated nonhuman primates

Mitochondrial dysfunction and oxidative stress in Parkinson's disease

Sensitive and reliable proarrhythmia in vivo animal models for predicting drug-induced torsades de pointes in patients with remodelled hearts

Cerebral ischemia elicits aberration in myocardium contractile function and intracellular calcium handling

Abnormal termination of Ca2 + release is a common defect of RyR2 mutations associated with cardiomyopathies

The cardiac valve interstitial cell

Upregulation of coronary endothelial P-selectin in a monkey heart ischemia reperfusion model

Aged vervet monkeys developing transthyretin amyloidosis with the human disease-causing Ile122 allele: A valid pathological model of the human disease

Protein kinase C mRNA and protein expressions in hypobaric hypoxia-induced cardiac hypertrophy in rats

Heart disease is common in humans and chimpanzees, but is caused by different pathological processes

From viral infection to pulmonary arterial hypertension: A role for viral proteins?

Mass-spectrometry-based draft of the human proteome

Proteomics analysis of human pericardial fluid

Electrocardiographic characterization of rhesus monkey model of ischemic myocardial infarction induced by left anterior descending artery ligation

Ca2 + -regulated-cAMP/PKA signaling in cardiac pacemaker cells links ATP supply to demand

Role of ryanodine receptor as a Ca(2)( + ) regulatory center in normal and failing hearts

Repair of infarcted myocardium mediated by transplanted bone marrow-derived CD34 + stem cells in a nonhuman primate model

Regionresolved proteomics profiling of monkey heart

We only displayed the five major heart-related disorders for one gene when exceeding this limit. HU ET AL. | 13731