Research Paper Volume 9, Issue 11 pp 2316—2333

Identification of new genetic variants of HLA-DQB1 associated with human longevity and lipid homeostasis—a cross-sectional study in a Chinese population

Fan Yang1,2, , Liang Sun1, , Xiaoquan Zhu1, , Jing Han1,2, , Yi Zeng3,4, , Chao Nie5, , Huiping Yuan1, , Xiaoling Li1,2, , Xiaohong Shi1, , Yige Yang1, , Caiyou Hu6, , Zeping Lv6, , Zezhi Huang7, , Chenguang Zheng8, , Siying Liang9, , Jin Huang10, , Gang Wan11, , Keyan Qi1, , Bin Qin1, , Suyan Cao1, , Xin Zhao1, , Yongqiang Zhang1, , Ze Yang1,2, ,

  • 1 The MOH Key Laboratory of Geriatrics, Beijing Hospital, National Center of Gerontology, Beijing, China
  • 2 Graduate School of Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China
  • 3 Center for the Study of Aging and Human Development and Geriatrics Division, Medical School of Duke University, Durham, NC 27708, USA
  • 4 Center for Healthy Aging and Development Studies, National School of Development, Peking University, Beijing, China
  • 5 BGI-Shenzhen, Shenzhen, Guangdong, China
  • 6 Jiangbing Hospital, Guangxi Zhuang Autonomous Region, Nanning, Guangxi, China
  • 7 Office of Longevity Cultural, People's Government of Yongfu County, Yongfu, Guangxi, China
  • 8 Birth Defects Prevention and Control Research Institute, Guangxi Zhuang Autonomous Region Women and Children Health Care Hospital, Nanning, Guangxi, China
  • 9 Genetic Testing Center Qingdao Women and Children’s Hospital, Qingdao University, Qingdao, China
  • 10 Department of Obstetrics and Gynecology, Aviation General Hospital of China Medical University, Beijing, China
  • 11 Capital Medical University Affiliated Beijing Ditan Hospital, Beijing, China

Received: September 1, 2017       Accepted: November 2, 2017       Published: November 10, 2017       Corrected online: March 15, 2018      

https://doi.org/10.18632/aging.101323
How to Cite
This article has been corrected. See Correction. Aging (Albany NY). 2018; 10:515-515 . https://doi.org/10.18632/aging.101400  PMID: 29615541

Copyright: © 2017 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Healthy longevity has been an unremitting pursuit of human, but its genetic and the environment causes are still unclear. As longevity population is a good healthy aging model for understanding how the body begin aging and the process of aging, and plasma lipids metabolism and balance is a very important to life maintain and physiologic functional turnover. It is important to explore how the effect of genetic variants associated long-life individuals on lipids metabolism and balance. Therefore, we developed a comparative study based population which contains 2816 longevity and 2819 control. Through whole-exome sequencing and sanger sequencing genotypes, we identified four new single nucleotide polymorphisms of HLA-DQB1(major histocompatibility complex, class II, DQ beta 1), rs41542812 rs1049107 rs1049100 rs3891176(Prange=0.048-2.811×10-8 for allele frequencies), associated with longevity in Chinese Longevity Cohort. Further, by analysis of the longevity-variants linked to blood lipids, we identified HLA-DQB1 rs1049107, T-carriers (PHDL=0.006, OR: 11.277; PTG=9.095×10-7, OR: 0.025; PLDL/HDL=0.047, OR: 1.901) and HLA-DQB1 rs1049100, T-carriers (PTG=1.799×10-6, OR: 0.028) associated with lipid homeostasis in long lived individuals. Our finding showed that longevity and lipid homeostasis were associated with HLA-DQB1 and suggested that immune gene variants could act on both new function of maintaining the homeostasis and anti-aging in longevity.

Introduction

For thousands of years, humans have pursued increased longevity. However, increasing the longevity of the human population has been unlikely. In recent years, with the increase in social economic levels and an ageing population, the likelihood of an increase in human longevity is gradually increasing and is becoming the focus of people’s attention. As we know, human longevity is determined by both genetic and environmental factors, but its genetics and the environment causes, particularly their interaction mechanism, are still not clear. Therefore, it is necessary to develop the study and obtain supportive data for us to be able to provide inspiration and direction for future research in the field of longevity and anti-aging.

A long-lived population provides a good healthy ageing model for understanding how the body begins to age and the process of ageing. Several previous studies, both cross-sectional and follow-up, investigated long-lived subjects who usually manage to delay major age-related diseases, such as metabolic-related diseases (cardiovascular disease, diabetes, and neurodegeneration disorders); in general, these age-related diseases had high morbidity and mortality, but most centenarians escape these diseases [13] and successfully exhibit a healthy ageing state.

The common view is that human longevity is determined by genetic factors. Several Genome-wide association studies (GWAS) in European, North American and Chinese individuals identified some loci associated with longevity, i.e., TOMM40/APOE/APOC1 loci, 5q33.3, IL6 and ANKRD20A9P in recent years [46]. With next-generation sequencing (NGS), we have an efficient way to discover the potential causal loci in the search for sequence variants. NGS technology has facilitated the identification of variants and shown that HLA-DQB1*05 and HLA-DQB1*03 were associated with longevity in Japanese individuals [7,8]. Previous studies of longevity genetics found that with increasing age, the contribution of heredity increased [9].

Along with ageing or senescence, the internal environment function of each hierarchy and each system is low and disordered, but there is still a lack of explicit evidence to define the genetically encoded program of ageing that functions to maintain the homeostasis of the human body, i.e., the plasma lipoproteins, in vivo.

The metabolism and balance of plasma lipids are very important to life and physiologic functional turnover. They are also the representative biomarker for cholesterol metabolism and lipid-related diseases such as cardiovascular and cerebrovascular diseases. A number of studies have shown that the plasma lipid levels were controlled by genetic factors and identified some associated variants of lipid levels, both common and rare [1012]; for example, HLA-DQB1 variations were associated with plasma lipid balance in coronary heart disease risk [13], but there are no reports on genetic variants and plasma lipids associated with healthy longevity.

