Risk of celiac disease autoimmunity is modified by interactions between CD247 and environmental exposures

Season of birth, viral infections, HLA haplogenotypes and non-HLA variants are implicated in the development of celiac disease and celiac disease autoimmunity, suggesting a combined role of genes and environmental exposures. The aim of the study was to further decipher the biological pathways conveying the season of birth effect in celiac disease autoimmunity to gain novel insights into the early pathogenesis of celiac disease. Interactions between season of birth, genetics, and early-life environmental factors on the risk of celiac autoimmunity were investigated in the multicenter TEDDY birth cohort study. Altogether 6523 genetically predisposed children were enrolled to long-term follow-up with prospective sampling and data collection at six research centers in the USA, Germany, Sweden and Finland. Celiac disease autoimmunity was defined as positive tissue transglutaminase antibodies in two consecutive serum samples. There was a significant season of birth effect on the risk of celiac autoimmunity. The effect was dependent on polymorphisms in CD247 gene encoding for CD3ζ chain of TCR-CD3 complex. In particular, children with major alleles for SNP rs864537A > G, in CD247 (AA genotype) had an excess risk of celiac autoimmunity when born March-August as compared to other months. The interaction of CD247 with season of birth on autoimmunity risk was accompanied by interactions with febrile infections between the ages of 3-6 months. Considering the important role of TCR-CD3 complex in the adaptive immune response and our findings here, CD247 variants and their possible effect of subgroups in autoimmunity development could be of interest in the design of future gene-environment studies of celiac disease. ClinicalTrials.gov Identifier: NCT00279318.

Revealing novel genomic insights and therapeutic targets for juvenile idiopathic arthritis through omics

BACKGROUND

The genetic architecture of JIA remains only partially comprehended. There is a clear imperative for continued endeavours to uncover insights into the underlying causes of JIA.

METHODS

This study encompassed a comprehensive spectrum of endeavours, including conducting a JIA genome-wide association study (GWAS) meta-analysis that incorporated data from 4550 JIA cases and 18 446 controls. We employed in silico and genome-editing approaches to prioritizing target genes. To investigate pleiotropic effects, we conducted phenome-wide association studies. Cell-type enrichment analyses were performed by integrating bulk and single-cell sequencing data. Finally, we delved into potential druggable targets for JIA.

RESULTS

Fourteen genome-wide significant non-HLA loci were identified, including four novel loci, each exhibiting pleiotropic associations with other autoimmune diseases or musculoskeletal traits. We uncovered strong genetic correlation between JIA and BMD traits at 52 genomic regions, including three GWAS loci for JIA. Candidate genes with immune functions were captured by in silico analyses at each novel locus, with additional findings identified through our experimental approach. Cell-type enrichment analysis revealed 21 specific immune cell types crucial for the affected organs in JIA, indicating their potential contribution to the disease. Finally, 24 known or candidate druggable target genes were prioritized.

CONCLUSIONS

Our identification of four novel JIA-associated genes, CD247, RHOH, COLEC10 and IRF8, broadens the novel potential drug repositioning opportunities. We established a new genetic link between COLEC10, TNFRSF11B and JIA/BMD. Additionally, the identification of RHOH underscores its role in positive thymocyte selection, thereby illuminating a critical facet of JIA's underlying biological mechanisms.

Immunogenetics of Systemic Sclerosis

Systemic sclerosis (SSc) is a rare autoimmune connective tissue disorder characterized by massive fibrosis, vascular damage, and immune imbalance. Advances in rheumatology and immunology over the past two decades have led to a redefinition of systemic sclerosis, shifting from its initial perception as primarily a "hyperfibrotic" state towards a recognition of systemic sclerosis as an immune-mediated disease. Consequently, the search for genetic markers has transitioned from focusing on fibrotic mechanisms to exploring immune regulatory pathways. Immunogenetics, an emerging field at the intersection of immunology, molecular biology, and genetics has provided valuable insights into inherited factors that influence immunity. Data from genetic studies conducted thus far indicate that alterations in genetic messages can significantly impact disease risk and progression. While certain genetic variations may confer protective effects, others may exacerbate disease susceptibility. This paper presents a comprehensive review of the most relevant genetic changes that influence both the risk and course of systemic sclerosis. Special emphasis is placed on factors regulating the immune response, recognizing their pivotal role in the pathogenesis of the disease.

Identification of biomarkers for abdominal aortic aneurysm in Behçet's disease via mendelian randomization and integrated bioinformatics analyses

