Screening and identification of key biomarkers associated with endometriosis using bioinformatics and next-generation sequencing data analysis

Background

Endometriosis is a common cause of endometrial-type mucosa outside the uterine cavity with symptoms such as painful periods, chronic pelvic pain, pain with intercourse and infertility. However, the early diagnosis of endometriosis is still restricted. The purpose of this investigation is to identify and validate the key biomarkers of endometriosis.

Methods

Next-generation sequencing dataset GSE243039 was obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) between endometriosis and normal control samples were identified. After screening of DEGs, gene ontology (GO) and REACTOME pathway enrichment analyses were performed. Furthermore, a protein–protein interaction (PPI) network was constructed and modules were analyzed using the Human Integrated Protein–Protein Interaction rEference database and Cytoscape software, and hub genes were identified. Subsequently, a network between miRNAs and hub genes, and network between TFs and hub genes were constructed using the miRNet and NetworkAnalyst tool, and possible key miRNAs and TFs were predicted. Finally, receiver operating characteristic curve analysis was used to validate the hub genes.

Results

A total of 958 DEGs, including 479 upregulated genes and 479 downregulated genes, were screened between endometriosis and normal control samples. GO and REACTOME pathway enrichment analyses of the 958 DEGs showed that they were mainly involved in multicellular organismal process, developmental process, signaling by GPCR and muscle contraction. Further analysis of the PPI network and modules identified 10 hub genes, including vcam1, snca, prkcb, adrb2, foxq1, mdfi, actbl2, prkd1, dapk1 and actc1. Possible target miRNAs, including hsa-mir-3143 and hsa-mir-2110, and target TFs, including tcf3 (transcription factor 3) and clock (clock circadian regulator), were predicted by constructing a miRNA-hub gene regulatory network and TF-hub gene regulatory network.

Conclusions

This investigation used bioinformatics techniques to explore the potential and novel biomarkers. These biomarkers might provide new ideas and methods for the early diagnosis, treatment and monitoring of endometriosis.

Fifty years of HLA-associated type 1 diabetes risk: history, current knowledge, and future directions

More than 50 years have elapsed since the association of human leukocyte antigens (HLA) with type 1 diabetes (T1D) was first reported. Since then, methods for identification of HLA have progressed from cell based to DNA based, and the number of recognized HLA variants has grown from a few to tens of thousands. Current genotyping methodology allows for exact identification of all HLA-encoding genes in an individual's genome, with statistical analysis methods evolving to digest the enormous amount of data that can be produced at an astonishing rate. The HLA region of the genome has been repeatedly shown to be the most important genetic risk factor for T1D, and the original reported associations have been replicated, refined, and expanded. Even with the remarkable progress through 50 years and over 5,000 reports, a comprehensive understanding of all effects of HLA on T1D remains elusive. This report represents a summary of the field as it evolved and as it stands now, enumerating many past and present challenges, and suggests possible paradigm shifts for moving forward with future studies in hopes of finally understanding all the ways in which HLA influences the pathophysiology of T1D.

Complete HLA genotyping of type 1 diabetes patients and controls from Mali reveals both expected and novel disease associations

