key: cord-0746639-upk6jqll
authors: Khor, Seik-Soon; Omae, Yosuke; Takeuchi, Junko S.; Fukunaga, Ami; Yamamoto, Shohei; Tanaka, Akihito; Matsuda, Kouki; Kimura, Moto; Maeda, Kenji; Ueda, Gohzoh; Mizoue, Tetsuya; Ujiie, Mugen; Mitsuya, Hiroaki; Ohmagari, Norio; Sugiura, Wataru; Tokunaga, Katsushi
title: An Association Study of HLA with the Kinetics of SARS-CoV-2 Spike Specific IgG Antibody Responses to BNT162b2 mRNA Vaccine
date: 2022-04-05
journal: Vaccines (Basel)
DOI: 10.3390/vaccines10040563
sha: 30802da73c7bcb7dcf140741ed26fce730154946
doc_id: 746639
cord_uid: upk6jqll

BNT162b2, an mRNA-based SARS-CoV-2 vaccine (Pfizer-BioNTech, New York, NY, USA), is one of the most effective COVID-19 vaccines and has been approved by more than 130 countries worldwide. However, several studies have reported that the COVID-19 vaccine shows high interpersonal variability in terms of humoral and cellular responses, such as those with respect to SARS-CoV-2 spike protein immunoglobulin (Ig)G, IgA, IgM, neutralizing antibodies, and CD4(+) and CD8(+) T cells. The objective of this study is to investigate the kinetic changes in anti-SARS-CoV-2 spike IgG (IgG-S) profiles and adverse reactions and their associations with HLA profiles (HLA-A, -C, -B, -DRB1, -DQA1, -DQB1, -DPA1 and -DPB1) among 100 hospital workers from the Center Hospital of the National Center for Global Health and Medicine (NCGM), Tokyo, Japan. DQA1*03:03:01 (p = 0.017; Odd ratio (OR) 2.80, 95%confidence interval (CI) 1.05–7.25) was significantly associated with higher IgG-S production after two doses of BNT162b2, while DQB1*06:01:01:01 (p = 0.028, OR 0.27, 95%CI 0.05–0.94) was significantly associated with IgG-S declines after two doses of BNT162b2. No HLA alleles were significantly associated with either local symptoms or fever. However, C*12:02:02 (p = 0.058; OR 0.42, 95%CI 0.15–1.16), B*52:01:01 (p = 0.031; OR 0.38, 95%CI 0.14–1.03), DQA1*03:02:01 (p = 0.028; OR 0.39, 95%CI 0.15–1.00) and DPB1*02:01:02 (p = 0.024; OR 0.45, 95%CI 0.21–0.97) appeared significantly associated with protection against systemic symptoms after two doses of BNT162b2 vaccination. Further studies with larger sample sizes are clearly warranted to determine HLA allele associations with the production and long-term sustainability of IgG-S after COVID-19 vaccination.

The concept of messenger RNA (mRNA) vaccination stems from 1987, when Robert Malone confirmed that human cells can absorb cationic liposomes containing mRNA and can create proteins from those mRNA sequences [1] . However, over the years, academic laboratories and companies working on mRNA have come to the consensus that mRNA is too prone to degradation to be used effectively as a drug or a vaccine [2] . From the end of the 20th century, research into mRNA vaccines mainly focused on influenza disease [3] and cancer [4] , with testing in animal models yielding satisfactory results. Such results inspired CureVac, BioNTech, and Moderna to focus on transforming mRNA into a drug platform [5] [6] [7] . In 2005, the Karikó and Weissman team performed a landmark experiment in which an mRNA vaccine was used to successfully suppress RNA recognition by Toll-like receptors (TLRs) using modified pseudouridine, a uridine analog [8] . All this previous research contributed to the rapid development of coronavirus disease 2019 (COVID- 19) mRNA vaccines within days of the genome for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) being released online [9] .

BNT162b2, an mRNA-based SARS-CoV-2 vaccine (Pfizer-BioNTech), is one of the most effective COVID-19 vaccines [10] and has been approved by more than 130 countries worldwide. However, several studies have reported that the COVID-19 vaccine shows high interpersonal variability in terms of humoral and cellular responses, such as those with respect to anti-SARS-CoV-2 spike immunoglobulin (Ig)G [11] [12] [13] [14] [15] [16] [17] , IgA [15, 16] , IgM [16] , neutralizing antibodies [18] [19] [20] , and CD4 + and CD8 + T cells [18, 19, 21] . The long-term sustainability of SARS-CoV-2 neutralizing antibodies is crucial for the protection against infection and hospitalization.

