Skip to main content

Plasma proteomic and autoantibody profiles reveal the proteomic characteristics involved in longevity families in Bama, China



Chinese Bama Yao Autonomous County is a well-known longevity region in the world. In the past 30 years, population and genome studies were undertaken to investigate the secret of longevity and showed that longevity is the result of a combination of multiple factors, such as genetic, environmental and other causes. In this study, characteristics of the blood plasma proteomic and autoantibody profiles of people from Bama longevity family were investigated.


Sixty-six plasma donors from Chinese Bama longevity area were recruited in this study. Thirty-three offsprings of longevous families were selected as case studies (Longevous group) and 33 ABO (blood type), age, and gender-matched subjects from non-longevous families were selected as controls (Normal group). Each group contains 3 biological replicates. Tandem mass tag-based proteomic technique was used to investigate the differentially expressed plasma proteins between the two groups. The auto-reactive IgG antibody profiles of the 3 pooled samples in each group were revealed by human proteome microarrays with 17,000 recombinant human proteins.


Firstly, 525 plasma proteins were quantified and 12 proteins were discovered differentially expressed between the two groups. Secondly, more than 500 proteins were recognized by plasma antibodies, 14 proteins ware differentially reacted with the autoantibodies in the two groups. Bioinformatics analysis showed some of the differential proteins and targeted autoantigens were involved in cancer, cardiovascular disease and immunity.


Proteomic and autoantibody profiles varied between the offspring of longevous and normal families which are from the same area and shared the same environmental factors. The identified differences were reported to be involved in several physiological and pathological pathways. The identified proteins will contribute to a better understanding of the proteomic characteristics of people from Bama longevous area and a revelation of the molecular mechanisms of longevity.


Human longevity is a complex phenotype with a significant familial component [1, 2]. Genome studies were taken and indicates that gene polymorphisms contribute to human longevity in the past years [3,4,5]. However, only examining the genome profiles does not give an integrated view, as proteome may not be totally accurately predicted by gene profiles due to several factors [6]. Furthermore, proteins are major components of the tissues and cells with diverse regulatory, enzymatic and structural functions [7, 8]. In the last 20 years, proteomics analysis has been widely used to find biomarkers, proteins that change in concentration or state in associations with a specific biological process or disease [9,10,11,12]. Blood plasma, which include thousands of proteins, is a predominant sample used for human proteomics study.

Bama Yao Autonomous County, a remote and mountainous county located in the northwest of Guangxi Zhuang Autonomous Region (Guangxi province) of China, is known for its longevity [13, 14]. The longevity township standard set by the United Nation is 7 centenarians among 100,000 people [13]. The population size and centenarian rate in Bama County are around 240,000 and 30/100,000 respectively, according to the National Population Census of China in the past decades [14, 15]. In 2003, Bama was awarded the certificate of The World’s Fifth Longevity Township by the International Natural Medical Association [13]. In the past decades, the people of Bama have been utilized for aging/longevity cohort [16,17,18,19].

To investigate the secret of Bama-induced longevity at the proteome levels, in this study, 66 plasma donors from Bama longevity hotspot which shared the same environmental factors were recruited and divided into two groups, offsprings of longevous families and offsprings of non-longevous families. Tandem mass tag (TMT)-based approach was used to systematically investigate the expression of proteins in the plasma of the two groups. The differentially expressed proteins (DEPs) were revealed by bioinformatics analysis.

In most proteomic profile studies, the high abundance proteins in plasma (include albumin, immunoglobulin G, transferrin, etc.) were removed to greatly improve detection of lower abundance protein [10, 20,21,22,23]. However, IgG antibodies are an important part of the immune system, which plays a vital role in physiological conditions in healthy individuals. Autoantibody or self-reactive antibody refers to antibody that are present in the blood plasma of individuals in the absence of deliberate immunization with the target antigen, it reacts with human own proteins or tissues [24]. It is generally believed that autoantibodies are critical for many physiological and pathological processes [25, 26]. In this study, the HuProt arrays, comprised of ~ 17,000 individually purified human proteins, were also employed to investigate the IgG autoantibody repertoires of the two groups of people from Bama.


Sample preparation

A questionnaire was sent to healthy plasma donors in the Bama longevity area, and the number of long-lived immediate family members (≥ 85 years old) within three generations in the donors were investigated. A donor with at least three long-lived immediate family members was defined as the offspring of longevous family, and donors without an immediate family member ≥ 85 years old within three generations were offsprings of non-longevity families.

In total, 33 offsprings of longevous families were recruited as case studies (Longevous group). And also 33 ABO-, age- and gender-matched offsprings of non-longevity families were recruited as controls (Normal group). The inclusion criteria of the participators were recruited as previously described and apheresis plasma of the participators was used in this study [27]. Briefly, all participators were over 18 years of age, healthy and unrelated. People who had prior history of thrombus or hemorrhage, usage of oral anticoagulation therapy, hepatic disease, HIV infection, pregnancy, diabetes, cardiopathy, hypertension, renal insufficiency, respiratory diseases and others were excluded from this study. Pooled plasma samples were generated by combining equal volumes of each 11 individual plasma samples from each group. The basic information of the volunteers was shown in Table 1. This study was approved by the Ethics Committee of the Institute of Blood Transfusion (No. 201716). In accordance with the Declaration of Helsinki, all participants gave informed consent prior to their entering the study.

Table 1 The basic information of the study cohort

TMT-based quantitative proteomic analysis

Tandem mass tagging labeling and HPLC fractionation

