Impact Factor 4.716 | CiteScore 4.71
More on impact ›

Original Research ARTICLE

Front. Immunol., 16 April 2018 | https://doi.org/10.3389/fimmu.2018.00777

Impact of Interferon-α Receptor-1 Promoter Polymorphisms on the Transcriptome of the Hepatitis B Virus-Associated Hepatocellular Carcinoma

imageTimokratis Karamitros1,2*, imageGeorge Papatheodoridis3, imageDimitrios Paraskevis4, imageAngelos Hatzakis4, imageJean L. Mbisa5, imageUrania Georgopoulou6, imagePaul Klenerman7 and imageGkikas Magiorkinis1,4*
  • 1Department of Zoology, University of Oxford, Oxford, United Kingdom
  • 2Department of Microbiology, Public Health Laboratories, Hellenic Pasteur Institute, Athens, Greece
  • 3Academic Department of Gastroenterology, Laiko General Hospital, Medical School of National and Kapodistrian University of Athens, Athens, Greece
  • 4Department of Hygiene and Epidemiology and Medical Statistics, Medical School of National and Kapodistrian University of Athens, Athens, Greece
  • 5Virus Reference Department, Public Health England, London, United Kingdom
  • 6Department of Microbiology, Molecular Virology Laboratory, Hellenic Pasteur Institute, Athens, Greece
  • 7Peter Medawar Building for Pathogen Research and Translational Gastroenterology Unit, University of Oxford, Oxford, United Kingdom

Background and aims: Genetic polymorphisms within the promoter of interferon-α receptor type-1 (IFNAR1) have been associated with the susceptibility to and the outcome of chronic hepatitis B virus (HBV) infection. However, the impact of these polymorphisms in the transcriptome of the HBV-associated hepatocellular carcinoma (HCC) remains largely unexplored.

Methods: Using whole-genome and exome sequencing data from The Cancer Genome Atlas project, we characterized three single-nucleotide polymorphisms (SNPs: −568G/C, −408C/T, −3C/T) and one variable number tandem repeat [VNTR: −77(GT)n] within the IFNAR1 promoter sequence in 49 HCC patients. RNAseq data from 10 genotyped HCC samples were grouped according to their −77VNTR or −3SNP genotype to evaluate the impact of these polymorphisms on the differential expression on the HCC transcriptome.

Results: There is a fourfold higher impact of the −77VNTR on the HCC transcriptome compared to the −3SNP (q < 0.1, p < 0.001). The expression of the primary IFNAR1 transcript is not affected by these polymorphisms but a secondary, HCC-specific transcript is expressed only in homozygous −77VNTR ≤8/≤8(GT)n samples (p < 0.05). At the same time, patients carrying at least one −77VNTR >8(GT) allele, presented a strong upregulation of the fibronectin-1 (FN-1) gene, which has been associated with the development of HCC. Gene Ontology and pathway enrichment analysis of the differentially expressed genes revealed a strong disruption of the PI3K–AKT signaling pathway, which can be partially triggered by the extracellular matrix FN-1.

Conclusion: The IFNAR-1 promoter polymorphisms are not involved in the expression levels of the main IFNAR-1 transcript. The −77VNTR has a regulatory role on the expression of a secondary, truncated, HCC-specific transcript, which in turn coincides with disruptions in cancer-associated pathways and in FN-1 expression modifications.

Introduction

Worldwide, more than two billion people have been infected with hepatitis B virus (HBV) and approximately 250 million individuals are chronically infected (1). Infected patients can be inactive chronic HBV carriers (IC) (eAg-negative, eAb-positive, with low levels of HBV DNA and no evidence of liver inflammation) or present with the progressive chronic hepatitis B (CHB) (25). HBV is responsible for >50% of hepatocellular carcinomas (HCCs) worldwide (2, 6, 7). On the other hand, HBeAg-negative ICs have a more benign prognosis with very low risk of cirrhosis or HCC, as indicated by long-term follow-up studies (810).

The genetic profile of the patient plays a substantial role in the clinical outcome of HBV infection (1113). The virus also modulates cellular mechanisms and signal pathways during the course of the infection. For example, the persistent expression of the HBV x antigen (HBxAg) is correlated with the development of fibrosis and cirrhosis during CHB, as it can activate fibronectin-1 (FN-1) gene, through the induction of the nuclear factor kappa B (NF-kappa B or NFkB) (14). FN is an omnipresent extracellular matrix glycoprotein. Plasma FN and cellular FN have distinct properties and roles in the strictly regulated mechanism of tissue repair (15). FN is a very important component of ECM and any dysfunction in the fibrinogenesis mechanisms can lead to the development of fibrotic disease (16).

Interferons (IFN)-α/β are cytokines involved in both innate and adaptive immune responses (17, 18), thus play a pivotal role as cancer defense mechanisms (19). Interferon-induced signal transduction pathways represent a fine-tuned network of interactions, which are triggered upon binding of IFN on its receptor (IFNR) (2022). IFNs are shown to regulate the transcription levels of more than 2,000 genes, which compose the Interferome, a gene-network created by integrating information collected from high-throughput experiments (23). As IFNs are anticancer and antiviral cytokines, it is expected that the genetic profiles of the genes involved in these signal transduction pathways (e.g., IFNR) will have an impact on the susceptibility of the patients to cancer and, when HBV is involved, to CHB and HBV-related HCC.