COVID-19 mRNA vaccines work by introducing the SARS-CoV-2 spike protein (Sprotein) mRNA sequence and subsequently the translated protein to the human immune system. This process first prompts an innate (non-specific) immune response by triggering immune cells such as natural killer cells, macrophages, neutrophils, dendritic cells, mast cells, basophils, and eosinophils. The second line of immune response is the adaptive immune response, in which T and B lymphocytes are activated to generate long-term specific immune responses to the SARS-CoV-2 virus.

Human leukocyte antigen (HLA) allelic variation may be associated with vaccine efficacy according to trials of the hepatitis B vaccine [22, 23] , influenza vaccine [24] , and HIV-1 vaccine [25] . However, data on possible associations with the kinetics of post-COVID-19 vaccination antibodies remain scarce. Here, we designed a longitudinal study aimed at investigating the association of host HLA polymorphisms with kinetic changes in IgG and adverse reactions during and after BNT162b2 vaccination among hospital workers at a national medical institution in Japan.

The National Center for Global Health and Medicine (NCGM) is a Japanese governmentdesignated medical center for the treatment of COVID-19 patients. We recruited and monitored a total of 100 hospital workers ≥20 years of age from the period between March and June 2021. Details of the study design are explained in previous studies [26, 27] . All participants were vaccinated with two doses of the BNT162b2 mRNA-based SARS-CoV-2 vaccine (Pfizer-BioNTech) according to the standard protocol (two doses of 30 µg, administered 3 weeks apart) [18, 28] . All participants received the first dose of the BNT162b2 vaccine in March 2021 and the second dose was administered 21 days after the first dose. Blood was drawn from volunteers on day 1 (immediately after the first dose). Samples were then collected on day 15, day 29 (7 days after the second dose), and day 61. None of the participants had any history of COVID-19 or administration of immunosuppressive medications. The study protocol was approved by the Ethics Committee of the NCGM, Japan (approval number: NCGM-A-004175-00). Written informed consent was obtained from all participants prior to enrollment.

Participants were requested to complete a questionnaire regarding the adverse reactions experienced after each dose of the vaccine. Based on the results, we classified side effects into local symptoms (pain at the injection site, swelling, and redness) and systemic symptoms (fever, dizziness, fatigue, headache, chills, vomiting, diarrhea, muscle pain, and joint pain) with reference to the US Food and Drug Administration guidance. Symptoms were further graded into 4 categories: grade 0, no adverse effect; grade 1, adverse effect without interfering with daily activities; grade 2, adverse effect with some degree of interference with daily activities; grade 3, adverse effect with considerable interference with daily activities. Redness or swelling at the side of infection was graded based on the size of the redness or swelling as grade 0, 0-2.0 cm; grade 1, 2.1-5.0 cm; grade 2, 5.1-10.0 cm; grade 3, >10.0 cm.

Antibodies against the SARS-CoV-2 receptor-binding domain (RBD) of the S1 subunit of the spike protein (IgG-S) (AdviseDx SARS-CoV-2 IgG II; Abbott, Chicago, IL, USA) and SARS-CoV-2 nucleocapsid (IgG-N) (Abbott ARCHITECT ® SARS-CoV-2 anti-N IgG; Abbott) were measured at each collection time point. For IgG-S, results > 50 AU/mL (as the cutoff set by the manufacturer) were considered indicative of seropositivity. For surrogate measurement of neutralizing antibody, a threshold of 4160 AU/mL was applied, as this threshold corresponds to a 95% probability of obtaining a positive result from plaque reduction neutralization test (PRNT) for SARS-CoV-2 at 1:250 dilution [14] . For IgG-N, results above the index value of 1.40 (as the cutoff set by the manufacturer) were considered indicative of seropositivity.

DNA specimens were extracted from peripheral blood mononuclear cells using QIAamp DNA Blood Maxi Kit (QIAGEN, Hilden, Germany) following the manufacturer protocol. The quality of the extracted DNA was checked using Qubit™ dsDNA HS and BR Assay Kits (Invitrogen, Waltham, MA, USA). NGS HLA genotyping was performed using AllType™ NGS Assays (One Lambda, West Hills, CA, USA) on an Ion GeneStudio S5 sequencing system (Thermo Fisher Scientific, Waltham, MA, USA). HLA-A, -C, -B, DRB1, DQA1, DQB1, DPA1, DPB1 targeted gene amplification, HLA library preparation, HLA template preparation, and HLA library loading onto an ion 530v1 chip in an Ion Chef library preparation robot (Thermo Fisher Scientific) and final sequencing in an Ion GeneStudio S5 sequencer (Thermo Fisher Scientific) were performed following the instructions from the vendor.

