Polymorphic variants in exon 8 at the 3' UTR of the HLA-G gene are associated with septic shock in critically ill patients

Genome-wide identification and functional analysis of dysregulated alternative splicing profiles in sepsis

BACKGROUND

An increasing body of evidence now shows that the long-term mortality of patients with sepsis are associated with various sepsis-related immune cell defects. Alternative splicing (AS), as a sepsis-related immune cell defect, is considered as a potential immunomodulatory therapy target to improve patient outcomes. However, our understanding of the role AS plays in sepsis is currently insufficient.

AIM

This study investigated possible associations between AS and the gene regulatory networks affecting immune cells. We also investigated apoptosis and AS functionality in sepsis pathophysiology.

METHODS

In this study, we assessed publicly available mRNA-seq data that was obtained from the NCBI GEO dataset (GSE154918), which included a healthy group (HLTY), a mild infection group (INF1), asepsis group (Seps), and a septic shock group (Shock). A total of 79 samples (excluding significant outliers) were identified by a poly-A capture method to generate RNA-seq data. The variable splicing events and highly correlated RNA binding protein (RBP) genes in each group were then systematically analyzed.

RESULTS

For the first time, we used systematic RNA-seq analysis of sepsis-related AS and identified 1505 variable AS events that differed significantly (p <= 0.01) across the four groups. In the sepsis group, the genes related to significant AS events, such as, SHISA5 and IFI27, were mostly enriched in the cell apoptosis pathway. Furthermore, we identified differential splicing patterns within each of the four groups. Significant differences in the expression of RNA Binding Protein(RBP) genes were observed between the control group and the sepsis group. RBP gene expression was highly correlated with variant splicing events in sepsis, as determined by co-expression analysis; The expression of DDX24, CBFA2T2, NOP, ILF3, DNMT1, FTO, PPRC1, NOLC1 RBPs were significant reduced in sepsis compared to the healthy group. Finally, we constructed an RBP-AS functional network.

CONCLUSION

Analysis indicated that the RBP-AS functional network serves as a critical post-transcriptional mechanism that regulates the development of sepsis. AS dysregulation is associated with alterations in the regulatory gene expression network that is involved in sepsis. Therefore, the RBP-AS expression network could be useful in refining biomarker predictions in the development of new therapeutic targets for the pathogenesis of sepsis.

Evidence for Epistatic Interaction between HLA-G and LILRB1 in the Pathogenesis of Nonsegmental Vitiligo

Vitiligo is the most frequent cause of depigmentation worldwide. Genetic association studies have discovered about 50 loci associated with disease, many with immunological functions. Among them is HLA-G, which modulates immunity by interacting with specific inhibitory receptors, mainly LILRB1 and LILRB2. Here we investigated the LILRB1 and LILRB2 association with vitiligo risk and evaluated the possible role of interactions between HLA-G and its receptors in this pathogenesis. We tested the association of the polymorphisms of HLA-G, LILRB1, and LILRB2 with vitiligo using logistic regression along with adjustment by ancestry. Further, methods based on the multifactor dimensionality reduction (MDR) approach (MDR v.3.0.2, GMDR v.0.9, and MB-MDR) were used to detect potential epistatic interactions between polymorphisms from the three genes. An interaction involving rs9380142 and rs2114511 polymorphisms was identified by all methods used. The polymorphism rs9380142 is an HLA-G 3'UTR variant (+3187) with a well-established role in mRNA stability. The polymorphism rs2114511 is located in the exonic region of LILRB1. Although no association involving this SNP has been reported, ChIP-Seq experiments have identified this position as an EBF1 binding site. These results highlight the role of an epistatic interaction between HLA-G and LILRB1 in vitiligo pathogenesis.

Non-classical human leukocyte antigen class I in Tunisian children with autism

Autism spectrum disorders (ASD) are one of the most common childhood morbidities characterized by deficits in communication and social skills. Increasing evidence has suggested associations between immune genes located in the human leukocyte antigen (HLA) complex and etiology of autism.

