Genetic Variants in TNFSF4 and TNFSF8 Are Associated With the Risk of HCV Infection Among Chinese High-Risk Population

Background: The tumor necrosis factor superfamily (TNFSF) and TNF receptor superfamily (TNFRSF) play important roles in the immune responses to infections. The aim of this study was to determine the impact of single nucleotide polymorphisms (SNPs) of several TNFSF/TNFRSF genes on the risk of hepatitis C virus (HCV) infection in the Chinese high-risk population.

Methods: The TNFSF4-rs1234313, TNFSF4-rs7514229, TNFSF8-rs3181366, TNFSF8-rs2295800, TNFRSF8-rs2298209, and TNFRSF8-rs2230625 SNPs were genotyped in 2309 uninfected controls, 597 subjects with spontaneous HCV clearance and 784 patients with persistent HCV infection using the TaqMan-MGB assay. The putative functions of the positive SNPs were determined using online bioinformatics tools.

Results: After adjusting for gender, age, high-risk population, alanine transaminase (ALT), aspartate aminotransferase (AST), IL28B-rs12979860 and rs8099917 genotypes, the non-conditional logistic regression showed that rs7514229-T, rs3181366-T, and rs2295800-C were associated with an increased risk of HCV infection (all P_FDR < 0.05). Combined analysis of rs7514229-T and rs3181366-T risk alleles showed that the subjects carrying 2–4 risk alleles were more susceptible to HCV infection compared with those lacking any risk allele (all P < 0.001). Furthermore, the risk of HCV infection increased with the number of risk alleles (P_trend < 0.001). In silico analysis showed that rs7514229, rs3181366, and rs2295800 polymorphisms may affect the transcription of mRNA by regulating miRNA binding, TF binding, and promoter activation, respectively, which may have biological consequences.

Conclusion: TNFSF4-rs7514229, TNFSF8-rs3181366, and TNFSF8-rs2295800 are associated with increased risk of HCV infection in the Chinese high-risk population.

Introduction

Hepatitis C is the result of hepatitis C virus (HCV) infection and currently afflicts around 185 million people worldwide, of which 71 million are chronically infected (WHO, 2018; Spearman et al., 2019). Approximately 55–85% of the infected patients subsequently develop chronic hepatitis C (CHC) due to lack of viral clearance, which can lead to decompensated cirrhosis, hepatocellular carcinoma (HCC), and even death (Chinese Society of Hepatology and Chinese Society of Infectious Disease, 2019). WHO and other international organizations have pledged to eliminate HCV infection by 2030, which is incumbent on the development of effective vaccines and therapeutics (WHO, 2018). Although direct-acting antiviral drugs (DAAs) have been effective against HCV, only 1.3% of the patients in China have received these drugs owing to their high costs (Wei et al., 2016; Li and Chung, 2019). Furthermore, the attempts to develop an HCV vaccine have been largely unsuccessful due to the high genetic variability of this virus. Therefore, the molecular mechanisms underlying HCV infection, especially virus–host interactions, need further elucidation to circumvent the aforementioned limitations (Heim et al., 2016).

Studies show that genetic factors are a major determinant of the host response to HCV infection (Matsuura and Tanaka, 2018). Genetic variants of immune-related genes such as IL28B (Thomas et al., 2009), HLA-DQB1 (Lee et al., 2018), and IFN (Welzel et al., 2009) have been implied in HCV infection. Tumor necrosis factor superfamily/tumor necrosis factor receptor superfamily (TNFSF/TNFRSF) proteins are expressed in immune cells, and are frequently activated or dysregulated in inflammatory diseases such as inflammatory bowel disease (IBD; De Voogd et al., 2016), systemic lupus erythematosus (SLE; Jackson and Davidson, 2019), rheumatoid arthritis (RA; Croft and Siegel, 2017), Crohn’s disease (CD; Hong et al., 2016), and hepatitis (Shin et al., 2016). Furthermore, genetic variants of TNFRSF1A (Yue et al., 2020), TNFRSF5 (Tian et al., 2019), TNFSF6 (Huang et al., 2020), and TNFRSF11B (Wu et al., 2019) influence the immune response to HCV infection. TNFSF4, TNFSF8, and their respective receptors also mediate the immune response in the manner similar to of sOX40L and sCD30L (Späth et al., 2017; Qin et al., 2018). However, it is still unclear whether polymorphisms in TNFSF4/TNFRSF4 and TNFSF8/TNFRSF8 have an effect on the host response to HCV infection.

In this study, we screened six single nucleotide polymorphisms (SNPs) of TNFSF4/TNFRSF4 and TNFSF8/TNFRSF8 in a Chinese cohort at high risk of HCV infection to determine their potential role in both HCV infection and CHC.

