Indonesians Human Leukocyte Antigen (HLA) Distributions and Correlations with Global Diseases

Design of multi-epitope-based therapeutic vaccine candidates from HBc and HBx proteins of hepatitis B virus using reverse vaccinology and immunoinformatics approaches

The major problem in cases of chronic hepatitis B (CHB) is the failure of the patient's immune response to eliminate the covalently closed circular DNA (cccDNA) minichromosome of hepatitis B virus (HBV). Epigenetic regulation involving the HBV core protein (HBc) and HBV X protein (HBx) influences the transcription and stability of the cccDNA minichromosome. The HBc and/or HBx-based therapeutic vaccines that have been developed cannot accommodate differences between HBV genotypes. This research aims to design a therapeutic vaccine candidate based on the multi-epitope of HBc and HBx using reverse vaccinology (RV) and immunoinformatics approach. HBc and HBx sequences from 10 HBV genotypes were obtained from the NCBI Entrez Protein database. Epitopes were predicted from consensus sequences, which consisted of 13,610 HBc sequences and 12,333 HBx sequences. The study identified four cytotoxic T lymphocyte epitopes, two helper T lymphocyte epitopes, and five linear B lymphocyte that met the inclusion criteria. The vaccine candidate designed using cholera toxin subunit B and pan HLA DR-binding epitope adjuvants was predicted to be safe, antigenic, stable, and has a global population coverage of 99.43%. Molecular docking and molecular dynamics simulations demonstrated that the vaccine candidate could stably bind to B cell receptor, cytotoxic T cell receptor, and TLR4 for 100 ns. Immune response simulation indicated that it can induce antibody production and the proliferation of B and T cells. It can be concluded that RV and immunoinformatics successfully facilitated the design of a multi-epitope therapeutic vaccine candidate for CHB.

Graves' disease and inflammatory bowel disease: A bidirectional Mendelian randomization

CONTEXT

Both Graves' disease (GD) and inflammatory bowel disease (IBD) are common autoimmune diseases that severely damage patients' quality of life. Previous epidemiological studies have suggested associations between GD and IBD. However, whether a causal relationship exists between these two diseases remains unknown.

OBJECTIVE

To infer a causal relationship between GD and IBD using bidirectional two-sample Mendelian randomization(MR).

METHODS

We performed bidirectional two-sample MR to infer a causal relationship between GD and IBD using GWAS summary data obtained from Biobank Japan (BBJ) and the International Inflammatory Bowel Disease Genetic Consortium (IIBDGC). Several methods (random-effect inverse variance weighted, weighted median, MR‒Egger regression, and MR-PRESSO) were used to ensure the robustness of the causal effect. Heterogeneity was measured based on Cochran's Q value. Horizontal pleiotropy was evaluated by MR‒Egger regression and leave-one-out analysis.

RESULTS

Genetically predicted IBD may increase the risk of GD by 24% (OR 1.24, 95% CI 1.01-1.52, p = 0.041). Crohn's disease (CD) may increase the risk of GD, whereas ulcerative colitis (UC) may prevent patients from developing GD. Conversely, genetically predicted GD may slightly increase the risk of CD, although evidence indicating that the presence of GD increased the risk of UC or IBD was lacking. Outlier-corrected results were consistent with raw causal estimates.

CONCLUSIONS

Our study revealed a potentially higher comorbidity rate for GD and CD. However, UC might represent a protective factor for GD. The underlying mechanism and potential common pathways await discovery.

Immunoinformatics Analysis of Citrullinated Antigen as Potential Multi-peptide Lung Cancer Vaccine Candidates for Indonesian Population