Therefore, we postulated that longer human survivors, who live to a median age of over 90 years, have a higher number of longevity genetic variants; briefly, ageing and longevity genetic variants serve as important biomarkers to represent the normal plasma lipid homeostasis in the internal environment of extremely elderly individuals. Therefore, based on our longevity cohort study in southern China, we developed a comparative study using a case-control design to detect the relationship between longevity-associated variants and plasma lipid levels in our long-lived individuals in China. Our study aims to suggest that a healthy lipid balance linked to genetic variants in long-lived individuals could play an important part in human longevity, and our data will provide new knowledge to help us understand human longevity and the possible potential applications for anti-aging primary care in social communities of the elderly.

Results

Baseline information in longevity and control subjects

There were significant differences in HDL-c (P=0.003), LDL-c (P=0.010), TG (P=3.052×10-4) and TC (P=0.012) between centenarians and controls. We detected that HDL-c (P=1.461×10-7), LDL-c (P=2.128×10-4) and TC (P=8.181×10-6) were significantly different between nonagenarians and controls. A comparison of longevity (centenarians and nonagenarians) in samples and controls, revealed that significant differences between HDL-c (P=3.438×10-9), LDL-c (P=1.269×10-5), TG (P=0.035) and TC (P=6.522×10-7) had been found. (Supplementary Table 1)

Identifying the new longevity-associated genetic variations

Whole-exome sequence screening in 100 longevity subjects

By whole exome sequence in 100 longevity subjects involving 74 nonagenarians and 26 centenarians, we compared both genomic differences between long-lived individuals and our genotype database of general Chinese (age <60 years, n=1000), and we primarily identified 2171 possible longevity-associated variants.

Bioinformatic and arrangement analysis

Aligning with GenBank, the analysis of the collection and arrangement of bioinformation based on the LongevityMap database [14] reported possible longevity variants, and we selected 17 genes with 26 variants as our longevity-associated candidate genes (SSupplementary Table 2). A flow chart of the consecutive analysis steps is depicted in Figure 1.

A flow chart of the consecutive analysis steps.

Figure 1. A flow chart of the consecutive analysis steps.

Identified new longevity-associated variants in HLA-DQB1

