Next Article in Journal
A Highly Conserved Epitope (RNNQIPQDF) of Porcine teschovirus Induced a Group-Specific Antiserum: A Bioinformatics-Predicted Model with Pan-PTV Potential
Next Article in Special Issue
Th1-Polarized, Dengue Virus-Activated Human Mast Cells Induce Endothelial Transcriptional Activation and Permeability
Previous Article in Journal
Multiplexed tat-Targeting CRISPR-Cas9 Protects T Cells from Acute HIV-1 Infection with Inhibition of Viral Escape
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Polymorphisms in Processing and Antigen Presentation-Related Genes and Their Association with Host Susceptibility to Influenza A/H1N1 2009 Pandemic in a Mexican Mestizo Population

1
HLA Laboratory, Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas, Mexico City 14080, Mexico
2
Tobacco Smoking and COPD Research Department, Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas, Mexico City 14080, Mexico
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2020 / Revised: 21 October 2020 / Accepted: 24 October 2020 / Published: 29 October 2020
(This article belongs to the Special Issue Host Factors in Viral Infections)

Abstract

:
(1) Background: The influenza A/H1N1 pdm09 virus rapidly spread throughout the world. Despite the inflammatory and virus-degradation pathways described in the pathogenesis of influenza A virus (IAV) infection, little is known about the role of the single nucleotide polymorphisms (SNPs) in the genes involved in the processing and antigenic presentation-related mechanisms. (2) Methods: In this case-control study, we evaluated 17 SNPs in five genes (TAP1, TAP2, TAPBP, PSMB8, and PSMB9). One hundred and twenty-eight patients with influenza A/H1N1 infection (INF-P) and 111 healthy contacts (HC) were included; all of them are Mexican mestizo. (3) Results: In allele and genotype comparison, the rs241433/C allele (TAP2), as well as AG haplotype (rs3763365 and rs4148882), are associated with reduced risk for influenza A/H1N1 infection (p < 0.05). On the other hand, the rs2071888G allele (TAPBP) and GG haplotype (rs3763365 and rs9276810) are associated with a higher risk for influenza A/H1N1 infection. In addition, after adjustment for covariates, the association to a reduced risk for influenza A/H1N1 infection remains with rs241433/C allele (p < 0.0001, OR = 0.24, 95% CI = 0.13–0.43), and the association with TAPBP is also maintained with the G allele (p = 0.0095, OR = 1.89, 95% CI = 1.17–3.06) and GG genotype models (p < 0.05, OR = 2.18, 95% CI = 1.27–3.74). (4) Conclusion: The rs241433/C allele and AC genotype (TAP2) and the AG haplotype are associated with a reduced risk for influenza A/H1N1 infection. In addition, the rs2071888/G allele and GG genotype (TAPBP) and the GG haplotype are associated with a higher risk for developing influenza A/H1N1 infection in a Mexican mestizo population.

Graphical Abstract

1. Introduction

In April 2009, the Mexican Secretary of Health reported the emergence of the influenza A/H1N1 pdm09 virus, which shares molecular characteristics of swine, avian, and human influenza viruses, and rapidly spread throughout the world [1,2]. In Mexico, the lethality rate was estimated to be 5% for confirmed cases of influenza A/H1N1, with several comorbidities such as immunosuppression, pre-existing respiratory and cardiovascular diseases, diabetes, tobacco smoking, and obesity, some of the most important predictors of the disease severity [2,3,4]. The most frequent signs and symptoms are fever, dyspnea, cough, myalgias, arthralgias, headache, and dyspnea. The most frequent adverse outcomes are mechanical ventilation, intensive unit care (ICU) admission, and death [1,5].
Additionally, previous studies have described several clinical and laboratory predictors of mortality or disease outcome. These parameters include sex (male), delayed medical attention and treatment onset, intensive care unit admission, low lymphocyte and platelet count, high creatine phosphokinase (CPK), serum creatine, blood urea nitrogen (BUN), and lactate dehydrogenase (LDH) [6,7,8].
Previous studies have determined that one of the most critical pathological mechanisms in this viral infection is systemic dysregulation of the inflammatory response correlated with disease severity [9,10,11]. In this sense, genetic variants in proinflammatory genes have been associated with host susceptibility for influenza A/H1N1 infection. For example, single nucleotide polymorphisms (SNPs) in TNF, LTA, IL1B, and IL6 are associated with a higher risk for influenza A/H1N1 infection in the Mexican population [11,12].
However, antigen presentation by influenza-infected cells through Human leukocyte antigens (HLA) class I is necessary to mount an effective antiviral immune response. It has been described that cytotoxic T lymphocytes (CTL) kill influenza-infected cells upon recognition of diverse HLA class I peptide complexes at the cell surface, including short viral peptides of 8–12 aa in length [13,14]. The assembly of these peptides occurs in the endoplasmic reticulum (ER) of infected cells and involves the transporter associated with antigen processing (TAP) and Tapasin (encoded by TAPBP gene), as well as ER chaperones. TAP is a critical factor comprising two subunits, TAP1 and TAP2 (encoded by TAP1 and TAP2 genes, respectively), required for translocation of peptides from the cytosol into the ER. On the other hand, Tapasin is another critical cofactor required for the assembly of HLA class I with exogenous peptides, which were obtained by intracellular degradation through proteasome, formed by proteasome subunit beta type 8 (PSMB8, encoded by PSMB8 gene) and 9 (PSMB9, encoded by PSMB9 gene) (previously known as large multifunctional peptidase 7 (LMP7) and 2 (LMP2), respectively) subunits [15,16,17,18].
Various studies have demonstrated that different HLA class I alleles are associated with a higher risk for influenza A/H1N1 infection in several populations [19,20]. Moreover, the role of SNPs in processing and antigen presentation-related genes has been studied in hypersensitivity pneumonitis (HP) [21], ankylosing spondylitis [22], and human papillomavirus (HPV) infection [23,24], but not in influenza A/H1N1 infection.
For all mentioned above, the study aimed to evaluate the association between processing and antigen presentation-related gene (TAP1, TAP2, TAPBP, PSMB8, and PSMB9) SNPs with the genetic susceptibility for influenza A/H1N1 infection in a Mexican mestizo population.

2. Materials and Methods

A retrospective case-control study was performed between May 2009 and March 2012 using nasopharyngeal swab samples from patients hospitalized with clinical symptoms of influenza and confirmed diagnosis by PCR techniques as previously stated [19].

2.1. Ethics Statements

This study was reviewed and approved by the Institutional Committees for Investigation, Ethics in Research, and Biosecurity of the Instituto Nacional de Enfermedades Respiratorias Ismael Cosío Villegas (INER) (approbation number: B05–10, approved on 2 March 2010). All participants were informed about the protocol’s aim after being given a detailed description of the study and invited to participate as volunteers. All of them signed an informed consent paper and supplied a privacy statement describing the legal protection of personal data; both documents approved by the institutional Research and Ethics in Research Committees.
A trained and responsible physician informed the patient about a possible diagnosis of influenza. All participants were invited to participate voluntarily by a trained researcher. All participants received medical care even if they chose not to participate in the research protocol. A clinical investigation was conducted according to the principles expressed in the Declaration of Helsinki. All experiments were performed following the relevant guidelines and regulations. The STREGA (STrengthening the REporting of Genetic Association) guidelines were considered to design this genetic association study [25].