Demultiplexing of barcodes and base-calling was carried out using Torrent Suite version 5.8.0 (Thermo Fisher Scientific). Raw fastq reads were extracted using the FileExporter function in Torrent Suite version 5.8.0. HLA genotype assignments were undertaken using HLATypeStream Visual (TSV v2.0; One Lambda, West Hills, CA, USA) and NGSengine ® (v2.22.0.22581; GenDX, Utrecht, The Netherlands). The default analysis parameters and healthy metrics thresholds were applied for TSV v2.0, while we applied the "ignore regions" function in NGSengine ® to eliminate known sequencing error sites in the ion S5 system. Four-field resolution HLA alleles were obtained for HLA-A, -C, -B, DRB1, DQA1, DQB1, and DPA1. Novel HLA alleles were subjected to Pacbio Sequel sequencing in collaboration with the HU Group Research Institute (Tokyo, Japan). Pacbio subreads and consensus reads were obtained from SMRT Link software and HLA calling was performed using NGSengine ® software (v2.22.0.22581; GenDX).

Correlation coefficient tests and Kruskal-Wallis tests were applied to describe the correlation between IgG-S titers with age, sex, BMI of participants, and adverse reactions after two doses of vaccination. IgG-S levels in participants plateaued on day 29 (7 days after the second dose), and were further classified into low responders (top 25th percentile of the IgG-S distribution) and high responders (bottom < 25th percentile of the IgG-S distribution). Strong decliners (top 25th percentile of ∆) and weak decliners (bottom 75th percentile of ∆) were defined as ∆ of IgG-S levels between day 29 (7 days after the second dose) and day 61 (39 days after the second vaccination).

Case-control HLA allele association tests, HLA haplotype estimations, and casecontrol HLA haplotype association tests were prepared and analyzed using the Bridging ImmunoGenomics Data Analysis Workflow Gaps (BIGDAWG) R package [29] .

A total of 100 hospital workers at the NCGM were recruited, comprising 32 men and 68 women. Participants were further classified into five age strata to observe the changes in IgG-S in relation to age: 20-29 years old; 30-39 years old; 40-49 years old; 50-59 years old; ≥ 60 years old ( Table 1 ). The BMI of participants ranged from 17 kg/m 2 to 38 kg/m 2 , with the highest BMI observed in the group ≥60 years old (mean ± standard deviation, 26 ± 6.30 kg/m 2 ). 

All participants tested negative for IgG-S, IgG-N, and IgM on day 1 (before BNT162b2 vaccination). IgG-S levels linearly increased after the first dose of BNT162b2 vaccination, plateauing on day 15. After the second dose of BNT162b2 vaccination, the IgG-S level peaked on day 29 (7 days after the second vaccination). Younger age (especially under 59 years old) (Kruskal-Wallis test, p = 0.024) (Supplementary Figure S1a ) and female sex (Kruskal-Wallis test, p = 0.034) (Supplementary Figure S1b ) were significantly associated with higher IgG-S production after two doses of BNT162b2 vaccinations (Supplementary Figure S1a-c). For the surrogate measurement of neutralizing antibody activities, we examined for IgG-S level > 4160 AU/mL, which corresponded to a 95% probability of obtaining a PRNT ID 50 (estimated number of virus particles required to produce infection in 50% of normal adult humans) at 1:250 dilution [14] . After 15 days from the first dose of BNT162b2 vaccination, none of the participants had reached a sufficient level of surrogate neutralizing activities, although 96% of participants obtained a sufficient level of surrogate neutralizing activities against SARS-CoV-2 by day 29 (7 days after the second dose). However, the proportion showing a sufficient level of surrogate neutralizing activities decreased to 81% on day 61.

On day 61 (39 days after the second vaccination), we observed a significant correlation between age (r = 0.325, p = 0.0012) (Supplementary Figure S2a) and decreased IgG-S levels, but no associations with sex (p = 0.0621) (Supplementary Figure S2b) or BMI (r = 0.173, p = 0.1065) (Supplementary Figure S2c) were apparent for participants. However, although it was not statistically significant, strong decliners (∆ ≥20,000 AU/mL) were observed to be mainly concentrated among the female group (n = 15, 23%), compared with the male group (n = 3, 9.4%).