Through an association study based on genotyping longevity and a local younger population, we identified four genetic variants in HLA-DQB1 as longevity-associated gene variants. Four variants in HLA-DQB1 showed significant difference in allele and genotype frequencies between longevity and controls. These variants include rs41542812, rs1049107, rs1049100 and rs3891176, and this result was detected in a sample set what contained 518 longevities and 277 controls, rs41542812 (Pallele=0.048, Pgenotype=0.003); rs1049107 (Pallele=0.002, Pgenotype=0.007); rs1049100 (Pallele=4.731×10-4, Pgenotype=0.002); and rs3891176 (Pallele=2.811×10-8, Pgenotype=4.033×10-8 (Table 1).

Table 1. Association of variants and longevity.

Major alleleMinor alleleMajor homoHetroMinor home
Case/ControlCase/ControlPOR95%CICase/ControlCase/ControlCase/ControlP
Longevity vs. Controlsrs41542812114/51786/4990.0481.419(1.001-2.011)364/22956/5129/50.003
rs1049107199/86653/4400.0021.5591.178-2.064278/19097/6051/130.007
rs1049100205/85647/4414.731×10-41.6441.242-2.176275/19097/6154/120.002
rs3891176248/84600/4402.811×10-82.1651.642-2.854253/18694/6877/84.033×10-8
Nonagenarians vs. Controlsrs4154281248/51262/4990.0061.7931.176 -2.732123/22916/5116/56.652×10-5
rs104910785/86203/4401.041×10-52.1421.521-3.01890/19060/6031/135.249×10-7
rs104910089/85199/4419.297×10-72.3201.650-3.26488/19023/6133/129.258×10-8
rs3891176151/84141/4403.832×10-275.6104.044-7.78263/18625/6858/86.193×10-21
Centenarians vs. Controlsrs4154281266/51524/4990.2871.232838-1.812242/22940/5113/50.079
rs1049107114/86450/4400.0011.2960.951-1.766188/19074/6020/130.317
rs1049100116/85448/4410.0611.3430.986-1.830187/19074/6121/120.215
rs3891176107/84449/4400.1661.2480.911-1.710190/18669/6819/80.131

Identified longevity-associated variation with ageing

We detected that rs41542812 (Pallele=0.006, Pgenotype=6.652×10-5), rs1049107 (Pallele=1.041×10-5, Pgenotype=5.249×10-7), rs1049100 (Pallele=9.297×10-7, Pgenotype=9.258×10-8) and rs3891176 (Pallele=3.832×10-27, Pgenotype=6.193×10-21) were associated with nonagenarians. Thus, rs1049107 (Pallele=0.001, Pgenotype=0.317) was associated with centenarians (Table 1; Figure 2).

Association analysis identified HLA-DQB1 as longevity-associated gene. (A) Manhattan map of whole-exome sequencing; (B) Four variants in HLA-DQB1; (C) Linkage disequilibrium analysis of four variants. C1 LD map of centenarians; C1 LD map of nonagenarians; C3: LD map of longevity; C4: LD map of controls.

Figure 2. Association analysis identified HLA-DQB1 as longevity-associated gene. (A) Manhattan map of whole-exome sequencing; (B) Four variants in HLA-DQB1; (C) Linkage disequilibrium analysis of four variants. C1 LD map of centenarians; C1 LD map of nonagenarians; C3: LD map of longevity; C4: LD map of controls.

Identified longevity-associated haplotypes

The link disequilibrium analysis detected that there was a block, formed by rs41542812, rs1049107 and rs1049100, on HLA-DQB1 (Figure 2). There were three haplotypes with frequencies >0.03. Compared with the CCC haplotype, the GTT and CTT haplotypes could increase the chance of longevity (P=3.996×109, OR: 4.367, 95%CI: 2.608-7.313; P=1.812×105, OR: 3.677, 95%CI: 1.970-6.865). Comparing between centenarians and controls, the GTT and CTT haplotypes were associated with centenarians (P=3.327×108, OR:6.484, 95%CI:3.128-13.440; and P=0.025, OR:3.373, 95%CI: 1.234-9.216, respectively). In nonagenarians, the GTT and CTT haplotypes would increase the chance to reach a nonagenarian age (P=1.921×106, OR:3.721, 95%CI: 2.116-6.543; and P=3.604×105, OR:3.770, 95%CI: 1.945-7.308, respectively; Supplementary Table 3).

Relationship of longevity-associated variants and plasma lipid homeostasis

For rs1049107, Supplementary Table 4 shows the proportion of participants with different plasma lipid levels, normal or abnormal, and LDL-c/HDL-c ratios ≤2 and those >2, according to the polymorphism genotypes. The rs1049107 T allele was associated with a lower LDL-c/HDL-c ratio (P=0.047, OR: 1.901) in the longevity group.

In HDL-c subgroups of longevity samples, the rs1049107 TT genotype samples and T allele carriers showed a better HDL-c level than did the rs1049107 CC genotype samples (P=0.007, OR: 10.636; P=0.006, OR: 11.277). The normal HDL-c subgroup contained a higher rs1049107 T allele proportion (P=1.458×10-4, OR: 10.689), but this subgroup also had a lower rs3891176 C allele ratio than the major allele (P=0.013, OR: 0.360).

The combined rs1049107 genotype (TT+CT) in the control group, assuming a dominant effect of the T allele, was associated with a higher LDL-c level (P=1.243×10-7, OR: 0.098). In longevity group rs1049107 CT and rs1049100 CT, the genotype was associated with a higher TG level (P=9.095×10-7, OR: 0.025; P=1.799×10-6, 0.028), but it was opposite in the controls (P=0.004, OR: 11.667; P=0.004, OR: 11.667) compared with major homozygotes.

Overall, the dominant effect of the rs1049107 T allele was associated with a normal HDL-c level (P=0.001, OR: 72.757,95%CI: 4.378-984.360) and a higher TG level (P=0.035, OR: 0.312, 95%CI: 0.105-0.924) in the longevity group and with a higher LDL-c level in controls (P=0.050, OR: 0.383, 95%CI: 0.147-0.999) by a logistic regression analysis.

For the rs3891176 CA genotype, the dominant effect of the variant C allele and rs3891176 C allele were associated with a lower LDL-c/HDL-c ratio (Supplementary Table 3) by a logistic regression analysis, and the dominant effect of the rs3891176 C allele was associated with a lower LDL-c/HDL-c ratio (P=0.034, OR: 2.402, 95%CI: 1.068-5.402). Furthermore, in the longevity group, we detected that major homozygotes of rs3891176 had significantly higher HDL-c levels compared to minor allele carriers (P=0.044) (Table 2).

Table 2. Allele and genotype distribution analysis in different plasma lipids level group.

LongevityControl
HDL>1.04HDL<=1.04POR95%CIHDL>1.04HDL<=1.04POR95%CI
rs41542812CC5912ref9829ref
GC310.5410.6100.058-6.3771480.1750.5180.198-1.356
GG2010.2854.0680.497-33.288301.0001.2120.13-11.262
GG+GC2320.2782.3390.485-11.2721780.3290.6290.246-1.605
CC+GC62130.2910.2380.029-1.939112371.0000.7430.081-6.858
C12125ref21066ref
G4330.0762.9610.851-10.3062080.5840.7860.331-1.867
rs1049107CC4413ref8024ref
CT200.9540.9330.090-9.70225110.3720.6820.293-1.584
TT3610.00710.6361.327-85.232401.0001.5430.172-13.835
TT+CT3810.00611.2271.403-89.84629110.5790.7910.345-1.815
CC+CT46130.0090.0980.012-0.787105351.0000.5890.067-5.210
C9026ref18559ref
T7421.458×10-410.6892.456-46.52333110.9070.9570.455-2.011
rs1049100CC449ref8024ref
CT110.3330.2050.012-3.58325110.3720.6820.293-1.584
TT3740.3141.8920.539-6.646401.0001.5430.172-13.835
TT+CT3850.4601.5550.479-5.0429110.5790.7910.345-1.815
CC+CT45100.2470.4860.141-1.678105351.0000.5890.067-5.210
C8919ref18559ref
T7590.1801.7790.760-4.16433110.9070.9570.455-2.011
rs3891176AA577ref7325ref
CA1111.0001.3510.151-12.12870.5131.3700.533-3.523
CC1460.0700.2870.083-0.987700.4552.8490.34-23.893
CC+CA2570.2190.4390.139-1.3833570.2531.7120.676-4.34
AA+CA6880.0683.6431.092-12.152101320.6860.3860.047-3.205
A12515ref17457ref
C39130.0130.3600.158-0.8224270.1161.9660.837-4.618
LDL<3.12LDL>=3.12POR95%CILDL<3.12LDL>=3.12POR95%CI
rs41542812CC4427ref10522ref
GC400.4063.1110.345-28.0291931.0001.3270.361-4.876
GG1560.4271.5340.531-4.433301.0000.8680.093-8.13
GG+GC1960.2041.9430.690-5.4722230.5121.5370.423-5.587
CC+GC48270.5260.7110.247-2.04812425
C9254ref22947ref
G34120.1751.6630.794-3.4822530.5921.7100.496-5.898
rs1049107CC3423ref8915ref
CT201.0002.0570.202-20.9752880.2760.5900.226-1.537
TT27100.1861.8260.744-4.483310.4790.5060.049-5.188
TT+CT29100.1361.9620.804-4.7893190.2440.5810.231-1.46
CC+CT36230.2300.5800.237-1.41811723
C7046ref20638ref
T56200.0571.8400.978-3.46134100.6261.2190.55-2.703
rs1049100CC3221ref8915ref
CT201.0002.0000.195-20.4862880.2760.5900.226-1.537
TT29120.2971.5860.665-3.783310.4790.5060.049-5.188
TT+CT31120.2291.6950.714-4.0243190.2440.5810.231-1.46
CC+CT34210.3630.6700.282-1.59111723
C6642ref20638ref
T60240.1351.5910.863-2.93234100.6261.2190.55-2.703
rs3891176AA4321ref7919ref
CA750.7410.6840.194-2.4123230.1392.5650.71-9.273
CC1370.8560.9070.315-2.609430.1580.3210.066-1.555
CC+CA20120.6490.8140.336-1.9743660.4701.4430.531-3.918
AA+CA50261.0001.0360.368-2.91211122
A9347ref19041ref
C33190.7000.8780.452-1.7064090.9180.9590.432-2.13
TG<1.7TG>=1.7POR95%CITG<1.7TG>=1.7POR95%CI
rs41542812CC5813ref8938ref
GC310.7380.6720.065-6.9931390.3060.6170.243-1.565
GG1740.9390.9530.275-3.306301.0001.7330.188-16.012
GG+GC2051.0000.8970.284-2.8311690.5310.7510.305-1.847
CC+GC61140.2870.2560.031-2.09010247
C11927ref19185ref
G3790.8710.9330.403-2.161990.8830.9390.408-2.162
rs1049107CC507ref7034ref
CT2119.095×10-70.0250.005-0.1402300.00411.6671.516-89.79
TT2600.1534.2350.503-35.6594130.0010.1490.045-0.493
TT+CT28110.0500.3560.124-1.02327130.9821.0090.463-2.197
CC+CT52180.0100.1030.013-0.8149334
C10225ref16368ref
T54110.6431.2030.550-2.6331260.0200.4970.275-0.9
rs1049100CC467ref7034ref
CT2111.799×10-60.0280.005-0.1522300.00411.6671.516-89.79
TT3000.1465.2770.628-44.3024130.0010.1490.045-0.493
TT+CT32110.1220.4430.155-1.26427130.9821.0090.463-2.197
CC+CT48180.0030.0830.011-0.6539334
C9425ref16368ref
T62110.3061.4990.688-3.26431260.0200.4970.275-0.9
rs3891176AA4915ref6533ref
CA1110.4423.3670.401-28.2582690.3851.4670.617-3.487
CC1820.3382.7550.573-13.259521.0001.2690.234-6.896
CC+CA2930.0962.9590.789-11.09931110.3821.4310.639-3.201
AA+CA60160.3470.4170.087-1.9869142
A10931ref15675ref
C4750.0482.6730.979-7.336130.4161.3310.667-2.658
TC<5.2TC>=5.2TC<5.2TC>=5.2
rs41542812CC4526ref11314ref
GC400.6612.4110.268-21.6611930.7190.7850.206-2.992
GG1560.4961.4440.499-4.181300.4760.5260.055-5.026
GG+GC1960.2501.8300.649-5.1612231.0000.9090.241-3.428
CC+GC49260.6000.7540.261-2.17513217
C9452ref24531ref
G34120.2321.5670.748-3.2862531.0001.0540.301-3.697
rs1049107CC3621ref9410ref
CT201.0001.7840.175-18.2223060.3600.5320.178-1.586
TT26110.4771.3790.568-3.347400.5020.5790.062-5.416
TT+CT28110.3781.4850.615-3.5833460.3810.6030.204-1.785
CC+CT38210.5530.7660.316-1.853124160.5471.4710.162-13.352
C7442ref21826ref
T54220.29735.0000.747-2.63860.6020.7550.291-1.957
rs1049100CC3419ref9410ref
CT201.0001.7140.167-17.63060.3600.5320.178-1.586
TT28130.6741.2040.507-2.858400.5020.5790.062-5.416
CC+CT36190.7700.8800.372-2.081124160.5471.4710.162-13.352
TT+CT30130.5621.2900.546-3.0463460.3810.6030.204-1.785
C7038ref21826ref
T58260.5371.2110.659-2.2253860.6020.7550.291-1.957
rs3891176AA4123ref8612ref
CA841.0001.1220.304-4.1353320.3542.3020.489-10.845
CC1550.3651.6830.542-5.229611.0000.8370.093-7.567
CC+CA2390.4441.4340.569-3.6133930.5531.8140.484-6.794
AA+CA49270.5520.7260.252-2.089119140.5561.4170.159-12.636
A9050ref20526ref
C38140.2511.5080.746-3.0474540.5251.4270.474-4.291
LDL/HDL≤2LDL/HDL>2POR95%CILDL/HDLPOR95%CI
rs41542812CC4724ref7156ref
GC311.0001.5320.151-15.52612100.9060.9460.381-2.35
GG1470.9681.0210.364-2.866210.7101.5770.139-17.845
GG+CC1780.8691.0850.410-2.87314110.9931.0040.423-2.381
CC+GC50251.0001.0000.358-2.79183660.7050.6290.056-7.086
C9749ref154122ref
G31150.9051.0440.515-2.11416120.8911.0560.482-2.316
rs1049107CC3522ref6143ref
CT111.0000.6290.037-10.57316200.1400.5640.263-1.211
TT2890.1501.9560.778-4.912221.0000.7050.096-5.201
TT+CT29100.1861.8230.745-4.46118220.1400.5770.277-1.203
CC+CT36230.2080.5030.201-1.25777531.0001.4530.198-10.638
C7145ref138106ref
T57190.0471.9011.003-3.60420240.1730.6400.336-1.22
rs1049100CC3419ref6143ref
CT111.0000.5590.033-9.45116200.1400.5640.263-1.211
TT29120.5011.3500.562-3.244221.0000.7050.096-5.201
TT+CT30130.5621.2900.546-3.04618220.1400.5770.277-1.203
CC+CT35200.4660.7240.304-1.72677531.0001.4530.198-10.638
C6939ref138106ref
T59250.3551.3340.724-2.45720240.1730.6400.336-1.22
rs3891176AA4420ref4949ref
CA660.2090.4550.130-1.5852690.0132.8891.228-6.794
CC1460.9161.0610.356-3.163430.7151.3330.283-6.272
CC+CA20120.5400.7580.311-1.84430120.0192.5001.149-5.442
AA+CA50260.7950.8240.283-2.396756810.8270.179-3.830
A9446ref124107ref
C34180.8180.9240.472-1.80934150.0441.9561.011-3.785

Analysis of plasma-lipid homeostasis and longevity-associated haplotypes

To analyse the effect of haplotypes on plasma lipids, we divided samples between plasma lipid levels and age. We found that rs41542812, rs1049107 and rs1049100 would form a block in HLA-DQB1 and that the CCC haplotype, which was formed by major alleles, seemed to be a reference. In the centenarian group, we obtained none of the significant haplotypes for four plasma lipid indices. In nonagenarians, the block contained three haplotypes, and the CTT haplotype was associated with a higher LDL-c/HDL-c ratio (P=0.037, OR: 0.301, 95%CI: 0.093-0.969). In the longevity group, the GTT and CTT haplotypes were associated with higher HDL-c levels (P=0.027, OR:4.750, 95%CI: 1.057-21.337; P=0.045, OR:6.667, 95%CI: 0.857-51.883). The CTT haplotype was associated with a higher TG level (P=1.096×10-4=, OR: 0.178, 95%CI: 0.070-0.451), and the CTT haplotype was associated with LDL-c/HDL-c <=2 (P=0.027, OR: 3.403, 95%CI: 1.098-10.544; Supplementary Table 5).

Interaction between longevity-associated variants and internal milieu lipid levels

An interaction analysis between longevity associated variants and internal milieu lipid levels on age detected significant results in a corrected model (R2=0.375, R2Adjusted=0.332). Significant values of rs14917 & rs1491399 (P=0.001), rs1491 & LDL-c/HDL-c (P=0.016), rs927299 & LDL-c/HDL-c (P=0.028), and rs1491 & rs927299 & LDL-c/HDL-c (P=0.031) were detected in this model.

Discussion

Longevity-associated variants and haplotypes identified

Because of the low prevalence of delayed age-related disease, the longevity population was an ideal population to be considered a “real negative control” in human successful healthy ageing and age-related disease. Genetic research about longevity had identify some longevity-associated factors, such as FOXO3A [15,16], SIRT1 [17], APOE [18,19], IL6 and ANKRD20A9P [20]. There had been multiple genetic studies of Chinese longevity, and FOXO3 [21], IGFBP-3 [22], CETP [23], and SIRT1 [24]. Studies in Sardinian centenarians showed the association between HLA-DQB1 and longevity [25,26]. The association between longevity and HLA-DQB1*05 or HLA-DQB1*03 was identified in a Japanese population [7,8]. However, there was no association study about SNPs in HLA-DQB1 and longevity. In this study, we identified four new SNPs in HLA-DQB1, rs41542812 (DQB1*05), rs1049107 (DQB1*03), rs1049100 (DQB1*03) and rs3891176 (DQB1*02), associated with longevity in our cohort (Table 1). Further, we identified that rs41542812, rs1049107, and rs1049100 were closely linked to disequilibrium in one block; that is, the haplotypes GTT and CTT, which significantly increased the chance of longevity (P=3.996×109, OR: 4.367, 95%CI: 2.608-7.313; P=1.812×105, OR: 3.677, 95%CI: 1.970-6.865) to our knowledge have also never been reported before.

In these new longevity-associated variants, we noted through the use of stratified analysis that only rs1049107 was a centenarian-associated variant in long-lived individuals, but the other three variants were associated with both centenarians and nonagenarians (Table 1). This result suggested that the longevity-associated genetic variance seems to increase with ageing, so we observed the variance only in centenarians but not nonagenarians.

Longevity variants associated with lipid homeostasis in LLIs

The main cholesterol-carrying lipoproteins are LDL-c and HDL-c. LDL-c is always considered the bad lipoprotein and can increase CVD and metabolic syndrome risk. In contrast, HDL-c is considered the good lipoprotein and can decrease metabolic-related disease and CVD risk. The importance of the balance of plasma lipids for maintaining health is obvious. The LDL-c/HDL-c ratio is increasingly being considered a valuable tool in CVD-related clinical research. Many variants have been identified as potentially affecting human lipids, as have some longevity-associated variants. The changes in lipid levels in our cohort were consistent with those in previous studies [2729]. However, evidence for the balance of plasma lipids being controlled by longevity-associated loci is still lacking.

Furthermore, we identified HLA-DQB1 rs1049107, T-carriers (PHDL=0.006, OR: 11.277; PTG=9.095×10-7, OR: 0.025; PLDL/HDL=0.047, OR: 1.901) and HLA-DQB1 rs1049100, T-carriers (PTG=1.799×10-6, OR: 0.028) associated with lipid homeostasis in LLIs.

The interaction analysis of longevity variants and lipid levels

The interaction analysis between lipids and genotype results show that HLA-DQB1 longevity variant alleles or haplotype carriers have a trend of better blood lipid levels, thus benefiting healthy longevity in LLIs (Table 3; Supplementary Figure 1). The interaction analysis results suggest that with longevity-associated variants enriched in the longevity population, the plasma lipid levels were better controlled than by adjusting age. They also suggested that it is possible for longevity variants to possess functions to regulate the homeostasis of plasma lipids in the internal environment of LLIs.

Table 3. The effect of variants and lipid levels interaction on long-life (years) in all subjects.

SourceType III Sum of SquaresdfMean SquareFP
Corrected Model104143.199323254.4758.738<0.001
8045.24318045.24321.602<0.001
rs14917*rs14913993.12413993.12410.5470.001
rs1491*LDL/HDL2164.92212164.9225.8130.016
rs927299*LDL/HDL2692.90821346.4543.6150.028
rs1491*rs927299*LDL/HDL6703.2108837.9012.1360.031
Error173924.649467372.430
Total2559434.000500
Corrected Total278067.848499

Functional analysis of the new longevity-associated variants in HLA-DQB1

HLA-DQB1 contained 5 exons. Exon 1 encodes the leader peptide, exons 2 and 3 encode the two extracellular protein domains, exon 4 encodes the transmembrane domain, and exon 5 encodes the cytoplasmic tail. All four novel variants in HLA-DQB1 are functional missense mutations. I These four variants, rs41542812 (Gln158His), rs1049107 (Gly157Ser) and rs1049100 (Val148Ile), lays in β chain of HLA-DQB1. Variety in β chain may change the binding site for special antigen, which plays a key role in maintaining internal homeostasis in vivo, such as lipid homeostasis. In addition, rs3891176 (Ala29Ser), located in exon 1, may change the function of the leader peptide and transform the signal that HLA-DQB1 received (Figure 3). Therefore, we first identified the association between longevity and these four variants. oxLDL (oxidized low-density lipoprotein) or LPS bands on HLA-DQB1 of T cell, and stimulate T cell producing cytokines. Cytokines effect on liver cell, and increase or decrease fatty acid metabolism. OxLDL could also bands on LDLR (low-density lipoprotein receptor) or HLA-DQB1 of liver cell. The signal could also change fatty acid metabolism and HDL-c producing of liver cell. With this possible process, plasma lipid homeostasis could affect by HLA-DQB1 variants (Figure 3).

Based on the resource from relevant information on HLA-DQB1 and homeostasis/metabolism phenotype, we draw the possible mechanic or interactive path way.

Figure 3. Based on the resource from relevant information on HLA-DQB1 and homeostasis/metabolism phenotype, we draw the possible mechanic or interactive path way.

We speculated that the HLA-DQB1 variants decreased the expression and/or antigen binding function of the HLA Class II antigen DQ protein β chain on the Antigen Presenting Cell (i.e., the macrophage), reduces the cytokines released by T-lymphocytes, thus reducing liver cell synthesis and release of cholesterol, to maintain the balance of cholesterol metabolism in vivo. This theory remains to be studied more in depth to obtain the corresponding evidence.

Findings in this research suggest a possible gene effecting model. HLA-DQB1 affects the longevity through a certain biological pathway, and these pathways could influence the balance of plasma lipids through unknown biological process. Longevity and the balance of plasma lipids could affect each other. A longevity gene could help keep plasma lipids balanced, and a better balance of plasma lipids could promote homeostasis and healthy longevity. However, the mechanisms of longevity and plasma lipid levels are still not clear. Furthermore, a gene associated with longevity and plasma lipid levels could be identified in multiple large longevity populations, and rare or low-frequency variants must receive more attention. The mechanism could be studied in cell or animal models by using new technology, such as Crispr/Cas 9. It would be useful for us to understand the mechanism of longevity and help the wider community achieve healthy ageing.

Methods

Subjects

The current study was based on the Longevity and Health of Aging Population (LHAP) study conducted in Bama County, Guangxi, China, in 2008 and the Chinese Longitudinal Healthy Longevity Survey (CLHLS) during 1997-2015 (Supplementary Table 6).

To investigate the association between genotype and plasma lipid levels, there were 100 long-lived individuals (aged 90-107 years, 32 males, 68 females), and 172 controls (aged 40-65 years, 75 males, 97 females). The longevity group contained 26 centenarians (5 males, 21 females) and 74 nonagenarians (27 males, 47 females). All control subjects lacked longevity history (no lineal family members aged above 85 for three generations). The Ethics Committee of Beijing Hospital, Ministry of Health, approved the study protocol. All participants were informed and provided informed consent in writing. All clinical investigations were conducted according to the principles of the Declarations of Helsinki. Laboratory parameters were recorded, including total cholesterol (TC), triglyceride (TG), high density lipoprotein-cholesterol (HDL-c) and low-density lipoprotein-cholesterol (LDL-c).

Exome target sequencing

We sequenced the segment of exomes using a targeted NGS approach to analyse SNVs (single nucleotide variations) and small insertions and/or deletions (Indels). The design parameters were as follows: 1) the bait length was 120 bp; 2) the bait tiling frequency was 2×; and 3) the human reference genome (GRch37) was used to search the baits. Genomic DNA (3 μg) was used to construct DNA libraries that contained index sequences for identifying DNA samples. The targeted genomic segment was captured with the Agilent Capture kit 39M (Agilent Technologies, Santa Clara, CA, USA) by NGS on the Illumina HiSeq 2500 platform to identify mutations in 100 Bama longevity subjects using 100-bp paired-end reads. Base sequences were aligned to the reference genome (GRch37), and variations SNVs and Indels) were called using the Samtools software and annotated by comparing them with the 1000 Genomes database. The average coverage of each base in the targeted region was 50×, which resulted in ≥90% of the targeted bases being sufficiently covered for variant calling (≥10×). The average coverage was ≥95% at the NGS stage. We compared the data from our 100 LLIs to publicly available genotype data from a CHS (Chinese Han South) population (n = 97) from the 1000 Genome project phase I. We found 115327 SNVs, in total. Of the 115327 SNVs, 84914 were consistent with those in the public SNP database, and 30413 were previously unknown.