In this study, we investigated whether the non-classical class I HLA-G, -E, and -F polymorphisms are associated with genetic predisposition to autism in Tunisia. We aimed to find a correlation between HLA-G genotypes and soluble HLA-G (sHLA-G) levels. We have analyzed the HLA-G, -E, and -F genotypes of 15 autistic children and their parents. DNA typing of HLA class I genes was performed using PCR-SSP and PCR-RFLP methods. Also, we evaluated the serum levels of HLA-G (1 and 5) by a validated ELISA technique in autistic probands and their parents.

No association was found between any polymorphism and autism in the study subjects. Additionally, we found no correlation between sHLA-G1 and sHLA-G5 and autism. Also, no significant difference in sHLA-G testing in parents and offspring was found. However, parents carrying [GG] genotype presented a higher sHLA-G levels than those carrying ([CC]+[GC]) genotypes (p = 0.037).

From this preliminary study, we conclude that the investigated polymorphisms of HLA-G, -E, and -F genes did not lead to autism susceptibility in Tunisian children. However, the CGTIGA haplotype was found to be associated with the disease.

Genetic association between HLA-G 14-bp polymorphism and diseases: A systematic review and meta-analysis

BACKGROUND

HLA-G is an immune checkpoint molecule. Since a differential molecule expression has been reported even for healthy individuals, many studies have focused on polymorphisms at HLA-G regulatory regions, particularly the 3' untranslated region (3'UTR). The presence/absence of a 14-bp sequence was the first polymorphism described and it is the most studied in association between HLA-G and disorders.

METHODS

In this study, we performed a systematic review and meta-analysis of all association studies published regarding the HLA-G 14-bp.

RESULTS

We verified association between 14-bp alleles and diseases in the following situations: (1) presence of 14-bp (insertion) conferred susceptibility to preeclampsia (child alleles evaluated) and systemic lupus erythematosus (OR = 1.42; 95%CI = 1.04-1.93; p = 0.026 and OR = 1.13; 95%CI = 1.01-1.27, p = 0.028); (2) 14-bp absence (deletion) was associated with increased risk to breast cancer (OR = 1.23; 95%CI = 1.06-1.43; p = 0.006) and human Cytomegalovirus infection (OR = 2.06; 95%CI = 1.60-2.64; p < 0.0001); and (3) a risk association was observed between the group of reproductive disorders and the 14-bp insertion (OR = 1.12; 95%CI = 1.01-1.24; p = 0.034).

CONCLUSIONS

Considering that others 14-bp associations were inconclusive and that other variation sites observed at HLA-G 3'UTR exhibit a proven role on post-transcriptional regulation of HLA-G expression, the complete 3'UTR segment should be analyzed in terms of disease susceptibility, instead of a single polymorphism.

Association of HLA-G 3'UTR polymorphisms and haplotypes with severe sepsis in a Brazilian population

BACKGROUND

The human leukocyte antigen G (HLA-G) is a molecule involved in immune system modulation, acting in the maintenance of a state of immune tolerance. Some polymorphisms in the HLA-G gene 3' untranslated region (3'UTR) were associated to distinct levels of HLA-G expression and to sepsis development. In the present study, haplotypes and polymorphisms of the HLA-G 3'UTR were analyzed in Brazilian septic patients.

METHODS

The HLA-G 3'UTR was amplified by PCR, sequenced and eight polymorphisms were genotyped (the 14bp insertion/deletion, +3003T/C, +3010C/G, +3027A/C, +3035C/T, +3142G/C, +3187A/G and+3196C/G) in DNA samples from septic patients (with severe sepsis or septic shock) and controls. The haplotypes were inferred and association tests were performed through Chi square test and binary logistic regression.

RESULTS

