|Year : 2017 | Volume
| Issue : 14 | Page : 1677-1682
Variants in the Promoter Region of HLA-DQA1 were Associated with Idiopathic Membranous Nephropathy in a Chinese Han Population
Xiao-Song Qin1, Jian-Hua Liu1, Guan-Ting Lyu2, Meng-Le Peng3, Fu-Ning Yang4, Dong-Chun Qin3, Yong-Zhe Li4, Yong Liu1
1 Department of Medical Laboratory, Shengjing Hospital of China Medical University, Shenyang, Liaoning 110011, China
2 Department of Blood Transfusion, Tangdu Hospital, The Fourth Military Medical University, Xi'an, Shaanxi 710038, China
3 Department of Medical Laboratory, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan 450052, China
4 Department of Rheumatology and Clinical Immunology, Peking Union Medical College Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Key Laboratory of Rheumatology and Clinical Immunology, Ministry of Education, Beijing 100005, China
|Date of Submission||04-Jan-2017|
|Date of Web Publication||7-Jul-2017|
Department of Medical Laboratory, Shengjing Hospital of China Medical University, 36# Sanhao Street, Heping District, Shenyang, Liaoning 110004
Source of Support: None, Conflict of Interest: None
Background: Idiopathic membranous nephropathy (IMN) is an autoimmune disease and the leading cause of adult nephritic syndrome. HLA-DQA1 had been identified to be associated with IMN in Europeans and the result was replicated in Chinese Han population. In this study, six single nucleotide polymorphisms (SNPs) in the promoter of HLA-DQA1 and other two SNPs with IgA nephropathy were included for the association analysis.
Methods: The SNPs were genotyped in 509 patients and 601 controls by the MassArray iPLEX. The quantification of anti-phospholipase A2 receptor (PLA2R) antibodies in sera of IMN patients was performed by anti-PLA2R ELISA (IgG) kit.
Results: After analysis, four SNPs were significantly associated with IMN, with rs2187668 and rs28383345 as the top two signals (P = 8.42×10-5 and 2.48×10-5, respectively). Even under dominant model, the two SNPs were still significantly associated with IMN (P = 3.50×10-3 for rs28383345 and P = 6.55×10-5 for rs2187668). After conditional study with rs2187668, rs28383345 was the only variant significantly correlated with IMN after Bonferroni correction (P = 0.016). The minor alleles of the two SNPs were also mutually exclusive in our cohort. This indicated that the two SNPs were independently associated with IMN in Chinese Han population. Levels of anti-PLA2R autoantibodies were correlated with the genotypes of the two SNPs, but not significantly (P>0.05).
Conclusions: Our results revealed that a novel independent variant in the promoter of HLA-DQA1 was associated with IMN in Chinese Han population. The locus possessed regulatory role according to the data of RegulomeDB. The exact role of the SNPs on the expression of HLA-DQA1 needs further investigation.
Keywords: Autoimmune Disease; HLA-DQA1; Idiopathic Membranous Nephropathy; Promoter; Chinese
|How to cite this article:|
Qin XS, Liu JH, Lyu GT, Peng ML, Yang FN, Qin DC, Li YZ, Liu Y. Variants in the Promoter Region of HLA-DQA1 were Associated with Idiopathic Membranous Nephropathy in a Chinese Han Population. Chin Med J 2017;130:1677-82
|How to cite this URL:|
Qin XS, Liu JH, Lyu GT, Peng ML, Yang FN, Qin DC, Li YZ, Liu Y. Variants in the Promoter Region of HLA-DQA1 were Associated with Idiopathic Membranous Nephropathy in a Chinese Han Population. Chin Med J [serial online] 2017 [cited 2018 Apr 26];130:1677-82. Available from: http://www.cmj.org/text.asp?2017/130/14/1677/209884
| Introduction|| |
Idiopathic membranous nephropathy (IMN) is considered as an autoimmune disease, which is the leading cause of adult nephritic syndrome., The incidence of IMN was increased significantly from 6.48% to 22.79% during 1997–2011; especially in 2009–2011, IMN has become the second common primary glomerular diseases after IgA nephropathy (IgAN) in Chinese. In recent years, phospholipase A2 receptor 1 (PLA2R1) and thrombospondin type-1 domain-containing 7A were identified as the two major antigens of IMN in adults. Immune deposits containing antibodies and their targeted antigens were formed on the outer layer of the glomerular basement membrane, which would evoke the complement activation cascade resulting in podocyte injury and urinary protein loss. Although IMN is not a typical hereditary disease, some researches had confirmed the importance of genetic factors in its pathogenesis.,, Especially, a genome-wide association study (GWAS) for Caucasians identified HLA-DQA1 and PLA2R1 as the two most closely associated signals with IMN.
HLA-DQA1 is located on chromosome 6p21, belonging to the group of major histocompatibility complex (MHC) Class II genes. A major function of MHC Class II molecule is to present antigens derived from extracellular proteins for recognition by T-cells. Since 1970s, genes in the HLA locus have been implicated to be associated with IMN.,HLA-DR2, DR3, B8, and B18 were reported to be associated with IMN in Japanese and Caucasians.,, To date, most studies about association between HLA variants and risk of IMN have been performed in Westerners; only a few reports for Chinese Han population are available., In later studies, merely, the proxy single nucleotide polymorphism (SNP) rs2187668 of HLA-DQA1 was detected in different cohorts of IMN patients.
In this study, we detected SNPs in the promoter region of HLA-DQA1 in a relative large Chinese Han population to evaluate whether polymorphisms were associated with the susceptibility to IMN. After curating the clinical information, we found that about 2.95% (15/509) IMNs were coupled with mild mesangial proliferative IgAN. Therefore, two other SNPs of the HLA region were also included in our analysis.
Except the previously reported variant rs2187668, the SNPs rs28383345, rs36173887, and rs1794275 were also significantly associated with IMN [Table 1]. Conditional analysis indicated that rs2187668 and rs28383345 were independently associated with the susceptibility of IMN. The minor alleles of these SNPs were correlated with levels of anti-PLA2R in sera, but not significantly. Through this study, a novel variant highly linked with IMN was identified in our cohort of Chinese Han population. The exact functional impact of this SNP on HLA-DQA1 still needs further study.
|Table 1: Association of HLA-DQA1 variants with IMN in Chinese Han population|
Click here to view
| Methods|| |
This study was approved by the Ethics Committee of Shengjing Hospital of China Medical University. All methods were carried out in accordance with the approved protocol guidelines of the Shengjing Hospital. Written informed consents were obtained from all participates in this study.
A total of 509 patients were enrolled from five hospitals (Shengjing Hospital of China Medical University, Peking Union Medical College Hospital, First Affiliated Hospital of Zhengzhou University, Second Hospital of Jilin University, and Yan'an People's Hospital) in different geographical areas of China between December 2014 and April 2016. The diagnosis of IMN was confirmed by renal biopsy. Cases of secondary membranous nephropathy were removed, which included membranous nephropathies induced by other autoimmune diseases, diabetes, infections of hepatitis B virus, malaria, and other pathogens or malignancies. The age- and gender-matched healthy individuals (n = 601) were recruited from Peking Union Medical College Hospital. All participants were genetically unrelated Han Chinese, determined by self-report. The clinical features of participants are listed in [Table 2].
Selection of single nucleotide polymorphisms
The 5 kb region (chr6:32601004–32606003) around the promoter of HLA-DQA1 was analyzed using RegulomeDB version 1.1 to identify DNA features and regulatory elements. Only SNPs with a score <3 were selected for subsequent screening. Next, SNPs with an MAF >1% were selected based on the genotyping data of Chinese Han population from 1000 Genome Projects. Two variants associated with IgAN were also included for our analysis to verify whether there existed a common pathway shared by IMN and IgAN. Totally, eight SNPs were selected for the association study. Among them, SNP rs2187668 has been identified to be significantly associated with IMN in the previous GWAS research.
DNA extraction and genotyping
Genomic DNAs of the participants were extracted from the peripheral blood using TIANamp Blood DNA Kit (Tiangen Biotech, Beijing, China) according to the manufacturer's instructions.
SNP genotyping was performed using the MassArray system (Sequenom, San Diego, CA, USA). The data were analyzed by the MassArray Typer software version 4.0 (Sequenom, San Diego, CA, USA). Assay Designer version 4.0 was used to design primers for multiplexed PCR and single-base extension reactions. Briefly, about 10 ng of genomic DNA was used for genotyping by MassArray iPLEX system at Beijing DNALead Co. Ltd., according to the manufacturer's instructions. The reactions were desalted and transferred to a 384-element SpectroCHIP array. Allele identification was obtained using MALDI-TOF mass spectrometry. The mass spectrograms and genotype data were resolved by the MassArray Typer software version 4.0 (Sequenom).
Detection of anti-phospholipase A2 receptor autoantibodies
The quantification of anti-PLA2R antibodies in sera of IMN patients was performed by anti-PLA2R ELISA (IgG) kit (Euroimmun, Luebeck, Germany) according to the manufacturer's instructions. The cutoff value for positive result was set to 20 U/ml for anti-PLA2R antibodies referred by the manufacturer's handbook.
The distributions of age, items for routine test, and levels for anti-PLA2R were analyzed through descriptive statistics of SPSS version 13.0 (SPSS Inc., USA). Values for mean and standard deviation were presented.
The statistical analysis was primarily carried out as previously reported., A PLINK tool set (Broad Institute, MA, USA) was used to conduct the association analysis. The Chi-square test was applied to compare allele frequencies between case and control groups under dominant and recessive models. The sex and age of samples were used as covariates. The odds ratio (OR) and 95% confidence interval were also calculated. SNPs with a P < 0.05 were taken as significantly associated with IMN. Conditional logistic regression analyses were also undertaken by including the most significant SNP as a covariate in calculation.
| Results|| |
Variants in the promoter region of HLA-DQA1 were strongly associated with idiopathic membranous nephropathy in Chinese Han population
As for our case group, the cases were predominantly male (male/female = 1.97:1), with a mean age of 48.1 ± 14.1 years. In the group of healthy controls, there were 403 men and 198 women (men/women = 2.03:1), with a mean age of 48.6 ± 14.0 years. No significant differences were observed in sex and age distributions (P > 0.05) [Table 2]. In the case group, the positivity for anti- PLA2R autoantibodies was about 68.9% among the first-visits. The frequency dropped to 22.5% for those treated. As for the renal pathology, more than half of the cases (55.92%) were at stage 2 based on biopsy, followed by 16.37% in stage 1, 3.02% in stage 3.
To investigate the promoter region of HLA-DQA1 with IMN in Chinese Han population, six SNPs were detected in 509 IMN patients and 601 healthy controls. Unexpectedly, two SNPs (rs72848263 and rs15222936) had no polymorphism in our cohort, which were inconsistent with results of the 1000 Genomes Project for Chinese Han population. Among the 103 individuals of Chinese Han population in Beijing, the MAFs for rs72848263 and rs15222936 were 0.238 and 0.039, respectively. As for the remaining four SNPs in the promoter region of HLA-DQA1, all of them were significantly associated with IMN based on the criteria with a P < 0.05 by Cochran–Armitage Trend Test. Among these SNPs, rs28383345 and rs2187668 were the top two signals associated with IMN (P = 2.48 × 10−5 for rs28383345 and P = 8.42 × 10−5 for rs2187668) [Table 1]. After conditional analysis, rs28383345 was still positively linked with IMN after Bonferroni correction (PBonf = 0.016) [Table 3].
Associations of single nucleotide polymorphisms with idiopathic membranous nephropathy under dominant and recessive models
The associations of these SNPs with IMN were also analyzed under dominant and recessive models. Except the two SNPs without polymorphism, four SNPs of HLA-DQA1 were significantly associated with IMN under dominant model, still with rs28383345 (P = 3.50 × 10−3) and rs2187668 (P = 6.55 × 10−5) as the top two signals [Table 4]. It was consistent with the results under allelic model [Table 1]. As for the recessive model, rs1794275 was marginally associated with IMN in our cohort (P = 0.043). However, since there were no homozygous genotypes of minor alleles for rs36173887, rs28383345, and 2187668 in our cohort, the associations with IMN of these SNPs under recessive models could not be calculated.
IgA nephropathy-associated variant linked marginally with idiopathic membranous nephropathy
Two previously reported IgAN-associated SNPs of HLA in Chinese Han population were also detected in our IMN cohort. Rs2523946 in the HLA-A region had negative relationship with the susceptibility of IMN (P = 0.429). Only rs1794275 was marginally associated with IMN in the current cohort (P = 0.013) [Table 1].
Most of the selected variants were functional
The functional characters of SNPs were explored by RegulomeDB [Table 1]. Rs2523946 was a reported cis-eQTL for HLA-A in lymphoblastoids and fibroblasts. Rs72848263, rs114929610, rs36173887, rs115222936, and rs28383345 could affect the binding motifs of transcription factors determined by footprinting in RegulomeDB. Interestingly, the previously reported IMN-associated variant, rs2187668, was a cis-eQTL for HLA-DQA1, HLA-DQB1, HLA-DRB1, and HLA-A in lymphoblastoid and monocytes. Rs28383345 was located at -2 site of the Open reading frame (ORF) of HLA-DQA1, which was the motif of the Kozak sequence influencing translational efficacy.
Rs2187668 and rs28383345 were correlated, but not significantly, with anti-phospholipase A2 receptor levels
For IMN individuals, anti-PLA2R antibodies were positive in 69.76% of patients among the first visits. After comparing anti-PLA2R level with the genotypes of rs2187668 and rs28383345, the minor alleles of these two SNPs were not significantly associated with the levels of anti-PLA2R (P = 0.163 for rs2187668, P = 0.184 for rs28383345). However, the minor allele T of rs2187668 had the higher level of anti-PLA2R than that of the major allele C (195.90 ± 62.58 vs. 142.30 ± 22.48). Minor allele A of rs28383345 was with lower level of anti-PLA2R (73.03 ± 19.35 vs. 104.20 ± 13.65) [Figure 1] and [Table 5].
|Figure 1: Correlation between genotypes of single nucleotide polymorphisms and levels anti-phospholipase A2 receptor. PLA2R: Phospholipase A2 receptor 1.|
Click here to view
| Discussion|| |
The variant rs2187668 in HLA-DQA1 had been reported to be associated with the susceptibility of IMN;,, however, the causal SNPs of this gene had not been located. To find out if the variants in the promoter region of HLA-DQA1 were associated with IMN, four SNPs in the promoter of HLA-DQA1 and the proxy SNP rs2187668 were genotyped in a cohort of IMN (509 patients and 601 healthy controls).
Except the proxy SNP rs2187668, other four SNPs were also associated significantly with IMN [Table 1]. Even under dominant model, these SNPs were still linked with the susceptibility of the disease [Table 4]. After conditioning analysis with rs2187668, only rs28383345 was significantly associated the susceptibility of IMN (PBonf = 0.01552). Interestingly, the risk allele T of rs2187668 was contributive to the disease (OR = 1.956 [1.393–2.746]), the risk allele A of rs28383345 protective against the susceptibility of IMN (OR = 0.418 [0.276–0.634]). Besides, in our sample cohort, the risk alleles of these two SNPs were mutually exclusive, indicating that these two SNPs were independently associated with IMN.
Very impressively, after curating the clinical information of our project, about 2.95% (15/509) IMNs were coupled with mild mesangial proliferative IgAN in our disease group. It is inferred that variants associated with IgAN might be also linked with the susceptibility of IMN. In this study, two previously reported IgAN-associated SNPs of the MHC gene were discovered to be linked with the pathogenesis of IMN. Rs2523946 in the HLA-A had a negative relationship with the susceptibility of IMN (P = 0.429). Only rs1794275 in the HLA-DQB1 was marginally associated with IMN in the current cohort (P = 0.013). However, these two variants were strongly associated with IgAN in Chinese (P = 1.74 × 10−11 for rs2523946 and P = 3.43 × 10−13 for rs1794275). This implied that the mechanisms for IgAN might be different from that for IMN at least in the HLA although these two diseases shared a common pathway for the pathogenic etiology.
According to eQTL and chromatin modification signatures of RegulomeDB, most of the selected variants possessed regulatory function. Rs2187668 was a cis-eQTL, affecting the expression of HLA-DQA1, HLA-DQB1, HLA-DRB1, and HLA-A in lymphoblastoid and monocytes. Rs28383345 was at the -2 site of the ORF of HLA-DQA1. This site was at the motif of the well-known Kozak sequence, which could affect the translational efficacy of protein coding genes. As Kozak sequence plays a major role in the initiation of the translation process in eukaryotic cells, the minor allele of rs28383345 might be able to affect the translation of HLA-DQA1.
Since anti-PLA2R1 was the most prevalent autoantibodies in IMN, with a prevalence at about 75% in IMN,, the relationship between variants of HLA-DQA1 and levels of anti-PLA2R was explored in our cohort. Although with a suggestive trend, the association of two SNPs was not significantly correlated with anti-PLA2R levels (P = 0.163 for rs2187668 and P = 0.184 for rs28383345) [Table 5]. This might be due to the inadequate number of first-visit nontreated IMN patients with anti-PLA2R detected.
In this study, we assessed, for the first time, the putative association between variants in the regulatory region of HLA-DQA1 and the pathogenesis of IMN. We applied RegulomeDB to identify variants in regulatory elements of HLA-DQA1 with a score <3a. The motifs of these categories were verified by footprinting and worthy of further experimental studies. Though we identified novel polymorphisms associated with IMN, there were still existed limitations for this study. First, the signals should be detected in a new cohort with larger samples to verify the relationship between variants and IMN. Second, the interplay between regulatory variants and the missense polymorphisms of HLA-DQA1 should be studied; it might be highly possible to identify haplotypes more strongly associated with IMN.
From this study, we selected SNPs in the promoter region of HLA-DQA1 through RegulomeDB and analyzed their association with IMNs in a cohort of Chinese Han individuals. Except the previously reported variant associated with IMN, we also found a novel SNP, rs28383345 in the motif of Kozak sequence, associated with the disease. The functional impact of this SNP on the expression of HLA-DQA1 is undergoing.
We thank all participants recruited for this study.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Beck LH Jr., Salant DJ. Membranous nephropathy: From models to man. J Clin Invest 2014;124:2307-14. doi: 10.1172/JCI72270.
Wang XP, Hu ZX, Guo DY, Tao Y. Remission of refractory membranous nephropathy by low-dose rituximab: A case report. Chin Med J 2016;129:871-3. doi: 10.4103/0366-6999.178967.
] [Full text]
Pan X, Xu J, Ren H, Zhang W, Xu Y, Shen P, et al.
Changing spectrum of biopsy-proven primary glomerular diseases over the past 15 years: A single-center study in China. Contrib Nephrol 2013;181:22-30. doi: 10.1159/000348638.
Beck LH Jr., Bonegio RG, Lambeau G, Beck DM, Powell DW, Cummins TD, et al.
M-type phospholipase A2 receptor as target antigen in idiopathic membranous nephropathy. N Engl J Med 2009;361:11-21. doi: 10.1056/NEJMoa0810457.
Ronco P, Debiec H. Pathophysiological advances in membranous nephropathy: Time for a shift in patient's care. Lancet 2015;385:1983-92. doi: 10.1016/S0140-6736(15)60731-0.
Vaughan RW, Demaine AG, Welsh KI. A DQA1 allele is strongly associated with idiopathic membranous nephropathy. Tissue Antigens 1989;34:261-9.
Bantis C, Heering PJ, Aker S, Siekierka M, Kuhr N, Grabensee B, et al.
Tumor necrosis factor-alpha gene G-308A polymorphism is a risk factor for the development of membranous glomerulonephritis. Am J Nephrol 2006;26:12-5. doi: 10.1159/000090706.
Kim S, Chin HJ, Na KY, Kim S, Oh J, Chung W, et al.
Single nucleotide polymorphisms in the phospholipase A2 receptor gene are associated with genetic susceptibility to idiopathic membranous nephropathy. Nephron Clin Pract 2011;117:c253-8. doi: 10.1159/000320194.
Stanescu HC, Arcos-Burgos M, Medlar A, Bockenhauer D, Kottgen A, Dragomirescu L, et al.
Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy. N Engl J Med 2011;364:616-26. doi: 10.1056/NEJMoa1009742.
Schmidt H, Williamson D, Ashley-Koch A. HLA-DR15 haplotype and multiple sclerosis: A HuGE review. Am J Epidemiol 2007;165:1097-109. doi: 10.1093/aje/kwk118.
Klouda PT, Manos J, Acheson EJ, Dyer PA, Goldby FS, Harris R, et al.
Strong association between idiopathic membranous nephropathy and HLA-DRW3. Lancet 1979;2:770-1.
Hiki Y, Kobayashi Y, Itoh I, Kashiwagi N. Strong association of HLA-DR2 and MT1 with idiopathic membranous nephropathy in Japan. Kidney Int 1984;25:953-7.
Le Petit JC, Laurent B, Berthoux FC. HLA-DR3 and idiopathic membranous nephritis (IMN) association. Tissue Antigens 1982;20:227-8.
Tomura S, Kashiwabara H, Tuchida H, Shishido H, Sakurai S, Miyajima T, et al.
Strong association of idiopathic membranous nephropathy with HLA-DR2 and MT1 in Japanese. Nephron 1984;36:242-5.
Lv J, Hou W, Zhou X, Liu G, Zhou F, Zhao N, et al.
Interaction between PLA2R1 and HLA-DQA1 variants associates with anti-PLA2R antibodies and membranous nephropathy. J Am Soc Nephrol 2013;24:1323-9. doi: 10.1681/ASN.2012080771.
Cui G, Zhang L, Xu Y, Cianflone K, Ding H, Wang DW. Development of a high resolution melting method for genotyping of risk HLA-DQA1 and PLA2R1 alleles and ethnic distribution of these risk alleles. Gene 2013;514:125-30. doi: 10.1016/j.gene.2012.11.004.
Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, et al.
Annotation of functional variation in personal genomes using RegulomeDB. Genome Res 2012;22:1790-7. doi: 10.1101/gr.137323.112.
Li P, Lu G, Cui Y, Wu Z, Chen S, Li J, et al.
Association of IL12A Expression Quantitative Trait Loci (eQTL) With primary biliary cirrhosis in a Chinese Han population. Medicine (Baltimore) 2016;95:e3665. doi: 10.1097/MD.0000000000003665.
Li P, Lu G, Wang L, Cui Y, Wu Z, Chen S, et al.
A rare nonsynonymous variant in the lipid metabolic gene HELZ2 related to primary biliary cirrhosis in Chinese Han. Allergy Asthma Clin Immunol 2016;12:14. doi: 10.1186/s13223-016-0120-6.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al.
PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559-75. doi: 10.1086/519795.
Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J, et al.
Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 2010;464:773-7. doi: 10.1038/nature08903.
Zeller T, Wild P, Szymczak S, Rotival M, Schillert A, Castagne R, et al.
Genetics and beyond – The transcriptome of human monocytes and disease susceptibility. PLoS One 2010;5:e10693. doi: 10.1371/journal.pone.0010693.
Bullich G, Ballarín J, Oliver A, Ayasreh N, Silva I, Santín S, et al.
HLA-DQA1 and PLA2R1 polymorphisms and risk of idiopathic membranous nephropathy. Clin J Am Soc Nephrol 2014;9:335-43. doi: 10.2215/CJN.05310513.
Kozak M. An analysis of 5'-noncoding sequences from 699 vertebrate messenger RNAs. Nucleic Acids Res 1987;15:8125-48.
Gopalakrishnan N, Abeesh P, Dineshkumar T, Murugananth S, Sakthirajan R, Raman GS, et al.
Prevalence of serum anti M-type phospholipase A2 receptor antibody in primary membranous nephropathy: A single center experience. Indian J Nephrol 2016;26:257-61. doi: 10.4103/0971-4065.160334.
] [Full text]
Xu NX, Xie QH, Sun ZX, Wang J, Li Y, Wang L, et al.
Renal phospholipase A2 receptor and the clinical features of idiopathic membranous nephropathy. Chin Med J 2017;130:892-8. doi: 10.4103/0366-6999.204096.
] [Full text]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5]