Behçet's disease (BD) is a complex autoimmune disorder impacting several organ systems. Although the involvement of abdominal aortic aneurysm (AAA) in BD is rare, it can be associated with severe consequences. In the present study, we identified diagnostic biomarkers in patients with BD having AAA. Mendelian randomization (MR) analysis was initially used to explore the potential causal association between BD and AAA. The Limma package, WGCNA, PPI and machine learning algorithms were employed to identify potential diagnostic genes. A receiver operating characteristic curve (ROC) for the nomogram was constructed to ascertain the diagnostic value of AAA in patients with BD. Finally, immune cell infiltration analyses and single-sample gene set enrichment analysis (ssGSEA) were conducted. The MR analysis indicated a suggestive association between BD and the risk of AAA (odds ratio [OR]: 1.0384, 95% confidence interval [CI]: 1.0081-1.0696, p = 0.0126). Three hub genes (CD247, CD2 and CCR7) were identified using the integrated bioinformatics analyses, which were subsequently utilised to construct a nomogram (area under the curve [AUC]: 0.982, 95% CI: 0.944-1.000). Finally, the immune cell infiltration assay revealed that dysregulation immune cells were positively correlated with the three hub genes. Our MR analyses revealed a higher susceptibility of patients with BD to AAA. We used a systematic approach to identify three potential hub genes (CD247, CD2 and CCR7) and developed a nomogram to assist in the diagnosis of AAA among patients with BD. In addition, immune cell infiltration analysis indicated the dysregulation in immune cell proportions.

Screening of core genes prognostic for sepsis and construction of a ceRNA regulatory network

OBJECTIVE

To screen out core genes potentially prognostic for sepsis and construct a competing endogenous RNA (ceRNA) regulatory network.

METHODS

Subjects included in this project were 23 sepsis patients and 10 healthy people. RNA-seq for lncRNA, miRNA and mRNA was performed in the peripheral blood samples. Differentially expressed RNAs (DER) were screened out for further analysis. GO annotation and GSEA functional clustering were performed to view the functional enrichment of DEmRNAs. Core genes of prognostic significance were screened out with the weighted correlation network analysis (WGCNA). Meta-analysis and Survival analysis was devised in different microarray datasets. RT-qPCR was conducted to validate these core genes. A ceRNA network was accordingly constructed according to the correlation analysis and molecular interaction prediction.

RESULTS

RNA-seq and differential analysis screened out 1,044 DEmRNAs, 66 DEmiRNAs and 155 DElncRNAs. The GO and GSEA analysis revealed that DEmRNAs are mainly involved in inflammatory response, immune regulation, neutrophil activation. WGCNA revealed 4 potential core genes, including CD247, IL-2Rβ, TGF-βR3 and IL-1R2. In vitro cellular experiment showed up-regulated expression of IL-1R2 while down-regulated of CD247, IL-2Rβ, TGF-βR3 in sepsis patients. Correspondingly, a ceRNA regulatory network was build based on the core genes, and multiple lncRNAs and miRNAs were identified to have a potential regulatory role in sepsis.

CONCLUSION

This study identified four core genes, including CD247, IL-1R2, IL-2Rβ and TGF-βR3, with potential to be novel biomarkers for the prognosis of sepsis. In the meantime, a ceRNA network was constructed aiming to guide further study on prognostic mechanism in sepsis.

Identification of Immune-Related Hub Genes in Thymoma: Defects in CD247 and Characteristics of Paraneoplastic Syndrome

Background: Thymomas (Ts) and thymic carcinomas (TCs) are rare primary tumors of the mediastinum. Paraneoplastic syndrome (PNS) is an important feature of thymoma, which presents great challenges to clinicians.

Methods: The present study uses the weighted gene co-expression network analysis (WGCNA) to identify possible immunologic mechanisms of thymoma. RNA sequencing data from thymoma samples were downloaded from the TCGA. Core genes were taken from the module that is closely related to the WHO’s stage of classification. Enhanced analysis using the online database “Metascape” and an overall survival (OS) analysis were carried out via the Kaplan–Meier method. The hub genes were obtained from the protein–protein interaction (PPI) network. In addition, we jointly analyzed multiple sets of PNS data related to thymomas from other sources to verify the correlation between thymomas and PNS. The impact of hub genes on the prognosis of PNS was evaluated via the ROC curve, with simultaneous analysis of immune infiltration by CIBERSORT.

Findings: The 14 immune hub genes closely related to thymomas were found to be jointly involved in the T-cell receptor signaling pathway. Compared to the normal thymus and type B1/B2 thymoma, there is a lower number of T-cells in type A/B3 thymoma and thymic carcinoma. The expression of genes related to the T-cell receptor signaling pathway appeared defective. The low expression of CD247 and the decrease in the number of mature T-cells are common features among thymomas, specific pulmonary fibrosis, rheumatoid arthritis, and systemic lupus erythematosus.