Materials and Methods

Study Population

A total of 3,976 subjects were recruited, including 816 hemodialysis (HD) patients from nine hospitals across southern China, 1,848 paid blood donors (PBD) from 20 villages within the Jiangsu Province, and 1312 people who use drugs (PWUD) from detoxification centers in Nanjing and Yixing City. In the 1990s, due to China’s poor medical level, HDs were usually infected with HCV through blood transmission. Besides, from 1980 to 1990, PBDs were monetarily compensated, and subsequently numerous donors were found to be infected with HCV. The plasma of some paid donors was separated and collected by plasmapheresis, and the other blood components that contained cross-contamination were returned to the donor. For PWUD, it was more likely to be the cross-use of contaminated needles and unsafe sex. The most common modes of transmission of HCV infection are blood transmission, sexual transmission, and mother-to-child transmission. As a result, we recruited HD, PBD, and PWUD as our research objects.

Participants were excluded for the following reasons: (1) failing to collect the blood or whose infection outcome undetermined; (2) HBV or HIV co-infection; (3) history of antiretroviral therapy; (4) age < 18 years and age > 80 years; (5) liver cirrhosis and other liver diseases; and (6) history of cancer or malignant tumor. The detailed flow diagram for recruiting participants is shown in Figure 1.

FIGURE 1

Figure 1. Flow diagram for inclusion of participants.

According to the diagnosis criteria of the “Guidelines for the prevention and treatment of hepatitis C (2019 version) (Chinese Society of Hepatology and Chinese Society of Infectious Disease, 2019),” the subjects were divided into the following groups: (1) uninfected controls—anti-HCV seronegative and HCV RNA seronegative, (2) spontaneous clearance—anti-HCV seropositive and HCV RNA seronegative, and (3) persistent infection—anti-HCV seropositive and HCV RNA seropositive. In addition, the spontaneous clearance and persistent infection groups were classified as HCV-infected.

Serological Testing

After individual surveys using questionnaires, 5 ml venous blood was collected from each subject for further testing. HCV antibodies and RNA load were detected within 4 h, and the remaining blood was used for DNA extraction. The presence of anti-HCV antibody (anti-HCV) was tested by ELISA (Shanghai Abbott Laboratories, Shanghai, China), and HCV-RNA was isolated using Trizol LS reagent (Takara Biotech, Tokyo, Japan). Genomic DNA of leukocytes was extracted by the standard phenolic chloroform extraction method involving proteinase K digestion, phenol-chloroform purification, and ethanol precipitation. All DNA samples were stored at −20°C.

SNP Selection and Genotyping Assays

The linkage disequilibrium (LD) data of the HapMap Phase II CHB (Chinese in Beijing) obtained from the 1000 Genomes Project database¹ was imported into the Haploview software (version 4.2; Broad Institute, Cambridge, MA, United States). The candidate SNPs were screened with minor allele frequency (MAF) ≥ 0.05 and correlation coefficient r² ≥ 0.8 as the thresholds. In addition, the sequences 2,000 bp upstream and downstream of the transcription initiation sites of TNFSF4/TNFRSF4 and TNFSF8/TNFRSF8 genes were also included in the analysis. The putative functional SNPs were then further screened using the RegulomeDB database² and UCSC Genome Browser.³ Previously published SNPs associated with immunological or infectious diseases were also retrieved. Six SNPs (TNFSF4-rs1234313, rs7514229; TNFSF8-rs3181366, rs2295800; TNFRSF4-rs2298209; and TNFRSF8-rs2230625) were finally selected for further analysis. The candidate SNPs were genotyped using TaqMan real-time PCR assay in the LightCycler 480 II Real-Time PCR System (Roche Diagnostics, Mannheim, Germany). The primer and probe sequences are shown in Supplementary Table 1. The reaction parameters consisted of preheating at 50°C for 2 min and pre-denaturation at 95°C for 10 min, followed by 45 cycles of denaturation at 95°C for 15 s and annealing, and extension at 60°C for 1 min. The success rate for each SNP was above 95%. The experiment was repeated with randomly selected 10% of the samples, with a consistency rate of 100%. Genotyping was performed in a manner blinded to the clinical data, and two blank controls were set up for each 384-well format for quality control.

In silico Analysis

The SNPs were functionally annotated using the SNPinfo website.⁴ The RegulomeDB online database⁵ was used to determine the regulatory role of the SNPs on the basis of RegulomeDB scores (Supplementary Table 2). The RNA Web Servers⁶ based on the latest Vienna RNA Package (version 2.4.16) was used to predict the secondary structures of single stranded RNA sequences and obtaining the minimum free energy (MFE), and the potential biological function was annotated using the UCSC Genome Bioinformatics website.⁷ The H3K27Ac histone marker data of seven cell lines (GM12878, H1-hESC, HSMM, HUVEC, K562, NHEK, and NHLF) was also analyzed.

