Article Text

Download PDFPDF

Response to: ‘Questions on ‘Sequencing of the MHC region defines HLA-DQA1 as the major genetic risk for seropositive rheumatoid arthritis in Han Chinese population’ by Guo et al’ by Regueiro and Gonzalez
Free
  1. Jianping Guo1,
  2. Tao Zhang2,3,
  3. Xiaowei Li3,4,
  4. Hongzhi Cao5,6,7,
  5. Zhanguo Li1
  1. 1Department of Rheumatology and Immunology, Peking University People's Hospital, Beijing, China
  2. 2Institute of Precision Health Research, Beijing genomics institute (BGI)-Shenzhen, Shenzhen, China
  3. 3Institute of Precision Health Research, China National GeneBank-Shenzhen, BGI-Shenzhen, Shenzhen, China
  4. 4Technology Service Center, Beijing genomics institute (BGI)-shenzhen, shenzhen, China
  5. 5Division of Immunology, Shenzhen Digital Life Institute, Shenzhen, China
  6. 6Division of Artificial Intelligence and Bioinformatics, iCarbonX, ShenZhen, China
  7. 7Division of Birth Cohort Study, Guangzhou Women and Children’s Medical Center, Guangzhou Medical University, Guangzhou, China
  1. Correspondence to Professor Jianping Guo, Department of Rheumatology and Immunology, Peking University People's Hospital, Beijing 100044, China; jianping.guo{at}bjmu.edu.cn

Statistics from Altmetric.com

We appreciate Dr Gonzalez’s interest and comments on our recent publication ‘Sequencing of the major histocompatibility complex (MHC) region defines human leukocyte antigen (HLA)-DQA1 as the major genetic risk for seropositive rheumatoid arthritis in Han Chinese population’.1 2 Dr Gonzalez’s comments provide us with an opportunity to clarify and discuss the frequencies of amino-acids at position DQα1:160 and the protective association of DRβ1:37N in rheumatoid arthritis (RA), and to improve our study.

One of the concerns Dr Gonzalez expressed is the frequencies of three amino-acids (Asp (D), Ala (A) and Ser (S)) at DQα1:160 do not sum up to 1.0, that is, 0.20, 0.22 and ≈0.02 in controls; 0.36, 0.37 and ≈0.01 in RA patients. The explanation is that the ‘minor frequencies’ was set as default for all variants in PLINK. The original frequencies of three amino-acids (D, A and S) were 0.20, 0.78 (1–0.22) and ≈0.02 (the sum is 1.00) in controls and 0.36, 0.63 (1–0.37) and ≈0.01 (the sum is 1.00) in cases, respectively. In our original paper the frequencies of DQα1:160A inhealthy controls were similar to those reported in other Asian studies.3 These results do not affect the calculation of p value, but do affect the odd ratio (OR) calculation. Indeed, by the omnibus test DQα1:160A showed a protective effect (OR=0.46, p=2.72 x 10-35, online supplementary table 10 in Guo et al).2 We appreciate Dr Gonzalez et al for this important point and have made a correction for our publication, in which all variants have been presented according to original frequencies instead of minor frequencies.4

Regarding the protective effect of DRβ1:37N, although the identified amino-acid DRβ1:37N did not show any significant association in univariate regression analysis, it reached second strong statistical significance after conditioning on DQα1:160D in both discovery and validation stages, indicating an independent association. This phenomenon could be potentially explained by the Simpson's paradox, a striking observation that an association between two variables at the population-level might increase or decrease in quantity, or even change direction within the subgroups, depending on the set of variables being controlled,5 6 and has been reported in several genetic association studies.7 8 Notably, the DRβ1:11D also showed an independent protective effect and was in high linkage disequilibrium (LD) with DRβ1:37N (r2=0.62; online supplementary tables 8 and 9 in Guo et al).

Regarding other DRB1 variants, as the author indicated, by omnibus test we replicated the findings reported in previous studies,9–12 including the position 11 and 13 at DRβ1, and the allele DRB1*04:05. However, our study focused on single nucleotide polymorphisms (SNPs), classical HLA alleles and the individual amino-acid variants rather than amino-acid positions, because a particular amino-acid(s) may have potential biological function(s). Furthermore, different amino-acids at same position may insert different functions.13 Taking this into consideration, DQα1:160D remained the top association in omnibus test (OR=2.30, p=1.82 x 10−38) (online supplementary table 10 in Guo et al). Furthermore, consistent with our findings, Hirata et al14 have also reported that one of DQα1:160D encoding allele DQA1*0303 was a strong risk for susceptibility to RA in Japanese population (OR=2.65, p=2.0 × 10−173, shown in table 1 in Hirata et al).

References

Footnotes

  • Handling editor Josef S Smolen

  • Contributors JG undertook the primary duties in writing and editing the manuscript. TZ, XWL and HC undertook the duties in checking the primary data. ZL oversaw the manuscript writing and edited the manuscript.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Patient consent for publication Not required.

  • Provenance and peer review Commissioned; internally peer reviewed.

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Linked Articles