Page 42 - Read Online
P. 42

Page 4 of 10                                                   Lu et al. Hepatoma Res 2018;4:21  I  http://dx.doi.org/10.20517/2394-5079.2018.44


                Staphylococcus aureus infection             KEGG          3.410    1.3766   52      504
                Serotonergic synapse                        KEGG          -1.854   1.3733   73      3128
                N-glycan antennae elongation in the medial/trans-Golgi   REACTOME  2.122   1.3729   14  396
                Integrins in angiogenesis                   NCI/NATURE    -1.929   1.3622   74      2110
                Tandem pore domain potassium channels       REACTOME      2.299    1.3589   4       206
                Fatty acid elongation in mitochondria       REACTOME      2.226    1.3579   12      170
                IL5-mediated signaling events               NCI/NATURE    2.080    1.3568   12      304
                Antigen processing and presentation         KEGG          3.506    1.3397   65      400
                Asthma                                      KEGG          3.713    1.3246   31      200
                Neurotransmitter release cycle              REACTOME      1.846    1.3233   9       326
                Classical complement pathway                BIOCARTA      2.682    1.3164   12      239
                Antigen processing and presentation         BIOCARTA      2.938    1.2857   9       52
                Interferon gamma signaling                  REACTOME      3.080    1.0558   61      2598
                Antigen processing-cross presentation       REACTOME      2.187    1.0371   59      1962


               the unit of analysis and accounts for interactions among features within the network. In this analysis “antigen
               processing and presentation” was identified as having significant differences in variability in a population
               of Korean HBV associate HCC cases and controls. Consistent with the results of the enrichment analysis,
               re-analysis of this dataset with an extended set of 1200 pathways again identified “antigen processing and
               presentation”, but also “interferon gamma signaling”, “TCR signaling”, and “T cell receptor signaling pathway”
               [Table 1] suggesting that immune response may be a key driver of HCC susceptibility.



               THE ROLE OF ANTIGEN PROCESSING AND PRESENTATION IN HCC
               To assess what might be the key factors within “antigen processing and presentation”, we performed analysis
               utilizing a modified version of PoDA using the Korean HCC dataset. In this analysis, all 400 of the SNPs
               genotyped in the data set for the 65 genes in the pathway were contrasted in the cases and controls. After
               assessing significance of the odds ratio for the entire set of SNPs, each individual SNP was removed one at
               a time from the dataset and the significance was re-assessed. The SNP which least affected the significance
               of the odds ratio was then removed and the process was repeated. SNPs were progressively removed in this
               “stepdown” procedure until the significance of the odds ratio was no longer improved. Interestingly, it was
               observed that initial removal of SNPs substantially improved significance of the difference between cases and
               controls. When stepdown was completed, a total of 49 SNPs in 26 genes were observed [Table 2].


               While the genes identified included key genes seen in the GWAS catalog, specifically members of HLA class
               II, other genes associated with antigen processing were also observed [Figure 1]. The design of Genome-
               wide association studies does not permit the specific etiologic effects of the variation. By design, the variation
               used in the studies is not chosen for function, but instead the ability to test differences between populations.
               The high linkage disequilibrium observed between variations in humans further complicates the capacity to
               interpret the molecular mechanisms of action.


               Nevertheless, this study identifies variation of genes of potential significance in etiology. Of particular interest
               are the proteasome (HSPA2, HSPA4, HSPA5 HSP90AB1), endoplasmic reticulum TAP1, TAP2, CANX), and
               exosome (LGMN) genes associated with the processing of antigens so that they may be presented by HLA
               loci. The pathway also identifies genes on the surface of immune cells - NK cells (KIR2DL3, KIR2DL4, and
               KIR2DL5) and CD4 T cells (CD4) that may compromise immune surveillance and regulation.


                                                                                                       [18]
               It is possible to examine the intra-pathway associations of the variants. Using the analytic tool PLINK ,
                                            2
               one can estimate the association (r ) between loci in cases and controls [Table 3]. As expected by the PoDA
               analysis, variants within the pathways are associated with one another. Both variants within loci and between
               loci are observed to be associated. Interestingly, the magnitude of associations differs between cases and
   37   38   39   40   41   42   43   44   45   46   47