Statistical Analysis

Deviations from Hardy–Weinberg equilibrium (HWE) for each SNP among the controls were analyzed with χ² test. Differences in demographic characteristics were described by mean ± SD or count (proportion) and compared by one-way ANOVA (for continuous variables) or the χ² test (for categorical variables). The association of the selected SNPs with HCV susceptibility and outcomes was estimated by constructing logistic regression models with age, gender, high-risk population, ALT, AST, IL28B-rs12979860, and IL28B-rs8099917 and each SNP as the covariates. The ORs with 95% CIs were calculated using co-dominant model, dominant model, recessive model, and additive model. All statistical analyses were performed by STATA 15.0 software (STATA Corp, College Station, TX, United States), and P < 0.05 was considered statistically significant. False discovery rate (FDR) correction was used to analyze the genotype distribution among the different groups.

Results

Basic Characteristics of Participants

As shown in Table 1, there were no significant differences in the distribution of the IL28B-rs12979860 genotypes among the three groups. In contrast, age, gender, ALT level, AST level, high-risk population, HCV genotype, and IL28B-rs8099917 (all P < 0.001) showed significant differences. The allele frequencies of five SNPs in the healthy uninfected controls were in accordance with the HWE (all P > 0.05), and only rs2298209 showed deviation (Supplementary Table 1).

TABLE 1

Table 1. Demographical and clinical characteristics of the HCV control, spontaneous clearance, and persistent infection populations.

Association Between Selected Gene Polymorphisms and HCV Infection Outcomes

To determine the association between the SNPs and the risk of HCV infection, the patients with natural clearance and persistent infection were grouped as HCV-infected and compared with the uninfected group. After adjusting for gender, age, high-risk population, ALT, AST, IL28B-rs12979860, and IL28B-rs8099917, the different genetic models showed that TNFSF4-rs7514229 (dominant model: adjusted OR 1.37, 95% CI 1.13–1.66, P = 0.001; recessive model: adjusted OR 6.40, 95% CI 3.70–11.05, P < 0.001; additive model: adjusted OR 1.52, 95% CI 1.29–1.78, P < 0.001), TNFSF8-rs3181366 (dominant model: adjusted OR 1.20, 95% CI 1.03–1.40, P = 0.018; recessive model: adjusted OR 1.50, 95% CI 1.15–1.97, P = 0.003; additive model: adjusted OR 1.21, 95% CI 1.07–1.36, P = 0.002) and TNFSF8-rs2295800 (dominant model: adjusted OR 1.22, 95% CI 1.04–1.42, P = 0.012) were significantly associated with the susceptibility to HCV infection, while no significant correlation was seen between TNFSF4-rs1234313, TNFRSF8-rs2230625, and the risk of HCV infection (Table 2). After correcting for multiple comparison, the difference between rs7514229-T, rs3181366-T, rs2295800-C, and HCV infection were still significant (all P_FDR < 0.05, Supplementary Table 3). In addition, there was no significant correlation between these SNPs and HCV infection chronicity when comparing the spontaneous clearance and persistent infection groups (all P > 0.05).

TABLE 2

Table 2. Distribution of TNFSF/TNFRSF genes among the uninfected control, spontaneous clearance, and persistent infection groups.

Independent Analysis and Combined Analysis

To further analyze the combined effect of rs7514229, rs3181366, and rs2295800 on HCV infection, we performed independent tests on the three SNPs. As shown in Supplementary Table 4, after adjusting for rs3181366 (P = 0.189) or rs3181366 and rs7514229 (P = 0.371), the effect of rs2295800 was significantly reduced. Therefore, the combined effects of the risk alleles rs7514229-T and rs3181366-T were assessed. As shown in Table 3, the presence of 2–4 risk alleles were linked to an increased risk of HCV infection compared lack of any risk allele (all P < 0.001). Furthermore, the risk of HCV infection increased with the number of risk alleles (P_trend < 0.001). Analysis of the combined risk genotypes (rs7514229-TT and rs3181366-TT) suggested that subjects with 1–2 risk genotypes were more susceptible (all P < 0.05), and both SNPs indicated a higher risk of HCV infection after the Cochran–Armitage trend test (P_trend = 0.010) (Supplementary Table 5).

TABLE 3

Table 3. The combined effects of risk alleles (rs7514229-T and rs3181366-T) on the risk of HCV infection.

Stratified Analysis and Interaction Analysis of Positive SNPs