Selection of variations

We selected candidate variation from the target sequence dataset using four continuous steps.

1) To minimize error rates and the impact of random sampling variance, at least three reads had to be obtained from both the forward and reverse DNA strands (double-stranded coverage) for all variants that were kept for association tests.

2) Variants located at repeat sequences, including STR (short tandem repeat) or single nucleotide repeat expansions, may be falsely called; thus, we first excluded these possible artificial variants.

3) The Longevity Gene database was downloaded, and all significant genes were screened. All variants of these significant genes were called and compared with the CHS population in the 1000 Genome project to obtain candidate variants with MAF (minor allele frequency) >10% in the 100 Bama longevity subjects compared to those in one public database. Seventeen genes that included 26 variants were selected as candidate longevity-associated genes. Through genotyping based on longevity and local younger populations, four genetic variants in HLA-DQB1 were identified as longevity-associated gene variants (Supplementary Table 1).

Genotyping and quality control for genotyping

The Sanger Sequencing method was used in the genotyping case-control study. To perform genotyping quality control, all associated allele carriers, and 10% of cases and 10% of controls who carried non-risk alleles, were re-genotyped by Sanger sequencing. Sequencing primers are shown in Supplementary Table 7.

Statistical analysis

Genotypes were evaluated for departure from the Hardy–Weinberg equilibrium (HWE) in cases and controls using the chi-squared goodness test. Variants with P < 0.05 were considered to deviate from the HWE. The MAFs of the variants were used as high-risk allele frequencies, and the defined type I error rate was 0.05. The genotype frequencies of the CHS population from the 1000 Genome database were used as the references for selecting candidate SNVs. A t-test was used to analyse the mean difference between the groups. A chi-squared test was used to establish the differences in the distribution of genotype and allele frequencies between the cases and controls. The odds ratio (OR) was used to estimate the strength of association between variables, with the OR and 95% confidence intervals (95%CI). The ORs and 95% CIs were calculated using the SPSS V17.0 software (SPSS; SAS Institute, Cary, NC, USA). A two-sided P value < 0.05 was considered statistically significant.