The+3027AC genotype was associated asa risk factor to sepsis development (OR 3.17, P_Bonferroni 0.048). Further, the presence of the UTR-7 haplotype (OR 2.97, P_Bonferroni 0.018), and of 14bp-Ins_+3142G_+3187A haplotype (OR 2.39, P_Bonferroni 0.045) were associated with sepsis, conferring susceptibility.

CONCLUSION

Our data confirm an important role of HLA-G 3'UTR polymorphisms in the development of severe forms of sepsis (severe sepsis and septic shock). The genotyping of HLA-G genetic variants and haplotypes could be useful as a prediction tool of increased risk to severe sepsis.

Correlations between ACE single nucleotide polymorphisms and prognosis of patients with septic shock

The aim of the present study is to investigate association between septic shock (SS) and angiotensin I-converting enzyme (ACE) single nucleotide polymorphisms (SNPs). From October 2009 to December 2016, 238 SS patients and 242 healthy individuals were selected for our study. ACE activity was detected, ACE rs4291 and rs4646994 polymorphisms were detected using PCR-restriction fragment length polymorphism (PCR-RFLP). The Kaplan–Meier survival curve was employed to evaluate the association between ACE SNPs and patients’ survival and univariate and multivariate analyses to estimate risk factors for SS. ACE activity in the case group was increased in comparison with the control group. Allele and genotype frequencies of rs4291 and rs4646994 were different between the case and control groups. The TT genotype frequency of the rs4291 polymorphisms and the DD genotype of the rs4646994 polymorphisms of the case group were higher than those in the control group. The AT and TT genotypes indicated a significant elevation of ACE activity than the AA genotype, while a significant decline was found in the DI and II genotypes in comparison with the DI genotype. Patients with TT or DD genotypes had increased fatality rate within 7 and 30 days when compared with those with non-TT or non-DD genotypes. Lower sepsis-related organ failure assessment (SOFA) scores, rs4291, serum ACE and rs4646994 were all considered as risky factors for SS patients. The study demonstrates that TT genotype of rs4291 or DD genotype of rs4646994 may be indicative of a higher risk of SS and a poorer prognosis in SS patients.

Modulation of LILRB2 protein and mRNA expressions in septic shock patients and after ex vivo lipopolysaccharide stimulation

Septic patients develop immune dysfunctions, the intensities and durations of which are associated with deleterious outcomes. LILRB2 (leukocyte immunoglobulin-like receptors subfamily B, member 2), an inhibitory member of the LILR family of receptors, is known for its immunoregulatory properties. In a microarray study, we identified LILRB2 as an upregulated gene in septic shock patients. On monocytes primed with LPS ex vivo, LILRB2 mRNA and protein expressions were dose-dependently downregulated and subsequently highly upregulated versus non-stimulated cells. This is concordant with clinical data, since both LILRB2 mRNA and protein expressions were significantly increased in septic shock patients at day 3. In a cohort of more than 700 patients, only after septic shock were LILRB2 mRNA levels increased compared with non-infected or less severely infected patients. This was preceded by a phase of downregulated mRNA expression during the first hours after septic shock. Interestingly, the intensity of this decrease was associated with increased risk of death after septic shock. LILRB2 protein and mRNA expressions are deregulated on monocytes after septic shock and this can be reproduced ex vivo after LPS challenge. Considering LILRB2 inhibitory properties, we can hypothesize that LILRB2 may participate in the altered immune response after septic shock.

Reduced frequency of two activating KIR genes in patients with sepsis