IFN-α is administrated as part of the first-line therapies for CHB (2, 24). The compatible receptor (IFNAR) consists of two subunits IFNAR-1 and IFNAR-2 (25). A number of polymorphisms in the promoter of IFNAR-1 gene and also in its coding sequence have been described (2629). Those observed in the promoter—three SNPs and one variable number tandem repeat (VNTR) at positions −568, −408, −3 and −77, respectively—are believed to affect the expression of the receptor and are associated with the clinical outcome of HBV infection (27, 30, 31).

We have previously shown that the same polymorphisms are associated with the clinical phase of HBeAg-negative chronic HBV infection in Caucasians. Briefly, patients with genotypes −568GC/CC, −408CT, and −3CT and patients with less than 8(GT) repeats in the −77VNTR were more frequent among inactive carriers (IC) vs. patients with HBeAg-negative CHB (32). Notably, the (GT)-repeats in the −77VNTR were strongly associated with the clinical outcome of the patients in our study; homozygotes carrying both alleles with ≤8(GT)-repeats, were more likely to be ICs, compared to those carrying both alleles with >8(GT) repeats (OR = 7.14, p ≤ 0.001).

In this study, we use whole-genome sequencing (WGS), whole-exome sequencing (WXS), and RNAseq data derived from The Cancer Genome Atlas (TCGA) project in combination with a bioinformatics pipeline, which efficiently and accurately genotypes the −77VNTR. We estimate the distribution of the genetic variations of the four IFNAR-1 promoter polymorphisms and we investigate the impact of −77VNTR and −3SNP polymorphisms on the regulation of HCC transcriptome and interferome. IFNAR-1 is transcribed into three transcripts, ENST00000270139 (which is translated into the receptor), ENST00000493503, and ENST00000442071; performing differential expression analysis in patients grouped according to their −77VNTR genotype, we show that −77VNTR modifies the transcriptional profile of IFNAR-1 gene by controlling the expression of the last one, which comprises only a fibronectin III domain and is translated into a truncated form of the receptor. The expression of this transcript coincides with significantly lower expression of FN-1. Finally, the PI3K–AKT signaling pathway, which is partially triggered by FN-1, was significantly enriched with differentially expressed genes. Our results indicate that this secondary transcript could potentially act as a regulatory element, but its functional role remains largely unknown.

Patients and Methods

Study Design

We examined four previously described polymorphisms within the IFNAR-1 promoter (27, 30, 31), from now on referred to as interferon-α receptor promoter polymorphisms (IFNARPPs). They include three single-nucleotide polymorphisms (SNPs) at positions −568G/C, −408C/T, and −3C/T and a variable number tandem repeat of the binucleotide GT, −77VNTR(GT)n.

We used 49 (15 WGS and 34 WXS) TCGA “liver hepatocellular carcinoma (LIHC)” samples (project ID 10464). We filtered all the available (up to March 2016) HCC samples with positive HBV surface antigen. We performed genotyping of the four IFNARPPs, using the pipeline described below. We compared the frequencies of the revealed genotypes with those of a group of 92 chronically infected HBeAg-negative IC as defined in our previous study (32) as they present very low risk of cirrhosis and HCC (810). Briefly, IC state was defined for patients with persistently normal ALT values under strict follow-up and maximum HBV DNA levels ≤20,000 IU/mL. None of these inactive carriers had cirrhosis (33).

To assess the impact of the −77VNTR and the −3SNP on the expression profiles in HBV-associated hepatocellular carcinoma, we used tumor RNAseq data available for 10 TCGA HCC genotyped samples after classifying them according to their genotype (for −77VNTR: ≤8/≤8 vs. ≤8/>8 or >8/>8(GT) repeats, and for −3SNP: CC vs. CT or TT) (Table 1).

TABLE 1
www.frontiersin.org

Table 1. RNAseq samples analyzed in this study.

Bioinformatics

The .bam files were transformed into paired-end .fastq files using Samtools-bam2fq (34) and were locally mapped against two artificial chromosomes using Bowtie2 (35) in—very sensitive local mode (Figure 1). The artificial reference chromosomes were created after splitting the IFNAR1 promoter sequence (genebank: X60459.1) at the −77VNTR, leaving 3 GT repeats at each hanging end, to avoid non-specific mapping of reads due to the low complexity of the microsatellite repeats. We extracted the mapped reads and used them to de novo assemble the alleles of IFNAR1 promoter for each TCGA patient, using MIRA (36). We further analyzed the de novo assembled contigs using R DNAstrings package to count the GT repeats. We used Samtools mpileup and bcftools (34) to call the variations of the three SNPs from the mapping alignments. The genotyping results were confirmed and quality-controlled by visual inspection of the mapping alignments of 15 random samples, using IGV (37).

FIGURE 1
www.frontiersin.org

Figure 1. Schematic representation of the computational pipeline applied to whole-genome sequencing (WGS) and WXS The Cancer Genome Atlas (TCGA) data (in blue) for the genotyping of the IFNARPPs. IFNAR1 promoter sequence (genebank: X60459.1) was used to for the construction of two pseudo-chromosomes, on which the reads were aligned. The mapped reads were extracted and used for the de novo assembly of the two alleles of each sample. This method allows the assembly of alleles varying substantially in length compared to the reference.

We performed the RNAseq analysis using Kallisto (38) to map the reads against the Human Transcriptome reference (v.GRCh38.rel79) and to calculate the transcripts abundances. We analyzed the impact of the IFNARPPs on the transcriptional landscape of the interferon-associated genes using the “Interferome” database (23) on a subset our whole-transcriptome results. We used Sleuth (39) and R-base functions to interpret and visualize the RNA-seq analysis results. We performed Gene Ontology (GO) and KEGG pathway enrichment analysis using the differentially expressed genes (p < 0.001).

Statistical Analysis

We used the t-test, the Fisher’s exact test, and the z-test to evaluate the association of the demographic and genetic characteristics of the patients with the disease outcome. We used RStudio v0.99.446 for R v3.2.3 programming language for all statistical computing and graphics.

Results

RNAseq Differential Expression Design

We tested for differentially expressed transcripts in 10 HCC RNAseq samples after defining their −77VNTR and −3SNP genotypes. We grouped and compared them according to their −77VNTR genotype [samples A–E: >8/≤8 or >8/>8(GT)n vs. F–J: ≤8/≤8(GT)n] and according to their −3SNP genotype (samples A, B, C, F, G: “CT” vs. D, E, H, I, J: “CC”) (Table 1). The quality-control assessment of the RNAseq analysis is summarized in Figure S1 in Supplementary Material.

The −77 VNTR and −3 SNP Polymorphisms and the IFNAR1 Transcription Profile

IFNAR1 gene generates three different transcripts, but only one of them is translated into the functional receptor protein subunit (transcript “001”: ENST00000270139). Of the two remaining transcripts, one is not translated (transcript “002”: ENST00000493503) and the other (transcript “003”: ENST00000442071) produces a truncated (136 aa long) isoform of the receptor, incorporating only a fibronectin type III domain. We found that −77 VNTR and −3 SNP polymorphisms do not significantly affect the expression levels of the primary transcript “001” and the secondary transcript “002.”

Interestingly, the expression of the secondary transcript “003” was detectable only in homozygous patients carrying both alleles with ≤8(GT) repeats, while it was absent in patients carrying at least one allele with >8(GT) repeats (2.47 vs. 0.01 TPM, respectively) (Figure 2). We further tested 10 TCGA normal tissue (liver) RNAseq samples for the expression of this transcript. We found them all negative except from one, which presented only basal expression (0.02 TPM). Thus, transcript “003” was HCC-specific.

FIGURE 2
www.frontiersin.org

Figure 2. Analysis of IFNAR1 transcripts: (A) structure of IFNAR-1 gene. Alternative splicing and transcription initiation sites produce three distinct transcripts that differ in length. Only transcripts 001 and 003 are translated into proteins, which share only the C-terminal part and a fibronectin type III domain. (B): RNAseq-derived expression levels of the three transcripts from 10 hepatocellular carcinoma (HCC) samples grouped based on their −77VNTR(GT)n genotype. The levels of transcripts 001 (responsible for the production of the receptor) and 002 do not differ significantly between the genotypes. Transcript 003 is expressed only in samples with <8/<8 (GT)n repeats in the −77VNTR polymorphism. The genotypes of HCC RNA samples (A–J) are described in Table 1.

Impact of −77 VNTR and −3 SNP on the HCC Transcriptome and Interferome