None of the participants were seropositive for IgG-N across the five survey time points, indicating no prior SARS-CoV-2 infection. No infection was observed during or after the BNT162b2 vaccination.

We classified the self-reported body temperature into four categories: T1, <37. (Table 2 ). Haplotype analysis of HLA-DQA1 and surrounding HLA genes did not reveal any significant associations with HLA haplotypes (DQA1*03:03:01-DQB1*04:01:01, p = 0.072, OR 2.45, 95%CI 0.92-6.52 and DRB1*04:05:01-DQA1*03:03:01-DQB1*04:01:01, p = 0.072, OR 2.45, 95%CI 0.92-6.52) ( Table 3) . Table 2 . Associational analysis of HLA-A, -C, -B, -DRB1, -DQB1, -DOA1 and -DPB1 alleles in 73 IgG-S low responders and 24 IgG-S high responders. Significant HLA alleles are highlighted. Abbreviations: OR, odds ratio; 95%CI, 95% confidence interval; binned, rare HLA alleles with expected count <5 are combined into a common class.

Two-Field Allele Figure S2a ). DQB1*06:01:01:01 (p = 0.028, OR 0.27, 95%CI 0.05-0.94) was significantly associated with IgG-S declines after two doses of BNT162b2 ( Table 5 ). Due to the limited sample size, none of the abovementioned HLA alleles remained significant associated with IgG-S levels after correction for multiple comparisons, but a potential association between individual HLA alleles and responsiveness to BNT162b2 vaccination was still suggested.

No HLA alleles were significantly associated with either local symptoms (Supplementary Table S1 ) or fever (Supplementary Table S2 Table S3 ). 

In this study, we described the kinetic changes in IgG-S profiles and adverse reactions and their associations with HLA profiles among hospital workers from a national medical institution in Japan. BNT162b2 was designed to present the entire spike glycoprotein of SARS-CoV-2 as a target to elicit neutralizing antibodies and so block viral entry through the angiotensin-converting enzyme 2 (ACE2) cellular receptor, particularly the viral sequences recognizing the receptor-binding domain [28] . A structural study of SARS-CoV-2 proteins identified the S-protein as a major inducer of protective immunity [30] . Assessment of changes in the kinetics of antibodies, adverse reactions after BNT162b2 vaccination, and possible associations between interpersonal differences and variations in HLA genes are thus important.

All participants in the present study displayed seroconversion on day 29 (7 days after the second dose) and displayed an increase in IgG-S against SARS-CoV-2. The youngest age group and female participants produced higher levels of IgG-S after two doses of BNT162b2 vaccinations, consistent with findings from previous studies [31, 32] . However, the IgG-S antibody response appears both time-and age-dependent (Supplementary Figure S2) [11, 13, 16, 33] . In this study, the female group produced significantly more IgG-S at the early stage of our longitudinal study (Supplementary Figure S1b ), but we also observed a strong IgG-S decline (∆20,000 AU/mL) in the female group in general (Supplementary Figure S2b ) with largest declines observed in 20-29 y.o. female. DQA1*03:03:01 (p = 0.017; OR 2.80, 95%CI 1.05-7.25) was significantly associated with strong responders for IgG-S, Interestingly, DQA1*03:03:01 is the second most common DQA1 allele in the Japanese population, with a frequency of 16.50% [34] .

The long-term sustainability of IgG-S after COVID-19 vaccination is central to the protection afforded against SARS-CoV-2. This study did not identify any HLA alleles associated with strong declines in IgG-S levels. Instead, we identified an association between DQB1*06:01:01:01 (p = 0.028; OR 0.27, 95%CI 0.05-0.94) and protection against IgG-S declines.

The severity of systemic adverse symptoms after two doses of BNT162b2 vaccination was associated with the level of IgG-S production, consistent with previous research [35] . We observed that carriers of C*12:02:02, B*52:01:01, DQA1*03:02:01, and DPB1*02:01:02 showed significant protection against systemic symptoms.