Advances in CD247

CD247, which is also known as CD3ζ, CD3H, CD3Q, CD3Z, IMD25, T3Z, and TCRZ, encodes CD3ζ protein, which is expressed primarily in natural killer (NK) and T cells. Since the discovery of the ζ peptide in 1986, it has been continuously investigated. In this paper, we review the composition, molecular mechanisms and regulatory factors of CD247 expression in T cells; and review the autoimmune diseases, tumors and inflammatory diseases associated with CD247, providing a detailed and comprehensive reference for further research on the mechanism of CD247 and related diseases.

Analysis of Sepsis Markers and Pathogenesis Based on Gene Differential Expression and Protein Interaction Network

Objective

The purpose of the present study is to screen the hub genes associated with sepsis, comprehensively understand the occurrence and progress mechanism of sepsis, and provide new targets for clinical diagnosis and treatment of sepsis.

Methods

The microarray data of GSE9692 and GSE95233 were downloaded from the Gene Expression Omnibus (GEO) database. The dataset GSE9692 contained 29 children with sepsis and 16 healthy children, while the dataset GSE95233 included 102 septic subjects and 22 healthy volunteers. Differentially expressed genes (DEGs) were screened by GEO2R online analysis. The DAVID database was applied to conduct functional enrichment analysis of the DEGs. The STRING database was adopted to acquire protein-protein interaction (PPI) networks.

Results

We identified 286 DEGs (217 upregulated DEGs and 69 downregulated DEGs) in the dataset GSE9692 and 357 DEGs (236 upregulated DEGs and 121 downregulated DEGs) in the dataset GSE95233. After the intersection of DEGs of the two datasets, a total of 98 co-DEGs were obtained. DEGs associated with sepsis were involved in inflammatory responses such as T cell activation, leukocyte cell-cell adhesion, leukocyte-mediated immunity, cytokine production, immune effector process, lymphocyte-mediated immunity, defense response to fungus, and lymphocyte-mediated immunity. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis suggested that sepsis was connected to bacterial and viral infections. Through PPI network analysis, we screened the most important hub genes, including ITK, CD247, MMP9, CD3D, MMP8, KLRK1, and GZMK.

Conclusions

In conclusion, the present study revealed an unbalanced immune response at the transcriptome level of sepsis and identified genes for potential biomarkers of sepsis, such as ITK, CD247, MMP9, CD3D, MMP8, KLRK1, and GZMK.

A Genome-Wide Association Study Identifies Novel Susceptibility loci in Chronic Chagas Cardiomyopathy

BACKGROUND

Chagas disease is an infectious disease caused by the parasite Trypanosoma cruzi and is endemic from Latin American countries. The goal of our study was to identify novel genetic loci associated with chronic Chagas cardiomyopathy development in Chagas disease patients from different Latin American populations.

METHODS

We performed a cross-sectional, nested case-control study including 3 sample collections from Colombia, Argentina, and Bolivia. Samples were genotyped to conduct a genome-wide association study (GWAS). These results were meta-analyzed with summary statistic data from Brazil, gathering a total of 3413 Chagas disease patients. To identify the functional impact of the associated variant and its proxies, we performed an in silico analysis of this region.

RESULTS

The meta-analysis revealed a novel genome-wide statistically significant association with chronic Chagas cardiomyopathy development in rs2458298 (OR = 0.90, 95%CI = 0.87-0.94, P-value = 3.27 × 10-08), nearby the SAC3D1 gene. In addition, further in silico analyses displayed functional relationships between the associated variant and the SNX15, BAFT2, and FERMT3 genes, related to cardiovascular traits.

CONCLUSIONS

Our findings support the role of the host genetic factors in the susceptibility to the development of the chronic cardiac form of this neglected disease.

Alternative Splicing: A New Cause and Potential Therapeutic Target in Autoimmune Disease

Alternative splicing (AS) is a complex coordinated transcriptional regulatory mechanism. It affects nearly 95% of all protein-coding genes and occurs in nearly all human organs. Aberrant alternative splicing can lead to various neurological diseases and cancers and is responsible for aging, infection, inflammation, immune and metabolic disorders, and so on. Though aberrant alternative splicing events and their regulatory mechanisms are widely recognized, the association between autoimmune disease and alternative splicing has not been extensively examined. Autoimmune diseases are characterized by the loss of tolerance of the immune system towards self-antigens and organ-specific or systemic inflammation and subsequent tissue damage. In the present review, we summarized the most recent reports on splicing events that occur in the immunopathogenesis of systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) and attempted to clarify the role that splicing events play in regulating autoimmune disease progression. We also identified the changes that occur in splicing factor expression. The foregoing information might improve our understanding of autoimmune diseases and help develop new diagnostic and therapeutic tools for them.