Natural killer (NK) cell activity is regulated by activating and inhibitory signals transduced by killer cell immunoglobulin-like receptors (KIR). Diversity in KIR gene repertoire among individuals may affect disease outcome. Sepsis development and severity may be influenced by genetic factors affecting the immune response. Here, we examined sixteen KIR genes and their human leucocyte antigen (HLA) class I ligands in critical patients, aiming to identify patterns that could be associated with sepsis. Male and female patients (ages ranging between 14 and 94years-old) were included. DNA samples from 211 patients with sepsis and 60 controls (critical care patients with no sepsis) collected between 2004 and 2010 were included and genotyped for KIR genes using the polymerase chain reaction method with sequence-specific oligonucleotide (PCR-SSO), and for HLA genes using the polymerase chain reaction method with sequence-specific primers (PCR-SSP). The frequencies of activating KIR2DS1 and KIR3DS1 in sepsis patients when compared to controls were 41.23% versus 55.00% and 36.49% versus 51.67% (p=0.077 and 0.037 respectively before Bonferroni correction). These results indicate that activating KIR genes 2DS1 and 3DS1 may more prevalent in critical patients without sepsis than in patients with sepsis, suggesting a potential protective role of activating KIR genes in sepsis.

Fc Gamma Receptor IIA (CD32A) R131 Polymorphism as a Marker of Genetic Susceptibility to Sepsis

Sepsis is a devastating disease that can affect humans at any time between neonates and the elderly and is associated with mortality rates that range from 30 to 80%. Despite intensive efforts, its treatment has remained the same over the last few decades. Fc receptors regulate multiple immune responses and have been investigated in diverse complex diseases. FcγRIIA (CD32A) is an immunoreceptor, tyrosine-based activation motif-bearing receptor that binds immunoglobulin G and C-reactive protein, important opsonins in host defense. We conducted a study of 702 patients (184 healthy individuals, 171 non-infected critically ill patients, and 347 sepsis patients) to investigate if genetic polymorphisms in the CD32A coding region affect the risk of septic shock. All individuals were genotyped for a variant at position 131 of the FcγRIIA gene. We found that allele G, associated with the R131 genotype, was significantly more frequent in septic patients than in the other groups (p = 0.05). Our data indicate that FcγRIIA genotyping can be used as a marker of genetic susceptibility to sepsis.

Immunogenetics of prostate cancer and benign hyperplasia--the potential use of an HLA-G variant as a tag SNP for prostate cancer risk

Human leukocyte antigen G (HLA-G) is an immunomodulatory molecule with important roles both physiologically as well as an escape mechanism of cancer cells. In this study, we evaluated the impact of eight polymorphisms at the 3' untranslated region (3'UTR) of the HLA-G gene in the development of prostate cancer (PCa) and benign prostatic hyperplasia (BPH). A total of 468 DNA samples of Brazilian men predominantly Euro-descendant with PCa (N = 187), BPH (N = 152) and healthy control individuals (N = 129) were evaluated. The HLA-G 3'UTR region was amplified by polymerase chain reaction (PCR), sequenced and genotyped to identify the 14 bp insertion/deletion (rs371194629), +3003T/C (rs1707), +3010C/G (rs1710), +3027A/C (rs17179101), +3035C/T (rs17179108), +3142G/C (rs1063320), +3187A/G (rs9380142) and +3196C/G (rs1610696) polymorphisms. Regression logistic and chi-square tests were performed to verify the influence of single nucleotide polymorphisms (SNPs) in PCa and/or BPH susceptibility, as well as in PCa progression (clinicopathological status). Our data showed the UTR-4 haplotype as a risk factor to PCa in comparison with control [odds ratio (OR) 2.35, 95% confidence interval (CI) 1.39-3.96, P adjusted = 0.003) and BPH groups (OR 1.82, 95% CI 1.15-2.86, P adjusted = 0.030). Further, the 'non-14bp Ins_ + 3142G_+3187A' haplotype (OR 1.56, 95% CI 1.10-2.20, P adjusted = 0.036), the +3003CT genotype (OR 4.44, 95% CI 1.33-4.50, P adjusted = 0.032) and the +3003C allele (OR 2.33, 95% CI 1.38-3.92, P adjusted = 0.016) also conferred susceptibility to PCa. Our data suggest an important influence of HLA-G 3'UTR polymorphisms in PCa susceptibility and support the use of the +3003 variant as a tag SNP for PCa risk.

The role of genetics and antibodies in sepsis