The high abundance proteins that account for approximately 80% of the total plasma protein concentration were removed from the pooled plasma samples, using a ProteoMiner™ Protein Enrichment Introductory Large-Capacity Kit (Bio-Rad, Richmond, USA) according to the manufacturer’s instructions. After being reduced with 5 mM DL-dithiothreitol for 30 min at 56 °C and alkylated with 11 mM iodoacetamide for 15 min in darkness at room temperature, the protein sample were digested with trypsin at a trypsin/protein mass ratio of 1:50 overnight in the first stage and 1:100 for 4 h in the second stage. Approximately 100 µg protein for each sample was digested with trypsin for the following experiments. Tryptic peptides were desalted by Strata X C18 SPE column (Phenomenex, CA, USA) and labeled with a 6-plex TMT kit (Thermo Fisher Scientific, CA, USA) according to manufacturer’s instructions. The 6 labeled samples then fractionated by high pH reverse-phase HPLC using 300Extend C18 column (5 µm particle size, 4.6 mm ID, 250 mm length, Agilent, CA, USA). Peptides were eluted with a gradient of acetonitrile (8% to 32%) in 10 mM ammonium bicarbonate (pH 9.0) over 60 min. At last, 60 fractions were collected, and they were combined into 18 fractions and dried by vacuum centrifuging.

LC–MS/MS and bioinformatics analysis

Samples were measured using LC–MS instrumentation consisting of an EASYnLC 1000 ultra-high-pressure system coupled via a nano-electrospray ion source to a Q Exactive™ Plus (Thermo Fisher Scientific). Purified peptides were separated on a reversed-phase analytical column (150 mm length, 75 μm ID). For each LC–MS/MS analysis, the electrospray voltage applied was 2.0 kV. Intact peptides were detected in the orbitrap at a resolution of 70,000, and ion fragments were detected in the orbitrap at a resolution of 17,500. In the MS survey scan, a data-dependent mode with an automatic alteration (1 MS scan followed by 20 MS/MS scans) was used for the top 20 precursor ions above a threshold ion count of 5 × 104 with 30 s dynamic exclusion. Automatic gain control was used to prevent overfilling of the orbitrap; 5 × 104 ions were accumulated for generation of MS/MS spectra. For MS scans, the m/z scan range was 350 to 1800, and the fixed first mass was set as 100 m/z.