Data regarding long-term antibody kinetics in vaccinated subjects remain scarce. In a population of 33 healthy adults having received the Moderna mRNA-1273 vaccine and followed up for 209 days, the estimated half-life of the antibody response was 52 days (95%CI 46-58 days), using an exponential decay model [13] . In a cohort of 188 unvaccinated COVID-19 patients (mostly not hospitalized: 174/188) who were followed up for as long as 8 months, the antibody half-life was 83 days (95%CI 62-126 days) [14] . In this study, we estimated that IgG-S levels would drop below 4160 AU/mL [14] after 101 days, suggesting a need for booster shots to sustain adequate levels of neutralizing antibodies.

Research into the associations between HLA alleles and IgG-S antibody levels after vaccination remains limited, but a small Italian cohort (n = 56) found no associations between HLA alleles and IgG-S levels after BNT162b2 vaccination.

In the present study, we identified several HLA alleles showing potential associations with the kinetics of IgG-S antibody levels. Further studies with larger sample sizes are clearly warranted to determine HLA allele associations with the production and long-term sustainability of IgG-S antibody levels after COVID-19 vaccination.

Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/vaccines10040563/s1. Figure S1 : SARS-CoV-2 spike-specific IgG production on day 29, stratified by age and sex: (a) SARS-CoV-2 spike-specific IgG production on day 29, stratified by 5 age groups: 20-29 years old (y.o), 30-39 y.o., 40-49 y.o. 50-59 y.o, and over 60 y.o.; (b) SARS-CoV-2 spike-specific IgG production on day 29, stratified by sex; (c) SARS-CoV-2 spike-specific IgG production on day 29, stratified by 10 groups including the categories listed in (a) and (b). Violin plots show median and 95% confidence interval (CI), and grey dotted line indicates IgG-S level of 4160 AU/mL which corresponded to a 95% probability of obtaining a PRNT ID50 (estimated number of virus particles required to produce infection in 50% of normal adult humans) at 1:250 dilution. Figure S2 : SARS-CoV-2 spike-specific IgG decline on day 61, stratified by age, sex, and BMI of participants: (a) Scatter plot showing the correlation between age and the change in IgG-S between day 29 and day 61. Correlation coefficients and p-values are calculated. Color intensity of the heatmap indicates the accumulation of data points in a specific area; (b) violin plot showing the correlation between sex and the change in IgG-S between day 29 and day 61. Violin plots show median and 95% confidence interval (CI), and grey dotted line indicates IgG-S level of 4160 AU/mL which corresponded to a 95% probability of obtaining a PRNT ID50 (estimated number of virus particles required to produce infection in 50% of normal adult humans) at 1:250 dilution; (c) scatter plot showing the correlation between the BMI of participants and the change in IgG-S between day 29 and day 61. Correlation coefficients and p-values are calculated. Color intensity of the heatmap indicates the accumulation of data points in a specific area. Figure S3 : Association between SARS-CoV-2 spike-specific IgG production on day 29 and adverse reactions after two doses of BNT162b2 vaccination: (a) violin plot shows the distribution of IgG-S production on day 29 stratified by grades of systemic symptoms experienced by participants after two doses of BNT162b2 vaccination; (b) violin plot shows the distribution of IgG-S production on day 29 stratified by fever experienced by participants after two doses of BNT162b2 vaccination. T1: <37.5 • C; T2: 37.5-37.9 • C; T3: 38.0-38.4 • C; and T4: 38.5-38.9 • C; (c) violin plot shows the distribution of IgG-S production on day 29 stratified by the grade of local symptoms experienced by participants after two doses of BNT162b2 vaccination. Violin plots show median and 95% confidence interval (CI), and grey dotted line indicates IgG-S level of 4160 AU/mL which corresponded to a 95% probability of obtaining a PRNT ID50 (estimated number of virus particles required to produce infection in 50% of normal adult humans) at 1:250 dilution. Table S1 . Associational analysis of HLA-A, -C, -B, -DRB1, -DQB1, -DOA1 and -DPB1 alleles in 55 participants with no local adverse effect versus 33 participants with local adverse effect; Table S2 Associational analysis of HLA-A, -C, -B, -DRB1, -DQB1, -DOA1 and -DPB1 alleles in 70 participants with fever versus 18 participants with fever more than 38 • C. Table S3 

Cationic liposome-mediated RNA transfection

mRNA vaccines-A new era in vaccinology

Characterization of a messenger RNA polynucleotide vaccine vector

Dendritic cells pulsed with RNA are potent antigen-presenting cells in vitro and in vivo

Spontaneous cellular uptake of exogenous messenger RNA in vivo is nucleic acid-specific, saturable and ion dependent

Phosphate-enhanced transfection of cationic lipid-complexed mRNA and plasmid DNA