HLA genotyping was performed on 99 type 1 diabetes (T1D) patients and 200 controls from Mali. Next-generation sequencing of the classical HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQA1, -DQB1, -DPA1, and -DPB1 loci revealed strong T1D association for all loci except HLA-C and -DPA1. Class II association is stronger than class I association, with most observed associations predisposing or protective as expected based on previous studies. For example, HLA-DRB1*03:01, HLA-DRB1*09:01, and HLA-DRB1*04:05 predispose for T1D, whereas HLA-DRB1*15:03 is protective. HLA-DPB1*04:02 (OR = 12.73, p = 2.92 × 10-05 ) and HLA-B*27:05 (OR = 21.36, p = 3.72 × 10-05 ) appear highly predisposing, although previous studies involving multiple populations have reported HLA-DPB1*04:02 as T1D-protective and HLA-B*27:05 as neutral. This result may reflect the linkage disequilibrium between alleles on the extended HLA-A*24:02~HLA-B*27:05~HLA-C*02:02~HLA-DRB1*04:05~HLA-DRB4*01:03~HLA-DQB1*02:02~HLA-DQA1*02:01~HLA-DPB1*04:02~HLA-DPA1*01:03 haplotype in this population rather than an effect of either allele itself. Individual amino acid (AA) analyses are consistent with most T1D association attributable to HLA class II rather than class I in this data set. AA-level analyses reveal previously undescribed differences of the HLA-C locus from the HLA-A and HLA-B loci, with more polymorphic positions, spanning a larger portion of the gene. This may reflect additional mechanisms for HLA-C to influence T1D risk, for example, through expression differences or through its role as the dominant ligand for killer cell immunoglobulin-like receptors (KIR). Comparison of these data to those from larger studies and on other populations may facilitate T1D prediction and help elucidate elusive mechanisms of how HLA contributes to T1D risk and autoimmunity.

New Evidence in the Pathogenesis of Celiac Disease and Type 1 Diabetes Mellitus: A Systematic Review

Celiac disease (CD) and type 1 diabetes mellitus (T1DM) are autoimmune diseases that coexist frequently. These illnesses share a common genetic background. This study aims to review the different pathophysiologic mechanisms that have been studied about the coexistence of CD and T1DM, to contrast them, and to summarize their specific role in these autoimmune diseases. We conducted a systematic review following the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) checklist and used the Medical Subject Headings (MeSH) search strategy to obtain relevant articles. We found 585 papers which were reduced to 355 after removing duplicates. Later, the filters and inclusion/exclusion criteria were applied which ended the search with 78 articles. Finally, we reviewed the articles that contained information about the pathogenesis of CD and T1DM, their coexistence, and how the pathogenesis impacts clinical outcomes. The reviewed studies strongly conclude that the presence of human leukocyte antigen (HLA) genes DQ2 and DQ8 are high-risk for developing the coexistence of CD and T1DM. We found that killer immunoglobulin-like receptor (KIR) genes, enterovirus infection in gut cells, and gut microbiota dysbiosis with the predominance of Bacteroides spp. also play a role in the pathogenesis and development of symptoms of CD in patients with the previous diagnosis of T1DM. CD4+ and CD8+ cell levels vary among patients and studies, consequently, more study on this topic is needed.

Association of KIR gene polymorphisms with Type 1 Diabetes: a meta-analysis

Purpose

Type 1 Diabetes (T1D) is a T cell-mediated disease, in which autoimmune destruction of insulin-producing β-cells in pancreatic islets occurs. In recent decades, the role of Killer cell immunoglobulin-like receptor (KIR) gene polymorphisms in susceptibility to T1D has been demonstrated in an increased number of studies. Nonetheless, inconsistency has been observed in the results of performed association studies. To comprehensively clarify the association among KIR gene polymorphisms and the risk of T1D, this meta-analysis on the previously published association studies was carried out via incorporating multiple research.

Methods

No publication has been recorded from Nov 2017 until July 2020 about the KIR genes and T1D. The PubMed/MEDLINE and Scopus databases were systematically searched up to November 2017 to identify investigations on the impact of the polymorphisms of KIR genes on susceptibility to T1D. The odds ratio (OR) with a 95% confidence interval (95% CI) was calculated. Funnel plot and Egger test were used to assess the publication bias. Thirteen qualified published case-control articles were found for evaluation in this meta-analysis.

Results

Our results show statistical significance between the genetic variations in the KIR2DL1 (OR = 0.42, 95% CI = 0.23-0.77; P = 0.005), KIR2DL2 (OR = 1.15, 95% CI = 1.00-1.32; P = 0.048), and KIR2DL5 (OR = 0.86, 95% CI = 0.75-0.98; P = 0.03) with susceptibility to T1D.