Non-small-cell lung cancer (NSCLC) is the most common lung cancer which has the highest mortality rate in Indonesia. One of the trends in treating cancer is by utilizing peptide vaccines, an immunotherapeutic approach that aims to stimulate the cell-mediated adaptive immune system to recognize cancer-associated peptides. Currently, no peptide vaccines are available in the market for NSCLC treatment. Therefore, this project aims to develop a multi-epitope peptide-based vaccine for NSCLC utilizing citrullinated peptides. Citrullination is a post-translational modification that occurs in cancer cells during autophagy that functions to induce immune responses towards modified self-epitopes such as tumor cells, through activation of PAD enzymes within the APC and target cells. It was found that introducing a common citrullinated neo-antigen peptide such as vimentin and enolase to the immune system could stimulate a higher specific CD4⁺ T cell response against NSCLC. Moreover, carcinoembryonic antigen (CEA), an antigen that is highly expressed in cancer cells, is also added to increase the vaccine’s specificity and to mobilize both CD4⁺ and CD8⁺ T cells. These antigens bind strongly to the MHC Class II alleles such as HLA-DRB1*07:01 and HLA-DRB*11:01, which are predominant alleles in Indonesian populations. Through in silico approach, the peptides generated from CEA, citrullinated vimentin and enolase, were analyzed for their MHC binding strength, immunogenicity, ability to induce IFNγ response, and population coverage. It is expected that the immunodominant antigens presentation is able to induce a potent immune response in NSCLC patients in Indonesia, resulting in tumor eradication.

Construction, expression, and in vitro assembly of virus-like particles of L1 protein of human papillomavirus type 52 in Escherichia coli BL21 DE3

Background

A major discovery in human etiology recognized that cervical cancer is a consequence of an infection caused by some mucosatropic types of human papillomavirus (HPV). Since L1 protein of HPV is able to induce the formation of neutralizing antibodies, it becomes a protein target to develop HPV vaccines. Therefore, this study aims to obtain and analyze the expression of HPV subunit recombinant protein, namely L1 HPV 52 in E. coli BL21 DE3. The raw material used was L1 HPV 52 protein, while the synthetic gene, which is measured at 1473 bp in pD451-MR plasmid, was codon-optimized (ATUM) and successfully integrated into 5643 base pairs (bps) of pETSUMO. Bioinformatic studies were also conducted to analyze B cell epitope, T cell epitope, and immunogenicity prediction for L1HPV52 protein.

Results

The pETSUMO-L1HPV52 construct was successfully obtained in a correct ligation size when it was cut with EcoRI. Digestion by EcoRI revealed a size of 5953 and 1160 bps for both TA cloning petSUMO vector and gene of interest, respectively. Furthermore, the right direction of construct pETSUMO-L1HPV52 was proven by PCR techniques using specific primer pairs then followed by sequencing, which shows 147 base pairs. Characterization of L1 HPV 52 by SDS-PAGE analysis confirms the presence of a protein band at a size of ~55 kDa with 6.12 mg/L of total protein concentration. Observation under by transmission electron microscope demonstrates the formation of VLP-L1 at a size between 30 and 40 nm in assembly buffer under the condition of pH 5.4. Based on bioinformatics studies, we found that there are three B cell epitopes (GFPDTSFYNPET, DYLQMASEPY, KEKFSADLDQFP) and four T cell epitopes (YLQMASEPY, PYGDSLFFF, DSLFFFLRR, MFVRHFFNR). Moreover, an immunogenicity study shows that among all the T cell epitopes, the one that has the highest affinity value is DSLFFFLRR for Indonesian HLAs.

Conclusion

Regarding the achievement on successful formation of L1 HPV52-VLPs, followed by some possibilities found from bioinformatics studies, this study suggests promising results for future development of L1 HPV type 52 vaccine in Indonesia.

HLA-DQA1 and HLA-DQB1 Gene Polymorphism in Indonesian Children with Type I Diabetes Mellitus

BACKGROUND

More than 40 genes influence the progression of type 1 diabetes mellitus (T1DM), including human leukocyte antigen (HLA) alleles. Different HLA genotype patterns result in diverse rates of T1DM development. HLA class II DR, DQ, and DP vary among different populations and ethnicities. Data on HLA polymorphism in T1DM in Indonesia are lacking. Therefore, this study was designed to evaluate the gene polymorphism of HLA-DQA1 and HLA-DQB1 in Indonesian children with T1DM.

PATIENTS AND METHODS

In this study, 31 patients with T1DM and 31 controls were enrolled from April 2020 to April 2021. This study was conducted at Dr. Soetomo Hospital, Indonesia. We evaluated the gene polymorphism of HLA-DQA1 and HLA-DQB1 using polymerase chain reaction-restriction fragment length polymorphism. The primers used were as follows: for HLA-DQA1, DQAS34: 5'-GGTGTAAACTTGTACCAG-3' (forward) and DQAA261: 5'-ATTGGTAGCAGCGGTAGA-3' (reverse); for HLA-DQB1, DQBS43: 5'-TGCTACT- TCACCAA(C/T)GGG-3' (forward) and DQBA249: 5'-GTAGTTGTGTCTGCA (C/T)AC-3' (reverse).