We classified 10 HCC TCGA samples according to their genotype, ≤8/≤8 vs. ≤8/>8 or >8/>8(GT) repeats for the −77VNTR, and CC vs. CT or TT for the −3SNP. There were significant changes in the transcription profiles between the groups tested. In detail, for the −77VNTR grouping, there were 246 differentially expressed genes (p < 0.001), while for the −3SNP grouping, only 57 genes showed significant changes in their expression levels (p < 0.001) (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Impact of (A) −77 VNTR (GT)n and (B) −3SNP C/T on the hepatocellular carcinoma transcriptome. Transcripts with statistically significant change in their expression are in red. Positive and negative normalized beta values correspond to up- or downregulation, respectively. There are approximately four times more significantly differentially expressed transcripts in the −77 VNTR (GT)n compared to the −3SNP C/T genotype group.

Focusing on the interferon-related genes, we created a subset of the transcriptomics results according to transcript IDs found in the Interferome database (23). The −77VNTR polymorphism had a fourfold higher impact on the Interferome compared to the −3SNP, as 46 and 11 transcripts, respectively, are either up- or downregulated (p < 0.001). Notably, among the most significantly differentially expressed genes was FN-1, which presented a 4.49-fold lower expression in ≤8/≤8(GT)n homozygous patients (Figure 4; Table S2 in Supplementary Material).

FIGURE 4
www.frontiersin.org

Figure 4. Impact of (A) −77 VNTR (GT)n and (B) −3SNP C/T on the hepatocellular carcinoma interferome. Only statistically significant differences in the expression levels of the transcripts are presented. Affected transcripts, the respective genes, and the transcripts’ product status are on the left.

Impact of −77 VNTR on Signaling Pathways

Performing GO and KEGG pathway enrichment analysis based on the −77VNTR grouping we found six significantly enriched pathways (Table 2). PI3K–AKT signaling pathway was ranked first, with the higher number and proportion of differentially expressed genes involved (14 genes, p < 0.05). The super-family “pathways in cancer” was also found significantly enriched with 15 differentially expressed genes, but since PI3K–AKT pathway was the main contributor, the family was excluded from Table 2. Intriguingly, PI3K–AKT signaling pathway can be triggered by FN-1 (Figure 5).

TABLE 2
www.frontiersin.org

Table 2. Pathway enrichment analysis.

FIGURE 5
www.frontiersin.org

Figure 5. The PI3K–AKT signaling pathway is significantly enriched with differentially expressed genes (magenta). Insulin, ErbB, mTOR, and VEGF signaling sub-pathways triggered by PDK1 have been removed for simplicity.

IFNARPPs and HBV-Associated Hepatocellular Carcinoma

We genotyped the four IFNARPPs in 49 HBV-associated HCC samples, 15 WGS and 34 WXS. We were able to confirm the existence of the −568G allele [previously associated with more pathogenic disease states (26, 32)] in 36 out of 49 HCC samples, as the read coverage was occasionally reduced over the −568SNP in the WXS dataset. We compared their distributions of the polymorphisms with those described previously for 92 ICs (32). The demographic characteristics of the 49 HCC and the 92 IC patients are summarized in Table S1 in Supplementary Material. There was no statistically significant difference in frequency of the −568SNP genotypes between the IC and HCC groups (Table 3). Given the previously described linkage between alleles −408C/T and −3C/T (−408C to −3C and −408T to −3T) (27) these were analyzed together. There was no statistically significant difference in the prevalence of the −408/−3 SNP polymorphisms in the IC and the HCC groups (Table 3). The −408/−3 TT genotype was not identified in any of the TCGA samples tested. Notably, the number of GT repeats at the −77VNTR of IFNAR-1 promoter was associated with the disease status. Heterozygous ≤8/>8 (GT)n patients are less likely to be IC compared to homozygotes carrying both alleles with ≤8 GT repeats (OR = 0.41, 95% CI: 0.19, 0.90). The same trend was observed when >8/>8 and ≤8/>8 (GT)n genotypes were grouped together (OR = 0.40, 95% CI: 0.19, 0.84) (Table 3; Figure S2 in Supplementary Material).

TABLE 3
www.frontiersin.org

Table 3. Genotypes of patient samples in relation to clinical status.

Discussion

Implementing a fully automated computational pipeline, we characterized the genetic profile of the four polymorphisms located in the IFNAR1 promoter region in 49 HBV-associated HCC samples derived from the TCGA project.

Whole-genome sequencing data have enormous size and they are usually available as ready-to-use alignments (.bam files), which are binary and compressed. Variation-calling from bam files can be performed with several routine pipelines, each characterized by different biases toward specific types of SNP and in/del genotyping errors (40). For example, reads that differ significantly in length can be ignored and remain unmapped. This is highly likely in VNTRs with heterozygous genotypes in loci where the two alleles differ substantially in length compared to the reference sequence (e.g., −77VNTR(GT)5/15), which in combination with low read coverage and read length (quite common in WGS data) can lead to inaccurate variant calling. Our pipeline solves this problem, by using an artificial reference of two pseudo-chromosomes, derived from the original sequence, segregated exactly on the VNTR locus, leaving only 3(GT)-repeats at each segregation edge. By using these two pseudo-chromosomes, we were able to specifically mine all the reads corresponding to the IFNAR1 promoter region—originally mapped or unmapped—that in turn precisely reconstructed the two alleles of each sample. Our method provides biologically meaningful results as only one (homozygous) or two (heterozygous) contigs were generated for all the samples tested. There were not any non-specific contigs generated, indicating that the mining of the reads was highly specific.

Since only 15 TCGA WGS samples met our selection criteria, we used WXS data as well. Although WXS libraries are optimized for mRNA sequencing, non-coding upstream and downstream sequences are randomly captured during the target enrichment process (41, 42). Thus, we were able to accurately genotype the proximate to the first IFNAR1 exon polymorphisms, with the exemption of −568SNP, which could only be genotyped in 36 out of the 49 HCC TCGA samples. −568SNP is the most distant one relatively to the start of the transcription, and this resulted in reduced read coverage in 11 WXS samples. The limited number of observations in the HCC group did not allow for statistically significant results with regard to −568 SNP and disease status, although a relative difference in the genotypic distributions was observed (Figure S2 in Supplementary Material). Further studies are needed to confirm the significance of this SNP, as it has been previously associated with the development of CHB and the spontaneous recovery after acute HBV infection (27, 32).

Our results suggest that the number of GT repeats in the −77VNTR is associated with the state of the disease, as alleles with more than 8(GT)n are more frequent among HCC samples. (Table 3; Figure S2 in Supplementary Material). The finding that more (GT) repeats are associated with a more severe outcome of the disease is in concordance with our previous observations where alleles carrying ≤8(GT) repeats in the −77VNTR, were associated with the IC phase, whereas >8(GT) repeats were more frequent in HBeAg-negative CHB patients (32). Zhou et al. studying CHB patients and spontaneously resolved cases after acute HBV infection showed that alleles carrying ≤9(GT) repeats were more frequent in CHB patients, while >9(GT) repeats were associated with spontaneous clearance after acute HBV infection (27). Moreover, this result was linked with the rest of the polymorphisms examined in haplotypes, [≤9(GT)n, −568C, −408/−3T] for the CHB patients and >9(GT)n, −568G, −408C, −3C in cases of spontaneous clearance after acute HBV infection. In our previous study, alleles −568C and −408/−3T were associated with the IC phase, whereas −568G and −408/−3C were more frequent in CHB patients. In this study, although genotypes −568CG and −408/−3CT show a trend of higher frequency in the HCC samples, compared to genotypes −568CC and −408/−3CC (Figure S2 in Supplementary Material), they were not associated with disease, as the limited availability of HCC samples and the lack of coverage over −568SNP in WXS samples led to lack of statistical power for these comparisons. Notably, genotype −408/−3TT was totally absent in the HCC dataset. Studies in larger datasets may discern whether alleles −568G and −408/−3T and the homozygous genotypes GG and TT are linked to HCC.

The effective expression of IFNAR-1 is essential for the linkage with IFN and the triggering of the downstream signal pathways (43). Even minor modifications in the receptor structure can impair its normal function and alter its antiviral and antiproliferative properties (44, 45). Moreover, the expression levels of both the receptor subunits IFNAR1 and IFNAR2 are associated with the INF-β treatment outcome in multiple sclerosis patients (46). Hepatitis C patients with 5GT repeats in the −77VNTR (genotypes GT5/5/GT5/14) have been reported to be better responders to IFN-based therapy (47). On the other hand, experiments based on luciferase reporter plasmids have suggested that the promoter of the IFNAR-1 is not affected by the −568SNP and the −77VNTR (26, 27), while in another study its activity is significantly affected by −3SNP (48). In detail, the promoter activity was reduced for −3T plasmids through reduced binding affinity to HMGB1, a factor that was suggested to bind to the −3 element to regulate the transcription levels of IFNAR1. Similarly, reduced levels of IFNAR1 expression in HBV patients with C > T substitution at the −3 position of the IFNAR1 promoter was reported in the same study (48). In another study, an important role of the −3CT SNP was evident regarding the regulation of the transcription factor High Mobility Group B protein 1 (HMGB-1); the C > T transition was shown to reduce the binding affinity of HMGB-1 to the IFNAR-1 promoter sequence, thus lead to reduced expression of IFNAR-1 (48).

In this study, we conducted a thorough RNAseq differential expression analysis on HCC samples grouped according to their −77VNTR and −3SNP genotype. In agreement with our previous findings, where a significant role of the −77VNTR genotypes had arisen with respect to the clinical course of HBV infection, here we report a more significant involvement of −77VNTR in the modification of the HCC transcriptome and interferome, as more genes were found to be significantly differentially expressed compared to those affected by the −3SNP. We also shed light on the controversial question about the impact of these promoter polymorphisms on the expression of IFNAR-1 gene. The expression levels of the receptor are associated with the clinical outcome of chronically HBV-infected patients (49). We conclude that the expression of the major IFNAR-1 transcript, which is responsible for the production of the functional receptor, is not significantly affected by these polymorphisms (Figure 2).

At the same time, a secondary IFNAR1 transcript, which is translated into a truncated form of the receptor, is only produced at detectable levels in samples with ≤8/≤8 (GT)n genotype in the −77VNTR, independently from the −3SNP genotype. This transcript is not detectable in normal liver tissue samples, thus appears to be HCC specific. Secondary IFNAR1 transcripts have been reported in tumor cell lines (50) and here, for the first time, we report that −77VNTR genotype is crucial for their production, which apparently leads to the transcriptional remodeling of IFNAR1 gene.

Grouping the samples according to their −77VNTR genotype revealed that one of the most differentially expressed genes was the extracellular matrix FN-1. Patients carrying at least one −77VNTR >8(GT) allele, presented a strong upregulation of the FN-1 gene, while FN-1 was significantly downregulated in ≤8/≤8(GT)-repeats carrying patients, a phenomenon coinciding with the expression of the secondary IFNAR-1 transcript, which implies its potentially protective role against the development of HCC. This transcript comprises only a fibronectin III domain; FN fragments and modules can inhibit FN-matrix assembly by competing for FN-assembly sites, which could act as a feedback system to regulate FN levels on the cell surface (51). Specifically, the first type III domain of the FN molecule is important for the matrix assembly (52, 53), while even small fragments derived from this module regulate FN polymerization, inducing it at moderate concentrations but inhibiting it at high concentrations (54, 55). FN-1 is involved in HCC as it can be upregulated by HBxAg in an NFkB-dependent way (14), while is generally over-expressed in several cancers [reviewed in Ref. (56)]. NFkB was not found differentially expressed in our dataset (p = 0.16) but is linked to PI3K-AKT, which was the most significantly enriched signaling pathway, down to the MDM2 proto-oncogene and the apoptosis regulator Bcl-2 (Figure 5). FN-1 can trigger this pathway and transduce multiple intracellular signals that control cell cycle. Our data suggest that apart from NFkB, complementary mechanisms, partially controlled by FN-1 and involved in cell cycle regulation might exist.

This study is based on data mining and in silico analyses, thus our findings here generate a strong hypothesis about links between the expression of the truncated IFNAR-1 transcript either with the PI3K–AKT signaling pathway and/or with the HCC development. To test this hypothesis, further wet-lab studies will be needed involving knock-down of the transcript, transfection with a plasmid that will produce this transcript and quantitative analysis of the respective transcripts and proteins would confirm these observations and would shed light in these complex regulatory networks.

In conclusion, our results suggest very minimal (if any) involvement of the IFNAR-1 promoter polymorphisms in the expression levels of the IFNAR-1 major transcript but at the same time raise a potential and intriguing role for the −77VNTR regarding the regulation of downstream genes. Our study shows that the majority of modifications of the Interferome coincided with the production of the truncated IFNAR-1 transcript. Thus, further study of this truncated transcript could clarify the mechanistic features of the combined antiviral and anticancer roles of IFNAR-1.

Author Contributions

TK designed and conducted the analyses, evaluated the results, and wrote the manuscript. GP, DP, AH, JM, UG, PK, and GM wrote and revised the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

The work was supported the Medical Research Council, UK (Project Reference: MR/K010565/1).

Supplementary Material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/articles/10.3389/fimmu.2018.00777/full#supplementary-material.

References

1. Schweitzer A, Horn J, Mikolajczyk RT, Krause G, Ott JJ. Estimations of worldwide prevalence of chronic hepatitis B virus infection: a systematic review of data published between 1965 and 2013. Lancet (2015) 386(10003):1546–55. doi:10.1016/S0140-6736(15)61412-X

PubMed Abstract | CrossRef Full Text | Google Scholar

2. European Association for the Study of the Liver. EASL clinical practice guidelines: management of chronic hepatitis B virus infection. J Hepatol (2012) 57:167–85. doi:10.1016/j.jhep.2012.02.010

CrossRef Full Text | Google Scholar

3. Funk ML, Rosenberg DM, Lok AS. World-wide epidemiology of HBeAg-negative chronic hepatitis B and associated precore and core promoter variants. J Viral Hepat (2002) 9(1):52–61. doi:10.1046/j.1365-2893.2002.00304.x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Papatheodoridis GV, Manolakopoulos S, Liaw YF, Lok A. Follow-up and indications for liver biopsy in HBeAg-negative chronic hepatitis B virus infection with persistently normal ALT: a systematic review. J Hepatol (2012) 57(1):196–202. doi:10.1016/j.jhep.2011.11.030

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hadziyannis SJ, Papatheodoridis GV. Hepatitis B e antigen-negative chronic hepatitis B: natural history and treatment. Semin Liver Dis (2006) 26(2):130–41. doi:10.1055/s-2006-939751

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ribes J, Cleries R, Esteban L, Moreno V, Bosch FX. The influence of alcohol consumption and hepatitis B and C infections on the risk of liver cancer in Europe. J Hepatol (2008) 49(2):233–42. doi:10.1016/j.jhep.2008.04.016

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bosetti C, Levi F, Boffetta P, Lucchini F, Negri E, La Vecchia C. Trends in mortality from hepatocellular carcinoma in Europe, 1980-2004. Hepatology (2008) 48(1):137–45. doi:10.1002/hep.22312

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Bellentani S, Scaglioni F, Ciccia S, Bedogni G, Tiribelli C. HCV, HBV and alcohol – the dionysos study. Dig Dis (2010) 28(6):799–801. doi:10.1159/000324288

CrossRef Full Text | Google Scholar

9. Hsu YS, Chien RN, Yeh CT, Sheen IS, Chiou HY, Chu CM, et al. Long-term outcome after spontaneous HBeAg seroconversion in patients with chronic hepatitis B. Hepatology (2002) 35(6):1522–7. doi:10.1053/jhep.2002.33638

PubMed Abstract | CrossRef Full Text | Google Scholar

10. de Franchis R, Meucci G, Vecchi M, Tatarella M, Colombo M, Del Ninno E, et al. The natural history of asymptomatic hepatitis B surface antigen carriers. Ann Intern Med (1993) 118(3):191–4. doi:10.7326/0003-4819-118-3-199302010-00006

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jung MC, Pape GR. Immunology of hepatitis B infection. Lancet Infect Dis (2002) 2(1):43–50. doi:10.1016/S1473-3099(01)00172-4

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Jiang YG, Wang YM, Liu TH, Liu J. Association between HLA class II gene and susceptibility or resistance to chronic hepatitis B. World J Gastroenterol (2003) 9(10):2221–5. doi:10.3748/wjg.v9.i10.2221

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Li N, Gao Y-F, Zhang T-C, Chen P, Li X, Su F. Relationship between interleukin 18 polymorphisms and susceptibility to chronic hepatitis B virus infection. World J Hepatol (2012) 4(3):105–9. doi:10.4254/wjh.v4.i3.105

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Norton PA, Reis HM, Prince S, Larkin J, Pan J, Liu J, et al. Activation of fibronectin gene expression by hepatitis B virus x antigen. J Viral Hepat (2004) 11(4):332–41. doi:10.1111/j.1365-2893.2004.00555.x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. To WS, Midwood KS. Plasma and cellular fibronectin: distinct and independent functions during tissue repair. Fibrogenesis Tissue Repair (2011) 4:21. doi:10.1186/1755-1536-4-21

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Liu X-Y, Liu R-X, Hou F, Cui L-J, Li C-Y, Chi C, et al. Fibronectin expression is critical for liver fibrogenesis in vivo and in vitro. Mol Med Rep (2016) 14(4):3669–75. doi:10.3892/mmr.2016.5673

CrossRef Full Text | Google Scholar

17. Shi M, Fu J, Shi F, Zhang B, Tang Z, Jin L, et al. Transfusion of autologous cytokine-induced killer cells inhibits viral replication in patients with chronic hepatitis B virus infection. Clin Immunol (2009) 132(1):43–54. doi:10.1016/j.clim.2009.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Pham TN, Mercer SE, Michalak TI. Chronic hepatitis C and persistent occult hepatitis C virus infection are characterized by distinct immune cell cytokine expression profiles. J Viral Hepat (2009) 16(8):547–56. doi:10.1111/j.1365-2893.2009.01092.x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Parker BS, Rautela J, Hertzog PJ. Antitumour actions of interferons: implications for cancer therapy. Nat Rev Cancer (2016) 16(3):131–44. doi:10.1038/nrc.2016.14

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Domanski P, Colamonici OR. The type-I interferon receptor. the long and short of it. Cytokine Growth Factor Rev (1996) 7(2):143–51. doi:10.1016/1359-6101(96)00017-2

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Mogensen KE, Lewerenz M, Reboul J, Lutfalla G, Uze G. The type I interferon receptor: structure, function, and evolution of a family business. J Interferon Cytokine Res (1999) 19(10):1069–98. doi:10.1089/107999099313019

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Langer JA, Cutrone EC, Kotenko S. The Class II cytokine receptor (CRF2) family: overview and patterns of receptor-ligand interactions. Cytokine Growth Factor Rev (2004) 15(1):33–48. doi:10.1016/j.cytogfr.2003.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rusinova I, Forster S, Yu S, Kannan A, Masse M, Cumming H, et al. Interferome v2.0: an updated database of annotated interferon-regulated genes. Nucleic Acids Res (2013) 41(Database issue):D1040–6. doi:10.1093/nar/gks1215

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Liaw Y-F, Kao J-H, Piratvisuth T, Chan H, Chien R-N, Liu C-J, et al. Asian-Pacific consensus statement on the management of chronic hepatitis B: a 2012 update. Hepatol Int (2012) 6(3):531–61. doi:10.1007/s12072-012-9365-4

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Pestka S, Krause CD, Walter MR. Interferons, interferon-like cytokines, and their receptors. Immunol Rev (2004) 202:8–32. doi:10.1111/j.0105-2896.2004.00204.x

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Zhou J, Smith DK, Lu L, Poon VK, Ng F, Chen DQ, et al. A non-synonymous single nucleotide polymorphism in IFNAR1 affects susceptibility to chronic hepatitis B virus infection. J Viral Hepat (2009) 16(1):45–52. doi:10.1111/j.1365-2893.2008.01040.x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Zhou J, Lu L, Yuen MF, Lam TW, Chung CP, Lam CL, et al. Polymorphisms of type I interferon receptor 1 promoter and their effects on chronic hepatitis B virus infection. J Hepatol (2007) 46(2):198–205. doi:10.1016/j.jhep.2006.08.017

PubMed Abstract | CrossRef Full Text | Google Scholar

28. He XX, Chang Y, Jiang HJ, Tang F, Meng FY, Xie QH, et al. Persistent effect of IFNAR-1 genetic polymorphism on the long-term pathogenesis of chronic HBV infection. Viral Immunol (2010) 23(3):251–7. doi:10.1089/vim.2009.0102

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Song le H, Xuan NT, Toan NL, Binh VQ, Boldt AB, Kremsner PG, et al. Association of two variants of the interferon-alpha receptor-1 gene with the presentation of hepatitis B virus infection. Eur Cytokine Netw (2008) 19(4):204–10. doi:10.1684/ecn.2008.0137

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Muldoon J, Uriel A, Khoo S, Ollier WE, Hajeer AH. Novel IFN-alpha receptor promoter polymorphisms. Genes Immun (2001) 2(3):159–60. doi:10.1038/sj.gene.6363757

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Aucan C, Walley AJ, Hennig BJ, Fitness J, Frodsham A, Zhang L, et al. Interferon-alpha receptor-1 (IFNAR1) variants are associated with protection against cerebral malaria in the Gambia. Genes Immun (2003) 4(4):275–82. doi:10.1038/sj.gene.6363962

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Karamitros T, Papatheodoridis G, Dimopoulou E, Papageorgiou MV, Paraskevis D, Magiorkinis G, et al. The interferon receptor-1 promoter polymorphisms affect the outcome of caucasians with HBeAg-negative chronic HBV infection. Liver Int (2015) 35(12):2506–13. doi:10.1111/liv.12859

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Papatheodoridis GV, Manolakopoulos S, Margariti A, Papageorgiou MV, Kranidioti H, Katoglou A, et al. The usefulness of transient elastography in the assessment of patients with HBeAg-negative chronic hepatitis B virus infection. J Viral Hepat (2014) 21(7):517–24. doi:10.1111/jvh.12176

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics (2009) 25(16):2078–9. doi:10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods (2012) 9(4):357–9. doi:10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. Computer Science and Biology Proceedings of the German Conference on Bioinformatics. Hannover (1999).

Google Scholar

37. Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform (2013) 14(2):178–92. doi:10.1093/bib/bbs017

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Bray NL, Pimentel H, Melsted P, Pachter L. Erratum: near-optimal probabilistic RNA-seq quantification. Nat Biotechnol (2016) 34(8):888. doi:10.1038/nbt.3519

CrossRef Full Text | Google Scholar

39. Pimentel H, Bray NL, Puente S, Melsted P, Pachter L. Differential analysis of RNA-seq incorporating quantification uncertainty. Nature Methods (2017)14: 687.

Google Scholar

40. Hwang S, Kim E, Lee I, Marcotte EM. Systematic comparison of variant calling pipelines using gold standard personal exome variants. Sci Rep (2015) 5:17875. doi:10.1038/srep17875

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Karamitros T, Magiorkinis G. A novel method for the multiplexed target enrichment of MinION next generation sequencing libraries using PCR-generated baits. Nucleic Acids Res (2015) 43(22):e152. doi:10.1093/nar/gkv773

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Tsangaras K, Wales N, Sicheritz-Ponten T, Rasmussen S, Michaux J, Ishida Y, et al. Hybridization capture using short PCR products enriches small genomes by capturing flanking sequences (CapFlank). PLoS One (2014) 9(10):e109101. doi:10.1371/journal.pone.0109101

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Massirer KB, Hirata MH, Silva AE, Ferraz ML, Nguyen NY, Hirata RD. Interferon-alpha receptor 1 mRNA expression in peripheral blood mononuclear cells is associated with response to interferon-alpha therapy of patients with chronic hepatitis C. Braz J Med Biol Res (2004) 37(5):643–7. doi:10.1590/S0100-879X2004000500003

PubMed Abstract | CrossRef Full Text | Google Scholar

44. de Weerd NA, Matthews AY, Pattie PR, Bourke NM, Lim SS, Vivian JP, et al. A hot spot on interferon alpha/beta receptor subunit 1 (IFNAR1) underpins its interaction with interferon-beta and dictates signaling. J Biol Chem (2017) 292(18):7554–65. doi:10.1074/jbc.M116.773788

CrossRef Full Text | Google Scholar

45. Zhang G, deWeerd NA, Stifter SA, Liu L, Zhou B, Wang W, et al. A proline deletion in IFNAR1 impairs IFN-signaling and underlies increased resistance to tuberculosis in humans. Nat Commun (2018) 9(1):85. doi:10.1038/s41467-017-02611-z

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Karam RA, Rezk NA, Amer MM, Fathy HA. Immune response genes receptors expression and polymorphisms in relation to multiple sclerosis susceptibility and response to INF-beta therapy. IUBMB Life (2016) 68(9):727–34. doi:10.1002/iub.1530

CrossRef Full Text | Google Scholar

47. Matsuyama N, Mishiro S, Sugimoto M, Furuichi Y, Hashimoto M, Hijikata M, et al. The dinucleotide microsatellite polymorphism of the IFNAR1 gene promoter correlates with responsiveness of hepatitis C patients to interferon. Hepatol Res (2003) 25(3):221–5. doi:10.1016/S1386-6346(02)00269-3

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zhou J, Huang JD, Poon VK, Chen DQ, Chan CC, Ng F, et al. Functional dissection of an IFN-alpha/beta receptor 1 promoter variant that confers higher risk to chronic hepatitis B virus infection. J Hepatol (2009) 51(2):322–32. doi:10.1016/j.jhep.2009.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Santos JC, de Deus DMV, de Moura IMF, Lopes EP, Alves MRB, Coêlho MRCD. Association between the IFNA1 (-2C&#x2192;T) polymorphism and increased IFNAR1 gene expression levels in chronic hepatitis B infection. Intervirology (2015) 58(6):393–402. doi:10.1159/000444365

CrossRef Full Text | Google Scholar

50. de Weerd NA, Samarajiwa SA, Hertzog PJ. Type I interferon receptors: biochemistry and biological functions. J Biol Chem (2007) 282(28):20053–7. doi:10.1074/jbc.R700006200

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Bourdoulous S, Orend G, MacKenna DA, Pasqualini R, Ruoslahti E. Fibronectin matrix regulates activation of RHO and CDC42 GTPases and cell cycle progression. J Cell Biol (1998) 143(1):267–76. doi:10.1083/jcb.143.1.267

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Aguirre KM, McCormick RJ, Schwarzbauer JE. Fibronectin self-association is mediated by complementary sites within the amino-terminal one-third of the molecule. J Biol Chem (1994) 269(45):27863–8.

PubMed Abstract | Google Scholar

53. Hocking DC, Sottile J, McKeown-Longo PJ. Fibronectin’s III-1 module contains a conformation-dependent binding site for the amino-terminal region of fibronectin. J Biol Chem (1994) 269(29):19183–7.

PubMed Abstract | Google Scholar

54. Morla A, Ruoslahti E. A fibronectin self-assembly site involved in fibronectin matrix assembly: reconstruction in a synthetic peptide. J Cell Biol (1992) 118(2):421–9. doi:10.1083/jcb.118.2.421

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Morla A, Zhang Z, Ruoslahti E. Superfibronectin is a functionally distinct form of fibronectin. Nature (1994) 367(6459):193–6. doi:10.1038/367193a0

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Wang JP, Hielscher A. Fibronectin: how its aberrant expression in tumors may improve therapeutic targeting. J Cancer (2017) 8(4):674–82. doi:10.7150/jca.16901

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hepatitis virus B, interferon-α receptor type-1 promoter, hepatocellular carcinoma, interferon-α receptor, polymorphism, transcriptome, RNAseq, the Cancer Genome Atlas project

Citation: Karamitros T, Papatheodoridis G, Paraskevis D, Hatzakis A, Mbisa JL, Georgopoulou U, Klenerman P and Magiorkinis G (2018) Impact of Interferon-α Receptor-1 Promoter Polymorphisms on the Transcriptome of the Hepatitis B Virus-Associated Hepatocellular Carcinoma. Front. Immunol. 9:777. doi: 10.3389/fimmu.2018.00777

Received: 29 September 2017; Accepted: 28 March 2018;
Published: 16 April 2018

Edited by:

Linda F. Van Dyk, University of Colorado Denver, United States

Reviewed by:

Hridayesh Prakash, All India Institute of Medical Sciences, India
Sandip Chakraborty, College of Veterinary Sciences and Animal Husbandry, India

Copyright: © 2018 Karamitros, Papatheodoridis, Paraskevis, Hatzakis, Mbisa, Georgopoulou, Klenerman and Magiorkinis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Timokratis Karamitros, tkaram@pasteur.gr;
Gkikas Magiorkinis, gmagi@med.uoa.gr