Conclusions

This meta-analysis study provides reliable evidence that KIR gene polymorphisms may contribute to T1D risk. KIR 2DL1 and 2DL5 genes might be considered as a protective factor for T1D, while 2DL2 seemed to be a susceptibility factor.

Epistatic interactions between killer immunoglobulin-like receptors and human leukocyte antigen ligands are associated with ankylosing spondylitis

The killer immunoglobulin-like receptors (KIRs), found predominantly on the surface of natural killer (NK) cells and some T-cells, are a collection of highly polymorphic activating and inhibitory receptors with variable specificity for class I human leukocyte antigen (HLA) ligands. Fifteen KIR genes are inherited in haplotypes of diverse gene content across the human population, and the repertoire of independently inherited KIR and HLA alleles is known to alter risk for immune-mediated and infectious disease by shifting the threshold of lymphocyte activation. We have conducted the largest disease-association study of KIR-HLA epistasis to date, enabled by the imputation of KIR gene and HLA allele dosages from genotype data for 12,214 healthy controls and 8,107 individuals with the HLA-B*27-associated immune-mediated arthritis, ankylosing spondylitis (AS). We identified epistatic interactions between KIR genes and their ligands (at both HLA subtype and allele resolution) that increase risk of disease, replicating analyses in a semi-independent cohort of 3,497 cases and 14,844 controls. We further confirmed that the strong AS-association with a pathogenic variant in the endoplasmic reticulum aminopeptidase gene ERAP1, known to alter the HLA-B*27 presented peptidome, is not modified by carriage of the canonical HLA-B receptor KIR3DL1/S1. Overall, our data suggests that AS risk is modified by the complement of KIRs and HLA ligands inherited, beyond the influence of HLA-B*27 alone, which collectively alter the proinflammatory capacity of KIR-expressing lymphocytes to contribute to disease immunopathogenesis.

Cells of the immune system utilise various cell-surface receptors to differentiate between healthy and infected or malignant cells, enabling targeted inflammatory responses while minimising damage to self-tissue. In instances where the immune system fails to correctly differentiate healthy from diseased tissue, or inflammatory activity is poorly regulated, autoimmune or autoinflammatory conditions can develop. Here we have investigated a possible role for a class of immune-cell activating and inhibitory receptors in the pathogenesis of ankylosing spondylitis (AS), a common but poorly understood inflammatory arthritis in which the immune system causes severe damage to the joints of the pelvis and spine. Using genetic information from 12,214 healthy controls and 8,107 individuals with AS we were able to identify combinations of independently inherited immune cell receptors and their ligands that increase or decrease an individual’s risk of disease. This research provides new insight into the nature of co-inherited genetic factors that may collectively alter the proinflammatory capacity of immune cells, contributing to the immunopathogenesis of immune-mediated diseases.

Incidence of Type 1 Diabetes Is Not Increasing in a Population-Based Cohort in Olmsted County, Minnesota, USA

OBJECTIVES

To investigate the recent incidence of T1D in a US Midwestern county to determine whether this increase has been sustained and compare it with the incidence of celiac disease (CD) and also investigate the prevalence of CD, an associated autoimmune disease, within the cohort.

PATIENTS AND METHODS

A broad search strategy was used to identify all incident cases of T1D in Olmsted County, Minnesota, between January 1, 1994, and December 31, 2010, using the Rochester Epidemiology Project. Diagnosis and residency status were confirmed through the medical record. Incidence rates were directly standardized to the 2010 US population. Poisson regression was used to test for a change in incidence rate. Clinical charts were reviewed to confirm case status.

RESULTS