RESULTS

The most common HLA-DQA1 subtype in the T1DM group was 0101/0102 accounting for 67.6%, and 01/03 and 02/03 were found in the T1DM group only. Meanwhile, the most common HLA-DQB1 subtype in the T1DM group was 0301, accounting for 54.8%. Most subjects in this study were Javanese.

CONCLUSION

HLA-DQA1 0101/0102 and HLA-DQB1 0301 were commonly found in Indonesian children with T1DM.

Evaluation of Human Leukocyte Antigen Class I and Class II in End-Stage Renal Disease Occurrence in Indonesian Transplantation Patients

Background

Genetic studies of end-stage renal disease (ESRD), including those of human leukocyte antigen (HLA) genes, have been reported in several populations but have not yet been evaluated in Indonesia. Some studies have reported that these genes had a substantial role in ESRD. This study aims to analyze the association between HLA genes and ESRD within the Indonesian community.

Method

A retrospective study to investigate HLA class I and II alleles to find out the distribution of HLA-A, -B, -C, -DPB1, -DQB1, and -DRB1 in renal transplant recipients and to ascertain their role in susceptibility to ESRD was performed on totally 149 subjects, consisting of 69 ESRD patients and 80 healthy controls. HLA typing was determined using Luminex techniques. The allele and haplotype frequencies were compared between ESRD patients and controls.

Result

High-frequency alleles were HLA-A∗24 (43.6%), B∗15 (38.2%), C∗08 (30.8%), DRB1∗12 (47.3%), DQB1∗03 (50.6%), and DPB1∗13 (22.5%). HLA-A∗24 (p=0.01) and HLA-B∗35 (p=0.02) were associated with a protective effect, with OR 0.537 (95%CI 0.34–0.86) and 0.316 (95%CI 0.11–0.88), respectively. There were some two-locus haplotypes associated with susceptibility to ESRD, such as B∗15-DRB1∗12, B∗13-DRB1∗15, A∗02-B∗15, A∗02-C∗08, and B∗13-DQB1∗05. HLA-A∗02-B∗15-DRB1∗12 and A∗24-B∗13-DRB1∗15 appear to be associated with susceptibility to ESRD.

Conclusion

The allele groups of HLA-A∗24 and HLA-B∗35 are associated with protection from ESRD. Meanwhile, HLA-B∗13-DRB1∗15 and A∗24-B∗13-DRB1∗15 are the most frequent HLAs associated with ESRD in two-locus and three-locus haplotype, respectively.

Association Study between Psoriatic Arthritis and Killer Immunoglobulin-Like Receptor (KIR) Genes: A Meta-Analysis

Background: Psoriatic Arthritis (PsA) is a seronegative spondyloarthropathy frequently associated with psoriasis. Studies have shown different members of the KIR (Killer Immunoglobulin-like Receptor) gene family may act as potential susceptibility factors; however, data have been inconsistent or with a reduced sample size. Therefore, the objective of this investigation was to determine associations between KIR genes and PsA susceptibility a meta-analysis approach. Methods: We performed a systemic search on PubMed, Scopus, and Web of Science to identify association studies linking KIR genes with PsA susceptibility. The search cut-off was May 2019. Odds Ratio (OR), 95% Confidence Intervals (95% CI), and forest plots were obtained for each KIR gene. Publication bias was evaluated by Begg and Egger linear regression tests. Results: Five articles were included in this meta-analysis. The KIR2DL2, 2DS1, 2DS2, and 2DS3 genes were positively associated with susceptibility to PsA (OR = 1.269, p = .003; OR = 1.392, p < .001; OR = 1.279, p = .002; and OR = 1.230, p = .038, respectively). In Caucasians, positive association with susceptibility to PsA were maintained by KIR2DL2, 2DS1, and 2DS2 genes (OR = 1.257, p = .005; OR = 1.535, p = .003; and OR = 1.267, p = .004, respectively). Conclusion: These associations suggest that KIR2DL2, 2DS1, 2DS2, and 2DS3 genes are susceptibility factors for PsA.