Screening of key genes related to the prognosis of mouse sepsis

Sepsis is a common clinical disease with high mortality, and patients with sepsis have varied prognoses. Researchers need to explore the underlying mechanisms that determine the prognosis of sepsis. Hence, a mouse model was used to evaluate new potential prognostic markers of sepsis. Mice were randomly divided into low-dose group (n=3, lipopolysaccharides [LPS], 20 mg/kg) and high-dose group (n=3; LPS, 40 mg/kg). Total RNA was extracted from the peripheral blood of mice, and samples were then subjected to RNA sequencing. When complete data were normalized, the high-dose group and low-dose group were screened for differentially expressed genes (DEGs, log2FC ≥ 1 and q value ≤ 0.05). DEGs were analyzed by gene ontology enrichment, and potential core genes were screened using protein–protein interaction (PPI) network and weighted gene co-expression network analysis (WGCNA). Moreover, the survival data in GSE65682 were used to observe the correlation between core genes and prognosis. A total of 967 DEGs were identified in the low-dose group, of which 390 were up-regulated and 577 were down-regulated. These genes were mainly enriched in white blood cell activation, lymphocyte activation, immune system response etc. LCK, ZAP70, ITK, CD247, and DOCK2 were found at the core of PPI network, while WGCNA found that interferon-inducible protein 35 (IFI35), ITGB3, and mediator complex subunit 25 (MED25) may be potential core genes. It was demonstrated that CD247, DOCK2, IFI35, ITK, and LCK core genes were positively correlated with prognosis based on GSE65682. CD247, DOCK2, IFI35, ITK, LCK, and MED25 might be important targets affecting the prognosis of sepsis.

Marek's Disease Virus Regulates the Ubiquitylome of Chicken CD4+ T Cells to Promote Tumorigenesis

Ubiquitination and deubiquitination of cellular proteins are reciprocal reactions catalyzed by ubiquitination-related enzymes and deubiquitinase (DUB) which regulate almost all cellular processes. Marek's disease virus (MDV) encodes a viral DUB that plays an important role in the MDV pathogenicity. Chicken CD4+ T-cell lymphoma induced by MDV is a key contributor to multiple visceral tumors and immunosuppression of chickens with Marek's disease (MD). However, alterations in the ubiquitylome of MDV-induced T lymphoma cells are still unclear. In this study, a specific antibody against K-ε-GG was used to isolate ubiquitinated peptides from CD4+ T cells and MD T lymphoma cells. Mass spectrometry was used to compare and analyze alterations in the ubiquitylome. Our results showed that the ubiquitination of 717 and 778 proteins was significantly up- and downregulated, respectively, in T lymphoma cells. MDV up- and downregulated ubiquitination of a similar percentage of proteins. The ubiquitination of transferases, especially serine/threonine kinases, was the main regulatory target of MDV. Compared with CD4+ T cells of the control group, MDV mainly altered the ubiquitylome associated with the signal transduction, immune system, cancer, and infectious disease pathways in T lymphoma cells. In these pathways, the ubiquitination of CDK1, IL-18, PRKCB, ETV6, and EST1 proteins was significantly up- or downregulated as shown by immunoblotting. The current study revealed that the MDV infection could exert a significant influence on the ubiquitylome of CD4+ T cells.

Association of CD247 (CD3ζ) gene polymorphisms with T1D and AITD in the population of northern Sweden

Background

T1D and AITD are autoimmune disorders commonly occurring in the same family and even in the same individual. The genetic contribution to these disorders is complex making uncovering of susceptibility genes very challenging. The general aim of this study was to identify loci and genes contributing to T1D/AITD susceptibility. Our strategy was to perform linkage and association studies in the relatively genetically homogenous population of northern Sweden. We performed a GWLS to find genomic regions linked to T1D/AITD in families from northern Sweden and we performed an association study in the families to test for association between T1D/AITD and variants in previously published candidate genes as well as a novel candidate gene, CD247.

Methods

DNA prepared from 459 individuals was used to perform a linkage and an association study. The ABI PRISM Linkage Mapping Set v2.5MD10 was employed for an initial 10-cM GWLS, and additional markers were added for fine mapping. Merlin was used for linkage calculations. For the association analysis, a GoldenGate Custom Panel from Illumina containing 79 SNPs of interest was used and FBAT was used for association calculations.

Results