There were 233 incident cases of T1D. Directly adjusting for age and sex with respect to the 2010 US white population, the overall annual incidence of T1D was 9.2 (95% CI, 8.0-10.4) per 100,000 people per year among all ages and 19.9 (95% CI, 16.6-23.2) per 100,000 people per year for those younger than 20 years. There was no significant increase in the incidence of T1D over time (P=.45). Despite the overall stability in annual incidence, there was an initial increasing trend followed by a plateau. Of the 109 patients with T1D (47%) tested for CD, 12% (13) had biopsy-proven CD.

CONCLUSION

The incidence of T1D has stopped increasing in Olmsted County, Minnesota, in the most recent decade. Further studies are needed to confirm this finding and explore reasons for this plateau.

High-risk genotypes HLA-DR3-DQ2/DR3-DQ2 and DR3-DQ2/DR4-DQ8 in co-occurrence of type 1 diabetes and celiac disease

Shared susceptibility alleles in the HLA region contribute to the co-existence of type 1 diabetes (T1D) and celiac disease (CD). The aim of our study was to identify HLA genotype variations that influence co-occurrence of T1D and CD (T1D + CD) and the order of their onset. Totally 244 patients, 67 with T1D, 68 with CD and 69 with T1D + CD, (split into "T1D first" and "CD first"), were analyzed. Control group consisted of 130 healthy unrelated individuals. Two-tailed Fisher's exact test was used for statistical analysis. The genetic background of Slovenian CD patients resembled more northern than southern European populations with DR3-DQ2/DR3-DQ2 (odds ratio [OR] = 19.68) conferring the highest risk. The T1D + CD was associated with DR3-DQ2/DR3-DQ2 (OR = 45.53) and even more with DR3-DQ2/DR4-DQ8 (OR = 93.76). DR3-DQ2/DR7-DQ2 played a neutral role in susceptibility for T1D + CD. The order of the onset of T1D or CD in patients with co-occurring diseases was not influenced by HLA risk genotype profile. DR3-DQ2/DR3-DQ2 was associated with an increased risk for developing CD in patients with T1D, whereas patients with CD carrying DR3-DQ2/DR4-DQ8 were at higher risk for developing T1D. In addition to other genetic factors including HLA class I alleles present on DR3-DQ2 extended haplotype, the second extended haplotype may moderate the risk for T1D + CD conferred by DR3-DQ2. Our results suggested that individuals carrying high-risk genotypes DR3-DQ2/DR3-DQ2 or DR3-DQ2/DR4-DQ8 would more likely develop both T1D and CD than either disease alone.

Significant association of the KIR2DL3/HLA-C1 genotype with susceptibility to Crohn's disease

We aimed to analyze the possible association of KIR/HLA-C genotypes with the susceptibility to Crohn's disease (CD) in a Spanish population. A total of 125 patients with CD and 339 healthy controls were selected for this study. KIR and HLA-C typing were developed by sequence-specific oligonucleotide probing. We found that the centromeric A/A genotype and HLA-C1 combination was significantly increased in CD patients (P<10(-3)). The KIR2DL3/2DL3 genotype was significantly increased in CD patients (P<0.0005). Moreover, we also observed a highly significant increase of KIR2DL3-HLA-C1 homozygosis in CD patients (P<0.0005). Our results confirm the relevance of the KIR2DL2/KIR2DL3 genes and their interaction with HLA-C to CD. We show that the contribution of the KIR genes to CD susceptibility extends beyond the association with individual KIRs, with an imbalance between activating and inhibitory KIR genes seeming to influence the susceptibility to CD.

Contribution of KIR genes, HLA class I ligands, and KIR/HLA class I ligand combinations on the genetic predisposition to celiac disease and coexisting celiac disease and type 1 diabetes mellitus

BACKGOUND AND AIM

There are some common genetic features between celiac disease (CD) and diabetes mellitus type 1 (DM). However, the genetic risk factors have not been fully clarified for CD and the co-occurrence of CD and DM. KIR (killer immunoglobulin-like receptor) genes regulate the cytolitic activity of NK-cells and T lymphocytes. The aim of this study is to evaluate the contribution of KIR genes, KIR ligands, and combinations of KIR/ KIR ligands on the genetic predisposition to CD and co-occurrence of CD and DM.