2.2. Diagnosis

All patients underwent the QuickVue Influenza A + B test (Quidel, San Diego, CA, USA). Additionally, peripheral blood from patients and asymptomatic healthy contacts (HC) was collected; these interventions adhered to the guidelines and recommendations of the Centers for Disease Control and Prevention of the United States (CDC) and the World Health Organization (WHO). Subsequently, the RespiFinder assay was done, and the presence of the influenza A/H1N1 virus was confirmed.

2.3. Population

Our study’s subjects were evaluated according to a clinical questionnaire, which contained questions regarding influenza vaccination. The inclusion period was from May 2009 to March 2012, so the subjects examined during the study belonged to three different influenza waves in Mexico. A total of 261 subjects with clinical symptoms and signs of influenza were enrolled, of which 128 were positive (INF-P) and 133 were negative for influenza A/H1N1 infection. On admission to the hospital, all patients with suspected influenza were treated with oseltamivir. In this study, most patients sought medical attention 3–4 days after the onset of symptoms. All clinical data were collected retrospectively; nonetheless, blood samples and nasopharyngeal swabs were collected at the hospital during the first 24 h after admission. We also included the variable “died” during the patients’ hospitalization to perform a mortality analysis.
In the healthy contacts (HC) group, 111 biologically unrelated, intrahousehold contacts were enrolled who were asymptomatic and in close, continuous contact with influenza A/H1N1 patients before and during hospital admission. None of the HCs developed symptoms or was hospitalized for influenza virus infection. Peripheral blood from these individuals was tested for anti-A/H1N1 antibody titers by hemagglutination inhibition (HI) assay according to the method described by Julkunen I, 1985 [18] to determine exposure to the virus. All subjects included in this group had HI titers greater than 1:16, indicating direct contact with the virus. According to the samples’ availability in the biobank generated since 2009 in the HLA Laboratory, the participants were selected.

2.4. DNA Extraction

Firstly, we obtained 15 mL of peripheral blood via venipuncture in two EDTA tubes (as anticoagulant) from all 472 subjects. After that, the DNA extraction was performed using the commercial BDtract Genomic DNA isolation kit (Maxim Biotech, San Francisco, CA, USA). The DNA was quantified by UV absorption spectrophotometry at the 260-nm wavelength employing a NanoDrop 2000 device (Thermo Scientific, Wilmington, DE, USA). Contamination with organic compounds and proteins was determined by measuring the ratio absorbance at 260/280. Samples were considered of good quality when the ratio was ~1.8. All samples were adjusted to 15 µg/µL for subsequent genotyping.

2.5. SNP Selection

Four SNPs (rs1057149, rs2127679, rs4148882, and 41561219) in TAP1; four (rs13501, rs241433, rs241441, and rs2071544) in TAP2; two (rs2071888 and rs2282851) in TAPBP; four (rs2071542, rs2071543, rs3763365, and rs9276810) in PSMB8; and three (rs17587, rs241418, and rs2071534) in PSMB9 genes were selected based on a bibliographic search in PubMed (NCBI), identifying polymorphisms previously associated with other virus infectious and respiratory diseases. The molecular characteristics of the SNPs evaluated are described in Supplementary Table S1 (frequencies were obtained from the 1000 genomes project [26]).

2.6. SNP Genotyping

The allele discrimination was performed using commercial TaqMan probes (Applied Biosystems, San Francisco CA, USA), employing qPCR in a 7300 Real-Time PCR System (Applied Biosystems/Thermo Fisher Scientific Inc., Singapore), and the analysis performed by sequence detection software version 1.4 (Applied Biosystems, San Francisco, CA, USA). Four controls without template (contamination controls) were included for each genotyping plate, and 5% of the genotyped in duplicate as controls for allele assignment.

2.7. Bioinformatic Analysis and in Silico Analysis

HaploReg v4.1 [27] was consulted to determine the associated SNPs’ putative biologic effect. In addition, with the RNAfold web server [26], we introduce the section’s sequence that includes the SNPs to elucidate for changes in the secondary RNA structures. We also conducted an in silico analysis with the associated SNPs to evaluate the possible function or biological consequences across diverse software, considering if they are coding or noncoding SNPs. Supplementary Figure S1 shows the algorithm used for the in silico analysis of the associated SNPs.

2.8. Hardy–Weinberg Equilibrium and Haplotypes

The Hardy–Weinberg Equilibrium (HWE) was assessed using SNPStats online software for the 17 SNPs evaluated [28]. The haplotype analysis was performed using Haploview software v. 4.2 [29], using the Gabriel et al. [30] established criteria.

2.9. Statistical Analysis

Clinical quantitative variables were analyzed using SPSS for Mac Os, v20.0 (SPSS software, IBM, Chicago, IL, USA). Kolmogorov–Smirnov normality test was used, and parametric or nonparametric tests were used as appropriate. Genotype analysis and genetic association models were carried out with Pearson’s chi-squared and Fisher’s exact tests using Epi Info v. 7.1 [31] and Epi Dat v.3.1 software, and 2 × N and 2 × 2 contingency tables were made to estimate the genetic association for Influenza A/H1N1 infection and mortality. A p-value < 0.05 was considered statistically significant. Bonferroni correction (multiple testing) was applied to adjust the significant values of the genetic association results. Comparisons were made between case and control groups. After that, a logistic regression analysis was performed to adjust for potential confounding variables (sex and body mass index (BMI)) using Plink v. 1.07 [32].

3. Results

3.1. Demographic Variables in Case and Control Groups

One hundred and twenty-eight patients with influenza A/H1N1 infection (INF-P) and 111 HCs were included. When we compared demographic data between groups, we found that there were more males in the INF-P group, and they had a higher body mass index (BMI) compared with HCs (p < 0.05) (Table 1). Additionally, we found that the most common comorbidity was hypertension (13.28%), followed by asthma, type 2 diabetes mellitus, and COPD. Besides, the most frequent symptom was fever (75.78%), followed by dyspnea (72.66%) and cough (53.13%). Finally, pneumonia was the most frequent complication (59.38%). The noncomparable clinical features of the INF-P group are shown in Supplementary Table S2.

3.2. Demographic and Clinical Variables in Survivors and Nonsurvivors

Additionally, we performed an intracase analysis to evaluate which factors are associated with a higher risk of mortality in those patients with influenza A/H1N1 virus infection. Interestingly, we found that nonsurvivors subjects presented a lower platelet count, and higher glucose levels, CPK, and LDH than survivors subjects (p < 0.05). Most of the nonsurvivors subjects were in ICU, compared with survivors group subjects (p < 0.05). Complete results are shown in Appendix A, Table A1.

3.3. Hardy–Weinberg Equilibrium

Six polymorphisms evaluated did not meet the HWE (four in TAP2—rs13501, rs241433, rs241441, rs2071544; and two in TAPBP—rs2071888 and rs2282851) for the HC group. However, allele and genotype comparisons were carried out in all SNPs selected for the study. The results are shown in Supplementary Table S3.