Conclusions

In this research, we identified four longevity-associated variants in HAL-DQB1, and these variants were associated with different kinds of plasma lipids. Our finding shows that longevity and plasma lipid levels were affected by multiple genes and that internal lipid homeostasis promoted healthy longevity. Further similar studies should be developed in human populations, especially investigating human lifespan and healthy longevity.

Supplementary Materials

Supplementary File

Author Contributions

F.Y., XQ. Z, and Ze. Y. designed the study. F. Y. contributed to sample management, exome sequencing and performed the genotyping. J. H., XL. L, L. S, YG. Y, XH. S. YG. Y performed bioinformatics and the statistical analysis. Yi. Z and C. N. contributed CLHLS database. CY. H., ZP. L, ZZ. H, CG. Z, B. Q, SY. C, X. Z, and YQ Z contributed Bama subjects. All authors revised and gave final approval of this manuscript.

Conflicts of Interest

The authors declare no competing financial interests. All financial and material support for this research has no potential conflicts.

Funding

This work was funded by the Natural Science Foundation of Guangxi and China (81400790, 81061120527, 81370445, 81472408, 81571385, 2014GXNSFDA118028, GKZ0991198), funding from National Department Public Benefit Research Foundation by Ministry of Health P. R. China (201302008), and 12th 5 Year National Program from the Chinese Ministry of Scientific Technology (2012BAI10B01, 2015BAI06B03).