MATERIAL AND METHODS

Forty six patients with CD (n = 46), 20 patients with CD+DM (n = 20), and 60 healthy controls (n = 60) were included in this study. KIR genes and KIR ligands were investigated with PCR-SSOP and PCR-SSP in all subjects, respectively.

RESULTS

This study showed that while the telomeric KIR genes (2DS5 and 3DS1), and combinations of 3DS1+HLA-BBw4-Thrand 3DS1+HLA-BBw4-Iso- (p < 0.001, p < 0.001, p < 0.001, and p < 0.001, respectively) were observed more frequently in patients with CD than in controls, the 2DS5, 3DS1 KIR genes, C1 ligand, and combinations of 3DS1+HLA-BBw4-Thr- and 3DS1+HLA-BBw4-Iso- (p = 0.002, p = 0.004, p = 0.036, p < 0.001, and p = 0.007, respectively) were observed more frequently in patients with CD+DM than in controls.

CONCLUSIONS

The results of this study indicated that some KIR genes, KIR ligands, and KIR/KIR ligand interactions may be responsible for a predisposition to CD and the coexistence of CD and DM. For development of coexisting CD and DM, the 2DS5 and 3DS1 genes, C1 ligand, and combinations of 3DS1+HLA-BBw4-Thr- and 3DS1+HLA-BBw4-Iso- were found to be risk factors.

KIR and HLA under pressure: evidences of coevolution across worldwide populations

KIR (killer cell immunoglobulin-like receptors) and HLA (human leukocyte antigens) are two distinct gene families with remarkable importance for human immune responses. The recognition of HLA molecules by activating and inhibitory KIR promotes a balance of signals that regulates NK cell function and is especially important for the innate defense against pathogens and early placentation. There is no documented gametic association between these two gene families and no evidence of common regulation. However, due to the critical role of KIR recognition for immunity and reproduction, the possibility of KIR-HLA combinations being under selective pressure is not surprising. In this manuscript, we first summarize the HLA-KIR system, the HLA molecules that are the putative ligands for KIR, and then we review the evidences that suggest these two gene families are coevolving as an integrated system.

Fifty-Five Years of Pediatric Endocrinology and 50 Years of the Department of Pediatric Endocrinology, Diabetes and Metabolic Diseases in Slovenia

Paediatric endocrinology started its independent development early in the general development of this specialty, with a strong focus on research and clinical excellence. Slovenian paediatric endocrinology was an integral part of the European paediatric endocrinology from its beginnings and a founding member of the first ‘International Study Group for Diabetes in Children and Adolescents’. After the pioneering work of Prof. Lev Matajc, Prof. Ciril Kržišnik firmly integrated the Department of Pediatric Endocrinology, Diabetes and Metabolic Diseases at the University Children’s Hospital in Ljubljana in the international scientific community. In the last decade, the department participates in cutting-edge research and provides clinical services at highest international standards.

Type 1 Diabetes in the Young: Organization of Two National Centers in Israel and Slovenia

Type 1 diabetes is a chronic autoimmune disease that affects mainly young people. In the last 50 years, a steady increase of the T1D incidence in the young is reported worldwide, with an average 4 % increase annually. In addition, the mean age at the diagnosis is decreasing. Studies show that good metabolic control is important not only for delaying the chronic complications of diabetes but also for improving the quality of life of patients and their families. Continuous education, together with modern technology, is crucial in achieving these goals. Longitudinal data on glycated hemoglobin (HbA1c), along with the data on severe hypoglycemia and severe diabetic ketoacidosis, can describe the quality of care in a defined population. Two national reference diabetes centres taking care of children, adolescents and young adults with diabetes in Israel and Slovenia are described.