MS/MS data was analyzed by MaxQuant software with integrated Andromeda search engine (Version The mass tolerance for precursor ions was set as 20 ppm in First search and 5 ppm in Main search, and that for fragment ions was set as 0.02 Da. False discovery rate (FDR) was adjusted to < 1% and minimum score for peptides was set > 40. Only unique peptides were used for protein quantification. And for protein quantification method, TMT 6-plex was selected in Mascot. FDR was adjusted to < 1% at protein, peptide and peptide spectrum match (PSM) level. Gene Ontology (GO) annotation (the biological process, cellular component and molecular function) of the differentially expressed proteins were derived from the UniProt-GOA database ( PSORT/PSORT II and SubLoc ( were used to predict subcellular localization of all identified DEPs. A p value < 0.05 was used as the threshold to determine the significant enrichments of GO annotation.

Plasma autoantibody profiling using the human proteome microarrays

Human proteome microarrays

The HuProt™ V3.1 microarray (CDI Laboratories, Inc. Mayaguez, USA) composed of about 17,000 recombinant human proteins was used in this study. All of the recombinant human proteins were generated by the Saccharomyces cerevisiae expression system and carried an N-terminal glutathione S-transferase (GST) tag. Each protein was spotted in duplicate [11, 28]. The microarray was firstly blocked with blocking buffer (5% BSA in 0.1 mol/l TBS-T, pH = 7.5) at room temperature for 1.5 h. Then it was washed with with TBS-T and then Milli-Q water for 10 min respectively. The six pooled plasma samples of two groups were diluted 1: 2000 in 0.01 mol/l Tris-Buffered saline containing 0.1% (v/v) Tween 20 detergent (TBS-T, pH = 7.5) with 5% BSA. Five ml of diluted plasma were added and incubated with microarray under a glass coverslip at room temperature for 1 h. After rinsing with TBS-T and then Milli-Q water, Cy3-conjugated anti-human IgG diluted 1:2000 in TBS-T was added to the washed microarray and incubated in darkness at room temperature for 1 h. The microarray was then washed again as described above. After drying by centrifugation (ChipMate PMC-082, Tomy, Japan), microarray was subjected for scanning with the Axon GenePix 4000B microarray scanner (Molecular Devices, CA, USA) using an excitation wavelength of 523 nm and an emission filter specifically for cy3 fluorescent dye. The optical signal was detected by a photomultipliers tube (PMT).

Human proteome microarray data analysis

Human proteome microarray data analysis GenePix Pro 6.0 software (Molecular Devices, CA, USA) was used to obtain microarray signal intensity. The foreground median intensity (F median) and background median intensity (B median) for each spot on the microarray images were acquired. To quantify the signal intensity for each spot, we calculated the signal intensity for each spot, which was defined as the foreground median intensity divided by its local background median intensity (F median/B median). The intra-array signal intensity normalization and positive hits identification were taken as described in the previous reports [29, 30]. Briefly, the signal intensities of the spots within a protein microarray were normalized by setting the median intensity of that microarray equal to 1. To identify autoantibodies that bind to a protein on microarray (positive hits), an intensity cutoff value needed to be assigned for each microarray. We used I to denote the normalized intensity of a spot (control spots omitted) on the protein microarray. And Meani denote the average of I (control spots omitted) on a protein microarray. The standard deviation (SD) away from the mean of the signal intensities for all the spots in a microarray was calculated. A cutoff was defined as Meani + 6SD, and spots producing a signal intensity greater than the cutoff were identified as “positive hits.” Moreover, each protein was printed in duplicate on a microarray, proteins and controls were considered as a positive feature only when both of their duplicate spots were simultaneously judged as a positive hit. T-test was chosen to assess the differential significance for each protein between the two groups based on signal value. The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 were used for the functional category analyzation of the identified autoantigens, subcellular localization prediction and GO annotation of the different autoantigens between groups.


Differentially expressed proteins revealed by quantitative proteomic analysis

A total of 598 proteins were identified using MaxQuant, among which 525 proteins were quantified at more than a 95% CI, and proteins with no less than two unique peptides were considered to be a positive identification (see Additional file 1). For comparison between the longevous and normal groups, a protein featuring a fold change of > 1.2 or < 0.83 and a p value of < 0.05 was regarded as DEP. The identified DEPs were showed by a volcano plot and a hierarchical clustering heat map (Fig. 1a, b), and described in Table 2. 17 DEPs were identified based on the criteria, which were corresponding to 12 different proteins, as some results corresponded to different type or different isoforms of the same protein (DEPs No. 3–7, No. 10 and 11). Among these 12 proteins, 5 were significantly down-regulated (LPA, DEFA3, immunoglobulin, BLVRB and PDIA3), 7 were accumulated (TPM, GBA, DSG2, SRC, CHGA, ITGB3 and TAGLN2) in the longevous group A than normal group.

Fig. 1

Differentially expressed proteins in the longevous and normal groups. TMT-based proteomic technique was used to study the protein characteristics in longevity families. Differentially expressed proteins and were defined as having fold changes > 1.2 (> 1.20 increased or < 0.83 decreased). A two tail Student’s T-test was performed and a p-value < 0.05 was considered to be statistically significant. a The volcano plot shows the up- (red) or down-regulated (green) proteins between the two groups. b Hierarchical clustering of differentially expressed proteins. The color scale bar locates in the bottom, and green and red indicate decreased and increased levels of the identified poteins, respectively. A1, A2 and A3, 3 replicates of the longevous group; B1, B2 and B3, 3 replicates of the normal group. c Biological processes, d cellular components and e molecular functions of GO annotation based function classification of the differentially expressed proteins

Table 2 Differentially expressed proteins in the longevous and normal groups detected by TMT-based proteomic technique

For an overview of the aforementioned DEPs, ontological functions were explored using UniProt-GOA. The ontology of biological processes indicated that these proteins were involved in 16 biological processes including single-organism process, response to stimulus, signaling, multicellular organismal process, multi-organism process, cellular process, cellular component organization or biogenesis, biological regulation, developmental process, reproduction, immune system process, biological adhesion, locomotion, localization, metabolic process and cell killing (Fig. 1c). For cellular components, Fig. 1d illustrated that the proteins were main membrane, cell, extracellular region, organelle and cell junction components. The ontology of molecular functions suggested that the represented function of the DEPs were binding and catalytic activity (Fig. 1e).

Plasma autoantibody identification by human proteome microarray

The HuProt arrays were employed to profile the immune-reactive autoantigens of the two groups of people from Bama longevity area. The 6 samples were individually incubated on the HuProt arrays, followed by incubation with a fluorescently labeled secondary anti-human IgG antibody to identify immune-reactive autoantigens. Using a relatively stringent cutoff value, human protein positively reacted to each sample were identified as the autoantigens that reacted with the autoantibodies in the 6 pooled plasma samples of the two groups. Finally, 652, 695 and 652 proteins on the human proteome microarray were recognized by the IgG in the 3 samples of the longevous group respectively. Among these proteins, 553 were recognized by all the 3 samples (Fig. 2a and Additional file 2). For the normal group, 668, 585 and 622 autoantigens were identified that directly reacted with the 3 samples respectively, and 508 proteins were recognized by all the samples of the normal group (Fig. 2b and Additional file 2).

Fig. 2

Venn diagram of the number of the immune-reactive autoantigens in the Longevous and Normal groups identified by human proteome microarray. a Longevous group, A1, A2 and A3, 3 replicates of offsprings of Bama longevous families; b Normal group. B1, B2 and B3, 3 replicates of offsprings of Bama non-longevous families

Functional category of the aforementioned 553 (in the longevous group) and 508 (in the normal group) autoantigens were then analyzed. The top 10 of the count number of terms were displayed by pie chat (Fig. 3). The result indicated that the functional category of the antigens recognized by the two groups are almost the same: more than 200 identified antigens were involved in polymorphism (284/258, 13%), alternative splicing (281/258, 13%) and phosphoprotein (278/251, 13%). And other significant function groups included response to cytoplasm (9–10%), nucleus (7%), acetylation (6–7%), coiled coil (5%), Ubl conjugation (3%), nucleotide-binding (2–3%) and transferase (2%), etc.

Fig. 3

Functional category analyzation of the indentified proteins (autoantigens) recognized by plasma autoantibodies, using human proteome microarrays. a 553 autoantigens recognized by all the 3 samples of longevous group; b 508 autoantigens recognized by all the 3 samples of normal group

For the determination of the specifically and differentially expressed autoantibodies in the offsprings from longevous and non-longevous families from Bama, all antigens recognized by IgG in plasma samples of the longevous group were compared with that of the normal group. A protein recognized in all three samples in one group and not recognized by any sample in another group was defined as a specific autoantigen. And an antibody with its antigen–antibody reactive signal intensity of featuring a fold change of > 1.2 or < 0.83 (p < 0.05) between the two groups was regarded as a differentially expressed autoantibody. The differentially recognized autoantigens were showed by a volcano plot and a hierarchical clustering heat map (Fig. 4a, b), and described in Table 3. 14 autoantigens were identified, including 6 autoantigens specifically reacted with the sample of the longevous group (CASC4, DPY30, PMVK, HSPA2, IFIT3, SCRN2), 1 specifically reacted with the sample of the normal group (ZNF207), and 7 were differentially recognized by the two groups (ARL2BP, GCK, ING3, NXPH2, RUVBL2, TPD52L1, ORMDL2). Furthermore, 12 of the 14 identified antigens were specifically or more intensively reacted with antibodies in the longevous group.

Fig. 4

Identified proteins differentially recognized by plasma autoantibodies between the longevous and normal groups. Human proteome microarrays were used to study the autoantibody characteristics in longevity families. A protein recognized in all three samples in one group and not recognized by any sample in another group was defined as a specifically autoantigen that represented a specifically expressed autoantibody. And an antibody with its antigen–antibody reactive signal intensity of featuring a fold change of > 1.2 (> 1.20 increased or < 0.83 decreased) between the two groups was regarded as a differentially expressed autoantibody. A two tail Student’s T-test was performed and a p-value < 0.05 was considered to be statistically significant. a The volcano plot shows the up- (red) or down-regulated (green) autoantibodies between the two groups. b Hierarchical clustering of specifically and differentially expressed autoantibodies. The color scale bar locates in the bottom, and green and red indicate decreased and increased levels of the identified autoantibodies, respectively. A1, A2 and A3, 3 replicates of the longevous group; B1, B2 and B3, 3 replicates of the normal group. c Biological processes, d cellular components and e molecular functions of GO annotation based function classification of the specifically and differentially expressed autoantibodies

Table 3 Identified proteins (autoantigens) differentially recognized by plasma autoantibodies of the longevous and normal groups, using human proteome microarrays

GO annotation of these specific and differential autoantigens were then explored. The ontology of biological process indicated that these proteins were main involved in 6 processes (Fig. 4c): ING3 and RUVBL2 were participated in the processes of histone H2A acetylation, histone H4 acetylation and regulation of growth, GCK and ARL2BP were related to the regulation of insulin secretion, ING3, RUVBL2, ZNF207 and RUVBL2 were involved in transcription or regulation of transcription. The cellular components ontology showed the identified antigens were located in cytoplasm, nucleus, cytosol, nucleoplasm, extracellular exosome, mitochondrion, nucleolus, membrane, integral component of membrane, Golgi apparatus, etc. (Fig. 4d). The main molecular functions of these antigens were binding, including identical protein binding, ATP binding, zinc ion binding, unfolded protein binding, etc. In addition, DPY30 and TPD52L1 had protein homodimerization activity (Fig. 4e).


Human longevity was considered to be the result of numerous interacting factors including genetic, environmental and behavioural components [31,32,33]. Several studies suggested that about 20–50% of the variation in human lifespan is accounted for by genetic factors [2, 34, 35]. A moderate clustering of extreme longevity within families was also reported [36, 37]. In the last 30 years, genome and transcriptome studies were employed and many candidate genes have been investigated for putative associations with human survival or longevity. The main categories of the genes are involved in metabolic and immune systems [5]. But little is known about the proteomic profiles involved in longevity. In addition to the whole genome association studies, proteomic analysis can be taken to investigate the molecular mechanism of longevity more accurately, as proteomic profiles do not accurately predicted by transcriptome profiles [6, 38].

For a further investigation of the potential mechanisms of longevity at the protein level, people from a well-known longevous area (Bama, China) were recruited to participate in this study. Offsprings of longevous families were recruited as case studies and ABO-, age- and gender-matched subjects from non-longevous families were selected as controls. TMT-based proteomic technique and human proteome microarrays were used to study the protein and autoantibody characteristics in longevity families. Differentially expressed proteins and autoantibodies were defined as having fold changes > 1.2 and < 0.83, respectively. In addition, a two tail Student’s T-test was performed and a p value < 0.05 with no strategy of correction for multiple testing was considered to be statistically significant based on the previous proteomic and microarray studies [30, 39,40,41,42].

Differentially expressed proteins revealed by quantitative proteomic analysis

By using TMT-based quantitative proteomic analysis, a set of low-abundance plasma proteins were quantified and 12 proteins were found that differentially expressed between the longevous and normal groups. Identified differences indicate involvement of several physiological pathways, such as immunity, metabolism, cell adhesion and signaling. Lipoprotein(a) (LPA) is known as a risk factor for atherosclerotic diseases such as coronary heart disease and stroke [43, 44]. Genetic and epidemiologic studies also showed that individuals without LPA or with very low LPA levels seem to be healthy [45], genes related to cardiovascular health may be implicated in exceptional longevity [5]. Interestingly, LPA was found to be down-regulated most in offsprings from longevous families in our results. Two components involved in immunity were found to be down-regulated in the longevous group : defensin 3 (DEFA3) is expressed primarily in neutrophils as well as in NK cells and certain T-lymphocyte subsets, plays a role in phagocyte-mediated host defense and is considered as a part of the innate immune response [46]. Protein disulfide-isomerase (PDIA3) is part of the major histocompatibility complex (MHC) class I peptide loading complex, which is essential for formation of the final antigen conformation and export from the endoplasmic reticulum to the cell surface [47]. PDIA3 also involved in cytokine-dependent signal transduction, overexpression of it is reported linked to renal fibrosis [48, 49]. Biliverdin reductase B (BLVRB), a NADPH-dependent oxidoreductase, was also low expressed in the offsprings of longevous families. BLVRB is considered to be a critical regulator of cellular redox, as it can physiologically regenerates bilirubin in a catalytic cycle [50].

Several proteins were found to be up-regulated in offsprings from longevous families. Integrin beta 3 (ITGB3, CD61) are known to participate in cell adhesion as well as cell-surface-mediated signaling [51, 52], chromogranin A (CHGA) is widely expressed by neuroendocrine cells, induces and promotes the generation of secretory granules [53]. Proto-oncogene tyrosine-protein kinase Src may play a role in the regulation of embryonic development and cell growth, an elevated level of activity of this protein is suggested to be linked to cancer progression by promoting other signals [54]. Desmoglein-2 (DSG2) is localized to desmosome structures at regions of cell–cell contact and functions to structurally adhere adjacent cells together. Tropomyosin (TMP) is involved in the cytoskeleton of non-muscle cells.

Several proteomic studies have been reported to reveal the age-related protein changes in plasma, cerebral and serum [55,56,57,58]. One study that was conducted in plasma of women enrolled in the TwinsUK study. SOMAscan protein data of 206 female twins were analyzed. Thirteen age‐associated proteins were discovered. Some of them were previously been described to related to cardiovascular disease, apoptosis and aging-associated neuroinflammation [55]. Another SOMAscan proteomic study identified 217 age‐associated plasma proteins,which were main involved in blood coagulation, chemokine and inflammatory pathways, axon guidance, peptidase activity, and apoptosis [56]. 12 of the 13 age-associated proteins reported in TwinsUK study were confirmed to be associated with age in the latter study. In our study, none of the differentially expressed plasma proteins between the longevity and normal group were found to be consistent with the results of the two studies. It is difficult to directly compare the results from the two studies and our work. Because in this study two groups of age-matched individuals were enrolled. And the subjects were grouped based on family longevity, not age.

Plasma autoantibody identification by human proteome microarray

Human plasma contains a large number of autoantibodies that recognize self proteins in healthy individuals, and the autoreactive repertoires are predominantly selected during fetal life [25]. Several functions have been proposed for autoantibody, such as defense against infection, clearance of aging cells, anti-tumoral surveillance and anti-inflammatory activity [25, 26]. In the second part of our study, by using human proteome microarray, more than 500 proteins were found to be recognized by plasma autoantibodies. 14 proteins ware differentially reacted with the autoantibodies in the longevous and normal groups. In addition, we found 12 proteins have higher affinity with antibodies or specifically react with antibodies in plasma of longevous group. Among the 12 proteins, 5 are reported to be related to cancer: CASC4 is a cancer susceptibility candidate, reduction of the expression of CASC4 has been shown to suppress the proliferation of gastric cancer [59]; tumor protein D52-like 1 (TPD52L1) is a cell cycle-regulated protein and maximally expressed at the G2-M transition in breast cancer cells [60]; secernin 2 (SCRN2) is found to be a tumor-associated antigen in gastric cancer [61]; inhibitor of growth family member 3 (ING3) is a member of the ING tumor suppressor family, overexpression of this gene has been shown to inhibit cell growth and induce apoptosis [62]; RuvB like AAA ATPase 2 (RUVBL2) plays roles in transcriptional regulation and DNA repair, it also concerns cancer-related processes via interaction with oncogenic transcription factors [63, 64]. Other autoantigens were also involved in some important biological processes: glucokinase (GCK) and phosphomevalonate kinase (PMVK) are metabolic enzymes that play important role in the carbohydrate metabolism and mevalonate pathway respectively. Histone methyltransferase complex regulatory subunit (DPY30) is located on the nucleus and acted as a transcription regulator [65]. Heat shock protein (HSPA2) protects cells from oxidative stress and inhibits apoptosis [66]. And interferon induced protein (IFIT3) takes part in cell signaling associated with the immune system [67]. In addition, only two proteins, a Zinc finger protein (ZNF207) and an endoplasmic reticulum-localized transmembrane protein (ORMDL2), were found specifically and more intensively reacted with IgG in the normal group.


By using TMT-based proteomic technique and human proteome microarray, the differences of proteomics and autoantibody profiles between offsprings from longevity and non-longevity families in Chinese Bama longevous area were investgated. As a result, 12 proteins and 14 autoantibodies were discovered differentially expressed between the longevous and normal groups. Most of the differential proteins and targeted autoantigens were reported playing critical roles in some physiological and pathological processes. These findings will be contributed to a better understanding of the proteomic characteristics of people from Bama longevous area and a revelation of the mystery of longevity. Further studies are needed to study the effect and molecule mechanism of the differential expressed proteins and autoantibodies on longevity.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.



differentially expressed protein


false discovery rate


gene ontology


glutathione S-transferase


immunoglobulin G


major histocompatibility complex


tandem mass tag


  1. 1.

    Gudmundsson H, Gudbjartsson DF, Kong A, Gudbjartsson H, Frigge M, Gulcher JR, et al. Inheritance of human longevity in Iceland (reprinted from European journal of human genetics, vol 8, pg 743–749, 2000). Eur J Hum Genet. 2017;25:S45–50.

    Article  Google Scholar 

  2. 2.

    Christensen K, Johnson TE, Vaupel JW. The quest for genetic determinants of human longevity: challenges and insights. Nat Rev Genet. 2006;7(6):436.

    CAS  Article  Google Scholar 

  3. 3.

    Willcox BJ, Donlon TA, He Q, Chen R, Grove JS, Yano K, et al. FOXO3A genotype is strongly associated with human longevity. Proc Natl Acad Sci. 2008;105(37):13987–92.

    CAS  Article  Google Scholar 

  4. 4.

    Sebastiani P, Gurinovich A, Bae H, Andersen S, Malovini A, Atzmon G, et al. Four genome-wide association studies identify new extreme longevity variants. J Gerontol Ser A Biomed Sci Med Sci. 2017;72(11):1453–64.

    CAS  Article  Google Scholar 

  5. 5.

    Revelas M, Thalamuthu A, Oldmeadow C, Evans T-J, Armstrong NJ, Kwok JB, et al. Review and meta-analysis of genetic polymorphisms associated with exceptional human longevity. Mech Ageing Dev. 2018.

    Article  PubMed  Google Scholar 

  6. 6.

    Poverennaya EV, Ilgisonis EV, Ponomarenko EA, Kopylov AT, Zgoda VG, Radko SP, et al. Why are the correlations between mRNA and protein levels so low among the 275 predicted protein-coding genes on human chromosome 18? J Proteome Res. 2017;16(12):4311–8.

    CAS  Article  Google Scholar 

  7. 7.

    Boeckmann B, Blatter MC, Famiglietti L, Hinz U, Lane L, Roechert B, et al. Protein variety and functional diversity: Swiss-Prot annotation in its biological context. C R Biol. 2005;328(10–11):882–99.

    CAS  Article  Google Scholar 

  8. 8.

    Jakob U, Kriwacki R, Uversky VN. Conditionally and transiently disordered proteins: awakening cryptic disorder to regulate protein function. Chem Rev. 2014;114(13):6779–805.

    CAS  Article  Google Scholar 

  9. 9.

    Savitski MM, Reinhard FB, Franken H, Werner T, Savitski MF, Eberhard D, et al. Tracking cancer drugs in living cells by thermal profiling of the proteome. Science. 2014;346(6205):1255784.

    Article  Google Scholar 

  10. 10.

    Geyer PE, Kulak NA, Pichler G, Holdt LM, Teupser D, Mann M. Plasma proteome profiling to assess human health and disease. Cell Syst. 2016;2(3):185–95.

    CAS  Article  Google Scholar 

  11. 11.

    Hu C, Pan J, Song G, Wen X, Wu Z, Chen S, et al. Identification of novel biomarkers for behcet disease diagnosis using HuProt array approach. Mol Cell Proteom. 2016;M116:061002.

    Google Scholar 

  12. 12.

    Ngo D, Sinha S, Shen D, Kuhn EW, Keyes MJ, Shi X, et al. Aptamer-based proteomic profiling reveals novel candidate biomarkers and pathways in cardiovascular disease clinical perspective. Circulation. 2016;134(4):270–85.

    CAS  Article  Google Scholar 

  13. 13.

    Liu Y, Nie L, Liao B. The recreational value of Bama in China: one of the five world’s longevity townships. Bus Manag Res. 2012;1(4):141.

    Article  Google Scholar 

  14. 14.

    Xiao Z, Xu Q, Yuan Y. Solving the mystery of the status and longevity of centenarians in Bama. Chin J Popul Sci. 1996;8:385–94.

    CAS  PubMed  Google Scholar 

  15. 15.

    Liang G, Qin J, He M, Li X, Shi L, Huang Q, et al. Genetic polymorphisms of PON1 in longevous people in Bama county of Guangxi province. Chin J Gerontol. 2012;32(7):1338–40.

    CAS  Google Scholar 

  16. 16.

    Ge L, Wu H-Y, Pan S-L, Huang L, Sun P, Liang Q-H, et al. COMT Val158Met polymorphism is associated with blood pressure and lipid levels in general families of Bama longevous area in China. Int J Clin Exp Pathol. 2015;8(11):15055.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Pan S, Wang F, Lu Z, Liu C, Hu C, Luo H, et al. Cholesteryl ester transfer protein Taq IB polymorphism and its association with serum lipid levels and longevity in Chinese Bama Zhuang population. Lipids Health Dis. 2012;11(1):26.

    CAS  Article  Google Scholar 

  18. 18.

    Yang X, Wang X, Yao H, Deng J, Jiang Q, Guo Y, et al. Mitochondrial DNA polymorphisms are associated with the longevity in the Guangxi Bama population of China. Mol Biol Rep. 2012;39(9):9123–31.

    CAS  Article  Google Scholar 

  19. 19.

    Luo C-Y, Liu C-W, Ge L, Pang G-F, Yang M, Hu C-Y, et al. PPARD + 294C overrepresentation in general and long-lived population in China Bama longevity area and unique relationships between PPARD + 294T/C polymorphism and serum lipid profiles. Lipids Health Dis. 2015;14(1):17.

    Article  Google Scholar 

  20. 20.

    Raida M, Schulz-Knappe P, Heine G, Forssmann W-G. Liquid chromatography and electrospray mass spectrometric mapping of peptides from human plasma filtrate. J Am Soc Mass Spectrom. 1999;10(1):45–54.

    CAS  Article  Google Scholar 

  21. 21.

    Georgiou HM, Rice GE, Baker MS. Proteomic analysis of human plasma: failure of centrifugal ultrafiltration to remove albumin and other high molecular weight proteins. Proteom Int Ed. 2001;1(12):1503–6.

    CAS  Article  Google Scholar 

  22. 22.

    Echan LA, Tang HY, Ali-Khan N, Lee K, Speicher DW. Depletion of multiple high-abundance proteins improves protein profiling capacities of human serum and plasma. Proteomics. 2005;5(13):3292–303.

    CAS  Article  Google Scholar 

  23. 23.

    Muenchhoff J, Poljak A, Thalamuthu A, Gupta VB, Chatterjee P, Raftery M, et al. Changes in the plasma proteome at asymptomatic and symptomatic stages of autosomal dominant Alzheimer’s disease. Sci Rep. 2016;6:29078.

    CAS  Article  Google Scholar 

  24. 24.

    Coutinho A, Kazatchkine MD, Avrameas S. Natural autoantibodies. Curr Opin Immunol. 1995;7(6):812–8.

    CAS  Article  Google Scholar 

  25. 25.

    Lacroix-Desmazes S, Kaveri SV, Mouthon L, Ayouba A, Malanchère E, Coutinho A, et al. Self-reactive antibodies (natural autoantibodies) in healthy individuals. J Immunol Methods. 1998;216(1–2):117–37.

    CAS  Article  Google Scholar 

  26. 26.

    Siloşi I, Siloşi CA, Boldeanu MV, Cojocaru M, Biciuşcă V, Avrămescu CS, et al. The role of autoantibodies in health and disease. Rom J Morphol Embryol. 2016;57(Supplement 2):633–8.

    PubMed  Google Scholar 

  27. 27.

    Wang Z, Dou M, Du X, Ma L, Sun P, Cao H, et al. Influences of ABO blood group, age and gender on plasma coagulation factor VIII, fibrinogen, von Willebrand factor and ADAMTS13 levels in a Chinese population. PeerJ. 2017;5:e3156.

    Article  Google Scholar 

  28. 28.

    Yang L, Wang J, Li J, Zhang H, Guo S, Yan M, et al. Identification of serum biomarkers for gastric cancer diagnosis using a human proteome microarray. Mol Cell Proteom. 2015;M115:051250.

    Google Scholar 

  29. 29.

    Hu S, Xie Z, Onishi A, Yu X, Jiang L, Lin J, et al. Profiling the human protein-DNA interactome reveals ERK2 as a transcriptional repressor of interferon signaling. Cell. 2009;139(3):610–22.

    CAS  Article  Google Scholar 

  30. 30.

    Hu C, Huang W, Chen H, Song G, Li P, Shan Q, et al. Autoantibody profiling on human proteome microarray for biomarker discovery in cerebrospinal fluid and sera of neuropsychiatric lupus. PLoS One. 2015;10(5):e0126643.

    Article  Google Scholar 

  31. 31.

    Christensen K, Vaupel JW. Determinants of longevity: genetic, environmental and medical factors. J Intern Med. 1996;240(6):333–41.

    CAS  Article  Google Scholar 

  32. 32.

    Passarino G, De Rango F, Montesanto A. Human longevity: genetics or lifestyle? It takes two to tango. Immun Ageing. 2016;13(1):12.

    Article  Google Scholar 

  33. 33.

    Dato S, Rose G, Crocco P, Monti D, Garagnani P, Franceschi C, et al. The genetics of human longevity: an intricacy of genes, environment, culture and microbiome. Mech Ageing Dev. 2017;165:147–55.

    CAS  Article  Google Scholar 

  34. 34.

    Iachine I, Skytthe A, Vaupel JW, McGue M, Koskenvuo M, Kaprio J, et al. Genetic influence on human lifespan and longevity. Hum Genet. 2006;119(3):312.

    Article  Google Scholar 

  35. 35.

    Bishop NA, Guarente L. Genetic links between diet and lifespan: shared mechanisms from yeast to humans. Nat Rev Genet. 2007;8(11):835.

    CAS  Article  Google Scholar 

  36. 36.

    Kerber RA, O’Brien E, Smith KR, Cawthon RM. Familial excess longevity in Utah genealogies. J Gerontol Ser A Biol Sci Med Sci. 2001;56(3):B130–9.

    CAS  Article  Google Scholar 

  37. 37.

    Perls T, Terry D. Understanding the determinants of exceptional longevity. Ann Intern Med. 2003;139(5_Part_2):445–9.

    Article  Google Scholar 

  38. 38.

    Greenbaum D, Colangelo C, Williams K, Gerstein M. Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biol. 2003;4(9):117.

    Article  Google Scholar 

  39. 39.

    Hu CJ, Song G, Huang W, Liu GZ, Deng CW, Zeng HP, et al. Identification of new autoantigens for primary biliary cirrhosis using human proteome microarrays. Mol Cell Proteom. 2012;11(9):669–80.

    CAS  Article  Google Scholar 

  40. 40.

    Savitski MM, Mathieson T, Zinn N, Sweetman G, Doce C, Becher I, et al. Measuring and managing ratio compression for accurate iTRAQ/TMT quantification. J Proteome Res. 2013;12(8):3586–98.

    CAS  Article  Google Scholar 

  41. 41.

    Shi Q, Chen LN, Zhang BY, Xiao K, Zhou W, Chen C, et al. Proteomics analyses for the global proteins in the brain tissues of different human prion diseases. Mol Cell Proteom. 2015;14(4):854–69.

    CAS  Article  Google Scholar 

  42. 42.

    Ferreira E, Shaw DM, Oddo S. Identification of learning-induced changes in protein networks in the hippocampi of a mouse model of Alzheimer’s disease. Transl Psychiatry. 2016;6(7):e849.

    CAS  Article  Google Scholar 

  43. 43.

    Nordestgaard BG, Chapman MJ, Ray K, Borén J, Andreotti F, Watts GF, et al. Lipoprotein (a) as a cardiovascular risk factor: current status. Eur Heart J. 2010;31(23):2844–53.

    CAS  Article  Google Scholar 

  44. 44.

    Kamstrup PR, Tybjærg-Hansen A, Nordestgaard BG. Lipoprotein (a) and risk of myocardial infarction–genetic epidemiologic evidence of causality. Scand J Clin Lab Invest. 2011;71(2):87–93.

    CAS  Article  Google Scholar 

  45. 45.

    Pati N, Rouf A, Pati U. Simultaneous mutations (A/G − 418 and C/T − 384) in the Apo (a) promoter of individuals with low Lp (a) levels. Mol Genet Metab. 2000;69(2):165–7.

    CAS  Article  Google Scholar 

  46. 46.

    Barrera GJ, Sanchez G, Gonzalez JE. Trefoil factor 3 isolated from human breast milk downregulates cytokines (IL8 and IL6) and promotes human beta defensin (hBD2 and hBD4) expression in intestinal epithelial cells HT-29. Bosn J Basic Med Sci. 2012;12(4):256.

    CAS  Article  Google Scholar 

  47. 47.

    Dong G, Wearsch PA, Peaper DR, Cresswell P, Reinisch KM. Insights into MHC class I peptide loading from the structure of the tapasin-ERp57 thiol oxidoreductase heterodimer. Immunity. 2009;30(1):21–32.

    Article  Google Scholar 

  48. 48.

    Santana-Codina N, Carretero R, Sanz-Pamplona R, Cabrera T, Guney E, Oliva B, et al. A transcriptome-proteome integrated network identifies ERp57 as a hub that mediates bone metastasis. Mol Cell Proteom. 2013;M112:022772.

    Google Scholar 

  49. 49.

    Dihazi H, Dihazi GH, Bibi A, Eltoweissy M, Mueller CA, Asif AR, et al. Correction: secretion of ERP57 is important for extracellular matrix accumulation and progression of renal fibrosis, and is an early sign of disease onset. J Cell Sci. 2018;131(13):jcs219014.

    Article  Google Scholar 

  50. 50.

    Atukeren P, Oner S, Baran O, Kemerdere R, Eren B, Cakatay U, et al. Oxidant and anti-oxidant status in common brain tumors: correlation to TP53 and human biliverdin reductase. Clin Neurol Neurosurg. 2017;158:72–6.

    Article  Google Scholar 

  51. 51.

    Columbo M, Bochner BS. Human skin mast cells adhere to vitronectin via the αvβ3 integrin receptor (CD51/CD61). J Allergy Clin Immunol. 2001;107(3):554.

    CAS  Article  Google Scholar 

  52. 52.

    Jakus Z, Fodor S, Abram CL, Lowell CA, Mócsai A. Immunoreceptor-like signaling by β2 and β3 integrins. Trends Cell Biol. 2007;17(10):493–501.

    CAS  Article  Google Scholar 

  53. 53.

    Campana D, Nori F, Piscitelli L, Morselli-Labate AM, Pezzilli R, Corinaldesi R, et al. Chromogranin A: is it a useful marker of neuroendocrine tumors? J Clin Oncol. 2007;25(15):1967–73.

    CAS  Article  Google Scholar 

  54. 54.

    Wheeler DL, Iida M, Dunn EF. The role of Src in solid tumors. Oncologist. 2009;14(7):667–78.

    CAS  Article  Google Scholar 

  55. 55.

    Menni C, Kiddle SJ, Mangino M, Vinuela A, Psatha M, Steves C, et al. Circulating proteomic signatures of chronological age. J Gerontol A Biol Sci Med Sci. 2015;70(7):809–16.

    CAS  Article  Google Scholar 

  56. 56.

    Tanaka T, Biancotto A, Moaddel R, Moore AZ, Gonzalez-Freire M, Aon MA, et al. Plasma proteomic signature of age in healthy humans. Aging Cell. 2018;17(5):e12799.

    Article  Google Scholar 

  57. 57.

    Baird AL, Westwood S, Lovestone S. Blood-based proteomic biomarkers of alzheimer’s disease pathology. Front Neurol. 2015;6:236.

    Article  Google Scholar 

  58. 58.

    Di Narzo AF, Telesco SE, Brodmerkel C, Argmann C, Peters LA. High-throughput characterization of blood serum proteomics of ibd patients with respect to aging and genetic factors. PLos Genet. 2017;13(1):e1006565.

    Article  Google Scholar 

  59. 59.

    Xiao D, Yun Z, Wu T, Zhang J, Shao J. Helicobacter pylori infection associated miR-451 suppressed the proliferation of gastric cancer by targeting CASC4. Int J Clin Exp Pathol. 2016;9(12):12390–8.

    CAS  Google Scholar 

  60. 60.

    Boutros R, Byrne JA. D53 (TPD52L1) is a cell cycle-regulated protein maximally expressed at the G2-M transition in breast cancer cells. Exp Cell Res. 2005;310(1):152–65.

    CAS  Article  Google Scholar 

  61. 61.

    Suda T, Tsunoda T, Uchida N, Watanabe T, Hasegawa S, Satoh S, et al. Identification of secernin 1 as a novel immunotherapy target for gastric cancer using the expression profiles of cDNA microarray. Cancer Sci. 2006;97(5):411–9.

    CAS  Article  Google Scholar 

  62. 62.

    Gunduz M, Ouchida M, Fukushima K, Ito S, Jitsumori Y, Nakashima T, et al. Allelic loss and reduced expression of the ING3, a candidate tumor suppressor gene at 7q31, in human head and neck cancers. Oncogene. 2002;21(28):4462–70.

    CAS  Article  Google Scholar 

  63. 63.

    Huber O, Ménard L, Haurie V, Nicou A, Taras D, Rosenbaum J. Pontin and reptin, two related ATPases with multiple roles in cancer. Cancer Res. 2008;68(17):6873–6.

    CAS  Article  Google Scholar 

  64. 64.

    Lakomek K, Stoehr G, Tosi A, Schmailzl M, Hopfner K-P. Structural basis for dodecameric assembly states and conformational plasticity of the full-length AAA + ATPases Rvb1 Rvb2. Structure. 2015;23(3):483–95.

    CAS  Article  Google Scholar 

  65. 65.

    Jiang H, Shukla A, Wang X, Chen W-Y, Bernstein BE, Roeder RG. Role for Dpy-30 in ES cell-fate specification by regulation of H3K4 methylation within bivalent domains. Cell. 2011;144(4):513–25.

    CAS  Article  Google Scholar 

  66. 66.

    Beere HM, Wolf BB, Cain K, Mosser DD, Mahboubi A, Kuwana T, et al. Heat-shock protein 70 inhibits apoptosis by preventing recruitment of procaspase-9 to the Apaf-1 apoptosome. Nat Cell Biol. 2000;2(8):469.

    CAS  Article  Google Scholar 

  67. 67.

    Hsu Y, Shi S, Wu W, Ho L, Lai J. Protective roles of interferon-induced protein with tetratricopeptide repeats 3 (IFIT3) in dengue virus infection of human lung epithelial cells. PLoS One. 2013;8(11):e79518.

    CAS  Article  Google Scholar 

Download references


Not applicable.


This study was supported by the CAMS Innovation Fund for Medical Sciences (CIFMS, No. 2017-I2M-3-021), the CAMS Project for Reform of Scientific Research System (2017-TZGG-16), and the Sichuan Provincial Science and Technology Project (Nos. 2017JY0064 and 2019YJ0280).

Author information




SY, ZW and CL conceived the study and participated in its design. JX, XZ, YL, ZS, LC and YH collected and classified the plasma samples. ZW, ML, RZ, XD and HC carried out the TMT- proteomic experiments. SY, FJ.L, PJ and FZ.L carried out the proteome microarray experiments. SY, ZW and ML analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zongkui Wang or Changqing Li.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the Institute of Blood Transfusion (No.201716). In accordance with the Declaration of Helsinki, all participants gave informed consent prior to their entering the study.

Consent for publication

All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ye, S., Ma, L., Zhang, R. et al. Plasma proteomic and autoantibody profiles reveal the proteomic characteristics involved in longevity families in Bama, China. Clin Proteom 16, 22 (2019).

Download citation


  • Bama longevity area
  • Plasma proteomics
  • Autoantibody profiles
  • Tandem mass tag
  • Human proteome microarray