3.4. Allele and Genotype Frequencies

We compared allele and genotype frequencies between INF-P and HC groups for the 17 SNPs selected for the study. In Table 2, the statistically significant results are shown; whereas in Supplementary Table S4, the allele and genotype frequencies for nonsignificative SNPs are shown. For TAP1, PSMB8, and PSMB9 SNPs, we did not find statistically significant differences in allele or genotype frequencies between INF-P and HC groups. However, the TAP2 rs241433/AC genotype and C allele are associated with a reduced risk for influenza A/H1N1 infection (p < 0.001, OR = 0.24, 95% CI = 0.14–0.22; OR = 0.44, 95% CI = 0.29–0.66, respectively), and the association is maintained after Bonferroni correction (p < 0.001).
On the other hand, TAPBP rs2071888/GG genotype and G allele are associated with higher risk for influenza A/H1N1 infection (p < 0.001, OR = 5.64, 95% CI = 2.00–15.86; p = 0.009, OR = 1.62, 95% CI = 1.13–2.34), and only the genotype association is maintained after Bonferroni correction (p = 0.014). In addition, the TAPBP rs2282851/T allele is associated with reduced risk for Influenza A/H1N1 infection (p = 0.048, OR = 0.69, 95% CI = 0.40–0.99). However, the association is not maintained after Bonferroni correction (p > 0.05).

3.5. Allele and Genotype Frequencies from Survivors and Nonsurvivors

We compared allele and genotype frequencies of the associated SNPs between Survivors and Nonsurvivors from the INF-P group. We did not find statistically significant differences in allele or genotype frequencies between groups. Complete results are shown in Appendix A, Table A2.

3.6. Logistic Regression Analysis

A logistic regression analysis adjusting by possible confounding covariates (sex and BMI) was carried out, which results in significant differences between INF-P and HC groups (Table 3). Interestingly, we found that the association with a reduced risk for Influenza A/H1N1 infection remains significant with the rs241433/C allele (p < 0.0001, OR = 0.24, 95% CI = 0.13–0.43). However, genotype analysis was not possible to perform with this SNP because we did not find minor allele homozygous in our population.
On the other hand, association to higher risk for Influenza A H1N1 infection with rs2071888/G allele (p = 0.0095, OR = 1.89, 95% CI = 1.17–3.06) and genotype (p = 0.005, OR = 2.18, 95% CI = 1.27–3.74) are also maintained.

3.7. Logistic Regression Analysis in Survivors and Nonsurvivors

After allele and genotype frequencies comparison between survivors and nonsurvivors, we performed a logistic regression analysis to adjust for possible confounding variables. However, we did not find statistically significant differences between groups (p > 0.05, table not shown).

3.8. Bioinformatic Analysis

We found that the rs241433/A allele has reduced affinity with the myeloid and B-cell development regulatory unwinding protein 1 (UP1). On the other hand, the rs2071888/C, affinity with hepatocyte nuclear factor 4 (HNF4) protein, a nuclear receptor, and transcription factor are increased (Table 4).
Moreover, no differences in predicted secondary messenger ribonucleotide acid (mRNA) structures were found in rs241433. However, for rs2071888, we can observe differences in the number of loops with rs2071888/C allele (Appendix B, Figure A1D).

3.9. In Silico Analysis