Our study revealed linkage to two previously identified chromosomal regions, 4q25 and 6p22, as well as to a novel chromosomal region, 1q23. The association study replicated association to PTPN22, HLA-DRB1, INS, IFIH1, CTLA4 and C12orf30. Evidence in favor of association was also found for SNPs in the novel susceptibility gene CD247.

Conclusions

Several risk loci for T1D/AITD identified in published association studies were replicated in a family material, of modest size, from northern Sweden. This provides evidence that these loci confer disease susceptibility in this population and emphasizes that small to intermediate sized family studies in this population can be used in a cost-effective manner for the search of genes involved in complex diseases. The linkage study revealed a chromosomal region in which a novel T1D/AITD susceptibility gene, CD247, is located. The association study showed association between T1D/AITD and several variants in this gene. These results suggests that common susceptibility genes act in concert with variants of CD247 to generate genetic risk for T1D/AITD in this population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-016-0333-z) contains supplementary material, which is available to authorized users.

A validated gene regulatory network and GWAS identifies early regulators of T cell-associated diseases

Early regulators of disease may increase understanding of disease mechanisms and serve as markers for presymptomatic diagnosis and treatment. However, early regulators are difficult to identify because patients generally present after they are symptomatic. We hypothesized that early regulators of T cell-associated diseases could be found by identifying upstream transcription factors (TFs) in T cell differentiation and by prioritizing hub TFs that were enriched for disease-associated polymorphisms. A gene regulatory network (GRN) was constructed by time series profiling of the transcriptomes and methylomes of human CD4(+) T cells during in vitro differentiation into four helper T cell lineages, in combination with sequence-based TF binding predictions. The TFs GATA3, MAF, and MYB were identified as early regulators and validated by ChIP-seq (chromatin immunoprecipitation sequencing) and small interfering RNA knockdowns. Differential mRNA expression of the TFs and their targets in T cell-associated diseases supports their clinical relevance. To directly test if the TFs were altered early in disease, T cells from patients with two T cell-mediated diseases, multiple sclerosis and seasonal allergic rhinitis, were analyzed. Strikingly, the TFs were differentially expressed during asymptomatic stages of both diseases, whereas their targets showed altered expression during symptomatic stages. This analytical strategy to identify early regulators of disease by combining GRNs with genome-wide association studies may be generally applicable for functional and clinical studies of early disease development.

Brief Exercises Affect Gene Expression in Circulating Monocytes

We aimed to give a systematic hypothesis on the functions of exercise on circulating monocytes by identifying a discrete set of genes in circulating monocytes that were altered by exercise. The microarray expression profile of GSE51835 was downloaded from gene expression omnibus (GEO) database for the identification of differentially expressed genes (DEGs) using limma and affy packages in R language. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for DEGs, followed by the construction of co-expression network and protein-protein interaction (PPI) network. The top 10 nodes in PPI network were screened, and subnetwork was constructed for the key genes identification. Totally, 35 DEGs, including 2 upregulated genes and 33 downregulated genes, were identified. The enriched GO terms were mainly linked to immune response and defence response, and the enriched KEGG pathways were mainly associated with natural killer cell-mediated cytotoxicity and graft-versus-host disease. Dual-specificity phosphatase 2 (DUSP2) was identified as a key node in the co-expression network. In the PPI network, CD247 module (CD247), chemokine (C-X-C motif) receptor 4 (CXCR4), granzyme B (GZMB) and perforin 1 (PRF1) were identified as key nodes. An important interaction, GZMB/PRF1, was detected. Five key genes, including DUSP2, CD247, CXCR4, GZMB and PRF1, and an interaction of GZMB/PRF1, were significant factors in the immune processes of circulating monocytes, which might be regulated by brief exercises, leading to the enhancement of immune function.

Genetics and epigenetics of systemic lupus erythematosus

Genetics unquestionably contributes to systemic lupus erythematosus (SLE) predisposition, progression and outcome. Nevertheless, single-gene defects causing lupus-like phenotypes have been infrequently documented. The majority of the identified genetic SLE risk factors are, therefore, common variants, responsible for a small effect on the global risk. Recently, genome wide association studies led to the identification of a growing number of gene variants associated with SLE susceptibility, particular disease phenotypes, and antibody profiles. Further studies addressed the biological effects of these variants. In addition, the role of epigenetics has recently been revealed. These combined efforts contributed to a better understanding of SLE pathogenesis and to the characterization of clinically relevant pathways. In this review, we describe SLE-associated single-gene defects, common variants, and epigenetic changes. We also discuss the limitations of current methods and the challenges that we still have to face in order to incorporate genomic and epigenomic data into clinical practice.