During the course of sepsis when immunosuppression predominates, the concentrations of circulating immunoglobulins (IGs) are decreased and this is associated with adverse outcomes. The production of IGs as response to invasive bacterial pathogens takes place through a complex pathway starting from the recognition of the antigen (Ag) by innate immune cells that process and present Ags to T cells. The orchestration of T-helper (Th) lymphocyte responses directs specific B cells and ends with the production of IGs by plasma cells. All molecules implicated in this process are encoded by genes bearing single nucleotide polymorphisms (SNPs). Meta-analysis of case-control studies have shown that the carriage of minor frequency SNPs of CD14, TLR2 and TNF is associated with increased sepsis risk. The ambiguity of results of clinical trials studying the clinical efficacy of exogenous IG administration in sepsis suggests that efficacy of treatment should be considered after adjustment for SNPs of all implicated genes in the pathway of IG production.

Nonclassical human leukocyte antigen (HLA-G, HLA-E, and HLA-F) in coronary artery disease

AIMS

Several evidences suggest the association between the evolution of coronary artery disease (CAD) and the development of coronary syndrome that is often associated with disrupted plaque and partial or complete thrombosis of the related artery. Because of the inflammatory nature of CAD, we investigated the human leukocyte antigen (HLA)-G, HLA-E, and HLA-F genetic polymorphisms within CAD patients and evaluated their potential association with this disease in Tunisian population.

METHODS

Different polymorphisms in HLA-G (14-bp Insertion/Deletion, +3142C/G), HLA-E (HLA-E*01:01/01:03 A/G), HLA-F (HLA-F*01:02 T/C, 01:03 C/T, 01:04 A/C) genes were typed using different laboratory techniques in a cohort of 89 CAD patients and 84 controls.

RESULTS

A significant association was reported between the HLA-G +3142 G allele (OR=1.64, 95% CI=1.05-2.56, p=0.02) and increased risk of CAD. No association was found for the other studied polymorphisms. When we considered the haplotypes, we found TDELCA and TDELGG haplotypes associated to CAD with p=0.008 and p=0.030, respectively, suggesting the potential interaction between HLA-G and HLA-E genes.

CONCLUSIONS

Our findings indicated that the HLA-G +3142C/G polymorphism and TDELCA and TDELGG haplotypes can harbour a reliable diagnosis value for the risk of CAD development suggesting that HLA-G, -E and -F molecules might be involved in the pathogenesis of the disease. However, further studies are necessary to confirm our results.

A combination of single nucleotide polymorphisms in the 3'untranslated region of HLA-G is associated with preeclampsia

Reduced expression of human leukocyte antigen-G (HLA-G) has been linked to onset of preeclampsia. Associations have also been reported between preeclampsia and single nucleotide polymorphisms (SNP) in the 3'-untranslated region (UTR) of the HLA-G gene. However, there are conflicting results between studies. This studied examined whether a SNP, by itself or in combination with other SNPs, in the 3'UTR of the HLA-G gene is associated with an increased risk of preeclampsia. Placenta samples were obtained from 47 preeclamptic and 68 control cases. DNA was extracted, and the 3'UTR was sequenced and analyzed for nine polymorphisms using different genetic models of inheritance. Four of these polymorphisms have never been analyzed for an association with preeclampsia. Disputing existing reports, preeclamptic cases were suggestively associated with a G/G-genotype at SNP +3187 (p<0.05). Several SNP combinations were more prevalent in preeclampsia cases. Following corrections for multiple testing, one SNP combination (+3027C/C and +3187G/G) was significantly more prevalent in preeclampsia cases using co-dominant, additive, and dominant models (p<0.001). Taken together with the current literature, the data suggests that HLA-G 3'UTR SNP-pair associations, and not individual SNPs, could be useful in a predictive test for the susceptibility to preeclampsia.

Insights into HLA-G Genetics Provided by Worldwide Haplotype Diversity