The rs2071888 (C/G) change from medium size and polar amino acid (threonine) change to large size, and basis (arginine) could alter the polarity of the protein structure (https://web.expasy.org/variant_pages/VAR_010253.html, https://www.uniprot.org/uniprot/O15533). The arginine presence increases protein stability (Gibbs free energy in 260 position = 0.31 kcal/mol at 25 °C and pH 7) (https://folding.biofold.org/cgi-bin/i-mutant2.0.cgi).
In silico analysis for rs241433 did not show sites of alternative splicing (http://wangcomputing.com/assp/, http://www.cbs.dtu.dk/services/NetGene2/) or micro-RNA (miRNAs) generation (http://www.mirbase.org/, https://tanuki.ibisc.univ-evry.fr/evryrna/mirnafold/mirnafold_form).

3.10. Linkage Disequilibrium (LD) and Haplotype Analysis

The haplotype analysis was carried out to determine its association with influenza A/H1N1 infection susceptibility and the LD between the same gene’s SNPs. This analysis included the 11 SNPs evaluated that met the HWE, located in TAP1, PSMB8, and PSMB9 genes, comparing INF-P versus HC groups.
Haplotypes and their frequencies are summarized in Figure 1. Figure 1A shows that the haplotype shaped by (both in PSMB8 gene) is in high LD (r2 = 61). Moreover, according to the frequencies of this haplotype in both groups, we found a significative difference between INF-P vs. HC with GG haplotype (conformed by the common allele of two SNPs, p = 0.038, OR = 1.95, 95% CI = 1.04–3.64).
On the other hand, the haplotype shaped by rs3763365 and rs4148882 is not in high LD (r2 = 49). However, when we compared frequencies of this haplotype in both groups, we found a reduced risk for influenza A/H1N1 infection with AG haplotype (minor alleles of both SNPs, p = 0.046, OR = 0.67, 95% CI = 0.45–0.99).

4. Discussion

Clinical features of influenza A/H1N1 pdm09 infection associated with higher risk (hospitalization, disease severity, and mortality) have been widely investigated. Our study found that patients with influenza A/H1N1 infection were predominantly smoker men with higher BMI. These findings agree with previous reports, where they are added to the mentioned risk factors—pre-existing pulmonary diseases, diabetes, cardiovascular diseases, immunosuppression, and pregnancy are included [3,4,33].
Additionally, previous studies in the Mexican population have reported that intensive care unit admission; low lymphocyte and platelet count; and high CPK, BUN, and LDH are associated with a higher risk for influenza A/H1N1 virus infection mortality [6,7,8], coinciding with our results.
The HWE was assessed for all the SNPs evaluated. Unfortunately, six of the 17 selected SNPs (in TAP2 and TAPBP) did not meet it (p > 0.05). However, we decided to continue with the subsequent analyses due to the SNPs position at the gene structure and their possible biological implications. After allele and genotype frequencies and adjustment results by covariates (sex and BMI), we found that rs2071888/GG genotype in TAPBP is associated with increased risk for influenza A/H1N1 infection. Conversely, rs241433/AG genotype in TAP2 is associated with reduced risk for influenza A/H1N1 infection. Nevertheless, both SNPs did not meet the HWE. Although this could seem a limitation, some authors have described that mestizo populations (such as Mexican) have a rich genetic variability, a product of the years of genetic recombination between ancestral populations (Amerindian, Caucasian, African, and Asian descendants), which could result in some SNPs having a different behavior from those described in the scientific literature [34,35].
TAP complex (shaped by TAP1 and TAP2 subunits) is a necessary transporter for translocating peptides, which were obtained via proteasome degradation, from the cytosol to ER, where HLA class I assembly takes place. Previous studies have demonstrated that in mice and humans, the TAP complex can play a crucial role in the translocation of viral epitopes to ER [36,37], and mutations in TAP1 or TAP2 result in structural changes that can alter antigen recognition and presentation. In this sense, SNPs in TAP2 have been previously associated with a higher risk for ankylosing spondylitis [22] and cervical carcinoma associated with HPV [24]. It has been proposed that the genetic variation in this gene may cause functional alterations, possibly by affecting the ATP-binding capacity and transport efficiency of TAP, as well as the expression and structural stability of the mRNA and protein [38]. The rs241433 has previously only been associated with a higher risk for rheumatoid arthritis [39]. Our work is the first in which this SNP is associated with reduced risk for influenza A/H1N1 infection, and more functional and genetic studies are required to establish its role in the pathogenesis of this infectious disease. In the bioinformatic analysis, the rs241433/A allele has reduced affinity with the myeloid and B-cell development regulatory protein UP1; and on the other hand, the C allele has a high affinity for this protein. Additionally, no differences in predicted secondary mRNA structures were found with the rs241433. The UP1 protein (also known as heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1)) plays a crucial role in the regulation of alternative splicing, and multiple alternatively spliced transcript variants have been found [40]. However, there are no previous studies in literature where the biological and pathological influence has been reported. In silico analysis for rs241433 shows no sites of alternative splicing or miRNA generation. However, functional studies are required to establish the specific biological participation of this SNP.
Tapasin is encoded by the TAPBP gene. This protein plays an essential function in the peptides assembly previously translocated by TAP complex to HLA class I into the ER of cells [17]. It has been reported that some viral products compromise the antigen presentation by interacting with the Tapasin protein [41]. For example, the molluscum contagiosum virus encodes MC80 protein, which generates ER and Tapasin degradation, accompanied by a TAP loss, limiting HLA class I access to cytosolic peptides and a deficient antigen presentation [42]. Moreover, in a Chinese population, the rs9277972 was associated with an increased risk for hepatitis C virus [43]. On the other hand, the rs2071888 (SNP-associated in our study) has only been linked with gastrointestinal stromal tumors [44] and aspirin-exacerbated respiratory disease [45].
There are no previous studies where rs2071888 has been associated with influenza A/H1N1 or another viral infection. The rs2071888 is a missense variant in the TAPBP exon 4, which generates a Threonine > Arginine change, which could alter the polarity of the protein structure, thereby modulating interactions with other molecules such as the HLA class I, leading to an altered antigen presentation and consequently, deficient infected cells clearance by CD8+ T cells, higher viral replication, and disease development [45]. In this sense, we conducted a bioinformatics analysis where we found that the rs2071888/G allele decreases the affinity for a transcriptional factor, affecting the TAPBP expression, probably leading to a deficient antigen presentation. In addition, this SNP modifies the secondary structure of the mRNA, which could alter the protein function. Interestingly, in the in silico analysis, we corroborate that—as we described before—the rs2071888/G allele change from medium size and polar amino acid (threonine) change to large size and basis (arginine), which modifies the polarity of the protein structure [46,47].
We did not find significant associations with the three SNPs associated with influenza A/H1N1 infection in the mortality association analysis, possibly due to the reduced sample size. Other studies have evaluated different SNPs with disease severity (where mortality is included as a severe characteristic). García-Ramírez et al. [11] described the rs361525 in TNF associated with greater disease severity in the same Mexican population. In addition, Chen and colleagues [48] found that the IFITM3 rs12252 SNP is associated with an increased risk of severe and mild influenza in Asian and Caucasian populations.
Furthermore, in the haplotype analysis, we found an increased frequency of the GG haplotype in the INF-P group, shaped by the two common alleles of rs3763365 and rs9276810 (both in the PSMB8 gene), as well as a higher frequency of the AG haplotype, shaped by minor alleles of rs3763365 (PSMB8) and rs4148882 (TAP1) in HC group. Previously, it was established that catalytic proteasome subunits—PSMB8 and PSMB9—were closely linked to the TAP1 and TAP2 genes, suggesting a critical role in processing and antigen presentation [41,49]. To the best of our knowledge, to date, no studies have determined whether the products of different PSMB alleles differ concerning the peptides produced. A previous study demonstrated that an SNP (rs2071543) in PSMB8 is associated with a higher risk for HP [21]. Additionally, a study where rs2071543 was evaluated in cervical cancer showed that a common allele is associated with a higher risk for developing the oncological disease [38], which is in the same sense as our haplotype results, suggesting an essential role of the proteasome in viral clearance and antigen presentation.
Other studies have evaluated single-nucleotide variants in diverse genes. For example, Zúñiga et al. [50] found four polymorphisms associated with severe pneumonia in Influenza A/H1N1, including SNPs in the Fc fragment of immunoglobulin (Ig) G; low-affinity IIA; receptor (FCGR2A), RPA interacting protein (RPAIN); and complement component 1, q subcomponent binding protein (C1QBP) genes. In this sense, Chatzopoulou and colleagues [51] also found positive associations with SNPs in complement-associated genes (FCGR2A, C1QBP, CD55).
Our study is not free of limitations. Firstly, we included a relatively small sample size for both cases and control; however, our study was done with samples from our laboratory’s biobank and exclusively with samples from patients with a confirmed influenza A/H1N1 diagnosis. Secondly, we carried out a mortality analysis; however, we have a reduced sample size. Thirdly, we did not measure the expression for any of the genes evaluated. Despite the associated SNPs not meet HWE, we conducted bioinformatics and in silico analysis, where we found data supporting our genetic-association study.

5. Conclusions

In summary, patients with influenza A/H1N1 infection are predominantly men, obese, and smokers. The rs241433 AC genotype from TAP2 and the AG haplotype (rs3763365 and rs4148882) are associated with a reduced risk for influenza A/H1N1 infection. The GG genotype from TAPBP rs2071888 and the GG haplotype (rs3763365 and rs9276810) are associated with a higher risk for developing influenza A/H1N1 infection in a Mexican mestizo population. On the other hand, these SNPs are not associated with a higher risk of mortality. More association and functional studies and replication in other populations are required to establish these SNPs’ participation in the influenza A/H1N1 infection pathogenesis.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/11/1224/s1: Supplementary Figure S1: Algorithm used for in silico analysis. Supplementary Table S1: Molecular characteristics of the evaluated SNPs. Supplementary Table S2: Noncomparable clinical features of INF-P group. Supplementary Table S3: Hardy–Weinberg equilibrium analysis. Supplementary Table S4: Allele and genotype frequencies from Case and Control groups of the non-significative SNPs evaluated.

Author Contributions

Conceptualization, G.P.-R. and R.F.-V.; data curation, M.A.P.-G. and R.F.-V.; formal analysis, M.A.P.-G., E.A.-O., and R.F.-V.; funding acquisition, A.R.-V. and R.F.-V.; investigation, M.A.P.-G., A.R.-C., and N.E.B.-F.; methodology, A.R.-C. and N.E.B.-F.; project administration, G.P.-R. and R.F.-V.; resources, A.R.-V.; software, E.A.-O.; supervision, A.R.-V. and R.F.-V.; validation, M.A.P.-G., G.P.-R., and R.F.-V.; visualization, A.R.-V. and R.F.-V.; writing—original draft, M.A.P.-G. and R.F.-V.; writing—review & editing, M.A.P.G. and R.F.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Consejo Nacional de Ciencia y Tecnología from México (Grant nos. SALUD-2009-C02-126699 and SALUD-2009-126731 to Ramcés Falfán-Valencia and Alejandra Ramírez-Venegas, respectively). Also supported by the allocated budget to research (RFV-HLA Laboratory) from the Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas (INER).

Acknowledgments

The authors wish to thank Instituto Nacional de Enfermedades Respiratorias Ismael Cosio Villegas for the facilities to carry out this research. The authors acknowledge the support received from physicians and technicians at INER to confirm the study participants’ diagnosis, data acquisition, and clinical care.

Conflicts of Interest

The authors declare no competing interests for this project.

Appendix A

Table A1. Demographic and clinical data of the Survivors and Non-survivors patients from INF-P group.
Table A1. Demographic and clinical data of the Survivors and Non-survivors patients from INF-P group.
VariablesSurvivors = 109Non-Survivors = 19p-Value
Age39 (20–89)44 (30–56)0.93
Sex, male (%)64 (58.72)13 (68.42)0.43
BMI30.9 (17.4–52.7)32.2 (24.4–44.6)0.99
Tobacco smoking (%)49 (44.95)8 (42.11)0.82
Comorbidities
Hypertension (%)17 (15.6)00.07
Asthma (%)9 (8.26)1 (5.26)0.65
Diabetes (%)8 (7.34)2 (10.53)0.64
COPD (%)4 (3.67)01
Symptomatology
Fever (%)85 (77.68)12 (63.16)0.16
Dyspnea (%)79 (72.48)14 (73.68)1
Cough (%)58 (53.21)10 (52.63)0.96
Rhinorrhea (%)38 (34.86)4 (21.05)0.29
Nasal congestion (%)15 (13.76)00.12
Clinical variables
Leukocytes9.5 (2.7–15)10.7 (1.9–10.9)0.87
Platelets208 (261–551)153 (62–233)0.006
Glucose107 (43–341)122 (113–189)0.007
Urea26.5 (8–178)38 (5–63)0.89
Creatinine0.8 (0.25–10.1)0.9 (0.75–1.35)0.25
CPK226 (11–1438)670 (20–7241)0.012
LDH481 (104–1481)1082 (524–1477)<0.001
AST45 (16–184)69 (15–252)0.22
ALT38 (10–174)37 (12–124)0.63
Complications
Pneumonia (%)67 (61.47)9 (47.37)0.25
ICU (%)32 (29.35)14 (73.68)<0.001
ARDS (%)29 (26.61)9 (47.37)0.06
Table A2. Allele and genotype frequencies between Non-survivors and Survivors.
Table A2. Allele and genotype frequencies between Non-survivors and Survivors.
GeneNon-SurvivorsSurvivorsp-ValueOR95% CI
nF (%)nF (%)
TAP2rs241433
Genotypes
AA1157.896458.720.950.970.36–2.60
AC842.114541.281.030.39–2.78
CC0000
19100109100
Alleles
A307917379.360.950.980.42–2.27
C8214520.641.030.44–2.39
TAPBPrs2071888
Genotypes
CC631.581816.670.111
CG736.846762.040.310.09–1.05
GG631.582321.301.280.35–4.64
19100108100
Alleles
C1950103480.791.090.55–2.19
G1950113520.910.46–1.82
TAPBPrs2282851
Genotypes
CC736.842926.850.1301
CT736.846560.190.450.14–1.38
TT526.321412.961.480.40–5.50
19100108100
Alleles
C215512356.940.850.930.47–1.87
T17459343.061.070.53–2.14

Appendix B

Figure A1. mRNA secondary structures for rs241433 (A,B) and rs2071888 (C,D). Blue box: Complete conformational change related to the nitrogenous base; Red box: Nitrogenous base change.
Figure A1. mRNA secondary structures for rs241433 (A,B) and rs2071888 (C,D). Blue box: Complete conformational change related to the nitrogenous base; Red box: Nitrogenous base change.
Viruses 12 01224 g0a1

References

  1. Perez-Padilla, R.; De La Rosa-Zamboni, D.; De Leon, S.P.; Hernandez, M.; Quiñones-Falconi, F.; Bautista, E.; Ramirez-Venegas, A.; Rojas-Serrano, J.; Ormsby, C.E.; Corrales, A.; et al. Pneumonia and Respiratory Failure from Swine-Origin Influenza A (H1N1) in Mexico. N. Engl. J. Med. 2009, 361, 680–689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Echevarría-Zuno, S.; Mejía-Aranguré, J.M.; Mar-Obeso, A.J.; Grajales-Muñiz, C.; Robles-Pérez, E.; González-León, M.; Ortega-Alvarez, M.C.; Gonzalez-Bonilla, C.; Rascón-Pacheco, R.A.; Borja-Aburto, V.H. Infection and death from influenza A H1N1 virus in Mexico: A retrospective analysis. Lancet 2009, 374, 2072–2079. [Google Scholar] [CrossRef]
  3. Ward, K.A.; Spokes, P.J.; McAnulty, J.M. Case-control Study of Risk Factors for Hospitalization Caused by Pandemic (H1N1) 2009. Emerg. Infect. Dis. 2011, 17, 1409–1416. [Google Scholar] [CrossRef] [PubMed]
  4. Fezeu, L.; Julia, C.; Henegar, A.; Bitu, J.; Hu, F.B.; Grobbee, D.E.; Kengne, A.P.; Hercberg, S.; Czernichow, S. Obesity is associated with higher risk of intensive care unit admission and death in influenza A (H1N1) patients: A systematic review and meta-analysis. Obes. Rev. 2011, 12, 653–659. Available online: https://pubmed.ncbi.nlm.nih.gov/21457180/ (accessed on 12 July 2020). [CrossRef] [PubMed]
  5. Chowell, G.; Echevarría-Zuno, S.; Viboud, C.; Simonsen, L.; Tamerius, J.; Miller, M.A.; Borja-Aburto, V.H. Characterizing the Epidemiology of the 2009 Influenza A/H1N1 Pandemic in Mexico. PLoS Med. 2011, 8, e1000436. [Google Scholar] [CrossRef] [PubMed]
  6. Carlos, J.; Salvador, N.; Marcela, L.; Garduño, B.; Molina, H.L.; Serrano, A.T. Factores predictores de defunción en pacientes con neumonía por influenza a H1N1. Med. Interna México 2011, 27. [Google Scholar]
  7. Shi, S.J.; Li, H.; Liu, M.; Liu, Y.M.; Zhou, F.; Liu, B.; Qu, J.X.; Cao, B. Mortality prediction to hospitalized patients with influenza pneumonia: PO2/FiO2combined lymphocyte count is the answer. Clin. Respir. J. 2015, 11, 352–360. [Google Scholar] [CrossRef] [PubMed]
  8. Mata-Marín, L.A.; Mata-Marín, J.A.; Vásquez-Mota, V.C.; Arroyo-Anduiza, C.I.; Gaytán-Martínez, J.; Manjarrez-Téllez, B.; Ochoa-Carrera, L.A.; Sandoval-Ramirez, J.L. Risk factors associated with mortality in patients infected with influenza A/H1N1 in Mexico. BMC Res. Notes 2015, 8, 432. [Google Scholar] [CrossRef] [Green Version]
  9. Bermejo-Martin, J.F.; De Lejarazu, R.O.; Pumarola, T.; Rello, J.; Almansa, R.; Ramírez, P.; Martin-Loeches, I.; Varillas, D.; Gallegos, M.C.; Serón, C.; et al. Th1 and Th17 hypercytokinemia as early host response signature in severe pandemic influenza. Crit. Care 2009, 13, R201. [Google Scholar] [CrossRef] [Green Version]
  10. To, K.K.W.; Hung, I.F.N.; Li, I.W.S.; Lee, K.; Koo, C.; Yan, W.; Liu, R.; Ho, K.; Chu, K.; Watt, C.; et al. Delayed Clearance of Viral Load and Marked Cytokine Activation in Severe Cases of Pandemic H1N1 2009 Influenza Virus Infection. Clin. Infect. Dis. 2010, 50, 850–859. [Google Scholar] [CrossRef] [Green Version]
  11. Garcia-Ramirez, R.A.; Ramirez-Venegas, A.; Quintana-Carrillo, R.; Camarena, Á.E.; Falfán-Valencia, R.; Mejía-Aranguré, J.M. TNF, IL6, and IL1B Polymorphisms Are Associated with Severe Influenza A (H1N1) Virus Infection in the Mexican Population. PLOS ONE 2015, 10, e0144832. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ramirez-Venegas, A.; Gonz lez-Bonilla, C.; Borja-Aburto, V. Pandemic influenza A/H1N1 virus infection and TNF, LTA, IL1B, IL6, IL8, and CCL polymorphisms in Mexican population: A case-control study. BMC Infect. Dis. 2012, 12, 299–307. Available online: http://0-www-biomedcentral-com.brum.beds.ac.uk/1471-2334/12/299 (accessed on 28 July 2020).
  13. Yewdell, J.W.; Bennink, J.R.; Hosaka, Y. Cells process exogenous proteins for recognition by cytotoxic T lymphocytes. Science 1988, 239, 637–640. [Google Scholar] [CrossRef]
  14. Wahl, A.; Schafer, F.; Bardet, W.; Buchli, R.; Air, G.M.; Hildebrand, W.H. HLA class I molecules consistently present internal influenza epitopes. Proc. Natl. Acad. Sci. USA 2009, 106, 540–545. Available online: www.pnas.org/cgi/content/full/ (accessed on 28 July 2020). [CrossRef] [Green Version]
  15. Ortmann, B.; Copeman, J.; Lehner, P.J.; Sadasivan, B.; Herberg, J.A.; Grandea, A.G.; Riddell, S.R.; Tampé, R.; Spies, T.; Trowsdale, J.; et al. A Critical Role for Tapasin in the Assembly and Function of Multimeric MHC Class I-TAP Complexes. Science 1997, 277, 1306–1309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Li, S.; Sjögren, H.-O.; Hellman, U.; Pettersson, R.F.; Wang, P. Cloning and functional characterization of a subunit of the transporter associated with antigen processing. Proc. Natl. Acad. Sci. USA 1997, 94, 8708–8713. Available online: www.pnas.org (accessed on 28 July 2020). [CrossRef] [Green Version]
  17. Rizvi, S.M.; Raghavan, M. Mechanisms of function of tapasin, a critical major histocompatibility complex class I assembly factor. Traffic 2009, 11, 332–347. [Google Scholar] [CrossRef] [PubMed]
  18. Zhou, X.; Momburgu, F.; Liuu, T.; Motalu, U.M.A.; Jondalu, M.; Hämmerlingu, G.J.; Ljunggren, H.-G.; Momburg, F.; Liu, T.; Motal, U.M.A.; et al. Presentation of viral antigens restricted by H-2Kb, Db or Kd in proteasome subunit LMP2-and LMP7-deficient cells. Eur. J. Immunol. 1994, 24, 1863–1868. [Google Scholar] [CrossRef]
  19. Falfán-Valencia, R.; Narayanankutty, A.; Reséndiz-Hernández, J.M.; Pérez-Rubio, G.; Ramírez-Venegas, A.; Nava-Quiroz, K.J.; Bautista-Félix, N.E.; Vargas-Alarcón, G.; Castillejos-López, M.D.J.; Hernández, A. An Increased Frequency in HLA Class I Alleles and Haplotypes Suggests Genetic Susceptibility to Influenza A (H1N1) 2009 Pandemic: A Case-Control Study. J. Immunol. Res. 2018, 2018, 1–12. [Google Scholar] [CrossRef] [PubMed]
  20. Dutta, M.; Dutta, P.; Medhi, S.; Borkakoty, B.; Biswas, D. Polymorphism of HLA class I and class II alleles in influenza A(H1N1)pdm09 virus infected population of Assam, Northeast India. J. Med Virol. 2018, 90, 854–860. [Google Scholar] [CrossRef] [PubMed]
  21. Camarena, A.; Aquino-Gálvez, A.; Falfán-Valencia, R.; Sanchez, G.; Montaño, M.; Ramos, C.; Juárez, A.; De Alba, C.G.; Granados, J.; Selman, M. PSMB8 (LMP7) but not PSMB9 (LMP2) gene polymorphisms are associated to pigeon breeder’s hypersensitivity pneumonitis. Respir. Med. 2010, 104, 889–894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Feng, M.; Yin, B.; Shen, T.; Ma, Q.; Liu, L.; Zheng, J.; Zhao, Y.; Qian, K.; Liu, D. TAP1 and TAP2 polymorphisms associated with ankylosing spondylitis in genetically homogenous Chinese Han population. Hum. Immunol. 2009, 70, 257–261. [Google Scholar] [CrossRef] [PubMed]
  23. Mehta, A.M.; Jordanova, E.S.; Corver, W.E.; Ven Wezel, T.; Uh, H.W.; Kenter, G.G. Single nucleotide polymorphisms in antigen processing machinery component ERAPI significantly associate with clinical outcome in cervical carcinoma. Genes Chromosomes Cancer 2009, 48, 410–418. Available online: https://pubmed.ncbi.nlm.nih.gov/19202550/ (accessed on 28 July 2020). [CrossRef] [PubMed]
  24. Deshpande, A.; Wheeler, C.M.; Hunt, W.C.; Peyton, C.L.; White, P.S.; Valdez, Y.E.; Nolan, J.P. Variation in HLA Class I Antigen-Processing Genes and Susceptibility to Human Papillomavirus Type 16–Associated Cervical Cancer. J. Infect. Dis. 2008, 197, 371–381. [Google Scholar] [CrossRef]
  25. Little, J.; Higgins, J.P.T.; Ioannidis, J.P.; Moher, D.; Gagnon, F.; Von Elm, E.; Khoury, M.J.; Cohen, B.; Davey-Smith, G.; Grimshaw, J.; et al. STrengthening the REporting of Genetic Association studies (STREGA)—An extension of the STROBE statement. Eur. J. Clin. Investig. 2009, 39, 247–266. [Google Scholar] [CrossRef] [Green Version]
  26. Ensembl Genome Browser 89. Available online: http://www.ensembl.org/index.html (accessed on 4 August 2017).
  27. Ward, L.D.; Kellis, M. HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [Green Version]
  28. SNPStats: A Web Tool for the Analysis of Association Studies—PubMed. Available online: https://pubmed.ncbi.nlm.nih.gov/16720584/ (accessed on 23 June 2020).
  29. Haploview: Analysis and Visualization of LD and Haplotype Maps—PubMed. Available online: https://pubmed.ncbi.nlm.nih.gov/15297300/ (accessed on 23 June 2020).
  30. Gabriel, S.; Schaffner, S.F.; Nguyen, H.; Moore, J.M.; Roy, J.; Blumenstiel, B.; Higgins, J.; DeFelice, M.; Lochner, A.; Faggart, M.; et al. The Structure of Haplotype Blocks in the Human Genome. Science 2002, 296, 2225–2229. [Google Scholar] [CrossRef] [Green Version]
  31. Dean, A.G.; Arner, T.G.; Sunki, G.G.; Friedman, R.; Lantinga, M.; Sangam, S.; Zubieta, J.C.; Sullivan, K.M.; Brendel, K.A.; Gao, Z.; et al. Epi Info, A Database and Statistics Program for Public Health Professionals; CDC: Atlanta, GA, USA, 2004.
  32. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; De Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  33. Lenzi, L.; De Mello, Â.M.; Da Silva, L.R.; Grochocki, M.H.C.; Pontarolo, R. Influenza pandêmica A (H1N1) 2009: fatores de risco para o internamento. J. Bras. Pneumol. 2012, 38, 57–65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ambrocio-Ortiz, E.; Pérez-Rubio, G.; Ramírez-Venegas, A.; Hernández-Zenteno, R.; Del Angel-Pablo, A.D.; Pérez-Rodríguez, M.E.; Salazar, A.M.; Abarca-Rojano, E.; Falfán-Valencia, R. Effect of SNPs in HSP Family Genes, Variation in the mRNA and Intracellular Hsp Levels in COPD Secondary to Tobacco Smoking and Biomass-Burning Smoke. Front. Genet. 2020, 10, 1307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Namipashaki, A.; Razaghi-Moghadam, Z.; Ansari-Pour, N. The essentiality of reporting Hardy-Weinberg equilibrium calculations in population-based genetic association studies. Cell J. Royan Inst. 2015, 17, 187–192. [Google Scholar]
  36. Braud, V.M.; McMichael, A.J.; Cerundolo, V. Differential processing of influenza nucleoprotein in human and mouse cells. Eur. J. Immunol. 1998, 28, 625–635. Available online: https://0-onlinelibrary-wiley-com.brum.beds.ac.uk/doi/full/10.1002/%28SICI%291521-4141%28199802%2928%3A02%3C625%3A%3AAID-IMMU625%3E3.0.CO%3B2-I (accessed on 29 July 2020). [CrossRef]
  37. Ho, A.W.S.; Prabhu, N.; Betts, R.J.; Ge, M.Q.; Dai, X.; Hutchinson, P.E.; Lew, F.C.; Wong, K.L.; Hanson, B.J.; Macary, P.A.; et al. Lung CD103 + Dendritic Cells Efficiently Transport Influenza Virus to the Lymph Node and Load Viral Antigen onto MHC Class I for Presentation to CD8 T Cells. J. Immunol. 2011, 187, 6011–6021. [Google Scholar] [CrossRef] [Green Version]
  38. Mehta, A.M.; Jordanova, E.S.; Van Wezel, T.; Uh, H.-W.; Corver, W.E.; Kwappenberg, K.M.C.; Verduijn, W.; Kenter, G.G.; Van Der Burg, S.H.; Fleuren, G.J. Genetic variation of antigen processing machinery components and association with cervical carcinoma. Genes Chromosomes Cancer 2007, 46, 577–586. [Google Scholar] [CrossRef] [PubMed]
  39. Lee, H.-S.; Lee, A.T.; Criswell, L.A.; Seldin, M.F.; Amos, C.I.; Carulli, J.P.; Navarrete, C.; Remmers, E.F.; Kastner, D.L.; Plenge, R.M.; et al. Several Regions in the Major Histocompatibility Complex Confer Risk for Anti-CCP-Antibody Positive Rheumatoid Arthritis, Independent of the DRB1 Locus. Mol. Med. 2008, 14, 293–300. [Google Scholar] [CrossRef]
  40. HNRNPA1 Heterogeneous Nuclear Ribonucleoprotein A1 [Homo Sapiens (Human)]—Gene—NCBI. [cited]. Available online: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/gene/3178 (accessed on 9 September 2020).
  41. Kelly, A.; Trowsdale, J. Genetics of antigen processing and presentation. Immunogenetics 2018, 71, 161–170. [Google Scholar] [CrossRef] [Green Version]
  42. Harvey, I.B.; Wang, X.; Fremont, D.H. Molluscum contagiosum virus MC80 sabotages MHC-I antigen presentation by targeting tapasin for ER-associated degradation. PLoS Pathog. 2019, 15, e1007711. [Google Scholar] [CrossRef] [PubMed]
  43. Huang, P.; Zhang, Y.; Lu, X.; Xu, Y.; Wang, J.; Zhang, Y.; Yu, R.; Su, J. Association of Polymorphisms in HLA Antigen Presentation-Related Genes with the Outcomes of HCV Infection. PLoS ONE 2015, 10, e0123513. [Google Scholar] [CrossRef] [Green Version]
  44. O’Brien, K.M.; Orlow, I.; Antonescu, C.R.; Ballman, K.V.; McCall, L.; DeMatteo, R.P.; Engel, L.S. Gastrointestinal stromal tumors: A case-only analysis of single nucleotide polymorphisms and somatic mutations. Clin. Sarcoma Res. 2013, 3, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Cho, S.-H.; Park, J.-S.; Park, B.L.; Bae, D.-J.; Uh, S.-T.; Kim, M.-K.; Choi, I.S.; Shin, H.D.; Park, C.-S. Association analysis of tapasin polymorphisms with aspirin-exacerbated respiratory disease in asthmatics. Pharm. Genom. 2013, 23, 341–348. [Google Scholar] [CrossRef]
  46. TAPBP—Tapasin Precursor—Homo Sapiens (Human)—TAPBP Gene & Protein. Available online: https://www.uniprot.org/uniprot/O15533 (accessed on 9 September 2020).
  47. UniProtKB/SwissProt Variant VAR_010253. Available online: https://web.expasy.org/variant_pages/VAR_010253.html (accessed on 9 September 2020).
  48. Chen, T.; Xiao, M.-F.; Yang, J.; Chen, Y.K.; Bai, T.; Tang, X.J.; Shu, Y. Association between rs12252 and influenza susceptibility and severity: An updated meta-analysis. Epidemiol. Infect. 2018, 147, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Glynne, R.; Powis, S.H.; Beck, S.; Kelly, A.; Kerr, L.-A.; Trowsdale, J. A proteasome-related gene between the two ABC transporter loci in the class II region of the human MHC. Nature 1991, 353, 357–360. [Google Scholar] [CrossRef] [PubMed]
  50. Zúñiga, J.; Buendia-Roldan, I.; Zhao, Y.; Jiménez, L.; Torres, D.; Romo, J.; Ramirez, G.; Cruz, A.; Vargas-Alarcón, G.; Sheu, C.-C.; et al. Genetic variants associated with severe pneumonia in A/H1N1 influenza infection. Eur. Respir. J. 2011, 39, 604–610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Chatzopoulou, F.; Gioula, G.; Kioumis, I.; Chatzidimitriou, D.; Exindari, M. Identification of complement-related host genetic risk factors associated with influenza A(H1N1)pdm09 outcome: Challenges ahead. Med. Microbiol. Immunol. 2019, 208, 631–640. [Google Scholar] [CrossRef]
Figure 1. Haplotype analysis. (A) Haplotype conformed by the common allele of rs3763365 and rs9276810; (B) Haplotype conformed by minor alleles of rsr3763365 and rs41488882. INF-P—Patients with influenza A/H1N1 infection; HC—Healthy contacts. p-value < 0.05 was considered as significative.
Figure 1. Haplotype analysis. (A) Haplotype conformed by the common allele of rs3763365 and rs9276810; (B) Haplotype conformed by minor alleles of rsr3763365 and rs41488882. INF-P—Patients with influenza A/H1N1 infection; HC—Healthy contacts. p-value < 0.05 was considered as significative.
Viruses 12 01224 g001
Table 1. Demographic and clinical data of the influenza A/H1N1 infection (INF-P) and healthy contact (HC) groups.
Table 1. Demographic and clinical data of the influenza A/H1N1 infection (INF-P) and healthy contact (HC) groups.
VariablesINF-P = 128HC = 111p-Value
Age41 (19–86)44 (19–82)0.986
Sex, male (%)77 (60.16)46 (41.44)0.003
BMI30.5 (17.4–52.8)26.8 (16–44.8)<0.001
Tobacco smoking, %57 (44.53)39 (35.14)0.139
All values are expressed as median and minimum–maximum values. We used the Mann–Whitney U test and Fisher Exact test. p-value < 0.05 was considered as significant.
Table 2. Allele and genotype frequencies between cases and controls.
Table 2. Allele and genotype frequencies between cases and controls.
GeneINF-PHCp-Valuep-Value CorrectedOR95% CI
nF (%)nF (%)
TAP2rs241433
Genotypes
AA7558.592825.23<0.0001<0.00014.192.41–7.30
AC5341.418374.770.240.14–0.42
CC0000
128100111100
Alleles
A2037913962.61<0.0001<0.00012.291.52–3.44
C53218337.390.440.29–0.66
TAPBPrs2071888
Genotypes
CC2418.902825.450.00080.0141
CG7458.277669.091.140.60–2.14
GG2922.8365.455.642.00–15.86
127100110100
Alleles
C12248.03132600.0090.1530.620.43–0.89
G13251.9788401.621.13–2.34
TAPBPrs2282851
Genotypes
CC3629.031513.640.0170.2891
CT6955.657568.180.380.19–0.76
TT1915.322018.180.400.17–0.94
124100110100
Alleles
C1415710547.730.0480.8161.441.00–2.08
T1074311552.270.690.40–0.99
INF-P—Patients with influenza A/H1N1 infection; HC—Healthy contacts. p-value < 0.05 was considered as significant.
Table 3. Logistic regression analysis for single nucleotide polymorphisms (SNPs) associated with influenza A virus (IAV) infection.
Table 3. Logistic regression analysis for single nucleotide polymorphisms (SNPs) associated with influenza A virus (IAV) infection.
Alleles
ChrGeneSNPA1Testp-ValueOR95% CI
6TAP2rs241433CAdd<0.00010.2410.133–0.434
Sex0.039791.8281.029–3.247
BMI<0.00011.1021.05–1.157
TAPBPrs2071888GAdd0.00951.8911.168–3.061
Sex0.00562.1841.257–3.794
BMI0.00011.0971.046–1.15
Genotypes
ChrGeneSNPA1Testp-ValueOR95% CI
6TAP2rs241433CCAddNANANA
DomdevNANANA
SexNANANA
BMINANANA
Geno_2DFNANANA
TAPBPrs2071888GGAdd0.0052.1841.274–3.745
Domdev0.0750.56120.297–1.059
Sex0.0092.1021.205–3.668
BMI0.00021.0931.043–1.145
Geno_2DF0.013NANA
Chr—Chromosome; SNP—Single Nucleotide Polymorphism; BP—Base pair location; A1—minor allele for each SNP; OR—Odds ratio; 95% CI—Confidence interval; Add—Additive effect model; Domdev—Dominance deviation; BMI—Body mass index; Geno_2DF—Genotypic 2 degree-freedom. p-value < 0.05 was considered as significant.
Table 4. Sequences and affinity for different proteins.
Table 4. Sequences and affinity for different proteins.
SNP AlleleSequenceAbsolute Affinity ValueProtein
rs241433/CGGATATACAGTCCCTTCTCCTACCATACAG12.4UP1
rs241433/AGGATATACAGTCCATTCTCCTACCATACAG0.4
rs2071888/GTTGAACTGTAGGCAGCC0.8HNF4
rs2071888/CTTGAACTCTAGGCAGCC10.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ponce-Gallegos, M.A.; Ruiz-Celis, A.; Ambrocio-Ortiz, E.; Pérez-Rubio, G.; Ramírez-Venegas, A.; Bautista-Félix, N.E.; Falfán-Valencia, R. Polymorphisms in Processing and Antigen Presentation-Related Genes and Their Association with Host Susceptibility to Influenza A/H1N1 2009 Pandemic in a Mexican Mestizo Population. Viruses 2020, 12, 1224. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111224

AMA Style

Ponce-Gallegos MA, Ruiz-Celis A, Ambrocio-Ortiz E, Pérez-Rubio G, Ramírez-Venegas A, Bautista-Félix NE, Falfán-Valencia R. Polymorphisms in Processing and Antigen Presentation-Related Genes and Their Association with Host Susceptibility to Influenza A/H1N1 2009 Pandemic in a Mexican Mestizo Population. Viruses. 2020; 12(11):1224. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111224

Chicago/Turabian Style

Ponce-Gallegos, Marco Antonio, Aseneth Ruiz-Celis, Enrique Ambrocio-Ortiz, Gloria Pérez-Rubio, Alejandra Ramírez-Venegas, Nora E. Bautista-Félix, and Ramcés Falfán-Valencia. 2020. "Polymorphisms in Processing and Antigen Presentation-Related Genes and Their Association with Host Susceptibility to Influenza A/H1N1 2009 Pandemic in a Mexican Mestizo Population" Viruses 12, no. 11: 1224. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111224

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop