High resolution HLA-DRB1 analysis and shared molecular amino acid signature of DRβ1 molecules in Occult hepatitis B infection

To investigate the association of human leukocyte antigen (HLA)-DRB1 alleles and the variations of polymorphic amino acid changes in DRβ1 chain in Shaanxi Han population with Occult hepatitis B infection (OBI). High-resolution HLA-DRB1 genotyping was performed in 107 OBI carriers and 280 normal controls. Sequence information was used to assign which amino acids were encoded at all polymorphic positions. Three-dimensional modeling was performed to explore the effect of the key residues on the HLA-DRB1 molecule. Strong susceptible association for allele DRB1*07:01 was observed in OBI carriers. The amino acid variation at HLA-DRβ1 molecule revealed susceptible associations for residues Gln4β, Val57β(P9), Ser60β(P9) and Val78β(P4), the amino acids Arg4β, Asp57β(P9), Tyr60β(P9) and Tyr78β(P4) showed protective associations. Alleles DRB1*07:01 showed strong susceptible associations in OBI carriers. The amino acid variations in DRβ molecules revealed significant molecular markers for susceptibility and protection from OBI in Shaanxi Han population.


Introduction
Occult hepatitis B virus infection (OBI) is a special type of HBV infection. When the infection is excluded during the window period, the HBsAg is negative in the serum, while HBV DNA is positive (the concentration of HBV DNA in the serum, or liver tissue, lymphatic system cells < 200 IU/ml) [1,2]. The prevalence of OBI among blood donors in different parts is quite variable depending on the HBV endemicity. The OBI prevalence rates of 1:7517, 1:8209 and 1:9819 were observed among blood donors in China, South African and Europeans [3][4][5]. The prevalence of OBI is higher in chronic liver disease and may be as high as 40 to 75% in HBsAg-negative hepatocellular carcinoma [6]. OBI may impact in several different clinical contexts, including the possible transmission of the infection [7,8], the risk of reactivation [9], the contribution to liver disease progression and to the development of hepatocellular carcinoma [10].
The mechanisms of OBI have not been completely understood. Studies have showed that HBV replication and transcription have been strongly suppressed by the host immune response which consists mainly of anti-viral cytotoxic T cells [11,12]. Furthermore, OBI carriers with anti-HBc show a typical T cells memory response [13].

Open Access
*Correspondence: qijun0802@163.com Human leukocyte antigen (HLA) plays an important role in anti-viral immunity [14]. HLA class II molecules recognize and bind to viral antigen peptides, then present them to CD4 + T cells. As a result, B cells produce antiviral antibodies. In addition, the HLA-DRB1 molecule is a heterodimer, composed of an α and a β chain; the heterodimer forms a peptide-binding groove consisting of nine peptide residue positions (P1-P9), of which positions P1, P4, P6, P7 and P9 constitute peptide anchoring pockets, for the various candidate antigenic peptides. The peptide-binding groove of HLA-DRB1 molecular accommodates 13-18 amino acid peptides. HLA molecular differing only at 1-2 amino acids showed different peptide binding preference. Lacking of evaluating the relationship of HLA-DRB1 alleles with OBI in Chinese Han population. The association between HLA-DRB1 amino acid signatures of the peptide-binding pockets and OBI infections has not been characterized previously.
In this study, we analyzed the risk HLA-DRB1 alleles in a case-control study in the Shaanxi Han populations associated with OBI. We also reported differences in the amino acid signatures of peptide-binding pockets of the HLA-DRB1 molecules in OBI carriers as compared with controls. Determination of the structural of DRB1 molecules associated with OBI help in identifying the disease mechanism and deepening understanding of OBI etiology.

Materials and methods
Subjects 178,941 healthy blood donors were collected from the Shaanxi Blood Center during the period from 1 January 2013 to 31 December 2013. All blood donors have no underlying diseases. Samples that had confirmatory HBsAg (−)/HBV DNA (+) were detected anti-HBs, HBeAg, anti-HBe and anti-HBc. Samples with HBsAg (−)/HBV DNA (+) and serologic markers (+) were collected. Samples with HBsAg (−)/HBV DNA (+) and serologic markers (−) were tested again after six months and only markers that were negative were collected.
A total of 107 blood samples which showed negative for anti-hepatitis C virus (HCV), anti-human immunodeficiency virus (HIV) and HBsAg while positive for HBV NDA were collected from the OBI infection volunteers during January to December in 2013. 280 healthy control subjects were involved from voluntary blood donors whose bloods were negative for anti-HCV and anti-HIV, HBsAg and HBV DNA.

HLA-DRB1 data
High-resolution HLA-DRB1 genotypes were available from a previous study [15]. High-resolution DRB1 sequencing was completed for all study participants using an AlleleAEQR HLA-DRB1 reagent kit and protocol (One Lambda, INC., Canoga Park, CA, USA) and an ABI 3730xl automated sequencer (Applied Biosystems, Foster City, CA). Four-digit DRB1 genotypes were assigned using Assign software (ThermoFisher, Waltham, USA).

Sequence's alignment and 3D protein structure modeling of HLA-DRB1 molecules
Sequences of DRB1 alleles were obtained from the IMGT/HLA database release 3.39.0, January 2020 release (http:// raw. githu buser conte nt. com/ ANHIG/ IMGTH LA/ Latest/ align ments/ DRB1_ pro. txt). Each individual was assigned two amino acids for each polymorphic residue. HLA-DRB1 proteins of known structure suitable as modeling templates were identified in the Protein Data Bank (PDB; http:// www. rcsb. org/ pdb/) and evaluated for structural quality. The binding pockets were identified based on the annotation from Stern et al. [16]. The structure model of OBI-susceptible HLA-DRB1 molecules were constructed based on the x-ray diffraction structure of HLA-DR1 in complex with an endogenous peptide [PDB ID: 1AQD].

Statistical analysis
The differences in allele frequencies of HLA-DRB1 between OBI carriers and controls were compared by chi-squared statistics. Odds ratio (OR) and 95% confidence intervals (95%CI) were calculated for each allele. All amino acid polymorphic positions derived from exon 2 of HLA-DRB1 were tested for their involvement as risk factors for the development of OBI. The association of amino acid residue carrier frequencies were analyzed using Fisher' s exact test with 2 × 2 contingency table. Adjustment for multiple comparisons were conducted with Bonferroni method, corrected P (Pc) values were calculated by multiplying the P value by the number of tested amino acid positions [17]. All the statistical analyses were carried out using the SPSS (SPSS Inc, Chicago, IL) and Graphad Prism (Graphad software Inc, San Diego, CA) packages.

Clinical characteristics of subjects
A total of 387 blood donor (107 OBI carriers and 280 normal controls) were enrolled in this study. General information about the subject was provided in Table 1. The sex, age, address and type of donor distributions showed no statistical difference between the OBI carriers and normal controls. All participants were Chinese Han living in the northwest of China. Among 107 cases of OBI, 104 (97.20%) were positive for either anti-HBs, anti-HBc, anti-HBe, or both. 48 (44.86%) of the OBI cases were positive for both anti-HBs and anti-HBc, while 6 (5.61%) were positive for only anti-HBs, and 37 (34.58%) were positive for only anti-HBc. 3 (2.80%) of the 107 cases of OBI were negative for all serological markers.

Susceptible of HLA-DRB1 amino acid residues associations with OBI
Unequivocal sequences were obtained in all 107 OBI carriers and in 280 controls. Sequencing analysis of the polymorphic exon 2 of the DRB1 gene showed 31 polymorphic amino acid positions. HLA-DRB1 alleles with incidence of over 2% in controls or patients with OBI were shown in Fig. 2.
Significant differences in the frequencies of amino acid between OBI carriers and controls were observed at 13 of the 31 polymorphic amino acids, at positions 4, 11, 13, 14, 25, 28, 30, 37, 57, 60, 73, 74 and 78. The results of amino acid variation for each position are provided in Table 4

Specific amino acid residues of HLA-DRB1 mediate risk and protection in OBI
Four amino acids at the same four positions showed the strongest protective effects (Table 5 Tables 4 and 5. The amino acid residue on position 57 participates in the formation of a α-helix in the ninth peptide-binding pocket of HLA-DRβ1 chain [18]. Among the four possible amino acid residues, Alanine(A), Aspartic acid(D), Serine(S), Valine(V), at this position, Valine(V) showed a risk effect for OBI (P = 0.0002, OR = 1.866, 95% CI 1.357-2.566). Aspartic acid(D) showed a protective effect (P = 0.0359, OR = 0.7027, 95% CI 0.5102-0.9678), the other two residues had no significant effect.

Analysis of genotype combination of HLA-DRB1 epitope
The analysis for homozygous and heterozygous genotype combination of HLA-DRB1 residues was presented (

Discussion
The present study examined HLA-DRB1 alleles and amino acid residues in OBI carriers and control group in Shaanxi Han population.     [26], GMLPVCPLIPGSTTTNTG [27], TSLNFLGGSPVC [27], TTNTGPCKTCTTPAQG [27], WASVRFSWLSLL-VPF [28] have been defined and identified in the Immune epitope database and analysis resource (IEDB) online server (http:// www. iedb. org). The expression and type of HBV related epitopes in OBI have not been reported. A detailed understanding of HBV related epitope in OBI may potentially reveal the pathogenesis of OBI. Determination of the structural of the DRB1 molecules associated with OBI could help in identifying the disease mechanism. Many studies reported that rheumatoid arthritis (RA) susceptibility is linked to amino acid residues within DRB1 peptide-binding groove. DRB1 alleles with a conserved sequence at 70-74 amino acid residues (the shared epitope) in the third hyper-variable region of DRB1 molecular are involved in RA etiology and pathogenesis [29][30][31]. Menconi et al. have identified DRβ-Tyr-26, DRβ-Leu-67, DRβ-Lys-71, and DRβ-Arg-74 are strongly associated with type 1A diabetes and autoimmune thyroid disease [32]. Shared amino acids in the peptide-binding pocket have been demonstrated in nasopharyngeal carcinoma [33], primary biliary cholangitis [34], Parkinson's disease [35]. Sakai A et al. have reported that HLA-DRβ1 position 26 and HLA-DPβ1 position 84-87 are independently associated with anti-HBs production against HBV vaccine [36].
To decipher the precise structural elements mediating the HLA-DRB1 associations with OBI, we examined the contribution of individual amino acids of the HLA-DRB1 molecule to disease risk. The amino-acid signatures for Three-dimensional ribbon models of the HLA proteins associated with OBI. The protein structures of HLA-DRB1 are based on Protein Data Bank (PDB) entries 1AQD. Antigen peptides are shown in blue which was previously shown to bind to HLA-DRB1*01:01. DR α chain is showed by green cyan, while β chain is showed by orange. The four critical residues Gln4β, Val57β(P9), Ser60β(P9) and Val78β(P4) are shown by yellow arrows. The residues in the binding groove of the HLA were selected by proximity (2.45 Å) to the peptide the peptide-binding pockets of DRβ-chain reveal specific molecular change for predisposition and protection from the disease. The data demonstrate that susceptible associations related to residues Gln 4β , Val 57β (P9), Ser 60β (P9) and Val 78β (P4), whereas the amino acids Arg 4β , Asp 57β (P9), Tyr 60β (P9) and Tyr 78β (P4) showed protective associations.
The immune synapse (IS) is initiated by T cell receptor (TCR) interaction with peptide-major histocompatibility complex (pMHC) to form a specificity complex. MHCpeptide monomers bind but do not activate T cells. A dimer of MHC-peptide complexes is necessary and sufficient for initiation of T cell activation. HLA-DR crystals form homodimers of ab heterodimers, a structure that in vitro leads to improved sensitivity on the part of CD4 + T cells [16,39]. The amino acid residue of DRB1 at position 4 is located on the outside of the antigen-binding groove. Arg 4β (a positively charged hydrophilic amino acid) is substituted for Gln 4β (common neutral amino acid). Gln 4β alleles might form an IS faster, precisely because of the location of this residue, near the aminoterminal domain of the β-chain, free to associate its neighboring chains, while not so easy to do in the case of Arg 4β molecules. In the next step research work, we will investigate the properties of HLA-DR7 as well as other β4Arg + alleles in their homozygous state (expression levels, molecular stability). The amino acid residues at position 57 and 60 influences the binding of antigen side chains associated with the P9 pocket. Previous studies have been suggested that the antigenic peptides bound to HLA-DR is influenced by interactions between pocket 6 and pocket 9 [40]. β57 is within pocket 9, with Asp/non-Asp occupancy at this position determining the type of anchor residue at this pocket (mostly small aliphatic vs acidic) [16,41]. The amino acid residue at 60 changes from Tyr to Ser, determining the type of anchor residue at this pocket (aromatic vs small) [42]. The amino acid residue at position 78 located in pocket 4 of the DR peptide-binding cleft [43], The P4 pocket of peptide binding groove can influence and/or modifying the peptide-binding repertoire of DR molecules. The exact mechanism by which Val 78β confers susceptibility to OBI is unknown. The change from the common neutral amino acid (Tyr) to hydrophobic amino acid (Val), the amino acid change possibility modifies the structure of the peptide binding cleft and altering the peptide binding properties.
Inferring polymorphic HLA-DRB1 amino acid positions from the alleles allows for interrogation of specific residues within the DRB1 binding pocket for association with OBI. We identified HLA-DRB1 amino acid positions 4, 57, 60 and 78 associated OBI, which have been associated with autoimmune disease: position 57 with Vogt-Kusanagi-Harada's disease [44], positions 60 with multiple sclerosis [45,46], position 78 with primary biliary cirrhosis [47], indicating the auto immunogenic potential of these HLA-DRB1 amino acid positions susceptible to OBI. The inheritance of two copies (homozygous genotype) and one copy (heterozygous genotype) of amino acid conferred a higher probability of risk and protect for OBI. The amino acids Arg 4β , Asp 57β , Tyr 60β and Tyr 78β show protective association with OBI, while the frequency of amino acid genotype Arg 4β /Arg 4β , Asp 57β /Asp 57β , Tyr 60β /Tyr 60β , Tyr 78β / Tyr 78β in control group were higher than OBI carriers, also reveal protective. The frequency of heterozygous genotype Gln 4β /Arg 4β , Asp 57β /Val 57β , Val 78β /Tyr 78β were higher in OBI carriers. Based on this result, we infer that certain types of amino acid affect the binding to antigenic peptides, resulting in the correlation between certain types of HLA-DRB1 and OBI. In doing so, we are able to identify the specific amino acid and elucidate structural features mediating both risk and protection in OBI, may the integrated influence of synergic effects, counteraction or offsetting effects. The susceptibility effects observed for Gln 4β , Val 57β , Ser 60β , Val 78β can be attributed to association of these variants with allele of HLA-DRB1*07:01, aligning with our initial observation of susceptibility mediated by DRB1*07:01. The protective effects observed for Arg 4β , Asp 57β /Ser 57β , Tyr 60β and Tyr 78β can be attributed to association of these variants with allele of HLA-DRB1*01:01, HLA-DRB1*08:03, and HLA-DRB1*15:01, aligning with our initial observation of protective mediated by these alleles. All 4 residues of DRB1 are located either under (Arg/Gln 4β under the P9 pocket) or within the peptide-binding groove (β57, β60, β78), suggesting that specific patterns of these residues recognize specific peptides from viral antigens as the functional mechanism.
There are three primary limitations to this study. First, immune evasion caused by mutations has not yet been done in this study. Second, the 107 OBI samples were collected from 178,941 blood donors, 103 of them were serologically positive and 4 were serologically negative OBI. The differences in HLA-DRB1 alleles and amino acids between serologically positive and negative OBI have not been compared. And finally, the HBV infection patients have not been studied in this study. Next, our study group will carry out this part of research.
In summary, we find a clear role for HLA-DRB1 in predisposition and protection in OBI carriers. The association of specific combinations of amino acids that participate in critical peptide-binding pockets of the DRB1 molecule implies antigen presentation as a key factor in HBV and other viral infectious diseases.

Conclusion
Our key finding are as follows: (1) HLA-DRB1*07:01 significantly associates with OBI, which agrees with several studies about chronic HBV persistence. (2) Determined and analyzed the amino acid sequence encoded by the DRB1 allele in relation to OBI. We can identify the specific amino acids and elucidate structural features mediating risk and protection to OBI. The amino acid variation at HLA-DRβ reveal susceptible associations for residues Gln 4β , Val 57β (P9), Ser 60β (P9) and Val 78β (P4), whereas the amino acids Arg 4β , Asp 57β (P9), Tyr 60β (P9) and Tyr 78β (P4) show protective associations.