Overexpression of urokinase receptor in mammalian cells following administration of the in vitro transcribed encoding mRNA

Suppression of RNA recognition by Toll-like receptors: The impact of nucleoside modification and the evolutionary origin of RNA

A new coronavirus associated with human respiratory disease in China

BNT162b2 mRNA Covid-19 Vaccine in a Nationwide Mass Vaccination Setting

Antibody titres decline 3-month post-vaccination with BNT162b2

Induction of immune response after SARS-CoV-2 mRNA BNT162b2 vaccination in healthcare workers

Antibody responses to BNT162b2 mRNA COVID-19 vaccine in 2015 healthcare workers in a single tertiary referral hospital in Japan

Antibody responses to the BNT162b2 mRNA vaccine in individuals previously infected with SARS-CoV-2

Human IgG and IgA responses to COVID-19 mRNA vaccines

BNT162b2 mRNA SARS-CoV-2 Vaccine Elicits High Avidity and Neutralizing Antibodies in Healthcare Workers

HLA Does Not Impact on Short-Medium-Term Antibody Response to Preventive Anti-SARS-Cov-2 Vaccine

BNT162b2 vaccine induces neutralizing antibodies and poly-specific T cells in humans

Early T cell and binding antibody responses are associated with COVID-19 RNA vaccine efficacy onset

Correlates of Neutralizing/SARS-CoV-2-S1-binding Antibody Response with Adverse Effects and Immune Kinetics in BNT162b2-Vaccinated Individuals

T-cell responses as a correlate of COVID-19 vaccination. A pilot study in Health Care Workers

Key HLA-DRB1-DQB1 haplotypes and role of the BTNL2 gene for response to a hepatitis B vaccine

The effect of HLA on immunological response to hepatitis B vaccine in healthy people: A meta-analysis

Correlation between human leukocyte antigen class II alleles and HAI titers detected post-influenza vaccination

Analysis of HLA A*02 association with vaccine efficacy in the RV144 HIV-1 vaccine trial

SARS-CoV-2 specific T cell and humoral immune responses upon vaccination with BNT162b2

Association between reactogenicity and SARS-CoV-2 antibodies after the second dose of the BNT162b2 COVID-19 vaccine

Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine

Bridging ImmunoGenomic Data Analysis Workflow Gaps (BIGDAWG): An integrated case-control analysis pipeline

Contributions of the structural proteins of severe acute respiratory syndrome coronavirus to protective immunity

Age-dependent immune response to the Biontech/Pfizer BNT162b2 COVID-19 vaccination

Age-dependent and gender-dependent antibody responses against SARS-CoV-2 in health workers and octogenarians after vaccination with the BNT162b2 mRNA vaccine

Anti-SARS-CoV-2 Receptor-Binding Domain Total Antibodies Response in Seropositive and Seronegative Healthcare Workers Undergoing COVID-19 mRNA BNT162b2 Vaccination

Determination of HLA-A, -C, -B, -DRB1 allele and haplotype frequency in Japanese population based on family study

Association between Immunoglobulin G Levels and Adverse Effects Following Vaccination with the BNT162b2 Vaccine among Japanese Healthcare Workers

We are grateful to the members of the working group for this study (Yoshimi Shigemori, Ayumi Nakayama, Yusuke Oshiro, Natsumi Inamura, Haruka Osawa, Maki Konishi, Azusa Kamikawa, and Yumiko Kito) for their support. We thank Ikue Ito and Tetsuya Sato from HU Group Research Institute for their technical support for Pacbio Sequel sequencing.

All authors (except Gohzoh Ueda) declare no conflict of interest. Gohzoh Ueda is one of the employees of Abbott Japan, which provided the antibody assay reagents and funding for the present study. The role of the funder/sponsor is described above.

Funding: This research is supported by Japan Agency for Medical Research and Development (AMED) under Grant Number JP20kk0205012 and JP20fk0108104, the NCGM Intramural Research Fund 20A2002D, 21A006, 19K059 and 2020-B-09) and Abbott Japan (Grant Number 20C050). The funders did not play any role in the design and conduct of the study; collection, management, analysis, or interpretation of the data; preparation, review, or approval of the manuscript; the decision to submit the manuscript for publication.

The study protocol was approved by the Ethics Committee of the NCGM, Japan (approval number: NCGM-A-004175-00).Informed Consent Statement: Written informed consent was obtained from all participants prior to enrollment.

The data presented in this study are available on request from the corresponding author.