A stratified analysis of the positive SNPs was next performed in the age, gender, ALT level, AST level, and high-risk subgroups. The additive model showed that compared with the rs7514229-G allele, a significantly increased risk of HCV infection was attributed to rs7514229-T in the subgroups of age (all P < 0.05), female (adjusted OR 1.77, 95% CI 1.42–2.21, P < 0.001), ALT ≤ 40 U/L (adjusted OR 1.62, 95% CI 1.36–1.92, P < 0.001), AST ≤ 40 U/L (adjusted OR 1.45, 95% CI 1.23–1.72, P < 0.001), HD (adjusted OR 1.61, 95% CI 1.10–2.34, P = 0.013) and PBD (adjusted OR 1.57, 95% CI 1.25–1.99, P < 0.001) (Table 4). No significant heterogeneity was found among age, AST, and high-risk population stratifications (all P > 0.05), while the difference between gender (P = 0.028) and ALT level (P = 0.011) was evident (Supplementary Table 6). Furthermore, rs3181366-T was associated with a significantly increased risk of HCV infection in the subgroups of age <50 years (adjusted OR 1.34, 95% CI 1.12–1.61, P = 0.001), male (adjusted OR 1.31, 95% CI 1.10–1.56, P = 0.002), ALT ≤ 40 U/L (adjusted OR 1.19, 95% CI 1.05–1.35, P = 0.008), AST ≤ 40 U/L (adjusted OR 1.20, 95% CI 1.06–1.36, P = 0.003) and PWUD (adjusted OR 1.64, 95% CI 1.31–2.04, P < 0.001) as opposed to the rs3181366-C allele. No significant heterogeneity was found among age, gender, ALT, and AST stratifications (All P > 0.05), while that between high-risk population (P = 0.007) stratification was evident (Supplementary Table 7).

TABLE 4

Table 4. Stratified analysis of rs7514229 and rs3181366 between HCV-infected and uninfected control groups.

To take into account the heterogeneity between SNPs and stratification, the interaction between the two SNPs and other factors in the context of HCV infection risk was analyzed (Table 5). Multiplicative interactions between rs7514229 genotypes and gender (P_interaction = 0.034) were assessed, and a significantly higher risk of HCV infection was observed for females (all P < 0.001) and males with TT genotype (P < 0.001) compared with males with the GG genotype. For rs3181366 genotypes and high-risk population (P_interaction = 0.016), the risk of HCV infection was significantly higher for PBD and PWUD compared with HD with CC genotype (all P < 0.05).

TABLE 5

Table 5. Interaction analysis between two SNPs and other factors on the risk of HCV infection.

In silico Analysis of Positive SNP Function

Rs7514229 is located in the three prime untranslated regions (3′UTR) of TNFSF4. Based on SNPinfo,⁸ rs7514229 is a putative microRNA-binding (miRNA-binding) site. To further analyze the effects of mutations on miRNA and transcriptional regulation, RNAfold⁹ was used to predict the secondary structure of mRNA and calculate the MFE of the centroid structure (one with minimal base pair distance). The secondary structure of the mRNA with mutant T allele differed from that of the wild G allele (Figure 2) and had a higher MFE (−16.60 kcal/mol vs. −17.30 kcal/mol).

FIGURE 2