Human leukocyte antigen G (HLA-G) belongs to the family of non-classical HLA class I genes, located within the major histocompatibility complex (MHC). HLA-G has been the target of most recent research regarding the function of class I non-classical genes. The main features that distinguish HLA-G from classical class I genes are (a) limited protein variability, (b) alternative splicing generating several membrane bound and soluble isoforms, (c) short cytoplasmic tail, (d) modulation of immune response (immune tolerance), and (e) restricted expression to certain tissues. In the present work, we describe the HLA-G gene structure and address the HLA-G variability and haplotype diversity among several populations around the world, considering each of its major segments [promoter, coding, and 3′ untranslated region (UTR)]. For this purpose, we developed a pipeline to reevaluate the 1000Genomes data and recover miscalled or missing genotypes and haplotypes. It became clear that the overall structure of the HLA-G molecule has been maintained during the evolutionary process and that most of the variation sites found in the HLA-G coding region are either coding synonymous or intronic mutations. In addition, only a few frequent and divergent extended haplotypes are found when the promoter, coding, and 3′UTRs are evaluated together. The divergence is particularly evident for the regulatory regions. The population comparisons confirmed that most of the HLA-G variability has originated before human dispersion from Africa and that the allele and haplotype frequencies have probably been shaped by strong selective pressures.

Energy crisis: the role of oxidative phosphorylation in acute inflammation and sepsis

Mitochondrial dysfunction is increasingly recognized as an accomplice in most of the common human diseases including cancer, neurodegeneration, diabetes, ischemia/reperfusion injury as seen in myocardial infarction and stroke, and sepsis. Inflammatory conditions, both acute and chronic, have recently been shown to affect mitochondrial function. We here discuss the role of oxidative phosphorylation (OxPhos), focusing on acute inflammatory conditions, in particular sepsis and experimental sepsis models. We discuss mitochondrial alterations, specifically the suppression of oxidative metabolism and the role of mitochondrial reactive oxygen species in disease pathology. Several signaling pathways including metabolic, proliferative, and cytokine signaling affect mitochondrial function and appear to be important in inflammatory disease conditions. Cytochrome c oxidase (COX) and cytochrome c, the latter of which plays a central role in apoptosis in addition to mitochondrial respiration, serve as examples for the entire OxPhos system since they have been studied in more detail with respect to cell signaling. We propose a model in which inflammatory signaling leads to changes in the phosphorylation state of mitochondrial proteins, including Tyr304 phosphorylation of COX catalytic subunit I. This results in an inhibition of OxPhos, a reduction of the mitochondrial membrane potential, and consequently a lack of energy, which can cause organ failure and death as seen in septic patients.

Upregulation of the inhibitory receptor ILT4 in monocytes from septic patients

Sepsis-induced immune dysfunction is a complex phenomenon that involves both innate and adaptive responses. Upregulation of the inhibitor receptor named immunoglobulin like transcript 4 (ILT4) is crucial to the tolerogenic function of monocytes. Here, ILT4 expression, endotoxin-induced IL-12 and IL-10 production and CD86 expression were investigated in circulating monocytes from 16 patients with severe sepsis and 16 age and sex matched controls. We found that monocytes from patients with severe sepsis express significantly higher levels of ILT4 than monocytes from controls. Upregulation of ILT4 expression appeared to be induced by soluble factors present in the serum of septic patients and directly correlated with the degree of organ dysfunction. ILT4(+) monocytes from septic patients also displayed an alteration in the cytokine response to endotoxin stimulation characterized by reduced IL-12 production and increased IL-10 production, and a reduced expression of the costimulatory molecule CD86. In conclusion, the increased ILT4 expression and IL-10 production and the decreased CD86 expression and IL-12 production indicate that during sepsis monocytes undergo substantial modulation of the surface and cytokine phenotype. These phenotypic changes may interfere with the antigen presenting cell activity of monocytes, which may contribute to the impairment of adaptive immune responses that takes place during sepsis.