References

  • 1. Franceschi C, Bonafè M. Centenarians as a model for healthy aging. Biochem Soc Trans. 2003; 31:457–61. https://doi.org/10.1042/bst0310457 [PubMed]
  • 2. Hitt R, Young-Xu Y, Silver M, Perls T. Centenarians: the older you get, the healthier you have been. Lancet. 1999; 354:652. https://doi.org/10.1016/S0140-6736(99)01987-X [PubMed]
  • 3. Terry DF, Wilcox M, McCormick MA, Lawler E, Perls TT. Cardiovascular advantages among the offspring of centenarians. J Gerontol A Biol Sci Med Sci. 2003; 58:M425–31. https://doi.org/10.1093/gerona/58.5.M425 [PubMed]
  • 4. Sebastiani P, Solovieff N, Dewan AT, Walsh KM, Puca A, Hartley SW, Melista E, Andersen S, Dworkis DA, Wilk JB, Myers RH, Steinberg MH, Montano M, et al. Genetic signatures of exceptional longevity in humans. PLoS One. 2012; 7:e29848. https://doi.org/10.1371/journal.pone.0029848 [PubMed]
  • 5. Deelen J, Beekman M, Uh HW, Broer L, Ayers KL, Tan Q, Kamatani Y, Bennet AM, Tamm R, Trompet S, Guðbjartsson DF, Flachsbart F, Rose G, et al. Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age. Hum Mol Genet. 2014; 23:4420–32. https://doi.org/10.1093/hmg/ddu139 [PubMed]
  • 6. Newman AB, Walter S, Lunetta KL, Garcia ME, Slagboom PE, Christensen K, Arnold AM, Aspelund T, Aulchenko YS, Benjamin EJ, Christiansen L, D’Agostino RBSr, Fitzpatrick AL, et al. A meta-analysis of four genome-wide association studies of survival to age 90 years or older: the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium. J Gerontol A Biol Sci Med Sci. 2010; 65:478–87. https://doi.org/10.1093/gerona/glq028 [PubMed]
  • 7. Akisaka M, Suzuki M, Inoko H. Molecular genetic studies on DNA polymorphism of the HLA class II genes associated with human longevity. Tissue Antigens. 1997; 50:489–93. https://doi.org/10.1111/j.1399-0039.1997.tb02904.x [PubMed]
  • 8. Akisaka M, Suzuki M. [Okinawa Longevity Study. Molecular genetic analysis of HLA genes in the very old]. Nippon Ronen Igakkai Zasshi. 1998; 35:294–98. https://doi.org/10.3143/geriatrics.35.294 [PubMed]
  • 9. Murabito JM, Yuan R, Lunetta KL. The search for longevity and healthy aging genes: insights from epidemiological studies and samples of long-lived individuals. J Gerontol A Biol Sci Med Sci. 2012; 67:470–79. https://doi.org/10.1093/gerona/gls089 [PubMed]
  • 10. Wellcome Trust Case Control C, and Wellcome Trust Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007; 447:661–78. https://doi.org/10.1038/nature05911 [PubMed]
  • 11. Do R, Willer CJ, Schmidt EM, Sengupta S, Gao C, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, Beckmann JS, et al. Common variants associated with plasma triglycerides and risk for coronary artery disease. Nat Genet. 2013; 45:1345–52. https://doi.org/10.1038/ng.2795 [PubMed]
  • 12. Lange LA, Hu Y, Zhang H, Xue C, Schmidt EM, Tang ZZ, Bizon C, Lange EM, Smith JD, Turner EH, Jun G, Kang HM, Peloso G, et al, and NHLBI Grand Opportunity Exome Sequencing Project. Whole-exome sequencing identifies rare and low-frequency coding variants associated with LDL cholesterol. Am J Hum Genet. 2014; 94:233–45. https://doi.org/10.1016/j.ajhg.2014.01.010 [PubMed]
  • 13. Giger JN, Strickland OL, Weaver M, Taylor H, Acton RT. Genetic predictors of coronary heart disease risk factors in premenopausal African-American women. Ethn Dis. 2005; 15:221–32. [PubMed]
  • 14. Budovsky A, Craig T, Wang J, Tacutu R, Csordas A, Lourenço J, Fraifeld VE, de Magalhães JP. LongevityMap: a database of human genetic variants associated with longevity. Trends Genet. 2013; 29:559–60. https://doi.org/10.1016/j.tig.2013.08.003 [PubMed]
  • 15. Willcox BJ, Donlon TA, He Q, Chen R, Grove JS, Yano K, Masaki KH, Willcox DC, Rodriguez B, Curb JD. FOXO3A genotype is strongly associated with human longevity. Proc Natl Acad Sci USA. 2008; 105:13987–92. https://doi.org/10.1073/pnas.0801030105 [PubMed]
  • 16. Broer L, Buchman AS, Deelen J, Evans DS, Faul JD, Lunetta KL, Sebastiani P, Smith JA, Smith AV, Tanaka T, Yu L, Arnold AM, Aspelund T, et al. GWAS of longevity in CHARGE consortium confirms APOE and FOXO3 candidacy. J Gerontol A Biol Sci Med Sci. 2015; 70:110–18. https://doi.org/10.1093/gerona/glu166 [PubMed]
  • 17. Kilic U, Gok O, Erenberk U, Dundaroz MR, Torun E, Kucukardali Y, Elibol-Can B, Uysal O, Dundar T. A remarkable age-related increase in SIRT1 protein expression against oxidative stress in elderly: SIRT1 gene variants and longevity in human. PLoS One. 2015; 10:e0117954. https://doi.org/10.1371/journal.pone.0117954 [PubMed]
  • 18. van Bockxmeer FM. ApoE and ACE genes: impact on human longevity. Nat Genet. 1994; 6:4–5. https://doi.org/10.1038/ng0194-4 [PubMed]
  • 19. Schächter F, Faure-Delanef L, Guénot F, Rouger H, Froguel P, Lesueur-Ginot L, Cohen D. Genetic associations with human longevity at the APOE and ACE loci. Nat Genet. 1994; 6:29–32. https://doi.org/10.1038/ng0194-29 [PubMed]
  • 20. Zeng Y, Nie C, Min J, Liu X, Li M, Chen H, Xu H, Wang M, Ni T, Li Y, Yan H, Zhang JP, Song C, et al. Novel loci and pathways significantly associated with longevity. Sci Rep. 2016; 6:21243. https://doi.org/10.1038/srep21243 [PubMed]
  • 21. Sun L, Hu C, Zheng C, Qian Y, Liang Q, Lv Z, Huang Z, Qi K, Gong H, Zhang Z, Huang J, Zhou Q, Yang Z. FOXO3 variants are beneficial for longevity in Southern Chinese living in the Red River Basin: A case-control study and meta-analysis. Sci Rep. 2015; 5:9852. https://doi.org/10.1038/srep09852 [PubMed]
  • 22. He YH, Lu X, Yang LQ, Xu LY, Kong QP. Association of the insulin-like growth factor binding protein 3 (IGFBP-3) polymorphism with longevity in Chinese nonagenarians and centenarians. Aging (Albany NY). 2014; 6:944–56. https://doi.org/10.18632/aging.100703 [PubMed]
  • 23. Sun L, Hu CY, Shi XH, Zheng CG, Huang ZZ, Lv ZP, Huang J, Wan G, Qi KY, Liang SY, Zhou L, Yang Z. Trans-ethnical shift of the risk genotype in the CETP I405V with longevity: a Chinese case-control study and meta-analysis. PLoS One. 2013; 8:e72537. https://doi.org/10.1371/journal.pone.0072537 [PubMed]
  • 24. Huang J, Sun L, Liu M, Zhou L, Lv ZP, Hu CY, Huang ZZ, Zheng CG, Zhou L, Yang Z. [Association between SIRT1 gene polymorphisms and longevity of populations from Yongfu region of Guangxi]. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2013; 30:55–59. Association between SIRT1 gene polymorphisms and longevity of populations from Yongfu region of Guangxi [PubMed]
  • 25. Scola L, Lio D, Candore G, Forte GI, Crivello A, Colonna-Romano G, Pes MG, Carru C, Ferrucci L, Deiana L, Baggio G, Franceschi C, Caruso C. Analysis of HLA-DRB1, DQA1, DQB1 haplotypes in Sardinian centenarians. Exp Gerontol. 2008; 43:114–18. https://doi.org/10.1016/j.exger.2007.06.007 [PubMed]
  • 26. Scola L, Lio D, Crivello A, Candore G, Forte GI, Colonna-Romano G, Pes MG, Carru C, Ferrucci L, Deiana L, Baggio G, Franceschi C, Caruso C. Analysis of HLA-DQA, HLA-DQB frequencies in a group of Sardinian centenarians. Rejuvenation Res. 2006; 9:157–60. https://doi.org/10.1089/rej.2006.9.157 [PubMed]
  • 27. Feng J, Zhang J, Liu M, Wan G, Qi K, Zheng C, Lv Z, Hu C, Zeng Y, Gregory SG, Yang Z. Association of mtDNA haplogroup F with healthy longevity in the female Chuang population, China. Exp Gerontol. 2011; 46:987–93. https://doi.org/10.1016/j.exger.2011.09.001 [PubMed]
  • 28. Vaarhorst AA, Beekman M, Suchiman EH, van Heemst D, Houwing-Duistermaat JJ, Westendorp RG, Slagboom PE, Heijmans BT, and Leiden Longevity Study (LLS) Group. Lipid metabolism in long-lived families: the Leiden Longevity Study. Age (Dordr). 2011; 33:219–27. https://doi.org/10.1007/s11357-010-9172-6 [PubMed]
  • 29. Heijmans BT, Beekman M, Houwing-Duistermaat JJ, Cobain MR, Powell J, Blauw GJ, van der Ouderaa F, Westendorp RG, Slagboom PE. Lipoprotein particle profiles mark familial and sporadic human longevity. PLoS Med. 2006; 3:e495. https://doi.org/10.1371/journal.pmed.0030495 [PubMed]