Figure 2. The influence of rs7514229 variants on TNFSF4 mRNA centroid secondary structures. Changes in the local structure were illustrated by the RNAfold Web Server. The arrows indicate the location of the mutation (50 bases upstream and 50 bases downstream of the mutation) shown in bold type. The underlined sequence has overlapping nucleotide letters. TNFSF4 rs7514229 changes the local structure and the minimum free energy of the mRNA centroid secondary structure (one with the smallest base pair distance) from −17.30 kcal/mol (A) to −16.60 kcal/mol (B). The wild-type and mutant-type sequences are also listed (available at http://rna.tbi.univie.ac.at//cgi-bin/RNAWebSuite/RNAfold.cgi).

Rs3181366 is located on the intron of TNFSF8. Its RegulomeDB score¹⁰ is 5, indicating potential functions like transcription factor (TF) binding or DNase peak (Supplementary Table 2). Rs2295800 is also located on the intron of TNFSF8, and its RegulomeDB score of 5 is suggestive of similar regulatory functions (Supplementary Table 2). Based on ENCODE and UCSC genome browser (see footnote 7), we found that rs2295800 was located on the highest peak of the histone H3 acetylated lysine 27 (H3K27Ac) histone marker, which was confirmed by the enrichment of H3K27Ac via ChIP-seq assay (Figure 3). The acetylation of lysine 27 may enhance transcription by blocking the spread of the repressive methylated H3K27Me3. Thus, rs2295800 polymorphism may affect the transcription of mRNA by affecting promoter activation, which may translate to disease susceptibility.

FIGURE 3

Figure 3. Functional annotation for rs2295800 using ENCODE data from UCSC genome browser. Transcription Factor ChIP-seq Clusters show binding regions of transcription factors and other regulatory proteins. Based on ENCODE project and the UCSC genome browser data, rs2295800 is located on the highest peak of the histone H3 acetylated lysine 27 (H3K27Ac) histone marker and this track shows enrichment of the H3K27Ac histone mark across the genome as determined by a ChIP-seq assay. The light blue line indicates the position of SNP rs2295800 (available at http://genome.ucsc.edu/).

Discussion

Our results show that TNFSF4 rs7514229-T, TNFSF8 rs3181366-T, and TNFSF8 rs2295800-C are associated with an increased risk of HCV infection in the Chinese high-risk population. Furthermore, the presence of both rs7514229 and rs3181366 is significantly linked to a higher risk of HCV infection, and the risk increases with the number of risk alleles or genotypes. In silico analysis further showed that rs7514229, rs3181366, and rs2295800 polymorphisms may affect the transcription of mRNA by affecting miRNA binding, TF binding, and promoter activation, respectively, and thus mediate disease susceptibility.

TNFSF4, also known as OX40L, could be capable of interacting with its receptor on the late proliferation and sustained activation of T lymphocytes by extending the half-life of the cytokine mRNA (Vogel et al., 2013). Previous studies have shown that genetic variants of TNFSF4/TNFRSF4 are associated with immune disorders such as Behcet’s Disease (Lu et al., 2016), autoimmune thyroid diseases (AITDs; Song et al., 2016), and SLEs (Cortini et al., 2017). We found that TNFSF4 rs7514229-T (the mutant allele) was linked to an increased risk of HCV infection, and this effect was more evident in the age, female, lower ALT level (≤ 40 U/L), lower AST level (≤ 40 U/L), HD, and PBD subgroups. In addition, compared with male with rs7514229 GG genotype, a significant increased risk of HCV infection was observed for those who are female with GT/TT. Generally, being female is regarded as the common protective factor for hepatitis C because of the estradiol-related, more effective immune response (Fish, 2008). Based on bioinformatics analysis, we hypothesized that rs7514229 polymorphism may affect mRNA transcription by affecting the binding of miRNA, resulting in structural changes in the former that may regulate disease susceptibility. This hypothesis will have to be validated with functional studies on cellular models.

TNFSF8, also known as CD30L, interacts with its receptor on effector or memory T helper cells following activation by neutrophils, CD4⁺ T, and antigen-presenting cells, eventually mediating inflammatory diseases like IBD (Sun et al., 2008), RA (Barbieri et al., 2015), and CD (Hong et al., 2016). Some studies have also reported an association between the rs3181366 polymorphism and lung cancer (Wei et al., 2011) and myeloma bone disease (Durie et al., 2009). We observed that the TNFSF8 rs3181366-T mutant allele was linked to an elevated risk of HCV infection. Moreover, the effect of rs3181366-T was prominent in <50 years of age, male, lower ALT level (≤40 U/L), lower AST level (≤40 U/L), and PWUD subgroups. The mutant TNFSF8 rs2295800-C allele was also identified as a susceptibility locus for HCV infection, especially in the >50 years of age, lower ALT level (≤40 U/L), and lower AST level (≤40 U/L) subgroups. Based on the in silico analysis, the rs3181366 polymorphism may affect mRNA transcription by affecting TF binding and inducing structural changes with biological consequences. Rs2295800 was located on the highest peak of the H3K27Ac histone marker and may enhance transcription by blocking the spread of the repressive histone mark H3K27Me3. Thus, rs2295800 likely affects promoter activation and transcription.

The rs7514229 and rs3181366 SNPs showed independent effects on the risk of HCV infection, and the combined analysis of the rs7514229-T and rs3181366-T risk alleles suggested that subjects carrying two or more risk alleles were more susceptible to HCV infection. Furthermore, the risk increased with the number of risk alleles. Therefore, we hypothesized that genetic variants of TNFSF4 and TNFSF8 may have synergistic effects during the course of HCV infection. Since TNFSF8, TNFSF4, and their receptors play key roles in the differentiation and expansion of Th17 cells (Sun et al., 2010; Zhang et al., 2010), rs7514229 and rs3181366 may influence the outcome of HCV infection by affecting the Th17 population.

Although the results of this study are reliable and representative due to the reasonable design and large sample size, we must acknowledge some potential limitations. First, only six SNPs were selected and genotyped, which may be insufficient to fully analyze the relationship between TNFSF4 and TNFSF8 polymorphisms and HCV infection outcomes. However, the selection of candidate SNPs was based on strict and reasonable criteria. Also, this study was indeed a continuation of the later studies; the differences were the increase of the sample size of the population and the change of the pathway genes, and that FDR correction was used to solve the multiple comparisons problem of multiple SNPs in this study. In addition, this study lacked information on the prevalence of HCV genotypes, which may also affect the outcomes of HCV infection. Nevertheless, a previous study reported that HCV genotype 1 was the most common genotype in the Chinese population (Chen et al., 2017). Finally, the predicted biological functions of SNPs will need experimental validation.

Taken together, TNFSF4 rs7514229-T, TNFSF8 rs3181366-T, and TNFSF8-2295800-C are linked to an increased risk of HCV infection among the Chinese high-risk population. Our findings provide new insights into HCV screening and prevention, as well as vaccine development.

Data Availability Statement

The datasets presented in this article are not readily available because the principle of confidentiality and Chinese relevant policies. Requests to access the datasets should be directed to Department of Science and Technology, Nanjing Medical University, kejichu@njmu.edu.cn.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the Eastern Theater Command Centers for Disease Control and Prevention, Nanjing, China. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ZF, WC, PH, and MY designed and organized the study and supervised the whole project. ZF, JS, HX, ZG, CD, and JZ contributed to field survey, data collection, laboratory detection, and quality control. ZF, HF, CS, and PH performed the data cleansing and statistical analysis. WC, CW, YZ, and MY provided analysis tools and performed data interpretation. ZF, WC, PH, and MY wrote and critically revised the manuscript. All authors made substantial contributions to editing and drafting of the manuscript and read and approved the final manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (81773499), the Science Foundation for Distinguished Young Scholars of Jiangsu Province (BK20190106), Jiangsu Program for Young Medical Talents (QNRC2016616), Top Talents of Excellent young and middle aged talents support program of Jiangsu Province Hospital (YNRCZN007), Clinical Research Center for emerging respiratory infectious diseases (HS2020002), and the Key Project of Yunnan Province Applied Basic Research Program (2019FA005).

Conflict of Interest

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.

Supplementary Material

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

Supplementary Table 1 | Sequences of probes and primers specific for the TNFSF4/TNFRSF4, TNFSF8/TNFRSF8, and IL28B SNPs.

Supplementary Table 2 | RegulomeDB scores.

Supplementary Table 3 | Association between the TNFSF/TNFRSF SNPs and HCV infection outcomes in the co-dominant, dominant, recessive, and additive models of multivariable analyses.

Supplementary Table 4 | Independent analysis of SNPs associated with the risk of HCV infection.

Supplementary Table 5 | The combined effects of risk genotypes (rs7514229-TT and rs3181366-TT) on the risk of HCV infection.

Supplementary Table 6 | Stratified analyses of the association between rs7514229 and the risk of HCV infection.

Supplementary Table 7 | Stratified analyses of the association between rs3181366 and the risk of HCV infection.

Footnotes

1. ^ http://www.1000genomes.org/
2. ^ http://www.regulomedb.org/
3. ^ https://genome.ucsc.edu/
4. ^ https://snpinfo.niehs.nih.gov/cgi-bin/snpinfo/snpfunc.cgi
5. ^ http://www.regulomedb.org/
6. ^ http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi
7. ^ http://genome.ucsc.edu/
8. ^ https://snpinfo.niehs.nih.gov/cgi-bin/snpinfo/snpfunc.cgi
9. ^ http://rna.tbi.univie.ac.at//cgi-bin/RNAWebSuite/RNAfold.cgi
10. ^ https://www.regulomedb.org/regulome-search/

References

Barbieri, A., Dolcino, M., Tinazzi, E., Rigo, A., Argentino, G., Patuzzo, G., et al. (2015). Characterization of CD30/CD30L(+) Cells in Peripheral Blood and Synovial Fluid of Patients with Rheumatoid Arthritis. J. Immunol. Res. 2015:729654.

Google Scholar

Chen, Y., Yu, C., Yin, X., Guo, X., Wu, S., and Hou, J. (2017). Hepatitis C virus genotypes and subtypes circulating in Mainland China. Emerg. Microb. Infect. 6:e95.

Google Scholar

Chinese Society of Hepatology and Chinese Society of Infectious Disease. (2019). CMA Guidelines for the prevention and treatment of hepatitis C (2019 version). J. Clin. Hepatol. 35, 2670–2686.

Google Scholar

Cortini, A., Ellinghaus, U., Malik, T. H., Cunninghame, G. D. S., Botto, M., and Vyse, T. J. (2017). B cell OX40L supports T follicular helper cell development and contributes to SLE pathogenesis. Ann. Rheum. Dis. 76, 2095–2103. doi: 10.1136/annrheumdis-2017-211499

PubMed Abstract | CrossRef Full Text | Google Scholar

Croft, M., and Siegel, R. M. (2017). Beyond TNF: TNF superfamily cytokines as targets for the treatment of rheumatic diseases. Nat. Rev. Rheumatol. 13, 217–233. doi: 10.1038/nrrheum.2017.22

PubMed Abstract | CrossRef Full Text | Google Scholar

De Voogd, F. A., Gearry, R. B., and Day, A. S. (2016). Osteoprotegerin: A novel biomarker for inflammatory bowel disease and gastrointestinal carcinoma. J. Gastroenterol. Hepatol. 31, 1386–1392. doi: 10.1111/jgh.13324

PubMed Abstract | CrossRef Full Text | Google Scholar

Durie, B. G., Van Ness, B., Ramos, C., Stephens, O., Haznadar, M., Hoering, A., et al. (2009). Genetic polymorphisms of EPHX1, Gsk3beta, TNFSF8 and myeloma cell DKK-1 expression linked to bone disease in myeloma. Leukemia 23, 1913–1919. doi: 10.1038/leu.2009.129

PubMed Abstract | CrossRef Full Text | Google Scholar

Fish, E. N. (2008). The X-files in immunity: sex-based differences predispose immune responses. Nat. Rev. Immunol. 8, 737–744. doi: 10.1038/nri2394

PubMed Abstract | CrossRef Full Text | Google Scholar

Heim, M. H., Bochud, P. Y., and George, J. (2016). Host - hepatitis C viral interactions: The role of genetics. J. Hepatol. 65(1 Suppl.), S22–S32.

Google Scholar

Hong, S. N., Park, C., Park, S. J., Lee, C. K., Ye, B. D., Kim, Y. S., et al. (2016). Deep resequencing of 131 Crohn’s disease associated genes in pooled DNA confirmed three reported variants and identified eight novel variants. Gut 65, 788–796. doi: 10.1136/gutjnl-2014-308617

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, P., Wang, C. H., Zhuo, L. Y., Xia, X. S., Yang, S., Zhang, J. W., et al. (2020). Polymorphisms rs763110 in FASL is linked to hepatitis C virus infection among high-risk populations. Br. J. Biomed. Sci. 77, 112–117. doi: 10.1080/09674845.2020.1747182

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, S. W., and Davidson, A. (2019). BAFF inhibition in SLE-Is tolerance restored? Immunol. Rev. 292, 102–119. doi: 10.1111/imr.12810

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. H., Huang, Y. H., Chen, H. Y., Khor, S. S., Chang, Y. H., Lin, Y. J., et al. (2018). Human leukocyte antigen variants and risk of hepatocellular carcinoma modified by hepatitis C virus genotypes: A genome-wide association study. Hepatol. 67, 651–661. doi: 10.1002/hep.29531

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D. K., and Chung, R. T. (2019). Overview of Direct-Acting Antiviral Drugs and Drug Resistance of Hepatitis C Virus. Methods Mol. Biol. 1911, 3–32. doi: 10.1007/978-1-4939-8976-8_1

CrossRef Full Text | Google Scholar

Lu, S., Song, S., Hou, S., and Yang, P. (2016). Association of TNFSF4 Polymorphisms with Vogt-Koyanagi-Harada and Behcet’s Disease in Han Chinese. Sci. Rep. 6:37257.

Google Scholar

Matsuura, K., and Tanaka, Y. (2018). Host genetic variations associated with disease progression in chronic hepatitis C virus infection. Hepatol. Res. 48, 127–133. doi: 10.1111/hepr.13042

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, W., Zhao, C., Yang, Z., Zhu, S., and Li, M. S. Y. H. (2018). Research advances in the role of OX40/OX40L in the immune pathogenesis of chronic hepatitis. J. Clin. Hepatol. 34, 1103–1106.

Google Scholar

Shin, G. C., Kang, H. S., Lee, A. R., and Kim, K. H. (2016). Hepatitis B virus-triggered autophagy targets TNFRSF10B/death receptor 5 for degradation to limit TNFSF10/TRAIL response. Autophagy 12, 2451–2466. doi: 10.1080/15548627.2016.1239002

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, R. H., Wang, Q., Yao, Q. M., Shao, X. Q., Li, L., Wang, W., et al. (2016). TNFSF4 Gene Variations Are Related to Early-Onset Autoimmune Thyroid Diseases and Hypothyroidism of Hashimoto’s Thyroiditis. Int. J. Mol. Sci. 2016:17.

Google Scholar

Späth, F., Wibom, C., Krop, E. J., Johansson, A. S., Bergdahl, I. A., Vermeulen, R., et al. (2017). Biomarker Dynamics in B-cell Lymphoma: A Longitudinal Prospective Study of Plasma Samples Up to 25 Years before Diagnosis. Cancer Res. 77, 1408–1415. doi: 10.1158/0008-5472.can-16-2345

PubMed Abstract | CrossRef Full Text | Google Scholar

Spearman, C. W., Dusheiko, G. M., Hellard, M., and Sonderup, M. (2019). Hepatitis C. Lancet 394, 1451–1466. doi: 10.1016/S0140-6736(19)32320-7

CrossRef Full Text | Google Scholar

Sun, X., Somada, S., Shibata, K., Muta, H., Yamada, H., Yoshihara, H., et al. (2008). A critical role of CD30 ligand/CD30 in controlling inflammatory bowel diseases in mice. Gastroenterology 134, 447–458. doi: 10.1053/j.gastro.2007.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, X., Yamada, H., Shibata, K., Muta, H., Tani, K., Podack, E. R., et al. (2010). CD30 ligand/CD30 plays a critical role in Th17 differentiation in mice. J. Immunol. 185, 2222–2230. doi: 10.4049/jimmunol.1000024

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, D. L., Thio, C. L., Martin, M. P., Qi, Y., Ge, D., O’Huigin, C., et al. (2009). Genetic variation in IL28B and spontaneous clearance of hepatitis C virus. Nature 461, 798–801.

Google Scholar

Tian, T., Huang, P., Wu, J., Wang, C., Fan, H., Zhang, Y., et al. (2019). CD40 polymorphisms were associated with HCV infection susceptibility among Chinese population. BMC Infect. Dis. 19:840. doi: 10.1186/s12879-019-4482-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogel, K. U., Edelmann, S. L., Jeltsch, K. M., Bertossi, A., Heger, K., Heinz, G. A., et al. (2013). Roquin paralogs 1 and 2 redundantly repress the Icos and Ox40 costimulator mRNAs and control follicular helper T cell differentiation. Immunity 38, 655–668. doi: 10.1016/j.immuni.2012.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, L., Li, J., Yang, X., Wang, G., Feng, B., Hou, J., et al. (2016). Nationwide survey of specialist knowledge on current standard of care (Peg-IFN/RBV) and barriers of care in chronic hepatitis C patients in China. J. Gastroenterol. Hepatol. 31, 1995–2003. doi: 10.1111/jgh.13399

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, S., Niu, J., Zhao, H., Liu, Z., Wang, L. E., Han, Y., et al. (2011). Association of a novel functional promoter variant (rs2075533 C>T) in the apoptosis gene TNFSF8 with risk of lung cancer–a finding from Texas lung cancer genome-wide association study. Carcinogenesis 32, 507–515. doi: 10.1093/carcin/bgr014

PubMed Abstract | CrossRef Full Text | Google Scholar

Welzel, T. M., Morgan, T. R., Bonkovsky, H. L., Naishadham, D., Pfeiffer, R. M., Wright, E. C., et al. (2009). Variants in interferon-alpha pathway genes and response to pegylated interferon-Alpha2a plus ribavirin for treatment of chronic hepatitis C virus infection in the hepatitis C antiviral long-term treatment against cirrhosis trial. Hepatology 49, 1847–1858. doi: 10.1002/hep.22877

PubMed Abstract | CrossRef Full Text | Google Scholar

WHO (2018). World Health Organization. Global hepatitis report. Geneva: World Health Organization.

Google Scholar

Wu, J. J., Huang, P., Yue, M., Wang, C. H., Wu, C., Shao, J. G., et al. (2019). Association between TNFRSF11A and TNFRSF11B gene polymorphisms and the outcome of hepatitis C virus infection. Zhonghua Liu Xing Bing Xue Za Zhi. 40, 1291–1295.

Google Scholar

Yue, M., Huang, P., Wang, C., Fan, H., Tian, T., Wu, J., et al. (2020). Genetic Variation on TNF/LTA and TNFRSF1A Genes is Associated with Outcomes of Hepatitis C Virus Infection. Immunol. Invest. 2020, 1–11. doi: 10.1080/08820139.2019.1708384

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Zhong, W., Hinrichs, D., Wu, X., Weinberg, A., Hall, M., et al. (2010). Activation of OX40 augments Th17 cytokine expression and antigen-specific uveitis. Am. J. Pathol. 177, 2912–2920. doi: 10.2353/ajpath.2010.100353

PubMed Abstract | CrossRef Full Text | Google Scholar